CN114236301B - Non-invasive industrial load event detection method and system - Google Patents

Non-invasive industrial load event detection method and system Download PDF

Info

Publication number
CN114236301B
CN114236301B CN202210183139.9A CN202210183139A CN114236301B CN 114236301 B CN114236301 B CN 114236301B CN 202210183139 A CN202210183139 A CN 202210183139A CN 114236301 B CN114236301 B CN 114236301B
Authority
CN
China
Prior art keywords
load
data
adjacent
event
data segment
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
CN202210183139.9A
Other languages
Chinese (zh)
Other versions
CN114236301A (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202210183139.9A priority Critical patent/CN114236301B/en
Publication of CN114236301A publication Critical patent/CN114236301A/en
Application granted granted Critical
Publication of CN114236301B publication Critical patent/CN114236301B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R21/00Arrangements for measuring electric power or power factor
    • G01R21/06Arrangements for measuring electric power or power factor by measuring current and voltage

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a non-invasive industrial load event detection method and a system, which relate to the technical field of intelligent power consumption management, and the method comprises the following steps: synchronously acquiring electricity utilization data through a non-invasive acquisition device to form a first electricity utilization data set; intercepting N adjacent first data segments from the first power data set, and calculating the average power difference between any two adjacent first data segments; determining whether a load state change event occurs according to the average power difference value corresponding to a plurality of adjacent two first data segments; constructing a second electrical data set according to the load state change event obtained by detection; intercepting M adjacent second data segments from the second electrical data set, and calculating the intersection point number difference between any two adjacent second data segments; and determining whether a load mode change event occurs according to the difference value of the number of the intersection points between a plurality of adjacent two second data segments, thereby realizing the type detection of the industrial load event.

Description

Non-invasive industrial load event detection method and system
Technical Field
The invention relates to the technical field of intelligent power consumption management, in particular to a non-invasive industrial load event detection method and system.
Background
The targets of carbon peak reaching and carbon neutralization provide new requirements for low carbon and energy saving for industrial development, and industrial users are promoted to continuously improve the production energy efficiency level. The industrial production consumes huge electricity, and the electricity consumption data of the load can reflect the actual situation of the production operation of the factory to a great extent. Therefore, in industrial production, it is necessary to monitor load electricity consumption information, further obtain the production condition of a user and restore the load action, so that the industrial user can be helped to improve the production efficiency and realize low-carbon transformation.
The traditional intrusive load detection method directly collects independent data of the operation load and then carries out detection based on the collected data, although the obtained load information is real and reliable, the required equipment quantity is large, the cost is high, the influence of installation and maintenance on load operation is large, and the traditional intrusive load detection method is not suitable for industrial load event detection.
Disclosure of Invention
The invention aims to provide a non-invasive industrial load event detection method and a non-invasive industrial load event detection system, so as to effectively identify the type of a load event in industrial production.
To achieve the above object, the present invention provides a non-intrusive industrial load event detection method, comprising:
synchronously acquiring electricity utilization data from an electric power metering position of an industrial user through a non-invasive acquisition device to form a first electricity utilization data set;
intercepting N adjacent first data segments from the first power data set, and calculating the average power difference between any two adjacent first data segments; wherein N is a positive integer greater than 2;
determining whether a load state change event occurs according to the average power difference value corresponding to a plurality of adjacent two first data segments;
constructing a second electrical data set according to the load state change event obtained by detection;
intercepting M adjacent second data segments from the second electrical data set, and calculating the intersection number difference between any two adjacent second data segments; wherein M is a positive integer greater than 2;
and determining whether a load mode change event occurs according to the difference of the number of the intersection points between a plurality of adjacent two second data segments.
Optionally, the method further comprises:
and determining the load action corresponding to the load state change event and/or the load mode change event.
Optionally, the intercepting N adjacent first data segments from the first electrical data set, and calculating an average power difference between any two adjacent first data segments specifically includes:
taking a collecting point corresponding to the change of the voltage value from negative to positive as an initial point of a voltage cycle;
intercepting N adjacent first data segments from the first electrical data set by taking an initial point as a critical point; each of the first data segments includesN 1Electricity consumption data collected in each voltage period; wherein,N 1is a positive integer greater than 1;
calculating the load power consumption of each voltage period in each first data segment according to the power consumption data acquired by each voltage period in each first data segment;
calculating the average power consumption value of the load in each first data segment according to the power consumption of the load in each voltage period in each first data segment;
and calculating the average power difference value between two adjacent first data segments according to the load average power utilization values of the two adjacent first data segments.
Optionally, the determining whether a load state change event occurs according to the average power difference corresponding to the two adjacent first data segments specifically includes:
step S31: calculating load power consumption variance corresponding to each first data segment according to the load power consumption of each voltage period in each first data segment and the load average power consumption value in each first data segment;
step S32: judging the average power difference deltaP X Whether or not greater than
Figure 100002_DEST_PATH_IMAGE001
(ii) a If ΔP X Is greater than
Figure 831791DEST_PATH_IMAGE001
If yes, the event is a load state change event, and step S33 is executed; if ΔP X Is less than or equal to
Figure 709748DEST_PATH_IMAGE001
If yes, let X = X +1, determine whether X is less than or equal to N, if X is less than or equal to N, execute "step S32", if X is greater than N, end; wherein,
Figure 420084DEST_PATH_IMAGE002
are all rational numbers greater than 0 and less than 1,P' X the power consumption variance of the load corresponding to the Xth first data segment is shown,P X-1representing the average power consumption value of the load in the X-1 th first data segment;
step S33: it is determined whether the first event end condition is satisfied, and if the first event end condition is satisfied, it indicates that the load-state variation event has ended, and if the first event end condition is not satisfied, the process proceeds to "step S33".
Optionally, the first event ending condition is specifically:
Figure 100002_DEST_PATH_IMAGE003
wherein, X1And X2Indicating two non-adjacent first data segment sequence numbers, X1Less than X2P' nThe power consumption variance of the load corresponding to the nth first data segment is shown,P n-1represents the average power consumption value of the load in the (n-1) th first data segment,K 1denotes a set integer,. DELTAP nRepresents the average power difference between the nth first data segment and the (n-1) th first data segment.
Optionally, the determining whether a load pattern change event occurs according to a difference between the number of intersections between two adjacent second data segments specifically includes:
step S61: the difference value delta of the number of the intersection points is judgedS YWhether or not greater than
Figure 974562DEST_PATH_IMAGE004
(ii) a If the number of intersections is different by ΔS YIs greater than
Figure 100002_DEST_PATH_IMAGE005
If yes, the event is a load mode change event, and step S62 is executed; if the number of intersections is different by ΔS YIs less than or equal to
Figure 75242DEST_PATH_IMAGE005
If yes, let Y = Y +1, determine whether Y is less than or equal to M, if Y is less than or equal to M, return to "step S61", if Y is greater than M, end; wherein,
Figure 543176DEST_PATH_IMAGE006
is a rational number greater than 0 and less than 1,S Y-1the number of the intersection points corresponding to the Y-1 th second data segment is represented;
step S62: judging whether a second event ending condition is met, if so, indicating that the load mode change event is ended; if the second event end condition is not satisfied, "step S62" is continuously performed.
Optionally, the second event ending condition is specifically:
Figure 100002_DEST_PATH_IMAGE007
wherein, Y1And Y2Indicating two non-adjacent second data segment sequence numbers, Y1Less than Y2S n-1Indicates the number of the intersections corresponding to the (n-1) th second data segment,K 2 Denotes a set integer,. DELTAS nIndicating the difference in the number of intersections between the nth second data segment and the (n-1) th second data segment.
The present invention also provides a non-intrusive industrial load event detection system, the system comprising:
the first electricity data set construction module is used for synchronously acquiring electricity data from the electric power metering position of an industrial user through a non-invasive acquisition device to form a first electricity data set;
the average power difference calculation module is used for intercepting N adjacent first data segments from the first power data set and calculating the average power difference between any two adjacent first data segments; wherein N is a positive integer greater than 2;
the load state change event determining module is used for determining whether a load state change event occurs according to the average power difference value corresponding to the two adjacent first data segments;
the second electrical data set building module is used for building a second electrical data set according to the load state change event obtained by detection;
the intersection number difference calculation module is used for intercepting M adjacent second data segments from the second electrical data set and calculating the intersection number difference between any two adjacent second data segments; wherein M is a positive integer greater than 2;
and the load mode change event determining module is used for determining whether a load mode change event occurs according to the difference value of the number of the intersection points between the plurality of adjacent second data segments.
Optionally, the average power difference calculation module specifically includes:
the initial point determining unit is used for taking the acquisition point corresponding to the voltage value from negative to positive as the initial point of the voltage period;
the first interception unit is used for intercepting N adjacent first data sections from the first electric data set by taking an initial point as a critical point; each of the first data segments includesN 1Electricity consumption data collected in each voltage period; wherein N is1Is a positive integer greater than 1;
the first load electricity consumption power calculation unit is used for calculating the load electricity consumption power of each voltage period in each first data segment according to the electricity consumption data acquired by each voltage period in each first data segment;
the first load average power consumption value calculation unit is used for calculating the load average power consumption value in each first data segment according to the load power consumption of each voltage period in each first data segment;
and the average power difference calculation unit is used for calculating the average power difference between the two adjacent first data segments according to the load average power utilization values of the two adjacent first data segments.
Optionally, the load state change event determining module specifically includes:
the load electricity consumption power variance calculating unit is used for calculating load electricity consumption power variances corresponding to the first data sections according to the load electricity consumption power of each voltage period in each first data section and the load average electricity consumption power value in each first data section;
a first judging unit for judging the average power difference DeltaP X Whether or not greater than
Figure 780122DEST_PATH_IMAGE001
(ii) a If ΔP X Is greater than
Figure 974343DEST_PATH_IMAGE001
If yes, the event is a load state change event, and a second judgment unit is executed; if ΔP X Is less than or equal to
Figure 359057DEST_PATH_IMAGE001
If yes, let X = X +1, determine whether X is less than or equal to N, if X is less than or equal to N, execute "first determination unit", if X is greater than N, end; wherein,
Figure 641134DEST_PATH_IMAGE008
are all rational numbers greater than 0 and less than 1,P' X the power consumption variance of the load corresponding to the Xth first data segment is shown,P X-1 representing the average power consumption value of the load in the X-1 th first data segment;
and a second judging unit for judging whether the first event ending condition is satisfied, if the first event ending condition is satisfied, indicating that the load state change event has ended, and if the first event ending condition is not satisfied, continuing to execute the second judging unit.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
according to the invention, the non-invasive acquisition device is used for synchronously acquiring the electricity consumption data from the electric power metering position of the industrial user, and then the type of the industrial load event is effectively detected based on the acquired electricity consumption data, so that the problems of complex data acquisition and high cost in the traditional method are solved, and the industrial load event can be detected.
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 flow chart of a non-intrusive industrial load event detection method of the present invention;
FIG. 2 is a block diagram of a non-intrusive industrial load event detection system of the present invention;
FIG. 3 is a schematic diagram of a non-intrusive industrial load event detection architecture in accordance with the present invention;
FIG. 4 shows the detection of the present inventiont 1The state change event at the moment and its extracted current and voltage waveforms;
FIG. 5 shows the detection of the present inventiont 2The pattern change event at the time and the power waveform extracted therefrom.
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 non-invasive industrial load event detection method and a non-invasive industrial load event detection system, so as to effectively identify the type of a load event in industrial production.
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.
Example 1
In the following embodiments, in order to distinguish the difference between the data pieces in step S2 and step S5, the data piece intercepted in step S2 is referred to as a first data piece, and the data piece intercepted in step S5 is referred to as a second data piece. In addition, the load events in the present invention include load state change events and load pattern change events; recording the event of the slow and large change of the power signal caused by the load switching action as a load state change event; and recording the event of regular change of the power signal with unchanged power before and after caused by the load change working mode as a load mode change event.
As shown in fig. 1, the present invention discloses a non-intrusive industrial load event detection method, comprising:
step S1: and synchronously acquiring electricity utilization data from the electric power metering position of the industrial user through a non-invasive acquisition device to form a first electricity utilization data set.
Step S2: intercepting N adjacent first data segments from the first power data set, and calculating the average power difference between any two adjacent first data segments; wherein N is a positive integer greater than 2.
Step S3: and determining whether a load state change event occurs according to the average power difference value corresponding to the plurality of adjacent two first data segments.
Step S4: and constructing a second electrical data set according to the detected load state change event.
Step S5: intercepting M adjacent second data segments from the second electrical data set, and calculating the intersection number difference between any two adjacent second data segments; wherein M is a positive integer greater than 2.
Step S6: and determining whether a load mode change event occurs according to the difference of the number of the intersection points between a plurality of adjacent two second data segments.
The individual steps are discussed in detail below:
step S1: synchronously acquiring electricity utilization data from an electric power metering position of an industrial user through a non-invasive acquisition device to form a first electricity utilization data set; the electricity consumption data includes: voltage and current. In this embodiment, the power consumption data is collected by a high-frequency collection method.
Step S2: intercepting N adjacent first data segments from the first power data set, and calculating an average power difference value between any two adjacent first data segments, specifically comprising:
step S21: taking a collecting point corresponding to the voltage value from negative to positive as an initial point of a voltage cycle, wherein a specific judgment formula is as follows:
Figure DEST_PATH_IMAGE009
(1)
wherein,U A,z representing the voltage collected at the z-th sample point for phase a.
Step S22: intercepting N adjacent first data segments from the first electrical data set by taking an initial point as a critical point; each of the first data segments includesN 1Electricity consumption data collected in each voltage period; in the present embodiment, the first and second electrodes are,N 1and 500 is taken.
Step S23: calculating the load power consumption of each voltage period in each first data segment according to the power consumption data acquired by each voltage period in each first data segment, wherein the specific calculation formula is as follows:
Figure 919537DEST_PATH_IMAGE010
(2)
wherein,P X,i the power consumption of the load of the ith voltage period in the Xth first data segment is represented,Tthe total number of points collected for each voltage cycle,U A,z representing the voltage acquired at the z-th sample point of phase a,I A,z represents the current collected by the z-th sampling point of the phase A, and X belongs to [1, N ]],i∈[1,N 1]。
Step S24: calculating the average power consumption value of the load in each first data segment according to the power consumption of the load in each voltage period in each first data segment, wherein the specific calculation formula is as follows:
Figure DEST_PATH_IMAGE011
(3)
wherein,P X represents the average power consumption value of the load in the Xth first data segment,N 1representing the total number of voltage cycles contained in each first data segment,P X,i and representing the electric power used by the load in the ith voltage period in the Xth first data segment.
Step S25: calculating the average power difference value between two adjacent first data segments according to the load average power utilization values of the two adjacent first data segments, wherein the specific calculation formula is as follows:
Figure 753501DEST_PATH_IMAGE012
(4)
wherein,P X represents the average power utilization value of the load in the Xth first data segment,P X-1 represents the average power consumption value of the load in the X-1 th first data segment, deltaP X Represents the average power difference between the xth first data segment and the xth-1 first data segment.
Step S3: determining whether a load state change event occurs according to an average power difference value corresponding to a plurality of adjacent two first data segments, specifically including:
step S31: calculating the load power consumption variance corresponding to each first data segment according to the load power consumption of each voltage period in each first data segment and the load average power consumption value in each first data segment, wherein the specific formula is as follows:
Figure DEST_PATH_IMAGE013
(5)
wherein,P'xrepresents the power variance of the load corresponding to the Xth first data segment,P X represents the average power consumption value of the load in the Xth first data segment,N 1which represents the total number of voltage cycles,P X,i and representing the electric power used by the load in the ith voltage period in the Xth first data segment.
Step S32: judging the average power difference deltaP X Whether or not greater than
Figure 947547DEST_PATH_IMAGE014
(ii) a If ΔP X Is greater than
Figure DEST_PATH_IMAGE015
If yes, the event is a load state change event, andexecute "step S33"; if ΔP X Is less than or equal to
Figure 689107DEST_PATH_IMAGE015
If yes, let X = X +1, determine whether X is less than or equal to N, if X is less than or equal to N, execute "step S32", if X is greater than N, end; wherein,
Figure 228542DEST_PATH_IMAGE016
are all rational numbers greater than 0 and less than 1,P' X the power consumption variance of the load corresponding to the Xth first data segment is shown,P X-1 representing the average power consumption value of the load in the X-1 th first data segment; in the present embodiment, the first and second electrodes are,
Figure DEST_PATH_IMAGE017
the content of the organic acid was 0.02,
Figure 295724DEST_PATH_IMAGE018
is 0.3.
Step S33: it is determined whether the first event end condition is satisfied, and if the first event end condition is satisfied, it indicates that the load-state variation event has ended, and if the first event end condition is not satisfied, the process proceeds to "step S33".
The first event ending condition is specifically as follows:
Figure 671342DEST_PATH_IMAGE003
(6)
wherein, X1And X2Indicating two non-adjacent first data segment sequence numbers, X1Less than X2
Figure 278909DEST_PATH_IMAGE016
Are all rational numbers greater than 0 and less than 1,P' nrepresents the power variance, P, of the load power corresponding to the nth first data segmentn-1Represents the average power consumption value of the load in the (n-1) th first data segment,K 1denotes a set integer,. DELTAP nRepresents the average power difference between the nth first data segment and the (n-1) th first data segment. In the present embodiment, the first and second electrodes are,K 1is 20.
Step S4: and constructing a second electrical data set according to the detected load state change event. The second electrical data set is current data and voltage data between two adjacent load-state-change events.
Step S5: intercepting M adjacent second data segments from the second electrical data set, and calculating the intersection number difference between any two adjacent second data segments, specifically comprising:
step S51: intercepting M adjacent second data segments from the second electrical data set by taking the initial point as a critical point; each of the second data segments includesN 2Electricity consumption data collected in each voltage period; in this exampleN 2And 300 is taken.
Step S52: and calculating the load power consumption of each voltage period in each second data segment according to the power consumption data acquired by each voltage period in each second data segment, wherein the specific calculation formula is as follows:
Figure DEST_PATH_IMAGE019
(7)
wherein,P Y,j the electric power of the load of the jth voltage period in the Yth second data segment is represented,Tfor the total number of points collected for each voltage cycle,U A,z representing the voltage acquired at the z-th sample point of phase a,I A,z represents the current collected by the z-th sampling point of the phase A, and Y belongs to [1, M ]],j∈[1,N 2]。
Step S53: calculating the average power consumption value of the load in each second data segment according to the power consumption of the load in each voltage period in each second data segment, wherein the specific calculation formula is as follows:
Figure 876113DEST_PATH_IMAGE020
(8)
wherein,P Y representing the load average power consumption value in the Yth second data segment,N 2indicating the total number of voltage cycles contained within each second data segment,P Y,j and the electric power consumption of the load in the jth voltage period in the Yth second data segment is shown.
Step S54: calculating the number of intersection points corresponding to each second data segment according to the load power consumption of each voltage period in each second data segment and the load average power consumption value in each second data segment, wherein the specific formula is as follows:
Figure DEST_PATH_IMAGE021
(9)
wherein,
Figure 379775DEST_PATH_IMAGE022
it is indicated that the set-up parameters,P Yrepresenting the load average power consumption value in the Yth second data segment,N 2representing the total number of voltage periods contained in each second data segment,P Y,j the electric power of the load of the jth voltage period in the Yth second data segment is represented,S Ythe number of the intersection points corresponding to the Y-th second data segment is shown, in this embodiment, the number of the intersection points is the number of the intersection points of the load power curve and the power average value,
Figure 242689DEST_PATH_IMAGE022
take 5W.
Step S55: calculating the difference value of the number of the intersections between two adjacent second data segments according to the number of the intersections corresponding to the two adjacent second data segments, wherein the specific formula is as follows:
S Y=S Y-S Y-1(10)
wherein, ΔS YIndicating the difference in the number of intersections between the Y-th second data segment and the Y-1 st second data segment,S Yindicating the number of intersections corresponding to the Y-th second data segment.
Step S6: determining whether a load mode change event occurs according to a difference of the number of intersections between a plurality of adjacent two second data segments, specifically comprising:
step S61: judging the difference value delta of the number of the intersection pointsS YWhether or not greater than
Figure 653948DEST_PATH_IMAGE005
(ii) a If the number of intersections is different by ΔS YIs greater than
Figure 918707DEST_PATH_IMAGE005
If yes, the event is a load mode change event, and step S62 is executed; if the number of intersections is different by ΔS YIs less than or equal to
Figure 661447DEST_PATH_IMAGE005
If yes, let Y = Y +1, determine whether Y is less than or equal to M, if Y is less than or equal to M, return to "step S61", if Y is greater than M, end; wherein,
Figure 11657DEST_PATH_IMAGE006
is a rational number greater than 0 and less than 1,S Y-1indicates the number of intersections, Delta, corresponding to the Y-1 th second data segmentS YRepresenting the difference of the number of the intersection points between the Yth second data segment and the Y-1 th second data segment; in the present embodiment of the present invention,
Figure 242918DEST_PATH_IMAGE006
take 0.5.
Step S62: judging whether a second event ending condition is met, if so, indicating that the load mode change event is ended; if the second event end condition is not satisfied, "step S62" is continuously performed.
The second event end condition is specifically:
Figure DEST_PATH_IMAGE023
(11)
wherein, Y1And Y2Indicating two non-adjacent second data segment sequence numbers, Y1Less than Y2
Figure 532823DEST_PATH_IMAGE006
Is a rational number greater than 0 and less than 1,S n-1indicates the number of the intersections corresponding to the (n-1) th second data segment,K 2 denotes a set integer,. DELTAS nIndicating the difference in the number of intersections between the nth second data segment and the (n-1) th second data segment. In the present embodiment, the first and second electrodes are,
Figure 925758DEST_PATH_IMAGE006
taking out the mixture of 0.5 percent,K 2and 5, taking.
Step S7: determining a load action corresponding to the load state change event and/or the load mode change event, specifically comprising:
step S71: determining a load action corresponding to the load state change event, specifically comprising:
step S711: and extracting a steady-state current waveform before the load state change event, a steady-state current waveform after the load state change event and a voltage waveform corresponding to each steady-state current waveform to obtain a steady-state waveform of the event.
The specific decision formula of the steady-state current waveform is as follows:
Figure 746953DEST_PATH_IMAGE024
(12)
wherein, X1And X2Indicating two non-adjacent first data segment sequence numbers, X1Less than X2
Figure DEST_PATH_IMAGE025
Indicating a set steady-state current waveform decision threshold,K 3denotes a set integer,. DELTAP nRepresenting the average power difference between the nth first data segment and the (n-1) th first data segment; in the present embodiment, the first and second electrodes are,
Figure 719588DEST_PATH_IMAGE025
is a 500-watt water-cooling water heater,K 3is 2.
The specific formula of the steady state current waveform and voltage waveform of the event is as follows:
Figure 473786DEST_PATH_IMAGE026
(13)
Figure 37623DEST_PATH_IMAGE027
(14)
wherein,I o,r a steady state current waveform indicative of a detected state change event,I r representing a steady state current waveform following the load condition varying event,I r-1 representing a steady state current waveform prior to the load condition alteration event,U o,r a voltage waveform indicative of the detected state change event,U r a voltage waveform indicative of a post-load condition alteration event,U r-1 a voltage waveform representing a voltage waveform prior to the load condition alteration event.
Step S712: determining the electrical characteristics corresponding to the load state change event according to the steady-state waveform of the event corresponding to the load state change event, wherein the specific formula is as follows:
Figure 628004DEST_PATH_IMAGE028
(15)
Figure DEST_PATH_IMAGE029
(16)
wherein,I s for the effective value of the event current,I o,r,z as an event currentI o,r The sampled current at the z-th sample point,P p in order to be the power of the event,U o,r,z is an event voltageU o,r The sampled voltage at the z-th sample point,Tfor the sampling period, in this embodimentTIs 100.
Step S713: and comparing the electrical characteristics corresponding to the load state change event with reference information, and determining the load action corresponding to the load state change event. The reference information is at least one of a parameter of the load device and production control information.
Step S72: determining a load action corresponding to the load mode change event, specifically comprising:
step S721: and extracting voltage and current waveforms of the whole process of the load mode variation event to obtain a transient waveform of the event.
The specific formula of the transient waveform of an event is:
Figure 388018DEST_PATH_IMAGE030
(17)
wherein,P v,A,i P v,B,i P v,C,i respectively representing the phase A, B and C during the load mode change eventiThe value of the extracted power for each voltage period,i∈[1,TN2L];I v,z,A I v,z,B I v,z,C respectively representing the A, B and C phases during the load mode change eventiThe z-th sampling current of each voltage period;U v,z,A U v,z,B U v,z,C respectively representing the A, B and C phases during the load mode change eventiThe z-th sampling voltage of each voltage period; t represents the number of sampling points in the voltage period; l represents the number of the second data segment in which the load pattern variation event continues;I v,h,A I v,h,B I v,h,C respectively representing the h-th sampling current of the phases A, B and C in the process of the load mode change event;U v,h,A U v,h,B U v,h,C respectively representing the h-th sampling voltages of the phases A, B and C in the process of the load mode change event; z is equal to [1, T ]]Indicating phase A, B, C during load mode change eventiThe z-th sampling point within a voltage period;h∈[1,TN2L]representing the h-th sampling point of A, B and C phases in the process of the load mode change event; in this embodiment, T is 100.
Step S722: determining the time characteristics corresponding to the load mode change event according to the transient waveform of the load mode change event, wherein the specific formula is as follows:
T v =TN 2 L(18)
wherein,T v indicating the duration of the load pattern change event,Tthe number of sampling points representing the voltage period,Lthe number of the second data segment indicating the continuation of the load pattern variation event; in the present embodiment, the first and second electrodes are, T is 100.
Step S723: and comparing the time characteristics corresponding to the load mode change event with reference information, and determining the load action corresponding to the load mode change event. The reference information is at least one of a parameter of the load device and production control information.
Example 2
As shown in fig. 2, the present invention also provides a non-intrusive industrial load event detection system, the system comprising:
the first electricity data set constructing module 201 is configured to synchronously acquire electricity data from an electric power metering place of an industrial user through a non-invasive acquisition device, so as to form a first electricity data set.
An average power difference calculation module 202, configured to intercept N adjacent first data segments from the first electrical data set, and calculate an average power difference between any two adjacent first data segments; wherein N is a positive integer greater than 2.
The load state change event determining module 203 is configured to determine whether a load state change event occurs according to an average power difference value corresponding to a plurality of adjacent two first data segments.
And a second electrical data set constructing module 204, configured to construct a second electrical data set according to the detected load state change event.
An intersection number difference calculation module 205, configured to intercept M adjacent second data segments from the second electrical data set, and calculate an intersection number difference between any two adjacent second data segments; wherein M is a positive integer greater than 2.
The load pattern change event determining module 206 is configured to determine whether a load pattern change event occurs according to a difference between the number of intersections between two adjacent second data segments.
As an optional implementation manner, the average power difference calculation module 202 according to the present invention specifically includes:
and the initial point determining unit is used for taking the acquisition point corresponding to the voltage value changed from negative to positive as the initial point of the voltage period.
The first interception unit is used for intercepting N adjacent first data sections from the first electric data set by taking an initial point as a critical point; each of the first data segments includesN 1Electricity consumption data collected in each voltage period; wherein,N 1is a positive integer greater than 1.
And the first load electricity power calculation unit is used for calculating the load electricity power of each voltage period in each first data segment according to the electricity data acquired by each voltage period in each first data segment.
And the first load average power consumption power value calculating unit is used for calculating the load average power consumption power value in each first data segment according to the load power consumption power of each voltage period in each first data segment.
And the average power difference calculation unit is used for calculating the average power difference between the two adjacent first data segments according to the load average power utilization values of the two adjacent first data segments.
As an optional implementation manner, the load state change event determining module 203 of the present invention specifically includes:
the load electricity consumption power variance calculating unit is used for calculating load electricity consumption power variances corresponding to the first data sections according to the load electricity consumption power of each voltage period in each first data section and the load average electricity consumption power value in each first data section;
a first judging unit for judging the average power difference DeltaP X Whether or not greater than
Figure 731144DEST_PATH_IMAGE001
(ii) a If a is equalP X Is greater than
Figure 731461DEST_PATH_IMAGE001
If yes, the event is a load state change event, and a second judgment unit is executed; if ΔP X Is less than or equal to
Figure 543559DEST_PATH_IMAGE001
If yes, let X = X +1, determine whether X is less than or equal to N, if X is less than or equal to N, execute "first determination unit", if X is greater than N, end; wherein,
Figure 169581DEST_PATH_IMAGE008
are all rational numbers greater than 0 and less than 1,P' X represents the power variance of the load corresponding to the Xth first data segment,P X-1 representing the average power consumption value of the load in the X-1 th first data segment;
and a second judging unit for judging whether the first event ending condition is satisfied, if the first event ending condition is satisfied, indicating that the load state change event has ended, and if the first event ending condition is not satisfied, continuing to execute the second judging unit.
As an optional implementation manner, the load mode change event determining module 206 of the present invention specifically includes:
a third judging unit for judging the difference value Delta of the number of the intersection pointsS YWhether or not greater than
Figure 383525DEST_PATH_IMAGE005
(ii) a If the number of intersections is different ΔS YIs greater than
Figure 556010DEST_PATH_IMAGE005
If yes, the event is a load mode change event, and a fourth judgment unit is executed; if the number of intersections is different by ΔS YIs less than or equal to
Figure 120983DEST_PATH_IMAGE005
If Y = Y +1, judging whether Y is less than or equal to M, if Y is less than or equal to M, returning to a third judgment unit, and if Y is greater than M, ending; wherein,
Figure 35850DEST_PATH_IMAGE006
is a rational number greater than 0 and less than 1,S Y-1indicates the number of intersections, Delta, corresponding to the Y-1 th second data segmentS YIndicating the difference in the number of intersections between the Y-th second data segment and the Y-1 st second data segment.
A fourth judging unit, configured to judge whether a second event end condition is satisfied, and if the second event end condition is satisfied, the fourth judging unit indicates that the load pattern change event has ended; if the second event end condition is not satisfied, the "fourth judgment unit" is continuously executed.
As an optional implementation, the system of the present invention further includes:
and the load action determining module is used for determining the load action corresponding to the load state change event and/or the load mode change event.
As an optional implementation manner, the load action determining module of the present invention specifically includes:
and a first load operation determination unit for determining a load operation corresponding to the load state change event.
And a second load operation determination unit for determining a load operation corresponding to the load pattern change event.
Example 3
As shown in FIG. 3, a non-invasive acquisition device is installed at the low-voltage side power metering position of a certain steel mill, voltage waveform and current waveform signal data are acquired, the sampling rate is 5000Hz, and the acquisition time is 60 minutes.
Experiments were performed according to the method disclosed in the present invention:
1) according to the type of the production process classification event of the user, taking the steel mill of the embodiment as an example, according to the production process, except for an iron-making blast furnace and a matched load thereof, the steel mill has the condition change of opening and closing, wherein loads such as a hot rolling mill, a cold rolling mill, a coiler, an uncoiler, a flying shear, a continuous casting machine and the like repeatedly change between the conditions of outputting force to the workpiece due to the existence of intervals among the workpieces, and a mode change event occurs.
2) The total current and voltage waveforms of the industrial user are collected.
3) And calculating the average power consumption of each first data segment.
4) Detection and extraction of State Change events by transient waveform features, FIG. 4 is a diagram of detection of State Change events according to the present inventiont 1The state change event at the moment and its extracted current and voltage waveforms.
5) The acquisition signal is segmented with a state change event.
6) Calculating the number of the intersection points of the power curve of each second data segment and the mean value transverse line thereofS Y
7) Detection of pattern change events by transient waveform characteristics, FIG. 5 illustrates detection of pattern change events according to the present inventiont 2The pattern change event at the time and the power waveform extracted therefrom.
8) Calculating the electrical and time characteristics of the event, comparing with the reference information in the embodiment 1, identifying the load type of the event, and restoring the load action.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
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 (10)

