CN110492493A - A kind of power system reactive power compensation method for optimizing configuration - Google Patents
A kind of power system reactive power compensation method for optimizing configuration Download PDFInfo
- Publication number
- CN110492493A CN110492493A CN201910756767.XA CN201910756767A CN110492493A CN 110492493 A CN110492493 A CN 110492493A CN 201910756767 A CN201910756767 A CN 201910756767A CN 110492493 A CN110492493 A CN 110492493A
- Authority
- CN
- China
- Prior art keywords
- node
- reactive
- compensation
- alternate
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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/18—Arrangements for adjusting, eliminating or compensating reactive power in 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
- H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
- H02J3/1871—Methods for planning installation of shunt reactive power compensators
-
- 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 present invention proposes a kind of power system reactive power compensation method for optimizing configuration, belongs to power system reactive power compensation field of configuration.This method initially sets up the electromechanical transient simulation model of electric system to be optimized, and obtains the flow data and short circuit current of each node in electric system to be optimized；After node-classification, node load rate is considered, obtain the initial alternate node of configuration reactive compensation, consider that the effect in reactive voltage source and more feed-in interaction factors are polymerize, finally obtain the reactive compensation alternate node of electric system to be optimized；Configuration reactive capability needed for each reactive compensation alternate node is optimized using genetic algorithm, obtains the reactive capability Optimization Compensation result of each reactive compensation alternate node.The present invention is easily achieved and accuracy is high, is a kind of method of practical power system reactive power compensation configuration optimization.
Description
Technical field
The present invention relates to a kind of power system reactive power compensation method for optimizing configuration, belong to power system reactive power compensation configuration neck
Domain.
Background technique
Since resource allocation is uneven, region disparate development, straight through the trans-regional friendship of interior development after a period of time
After stream transmission of electricity, more multi-feed HVDC loading zone is formd.In AC/DC Power System, topological structure constantly changes
Alternating current-direct current inverter, either rectification or inverter operation, become AC system because it need to consume a large amount of reactive power
Nonlinear load.I.e. for AC system, DC transmission system is a load or burden without work.And not by AC system voltage
It can restore or persistently reduce may cause that multiple-circuit line transmission system is caused to recur the even locking of multiple commutation failure.Its
Geological Problems are that System Reactive Power support is mismatched with reactive requirement.Therefore, being equipped with a certain amount of reactive compensation in the power system is
It is very necessary.
There is research to optimize using modified particle swarm optiziation to alternating current-direct current electric network reactive compensation at present, in mentioned method
Reconnaissance is carried out to avoid overcompensation using multiple Load flow calculation, but increases certain calculation amount；In addition, in terms of capacity configuration
It may not be able to realize the effect of global compensation optimization by the way of being grouped optimizing.
However, AC/DC Power System interior joint is numerous, in order to reasonably be equipped with reactive power compensator, it is therefore desirable to a kind of
The method of effective reactive compensation configuration optimization, node selection and capacity optimization including idle configuration.
Summary of the invention
The purpose of the present invention is the shortcoming to overcome prior art, propose that a kind of power system reactive power compensation configuration is excellent
Change method.On the one hand this method polymerization passes through load factor primary dcreening operation alternate node after, using more feed-in interaction factors (MIIF)
The strong node of reciprocation is polymerize in alternate node, while being considered reactive voltage source-generator and being combined with MIIF
Effect, choose in polymerization with the strongest node of other node reciprocations alternately node, to system node network into
Row is shunk, and improves the efficiency of alternate node screening；On the other hand overcome optimization of the tradition to the idle configuration capacity of single node
Method is uniformly processed the optimization of global alternate node using genetic algorithm, and carries out multiple Load flow calculation, to obtain more
For the offset of optimization, to solve the problems, such as that AC/DC Power System wattles power economic equivalent (Voltage-stabilizing Problems) provides a kind of row
Effective method.
The present invention proposes a kind of power system reactive power compensation method for optimizing configuration, which comprises the following steps:
1) it is based on electromechanical transient simulation program, the electromechanical transient for establishing electric system to be optimized using lumped parameter model is imitative
True mode；
2) Load flow calculation and short circuit current meter are carried out to the simulation model that step 1) is established by electromechanical transient program respectively
It calculates, obtains the flow data of each node in electric system to be optimized, include: each node voltage real part ei, each node voltage
Imaginary part fi, each node power and short circuit current；
3) the reactive compensation alternate node of electric system to be optimized is chosen；Specific step is as follows:
The node in electric system to be optimized 3-1) is divided into load bus, generator node and inverter node, wherein
Load bus and generator node form AC system node, and inverter node is direct current system node；
3-2) calculate separately each node of load factor and direct current system of each load bus in AC system node
Load factor；
Wherein in AC system node the load factor of load bus calculation expression are as follows:
Wherein, SiFor the apparent energy of node i, QliFor node i reactive load power, l indicates load；ΔQdp,i|ΔUiTable
Show that disturbance posterior nodal point i voltage occurs at node i occurs Δ UiVariation when correspond to the variable quantity of reactive power；
The calculation expression of DC node load factor is as follows:
3-3) according to the calculated result of step 3-2), each AC system load bus and DC node are determined:
If the load factor of any node is higher than system Rate of average load and is somebody's turn to do in AC system load bus or DC node
Node load rate ranking the similar node load rate of system preceding one third, then using the node as configuration reactive compensation just
Beginning alternate node；
3-4) all initial alternate nodes that step 3-3) is obtained are screened；Specific step is as follows:
3-4-1) calculate more feed-in interaction factors of other all nodes in each initial alternate node and system
MIIFij；
The wherein MIIF between node i and node jijCalculation expression are as follows:
Wherein, ZijIndicate the mutual impedance between node i and node j, ZiiIndicate the self-impedance of node i；
3-4-2) using the calculated result of step 3-4-1), by all MIIFij> 0.3 initial alternate node is polymerize；
It is chosen in the node set of each polymerization and is all larger than 0.3 with other all more feed-in interaction factor values of node in the set
Node screened after initial alternate node；
3-5) when node i is generator node, the more of other all nodes in each generator node and system are calculated
Feed-in interaction factor MIIFij；By MIIF in calculated resultij> 0.3 corresponding node j is polymerize, if what polymerization obtained
As a result there are nodes to belong to step 3-4 in) initial alternate node after screening, then by the node from the initial alternative section after screening
It is deleted in point, finally obtains the number of reactive compensation alternate node number m and each alternate node；
4) configuration reactive capability needed for each reactive compensation alternate node obtained using genetic algorithm to step 3) is carried out
Optimization, obtains the reactive capability Optimization Compensation result of each reactive compensation alternate node；Specific step is as follows:
4-1) obtain parameters of electric power system and variable to be optimized, comprising: the number of all nodes, the voltage of each node are real
Portion eiWith imaginary part fi, connection relationship between node power and node；
4-2) by configuration reactive capability coding needed for each reactive compensation alternate node of electric system to be optimized and to system
Trend is calculated, according to reactive compensation alternate node reactive power value Q each in the result acquisition system of step 4-1)i[k]As this
The initial reactive power value of node, and assignment qi[k]=Qi[k]；Enabling each alternate node compensation rate capacity is QCi[k]And assign initial value QCi[k]
=0, update the reactive power Q of alternate nodei[k]=Qi[k]+QCi[k], wherein Qi[k]Indicate k-th of reactive compensation alternate node
Reactive power value, QCi[k]Indicate that the reactive compensation capacity of k-th of reactive compensation alternate node, k indicate alternate node number, k=
1,2…m；
Initial population kind 4-3) is randomly generated；
Calculating 4-4) is iterated to system load flow, the system interior joint voltage and node power after updating reactive compensation；
4-5) determine system optimization target, with the minimum objective function of voltage deviation:
F=min (fΔV)
Wherein,
Wherein, Δ V indicates voltage variety, V0For voltage reference value, ViFor the voltage value of node i, VimaxWith ViminRespectively
Indicate the voltage magnitude bound of node i, i indicates that node serial number, i=1,2 ... N, N are system node sum；
Constraint condition includes:
Equality constraint:
Wherein, Pi SP、Qi SPThe respectively active power of PQ node i, reactive power；ΔPi、ΔQiRespectively node i is active
The departure of power and reactive power；GijAnd BijIt is the real and imaginary parts of admittance matrix the i-th row jth column element respectively；θijFor section
Phase difference of voltage between point i and j；QCiFor capacity of reactive power compensation device at node i；Pdi、QdiRespectively alternating current-direct current connecting node
Active power, reactive power at i；
Inequality constraints:
Vimin≤Vi≤Vimax,i∈N
QCimin≤QCi≤QCimax,i∈NC
Wherein, NCNumber is compensated for Reactive Power Device；VimaxWith ViminRespectively indicate the voltage magnitude upper limit and lower limit of node i；
QCimax、QCiminRespectively indicate the upper and lower bound of compensation capacity；
The objective function of step 4-5) 4-6) is converted into fitness function, and successively selected according to genetic algorithm,
Intersect, mutation process, obtains the fitness value of each reactive compensation alternate node and the compensation capacity correction value Δ of the node
qCi[k], each node voltage real part correction amount eiWith imaginary part correction amount fi；The fitness function is objective function
It is reciprocal；
4-7) real and imaginary parts of all node voltages are modified respectively:
ei=ei+Δei, fi=fi+Δfi
4-8) judge whether the fitness value of each reactive compensation alternate node is greater than the inverse of target function value:
If it is, calculating the reactive compensation amount Δ Q of each reactive compensation alternate nodeCi[k]=Qi[k]-qi[k], export Δ
QCi[k]As the Optimization Compensation for configuring reactive capability needed for each reactive compensation alternate node as a result, optimization terminates；
If it is not, then according to the compensation capacity correction value Δ q of each reactive compensation alternate nodeCi[k], update QCi[k]=
QCi[k]+ΔqCi[k], then return to step 4-4), it is adapted to until the fitness value of all reactive compensation alternate nodes is greater than
Angle value index calculates the reactive compensation amount Δ Q of each reactive compensation alternate nodeCi[k]=Qi[k]-qi[k], export Δ QCi[k]As every
The Optimization Compensation of configuration reactive capability needed for a reactive compensation alternate node is as a result, optimization terminates.
The features of the present invention and beneficial effect are:
This method by node load rate in AC system node load bus and direct current system node carry out it is preliminary
Screening, fully considers its idle dynamic change for AC system node；It will using more feed-in interaction factors (MIIF)
The big node of mutual reciprocation is polymerize in system, i.e., if there is between some node and other several nodes
MIIF it is bigger, then can choose the point alternately node, not need to compensate corresponding each node；In addition, by nothing
The potential spare value of function voltage source is combined with more feed-in interaction factors, will be existed between reactive voltage source node strong
Interactive node is polymerize, and such aggregation is excluded within the scope of alternate node.By considering node load
Rate and it is polymerize by MIIF, on the one hand reduces the range and configuration quantity of idle configuration reconnaissance in system, on the other hand
It can play the role of effective pinch system meshed network, the efficiency of lifting section point selection；Overcome tradition idle to single node
The optimization method of configuration capacity is uniformly processed the optimization of global alternate node using genetic algorithm, and carries out repeatedly tide
Stream calculation, to obtain the offset more optimized.This method is to solve the problems, such as AC/DC Power System wattles power economic equivalent (voltage
Stable problem) provide a kind of effective method.
Load flow calculation that this method refers to, calculation of short-circuit current, genetic algorithm have been maturation method, it is easy to accomplish and it is quasi-
Exactness is high, is a kind of method of practical power system reactive power compensation configuration optimization.
Detailed description of the invention
Fig. 1 is the overview flow chart of the method for the present invention.
Specific embodiment
A kind of power system reactive power compensation method for optimizing configuration proposed by the present invention, in the following with reference to the drawings and specific embodiments
It is as follows that the present invention is described in more detail.
A kind of power system reactive power compensation method for optimizing configuration proposed by the present invention, overall flow as shown in Figure 1, include with
Lower step:
1) it is based on electromechanical transient simulation program (such as BPA), electric system to be optimized is established using lumped parameter model
Electromechanical transient simulation model；
2) Load flow calculation and short circuit current meter are carried out to the simulation model that step 1) is established by electromechanical transient program respectively
It calculates, obtains the flow data of each node in electric system to be optimized, include: each node voltage real part ei, each node voltage
Imaginary part fi, each node power and short circuit current；
3) the reactive compensation alternate node of electric system to be optimized is chosen；Specific step is as follows:
3-1) node in electric system can be divided into load bus, generator node and inverter node, by power train
System node is classified according to ac and dc systems, is divided into AC system node (load bus and generator node) and direct current system
Node (inverter node)；
3-2) calculate separately each node of load factor and direct current system of each load bus in AC system node
Load factor；
Wherein in AC system node the load factor of load bus calculation expression are as follows:
Wherein, SiFor the apparent energy of node i, QliFor node i reactive load power (l indicates load)；ΔQdp,i|ΔUiTable
Show that disturbance posterior nodal point i voltage occurs at node i occurs Δ UiVariation when correspond to the variable quantity of reactive power；
The calculation expression of DC node load factor is as follows:
3-3) according to the calculated result of step 3-2), each AC system load bus and DC node are sentenced respectively
It is fixed:
The load factor of load bus or direct current system node is higher (if the node load rate is flat higher than system in AC system
The preceding one third of equal load factor and load factor ranking in the similar node load rate of system), the probability that the node demand is idle is got over
(probability that voltage stability problem occurs is bigger) greatly, then the node is considered as the weak link of electric system, may be selected as matching
Set the initial alternate node of reactive compensation；
3-4) the initial alternate node obtained to step 3-3) utilizes more feed-in interaction factor (MIIFij) calculate just
Reciprocation size in beginning alternate node in each node and system between other all nodes, in this, as AC and DC transmission
Reactive compensation is layouted selection factor in system, and the strong node of more feed-in reciprocations is polymerize by consideration (works as MIIFij>0.3
When, then it is assumed that the reciprocation between node i and node j is strong).It is chosen in the node set of each polymerization and other in set
There is strong interactive node as the reactive compensation alternate node after screening and (work as MIIF in nodeijWhen > 0.3, then it is assumed that section
Reciprocation between point i and node j is strong), it is shunk with treating the meshed network of optimization system, while being conducive to be lifted at
Computational efficiency in reactive compensation capacity optimization process.
The wherein MIIF between node i and node jijCalculation expression are as follows:
Wherein, ZijIndicate the mutual impedance between node i and node j, ZiiIndicate the self-impedance of node i.Work as MIIFij>0.3
When, then it is assumed that the reciprocation between node i and node j is strong.
3-5) synchronous generator is important reactive power source in system, and therefore, reactive compensation matches and purchases in consideration system
When selecting node, more feed-in interaction factors of generator node are also calculated；It, will be with hair when node i is generator node
There are strong reciprocation (MIIF for motor nodeij> 0.3) node j is polymerize, and the result and step 3-4 that polymerization is obtained)
Result lap after screening is deleted, and the final reactive compensation alternate node number for needing to configure reactive power is obtained
The number of m and each alternate node；
4) configuration reactive capability needed for each reactive compensation alternate node obtained using genetic algorithm to step 3) is carried out
Optimization, obtains the reactive capability Optimization Compensation result of each reactive compensation alternate node；Specific step is as follows:
4-1) obtain parameters of electric power system and variable to be optimized, comprising: the voltage of the number of all nodes, each node
(including real part ei, imaginary part fi), the connection relationship between node power (including reactive power and active power) and node；
4-2) each reactive compensation alternate node reactive capability of electric system to be optimized offset that may be present is carried out
(genetic algorithm needs for the solution to study a question to be encoded into the form of character string to coding, the corresponding system to be optimized of " problem " herein
The reactive compensation capacity that may be present of each alternate node) and to system-computed trend, system is obtained according to the result of step 4-1)
Each reactive compensation alternate node reactive power value Q in systemi[k]As the initial reactive power value of the node, and assignment qi[k]=
Qi[k]；Enabling each alternate node compensation rate capacity is QCi[k]And assign initial value QCi[k]=0, the reactive power Q of alternate nodei[k]=Qi[k]
+QCi[k], wherein Qi[k]Indicate the reactive power value of k-th of reactive compensation alternate node, QCi[k]Indicate that k-th of reactive compensation is alternative
The reactive compensation capacity of node, k indicate alternate node number, k=1,2 ... m；
Initial population kind 4-3) is randomly generated (size and number of group's kind can be configured according to the scale of system)；
Iterative calculation 4-4) is re-started to system load flow, the system interior joint voltage and node function after updating reactive compensation
Rate；
4-5) determine system optimization target: in the case where meeting various constraint conditions, with the smallest idle investment, to greatest extent
Voltage stabilization degree is improved, improves quality of voltage, reduce grid loss in ground.This sentences the minimum target letter of voltage deviation
Number:
F=min (fΔV)
Wherein,
Wherein, Δ V indicates that voltage variety, N are electric system node total number, V0For voltage reference value, ViFor node i
Voltage value, VimaxWith ViminThe voltage magnitude bound of node i is respectively indicated, i indicates node serial number, and i=1,2 ... N, N are system
Node total number.
The constraint condition are as follows:
Equality constraint (node power constraint):
Wherein, Pi SP、Qi SPThe respectively active power of PQ node i, reactive power；ΔPi、ΔQiRespectively node i is active
The departure of power and reactive power；GijAnd BijIt is the real and imaginary parts of admittance matrix the i-th row jth column element respectively；θijFor section
Phase difference of voltage between point i and j；QCiFor capacity of reactive power compensation device at node i；Pdi、QdiRespectively alternating current-direct current connecting node
Active power, reactive power at i；
Inequality constraints:
Vimin≤Vi≤Vimax,i∈N
QCimin≤QCi≤QCimax,i∈NC
Wherein, NCNumber is compensated for Reactive Power Device；VimaxWith ViminRespectively indicate the voltage magnitude upper limit and lower limit of node i
(0.95~1.05)；QCimax、QCiminRespectively indicate the upper and lower bound of compensation capacity.
4-6) objective function of step 4-5) is converted into fitness function (inverse of objective function), and is calculated according to heredity
Method successively selected, intersected, mutation process, obtains the fitness value of each reactive compensation alternate node and the benefit of the node
Repay cubage correction value Δ qCi[k], each node voltage real part correction amount eiWith imaginary part correction amount fi；
4-7) real and imaginary parts of all node voltages are modified respectively: ei=ei+Δei, fi=fi+Δfi, eiWith
fiThe respectively real and imaginary parts of the i-th row of node voltage matrix；
4-8) judge whether the fitness value of each reactive compensation alternate node is greater than set fitness value index (mesh
The inverse of offer of tender numerical value):
If it is, the reactive compensation amount Δ Q of each reactive compensation alternate node can be calculatedCi[k]=Qi[k]-qi[k], output
ΔQCi[k]As the Optimization Compensation for configuring reactive capability needed for each reactive compensation alternate node as a result, optimization terminates；If
It is no, update QCi[k]=QCi[k]+ΔqCi[k], then return to step 4-4), until all reactive compensation alternate nodes is suitable
It answers angle value to be greater than fitness value index, calculates the reactive compensation amount Δ Q of each reactive compensation alternate nodeCi[k]=Qi[k]-qi[k], defeated
Δ Q outCi[k]As the Optimization Compensation for configuring reactive capability needed for each reactive compensation alternate node as a result, optimization terminates.
Claims (1)
1. a kind of power system reactive power compensation method for optimizing configuration, which comprises the following steps:
1) it is based on electromechanical transient simulation program, the electromechanical transient simulation mould of electric system to be optimized is established using lumped parameter model
Type；
2) Load flow calculation and calculation of short-circuit current are carried out to the simulation model that step 1) is established by electromechanical transient program respectively, obtained
The flow data for taking each node in electric system to be optimized includes: each node voltage real part ei, each node voltage imaginary part fi、
Each node power and short circuit current；
3) the reactive compensation alternate node of electric system to be optimized is chosen；Specific step is as follows:
The node in electric system to be optimized 3-1) is divided into load bus, generator node and inverter node, wherein load
Node and generator node form AC system node, and inverter node is direct current system node；
3-2) calculate separately the load of each node of load factor and direct current system of each load bus in AC system node
Rate；
Wherein in AC system node the load factor of load bus calculation expression are as follows:
Wherein, SiFor the apparent energy of node i, QliFor node i reactive load power, l indicates load；Indicate node i
Place occurs disturbance posterior nodal point i voltage and Δ U occursiVariation when correspond to the variable quantity of reactive power；
The calculation expression of DC node load factor is as follows:
3-3) according to the calculated result of step 3-2), each AC system load bus and DC node are determined:
If the load factor of any node is higher than system Rate of average load and the node in AC system load bus or DC node
Load factor ranking the similar node load rate of system preceding one third, then using the node as configuration reactive compensation it is initial standby
Select node；
3-4) all initial alternate nodes that step 3-3) is obtained are screened；Specific step is as follows:
3-4-1) calculate more feed-in interaction factor MIIF of other all nodes in each initial alternate node and systemij；
The wherein MIIF between node i and node jijCalculation expression are as follows:
Wherein, ZijIndicate the mutual impedance between node i and node j, ZiiIndicate the self-impedance of node i；
3-4-2) using the calculated result of step 3-4-1), by all MIIFij> 0.3 initial alternate node is polymerize；Every
Selection is all larger than 0.3 section with other all more feed-in interaction factor values of node in the set in the node set of a polymerization
Put the initial alternate node after being screened；
3-5) when node i is generator node, more feed-ins of other all nodes in each generator node and system are calculated
Interaction factor MIIFij；By MIIF in calculated resultij> 0.3 corresponding node j is polymerize, if the result that polymerization obtains
It is middle that there are nodes to belong to step 3-4) screening after initial alternate node, then by the node from the initial alternate node after screening
It deletes, finally obtains the number of reactive compensation alternate node number m and each alternate node；
4) configuration reactive capability needed for each reactive compensation alternate node obtained using genetic algorithm to step 3) is optimized,
Obtain the reactive capability Optimization Compensation result of each reactive compensation alternate node；Specific step is as follows:
4-1) obtain parameters of electric power system and variable to be optimized, comprising: the voltage real part e of the number of all nodes, each nodei
With imaginary part fi, connection relationship between node power and node；
4-2) by configuration reactive capability coding needed for each reactive compensation alternate node of electric system to be optimized and to system-computed
Trend, according to reactive compensation alternate node reactive power value Q each in the result acquisition system of step 4-1)i[k]As the node
Initial reactive power value, and assignment qi[k]=Qi[k]；Enabling each alternate node compensation rate capacity is QCi[k]And assign initial value QCi[k]=0,
Update the reactive power Q of alternate nodei[k]=Qi[k]+QCi[k], wherein Qi[k]Indicate the idle of k-th of reactive compensation alternate node
Performance number, QCi[k]Indicate that the reactive compensation capacity of k-th of reactive compensation alternate node, k indicate that alternate node is numbered, k=1,2 ...
m；
Initial population kind 4-3) is randomly generated；
Calculating 4-4) is iterated to system load flow, the system interior joint voltage and node power after updating reactive compensation；
4-5) determine system optimization target, with the minimum objective function of voltage deviation:
Wherein,
Wherein, Δ V indicates voltage variety, V0For voltage reference value, ViFor the voltage value of node i, VimaxWith ViminIt respectively indicates
The voltage magnitude bound of node i, i indicate that node serial number, i=1,2 ... N, N are system node sum；
Constraint condition includes:
Equality constraint:
Wherein, Pi SP、The respectively active power of PQ node i, reactive power；ΔPi、ΔQiRespectively node i active power
With the departure of reactive power；GijAnd BijIt is the real and imaginary parts of admittance matrix the i-th row jth column element respectively；θijFor node i
Phase difference of voltage between j；QCiFor capacity of reactive power compensation device at node i；Pdi、QdiAt respectively alternating current-direct current connecting node i
Active power, reactive power；
Inequality constraints:
Vimin≤Vi≤Vimax,i∈N
QCimin≤QCi≤QCimax,i∈NC
Wherein, NCNumber is compensated for Reactive Power Device；VimaxWith ViminRespectively indicate the voltage magnitude upper limit and lower limit of node i；QCimax、
QCiminRespectively indicate the upper and lower bound of compensation capacity；
The objective function of step 4-5) 4-6) is converted into fitness function, and successively selected, intersected according to genetic algorithm,
Mutation process obtains the fitness value of each reactive compensation alternate node and the compensation capacity correction value Δ q of the nodeCi[k]、
The voltage real part correction amount e of each nodeiWith imaginary part correction amount fi；The fitness function is the inverse of objective function；
4-7) real and imaginary parts of all node voltages are modified respectively:
ei=ei+Δei, fi=fi+Δfi
4-8) judge whether the fitness value of each reactive compensation alternate node is greater than the inverse of target function value:
If it is, calculating the reactive compensation amount Δ Q of each reactive compensation alternate nodeCi[k]=Qi[k]-qi[k], export Δ QCi[k]Make
To configure the Optimization Compensation of reactive capability needed for each reactive compensation alternate node as a result, optimization terminates；If it is not, then according to every
The compensation capacity correction value Δ q of a reactive compensation alternate nodeCi[k], update QCi[k]=QCi[k]+ΔqCi[k], then return to
Step 4-4), until the fitness value of all reactive compensation alternate nodes is greater than fitness value index, it is standby to calculate each reactive compensation
Select the reactive compensation amount Δ Q of nodeCi[k]=Qi[k]-qi[k], export Δ QCi[k]Match as needed for each reactive compensation alternate node
The Optimization Compensation of reactive capability is set as a result, optimization terminates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910756767.XA CN110492493B (en) | 2019-08-16 | 2019-08-16 | Reactive compensation configuration optimization method for power system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910756767.XA CN110492493B (en) | 2019-08-16 | 2019-08-16 | Reactive compensation configuration optimization method for power system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110492493A true CN110492493A (en) | 2019-11-22 |
| CN110492493B CN110492493B (en) | 2020-10-27 |
Family
ID=68551288
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910756767.XA Expired - Fee Related CN110492493B (en) | 2019-08-16 | 2019-08-16 | Reactive compensation configuration optimization method for power system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110492493B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112887345A (en) * | 2019-11-29 | 2021-06-01 | 上海交通大学 | Node load balancing scheduling method for edge computing environment |
| CN118157150A (en) * | 2024-05-09 | 2024-06-07 | 江苏多端科技有限公司 | Voltage optimization analysis method and system based on reactive resource aggregation technology of power distribution network |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20120028831A (en) * | 2010-09-15 | 2012-03-23 | 이광직 | The apparatus and method of distribution transformer attached a condenser controll moudle with switching |
| CN105610168A (en) * | 2016-02-03 | 2016-05-25 | 南方电网科学研究院有限责任公司 | Application method of node voltage stability index of multi-feed-in alternating current and direct current system |
| CN105914755A (en) * | 2016-04-21 | 2016-08-31 | 广州供电局有限公司 | Electric power system dynamic reactive power optimization configuration method and system |
| CN105939020A (en) * | 2016-05-14 | 2016-09-14 | 国家电网公司 | Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability |
| CN109412169A (en) * | 2018-09-05 | 2019-03-01 | 贵州电网有限责任公司凯里供电局 | Based on the reactive power compensation ability encouraged by force under mode and the anti-freeze method of power distribution network |
-
2019
- 2019-08-16 CN CN201910756767.XA patent/CN110492493B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20120028831A (en) * | 2010-09-15 | 2012-03-23 | 이광직 | The apparatus and method of distribution transformer attached a condenser controll moudle with switching |
| CN105610168A (en) * | 2016-02-03 | 2016-05-25 | 南方电网科学研究院有限责任公司 | Application method of node voltage stability index of multi-feed-in alternating current and direct current system |
| CN105914755A (en) * | 2016-04-21 | 2016-08-31 | 广州供电局有限公司 | Electric power system dynamic reactive power optimization configuration method and system |
| CN105939020A (en) * | 2016-05-14 | 2016-09-14 | 国家电网公司 | Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability |
| CN109412169A (en) * | 2018-09-05 | 2019-03-01 | 贵州电网有限责任公司凯里供电局 | Based on the reactive power compensation ability encouraged by force under mode and the anti-freeze method of power distribution network |
Non-Patent Citations (2)
| Title |
|---|
| MIN TANG ET AL.: "Optimal Placement of Dynamic Reactive Power Compensation Device for Improving Immunity to Commutation Failure in Muti-Infeed HVDC System", 《IEEE》 * |
| 彭思华 等: "配电网无功配置基态分析规划法", 《电网技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112887345A (en) * | 2019-11-29 | 2021-06-01 | 上海交通大学 | Node load balancing scheduling method for edge computing environment |
| CN118157150A (en) * | 2024-05-09 | 2024-06-07 | 江苏多端科技有限公司 | Voltage optimization analysis method and system based on reactive resource aggregation technology of power distribution network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110492493B (en) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110729765A (en) | Distribution network flexibility evaluation index system considering SOP and optimal scheduling method | |
| CN102280889B (en) | Method for reactive power optimization of electric power system on basis of clone-particle swarm hybrid algorithm | |
| CN105226708B (en) | A kind of large scale electric network operation plan equivalence coordination optimizing method | |
| CN111049173B (en) | Self-organizing droop control method for multi-terminal direct-current distribution network | |
| CN110601198B (en) | Hybrid micro-grid optimized operation method considering harmonic and voltage unbalance constraints | |
| WO2019056753A1 (en) | Dynamic equivalent modeling method for distributed photovoltaic power station cluster | |
| CN107947175B (en) | Micro-grid economic dispatching method based on distributed network control | |
| CN106374513B (en) | A kind of more microgrid dominant eigenvalues optimization methods based on leader-followers games | |
| CN106058858B (en) | Power distribution network optimization method and device | |
| Zhao et al. | Robust distributed coordination of parallel restored subsystems in wind power penetrated transmission system | |
| CN109494727A (en) | Consider the active and idle coordination optimization operation method of power distribution network of demand response | |
| CN110492493A (en) | A kind of power system reactive power compensation method for optimizing configuration | |
| CN109034587A (en) | A kind of active distribution system Optimization Scheduling for coordinating a variety of controllables | |
| CN107546773A (en) | A kind of more micro-capacitance sensor dynamic networking methods in region based on graph theory | |
| CN108376984B (en) | Automatic adjustment method for node voltage of power system | |
| CN112598159A (en) | Island fusion recovery method based on multi-source network reconstruction | |
| CN109377020B (en) | Power transmission network planning method considering load transfer capacity of power distribution network | |
| CN109378819A (en) | The power system voltage partition method and system verified based on spectral clustering and modularity | |
| CN117039882A (en) | Resource aggregation regulation and control method and system based on binary consistency algorithm | |
| CN116054213A (en) | Distributed energy storage group cooperative voltage control method based on cluster division | |
| CN108933448A (en) | A kind of low and medium voltage distribution network control method for coordinating and system containing photo-voltaic power supply | |
| CN111724064A (en) | Energy-storage-containing power distribution network planning method based on improved immune algorithm | |
| CN108376997B (en) | Active power distribution network island division method considering distributed power supply uncertainty | |
| Ellithy et al. | Optimal shunt capacitors allocation in distribution networks using genetic algorithm-practical case study | |
| CN113346501B (en) | Power distribution network voltage optimization method and system based on brainstorming algorithm |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201027 Termination date: 20210816 |