CN107317345B - Method for controlling frequency of isolated power grid by participation of electrolysis load - Google Patents

Method for controlling frequency of isolated power grid by participation of electrolysis load Download PDF

Info

Publication number
CN107317345B
CN107317345B CN201710681615.9A CN201710681615A CN107317345B CN 107317345 B CN107317345 B CN 107317345B CN 201710681615 A CN201710681615 A CN 201710681615A CN 107317345 B CN107317345 B CN 107317345B
Authority
CN
China
Prior art keywords
unit
load
electrolytic aluminum
power
formula
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
CN201710681615.9A
Other languages
Chinese (zh)
Other versions
CN107317345A (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.)
State Grid Corp of China SGCC
Wuhan University WHU
State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Wuhan University WHU
State Grid Hubei 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 State Grid Corp of China SGCC, Wuhan University WHU, State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201710681615.9A priority Critical patent/CN107317345B/en
Publication of CN107317345A publication Critical patent/CN107317345A/en
Application granted granted Critical
Publication of CN107317345B publication Critical patent/CN107317345B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/382
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/388Islanding, i.e. disconnection of local power supply from the network
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

The invention belongs to the technical field of power system operation and control, and particularly relates to a method for enabling electrolysis loads to participate in isolated power grid frequency control, wherein system frequency adjustment is realized through an electrolytic aluminum load control system, and the electrolytic aluminum load control system comprises a network control unit NCU, a dead zone unit, a proportional amplification unit, a WASHOUT unit and an amplitude limiting unit; the proportion amplification unit is connected with the WASHOUT unit in parallel, one end of the proportion amplification unit is sequentially connected with the dead zone unit and the network control unit NCU, and the other end of the proportion amplification unit is connected with the amplitude limiting unit; the network control unit NCU is connected with the isolated power grid system, and the amplitude limiting unit is connected with the WAMS control master station; the method specifically comprises the following steps: 1. establishing an isolated power grid system model; 2. designing control logics of the electrolytic aluminum load controller under different power disturbance modes; 3. and designing parameters of the electrolytic aluminum load controller under different power disturbance modes. The method utilizes the heat energy storage characteristic of the electrolytic aluminum load to realize effective power regulation of the electrolytic aluminum load and maintain the frequency stability of the system.

Description

