CN108872665B - Load switch event detection method and system based on minimum information criterion - Google Patents

Load switch event detection method and system based on minimum information criterion Download PDF

Info

Publication number
CN108872665B
CN108872665B CN201811097516.7A CN201811097516A CN108872665B CN 108872665 B CN108872665 B CN 108872665B CN 201811097516 A CN201811097516 A CN 201811097516A CN 108872665 B CN108872665 B CN 108872665B
Authority
CN
China
Prior art keywords
minimum information
information criterion
max
time period
time
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.)
Expired - Fee Related
Application number
CN201811097516.7A
Other languages
Chinese (zh)
Other versions
CN108872665A (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.)
Guangdong University of Petrochemical Technology
Original Assignee
Guangdong University of Petrochemical Technology
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 Guangdong University of Petrochemical Technology filed Critical Guangdong University of Petrochemical Technology
Priority to CN201811097516.7A priority Critical patent/CN108872665B/en
Publication of CN108872665A publication Critical patent/CN108872665A/en
Application granted granted Critical
Publication of CN108872665B publication Critical patent/CN108872665B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R11/00—Electromechanical arrangements for measuring time integral of electric power or current, e.g. of consumption
    • G01R11/48—Meters specially adapted for measuring real or reactive components; Meters specially adapted for measuring apparent energy

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a load switch event detection method and system based on a minimum information criterion. The method comprises the following steps: establishing a switching event detection time sequence; obtaining the active power of each moment in the time sequence, and calculating the active power difference value of two adjacent moments; determining a first time period of occurrence of a load switch event according to the active power difference value; calculating each power difference value in a first time period by adopting a minimum information criterion, and determining a time point corresponding to the maximum value of the minimum information criterion in the first time period as a first time point; selecting a plurality of moments before and after the first time point respectively, and determining a second time period of the load switch event; and calculating each power difference value in the second time period by adopting a minimum information criterion, and determining a time point corresponding to the maximum value of the minimum information criterion in the second time period as the moment when the load switch event occurs. The method and the system can eliminate the influence of noise to a great extent and can effectively resist the noise.

Description