1. A method of non-intrusive industrial load event detection, the method comprising:
synchronously acquiring electricity utilization data from an electric power metering position of an industrial user through a non-invasive acquisition device to form a first electricity utilization data set;
intercepting N adjacent first data segments from the first power data set, and calculating the average power difference between any two adjacent first data segments; wherein N is a positive integer greater than 2;
determining whether a load state change event occurs according to the average power difference value corresponding to a plurality of adjacent two first data segments;
constructing a second electrical data set according to the load state change event obtained by detection;
intercepting M adjacent second data segments from the second electrical data set, and calculating the intersection number difference between any two adjacent second data segments; wherein M is a positive integer greater than 2;
and determining whether a load mode change event occurs according to the difference of the number of the intersection points between a plurality of adjacent two second data segments.
2. The method of non-intrusive industrial load event detection as defined in claim 1, further comprising:
and determining the load action corresponding to the load state change event and/or the load mode change event.
3. The method according to claim 1, wherein the step of extracting N adjacent first data segments from the first electrical data set and calculating the average power difference between any two adjacent first data segments comprises:
taking a collecting point corresponding to the change of the voltage value from negative to positive as an initial point of a voltage cycle;
intercepting N adjacent first data segments from the first electrical data set by taking an initial point as a critical point; each of the first data segments includesN 1Electricity consumption data collected in each voltage period; wherein,N 1is a positive integer greater than 1;
calculating the load power consumption of each voltage period in each first data segment according to the power consumption data acquired by each voltage period in each first data segment;
calculating the average power consumption value of the load in each first data segment according to the power consumption of the load in each voltage period in each first data segment;
and calculating the average power difference value between two adjacent first data segments according to the load average power utilization values of the two adjacent first data segments.
4. The method according to claim 3, wherein the determining whether the load state change event occurs according to the average power difference corresponding to two adjacent first data segments specifically comprises:
step S31: calculating load power consumption variance corresponding to each first data segment according to the load power consumption of each voltage period in each first data segment and the load average power consumption value in each first data segment;
step S32: judging the average power difference deltaP X Whether or not greater than
Figure DEST_PATH_IMAGE001
(ii) a If ΔP X Is greater than
Figure 170040DEST_PATH_IMAGE001
If yes, the event is a load state change event, and step S33 is executed; if ΔP X Is less than or equal to
Figure 352760DEST_PATH_IMAGE001
If yes, let X = X +1, determine whether X is less than or equal to N, if X is less than or equal to N, execute "step S32", if X is greater than N, end; wherein,
Figure DEST_PATH_IMAGE002
are all rational numbers greater than 0 and less than 1,P' X the power consumption variance of the load corresponding to the Xth first data segment is shown,P X-1 representing the average power consumption value of the load in the X-1 th first data segment;
step S33: it is determined whether the first event end condition is satisfied, and if the first event end condition is satisfied, it indicates that the load-state variation event has ended, and if the first event end condition is not satisfied, the process proceeds to "step S33".
5. The method according to claim 4, wherein the first event ending condition is:
Figure DEST_PATH_IMAGE003
wherein, X1And X2Indicating two non-adjacent first data segment sequence numbers, X1Less than X2P' nThe power consumption variance of the load corresponding to the nth first data segment is shown,P n-1represents the average power consumption value of the load in the (n-1) th first data segment,K 1denotes a set integer,. DELTAP nRepresents the average power difference between the nth first data segment and the (n-1) th first data segment.
6. The method according to claim 1, wherein the determining whether the load pattern change event occurs according to the difference between the number of intersections between two adjacent second data segments specifically comprises:
step S61: the difference value delta of the number of the intersection points is judgedS YWhether or not greater than
Figure DEST_PATH_IMAGE004
(ii) a If the number of intersections is different by ΔS YIs greater than
Figure 582097DEST_PATH_IMAGE004
If yes, the event is a load mode change event, and step S62 is executed; if the number of intersections is different by ΔS YIs less than or equal to
Figure 242886DEST_PATH_IMAGE004
If yes, let Y = Y +1, determine whether Y is less than or equal to M, if Y is less than or equal to M, return to "step S61", if Y is greater than M, end; wherein,
Figure DEST_PATH_IMAGE005
is a rational number greater than 0 and less than 1,S Y-1the number of the intersection points corresponding to the Y-1 th second data segment is represented;
step S62: judging whether a second event ending condition is met, if so, indicating that the load mode change event is ended; if the second event end condition is not satisfied, "step S62" is continuously performed.
7. The method according to claim 6, wherein the second event ending condition is:
Figure DEST_PATH_IMAGE006
wherein, Y1And Y2Indicating two non-adjacent second data segment sequence numbers, Y1Less than Y2S n-1Indicates the number of the intersection points corresponding to the (n-1) th second data segmentNumber of,K 2 Denotes a set integer,. DELTAS nIndicating the difference in the number of intersections between the nth second data segment and the (n-1) th second data segment.
8. A non-intrusive industrial load event detection system, the system comprising:
the first electricity data set construction module is used for synchronously acquiring electricity data from the electric power metering position of an industrial user through a non-invasive acquisition device to form a first electricity data set;
the average power difference calculation module is used for intercepting N adjacent first data segments from the first power data set and calculating the average power difference between any two adjacent first data segments; wherein N is a positive integer greater than 2;
the load state change event determining module is used for determining whether a load state change event occurs according to the average power difference value corresponding to the two adjacent first data segments;
the second electrical data set building module is used for building a second electrical data set according to the load state change event obtained by detection;
the intersection number difference calculation module is used for intercepting M adjacent second data segments from the second electrical data set and calculating the intersection number difference between any two adjacent second data segments; wherein M is a positive integer greater than 2;
and the load mode change event determining module is used for determining whether a load mode change event occurs according to the difference value of the number of the intersection points between the plurality of adjacent second data segments.
9. The system of claim 8, wherein the mean power difference calculation module comprises:
the initial point determining unit is used for taking the acquisition point corresponding to the voltage value from negative to positive as the initial point of the voltage cycle;
a first cut-off unit for cutting off the first light beam from the initial point as the critical pointIntercepting N adjacent first data segments on the electricity utilization data set; each of the first data segments includesN 1Electricity consumption data collected in each voltage period; wherein N is1Is a positive integer greater than 1;
the first load electricity consumption power calculation unit is used for calculating the load electricity consumption power of each voltage period in each first data segment according to the electricity consumption data acquired by each voltage period in each first data segment;
the first load average power consumption value calculation unit is used for calculating the load average power consumption value in each first data segment according to the load power consumption of each voltage period in each first data segment;
and the average power difference calculation unit is used for calculating the average power difference between the two adjacent first data segments according to the load average power utilization values of the two adjacent first data segments.
10. The system of claim 9, wherein the load-state-change event determination module comprises:
the load electricity consumption power variance calculating unit is used for calculating load electricity consumption power variances corresponding to the first data sections according to the load electricity consumption power of each voltage period in each first data section and the load average electricity consumption power value in each first data section;
a first judgment unit for judging the average power difference value DeltaP X Whether or not greater than
Figure 181892DEST_PATH_IMAGE001
(ii) a If ΔP X Is greater than
Figure 332250DEST_PATH_IMAGE001
If yes, the event is a load state change event, and a second judgment unit is executed; if a is equalP X Is less than or equal to
Figure 864863DEST_PATH_IMAGE001
If yes, let X = X +1, determine whether X is less than or equal to N, if X is less than or equal to N, execute "first determination unit", if X is greater than N, end; wherein,
Figure DEST_PATH_IMAGE007
are all rational numbers greater than 0 and less than 1,P' X the power consumption variance of the load corresponding to the Xth first data segment is shown,P X-1 representing the average power consumption value of the load in the X-1 th first data segment;
and a second judging unit for judging whether the first event ending condition is satisfied, if the first event ending condition is satisfied, indicating that the load state change event has ended, and if the first event ending condition is not satisfied, continuing to execute the second judging unit.
CN202210183139.9A 2022-02-28 2022-02-28 Non-invasive industrial load event detection method and system Active CN114236301B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210183139.9A CN114236301B (en) 2022-02-28 2022-02-28 Non-invasive industrial load event detection method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210183139.9A CN114236301B (en) 2022-02-28 2022-02-28 Non-invasive industrial load event detection method and system

