CN103729537A - Extension cone method capable of sensing real-time situation of power distribution network - Google Patents
Extension cone method capable of sensing real-time situation of power distribution network Download PDFInfo
- Publication number
- CN103729537A CN103729537A CN201310412471.9A CN201310412471A CN103729537A CN 103729537 A CN103729537 A CN 103729537A CN 201310412471 A CN201310412471 A CN 201310412471A CN 103729537 A CN103729537 A CN 103729537A
- Authority
- CN
- China
- Prior art keywords
- node
- distribution network
- power
- power distribution
- real
- 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.)
- Pending
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses an extension cone method capable of sensing the real-time situation of a power distribution network. According to the extension cone method, an extension cone planning model is established, real-time trend results such as the voltage amplitude value and the phase angle of nodes of the whole network are calculated, and the purpose of sensing the real-time situation of the power distribution network is achieved. The extension cone method comprises the following steps that first, net rack data are read, the net rack data comprise the connection relation of the nodes and branch circuits, real-time data returned by points which are provided with measuring devices and historical loads of points which are not provided with the measuring devices; second, a planning mathematical model is established, the planning mathematical model comprises an objective function model, a trend equation equality constraint, an inequality constraint of node voltages and load nodes; third, secondary cone planning processing is conducted on the trend equation equality constraint, the node voltages, node active power and node reactive power respectively, and state estimation models of an objective function, the equality constraint and the inequality constraint are obtained; fourth, the state estimation models are solved, so that the voltage amplitude of the nodes of the power distribution network and the trend information of the phase angle are obtained, and the purpose for sensing the real-time situation of the power distribution network is achieved.
Description
Technical field
The present invention relates to Steady-State Analysis of Power System technical field, relate in particular to the computing method of a kind of expansion of the employing for perception power distribution network real-time situation cone planning.
Background technology
Along with further developing of intelligent grid, distributed energy progressively accesses in power distribution network, for the operation of power distribution network and planning bring huge challenge, therefore, find a kind of exact method for perception power distribution network real-time situation and seem particularly important.
In electric power networks, power distribution network measuring equipment need to be arranged along the line at feeder line, site environment complexity, and maintenance cost is higher.In addition, the grid structure of power distribution network changes very fast, newly-built many with reconstruction circuit, needs timely install and adjust measuring equipment, due to cost and workload restriction, in fact all only at part key node, measuring equipment is installed.Therefore just determined that power distribution network has the low feature of measurement amount redundance, the more accurate service data that how to obtain other node is a great problem of perception power distribution network real-time situation always.The method of traditional perception power distribution network situation has substitutional resistance method, rms current method, and the model that these two kinds of methods adopt is fairly simple, does not use detailed service data, and a large amount of metric data can not get application, so computational accuracy is very poor.
Open day is on 01 27th, 2010, publication number is that the patent documentation of CN101635456A discloses the technical scheme that name is called a kind of transmission and distribution network united state method of estimation, it includes following steps, first search for power distribution network, each feeder line load is equivalent to each outlet switch, according to each outlet switch of power distribution network, power transmission network is carried out to above state estimation at least three times again, then calculate the symmetrical three-phase outlet load value of each outlet switch of power distribution network and the asymmetric outlet load value of three-phase of each outlet switch of power distribution network, and then calculate the equivalent binding occurrence of all switches on power transmission network bus, and using result of calculation as the bus of power transmission network state estimation next time, inject load unbalanced amount and participate in power transmission network state estimation next time, until it converges to setting binding occurrence.The weak point of this technical scheme is, the one, and need to repeatedly estimate just and can draw estimated value, so real-time is bad; The 2nd, can only reduce the upper trend error in transmission and distribution network boundary in traditional evaluation method, can not improve evaluation of integrals precision in state of electric distribution network estimation.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of expansion cone method of perception power distribution network real-time situation, by foundation, expand cone plan model, calculate the Real-time Power Flow results such as the voltage magnitude of the whole network node and phase angle, having reached the object of perception power distribution network real-time situation, is a kind of computing method of comparatively accurate perception power distribution network real-time situation.
The present invention is directed to prior art problem and mainly by following technical proposals, solved, a kind of expansion cone method of perception power distribution network real-time situation, comprises the following steps:
Step 1, reads in rack data, comprises the annexation of node and branch road, the real time data that the point of measuring equipment passes back, the historical load that the point of measuring equipment is not installed are installed, and the point that measuring equipment is not installed also claims pseudo-gauge point;
Step 2, sets up mathematics for programming model, comprises the inequality constrain of objective function model, power flow equation equality constraint, node voltage and load bus, and objective function model is
,
Wherein,
for state variable,
,
for the amplitude of node voltage,
for node voltage phase angle,
for amount measurement
weight,
for measurement amount
measuring value,
for measuring function, power flow equation equality constraint is
,
Wherein,
,
for the set of all nodes,
be respectively
the first and end voltage magnitude of period branch road,
for
the first and end voltage phase angle of period branch road poor,
,
meritorious, idle the exerting oneself that represent a node,
,
represent load bus
meritorious and reactive power,
,
for node
between real part and the imaginary part of transadmittance, node voltage inequality constrain is to set node voltage near ratings, expression formula is
,
Load bus inequality constrain comprises the inequality constrain of node active power and the inequality constrain of node reactive power, the inequality constrain of node active power is that the historical load data by node obtain average load, be multiplied by respectively predetermined coefficients again and obtain upper and lower limit data, expression formula is
,
Wherein,
for the set of load bus, the inequality constrain of node reactive power is to obtain average reactive power by historical load data, then is multiplied by respectively predetermined coefficients and obtains upper and lower limit data, and expression formula is
;
Step 3, respectively power flow equation equality constraint and node voltage, node active power, node reactive power are carried out to second-order coneprogram processing, draw the state aware model that objective function model, equality constraint model and inequality constrain form, objective function model is
,
Wherein,
, equality constraint model is
,
Inequality constrain is
;
Step 4, solving state appraising model, calculates the load of pseudo-gauge point, and then obtains power distribution network node voltage amplitude and phase angle Power Flow Information, realizes perception power distribution network real-time situation.
This method, first read in rack data, the annexation that comprises node and branch road, meritorious and the reactive power of real-time measurement point, node voltage upper, lower limit, the historical load curve of pseudo-gauge point, set up expansion cone plan model, be about to objective function, power flow equation constraint, gain merit/reactive power the constraint specification of node voltage and load bus becomes expansion second-order coneprogram model, thereby the non-zero entry of the gloomy battle array in sea in solution procedure is only concentrated on and rotated in quadric cone and trigonometric function equation, reduced computation complexity, processing to non-linear trend equation is more convenient, and there is equally the good characteristic that solves with model, and can guarantee the Global Optimality separated, after solving state appraising model, can try to achieve the meritorious load or burden without work of pseudo-gauge point, and then calculate the Real-time Power Flow results such as the voltage magnitude of the whole network node and phase angle, reach the object of perception power distribution network real-time situation.
Power distribution network in actual motion must meet power flow equation constraint, adds this constraint condition to be used for calculating the data message of pseudo-gauge point in state estimation, has improved rationality and the estimated accuracy estimated; The inequality constrain of node voltage and meritorious/reactive power can be clamped down in the reasonable scope departing from the quantity of state that actual value is far away, rationally bad data is processed, and has avoided the impact on all the other estimators.
Wherein second-order coneprogram is a kind of special convex programming, and it is widely used in the aspects such as control, finance, Combinatorial Optimization and engineering as an important branch in mathematical programming field.Second-order coneprogram asks the minimum value of a linear objective function on the cartesian product of limited quadric cone and the common factor of affine subspace, and its constraint is non-linear and protruding, belongs to convex programming problem.Its canonical form is as follows:
Wherein, variable
, constant
,
,
,
be
the quadric cone of dimension.
Common quadric cone has following two kinds:
Quadric cone:
Rotation quadric cone:
Second-order coneprogram is between linear programming and Semidefinite Programming, and three has much similar character, as duality theory, and available interior some Algorithm for Solving etc.
According to the feature of traditional nonlinear power system model, by nonlinear programming to the conversion of second-order coneprogram mainly by introducing new variable, original trend equality constraint is become to the combination of linear equality constraints, the constraint of rotation quadric cone and arc tangent constraint.
The variable-definition of wherein introducing is as follows:
By
,
definition known
,
, be the equation number in balance equation group and variable number, for every branch road, draw following statement:
These two approximately intrafascicular, equation (4) does not meet strict rotation quadric cone form, equation (5) is nonlinear arctan function, after model deformation, does not meet strict second-order coneprogram form, therefore be called expansion cone plan model.
Expansion cone plan model after distortion can only concentrate on the non-zero entry of the gloomy battle array in sea in solution procedure and rotate in quadric cone and trigonometric function equation, reduced computation complexity, processing to non-linear trend equation is more convenient, and has equally the good characteristic that solves with former Nonlinear programming Model.
As preferably, this expansion cone of state aware model in adopting described in some cone Algorithm for Solving step 3.That interior some cone algorithm has is insensitive to initial value, convergence rapidly, the advantage such as strong robustness.After solving, can calculate the load of pseudo-gauge point, make it more to approach actual value, and obtain the system load flow results such as each node voltage amplitude of power distribution network and phase angle.
As preferably, in step 2, node voltage inequality constrain is to set node voltage near ratings, in the supply voltage deviation of 20kV and following power distribution network, is nominal voltage
, expression formula is
.The higher limit of node voltage is that 7% of node nominal voltage adds node nominal voltage, node voltage lower limit be node nominal voltage cut node nominal voltage 7%.
As preferably, in step 2, the inequality constrain of node active power is that historical load data by node obtain average active power, then presses
deviation obtains active power upper and lower limit data, and expression formula is
, wherein,
for the set of load bus, the inequality constrain of node reactive power is to obtain average reactive power by historical load data, then presses
deviation obtains reactive power upper and lower limit data, and expression formula is
.The inequality constrain of node active power is that the historical load data by node obtain average load, by average load+20% deviation, obtains active power higher limit, by average load-20% deviation, obtains active power lower limit; The inequality constrain of node reactive power is to obtain average reactive power by historical load data, by reactive power+20% deviation, obtains reactive power higher limit, by reactive power-20% deviation, obtains reactive power lower limit.
The beneficial effect that the present invention brings is, in Steady-State Analysis of Power System, adopt the expansion cone plan model after distortion that the non-zero entry of the gloomy battle array in sea in solution procedure is only concentrated on and rotated in quadric cone and trigonometric function equation, reduced computation complexity, processing to non-linear trend equation is more convenient, and compare with former Nonlinear programming Model and there is equally the good characteristic that solves, thereby conveniently calculate the Real-time Power Flow results such as the voltage magnitude of the whole network node and phase angle, simultaneously owing to adopting plan model to calculate the data message of pseudo-gauge point, therefore rationality and the estimated accuracy of result of calculation have been improved, reached the object of accurate perception power distribution network real-time situation.
Accompanying drawing explanation
Fig. 1 is operational flowchart of the present invention.
Embodiment
Below by embodiment, and by reference to the accompanying drawings, the technical scheme of invention is further described in detail.
Embodiment: as shown in Figure 1, the present invention is to provide a kind of expansion cone method of perception power distribution network real-time situation, is a kind of distribution power system load flow calculation method based on state estimation, and concrete steps are as follows:
Step 1: reading data, comprises the annexation of node and branch road, the real time data that installation measuring point is passed back, the historical load that measurement mounted point (also claiming pseudo-gauge point) is not installed.
Step 2: set up mathematics for programming model, comprise the inequality constrain of objective function, power flow equation equality constraint, voltage and load;
Objective function is:
Wherein,
for state variable,
,
for node voltage phase angle;
for amount measurement
weight;
for measurement amount
measuring value;
for measuring function.
Power flow equation equality constraint is:
（7）
Wherein,
,
set for all nodes;
be respectively
the first and end voltage magnitude of period branch road;
for
the first and end voltage phase angle of period branch road poor;
,
meritorious, idle the exerting oneself that represent a node;
,
represent load bus
meritorious and reactive power;
,
for node
between real part and the imaginary part of transadmittance.
Inequality constrain comprises following a few part:
Node voltage inequality constrain, sets node voltage near ratings, and the supply voltage deviation of general 20kV and following power distribution network is nominal voltage
, expression formula is:
The inequality constrain of load bus active power, is that the historical load data by node obtain average load, is multiplied by respectively the data that certain coefficient obtains upper and lower limit, generally gets
, expression formula is:
Wherein,
set for load bus.
The expression formula of load bus reactive power inequality constrain is:
Above-mentioned model description be a nonlinear programming problem with equation and inequality constrain.
Step 3: respectively power flow equation equality constraint and node voltage, node active power, node reactive power are carried out to second-order coneprogram processing, obtaining based on expansion cone plan model is state aware model, be intended to nonlinear programming problem to be converted into undemanding convex programming problem, making in solution procedure the non-zero entry of the gloomy battle array in sea only concentrate on rotates in quadric cone and trigonometric function equation, to reduce computation complexity.
Introducing variable-definition is as follows:
Thereby power flow equation equality constraint is become:
And two due to the equality constraint of introducing variable and increasing:
Thereby final mask comprises:
Objective function model is:
(14)
Equality constraint model is
(15)
Inequality constrain is:
(16)
Step 4: in adopting, some cone algorithm solves (14-16) this expansion cone.That interior some cone algorithm has is insensitive to initial value, convergence rapidly, the advantage such as strong robustness.After solving, can calculate the load of pseudo-gauge point, make it more to approach actual value, and obtain the system load flow results such as each node voltage amplitude of power distribution network and phase angle, reach the object of accurate perception power distribution network situation.
Claims (4)
1. the expansion cone method of a perception power distribution network real-time situation, is characterized in that comprising the following steps:
Step 1, reads in rack data, comprises the annexation of node and branch road, the real time data that the point of measuring equipment passes back, the historical load that the point of measuring equipment is not installed are installed, and the point that measuring equipment is not installed also claims pseudo-gauge point;
Step 2, sets up mathematics for programming model, comprises the inequality constrain of objective function model, power flow equation equality constraint, node voltage and load bus, and described objective function model is
,
Wherein,
for state variable,
,
for the amplitude of node voltage,
for node voltage phase angle,
for amount measurement
weight,
for measuring
measuring value,
for measuring function, described power flow equation equality constraint is
,
Wherein,
,
for the set of all nodes,
be respectively
the first and end voltage magnitude of period branch road,
for
the first and end voltage phase angle of period branch road poor,
,
meritorious, idle the exerting oneself that represent a node,
,
represent load bus
meritorious and reactive power,
,
for node
between real part and the imaginary part of transadmittance, described node voltage inequality constrain is to set node voltage near ratings, expression formula is
,
Described load bus inequality constrain comprises the inequality constrain of node active power and the inequality constrain of node reactive power, the inequality constrain of described node active power is that the historical load data by node obtain average load, be multiplied by respectively predetermined coefficients again and obtain upper and lower limit data, expression formula is
,
Wherein,
for the set of load bus, the inequality constrain of described node reactive power is to obtain average reactive power by historical load data, then is multiplied by respectively predetermined coefficients and obtains upper and lower limit data, and expression formula is
;
Step 3, respectively power flow equation equality constraint and node voltage, node active power, node reactive power are carried out to second-order coneprogram processing, draw the state aware model that objective function model, equality constraint model and inequality constrain form, described objective function model is
,
Wherein,
, described equality constraint model is
,
Described inequality constrain is
;
Step 4, solves described state estimation model, calculates the load of pseudo-gauge point, and then obtains power distribution network node voltage amplitude and phase angle Power Flow Information, realizes perception power distribution network real-time situation.
2. a kind of expansion cone method of perception power distribution network real-time situation according to claim 1, is characterized in that, this expansion cone of state aware model in adopting described in some cone Algorithm for Solving step 3.
3. a kind of expansion cone method of perception power distribution network real-time situation according to claim 1, it is characterized in that, the inequality constrain of node voltage described in step 2 is to set node voltage near ratings, in the supply voltage deviation of 20kV and following power distribution network, is nominal voltage
, expression formula is
.
4. according to the expansion cone method of a kind of perception power distribution network real-time situation described in claim 1 or 2 or 3, it is characterized in that, the inequality constrain of the active power of node described in step 2 is that the historical load data by node obtain average active power, then presses
deviation obtains active power upper and lower limit data, and expression formula is
, wherein,
for the set of load bus, the inequality constrain of described node reactive power is to obtain average reactive power by historical load data, then presses
deviation obtains reactive power upper and lower limit data, and expression formula is
.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310412471.9A CN103729537A (en) | 2013-09-12 | 2013-09-12 | Extension cone method capable of sensing real-time situation of power distribution network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310412471.9A CN103729537A (en) | 2013-09-12 | 2013-09-12 | Extension cone method capable of sensing real-time situation of power distribution network |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103729537A true CN103729537A (en) | 2014-04-16 |
Family
ID=50453610
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310412471.9A Pending CN103729537A (en) | 2013-09-12 | 2013-09-12 | Extension cone method capable of sensing real-time situation of power distribution network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103729537A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930924A (en) * | 2016-04-15 | 2016-09-07 | 中国电力科学研究院 | Power distribution network situation sensing method based on complex event processing technology and decision tree |
| CN105956760A (en) * | 2016-04-27 | 2016-09-21 | 河海大学 | Intelligent power distribution network situation perception method Based on multivariate spatio-temporal information modeling |
| CN107230982A (en) * | 2017-07-07 | 2017-10-03 | 广西大学 | A kind of micro-capacitance sensor linearizes tidal current computing method |
| CN118249519A (en) * | 2024-05-28 | 2024-06-25 | 国网江西省电力有限公司电力科学研究院 | State sensing method and system for distributed new energy access active power distribution network |
-
2013
- 2013-09-12 CN CN201310412471.9A patent/CN103729537A/en active Pending
Non-Patent Citations (7)
| Title |
|---|
| RABIH A. JABR等: ""Optimal Power Flow Using an Extended Conic Quadratic Formulation"", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
| 帅军庆: "《特大型电网高级调度中心关键技术》", 30 June 2010, 中国电力出版社 * |
| 李宏仲等: "《地区电网无功补偿与电压无功控制》", 30 April 2012, 机械工业出版社 * |
| 王继鑫: ""配电网状态估计的研究及应用"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 白晓清等: ""求解最优潮流问题的内点半定规划法"", 《中国电机工程学报》 * |
| 范金月等: ""最优潮流的扩展二次锥规划模型"", 《豆丁网，网址：HTTP://WWW.DOCIN.COM/P-515303396.HTML》 * |
| 邓超志: ""配电网状态估计及其应用技术研究与仿真试验"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930924A (en) * | 2016-04-15 | 2016-09-07 | 中国电力科学研究院 | Power distribution network situation sensing method based on complex event processing technology and decision tree |
| CN105956760A (en) * | 2016-04-27 | 2016-09-21 | 河海大学 | Intelligent power distribution network situation perception method Based on multivariate spatio-temporal information modeling |
| CN105956760B (en) * | 2016-04-27 | 2019-12-06 | 河海大学 | Intelligent power distribution network situation perception method based on multivariate time-space information modeling |
| CN107230982A (en) * | 2017-07-07 | 2017-10-03 | 广西大学 | A kind of micro-capacitance sensor linearizes tidal current computing method |
| CN107230982B (en) * | 2017-07-07 | 2021-02-02 | 广西大学 | Micro-grid linear load flow calculation method |
| CN118249519A (en) * | 2024-05-28 | 2024-06-25 | 国网江西省电力有限公司电力科学研究院 | State sensing method and system for distributed new energy access active power distribution network |
| CN118249519B (en) * | 2024-05-28 | 2024-09-10 | 国网江西省电力有限公司电力科学研究院 | State sensing method and system for distributed new energy access active power distribution network |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Beerten et al. | Analysis of power sharing and voltage deviations in droop-controlled DC grids | |
| Yang et al. | LMP revisited: A linear model for the loss-embedded LMP | |
| Liu et al. | State estimation of three-phase four-conductor distribution systems with real-time data from selective smart meters | |
| Dhople et al. | Analysis of power system dynamics subject to stochastic power injections | |
| Chen et al. | Optimal meter placement for distribution network state estimation: A circuit representation based MILP approach | |
| CN106532711B (en) | Change the Newton load flow calculation method of Jacobian matrix with iteration and node type | |
| Dörfler et al. | Novel insights into lossless AC and DC power flow | |
| CN102723721A (en) | Power system reactive power optimization method based on individual optimal position self-adaptive variation disturbance particle swarm algorithm | |
| CN107577870A (en) | The distribution network voltage power sensitivity robust estimation method measured based on synchronized phasor | |
| CN103729537A (en) | Extension cone method capable of sensing real-time situation of power distribution network | |
| Oluic et al. | Methodology for rotor angle transient stability assessment in parameter space | |
| CN104156542B (en) | It is a kind of based on the active distribution system Simulation of stability method implicitly projected | |
| CN102420427A (en) | Method for calculating voltage stability margin of regional power grid by considering outer network equivalence | |
| CN104218569A (en) | Evaluative analysis method for static security check of large-scaled power grid | |
| CN104993494A (en) | Motor simulator based on four-quadrant power electronic converter and method | |
| CN106026086A (en) | Power grid operation state dynamic estimation method | |
| Cheng et al. | An online feedback-based linearized power flow model for unbalanced distribution networks | |
| CN102522742B (en) | Estimation method of external network Thevenin equivalent parameters based on single point measurement information | |
| CN106410811B (en) | Iteration small impedance branches endpoint changes the tidal current computing method of Jacobian matrix for the first time | |
| CN103246207B (en) | A kind of online power-less optimized controlling method based on real-time emulation system | |
| Joseph et al. | Unbiased minimum variance filter-based generator state estimation using PMU measurements for unknown generator input | |
| Selim et al. | Efficient time series simulation of distribution systems with voltage regulation and PV penetration | |
| Jabr | Power flow based volt/var optimization under uncertainty | |
| Machowski et al. | State-variable control of shunt FACTS devices using phasor measurements | |
| Ren et al. | Steady-state Security Region Boundary Modification Model: A Hybrid Physical Model-Driven and Data-Driven Approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140416 |