Method for controlling frequency of isolated power grid by participation of electrolysis load
Technical Field
The invention belongs to the technical field of operation and control of power systems, and particularly relates to a method for controlling the frequency of an isolated power grid by electrolytic loads.
Background
The mode that renewable energy is consumed on the spot by forming an isolated power grid by using high energy consumption load, thermal power and renewable energy has attracted wide attention, however, due to the randomness and volatility of the renewable energy and the defects of small inertia, low reserve capacity and the like of the isolated power grid, along with the gradual increase of the proportion of the renewable energy in the isolated power grid, the frequency stability of the isolated power grid containing high-permeability wind power faces a great challenge.
For the problem of safe and stable operation of an isolated power grid containing high-permeability renewable energy, one method is to equip a local power grid with an energy storage device, however, for the local power grid containing centralized renewable energy, the capacity of the local power grid is often as high as million kilowatts, and the capacity of the conventional energy storage device is difficult to match.
The other method is to develop local load control capability to enable the load to participate in system frequency adjustment, and corresponding research is available for direct load control based on commercial loads such as HVAC, electric pumps and the like at home and abroad. However, the commercial load and the residential load are small in single capacity and are relatively dispersed, and it is difficult to perform large-scale aggregation to stabilize power fluctuation of the centralized renewable energy source.
The research of the participation of the electrolytic aluminum load in the system frequency control draws great attention, however, the existing control method only aims at the instant power disturbance of an isolated power grid, and the problem of continuous power disturbance of the isolated power grid is difficult to solve.
Disclosure of Invention
The invention aims to provide a control structure based on participation of electrolytic aluminum in frequency control, and designs control logics of a controller and a calculation method of control parameters under different power disturbance scenes.
In order to achieve the purpose, the invention adopts the technical scheme that: a method for controlling electrolytic load to participate in isolated power grid frequency realizes system frequency regulation through an electrolytic aluminum load control system, wherein the electrolytic aluminum load control system comprises a network control unit NCU, a dead zone unit, a proportional amplification unit, a WASHOUT unit and an amplitude limiting unit; the proportion amplification unit is connected with the WASHOUT unit in parallel, one end of the proportion amplification unit is sequentially connected with the dead zone unit and the network control unit NCU, and the other end of the proportion amplification unit is connected with the amplitude limiting unit; the network control unit NCU is connected with the isolated power grid system, and the amplitude limiting unit is connected with the WAMS control master station; the method specifically comprises the following steps:
step 1, establishing an isolated power grid system model;
step 2, designing control logics of the electrolytic aluminum load controller under different power disturbance modes;
and 3, designing parameters of the electrolytic aluminum load controller in different power disturbance modes.
In the method for participating in isolated grid frequency control of electrolysis loads, the implementation of step 1 comprises the following steps:
step 1.1, obtaining an electrolytic aluminum load power-voltage mathematical model by establishing an electrolytic aluminum load topological structure; the method comprises the following specific steps:
1) the electrolytic aluminum load topological structure is equivalent to that a back electromotive force E is connected with an equivalent resistor R in series, and is connected with an alternating current bus line through a rectifier bridge, a saturable reactor and an on-load tap changer;
2) active power P of electrolytic aluminium loadASLAnd a DC side voltage VBThe relationship is as follows:
Figure BDA0001375635400000021
3) electrolytic aluminum load DC side bus voltage VBBus voltage V on the AC side of loadAHIs as in formula (2);
Figure BDA0001375635400000022
(2) in the formula VSRThe voltage drop of the saturable reactor is shown, and k is the transformation ratio of the on-load tap changing transformer;
4) changing the voltage of the high voltage side of the load alternating current bus by changing the voltage of the generator terminal; obtaining a voltage change Δ V of the generator terminal according to a voltage sensitivity methodGVoltage delta V of high voltage side of AC bus with loadAHThe relationship of (a) is as in formula (3);
ΔVAH=Ksens·ΔVG(3)
(3) in the formula, KsensThe voltage sensitivity coefficient of the bus i to the bus j is obtained;
5) the generator terminal voltage variation delta V can be obtained by the formulas (2) and (3)GActive power variation delta P with loadASLThe corresponding relationship of (A) is:
ΔPASL=KASLΔVG(4)
(4) in the formula, KASLThe load comprehensive proportionality coefficient;
step 1.2, acquiring electrical parameters of an isolated power grid system, and establishing a mathematical model of the isolated power grid system, wherein the mathematical model comprises a thermal power unit excitation system, a thermal power unit speed regulator, a thermal power unit, a transmission line, a transformer, a double-fed wind driven generator and an electrolytic aluminum load;
obtaining parameters of an inertia link by using least square curve fitting, wherein the parameters are as shown in a formula (5);
ΔPASL(t)=-0.054+0.054·e-2.782t(5)
written in the form of a transfer function, as in equation (6);
Figure BDA0001375635400000031
(6) in the formula, KASLIs the load amplification factor of the electrolytic aluminum; t isASLThe time constant of active power of electrolytic aluminum load.
In the method for participating in isolated grid frequency control of electrolysis loads, the step 2 is realized by:
step 2.1, establishing an electrolytic aluminum load closed-loop control system;
the electrolytic aluminum load closed-loop control system comprises a conventional WAMS system, an additional WAMS control master station, a downlink channel and a network control unit NCU; the conventional WAMS system comprises a PMU substation, an upstream channel and a WAMS main station; the WAMS main station communicates with the PMU substation through a TCP/IP protocol to realize system data acquisition and state monitoring; the WAMS control main station communicates with the PMU substation by adopting a UDP protocol, and directly acquires the state quantity of the PMU substation; the WAMS controls the main station to calculate control parameters in real time and sends the control parameters to a network control unit NCU;
the network control unit NCU is arranged on the side of the generator excitation cabinet to realize the control of the active power of the electrolytic aluminum load; the output port of the NCU is accessed to the reference voltage superposition point of the generator excitation system, and the acquired frequency signal is used as a feedback signal to be input; when power disturbance occurs in the system, the frequency deviation acts on a reference voltage superposition point of an excitation system of the thermal power generating unit through a Network Control Unit (NCU) to change the excitation voltage, so that the voltage of the alternating current side of the electrolytic aluminum load is changed, and the active power of the electrolytic aluminum load is adjusted to stabilize the unbalanced power of the system;
step 2.2, designing the control logic of the electrolytic aluminum load controller;
I) monitoring system power disturbance according to PMU substation information, and judging the system state;
II) when the system state is instantaneous power disturbance or continuous power fluctuation, the WAMS controls the master station to issue control parameters;
III) collecting local frequency signals, and judging the communication state of the WAMS control master station;
IV) if the communication state of the WAMS control master station is normal, receiving the WAMS master station parameters, and updating the NCU parameters of the network control unit to control the excitation voltage of the generator;
v) if the communication state of the main station is abnormal, the excitation voltage of the generator is controlled by locally storing control parameters by the network control unit NCU.
In the method for participating in isolated grid frequency control of electrolysis loads, the step 3 is realized by:
step 3.1, designing parameters of a proportional amplification unit of the electrolytic aluminum load controller;
the proportion amplifying unit provides steady-state power support for the system, and the power disturbance quantity delta P generated by the system comprises the load adjustment quantity delta P of the electrolytic aluminumASLregPrimary regulating quantity delta P of thermal power generating unitGreg(ii) a The primary regulating quantity and the system frequency variable quantity of the thermal power generating unit satisfy the following formula:
Figure BDA0001375635400000041
(7) in the formula: f. ofNTo a rated frequency, fregIs a steady-state value of the system frequency after primary frequency modulation, PGNRated power of the generator set, RGThe difference adjustment coefficient of the thermal power generating unit is obtained;
the relationship between the electrolytic aluminum load adjustment amount and the system frequency variation satisfies the following formula (8):
Figure BDA0001375635400000051
(8) in the formula: Δ fASLLbFor the dead band range of the load controller, PASLNRated power for electrolytic aluminum load, RASLThe load equivalent adjustment coefficient of the electrolytic aluminum is obtained;
the power regulation amount borne by the electrolytic aluminum load is as follows:
ΔPASLreg=ΔP-ΔPGreg(9)
(9) in the formula, the system power disturbance quantity delta P is calculated in two situations;
a. under the condition of instantaneous power disturbance, the system power disturbance quantity is delta PSThen, there is formula (10):
ΔP=ΔPS(10)
b. counting the maximum fluctuation quantity delta p of the wind power for a period of time under continuous power fluctuationwindmax(ii) a The fluctuation amount is used as a system power disturbance amount delta P, and an amplification factor K is calculatedPThen, formula (11):
ΔP=ΔPwindmax(11)
the load equivalent adjustment coefficient of the electrolytic aluminum is obtained by the following equations (9), (10) and (11):
for an electrolytic aluminum load controller acting on a generator excitation system, equation (13) is satisfied:
when the electrolytic aluminum load active power changesThe chemical quantity is delta PASLregIn time, the variation of the terminal voltage of the thermal power generating unit is as follows:
Figure BDA0001375635400000054
the proportional amplification factor K is obtained from the formula (13) and the formula (14)P:
Figure BDA0001375635400000055
Step 3.2, designing parameters of a WASHOUT unit of the electrolytic aluminum load controller;
calculating WASHOUT unit parameters through an SFR model;
WASHOUT unit parameter TDFor high-pass filtering, generally selecting 6-7 s;
Δ f(s) and Δ PLThe relationship of(s) is:
Figure BDA0001375635400000061
(16) in the formula:
Figure BDA0001375635400000062
Figure BDA0001375635400000063
(17) in the formulas (18), M is the equivalent inertia coefficient of the system, D is the load damping coefficient, and TD、TR、FH、KMThe time constant of the generator is R is the droop coefficient of the speed regulator of the generator;
ignoring the pole-zero far from the imaginary axis, equation (16) can be converted to:
Figure BDA0001375635400000064
(19) in the formula, the parameter z2,ζ,ωndAll of which can pass through SFR modelCalculating parameters to obtain;
the lagrange transform of the output signal of the SFR model under the unit step function is obtained from equation (19):
Figure BDA0001375635400000065
in formula (20):
Figure BDA0001375635400000066
performing reverse Laplace transformation on the formula (20) to obtain
(21) In the formula:
let t bezAt this time, the system frequency offset is maximum, and the slope of the frequency change amount at this time is zero, which can be expressed by equation (22):
Figure BDA0001375635400000071
the maximum offset occurrence time t of the system frequency can be obtained by the formula (22)zBy adjusting WASHOUT unit KDParameter, maximum frequency offset of system Δ fdevmaxIs maintained at Δ fmaxAccordingly, formula (23) can be obtained:
the WASHOUT unit coefficient K can be obtained from the expressions (22) - (23)D
The invention has the beneficial effects that:
1. the load controller fully considers the complex working conditions of the actual site, and the industrial application of the load controller can be realized.
2. The method can effectively identify the instantaneous power disturbance and the continuous power disturbance of an isolated power grid system, automatically use different control modes for different disturbances, realize effective power regulation of electrolytic aluminum load and maintain the stable frequency of the system.
3. On-site tests have been carried out on the industrial site, and the test results prove the effectiveness of the control method of the invention.
Drawings
FIG. 1 is a mathematical model of the load of electrolytic aluminum created according to one embodiment of the present invention in conjunction with the principles of the electrolytic aluminum process;
FIG. 2 is a dynamic power fitting curve of the load of the electrolytic aluminum, which is established by the measured data according to one embodiment of the invention;
FIG. 3 is a diagram of the topology of an electrolytic aluminum load controller in accordance with one embodiment of the present invention;
FIG. 4 is a control logic diagram of an electrolytic aluminum load controller in accordance with an embodiment of the present invention;
FIG. 5(a) is a wind power graph under a closed loop control test considering a rapid drop in wind power according to an embodiment of the present invention;
FIG. 5(b) is a system frequency variation graph under a closed-loop control test considering a rapid drop in wind power according to an embodiment of the present invention;
FIG. 6(a) is a graph of voltage variation of generator terminal under a closed-loop control test considering rapid decrease of wind power according to an embodiment of the present invention;
FIG. 6(b) is a graph of the voltage variation of the load bus under a closed-loop control test considering the rapid decrease of the wind power according to an embodiment of the present invention;
FIG. 7(a) is a graph of active power change of electrolytic aluminum in a closed loop control test considering rapid wind power droop in accordance with an embodiment of the present invention;
FIG. 7(b) is a graph of the active power change of the generator under a closed-loop control test considering the rapid wind power drop according to an embodiment of the present invention;
FIG. 8 is a graph of wind power change under a closed loop control test with consideration of wind farm removal in accordance with an embodiment of the present invention;
FIG. 9(a) is a graph of the frequency change of a loaded controller system under a closed loop control test with consideration of wind farm shedding for one embodiment of the present invention;
FIG. 9(b) is a graph of the frequency change of a system without a load controller under a closed loop control test with consideration of wind farm shedding in accordance with an embodiment of the present invention;
FIG. 10(a) is a graph of active power change for an electrolytic aluminum load with a load controller under a closed loop control test with consideration of wind farm removal in accordance with an embodiment of the present invention;
FIG. 10(b) is a graph of the change in active power without the load controller electrolytic aluminum load under a closed loop control test with consideration of wind farm removal, in accordance with an embodiment of the present invention.
Detailed Description
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
The embodiment is realized by the following technical scheme that the method for controlling the frequency of the isolated power grid by the electrolytic load is characterized in that the system frequency is adjusted by an electrolytic aluminum load control system, wherein the electrolytic aluminum load control system comprises a network control unit NCU, a dead zone unit, a proportional amplification unit, a WASHOUT unit and an amplitude limiting unit; the proportion amplification unit is connected with the WASHOUT unit in parallel, one end of the proportion amplification unit is sequentially connected with the dead zone unit and the network control unit NCU, and the other end of the proportion amplification unit is connected with the amplitude limiting unit; the network control unit NCU is connected with the isolated power grid system, and the amplitude limiting unit is connected with the WAMS control master station; the method specifically comprises the following steps:
step 1, establishing an isolated power grid system model;
step 2, designing control logics of the electrolytic aluminum load controller under different power disturbance modes;
and 3, designing parameters of the electrolytic aluminum load controller in different power disturbance modes.
Further, the implementation of step 1 includes:
step 1.1, obtaining an electrolytic aluminum load power-voltage mathematical model by establishing an electrolytic aluminum load topological structure; the method comprises the following specific steps:
1) the electrolytic aluminum load topological structure is equivalent to that a back electromotive force E is connected with an equivalent resistor R in series, and is connected with an alternating current bus line through a rectifier bridge, a saturable reactor and an on-load tap changer;
2) active power P of electrolytic aluminium loadASLAnd a DC side voltage VBThe relationship is as follows:
Figure BDA0001375635400000091
3) electrolytic aluminum load DC side bus voltage VBBus voltage V on the AC side of loadAHIs as in formula (2);
Figure BDA0001375635400000092
(2) in the formula VSRThe voltage drop of the saturable reactor is shown, and k is the transformation ratio of the on-load tap changing transformer;
4) changing the voltage of the high voltage side of the load alternating current bus by changing the voltage of the generator terminal; obtaining a voltage change Δ V of the generator terminal according to a voltage sensitivity methodGVoltage delta V of high voltage side of AC bus with loadAHThe relationship of (a) is as in formula (3);
ΔVAH=Ksens·ΔVG(3)
(3) in the formula, KsensThe voltage sensitivity coefficient of the bus i to the bus j is obtained;
5) the generator terminal voltage variation delta V can be obtained by the formulas (2) and (3)GActive power variation delta P with loadASLThe corresponding relationship of (A) is:
ΔPASL=KASLΔVG(4)
(4) in the formula, KASLThe load comprehensive proportionality coefficient;
step 1.2, acquiring electrical parameters of an isolated power grid system, and establishing a mathematical model of the isolated power grid system, wherein the mathematical model comprises a thermal power unit excitation system, a thermal power unit speed regulator, a thermal power unit, a transmission line, a transformer, a double-fed wind driven generator and an electrolytic aluminum load;
obtaining parameters of an inertia link by using least square curve fitting, wherein the parameters are as shown in a formula (5);
ΔPASL(t)=-0.054+0.054·e-2.782t(5)
written in the form of a transfer function, as in equation (6);
Figure BDA0001375635400000101
(6) in the formula, KASLIs the load amplification factor of the electrolytic aluminum; t isASLThe time constant of active power of electrolytic aluminum load.
Further, the implementation of step 2 includes:
step 2.1, establishing an electrolytic aluminum load closed-loop control system;
the electrolytic aluminum load closed-loop control system comprises a conventional WAMS system, an additional WAMS control master station, a downlink channel and a network control unit NCU; the conventional WAMS system comprises a PMU substation, an upstream channel and a WAMS main station; the WAMS main station communicates with the PMU substation through a TCP/IP protocol to realize system data acquisition and state monitoring; the WAMS control main station communicates with the PMU substation by adopting a UDP protocol, and directly acquires the state quantity of the PMU substation; the WAMS controls the main station to calculate control parameters in real time and sends the control parameters to a network control unit NCU;
the network control unit NCU is arranged on the side of the generator excitation cabinet to realize the control of the active power of the electrolytic aluminum load; the output port of the NCU is accessed to the reference voltage superposition point of the generator excitation system, and the acquired frequency signal is used as a feedback signal to be input; when power disturbance occurs in the system, the frequency deviation acts on a reference voltage superposition point of an excitation system of the thermal power generating unit through a Network Control Unit (NCU) to change the excitation voltage, so that the voltage of the alternating current side of the electrolytic aluminum load is changed, and the active power of the electrolytic aluminum load is adjusted to stabilize the unbalanced power of the system;
step 2.2, designing the control logic of the electrolytic aluminum load controller;
I) monitoring system power disturbance according to PMU substation information, and judging the system state;
II) when the system state is instantaneous power disturbance or continuous power fluctuation, the WAMS controls the master station to issue control parameters;
III) collecting local frequency signals, and judging the communication state of the WAMS control master station;
IV) if the communication state of the WAMS control master station is normal, receiving the WAMS master station parameters, and updating the NCU parameters of the network control unit to control the excitation voltage of the generator;
v) if the communication state of the main station is abnormal, the excitation voltage of the generator is controlled by locally storing control parameters by the network control unit NCU.
Further, the implementation of step 3 includes:
step 3.1, designing parameters of a proportional amplification unit of the electrolytic aluminum load controller;
the proportion amplifying unit provides steady-state power support for the system, and the power disturbance quantity delta P generated by the system comprises the load adjustment quantity delta P of the electrolytic aluminumASLregPrimary regulating quantity delta P of thermal power generating unitGreg(ii) a The primary regulating quantity and the system frequency variable quantity of the thermal power generating unit satisfy the following formula:
Figure BDA0001375635400000121
(7) in the formula: f. ofNTo a rated frequency, fregIs a steady-state value of the system frequency after primary frequency modulation, PGNRated power of the generator set, RGThe difference adjustment coefficient of the thermal power generating unit is obtained;
the relationship between the electrolytic aluminum load adjustment amount and the system frequency variation satisfies the following formula (8):
Figure BDA0001375635400000122
(8) in the formula: Δ fASLLbFor the dead band range of the load controller, PASLNRated power for electrolytic aluminum load, RASLThe load equivalent adjustment coefficient of the electrolytic aluminum is obtained;
the power regulation amount borne by the electrolytic aluminum load is as follows:
ΔPASLreg=ΔP-ΔPGreg(9)
(9) in the formula, the system power disturbance quantity delta P is calculated in two situations;
a. under the condition of instantaneous power disturbance, the system power disturbance quantity is delta PSThen, there is formula (10):
ΔP=ΔPS(10)
b. under continuous power fluctuation, counting the maximum fluctuation quantity delta pwindmax of the wind power for a period of time; the fluctuation amount is used as a system power disturbance amount delta P, and an amplification factor K is calculatedPThen, formula (11):
ΔP=ΔPwindmax(11)
the load equivalent adjustment coefficient of the electrolytic aluminum is obtained by the following equations (9), (10) and (11):
Figure BDA0001375635400000123
for an electrolytic aluminum load controller acting on a generator excitation system, equation (13) is satisfied:
when the active power variation of the electrolytic aluminum load is delta PASLregIn time, the variation of the terminal voltage of the thermal power generating unit is as follows:
Figure BDA0001375635400000125
the proportional amplification factor K is obtained from the formula (13) and the formula (14)P:
Figure BDA0001375635400000131
Step 3.2, designing parameters of a WASHOUT unit of the electrolytic aluminum load controller;
calculating WASHOUT unit parameters through an SFR model;
WASHOUT unit parameter TDFor high-pass filtering, generally selecting 6-7 s;
Δ f(s) and Δ PLThe relationship of(s) is:
Figure BDA0001375635400000132
(16) in the formula:
Figure BDA0001375635400000133
Figure BDA0001375635400000134
(17) in the formulas (18), M is the equivalent inertia coefficient of the system, D is the load damping coefficient, and TD、TR、FH、KMThe time constant of the generator is R is the droop coefficient of the speed regulator of the generator;
ignoring the pole-zero far from the imaginary axis, equation (16) can be converted to:
(19) in the formula, the parameter z2,ζ,ωndAll can be obtained by the parameter calculation of the SFR model;
the lagrange transform of the output signal of the SFR model under the unit step function is obtained from equation (19):
Figure BDA0001375635400000136
in formula (20):
Figure BDA0001375635400000137
performing reverse Laplace transformation on the formula (20) to obtain
Figure BDA0001375635400000141
(21) In the formula:
let t bezAt this time, the system frequency offset is maximum, and the slope of the frequency change amount at this time is zero, which can be expressed by equation (22):
Figure BDA0001375635400000143
the maximum offset occurrence time t of the system frequency can be obtained by the formula (22)zBy adjusting WASHOUT unit KDParameter, maximum frequency offset of system Δ fdevmaxIs maintained at Δ fmaxAccordingly, formula (23) can be obtained:
Figure BDA0001375635400000144
the WASHOUT unit coefficient K can be obtained from the expressions (22) - (23)D
In specific implementation, the method for controlling the frequency of the isolated power grid by the electrolytic load comprises the following steps:
s1, modeling an isolated power grid system: obtaining an electrolytic aluminum load topological structure, and establishing an electrolytic aluminum load power-voltage mathematical model; obtaining electrical parameters of an isolated power grid, establishing a mathematical model of an isolated power grid system, wherein the mathematical model comprises a thermal power unit excitation system, a thermal power unit speed regulator, a thermal power unit, a transmission line, a transformer, a double-fed wind driven generator and an electrolytic aluminum load, and establishing a simulation model of the isolated power grid system in a real-time digital simulation system (RTDS).
S2, designing a topological structure of the electrolytic aluminum load controller. And designing a topological structure of the controller according to the transient process power support and the steady process power support by combining the wind power fluctuation characteristic of the isolated power grid. The topological structure of the controller is divided into a proportion amplifying part, and power support is mainly provided in a steady-state process; the wasbout portion, mainly provides power support during transients.
And S3, designing action logics of the electrolytic aluminum load controller under different disturbance conditions. The electrolytic aluminum load controller can automatically identify the power disturbance of the system, and the power disturbance is divided into instant power disturbance and continuous power disturbance. Under different power disturbance modes, the control parameters of the electrolytic aluminum load controller are different. Meanwhile, the electrolytic aluminum load controller has the functions of local control and online control. When the communication is normal, the electrolytic aluminum load controller adopts an online control function; when the communication is abnormal, the electrolytic aluminum load controller adopts a local control mode.
And S4, designing parameters of the electrolytic aluminum load controller under different power disturbance modes. When the system has instantaneous power disturbance, the proportion amplifying part of the electrolytic aluminum load controller and the WASHOUT part act together, and the adjusting power born by the thermal power generating unit and the electrolytic aluminum load controller is calculated according to the unbalanced power generated by the wide area measurement system. And then calculating a parameter KP of the proportional amplification part. And calculating the parameters of the WASHOUT part according to the maximum offset of the system frequency greater than a certain limit value. When the system has continuous power fluctuation, the parameter of the proportional amplification part is calculated according to historical wind power data, and the WASHOUT part does not act in the mode.
And S5, respectively carrying out comparison tests on the designed electrolytic aluminum load controller on an RTDS simulation platform and an industrial field, and checking the control effect of the load controller under the instantaneous disturbance of the wind power and the continuous fluctuation of the wind power.
The following is further explained in conjunction with the accompanying drawings.
Modeling of electrolytic aluminum load
The electrolytic aluminum load is a heat storage type load, the power regulation speed is high, and dynamic power support can be provided for an isolated power grid system. The electrolytic aluminum load can be equivalent to that the back electromotive force E is connected with the equivalent resistor R in series, and is connected with an alternating current bus line through a rectifier bridge, a saturable reactor and an on-load tap changing transformer, and an equivalent circuit diagram is shown in figure 1.
Active power P of electrolytic aluminium loadASLAnd a DC side voltage VBThe relationship is as follows:
Figure BDA0001375635400000151
electrolytic aluminum load DC side bus voltage VBBus voltage on the AC side of loadVAHIs as in formula (2):
Figure BDA0001375635400000161
in the formula VSRAnd k is the voltage drop of the saturable reactor, and k is the transformation ratio of the on-load tap changing transformer.
Because each bus of the isolated power grid is small in electrical distance, the voltage of the high-voltage side of the load alternating-current bus can be changed by changing the voltage of the generator terminal. According to the voltage sensitivity method, the voltage change quantity delta V of the generator terminal can be obtainedGVoltage delta V of high voltage side of AC bus with loadAHThe relationship is as in formula (3).
ΔVAH=Ksens·ΔVG(3)
Wherein, KsensIs the voltage sensitivity coefficient of bus i to bus j.
The generator terminal voltage variation Δ V can be obtained by the equations (2) and (3)GActive power variation delta P with loadASLThe corresponding relationship of (1).
ΔPASL=KASLΔVG(4)
Wherein, KASLThe load comprehensive proportionality coefficient is obtained.
In FIG. 2, the thick line is the change curve of the active power of the electrolytic aluminum load when the voltage of the generator terminal measured on site has a-0.15 p.u. step. As can be seen from the figure, the dynamic response of the electrolytic aluminum load has an obvious inertia link. A first-order inertia link is adopted to describe the dynamic response of the electrolytic aluminum load, and each parameter of the inertia link is obtained by using least square curve fitting, as shown in formula (5), and a thin line in figure 2 is a fitting curve.
ΔPASL(t)=-0.054+0.054·e-2.782t(5)
Written in the form of a transfer function, as shown in equation (6).
Figure BDA0001375635400000162
In the formula, KASLFor load amplification of electrolytic aluminiumA coefficient; t isASLThe time constant of active power of electrolytic aluminum load.
Topological structure design of electrolytic aluminum load controller
The long-term frequent participation of the electrolytic aluminum load in the system frequency adjustment has certain influence on the production efficiency of the electrolytic aluminum, and the control range of the electrolytic aluminum load needs to be limited. In addition, as can be seen from fig. 2, the electrolytic aluminum load dynamic response is faster than the primary frequency modulation of the thermal power generating unit, and more power support should be provided during the transient of the frequency modulation. The logic block diagram of the electrolytic aluminum load controller is shown in figure 3 according to the dynamic characteristics of the electrolytic aluminum load. The load controller comprises a proportional amplification unit, a WASHOUT unit, a dead zone unit and a limiting unit.
The dead zone unit realizes the locking of the electrolytic aluminum load controller, when the system is in a normal operation condition, the frequency deviation of the system is small, and the dead zone unit locking load controller is adjusted only by primary frequency modulation of the thermal power generating unit.
The WASHOUT unit presents a high-pass filtering characteristic and provides transient power support only when the system has an emergency condition, which causes the frequency to change violently. The wasbout cell effect diminishes as the frequency levels back off. Controlling parameter K by WASHOUT unitDThe maximum variation of the system frequency can be reduced, and the system frequency is prevented from being reduced too low.
The function of the proportional amplification unit is similar to that of a primary frequency modulation coefficient R of the thermal power generating unit, and the proportional amplification unit and the primary frequency modulation coefficient R of the thermal power generating unit jointly provide steady-state power support for the system. By controlling the proportional amplification unit, the primary frequency modulation spare capacity of the system can be improved.
Control logic of electrolytic aluminum load controller
The following is an electrolytic aluminum load closed loop control system applied to the site. The closed-loop control system is additionally provided with a WAMS control main station, a downlink channel and a network control unit NCU on the basis of a conventional WAMS system (comprising a PMU substation, an uplink channel and a WAMS main station).
The conventional WAMS main station and the PMU realize the functions of system data acquisition and state monitoring through TCP/IP protocol communication. The WAMS control main station communicates with the PMU by adopting a UDP protocol, and directly obtains the state quantity of the PMU substation. The WAMS control main station comprises a high-level application program, calculates control parameters in real time according to the system state and sends the control parameters to the network control unit NCU.
And the network control unit NCU is arranged on the side of the generator excitation cabinet to realize the control of the active power of the electrolytic aluminum load. The network control unit NCU has the functions of frequency signal acquisition and real-time communication with the WAMS control master station, has independent logical operation capacity and can convert digital quantity into physical quantity for output. And the output port of the NCU is accessed to the reference voltage superposition point of the generator excitation system, and the acquired frequency signal is used as a feedback signal to be input. When power disturbance occurs in the system, the frequency deviation acts on a reference voltage superposition point of an excitation system of the thermal power generating unit through the network control unit NCU, so that the excitation voltage is changed, the voltage of the alternating current side of the electrolytic aluminum load is changed, and the active power of the electrolytic aluminum load is adjusted to stabilize the unbalanced power of the system.
In order to improve the reliability of the closed-loop control system, a control mode which takes local control as a main mode and takes on-line control as an auxiliary mode is adopted, and a local device network control unit NCU comprises a main control logic and can still independently complete local closed-loop control under the condition of communication interruption with a WAMS control main station. And the WAMS control master station modifies the control parameters of the NCU of the network control unit on line according to the running state of the system, and optimizes the effect of the closed-loop control system. The control logic of the real-time closed-loop control system based on the electrolytic aluminum load is shown in FIG. 4.
Fourth, design of parameters of electrolytic aluminum load controller
4.1. Parameter design of proportional amplification unit of electrolytic aluminum load controller
The scale up unit provides steady state power support for the system. The power disturbance quantity delta P generated by the system is adjusted by the load of the electrolytic aluminumASLregPrimary regulation delta P of thermal power generating unitGregAnd (4) sharing the common charge. The primary regulating quantity and the system frequency variable quantity of the thermal power generating unit satisfy the following formula:
wherein: f. ofNTo a rated frequency, fregIs a steady-state value of the system frequency after primary frequency modulation, PGNRated power of the generator set, RGAnd the difference adjustment coefficient is the thermal power generating unit.
The relationship between the electrolytic aluminum load adjustment amount and the system frequency variation satisfies the following formula (8):
Figure BDA0001375635400000191
wherein: Δ fASLLbFor the dead band range of the load controller, PASLNRated power for electrolytic aluminum load, RASLThe load equivalent adjustment coefficient of the electrolytic aluminum is shown.
The power regulation amount borne by the electrolytic aluminum load is shown as the formula (9):
ΔPASLreg=ΔP-ΔPGreg(9)
the system power disturbance amount Δ P is calculated in two cases.
i) Under the condition of instantaneous power disturbance, the real-time monitoring system can be used for monitoring the power disturbance generated by the system to directly obtain the system power disturbance quantity delta P, and the system power disturbance quantity is recorded as delta P under the conditionSThen, there is formula (10):
ΔP=ΔPS(10)
ii) under the continuous fluctuation of the wind power, the fluctuation quantity of the wind power is difficult to measure and calculate in real time, and the maximum fluctuation quantity delta Pwindmax of the wind power is counted for a period of time on the basis of historical change data of the wind power. The fluctuation amount is used as a system power disturbance amount delta P, and an amplification factor K is calculatedPThen, formula (11):
ΔP=ΔPwindmax(11)
the load equivalent adjustment coefficient of the electrolytic aluminum is obtained by the following equations (9), (10) and (11):
Figure BDA0001375635400000192
for an electrolytic aluminum load controller acting on a generator excitation system, equation (13) is satisfied:
Figure BDA0001375635400000193
when the active power variation of the electrolytic aluminum load is delta PASLregIn time, the variation of the terminal voltage of the thermal power generating unit is as follows:
Figure BDA0001375635400000201
the proportional amplification factor K is obtained from the formula (13) and the formula (14)P:
Figure BDA0001375635400000202
4.2. Design of WASHOUT unit of electrolytic aluminum load control system
The wasbouut unit provides transient power support only when the system frequency changes rapidly, reducing the maximum frequency offset. The WASHOUT unit parameters may be calculated by an SFR model.
Assume a coherent response of all generators to system load changes and equate them to one machine. Further, for such isolated power systems, automatic generation control is not invested, and therefore the role of the secondary frequency modulation will not be considered in the SFR model. The prime movers of all generators in the isolated grid under study are reheat turbines. Due to the fact that the response time of the thermal power unit speed regulator is short, the SFR model ignores the dynamic response process of the thermal power unit speed regulator. WASHOUT unit parameter TDSelecting 6-7s for high-pass filtering according to engineering experience;
Δ f(s) and Δ PLThe relationship of(s) is shown by the following formula:
Figure BDA0001375635400000203
wherein:
Figure BDA0001375635400000204
Figure BDA0001375635400000205
wherein M is the equivalent inertia coefficient of the system, D is the load damping coefficient, and TD、TR、FH、KMAnd R is the droop coefficient of the speed regulator of the generator.
Since the system has a pole-zero far from the imaginary axis, neglecting the pole-zero far from the imaginary axis, equation (16) can be expressed as follows:
Figure BDA0001375635400000211
in the formula, the parameter z2,ζ,ωndAll can be obtained by the parameter calculation of the SFR model.
The lagrange transform of the output signal of the SFR model under the unit step function is obtained from equation (19):
Figure BDA0001375635400000212
wherein:
Figure BDA0001375635400000213
performing reverse Laplace transformation on the formula (20) to obtain
Figure BDA0001375635400000214
Wherein:
let t bezAt this time, the system frequency offset is maximum, and the slope of the frequency change amount at this time is zero, which can be expressed by equation (22):
Figure BDA0001375635400000216
the maximum offset occurrence time t of the system frequency can be obtained by the formula (22)zBy adjusting WASHOUT unit KDParameter, maximum frequency offset of system Δ fdevmaxIs maintained at Δ fmaxAccordingly, formula (23) can be obtained:
Figure BDA0001375635400000217
the WASHOUT unit coefficient K can be obtained from the expressions (22) - (23)D
Five examples and simulations
The following research is conducted based on a certain actual isolated grid in China.
Example 1: considering a closed-loop control test of the rapid reduction of the wind power;
the wind farm output at the initial moment was 61 MW. The power of the wind power plant is manually intervened, the active output of the wind power plant is reduced by 30MW within 30 seconds, and the active output curve of the wind power plant is shown in fig. 5 (a). When the closed loop control system is engaged, the system frequency curve is shown in fig. 5 (b). The generator terminal voltage change curve is shown by a solid line in FIG. 6(a), and the electrolytic aluminum load bus voltage change curve is shown by a solid line in FIG. 6 (b). The active power of the electrolytic aluminum load is shown by a solid line in FIG. 7 (a). The active power of the thermal power generating unit is shown by a solid line in fig. 7 (b). Because the thermal power generating unit and the electrolytic aluminum load jointly stabilize the disturbance of wind power, the system frequency is finally stabilized at 49.81Hz, and the system keeps stable operation.
As a comparative test, when the closed-loop control system is not put into operation, the response curves of the system frequency, the generator terminal voltage, the load bus voltage, the load active power and the thermal power unit active power are shown by dotted lines in fig. 5(b), 6(a) (b), 7(a) (b). Comparing the frequency response curves represented by the solid line and the dotted line in fig. 5(b), it can be seen that the proposed control method can effectively control the electrolytic aluminum load to participate in the frequency control of the isolated power grid, share the primary frequency modulation pressure of the thermal power generating unit, and significantly improve the frequency response characteristics of the isolated power grid.
Example 2: considering a closed-loop control test for cutting off a wind power plant;
in order to verify the effect of participating in system frequency control when large disturbance of a control system occurs, a field test for cutting off a wind power plant is carried out. Considering the danger of the field test, the comparative test that the electrolytic aluminum load does not participate in the frequency control adopts RTDS-based hardware-in-loop simulation. And cutting off the active power output of the wind power plant before cutting off, and cutting off the wind power plant when t is 5s, as shown in FIG. 8. The sudden change of the wind power causes the system frequency to drop rapidly, and under the action of the closed-loop control system, the electrolytic aluminum load can actively and rapidly adjust the drop of the automatic response system frequency, so that active support is provided for the system, the active power of the system drops by 30MW within 3s, and sufficient active support can be provided in the transient process, as shown by the solid curve of fig. 10 (a). And then, the primary frequency modulation of the thermal power generating unit also plays a role gradually, part of power adjustment amount is born, and finally, under the common response of the thermal power generating unit and the electrolytic aluminum load, the system frequency is stabilized at about 49.9 Hz. And simultaneously, simulation tests under the same working condition are carried out on the RTDS-based hardware-in-loop simulation platform. The system frequency response curves of the hardware-in-loop platform simulation at the time of the control method input and output are respectively shown by dotted lines in fig. 9(a) and (b). As can be seen from fig. 9(a), under the condition that the boundary conditions and the fault settings are the same, the frequency response curve of the system obtained by simulation can more accurately simulate the frequency response curve actually measured on site. Without a closed loop control system, the electrolytic aluminum load power would not respond to the frequency change of the system, as shown by the dashed curve in FIG. 9 (b). Although the frequency of the system does not collapse under the action of the primary frequency modulation of the thermal power generating unit, compared with the situation that the control method is put into use, the frequency of the system is greatly deviated, and the frequency is reduced to 49.5Hz at the lowest, as shown in FIG. 9 (b). For large power disturbances, a too low frequency drop may cause the low frequency load shedding device to act.
In summary, the control method of the embodiment can provide reference for stability research and control strategy formulation and verification of a system similar to the isolated power grid.
It should be understood that parts of the specification not set forth in detail are well within the prior art.
Although specific embodiments of the present invention have been described above with reference to the accompanying drawings, it will be appreciated by those skilled in the art that these are merely illustrative and that various changes or modifications may be made to these embodiments without departing from the principles and spirit of the invention. The scope of the invention is only limited by the appended claims.

