CN110083863B - Static reactive compensator-based electromagnetic transient rapid simulation method - Google Patents

Static reactive compensator-based electromagnetic transient rapid simulation method Download PDF

Info

Publication number
CN110083863B
CN110083863B CN201910208307.3A CN201910208307A CN110083863B CN 110083863 B CN110083863 B CN 110083863B CN 201910208307 A CN201910208307 A CN 201910208307A CN 110083863 B CN110083863 B CN 110083863B
Authority
CN
China
Prior art keywords
thyristor
state
phase
svc
phases
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.)
Expired - Fee Related
Application number
CN201910208307.3A
Other languages
Chinese (zh)
Other versions
CN110083863A (en
Inventor
麦立昀
毛晓明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201910208307.3A priority Critical patent/CN110083863B/en
Publication of CN110083863A publication Critical patent/CN110083863A/en
Application granted granted Critical
Publication of CN110083863B publication Critical patent/CN110083863B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

The invention discloses an electromagnetic transient rapid simulation method based on a static var compensator, which comprises the following steps: acquiring basic parameters of an SVC device; the parameters to be acquired include: inductance L of TCR branch circuit, capacitance C of fixed capacitor FCsResistance R of filterfInductor LfAnd a capacitor CfExternal input voltage uA、uB、uCAnd a thyristor firing angle α; the method introduces a sine auxiliary variable, converts the state equation of the SVC from a non-homogeneous linear differential equation set into a homogeneous linear differential equation set, quickly obtains the solution of the equation set by a matrix exponential integration method, and meanwhile, the state equations under various working conditions are unified in form; compared with the traditional electromagnetic transient calculation in which the solution required system model changes at each time step, the method provided by the invention has the advantages that the solution required system model changes only when the on-off state of the thyristor changes, so that the calculation of the method provided by the invention is more efficient, and the working efficiency is greatly improved.

Description

Static reactive compensator-based electromagnetic transient rapid simulation method
Technical Field
The invention relates to the technical field of electrical simulation, in particular to an electromagnetic transient rapid simulation method based on a static var compensator.
Background
Static Var compensator svc (static Var compensator) is a power electronic device, and is usually composed of thyristor Controlled reactor tcr (thyristor Controlled reactor) and fixed capacitor fc (fixed capacitor). When the electric power system has a fault, the reactive power output by the SVC to the electric power system can be quickly and smoothly adjusted by quickly adjusting the trigger angle of the TCR in the SVC device, so that the purpose of adjusting the voltage of the electric power system is achieved. With the improvement of the performance of power electronic devices and the reduction of the manufacturing cost, the application of the SVC in power systems is more and more common.
Time domain simulation is an important means for researching the operating characteristics of the power system, and electromechanical transient simulation and electromagnetic transient simulation are two important contents of power system off-line simulation. For power electronic equipment, electromechanical transient simulation ignores the dynamic process inside the equipment, a quasi-steady-state modeling method is adopted, only the operation and control characteristics of the equipment under the positive sequence fundamental frequency are considered, and the method has great limitation on researching the influence of the asymmetric fault of the power grid on the equipment and the influence of the control method inside the equipment on the operation of a large power grid; the electromagnetic transient simulation carries out three-phase modeling on the power electronic equipment on an accurate circuit layer, has the advantages of accurate calculation, high reliability and the like, but has small calculation step length and low calculation speed.
In order to find a calculation method which is suitable for analyzing the electromechanical transient process of a large power grid and can fully consider the quick response of power electronic equipment, the electromechanical-electromagnetic transient hybrid simulation technology of a power system is developed. The hybrid simulation is used for performing electromagnetic transient simulation on the power electronic equipment by using a detailed device model, and performing electromechanical transient simulation on an external alternating current power grid, so that the rapid transient process of a specific device is reflected while the dynamic response of the large power grid is accurately obtained. When the infiltration level of power electronic equipment such as SVC is not high, the traditional hybrid simulation technology provides an effective means for researching the interaction between a large power grid and the power electronic equipment. With the scale enlargement of the power grid, more and more high-power electronic devices are connected into the power system. For the safety analysis of a large power grid, on the premise of ensuring the precision, it is necessary to develop a new method for electromagnetic transient modeling simulation of power electronic equipment such as SVC (static var compensator) with higher efficiency.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide an electromagnetic transient rapid simulation method based on a static var compensator, which can accurately reflect the dynamic behavior of SVC under normal, symmetrical and asymmetrical states of a system, greatly improve the simulation speed while ensuring the calculation precision and greatly improve the working efficiency.
The purpose of the invention is realized by the following technical scheme:
an electromagnetic transient rapid simulation method based on a static var compensator comprises the following steps:
acquiring basic parameters of an SVC device; the parameters to be acquired include: inductance L of TCR branch circuit, capacitance C of fixed capacitor FCsResistance R of filterfInductor LfAnd a capacitor CfExternal input voltage uA、uB、uCAnd a thyristor firing angle α;
because the phase voltages input into each phase line of the SVC are fundamental frequency sinusoidal quantities, and the generality is not lost, an expression of three-phase fundamental wave components of the connection bus voltage is shown as the following formula (1):
Figure BDA0001999709760000031
and introducing auxiliary variables as shown in the following formula (2):
Figure BDA0001999709760000032
step two, finding the initial value X of the state variable at the initial time, namely, at the time when t is 00
X0=[i1(0),i2(0),i3(0),uaf(0),iaf(0),ubf(0),ibf(0),ucf(0),icf(0),uA(0),uB(0),uC(0),vA(0),vB(0),vC(0)]T
Wherein:
Figure BDA0001999709760000033
step three, making k equal to 1;
step four, judging whether t is equal to 0 or not; if yes, turning to the step ten; if not, turning to the fifth step;
step five, t is t + [ delta ] t;
judging whether the on-off condition of the thyristor changes or not; if yes, go to step seven; if not, go to step eleven;
the method for judging the on-off condition of the thyristor specifically comprises the following steps: when one of the two thyristors corresponding to any one of the three branches of the TCR is switched on or off, the on-off state of the thyristor is considered to be changed; wherein, the conduction condition of the thyristor is as follows: when the thyristor acts with trigger pulse, and the polarity of the voltage at the two ends is positive; the turn-off conditions of the thyristor were: when two adjacent calculation time steps are different in sign of the inductive current, namely: i.e. ix(t-Δt)·ix(t) < 0, or when ix(t- Δ t) ≠ 0 and ixWhen (t) is 0, the thyristor is considered to be turned off at the time t, and the inductive current i at the time is juxtaposedx(t)=0;
Step seven, making the kth time interval end time
Figure BDA0001999709760000041
Step eight, k is k + 1;
step nine, let the k time interval start time
Figure BDA0001999709760000042
Step ten, forming a three-phase state equation
Figure BDA0001999709760000043
Wherein:
Figure BDA0001999709760000044
wherein, the matrix Z, F, W is a constant matrix, which is respectively:
Figure BDA0001999709760000051
Figure BDA0001999709760000052
Figure BDA0001999709760000053
wherein, the matrix LkThe matrix is changed along with the period and is formed according to the on-off condition of the TCR branch thyristor in the kth period;
Lk=[LAB,LBC,LCA]T
the three-phase TCR branch circuit thyristor on-off conditions comprise:
when the thyristor between the AB phases is in a conducting state, LAB=[1/L -1/L 0]When both thyristors between the AB phases are in the OFF state, LAB=[0 0 0];
When a thyristor is in a conducting state between BC phases, LBC=[0 1/L -1/L]When both thyristors between the BC phases are in the OFF state, LBC=[0 0 0];
When a thyristor is in a conducting state between the CA phases, LCA=[-1/L 0 1/L]When both thyristors between CA phases are in OFF state, LCA=[0 0 0];
Step eleven, equation of state
Figure BDA0001999709760000061
In that
Figure BDA0001999709760000062
The solution of time is
Figure BDA0001999709760000063
Step twelve, absorbing three-phase current Y ═ i by SVCa,ib,ic]TCan be determined by Y ═ BX, where matrix B is:
Figure BDA0001999709760000064
step thirteen, judging whether t is larger than tend(ii) a If yes, ending the program; if not, returning to the fourth step.
Compared with the prior art, the invention has the following beneficial effects:
(1) the method introduces a sine auxiliary variable, converts the state equation of the SVC from a non-homogeneous linear differential equation set into a homogeneous linear differential equation set, quickly obtains the solution of the equation set by a matrix exponential integration method, and meanwhile, the state equations under various working conditions are unified in form;
(2) compared with the traditional electromagnetic transient calculation in which the solution required system model changes at each time step, the method provided by the invention has the advantages that the solution required system model changes only when the on-off state of the thyristor changes, so that the calculation of the method provided by the invention is more efficient, and the working efficiency is greatly improved.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic three-phase connection diagram of a six-pulse TCR-FC type SVC device according to the present invention;
FIG. 3 is a schematic view of the firing angle α of the present invention;
FIG. 4 is a schematic diagram showing comparison of simulation results under symmetrical steady-state operation of the present invention;
FIG. 5 is a schematic diagram showing a comparison of simulation results when a three-phase power supply amplitude sudden change occurs in the present invention;
FIG. 6 is a diagram showing comparison of simulation results when the amplitude of the A-phase power supply suddenly changes according to the present invention;
FIG. 7 is a schematic diagram showing a comparison of simulation results under the condition of sudden change of the firing angle during symmetrical steady-state operation;
FIG. 8 is a diagram showing a comparison of simulation results when the amplitude of the A-phase power supply and the flip angle are suddenly changed at the same time.
Detailed Description
The present invention will be described in further detail with reference to examples and drawings, but the present invention is not limited thereto.
The invention takes three-phase fundamental frequency sinusoidal voltage as input of SVC, adopts a state equation of a state space method column writing system, converts a model from a non-homogeneous linear differential equation set into a homogeneous linear differential equation set through processing of each step, and the form of the homogeneous linear differential equation set can be changed when the on-off state of a thyristor is changed, so that the homogeneous linear differential equation set is divided into time intervals and can be solved through a matrix exponential integration method, thereby improving the calculation speed.
The form of the established homogeneous system of linear differential equations is as follows:
Figure BDA0001999709760000071
wherein X is a state variable of the system;
Figure BDA0001999709760000072
is a state matrix; k is a time period number; t is tk stIs the starting time of the kth time interval, and the initial value of the corresponding state variable is Xk st;tk edIs the end time of the kth time interval, and the corresponding state variable has the end value of Xk ed. The following relationship holds:
Figure BDA0001999709760000073
the solving idea and solving steps of the present invention are described as shown in the six-pulse TCR-FC type SVC device of FIG. 2. In FIG. 2, uaf、iaf、ubf、ibf、ucf、icfAnd the capacitance voltage and the inductance current of each phase of the filter branch circuit are represented.Let the simulation time period be t e [0, t ∈end]The calculation step is delta t (generally 50 mu s), and the three-phase input voltage u of the system is calculatedA、uB、uCThe inductance current i of the TCR branch is a known quantity1、i2、i3And SVC absorbed three-phase current ia、ib、icIs the amount to be requested.
As shown in fig. 1, a method for rapidly simulating an electromagnetic transient based on a static var compensator includes the following steps:
acquiring basic parameters of an SVC device; the parameters to be acquired include: inductance L of TCR branch circuit, capacitance C of fixed capacitor FCsResistance R of filterfInductor LfAnd a capacitor CfExternal input voltage uA、uB、uCAnd a thyristor firing angle α;
as shown in fig. 3, the firing angle α is the electrical angle from the TCR leg inductance voltage zero crossing to the time the firing pulse is emitted.
Because the phase voltages input into each phase line of the SVC are fundamental frequency sinusoidal quantities, and the generality is not lost, an expression of three-phase fundamental wave components of the connection bus voltage is shown as the following formula (1):
Figure BDA0001999709760000081
and introducing auxiliary variables as shown in the following formula (2):
Figure BDA0001999709760000082
in the above formulas (1) and (2), U is acquiredAm、UBm、UCm
Figure BDA0001999709760000083
And ω. Wherein, UAm、UBm、UCmAnd
Figure BDA0001999709760000084
respectively externally input three-phase voltage uA、uB、uCThe amplitude and phase angle of (d); ω is related to the frequency f of the system, ω 2 π f.
Step two, obtaining an initial value X of the state variable at the initial time, namely, the time when t is equal to 00
X0=[i1(0),i2(0),i3(0),uaf(0),iaf(0),ubf(0),ibf(0),ucf(0),icf(0),uA(0),uB(0),uC(0),vA(0),vB(0),vC(0)]T
Wherein:
Figure BDA0001999709760000091
step three, making k equal to 1;
step four, judging whether t is equal to 0 or not; if yes, turning to the step ten; if not, turning to the fifth step;
step five, t is t + [ delta ] t;
judging whether the on-off condition of the thyristor changes or not; if yes, go to step seven; if not, turning to the eleventh step;
the method for judging the on-off condition of the thyristor specifically comprises the following steps: when one of the two thyristors corresponding to any one of the three branches of the TCR is switched on or off, the on-off state of the thyristor is changed; wherein, the conduction condition of the thyristor is as follows: when the thyristor acts with trigger pulse, and the polarity of the voltage at the two ends is positive; the turn-off conditions of the thyristor were: when two adjacent calculation time steps are different in sign of the inductive current, namely: i all right anglex(t-Δt)·ix(t) < 0, or when ix(t- Δ t) ≠ 0 and ixWhen (t) is 0, the thyristor is considered to be turned off at the time t, and the inductive current i at the time is juxtaposedx(t) ═ 0; wherein, x can be 1, 2 or 3; i.e. ix(t) represents the current in any of the TCR legs.
Step seven, making the kth time interval end time
Figure BDA0001999709760000101
Step eight, k is k + 1;
step nine, let the k time interval start time
Figure BDA0001999709760000102
Step ten, forming a three-phase state equation
Figure BDA0001999709760000103
Wherein:
Figure BDA0001999709760000104
wherein, the matrix Z, F, W is a constant matrix, which is respectively:
Figure BDA0001999709760000105
Figure BDA0001999709760000111
Figure BDA0001999709760000112
wherein, the matrix LkThe matrix is changed along with the period and is formed according to the on-off condition of the TCR branch thyristor in the kth period;
Lk=[LAB,LBC,LCA]T
the three-phase TCR branch circuit thyristor on-off conditions comprise:
when the thyristor between the AB phases is in a conducting state, LAB=[1/L -1/L 0]When both thyristors between the AB phases are in the OFF state, LAB=[0 0 0];
When a thyristor is in a conducting state between BC phases, LBC=[0 1/L -1/L]When both thyristors between the BC phases are in the OFF state, LBC=[0 0 0];
When a thyristor is in a conducting state between the CA phases, LCA=[-1/L 0 1/L]When both thyristors between CA phases are in OFF state, LCA=[0 0 0];
Step eleven, equation of state
Figure BDA0001999709760000113
In that
Figure BDA0001999709760000114
The solution of time is
Figure BDA0001999709760000115
Step twelve, absorbing three-phase current Y ═ i by SVCa,ib,ic]TCan be determined by Y ═ BX, where matrix B is:
Figure BDA0001999709760000121
step thirteen, judging whether t is larger than tend(ii) a If yes, ending the program; if not, returning to the fourth step.
The following is a specific embodiment of the present invention:
taking the six-pulse TCR-FC type SVC device shown in fig. 2 as an example, the calculation result is compared with the calculation result of the present invention by using PSCAD/EMTDC electromagnetic transient simulation software. In the system, three-phase symmetrical voltage sources are input:
Figure BDA0001999709760000122
wherein U isAm=UBm=Ucm10kV, w 2 pi f 100 pi, inductance L0.0975H of TCR branch, filteringResistance R of the device branchf0.1 Ω, inductance Lf10H, capacitance CfThe capacitance C of the capacitor FC is fixed at 200 μ Fs3 μ F, flip angle α 100 °; the simulation calculation step length is 50 mus, and the simulation duration is 1 s.
(1) The system runs symmetrically and stably;
the simulation results are as follows: FIG. 4(a) shows the current i of the TCR branch1、i2、i3(b) gives the three-phase current i absorbed by the SVC from the systema、ib、icThe waveform of (a); the black dotted line is the calculation result of PSCAD/EMTDC electromagnetic transient software, and the black solid line is the calculation result of the invention.
(2) The system generates a three-phase voltage drop fault with the duration of 0.06s at a power supply end of 0.1s, and the voltage amplitude falls from 10kV to 7 kV;
the simulation results are as follows: FIG. 5(a) shows the current i of the TCR branch1、i2、i3(b) gives the three-phase current i absorbed by the SVC from the systema、ib、icThe waveform of (a); the black dotted line is the calculation result of PSCAD/EMTDC electromagnetic transient software, the black solid line is the calculation result of the invention, and the black dotted line is the time when the fault occurs.
(3) The system has an A-phase voltage dropping fault with the duration of 0.06s at a power supply end of 0.1s, and the voltage amplitude drops from 10kV to 7 kV;
the simulation results are as follows: FIG. 6(a) shows the current i of the TCR branch1、i2、i3(b) gives the three-phase current i absorbed by the SVC from the systema、ib、icThe waveform of (a); the black dotted line is the calculation result of PSCAD/EMTDC electromagnetic transient software, the black solid line is the calculation result of the invention, and the black dotted line is the time when the fault occurs.
(4) When the system operates symmetrically and stably, the triggering angle alpha is mutated at 0.1s and is changed from 100 degrees to 115 degrees;
the simulation results are as follows: FIG. 7(a) shows the current i for the TCR branch1、i2、i3(b) gives the three-phase current i absorbed by the SVC from the systema、ib、icThe waveform of (a); the black dotted line is the calculation result of PSCAD/EMTDC electromagnetic transient software, the black solid line is the calculation result of the invention, and the black dotted line is the moment of sudden change of the trigger angle.
(5) The system has an A-phase voltage drop fault with the duration of 0.06s at a power supply end of 0.1s, the voltage amplitude is dropped from 10kV to 7kV, and the trigger angle alpha is changed from 100 degrees to 115 degrees;
the simulation results are as follows: FIG. 8(a) shows the current i of the TCR branch1、i2、i3(b) gives the three-phase current i absorbed by the SVC from the systema、ib、icThe waveform of (a); the black dotted line is the calculation result of PSCAD/EMTDC electromagnetic transient software, the black solid line is the calculation result of the invention, and the black dotted line is the time when the fault occurs.
As can be seen from FIGS. 4 to 8, the result obtained by the method is well matched with the electromagnetic transient simulation result, and the effectiveness of the method is verified. Meanwhile, as can be seen from the following table 1, compared with the electromagnetic transient simulation program PSCAD, the simulation of the present invention is short in time consumption and has an obvious advantage in calculation speed.
TABLE 1 comparison of the invention with PSCAD simulation elapsed time
Figure BDA0001999709760000141
The invention introduces a sine auxiliary variable, converts the state equation of the SVC from a non-homogeneous linear differential equation set into a homogeneous linear differential equation set, quickly obtains the solution of the equation set by a matrix exponential integration method, and simultaneously, the state equations under various working conditions are unified in form; compared with the traditional electromagnetic transient calculation in which the solution required system model changes at each time step, the method provided by the invention has the advantages that the solution required system model changes only when the on-off state of the thyristor changes, so that the calculation of the method provided by the invention is more efficient, and the working efficiency is greatly improved.
The present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents and equivalents thereof, which are intended to be included in the scope of the present invention.