Load switch event detection method and system based on minimum information criterion
Technical Field
The invention relates to the technical field of load switch event detection, in particular to a load switch event detection method and system based on a minimum information criterion.
Background
The load switch event is an operation of turning on or off a power switch of a load or an electrical device. Load switch event detection is the most important step in energy decomposition, and energy decomposition refers to decomposing the power value read by an electric meter into the power value consumed by a single load. The current energy decomposition of the electric load is mainly divided into an invasive load decomposition method and a non-invasive load decomposition method. The non-invasive load decomposition method does not need to install monitoring equipment on internal electric equipment of the load, and can obtain the load information of each electric equipment only according to the total information of the electric load.
In the non-invasive load decomposition algorithm, the detection of the switching event of the electrical equipment is the most important link. The existing detection method can identify the switching event with a relatively large active power change value (such as 70W). Due to the noise in the system, the active power increase caused by some of the consumers with active power close to the threshold (i.e. 50W) at startup may be attenuated by the system noise and thus not accurately identified. Some household appliances (such as a crusher and a juicer) using a motor can generate serious impulse noise, the influence of common event detection algorithm noise is large, the detection performance is reduced sharply, especially, the noise has a very large influence on the confirmation of the event occurrence time, and a large error is often generated in the event occurrence time.
Disclosure of Invention
The invention aims to provide a load switch event detection method and system based on a minimum information criterion, which can effectively eliminate the influence of noise.
In order to achieve the purpose, the invention provides the following scheme:
a method of load switch event detection based on minimum information criteria, comprising:
acquiring the occurrence time of a marked load switch event, selecting a plurality of times with the same number before and after the occurrence time, and establishing a switch event detection time sequence according to the time sequence;
obtaining the active power of each moment in the time sequence, and calculating the active power difference value of two adjacent moments;
determining a first time period of occurrence of a load switch event according to the active power difference value;
calculating each power difference value in the first time period by adopting a minimum information criterion, comparing minimum information criterion values to obtain a maximum value of the minimum information criterion in the first time period, and determining a time point corresponding to the maximum value of the minimum information criterion in the first time period as a first time point;
selecting a plurality of moments before and after the first time point respectively, and determining a second time period of the load switch event;
calculating each power difference value in the second time period by adopting a minimum information criterion, comparing the minimum information criterion values to obtain the maximum value of the minimum information criterion of the second time period, and determining the time point corresponding to the maximum value of the minimum information criterion of the second time period as a second time point; the second time point is the moment when the load switch event occurs.
Optionally, determining a first time period of occurrence of a load switch event according to the active power difference specifically includes:
comparing the active power difference value to obtain the maximum value of the active power difference value;
determining the sequence number i of the moment of occurrence of the load switch event corresponding to the maximum active power difference value in the switch event detection time sequencemaxThen the first time period is [0, i ]max]。
Optionally, the calculating the power difference values in the first time period by using a minimum information criterion specifically includes:
to [0, i ]max-1]Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
AIC(i)＝ilog[var{ΔP(0:i)}]+(N-i-1)log[var{ΔP(i+1:imax)}]
wherein i is the serial number of the time sequence of the ith moment after the time sequence is arranged, AIC (i) is the minimum information criterion value with the serial number of i, and var { delta P (0: i) } is { delta P0,ΔP1,…,ΔPiThe variance of }, Δ PiIs the difference value of the active power corresponding to the serial number i and the active power corresponding to the serial number i-1, N is the total number of the active power in the time sequence, and var { delta P (i +1: i:)max) Is as
Figure BDA0001805870710000021
Variance of (1), Δ P0＝0，AIC(0)＝0，AIC(imax-1)＝0；
To imaxThe calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
AIC(imax)＝imaxlog[var{ΔP(0:imax)}]
wherein, var { Δ P (0: i)max) Is { Δ P }0,ΔP1,…,ΔPmaxThe variance of.
Optionally, a plurality of times are respectively selected before and after the first time point, and a second time period in which a load switch event occurs is determined, which specifically includes:
determining the serial number i of the time sequence of the first time point0 max；
At the serial number i0 maxFront selection i0The number of each of the serial numbers is,
Figure BDA0001805870710000031
Figure BDA0001805870710000032
represents rounding down;
at the serial number i0 maxPost selection of i1The number of each of the serial numbers is,
Figure BDA0001805870710000033
Figure BDA0001805870710000034
a second time period during which a load switch event occurs.
Optionally, calculating each power difference value in the second time period by using a minimum information criterion specifically includes:
to pair
Figure BDA0001805870710000035
Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA0001805870710000036
wherein,
Figure BDA0001805870710000037
AIC*(i*) Is a serial number i*Is determined by the minimum information criterion value of (c),
Figure BDA0001805870710000038
is composed of
Figure BDA0001805870710000039
The variance of (a) is determined,
Figure BDA00018058707100000310
is a serial number i*Corresponding active power and serial number i*-1 corresponding difference in active power,
Figure BDA00018058707100000311
is composed of
Figure BDA00018058707100000312
The variance of (a) is determined,
Figure BDA00018058707100000313
to pair
Figure BDA00018058707100000314
The calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA00018058707100000315
wherein,
Figure BDA00018058707100000316
is composed of
Figure BDA00018058707100000317
Variance of (1), Δ P0＝0。
The present invention also provides a load switch event detection system based on minimum information criteria, comprising:
the switching event detection time sequence generation module is used for acquiring the occurrence time of the marked load switching event, selecting a plurality of times with the same number before and after the occurrence time, and establishing a switching event detection time sequence according to the time sequence;
the active power difference value calculation module is used for acquiring active power at each moment in the time sequence and calculating the active power difference value of two adjacent moments;
the first time period generating module is used for determining a first time period of the occurrence of the load switch event according to the active power difference value;
a first time point generating module, configured to calculate each power difference value in the first time period by using a minimum information criterion, compare minimum information criterion values to obtain a minimum information criterion maximum value in the first time period, and determine a time point corresponding to the minimum information criterion maximum value in the first time period as a first time point;
the second time period generation module is used for selecting a plurality of moments before and after the first time point and determining a second time period of the load switch event;
the load switch event occurrence moment generation module is used for calculating each power difference value in the second time period by adopting a minimum information criterion, comparing minimum information criterion values to obtain a minimum information criterion maximum value of the second time period, and determining a time point corresponding to the minimum information criterion maximum value of the second time period as a second time point; the second time point is the moment when the load switch event occurs.
Optionally, the first time period generating module specifically includes:
the active power difference value comparison unit is used for comparing the active power difference value to obtain the maximum value of the active power difference value;
a first time period generating unit, configured to determine a serial number i in the switching event detection time sequence of the load switching event occurrence time corresponding to the maximum active power difference valuemaxThen the first time period is [0, i ]max]。
Optionally, the first time point generating module specifically includes:
a first minimum information criterion calculation unit for calculating a [0, imax-1]Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
AIC(i)＝ilog[var{ΔP(0:i)}]+(N-i-1)log[var{ΔP(i+1:imax)}]
wherein i is a serial number of a time sequence of the ith moment after the time sequence is arranged, and the serial number of AIC (i) isi, var { Δ P (0: i) } is { Δ P0,ΔP1,…,ΔPiThe variance of }, Δ PiIs the difference value of the active power corresponding to the serial number i and the active power corresponding to the serial number i-1, N is the total number of the active power in the time sequence, and var { delta P (i +1: i:)max) Is as
Figure BDA0001805870710000041
Variance of (1), Δ P0＝0，AIC(0)＝0，AIC(imax-1)＝0；
A second minimum information criterion calculation unit for calculating a second minimum information criterion for imaxThe calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
AIC(imax)＝imaxlog[var{ΔP(0:imax)}]
wherein, var { Δ P (0: i)max) Is { Δ P }0,ΔP1,…,ΔPmaxThe variance of.
Optionally, the second time period generating module specifically includes:
a sequence number determining unit for determining the sequence number i of the time sequence of the first time point0 max；
A first sequence number selection unit for selecting the sequence number i0 maxFront selection i0The number of each of the serial numbers is,
Figure BDA0001805870710000042
Figure BDA0001805870710000051
represents rounding down;
a second sequence number selection unit for selecting the sequence number i0 maxPost selection of i1The number of each of the serial numbers is,
Figure BDA0001805870710000052
a second time period generation unit for generating a second time period
Figure BDA0001805870710000053
A second time period is determined for the occurrence of the load switch event.
Optionally, the load switch event occurrence time generation module specifically includes:
a third minimum information criterion calculation unit for calculating
Figure BDA0001805870710000054
Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA0001805870710000055
wherein,
Figure BDA0001805870710000056
AIC*(i*) Is a serial number i*Is determined by the minimum information criterion value of (c),
Figure BDA0001805870710000057
is composed of
Figure BDA0001805870710000058
The variance of (a) is determined,
Figure BDA0001805870710000059
is a serial number i*Corresponding active power and serial number i*-1 corresponding difference in active power,
Figure BDA00018058707100000510
is composed of
Figure BDA00018058707100000511
The variance of (a) is determined,
Figure BDA00018058707100000512
a fourth minimum information criterion calculation unit for calculating
Figure BDA00018058707100000513
The calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA00018058707100000514
wherein,
Figure BDA00018058707100000515
is composed of
Figure BDA00018058707100000516
Variance of (1), Δ P0＝0。
Compared with the prior art, the invention has the beneficial effects that:
the invention provides a load switch event detection method and system based on a minimum information criterion, which can eliminate the influence of noise to a great extent by adopting the minimum information criterion for calculation, can effectively resist the noise, and improves the precision of load switch event detection by shortening the time period for determining the occurrence of the load switch event twice.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flowchart of a method for detecting a load switch event based on a minimum information criterion according to an embodiment of the present invention;
fig. 2 is a diagram of a load switch event detection system based on a minimum information criterion according to a second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a load switching event detection method based on a minimum information criterion, which can effectively eliminate the influence of noise.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
The first embodiment is as follows:
fig. 1 is a flowchart of a load switch event detection method based on a minimum information criterion according to an embodiment of the present invention, and as shown in fig. 1, the load switch event detection method based on the minimum information criterion includes:
step 101: the method comprises the steps of obtaining the occurrence time of a marked load switch event, selecting a plurality of times with the same number before and after the occurrence time, and establishing a switch event detection time sequence according to the time sequence.
Wherein, the marked load switch event occurrence time is obtained by adopting a switch event detection method in the prior art, and the switch event detection method in the prior art comprises the following steps:
step 1: and (3) calculating an absolute value of the difference between two adjacent active power data, judging whether the absolute value is greater than or equal to 30W, if so, executing the step 3, and otherwise, executing the step 2.
Step 2: and (4) after reading the active power data at the next moment, continuing to execute the step 1.
And step 3: the duration T of the occurrence of the event is increased by 1 second and the execution continues with step 4, with an initial value T equal to 0.
And 4, step 4: reading the active power data at the next moment and calculating delta Pt+1＝Pt+1-PtAnd determining Δ Pt+1And (5) whether the absolute value is greater than or equal to 30W, if so, executing the step 5, otherwise, returning to execute the step 6.
And 5: and reading the active power data at the next moment, and executing the step 3.
Step 6: obtaining the end time T + T of the event according to the duration T of the event, and calculating the change value delta P of the active power before and after the event occurst+T＝Pt+T-PtIf Δ Pt+TAnd if the absolute value is greater than or equal to 50W, executing the step 7, otherwise, determining that the absolute value is abnormal, and returning to execute the step 2.
And 7: and outputting a result: according to Δ Pt+TWhether the occurrence is a rising edge event or a falling edge event may be determined. If Δ Pt+TA positive, indicating an increase in active power, is determined as a rising edge event, typically caused by an electrical appliance being put into operation or a state change; otherwise, the active power is reduced, and the active power is determined as a falling edge event, which is generally caused by the fact that the electrical appliance is out of operation or the state changes. time T is the starting time of the event, and time T + T is the ending time of the event.
Step 102: and obtaining the active power of each moment in the time sequence, and calculating the active power difference value of two adjacent moments.
Step 103: and determining a first time period of the load switch event according to the active power difference value.
Comparing the active power difference value to obtain the maximum value of the active power difference value;
determining the sequence number i of the moment of occurrence of the load switch event corresponding to the maximum active power difference value in the switch event detection time sequencemaxThen the first time period is [0, i ]max]。
Step 104: calculating each power difference value in the first time period by adopting a minimum information criterion, which specifically comprises the following steps:
to [0, i ]max-1]Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
AIC(i)＝ilog[var{ΔP(0:i)}]+(N-i-1)log[var{ΔP(i+1:imax)}]
Figure BDA0001805870710000071
Figure BDA0001805870710000072
Figure BDA0001805870710000073
Figure BDA0001805870710000074
wherein i is the serial number of the time sequence of the ith moment after the time sequence is arranged, AIC (i) is the minimum information criterion value with the serial number of i, and var { delta P (0: i) } is { delta P0,ΔP1,…,ΔPiThe variance of }, Δ PiIs the difference value of the active power corresponding to the serial number i and the active power corresponding to the serial number i-1, N is the total number of the active power in the time sequence, and var { delta P (i +1: i:)max) Is as
Figure BDA0001805870710000075
Variance of (1), Δ P0＝0，AIC(0)＝0，AIC(imax-1)＝0。
To imaxThe calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
AIC(imax)＝imaxlog[var{ΔP(0:imax)}]
Figure BDA0001805870710000081
Figure BDA0001805870710000082
wherein, var { Δ P (0: i)max) Is { Δ P }0,ΔP1,…,ΔPmaxThe variance of.
And comparing the minimum information criterion values to obtain the maximum value of the minimum information criterion of the first time period, and determining the time point corresponding to the maximum value of the minimum information criterion of the first time period as the first time point.
Step 105: and selecting a plurality of moments before and after the first time point, and determining a second time period of the load switch event.
Determining the serial number i of the time sequence of the first time point0 max；
At the serial number i0 maxFront selection i0The number of each of the serial numbers is,
Figure BDA0001805870710000083
Figure BDA0001805870710000084
represents rounding down;
at the serial number i0 maxPost selection of i1The number of each of the serial numbers is,
Figure BDA0001805870710000085
Figure BDA0001805870710000086
a second time period during which a load switch event occurs.
Step 106: calculating each power difference value in the second time period by adopting a minimum information criterion, which specifically comprises the following steps:
to pair
Figure BDA0001805870710000087
Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA0001805870710000088
Figure BDA0001805870710000089
Figure BDA00018058707100000810
Figure BDA00018058707100000811
Figure BDA00018058707100000812
wherein,
Figure BDA00018058707100000813
AIC*(i*) Is a serial number i*Is determined by the minimum information criterion value of (c),
Figure BDA00018058707100000814
is composed of
Figure BDA00018058707100000815
The variance of (a) is determined,
Figure BDA00018058707100000816
is a serial number i*Corresponding active power and serial number i*-1 corresponding difference in active power,
Figure BDA00018058707100000817
is composed of
Figure BDA0001805870710000091
The variance of (a) is determined,
Figure BDA0001805870710000092
to pair
Figure BDA0001805870710000093
The calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA0001805870710000094
Figure BDA0001805870710000095
Figure BDA0001805870710000096
wherein,
Figure BDA0001805870710000097
is composed of
Figure BDA0001805870710000098
Variance of (1), Δ P0＝0。
Comparing the minimum information criterion value to obtain the maximum value of the minimum information criterion of a second time period, and determining a time point corresponding to the maximum value of the minimum information criterion of the second time period as a second time point; the second time point is the moment when the load switch event occurs.
Example two:
fig. 2 is a system structure diagram of a load switch event detection method based on a minimum information criterion in an embodiment of the present invention, and as shown in fig. 2, the load switch event detection system based on the minimum information criterion includes:
the switching event detection time sequence generating module 201 is configured to obtain a marked occurrence time of a load switching event, select a plurality of times with the same number before and after the occurrence time, and establish a switching event detection time sequence according to a time sequence.
And the active power difference value calculating module 202 is configured to obtain active power at each time in the time sequence, and calculate an active power difference value between two adjacent times.
The first time period generating module 203 is configured to determine a first time period of occurrence of a load switch event according to the active power difference, and specifically includes:
the active power difference value comparison unit is used for comparing the active power difference value to obtain the maximum value of the active power difference value;
a first time period generating unit, configured to determine a serial number of a load switch event occurrence time corresponding to the maximum active power difference in the switch event detection time sequenceimaxThen the first time period is [0, i ]max]。
A first time point generating module 204, configured to calculate each power difference value in the first time period by using a minimum information criterion, compare minimum information criterion values to obtain a maximum value of the minimum information criterion in the first time period, and determine a time point corresponding to the maximum value of the minimum information criterion in the first time period as the first time point.
The first time point generating module 204 specifically includes:
a first minimum information criterion calculation unit for calculating a [0, imax-1]Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
AIC(i)＝ilog[var{ΔP(0:i)}]+(N-i-1)log[var{ΔP(i+1:imax)}]
Figure BDA0001805870710000101
Figure BDA0001805870710000102
Figure BDA0001805870710000103
Figure BDA0001805870710000104
wherein i is the serial number of the time sequence of the ith moment after the time sequence is arranged, AIC (i) is the minimum information criterion value with the serial number of i, and var { delta P (0: i) } is { delta P0,ΔP1,…,ΔPiThe variance of }, Δ PiIs the difference value of the active power corresponding to the serial number i and the active power corresponding to the serial number i-1, N is the total number of the active power in the time sequence, and var { delta P (i +1: i:)max) Is as
Figure BDA0001805870710000105
Variance of (1), Δ P0＝0，AIC(0)＝0，AIC(imax-1)＝0；
A second minimum information criterion calculation unit for calculating a second minimum information criterion for imaxThe calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
AIC(imax)＝imaxlog[var{ΔP(0:imax)}]
Figure BDA0001805870710000106
Figure BDA0001805870710000107
wherein, var { Δ P (0: i)max) Is { Δ P }0,ΔP1,…,ΔPmaxThe variance of.
The second time period generating module 205 is configured to select multiple times before and after the first time point, and determine a second time period in which a load switch event occurs, where the second time period generating module specifically includes:
a sequence number determining unit for determining the sequence number i of the time sequence of the first time point0 max；
A first sequence number selection unit for selecting the sequence number i0 maxFront selection i0The number of each of the serial numbers is,
Figure BDA0001805870710000108
Figure BDA0001805870710000109
represents rounding down;
a second sequence number selection unit for selecting the sequence number i0 maxPost selection of i1The number of each of the serial numbers is,
Figure BDA0001805870710000111
a second time period generation unit for generating a second time period
Figure BDA0001805870710000112
A second time period is determined for the occurrence of the load switch event.
A load switch event occurrence time generation module 206, configured to calculate each power difference value in the second time period by using a minimum information criterion, compare minimum information criterion values to obtain a minimum information criterion maximum value in the second time period, and determine a time point corresponding to the minimum information criterion maximum value in the second time period as a second time point; the second time point is the moment when the load switch event occurs.
The load switch event occurrence time generation module 206 specifically includes:
a third minimum information criterion calculation unit for calculating
Figure BDA0001805870710000113
Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA0001805870710000114
Figure BDA0001805870710000115
Figure BDA0001805870710000116
Figure BDA0001805870710000117
Figure BDA0001805870710000118
wherein,
Figure BDA0001805870710000119
AIC*(i*) Is a serial number i*Is determined by the minimum information criterion value of (c),
Figure BDA00018058707100001110
is composed of
Figure BDA00018058707100001111
The variance of (a) is determined,
Figure BDA00018058707100001112
is a serial number i*Corresponding active power and serial number i*-1 corresponding difference in active power,
Figure BDA00018058707100001113
is composed of
Figure BDA00018058707100001114
The variance of (a) is determined,
Figure BDA00018058707100001115
a fourth minimum information criterion calculation unit for calculating
Figure BDA00018058707100001116
The calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
Figure BDA00018058707100001117
Figure BDA00018058707100001118
Figure BDA00018058707100001119
wherein,
Figure BDA00018058707100001120
is composed of
Figure BDA00018058707100001121
Variance of (1), Δ P0＝0。
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (2)

