CN106953348A - A kind of Power network fault calculation method containing back-to-back DC power transmission - Google Patents
A kind of Power network fault calculation method containing back-to-back DC power transmission Download PDFInfo
- Publication number
- CN106953348A CN106953348A CN201710285879.2A CN201710285879A CN106953348A CN 106953348 A CN106953348 A CN 106953348A CN 201710285879 A CN201710285879 A CN 201710285879A CN 106953348 A CN106953348 A CN 106953348A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- voltage
- current
- formula
- component
- 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/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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Emergency Protection Circuit Devices (AREA)
Abstract
The invention discloses a kind of Power network fault calculation method containing back-to-back (BTB) direct current transportation, comprise the steps:Input the systematic parameter and control parameter of back-to-back transmission system;Determine the exchange of failure side network topology structure and failure boundary, the network topology structure of non-faulting top-cross streaming system, the constraint equation of the straight-flow system of both sides;Set up the Inverter circuit switch function model of the straight-flow system at two endsWith voltage switch function model;Combination failure top-cross stream constraint equation, direct current constraint equation and transverter switch function solve failure side power frequency busbar voltage;Iterative calculation, exits calculating beyond maximum magnitude when variable error meets calculating requirement or iterations, exports BTB systems two ends AC network change of current bus power-frequency voltage result of calculation.Calculating speed of the present invention is fast, the degree of accuracy is high, strong applicability, has very strong engineering practical value in terms of the AC network calculation of fault and its setting value order of relay protection containing back-to-back DC power transmission.
Description
Technical field
The present invention relates to the research field of electric power system fault computational methods, and in particular to one kind contains back-to-back DC power transmission
Power network fault calculation method.
Background technology
As voltage class and construction scale are the western Shandong back-to-back DC engineering of the first in the world compared with similar projects
Put into operation, indicate that the application of the back-to-back engineering of China enters the new stage.As between bulk power grid asynchronous networking it is important
Mode, back-to-back DC can effectively reduce both sides and connect influencing each other between power network;However, compared with customary DC, the back of the body
Backrest direct current is due to without DC power transmission line so that the electrical link between its rectification and inversion is close, thus both sides AC system
Between there is certain failure interaction.I.e. when side AC network breaks down, pass through the friendship with back-to-back DC system
The interaction of stream-DC-AC, offside power network will be influenceed by certain.Therefore, power network event of the research containing back-to-back DC
Barrier computational methods will establish important foundation for the accident analysis of such power network.
Failure interaction and its accident analysis problem have obtained the extensive of experts and scholars all the time between ac and dc systemses
Concern, related research work has been carried out for straight-flow system modeling to fault calculation methods for transmission etc..On straight-flow system model,
There is research to establish the transverter model based on switch function;Further research work then considers DC control system
Failure response, establishes straight-flow system Equivalent Model;The calculation of fault of alternating current-direct current power network can be realized using these models.But on
Stating research, to be all based on customary DC transmission line of electricity longer, and only this side current conversion station is impacted during the AC network failure of side,
And opposite side current conversion station is substantially unaffected this precondition.It is electric due to back-to-back DC rectifier and inverter
Closely, the above-mentioned analysis method for considering failure side current conversion station is not particularly suited for back-to-back DC power transmission system for system.And mesh
The preceding research for back-to-back DC then focuses primarily upon analysis and optimization to engineering proposal research and control system, exchange
Side protective value analysis and the back-to-back DC commutation failure analysis of causes etc., there is no and be related to the alternating current-direct current containing back-to-back DC
Electric network fault is analyzed.Therefore, it is necessary to propose a kind of Power network fault calculation method containing back-to-back DC power transmission.
The content of the invention
It is a primary object of the present invention to fill up the blank of this research field, overcome the shortcoming and deficiency of prior art, carry
Go out a kind of Power network fault calculation method containing back-to-back DC power transmission.
In order to achieve the above object, the present invention uses following technical scheme:
A kind of Power network fault calculation method containing back-to-back DC power transmission, comprises the steps:
S1, the back-to-back transmission system of input systematic parameter and control parameter;
S2, form by the network topology structure and failure boundary condition of failure top-cross streaming system the constraint of failure top-cross streaming system
Equation, is formed the AC system constraint equation of the side by the network topology structure of non-faulting top-cross streaming system, by the direct current of both sides
The direct current constraint equation that control characteristic is set up between its DC side DC voltage and DC current respectively;
S3, failure side and non-faulting side change of current busbar voltage areWherein k is the secondary of interative computation
Number;Whole variables are put with initial value, i.e. failure side and change of current busbar voltage in non-faulting side is
S4, set up according to AC power-frequency voltage, DC side current dc component and initial trigger angle voltage x current switch
Function model, obtains the switch function of voltage, electric currentWith
S5, by failure side power frequency busbar voltageThe voltage of failure side straight-flow system is tried to achieve with reference to transverter switch function
DC componentAnd second harmonic componentSo as to try to achieve the fundamental component of DC side electric currentAnd second harmonic component
S6, with reference to transverter switch function and non-faulting top-cross stream constraint equation try to achieve non-faulting side power frequency busbar voltageAnd thus correct the fundamental component of DC sideAnd second harmonic component
S7, try to achieve the electric current that DC side is injected into the AC network of failure side
S8, the electric current in the AC network of injection failure sideThe exchange constraint equation of combination failure side can try to achieve exchange
Try to achieve failure side change of current busbar voltage
S9, judgementWithDifference andWithDifference whether be satisfied by convergence
Condition, if being unsatisfactory forSubstitute into respectivelyAnd S4~S9 steps are repeated, until
Restrain or reach that maximum iteration backed off after random is calculated, and export result of calculation.
It is preferred that, in step S2, the sequence network equation in the non-faulting top-cross stream constraint equation is as follows:
And then formula (2) is shown for failure boundary conditional equation.
Formula (1), (2) constitute failure top-cross stream constraint equation.In the formula of the above two,AndRespectively event
Hinder positive sequence, negative phase-sequence and the zero-sequence component of side change of current bus fundamental frequency voltages.For the electric current positive-sequence component of AC network equivalent power supply,Respectively DC side is injected into the Fundamental-frequency Current positive sequence and negative sequence component of AC network;WithRespectively
Positive-sequence component, negative sequence component and the zero-sequence component of fault current at failure;
For non-faulting side exchange constraint equation can be obtained according to the topological structure and principle of stacking of its AC network as
Lower expression formula:
In formula (3),ForAC network power frequency order impedance matrix during independent role,ForDuring independent role
AC network power frequency positive sequence admittance matrix.
For the direct current set up respectively between its DC side DC voltage and DC current according to the DC control characteristic of both sides
Constraint equation is expressed as follows:
Idi0=fi(Udi0) (4)
Idr0=fr(Udr0) (5)
In formula (4), (5), Udi0、Idi0、Udr0、Idr0The DC component of rectification side DC voltage and electric current is expressed as,
The DC voltage of inverter side and the DC component of electric current, subscript " i " and " r " wherein in following table represent inverter side and rectification respectively
Side;Ud0、Id0Then refer to the voltage DC component and current dc component of straight-flow system.
It is preferred that, in step S3, k is the number of times of interative computation, and when calculation error is unsatisfactory for requiring, k values Jia 1, is continued
Cycle calculations.
It is preferred that, in step S4, built according to AC power-frequency voltage, DC side current dc component and initial trigger angle
Vertical voltage x current switch function model, the calculating of the transverter voltage x current switch function power frequency order components is specific as follows:
A) according to the definition of Park Transformation, change of current bus can be converted to a, b, c three-phase voltage by positive and negative, residual voltage:
In formula (6),Respectively the fundamental frequency phase voltage of a, b, c phase of change of current bus, a=ej2π/3.By phase
Voltage and then try to achieve corresponding fundamental frequency line voltageWithAnd then obtain the skew of synchronizing voltage phase
WillWithThe α components and β components of commutation voltage are represented respectively, are counted by formula (7)
Calculate:
The phase of DC control system synchronizing voltage can be tried to achieve using the α components and β components of commutation voltage
WillPhase angle bePhase angle bePhase angle beIt can be calculated by following formula
The phase offset of synchronizing voltage
In formula (9),For ca phases and the phase offset of synchronizing voltage,Phase for ab phases and synchronizing voltage is inclined
Move,For bc phases and the phase offset of synchronizing voltage;
B) according to AC power-frequency voltageWithDC side current dc component Id0And initial trigger angle
α0, calculate converter valve turn on delay angle θmn, actual Trigger Angle αmnWith actual angle of overlap μmn;
Turn on delay angle θmnComputing formula be:
In formula (10), xy=ab, bc, ca, wherein a, b, c represent the phase in three-phase respectively, similarly hereinafter;
Actual Trigger Angle αxyComputing formula be:
In formula (10) and (11), all angles are with delayed for just, advanced is negative;
Angle of overlap μ during xy two-phase commutationsxyComputing formula be:
μxy=cos-1(cosαxy-2XrId0/Uxy1)-αxy (12)
In formula (12), XrFor the commutating reactance of straight-flow system transverter;
C) according to θxy、αxyAnd μxyThree-phase voltage current switching waveform is made, Fu is utilized by the three-phase voltage current waveform
In leaf-size class number derive each order component of voltage x current switch function:
In formula (13) and (14),Respectively k order components of three-phase voltage switch function, Respectively k order components of three-phase current switch function, k=0,1,2,3 ..., T is 2 π;
Set up the order components of voltage x current switch functionWith
In formula (15) and (16),For m-n phase of abc three-phase voltage switch functions
Amount;For its positive and negative sequence component;For n-m phasor of three-phase current switch function;For its positive and negative sequence component;M=0,1,2,3 ..., n=1,2,3 ..., a=ej2π/3。
It is preferred that, in step S5, the relational expression of DC voltage is tried to achieve by exchange change of current busbar voltage and voltage switch function
For:
In formula (17),For m phasor of DC voltage;Thus the voltage DC component of straight-flow system can be tried to achieveAnd second harmonic componentIt can divide with reference to direct current constraint equation (4), formula (5) in the hope of the fundamental frequency of DC side electric current
AmountAnd the m subharmonic currents of DC side can be tried to achieve by following formula:
Zd(m)The DC side equivalent harmonic wave impedance seen into for AC, wherein m=2 can be in the hope of the two of DC side electric current
Order harmonic components
It is preferred that, in step S6, the relational expression for trying to achieve alternating current by DC current and current switch function is:
For non-faulting side, n=1 then has:
With reference to non-faulting side exchange constraint equation (3) can in the hope of non-faulting side change of current bus power-frequency voltageIt is possible thereby to reference to the specific steps for setting up switch function, non-faulting side voltage x current switch function is corrected, and then
With reference to the fundamental component of calculation procedure amendment DC sideAnd second harmonic component
It is preferred that, in step S7, the electric current that DC side is injected into the AC network of failure side can be tried to achieve according to formula (19)
It is preferred that, in step S8, exchange can be tried to achieve according to formula (1), formula (2) and try to achieve failure side change of current busbar voltage
It is preferred that, in step S9, judgeWithDifference andWithDifference be
It is no to be satisfied by the condition of convergence, if being unsatisfactory forSubstitute into respectivelyUnder carrying out afterwards
Iterative calculation once, until meeting the condition of convergence or reaching that maximum iteration backed off after random is calculated, and exports result of calculation.
The present invention compared with prior art, has the following advantages that and beneficial effect:
1st, research object of the invention is the AC network containing back-to-back DC power transmission system, has been filled up containing back-to-back DC
The blank of research, is the stable operation of BTB systems and its adjusting for relay protection in terms of transmission system AC network calculation of fault
Cooperation is laid a good foundation.
2nd, present invention meter and the influence of failure top-cross direct current interphase interaction, at the same consider in side AC network therefore
By the AC-DC with back-to-back DC system-exchange interaction during barrier, one caused to non-faulting side power network is fixed
Ring.Push through journey by being passed back to and be iterated calculating to whole network, therefore the fault calculation methods for transmission can more realistic situation.
3rd, calculating speed of the present invention is fast, the degree of accuracy is high, in the AC network calculation of fault containing back-to-back DC power transmission system
And its there is very strong engineering practical value in terms of the setting value order of relay protection.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention;
Fig. 2 is back-to-back DC power transmission system pie graph of the present invention;
Fig. 3 (a) is the sequence diagrams of failure side AC network positive sequence after BTB systems AC network failure of the present invention;
Fig. 3 (b) is the sequence diagrams of failure top-cross stream electrical network negative phase-sequence after BTB systems AC network failure of the present invention;
Fig. 3 (c) is the sequence diagrams of failure side AC network zero sequence after BTB systems AC network failure of the present invention;
Fig. 3 (d) is the sequence diagrams of non-faulting side AC network positive sequence after BTB systems AC network failure of the present invention;
Fig. 3 (e) is the sequence diagrams of non-faulting top-cross stream electrical network negative phase-sequence after BTB systems AC network failure of the present invention;
Fig. 4 is straight-flow system control characteristic curve in the present embodiment.
Embodiment
With reference to embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention are not limited
In this.
Embodiment
As shown in figure 1, a kind of Power network fault calculation method containing back-to-back DC power transmission of the present invention, implements process such as
Under:
Step 1:Input the systematic parameter and control parameter of back-to-back transmission system.
Step 2:According to the structure and parameter of BTB system AC networks, for failure side AC network, with reference to AC system
Network topology structure and failure boundary condition can set up what is determined by AC systemWithAC system constraint equation;
For non-faulting side AC network, it can be set up with reference to the network structure of AC systemWithAC system constraint side
Journey;And the direct current constraint side set up respectively between its DC side DC voltage and DC current according to the DC control characteristic of both sides
Journey;Exchange constraint equation and direct current constraint equation specify that the function between AC voltage, electric current between DC voltage, electric current
Relation;Sequence network equation in the non-faulting top-cross stream constraint equation is as follows:
And failure boundary conditional equation is then as shown in formula (2):
Formula (1), (2) constitute failure top-cross stream constraint equation.In the formula of the above two,AndRespectively event
Hinder positive sequence, negative phase-sequence and the zero-sequence component of side change of current bus fundamental frequency voltages.For the electric current positive-sequence component of AC network equivalent power supply,Respectively DC side is injected into the Fundamental-frequency Current positive sequence and negative sequence component of AC network.WithRespectively
Positive-sequence component, negative sequence component and the zero-sequence component of fault current at failure.Can for the exchange constraint equation of non-faulting side
Following expression is obtained according to the topological structure and principle of stacking of its AC network:
In formula (3),ForAC network power frequency order impedance matrix during independent role,ForDuring independent role
AC network power frequency positive sequence admittance matrix.
For the direct current set up respectively between its DC side DC voltage and DC current according to the DC control characteristic of both sides
Constraint equation is expressed as follows:
Idi0=fi(Udi0) (4)
Idr0=fr(Udr0) (5)
In formula (4) and (5), Udi0、Idi0、Udr0、Idr0The DC component of rectification side DC voltage and electric current is expressed as,
The DC voltage of inverter side and the DC component of electric current, subscript " i " and " r " wherein in following table represent inverter side and rectification respectively
Side, is such as omitted, i.e. Ud0、Id0Then refer to the voltage DC component and current dc component of straight-flow system.
Step 3:K=0 is made, whole variables are put with initial value, including failure side and non-faulting side change of current busbar voltage
Step 4:Voltage x current is set up according to AC power-frequency voltage, DC side current dc component and initial trigger angle
Switch function model.More specifically, the calculating of the transverter voltage x current switch function power frequency order components is specific as follows:
A) according to the definition of Park Transformation, change of current bus can be converted to a, b, c three-phase voltage by positive and negative, residual voltage:
In formula (6),Respectively the fundamental frequency phase voltage of a, b, c phase of change of current bus, a=ej2π/3.Thus
Try to achieve corresponding fundamental frequency line voltageWithAnd then obtain the skew of synchronizing voltage phase
IfWithThe α components and β components of commutation voltage are represented respectively, are calculated by formula (7):
The phase of DC control system synchronizing voltage can be tried to achieve using the α components and β components of commutation voltage
IfPhase angle bePhase angle bePhase angle beIt can be counted by formula (9)
Calculate the phase offset of synchronizing voltage
In formula (9),For ca phases and the phase offset of synchronizing voltage,Phase for ab phases and synchronizing voltage is inclined
Move,For bc phases and the phase offset of synchronizing voltage;
B) according to AC power-frequency voltageWithDC side current dc component Id0And initial trigger angle
α0, calculate converter valve turn on delay angle θmn, actual Trigger Angle αmnWith actual angle of overlap μmn;
Turn on delay angle θmnComputing formula be:
In formula (10), xy=ab, bc, ca, wherein a, b, c represent the phase in three-phase respectively, similarly hereinafter;
Actual Trigger Angle αxyComputing formula be:
In formula (10) and (11), all angles are with delayed for just, advanced is negative;
Angle of overlap μ during xy two-phase commutationsxyComputing formula be:
μxy=cos-1(cosαxy-2XrId0/Uxy1)-αxy (12)
In formula (12), XrFor the commutating reactance of straight-flow system transverter;
C) according to θxy、αxyAnd μxyThree-phase voltage current switching waveform is made, Fu is utilized by the three-phase voltage current waveform
In leaf-size class number derive each order component of voltage x current switch function:
In formula (13) and (14),Respectively k order components of three-phase voltage switch function, Respectively k order components of three-phase current switch function, k=0,1,2,3 ..., T is 2 π;
Set up the order components of voltage x current switch functionWith
In formula (15) and (16),For m-n phase of abc three-phase voltage switch functions
Amount;For its positive and negative sequence component;For n-m phasor of three-phase current switch function;For its positive and negative sequence component;M=0,1,2,3 ..., n=1,2,3 ..., a=ej2π/3。
Step 5:The relational expression for trying to achieve DC voltage by exchange change of current busbar voltage and voltage switch function is:
In formula (17),For m phasor of DC voltage;Thus the voltage DC component of straight-flow system can be tried to achieveAnd second harmonic componentIt can divide with reference to direct current constraint equation (4), formula (5) in the hope of the fundamental frequency of DC side electric current
AmountAnd the m subharmonic currents of DC side can be tried to achieve by following formula:
Zd(m)The DC side equivalent harmonic wave impedance seen into for AC, makes the m=2 can be in the hope of the secondary of DC side electric current
Harmonic component
Step 6:The relational expression for trying to achieve alternating current by DC current and current switch function is:
For non-faulting side, n=1 is made, then is had:
With reference to non-faulting side exchange constraint equation (3) can in the hope of non-faulting side change of current bus power-frequency voltageIt is possible thereby to reference to the specific steps for setting up switch function, non-faulting side voltage x current switch function is corrected, and then
The fundamental component I of calculation procedure amendment DC sided(0)And second harmonic component
Step 7:The electric current that DC side is injected into the AC network of failure side can be tried to achieve according to formula (19)According to formula
(1), formula (2) can try to achieve exchange and try to achieve failure side change of current busbar voltage
Step 8:JudgeWithDifference andWithDifference whether be satisfied by receive
Hold back condition, if being unsatisfactory forSubstitute into respectivelyIteration next time is carried out afterwards
Calculate, until meeting the condition of convergence or reaching that maximum iteration backed off after random is calculated, and export result of calculation.
Fig. 2 is is back-to-back DC power transmission system pie graph in the present invention, back-to-back DC power transmission system as seen from the figure
Main element includes change of current bus, converter power transformer, rectifier and the inverter at two ends, while matching somebody with somebody on the ac bus at two ends
For more complete wave filter and reactive power compensator.Do not have transmission line of electricity between rectification side and inverter side, only pass through flat ripple
Reactor is connected.
The present embodiment establishes Lingbao BTB HVDC Engineering straight-flow system model based on PSCAD/EMTDC.Lingbao City
Back-to-back rated capacity is 360MW, and grade of rated voltage is 120kV, and rated current is 3000A, and smoothing reactor is
120mH, rectification side change of current busbar voltage is 330kV, and inverter side change of current busbar voltage is 220kV.For Fig. 2 BTB system knots
Composition, failure side and non-faulting side AC network sequence diagrams such as Fig. 3 (a) -3 are can obtain in a certain top-cross stream grid collapses
(e) show.In Fig. 3 (a) -3 (e),AndRespectively the positive sequence of failure side change of current bus fundamental frequency voltages, negative phase-sequence and
Zero-sequence component.For the electric current positive-sequence component of AC network equivalent power supply,Respectively DC side is injected into failure side
The Fundamental-frequency Current positive sequence and negative sequence component of AC network.WithThe positive-sequence component of fault current respectively at failure,
Negative sequence component and zero-sequence component.Respectively DC side be injected into failure side AC network Fundamental-frequency Current positive sequence and
Negative sequence component.In Fig. 3 (a) -3 (e), the parameter of BTB systems is:
Unit is kV; Unit Ω.
DC control system characteristic curve in the present embodiment is as shown in Figure 4.Can be with for the control function of different zones
Represented with piecewise function.Fine line A-I in figure is that rectifier runs on minimum trigger angle αminInversion side controller during control
Steady curve, heavy line A-Z is invertor operation in determining shut-off angle γ0The stable state fortune of rectification side controller during control
Row characteristic curve.Therefore, the constraint equation expression of the DC voltage of straight-flow system and electric current is as follows:
Rectification side:
The exchange constraint equation and direct current constraint equation of system can be set up according to above-mentioned parameter.K=0 is made, sets complete
Portion's variable initial value is steady-state operation value, i.e., DC voltage is Ud=1.0, DC current is Id=1.0 (perunit values).
Progressively calculated according to step 4~step 9, step-up error is no more than for front and rear change of current busbar voltage difference twice
0.001, the result of calculation of this computational methods is calculated by successive ignition.
Example one:It is set to rectification side AC network change of current bus and occurs singlephase earth fault, fault resstance is respectively 30
Ω、50Ω;Two-phase short-circuit fault, fault resstance is respectively 30 Ω, 50 Ω;Line to line fault earth fault, fault resstance is respectively
30 Ω, 70 Ω, calculate positive sequence, negative phase-sequence and the residual voltage of stable state change of current bus after failure, this algorithm result of calculation and PSCAD/
The simulation result of EMTDC simulation softwares is as shown in table 1.
Lingbao City's back-to-back DC model rectification side unbalanced fault change of current busbar voltage of table 1 is calculated and simulation result
Example two:It is set to inverter side AC network change of current bus and occurs singlephase earth fault, fault resstance is respectively 50
Ω、70Ω;Two-phase short-circuit fault, fault resstance is respectively 50 Ω, 70 Ω;Line to line fault earth fault, fault resstance is respectively
50 Ω, 70 Ω, calculate positive sequence, negative phase-sequence and the residual voltage of stable state change of current bus after failure, this algorithm result of calculation and PSCAD/
The simulation result of EMTDC simulation softwares is as shown in table 2.
Lingbao City's back-to-back DC model inverter side unbalanced fault change of current busbar voltage of table 2 is calculated and simulation result
Knowable to upper Tables 1 and 2, this algorithm is different in not homonymy failure, homonymy different type failure and same fault
The equal very little of error between result of calculation and simulation model simulation value under transition resistance, thus illustrates that this computational methods is calculated accurate
Degree is high, can be suitably used for various short troubles, and calculating speed is fast, is a kind of stronger electricity containing back-to-back DC power transmission of practicality
Net fault calculation methods for transmission.
Above-described embodiment is preferably embodiment, but embodiments of the present invention are not by above-described embodiment of the invention
Limitation, other any Spirit Essences without departing from the present invention and the change made under principle, modification, replacement, combine, simplification,
Equivalent substitute mode is should be, is included within protection scope of the present invention.
Claims (9)
1. a kind of Power network fault calculation method containing back-to-back DC power transmission, it is characterised in that comprise the steps:
S1, the back-to-back transmission system of input systematic parameter and control parameter;
S2, failure top-cross streaming system constraint side formed by the network topology structure and failure boundary condition of failure top-cross streaming system
Journey, is formed the AC system constraint equation of the side by the network topology structure of non-faulting top-cross streaming system, by the direct current control of both sides
The direct current constraint equation that characteristic processed is set up between its DC side DC voltage and DC current respectively;
S3, failure side and non-faulting side change of current busbar voltage areWherein k is the number of times of interative computation;To complete
Portion's variable puts initial value, i.e. failure side and change of current busbar voltage in non-faulting side is
S4, voltage x current switch function set up according to AC power-frequency voltage, DC side current dc component and initial trigger angle
Model, obtains the switch function of voltage, electric currentWith
S5, by failure side power frequency busbar voltageThe voltage DC of failure side straight-flow system is tried to achieve with reference to transverter switch function
ComponentAnd second harmonic componentSo as to try to achieve the fundamental component of DC side electric currentAnd second harmonic component
S6, with reference to transverter switch function and non-faulting top-cross stream constraint equation try to achieve non-faulting side power frequency busbar voltageAnd thus correct the fundamental component of DC sideAnd second harmonic component
S7, try to achieve the electric current that DC side is injected into the AC network of failure side
S8, the electric current in the AC network of injection failure sideThe exchange constraint equation of combination failure side can try to achieve exchange and try to achieve
Failure side change of current busbar voltage
S9, judgementWithDifference andWithDifference whether be satisfied by the condition of convergence,
If being unsatisfactory forSubstitute into respectivelyAnd repeat S4~S9 step, until convergence or
Reach that maximum iteration backed off after random is calculated, and export result of calculation.
2. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 1, it is characterised in that step
In rapid S2, the sequence network equation in the non-faulting top-cross stream constraint equation is as follows:
And failure boundary conditional equation is then as shown in formula (2):
Formula (1), (2) constitute failure top-cross stream constraint equation;In the formula of the above two,AndRespectively failure side
Positive sequence, negative phase-sequence and the zero-sequence component of change of current bus fundamental frequency voltages;For the electric current positive-sequence component of AC network equivalent power supply,Respectively DC side is injected into the Fundamental-frequency Current positive sequence and negative sequence component of AC network;WithRespectively
Positive-sequence component, negative sequence component and the zero-sequence component of fault current at failure;
It can be obtained for the exchange constraint equation of non-faulting side according to the topological structure and principle of stacking of its AC network such as following table
Up to formula:
In formula (3),ForAC network power frequency order impedance matrix during independent role,ForExchanged during independent role
Power network power frequency positive sequence admittance matrix;
For setting up the constraint of the direct current between its DC side DC voltage and DC current respectively according to the DC control characteristic of both sides
Equation is expressed as follows:
Idi0=fi(Udi0) (4)
Idr0=fr(Udr0) (5)
With in (4), (5), Udi0、Idi0、Udr0、Idr0It is expressed as the DC component of rectification side DC voltage and electric current, inversion
The DC voltage of side and the DC component of electric current, subscript " i " and " r " wherein in following table represent inverter side and rectification side respectively;
Ud0、Id0Then refer to the voltage DC component and current dc component of straight-flow system.
3. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 1, it is characterised in that step
In rapid S3, k is the number of times of interative computation, and when calculation error is unsatisfactory for requiring, k values Jia 1, continue cycling through calculating.
4. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 1, it is characterised in that step
In rapid S4, voltage x current switch function is set up according to AC power-frequency voltage, DC side current dc component and initial trigger angle
Model, the calculating of the transverter voltage x current switch function power frequency order components is specific as follows:
A) according to the definition of Park Transformation, change of current bus can be converted to a, b, c three-phase voltage by positive and negative, residual voltage:
In formula (6),Respectively the fundamental frequency phase voltage of a, b, c phase of change of current bus, a=ej2π/3;By phase voltage
And then try to achieve corresponding fundamental frequency line voltageWithAnd then obtain the skew of synchronizing voltage phase
WillWithThe α components and β components of commutation voltage are represented respectively, are calculated by following formula (7):
The phase of DC control system synchronizing voltage can be tried to achieve using the α components and β components of commutation voltage
WillPhase Angle Table be shown asPhase Angle Table be shown asPhase Angle Table be shown asBy formula (9)
The phase offset of synchronizing voltage can be calculated
In formula (9),For ca phases and the phase offset of synchronizing voltage,For ab phases and the phase offset of synchronizing voltage,For bc phases and the phase offset of synchronizing voltage;
B) according to AC power-frequency voltageWithDC side current dc component Id0And initial trigger angle α0, meter
Calculate converter valve turn on delay angle θmn, actual Trigger Angle αmnWith actual angle of overlap μmn;
Turn on delay angle θmnComputing formula be:
In formula (10), xy=ab, bc, ca, wherein a, b, c represent the phase in three-phase respectively, similarly hereinafter;
Actual Trigger Angle αxyComputing formula be:
In formula (10) and (11), all angles are with delayed for just, advanced is negative;
Angle of overlap μ during xy two-phase commutationsxyComputing formula be:
μxy=cos-1(cosαxy-2XrId0/Uxy1)-αxy (12)
In formula (12), XrFor the commutating reactance of straight-flow system transverter;
C) according to θxy、αxyAnd μxyThree-phase voltage current switching waveform is made, Fourier is utilized by the three-phase voltage current waveform
Series derives each order component of voltage x current switch function:
In formula (13) and (14),Respectively k order components of three-phase voltage switch function, Respectively k order components of three-phase current switch function, k=0,1,2,3 ..., T is 2 π;
Set up the order components of voltage x current switch functionWith
In formula,For m-n phasor of abc three-phase voltage switch functions;For it just,
Negative sequence component;For n-m phasor of three-phase current switch function;For its positive and negative sequence
Component;M=0,1,2,3 ..., n=1,2,3 ..., a=ej2π/3。
5. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 1, it is characterised in that
In step S5, the relational expression for trying to achieve DC voltage by exchange change of current busbar voltage and voltage switch function is:
In formula (17),For m phasor of DC voltage;Thus the voltage DC component of straight-flow system can be tried to achieveWith
Second harmonic componentWith reference to direct current constraint equation (4), formula (5) can in the hope of DC side electric current fundamental component
And the m subharmonic currents of DC side can be tried to achieve by following formula:
Zd(m)The DC side equivalent harmonic wave impedance seen into for AC, wherein m=2 can be in the hope of the secondary humorous of DC side electric current
Wave component
6. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 1, it is characterised in that
In step S6, the relational expression for trying to achieve alternating current by DC current and current switch function is:
For non-faulting side, n=1 then has:
With reference to non-faulting side exchange constraint equation (3) can in the hope of non-faulting side change of current bus power-frequency voltage
It is possible thereby to reference to the specific steps for setting up switch function, amendment non-faulting side voltage x current switch function, and then combine calculating
The fundamental component of step amendment DC sideAnd second harmonic component
7. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 6, it is characterised in that
In step S7, the electric current that DC side is injected into the AC network of failure side can be tried to achieve according to formula (19)
8. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 2, it is characterised in that
In step S8, exchange can be tried to achieve according to formula (1), formula (2) and try to achieve failure side change of current busbar voltage
9. a kind of Power network fault calculation method containing back-to-back DC power transmission according to claim 1, it is characterised in that
In step S9, judgeWithDifference andWithDifference whether be satisfied by restrain bar
Part, if being unsatisfactory forSubstitute into respectivelyIterative calculation next time is carried out afterwards,
Until meeting the condition of convergence or reaching that maximum iteration backed off after random is calculated, and export result of calculation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710285879.2A CN106953348A (en) | 2017-04-27 | 2017-04-27 | A kind of Power network fault calculation method containing back-to-back DC power transmission |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710285879.2A CN106953348A (en) | 2017-04-27 | 2017-04-27 | A kind of Power network fault calculation method containing back-to-back DC power transmission |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106953348A true CN106953348A (en) | 2017-07-14 |
Family
ID=59477792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710285879.2A Pending CN106953348A (en) | 2017-04-27 | 2017-04-27 | A kind of Power network fault calculation method containing back-to-back DC power transmission |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106953348A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110011346A (en) * | 2019-03-19 | 2019-07-12 | 国家电网有限公司 | A kind of interactional commutation failure methods of risk assessment of consideration Inverter Station |
CN110441653A (en) * | 2019-07-26 | 2019-11-12 | 南京工程学院 | The feature extraction method of discrimination of the exchange side fault type of alternating current-direct current mixing power grid |
CN110569576A (en) * | 2019-08-22 | 2019-12-13 | 天津大学 | alternating current system fault steady-state analysis method containing direct current feed-in |
CN112994072A (en) * | 2021-04-27 | 2021-06-18 | 华北电力大学 | Two-port equivalent modeling method suitable for back-to-back VSC fault transient calculation |
CN114374208A (en) * | 2022-01-13 | 2022-04-19 | 广东电网有限责任公司 | Subsynchronous oscillation prejudging method and subsynchronous oscillation prejudging device for direct-drive fan |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140009980A1 (en) * | 2011-11-11 | 2014-01-09 | Varentec, Inc. | Power flow controller with a fractionally rated back-to-back converter |
CN104617576A (en) * | 2015-01-29 | 2015-05-13 | 国网宁夏电力公司 | Multi-DC-feed (multi-direct current-feed) AC (alternate current) grid fault calculating method taking DC control characteristics into consideration |
-
2017
- 2017-04-27 CN CN201710285879.2A patent/CN106953348A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140009980A1 (en) * | 2011-11-11 | 2014-01-09 | Varentec, Inc. | Power flow controller with a fractionally rated back-to-back converter |
CN104617576A (en) * | 2015-01-29 | 2015-05-13 | 国网宁夏电力公司 | Multi-DC-feed (multi-direct current-feed) AC (alternate current) grid fault calculating method taking DC control characteristics into consideration |
Non-Patent Citations (2)
Title |
---|
刘俊磊等: ""计及直流控制特性的直流系统等值模型及其谐波计算"", 《电力系统自动化》 * |
麦国浩等: ""背靠背直流输电的交流电网故障计算模型研究"", 《广东电力》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110011346A (en) * | 2019-03-19 | 2019-07-12 | 国家电网有限公司 | A kind of interactional commutation failure methods of risk assessment of consideration Inverter Station |
CN110441653A (en) * | 2019-07-26 | 2019-11-12 | 南京工程学院 | The feature extraction method of discrimination of the exchange side fault type of alternating current-direct current mixing power grid |
CN110569576A (en) * | 2019-08-22 | 2019-12-13 | 天津大学 | alternating current system fault steady-state analysis method containing direct current feed-in |
CN110569576B (en) * | 2019-08-22 | 2023-04-18 | 天津大学 | Alternating current system fault steady-state analysis method containing direct current feed-in |
CN112994072A (en) * | 2021-04-27 | 2021-06-18 | 华北电力大学 | Two-port equivalent modeling method suitable for back-to-back VSC fault transient calculation |
CN112994072B (en) * | 2021-04-27 | 2022-08-19 | 华北电力大学 | Two-port equivalent modeling method suitable for back-to-back VSC fault transient calculation |
CN114374208A (en) * | 2022-01-13 | 2022-04-19 | 广东电网有限责任公司 | Subsynchronous oscillation prejudging method and subsynchronous oscillation prejudging device for direct-drive fan |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106953348A (en) | A kind of Power network fault calculation method containing back-to-back DC power transmission | |
CN104811067B (en) | PR (proportional resonant) controller-based NMC-HVDC (modular multilevel converter-high voltage direct current) circulating current suppression method | |
CN104617576B (en) | Multi-DC-feed (multi-direct current-feed) AC (alternate current) grid fault calculating method taking DC control characteristics into consideration | |
CN101572408B (en) | Method for improving simulation capacity of HVDC transmission system | |
CN107171584A (en) | The model predictive control method and device of NPC three-phase tri-level combining inverters | |
CN106340907A (en) | Power system security and stability control strategy determining method and device | |
CN104993494B (en) | Motor simulator based on four-quadrant power electronic converter and method | |
CN109256970B (en) | MMC-MTDC transmission system monopolar grounding fault current calculation method | |
CN104866689B (en) | A kind of unbalanced fault straight-flow system quasi steady state model construction method | |
CN104333030B (en) | A kind of analysis method of multi-infeed HVDC interaction factor | |
CN105958485A (en) | Power flow calculation method for flexible interconnecting alternating current-direct current hybrid power distribution network | |
CN106066918A (en) | Based on homotopy containing distributed power source and the short-circuit current calculation method of nonlinear-load | |
CN102361329A (en) | Modeling method for performing dynamic characteristic research on hybrid alternating current/direct current (AC/DC) transmission system | |
Sybille et al. | Simulation of FACTS controllers using the MATLAB power system blockset and hypersim real-time simulator | |
CN102508098A (en) | Quick distinguishing method for zero-sequence directional elements of alternating-current and direct-current serial-parallel power grid | |
CN104407195B (en) | A kind of unbalanced fault short-circuit current detection method of alternating current-direct current combined hybrid system | |
CN103793755A (en) | Method for obtaining optimal power flow of mixed direct-current transmission system based on prime-dual interior point method | |
CN105095590B (en) | A kind of modeling method of the electromechanical transient simulation system based on three sequence equivalent impedances | |
CN110556853A (en) | calculation method and system for providing initial value for electromagnetic transient simulation | |
CN102801152B (en) | Electromechanical transient model for cascade connection type multi-terminal direct current power transmission system and inter-station coordination method | |
CN103715704A (en) | Micro electrical network common bus voltage imbalance inhibition method | |
CN105514972B (en) | The PSCAD modelings of grid-connected converter and emulation mode during unbalanced grid faults | |
Huang et al. | A calculation method for maximum steady-state loss of VSC-HVDC grounding resistor in Luxi back-back HVDC interconnector | |
CN103560494B (en) | The short circuit current acquisition methods that distribution protection is adjusted | |
CN102593876A (en) | Continuous power flow algorithm of large-scale alternating current and direct current interconnected electric power system |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170714 |
|
WD01 | Invention patent application deemed withdrawn after publication |