CN113361831B - Non-invasive load identification electric quantity decomposition method and system based on priority distribution - Google Patents

Non-invasive load identification electric quantity decomposition method and system based on priority distribution Download PDF

Info

Publication number
CN113361831B
CN113361831B CN202110911018.7A CN202110911018A CN113361831B CN 113361831 B CN113361831 B CN 113361831B CN 202110911018 A CN202110911018 A CN 202110911018A CN 113361831 B CN113361831 B CN 113361831B
Authority
CN
China
Prior art keywords
load
priority
edge
increment
electric quantity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110911018.7A
Other languages
Chinese (zh)
Other versions
CN113361831A (en
Inventor
黄奇峰
方凯杰
左强
杨世海
黄艺璇
刘恬畅
程含渺
陈铭明
李波
陆婋泉
曹晓冬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Jiangsu Electric Power Co ltd Marketing Service Center
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Jiangsu Electric Power Co ltd Marketing Service Center
State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Jiangsu Electric Power Co ltd Marketing Service Center, State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Jiangsu Electric Power Co ltd Marketing Service Center
Priority to CN202110911018.7A priority Critical patent/CN113361831B/en
Publication of CN113361831A publication Critical patent/CN113361831A/en
Application granted granted Critical
Publication of CN113361831B publication Critical patent/CN113361831B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2310/00The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/70Load identification

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Power Engineering (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The non-invasive load identification electric quantity decomposition method and system based on priority distribution comprise load start and stop edge detection, load start and stop edge classification, first priority load electric quantity decomposition and second priority load electric quantity decomposition; according to the high-frequency sampling information collected and calculated by the non-invasive load identification equipment, detecting the start and stop edges of the load by adopting a maximum likelihood ratio algorithm, and acquiring the characteristic information of active power, reactive power and harmonic waves at the edges; then, calculating Euclidean distance of edge characteristic data, and classifying load edges so as to realize load type identification; finally, performing coulombic decomposition according to the two-step priority. The invention executes different electric quantity decomposition strategies based on the priority sequence, solves the problem of inaccurate electric quantity identification of power fluctuation characteristics of the air conditioner and the like, and can realize accurate electric quantity decomposition of a plurality of electric appliances in a simultaneous working state.

Description

Non-invasive load identification electric quantity decomposition method and system based on priority distribution
Technical Field
The invention relates to the technical field of non-intrusive load identification, in particular to a non-intrusive load identification electric quantity decomposition method and system based on priority distribution.
Background
For a long time, the information sensing capability of the equivalent equipment of the intelligent electric energy meter at the user side is weak, accurate fine-grained load data and analysis service thereof cannot be provided, and the bidirectional information interaction between an electric power company and residential users is difficult to support.
The non-invasive load identification technology is a new technology for monitoring the electricity consumption state of residents in real time, and the high-frequency sampling data of equipment such as a distributed terminal and an intelligent electric energy meter are utilized, and signal analysis and a machine learning algorithm are applied, so that the electricity consumption information of the residents, such as the start-stop time, the energy consumption level, the use rule and the like of an electric appliance, can be acquired more finely. The technology can realize the accurate identification of the composition and energy consumption of residential electric equipment in all time periods, thereby providing value-added services such as energy-saving diagnosis, identification of potential safety hazards of electricity consumption, service life prediction of electric appliances, utilization rate statistics of electric appliances and the like for customers.
Under the background, the non-invasive load identification technology is widely concerned by researchers at home and abroad, and is always in high research heat. In the prior art, a non-invasive load identification process includes data acquisition, data preprocessing, event detection, feature extraction, load identification and other steps. The data preprocessing comprises the substeps of abnormal value elimination, power value calculation and the like, and the load identification comprises the substeps of load identification, power decomposition and the like. The research literature and patents of the load identification technology direction are numerous, and the Chinese patent application (CN 112039059A) 'a long transient load event detection method based on power step continuity judgment', the Chinese patent application (CN 111191908A) 'a load event detection method, a system and a storage medium' and the Chinese patent application (CN 104483575A) 'a self-adaptive load event detection method for non-invasive power monitoring' mainly concentrate on load event detection; the Chinese patent application (CN 112327046A) is a non-invasive load monitoring method based on fuzzy clustering and support vector regression, the Chinese patent application (CN 111949707A) is a non-invasive load decomposition method based on a hidden Markov model of an image field, the Chinese patent application (CN 111830347A) is a two-stage non-invasive load monitoring method based on events, and the Chinese patent application (CN 111553444A) is a load identification method based on non-invasive load terminal data, and the load type identification is concerned; in the terminal device and system direction, the representative technology is chinese utility model patent (CN 210781189U) "hardware device for collecting user electric meter load identification data", chinese utility model patent (CN 208297577U) "a user load identification module" capable of being embedded in an electric meter ", and chinese utility model patent (CN 204992779U)" a non-invasive electric load identification system ". However, there are few studies on non-invasive load electricity decomposition, and in the prior art, the chinese patent application (CN 111753968A) "non-invasive load monitoring intelligent ammeter and electricity decomposition method" adopts a sequence-to-point electricity decomposition model with an attention mechanism added, and performs fine total electricity decomposition by self-learning weights, whereas the chinese patent application (CN 109840691A) constructs a itemized electricity decomposition estimation method based on a deep neural network.
Whether the deep neural network model or the sequence-to-point electricity decomposition model is adopted, a large amount of training sample data required by the model is difficult to obtain in practical application. The difference of the types of the electric appliances among different families is huge, the situations of updating iteration and the like of the electric appliances of the same type cause the quantity of training samples to be increased sharply, otherwise, the identification accuracy is reduced. Meanwhile, the model needs to go through a long training and learning process, is time-consuming and is not suitable for a localized load identification device for real-time acquisition and real-time processing and analysis. Therefore, it is desirable to provide a non-intrusive load identification power decomposition method for distributed devices.
Disclosure of Invention
In order to solve the defects in the prior art, the invention aims to provide a non-invasive load identification electric quantity decomposition method and a non-invasive load identification electric quantity decomposition system based on priority distribution, which are used for performing decomposition calculation on the electric quantity of a load to be identified based on a multi-priority distribution method on the premise of completing quick detection and accurate matching of a load edge, so that the electric quantity decomposition precision in the non-invasive load identification technology is improved, and the influence of fluctuating load on electric quantity decomposition is avoided.
The invention adopts the following technical scheme.
The non-intrusive load identification electric quantity decomposition method based on priority distribution comprises the following steps:
step 1, carrying out high-frequency sampling on current and voltage of an in-house bus, and carrying out fast Fourier transform on high-frequency sampling values of the current and the voltage so as to obtain active power, reactive power, direct current component and harmonic current component of a power frequency period;
step 2, detecting a load edge by using the active power of a power frequency cycle and adopting a maximum likelihood ratio algorithm, and synchronously calculating the active power increment, the reactive power increment, the increment of a direct current component and the increment of a harmonic current component of the load edge;
step 3, extracting characteristic vectors from the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge; identifying the type of the load edge by using the characteristic vector and each electric appliance characteristic vector in the electric appliance characteristic library through an Euclidean distance judgment method;
step 4, according to different types of load edges, performing multi-priority distribution on various loads to obtain a first priority load and a second priority load; the first priority load is a stability load, and the second priority load is a volatility load;
step 5, for any load edge, matching the load shutdown edge with the load startup edge to form a load identification pair; and under any load identification pair, carrying out electric quantity decomposition on the first priority load to obtain the consumed electric quantity of the corresponding electric appliance, simultaneously obtaining the total electric quantity of the second priority load from starting to stopping, and subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
Preferably, in step 1, the current and voltage of the household bus are sampled at high frequency by using an electric energy meter or other terminal equipment, and the high-frequency sampled current and voltage are output by the metering core; wherein, the sampling frequency of the high-frequency sampling period is 1.6 kHz.
The step 1 comprises the following steps:
step 1.1, collecting current and voltage of the household bus based on a high-frequency sampling period, and respectively forming a current sampling sequence by utilizing high-frequency sampling values of the current and the voltage
Figure 100002_DEST_PATH_IMAGE001
And a voltage sampling sequence
Figure 100002_DEST_PATH_IMAGE002
Step 1.2, sampling sequence of current according to the following relation
Figure 526623DEST_PATH_IMAGE001
And a voltage sampling sequence
Figure 974922DEST_PATH_IMAGE002
Performing fast Fourier transform to obtain current amplitude result sequence
Figure 100002_DEST_PATH_IMAGE003
Current phase result sequence
Figure 100002_DEST_PATH_IMAGE004
Voltage amplitude result sequence
Figure 100002_DEST_PATH_IMAGE005
And a sequence of voltage phase results
Figure 100002_DEST_PATH_IMAGE006
Figure 100002_DEST_PATH_IMAGE007
In the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE008
is a fast fourier transform of a sequence of current samples,
Figure 100002_DEST_PATH_IMAGE009
fast fourier transform of a sequence of voltage samples;
step 1.3, calculating the active power, the reactive power, the direct current component and the second component of the power frequency period according to the following relational expressionskSubharmonic current component:
Figure 100002_DEST_PATH_IMAGE010
Figure 100002_DEST_PATH_IMAGE011
in the formula (I), the compound is shown in the specification,kin order to be the harmonic frequency, the frequency of the harmonic wave,Pis the active power of the power frequency period,QIs the reactive power of the power frequency period,
Figure 100002_DEST_PATH_IMAGE012
Is a direct current of a power frequency period,
Figure 100002_DEST_PATH_IMAGE013
is the first of power frequency cyclekThe effective value of the sub-harmonic current,
Figure 100002_DEST_PATH_IMAGE014
is as followskThe magnitude of the current at the sub-harmonic,
Figure 100002_DEST_PATH_IMAGE015
is as followskThe voltage amplitude of the sub-harmonics,
Figure 100002_DEST_PATH_IMAGE016
is as followskThe phase of the current of the sub-harmonic,
Figure 100002_DEST_PATH_IMAGE017
is as followskThe phase of the voltage of the sub-harmonic,
Figure 100002_DEST_PATH_IMAGE018
is the current amplitude of the fundamental wave.
Preferably, step 2 comprises:
step 2.1, arranging active power, reactive power, direct current component and harmonic current component of a power frequency period into an active power sequence according to a period time sequence
Figure 100002_DEST_PATH_IMAGE019
Reactive power sequence
Figure 100002_DEST_PATH_IMAGE020
Sequence of direct current components
Figure 100002_DEST_PATH_IMAGE021
And a firstkSequence of sub-harmonic current components
Figure 100002_DEST_PATH_IMAGE022
Step 2.2, for the mth period, constructing a sliding data window for the active power sequence; wherein the sliding data window comprises a front sliding data window
Figure 100002_DEST_PATH_IMAGE023
And a rear sliding data window
Figure 100002_DEST_PATH_IMAGE024
(ii) a Wherein, the front sliding data window is an active power sequence formed by 50 active powers before the m period, and the front sliding data window is formed by 50 active powers after the m periodA resulting active power sequence;
step 2.3, calculating the maximum likelihood ratio of the sliding data window at the mth period according to the following relational expression
Figure 100002_DEST_PATH_IMAGE025
Figure 100002_DEST_PATH_IMAGE026
In the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE027
is a function of the probability of the normal distribution,
Figure 100002_DEST_PATH_IMAGE028
for the mean value of the post-sliding data window,
Figure 100002_DEST_PATH_IMAGE029
for the post-sliding data window standard deviation,
Figure 100002_DEST_PATH_IMAGE030
is the mean of the front sliding window of data,
Figure 100002_DEST_PATH_IMAGE031
a front sliding data window standard deviation;
step 2.4, if the maximum likelihood ratio of the sliding data window at the mth period is larger than the set event threshold, judging that a load edge event exists at the mth period
Figure 100002_DEST_PATH_IMAGE032
(ii) a Wherein the set event threshold value is 0.5;
step 2.5, for the presence of load edge events
Figure 800664DEST_PATH_IMAGE032
Period m of (a), if
Figure 100002_DEST_PATH_IMAGE033
Then the load edge event is judged to be a load starting edge event
Figure 100002_DEST_PATH_IMAGE034
(ii) a On the contrary, if
Figure 100002_DEST_PATH_IMAGE035
Then the load edge event is determined to be a load shutdown edge event
Figure 100002_DEST_PATH_IMAGE036
Step 2.6, respectively calculating the active power increment of the load edge occurring at the mth period according to the following relational expression
Figure 100002_DEST_PATH_IMAGE037
Reactive power increment
Figure 100002_DEST_PATH_IMAGE038
Delta of DC current component
Figure 100002_DEST_PATH_IMAGE039
And a firstkIncrement of subharmonic current component
Figure 100002_DEST_PATH_IMAGE040
Figure 100002_DEST_PATH_IMAGE041
In the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE042
active power for the x-th cycle,
Figure 100002_DEST_PATH_IMAGE043
Is the x-th period reactive power,
Figure 100002_DEST_PATH_IMAGE044
Is the x-th cycle direct currentThe components of the first and second images are,
Figure 100002_DEST_PATH_IMAGE045
is the x periodkThe sub-harmonic current effective value.
Preferably, in step 3, the feature vector is a set of 9-dimensional feature vectors, including an active power increment, a reactive power increment, an increment of a direct current component, and an increment of a 2-7 th harmonic current component of a load edge occurring at the mth cycle; the electrical appliance feature library comprises active power increment, reactive power increment, increment of direct current component and increment of 2-7 harmonic current component of each electrical appliance.
The step 3 comprises the following steps:
step 3.1, traversing and calculating all Euclidean distances between the characteristic vectors and all the electric appliance characteristic vectors in the electric appliance characteristic library, and solving a maximum Euclidean distance value in the characteristic vectors;
step 3.2, if the obtained maximum Euclidean distance value is larger than a set distance threshold value, taking the type of the electric appliance in the electric appliance feature library corresponding to the maximum Euclidean distance value as the type of the load edge corresponding to the feature vector; wherein the set distance threshold value is 200.
Preferably, in step 4, the different types of the load edges are different types of the electric appliances in the electric appliance feature library, and the types of the electric appliances comprise electric appliances with stable working states and electric appliances with fluctuating working states; the electric appliances with stable working states are classified into a first priority load, and the electric appliances with fluctuating working states are classified into a second priority load.
Preferably, step 5 comprises:
step 5.1, for any load edge, shutting down the load to edge event
Figure 100002_DEST_PATH_IMAGE046
With a load start edge event occurring before the time of the load stop edge event
Figure 100002_DEST_PATH_IMAGE047
Matching to obtain load identification pairs
Figure 100002_DEST_PATH_IMAGE048
(ii) a Wherein the load stops the edge event
Figure 836621DEST_PATH_IMAGE046
Take place inmAt cycle, load initiated edge events
Figure 385414DEST_PATH_IMAGE034
Take place innAt period of time, satisfy
Figure 100002_DEST_PATH_IMAGE049
And the load edges corresponding to the two events are of the same kind;
step 5.2, identify the pair in any load
Figure 286505DEST_PATH_IMAGE048
Next, the first priority load is decomposed to obtain the consumed power of the corresponding electrical appliance according to the following relation:
Figure 100002_DEST_PATH_IMAGE050
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE051
for identifying the pair in the load
Figure 615855DEST_PATH_IMAGE048
The first priority load below corresponds to the consumed power of the appliance,
Figure 96646DEST_PATH_IMAGE037
the first priority load corresponds to the active power increment of the electric appliance under the load shutdown edge event of the mth period,
Figure 100002_DEST_PATH_IMAGE052
under the condition of starting the edge event for the load of the nth period, the first priority load corresponds to the active power increment of the electric appliance;
step 5.3, identify the pair in any load
Figure 31104DEST_PATH_IMAGE048
Next, integrating the total power of the second priority load from the starting time to the ending time of the second priority load to obtain the total electric quantity;
and 5.4, subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
The non-intrusive load identification electric quantity decomposition system based on priority distribution comprises: the system comprises a high-frequency sampling module, a load edge detection module, a load edge identification module, a load priority distribution module and a load identification electric quantity decomposition module;
the high-frequency sampling module is used for carrying out high-frequency sampling on the current and the voltage of the household bus and carrying out fast Fourier transform on high-frequency sampling values of the current and the voltage so as to obtain active power, reactive power, direct current component and harmonic current component of a power frequency period;
the load edge detection module is used for detecting the load edge by utilizing the active power of the power frequency cycle and adopting a maximum likelihood ratio algorithm, and synchronously calculating the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge;
the load edge identification module is used for extracting a characteristic vector from an active power increment, a reactive power increment, an increment of a direct current component and an increment of a harmonic current component of a load edge; identifying the type of the load edge by using the characteristic vector and each electric appliance characteristic vector in the electric appliance characteristic library through an Euclidean distance judgment method;
the load priority distribution module is used for carrying out multi-priority distribution on various loads according to different types of load edges so as to obtain a first priority load and a second priority load; the first priority load is a stability load, and the second priority load is a volatility load;
the load identification electric quantity decomposition module is used for matching a load shutdown edge with a load start edge for any load edge to form a load identification pair; and under any load identification pair, carrying out electric quantity decomposition on the first priority load to obtain the consumed electric quantity of the corresponding electric appliance, simultaneously obtaining the total electric quantity of the second priority load from starting to stopping, and subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
Compared with the prior art, the load identification electric quantity decomposition strategy based on the multi-priority distribution method has the advantages that under the condition that the load edge is rapidly detected and accurately matched, non-air-conditioning electric appliances such as an electric water heater and an electric hair drier are brought into the first priority load, the electric quantity decomposition is carried out by adopting a trapezoidal method, and the electric quantity calculation is accurate and reliable by adopting the trapezoidal method due to the stable working state of the electric appliances; and aiming at the air-conditioning electric appliance, a second priority load is included, and the electric quantity is estimated by subtracting the stable non-air-conditioning electric quantity from the total electric quantity, so that the problem of inaccurate electric quantity decomposition caused by fluctuation of the working state of the air-conditioning electric appliance is effectively solved.
In addition, the method provided by the invention is easy to operate and realize in the methods of data acquisition, data preprocessing, event detection, feature extraction, load identification and the like, so that the method is more suitable for the localized load identification equipment with real-time acquisition and real-time processing analysis.
Drawings
FIG. 1 is a block flow diagram of a non-intrusive load identification power decomposition method based on priority assignment according to the present invention;
fig. 2 is a schematic structural diagram of a non-intrusive load identification electric quantity decomposition system based on priority assignment according to the present invention.
Detailed Description
The present application is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present application is not limited thereby.
Referring to fig. 1, a method for non-intrusive load identification power decomposition based on multi-priority distribution includes:
step 1, carrying out high-frequency sampling on current and voltage of an in-house bus, and carrying out fast Fourier transform on high-frequency sampling values of the current and the voltage so as to obtain active power, reactive power, direct current component and harmonic current component of a power frequency period.
In particular, the amount of the solvent to be used,
in the step 1, an electric energy meter or other terminal equipment is used for carrying out high-frequency sampling on current and voltage of a home bus, and a metering core outputs high-frequency sampling current and voltage; wherein, the sampling frequency of the high-frequency sampling period is 1.6 kHz.
It should be noted that, those skilled in the art can calibrate the optimal high-frequency sampling frequency according to the requirements of multiple experiments or data acquisition, and in the preferred embodiment, 1.6kHz is used as the high-frequency sampling frequency, which is a non-limiting preferred choice.
Further, step 1 comprises:
step 1.1, collecting current and voltage of the household bus based on a high-frequency sampling period, and respectively forming a current sampling sequence by utilizing high-frequency sampling values of the current and the voltage
Figure 821205DEST_PATH_IMAGE001
And a voltage sampling sequence
Figure 309955DEST_PATH_IMAGE002
Step 1.2, sampling sequence of current according to the following relation
Figure 578126DEST_PATH_IMAGE001
And a voltage sampling sequence
Figure 646051DEST_PATH_IMAGE002
Performing fast Fourier transform to obtain current amplitude result sequence
Figure 607054DEST_PATH_IMAGE003
Current phase result sequence
Figure 583100DEST_PATH_IMAGE004
Voltage amplitude result sequence
Figure 389382DEST_PATH_IMAGE005
And a sequence of voltage phase results
Figure 829591DEST_PATH_IMAGE006
Figure 899178DEST_PATH_IMAGE007
In the formula (I), the compound is shown in the specification,
Figure 909991DEST_PATH_IMAGE008
is a fast fourier transform of a sequence of current samples,
Figure 785543DEST_PATH_IMAGE009
is a fast fourier transform of a sequence of voltage samples.
Step 1.3, calculating the active power, the reactive power, the direct current component and the second component of the power frequency period according to the following relational expressionskSubharmonic current component:
Figure 814679DEST_PATH_IMAGE010
Figure 320746DEST_PATH_IMAGE011
in the formula (I), the compound is shown in the specification,kin order to be the harmonic frequency, the frequency of the harmonic wave,Pis the active power of the power frequency period,QIs the reactive power of the power frequency period,
Figure 802543DEST_PATH_IMAGE012
Is a direct current of a power frequency period,
Figure 232519DEST_PATH_IMAGE013
is the first of power frequency cyclekThe effective value of the sub-harmonic current,
Figure 319423DEST_PATH_IMAGE014
is as followskThe magnitude of the current at the sub-harmonic,
Figure 793130DEST_PATH_IMAGE015
is as followskThe voltage amplitude of the sub-harmonics,
Figure 27802DEST_PATH_IMAGE016
is as followskThe phase of the current of the sub-harmonic,
Figure 979578DEST_PATH_IMAGE017
is as followskThe phase of the voltage of the sub-harmonic,
Figure 186568DEST_PATH_IMAGE018
is the current amplitude of the fundamental wave.
And 2, detecting the load edge by using the power frequency cycle active power and adopting a maximum likelihood ratio algorithm, and synchronously calculating the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge.
Specifically, step 2 comprises:
step 2.1, arranging active power, reactive power, direct current component and harmonic current component of a power frequency period into an active power sequence according to a period time sequence
Figure 581909DEST_PATH_IMAGE019
Reactive power sequence
Figure 303877DEST_PATH_IMAGE020
Sequence of direct current components
Figure 59343DEST_PATH_IMAGE021
And a firstkSequence of sub-harmonic current components
Figure 855261DEST_PATH_IMAGE022
Step 2.2, for the mth period, constructing a sliding data window for the active power sequence; wherein the sliding data window comprises a front sliding data window
Figure 936350DEST_PATH_IMAGE023
And a rear sliding data window
Figure 633697DEST_PATH_IMAGE024
(ii) a The front sliding data window is an active power sequence formed by 50 active powers before the mth period, and the front sliding data window is an active power sequence formed by 50 active powers after the mth period.
It should be noted that, a person skilled in the art may calibrate the number of data in the active power sequence according to multiple experiments or data acquisition requirements, and the active power sequence composed of 50 active powers in the preferred embodiment is a non-limiting preferred choice.
Step 2.3, calculating the maximum likelihood ratio of the sliding data window at the mth period according to the following relational expression
Figure 192854DEST_PATH_IMAGE025
Figure DEST_PATH_IMAGE053
In the formula (I), the compound is shown in the specification,
Figure 640016DEST_PATH_IMAGE027
is a function of the probability of the normal distribution,
Figure 626427DEST_PATH_IMAGE028
for the mean value of the post-sliding data window,
Figure 260670DEST_PATH_IMAGE029
for the post-sliding data window standard deviation,
Figure 108672DEST_PATH_IMAGE030
is the mean of the front sliding window of data,
Figure 207078DEST_PATH_IMAGE031
is the standard deviation of the front sliding data window.
Step 2.4, if the maximum likelihood ratio of the sliding data window at the mth period is larger than the set event threshold, judging that a load edge event exists at the mth period
Figure 98810DEST_PATH_IMAGE032
(ii) a Wherein the set event threshold value is 0.5.
It should be noted that, in the preferred embodiment, the value of the event threshold is 0.5, which is a non-limiting preferred choice.
Step 2.5, for the presence of load edge events
Figure 220350DEST_PATH_IMAGE032
Period m of (a), if
Figure 386889DEST_PATH_IMAGE033
Then the load edge event is judged to be a load starting edge event
Figure 559376DEST_PATH_IMAGE034
(ii) a On the contrary, if
Figure 90851DEST_PATH_IMAGE035
Then the load edge event is determined to be a load shutdown edge event
Figure 230846DEST_PATH_IMAGE036
Step 2.6, respectively calculating the active power increment of the load edge occurring at the mth period according to the following relational expression
Figure 201076DEST_PATH_IMAGE037
Reactive power increment
Figure 742915DEST_PATH_IMAGE038
D.c. currentIncrement of component
Figure 445292DEST_PATH_IMAGE039
And a firstkIncrement of subharmonic current component
Figure 88894DEST_PATH_IMAGE040
Figure 331657DEST_PATH_IMAGE041
In the formula (I), the compound is shown in the specification,
Figure 931265DEST_PATH_IMAGE042
active power for the x-th cycle,
Figure 601281DEST_PATH_IMAGE043
Is the x-th period reactive power,
Figure 981447DEST_PATH_IMAGE044
Is the x-th period direct current component,
Figure 965583DEST_PATH_IMAGE045
is the x periodkThe sub-harmonic current effective value.
Step 3, extracting characteristic vectors from the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge; and identifying the type of the load edge by using the characteristic vector and the characteristic vectors of the electric appliances in the electric appliance characteristic library through an Euclidean distance judgment method.
Specifically, in step 3, the feature vector is a set of 9-dimensional feature vectors, including an active power increment, a reactive power increment, an increment of a direct current component, and an increment of a 2-7 th harmonic current component of a load edge occurring at the m-th cycle.
The electrical appliance feature library comprises active power increment, reactive power increment, increment of direct current component and increment of 2-7 harmonic current component of each electrical appliance.
Further, step 3 comprises:
and 3.1, traversing and calculating all Euclidean distances between the characteristic vectors and all the electric appliance characteristic vectors in the electric appliance characteristic library, and solving a maximum Euclidean distance value.
Step 3.2, if the obtained maximum Euclidean distance value is larger than a set distance threshold value, taking the type of the electric appliance in the electric appliance feature library corresponding to the maximum Euclidean distance value as the type of the load edge corresponding to the feature vector; wherein the set distance threshold value is 200.
It should be noted that, in the preferred embodiment, the distance threshold is 200, which is a non-limiting preferred choice.
Step 4, according to different types of load edges, performing multi-priority distribution on various loads to obtain a first priority load and a second priority load; the first priority load is a stability load, and the second priority load is a volatility load.
Specifically, in step 4, the different types of the load edges are different types of the respective appliances in the appliance feature library, and the types of the appliances include appliances with stable working states and appliances with fluctuating working states.
The electric appliances with stable working states are classified into a first priority load, and the electric appliances with fluctuating working states are classified into a second priority load.
In the preferred embodiment of the present invention, the first priority load includes but is not limited to non-air-conditioning appliances such as electric water heaters and hair dryers, and the second priority load is an air conditioner.
Step 5, for any load edge, matching the load shutdown edge with the load startup edge to form a load identification pair; and under any load identification pair, carrying out electric quantity decomposition on the first priority load to obtain the consumed electric quantity of the corresponding electric appliance, simultaneously obtaining the total electric quantity of the second priority load from starting to stopping, and subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
Further, step 5 comprises:
step 5.1, for any loadEdge, load off edge event
Figure 229818DEST_PATH_IMAGE046
With a load start edge event occurring before the time of the load stop edge event
Figure 70735DEST_PATH_IMAGE047
Matching to obtain load identification pairs
Figure 938197DEST_PATH_IMAGE048
(ii) a Wherein the load stops the edge event
Figure 522762DEST_PATH_IMAGE046
Take place inmAt cycle, load initiated edge events
Figure 644433DEST_PATH_IMAGE047
Take place innAt period of time, satisfy
Figure 859514DEST_PATH_IMAGE049
And the load edges for both events are of the same kind.
Step 5.2, identify the pair in any load
Figure 683113DEST_PATH_IMAGE048
Then, the power consumption of the corresponding electrical appliance is obtained by decomposing the power of the first priority load by using a trapezoidal method according to the following relation:
Figure 71369DEST_PATH_IMAGE050
in the formula (I), the compound is shown in the specification,
Figure 296814DEST_PATH_IMAGE051
for identifying the pair in the load
Figure 230266DEST_PATH_IMAGE048
First ofThe priority load corresponds to the consumed electric power of the electric appliance,
Figure 72320DEST_PATH_IMAGE037
the first priority load corresponds to the active power increment of the electric appliance under the load shutdown edge event of the mth period,
Figure 998688DEST_PATH_IMAGE052
under the condition of starting the edge event for the load of the nth period, the first priority load corresponds to the active power increment of the electric appliance;
step 5.3, identify the pair in any load
Figure 78639DEST_PATH_IMAGE048
And integrating the total power of the second priority load from the starting time to the ending time of the second priority load to obtain the total electric quantity.
And 5.4, subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
As shown in fig. 2, the non-intrusive load identification electric quantity decomposition system based on priority distribution includes: the system comprises a high-frequency sampling module, a load edge detection module, a load edge identification module, a load priority distribution module and a load identification electric quantity decomposition module.
The high-frequency sampling module is used for carrying out high-frequency sampling on the current and the voltage of the household bus and carrying out fast Fourier transform on high-frequency sampling values of the current and the voltage so as to obtain active power, reactive power, direct current component and harmonic current component of a power frequency period;
the load edge detection module is used for detecting the load edge by utilizing the active power of the power frequency cycle and adopting a maximum likelihood ratio algorithm, and synchronously calculating the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge;
the load edge identification module is used for extracting a characteristic vector from an active power increment, a reactive power increment, an increment of a direct current component and an increment of a harmonic current component of a load edge; identifying the type of the load edge by using the characteristic vector and each electric appliance characteristic vector in the electric appliance characteristic library through an Euclidean distance judgment method;
the load priority distribution module is used for carrying out multi-priority distribution on various loads according to different types of load edges so as to obtain a first priority load and a second priority load; the first priority load is a stability load, and the second priority load is a volatility load;
the load identification electric quantity decomposition module is used for matching a load shutdown edge with a load start edge for any load edge to form a load identification pair; and under any load identification pair, carrying out electric quantity decomposition on the first priority load to obtain the consumed electric quantity of the corresponding electric appliance, simultaneously obtaining the total electric quantity of the second priority load from starting to stopping, and subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
Compared with the prior art, the load identification electric quantity decomposition strategy based on the multi-priority distribution method has the advantages that under the condition that the load edge is rapidly detected and accurately matched, non-air-conditioning electric appliances such as an electric water heater and an electric hair drier are brought into the first priority load, the electric quantity decomposition is carried out by adopting a trapezoidal method, and the electric quantity calculation is accurate and reliable by adopting the trapezoidal method due to the stable working state of the electric appliances; and aiming at the air-conditioning electric appliance, a second priority load is included, and the electric quantity is estimated by subtracting the stable non-air-conditioning electric quantity from the total electric quantity, so that the problem of inaccurate electric quantity decomposition caused by fluctuation of the working state of the air-conditioning electric appliance is effectively solved.
In addition, the method provided by the invention is easy to operate and realize in the methods of data acquisition, data preprocessing, event detection, feature extraction, load identification and the like, so that the method is more suitable for the localized load identification equipment with real-time acquisition and real-time processing analysis.
The present applicant has described and illustrated embodiments of the present invention in detail with reference to the accompanying drawings, but it should be understood by those skilled in the art that the above embodiments are merely preferred embodiments of the present invention, and the detailed description is only for the purpose of helping the reader to better understand the spirit of the present invention, and not for limiting the scope of the present invention, and on the contrary, any improvement or modification made based on the spirit of the present invention should fall within the scope of the present invention.

Claims (9)

1. The non-intrusive load identification electric quantity decomposition method based on priority distribution is characterized in that,
the decomposition method comprises the following steps:
step 1, carrying out high-frequency sampling on current and voltage of an in-house bus, and carrying out fast Fourier transform on high-frequency sampling values of the current and the voltage so as to obtain active power, reactive power, direct current component and harmonic current component of a power frequency period;
step 2, detecting a load edge by using the active power of a power frequency cycle and adopting a maximum likelihood ratio algorithm, and synchronously calculating the active power increment, the reactive power increment, the increment of a direct current component and the increment of a harmonic current component of the load edge;
step 3, extracting characteristic vectors from the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge; identifying the type of the load edge by using the characteristic vector and each electric appliance characteristic vector in the electric appliance characteristic library through an Euclidean distance judgment method;
step 4, according to different types of load edges, performing multi-priority distribution on various loads to obtain a first priority load and a second priority load; the first priority load is a stability load, and the second priority load is a volatility load;
step 5, for any load edge, matching the load shutdown edge with the load startup edge to form a load identification pair; and under any load identification pair, carrying out electric quantity decomposition on the first priority load to obtain the consumed electric quantity of the corresponding electric appliance, simultaneously obtaining the total electric quantity of the second priority load from starting to stopping, and subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
2. The method of non-intrusive load recognition battery power decomposition based on priority assignment as defined in claim 1,
in the step 1, an electric energy meter or other terminal equipment is used for carrying out high-frequency sampling on current and voltage of a home bus, and a metering core outputs high-frequency sampling current and voltage; wherein, the sampling frequency of the high-frequency sampling period is 1.6 kHz.
3. The method of non-intrusive load recognition battery decomposition based on priority assignment as defined in claim 2,
the step 1 comprises the following steps:
step 1.1, collecting current and voltage of the household bus based on a high-frequency sampling period, and respectively forming a current sampling sequence by utilizing high-frequency sampling values of the current and the voltage
Figure DEST_PATH_IMAGE001
And a voltage sampling sequence
Figure DEST_PATH_IMAGE002
Step 1.2, sampling sequence of current according to the following relation
Figure 544018DEST_PATH_IMAGE001
And a voltage sampling sequence
Figure 539787DEST_PATH_IMAGE002
Performing fast Fourier transform to obtain current amplitude result sequence
Figure DEST_PATH_IMAGE003
Current phase result sequence
Figure DEST_PATH_IMAGE004
Voltage amplitude result sequence
Figure DEST_PATH_IMAGE005
And a sequence of voltage phase results
Figure DEST_PATH_IMAGE006
Figure DEST_PATH_IMAGE007
In the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE008
is a fast fourier transform of a sequence of current samples,
Figure DEST_PATH_IMAGE009
fast fourier transform of a sequence of voltage samples;
step 1.3, calculating the active power, the reactive power, the direct current component and the second component of the power frequency period according to the following relational expressionskSubharmonic current component:
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE011
in the formula (I), the compound is shown in the specification,kin order to be the harmonic frequency, the frequency of the harmonic wave,Pis the active power of the power frequency period,QIs the reactive power of the power frequency period,
Figure DEST_PATH_IMAGE012
Is a direct current of a power frequency period,
Figure DEST_PATH_IMAGE013
is the first of power frequency cyclekSubharmonic electricityThe effective value of the stream is,
Figure DEST_PATH_IMAGE014
is as followskThe magnitude of the current at the sub-harmonic,
Figure DEST_PATH_IMAGE015
is as followskThe voltage amplitude of the sub-harmonics,
Figure DEST_PATH_IMAGE016
is as followskThe phase of the current of the sub-harmonic,
Figure DEST_PATH_IMAGE017
is as followskThe phase of the voltage of the sub-harmonic,
Figure DEST_PATH_IMAGE018
is the current amplitude of the fundamental wave.
4. The method of non-intrusive load recognition battery power decomposition based on priority assignment as defined in claim 1,
the step 2 comprises the following steps:
step 2.1, arranging active power, reactive power, direct current component and harmonic current component of a power frequency period into an active power sequence according to a period time sequence
Figure DEST_PATH_IMAGE019
Reactive power sequence
Figure DEST_PATH_IMAGE020
Sequence of direct current components
Figure DEST_PATH_IMAGE021
And a firstkSequence of sub-harmonic current components
Figure DEST_PATH_IMAGE022
Step 2.2, forAt m periods, constructing a sliding data window for the active power sequence; wherein the sliding data window comprises a front sliding data window
Figure DEST_PATH_IMAGE023
And a rear sliding data window
Figure DEST_PATH_IMAGE024
(ii) a The front sliding data window is an active power sequence formed by 50 active powers before the mth period, and the front sliding data window is an active power sequence formed by 50 active powers after the mth period;
step 2.3, calculating the maximum likelihood ratio of the sliding data window at the mth period according to the following relational expression
Figure DEST_PATH_IMAGE025
Figure DEST_PATH_IMAGE026
In the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE027
is a function of the probability of the normal distribution,
Figure DEST_PATH_IMAGE028
for the mean value of the post-sliding data window,
Figure DEST_PATH_IMAGE029
for the post-sliding data window standard deviation,
Figure DEST_PATH_IMAGE030
is the mean of the front sliding window of data,
Figure DEST_PATH_IMAGE031
is the standard deviation of the front sliding data window;
step 2.4, if the sliding data window at the m-th periodIs greater than the set event threshold value, the load edge event is judged to exist at the mth period
Figure DEST_PATH_IMAGE032
(ii) a Wherein the set event threshold value is 0.5;
step 2.5, for the presence of load edge events
Figure 733571DEST_PATH_IMAGE032
Period m of (a), if
Figure DEST_PATH_IMAGE033
Then the load edge event is judged to be a load starting edge event
Figure DEST_PATH_IMAGE034
(ii) a On the contrary, if
Figure DEST_PATH_IMAGE035
Then the load edge event is determined to be a load shutdown edge event
Figure DEST_PATH_IMAGE036
Step 2.6, respectively calculating the active power increment of the load edge occurring at the mth period according to the following relational expression
Figure DEST_PATH_IMAGE037
Reactive power increment
Figure DEST_PATH_IMAGE038
Delta of DC current component
Figure DEST_PATH_IMAGE039
And a firstkIncrement of subharmonic current component
Figure DEST_PATH_IMAGE040
Figure DEST_PATH_IMAGE041
In the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE042
active power for the x-th cycle,
Figure DEST_PATH_IMAGE043
Is the x-th period reactive power,
Figure DEST_PATH_IMAGE044
Is the x-th period direct current component,
Figure DEST_PATH_IMAGE045
is the x periodkThe sub-harmonic current effective value.
5. The method of non-intrusive load recognition battery power decomposition based on priority assignment as defined in claim 1,
in step 3, the eigenvector is a group of 9-dimensional eigenvectors, including an active power increment, a reactive power increment, an increment of a direct current component and an increment of a 2-7 harmonic current component of a load edge occurring at the mth cycle;
the electric appliance feature library comprises active power increment, reactive power increment, increment of direct current component and increment of 2-7 harmonic current component of each electric appliance.
6. The method of non-intrusive load recognition battery decomposition based on priority assignment as defined in claim 5,
the step 3 comprises the following steps:
step 3.1, traversing and calculating all Euclidean distances between the characteristic vectors and all the electric appliance characteristic vectors in the electric appliance characteristic library, and solving a maximum Euclidean distance value in the characteristic vectors;
step 3.2, if the obtained maximum Euclidean distance value is larger than a set distance threshold value, taking the type of the electric appliance in the electric appliance feature library corresponding to the maximum Euclidean distance value as the type of the load edge corresponding to the feature vector; wherein the set distance threshold value is 200.
7. The method of non-intrusive load recognition battery power decomposition based on priority assignment as defined in claim 1,
in step 4, the different types of the load edges are different types of the electric appliances in the electric appliance feature library, and the types of the electric appliances comprise electric appliances with stable working states and electric appliances with fluctuating working states;
the electric appliances with stable working states are classified into a first priority load, and the electric appliances with fluctuating working states are classified into a second priority load.
8. The method of non-intrusive load recognition battery power decomposition based on priority assignment as defined in claim 7,
the step 5 comprises the following steps:
step 5.1, for any load edge, shutting down the load to edge event
Figure DEST_PATH_IMAGE046
With a load start edge event occurring before the time of the load stop edge event
Figure DEST_PATH_IMAGE047
Matching to obtain load identification pairs
Figure DEST_PATH_IMAGE048
(ii) a Wherein the load stops the edge event
Figure 909924DEST_PATH_IMAGE046
Take place inmAt cycle, load initiated edge events
Figure 255454DEST_PATH_IMAGE047
Take place innAt period of time, satisfy
Figure DEST_PATH_IMAGE049
And the load edges corresponding to the two events are of the same kind;
step 5.2, identify the pair in any load
Figure 890966DEST_PATH_IMAGE048
Next, the first priority load is decomposed to obtain the consumed power of the corresponding electrical appliance according to the following relation:
Figure DEST_PATH_IMAGE050
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE051
for identifying the pair in the load
Figure 751475DEST_PATH_IMAGE048
The first priority load below corresponds to the consumed power of the appliance,
Figure 730801DEST_PATH_IMAGE037
the first priority load corresponds to the active power increment of the electric appliance under the load shutdown edge event of the mth period,
Figure DEST_PATH_IMAGE052
under the condition of starting the edge event for the load of the nth period, the first priority load corresponds to the active power increment of the electric appliance;
step 5.3, identify the pair in any load
Figure 461997DEST_PATH_IMAGE048
Next, integrating the total power of the second priority load from the starting time to the ending time of the second priority load to obtain the total electric quantity;
and 5.4, subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
9. The non-intrusive load identification electric quantity decomposition system based on priority distribution realized by using the non-intrusive load identification electric quantity decomposition method based on priority distribution of any one of claims 1 to 8,
the system comprises a high-frequency sampling module, a load edge detection module, a load edge identification module, a load priority distribution module and a load identification electric quantity decomposition module;
the high-frequency sampling module is used for carrying out high-frequency sampling on the current and the voltage of the household bus and carrying out fast Fourier transform on high-frequency sampling values of the current and the voltage so as to obtain active power, reactive power, direct current component and harmonic current component of a power frequency period;
the load edge detection module is used for detecting the load edge by utilizing the active power of the power frequency cycle and adopting a maximum likelihood ratio algorithm, and synchronously calculating the active power increment, the reactive power increment, the increment of the direct current component and the increment of the harmonic current component of the load edge;
the load edge identification module is used for extracting a characteristic vector from an active power increment, a reactive power increment, an increment of a direct current component and an increment of a harmonic current component of a load edge; identifying the type of the load edge by using the characteristic vector and each electric appliance characteristic vector in the electric appliance characteristic library through an Euclidean distance judgment method;
the load priority distribution module is used for carrying out multi-priority distribution on various loads according to different types of load edges so as to obtain a first priority load and a second priority load; the first priority load is a stability load, and the second priority load is a volatility load;
the load identification electric quantity decomposition module is used for matching a load shutdown edge with a load start edge for any load edge to form a load identification pair; and under any load identification pair, carrying out electric quantity decomposition on the first priority load to obtain the consumed electric quantity of the corresponding electric appliance, simultaneously obtaining the total electric quantity of the second priority load from starting to stopping, and subtracting the consumed electric quantity of the electric appliance corresponding to the first priority load from the total electric quantity to obtain the consumed electric quantity of the electric appliance corresponding to the second priority load.
CN202110911018.7A 2021-08-10 2021-08-10 Non-invasive load identification electric quantity decomposition method and system based on priority distribution Active CN113361831B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110911018.7A CN113361831B (en) 2021-08-10 2021-08-10 Non-invasive load identification electric quantity decomposition method and system based on priority distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110911018.7A CN113361831B (en) 2021-08-10 2021-08-10 Non-invasive load identification electric quantity decomposition method and system based on priority distribution

Publications (2)

Publication Number Publication Date
CN113361831A CN113361831A (en) 2021-09-07
CN113361831B true CN113361831B (en) 2021-11-09

Family

ID=77540865

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110911018.7A Active CN113361831B (en) 2021-08-10 2021-08-10 Non-invasive load identification electric quantity decomposition method and system based on priority distribution

Country Status (1)

Country Link
CN (1) CN113361831B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113985266B (en) * 2021-10-28 2024-02-09 中国计量大学现代科技学院 Centralized load identification method for multiple intelligent electric meters
CN115180161B (en) * 2022-05-18 2023-03-21 北京锐士装备科技有限公司 Driving equipment and driving method for unmanned aerial vehicle load device
CN115439017B (en) * 2022-10-24 2023-04-18 国网江苏省电力有限公司营销服务中心 Non-intervention load identification event detection method and device in multi-event scene

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103001230B (en) * 2012-11-16 2014-10-15 天津大学 Non-invasive power load monitoring and decomposing current mode matching method
CN106532719B (en) * 2016-11-28 2019-06-18 国网江苏省电力公司苏州供电公司 Based on electric current second harmonic and idle non-frequency air-conditioning non-intrusion type discrimination method
CN110555369A (en) * 2019-07-16 2019-12-10 浙江工业大学 MLCDTL-based non-intrusive load identification method
CN112366699A (en) * 2020-11-03 2021-02-12 河海大学 Household energy double-layer optimization method for realizing interaction between power grid side and user side
CN112580471A (en) * 2020-12-11 2021-03-30 浙江大学 Non-invasive load identification method based on AdaBoost feature extraction and RNN model

Also Published As

Publication number Publication date
CN113361831A (en) 2021-09-07

Similar Documents

Publication Publication Date Title
CN113361831B (en) Non-invasive load identification electric quantity decomposition method and system based on priority distribution
CN110956220B (en) Non-invasive household appliance load identification method
CN106936129B (en) Power load identification method and system based on multi-feature fusion
CN108021736B (en) Load switching action monitoring method based on sliding window residual error model
CN111830347B (en) Two-stage non-invasive load monitoring method based on event
CN111564842B (en) Method for statistical estimation of medium-and-long-term energy consumption in non-invasive electric load monitoring
CN111027408A (en) Load identification method based on support vector machine and V-I curve characteristics
CN110488112A (en) Classification metering method non-intrusion type load identification and its realized based on recognition result
CN111209976B (en) Non-invasive load monitoring method and system
CN106093565A (en) A kind of electricity subentry measurement method and device based on steady state characteristic Waveform Matching
CN111932051A (en) Malicious behavior detection method based on non-invasive power terminal time sequence monitoring
CN113567794A (en) Electric bicycle indoor charging identification method and system based on dynamic time warping
CN110555787B (en) Non-invasive load identification method based on improved particle swarm optimization
Zhou et al. Neural network pattern recognition based non-intrusive load monitoring for a residential energy management system
CN114859169A (en) Intelligent identification method and system for distribution transformer outgoing line load and storage medium
CN103246814A (en) Personal electric device state identification method based on K-means modeling
Mei et al. Random forest based adaptive non-intrusive load identification
Liu et al. Nonintrusive load disaggregation combining with external attention mechanism and Seq2Piont
CN115687999A (en) Non-invasive load sensing method fusing multi-time scale electric appliance characteristics
CN114662576A (en) Non-invasive storage battery car charging detection method and system based on supervised classification
CN115146709A (en) Non-invasive load identification method facing surface end
CN114759558A (en) Non-invasive online rapid detection method for charging load of electric bicycle
CN112862163A (en) Short-term power load prediction system based on DTW algorithm
Du et al. Two-Step Non-Intrusive Load Monitoring Based on High-Frequency Current Data
CN113420728A (en) Non-invasive air conditioner load identification method and system integrating multi-time scale information

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210907

Assignee: Nanjing Shuxing Power Technology Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2023110000054

Denomination of invention: A Non Invasive Load Identification Method and System for Energy Decomposition Based on Priority Allocation

Granted publication date: 20211109

License type: Common License

Record date: 20230410

Application publication date: 20210907

Assignee: NANJING FEAT ELECTRONIC TECHNOLOGY CO.,LTD.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2023110000055

Denomination of invention: A Non Invasive Load Identification Method and System for Energy Decomposition Based on Priority Allocation

Granted publication date: 20211109

License type: Common License

Record date: 20230410

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210907

Assignee: JIANGSU FANGTIAN POWER TECHNOLOGY Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980000801

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240117

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210907

Assignee: WUHAN SAN FRAN ELECTRONICS Corp.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001384

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240125

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210907

Assignee: JIANGSU WIZPOWER ELECTRIC POWER TECHNOLOGY Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001489

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240129

Application publication date: 20210907

Assignee: Jiangsu suneng Power Technology Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001488

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240129

Application publication date: 20210907

Assignee: WASION Group Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001629

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240130

Application publication date: 20210907

Assignee: YANTAI DONGFANG WISDOM ELECTRIC Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001616

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240130

Application publication date: 20210907

Assignee: HENAN XUJI METERING Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001535

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240129

Application publication date: 20210907

Assignee: JIANGSU LINYANG ENERGY Co.,Ltd.

Assignor: State Grid Jiangsu Electric Power Co.,Ltd. Marketing Service Center

Contract record no.: X2024980001534

Denomination of invention: A non-invasive load identification method and system for electricity quantity decomposition based on priority allocation

Granted publication date: 20211109

License type: Common License

Record date: 20240129