1. A method for load switch event detection based on minimum information criteria, comprising:
acquiring the occurrence time of a marked load switch event, selecting a plurality of times with the same number before and after the occurrence time, and establishing a switch event detection time sequence according to the time sequence;
obtaining the active power of each moment in the time sequence, and calculating the active power difference value of two adjacent moments;
determining a first time period of occurrence of a load switch event according to the active power difference, specifically comprising:
comparing the active power difference value to obtain the maximum value of the active power difference value;
determining the sequence number i of the moment of occurrence of the load switch event corresponding to the maximum active power difference value in the switch event detection time sequencemaxThen the first time period is [0, i ]max]；
Calculating each power difference value in the first time period by adopting a minimum information criterion, comparing minimum information criterion values to obtain a maximum value of the minimum information criterion in the first time period, and determining a time point corresponding to the maximum value of the minimum information criterion in the first time period as a first time point;
the calculating each power difference value in the first time period by using a minimum information criterion specifically includes:
to [0, i ]max-1]Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
AIC(i)＝i log[var{ΔP(0:i)}]+(N-i-1)log[var{ΔP(i+1:imax)}]
wherein i is the serial number of the time sequence of the ith moment after the time sequence is arranged, AIC (i) is the minimum information criterion value with the serial number of i, and var { delta P (0: i) } is { delta P0,ΔP1,…,ΔPiThe variance of }, Δ PiIs the difference value of the active power corresponding to the serial number i and the active power corresponding to the serial number i-1, N is the total number of the active power in the time sequence, and var { delta P (i +1: i:)max) Is as
Figure FDA0002461455090000011
Variance of (1), Δ P0＝0，AIC(0)＝0，AIC(imax-1)＝0；
To imaxThe calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
AIC(imax)＝imaxlog[var{ΔP(0:imax)}]
wherein, var { Δ P (0: i)max) Is { Δ P }0,ΔP1,…,ΔPmaxThe variance of } in the mean square wave
Selecting a plurality of moments before and after the first time point, and determining a second time period of the load switch event, specifically comprising:
determining the serial number i of the time sequence of the first time point0 max；
At the serial number i0 maxFront selection i0The number of each of the serial numbers is,
Figure FDA0002461455090000012
Figure FDA0002461455090000013
represents rounding down;
at the serial number i0 maxPost selection of i1The number of each of the serial numbers is,
Figure FDA0002461455090000014
Figure FDA0002461455090000015
a second time period for a load switch event to occur;
calculating each power difference value in the second time period by adopting a minimum information criterion, comparing the minimum information criterion values to obtain the maximum value of the minimum information criterion of the second time period, and determining the time point corresponding to the maximum value of the minimum information criterion of the second time period as a second time point; the second time point is the moment when the load switch event occurs;
calculating each power difference value in the second time period by adopting a minimum information criterion, which specifically comprises the following steps:
to pair
Figure FDA0002461455090000021
Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
Figure FDA0002461455090000022
wherein,
Figure FDA0002461455090000023
AIC*(i*) Is a serial number i*Is determined by the minimum information criterion value of (c),
Figure FDA0002461455090000024
is composed of
Figure FDA0002461455090000025
The variance of (a) is determined,
Figure FDA0002461455090000026
is a serial number i*Corresponding active power and serial number i*-1 corresponding difference in active power,
Figure FDA0002461455090000027
is composed of
Figure FDA0002461455090000028
The variance of (a) is determined,
Figure FDA0002461455090000029
to pair
Figure FDA00024614550900000210
The calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
Figure FDA00024614550900000211
wherein,
Figure FDA00024614550900000212
is composed of
Figure FDA00024614550900000213
Variance of (1), Δ P0＝0。
2. A minimum information criteria based load switch event detection system, comprising:
the switching event detection time sequence generation module is used for acquiring the occurrence time of the marked load switching event, selecting a plurality of times with the same number before and after the occurrence time, and establishing a switching event detection time sequence according to the time sequence;
the active power difference value calculation module is used for acquiring active power at each moment in the time sequence and calculating the active power difference value of two adjacent moments;
the first time period generating module is used for determining a first time period of the occurrence of the load switch event according to the active power difference value; the first time period generating module specifically includes:
the active power difference value comparison unit is used for comparing the active power difference value to obtain the maximum value of the active power difference value;
a first time period generating unit, configured to determine a serial number i in the switching event detection time sequence of the load switching event occurrence time corresponding to the maximum active power difference valuemaxThen the first time period is [0, i ]max]；
A first time point generating module, configured to calculate each power difference value in the first time period by using a minimum information criterion, compare minimum information criterion values to obtain a minimum information criterion maximum value in the first time period, and determine a time point corresponding to the minimum information criterion maximum value in the first time period as a first time point; the first time point generating module specifically includes:
a first minimum information criterion calculation unit for calculating a [0, imax-1]Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
AIC(i)＝ilog[var{ΔP(0:i)}]+(N-i-1)log[var{ΔP(i+1:imax)}]
wherein i is the serial number of the time sequence of the ith moment after the time sequence is arranged, AIC (i) is the minimum information criterion value with the serial number of i, and var { delta P (0: i) } is { delta P0,ΔP1,…,ΔPiThe variance of }, Δ PiIs the difference value of the active power corresponding to the serial number i and the active power corresponding to the serial number i-1, N is the total number of the active power in the time sequence, and var { delta P (i +1: i:)max) Is as
Figure FDA0002461455090000031
Variance of (1), Δ P0＝0，AIC(0)＝0，AIC(imax-1)＝0；
A second minimum information criterion calculation unit for calculating a second minimum information criterion for imaxThe calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
AIC(imax)＝imaxlog[var{ΔP(0:imax)}]
wherein, var { Δ P (0: i)max) Is { Δ P }0,ΔP1,…,ΔPmaxThe variance of };
the second time period generation module is used for selecting a plurality of moments before and after the first time point and determining a second time period of the load switch event; the second time period generation module specifically includes:
a sequence number determining unit for determining the sequence number i of the time sequence of the first time point0 max；
A first sequence number selection unit for selecting the sequence number i0 maxFront selection i0The number of each of the serial numbers is,
Figure FDA0002461455090000032
Figure FDA0002461455090000033
represents rounding down;
a second sequence number selection unit for selecting the sequence number i0 maxPost selection of i1The number of each of the serial numbers is,
Figure FDA0002461455090000034
a second time period generation unit for generating a second time period
Figure FDA0002461455090000035
Determining a second time period for the occurrence of the load switch event;
the load switch event occurrence moment generation module is used for calculating each power difference value in the second time period by adopting a minimum information criterion, comparing minimum information criterion values to obtain a minimum information criterion maximum value of the second time period, and determining a time point corresponding to the minimum information criterion maximum value of the second time period as a second time point; the second time point is the moment when the load switch event occurs; the load switch event occurrence time generation module specifically includes:
a third minimum information criterion calculation unit for calculating
Figure FDA0002461455090000036
Each power difference value is calculated by adopting a minimum information criterion, and the specific formula is as follows:
Figure FDA0002461455090000037
wherein,
Figure FDA0002461455090000041
AIC*(i*) Is a serial number i*Is determined by the minimum information criterion value of (c),
Figure FDA0002461455090000042
is composed of
Figure FDA0002461455090000043
The variance of (a) is determined,
Figure FDA0002461455090000044
is a serial number i*Corresponding active power and serial number i*-1 corresponding difference in active power,
Figure FDA0002461455090000045
is composed of
Figure FDA0002461455090000046
The variance of (a) is determined,
Figure FDA0002461455090000047
a fourth minimum information criterion calculation unit for calculating
Figure FDA0002461455090000048
The calculation is carried out by adopting a minimum information criterion, and the specific formula is as follows:
Figure FDA0002461455090000049
wherein,
Figure FDA00024614550900000410
is composed of
Figure FDA00024614550900000411
Variance of (1), Δ P0＝0。
CN201811097516.7A 2018-09-20 2018-09-20 Load switch event detection method and system based on minimum information criterion Expired - Fee Related CN108872665B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811097516.7A CN108872665B (en) 2018-09-20 2018-09-20 Load switch event detection method and system based on minimum information criterion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811097516.7A CN108872665B (en) 2018-09-20 2018-09-20 Load switch event detection method and system based on minimum information criterion

