CN110011313A - A kind of flexible direct current power grid load flow calculation method and system - Google Patents
A kind of flexible direct current power grid load flow calculation method and system Download PDFInfo
- Publication number
- CN110011313A CN110011313A CN201910129301.7A CN201910129301A CN110011313A CN 110011313 A CN110011313 A CN 110011313A CN 201910129301 A CN201910129301 A CN 201910129301A CN 110011313 A CN110011313 A CN 110011313A
- Authority
- CN
- China
- Prior art keywords
- direct current
- node
- flexible direct
- power grid
- current 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.)
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Abstract
The invention discloses a kind of flexible direct current power grid load flow calculation methods, comprising: calculates the node admittance matrix of the flexible direct current power grid；The voltage initial value of the node of the flexible direct current power grid is set, and obtains maximum number of iterations；The error vector of converter station in the flexible direct current power grid is calculated, and judges whether the error vector meets the margin of error；If so, trend distribution is calculated according to node voltage, to export the result of Load flow calculation；If it is not, the node voltage on each polar curve road is then corrected, until the error vector meets the margin of error or times of revision reaches the maximum number of iterations；When reaching the maximum number of iterations, the result of Load flow calculation is exported as not coverage information.Invention additionally discloses a kind of flexible direct current electric network swim computing systems.Using the embodiment of the present invention, the trend distribution under system difference operational mode can be acquired with degree of precision, versatility is high.
Description
Technical field
The present invention relates to electric power system protection and control field more particularly to a kind of flexible direct current power grid load flow calculation methods
And system.
Background technique
Flexible DC transmission technology is can to turn off power electronic devices as core using voltage source converter and full-control type
New DC transmission mode is smart grid research and representative one of the key technology of application field.It has built at present
At flexible DC transmission engineering largely using the pseudo- bipolar mode of connection, research is suitable for the tide of pseudo- bipolar flexible DC grid
Flow calculation methodologies are the bases that control is optimized to trend.Existing puppet bipolar flexible DC grid tidal current computing method is only fitted
For running on the network system of symmetric pattern, i.e., tidal current computing method poor universality in the prior art is not suitable for symmetrical
The scene that mode and any bar route are removed by failure.
Summary of the invention
The purpose of the embodiment of the present invention is that a kind of flexible direct current power grid load flow calculation method and system are provided, it can be with more high-precision
Degree acquires the distribution of the trend under system difference operational mode, and versatility is high.
To achieve the above object, the embodiment of the invention provides a kind of flexible direct current power grid load flow calculation methods, comprising:
Calculate the node admittance matrix of the flexible direct current power grid；
The node voltage initial value of the flexible direct current power grid is set, and obtains maximum number of iterations；
The error vector of converter station in the flexible direct current power grid is calculated, and judges whether the error vector meets maximum
Allowable error；
If so, trend distribution is calculated according to node voltage, to export the result of Load flow calculation；If it is not, then solve it is refined can
Than matrix and the node voltage on each polar curve road is corrected, until the error vector meets the margin of error or times of revision reaches
The maximum number of iterations；
When reaching the maximum number of iterations, the result of Load flow calculation is exported as not coverage information.
Compared with prior art, flexible direct current power grid load flow calculation method disclosed by the invention, firstly, calculating flexible direct current
The node admittance matrix of power grid；Then, the error vector of converter station in flexible direct current power grid is calculated, and judges the error vector
Whether the margin of error is met；If satisfied, then export Load flow calculation as a result, if not satisfied, then solving Jacobian matrix simultaneously
The node voltage on each polar curve road is corrected, and when reaching the maximum number of iterations, the result for exporting Load flow calculation is not restrain
Information.It solves tidal current computing method poor universality in the prior art, is not suitable for symmetric pattern and any bar route for some reason
The problem of scene that barrier is removed, the trend distribution under system difference operational mode can be acquired with degree of precision, versatility is high.
As an improvement of the above scheme, the node admittance matrix for calculating the flexible direct current power grid meets following public
Formula:
Wherein, i and j indicates that node serial number, j ∈ i indicate that node i is associated with node j by route；K=p, n；Work as k=
When p, the corresponding electrical quantity of positive route is indicated, as k=n, indicate the corresponding electrical quantity of anode circuit；Indicate connection section
The conductance of the pole the k transmission line of electricity of point i and j.
As an improvement of the above scheme, the error vector for calculating converter station in the flexible direct current power grid, meet with
Lower formula:
Wherein, Δ WiFor the error vector of converter station i；K1And K2For control coefrficient；Work as K1=0, K2When=1, it is fixed to be equivalent to
DC voltage control；Work as K1=1, K2When=0, it is equivalent to and determines active power controller；K1≠ 0, K2When ≠ 0, it is equivalent to sagging control
System, K1And K2Respectively indicate the corresponding power coefficient of droop control device and voltage coefficient；Pi *Indicate that converter station i corresponds to controller
Active power reference value；Indicate that converter station i corresponds to the direct voltage reference value of controller；Indicate the change of current
The voltage-to-ground of the corresponding pole the k outlet of the i that stands, p are anode, and n is cathode；PiIndicate that converter station i is received active from AC system
Power.
As an improvement of the above scheme, described to solve Jacobian matrix and correct the node voltage on each polar curve road, meet with
Lower formula:
Wherein, Δ Up(k)For the correction amount vector of cathode node voltage to be solved；ΔUn(k)For negative pole node to be solved
The correction amount vector of voltage；Jp(k)And Jn(k)Respectively corresponding Jacobian matrix.
As an improvement of the above scheme, described that trend distribution is calculated according to node voltage, meet following formula:
Wherein,Indicate the load current of the pole the k transmission line of electricity of connecting node i and j.
As an improvement of the above scheme, the node voltage initial value of the setting flexible direct current power grid, specifically includes:
According to the node voltage initial value of the voltage rating of flexible direct current power grid setting anode and anode circuit.
To achieve the above object, the embodiment of the invention also provides a kind of flexible direct current electric network swim computing systems, comprising:
Node admittance matrix computing unit, for calculating the node admittance matrix of the flexible direct current power grid；
Node voltage initial value design unit, for setting the node voltage initial value of the flexible direct current power grid；
Maximum number of iterations acquiring unit, for obtaining maximum number of iterations；
Error vector computing unit, for calculating the error vector of converter station in the flexible direct current power grid；
Judging unit, for judging whether the error vector meets the margin of error；
Load flow calculation unit, for calculating tide according to node voltage when the error vector meets the margin of error
Flow distribution, to export the result of Load flow calculation；It is also used to when the error vector is unsatisfactory for the margin of error, solution is refined can
Than matrix and the node voltage on each polar curve road is corrected, until the error vector meets the margin of error or times of revision reaches
The maximum number of iterations；When reaching the maximum number of iterations, the result of Load flow calculation is exported as not coverage information.
Compared with prior art, flexible direct current electric network swim computing system disclosed by the invention, firstly, node admittance matrix
The node admittance matrix of computing unit calculating flexible direct current power grid；Then, error vector computing unit calculates flexible direct current power grid
The error vector of middle converter station, judging unit judge whether the error vector meets the margin of error；If satisfied, then trend
Computing unit output Load flow calculation as a result, if not satisfied, then Load flow calculation unit solves and Jacobian matrix and corrects each polar curve
The node voltage on road, and when reaching the maximum number of iterations, the result of Load flow calculation unit output Load flow calculation is not receive
Hold back information.Solve tidal current computing method poor universality in the prior art, be not suitable for symmetric pattern and any bar route because
The problem of scene that failure is removed, the trend distribution under system difference operational mode can be acquired with degree of precision, versatility is high.
As an improvement of the above scheme, the node admittance matrix for calculating the flexible direct current power grid meets following public
Formula:
Wherein, i and j indicates that node serial number, j ∈ i indicate that node i is associated with node j by route；K=p, n；Work as k=
When p, the corresponding electrical quantity of positive route is indicated, as k=n, indicate the corresponding electrical quantity of anode circuit；Indicate connection section
The conductance of the pole the k transmission line of electricity of point i and j.
As an improvement of the above scheme, the error vector for calculating converter station in the flexible direct current power grid, meet with
Lower formula:
Wherein, Δ WiFor the error vector of converter station i；K1And K2For control coefrficient；Work as K1=0, K2When=1, it is fixed to be equivalent to
DC voltage control；Work as K1=1, K2When=0, it is equivalent to and determines active power controller；K1≠ 0, K2When ≠ 0, it is equivalent to sagging control
System, K1And K2Respectively indicate the corresponding power coefficient of droop control device and voltage coefficient；Pi *Indicate that converter station i corresponds to controller
Active power reference value；Indicate that converter station i corresponds to the direct voltage reference value of controller；Indicate the change of current
The voltage-to-ground of the corresponding pole the k outlet of the i that stands, p are anode, and n is cathode；PiIndicate that converter station i is received active from AC system
Power.
As an improvement of the above scheme, the node voltage initial value design unit is specifically used for:
According to the node voltage initial value of the voltage rating of flexible direct current power grid setting anode and anode circuit.
Detailed description of the invention
Fig. 1 is a kind of flow chart of flexible direct current power grid load flow calculation method provided in an embodiment of the present invention；
Fig. 2 is a kind of pseudo- bipolar flexible DC grid typical structure schematic diagram provided in an embodiment of the present invention；
Fig. 3 is a kind of another flow chart of flexible direct current power grid load flow calculation method provided in an embodiment of the present invention；
Fig. 4 is a kind of structural block diagram of flexible direct current electric network swim computing system provided in an embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Embodiment one
It is a kind of flow chart of flexible direct current power grid load flow calculation method provided in an embodiment of the present invention referring to Fig. 1, Fig. 1；
Include:
S1, the node admittance matrix for calculating the flexible direct current power grid；
The node voltage initial value of S2, the setting flexible direct current power grid, and obtain maximum number of iterations；
S3, the error vector for calculating converter station in the flexible direct current power grid, and judge whether the error vector meets
The margin of error；
S4, it is distributed if so, calculating trend according to node voltage, to export the result of Load flow calculation；If it is not, then solving refined
Comparable matrix and the node voltage for correcting each polar curve road, until the error vector meets the margin of error or times of revision reaches
To the maximum number of iterations；
S5, when reaching the maximum number of iterations, the result for exporting Load flow calculation is not coverage information.
It is worth noting that the applicable object of the embodiment of the present invention is the flexible direct current power grid for using pseudo- bipolar wiring, and
Grounding point is arranged in the exchange side of converter station or DC side at one end, provides reference potential for DC equipment.Referring to fig. 2, Fig. 2 is this
A kind of pseudo- bipolar flexible DC grid typical structure schematic diagram that inventive embodiments provide；Specific converter station form is on the right side of figure
It is described in detail, DC line includes the two poles of the earth altogether: anode and cathode.In order to fix the above earth potential of the various equipment of DC side, select
The DC side of converter station 5 passes through big resistance eutral grounding.
Flexible direct current power grid load flow calculation method described in the embodiment of the present invention can acquire system difference with degree of precision
Trend distribution under operational mode；Wherein, the operational mode includes: that system operates normally, positive route is removed, negative line
Road is removed, positive route and anode circuit are removed isotype.
Specifically, in step sl, the node admittance matrix for calculating the flexible direct current power grid meets following public
Formula:
Wherein, i and j indicates that node serial number, j ∈ i indicate that node i is associated with node j by route；K=p, n；Work as k=
When p, the corresponding electrical quantity of positive route is indicated, as k=n, indicate the corresponding electrical quantity of anode circuit；Indicate connection section
The conductance of the pole the k transmission line of electricity of point i and j.
Specifically, in step s 2, the node voltage initial value of the setting flexible direct current power grid specifically includes: root
The node voltage initial value U of positive route is set according to the voltage rating of the flexible direct current power gridp(0)With the node electricity of anode circuit
Press initial value Un(0).Wherein, the voltage rating is determined according to engineering is practical, if DC voltage grade is ± 500kV, then positive
Node voltage initial value is set as 500kV, and negative pole node voltage initial value is set as -500kV.The maximum number of iterations KmaxGenerally
It rule of thumb sets, such as may be configured as 1000, if reaching maximum number of iterations KmaxThere are no the margin of error is met, then
Expression does not restrain.
Specifically, in step s3, the error vector for calculating converter station in the flexible direct current power grid meets following
Formula:
Wherein, Δ WiFor the error vector of converter station i；K1And K2For control coefrficient；Work as K1=0, K2When=1, it is fixed to be equivalent to
DC voltage control；Work as K1=1, K2When=0, it is equivalent to and determines active power controller；K1≠ 0, K2When ≠ 0, it is equivalent to sagging control
System, K1And K2Respectively indicate the corresponding power coefficient of droop control device and voltage coefficient；Pi *Indicate that converter station i corresponds to controller
Active power reference value；Indicate that converter station i corresponds to the direct voltage reference value of controller；Indicate the change of current
The voltage-to-ground of the corresponding pole the k outlet of the i that stands, p are anode, and n is cathode；PiIndicate that converter station i is received active from AC system
Power.
Specifically, in step s 4, as the error vector Δ W(k)When meeting margin of error ε (meet | Δ W(k)
| < ε), trend distribution is calculated according to node voltage, to export the result of Load flow calculation.At this point, described calculate according to node voltage
Trend distribution, meets following formula:
Wherein,Indicate the load current of the pole the k transmission line of electricity of connecting node i and j.
When the error vector is unsatisfactory for the margin of error, solves Jacobian matrix and correct the node on each polar curve road
Voltage, until the error vector meets the margin of error or times of revision reaches the maximum number of iterations.Preferably, institute
The node voltage for solving Jacobian matrix and correcting each polar curve road is stated, following formula is met:
Wherein, Δ Up(k)For the correction amount vector of cathode node voltage to be solved；ΔUn(k)For negative pole node to be solved
The correction amount vector of voltage；Jp(k)And Jn(k)Respectively corresponding Jacobian matrix.
The process for specifically solving the Jacobian matrix is as follows:
Wherein, G indicates the set of the converter station of exchange side or DC side ground connection.
Simultaneous formula (5)~(7) and formula (10)~(11), can acquire positive and negative electrode line node voltage increment,
Each polar curve circuit node voltage is updated on the basis of this and carries out next step iteration, until meeting the margin of error or reaching greatest iteration
Number.
Specifically, in step s 5, when reaching the maximum number of iterations, the result for exporting Load flow calculation is not restrain
Information.
Preferably, the process of above-mentioned steps S1~S5 can also refer to Fig. 3, and Fig. 3 is one kind provided in an embodiment of the present invention
Another flow chart of flexible direct current power grid load flow calculation method.
When it is implemented, firstly, calculating the node admittance matrix of flexible direct current power grid；Then, flexible direct current power grid is calculated
The error vector of middle converter station, and judge whether the error vector meets the margin of error；If satisfied, then exporting trend meter
It is calculating as a result, if not satisfied, then solve Jacobian matrix and correct the node voltage on each polar curve road, and change when reach the maximum
When generation number, the result of Load flow calculation is exported as not coverage information.
Compared with prior art, flexible direct current power grid load flow calculation method disclosed by the invention, solves in the prior art
Tidal current computing method poor universality, be not suitable for symmetric pattern and any bar route because failure be removed scene the problem of,
The trend distribution under system difference operational mode can be acquired with degree of precision, versatility is high.
Embodiment two
Referring to fig. 4, Fig. 4 is a kind of structural frames of flexible direct current electric network swim computing system provided in an embodiment of the present invention
Figure；Include:
Node admittance matrix computing unit 11, for calculating the node admittance matrix of the flexible direct current power grid；
Node voltage initial value design unit 12, for setting the node voltage initial value of the flexible direct current power grid；
Maximum number of iterations acquiring unit 13, for obtaining maximum number of iterations；
Error vector computing unit 14, for calculating the error vector of converter station in the flexible direct current power grid；
Judging unit 15, for judging whether the error vector meets the margin of error；
Load flow calculation unit 16, for being calculated according to node voltage when the error vector meets the margin of error
Trend distribution, to export the result of Load flow calculation；It is also used to when the error vector is unsatisfactory for the margin of error, solves refined
Comparable matrix and the node voltage for correcting each polar curve road, until the error vector meets the margin of error or times of revision reaches
To the maximum number of iterations；When reaching the maximum number of iterations, the result of Load flow calculation is exported as not coverage information.
Specifically, the node admittance matrix for calculating the flexible direct current power grid, meets following formula:
Wherein, i and j indicates that node serial number, j ∈ i indicate that node i is associated with node j by route；K=p, n；Work as k=
When p, the corresponding electrical quantity of positive route is indicated, as k=n, indicate the corresponding electrical quantity of anode circuit；Indicate connection section
The conductance of the pole the k transmission line of electricity of point i and j.
Specifically, the node voltage initial value design unit 12 is specifically used for: according to the specified of the flexible direct current power grid
Voltage sets the node voltage initial value U of positive routep(0)With the node voltage initial value U of anode circuitn(0).Wherein, the specified electricity
It is determined at all according to engineering is practical, if DC voltage grade is ± 500kV, then cathode node voltage initial value is set as 500kV, bears
Pole node voltage initial value is set as -500kV.The maximum number of iterations KmaxIt generally rule of thumb sets, such as may be configured as
1000, if reaching maximum number of iterations KmaxThere are no meet the margin of error, then it represents that does not restrain.
Specifically, the error vector for calculating converter station in the flexible direct current power grid, meets following formula:
Wherein, Δ WiFor the error vector of converter station i；K1And K2For control coefrficient；Work as K1=0, K2When=1, it is fixed to be equivalent to
DC voltage control；Work as K1=1, K2When=0, it is equivalent to and determines active power controller；K1≠ 0, K2When ≠ 0, it is equivalent to sagging control
System, K1And K2Respectively indicate the corresponding power coefficient of droop control device and voltage coefficient；Pi *Indicate that converter station i corresponds to controller
Active power reference value；Indicate that converter station i corresponds to the direct voltage reference value of controller；Indicate the change of current
The voltage-to-ground of the corresponding pole the k outlet of the i that stands, p are anode, and n is cathode；PiIndicate that converter station i is received active from AC system
Power.
Specifically, when the judging unit 15 determines the error vector Δ W(k)Meet (i.e. full when margin of error ε
Foot | Δ W(k)| < ε), the Load flow calculation unit 16 calculates trend according to node voltage and is distributed, to export the knot of Load flow calculation
Fruit.At this point, described calculate trend distribution according to node voltage, meet following formula:
Wherein,Indicate the load current of the pole the k transmission line of electricity of connecting node i and j.
When the judging unit 15 determines that the error vector is unsatisfactory for the margin of error, the Load flow calculation unit
16 solution Jacobian matrixes and the node voltage for correcting each polar curve road, until the error vector meets the margin of error or repairs
Positive number reaches the maximum number of iterations.Preferably, the node voltage for solving Jacobian matrix and correcting each polar curve road,
Meet following formula:
Wherein, Δ Up(k)For the correction amount vector of cathode node voltage to be solved；ΔUn(k)For negative pole node to be solved
The correction amount vector of voltage；Jp(k)And Jn(k)Respectively corresponding Jacobian matrix.
The process for specifically solving the Jacobian matrix is as follows:
Wherein, G indicates the set of the converter station of exchange side or DC side ground connection.
Simultaneous formula (5)~(7) and formula (10)~(11), can acquire positive and negative electrode line node voltage increment,
Each polar curve circuit node voltage is updated on the basis of this and carries out next step iteration, until meeting the margin of error or reaching greatest iteration
Number.
Specifically, the Load flow calculation unit 16 exports the result of Load flow calculation when reaching the maximum number of iterations
For not coverage information.
When it is implemented, firstly, node admittance matrix computing unit 11 calculates the node admittance matrix of flexible direct current power grid；
Then, error vector computing unit 14 calculates the error vector of converter station in flexible direct current power grid, and judging unit 15 judges described
Whether error vector meets the margin of error；If satisfied, then Load flow calculation unit 16 export Load flow calculation as a result, if discontented
Foot, then Load flow calculation unit 16 solves Jacobian matrix and corrects the node voltage on each polar curve road, and changes when reaching the maximum
When generation number, Load flow calculation unit 16 exports the result of Load flow calculation as not coverage information.
Compared with prior art, flexible direct current electric network swim computing system disclosed by the invention, solves in the prior art
Tidal current computing method poor universality, be not suitable for symmetric pattern and any bar route because failure be removed scene the problem of,
The trend distribution under system difference operational mode can be acquired with degree of precision, versatility is high.
The above is a preferred embodiment of the present invention, it is noted that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.
Claims (10)
1. a kind of flexible direct current power grid load flow calculation method characterized by comprising
Calculate the node admittance matrix of the flexible direct current power grid；
The node voltage initial value of the flexible direct current power grid is set, and obtains maximum number of iterations；
The error vector of converter station in the flexible direct current power grid is calculated, and it is maximum allowable to judge whether the error vector meets
Error；
If so, trend distribution is calculated according to node voltage, to export the result of Load flow calculation；If it is not, then solving Jacobean matrix
Battle array and the node voltage for correcting each polar curve road, until the error vector meets the margin of error or times of revision reaches described
Maximum number of iterations；
When reaching the maximum number of iterations, the result of Load flow calculation is exported as not coverage information.
2. flexible direct current power grid load flow calculation method as described in claim 1, which is characterized in that described to calculate the flexibility directly
The node admittance matrix of galvanic electricity net, meets following formula:
Wherein, i and j indicates that node serial number, j ∈ i indicate that node i is associated with node j by route；K=p, n；As k=p,
It indicates the corresponding electrical quantity of positive route, as k=n, indicates the corresponding electrical quantity of anode circuit；Indicate connecting node i with
The conductance of the pole the k transmission line of electricity of j.
3. flexible direct current power grid load flow calculation method as claimed in claim 2, which is characterized in that described to calculate the flexibility directly
The error vector of converter station in galvanic electricity net, meets following formula:
Wherein, △ WiFor the error vector of converter station i；K1And K2For control coefrficient；Work as K1=0, K2When=1, it is equivalent to and determines direct current
Voltage control；Work as K1=1, K2When=0, it is equivalent to and determines active power controller；K1≠ 0, K2When ≠ 0, it is equivalent to sagging control, K1With
K2Respectively indicate the corresponding power coefficient of droop control device and voltage coefficient；Pi *Indicate that converter station i corresponds to the wattful power of controller
Rate reference value；Indicate that converter station i corresponds to the direct voltage reference value of controller；(k=p, n) indicates that converter station i is corresponding
The voltage-to-ground of the pole k outlet, p are anode, and n is cathode；PiIndicate converter station i from the received active power of AC system.
4. flexible direct current power grid load flow calculation method as described in claim 1, which is characterized in that the solution Jacobian matrix
And the node voltage on each polar curve road is corrected, meet following formula:
Wherein, △ Up(k)For the correction amount vector of cathode node voltage to be solved；△Un(k)For negative pole node voltage to be solved
Correction amount vector；Jp(k)And Jn(k)Respectively corresponding Jacobian matrix.
5. flexible direct current power grid load flow calculation method as claimed in claim 3, which is characterized in that described according to node voltage meter
Trend distribution is calculated, following formula is met:
Wherein,Indicate the load current of the pole the k transmission line of electricity of connecting node i and j.
6. flexible direct current power grid load flow calculation method as described in claim 1, which is characterized in that the setting is described flexible straight
The node voltage initial value of galvanic electricity net, specifically includes:
According to the node voltage initial value of the voltage rating of flexible direct current power grid setting anode and anode circuit.
7. a kind of flexible direct current electric network swim computing system characterized by comprising
Node admittance matrix computing unit, for calculating the node admittance matrix of the flexible direct current power grid；
Node voltage initial value design unit, for setting the node voltage initial value of the flexible direct current power grid；
Maximum number of iterations acquiring unit, for obtaining maximum number of iterations；
Error vector computing unit, for calculating the error vector of converter station in the flexible direct current power grid；
Judging unit, for judging whether the error vector meets the margin of error；
Load flow calculation unit, for calculating trend point according to node voltage when the error vector meets the margin of error
Cloth, to export the result of Load flow calculation；It is also used to when the error vector is unsatisfactory for the margin of error, solves Jacobean matrix
Battle array and the node voltage for correcting each polar curve road, until the error vector meets the margin of error or times of revision reaches described
Maximum number of iterations；When reaching the maximum number of iterations, the result of Load flow calculation is exported as not coverage information.
8. flexible direct current electric network swim computing system as claimed in claim 7, which is characterized in that described to calculate the flexibility directly
The node admittance matrix of galvanic electricity net, meets following formula:
Wherein, i and j indicates that node serial number, j ∈ i indicate that node i is associated with node j by route；K=p, n；As k=p,
It indicates the corresponding electrical quantity of positive route, as k=n, indicates the corresponding electrical quantity of anode circuit；Indicate connecting node i with
The conductance of the pole the k transmission line of electricity of j.
9. flexible direct current electric network swim computing system as claimed in claim 8, which is characterized in that described to calculate the flexibility directly
The error vector of converter station in galvanic electricity net, meets following formula:
Wherein, △ WiFor the error vector of converter station i；K1And K2For control coefrficient；Work as K1=0, K2When=1, it is equivalent to and determines direct current
Voltage control；Work as K1=1, K2When=0, it is equivalent to and determines active power controller；K1≠ 0, K2When ≠ 0, it is equivalent to sagging control, K1With
K2Respectively indicate the corresponding power coefficient of droop control device and voltage coefficient；Pi *Indicate that converter station i corresponds to the wattful power of controller
Rate reference value；Indicate that converter station i corresponds to the direct voltage reference value of controller；(k=p, n) indicates that converter station i is corresponding
The voltage-to-ground of the pole k outlet, p are anode, and n is cathode；PiIndicate converter station i from the received active power of AC system.
10. flexible direct current electric network swim computing system as claimed in claim 7, which is characterized in that the node voltage initial value
Setup unit is specifically used for:
According to the node voltage initial value of the voltage rating of flexible direct current power grid setting anode and anode circuit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910129301.7A CN110011313A (en) | 2019-02-21 | 2019-02-21 | A kind of flexible direct current power grid load flow calculation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910129301.7A CN110011313A (en) | 2019-02-21 | 2019-02-21 | A kind of flexible direct current power grid load flow calculation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110011313A true CN110011313A (en) | 2019-07-12 |
Family
ID=67165764
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910129301.7A Pending CN110011313A (en) | 2019-02-21 | 2019-02-21 | A kind of flexible direct current power grid load flow calculation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110011313A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110492465A (en) * | 2019-07-29 | 2019-11-22 | 深圳供电局有限公司 | Tidal current computing method, device, computer equipment and readable storage medium storing program for executing |
| CN114819281A (en) * | 2022-03-29 | 2022-07-29 | 四川大学 | Method for optimizing coordinated power flow between flexible direct-current power grid stations |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106549384A (en) * | 2016-12-09 | 2017-03-29 | 国网江苏省电力公司经济技术研究院 | A kind of general tidal current computing method containing UPFC power systems |
| CN107769213A (en) * | 2017-10-31 | 2018-03-06 | 贵州电网有限责任公司电力科学研究院 | A kind of more transverters alternating current-direct current distribution power system load flow calculation method in parallel |
| CN109066693A (en) * | 2018-09-04 | 2018-12-21 | 河南理工大学 | A kind of flexiblesystem tidal current computing method with Higher Order of Convergence rate |
| CN109921406A (en) * | 2018-12-26 | 2019-06-21 | 天津大学 | True bipolar flexible DC grid is layered tidal current computing method |
-
2019
- 2019-02-21 CN CN201910129301.7A patent/CN110011313A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106549384A (en) * | 2016-12-09 | 2017-03-29 | 国网江苏省电力公司经济技术研究院 | A kind of general tidal current computing method containing UPFC power systems |
| CN107769213A (en) * | 2017-10-31 | 2018-03-06 | 贵州电网有限责任公司电力科学研究院 | A kind of more transverters alternating current-direct current distribution power system load flow calculation method in parallel |
| CN109066693A (en) * | 2018-09-04 | 2018-12-21 | 河南理工大学 | A kind of flexiblesystem tidal current computing method with Higher Order of Convergence rate |
| CN109921406A (en) * | 2018-12-26 | 2019-06-21 | 天津大学 | True bipolar flexible DC grid is layered tidal current computing method |
Non-Patent Citations (2)
| Title |
|---|
| YICHAO LIU等: ""Research on Power Flow Calculation and Optimization Method of Real Bipolar VSC-HVDC Grid under Asymmetrical Mode"", 《2018 IEEE POWER & ENERGY SOCIETY GENERAL MEETING (PESGM)》 * |
| 钱甜甜 等: ""基于节点电流关系的多端柔性直流输电系统潮流计算"", 《电网技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110492465A (en) * | 2019-07-29 | 2019-11-22 | 深圳供电局有限公司 | Tidal current computing method, device, computer equipment and readable storage medium storing program for executing |
| CN114819281A (en) * | 2022-03-29 | 2022-07-29 | 四川大学 | Method for optimizing coordinated power flow between flexible direct-current power grid stations |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105356481B (en) | A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes | |
| CN108429257B (en) | Distribution network frame planning method based on load point clustering partitioning and considering geographic factors | |
| US20140129195A1 (en) | Real time dynamic physics simulation device of flexible dc transmission system | |
| WO2022016622A1 (en) | Adaptive optimization and control method in event of failure of true bipolar flexible direct-current power transmission system | |
| CN102273039A (en) | Photovoltaic power plant having an offset voltage source controlling the DC potential at the inverter output | |
| CN105119267B (en) | A kind of power grid static voltage stability horizontal analysis method | |
| CN107508322B (en) | A kind of short-circuit fault of power system computational methods for considering photovoltaic electric station grid connection | |
| CN110489806B (en) | Electromagnetic transient modeling and calculating method for power grid comprising multi-voltage-source converter | |
| CN109638839B (en) | Load flow calculation method of bipolar flexible direct-current transmission system | |
| CN105119316A (en) | VSC-MTDC (Voltage source converter-multi-terminal high voltage direct current) voltage control method for offshore wind farm integration | |
| CN110011313A (en) | A kind of flexible direct current power grid load flow calculation method and system | |
| CN106786716B (en) | A kind of method and device for the degree that influences each other between calculating direct current | |
| CN105277773A (en) | Method for completely calculating grounding short-circuit current of transformation station | |
| CN108599227B (en) | MMC direct-current voltage balance control method for forming direct-current converter station based on MMC cascade connection | |
| CN108521141A (en) | It is a kind of meter and wind power plant voltage's distribiuting characteristic short-circuit current calculation method | |
| Jena et al. | Setting a fostered energy network by decarbonizing the grid: H ybridization, control, and future solutions upon storage | |
| CN108199381A (en) | A kind of cutting load place selection method for wind power integration system voltage stability | |
| CN105515034B (en) | A kind of Poewr control method of the bipolar MMC-HVDC systems in both ends | |
| CN106972526A (en) | The method that virtual plant carries out auto-control according to power supply capacity and network load | |
| CN107017658A (en) | The control method that virtual plant is self-regulated according to prediction power load | |
| CN111047346A (en) | Local area collaborative reactive compensation service value making method | |
| CN109599870A (en) | A kind of interconnection reactance quick calculation method inhibiting short circuit current | |
| CN110212558A (en) | A kind of addressing constant volume configuration method of distributed energy storage system | |
| CN109687499A (en) | A kind of voltage balancing control method of flexible direct current converter station series connection converter valve | |
| CN108988401B (en) | Method for evaluating power transmission capacity and constraint conditions of multi-channel alternating current-direct current series-parallel connection section |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210721 Address after: No. 11, Kexiang Road, Science City, Huangpu District, Guangzhou City, Guangdong Province Applicant after: China Southern Power Grid Research Institute Co.,Ltd. Applicant after: Tianjin University Address before: No. 11, Kexiang Road, Science City, Huangpu District, Guangzhou City, Guangdong Province Applicant before: China Southern Power Grid Research Institute Co.,Ltd. Applicant before: CSG POWER GRID TECHNOLOGY RESEARCH CENTER Applicant before: Tianjin University |
|
| TA01 | Transfer of patent application right | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190712 |
|
| RJ01 | Rejection of invention patent application after publication |