Publications (2)

Publication Number Publication Date
CN114236301A CN114236301A (en) 2022-03-25
CN114236301B true CN114236301B (en) 2022-05-27

Family

ID=80748204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210183139.9A Active CN114236301B (en) 2022-02-28 2022-02-28 Non-invasive industrial load event detection method and system

Country Status (1)

Country Link
CN (1) CN114236301B (en)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140132523A (en) * 2013-05-08 2014-11-18 삼성전자주식회사 Non-intrusive load monitoring apparatus and method
CN107302217A (en) * 2017-06-02 2017-10-27 广东中粤售电科技有限公司 Non-intrusion type load decomposition method and system based on V I tracks
CN107525964B (en) * 2017-10-23 2023-11-21 云南电网有限责任公司电力科学研究院 Non-invasive load identification method and device based on fusion decision
CN108899892A (en) * 2018-06-08 2018-11-27 上海电力学院 One kind being based on CNN non-intrusion type power load decomposition method
CN110569877A (en) * 2019-08-07 2019-12-13 武汉中原电子信息有限公司 Non-invasive load identification method and device and computing equipment
CN111830347B (en) * 2020-07-17 2021-03-19 四川大学 Two-stage non-invasive load monitoring method based on event
CN111736012B (en) * 2020-08-20 2021-01-05 国网浙江省电力有限公司湖州供电公司 Non-invasive power load identification method and system

