This invention belongs to a high voltage direct current transmission technology, in particular to a fault transient features optimization method and system of AC/DC system based on dissipated energy.
When the DC Line Protection detects a fault on the DC line, the rectifier station will enter the states of phase shift triggering and deionization, to wait for the fault to be cleared. Then, the DC restart process will be carried out by restarting the system. However, if the fault is not cleared during the process, the high DC pulsed current generated by restarting at the fault will cause a secondary impact to the system, which will cause the transient power angle of the rectifier side AC generators to swing greatly. In severe case, it may even lead to system instability.
Currently, the solutions to the above problems are mainly divided into two types:
The first type of measures uses DC current limiters or circuit breakers with current limiting capability to limit the restart DC pulsed current, thus reducing the strength of impact. A current limiter with power management capability was put forward, which can limit the fault current and control part of the power flows in AC/DC system. The superconducting fault current limiter exhibits zero resistance under steady-state operation and high resistance under fault pulsed current. Thus, the superconducting fault current limiter can limit the fault current spontaneously, which requires no need for the complicated cooperation between the control method and the relay protection. Circuit breakers with fault current limitation capability can cut off the fault current and limit the restart DC pulsed current. However, this method is difficult to be directly applied to power system due to its high cost and immature technology.
The second type of measures reduces the risk of instability during the fault process by adding control branches or regulating the control parameters to improve the transient stability of the system, and reduce the instability risk brought to the system during the fault. A nonlinear coordinated control strategy for multiple HVDC links was proposed, which can improve both the stability margin and power transfer capability of the system. Based on the Lyapunov's energy function, the operating features of the AC/DC system can be improved, and partially improved Lyapunov's functions can coordinate the power modulation of the MTDC or multiple point-to-point DC systems, which increases the domain of attraction and improves the stability margin of the generator power angle. The disadvantage of this method is that it is difficult to improve the transient features of the system greatly during restart.
In view of the above analysis, the embodiment of the invention provides a fault transient features optimization method and system of AC/DC system based on dissipated energy, which can be used to optimize the transient features of the system in the process of DC restart failure, so as to solve the problem that the pulsed current in DC system is too high and the rectifier side AC system generator power angle swing is too large in the process of HVDC system restart on fault.
On the one hand, the invention discloses a fault transient features optimization method of AC/DC system based on dissipated energy, the method including the following steps:
Obtaining the output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC system under the current set of control parameters when every time adjusting the control parameters of the rectifier, and obtain multiple sets of control parameters and their corresponding output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC system.
For each set of control parameters, calculating cumulative value of dynamic dissipated energy of the generator under the set of control parameters based on the output current and outlet bus voltage of the generator during the whole process of DC system fault recovery under the set of control parameters.
Selecting the set of control parameters that minimizes the cumulative value of dynamic dissipated energy of the generator as optimal control parameters in the process of fault transient features optimization.
On the basis of the above scheme, the invention makes the following improvements:
Further, the cumulative value of dynamic dissipated energy of the generator under the set of control parameters is obtained by performing the following operations:
Obtaining every commutation interval of Y/D bridge of inverter in the whole process of AC/DC system fault recovery under the set of control parameters.
Separately calculating the dynamic dissipated energy of the generator in each commutation interval.
Accumulating the dynamic dissipated energy of the generator in all commutation intervals and obtaining the cumulative value of dynamic dissipated energy of the generator under the set of control parameters.
Further, the rectifier adopts a constant current control mode.
Further, the cumulative value of dynamic dissipated energy of the generator WΣ is:
WΣ=ΣWG(n)(Kr, Tr, tP2Y, tP2Δ) (1)
Where, WG(n) represents the dynamic dissipated energy of the generator in the nth commutation interval; Kr, Tr is the control parameter of the rectifier in constant current control mode, where Kr is proportional coefficient of PI unit, Tr is the time constant of PI unit; tP2Y and tP2Δ respectively represent start time and end time of the commutation interval of Y/D bridge of inverter.
Further, the dynamic dissipated energy of the generator WG is:
W
G
=∫i
dG
du
qG
−∫i
qG
du
dG+∫(udGidG+uqGiqG)dθ (2)
Where idG and iqG respectively represent components of the output current of the generator in axis d and axis q; udG and uqG respectively represent components of generator outlet bus voltage on axis d and q, and θ is the included angle between axis q and axis x.
On the other hand, the invention also discloses a fault transient features optimization system of AC/DC system based on dissipated energy, including:
A data acquisition module is configured to obtain the output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC system under the current set of control parameters when adjusting the control parameters of the rectifier, and obtain multiple sets of control parameters and their corresponding output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC systems;
A cumulative value of dynamic dissipated energy acquisition module is configured to, for each set of control parameters, calculate the cumulative value of dynamic dissipated energy of the generator under the set of control parameters based on the output current and outlet bus voltage of the generator during the whole process of DC system fault recovery under the set of control parameters.
An optimized control parameter acquisition module is configured to select the set of control parameters that minimizes the cumulative value of dynamic dissipated energy of the generator as optimal control parameters in the process of fault transient features optimization.
On the basis of the above scheme, the invention also makes the following improvements:
Further, the cumulative value of dynamic dissipated energy of the generator under the set of control parameters is obtained by performing the following operations:
Obtaining every commutation interval of Y/D bridge in the whole process of AC/DC system fault recovery under the set of control parameters.
Separately calculating the dynamic dissipated energy of the generator in each commutation interval.
Accumulating the dynamic dissipated energy of the generator in all commutation intervals and obtaining the cumulative value of dynamic dissipated energy of the generator under the set of control parameters.
Further, the rectifier adopts a constant current control mode.
Further, the cumulative value of dynamic dissipated energy of the generator WΣ is:
WΣ=ΣWG(n)(Kr, Tr, tP2Y, tP2Δ) (3)
Where, WG(n) represents the dynamic dissipated energy of the generator in the nth commutation interval; Kr, Tr is the control parameter of the rectifier in constant current control mode, where Kr is proportional coefficient of PI unit, Tr is the time constant of PI unit; tP2Y and tP2Δ respectively represent start time and end time of the commutation interval of Y/D bridge of inverter.
Further, the dynamic dissipated energy of the generator WG is:
W
G
=∫i
dG
du
qG
−∫i
qG
du
dG+∫(udGiDG+uqGiqG)dθ (4)
Where, idG and iqG respectively represent components of the output current of the generator in axis d and axis q; udG and uqG respectively represent the components of generator outlet bus voltage on axis d and q, and θ is the included angle between axis q and axis x.
The fault transient features optimization method and system of AC/DC system based on dissipated energy provided by the invention have the following advantages:
The topology change of the rectifier during the restart process is considered in detail, which can realize the analysis of the impact of the converter conduction topology change on the system during the restart failure process, and realize the optimization of the transient features of the AC/DC hybrid system during the restart process.
It well solves the problems that the pulsed current of the DC system is too high and the power angle swing of the AC system generator on the rectifier side is too large during the restart failure of the HVDC transmission system.
In the invention, the above technical schemes can also be combined with each other to realize more optimal combination schemes. Other features and advantages of the invention will be described in a subsequent specification, and some of the advantages may become apparent from the specification or be known by the implementation of the invention. The purpose and other advantages of the invention can be realized and obtained from the contents specially indicated in the specification and appended drawings.
The drawings are only for the purpose of showing specific embodiments and are not considered to be a limitation of the disclosure. Throughout the drawings, the same reference symbols represent the same components.
The preferred embodiments of the present disclosure are described in detail with reference to the attached drawings. The accompanying drawings form a part of the present application and are used together with embodiments of the present disclosure to explain the principles of the disclosure and are not intended to limit the scope of the disclosure.
A specific embodiment of the invention discloses a fault transient features optimization method of AC/DC system based on dissipated energy. The flow chart is shown in
Step S1: Obtaining the output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC system under the current set of control parameters when every time adjusting the control parameters of the rectifier, and obtain multiple sets of control parameters and their corresponding output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC system.
Step S2: for each set of control parameters, calculating the cumulative value of dynamic dissipated energy of the generator under the set of control parameters based on the output current and outlet bus voltage of the generator during the whole process of DC system fault recovery under the set of control parameters.
Step S3: selecting the set of control parameters that minimizes the cumulative value of dynamic dissipated energy of the generator as optimal control parameters in the process of fault transient features optimization.
Compared with the prior art, the fault transient features optimization method of AC/DC system based on dissipated energy provided in this embodiment, the topology change of the rectifier in the restart process is considered in detail to realize the analysis of the impact of the converter conduction topology change on the system in the restart failure process, and the optimization of the transient features of the AC/DC hybrid system in the restart process. The problem that the DC system pulsed current is too high and the power angle swing of the AC system generator on the rectifier side is too large in the process of the restart failure of the HVDC transmission system is well solved.
Preferably, in this embodiment, the cumulative value of dynamic dissipated energy of the generator under the set of control parameter is obtained by performing the following operations:
Obtaining every commutation interval of Y/D bridge of inverter in the whole process of AC/DC system fault recovery under the set of control parameters.
Separately calculating dynamic dissipated energy of the generator in each commutation interval.
Accumulating the dynamic dissipated energy of the generator in all commutation intervals and obtaining the cumulative value of dynamic dissipated energy of the generator under the set of control parameters.
Specifically, the rectifier adopts a constant current control mode. At this point, the cumulative value of dynamic dissipated energy of the generator WΣ is:
WΣ=ΣWG(n)(Kr, Tr, tP2Y, tP2Δ)
Where, WG(n) represents the dynamic dissipated energy of the generator in the nth commutation interval; Kr, Tr is the control parameter of the rectifier in constant current control mode, where Kr is proportional coefficient of PI unit, Tr is time constant of PI unit; tP2Y and tP2Δ respectively represent start time and end time of the commutation interval of Y/D bridge of inverter.
The dynamic dissipated energy of the generator WG is:
W
G
=∫i
dG
du
qG
−∫i
qG
du
dG+∫(udGidG+uqGiqG)dθ
Where, idG and iqG respectively represent components of the output current of the generator in axis d and axis q; udG and uqG respectively represent components of generator outlet bus voltage in axis d and axis q.
In order to facilitate the person skilled in the art to better understand the implementation process of this embodiment, the scheme of this embodiment is introduced as follows:
Step 1. Considering the influence of phase shift triggering, deionization and restart stages in the whole transient process of the fault recovery stage on the generator, the dissipated energy model of the generator injection system including the rectifier control parameters is established.
In this embodiment, considering the influence of phase shift triggering, deionization and restart stages on the generator during the whole transient process of fault recovery stage, the dissipated energy model of generator injection system including rectifier control parameters is established.
The diagram of the rectifier side of AC/DC system (“AC/DC hybrid system”) is shown in
The wiring diagram of rectifier side converter is shown in
In DC system, the whole fault recovery process contains three stages: ‘phase shift triggering’, ‘deionization’ and ‘restart’. When a DC line fault is detected, the rectifier side will enter the ‘phase shift triggering’ stage first. In this stage, the control system will adjust the rectifier-side firing angle to 120°, so that the rectifier will turn to the inverter-operation-state and the fault current (idc) will quickly drop to 0. After the phase shifting, the firing angle of the rectifier station will be adjusted to 160°, so that the rectifier valves will stop commutating and the DC system will enter the ‘deionization’ stage to wait for the fault arc to be extinguished. After the deionization stage is completed, the control system will restore the real-time firing angle of the rectifier side to the initial angle (17° in the example system) and try to restart. If the fault is not cleared, restarting on the fault will lead to a high restart DC pulsed current, which will cause the DC protection to operate again, and then the above fault recovery process will repeat again. Since the switching of the rectifier operation mode in the above process, the traditional quasi-steady-state mean value model cannot depict the dynamic fault recovery process in details. This embodiment uses the detailed model of the rectifier shown in
During the DC phase shifting and restart process, the rectifier can commutate and conduct according to the normal trigger sequence; during the deionization process, the commutator valve maintains its conducting state and stops commutation. Thus, for the AC line and constant impedance load where the inductor L1 is located, there are:
Concerning the DC line, there are:
Concerning the AC filter, there are:
Concerning the 12-pulse rectifier, suppose Y bridge VT1, VT2 and D bridge VT1, VT6 are turned on:
Concerning the rectifier-side converter bus in the ‘restart’ stage:
Concerning the generator outlet bus:
In power system, the node current equation can be used to construct the dissipated energy function. The equivalent node current equation of the AC/DC system can be expressed as:
I
in
=YU
B
−I
G
+I
L (7)
where Y is the admittance matrix of the system; UB is the column vector of bus voltages; IG and IL are the column vectors of injection currents at generator nodes and load nodes. Integrating the imaginary part i*ij of the injection currents at generator nodes IG yields the dynamic dissipated energy of all generator injected to system:
where iij the current on the branch Lij, which is connected with generator node i; ui is the voltage at node i; subscripts x and y represent the variables in xy coordinate system. The relationship between the xy coordinate system and the dq coordinate system is shown in
According to
Similarly:
u
i=(ud sin θ+uq cos θ)+j(uq sin θ−ud cos θ) (10)
In equations (9) and (10), subscripts d and q represent the variables in dq coordinate system. By applying equations (9) and (10) to equation (8), the expression of the dynamic dissipated energy injected from the generator to the system can be obtained:
W=∫i
d
du
q
−∫i
q
du
d+∫(udid+uqiqdθ (11)
The AC/DC system rectifier-side system includes the generator, AC lines, loads, DC terminal and AC filters, thus according to (7) and the conservation law of energy:
W
G
=∫i
dG
du
qG
−∫i
qG
du
dG+∫(udGidG+uqGqGdθ=WHVDC+WLn+WL+WC (12)
Where, WHVDC, WLn,WC and WG are the dynamic dissipated energy of DC terminals, AC lines, loads, AC filters, and generator under AC/DC systems.
By combining equations (1)-(6), the voltage and current at the generator terminal can be calculated. Apply Park transformation to the voltage and current and then apply the results to equation (12), so that the dynamic dissipated energy injected from the generator to the system under the above conducting topology of the rectifier can be obtained:
Where, Um is the amplitude of generator outlet bus voltage; λ is the increasing coefficient of the voltage amplitude (generator actual outlet bus voltage divided by Um); ω is the rated speed of the generator; ωi is the actual speed of the generator; C1, C2 and C3 are constants generated in the integration operation. The specific expressions of other coefficients are as follows:
Where, C is the equivalent capacitance of the AC filter, and Rc is the equivalent resistance of the AC filter.
Due to the symmetry of the rectifier's conduction topology, the expression of the generator's dynamic dissipated energy under different conduction states is consistent with the form of equation (13).
Step 2: calculate the cumulative value of dynamic dissipated energy during the whole process of AC/DC system fault recovery.
Taking the Cigre HVDC transmission model as an example. Its rectifier adopts constant current control, and its control block diagram is shown in
The control equation of the rectifier is:
Where, αrec is the firing angle of the rectifier; Id_ref is the reference value of rectifier-side DC current; idc is the measured value of DC current; Kr and Tr are the PI parameters of rectifier-side constant current controller.
The intersection of phase voltages is taken as the absolute zero point, and taking the trigger pulse of TV2 in one cycle as the reference zero point, which is shown in
In normal operation state, the firing angles at tP2Y and tP2Δ are the same, i.e. α(tP2Y)=α(tP2Δ). Thus, the angle difference between P2Y and P2Δ is 30°, and the 1st commutation interval in one cycle period is
(with P2Y as the reference zero point). When a fault occurs on the DC line, affected by the constant current control, α(tP2Y)≠α(tP2Δ). Thus, the angle difference between P2Y and P2Δ will be 30°+αrec (tP2Δ)−αrec(tP2Y), and the 1st turn-on interval will become
According to equation (14), the rectifier firing angle αrec is a function of t. Hence, the expressions of αrec(t) in different turn-on interval of the converter valves are different. When VT1 and VT2 of Y bridge and VT1 and VT6 of D bridge are turned on, the turn-on interval is calculated as follows. Suppose this interval starts from tP2Y and ends at tP2Δ, thus, the firing angles at tP2Y and tP2Δ are:
The time length corresponding to this turn-on interval is:
Equation (17) is denoted as φ(Kr, Tr), a function of Kr and Tr.
In equation (17), the expressions of idc(t), tP2Y and tP2Δ are to be determined. By combining equations (1)-(6), the expression of DC current idc in this commutation interval can be obtained:
Where C4 is a constant generated during the integration operation. The relationship between the turn-on moments (tP2Y and tP2Δ) and the corresponding firing angles (αrec(tP2Y) and αrec(tP2Δ)) is:
According to equations (14), (15) and (18)-(20), tP2Y, tP2Δ and rectifier control parameters Kr and Tr form a transcendental equation containing implicit functions. Combining equation (17), the cumulative value of dissipated energy of the generator in the whole process of fault recovery' can be expressed as:
WΣ=ΣWG(n)(Kr,Tr,tP2Y,tP2Δ) (21)
Where-n represents the nth turn-on interval.
Step 3: To minimize the cumulative value of dissipated energy, solve the optimal control parameters of the rectifier and output them.
In order to reduce the impact of the dissipated energy on the sending-end system, WΣ should be made as small as possible. Denote the integrand function of WG in the nth turn-on interval as fn(t), and the upper and lower limits of the integration in the nth turn-on interval as φn(Kr, Tr) and φn+1(Kr, Tr). By combining equations (16) and (21), the partial derivative expressions of the cumulative value WΣ of the total system energy during the fault recovery period to Kr and Tr can be calculated:
According to equations (22) and (23), the conditions under which WΣ reaches the extreme values are:
When VT1 and VT2 of Y bridge and VT1 and VT6 of D bridge are turned on, fn[φn+1(Kr,Tr)] can be further expanded as
Equation (26) is the multivariate function of variables Kr, Tr, tP2Y and tP2Δ.
Further, combining equation (16), the partial derivative of φ(Kr, Tr) with respect to parameter Kr can be calculated:
By combining equations (14), (15), (19) and (20), the implicit functions of tP2Y, tP2Δ, Kr and Tr can be obtained:
Denote the implicit function in equation (28) as F(Kr, Tr, tP2Y, tP2Δ) and the implicit function in equation (29) as G(Kr, Tr, tP2Y, tP2Δ). According to the implicit function theorem of multi-variable functions:
As can be seen from equations (26)-(31), control parameters Kr and Tr are transcendental equations of the implicit functions of tP2Δ and tP2Y. By calculating equations (24) and (25), the system of the solutions of Kr and Tr can be obtained:
(Kr(i), Tr(i))={(Kr(1), Tr(1)),(Kr(2),Tr(2)), . . . ,(Kr(0),Tr(0))} (32)
By applying the elements in equation (32) to equation (21) and screening out the local minimum values of WΣ, the solution system of local minimum values WimΣ can be established:
WimΣ={W1mΣ|(K
After the control parameters are optimized, the cumulative value of dissipated energy of the system WiΣ should satisfy:
WiΣ<W0Σ (34)
Where, W0Σ is the accumulated energy of the system without parameter optimization. By combining equations (33) and (34), the system of solutions that satisfy the optimization conditions can be further screened:
WiΣ={W1Σ|(K
Then, the values of Kr and Tr corresponding to the local minimum values in WiΣ (denoted as WiΣmin) are selected as the optimal control parameters, which can minimize the impact exerted on the system during fault recovery.
Embodiment 2 of the invention discloses a fault transient features optimization system of AC/DC system based on dissipated energy. As shown in
A data acquisition module is configured to obtain the output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC system under the current set of control parameters when every time adjusting the control parameters of the rectifier, and obtain multiple sets of control parameters and their corresponding output current and outlet bus voltage of the generator during the whole process of fault recovery of the AC/DC systems;
A cumulative value of dynamic dissipated energy acquisition module is used to, for each set of control parameters, calculate cumulative the value of dynamic dissipated energy of the generator under the set of control parameters based on the output current and outlet bus voltage of the generator during the whole process of dc system fault recovery under the set of control parameters.
An optimized control parameter acquisition module is used to select the set of control parameters that minimizes the cumulative value of dynamic dissipated energy of the generator as optimal control parameters in the process of fault transient features optimization.
The specific implementation process of the second embodiment of the invention can be referred to the above embodiment of the method, and this embodiment is not repeated here.
As the principle of this embodiment is the same as that of the above embodiment of the method, the system also has the corresponding technical effect of the above embodiment of the method.
Take Cigre_Benchmark AC/DC transmission system as shown in
After the fault occurs, the system enters the first phase shilling stage at t=5.02 s. After the phase shifting, the system enters the deionization stage at 5.07 s. After deionization is completed, the system tries to restart for the first time at t=5.32 s. Because the fault is not cleared at this time, the system enters the second “phase shifting-deionization” at t=5.395 s after the restart failure, and tries the second restart at t=5.715 s. However, the fault is still not cleared at this time, so the system enters the third “phase shifting-deionization” and tries to restart again at t=6.16 s. The fault is rectified and the system restarts successfully for the third time.
Calculate the cumulative value of dynamic dissipated energy of the system. The system tries to restart for the first time at t=5.32 s. At this time, the restart fails because the line fault is not cleared. During this period, the restart on fault will generate high pulsed current in the DC system; and the higher the amplitude of the pulsed current the higher the threat to system stability and device security. By substituting the control parameters in Table 1 into the system, the DC current waveform under different parameters when the fault restarts can be obtained, as shown in
As can be seen from
According to table 1, when the control parameters adopt the standard values, the peak value of the restart DC pulsed current reaches as high as 1.57556 p.u. When the control parameters are adjusted to the optimal values (Kr=1.1074, Tr=8.31070), the peak value of the restart DC pulsed current decreases significantly (by 39.59%) to 0.95185 p.u.
At the same time, the transient power angle of the generator will also swing due to the high restart DC pulsed current. The larger the peak-valley difference in the swing curve of the power angle is, the more seriously the generator is disturbed, and the worse it is for the safe and stable operation of the system. After regulating the control parameters, the variation curves of the generator power angle during the 1st ‘restart at fault’ are shown in
According to
According to Table 2, when the control parameters are not regulated, the generator power angle swings violently, with the peak-valley difference reaching as high as 0.1180 p.u. When the control parameters are adjusted to the optimal values (Kr=1.1074, Tr=8.31070), the swing of generator power angle obviously softens. The peak-valley difference of the power angle decreases (by 75.08%) to 0.0294 p.u.
After the 1st restart fails, the system tries the 2nd restart at t=5.715 s after the 2nd ‘phase shifting-deionization’ process. The 2nd restart still fails because the fault is still not cleared. During the 2nd restart failure, the waveforms of the restart DC pulsed current under different parameter settings are shown in
As can be seen from
During the 2nd ‘restart at fault’, the variation curves of the generator power angle under different control parameters are shown in
According to
The person skilled in the art can understand that all or part of the process of realizing the above embodiments can be completed by instructing the relevant hardware through a computer program, and the program can be stored in a computer readable storage medium. Wherein, the computer readable storage medium for disk, optical disk, read-only storage memory or random storage memory.
The above is only a better specific implementation of the invention, but the protection scope of the invention is not limited to this, any change or replacement that can be easily thought of by a technical personnel familiar with the technical field within the scope of the technology disclosed by the invention shall be covered by the protection scope of the invention.
Number | Date | Country | Kind |
---|---|---|---|
202210190459.7 | Feb 2022 | CN | national |