Claims (1)

1. An electromagnetic transient rapid simulation method based on a static var compensator is characterized by comprising the following steps:
acquiring basic parameters of an SVC device; the parameters to be acquired include: inductance L of TCR branch circuit, capacitance C of fixed capacitor FCsResistance R of filterfInductor LfAnd a capacitor CfExternal input voltage uA、uB、uCAnd a thyristor firing angle α;
because the phase voltages input into each phase line of the SVC are fundamental frequency sinusoidal quantities, and the generality is not lost, an expression of three-phase fundamental wave components of the connection bus voltage is shown as the following formula (1):
Figure FDA0003608170350000011
wherein t represents time, UAm、UBm、UCmRespectively A, B, C phase voltage amplitudes,
Figure FDA0003608170350000012
Figure FDA0003608170350000013
initial phase angle, v, of A, B, C phase voltageA、vB、vCA, B, C auxiliary variables, respectively;
and introducing auxiliary variables as shown in the following formula (2):
Figure FDA0003608170350000014
step two, solving the initial value X of the state variable at the moment when t is equal to 00
X0=[i1(0),i2(0),i3(0),uaf(0),iaf(0),ubf(0),ibf(0),ucf(0),icf(0),uA(0),uB(0),uC(0),vA(0),vB(0),vC(0)]T
Wherein:
Figure FDA0003608170350000021
step three, making k equal to 1, wherein k represents a time interval serial number;
step four, judging whether t is equal to 0 or not; if yes, turning to the step ten; if not, turning to the fifth step;
step five, t is t + [ delta ] t;
step six, judging whether the on-off condition of the thyristor changes; if yes, go to step seven; if not, turning to the eleventh step;
the method for judging the on-off condition of the thyristor specifically comprises the following steps: when one of the two thyristors corresponding to any one of the three branches of the TCR is switched on or off, the on-off state of the thyristor is considered to be changed; wherein, the conduction condition of the thyristor is as follows: when a trigger pulse acts on the thyristor, the polarity of the voltage at two ends is positive; when one of the following two conditions is met, the thyristor is turned off at the moment t, and the inductive current i at the moment is juxtaposedx(t)=0:
Two adjacent calculation time steps have different signs of the inductive current, and the formula is as follows: i.e. ix(t-Δt)·ix(t)<0;
② the inductive current in two adjacent calculation time steps conforms to ix(t- Δ t) ≠ 0 and ix(t)=0。
Step seven, the kth time interval end time is ordered
Figure FDA0003608170350000031
Step eight, k is k + 1;
step nine, let the k time interval start time
Figure FDA0003608170350000032
Step ten, forming a three-phase state equation
Figure FDA0003608170350000033
Wherein:
Figure FDA0003608170350000034
Lka matrix representing a random transformation;
wherein, the matrix Z, F, W is a constant matrix, which is respectively:
Figure FDA0003608170350000035
Figure FDA0003608170350000036
Figure FDA0003608170350000041
wherein, the matrix LkThe matrix is changed along with the period and is formed according to the on-off condition of the TCR branch thyristor in the kth period;
Lk=[LAB,LBC,LCA]T
the three-phase TCR branch circuit thyristor on-off conditions comprise:
when the thyristor between the AB phases is in a conducting state, LAB=[1/L -1/L 0]When both thyristors between the AB phases are in the OFF state, LAB=[0 0 0];
When a thyristor is in a conducting state between BC phases, LBC=[0 1/L -1/L]When both thyristors between the BC phases are in the OFF state, LBC=[0 0 0];
When a thyristor is in a conducting state between the CA phases, LCA=[-1/L 0 1/L]When both thyristors between CA phases are in OFF state, LCA=[0 0 0];
Step eleven, equation of state
Figure FDA0003608170350000042
In that
Figure FDA0003608170350000043
The solution of time is
Figure FDA0003608170350000044
Step twelve, absorbing three-phase current Y ═ i by SVCa,ib,ic]TCan be determined by Y ═ BX, where matrix B is:
Figure FDA0003608170350000045
step thirteen, judging whether t is larger than tend(ii) a If yes, ending the program; if not, returning to the fourth step.
CN201910208307.3A 2019-03-19 2019-03-19 Static reactive compensator-based electromagnetic transient rapid simulation method Expired - Fee Related CN110083863B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910208307.3A CN110083863B (en) 2019-03-19 2019-03-19 Static reactive compensator-based electromagnetic transient rapid simulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910208307.3A CN110083863B (en) 2019-03-19 2019-03-19 Static reactive compensator-based electromagnetic transient rapid simulation method

