CN108448568B - Power distribution network hybrid state estimation method based on multiple time period measurement data - Google Patents
Power distribution network hybrid state estimation method based on multiple time period measurement data Download PDFInfo
- Publication number
- CN108448568B CN108448568B CN201810190874.6A CN201810190874A CN108448568B CN 108448568 B CN108448568 B CN 108448568B CN 201810190874 A CN201810190874 A CN 201810190874A CN 108448568 B CN108448568 B CN 108448568B
- Authority
- CN
- China
- Prior art keywords
- rtu
- state
- pmu
- pseudo
- data
- 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
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 86
- 239000011159 matrix material Substances 0.000 claims description 81
- 238000002347 injection Methods 0.000 claims description 28
- 239000007924 injection Substances 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 8
- 238000011835 investigation Methods 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 238000009434 installation Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 14
- 230000008859 change Effects 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
The invention discloses a power distribution network mixed state estimation method based on multiple time period measurement data. The proposed method can significantly shorten the period of estimation and provide real-time status information for advanced control applications like voltage control and system congestion control. The method is helpful for detecting the sudden change of the state. The robustness of the algorithm can be further optimized, and the influence of bad data and leverage points on the estimation result can be reduced. In addition, the details of the DKF process may be further optimized.
Description
Technical Field
The invention relates to a power distribution network hybrid state estimation method based on multiple time period measurement data.
Background
The concept of power system state estimation was proposed in the 70 s of the 20 th century. Since then, state estimation has become an indispensable element in EMS systems, which can provide the operator of the power system with the current real-time state information of the system, and control and development of high-level applications of the system are performed on the basis of system state perception. The use of these techniques maximizes the utilization of current system assets and provides safeguards for safe, stable operation of the system. In traditional research, the distribution network is regarded as a single passive network, and the flow of energy is unidirectional only as a link from the transmission network to the users. In addition, the influence of the traditional power distribution network fault on the stability of a large power grid is very limited, so that the attention degree is low. In recent years, with the application of a large number of renewable energy power sources in a power distribution network system, the increase of new loads and the development of demand side response technologies and policies, the functions of a power distribution network are more comprehensive, and the importance of the power distribution network in the energy system is higher and higher. As the importance of power distribution networks has increased, the measurement and communication infrastructure in power distribution network systems has also improved.
The observability of the system is satisfied, certain redundancy is the basis of state estimation research, calculation of observability guarantee can be solved theoretically, and the redundancy of information can improve the estimation quality and facilitate identification and detection of bad data. In a power distribution network system, in order to make up for the deficiency of measured information, many Pseudo measurement information (Pseudo Measurements) are added in the past state estimation research. Over the past few decades, a great deal of literature has been devoted to the establishment of pseudo-metrology information and weight determination. The early pseudo measurement information is derived from load survey information based on an electric charge metering system, several load characteristic curves with typical characteristics are obtained by combining the influences of the load, the season, the time period and other factors on the load, which are obtained in the long-term power system operation process, and the pseudo measurement information set in the calculation is further determined by combining the capacity of each distribution transformer in a power distribution network. The obtained pseudo measurement information has high uncertainty and often has large access to actual load information. For a modern power distribution network system, because the controllability of a user on loads is greatly enhanced, and the number of movable loads such as type electric vehicles is increased, the traditional pseudo measurement information cannot adapt to the operation management requirements of the modern power distribution network.
The energy internet system comprises a primary energy system, an electric power system and an information system, and emphasizes the maximization of comprehensive energy efficiency. From the perspective of the energy internet, the power distribution network system becomes an interface between a large power grid system and a primary energy system and a data system. The energy internet is an information System (Cyber-Physical System) and has an open interface, and power enterprises and various users participate together to ensure that more information and data are reliably processed. The change of the operating conditions of the power system has certain influence on the optimal operating state of the energy system. New measurement technologies, including synchrophasor measurement technologies and advanced measurement systems (AMI), offer great possibilities for better performing state estimation studies.
The synchronous phasor measurement technology provides an accurate time synchronization characteristic by a GPS signal, can accurately measure amplitude and phase information in a power system and has high information uploading frequency. PMUs (power management units) are currently used as terminals for synchrophasor measurement technologies, covering all substations of 500KV and higher voltage levels and most 220KV substations. With the technical development and the continuous reduction of the cost of the Micro-PMU, the PMU is expected to play an important role in the state estimation of the power distribution network. An RTU (Remote Terminal Unit) is used for monitoring, controlling and data acquisition. Has the functions of remote measurement, remote signaling, remote regulation and remote control. The AMI (advanced measurement infrastructure) comprises a smart meter, a data communication network and a data concentrator, and can acquire user or centralized load information and upload the information at certain time intervals (15 minutes or 30 minutes). The application of AMI technology enables the load information of the nodes not to be black boxes any more, and provides more powerful support for the pseudo-measurement information data. The current widespread use of smart meters on the consumer side has led to the power distribution network System being referred to from an Underdetermined System as an Overdetermined System.
Disclosure of Invention
The invention aims to solve the problems and provides a power distribution network mixed state estimation method based on various time period measurement data, which is used for researching a state estimation result with short time limit and high quality obtained by data fusion of various information based on the current situation of current power distribution network system measurement and an information system; on the other hand, by introducing reasonable technology, the problems of excessive investment on the infrastructure level and the maintenance difficulty caused by the excessive investment are avoided.
In order to achieve the purpose, the invention adopts the following technical scheme that:
the invention discloses a power distribution network hybrid state estimation method based on multiple time period measurement data, which comprises the following steps of:
(1) determining topological structure information of the power distribution network system and installation conditions of different measuring devices, including type and quantity information;
(2) determining the data acquisition frequency of each of PMU, RTU and AMI measuring devices in the power distribution network system according to specific manufacturer instructions and system state estimation needs;
(3) determining a system state variable used in calculation according to a topological structure of the system and a load investigation condition;
(4) carrying out load flow calculation according to the load investigation condition, and taking the calculation result as an initial value of a system state variable;
(5) determining an initial prediction model error matrix Q;
(6) predicting the state by using a Holt two-parameter exponential smoothing method to obtain a predicted state variablePrediction of a measured variableAnd predicting state variablesCovariance matrix of
(7) When the actual measurement data of the PMU is updated, whether bad data exist or the running state of the system changes is determined according to the difference value between the predicted value and the actual measurement value of the measurement variable;
if the above condition does not exist, using Kalman filter method to calculate and obtain the corrected state variableAnd estimating a covariance matrix of the state variables
(8) When the RTU measured data is updated, carrying out rectangular coordinate processing on the information directly measured by the RTU by using the result of the last linear dynamic state estimation;
(9) establishing a state estimation objective function according to actual measurement data of PMU and RTU and node injection current or load current equivalent to AMI last time, and calculating system state variable by using a linear least square method The estimated reference state vector is used as the dynamic state estimation of the next moment, and the linear prediction is carried out on the pseudo measurement data of the next RTU data updating moment;
(10) when AMI measurement data is updated, converting the information of load power or node injection power into equivalent node injection current or load current;
(11) and according to actual measurement data of the PMU, the RTU and the AMI, estimating the static linear state of the whole system to obtain a system state variable vector x.
Further, in the step (6), a predicted state variable is obtainedThe method specifically comprises the following steps:
using Holt two-parameter exponential smoothing method to process state variable under k timePerforming calculations with parameters α and β between 0 and 1;
Fk=α·(1+β)·I
predicting state variablesCovariance matrix ofThe method specifically comprises the following steps:
wherein,representing the predicted value of the state variable at time k,for the state variables after correction, matrix FkIs a transition matrix of state variables from time k-1 to time k, gkIs a control variable for calculating the predicted value of the state variable; sk、bkRespectively intermediate amounts;
a covariance matrix representing the predicted values of the state variables at time k,covariance matrix of the estimated state variable after correction.
Further, in the step (6), a predicted value of the measured variable is obtainedThe method specifically comprises the following steps:
the measured variables include: real and imaginary parts of the node voltage; the real part and the imaginary part of the branch current; injecting real and imaginary parts of current into the nodes;
and (c) calling a function H (x) to calculate a measurement variable, wherein H (x) is H x, and H is a Jacobian matrix:
a) for node voltages, the corresponding h (x) function is as follows:
b) for the branch current, the corresponding h (x) function is as follows:
c) for the node injection current, the addition and subtraction of several branch current phases are written according to the KCL definition.
Wherein, Ibr、IbxRespectively the real and imaginary parts, V, of the branch currenti r、Vi xVoltage real and imaginary parts, V, of node i, respectivelys r、Vs xThe real and imaginary parts, r, of the root node voltage, respectivelykIs the resistance of branch k, xkIs the reactance of branch k, /)iRepresenting a path from the root node to node i.
Further, in the step (7), the state change after the correction is calculatedMeasurement ofAnd estimating a covariance matrix of the state variablesThe method specifically comprises the following steps:
wherein R iszkAs a matrix of standard deviations of the measured information, KkReferred to as a Kalman gain matrix, is,in order to correct the state variable after the correction,a covariance matrix which is the state information variable after correction;representing the predicted value of the state variable at time k,a covariance matrix representing the predicted value of the state variable at the time k;
for measuring information zkIs estimated, corresponding toThe covariance matrix of (a) is:Hkis the Jacobian matrix below time k.
Further, in the step (7), determining whether there is bad data or the operation state of the system changes, specifically:
calculating variable vkNormalized form of (a):
wherein,for the predicted value of the measured variable, zkIs the measured value of the measured variable; v. ofk,iAn ith value representing a difference between an actually measured value and a predicted value of the measured variable,representing the difference between the actual and predicted values of the normalized measured variable, RkThe covariance matrix of the measured data and H are the size of a Jacobian matrix;
setting a threshold t, comparing the variablesAnd a threshold value t; if it is notAnd if the threshold value t is smaller than the threshold value t, the system is considered to be in a stable operation state, and the influence of bad data is not considered.
Further, in the step (8), the information directly measured by the RTU is processed by rectangular coordinate, specifically:
wherein,is the real part of the voltage estimated value obtained after rectangular coordinate processing,the imaginary part of the voltage estimation value is obtained after rectangular coordinate processing;respectively the real part and the imaginary part of the node voltage estimated at the last moment,respectively obtaining a real part and an imaginary part of the branch current estimated at the last moment; vm,iNode voltage amplitude, I, measured for RTUm,iAnd directly measuring the obtained branch current amplitude for the RTU.
Further, in the step (9), an objective function of state estimation is established, specifically:
J=(zpmu-hpmu)TWpmu(zpmu-hpmu)+(zrtu-hrtu)TWrtu(zrtu-hrtu)+(zpseudo-hpseudo)TWpseudo(zpseudo-hpseudo)
wherein z ispmuIs a measurement data vector, h, provided by the PMUpmuIs a calculated value corresponding to a PMU data vector, WpmuIs a weight matrix for PMU data; z is a radical ofrtuIs a measurement data vector, h, provided by the RTUrtuIs a calculated value, W, corresponding to a measured data vector provided by the RTUrtuIs a weight matrix corresponding to RTU data; z is a radical ofpseudoIs the last pseudo-measured data vector, hpseudoIs a pseudo-measurement vector, W, calculated from the current state variablespseudoA pseudo metrology data weight matrix.
Further, in the step (9), the linear prediction is performed on the pseudo measurement data at the next RTU data update time, and the specific implementation method is as follows:
wherein, tkFrom tj-1Time tjTaking each RTU data update as an interval between moments; z is a radical ofp,kRTU measurement data corresponding to time k, zp,jRTU measurement data corresponding to time j, TpThe data update period is measured for the RTU.
Further, in the step (10), the active load, reactive load and concentrated load variables of the user directly measured by the AMI are converted into a node injection power form, specifically:
wherein,is the real part of the equivalent current injected by the node,is the imaginary part, P, of the node injected equivalent currentinj,iIs node injection active power, Q, provided by AMIinj,iIs node injected reactive power, V, provided by AMIi rIs the estimated value of the real part of the node voltage, V, obtained at the last momenti xIs the estimated value of the imaginary part of the node voltage obtained at the last moment.
Further, in the step (11), according to measured data of the PMU, the RTU, and the AMI, a full-system static linear state estimation is performed to obtain a system state variable vector x, which specifically includes:
establishing an objective function of state estimation:
J=(zpmu-hpmu)TWpmu(zpmu-hpmu)+(zrtu-hrtu)TWrtu(zrtu-hrtu)+(zpseudo-hpseudo)TWpseudo(zpseudo-hpseudo)
wherein z ispmuIs a measurement data vector, h, provided by the PMUpmuIs a calculated value corresponding to a PMU data vector, WpmuIs a weight matrix for PMU data; z is a radical ofrtuIs a measurement data vector, h, provided by the RTUrtuIs a calculated value, W, corresponding to a measured data vector provided by the RTUrtuIs a weight matrix corresponding to RTU data; z is a radical ofpseudoIs a measurement data vector, h, provided by AMIpseudoIs a calculated value corresponding to a measured data vector provided by AMI, WpseudoA pseudo metrology data weight matrix.
The invention has the beneficial effects that:
the inventive method can significantly shorten the period of estimation and provide real-time status information for advanced control applications like voltage control and system congestion control. The method is helpful for detecting the sudden change of the state. The method can further reduce the influence of bad data and leverage points on the estimation result by optimizing the robustness of the algorithm. In addition, the details of the dkf (discrete Kalman filter) process can be further optimized.
Drawings
FIG. 1 is a schematic diagram of a hybrid power distribution network state estimation method according to the present invention;
FIG. 2 is a single line diagram of an IEEE-69 node system;
fig. 3 is the active injected power at node 38;
FIG. 4 is a graph of the voltage amplitude at node 9;
fig. 5 shows the voltage amplitude at node 5.
The specific implementation mode is as follows:
the invention will be further explained with reference to the drawings.
Kalman Filter principle
The Kalman filter is applied to the dynamic state estimation research of a linear system, and the whole process is divided into two parts: a prediction step and a filtering step. In the prediction link, the numerical value to be predicted and the corresponding covariance matrix are subjected to certain prediction,
wherein Fk-1The matrix connects the state variables of the system at two moments k-1 and k, qkRepresenting process errors in the prediction.
The second step is a filtering estimation part, which updates the state value by using the latest actually measured data.
Wherein the matrix KkIs a Gain Matrix (Gain Matrix). The estimation link of the Kalman filter is essentially a balance compromise of the predicted information and the actually measured information, and the judgment of the weight or the precision of different information is very important here.
The advantage of the Kalman filter is that it gives a reasonable interpretation to determine the weights that various information should have, and on the other hand the Kalman filter also provides more redundant information.
Holt’s Exponential Smoothing Method
Unlike the conventional linear dynamic system, the actual power system is a nonlinear system, and it is difficult to write the state variables (including node voltage vector and branch current vector) as the operating system, and the state variables are daes (differential and alternate equations) in a sequence of front and back time. Therefore, the prediction link needs to be supplemented according to the research in the aspects of statistics and linear regression.
The ARIMA (0,1,0) process model is applied to the research of partial dynamic state estimation, and the reason is that the data updating time interval of PMU is very small. This can be applied well under quasi-steady state condition, in order to capture the change characteristic of the state variable better, the invention uses the holt's linear exponential smoothing technique of two parameters. The method is characterized in that two parts of the current estimated level and the changing trend are considered simultaneously,
Sk+1=αxk+(1-α)(Sk+bk) (7)
bk+1=β(Sk+1-Sk)+(1-β)bk(8)
where the parameters a and β are directly settable, values between 0 and 1, equation (9) gives that the predicted value at time k +1 is given by two parts, i.e. the predicted estimated level and the trend of the change, equation (7) also gives that the estimated level S at time k +1 is divided into two parts, one being the final estimated result x at the previous time kkAnd the other is the predicted quantity S given at the previous momentk+bk。
Hybrid power distribution network state estimation scheme
As mentioned in the first section, there are three different measurement devices in the distribution network system, PMU, RTU and AMI. As shown in fig. 1, AMI in the present invention sets data update every 15 minutes, and if waiting for all the latest measurement information to appear, it will wait for every 15 minutes to perform static state estimation calculation, which does not meet the requirement of state sensing of modern power distribution systems. The RTU measured voltage amplitude and branch current amplitude information, the update frequency is to update data every 20 seconds. PMUs can upload 50 data per second as previously described. The whole process has three parts: only the time period when the PMU data is updated, the time when the PMU and RTU data are updated, and the time when the PMU, RTU and AMI data exist. For the first two cases, the system has a great problem in observability and data redundancy, and needs to add pseudo-metric information.
In the invention, the state estimation of the power distribution network based on branch current is selected, the state variable of the system is selected as the voltage of a root node (connected with a large power grid) and the real part and the imaginary part of each branch current, the calculation is carried out by utilizing a rectangular coordinate system, the linearization of the whole calculation process is convenient,
the relationship between the information measured by the PMUs and the state variables is linear and constant, and the corresponding Jacobian matrix does not need to be modified. For the current vector measurement itself, there is a relationship between 0 and 1 with the state variable, and the corresponding parameters between the voltage measurement information and the state variable are determined entirely by the topology of the system and the parameters of the line or other components.
The information measured by the RTU is a nonlinear functional relationship of state variables, and in order to make the calculation process more efficient, it is necessary to perform rectangular coordinate processing on the measured information. An ideal method is to fuse the amplitude information by using the linear dynamic estimation result to obtain the real part and the imaginary part of the electrical quantity.
In formula (11)Andrespectively, the node electric compaction estimated at the last momentThe real and imaginary parts of the branch current. The variable weights derived from equation (11) also need to be modified according to equation (11). Through rectangular coordinate processing, the Jacobian matrix element corresponding to the RTU measurement information is also fixed and unchangeable, and the content of the element is only 0 or 1.
Variables directly measured by AMI include active load and reactive load of users and concentrated load of a plurality of users, and are presented in a system in a node injection power form, the node injection power can be expressed into a linear function of a branch circuit (information of a Jacobian matrix needs to be updated) and can also be converted into a node injection power form, the latter method is selected in the invention to process according to a formula (12) and process uncertainty characteristics after conversion,
in the formula (12), the variable Pinj,i,And Qinj,iThe information is not available at every moment, and the corresponding node injection current information is obtained through the SE calculation at the previous moment when AMI does not perform data updating. After processing according to equation (12), the corresponding Jacobian matrix partial elements have a relationship only with the node-branch incidence matrix.
The reason for this is that: firstly, all elements of a Jacobian matrix are derived from parameters of 0,1 and a line, and the matrix is kept constant; the condition number of the Jacobian matrix H is reduced.
The update of AMI data is considered as a flag for one completion cycle. When new AMI data information (power injection information derived from the smart meter system) is available, the system completes the state estimation in accordance with the rigorous WLS calculations.
The detailed implementation process of the method for estimating the hybrid state of the power distribution network based on the measurement data of various time periods is described in detail below, and the method specifically comprises the following steps:
(1) determining topological structure information of the power distribution network system and installation conditions of different measuring devices, including type and quantity information;
(2) determining the data acquisition frequency of each of PMU, RTU and AMI measuring devices in the power distribution network system according to specific manufacturer instructions and system state estimation needs;
(3) determining a system state variable used in calculation according to a topological structure of the system and a load investigation condition;
(4) carrying out load flow calculation according to the load investigation condition, and taking the calculation result as an initial value of a system state variable;
(5) determining an initial prediction model error matrix Q;
the prediction model error matrix Q is a diagonal matrix with initial elements set to 1 × 10-6。
(6) Predicting the state by using a Holt two-parameter exponential smoothing method to obtain a predicted state variablePrediction of a measured variableAnd predicting state variablesCovariance matrix of
Using the Holt two-parameter exponential processing (Holt's two-parameter exponential smoothing) method, firstly, the state variable under the k time is processedCalculations were performed with parameters α and β between 0-1.
Fk=α·(1+β)·I
bk=β(Sk-Sk-1)+(1-β)·bk-1
Wherein,representing the predicted values of the state variables at time k, matrix FkIs the state variable transition matrix from time k-1 to time k, gkIs a control variable for calculating the predicted value of the state.
Measurement variable components: real and imaginary parts of the node voltage; the real part and the imaginary part of the branch current; the nodes inject the real and imaginary parts of the current.
The function H (x) is called to calculate the measured variable, and in the linear condition, the calculation methods of H (x) ═ H × x, the function H (x) and the Jacobian matrix H are explained as follows.
For node voltage, the h (x) function is as follows
The corresponding portion of Jacobian matrix H is as follows
For the branch current, the corresponding h (x) function is as follows
The corresponding portion of Jacobian matrix H is as follows
For the node injection current, no matter the element corresponding to the h (x) expression and the Jacobian matrix can be written into a form of addition and subtraction of several branch current phases according to the KCL definition, and details are not repeated.
(7) When the actual measurement data of the PMU is updated, whether bad data exist or the running state of the system changes is determined according to the difference value between the predicted value and the actual measurement value of the measurement variable;
if bad data exist or the operation state of the system changes, a proper method is selected for processing according to the data redundancy condition. Otherwise, calculating to obtain the corrected state variable by using a Kalman filter methodAnd estimating a covariance matrix of the state variables
The method for judging whether bad data exist or not or the operation state of the system changes comprises the following steps:
comparing the obtained variablesAnd a set threshold value t, which is set here by empirical values. If the value is less than the threshold value t, the system is considered to be in a (quasi-) stable operation state, and the influence of bad data is not considered.
Calculating the corrected state variableAnd estimating a covariance matrix of the state variablesThe method specifically comprises the following steps:
wherein R iszkAs a matrix of standard deviations of the measured information, KkReferred to as a Kalman gain matrix, is,is the state variable after correction.
(8) when the RTU measured data is updated, the result of the last linear dynamic state estimation is utilized to carry out rectangular coordinate processing on the information directly measured by the RTU, and the method specifically comprises the following steps:
wherein,is a rectangular seatThe real part of the voltage estimate obtained after the normalization process,the imaginary part of the voltage estimation value is obtained after rectangular coordinate processing;respectively the real part and the imaginary part of the node voltage estimated at the last moment,respectively obtaining a real part and an imaginary part of the branch current estimated at the last moment; vm,iNode voltage amplitude, I, measured for RTUm,iAnd directly measuring the obtained branch current amplitude for the RTU.
(9) Establishing a state estimation objective function according to actual measurement data of PMU and RTU and node injection current or load current equivalent to AMI last time, and calculating system state variable by using a linear least square method The estimated state vector is used as a reference state vector of the dynamic state estimation at the next moment, and the linear prediction is carried out on the pseudo measurement data (node injection current) at the next RTU data updating moment;
establishing an objective function of state estimation, specifically:
J=(zpmu-hpmu)TWpmu(zpmu-hpmu)+(zrtu-hrtu)TWrtu(zrtu-hrtu)+(zpseudo-hpseudo)TWpseudo(zpseudo-hpseudo)
wherein z ispmuIs a measurement data vector, h, provided by the PMUpmuIs a calculated value corresponding to a PMU data vector, WpmuIs a weight matrix for PMU data; z is a radical ofrtuIs RTMeasurement data vector, h, provided by UrtuIs a calculated value, W, corresponding to a measured data vector provided by the RTUrtuIs a weight matrix corresponding to RTU data; z is a radical ofpseudoIs the last pseudo-measured data vector, hpseudoIs a pseudo-measurement vector, W, calculated from the current state variablespseudoA pseudo metrology data weight matrix.
The method for linearly predicting the pseudo measurement data at the next RTU data updating moment comprises the following specific steps:
wherein, tkIs from tj-1Time tjOne time between the moments, with each RTU data update as an interval.
(10) When AMI measurement data is updated, converting the information of load power or node injection power into equivalent node injection current or load current; the method specifically comprises the following steps:
wherein,is the real part of the equivalent current injected by the node,is the imaginary part, P, of the node injected equivalent currentinj,iIs node injection active power, Q, provided by AMIinj,iIs node injected reactive power, V, provided by AMIi rIs the estimated value of the real part of the node voltage, V, obtained at the last momenti xIs the estimated value of the imaginary part of the node voltage obtained at the last moment.
(11) And according to actual measurement data of the PMU, the RTU and the AMI, estimating the static linear state of the whole system to obtain a system state variable vector x.
The specific implementation method comprises the following steps:
the objective function is established as follows:
Min J=(zpmu-hpmu)TWpmu(zpmu-hpmu)+(zrtu-hrtu)TWrtu(zrtu-hrtu)+(zpseudo-hpseudo)TWpseudo(zpseudo-hpseudo)
calculation using WLS (weighted least squares method);
Δzk=zmk-h(xk)
zmk=[zpmu,zrtu,zpseudo]T
h(xk)=[hpmu,hrtu,hpseudo]T
hpmusee step (6);
hrtuthe processing (voltage amplitude and branch current amplitude) of (1) is as follows:
representing the voltage amplitude, delta, of node iiIs the voltage phase of node i, /)iRepresenting a path from the root node to node i.
hpseudoSee step (6);
H=[Hpmu;Hrtu;Hpseudo]。
it should be noted that, the least square method is used in both the step (9) and the step (11), but the corresponding scenarios are different, and the sources of the used measurement information z are not completely the same. Step (9) corresponds to the updating of only PMU and RTU real-time data, and the node injection power or the node injection current is estimated through the previous estimation result; and (11) updating PMU, RTU and AMI data, and directly using real-time data for calculation.
Simulation case analysis
In this section, two simulation analyses were performed using an IEEE 69 node power distribution network system. The first simulation case is used to illustrate the ability of the branch current and DKF based state estimation method to track the state of the system. For better illustration, results obtained with a WLS-based method were also compared together. The second simulation case gives the above-mentioned result of the power distribution network comprehensive state estimation strategy.
A. Testing for traceability
Figure 2 gives a single line diagram of a 69 node system. The IEEE 69 node system is a fully radial distributed three-phase balanced system. The relevant parameters and basic power flow information can be consulted in the existing literature.
To test the ability of the DKF method to track system state variables with a limited number of PMUs, the results of the two methods are compared: the results obtained using the DKF method and the results obtained using the conventional WLS method. Both methods were compared under the same conditions (including the load case).
The measurement error of the node injection power variable in the system is considered to follow a normal distribution law with a mean value of 0. The standard deviation of the distribution is set in such a way that the maximum error does not exceed 20% of the measured value. The results of the active injected power at node 38 are given in fig. 3, and the voltage amplitude at node 9 is given in fig. 4.
The simulation shows 50 different time points, and the result shows that the DKF-based method can obtain a method which is not inferior to the WLS method within a certain time range. It has to be noted that in a real scenario, it is not possible for the system's operational engineer to get the system's node injection power information at every moment. The purpose of this simulation is only to demonstrate the effectiveness and potential applicability of the DKF-based method.
A. Case analysis under mixed measurement information
In this section, the simulation involved three measurement techniques with different time periods. As shown in fig. 1, the state estimation calculation based on DKF is performed every two seconds. Linear WLS (weighted Least Square) state estimation is performed every 20 seconds, using information including real-time measurement information from PMUs and RTUs and pseudo-measurement information for node injection power presentation. WLS calculations based on all actual measurement information were performed every 15 minutes.
The amplitude measurement error of the PMU is set to be less than 1% of the actual value, the measurement error of the phase can be set according to the maximum error of the time synchronization signal, and the time error of the GPS is within 1 microsecond according to the current technical standard. The measurement error of the RTU information is set to be less than 2% of the true value. The error setting for the AMI actual measurement is less than 5% of the true quantity. The above errors are all considered to follow a normal distribution.
Fig. 5 gives information on the magnitude of the voltage at the fifth node in the system. The line with greater undulations represents the estimated value and the substantially smooth line represents the true value. Compared with the real numerical value, the estimated result has an error within 0.2 percent, and meets the requirement of power distribution network state estimation.
Although the embodiments of the present invention have been described with reference to the accompanying drawings, it is not intended to limit the scope of the present invention, and it should be understood by those skilled in the art that various modifications and variations can be made without inventive efforts by those skilled in the art based on the technical solution of the present invention.
Claims (10)
1. The method for estimating the hybrid state of the power distribution network based on the measurement data of various time periods is characterized by comprising the following steps of:
(1) determining topological structure information of the power distribution network system and installation conditions of different measuring devices, including type and quantity information;
(2) determining the data acquisition frequency of each of PMU, RTU and AMI measuring devices in the power distribution network system according to specific manufacturer instructions and system state estimation needs;
(3) determining a system state variable used in calculation according to a topological structure of the system and a load investigation condition;
(4) carrying out load flow calculation according to the load investigation condition, and taking the calculation result as an initial value of a system state variable;
(5) determining an initial prediction model error matrix Q;
(6) predicting the state by using a Holt two-parameter exponential smoothing method to obtain a predicted state variablePrediction of a measured variableAnd predicted state variablesCovariance matrix of
(7) When the actual measurement data of the PMU is updated, whether bad data exist or the running state of the system changes is determined according to the difference value between the predicted value and the actual measurement value of the measurement variable;
if the above condition does not exist, using Kalman filter method to calculate and obtain the corrected state variableAnd estimating a covariance matrix of the state variables
(8) When the RTU measured data is updated, carrying out rectangular coordinate processing on the information directly measured by the RTU by using the result of the last linear dynamic state estimation;
(9) establishing a state estimation objective function according to actual measurement data of PMU and RTU and node injection current or load current equivalent to AMI last time, and calculating system state variable by using a linear least square method The estimated reference state vector is used as the dynamic state estimation of the next moment, and the linear prediction is carried out on the pseudo measurement data of the next RTU data updating moment;
(10) when AMI measurement data is updated, converting the information of load power or node injection power into equivalent node injection current or load current;
(11) and according to actual measurement data of the PMU, the RTU and the AMI, estimating the static linear state of the whole system to obtain a system state variable vector x.
2. The method for estimating hybrid state of distribution network based on multiple time period measurement data according to claim 1, wherein in the step (6), the predicted state variable is obtainedThe method specifically comprises the following steps:
using Holt two-parameter exponential smoothing method to predict state variable under k timePerforming calculations with parameters α and β between 0 and 1;
Fk=α·(1+β)·I
bk=β(Sk-Sk-1)+(1-β)·bk-1
predicting state variablesCovariance matrix ofThe method specifically comprises the following steps:
wherein,for the state variables after correction, matrix FkIs a transition matrix of state variables from time k-1 to time k, gkIs a control variable for calculating the predicted value of the state variable; sk、bkRespectively intermediate amounts;
3. The method for estimating hybrid state of distribution network based on multiple time period measurement data according to claim 1, wherein in the step (6), the predicted values of the measurement variables are obtainedThe method specifically comprises the following steps:
the measured variables include: real and imaginary parts of the node voltage; the real part and the imaginary part of the branch current; injecting real and imaginary parts of current into the nodes;
and (c) calling a function H (x) to calculate a measurement variable, wherein H (x) is H x, and H is a Jacobian matrix:
for node voltages, the corresponding h (x) function is as follows:
for the branch current, the corresponding h (x) function is as follows:
for node injection current, writing a mode of adding and subtracting a plurality of branch current phases according to KCL definition;
wherein,Ibr、Ibxrespectively the real and imaginary parts, V, of the branch currenti r、Vi xRespectively the real and imaginary parts of the voltage at node i,the real and imaginary parts, r, of the root node voltage, respectivelykIs the resistance of branch k, xkIs the reactance of branch k, /)iRepresenting a path from the root node to node i.
4. The method for estimating hybrid state of distribution network based on multiple time period measurement data according to claim 1, wherein in the step (7), the state variable after being corrected is calculatedAnd estimating a covariance matrix of the state variablesThe method specifically comprises the following steps:
wherein R iszkAs a matrix of standard deviations of the measured information, KkIn the form of a Kalman gain matrix, the matrix,in order to correct the state variable after the correction,a covariance matrix which is the state information variable after correction;representing the predicted value of the state variable at time k,a covariance matrix representing the predicted value of the state variable at the time k;
5. The method for estimating the hybrid state of the power distribution network based on the measurement data of multiple time periods as claimed in claim 1, wherein in the step (7), it is determined whether there is bad data or the operation state of the system changes, specifically:
calculating variable vkNormalized form of (a):
wherein,for measuring variablesPredicted value, zkIs the measured value of the measured variable; v. ofk,iAn ith value representing a difference between an actually measured value and a predicted value of the measured variable,representing the difference between the actual and predicted values of the normalized measured variable, RkThe covariance matrix of the measured data and H are the size of a Jacobian matrix;
6. The method for estimating the hybrid state of the power distribution network based on the measurement data of multiple time periods as claimed in claim 1, wherein in the step (8), the information directly measured by the RTU is processed by rectangular coordinate, specifically:
wherein,is the real part of the voltage estimated value obtained after rectangular coordinate processing,the imaginary part of the voltage estimation value is obtained after rectangular coordinate processing;respectively the real part and the imaginary part of the node voltage estimated at the last moment,respectively obtaining a real part and an imaginary part of the branch current estimated at the last moment; vm,iNode voltage amplitude, I, measured for RTUm,iAnd directly measuring the obtained branch current amplitude for the RTU.
7. The method for estimating the hybrid state of the power distribution network based on the measurement data of the multiple time periods as claimed in claim 1, wherein in the step (9), an objective function of state estimation is established, specifically:
J=(zpmu-hpmu)TWpmu(zpmu-hpmu)+(zrtu-hrtu)TWrtu(zrtu-hrtu)+(zpseudo-hpseudo)TWpseudo(zpseudo-hpseudo)
wherein z ispmuIs a measurement data vector, h, provided by the PMUpmuIs a calculated value corresponding to a PMU data vector, WpmuIs a weight matrix for PMU data; z is a radical ofrtuIs a measurement data vector, h, provided by the RTUrtuIs a calculated value, W, corresponding to a measured data vector provided by the RTUrtuIs a weight matrix corresponding to RTU data; z is a radical ofpseudoIs the last pseudo-measured data vector, hpseudoIs a pseudo-measurement vector, W, calculated from the current state variablespseudoIs a pseudo metrology data weight matrix.
8. The method for estimating the hybrid state of the power distribution network based on the measurement data of multiple time periods as claimed in claim 1, wherein in the step (9), the pseudo measurement data at the next RTU data update time is linearly predicted, and the specific implementation method is as follows:
wherein, tkFrom tj-1Time tjTaking each RTU data update as an interval between moments; z is a radical ofp,kRTU measurement data corresponding to time k, zp,jRTU measurement data corresponding to time j, TpThe data update period is measured for the RTU.
9. The method for estimating the hybrid state of the power distribution network based on the measurement data of multiple time periods as claimed in claim 1, wherein in the step (10), the user active, reactive load and concentrated load variables directly measured by the AMI are converted into a form of node injection power, specifically:
wherein,is the real part of the equivalent current injected by the node,is the imaginary part, P, of the node injected equivalent currentinj,iIs node injection active power, Q, provided by AMIinj,iIs node injected reactive power, V, provided by AMIi rIs the estimated value of the real part of the node voltage, V, obtained at the last momenti xIs the estimated value of the imaginary part of the node voltage obtained at the last moment.
10. The method for estimating hybrid state of power distribution network based on multiple time period measurement data according to claim 1, wherein in step (11), a system-wide static linear state estimation is performed according to measured data of PMU, RTU, and AMI to obtain a system state variable vector x, specifically:
establishing an objective function of state estimation:
J=(zpmu-hpmu)TWpmu(zpmu-hpmu)+(zrtu-hrtu)TWrtu(zrtu-hrtu)+(zpseudo-hpseudo)TWpseudo(zpseudo-hpseudo)
wherein z ispmuIs a measurement data vector, h, provided by the PMUpmuIs a calculated value corresponding to a PMU data vector, WpmuIs a weight matrix for PMU data; z is a radical ofrtuIs a measurement data vector, h, provided by the RTUrtuIs a calculated value, W, corresponding to a measured data vector provided by the RTUrtuIs a weight matrix corresponding to RTU data; z is a radical ofpseudoIs a measurement data vector, h, provided by AMIpseudoIs a calculated value corresponding to a measured data vector provided by AMI, WpseudoIs the AMI data weight matrix.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810190874.6A CN108448568B (en) | 2018-03-08 | 2018-03-08 | Power distribution network hybrid state estimation method based on multiple time period measurement data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810190874.6A CN108448568B (en) | 2018-03-08 | 2018-03-08 | Power distribution network hybrid state estimation method based on multiple time period measurement data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108448568A CN108448568A (en) | 2018-08-24 |
CN108448568B true CN108448568B (en) | 2020-03-27 |
Family
ID=63193832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810190874.6A Active CN108448568B (en) | 2018-03-08 | 2018-03-08 | Power distribution network hybrid state estimation method based on multiple time period measurement data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108448568B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109711662B (en) * | 2018-11-14 | 2021-03-09 | 华北电力大学 | Power grid robust state estimation method based on multi-source data fusion |
CN109586289B (en) * | 2018-12-13 | 2020-12-29 | 国网山东省电力公司青岛供电公司 | Power distribution network multi-time scale recursive dynamic state estimation method and system |
CN109829613A (en) * | 2018-12-26 | 2019-05-31 | 国网山东省电力公司淄博供电公司 | A kind of power distribution network multiple time scale model method for estimating state and system |
CN109754013B (en) * | 2018-12-31 | 2021-09-10 | 浙江大学 | Electric power system hybrid measurement fusion method based on unscented Kalman filtering |
CN109921415B (en) * | 2019-03-13 | 2022-08-02 | 东北大学 | Hybrid measurement-oriented online defense method for power grid malignant data injection attack |
CN111198891B (en) * | 2019-12-05 | 2023-02-07 | 国网甘肃省电力公司电力科学研究院 | Data source fusion method, electronic device and non-transitory computer readable storage medium |
CN111898816A (en) * | 2020-07-21 | 2020-11-06 | 南京信息工程大学 | Dynamic state estimation method for comprehensive energy system |
CN112083287B (en) * | 2020-09-08 | 2021-06-01 | 上海交通大学 | Power distribution network dynamic event positioning method and system based on characteristic value spectrum distribution model |
CN113488989B (en) * | 2021-06-09 | 2023-08-01 | 河海大学 | Local hybrid linear state estimation method based on stream calculation |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103606113A (en) * | 2013-12-02 | 2014-02-26 | 东北石油大学 | Static state estimation method for electrical power system based on PMU device |
CN106921156B (en) * | 2015-12-25 | 2019-11-01 | 中国电力科学研究院 | A kind of active distribution network method for estimating state based on more sampling period hybrid measurements |
CN106374497A (en) * | 2016-10-20 | 2017-02-01 | 天津大学 | Integrated state estimation method of middle-low-voltage hybrid distribution network based on AMI (Advanced Metering Infrastructure) |
CN106707061A (en) * | 2016-12-16 | 2017-05-24 | 湖南大学 | Hybrid measurement based power distribution network dynamic state estimation method |
-
2018
- 2018-03-08 CN CN201810190874.6A patent/CN108448568B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN108448568A (en) | 2018-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108448568B (en) | Power distribution network hybrid state estimation method based on multiple time period measurement data | |
Džafić et al. | Real time estimation of loads in radial and unsymmetrical three-phase distribution networks | |
Della Giustina et al. | Electrical distribution system state estimation: measurement issues and challenges | |
CN106921156B (en) | A kind of active distribution network method for estimating state based on more sampling period hybrid measurements | |
CN110880760A (en) | Low-voltage distribution network line loss and reactive compensation collaborative optimization method | |
CN110299762B (en) | PMU (phasor measurement Unit) quasi-real-time data-based active power distribution network robust estimation method | |
CN113078630B (en) | Low-voltage distribution network topology identification method based on real-time measurement data | |
Jabr et al. | Distribution management systems for smart grid: Architecture, work flows, and interoperability | |
CN106655152B (en) | A kind of State Estimation for Distribution Network measuring characteristic based on AMI | |
CN103199528A (en) | Status estimating and coordinating method of wide-area power system | |
Picallo et al. | A two-step distribution system state estimator with grid constraints and mixed measurements | |
CN106372440B (en) | A kind of adaptive robust state estimation method of the power distribution network of parallel computation and device | |
CN109858061B (en) | Power distribution network equivalence and simplification method for voltage power sensitivity estimation | |
Wang et al. | Research on time-sharing ZIP load modeling based on linear BP network | |
Wu et al. | Research on improvement of line loss algorithm based on three-phase unbalance degree | |
Liu et al. | Dynamic equivalent modeling for microgrids based on LSTM recurrent neural network | |
CN109586289A (en) | A kind of power distribution network Multiple Time Scales recurrence dynamic state estimator method and system | |
CN113468729B (en) | Method and system for measuring and calculating operation errors of electric power station metering device without calibration mode | |
CN116094060A (en) | Main and distribution network integrated state estimation method based on multi-source measurement | |
CN112103956B (en) | Distribution network state estimation method based on intelligent electric meter dynamic measurement point | |
US20220329101A1 (en) | Method and system for estimation and monitoring of distributed network conditions | |
CN112749489A (en) | Power grid model generation method based on splicing technology | |
CN111709612A (en) | Power distribution network state estimation method considering collected historical data | |
Yi et al. | A line loss calculation method for arbitrary partition of power grid based on measurement resource and matpower | |
CN103683278A (en) | Power distribution network load calibration method applied to on-line loop closing trend optimization analysis |
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 |