CN107732923A - The smoothing of discrete magnitude coherency function mixes secondary voltage control solving method with continuous quantity - Google Patents
The smoothing of discrete magnitude coherency function mixes secondary voltage control solving method with continuous quantity Download PDFInfo
- Publication number
- CN107732923A CN107732923A CN201710846406.5A CN201710846406A CN107732923A CN 107732923 A CN107732923 A CN 107732923A CN 201710846406 A CN201710846406 A CN 201710846406A CN 107732923 A CN107732923 A CN 107732923A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- discrete
- mtd
- voltage control
- 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
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
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
-
- 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/30—Reactive power compensation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of discrete magnitude coherency function smoothing to mix secondary voltage control solving method with continuous quantity, in the theoretical foundation of traditional secondary voltage control, aim to solve the problem that discrete magnitude smoothing mixes secondary voltage control method with continuous quantity, realize that optimal coordination calculates to discrete in region and continuous reactive source, this patent provides a kind of discrete magnitude coherency function smoothing and secondary voltage control method for solving is mixed with continuous quantity, extension secondary voltage control optimization is improved to coordinate, ensure the control of discrete in the region and optimized running of reactive source, effectively increase the economical of reactive power/voltage control in power network, accuracy and security.
Description
Technical field
The present invention relates to a kind of discrete magnitude coherency function smoothing to mix secondary voltage control solving method with continuous quantity, belongs to
Operation and control of electric power system technical field.
Background technology
Secondary voltage control is the Core Feature that automatism voltage control realizes that reactive Voltage Optimum distribution calculates in region, is
The coordination control of discrete magnitude and continuous quantity in region is realized, traditional way is solved by two-phase method, the first rank
Section carries out voltage correction in the out-of-limit situation of discrete magnitude direct correlation busbar voltage, and otherwise discrete magnitude is without regulation, second stage
Idle adjustment amount conversion correction is added into secondary voltage control object function and is iterated solution, calculates the control mesh of continuous quantity
Mark.This control mode has preferably decomposed discrete magnitude and continuous quantity coordination problem in region, has good control convergence
Property, the voltage-controlled demand of electric network reactive-load is disclosure satisfy that substantially.Applied however as extensive discrete magnitude in region and city is electric
The appearance of the electric network composition situations such as net reactive power support source deficiency, the region idle work optimization method of this more " rough " can not meet
Actual reactive power/voltage control demand, to electric network reactive-load voltage security control can exist not in time, irrational phenomenon occur, to electricity
Net safe operation causes to have a strong impact on.
The content of the invention
To solve the deficiencies in the prior art, it is an object of the invention to provide a kind of discrete magnitude coherency function smoothing and company
Continuous amount mixing secondary voltage control solving method, the improvement to traditional secondary voltage control, discrete magnitude and continuous quantity in region
Reactive voltage coordinates the amendment of mixed integer programming problem.
In order to realize above-mentioned target, the present invention adopts the following technical scheme that：
A kind of discrete magnitude coherency function smoothing mixes secondary voltage control solving method with continuous quantity, it is characterized in that, including
Following steps：
1) method of operation of discrete reactive apparatus is passed through into serialization function representation；
2) NCP functions being substituted using smoothing functions, introducing coherency function is smooth to the progress of Constraints condition to be approached,
Constraints condition is converted into the smooth Nonlinear System of Equations of series；
3) Nonlinear System of Equations for obtaining step 2) is brought into secondary voltage control model and is modified；
4) region reactive voltage target is calculated into optimal solution.
Further, particular content is in the step 1)：By the method for operation of discrete reactive apparatus be denoted as x (x-1)=
0, wherein,
Further, the discrete magnitude of the reactive apparatus discrete described in step 2) is when carrying out smoothing, according to mutual
Theoretical and coherency function is mended, Constraints condition is converted into the smooth Nonlinear System of Equations of series.
Further, resulting Nonlinear System of Equations is φ (x)=- μ ln (e-(x-1)/μ+e(x-1)/μ), wherein μ is smooth
Parameter.
Further, the secondary voltage control model in the step 3) is expressed as：
Object function V in formulapWith
Vp refMaincenter bus voltage value and voltage optimization desired value respectively in region；ΔQgFor continuous reactive source idle adjustment amount to
Amount, Δ QcFor the idle adjustment amount vector of discrete reactive source；For Δ QgTransposed vector,For Δ QcTransposed vector；
H is the weight coefficient in continuous reactive source, and r is the weight coefficient of discrete reactive source；CpgIt is continuous reactive source to maincenter busbar voltage
Sensitivity, CpcIt is idle to maincenter busbar voltage sensitivity for discrete reactive source；
Constrain Q in inequalitygFor continuous reactive source currently idle output, QcFor the current idle output of discrete reactive source；V pWithRespectively maincenter busbar voltage bound constraint,Q gWithRespectively the idle output in continuous reactive source bound constraint,Q cWithThe bound constraint of the respectively discrete idle output of reactive source；C'pcIt is idle to maincenter bus electricity for discrete magnitude reactive source
Press the famous value of sensitivity.
Further, the correction formula obtained in the step 3) is
The beneficial effect that the present invention is reached：This patent is in traditional secondary voltage control theoretical foundation, by condensing letter
It is several that continuous treatment is carried out to discrete magnitude, complex mixed integer programming problem is converted into nonlinear programming problem, filled
Divide and consider discrete magnitude, disposably solve discrete and continuous quantity coordination optimization object function, while discrete and continuous quantity is generated
Control targe, this processing mode disposably solve discrete magnitude and continuous quantity Harmonic Control, and its calculate time also compared with
To be rapid, the degree of convergence is also higher, is effectively improved the precision of region reactive power/voltage control, reduces idle network loss in region,
Ensure the safety and economic operation of electric network reactive-load voltage.
Embodiment
The invention will be further described below.Following examples are only used for the technical side for clearly illustrating the present invention
Case, and can not be limited the scope of the invention with this.
The present invention is the improvement to traditional secondary voltage control, and discrete magnitude and continuous quantity reactive voltage are coordinated in region
The amendment of mixed integer programming problem, the foundation of invention are complementary theory and coherency function method by discrete magnitude smoothing and region
Secondary voltage control function, is described as follows：
Because discrete reactive apparatus are discrete switch amount, can be designated as：
Its method of operation meets x (x-1)=0, and it is secondary that discrete variable characteristic can thus be changed into continuous quantity addition
In planning, the problem is set to be converted into continuous optimization problems, but for quadratically constrained quadratic programming problem, it is unsatisfactory for optimality
Condition, and this is also a challenge for being very difficult to solve in itself.
Discrete magnitude linear complementary problem is handled by NCP functions, but NCP functions are non-smooth functions, in order to solve
The problem, this patent are handled using smoothing, by introducing smoothing functions alternative functions, are further introduced into coherency function pair
The progress of Constraints condition is smooth to approach, as follows so as to which Constraints condition is converted into the smooth Nonlinear System of Equations of series
Represent：φ (x)=- μ ln (e-(x-1)/μ+e(x-1)/μ)。
Secondary voltage control is extended using the coherency function, traditional secondary voltage control can be corrected.
Conventional discrete amount and continuous quantity extension secondary voltage control model are expressed as below：
Object function V in formulapWith
Vp refMaincenter bus voltage value and voltage optimization desired value respectively in region；ΔQgWith Δ QcRespectively continuous reactive source and from
Dissipate the idle adjustment amount of reactive source；CpgAnd CpcRespectively continuous reactive source and discrete reactive source are idle to the spirit of maincenter busbar voltage
Sensitivity；Constrain Q in inequalityg、QcThe current idle output of respectively continuous and discrete reactive source；V pWith Q gWith Q cWith
Respectively maincenter busbar voltage, continuous reactive source is idle contributes, the constraint of discrete reactive source is idle output；C'pcFor discrete magnitude without
Work(source is idle to the famous value of maincenter busbar voltage sensitivity.
Discrete magnitude is smoothed willSecondary voltage control model is substituted into, then formula is modified to：
So secondary voltage control model is converted into Nonlinear programming Model from mixed-integer programming model, reduces and solves hardly possible
Degree, realize discrete, continuous control variable coordination control.
It can be seen that region secondary voltage control target is modified to nonlinear programming problem by this method, effectively to discrete magnitude
Coordination calculating is carried out with continuous quantity, region reactive voltage target is calculated into optimal solution.
Described above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, some improvement and deformation can also be made, these are improved and deformation
Also it should be regarded as protection scope of the present invention.
Claims (6)
1. a kind of discrete magnitude coherency function smoothing mixes secondary voltage control solving method with continuous quantity, it is characterized in that, including such as
Lower step：
1) method of operation of discrete reactive apparatus is passed through into serialization function representation；
2) NCP functions are substituted using smoothing functions, introducing coherency function is smooth to the progress of Constraints condition to be approached, will be mutual
Mend constraints and be converted into the smooth Nonlinear System of Equations of series；
3) Nonlinear System of Equations for obtaining step 2) is brought into secondary voltage control model and is modified；
4) region reactive voltage target is calculated into optimal solution.
2. a kind of discrete magnitude coherency function smoothing according to claim 1 mixes secondary voltage control solution with continuous quantity
Method, it is characterized in that, particular content is in the step 1)：The method of operation of discrete reactive apparatus is denoted as x (x-1)=0, its
In,
3. a kind of discrete magnitude coherency function smoothing according to claim 2 mixes secondary voltage control solution with continuous quantity
Method, it is characterized in that, the discrete magnitude of the reactive apparatus discrete described in step 2) is when carrying out smoothing, according to complementary theory
And coherency function, Constraints condition is converted into the smooth Nonlinear System of Equations of series.
4. a kind of discrete magnitude coherency function smoothing according to claim 3 mixes secondary voltage control solution with continuous quantity
Method, it is characterized in that, resulting Nonlinear System of Equations is φ (x)=- μ ln (e-(x-1)/μ+e(x-1)/μ), wherein μ is smoothing parameter.
5. a kind of discrete magnitude coherency function smoothing according to claim 1 mixes secondary voltage control solution with continuous quantity
Method, it is characterized in that, the secondary voltage control model in the step 3) is expressed as：
Object function V in formulapAnd Vp refPoint
Wei not maincenter bus voltage value and voltage optimization desired value in region；ΔQgFor the idle adjustment amount vector in continuous reactive source, Δ Qc
For the idle adjustment amount vector of discrete reactive source；For Δ QgTransposed vector,For Δ QcTransposed vector；H is continuous
The weight coefficient of reactive source, r are the weight coefficient of discrete reactive source；CpgIt is continuous reactive source to maincenter busbar voltage sensitivity,
CpcIt is idle to maincenter busbar voltage sensitivity for discrete reactive source；
Constrain Q in inequalitygFor continuous reactive source currently idle output, QcFor the current idle output of discrete reactive source；V pWithPoint
Not Wei maincenter busbar voltage bound constraint,Q gWithRespectively the idle output in continuous reactive source bound constraint,Q cWithThe bound constraint of the respectively discrete idle output of reactive source；C'pcIt is idle to the spirit of maincenter busbar voltage for discrete magnitude reactive source
The famous value of sensitivity.
6. a kind of discrete magnitude coherency function smoothing according to claim 5 mixes secondary voltage control solution with continuous quantity
Method, it is characterized in that, the correction formula obtained in the step 3) is
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mi>min</mi>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>&Delta;</mi>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mo>|</mo>
<mo>|</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>g</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;Q</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mi>&Delta;</mi>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<msubsup>
<mi>&Delta;Q</mi>
<mi>g</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>h&Delta;Q</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>&phi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>&gamma;</mi>
<mi>&Delta;</mi>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<munder>
<mi>V</mi>
<mo>&OverBar;</mo>
</munder>
<mi>p</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>g</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;Q</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mi>&Delta;</mi>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mi>p</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<munder>
<mi>Q</mi>
<mo>&OverBar;</mo>
</munder>
<mi>g</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;Q</mi>
<mi>g</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>g</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<mi>&Delta;</mi>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710846406.5A CN107732923A (en) | 2017-09-19 | 2017-09-19 | The smoothing of discrete magnitude coherency function mixes secondary voltage control solving method with continuous quantity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710846406.5A CN107732923A (en) | 2017-09-19 | 2017-09-19 | The smoothing of discrete magnitude coherency function mixes secondary voltage control solving method with continuous quantity |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107732923A true CN107732923A (en) | 2018-02-23 |
Family
ID=61207497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710846406.5A Pending CN107732923A (en) | 2017-09-19 | 2017-09-19 | The smoothing of discrete magnitude coherency function mixes secondary voltage control solving method with continuous quantity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107732923A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102664417A (en) * | 2012-04-26 | 2012-09-12 | 广东省电力调度中心 | Control method and device for secondary voltage |
| CN105633974A (en) * | 2016-03-03 | 2016-06-01 | 甘肃省电力公司风电技术中心 | Real-time regional voltage coordination control method |
| US20160204609A1 (en) * | 2015-01-12 | 2016-07-14 | Dominion Resources, Inc. | Systems and methods for volt-ampere reactive control and optimization |
-
2017
- 2017-09-19 CN CN201710846406.5A patent/CN107732923A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102664417A (en) * | 2012-04-26 | 2012-09-12 | 广东省电力调度中心 | Control method and device for secondary voltage |
| US20160204609A1 (en) * | 2015-01-12 | 2016-07-14 | Dominion Resources, Inc. | Systems and methods for volt-ampere reactive control and optimization |
| CN105633974A (en) * | 2016-03-03 | 2016-06-01 | 甘肃省电力公司风电技术中心 | Real-time regional voltage coordination control method |
Non-Patent Citations (1)
| Title |
|---|
| 韦园清: "基于Sigmoid和NCP函数的电力系统无功优化连续化方法", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9728974B2 (en) | Renewable energy site reactive power control | |
| CN103368173A (en) | Active power flow optimizing distribution method for alternating-direct current parallel system containing soft direct current power transmission | |
| CN104578857A (en) | Control method and device of photovoltaic power generation system and photovoltaic power generation system | |
| CN106817041A (en) | Power inverter, electricity generation system, control device and power transferring method | |
| Vega et al. | Scalar hadrons in AdS 5× S 5 | |
| CN105867126A (en) | Fractional order PI optimization control method of three-phase voltage source type inverter system | |
| CN107732923A (en) | The smoothing of discrete magnitude coherency function mixes secondary voltage control solving method with continuous quantity | |
| CN110962666A (en) | Electric vehicle charging equipment based on load virtual synchronous machine technology and control method | |
| CN103683329A (en) | Method for optimizing grid-connected unit excitation difference adjustment coefficient based on minimum overall network loss | |
| CN111864766A (en) | Control method for additional damping of direct current power transmission system based on angular frequency change rate | |
| CN111291985A (en) | Day-ahead power market clearing method and device based on conditional section modeling | |
| CN106992522A (en) | Voltage source converter power flow algorithm and active power distribution network tidal current computing method | |
| CN106099935B (en) | The droop control method and droop control device that busbar voltage deviation compensates automatically | |
| CN107069748A (en) | Using the dynamic electric voltage recovery device compensation control system and method for minimum current injection method | |
| CN111293710B (en) | Bidirectional charge and discharge control system of electric automobile and control method thereof | |
| CN106505582A (en) | A kind of dynamic reactive power voltage cooperative control method based on neural network forecast mode | |
| CN112803429B (en) | Reactive voltage coordination control method and system for energy storage power station and clean energy power station | |
| Ahmed et al. | Modeling and simulation of power electronic distribution transformer based on a three level converter | |
| CN108429266B (en) | Continuous and discrete reactive power source coordination control method based on dynamic coordination boundary threshold | |
| CN201174679Y (en) | Energy saving controller for loaded pressure regulating switch | |
| CN108199387A (en) | A kind of active distribution network reactive voltage progress control method based on expansion cone planning | |
| Monteiro et al. | High Power Factor Rectifier Using One Cycle Control Strategy | |
| CN111884227B (en) | UPFC model parameter adjustment method and system | |
| CN109626853A (en) | A kind of harmless treatment process of electrolytic manganese residues | |
| Landman | Proton collisional excitation in the ground/3/P terms of oxygen and sulfur isoelectronic ions |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180223 |