Publications (2)

Publication Number Publication Date
CN110083863A CN110083863A (en) 2019-08-02
CN110083863B true CN110083863B (en) 2022-07-05

Family

ID=67412579

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910208307.3A Expired - Fee Related CN110083863B (en) 2019-03-19 2019-03-19 Static reactive compensator-based electromagnetic transient rapid simulation method

Country Status (1)

Country Link
CN (1) CN110083863B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104820752A (en) * 2015-05-13 2015-08-05 中国南方电网有限责任公司电网技术研究中心 Transient response analytic calculation method for high-voltage direct current power transmission system
CN104852389A (en) * 2015-05-27 2015-08-19 西安工程大学 Static var compensator controlling device for improving system transient stability
CN105468864A (en) * 2015-12-14 2016-04-06 三峡大学 Electromagnetic transient numerical computation method of high-voltage power transmission line based on increment dimension precise integration
KR20180091491A (en) * 2017-02-07 2018-08-16 엘에스산전 주식회사 Static var compensator apparatus and operating method thereof
CN108536949A (en) * 2018-04-03 2018-09-14 上海交通大学 LCC-HVDC simulation models and parameter optimization method based on inductance/capacitance switch models

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6399462B2 (en) * 2014-01-25 2018-10-03 ▲コ▼▲トク▼科技服▲務▼(嘉▲興▼)有限公司Good Science And Technology Service (Jiaxing) Co.,Ltd. AC pressure regulating electronic switch

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104820752A (en) * 2015-05-13 2015-08-05 中国南方电网有限责任公司电网技术研究中心 Transient response analytic calculation method for high-voltage direct current power transmission system
CN104852389A (en) * 2015-05-27 2015-08-19 西安工程大学 Static var compensator controlling device for improving system transient stability
CN105468864A (en) * 2015-12-14 2016-04-06 三峡大学 Electromagnetic transient numerical computation method of high-voltage power transmission line based on increment dimension precise integration
KR20180091491A (en) * 2017-02-07 2018-08-16 엘에스산전 주식회사 Static var compensator apparatus and operating method thereof
CN108536949A (en) * 2018-04-03 2018-09-14 上海交通大学 LCC-HVDC simulation models and parameter optimization method based on inductance/capacitance switch models

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Analytical method for computing electromagnetic transient models of FACTS in hybrid simulation;Song Zou et al;《2016 IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC)》;20161231;第1316-1321页 *
HVDC电磁暂态仿真的算法研究;徐春明;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20111215;第1230-1239页 *
Methodology for rotor angle transient stability assessment in parameter space;Oluic M et al;《IEEE Trans. Power Syst》;20171231;第1202-1211页 *
Modelling and Operation Analysis of Electromechanical Transient- Electromagnetic Transient Hybrid System;Fang F et al;《Chemical Engineering Transactions》;20171231;第805-810页 *
含静止无功补偿器电力系统机电暂态和电磁暂态混合仿真;王立伟;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20040415;第C042-307页 *

