CN110086175B - Direct-current power fluctuation analog simulation method suitable for PSASP - Google Patents

Direct-current power fluctuation analog simulation method suitable for PSASP Download PDF

Info

Publication number
CN110086175B
CN110086175B CN201910279788.7A CN201910279788A CN110086175B CN 110086175 B CN110086175 B CN 110086175B CN 201910279788 A CN201910279788 A CN 201910279788A CN 110086175 B CN110086175 B CN 110086175B
Authority
CN
China
Prior art keywords
power
direct current
reactive
load
impact load
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.)
Active
Application number
CN201910279788.7A
Other languages
Chinese (zh)
Other versions
CN110086175A (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.)
Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
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 Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Priority to CN201910279788.7A priority Critical patent/CN110086175B/en
Publication of CN110086175A publication Critical patent/CN110086175A/en
Application granted granted Critical
Publication of CN110086175B publication Critical patent/CN110086175B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Abstract

The invention relates to the technical field of power system simulation, in particular to a direct current power fluctuation simulation method suitable for PSASP, which is characterized by comprising the following steps: calculating the power grid load flow in a normal operation state by adopting PSASP, and determining the direct current power of a direct current system; a constant load is adopted to replace a direct current system; adding impact load and setting an impact load power value, which specifically comprises the following steps: adding impact load at a converter bus of the inverter station, selecting a disturbance type as the impact load, selecting the added impact load, and setting an active value and a reactive value of the impact load according to the known direct current output power fluctuation characteristic data required by a user; setting a reactive compensation cutting sequence; and carrying out simulation calculation and outputting the direct current power required by the user. The invention can accurately simulate the DC output power fluctuation characteristic required by a user, and has simple and convenient operation.

Description