Publications (2)

Publication Number Publication Date
CN108872665A CN108872665A (en) 2018-11-23
CN108872665B true CN108872665B (en) 2020-07-03

Family

ID=64324406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811097516.7A Expired - Fee Related CN108872665B (en) 2018-09-20 2018-09-20 Load switch event detection method and system based on minimum information criterion

Country Status (1)

Country Link
CN (1) CN108872665B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110196354B (en) * 2019-04-23 2021-09-17 广东石油化工学院 Method and device for detecting switching event of load
CN110221119B (en) * 2019-06-25 2021-04-06 广东石油化工学院 Load switch event detection method and system based on power and akie fusion information
CN110244115B (en) * 2019-07-17 2020-12-04 华北电力大学 Load switch event detection method and system based on signal connectivity

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1731446A (en) * 2005-09-09 2006-02-08 清华大学 Comprehensive load prediction method based on minimum information loss
CN103513128A (en) * 2012-06-27 2014-01-15 财团法人工业技术研究院 Method and system for detecting load of electric appliance
KR20140098971A (en) * 2013-02-01 2014-08-11 전자부품연구원 Apparatus and Method for Test
CN104483575A (en) * 2014-12-22 2015-04-01 天津求实智源科技有限公司 Self-adaptive load event detection method for noninvasive power monitoring
CN106127303A (en) * 2016-06-15 2016-11-16 国网山东省电力公司菏泽供电公司 A kind of short-term load forecasting method towards multi-source data
CN106600074A (en) * 2016-12-28 2017-04-26 天津求实智源科技有限公司 DFHSMM-based non-intrusion type electric power load monitoring method and system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1731446A (en) * 2005-09-09 2006-02-08 清华大学 Comprehensive load prediction method based on minimum information loss
CN103513128A (en) * 2012-06-27 2014-01-15 财团法人工业技术研究院 Method and system for detecting load of electric appliance
KR20140098971A (en) * 2013-02-01 2014-08-11 전자부품연구원 Apparatus and Method for Test
CN104483575A (en) * 2014-12-22 2015-04-01 天津求实智源科技有限公司 Self-adaptive load event detection method for noninvasive power monitoring
CN106127303A (en) * 2016-06-15 2016-11-16 国网山东省电力公司菏泽供电公司 A kind of short-term load forecasting method towards multi-source data
CN106600074A (en) * 2016-12-28 2017-04-26 天津求实智源科技有限公司 DFHSMM-based non-intrusion type electric power load monitoring method and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A Fully Unsupervised Non-intrusive Load Monitoring Framework;Ruoxi Jia 等;《2015 IEEE International Conference on Smart Grid Communications (SmartGridComm): Data Management, Grid Analytics, and Dynamic Pricing》;20151231;872-878 *
Supervised Household’s Loads Pattern Recognition;Maher Azaza 等;《2016 IEEE Electrical Power and Energy Conference》;20161231;1-5 *
基于贝叶斯准则的非侵入式负荷监测方法;周晨轶 等;《浙江电力》;20180525;第37卷(第5期);7-11 *

