CN117060592A - Multi-domain cooperation-based power grid data synchronous calibration method and system - Google Patents
Multi-domain cooperation-based power grid data synchronous calibration method and system Download PDFInfo
- Publication number
- CN117060592A CN117060592A CN202311101913.8A CN202311101913A CN117060592A CN 117060592 A CN117060592 A CN 117060592A CN 202311101913 A CN202311101913 A CN 202311101913A CN 117060592 A CN117060592 A CN 117060592A
- Authority
- CN
- China
- Prior art keywords
- value
- data
- time
- correction step
- pmu
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 42
- 230000005540 biological transmission Effects 0.000 claims abstract description 131
- 238000012937 correction Methods 0.000 claims description 186
- 238000005259 measurement Methods 0.000 claims description 124
- 238000005457 optimization Methods 0.000 claims description 40
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000004891 communication Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000013500 data storage Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 235000008694 Humulus lupulus Nutrition 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J13/00—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
- H02J13/00001—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by the display of information or by user interaction, e.g. supervisory control and data acquisition systems [SCADA] or graphical user interfaces [GUI]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J13/00—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
- H02J13/00006—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L43/00—Arrangements for monitoring or testing data switching networks
- H04L43/08—Monitoring or testing based on specific metrics, e.g. QoS, energy consumption or environmental parameters
- H04L43/0852—Delays
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L67/00—Network arrangements or protocols for supporting network services or applications
- H04L67/01—Protocols
- H04L67/12—Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L7/00—Arrangements for synchronising receiver with transmitter
- H04L7/0016—Arrangements for synchronising receiver with transmitter correction of synchronization errors
- H04L7/0033—Correction by delay
Abstract
The invention discloses a multi-domain collaborative power grid data synchronous calibration method and system, which take network topology complexity, measured data uploading information quantity, signal to interference and noise ratio of a measured data transmission link and other information into consideration.
Description
Technical Field
The invention belongs to the field of novel power systems, and particularly relates to a multi-domain cooperation-based power grid data synchronous calibration method and system.
Background
With the development of a novel power system, a power grid operation mode of distributed resource regulation and control and source network charge storage is widely applied, a power grid terminal is explosive, terminal measurement data acquisition and interaction requirements are increasingly improved, an enterprise-level real-time measurement center is built, measurement data acquired by all link equipment of the novel power system are converged in real time, and the novel power system digital technology core foundation is realized. Therefore, the design of the synchronous calibration technology for the power grid data is urgently needed, the synchronous calibration is carried out on the data acquired from different equipment, the data are ensured to be on the same time section, the power grid data processing level is improved, and the normal operation of the power grid is ensured.
The patent of the invention 201110330295.5 proposes a section raw data alignment method based on double time scale delay estimation, which quantitatively estimates the delay time of each measuring point by adopting sampling sequences of two time scales according to the difference of data delay of different equipment transmission types, and predicts and corrects the current time discontinuous surface value of the measuring point by an autoregressive-sliding average model. The invention patent number 201810658052.6 proposes a method for quasi-synchronous starting and synchronous data acquisition of a harmonic on-line monitoring terminal under the condition of no time synchronization, which enables terminal equipment to actively send a starting signal to an adjacent terminal to realize synchronous starting of the adjacent terminal, and a measuring center realizes time section synchronization by comparing measured data with the starting signal. The traditional data synchronization calibration technology usually ignores the influence of a power grid structure on data transmission delay, and faces two problems when facing to the synchronous calibration of power grid data:
firstly, the patent with the invention number 201110330295.5 carries out time section synchronization by estimating the time delay of data transmission and forwarding of different devices, but the influence of fluctuation of link conditions on the estimated time delay is not considered, the estimation on the uncertainty of the time delay is lacked, the accuracy of data time synchronization is to be improved, and the time section deviation still exists. Therefore, how to design a more reasonable time slice correction step length determination method and optimize a data synchronization calibration method according to the arrival time of data uploaded by non-PMU equipment by using the network topology complexity, the measured data uploading information quantity, the historical signal to interference and noise ratio of a measured data transmission link, the link condition and other estimated time delay values and the uncertainty of the estimated time delay values is a urgent problem to be solved.
Secondly, the patent with the invention number 201810658052.6 realizes quasi-synchronous starting by sending a starting signal by adjacent equipment so as to perform time section synchronization, but the correlation of data uploaded by PMU equipment and non-PMU equipment is not considered, so that the data synchronization calibration result is inaccurate, and therefore, how to consider the correlation of active, reactive, voltage, current, switching value and other multi-domain data uploaded by the PMU equipment and the non-PMU equipment, calculate the time section deviation, optimize and correct the step length according to the uncertainty of the time delay estimated value and the time delay estimated value, and reduce the time section deviation is a urgent problem to be solved.
Aiming at the problems, the invention provides a multi-domain collaborative power grid data synchronous calibration technology and system based on network topology complexity, measurement data uploading information quantity, signal-to-interference-and-noise ratio of a measurement data transmission link and the like, a correction step length of a measured data section value is obtained based on a time delay estimated value and uncertainty of the time delay estimated value, synchronous calibration is carried out on data, and then time section deviation is obtained by solving correlation between the calibrated measured data and actual multi-domain data such as active, reactive, voltage, current, switching value and the like of PMU equipment, and correction step length optimization is carried out on acquired data without synchronous time marks.
Disclosure of Invention
The invention aims to solve the technical problems in the background technology and provides a multi-domain collaboration-based power grid data synchronous calibration method and system.
In order to solve the technical problems, the technical scheme of the invention is as follows:
a multi-domain collaboration-based grid data synchronization calibration method, the method comprising:
determining a time delay estimated value by using network topology complexity, measurement data uploading information quantity and historical signal-to-interference-and-noise ratio information of a measurement data transmission link, determining uncertainty of the transmission delay estimated value of the measurement data by using fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of the uploaded measurement data transmission link, and correcting the time delay estimated value to obtain a correction step length of a discontinuous face value when the measurement data is obtained;
according to the correction step length of the intermittent face value when measuring data, synchronously calibrating the data, and calculating the time section deviation between the measured data after synchronous calibration and the measured data actually transmitted by the PMU equipment;
and calculating an optimized correction step length of multi-domain data uploaded by non-PMU equipment based on the correction step length and the time section deviation, and updating a delay estimated value correction parameter and an uncertainty related parameter thereof based on the optimized correction step length.
Further, the determining the time delay estimation value specifically includes:
based on T moments, a set of moments is defined asIn an enterprise-level real-time measurement center, multiple sources of collected data exist at any one time, including collected data uploaded by PMU and non-PMU devices; consider that there are M non-PMU data collection devices in the grid that do not have synchronized time stamps, aggregated +.>Obtaining a time delay estimated value based on network topology complexity, measurement data uploading information quantity and historical signal to interference plus noise ratio information of a measurement data transmission link, and defining at time t, enterprise-level realityThe time measuring center receives the device r m After the measurement data is uploaded, the estimated time delay value is estimated to be tau m (t) expressed as:
wherein U is m (t) is the time t device r m Uploading measured data, wherein W (t) is the network topology complexity at the moment t, and the more the network topology is complex, the more the number of route hops is, the larger W (t) is, and the time delay estimated value tau is m The larger (t), B m (t) is the time t device r m Measuring data transmission bandwidth, SINR of data transmission link m (t) is the time t device r m The signal-to-interference-and-noise ratio of the data transmission link is measured.
Further, the determining the uncertainty of the estimated transmission delay value of the measurement data specifically includes:
based on r m Estimating r by uploading fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of transmission links of measurement data m Uncertainty delta of time delay estimated value of uploading measurement data m (t) expressed as:
wherein L is m (t) is time t, and r is obtained based on network topology complexity of time t-1, data transmission bandwidth of measurement data transmission link and signal-to-interference-and-noise ratio information of measurement data transmission link m Uploading link condition experience value of measurement data, the better the link condition, L m The larger (t), μ m (t) is t time r m Transmission link fluctuation mean weight, alpha, of uploading measured data m (t) is time t, r is obtained based on the empirical value of the link condition at the previous time t-1 m Uploading fluctuation condition of experience value of data transmission link condition, E m (t) is time t, based on the previous time t-1 r m Fluctuation obtained by uploading transmission link condition experience value of measurement dataMean value, beta m (t) is t time r m Transmission link condition fluctuation standard deviation weight, sigma, of uploading measured data m (t) is time t, based on the previous time t-1 r m Fluctuation standard deviation lambda obtained by uploading link condition experience value of measurement data m (t) is t time r m Maximum transmission link volatility weight for uploading measured data, wherein the formula is r m The larger the fluctuation mean value, variance and maximum value of the link condition experience value of the uploading measurement data, the uncertainty delta of the uploading time delay estimated value m The larger (t), L m (t)，α m (t)、E m (t) and sigma m (t) are respectively expressed as:
further, based on the uplink delay estimated value of the measurement data and the uncertainty of the transmission delay of the measurement data, the correction processing is performed on the delay estimated value to obtain a correction step length of the discontinuous face value when the measurement data is obtained, and the method specifically comprises the following steps:
correcting the time delay estimated value based on the obtained time delay estimated value and uncertainty of the time delay estimated value to obtain a correction step length of the measured data section value, synchronously calibrating the measured data section value through the correction step length, and defining a t-th moment measuring center pair r m The correction step length of the synchronous calibration of the uploaded measurement data is eta m (t) expressed as:
wherein omega u (t) is a delay estimation value correction parameter for adjusting the correction size of the delay estimation value;for event indicating variables, the link condition is better than the empirical condition at time t, i.e. the fluctuation of the empirical value of the link condition alpha m (t) > 0->On the contrary, let(s)>δ max The uncertainty threshold is used for measuring the uncertainty of the delay estimated value; the second term of the formula represents the deviation of the correction step from the delay estimate, when +. >When the link condition at time t is better than the experience condition, the correction step length should be smaller than the time delay estimated value, and the time delay estimated value tau m Subtracting the deviation on the basis of (t) to obtain a correction step length, otherwise, obtaining a time delay estimated value tau m Adding a bias to (t); at the same time, uncertainty delta of time delay estimated value m The larger the (t), the larger the deviation of the correction step from the original estimate; and based on the correction step length of the intermittent face value in the measurement data, carrying out synchronous calibration on the measurement data.
Further, calculating the time section deviation between the measurement data after synchronous calibration and the measurement data actually transmitted by the PMU device specifically includes:
for non-PMU devices r m The uploaded measurement data is estimated and calibrated, and a data set from multiple fields uploaded by the non-PMU equipment after calibration is defined as X m (t 0 )＝{x m,1 (t 0 ),...,x m,k (t 0 ),...,x m,K (t 0 ) K isThe number of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment, t 0 The non-PMU device measurement data obtained by calibration is uploaded at the moment x m,k (t 0 ) R represents m At t 0 K-th data uploaded at the moment; the presence of PMU devices around non-PMU devices defines that the PMU devices in the vicinity are at t 0 Time-of-day acquisition of uploaded data sets from multiple domains as X p (t 0 )＝{x p,1 (t 0 ),...,x p,k (t 0 ),...,x p,K (t 0 ) Definitions D m (t) is r after calibration m The calculation formula of the time section deviation between the uploaded measurement data and the actual measurement data uploaded by the PMU equipment is as follows:
Wherein MIC (X) m (t 0 -τ)；X p (t 0 ) X) represents X obtained based on maximum mutual information coefficient method m (t 0 - τ) and X p (t 0 ) Correlation between I (X) m (t 0 -τ)；X p (t 0 ) For X) m (t 0 - τ) and X p (t 0 ) The maximum mutual information value obtained after the meshing of the row a and the column B is carried out, B is a grid scale parameter, and the data quantity is taken to the power of 0.6; the formula means X m (t 0 - τ) and X p (t 0 ) The greater the correlation between them, the more likely the non-PMU device is at t 0 - τ time and PMU device at t 0 The greater the correlation of the time-point uploading data, the time section deviation is taken as tau at the moment.
Further, calculating an optimized correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment, wherein the optimized correction step length specifically comprises the following steps:
combined with time section deviation D m (t) time delay estimationValue τ m (t) uncertainty of the time-lapse estimation value delta m (t) performing correction step optimization on the fault values of the active, reactive, voltage, current and switching value multi-domain data uploaded by the non-PMU equipment, and then performing correction step eta 'on the fault values of the multi-domain measurement data uploaded by the non-PMU equipment' m The calculation formula of (t) is as follows:
this formula shows that the time section deviation D calculated by formula (5) m (t) is the difference between the time section between the non-PMU device uploading measurement data after calibration and the actual PMU device uploading measurement data, and the time section of the PMU device is known to be accurate, so that the t-th instant enterprise level real-time measurement center pair r obtained by the formula (4) m Correction step length eta for synchronous calibration of uploaded measurement data m (t) subtracting the time-section deviation D obtained by the formula (5) m (t) is the optimized correction step length eta 'of the active, reactive, voltage, current and switching value multi-domain data uploaded by the non-PMU equipment' m (t)。
Further, the updating the delay estimation value correction parameter and the uncertainty related parameter thereof includes: updating the time delay estimated value correction parameter, the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight parameter.
Further, the updating the delay estimation value correction parameter and the uncertainty related parameter thereof further includes:
based on correction step optimization of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment, updating a formula for obtaining uncertainty of a time delay estimated value, if the correction step is reduced after optimization, the uncertainty of the time delay estimated value should be correspondingly increased, if the correction step is increased after optimization, the uncertainty of the time delay estimated value should be correspondingly reduced, so that the correction step for synchronously calibrating multi-domain measurement data at the next moment is closer to an accurate value, and the fluctuation average value weight of a transmission link is updated Heavy mu m (t) transmission link volatility standard deviation weight beta m (t), maximum value weight lambda of transmission link fluctuation m (t) updating the calculation of uncertainty of the delay estimation value, then transmitting the link fluctuation mean weight mu m (t) transmission link volatility standard deviation weight beta m (t), maximum value weight lambda of transmission link fluctuation m (t) updating the formula:
wherein a is 1 ，b 1 Updating coefficients, a, for transmission link fluctuation mean weights 2 ，b 2 Updating coefficients, a, for transmission link volatility standard deviation weights 3 Updating coefficients for transmission link volatility maximum value weights; the formula is expressed as the ratio of the update of the transmission link fluctuation mean weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight to the optimized correction step length and the correction step length before optimizationRelated to; the smaller the optimized correction step length is, the larger the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight are, the larger the uncertainty of the time delay estimated value is, and the correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment at the time t+1 is reduced to be closer to the actual uploading time of the PMU equipment measured data;
Based on active, reactive, voltage, current, and switching value multi-domain data uploaded to non-PMU devicesCorrection step optimization is carried out, and the parameter omega is corrected for 1 time delay estimated value u (t) updating, if the optimized correction step length is reduced, the time delay estimation value correction parameter should be correspondingly increased, and if the optimized correction step length is increased, the time delay estimation value correction parameter should be correspondingly reduced, so that the correction step length for synchronously calibrating the multi-domain measurement data at the next moment is closer to an accurate value, and the time delay estimation value correction parameter omega u (t) updating the formula:
wherein a is 4 ，b 3 Correcting parameter updating coefficients for the time delay estimated value; the formula shows that the ratio of the update of the correction parameters of the time delay estimated value to the correction step length before optimization and the correction step length after optimizationIn the related process, the smaller the optimized correction step length is, the larger the correction parameter of the delay estimated value is, and the correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment at the time t+1 is reduced to be closer to the actual uploading time of the PMU equipment measurement data.
A multi-domain collaboration-based grid data synchronization calibration system, the system being applied to the method of any one of the above, the system comprising: the device comprises a communication module, a data storage module, a time delay estimation module, a time delay uncertainty estimation module, a correction step length calculation module, a time section deviation calculation module, a correction step length optimization module and a parameter updating module;
The communication module: the device is used for communicating the PMU and the non-PMU equipment to the real-time measurement center, and uploading measurement data acquired by the PMU and the non-PMU equipment to the real-time measurement center;
the data storage module: the PMU and non-PMU equipment are used for storing measurement data uploaded to the real-time measurement center by the PMU and non-PMU equipment;
the delay estimation module: estimating a transmission delay estimated value of the measurement data based on fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of the experience condition of the transmission link of the uploaded measurement data;
the delay uncertainty estimation module: estimating uncertainty of transmission delay of the measurement data based on fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of the transmission link;
the correction step length calculation module: correcting the delay estimated value based on the estimated value of the transmission delay and the uncertainty of the delay estimated value to obtain a correction step length of the discontinuous face value when measuring data;
the data synchronization calibration module: the method comprises the steps of performing synchronous calibration on data according to the correction step length of intermittent face values during data measurement;
the time section deviation calculation module is used for: the PMU device is used for calculating the time section deviation between the measurement data after synchronous calibration and the measurement data actually transmitted by the PMU device;
The correction step optimization module: based on the correction step length and the time section deviation, calculating an optimized correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment;
the parameter updating module: and updating the time delay estimated value correction parameter, the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight parameter based on the optimized correction step length.
A computer readable storage medium having stored therein computer executable instructions which when executed by a processor are adapted to carry out the method of any one of the preceding claims.
Compared with the prior art, the invention has the advantages that:
(1) The invention provides a correction step length determining method based on delay uncertainty estimation, which estimates an uplink delay estimated value of measured data based on network topology complexity, measured data uplink information quantity, signal to interference and noise ratio of a measured data transmission link and other information. In consideration of uncertainty of the delay estimated value caused by fluctuation of the link condition, the estimated value cannot be directly used for correcting the measured data section. Based on the information such as fluctuation mean value, fluctuation maximum value, fluctuation standard deviation and the like of the experience condition of the transmission link of the measurement data, the uncertainty of the delay estimated value is estimated, the delay estimated value is corrected, the correction step length of the section value of the measurement data is obtained, and the synchronous calibration of the section value of the measurement data is realized.
(2) The invention provides a correction step optimization method based on multi-domain synergy, which is characterized in that the correction step optimization is carried out on the section values of active, reactive, voltage, current, switching value and other multi-domain data uploaded by non-PMU equipment by comparing the calibrated measurement data with actual measurement data uploaded by PMU equipment to obtain time section deviation and combining network topology complexity, time delay estimated value and uncertainty of the time delay estimated value, and the correction step optimization is carried out on the section values of the time delay estimated value based on the correction step optimization result, and the correction parameters of the time delay estimated value are updated to realize the section data alignment of a large-scale power grid.
Drawings
FIG. 1, a multi-domain collaboration-based grid data synchronization calibration system;
FIG. 2 is a flow chart of synchronous calibration of power grid data based on multi-domain cooperation.
Detailed Description
The following describes specific embodiments of the present invention with reference to examples:
it should be noted that the structures, proportions, sizes and the like illustrated in the present specification are used for being understood and read by those skilled in the art in combination with the disclosure of the present invention, and are not intended to limit the applicable limitations of the present invention, and any structural modifications, proportional changes or size adjustments should still fall within the scope of the disclosure of the present invention without affecting the efficacy and achievement of the present invention.
Also, the terms such as "upper," "lower," "left," "right," "middle," and "a" and the like recited in the present specification are merely for descriptive purposes and are not intended to limit the scope of the invention, but are intended to provide relative positional changes or modifications without materially altering the technical context in which the invention may be practiced.
Example 1:
as shown in fig. 1, the embodiment provides a multi-domain collaboration-based power grid data synchronization calibration system, which comprises a power supply module, a communication module, a data storage module, a time delay estimation module, a time delay uncertainty estimation module, a correction step calculation module, a time section deviation calculation module, a correction step optimization module and a parameter update module, wherein the modules are as follows:
and a power supply module: the power supply unit supplies power to all modules in the power grid data synchronous calibration system.
And a communication module: and the communication from the PMU and the non-PMU equipment to the enterprise-level real-time measurement center is realized, and the measurement data collected by the PMU and the non-PMU equipment are uploaded to the enterprise-level real-time measurement center.
And a data storage module: and the device is responsible for storing measurement data uploaded to the enterprise-level real-time measurement center by the PMU and the non-PMU equipment.
And a time delay estimation module: based on the information such as the fluctuation mean value, the fluctuation maximum value, the fluctuation standard deviation and the like of the experience condition of the transmission link of the uploaded measurement data, the transmission delay of the measurement data is estimated.
A delay uncertainty estimation module: based on the information such as the fluctuation mean value, the fluctuation maximum value, the fluctuation standard deviation and the like of the experience condition of the transmission link, the uncertainty of the transmission delay of the measurement data is estimated.
And a correction step length calculation module: and correcting the delay estimated value based on the uncertainty of the delay estimated value and the delay estimated value to obtain a correction step length of the measured data fracture surface value.
And the data synchronization calibration module: and according to the correction step length of the intermittent face value when measuring the data, synchronously calibrating the data.
The time section deviation calculating module is used for: and calculating the time section deviation between the measured data after calibration and the measured data actually transmitted by the PMU equipment.
Correction step optimization module: and calculating the optimal correction step length of the multi-domain data such as active power, reactive power, voltage, current, switching value and the like uploaded by the non-PMU equipment based on the correction step length and the time section deviation.
Parameter updating module: based on the correction step length before and after optimization, updating the parameters such as the correction parameter of the time delay estimated value, the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight, the transmission link fluctuation maximum value weight and the like.
Example 2:
the embodiment provides a multi-domain collaboration-based power grid data synchronous calibration method, which is applied to the system of the embodiment 1, and as shown in fig. 2, by estimating the uncertainty of a time delay estimated value and a time delay estimated value, a correction step length of a measured data section value is obtained, and the correction step length is optimized based on a time section deviation, so that data synchronous calibration is realized, and the data is ensured to be on the same time section. The method provided by the embodiment mainly comprises a correction step length determining method based on delay uncertainty estimation and a correction step length optimizing method based on multi-domain cooperation, and specifically comprises the following steps:
1. Correction step length determining method based on time delay uncertainty estimation
1.1 estimating the measurement data upload delay estimate
The present embodiment considers T moments, the set of moments is defined asIn an enterprise-level real-time measurement center, there are multiple sources of collected data at any one time, including collected data uploaded by PMU, non-PMU devices. Consider that there are M non-PMU data collection devices in the grid that do not have synchronized time stamps, aggregated +.>Because the data uploaded by the non-PMU device does not have a synchronous time scale, the enterprise-level real-time measurement center cannot obtain the uploading time of the received measurement data of the non-PMU device, so that the time sections of the measurement data are difficult to align. The time section is the running state of the power system at a certain moment and comprises active, reactive, voltage, current, switching value and other multi-domain data of the running working conditions in the power system at the moment. Thus, enterprise level realityThe time measurement center performs synchronous calibration on the received measurement data by estimating the time delay of uploading the measurement data. In this embodiment, a time delay estimated value is obtained based on information such as network topology complexity, measurement data uploading information quantity, and historical signal-to-interference-and-noise ratio of a measurement data transmission link, and is defined at time t, and an enterprise-level real-time measurement center receives a device r m After the measurement data is uploaded, the estimated time delay value is estimated to be tau m (t) expressed as
Wherein U is m (t) is the time t device r m Uploading measured data, wherein W (t) is the network topology complexity at the moment t, and the more the network topology is complex, the more the number of route hops is, the larger W (t) is, and the time delay estimated value tau is m The larger (t), B m (t) is the time t device r m Measuring data transmission bandwidth, SINR of data transmission link m (t) is the time t device r m The signal-to-interference-and-noise ratio of the data transmission link is measured.
1.2 estimation of uncertainty of the upload delay estimation value
Based on the uncertainty of the delay estimated in step 1.1 due to fluctuation of the link condition, the uncertainty represents the uncertainty probability of the estimated delay value, the worse the link condition is, the more inaccurate the estimated time delay is, the larger the uncertainty of the delay estimated value is, and the time delay of the acquired data uploading is difficult to be estimated accurately, and the embodiment is based on r m Estimating r by uploading information such as fluctuation mean value, fluctuation maximum value, fluctuation standard deviation and the like of experience conditions of transmission links of measurement data m Uncertainty delta of time delay estimated value of uploading measurement data m (t) expressed as
Wherein L is m (t) is time t, and based on network topology complexity at time t-1, measuring data transmission bandwidth of data transmission link R obtained by measuring signal-interference-noise ratio and other information of data transmission link m Uploading link condition experience value of measurement data, the better the link condition, L m The larger (t). Mu (mu) m (t) is t time r m Transmission link fluctuation mean weight, alpha, of uploading measured data m (t) is time t, r is obtained based on the empirical value of the link condition at the previous time t-1 m Uploading fluctuation condition of experience value of data transmission link condition, E m (t) is time t, based on the previous time t-1 r m Fluctuation mean value beta obtained by uploading transmission link condition experience value of measurement data m (t) is t time r m Transmission link condition fluctuation standard deviation weight, sigma, of uploading measured data m (t) is time t, based on the previous time t-1 r m Fluctuation standard deviation lambda obtained by uploading link condition experience value of measurement data m (t) is t time r m Maximum transmission link volatility weight for uploading measured data, wherein the formula is r m The larger the fluctuation mean value, variance and maximum value of the link condition experience value of the uploading measurement data, the uncertainty delta of the uploading time delay estimated value m The larger (t), L m (t)，α m (t)、E m (t) and sigma m (t) are respectively expressed as
1.3 calculating correction step size of measured data section value
ConsiderSince the delay estimated value has uncertainty caused by fluctuation of the link condition and cannot be directly used for correcting the measured data fracture value, the delay estimated value is corrected based on the uncertainty of the delay estimated value obtained in the step 1.1 and the delay estimated value obtained in the step 1.2, the correction step length of the measured data fracture value is obtained, the measured data fracture value is synchronously calibrated through the correction step length, and the t moment measuring center pair r is defined m The correction step length of the synchronous calibration of the uploaded measurement data is eta m (t) expressed as
Wherein omega u And (t) is a delay estimated value correction parameter used for adjusting the correction size of the delay estimated value.For event indicating variables, the link condition is better than the empirical condition at time t, i.e. the fluctuation of the empirical value of the link condition alpha m (t) > 0->On the contrary, let(s)>δ max Is an uncertainty threshold and is used for measuring the uncertainty of the delay estimated value. The second term of the formula represents the deviation of the correction step from the delay estimate, when +.>When the link condition at time t is better than the experience condition, the correction step length should be smaller than the time delay estimated value, and the time delay estimated value tau m Subtracting the deviation on the basis of (t) to obtain a correction step length, otherwise, obtaining a time delay estimated value tau m Adding a bias to (t); at the same time, uncertainty delta of time delay estimated value m The larger (t), the larger the deviation of the correction step from the original estimate. Based on the amountAnd correcting step length of the intermittent face value during data measurement, and synchronously calibrating the measured data.
2. Correction step length optimization method based on multi-domain cooperation
In order to further calibrate the data synchronously, a correction step optimization method based on multi-domain cooperation is provided, the corrected data obtained in the step 2.1.3 is compared with the PMU uploading data to obtain the time section deviation, and the correction step and the updating parameters are adjusted based on the time section deviation.
2.1 obtaining the time section deviation
The measurement center calculates the non-PMU device r according to step 2.1 m The uploaded measurement data is estimated and calibrated, and a data set from multiple fields uploaded by the non-PMU equipment after calibration is defined as X m (t 0 )＝{x m,1 (t 0 ),...,x m,k (t 0 ),...,x m,K (t 0 ) K is the number of active, reactive, voltage, current, switching value and other multi-domain data uploaded by non-PMU equipment, t 0 The non-PMU device measurement data obtained by calibration is uploaded at the moment x m,k (t 0 ) R represents m At t 0 And uploading kth data at the moment. The presence of PMU devices around non-PMU devices defines that the PMU devices in the vicinity are at t 0 Time-of-day acquisition of uploaded data sets from multiple domains as X p (t 0 )＝{x p,1 (t 0 ),...,x p,k (t 0 ),...,x p,K (t 0 ) Definitions D m (t) is r after calibration m The time section deviation between the uploaded measurement data and the actual measurement data of the PMU equipment is calculated by the following formula
Wherein MIC (X) m (t 0 -τ)；X p (t 0 ) Representation based on maximumX obtained by mutual information coefficient method m (t 0 - τ) and X p (t 0 ) Correlation between I (X) m (t 0 -τ)；X p (t 0 ) For X) m (t 0 - τ) and X p (t 0 ) The maximum mutual information value obtained after the meshing of the a rows and the B columns is obtained, B is a grid scale parameter, and the data quantity is generally taken to the power of 0.6. The formula means X m (t 0 - τ) and X p (t 0 ) The greater the correlation between them, the more likely the non-PMU device is at t 0 - τ time and PMU device at t 0 The greater the correlation of the time-point uploading data, the time section deviation is taken as tau at the moment.
2.2 correction step optimization
Combining the time section deviation D obtained in the step 2.1 m (t) time delay estimated value tau obtained in step 1.1 m (t), and uncertainty delta of the delay estimation value obtained in the step 1.2 m (t) performing correction step optimization on the fault values of the multi-domain data such as active, reactive, voltage, current, switching value and the like uploaded by the non-PMU equipment, and then performing correction step eta 'on the fault values of the multi-domain measured data uploaded by the non-PMU equipment' m (t) the calculation formula is
This formula shows that the time section deviation D calculated by formula (5) m (t) is the difference between the time section between the non-PMU device uploading measurement data after calibration and the actual PMU device uploading measurement data, and the time section of the PMU device is known to be accurate, so that the t-th instant enterprise level real-time measurement center pair r obtained by the formula (4) m Correction step length eta for synchronous calibration of uploaded measurement data m (t) subtracting the time-section deviation D obtained by the formula (5) m (t) is the optimized correction step length eta 'of the multi-domain data such as active power, reactive power, voltage, current, switching value and the like uploaded by the non-PMU equipment' m (t)。
2.3 time delay estimation value and uncertainty related parameter update
Based on the correction step optimization of the multi-domain data such as active, reactive, voltage, current, switching value and the like uploaded by non-PMU equipment in 2.2, updating the formula for obtaining the uncertainty of the time delay estimated value in 1.2, if the correction step is reduced after the optimization, the uncertainty of the time delay estimated value should be correspondingly increased, if the correction step is increased after the optimization, the uncertainty of the time delay estimated value should be correspondingly reduced, so that the correction step for synchronously calibrating the multi-domain measurement data at the next moment is closer to the accurate value, and the fluctuation mean weight mu of the transmission link can be updated m (t) transmission link volatility standard deviation weight beta m (t), maximum value weight lambda of transmission link fluctuation m (t) updating the calculation of uncertainty of the delay estimation value, then transmitting the link fluctuation mean weight mu m (t) transmission link volatility standard deviation weight beta m (t), maximum value weight lambda of transmission link fluctuation m (t) updating the formula to
Wherein a is 1 ，b 1 Updating coefficients, a, for transmission link fluctuation mean weights 2 ，b 2 Updating coefficients, a, for transmission link volatility standard deviation weights 3 The coefficients are updated for the transmission link volatility maximum value weight. The formula shows that the ratio of the average fluctuation value weight of the transmission link, the standard deviation fluctuation weight of the transmission link and the maximum fluctuation value weight of the transmission link to the optimized correction step length and the correction step length before optimizationRelated to the following. The smaller the optimized correction step length is, the larger the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight are, the larger the uncertainty of the time delay estimated value is, and the correction step length of the multi-domain data such as active, reactive, voltage, current, switching value and the like uploaded by the non-PMU equipment at the time t+1 is reduced to be closer to the actual uploading time of the PMU equipment measured data, and vice versa.
Correction step optimization of the time delay estimated value correction parameter omega in 1.3 is based on correction step optimization of 2.2 on active, reactive, voltage, current, switching value and other multi-domain data uploaded by non-PMU equipment u (t) updating, if the optimized correction step length is reduced, the time delay estimation value correction parameter should be correspondingly increased, and if the optimized correction step length is increased, the time delay estimation value correction parameter should be correspondingly reduced, so that the correction step length for synchronously calibrating the multi-domain measurement data at the next moment is closer to an accurate value, and the time delay estimation value correction parameter omega u (t) updating the formula to
Wherein a is 4 ，b 3 And correcting parameter updating coefficients for the time delay estimated value. The formula shows that the ratio of the update of the correction parameters of the time delay estimated value to the correction step length before optimization and the correction step length after optimizationIn this regard, the smaller the optimized correction step size is, the larger the correction parameter of the delay estimation value is, and the correction step size of the multi-domain data such as active, reactive, voltage, current, switching value and the like uploaded by the non-PMU device at the time t+1 is reduced to be closer to the actual uploading time of the measured data of the PMU device, and vice versa.
Example 3:
those of ordinary skill in the art will appreciate that all or a portion of the steps of the various methods of the above embodiments may be performed by instructions, or by instructions controlling associated hardware, which may be stored in a computer-readable storage medium and loaded and executed by a processor.
To this end, an embodiment of the present invention provides a storage medium having stored therein a plurality of instructions that can be loaded by a processor to perform the steps of any of the multi-domain collaboration-based grid data synchronization calibration methods provided by the embodiments of the present invention.
For example, the instructions may perform the steps of:
a power grid data synchronous calibration method based on multi-domain cooperation comprises the following steps:
determining a time delay estimated value by using network topology complexity, measurement data uploading information quantity and historical signal-to-interference-and-noise ratio information of a measurement data transmission link, determining uncertainty of the transmission delay estimated value of the measurement data by using fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of the uploaded measurement data transmission link, and correcting the time delay estimated value to obtain a correction step length of a discontinuous face value when the measurement data is obtained;
according to the correction step length of the intermittent face value when measuring data, synchronously calibrating the data, and calculating the time section deviation between the measured data after synchronous calibration and the measured data actually transmitted by the PMU equipment;
and calculating an optimized correction step length of multi-domain data uploaded by non-PMU equipment based on the correction step length and the time section deviation, and updating a delay estimated value correction parameter and an uncertainty related parameter thereof based on the optimized correction step length.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While the preferred embodiments of the present invention have been described in detail, the present invention is not limited to the above embodiments, and various changes may be made without departing from the spirit of the present invention within the knowledge of those skilled in the art.
Many other changes and modifications may be made without departing from the spirit and scope of the invention. It is to be understood that the invention is not to be limited to the specific embodiments, but only by the scope of the appended claims.
Claims (10)
1. A multi-domain collaboration-based power grid data synchronous calibration method, the method comprising:
determining a time delay estimated value by using network topology complexity, measurement data uploading information quantity and historical signal-to-interference-and-noise ratio information of a measurement data transmission link, determining uncertainty of the transmission delay estimated value of the measurement data by using fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of the uploaded measurement data transmission link, and correcting the time delay estimated value to obtain a correction step length of a discontinuous face value when the measurement data is obtained;
According to the correction step length of the intermittent face value when measuring data, synchronously calibrating the data, and calculating the time section deviation between the measured data after synchronous calibration and the measured data actually transmitted by the PMU equipment;
and calculating an optimized correction step length of multi-domain data uploaded by non-PMU equipment based on the correction step length and the time section deviation, and updating a delay estimated value correction parameter and an uncertainty related parameter thereof based on the optimized correction step length.
2. The method for synchronously calibrating the power grid data based on multi-domain cooperation according to claim 1, wherein the determining the time delay estimated value specifically comprises:
based on T moments, a set of moments is defined asIn an enterprise-level real-time measurement center, multiple sources of collected data exist at any one time, including collected data uploaded by PMU and non-PMU devices; consider that there are M non-PMU data collection devices in the grid that do not have synchronized time stamps, aggregated +.>Obtaining a time delay estimated value based on network topology complexity, measured data uploading information quantity and historical signal-to-interference-and-noise ratio information of a measured data transmission link, and defining that at time t, an enterprise-level real-time measuring center receives equipment r m After the measurement data is uploaded, the estimated time delay value is estimated to be tau m (t) expressed as:
wherein U is m (t) is the time t device r m Uploading measured data, wherein W (t) is the network topology complexity at the moment t, and the more the network topology is complex, the more the number of route hops is, the larger W (t) is, and the time delay estimated value tau is m The larger (t), B m (t) is the time t device r m Measuring data transmission bandwidth, SINR of data transmission link m (t) is the time t device r m The signal-to-interference-and-noise ratio of the data transmission link is measured.
3. The method for synchronously calibrating power grid data based on multi-domain cooperation according to claim 1, wherein the determining the uncertainty of the estimated transmission delay value of the measurement data specifically comprises:
based on r m Estimating r by uploading fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of transmission links of measurement data m Uncertainty delta of time delay estimated value of uploading measurement data m (t) expressed as:
wherein L is m (t) is time t, and r is obtained based on network topology complexity of time t-1, data transmission bandwidth of measurement data transmission link and signal-to-interference-and-noise ratio information of measurement data transmission link m Uploading link condition experience value of measurement data, the better the link condition, L m The larger (t), μ m (t) is t time r m Transmission link fluctuation mean weight, alpha, of uploading measured data m (t) is time t, r is obtained based on the empirical value of the link condition at the previous time t-1 m Uploading fluctuation condition of experience value of data transmission link condition, E m (t) is time t, based on the previous time t-1 r m Fluctuation mean value beta obtained by uploading transmission link condition experience value of measurement data m (t) is t time r m Transmission link condition fluctuation standard deviation weight, sigma, of uploading measured data m (t) is time t, based on the previous time t-1 r m Fluctuation standard deviation lambda obtained by uploading link condition experience value of measurement data m (t) is t time r m Maximum transmission link volatility weight for uploading measured data, wherein the formula is r m Upload measurementThe larger the fluctuation mean value, variance and maximum value of the link condition experience value of the data, the uncertainty delta of the uploading delay estimated value m The larger (t), L m (t)，α m (t)、E m (t) and sigma m (t) are respectively expressed as:
4. the method for synchronously calibrating the power grid data based on multi-domain cooperation according to claim 1, wherein the method is characterized by carrying out correction processing on the time delay estimated value based on the uploading time delay estimated value of the measured data and the uncertainty of the transmission time delay of the measured data to obtain the correction step length of the intermittent face value when the measured data is obtained, and specifically comprises the following steps:
correcting the time delay estimated value based on the obtained time delay estimated value and uncertainty of the time delay estimated value to obtain a correction step length of the measured data section value, synchronously calibrating the measured data section value through the correction step length, and defining a t-th moment measuring center pair r m The correction step length of the synchronous calibration of the uploaded measurement data is eta m (t) expressed as:
wherein omega u (t) is a delay estimation value correction parameter for adjusting the correction size of the delay estimation value;for event indicating variables, the link condition is better than the empirical condition at time t, i.e. the fluctuation of the empirical value of the link condition alpha m When (t) > 0 is used,on the contrary, let(s)>δ max The uncertainty threshold is used for measuring the uncertainty of the delay estimated value; the second term of the formula represents the deviation of the correction step from the delay estimate, when +.>When the link condition at time t is better than the experience condition, the correction step length should be smaller than the time delay estimated value, and the time delay estimated value tau m Subtracting the deviation on the basis of (t) to obtain a correction step length, otherwise, obtaining a time delay estimated value tau m Adding a bias to (t); at the same time, uncertainty delta of time delay estimated value m The larger the (t), the larger the deviation of the correction step from the original estimate; and based on the correction step length of the intermittent face value in the measurement data, carrying out synchronous calibration on the measurement data.
5. The method for synchronously calibrating power grid data based on multi-domain cooperation according to claim 1, wherein calculating the time section deviation between the measured data after synchronous calibration and the measured data actually transmitted by the PMU device comprises:
For non-PMU devices r m The uploaded measurement data is estimated and calibrated, and a data set from multiple fields uploaded by the non-PMU equipment after calibration is defined as X m (t 0 )＝{x m,1 (t 0 ),...,x m,k (t 0 ),...,x m,K (t 0 ) K is the number of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment, t 0 The non-PMU device measurement data obtained by calibration is uploaded at the moment x m,k (t 0 ) R represents m At t 0 K-th data uploaded at the moment; the presence of PMU devices around non-PMU devices defines that the PMU devices in the vicinity are at t 0 Time-of-day acquisition of uploaded data sets from multiple domains as X p (t 0 )＝{x p,1 (t 0 ),...,x p,k (t 0 ),...,x p,K (t 0 ) Definitions D m (t) is r after calibration m The calculation formula of the time section deviation between the uploaded measurement data and the actual measurement data uploaded by the PMU equipment is as follows:
wherein MIC (X) m (t 0 -τ)；X p (t 0 ) X) represents X obtained based on maximum mutual information coefficient method m (t 0 - τ) and X p (t 0 ) Correlation between I (X) m (t 0 -τ)；X p (t 0 ) For X) m (t 0 - τ) and X p (t 0 ) The maximum mutual information value obtained after the meshing of the row a and the column B is carried out, B is a grid scale parameter, and the data quantity is taken to the power of 0.6; the formula means X m (t 0 - τ) and X p (t 0 ) The greater the correlation between them, the more likely the non-PMU device is at t 0 - τ time and PMU device at t 0 The greater the correlation of the time-point uploading data, the time section deviation is taken as tau at the moment.
6. The method for synchronously calibrating power grid data based on multi-domain cooperation according to claim 1, wherein the step of optimizing and correcting the active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment is calculated, and specifically comprises the following steps:
Combined with time section deviation D m (t) time delay estimation value tau m (t) uncertainty of the time-lapse estimation value delta m (t) performing correction step optimization on the fault values of the active, reactive, voltage, current and switching value multi-domain data uploaded by the non-PMU equipment, and then performing correction step eta 'on the fault values of the multi-domain measurement data uploaded by the non-PMU equipment' m The calculation formula of (t) is as follows:
this formula shows that the time section deviation D calculated by formula (5) m (t) is the difference between the time section between the non-PMU device uploading measurement data after calibration and the actual PMU device uploading measurement data, and the time section of the PMU device is known to be accurate, so that the t-th instant enterprise level real-time measurement center pair r obtained by the formula (4) m Correction step length eta for synchronous calibration of uploaded measurement data m (t) subtracting the time-section deviation D obtained by the formula (5) m (t) is the optimized correction step length eta 'of the active, reactive, voltage, current and switching value multi-domain data uploaded by the non-PMU equipment' m (t)。
7. The method for synchronously calibrating power grid data based on multi-domain cooperation according to claim 1, wherein the updating the delay estimation value correction parameter and the uncertainty related parameter thereof comprises the following steps: updating the time delay estimated value correction parameter, the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight parameter.
8. The method for synchronously calibrating power grid data based on multi-domain cooperation according to claim 7, wherein the updating the delay estimation value correction parameter and the uncertainty related parameter thereof further comprises:
based on correction step optimization of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment, updating a formula for obtaining uncertainty of a time delay estimated value, if the correction step is reduced after optimization, the uncertainty of the time delay estimated value should be correspondingly increased, if the correction step is increased after optimization, the uncertainty of the time delay estimated value should be correspondingly reduced, so that the correction step for synchronously calibrating multi-domain measurement data at the next moment is closer to an accurate value, and the fluctuation mean weight mu of a transmission link is updated m (t) transmission link volatility standard deviation weight beta m (t), maximum value weight lambda of transmission link fluctuation m (t) updating the calculation of uncertainty of the delay estimation value, then transmitting the link fluctuation mean weight mu m (t) transmission link volatility standard deviation weight beta m (t), maximum value weight lambda of transmission link fluctuation m (t) updating the formula:
wherein a is 1 ，b 1 Updating coefficients, a, for transmission link fluctuation mean weights 2 ，b 2 Updating coefficients, a, for transmission link volatility standard deviation weights 3 Updating coefficients for transmission link volatility maximum value weights; the formula is expressed as the ratio of the update of the transmission link fluctuation mean weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight to the optimized correction step length and the correction step length before optimizationRelated to; the smaller the optimized correction step length is, the larger the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight are, the larger the uncertainty of the time delay estimated value is, and the correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment at the time t+1 is reduced to be closer to the actual uploading time of the PMU equipment measured data;
correction step optimization of 1-time delay estimated value correction parameter omega based on active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment u (t) updating, if the optimized correction step length is reduced, the time delay estimation value correction parameter should be correspondingly increased, and if the optimized correction step length is increased, the time delay estimation value correction parameter should be correspondingly reduced, so that the correction step length for synchronously calibrating the multi-domain measurement data at the next moment is closer to an accurate value, and the time delay estimation value correction parameter omega u (t) updating the formula:
wherein a is 4 ，b 3 Correcting parameter updating coefficients for the time delay estimated value; the formula shows that the ratio of the update of the correction parameters of the time delay estimated value to the correction step length before optimization and the correction step length after optimizationIn the related process, the smaller the optimized correction step length is, the larger the correction parameter of the delay estimated value is, and the correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment at the time t+1 is reduced to be closer to the actual uploading time of the PMU equipment measurement data.
9. A multi-domain collaboration based grid data synchronization calibration system, characterized in that the system is applied to the method of any one of claims 1-8, the system comprising: the device comprises a communication module, a data storage module, a time delay estimation module, a time delay uncertainty estimation module, a correction step length calculation module, a time section deviation calculation module, a correction step length optimization module and a parameter updating module;
the communication module: the device is used for communicating the PMU and the non-PMU equipment to the real-time measurement center, and uploading measurement data acquired by the PMU and the non-PMU equipment to the real-time measurement center;
the data storage module: the PMU and non-PMU equipment are used for storing measurement data uploaded to the real-time measurement center by the PMU and non-PMU equipment;
The delay estimation module: estimating a transmission delay estimated value of the measurement data based on fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of the experience condition of the transmission link of the uploaded measurement data;
the delay uncertainty estimation module: estimating uncertainty of transmission delay of the measurement data based on fluctuation mean value, fluctuation maximum value and fluctuation standard deviation information of experience conditions of the transmission link;
the correction step length calculation module: correcting the delay estimated value based on the estimated value of the transmission delay and the uncertainty of the delay estimated value to obtain a correction step length of the discontinuous face value when measuring data;
the data synchronization calibration module: the method comprises the steps of performing synchronous calibration on data according to the correction step length of intermittent face values during data measurement;
the time section deviation calculation module is used for: the PMU device is used for calculating the time section deviation between the measurement data after synchronous calibration and the measurement data actually transmitted by the PMU device;
the correction step optimization module: based on the correction step length and the time section deviation, calculating an optimized correction step length of active, reactive, voltage, current and switching value multi-domain data uploaded by non-PMU equipment;
the parameter updating module: and updating the time delay estimated value correction parameter, the transmission link fluctuation mean value weight, the transmission link fluctuation standard deviation weight and the transmission link fluctuation maximum value weight parameter based on the optimized correction step length.
10. A computer readable storage medium having stored therein computer executable instructions which when executed by a processor are adapted to carry out the method of any one of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311101913.8A CN117060592B (en) | 2023-08-30 | 2023-08-30 | Multi-domain cooperation-based power grid data synchronous calibration method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311101913.8A CN117060592B (en) | 2023-08-30 | 2023-08-30 | Multi-domain cooperation-based power grid data synchronous calibration method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117060592A true CN117060592A (en) | 2023-11-14 |
| CN117060592B CN117060592B (en) | 2024-02-23 |
Family
ID=88662460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311101913.8A Active CN117060592B (en) | 2023-08-30 | 2023-08-30 | Multi-domain cooperation-based power grid data synchronous calibration method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117060592B (en) |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20080035116A (en) * | 2006-10-18 | 2008-04-23 | 주식회사 세니온 | Automation system and method for automatically detecting phase of electric power distribution line |
| CN102195354A (en) * | 2011-05-16 | 2011-09-21 | 中国电力科学研究院 | Large-grid wide-area asynchronous coordination control method |
| CN102361353A (en) * | 2011-10-26 | 2012-02-22 | 东北电网有限公司 | Method for aligning section raw data on basis of double time scale delay evaluation |
| CN103323742A (en) * | 2013-05-31 | 2013-09-25 | 华北电力大学 | System and method for grid fault identification based on random point distribution PMU |
| CN105914798A (en) * | 2016-04-08 | 2016-08-31 | 李仕乙 | Wind power field reactive power optimization operation control method based on whole network same time section |
| CN106383287A (en) * | 2016-12-05 | 2017-02-08 | 国网辽宁省电力有限公司丹东供电公司 | Method for improving power grid stable section measurement precision through improving data synchronism |
| CN107144760A (en) * | 2016-12-31 | 2017-09-08 | 国电南瑞科技股份有限公司 | It is a kind of to support electric network fault to perceive panorama Wave record method and plateform system with analysis |
| GB201810680D0 (en) * | 2018-06-28 | 2018-08-15 | Synaptec Ltd | Methods and apparatus for making a time-synchronised phasor measurement |
| CN109100614A (en) * | 2018-09-28 | 2018-12-28 | 东北大学 | A kind of transmission open acess system and method based on PMU device |
| CN110458330A (en) * | 2019-07-15 | 2019-11-15 | 北京科东电力控制系统有限责任公司 | A kind of cascade section off-the-line optimization method and system based on gray relative analysis method |
| CN110971492A (en) * | 2019-12-23 | 2020-04-07 | 国电南瑞科技股份有限公司 | RTU (remote terminal Unit) full-network simultaneous section generation method and system considering wide-area measurement time delay |
| CN113595058A (en) * | 2021-04-29 | 2021-11-02 | 湖南工学院 | Processing method of power distribution network state estimation measurement data |
| CN113962053A (en) * | 2021-11-08 | 2022-01-21 | 国网江苏省电力有限公司无锡供电分公司 | Power distribution network state evaluation method based on multi-section intelligent instrument data |
| CN114925973A (en) * | 2022-04-13 | 2022-08-19 | 国电南瑞科技股份有限公司 | Power grid running state obtaining method and system based on multi-source data integration |
| CN115184698A (en) * | 2022-07-11 | 2022-10-14 | 国网北京市电力公司 | Synchronous measurement system for power system |
| CN115498651A (en) * | 2022-08-22 | 2022-12-20 | 国网上海市电力公司 | Model prediction dynamic voltage cooperative control method based on PMU measurement data |
| CN116345564A (en) * | 2023-02-17 | 2023-06-27 | 国网山东省电力公司德州供电公司 | Multi-time-scale distributed collaborative optimization scheduling method and system for comprehensive energy system |
-
2023
- 2023-08-30 CN CN202311101913.8A patent/CN117060592B/en active Active
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20080035116A (en) * | 2006-10-18 | 2008-04-23 | 주식회사 세니온 | Automation system and method for automatically detecting phase of electric power distribution line |
| CN102195354A (en) * | 2011-05-16 | 2011-09-21 | 中国电力科学研究院 | Large-grid wide-area asynchronous coordination control method |
| CN102361353A (en) * | 2011-10-26 | 2012-02-22 | 东北电网有限公司 | Method for aligning section raw data on basis of double time scale delay evaluation |
| CN103323742A (en) * | 2013-05-31 | 2013-09-25 | 华北电力大学 | System and method for grid fault identification based on random point distribution PMU |
| CN105914798A (en) * | 2016-04-08 | 2016-08-31 | 李仕乙 | Wind power field reactive power optimization operation control method based on whole network same time section |
| CN106383287A (en) * | 2016-12-05 | 2017-02-08 | 国网辽宁省电力有限公司丹东供电公司 | Method for improving power grid stable section measurement precision through improving data synchronism |
| CN107144760A (en) * | 2016-12-31 | 2017-09-08 | 国电南瑞科技股份有限公司 | It is a kind of to support electric network fault to perceive panorama Wave record method and plateform system with analysis |
| GB201810680D0 (en) * | 2018-06-28 | 2018-08-15 | Synaptec Ltd | Methods and apparatus for making a time-synchronised phasor measurement |
| CN109100614A (en) * | 2018-09-28 | 2018-12-28 | 东北大学 | A kind of transmission open acess system and method based on PMU device |
| CN110458330A (en) * | 2019-07-15 | 2019-11-15 | 北京科东电力控制系统有限责任公司 | A kind of cascade section off-the-line optimization method and system based on gray relative analysis method |
| CN110971492A (en) * | 2019-12-23 | 2020-04-07 | 国电南瑞科技股份有限公司 | RTU (remote terminal Unit) full-network simultaneous section generation method and system considering wide-area measurement time delay |
| CN113595058A (en) * | 2021-04-29 | 2021-11-02 | 湖南工学院 | Processing method of power distribution network state estimation measurement data |
| CN113962053A (en) * | 2021-11-08 | 2022-01-21 | 国网江苏省电力有限公司无锡供电分公司 | Power distribution network state evaluation method based on multi-section intelligent instrument data |
| CN114925973A (en) * | 2022-04-13 | 2022-08-19 | 国电南瑞科技股份有限公司 | Power grid running state obtaining method and system based on multi-source data integration |
| CN115184698A (en) * | 2022-07-11 | 2022-10-14 | 国网北京市电力公司 | Synchronous measurement system for power system |
| CN115498651A (en) * | 2022-08-22 | 2022-12-20 | 国网上海市电力公司 | Model prediction dynamic voltage cooperative control method based on PMU measurement data |
| CN116345564A (en) * | 2023-02-17 | 2023-06-27 | 国网山东省电力公司德州供电公司 | Multi-time-scale distributed collaborative optimization scheduling method and system for comprehensive energy system |
Non-Patent Citations (3)
| Title |
|---|
| CHUNXIU LIU等: "State Estimation of Distribution Network Equipment Based on Genetic Algorithm", 《2023 IEEE 2ND INTERNATIONAL CONFERENCE ON ELECTRICAL ENGINEERING, BIG DATA AND ALGORITHMS (EEBDA)》, pages 929 - 932 * |
| 张武其等: "电力系统惯量评估研究框架与展望", 《中国电机工程学报》, vol. 41, no. 20, pages 6842 - 6852 * |
| 蒋振国等: "计及多扰动的同步发电机模型参数有效性评估", 《电网分析与研究》, vol. 48, no. 7, pages 100 - 105 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117060592B (en) | 2024-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6482699B2 (en) | Use of predictive means in process control systems with wireless or intermittent process measurements | |
| JP6356977B2 (en) | Improved Kalman filter in process control systems | |
| CN101675614A (en) | Method for synchronizing a clock of a network component with a clock of further network component and network component therefor | |
| CN102754369B (en) | For carrying out the method for time synchronized in a communication network | |
| CN103648108A (en) | Sensor network distributed consistency object state estimation method | |
| CN112713881B (en) | Synchronous clock maintaining system and method based on edge calculation | |
| CN111245670A (en) | Method and system for measuring real-time packet loss rate of link | |
| CN108303672A (en) | WLAN indoor positionings error correcting method based on location fingerprint and system | |
| CN117060592B (en) | Multi-domain cooperation-based power grid data synchronous calibration method and system | |
| WO2017109577A1 (en) | Self-organization and self-optimization method and apparatus for an ultra-dense network | |
| CN116974313B (en) | Drying chamber temperature control system suitable for lactobacillus thermal dry production | |
| CN105334734A (en) | Time delay and packet loss compensation method and device of data-based networked control system | |
| CN105607603B (en) | Networked control system and control method capable of compensating time delay and packet loss | |
| CN109788547B (en) | Adaptive time synchronization method for temperature compensation of low communication load | |
| CN105052045A (en) | Apparatus and method for controlling interference in a wireless communication system | |
| Zhang et al. | Two triggered information transmission algorithms for distributed moving horizon state estimation | |
| Hu et al. | Design and implementation of networked predictive control systems based on round trip time delay measurement | |
| Song et al. | Multisensor-based periodic estimation in sensor networks with transmission constraint and periodic mixed storage | |
| CN109525349B (en) | Distributed time synchronization method based on noise estimation and trust weighting | |
| Henriksson | Predictive control for wireless networked systems in process industry | |
| CN112702280A (en) | Network data delay control method | |
| CN107343314B (en) | Time synchronization method and device of wireless sensor network | |
| JPWO2009081712A1 (en) | Wireless communication apparatus and wireless communication method | |
| CN115128961B (en) | Self-triggering model prediction method and system based on LoRaWAN (Long Range network) | |
| Fan et al. | A servo design for slave clock in IEEE 1588 synchronization networks |
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 |