Direct-current power fluctuation analog simulation method suitable for PSASP
Technical Field
The invention relates to the technical field of power system simulation, in particular to a direct-current power fluctuation simulation method suitable for PSASP.
Background
High voltage direct current transmission (HVDC) has the advantages of large transmission capacity, small transmission loss, long transmission distance and the like, so the development is very rapid. Because the transmission capacity of the direct current transmission is large, the output power of the direct current transmission fluctuates greatly after the direct current transmission is subjected to large disturbance; therefore, research on alternating current and direct current interaction influence is necessary for safe and stable operation of a power grid.
The existing direct current system model with direct current fault simulation capability comprises an electromagnetic transient model and a quasi-steady state model. The electromagnetic transient model can simulate the dynamic characteristics of direct current in detail, but the simulation scale is small, so that the electromagnetic transient model is not suitable for simulation of a large-scale actual power grid; the quasi-steady-state model can completely reserve the scale of the alternating current system, so the quasi-steady-state model is widely applied to electromechanical transient simulation programs, but is not accurate enough; considering that the operation and scheduling of an actual power grid usually use a PSASP (power system analysis software package) as a platform, but the direct current system in the PSASP at present adopts a quasi-steady-state model, so that the output power characteristic of the direct current system under disturbance cannot be well simulated, and the simulation result of the direct current system has a large input and output with the requirement of a user; meanwhile, although the PSASP has a custom module, redefining the dc model takes into account various problems under static and transient states, which is too complicated.
Therefore, if the user has an accurate set of known data (such as measured data collected from WAMS or electromagnetic transient simulation data obtained from PSCAD), the dc output power fluctuation characteristics required by the user are obtained, but it is difficult to put the data into PSASP for simulation calculation.
Disclosure of Invention
The invention aims to provide a direct current power fluctuation simulation method suitable for PSASP, which can accurately simulate the direct current output power fluctuation characteristics required by users and is simple and convenient to operate.
In order to solve the technical problems, the technical scheme of the invention is as follows: a direct current power fluctuation simulation method suitable for PSASP comprises the following steps:
step 1: calculating the power grid load flow in a normal operation state by adopting PSASP, and determining the direct current power of a direct current system, wherein the method specifically comprises the following steps: in the invention, the existing AC/DC system calculation example adopts a PSASP (power system analysis software package) self-carried DC model and further determines the transmitting end active power P of a DC system 0R Sending end consuming reactive power Q 0R And receiving end active power P 0I Receiving end consuming reactive power Q 0I 。
It should be noted that the self-contained direct current model in the PSASP is a quasi-steady-state model, so that the accuracy is better during load flow calculation, and the defect is that the power response is inaccurate under disturbance, so that the self-contained PSASP model can be used during calculation of the load flow of the whole alternating current and direct current power grid;
step 2: the method adopts constant load to replace a direct current system, and specifically comprises the following steps: in an alternating current-direct current power grid, if the alternating current-direct current power grid normally operates and no fault occurs, a direct current system is equivalent to injecting active power P into a receiving end alternating current system for the receiving end alternating current system 0I And consumes reactive power Q 0I That is, it can be equivalent to a constant-power load under the converter bus, and the active power of the constant-power load is-P 0I (minus sign indicates active power sent to the converter bus) and reactive power is Q 0I (ii) a Similarly, for a sending-end alternating current system, a direct current system can be equivalent to an active power P 0R And the reactive power is Q 0R The load of (2);
therefore, a load is added in the PSASP element database, the node name of the load is selected as an inversion station current conversion bus and is an inversion side load, and meanwhile, the active power of the inversion side load is set to be P 0I With reactive power set to Q 0I Setting an upper limit and a lower limit; similarly, active power is added at a converter bus of the rectifier station to be P 0R And the reactive power is Q 0R The load of (1) is a rectification side load; so far, replacing a direct current system with a constant load is completed;
and 3, step 3: adding impact load and setting an impact load power value; the method specifically comprises the following steps:
step 3.1: adding impact load at a converter bus of the inverter station: in the transient process, the direct current power between the direct current system and the receiving end alternating current system comprises two parts, wherein the first part is the direct current initial power when no disturbance exists, the constant load is adopted to represent, and the active power P of the constant load 0 ＝﹣P 0I Reactive power Q 0 ＝Q 0I (ii) a The second part is a direct current power change amount during disturbance, and after the system is disturbed, the power change amount generated by the direct current system can be regarded as power impact at a converter bus of the inverter station, so that in order to research the influence of impact load on a receiving end alternating current power grid, an inverter side and a rectifier side are isolated, and the impact load is arranged at the converter bus of the inverter station to simulate direct current powerFluctuation of (2); it should be noted that the added impact load cannot affect the power flow of the alternating current and direct current power grid, so that the active value and the reactive value of the impact load are both 0, and the power flow of the original alternating current and direct current power grid cannot be changed during power flow calculation; similarly to the step 2, adding impact load with the active power of 0 and the reactive power of 0 at the position of a converter bus of the inverter station, and giving a range of the active power and the reactive power of the impact load, wherein the active power does not exceed +/-P 0 Specifically, is. + -. P 0I Reactive power not exceeding +/-Q 0 Specifically, is. + -. Q 0I (ii) a Step 3.2: and (3) setting the temporary stable operation, selecting the disturbance type as the impact load, and selecting the added impact load: when a user uses the transient stability calculation function of the PSASP to simulate and research an alternating current-direct current power grid in actual work, the power output by a direct current system can hardly meet the characteristics required by the user through fault setting; in order to simulate the fluctuation characteristic of the direct current output power under the fault, the disturbance type is selected as the impact load in the transient stability operation, and the added impact load name is selected;
step 3.3: setting the active value P (t) and the reactive value Q (t) of the impact load according to the known data with the DC output power fluctuation characteristic required by the user:
defining the active power output by the direct current system as P d (t), the active power output by the direct current system in the known data is P d ' (t) the positive direction is defined as that the direct current system transmits active power to the inversion station converter bus, and when no impact load exists, P d (t)＝-(-P 0 )＝P 0 (ii) a After impact load is added to a converter bus of the inverter station, P d (t)＝-(-P 0 +P(t))＝P 0 -P (t), so that the active power fluctuation of the dc system output satisfies P d ' (t), the active value calculation formula of the impact load is P (t) ═ P 0 -P d ' (t), calculating the active value of the impact load at each time point according to the formula, and setting the active value in an impact load setting interface;
defining the reactive power consumed by the inverter station of the direct current system as Q d (t), the reactive power consumed by the DC system inverter station in the known data is Q d ' (t) the positive direction is defined as the reverse direction of the DC systemThe converter bus of the transformer station absorbs idle work, and Q is obtained when no impact load is applied d (t)＝ Q 0 (ii) a When impact load is added at a converter bus of the inverter station, the impact load consumes reactive power Q d (t)＝Q 0 + Q (t), so that the reactive power consumed by the DC system satisfies Q d ' (t), the reactive value calculation formula of the impact load is Q (t) ═ Q d ’(t)-Q 0 Calculating the reactive value of the impact load at each time point according to the formula, and setting the reactive value in an impact load setting interface; and 4, step 4: setting a reactive compensation cutting sequence, specifically: defining the active power exchanged between the DC system and the receiving end AC system as P s Defining the reactive power exchanged between the DC system and the receiving end AC system as Q s And the positive direction is defined as that the direct current system transmits power to the receiving end alternating current system; the exchange active power being equal to the active power output by the DC system, i.e. P s ＝P d (t)＝P 0 -p (t); the reactive power exchanged includes the reactive power Q consumed by the DC system inverter station d (t) and reactive power Q provided by reactive power compensation equipment f (t), i.e. Q s ＝Q f (t)-Q d (t) wherein Q d (t)＝Q 0 + Q (t); therefore, if the reactive power of the direct current system at the direct current side switched to the receiving end alternating current system at the alternating current side needs to be simulated correctly, the reactive power compensation equipment cutting sequence needs to be set, the fault type is selected in the transient and stable operation, the 'circuit switching' is selected from the fault type, the current conversion bus where the reactive power equipment to be cut is located is selected, and the cutting time is set;
and 5: carrying out simulation calculation, and outputting direct current power required by a user, specifically: performing load flow calculation again, and performing transient simulation calculation according to the transient operation set in the step, wherein the power output at the converter bus of the inverter station in the calculation result is P of the active power exchanged between the direct current system and the receiving end alternating current system s The reactive power exchanged between the DC system and the receiving end AC system is Q s I.e. the dc power required by the user.
According to the scheme, in the step 1, the constant load is arranged at the position of the converter bus of the inverter stationThe upper and lower limits of the power are set to P 0I The upper and lower limits of the reactive power are set to Q 0I (ii) a The upper limit and the lower limit of the active power of the constant load at the position of a converter bus of the rectifier station are set to be P 0R The upper limit and the lower limit of the reactive power are set to be Q 0R (ii) a The upper limit and the lower limit of the constant load are both the same as the corresponding power value so as to ensure that the power of the load cannot change in the transient process;
compared with the prior art, the invention has the following beneficial effects: the invention provides a direct current power fluctuation simulation method suitable for PSASP (power system analysis software package) aiming at two problems that a direct current model in PSASP (power system analysis software package) cannot well simulate the output response characteristic of a direct current system under disturbance and a newly defined direct current model is too complex, wherein the application range specifically aims at researching the influence of the direct current system power fluctuation on a receiving end alternating current system at an alternating current side; the simulation method has practicability and universality, can be applied to simulation research of the PSASP on the AC/DC power grid, and provides higher accuracy for researching the dynamic behavior of the AC/DC system under large disturbance in the PSASP.
Drawings
FIG. 1 is a flow chart of a simulation method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the system after impact load is added;
FIG. 3a is an exemplary plot of constant load rectified side load in PSASP;
FIG. 3b is an exemplary diagram of the inverter side load for a constant load in PSASP;
FIG. 4 is an exemplary diagram of adding impact load to a PSASP;
FIG. 5 is an exemplary diagram of transient stability node disturbance data setup in PSASP;
FIG. 6a is a graph showing the fluctuation of the active power output from the DC system in the known data;
FIG. 6b is a graph of the fluctuation of reactive power consumed by the DC system in the known data;
FIG. 7 is an exemplary diagram of impact load power value settings in PSASP;
FIG. 8 is an exemplary diagram of a PSASP with reactive device cut-off;
fig. 9a is a comparison graph of active power exchanged between the dc system and the receiving ac system and a user's desired curve after the simulation method is applied in the PSASP;
fig. 9b is a comparison graph of reactive power exchanged between the dc system and the ac system at the receiving end and the curve assumed by the user after the simulation method is applied in the PSASP.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
Referring to fig. 1 to 9b, the present invention is a direct current power fluctuation simulation method suitable for a PSASP, which includes the following steps:
step 1: calculating the power flow of the power grid in a normal operation state, and determining the direct current power of the direct current system;
in the invention, the existing AC/DC system calculation example adopts a PSASP (power system analysis software package) self-carried DC model and further determines the transmitting end active power P of a DC system 0R Sending end consuming reactive power Q 0R And receiving end active power P 0I Receiving end consuming reactive power Q 0I 。
It should be noted that the self-contained direct current model in the PSASP is a quasi-steady-state model, so that the accuracy is better during load flow calculation, and the defect is that the power response is inaccurate under disturbance, so that the self-contained PSASP model can be used during the calculation of the load flow of the whole alternating current and direct current power grid.
Step 2: a constant load is adopted to replace a direct current system;
referring to fig. 2, in the ac/dc power grid, if the ac/dc power grid is operating normally and no fault occurs, the dc system is applied to the receiving endThe alternating current system is equivalent to injecting active power P into the receiving end alternating current system 0I And consumes reactive power Q 0I That is, it can be equivalent to a constant-power load under the converter bus, and the active power of the constant-power load is-P 0I (minus sign indicates active power sent to the converter bus) and reactive power is Q 0I (ii) a Similarly, for a sending-end alternating current system, a direct current system can be equivalent to an active power P 0R And the reactive power is Q 0R The load of (2);
therefore, a load is added in the PSASP element database, the node name of the load is selected as an inversion station current conversion bus and is an inversion side load, and meanwhile, the active power of the inversion side load is set to be P 0I With reactive power set to P 0I And respectively setting the upper limit and the lower limit to the same value, namely setting the upper limit and the lower limit of the active power of the constant load at the position of a converter bus of the inverter station to be P 0I The upper and lower limits of the reactive power are set to Q 0I To ensure that the power of the load does not change during transients; similarly, active power is added at a converter bus of the rectifier station to be P 0R And the reactive power is Q 0R The load of (1) is a rectification side load; the upper limit and the lower limit of the active power of the constant load at the position of a converter bus of the rectifier station are set to be P 0R The upper limit and the lower limit of the reactive power are set to be Q 0R (ii) a Meanwhile, the original direct current system is removed. So far, the adoption of a constant load instead of a direct current system has been completed.
And step 3: adding impact load and setting an impact load power value;
step 3.1: adding impact load at a converter bus of the inverter station: referring to fig. 2, in the transient process, the dc power between the dc system and the receiving ac system includes two parts, the first part is the dc initial power without disturbance, which is represented by the constant load, and the active power P of the constant load 0 ＝﹣P 0I Reactive power Q 0 ＝Q 0I (ii) a The second part is a direct current power change quantity during disturbance, and after the system is disturbed, the power change quantity generated by the direct current system can be regarded as the power at a converter bus of the inverter stationTherefore, in order to research the influence of impact load on a receiving end alternating current power grid, an inversion side and a rectification side are isolated, and the impact load is arranged at a converter bus of an inversion station to simulate the fluctuation of direct current power.
It should be noted that the added impact load cannot affect the power flow of the ac/dc power grid, so the active value and the reactive value of the impact load are both 0, and thus the power flow of the original ac/dc power grid is not changed during power flow calculation.
Similarly to the step 2, adding impact load with the active power of 0 and the reactive power of 0 at the position of a converter bus of the inverter station, and giving a range of the active power and the reactive power of the impact load, wherein the active power does not exceed +/-P 0 Specifically, is. + -. P 0I Reactive power not exceeding +/-Q 0 Specifically, is. + -. Q 0I 。
Step 3.2: the transient and stable operation setting specifically comprises the following steps: the disturbance type is selected as the impact load, and the added impact load is selected: when a user uses the transient stability calculation function of the PSASP to simulate and research an alternating current-direct current power grid in actual work, the power output by a direct current system can hardly meet the characteristics required by the user through fault setting; in order to simulate the fluctuation characteristic of the direct current output power under the fault, the disturbance type is selected as the impact load in the transient stability operation, and the added impact load name is selected;
step 3.3: setting the power value of the impact load, namely the active value P (t) and the reactive value Q (t) of the impact load according to the known data with the DC output power fluctuation characteristic required by the user:
referring to fig. 2, the active power output by the dc system is defined as P d (t), the active power output by the direct current system in the known data is P d ' (t) the positive direction is defined as that the direct current system transmits active power to the inversion station converter bus, and when no impact load exists, P d (t)＝-(- P 0 )＝P 0 (ii) a After impact load is added to a converter bus of the inverter station, P d (t)＝-(-P 0 +P(t))＝P 0 -P (t), so that the active power fluctuation of the DC system output satisfies P d ' (t), the active value calculation formula of the impact load is P (t) ═ P 0 -P d ’(t), calculating the active value of the impact load at each time point according to the formula, and setting the active value in an impact load setting interface.
Referring to fig. 2, the reactive power consumed by the inverter station of the dc system is defined as Q d (t), the reactive power consumed by the DC system inverter station in the known data is Q d ' (t) the positive direction is defined as that the DC system absorbs reactive power from the inversion station converter bus, and when no impact load is applied, Q d (t)＝Q 0 (ii) a When impact load is added at a converter bus of the inverter station, the impact load consumes reactive power Q d (t)＝ Q 0 + Q (t), so if the reactive power consumed by the inverter station of the DC system is required to satisfy Q d ' (t), the reactive value calculation formula of the impact load is Q (t) ═ Q d ’(t)-Q 0 And calculating the reactive value of the impact load at each time point according to the formula, and setting the reactive value in an impact load setting interface.
And 4, step 4: setting a reactive compensation cutting sequence:
referring to fig. 2, the active power exchanged between the dc system and the receiving ac system is defined as P s Defining the reactive power exchanged between the DC system and the receiving end AC system as Q s And the positive direction is defined as that the direct current system transmits power to the receiving end alternating current system; as can be seen, the exchange active power is equal to the active power output by the DC system, i.e., P s ＝P d (t)＝P 0 -p (t); the reactive power exchanged consists of two parts, reactive power Q consumed by the DC system d (t) and reactive power Q provided by the reactive power compensation equipment f (t), i.e. Q s ＝Q f (t)-Q d (t), i.e. Q s ＝Q f (t)-Q d (t) wherein Q d (t)＝Q 0 +Q(t)；
Therefore, if the reactive power of the direct current system at the direct current side to be exchanged to the receiving end alternating current system at the alternating current side needs to be simulated correctly, the cutting sequence of the reactive compensation equipment needs to be set; and selecting a fault type in the transient stability operation, selecting a 'circuit switching' in the fault type, selecting a current conversion bus where the reactive equipment to be switched off is located, and setting the switching-off time.
And 5: carrying out simulation calculation, and outputting direct current power required by a user:
because a constant load is adopted to replace a direct current system, load flow calculation needs to be carried out again, then transient simulation calculation is carried out according to transient stability operation set in the steps, the power output at a converter bus of the inverter station in the calculation result is direct current power required by a user, namely the active power exchanged between the direct current system and a receiving end alternating current system in the step 5 is P s Reactive power Q exchanged between DC system and receiving end AC system s 。
Fig. 3a to 9b show an embodiment of the dc power fluctuation simulation method suitable for the PSASP of the present invention, which specifically includes the following steps:
step 1: calculating the power flow of the power grid, and determining the direct current power of the direct current system;
the Henan power grid is a typical alternating current-direct current interconnected power grid, and a large amount of power is transmitted to the Henan power grid through direct current in the sky; taking the Henan power grid as an example, calculating the load flow in the PSASP, and calculating the active power P of the sending end of the day-middle direct current system in a certain running state 0R Is 3200MW, the transmitting end consumes the reactive power Q 0R At 1333Mvar, active power P of receiving end 0I 3100MW receiving end consumes reactive power Q 0I For 1880Mvar, the reactive compensation equipment of the inversion station comprises 8 groups, and the total compensated reactive power is about 2080 MW.
Step 2: a constant load is adopted to replace a direct current system;
referring to fig. 3a and 3b, fig. 3a is an exemplary graph of a commutation side load of a constant load in the PSASP, and fig. 3b is an exemplary graph of an inversion side load of a constant load in the PSASP; a constant load with active power of 3200MW and reactive power of 1333Mvar is set at a converter bus of the rectifier station, and the upper limit and the lower limit are set to be the same values respectively; a constant load with active power of-3100 MW and reactive power of 1880Mvar is set at a converter bus of the inverter station, and the upper limit and the lower limit are set to be the same values respectively; then, the original DC system is removed.
And step 3: adding impact load and setting an impact load power value;
step 3.1: adding impact load at a converter bus of the inverter station: referring to fig. 4, a surge load is added at the inverter station converter bus "midwim 500", the surge load is named "load midwim", both the active power value and the reactive power value are set to 0, the active power range is set to ± 31.0p.u. (± 3100MW), and the reactive power range is set to ± 18.8p.u. (± 1880 Mvar).
Step 3.2: and (3) setting a transient and stable operation: referring to fig. 5, the type of disturbance is selected as the impact load in the transient operation, and the impact load "load middle" added is selected.
Step 3.3: and setting the active value and the reactive value of the impact load according to the known data with the DC output power fluctuation characteristics required by a user:
after the direct current is subjected to single-pole blocking in the same day, the direct current active power P which has the direct current output power fluctuation characteristic data required by a user, namely the direct current active power P required by the user is known d ' (t) and consumed reactive power Q d ' (t) the fluctuation characteristics are shown in fig. 6a and 6b, and are based on the active value formula P (t) ═ P of the impact load 0 -P d ' (t) and reactive value of impact load equation Q (t) Q d ’(t)-Q 0 And the tidal current value P obtained in the step 1 0 And Q 0 The active value p (t) and the reactive value q (t) of the impact load are calculated and set in the PSASP, as shown in fig. 7.
And 4, step 4: setting a reactive compensation cutting sequence;
referring to fig. 8, in the present embodiment, four groups of reactive devices are cut off at 1.2s after the dc blocking in the day, and according to the exemplary diagram of fig. 8, the commutation bus where each group of reactive devices is located is selected, and the cutting-off time is determined.
And 5: carrying out simulation calculation and outputting direct current power required by a user;
after the impact load is set, carrying out load flow calculation again, carrying out transient simulation calculation according to the transient stability operation set in the step, and observing the tie line power between the direct current system and the receiving end alternating current system, namely the power output at the position of a converter bus of the inverter station; the connecting line between the direct current and the Henan power grid in the sky is shown as a dotted line in fig. 9a and fig. 9b, the solid line is an exchange power curve obtained by calculation of a user, and comparison between the solid line and the dotted line shows that the characteristic of direct current power fluctuation required by the user can be simulated accurately by adopting the embodiment, and the method is simple and practical to use and can be applied to simulation research on an alternating current-direct current power grid in PSASP.
The foregoing is a more detailed description of the present invention that is presented in conjunction with specific embodiments, and the practice of the invention is not to be considered limited to those descriptions. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (2)