Also Published As

Publication number Publication date
CN110083863A (en) 2019-08-02

Similar Documents

Publication Publication Date Title
CN102663174B (en) Simulation method of MMC (modular multilevel converter) and application thereof
CN106649927B (en) FPGA-based real-time simulation combined modeling method for power electronic element
CN106374767B (en) A kind of Modularized multi-level converter sub-module simulation model considering secondary circuit
CN202275142U (en) Back-to-back test device of a flexible direct-current transmission MMC valve steady-state operation test
CN104993715B (en) Modularization multi-level converter valve damages rapid evaluation algorithm
CN103605850A (en) MMC (modular multilevel converter) equivalent modeling method with module latching function
CN106950512A (en) A kind of energy accumulation current converter and off-network characteristic integrated detection system and method
CN108536949A (en) LCC-HVDC simulation models and parameter optimization method based on inductance/capacitance switch models
CN106886617A (en) A kind of multi tate electro-magnetic transient subnetting method containing many VSC
CN110083863B (en) Static reactive compensator-based electromagnetic transient rapid simulation method
CN110895644A (en) Simulation system for large-scale converter group grid connection and equivalent circuit of converter
CN110059346B (en) Electromagnetic transient rapid simulation method of controllable series compensation device
Shujun et al. Modeling for VSC-HVDC electromechanical transient based on dynamic phasor method
CN107257205A (en) A kind of MMC power models nonlinear characteristic simulation model
CN111898282B (en) Improved modularized multi-level converter Thevenin equivalent modeling method
CN113190993B (en) Electromagnetic transient rapid simulation modeling method and system for high-voltage high-power energy storage converter
CN103487702A (en) Small-power movable micro-grid connection detecting system
Guo et al. Analysis on voltage distribution performance of HVDC thyristor valves
Liu et al. Overview of modeling schemes for MMC converter in flexible DC distribution network
CN110046446A (en) A kind of operating status control method, device, medium and the equipment of converter
CN104267615A (en) MMC flexible direct-current system
Wang et al. Analysis and optimization of multi-DC current flow controllers with different configuration positions in M 2 TDC grid
Xing et al. Research on fault analysis technology for PV system integration with dc booster system based on RTDS simulator
CN110896229A (en) Simulation system for large-scale converter group grid connection
CN204967269U (en) Control of power quality device is synthesized to hybrid

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220705

CF01 Termination of patent right due to non-payment of annual fee