Also Published As

Publication number Publication date
CN108872665A (en) 2018-11-23

Similar Documents

Publication Publication Date Title
CN109324285B (en) Load switch event detection method and system
CN108872666B (en) Load switch event detection method and system
CN108872665B (en) Load switch event detection method and system based on minimum information criterion
CN109085423B (en) Load switch event detection method and system
CN102279342B (en) Cage asynchronous motor rotor broken-bar fault detection method based on multi-signal classification technology (MUSIC) and simulation annealing algorithm (SAA)
CN109307798A (en) A kind of power signal filtering method for switch events detection
CN110460305A (en) Fault detection method, device, readable medium and the electronic equipment of photovoltaic array
CN109239592A (en) Contact of breaker performance monitoring method, device and monitoring system
WO2022233153A1 (en) Arc fault detection method and apparatus, and direct-current electric appliance
CN111664934A (en) Transformer state vibration and sound detection signal filtering method and system using feature selection
CN117783955A (en) Transformer winding state detection method, device and equipment
CN110244115B (en) Load switch event detection method and system based on signal connectivity
CN110196354B (en) Method and device for detecting switching event of load
CN110221119B (en) Load switch event detection method and system based on power and akie fusion information
CN111830405A (en) Load switch event detection method and system by using frequency difference
CN112434567B (en) Power signal filtering method and system by using noise jitter property
CN112307986B (en) Load switch event detection method and system by utilizing Gaussian gradient
CN113805050A (en) Phase selection closing angle monitoring method and device, computer equipment and storage medium
CN110646691B (en) Transformer vibration sound signal filtering method and system by utilizing stretching transformation
CN110286289B (en) Filtering method for vibration and sound detection signal of transformer
CN110702981A (en) Load switch event detection method and system using classification tree
CN112257576B (en) Load switch event detection method and system using Maha distance measure
CN112180152A (en) Load switch event detection method and system by means of mean shift clustering
CN112180155A (en) Load switch event detection method and system using tight support set
CN112200061A (en) Power signal reconstruction method and system by using first-order smooth model

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200703

Termination date: 20210920