Claims (2)

1. A method for enabling electrolysis loads to participate in frequency control of an isolated power grid is characterized in that system frequency adjustment is achieved through an electrolytic aluminum load control system, and the electrolytic aluminum load control system comprises a network control unit NCU, a dead zone unit, a proportional amplification unit, a WASHOUT unit and an amplitude limiting unit; the proportion amplification unit is connected with the WASHOUT unit in parallel, one end of the proportion amplification unit is sequentially connected with the dead zone unit and the network control unit NCU, and the other end of the proportion amplification unit is connected with the amplitude limiting unit; the network control unit NCU is connected with the isolated power grid system, and the amplitude limiting unit is connected with the WAMS control master station; the method specifically comprises the following steps:
step 1, establishing an isolated power grid system model;
step 2, designing control logics of the electrolytic aluminum load controller under different power disturbance modes;
step 3, designing parameters of the electrolytic aluminum load controller in different power disturbance modes;
the implementation of step 1 comprises:
step 1.1, obtaining an electrolytic aluminum load power-voltage mathematical model by establishing an electrolytic aluminum load topological structure; the method comprises the following specific steps:
1) the electrolytic aluminum load topological structure is equivalent to that a back electromotive force E is connected with an equivalent resistor R in series, and is connected with an alternating current bus line through a rectifier bridge, a saturable reactor and an on-load tap changer;
2) active power P of electrolytic aluminium loadASLAnd a DC side voltage VBThe relationship is as follows:
3) voltage V on DC side of electrolytic aluminium loadBBus voltage on the AC side of loadVAHIs as in formula (2);
Figure FDA0002209079190000012
(2) in the formula VSRThe voltage drop of the saturable reactor is shown, and k is the transformation ratio of the on-load tap changing transformer;
4) changing the voltage of the high voltage side of the load alternating current bus by changing the voltage of the generator terminal; obtaining a voltage change Δ V of the generator terminal according to a voltage sensitivity methodGVoltage delta V of high voltage side of AC bus with loadAHThe relationship of (a) is as in formula (3);
ΔVAH=Ksens·ΔVG(3)
(3) in the formula, KsensThe voltage sensitivity coefficient of the bus i to the bus j is obtained;
5) the generator terminal voltage variation delta V can be obtained by the formulas (2) and (3)GActive power variation delta P with loadASLThe corresponding relationship of (A) is:
ΔPASL=KASLΔVG(4)
(4) in the formula, KASLThe load comprehensive proportionality coefficient;
step 1.2, acquiring electrical parameters of an isolated power grid system, and establishing a mathematical model of the isolated power grid system, wherein the mathematical model comprises a thermal power unit excitation system, a thermal power unit speed regulator, a thermal power unit, a transmission line, a transformer, a double-fed wind driven generator and an electrolytic aluminum load;
obtaining parameters of an inertia link by using least square curve fitting, wherein the parameters are as shown in a formula (5);
ΔPASL(t)=-0.054+0.054·e-2.782t(5)
written in the form of a transfer function, as in equation (6);
Figure FDA0002209079190000021
(6) in the formula, KASLThe load comprehensive proportionality coefficient; t isASLIs the time constant of active power of electrolytic aluminum load;
the implementation of step 2 comprises:
step 2.1, establishing an electrolytic aluminum load closed-loop control system;
the electrolytic aluminum load closed-loop control system comprises a conventional WAMS system, an additional WAMS control master station, a downlink channel and a network control unit NCU; the conventional WAMS system comprises a PMU substation, an upstream channel and a WAMS main station; the WAMS main station communicates with the PMU substation through a TCP/IP protocol to realize system data acquisition and state monitoring; the WAMS control main station communicates with the PMU substation by adopting a UDP protocol, and directly acquires the state quantity of the PMU substation; the WAMS controls the main station to calculate control parameters in real time and sends the control parameters to a network control unit NCU;
the network control unit NCU is arranged on the side of the generator excitation cabinet to realize the control of the active power of the electrolytic aluminum load; the output port of the NCU is accessed to the reference voltage superposition point of the generator excitation system, and the acquired frequency signal is used as a feedback signal to be input; when power disturbance occurs in the system, the frequency deviation acts on a reference voltage superposition point of an excitation system of the thermal power generating unit through a Network Control Unit (NCU) to change the excitation voltage, so that the voltage of the alternating current side of the electrolytic aluminum load is changed, and the active power of the electrolytic aluminum load is adjusted to stabilize the unbalanced power of the system;
step 2.2, designing the control logic of the electrolytic aluminum load controller;
I) monitoring system power disturbance according to PMU substation information, and judging the system state;
II) when the system state is instantaneous power disturbance or continuous power fluctuation, the WAMS controls the master station to issue control parameters;
III) collecting local frequency signals, and judging the communication state of the WAMS control master station;
IV) if the communication state of the WAMS control master station is normal, receiving the WAMS master station parameters, and updating the NCU parameters of the network control unit to control the excitation voltage of the generator;
v) if the communication state of the main station is abnormal, the excitation voltage of the generator is controlled by locally storing control parameters by the network control unit NCU.
2. The method for participating in isolated grid frequency control of electrolysis-type loads according to claim 1, wherein the step 3 is implemented by:
step 3.1, designing parameters of a proportional amplification unit of the electrolytic aluminum load controller;
the proportion amplifying unit provides steady-state power support for the system, and the power disturbance quantity delta P generated by the system comprises the load adjustment quantity delta P of the electrolytic aluminumASLregPrimary regulating quantity delta P of thermal power generating unitGreg(ii) a The primary regulating quantity and the system frequency variable quantity of the thermal power generating unit satisfy the following formula:
(7) in the formula: f. ofNTo a rated frequency, fregIs a steady-state value of the system frequency after primary frequency modulation, PGNRated power of the generator set, RGThe difference adjustment coefficient of the thermal power generating unit is obtained;
the relationship between the electrolytic aluminum load adjustment amount and the system frequency variation satisfies the following formula (8):
Figure FDA0002209079190000041
(8) in the formula: Δ fASLLbFor the dead band range of the load controller, PASLNRated power for electrolytic aluminum load, RASLThe load equivalent adjustment coefficient of the electrolytic aluminum is obtained;
the power regulation amount borne by the electrolytic aluminum load is as follows:
ΔPASLreg=ΔP-ΔPGreg(9)
(9) in the formula, the system power disturbance quantity delta P is calculated in two situations;
a. under the condition of instantaneous power disturbance, the system power disturbance quantity is delta PSThen, there is formula (10):
ΔP=ΔPS(10);
b. counting the maximum fluctuation quantity delta p of the wind power for a period of time under continuous power fluctuationwindmax(ii) a The fluctuation amount is used as a system power disturbance amount delta P, and an amplification factor K is calculatedPThen, formula (11):
ΔP=ΔPwindmax(11);
the load equivalent adjustment coefficient of the electrolytic aluminum is obtained by the following equations (9), (10) and (11):
Figure FDA0002209079190000042
for an electrolytic aluminum load controller acting on a generator excitation system, equation (13) is satisfied:
Figure FDA0002209079190000043
when the active power variation of the electrolytic aluminum load is delta PASLregIn time, the variation of the terminal voltage of the thermal power generating unit is as follows:
Figure FDA0002209079190000044
the proportional amplification factor K is obtained from the formula (13) and the formula (14)P:
Step 3.2, designing parameters of a WASHOUT unit of the electrolytic aluminum load controller;
calculating WASHOUT unit parameters through an SFR model;
WASHOUT unit parameter TDFor high-pass filtering, generally selecting 6-7 s;
Δ f(s) and Δ PLThe relationship of(s) is:
Figure FDA0002209079190000051
(16) in the formula:
Figure FDA0002209079190000052
Figure FDA0002209079190000053
(17) in the formulas (18), M is the equivalent inertia coefficient of the system, D is the load damping coefficient, and TD、TR、FH、KMThe time constant of the generator is R is the droop coefficient of the speed regulator of the generator;
ignoring the pole-zero far from the imaginary axis, equation (16) can be converted to:
Figure FDA0002209079190000054
(19) in the formula, the parameter z2,ζ,ωndAll can be obtained by the parameter calculation of the SFR model;
the lagrange transform of the output signal of the SFR model under the unit step function is obtained from equation (19):
Figure FDA0002209079190000055
in formula (20):
Figure FDA0002209079190000056
performing reverse Laplace transformation on the formula (20) to obtain
Figure FDA0002209079190000057
(21) In the formula:
Figure FDA0002209079190000061
let t bezAt this time, the system frequency offset is maximum, and the slope of the frequency change amount at this time is zero, which can be expressed by equation (22):
Figure FDA0002209079190000062
the maximum offset occurrence time t of the system frequency can be obtained by the formula (22)zBy adjusting WASHOUT unit KDParameter, maximum frequency offset of system Δ fdevmaxIs maintained at Δ fmaxAccordingly, formula (23) can be obtained:
the WASHOUT unit coefficient K can be obtained from the expressions (22) - (23)D
CN201710681615.9A 2017-08-10 2017-08-10 Method for controlling frequency of isolated power grid by participation of electrolysis load Active CN107317345B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710681615.9A CN107317345B (en) 2017-08-10 2017-08-10 Method for controlling frequency of isolated power grid by participation of electrolysis load

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710681615.9A CN107317345B (en) 2017-08-10 2017-08-10 Method for controlling frequency of isolated power grid by participation of electrolysis load