1. A direct current power fluctuation simulation method suitable for PSASP is characterized by comprising the following steps:
step 1: calculating the power grid load flow in a normal operation state by adopting PSASP, and determining the direct current power of a direct current system, wherein the method specifically comprises the following steps: calculating the load flow of the whole power grid by adopting a PSASP (Power System analysis software Package) self-contained direct current model, and determining the transmitting end active power P of a direct current system 0R Sending end consuming reactive power Q 0R Receiving end active power P 0I Reactive power Q consumed by receiving end 0I ；
Step 2: the method adopts constant load to replace a direct current system, and specifically comprises the following steps: adding a load in a PSASP element database, selecting a node name of the load as an inversion station converter bus, and setting the active power of the load to be-P 0I With reactive power set to Q 0I Setting an upper limit and a lower limit; similarly, active power is added at a converter bus of the rectifier station to be P 0R And the reactive power is Q 0R The load of (2); removing the original direct current system;
and step 3: adding impact load and setting an impact load power value; the method comprises the following specific steps:
step 3.1: adding impact load at a converter bus of the inverter station: in the transient process, the direct current power between the direct current system and the receiving end alternating current system comprises the direct current initial power when no disturbance exists and the direct current power change amount when disturbance exists, wherein the direct current initial power adopts constantLoad representation, active power P of constant load 0 ＝﹣P 0I Reactive power Q 0 ＝Q 0I (ii) a The direct current power change quantity is regarded as power impact at a converter bus of the inverter station, so that impact load with the active power of 0 and the reactive power of 0 is added at the converter bus of the inverter station, the active power and the reactive power of the impact load are set to be within a range, and the active power does not exceed +/-P 0 Reactive power not exceeding +/-Q 0 ；
Step 3.2: setting a transient and stable operation, selecting a disturbance type as an impact load, and selecting the added impact load;
step 3.3: setting the active value and the reactive value of the impact load according to the known data with the DC output power fluctuation characteristic required by a user:
defining the active power output by the direct current system as P d (t), the active power output by the direct current system in the known data is P d ' (t), the positive direction is defined as that the direct current system transmits active power to the inversion station converter bus, and the active power value calculation formula of the impact load is P (t) ═ P 0 -P d ' (t), calculating the active value of the impact load at each time point according to the formula, and setting the active value in an impact load setting interface;
defining the reactive power consumed by the inverter station of the direct current system as Q d (t), the reactive power consumed by the DC system inverter station in the known data is Q d ' (t), the positive direction is defined as that the direct current system absorbs reactive power from the inversion station converter bus, and the reactive power value calculation formula of the impact load is Q (t) ═ Q d ’(t)-Q 0 Calculating the reactive value of the impact load at each time point according to the formula, and setting the reactive value in an impact load setting interface;
and 4, step 4: setting a reactive compensation cutting sequence, specifically: defining the active power exchanged between the DC system and the receiving end AC system as P s Defining the reactive power exchanged between the DC system and the receiving end AC system as Q s And the positive direction is defined as that the direct current system transmits power to the receiving end alternating current system; the exchange active power being equal to the active power output by the DC system, i.e. P s ＝P d (t)＝P 0 -p (t); the reactive power exchanged includesReactive power Q consumed by inversion station of flow system d (t) and reactive power Q provided by reactive power compensation equipment f (t), i.e. Q s ＝Q f (t)-Q d (t) wherein Q d (t)＝Q 0 + Q (t); therefore, a reactive compensation cutting sequence needs to be set, a fault type is selected in transient and stable operation, a 'circuit switching' is selected from the fault type, a current conversion bus where reactive equipment to be cut is located is selected, and cutting time is set;
and 5: performing simulation calculation, and outputting the direct current power required by the user, specifically: performing load flow calculation again, and performing transient simulation calculation according to the set transient operation, wherein the power output at the inversion station conversion bus in the calculation result is P as the active power exchanged between the direct current system and the receiving end alternating current system s The reactive power exchanged between the DC system and the receiving end AC system is Q s I.e. the dc power required by the user.
2. The method of claim 1, wherein the method comprises: in the step 1, the upper and lower limits of the active power of the constant load at the position of the inversion station current conversion bus are set to be P 0I The upper and lower limits of the reactive power are set to Q 0I (ii) a The upper limit and the lower limit of the active power of the constant load at the position of a converter bus of the rectifier station are set to be P 0R The upper limit and the lower limit of the reactive power are set to be Q 0R 。
CN201910279788.7A 2019-04-09 2019-04-09 Direct-current power fluctuation analog simulation method suitable for PSASP Active CN110086175B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910279788.7A CN110086175B (en) 2019-04-09 2019-04-09 Direct-current power fluctuation analog simulation method suitable for PSASP

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910279788.7A CN110086175B (en) 2019-04-09 2019-04-09 Direct-current power fluctuation analog simulation method suitable for PSASP