Also Published As

Publication number Publication date
CN114236301A (en) 2022-03-25

Similar Documents

Publication Publication Date Title
CN110956220B (en) Non-invasive household appliance load identification method
Su et al. Feature selection of non-intrusive load monitoring system using STFT and wavelet transform
CN107025365A (en) A kind of non-intruding load discrimination method for user side
Chandra et al. Online voltage stability and load margin assessment using wide area measurements
CN111812451B (en) Phase current transient fault component-based distributed line selection method for power distribution network
CN106932686A (en) Power network exception method for quick based on minimax method of characteristic
CN102135560A (en) Disturbance identification method used during intrusion of lightning waves in transformer substation
CN106443433A (en) Monitoring system and monitoring method for circuit breaker state based on divide-shut brake coil current
Li et al. Fast event identification through subspace characterization of PMU data in power systems
CN107797028A (en) A kind of distribution arc grounding fault recognition method
CN110187242A (en) A kind of power distribution network abnormality detection and localization method based on spatio-temporal correlation analysis
CN114236301B (en) Non-invasive industrial load event detection method and system
Yi et al. A new event-detection method based on composite windows in nilm for industrial settings
CN102445660A (en) Gray Verhulst model-based prediction method of power angle of generator
CN109615160A (en) CVT electric voltage exception data analysing method
CN118150942A (en) Distribution network current ground fault positioning method, device, equipment and medium
CN102298115B (en) Method for online monitoring and analyzing adjustment performance of excitation system based on phasor measurement technology
CN117471157A (en) Non-invasive electric bicycle load charging event detection and extraction method
Cui et al. Anode effect prediction method based on local effect detection
CN111814325B (en) Single-variable alarm threshold optimization method for non-stationary operation of coal-fired power generation equipment
CN111146779B (en) Flexible safety control method and system for large power grid with successive faults of power equipment
CN115905891B (en) Power distribution network operation mode and key influence factor identification method based on PMU data
Sun et al. A Non-Intrusive Load Identification Method Based on AP Clustering
CN104549758B (en) Electrostatic field voltage control method and system for electrostatic dust collector
CN118353172B (en) On-line monitoring and controlling system for temperature and humidity of switch cabinet

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