Publications (2)

Publication Number Publication Date
CN107317345A CN107317345A (en) 2017-11-03
CN107317345B true CN107317345B (en) 2020-01-24

Family

ID=60175652

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710681615.9A Active CN107317345B (en) 2017-08-10 2017-08-10 Method for controlling frequency of isolated power grid by participation of electrolysis load

Country Status (1)

Country Link
CN (1) CN107317345B (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108565874B (en) * 2018-03-06 2020-02-14 合肥工业大学 Source-load cooperative frequency modulation method based on load frequency control model
CN108808688B (en) * 2018-07-04 2021-07-02 云南电网有限责任公司 Stable control load shedding method based on electrolytic aluminum control characteristics
CN110120676B (en) * 2019-06-18 2023-03-24 国网新疆电力有限公司电力科学研究院 MMC-HVDC converter station power control method and system based on simulation synchronous motor characteristics
CN110492523A (en) * 2019-07-11 2019-11-22 北京科东电力控制系统有限责任公司 The low-frequency load shedding optimization method of the high accounting system of new energy
CN110429619B (en) * 2019-08-31 2022-08-09 上海交通大学 Disaster recovery reaction method for demand side frequency emergency control of super-large-scale power grid
CN111258232B (en) * 2020-02-05 2023-05-12 国能信控互联技术有限公司 Fire-storage combined frequency modulation semi-physical simulation system
CN111864768A (en) * 2020-07-22 2020-10-30 山东大学 Control method and system for participation of electrolytic aluminum load in primary frequency modulation
CN112564128A (en) * 2020-11-17 2021-03-26 云南电网有限责任公司 Control system and method for electrolytic aluminum load participating in power grid frequency modulation
CN112417792B (en) * 2020-11-23 2023-05-26 云南电网有限责任公司 Electrolytic aluminum external characteristic modeling method based on controllable boundary of saturation reactor
CN112993987A (en) * 2021-03-05 2021-06-18 云南电网有限责任公司 Active adjustable capacity cooperative electrolytic aluminum load coordination power grid control method
CN113162065A (en) * 2021-03-11 2021-07-23 国网山东省电力公司 Power grid frequency modulation method and system based on cooperative coordination between production loads
CN113141015B (en) * 2021-03-18 2023-05-26 云南电网有限责任公司 Control method for electrolytic aluminum load to participate in frequency modulation of power grid at transmitting end
CN113258589A (en) * 2021-04-25 2021-08-13 国网山东省电力公司 Energy storage and electrolytic aluminum load polymerization frequency modulation method and device based on production operation
CN113765120A (en) * 2021-06-11 2021-12-07 中国电力科学研究院有限公司 Power grid frequency modulation method based on cooperative coordination of industrial park submerged arc furnace and polycrystalline silicon
CN114243731B (en) * 2022-02-23 2022-06-14 武汉大学 Controlled load damping factor controller and frequency modulation method adopting controller

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102593833A (en) * 2012-03-19 2012-07-18 上海电力学院 Method for designing frequency controller of independent micro-grid
CN103457297A (en) * 2013-07-12 2013-12-18 武汉大学 Closed-loop control method for achieving active power adjustment by introducing frequency deviation feedback to have control over voltage
CN105406496A (en) * 2015-12-15 2016-03-16 北京四方继保自动化股份有限公司 Isolated microgrid frequency modulation control method based on measured frequency response identification

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102593833A (en) * 2012-03-19 2012-07-18 上海电力学院 Method for designing frequency controller of independent micro-grid
CN103457297A (en) * 2013-07-12 2013-12-18 武汉大学 Closed-loop control method for achieving active power adjustment by introducing frequency deviation feedback to have control over voltage
CN105406496A (en) * 2015-12-15 2016-03-16 北京四方继保自动化股份有限公司 Isolated microgrid frequency modulation control method based on measured frequency response identification

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
An Industrial Verification of Frequency Regulation by Electrolytic Aluminum in an Isolated Power System;Yi Bao,etc;《2016 IEEE Power and Energy Society General Meeting》;20160717;第1-5页 *

Also Published As

Publication number Publication date
CN107317345A (en) 2017-11-03

Similar Documents

Publication Publication Date Title
CN107317345B (en) Method for controlling frequency of isolated power grid by participation of electrolysis load
US8058753B2 (en) Wide area transmission control of windfarms
CN102761133B (en) Method for controlling frequency modulation of micro-grid battery energy storage system based on fuzzy control
CN106786796B (en) Control method and system for wind power to participate in frequency modulation of power system
CN102280879B (en) Method and system for regulating power of large-scale energy storage power station of wind farm
CN105743107B (en) A kind of control method of electric power Isolated Network System frequency regulation
CN109103927B (en) PID controller parameter setting method for speed regulating system
CN103972899B (en) A kind of STATCOM access point voltage compensating method
CN108462212B (en) Control method of new energy power system in multi-source multi-regulation-control-domain operation mode
Bao et al. Field verification of frequency control by energy-intensive loads for isolated power systems with high penetration of wind power
CN101860044A (en) Method for cooperatively controlling reactive voltage of wind farm
CN111258232B (en) Fire-storage combined frequency modulation semi-physical simulation system
CN107658909A (en) A kind of low and medium voltage distribution network voltage lifting suppressing method of the access containing photovoltaic
CN103580569B (en) A kind of flexible extendible exciter control system
CN105048518A (en) Control method for photovoltaic participated isolated grid frequency modulation
CN110808602A (en) Improved additional frequency control method and system for multi-terminal flexible direct current power transmission system
CN103151795A (en) Scattered-type wind power plant reactive power optimization control method capable of reducing fan losses and system
CN106487024B (en) Wind power plant reactive compensation device and reactive replacement method and device of wind power generation set
CN105135409A (en) Supercritical unit boiler master controller control method based on primary frequency modulation action amplitude
CN104485670B (en) The control method of voltage sensitivity industrial load time-varying damping characteristic in island network
Dragosavac et al. Coordinated reactive power-voltage controller for multimachine power plant
CN108879721A (en) A kind of control method for frequency based on wind electricity digestion
CN105610200A (en) Synchronous coordinated control based full-power control method for thermal power plant
CN105717793A (en) Optimizing control method and system capable of ensuring motion direction of primary frequency modulation
CN109888797A (en) The control method for frequency of sending containing wind power integration based on Situation Awareness

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