Publications (2)

Publication Number Publication Date
CN110086175A CN110086175A (en) 2019-08-02
CN110086175B true CN110086175B (en) 2022-08-23

Family

ID=67414679

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910279788.7A Active CN110086175B (en) 2019-04-09 2019-04-09 Direct-current power fluctuation analog simulation method suitable for PSASP

Country Status (1)

Country Link
CN (1) CN110086175B (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103605871B (en) * 2013-12-03 2016-09-21 国家电网公司 A kind of Bulk power system simulation based on direct-current external characteristic matching analyzes method
CN108470091A (en) * 2018-02-24 2018-08-31 国网辽宁省电力有限公司 PSASP hybrid dynamic simulation methods based on the injection of adaptive impact load

Also Published As

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

Similar Documents

Publication Publication Date Title
CN102403916B (en) Design method of simulation accelerating circuit
CN102035203B (en) Fast calculation method of tie-line transient state stable limiting transmission power
SE540926C2 (en) Converter station power set point analysis system and method
Flatabo et al. A method for calculation of margins to voltage instability applied on the Norwegian system for maintaining required security level
CN110198044B (en) LCC-HVDC commutation failure immunity evaluation method and system considering transmission power limit
CN106202793B (en) A kind of mains frequency stability simulation method based on primary frequency modulation limitation
CN102403716A (en) Dynamic equalizing method for multi-infeed alternating/direct-current power grid
CN103972900A (en) Method for determining distribution of reactive compensation devices of multi-feed direct current transmission system based on voltage control sensitive factors
CN106814265A (en) A kind of photovoltaic DC-to-AC converter generating efficiency test system
CN104166753A (en) Method for measuring strength index of interaction between converter stations of multi-feed direct current system
CN110401184A (en) Multi-infeed DC receiving end power grid emergency control optimization method and system
CN106655195B (en) Active power distribution network high-frequency harmonic power flow calculation method
Li et al. An integrated online dynamic security assessment system for improved situational awareness and economic operation
CN104485661A (en) Tidal current automatic regulating method based on node type conversion
CN117081160B (en) Parallel-off-grid switching system for micro-grid
Shun et al. Influence of VSC HVDC on transient stability: Case study of the Belgian grid
CN110086175B (en) Direct-current power fluctuation analog simulation method suitable for PSASP
Hou et al. Measurement-based load modeling at distribution level with complete model structure
CN107370174A (en) A kind of HVDC transmission system simplifies modeling method
CN104993466A (en) Cascading fault fast dynamic simulation method applicable to alternating current-direct current power grid
Polster et al. Voltage stability monitoring using a modified thevenin impedance
CN109829178B (en) Transformer substation real-time simulation method based on model segmentation method
CN109327027A (en) Scissors crossing line fault check method based on PSD-BPA
CN104242354A (en) New energy resource centralized conveying running characteristic assessment method taking wind and light output relevance into consideration
CN108062449B (en) Power system external equivalence method based on PSD-BPA auxiliary switch

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