CN110880772A - Electricity selling company response power grid control method based on industrial park load aggregation - Google Patents

Electricity selling company response power grid control method based on industrial park load aggregation Download PDF

Info

Publication number
CN110880772A
CN110880772A CN201911088544.7A CN201911088544A CN110880772A CN 110880772 A CN110880772 A CN 110880772A CN 201911088544 A CN201911088544 A CN 201911088544A CN 110880772 A CN110880772 A CN 110880772A
Authority
CN
China
Prior art keywords
load
power
voltage
industrial park
regulation
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.)
Granted
Application number
CN201911088544.7A
Other languages
Chinese (zh)
Other versions
CN110880772B (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 CN201911088544.7A priority Critical patent/CN110880772B/en
Publication of CN110880772A publication Critical patent/CN110880772A/en
Application granted granted Critical
Publication of CN110880772B publication Critical patent/CN110880772B/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/24—Arrangements for preventing or reducing oscillations of power in networks
    • 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
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00—Climate change mitigation technologies for sector-wide applications
    • Y02P80/20—Climate change mitigation technologies for sector-wide applications using renewable energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Electrical Variables (AREA)

Abstract

The invention belongs to the power system operation and control technology, and particularly relates to an electric power selling company response power grid control method based on industrial park load aggregation, which models the power regulation characteristic of a typical industrial park load, provides a proportional aggregation mode in which a user reports a regulation and control range in priority by combining a load regulation willingness and a limit regulation range, and establishes an aggregation control strategy of the electric power selling company on an industrial park; considering load production benefits, establishing control cost models of various types of loads, and simultaneously providing a coordination control strategy for the power selling company to participate in demand response in multiple parks; aiming at the scene of wind power fluctuation of a large power grid, the stability of the power grid is maintained through load response power fluctuation, and the correctness of the provided control strategy is verified. For a power grid company, the evaluation indexes of low wind abandoning and light abandoning rates are facilitated to be realized; for a new energy power generation party, the new energy access rate is improved, and the total economic benefit of load-network-source is improved on the whole.

Description

Electricity selling company response power grid control method based on industrial park load aggregation
Technical Field
The invention belongs to the technical field of operation and control of power systems, and particularly relates to a power selling company response power grid control method based on load aggregation of an industrial park.
Background
With the gradual advance of energy strategies to large-scale new energy power generation and networking, the influence of uncertainty of new energy on safe and stable operation of a power grid becomes an urgent problem to be solved. In consideration of economy and feasibility, the novel flexible resource of load demand side response is applied, and the important significance is achieved in participating in peak regulation and frequency modulation of a power grid and absorbing new energy. The industrial load, especially the high energy consumption loads such as electrolysis load and arc load, is the main type, has the characteristics of high power consumption and stable power, and has great power regulation and control potential. For the networking form of the industrial park and the large power grid system, how to develop load demand response and formulating a coordination control method of various loads are problems to be solved urgently.
Meanwhile, the business of the electricity selling side is gradually expanded, and the electricity selling company also faces continuous challenges under market opportunity. The business mode of simply earning the selling electricity price difference is not enough, and the expansion of the business modes such as participation in power grid auxiliary service, construction of comprehensive energy service and the like is a core problem faced by electricity selling companies.
In conclusion, research on an aggregation strategy of the power selling company on the load of the industrial park to participate in demand response and auxiliary power grid stabilization of power fluctuation caused by new energy access are good modes of multi-party participation and multi-party mutual profit. The power selling company participates in the power grid auxiliary service through an aggregation control strategy, and business content is expanded; the industrial load can obtain certain economic compensation under the condition of permission of production conditions; the power grid company can guarantee the stable and reliable operation of the power grid at lower cost; the new energy power generation party improves the new energy access rate; the improvement of the total economic benefit of the load-net-source is realized on the whole.
Disclosure of Invention
The invention aims to provide a method for obtaining load aggregation characteristic delta P of an industrial park by aggregating loads of the industrial park by an electricity selling company∑,t-ntAnd control cost model F of industrial parkIP-ΔPΣReporting the adjustable capacity and the regulation and control cost of the power selling company to the power grid, issuing a power regulation instruction to each industrial park and different types of loads thereof by a bidding preferred coordination method according to the power grid regulation instruction, and finally realizing the method for stabilizing the power fluctuation.
In order to achieve the purpose, the invention adopts the technical scheme that the power selling company response power grid control method based on the load aggregation of the industrial park comprises the following steps:
step 1, carrying out load power normalization processing on an industrial park according to the power characteristics of electrolytic aluminum, a submerged arc furnace and a polycrystalline silicon load, and analyzing the power theoretical regulation boundary of each load;
step 2, a proportional type aggregation model which gives priority to the reporting regulation and control range of the user is provided in combination with the load regulation willingness and the limit regulation range, and an aggregation control strategy of the power selling company to the industrial park is established;
step 3, establishing a response strategy of the power selling company to the power grid control instruction according to the aggregation model, and realizing reasonable distribution and aggregation control of the load power of the industrial park;
step 4, analyzing the influence of the aggregation control strategy on the load production benefit, establishing a load control cost model, and establishing an industrial park load control cost model by combining the aggregation model;
and 5, based on the industrial park load control cost model, providing a bidding and preferential coordination strategy of the power selling company for multiple industrial parks, and realizing the coordinated distribution and control of the power of the multiple parks.
In the above method for controlling a power grid response of an electric power selling company based on load aggregation in an industrial park, the specific implementation of step 1 includes:
step 1.1, the load power characteristic of the electrolytic aluminum is as follows:
Figure BDA0002266180100000021
wherein, PALFor power of electrolytic aluminum, VBIs the DC bus voltage of the electrolyzer IdIs direct current of the electrolytic cell, RECIs the equivalent resistance of the electrolytic cell in series connection, and E is the equivalent potential of the electrolytic cell;
electrolytic aluminum DC bus voltage VBVoltage V of high-voltage bus connected with electrolytic aluminium loadAL-AHThe relationship of (1) is:
Figure BDA0002266180100000022
in the formula, LSRIs the equivalent value of the saturable reactor, VAL-AHIs a high voltage bus voltage, and omega is a voltage angular frequency;
considering that the electrolytic aluminum load is connected with an on-load tap changer, the transformation ratio value of the on-load tap changer is m1Stage, then change its transformation ratio kALCan realize m1Stage adjustment:
Figure BDA0002266180100000031
saturable reactor equivalent value L when production requirements are metSRHas a regulation range of [ LSR min,LSR max]Then, from the formula (2), the strain ratio k can be obtainedAL-iVoltage drop V of direct current bus of lower electrolytic aluminumBThe adjusting range is as follows:
Figure BDA0002266180100000032
Figure BDA0002266180100000033
the maximum adjustable power range of the electrolytic aluminum load obtained by the formula (1) is as follows:
Figure BDA0002266180100000034
step 1.2, the power characteristics of the submerged arc furnace load are as follows:
Figure BDA0002266180100000035
Figure BDA0002266180100000036
wherein, PSAF、QSAFRespectively active power and reactive power of the submerged arc furnace, USAFIs the voltage of the low-voltage side of the submerged arc furnace, RlineAnd XlineThe equivalent resistance and the equivalent reactance of the short net;
the equivalent impedance of the arc is defined by an arc static resistance R in a certain time break planearcAnd arc static reactance XarcTo characterize;
low-voltage side voltage U of submerged arc furnaceSAFThe voltage U of a high-voltage bus of the submerged arc furnaceSAF-AHObtained by a special transformer for a submerged arc furnace:
USAF＝USAF-AH/kSAF(9)
in the formula, kSAFThe transformer is a transformer special for the submerged arc furnace;
the arc impedance relational expression of the submerged arc furnace is obtained by fitting actual production data:
Xarc＝a1Rarc 2+a2Rarc+a3(10)
wherein, a1、a2、a3Is the arc impedance fitting coefficient, which is a constant;
when the power of the submerged arc furnace is adjusted by adopting a constant impedance voltage-adjusting mode, the transformation ratio of the transformer special for the submerged arc furnace load is set to be m in total2A stage; the voltage of the transformer is m in total by changing the transformation ratio of the transformer special for the submerged arc furnace2Stage regulation range:
Figure BDA0002266180100000041
arc static resistance RarcThe limiting range of (2):
Rarc min≤Rarc≤Rarc max(12)
upper limit of arc partial power factor
Figure BDA0002266180100000042
There is a lower limit to the arc segment power factor
Figure BDA0002266180100000043
Arc partial power factor
Figure BDA0002266180100000044
The limitations of (2) are:
Figure BDA0002266180100000045
namely, it is
Figure BDA0002266180100000046
Wherein the content of the first and second substances,
Figure BDA0002266180100000047
meanwhile, each smelting stage should meet the minimum power constraint:
PSAF,t≥PSAF,t min(15)
in the formula, PSAF,t minThe minimum power of the submerged arc furnace;
the rated voltage of the low-voltage side of the submerged arc furnace is set as USAF,NBy the formulas (7), (8), (12) and (15), when the constant voltage impedance power regulation method is adopted, the active power P of the submerged arc furnace in each smelting stageSAF,tAnd reactive power QSAF,tThe adjustment range of (a) is as follows:
c1≤PSAF,t≤c2(16)
d1≤QSAF,t≤d2(17)
wherein:
Figure BDA0002266180100000051
Figure BDA0002266180100000052
Figure BDA0002266180100000053
Figure BDA0002266180100000054
the power characteristic parameter of the ore furnace during rated operation is
Figure BDA0002266180100000055
The power regulation range at this time is converted into an equivalent low-voltage side voltage variation range, namely
Figure BDA0002266180100000056
From the equations (10) and (18), when the impedance-voltage coordinated regulation method is adopted, the voltage variation range of the low-voltage side of the submerged arc furnace load is as follows:
Figure BDA0002266180100000057
step 1.3, the polysilicon load power meets the following electrical relationship:
Figure BDA0002266180100000058
in the formula, PPCSFor polysilicon load AC total power, UvalLoading single phase electricity for polysiliconPressure, RPCSIs a polysilicon rod single-phase resistor;
the production process energy conversion of a polycrystalline silicon rod in the time delta t is as follows:
Figure BDA0002266180100000061
in the formula, PPCSFor polysilicon load AC total power, Δ Qout1V represents the amount of heat used to heat the reactant gas, and is given by the gas specific heat capacity formula1·Δt·s1·ρg·c·(Tx-Tg) Wherein v is1、s1、ρg、c、TgRespectively, the air inlet rate, the air inlet area, the mixed gas density, the mixed gas specific heat capacity and the air inlet temperature are constants; tx, surface temperature of silicon rod,. DELTA.Qout2And Δ Qout3Respectively, the heat quantity of the heat of the endothermic reaction maintained and dissipated through the base and the wall of the reduction furnace by thermal radiation
Figure BDA0002266180100000062
Corresponds to (Δ Q)out2+ΔQout3) Wherein, η, K, L, ToutThe reaction heat absorption ratio, the total heat transfer coefficient of the silicon rod and the mixed gas, the total length of the silicon rod, and the equivalent temperature of the chassis and the furnace wall surface are all constants; r is the radius of the polysilicon rod, and r can be regarded as a constant in a short time;
the power characteristic equation of the polysilicon load with radius r obtained from equations (20) and (21) is as follows:
Figure BDA0002266180100000063
in the formula, A, B, C, D, G, H is a fitting coefficient of the power characteristic of the polysilicon, is a constant and is obtained by fitting actual production rated operation data; i is polysilicon load single-phase current;
for the surface temperature T of the silicon rodxThe control range is as follows:
Tx min≤Tx≤Tx max(23)
when the temperature is less than or equal to 1000 ℃ TxCan ensure production at the temperature of less than or equal to 1100 ℃, and can ensure the production at Tx＝Tx，NThe temperature is the optimum temperature at 1080 ℃; when engaged, the polysilicon load generally engages in downward regulation of power, then there is Tx min＝1000℃,Tx max＝1080℃；
The flow rate of the cooling water is generally regulated, and the flow rate regulation rate of the cooling water is α
Figure BDA0002266180100000071
αmin≤α≤αmax(25)
Wherein, αmin＝90％，αmaxWhen the operation is rated, α is 100%;
the following equations (20) and (21) can be obtained:
Figure BDA0002266180100000072
from equations (23) and (26), the adjustment range of the polysilicon load power can be determined as follows:
Figure BDA0002266180100000073
Figure BDA0002266180100000074
effective value U of single-phase voltage of polysilicon loadvalThe adjusting range of (A) is as follows:
Figure BDA0002266180100000075
step 1.4, carrying out load normalization processing on the industrial park;
setting the number of electrolytic aluminum load, submerged arc furnace load and polysilicon load in the industrial park as N respectivelyAL、NSAF、NPCSGet itRated voltage of each load is a voltage reference value, and voltage of each load is U*The minimum and maximum limit values of each load voltage are respectively
Figure BDA0002266180100000078
Step 1.4.1. for electrolytic aluminum load i ═ {1,2, …, NALHas the following components:
PAL,i＝m1,i·U*2+m2,iU*(30)
Figure BDA0002266180100000076
in the formula, PAL,iThe power of the ith electrolytic aluminum load,
Figure BDA0002266180100000077
direct current bus voltage rating, R, for the ith cell for aluminum electrolysisEC,iThe electrolytic cell for the ith electrolytic aluminum is connected with an equivalent resistance in series, EiThe cell equivalent potential for the ith electrolytic aluminum;
DC bus voltage per unit value U of electrolytic aluminum load*The limit variation range of (2) is:
Figure BDA0002266180100000081
step 1.4.2. for ore furnace load i ═ NAL+1,NAL+2,…,NAL+NSAF}：
PSAF,i＝m1,i·U*2(33)
QSAF,i＝ni·U*2(34)
Figure BDA0002266180100000082
Figure BDA0002266180100000083
In the formula, PSAF,i、QSAF,iRespectively the load active power and the load reactive power of the ith submerged arc furnace,
Figure BDA0002266180100000084
is the low side voltage rating of the ith furnace,
Figure BDA0002266180100000085
for the power characteristic parameter, R, of the ith submerged arc furnace during rated operationline,iAnd Xline,iIs the equivalent resistance and the equivalent reactance, R, of the ith submerged arc furnace short netarc,iAnd Xarc,iThe resistance and the reactance of the arc static resistance of the ith submerged arc furnace are shown;
per unit value U of voltage at low voltage side of submerged arc furnace*The limit variation range of (2) is:
Figure BDA0002266180100000086
step 1.4.3. load i ═ N to polysiliconAL+NSAF+1,NAL+NSAF+2,…,NAL+NSAF+NPCS}：
PPCS,i＝m1,i·U*2(38)
Figure BDA0002266180100000087
In the formula, PPCS,iThe power is loaded for the ith poly-silicon,
Figure BDA0002266180100000088
for the ith polysilicon load single phase voltage rating, RPCS,iThe resistance is the single-phase resistance of the ith polycrystalline silicon rod;
the limit variation range of the per unit value of the polysilicon load single-phase voltage is as follows:
Figure BDA0002266180100000091
and 1.4.4. Total load power characteristic P of industrial park∑-U*Comprises the following steps:
P∑＝M·U*2+M'·U*(41)
wherein the content of the first and second substances,
Figure BDA0002266180100000092
m and M' are the quadratic coefficient and the first-order coefficient of the total load power characteristic, respectively, M1,iAnd m2,iThe quadratic coefficient and the first order coefficient for each type of load power characteristic are obtained by equations (31), (35), and (39).
In the above method for controlling a power grid response of an electric power selling company based on load aggregation in an industrial park, the implementation of step 2 includes the following steps:
step 2.1, giving priority to a proportional aggregation model of a user reporting regulation and control range:
the voltage variation range which is expected to participate in regulation is automatically reported by setting the load as
Figure BDA0002266180100000095
Voltage regulation dead zone of
Figure BDA0002266180100000096
Considering the voltage limit variation range of each type of load determined by the equations (6), (19) and (40), the voltage variation range in which the load actually participates in regulation is:
Figure BDA0002266180100000093
the voltage of each load can be adjusted to a practically adjustable minimum/maximum value of
Figure BDA0002266180100000094
The actual maximum allowable up/down adjustment amount of each load voltage is respectively as follows:
Figure BDA0002266180100000101
when power regulating quantity distribution is carried out, voltage is regulated among all loads in equal proportion according to the voltage regulation range;
when the system is initially set and is not involved in regulation, all loads are in a rated operation state, namely U * i,01 is ═ 1; the variation between the load voltage at time t and the initial voltage is recorded as Δ U* i,tI.e. by
ΔU* i,t＝U* i,t-U* i,0＝U* i,t-1 (44)
The load voltage adjustment amount at the adjacent time is:
U* i,t-U* i,t-1＝ΔU* i,t-ΔU* i,t-1(45)
namely, it is
Figure BDA0002266180100000102
In the formula, ntThe load adjusting parameters are proportional adjusting parameters, the adjusting parameters of each load are consistent during each adjustment, and the actual adjusting quantity is in positive correlation with the adjustable range of each load;
setting the power regulating quantity of the ith load in the park at the time t as delta P relative to the initial powerload,i,t＝Pload,i,t-Pload,i,0The total power regulation quantity of the load in the park is delta P relative to the initial power∑,t＝PΣ,t-P∑,0；
The aggregate characteristics of the respective loads can be obtained from the load regulation characteristics (30), (33), (38):
Figure BDA0002266180100000103
and also
Figure BDA0002266180100000111
The load polymerization characteristics of the industrial park were:
Figure BDA0002266180100000112
wherein the content of the first and second substances,
Figure BDA0002266180100000113
n is to betSubstituting 1 into the expressions (47) and (48) can respectively obtain the maximum upward adjustment capacity delta P of each load and industrial parkload,i,up max、ΔPΣ,up max(ii) a N is atThe maximum downward regulating capacity delta P of each load and industrial park can be respectively obtained by substituting the value of-1load,i,down max、ΔPΣ,down max；
Step 2.2, the power selling company responds to the aggregation algorithm of the power grid requirement;
step 2.2.1, an off-line polymerization stage;
step 2.2.1.1, solving the load power characteristic of the industrial park: respectively determining power characteristics of electrolytic aluminum, a submerged arc furnace and a polycrystalline silicon load according to formulas (1), (7), (8), (20) and (22) based on the measured data;
step 2.2.1.2. solving the total load power characteristic: obtaining the load power characteristic parameters of the expressions (31), (35), (36) and (39) and obtaining the total load power characteristic of the expression (41);
step 2.2.1.3, solving the voltage limit regulation range: the voltage limit variation range of each load is obtained by combining equations (4), (5), (19) and (29) with equations (32), (37) and (40)
Figure BDA0002266180100000114
Step 2.2.1.4, determining the actual voltage regulation range: recording the voltage regulation range and the regulation dead zone reported by each load independently, and solving the actual maximum upward regulation quantity of the voltage of each load by the formulas (42) and (43);
step 2.2.1.5. calculating the load polymerization characteristics and the maximum regulating capacity of the industrial park: the load of the industrial park is obtained from the formulae (48) and (49)Polymerization characteristics Δ P∑,t-nt(ii) a Let n betDetermining the maximum upward and downward regulating capacity DeltaP as 1 or-1Σ,up max、ΔPΣ,down max；
Step 2.2.2, an on-line calculation stage;
setting the power of the industrial park to be delta P at the previous moment of load controlΣ,t-1Initial adjustment time Δ PΣ,0When the maximum upward and downward regulation capacity of the industrial park is equal to 0, the maximum upward and downward regulation capacity of the industrial park at the current time t is delta PΣ,up,t max＝ΔPΣ,up max-ΔPΣ,t-1；
Step 2.2.3, reporting the adjustable capacity of the power grid at the current moment;
the electricity selling company adjusts the maximum upward and downward capacity delta P at the momentΣ,up,t max、ΔPΣ,down,t maxAnd uploading to a power grid dispatching center.
In the above method for controlling a power grid response of an electric power selling company based on load aggregation in an industrial park, the implementation of step 3 includes the following steps:
step 3.1, a load distribution principle;
when the power grid issues a control target power P to the power selling companyt netThen, the power selling company distributes power to each load in the industrial park;
setting the total power regulating quantity of the industrial park at t moment as delta P∑,tThe solution is as follows:
ΔP∑,t＝ΔPt net＝Pt net-Pt-1 net(50)
the power regulation quantity Delta P is obtained from the equation (48)∑,tCorresponding aggregate instruction nt：
Figure BDA0002266180100000121
From equations (44) and (46), the voltage control target of each load is:
Figure BDA0002266180100000122
the actual power of each load after participating in the regulation can be obtained according to the formulas (1), (7), (8) and the formulas (20) and (22) respectively;
step 3.2, a load control instruction is carried out;
step 3.2.1, controlling the load of the electrolytic aluminum;
the direct current bus voltage V of the electrolytic cellB,i,t＝VBN,i·U* i,tSubstituting into formula (2) can solve the saturated reactor equivalent value L of the electrolytic aluminum loadSR,iComprises the following steps:
Figure BDA0002266180100000131
the direct current bus voltage V of the electrolytic cellB,i,t＝VBN,i·U* i,tSubstituting the formula (1) into the load to obtain the corresponding rectified current value I of the loadd,i,tComprises the following steps:
Figure BDA0002266180100000132
step 3.2.2, a control instruction of the submerged arc furnace load;
using electrode current command value IrefTo control the lifting of the electrode; the electrode hydraulic lifting model comprises the following steps:
Figure BDA0002266180100000133
wherein L is the distance between the electrode and the furnace charge, is equal to the arc length of the electric arc, and L0As the initial position of the electrode, IrefIs an electrode current command value; v. ofupAnd vdownMaximum values representing the electrode rising and falling speeds, both constant, cupAnd cdownIs a proportionality coefficient between the current difference and the speed;
when changing the arc impedance, correspondingly adjusting I in the lifting model according to the following formularefAnd the electrode lifting can be realized:
Figure BDA0002266180100000134
when the voltage of the low-voltage side of the submerged arc furnace is
Figure BDA0002266180100000135
The electrode control command is then obtained by:
Figure BDA0002266180100000141
step 3.2.3, controlling the polysilicon load;
from the equation (20), it can be seen that the polysilicon load single-phase voltage U is changedvalThe adjustment of the load power of the polysilicon can be realized; controlling U by wave splicing technologyvalThe following are:
Figure BDA0002266180100000142
in the formula, the voltage angular frequency omega is constant, U1、U2To the splicing voltage, t1The wave splicing time;
when obtaining the polysilicon load single-phase voltage Uval,i,t＝UvalN,i·U* i,tThen, firstly, the wave splicing voltage U is determined1、U2: the value of the wave-splicing voltage is U1,U2E is (0,380,600,800,1500) V, and the fixed splicing wave voltage U is obtained1、U2Are respectively a target voltage Uval,i,tTwo adjacent voltage levels U1＜U2(ii) a Then put Uval,i,t、U1、U2Substitution formula (58) for calculating splicing wave voltage time t1(ii) a Will Uval,i,tSubstituting formula (26) and solving the surface temperature T of the silicon rodxCooling water flow rate α, wherein the surface temperature T of the silicon rod is maintained preferentially during regulationx＝1080℃。
In the above method for controlling a power grid response of an electric power selling company based on load aggregation in an industrial park, the implementation of step 4 includes the following steps:
step 4.1, controlling the influence of the strategy on the load production benefit;
according to different directions of the load participation of the industrial park in the upward/downward power regulation, the load control cost is divided into:
step 4.1.1, when load power regulating quantity delta PloadWhen less than 0, the loss of unit value FvCalculated from the following formula:
Fv＝(Fp-Fc)/CE(59)
wherein FvIs the loss of value per unit load, unit cell/kilowatt-hour, FpIs sold per unit load, unit/ton, FcProduction cost per unit load, unit/ton, CEThe unit of power consumption is the unit of output and kilowatt hour/ton;
step 4.1.2 when load power regulating quantity delta PloadAnd when the load is more than 0, modeling the loss of the overload operation based on the life cycle model:
the maintenance cost of the load is set as lambdaiUnit cell, rated maintenance life cycle τi,NThe unit h; when the load operation is set at time t, the service life changes to tau due to the overload operationi,tUnit h, the amount of change ρ in the hourly converted maintenance cost of the loadi,tUnit cell/h, calculated by:
Figure BDA0002266180100000151
recording the load power as PN,iIn kW unit, the reduced maintenance cost of the unit power consumption of the load is the maintenance unit price FreUnit cell/kwh, is solved by:
Figure BDA0002266180100000152
in summary, the control cost F of the unit electric quantity of the load0Unit cell/kwh, is:
Figure BDA0002266180100000153
step 4.2, controlling a cost model in the industrial park;
setting the control cost of each high energy consumption load in an industrial park as F0,iWherein i is 1,2, … NAL+NSAF+NPCSThe electrolytic aluminum, the submerged arc furnace and the polycrystalline silicon are sequentially numbered according to the load sequence;
the control cost of each load is related to the adjustment power by the following formula:
Figure BDA0002266180100000154
wherein Fload,i,tThe control cost for the ith load;
power regulation of industrial park
Figure BDA0002266180100000155
Let the total control cost of the industrial park be FIP,t：
Figure BDA0002266180100000156
Solving the control cost F of the industrial parkIP,tAnd adjusting the electric quantity delta PΣ,tFunctional relationship F ofIP,t-ΔPΣ,tIn combination with the load distribution and characteristic Δ Pload,i,t-ntSolving the control cost of the industrial park-polymerization model FIP,tNt, and then the model F of the power regulation cost of the industrial park is solvedIP,t-ΔPΣ,tNamely:
Figure BDA0002266180100000161
the power adjustment for the ith load on the industrial park is determined by:
Figure BDA0002266180100000162
from the equations (65) and (66), the polymerization model F, which is the control cost of the industrial park, can be determinedIP,tNt is as follows:
Figure BDA0002266180100000163
wherein the content of the first and second substances,
Figure BDA0002266180100000164
when n istWhen 1, the cost F of providing the maximum upward adjustable power for the industrial park can be obtainedIP,t max(ii) a When n istWhen-1, the cost F for providing the maximum downwardly adjustable power for the industrial park is obtainedIP,t min；
The power regulation cost model F of the industrial park can be obtained by the formulas (51) and (67)IP,t-ΔPΣ,tThe following formula:
Figure BDA0002266180100000165
wherein the content of the first and second substances,
Figure BDA0002266180100000171
in the above method for controlling a power grid response of an electric power selling company based on load aggregation of an industrial park, the implementation of step 5 includes the following steps:
step 5.1, a control cost model and a coordination control method of the power selling company;
the number of the management industrial parks of the electricity selling company is NIPThe industrial parks are sorted according to the minimum price of the times of power up regulation, and the serial numbers are recorded as u1 and u2 … uNIP(ii) a The times of price minimum sequencing of each industrial park are regulated downwards according to power, and the serial numbers are recorded as d1 and d2 … dNIPThus, the control cost model F of the electricity selling companyPSC-ΔPPSCComprises the following steps:
Figure BDA0002266180100000172
in the formula, FPSCIs the control cost, Δ P, of the electricity selling companyPSCIs the sum of the total regulated power of the power selling company, i.e. the regulated power of each industrial park, delta PΣ,up ui,max、ΔPΣ,down di,maxThe maximum upward regulating capacity of the industrial park numbered ui and the maximum downward regulating capacity of the industrial park numbered di are respectively;
according to a bidding preference strategy, the power of each industrial park is regulated in a sequence from low to high in control cost according to regulation requirements; the relationship between the total regulated power of the power selling company and the power of each park is as follows:
Figure BDA0002266180100000181
the total power of each industrial park can be determined by equation (70), and the aggregate demand n for each parktThe voltage target value of each load and the actual control command of each load can be obtained by the equation (51) and the equations (52), (54), (57), (58) and (26);
step 5.2, the coordination control algorithm of the power selling company to the multiple parks;
step 5.2.1, an off-line polymerization stage;
step 5.2.1.1. each park load reports the voltage regulation range independently
Figure BDA0002266180100000182
Calculating control cost by the formula (62);
step 5.2.1.2. evaluation of the formula (43) confirms the respective load control ranges, and calculation of the polymerization characteristics Δ P of the respective parks from the formula (48)Σ,t-ntCalculating a power regulation cost model F for each park from equation (67)IP-ΔPΣ；
Step 5.2.1.3, calculating the control cost model F of the power selling company by the formula (69)PSC-ΔPPSC；
Step 5.2.2, an on-line polymerization calculation stage;
solving the maximum upward/downward adjustable capacity delta P of each park at the current time t from the last control time t-1Σ,up,t ui,max＝ΔPΣ,up ui,max-ΔPΣ,t-1、ΔPΣ,down,t di,max＝ΔPΣ,down di,max-ΔPΣ,t-1And the control cost of the power selling company at the time t is FPSC,t＝FPSC(ΔPPSC-ΔPΣ,i)；
Step 5.2.2, reporting to the power grid; updating the adjustable capacity and cost model and reporting to the power grid;
step 5.2.3, receiving a power grid power regulation requirement;
step 5.2.3.1 receiving a grid power regulation demand Δ Pnet,tThen, the total power regulating quantity delta P of each industrial park is confirmed in sequence according to the formula (70) and the current adjustable capacity of each parkΣ,t；
Step 5.2.3.2, distributing the load power in each industrial park according to the step 3, and determining each control instruction;
and 5.2.3.3, recording and updating the current load state and waiting for the next adjustment.
The invention has the beneficial effects that: modeling the power regulation characteristic of the load of a typical industrial park based on the load production characteristic and the regulation and control means; synthesizing load regulation willingness and limit regulation range, providing a proportional aggregation model which gives priority to the regulation and control range reported by a user, and establishing an aggregation control strategy of an electricity selling company for an industrial park according to the proportional aggregation model to respond to the power grid frequency modulation requirement; considering load production benefits, establishing control cost models of various types of loads, and providing a coordination control strategy for the electricity selling company to respond to the demands of participation of multiple parks by combining the cost models and the aggregation models; aiming at the scene of wind power fluctuation of a large power grid, the stability of the power grid is maintained through load response power fluctuation, and the correctness of the provided control strategy is verified.
Research electric selling companies participate in demand response by an aggregation strategy of industrial park loads, and an auxiliary power grid stabilizes power fluctuation caused by new energy access, so that the electric selling companies are a good mode of multi-party participation and multi-party mutual profit. The power selling company participates in the power grid auxiliary service through an aggregation control strategy, and business content is expanded; the industrial load can obtain certain economic compensation under the condition of permission of production conditions, and the operational dilemma of the current high energy consumption load is solved; for a power grid company, the stability and reliability of a power grid can be guaranteed at a lower cost, and the evaluation indexes of lower 'wind abandoning' and 'light abandoning' rates are facilitated to be realized; for a new energy power generation party, the new energy access rate is improved, and more new energy power generation supplementary benefits are obtained; the improvement of the total economic benefit of the load-net-source is realized on the whole.
Drawings
FIG. 1 is an overall block diagram of an electricity vendor's response to a power grid aggregate control strategy according to one embodiment of the present invention;
FIG. 2 is a representative industrial park total load power characteristic of one embodiment of the present invention;
FIG. 3 is a model of the aggregate characteristics of a typical industrial park according to one embodiment of the present invention;
FIG. 4 is a wind power fluctuation plot for one minute for one area according to an embodiment of the present invention;
FIG. 5(a) is a graph of aggregate command changes in response to the windage fluctuations of FIG. 3 in accordance with an embodiment of the present invention;
FIG. 5(b) is a graph of the change in electrolytic aluminum load voltage in response to wind power fluctuations of FIG. 3 in accordance with an embodiment of the present invention;
FIG. 5(c) is a graph of the change in the load voltage of the submerged arc furnace in response to the wind power fluctuation of FIG. 3 according to an embodiment of the present invention;
FIG. 5(d) is a graph of the change in voltage for 75mm polysilicon in response to the wind electrical fluctuations of FIG. 3 in accordance with one embodiment of the present invention;
FIG. 6(a) is a graph illustrating a tie line fluctuation without using the aggregation control strategy for the wind power fluctuation of FIG. 3 according to an embodiment of the present invention;
FIG. 6(b) is a tie line fluctuation diagram of the wind power fluctuation of FIG. 3 using the aggregation control strategy according to an embodiment of the present invention;
FIG. 7(a) is a graph showing the change in load power of the electrolytic aluminum in response to the wind wave fluctuation of FIG. 3 in accordance with one embodiment of the present invention;
FIG. 7(b) is a graph showing the change of the load power of the submerged arc furnace in response to the windage wave fluctuation of FIG. 3 according to an embodiment of the present invention;
FIG. 7(c) is a graph of polysilicon load power variation in response to the wind electrical fluctuations of FIG. 3 in accordance with one embodiment of the present invention;
FIG. 8 is a power regulation penalty F for a typical industrial park of the present inventionIP-ΔPΣ；
FIG. 9 is a graph of cost control characteristics for three industrial parks in example two of the present invention;
FIG. 10 is a wind power fluctuation plot for one minute for one area according to the second embodiment of the present invention;
FIG. 11 is a cost control characteristic of the utility company for the load aggregation of the industrial park of example two, in accordance with the present invention;
fig. 12 is a graph of the control cost of the electricity selling company responding to the wind power fluctuation of fig. 9 in the second 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, and the power grid response control method for the power selling companies based on the load aggregation of the industrial park comprises the following steps of:
step 1, carrying out load power normalization processing on an industrial park according to the power characteristics of electrolytic aluminum, a submerged arc furnace and a polycrystalline silicon load, and analyzing the power theoretical regulation boundary of each load;
step 2, a proportional type aggregation model which gives priority to the reporting regulation and control range of the user is provided in combination with the load regulation willingness and the limit regulation range, and an aggregation control strategy of the power selling company to the industrial park is established;
step 3, establishing a response strategy of the power selling company to the power grid control instruction according to the aggregation model, and realizing reasonable distribution and matching control of the load power of the industrial park;
step 4, analyzing the influence of the control strategy on the load production benefit, establishing a load control cost model, and establishing an industrial park load control cost model by combining a polymerization model;
and 5, based on the industrial park control cost model, proposing a bidding and preferential coordination strategy of the power selling company for multiple industrial parks, and realizing the coordinated distribution and control of the power of the multiple parks.
The steps form an aggregation control strategy for responding to the power grid by the power selling company. By aggregating the industrial park loads for the electricity selling companies,obtaining the load polymerization characteristic Delta P of an industrial park∑,t-ntAnd control cost model F of industrial parkIP-ΔPΣAdjusting the capacity (maximum upward adjustment capacity Δ P)Σ,up maxMaximum turndown capacity Δ PΣ,down max) And the control cost of the electricity selling company (F)PSC-PPSC) Reporting to the power grid, and regulating the power (electrolytic aluminum rectification current value I) by a bidding preferred coordination method according to the power grid regulation instructiondElectrode current command value I of submerged arc furnacerefTime of splicing wave ω t of polysilicon1Wave-splicing voltage U1、U2) And the data are transmitted to each industrial park and loads of different types, and finally the stabilization of power fluctuation is realized. As shown in fig. 1.
In the specific implementation, in the step 1, the numbers of the electrolytic aluminum load, the submerged arc furnace load and the polycrystalline silicon load in the industrial park are respectively NAL、NSAF、NPCSTaking the rated voltage of each load as a voltage reference value, and recording the voltage quantity of each load as U*The minimum limit and the maximum limit of each load voltage are respectively recorded as
Figure BDA0002266180100000211
For electrolytic aluminum load (i ═ {1,2, …, N)AL}) there are:
PAL,i＝m1,i·U*2+m2,iU*(30)
Figure BDA0002266180100000221
in the formula, PAL,iThe power of the ith electrolytic aluminum load,
Figure BDA0002266180100000222
direct current bus voltage rating, R, for the ith cell for aluminum electrolysisEC,iThe electrolytic cell for the ith electrolytic aluminum is connected with an equivalent resistance in series, EiThe equivalent potential of the electrolytic cell for the ith electrolytic aluminum.
Loaded with electrolytic aluminiumPer unit value U of DC bus voltage*The limit variation range of (2) is:
Figure BDA0002266180100000223
to submerged arc furnace load (i ═ N)AL+1,NAL+2,…,NAL+NSAF})：
PSAF,i＝m1,i·U*2(33)
QSAF,i＝ni·U*2(34)
Figure BDA0002266180100000224
Figure BDA0002266180100000225
In the formula, PSAF,i、QSAF,iRespectively the load active power and the load reactive power of the ith submerged arc furnace,
Figure BDA0002266180100000226
is the low side voltage rating of the ith furnace,
Figure BDA0002266180100000227
for the power characteristic parameter, R, of the ith submerged arc furnace during rated operationline,iAnd Xline,iIs the equivalent resistance and the equivalent reactance, R, of the ith submerged arc furnace short netarc,iAnd Xarc,iThe arc static resistance and the arc static reactance of the ith submerged arc furnace.
Per unit value U of voltage at low voltage side of submerged arc furnace*The limit variation range of (2) is:
Figure BDA0002266180100000228
for polysilicon loading (i ═ N)AL+NSAF+1,NAL+NSAF+2,…,NAL+NSAF+NPCS})：
PPCS,i＝m1,i·U*2(38)
Figure BDA0002266180100000231
In the formula, PPCS,iThe power is loaded for the ith poly-silicon,
Figure BDA0002266180100000235
for the ith polysilicon load single phase voltage rating, RPCS,iThe resistance is the single-phase resistance of the ith polycrystalline silicon rod.
The limit variation range of the per unit value of the polysilicon load single-phase voltage is as follows:
Figure BDA0002266180100000232
the power regulation characteristics of typical industrial park load electrolytic aluminum, submerged arc furnace, polysilicon are first modeled. The industrial park comprises an electrolytic aluminum plant with a rated power of 410MW for a period of electrolytic aluminum load; a submerged arc furnace load with the rated power of 76 MW; a polysilicon factory has a load rating of 122 MW. The three types of typical load parameters are shown in tables 1,2 and 3 below.
TABLE 1 electrolytic aluminum load parameters
Figure BDA0002266180100000233
In the table, VAL-AHNHigh voltage busbar voltage rating, k, for electrolytic aluminum loadsALNFor on-load tap-changing transformer ratio rating, VBNIs the DC bus voltage rating of the cell, E is the cell equivalent potential, RECEquivalent resistance, L, for cells in seriesSRNThe equivalent value of the saturable reactor.
TABLE 2 submerged arc furnace load parameters
Figure BDA0002266180100000234
Figure BDA0002266180100000241
In the table, VSAF-AHNIs a voltage rated value k of a high-voltage bus of the submerged arc furnaceSAFNIs a transformer ratio rated value, U, special for a submerged arc furnaceSAFNFor the rated voltage value of the low-voltage side of the submerged arc furnace, Rline+jXlineFor equivalent reactance rating of short networks, RarcN+jXarcNIs the equivalent impedance rating of the arc, IrefNFor rated values of command values of electrode current, RarcIs an arc equivalent resistance value, Xarc＝a1Rarc 2+a2Rarc+a3Is the arc impedance relation of the submerged arc furnace, PSAF,t minIs the minimum power of the submerged arc furnace, cos psiarcIs the arc portion power factor.
Table 3 polysilicon loading parameters (6 each silicon rod radius r 5,10, …,75mm)
Figure BDA0002266180100000242
In table I2＝11.711r3+0.765·r2The single-phase current variation characteristic of the polysilicon load under the rated condition along with the radius,
Figure BDA0002266180100000243
is the characteristic of single-phase voltage variation with radius of the polysilicon load under rated conditions, PPCS33.942 r +2.217 is the characteristic that the total AC power of the polysilicon load changes with the radius under the rated condition, η is the ratio of reaction heat absorption, K is the total heat transfer coefficient of the silicon rod and the mixed gas, L is the total length of the silicon rod, T is the total length of the silicon rodoutIs equivalent temperature, T, of the surface of the base plate and the furnace wallXAnd TXnα and α for the surface temperature of the silicon rod and its nominal valueNThe rate and its nominal value are regulated for the cooling water flow rate.
The load power characteristic parameters are obtained by substituting the parameters into equations (31), (35), (36), and (39). And the total load power characteristic is obtained from the following equation, as shown in fig. 2.
P∑＝M·U*2+M'·U*(41)
Each load independently reports voltage change range participating in regulation
Figure BDA0002266180100000255
Voltage regulation dead band
Figure BDA0002266180100000256
As shown in table 4 below.
TABLE 4 load Voltage Limit Change constraints and autonomous reporting parameters
Figure BDA0002266180100000251
In the table, the number of the first and second,
Figure BDA0002266180100000257
is the minimum value of the limit of the load voltage,
Figure BDA00022661801000002510
the maximum value of the limit voltage of the load,
Figure BDA0002266180100000258
reporting the minimum value of the voltage for the load,
Figure BDA0002266180100000259
reporting the maximum value of the voltage for the load,
Figure BDA00022661801000002511
adjusting the dead zone for the load voltage; the polysilicon load data is the parameter r for 75mm silicon rod.
The voltage of each load is actually adjustable to a minimum/maximum value of
Figure BDA0002266180100000252
The actual maximum allowable up/down adjustment amount of each load voltage is respectively as follows:
Figure BDA0002266180100000253
therefore, the range of the actual up/down adjustable capacity for each load can be calculated from equations (42), (43), see table 5 below.
TABLE 5 actual Up/Down Adjustable Voltage, Power of the load
Figure BDA0002266180100000254
In the table, the number of the first and second,
Figure BDA00022661801000002512
the maximum upward adjustment is actually allowed for the load,
Figure BDA00022661801000002513
for the maximum downward adjustment actually allowed for the load, Δ Pload,i,up maxFor the actual upward adjustable capacity of the load, Δ Pload,i,down maxThe capacity is actually adjusted downwards for the load.
In step 2, it is assumed that each load is in a rated operation state at the beginning (when not participating in regulation), i.e. U * i,01. the variation between the load voltage at time t and the initial voltage is recorded as Δ U* i,tI.e. by
ΔU* i,t＝U* i,t-U* i,0＝U* i,t-1 (44)
The load voltage adjustment amount at the adjacent time is:
U* i,t-U* i,t-1＝ΔU* i,t-ΔU* i,t-1(45)
namely, it is
Figure BDA0002266180100000261
In the formula, ntFor proportional adjustment coefficients, each time of adjustmentThe adjusting coefficients of all the loads are consistent, and the actual adjusting quantity is in positive correlation with the adjustable range of each load.
Setting the power regulation quantity (relative to the initial power) of the ith load in the park at the time t as delta Pload,i,t＝Pload,i,t-Pload,i,0The total power regulation quantity (relative initial power) of the load of the park is delta P∑,t＝P∑,t-P∑,0。
The aggregate characteristics of the individual loads can be obtained from the load regulation characteristic equations (30), (33), (38):
Figure BDA0002266180100000262
and also
Figure BDA0002266180100000263
The load polymerization characteristics of the industrial park were:
Figure BDA0002266180100000264
wherein the content of the first and second substances,
Figure BDA0002266180100000271
the load polymerization characteristics DeltaP of the corresponding zone were obtained from the formulas (48) and (49)∑,t-ntAs shown in fig. 3.
Let n betDetermining the maximum upward and downward regulating capacity DeltaP of the industrial park as 1 or-1Σ,up max、ΔPΣ,down max22.3415MW and-80.3150 MW.
In step 3, for the above example, wind power fluctuation occurs in 1 minute on the contact line between the industrial park and the local power grid as shown in fig. 4. The electricity selling company utilizes the total power regulating quantity delta P of the industrial park∑,tIn response to the tie-line power fluctuations, i.e.
ΔP∑,t＝ΔPt net＝Pt net-Pt-1 net(50)
The power regulation quantity Delta P is obtained from the equation (48)∑,tCorresponding aggregate instruction nt：
Figure BDA0002266180100000272
From equations (44) and (46), the voltage control target of each load is:
Figure BDA0002266180100000273
the aggregate command at each time of the power selling company can be obtained by the equation (51), the target voltage of each load can be obtained by the equation (52), and the actual power of each load can be obtained by the equations (1), (7), (8), (20) and (22).
The aggregate command and the change in the load voltage in response to the wind electrical fluctuation in fig. 3 are shown in fig. 5(a), fig. (b), fig. (c), and fig. (d). The tie line fluctuation graphs before and after the aggregation control strategy are adopted are shown in fig. 6(a) and (b), the tie line power fluctuation at each moment is reduced to be near-0.3882 MW, and the stabilizing effect on wind power fluctuation is obvious.
The power change of each load when the aggregation control strategy is adopted to respond to the power fluctuation is shown in fig. 7(a), fig. 7(b) and fig. 7(c), and the power range with adjustable load is marked in the graph and is shown by a dotted line in the graph. It can be seen from the figure that each load strictly meets the requirement of autonomous reporting of the load at the control moment, the regulation capacity of the load is not exceeded, and the adjustable capacity of each load is fully utilized.
In step 4, the production benefits of each type of load are analyzed first. According to different directions of the load participation of the industrial park in the upward/downward power regulation, the load control cost is divided into:
① when load power adjustment amount delta PloadWhen the power is less than 0, the energy loss is caused by power regulation, the load yield is influenced, and the load value loss, the unit value loss F of the load value loss is causedvCalculated from the following formula:
Fv＝(Fp-Fc)/CE(59)
wherein FvValue loss per unit load (yuan/kilowatt-hour), FpSelling price per unit load (yuan/ton), FcProduction cost per unit load (yuan/ton), CEIs the power consumption per unit of production (kilowatt-hour/ton).
② when load power adjustment amount delta PloadWhen the energy is more than 0, although the input energy of the load is increased, the energy cannot correspondingly increase the load due to the capacity limit of the load equipment, namely, the electric energy cost of the load is additionally increased, namely the electricity price PrE(Yuan/kilowatt-hour).
In addition, since the overload operation under the high power condition increases the failure rate of the device and affects the life of the device, the loss of the overload operation will be modeled based on the life cycle model.
The maintenance cost of the load is set as lambdai(element) rated service life cycle of τi,N(h) In that respect When the load recording operation is carried out at the time t, the service life is changed to tau due to overload operationi,t(h) The amount of change ρ in the hourly converted maintenance cost of the loadi,t(m/h) can be calculated by:
Figure BDA0002266180100000281
recording the load power as PN,i(kW), the reduced maintenance cost per unit power consumption of the load is maintenance unit price Fre(yuan/kwh) can be solved by:
Figure BDA0002266180100000291
in summary, the control cost F of the unit electric quantity of the load0(yuan/kwh) is:
Figure BDA0002266180100000292
the economic parameters and control costs for a typical campus load control are shown in table 6 below.
TABLE 6 economic parameters for high energy consumption load control in typical parks
Figure BDA0002266180100000293
Power regulation cost model F of industrial parkIP,t-ΔPΣ,tThe following formula:
Figure BDA0002266180100000294
wherein the content of the first and second substances,
Figure BDA0002266180100000295
the power regulation cost F of the typical industrial park can be solved by substituting the parameters (62) and (67) in the table 6IP-ΔPΣAs shown in fig. 8.
In step 5, recording the number of the management industrial parks of the power selling company as NIPThe industrial parks are sorted according to the minimum price of the times of power up regulation, and the serial numbers are recorded as u1 and u2 … uNIP(ii) a The times of price minimum sequencing of each industrial park are regulated downwards according to power, and the serial numbers are recorded as d1 and d2 … dNIPThus, the control cost model F of the electricity selling companyPSC-ΔPPSCComprises the following steps:
Figure BDA0002266180100000301
wherein FPSCIs the control cost, Δ P, of the electricity selling companyPSCIs the sum of the total regulated power of the power selling company, i.e. the regulated power of each industrial park, delta PΣ,up ui,max、ΔPΣ,down di,maxThe maximum upward regulation capacity of the industrial park numbered ui and the maximum downward regulation capacity of the industrial park numbered di, respectively.
According to the bidding preference strategy, the power of each industrial park is regulated according to the regulation requirement in the order of controlling the cost from low to high. The relationship between the total regulated power of the power selling company and the power of each park is as follows:
Figure BDA0002266180100000302
three industrial parks are managed by an electricity selling company, and the high energy consumption load conditions of each park are as follows:
a first park: the first-stage electrolytic aluminum load is rated at 410 MW; the first-stage submerged arc furnace load is rated at 76 MW; the load of the first stage electrolytic aluminum is 122MW rated power, and the load is divided into groups according to the radius of the produced silicon rod, (r is 5,10, 15, … and 75mm), and each group comprises 6 silicon rods.
And a second park: the two-stage electrolytic aluminum load has the rated power of 410MW and 620MW respectively; the first-stage submerged arc furnace is loaded and rated at 76 MW.
And a third park: the two-stage submerged arc furnace is loaded, and the rated power is 76MW and 90MW respectively; the two-stage electrolytic aluminum load has rated power of 122MW and rated power of 203MW respectively, and is divided into groups according to the radius of the produced silicon rod, wherein r is 5,10, 15, … and 75mm, and each group comprises 6 silicon rods and 8 silicon rods.
The power adjustment cost F of each campus is solved by equation (67)IP-ΔPΣAs shown in fig. 9. The maximum upward regulating capacity of the industrial park is 22.335MW, and the maximum downward regulating capacity is-80.315 MW; the regulation capacity of the industrial park II is large, the maximum upward regulation capacity is 52.6303MW, the maximum downward regulation capacity is-136.4597 MW, but the control cost is relatively large; the maximum upward regulating capacity of the third industrial park is 13.3976MW, the maximum downward regulating capacity is-93.9124 MW, and the control cost is relatively small.
The tie-line fluctuation caused by wind power fluctuation in 1 minute on the tie-line of the ground and the local power grid is shown in fig. 10.
The power selling company adopts a multi-park coordination control strategy shown in a formula (70), namely when the power is regulated upwards, the park with low control cost is regulated preferentially, namely the park is regulated according to the sequence of the park I, the park II and the park III; when the power is adjusted downwards, the park with low control cost is adjusted preferentially, namely, the parks are adjusted according to the sequence of three parks, one park and two parks. Control cost model F capable of solving electricity selling companyPSC-ΔPPSCAs shown in fig. 11.
The wind power fluctuation shown in fig. 11 is stabilized, and the control cost of the power selling company at each moment in the control time of 60s is shown in fig. 12 (yuan/MWh). The total amount of power for power stabilization in 60s control time is 2.1765MWh, and the total control cost is 189.1082 yuan. The average control cost was 86.8869M/MWh.
In summary, the present embodiment provides a limit adjustment range of each load based on modeling the power adjustment characteristic of the load in a typical industrial park; in combination with the load regulation will, a proportional aggregation model which gives priority to the reporting regulation and control range of the user is provided, and thus an aggregation control strategy of the power selling company for the industrial park is established to respond to the power grid frequency modulation requirement; establishing an industrial park control cost model by quantifying the influence of various types of load control on production benefits; and (4) providing a coordination control strategy for the power selling company to respond to the participation demands of the plurality of parks by combining the cost and the aggregation model.
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 (6)

1. A power selling company response power grid control method based on industrial park load aggregation is characterized by comprising the following steps:
step 1, carrying out load power normalization processing on an industrial park according to the power characteristics of electrolytic aluminum, a submerged arc furnace and a polycrystalline silicon load, and analyzing the power theoretical regulation boundary of each load;
step 2, a proportional type aggregation model which gives priority to the reporting regulation and control range of the user is provided in combination with the load regulation willingness and the limit regulation range, and an aggregation control strategy of the power selling company to the industrial park is established;
step 3, establishing a response strategy of the power selling company to the power grid control instruction according to the aggregation model, and realizing reasonable distribution and aggregation control of the load power of the industrial park;
step 4, analyzing the influence of the aggregation control strategy on the load production benefit, establishing a load control cost model, and establishing an industrial park load control cost model by combining the aggregation model;
and 5, based on the industrial park load control cost model, providing a bidding and preferential coordination strategy of the power selling company for multiple industrial parks, and realizing the coordinated distribution and control of the power of the multiple parks.
2. The electric power selling company response power grid control method based on industrial park load aggregation as claimed in claim 1, wherein the concrete implementation of the step 1 comprises:
step 1.1, the load power characteristic of the electrolytic aluminum is as follows:
Figure FDA0002266180090000011
wherein, PALFor power of electrolytic aluminum, VBIs the DC bus voltage of the electrolyzer IdIs direct current of the electrolytic cell, RECIs the equivalent resistance of the electrolytic cell in series connection, and E is the equivalent potential of the electrolytic cell;
electrolytic aluminum DC bus voltage VBVoltage V of high-voltage bus connected with electrolytic aluminium loadAL-AHThe relationship of (1) is:
Figure FDA0002266180090000012
in the formula, LSRIs the equivalent value of the saturable reactor, VAL-AHIs a high voltage bus voltage, and omega is a voltage angular frequency;
considering that the electrolytic aluminum load is connected with an on-load tap changer, the transformation ratio value of the on-load tap changer is m1Stage, then change its transformation ratio kALCan realize m1Stage adjustment:
Figure FDA0002266180090000021
saturable reactor equivalent value L when production requirements are metSRHas a regulation range of [ LSR min,LSR max]Then, from the formula (2), the strain ratio k can be obtainedAL-iVoltage drop V of direct current bus of lower electrolytic aluminumBThe adjusting range is as follows:
Figure FDA0002266180090000022
Figure FDA0002266180090000023
the maximum adjustable power range of the electrolytic aluminum load obtained by the formula (1) is as follows:
Figure FDA0002266180090000024
step 1.2, the power characteristics of the submerged arc furnace load are as follows:
Figure FDA0002266180090000025
Figure FDA0002266180090000026
wherein, PSAF、QSAFRespectively active power and reactive power of the submerged arc furnace, USAFIs the voltage of the low-voltage side of the submerged arc furnace, RlineAnd XlineThe equivalent resistance and the equivalent reactance of the short net;
the equivalent impedance of the arc is defined by an arc static resistance R in a certain time break planearcAnd arc static reactance XarcTo characterize;
low-voltage side voltage U of submerged arc furnaceSAFThe voltage U of a high-voltage bus of the submerged arc furnaceSAF-AHObtained by a special transformer for a submerged arc furnace:
USAF＝USAF-AH/kSAF(9)
in the formula, kSAFThe transformer is a transformer special for the submerged arc furnace;
the arc impedance relational expression of the submerged arc furnace is obtained by fitting actual production data:
Xarc＝a1Rarc 2+a2Rarc+a3(10)
wherein, a1、a2、a3Is the arc impedance fitting coefficient, which is a constant;
when the power of the submerged arc furnace is adjusted by adopting a constant impedance voltage-adjusting mode, the transformation ratio of the transformer special for the submerged arc furnace load is set to be m in total2A stage; the voltage of the transformer is m in total by changing the transformation ratio of the transformer special for the submerged arc furnace2Stage regulation range:
Figure FDA0002266180090000031
arc static resistance RarcThe limiting range of (2):
Rarc min≤Rarc≤Rarc max(12)
upper limit of arc partial power factor
Figure FDA0002266180090000032
There is a lower limit to the arc segment power factor
Figure FDA0002266180090000033
Arc partial power factor
Figure FDA0002266180090000034
The limitations of (2) are:
Figure FDA0002266180090000035
namely, it is
Figure FDA0002266180090000036
Wherein the content of the first and second substances,
Figure FDA0002266180090000037
meanwhile, each smelting stage should meet the minimum power constraint:
PSAF,t≥PSAF,t min(15)
in the formula, PSAF,t minThe minimum power of the submerged arc furnace;
the rated voltage of the low-voltage side of the submerged arc furnace is set as USAF,NBy the formulas (7), (8), (12) and (15), when the constant voltage impedance power regulation method is adopted, the active power P of the submerged arc furnace in each smelting stageSAF,tAnd reactive power QSAF,tThe adjustment range of (a) is as follows:
c1≤PSAF,t≤c2(16)
d1≤QSAF,t≤d2(17)
wherein:
Figure FDA0002266180090000041
Figure FDA0002266180090000042
Figure FDA0002266180090000043
Figure FDA0002266180090000044
the power characteristic parameter of the ore furnace during rated operation is
Figure FDA0002266180090000045
Adjusting the power at the momentInto an equivalent low-side voltage variation range, i.e.
Figure FDA0002266180090000046
From the equations (10) and (18), when the impedance-voltage coordinated regulation method is adopted, the voltage variation range of the low-voltage side of the submerged arc furnace load is as follows:
Figure FDA0002266180090000047
step 1.3, the polysilicon load power meets the following electrical relationship:
Figure FDA0002266180090000048
in the formula, PPCSFor polysilicon load AC total power, UvalFor loading polycrystalline silicon with a single-phase voltage, RPCSIs a polysilicon rod single-phase resistor;
the production process energy conversion of a polycrystalline silicon rod in the time delta t is as follows:
Figure FDA0002266180090000051
in the formula, PPCSFor polysilicon load AC total power, Δ Qout1V represents the amount of heat used to heat the reactant gas, and is given by the gas specific heat capacity formula1·Δt·s1·ρg·c·(Tx-Tg) Wherein v is1、s1、ρg、c、TgRespectively, the air inlet rate, the air inlet area, the mixed gas density, the mixed gas specific heat capacity and the air inlet temperature are constants; t isxSurface temperature, Δ Q, of the silicon rodout2And Δ Qout3Respectively, the heat quantity of the heat of the endothermic reaction maintained and dissipated through the base and the wall of the reduction furnace by thermal radiation
Figure FDA0002266180090000052
Corresponds to (Δ Q)out2+ΔQout3) Wherein, η, K, L, ToutThe reaction heat absorption ratio, the total heat transfer coefficient of the silicon rod and the mixed gas, the total length of the silicon rod, and the equivalent temperature of the chassis and the furnace wall surface are all constants; r is the radius of the polysilicon rod, and r can be regarded as a constant in a short time;
the power characteristic equation of the polysilicon load with radius r obtained from equations (20) and (21) is as follows:
Figure FDA0002266180090000053
in the formula, A, B, C, D, G, H is a fitting coefficient of the power characteristic of the polysilicon, is a constant and is obtained by fitting actual production rated operation data; i is polysilicon load single-phase current;
for the surface temperature T of the silicon rodxThe control range is as follows:
Tx min≤Tx≤Tx max(23)
when the temperature is less than or equal to 1000 ℃ TxCan ensure production at the temperature of less than or equal to 1100 ℃, and can ensure the production at Tx＝Tx，NThe temperature is the optimum temperature at 1080 ℃; when engaged, the polysilicon load generally engages in downward regulation of power, then there is Tx min＝1000℃,Tx max＝1080℃；
The flow rate of the cooling water is generally regulated, and the flow rate regulation rate of the cooling water is α
Figure FDA0002266180090000054
αmin≤α≤αmax(25)
Wherein, αmin＝90％，αmaxWhen the operation is rated, α is 100%;
the following equations (20) and (21) can be obtained:
Figure FDA0002266180090000061
from equations (23) and (26), the adjustment range of the polysilicon load power can be determined as follows:
Figure FDA0002266180090000062
Figure FDA0002266180090000063
effective value U of single-phase voltage of polysilicon loadvalThe adjusting range of (A) is as follows:
Figure FDA0002266180090000064
step 1.4, carrying out load normalization processing on the industrial park;
setting the number of electrolytic aluminum load, submerged arc furnace load and polysilicon load in the industrial park as N respectivelyAL、NSAF、NPCSThe rated voltage of each load is taken as a voltage reference value, and the voltage quantity of each load is U*The minimum and maximum limit values of each load voltage are respectively
Figure FDA0002266180090000065
Step 1.4.1. for electrolytic aluminum load i ═ {1,2, …, NALHas the following components:
PAL,i＝m1,i·U*2+m2,iU*(30)
Figure FDA0002266180090000066
in the formula, PAL,iPower of the i-th electrolytic aluminum load, VBN,iDirect current bus voltage rating, R, for the ith cell for aluminum electrolysisEC,iThe electrolytic cell for the ith electrolytic aluminum is connected with an equivalent resistance in series, EiThe cell equivalent potential for the ith electrolytic aluminum;
DC bus voltage per unit value U of electrolytic aluminum load*The limit variation range of (2) is:
Figure FDA0002266180090000071
step 1.4.2. for ore furnace load i ═ NAL+1,NAL+2,…,NAL+NSAF}：
PSAF,i＝m1,i·U*2(33)
QSAF,i＝ni·U*2(34)
Figure FDA0002266180090000072
Figure FDA0002266180090000073
In the formula, PSAF,i、QSAF,iRespectively the load active power and the load reactive power of the ith submerged arc furnace,
Figure FDA0002266180090000074
is the low side voltage rating of the ith furnace,
Figure FDA0002266180090000075
for the power characteristic parameter, R, of the ith submerged arc furnace during rated operationline,iAnd Xline,iIs the equivalent resistance and the equivalent reactance, R, of the ith submerged arc furnace short netarc,iAnd Xarc,iThe resistance and the reactance of the arc static resistance of the ith submerged arc furnace are shown;
per unit value U of voltage at low voltage side of submerged arc furnace*The limit variation range of (2) is:
Figure FDA0002266180090000076
step 1.4.3. load i ═ N to polysiliconAL+NSAF+1,NAL+NSAF+2,…,NAL+NSAF+NPCS}：
PPCS,i＝m1,i·U*2(38)
Figure FDA0002266180090000077
In the formula, PPCS,iThe power is loaded for the ith poly-silicon,
Figure FDA0002266180090000078
for the ith polysilicon load single phase voltage rating, RPCS,iThe resistance is the single-phase resistance of the ith polycrystalline silicon rod;
the limit variation range of the per unit value of the polysilicon load single-phase voltage is as follows:
Figure FDA0002266180090000081
step 1.4.4. Total load Power characteristic P of Industrial parkΣ-U*Comprises the following steps:
PΣ＝M·U*2+M'·U*(41)
wherein the content of the first and second substances,
Figure FDA0002266180090000082
m and M' are the quadratic coefficient and the first-order coefficient of the total load power characteristic, respectively, M1,iAnd m2,iThe quadratic coefficient and the first order coefficient for each type of load power characteristic are obtained by equations (31), (35), and (39).
3. The electric power selling company response power grid control method based on industrial park load aggregation as claimed in claim 2, wherein the implementation of step 2 comprises the following steps:
step 2.1, giving priority to a proportional aggregation model of a user reporting regulation and control range:
let a burdenThe load reports the voltage change range hoped to participate in regulation independently as
Figure FDA0002266180090000083
Voltage regulation dead zone of
Figure FDA0002266180090000084
Considering the voltage limit variation range of each type of load determined by the equations (6), (19) and (40), the voltage variation range in which the load actually participates in regulation is:
Figure FDA0002266180090000085
the voltage of each load can be adjusted to a practically adjustable minimum/maximum value of
Figure FDA0002266180090000086
The actual maximum allowable up/down adjustment amount of each load voltage is respectively as follows:
Figure FDA0002266180090000091
when power regulating quantity distribution is carried out, voltage is regulated among all loads in equal proportion according to the voltage regulation range;
when the system is initially set and is not involved in regulation, all loads are in a rated operation state, namely U* i,01 is ═ 1; the variation between the load voltage at time t and the initial voltage is recorded as Δ U* i,tI.e. by
ΔU* i,t＝U* i,t-U* i,0＝U* i,t-1 (44)
The load voltage adjustment amount at the adjacent time is:
U* i,t-U* i,t-1＝ΔU* i,t-ΔU* i,t-1(45)
namely, it is
Figure FDA0002266180090000092
In the formula, ntThe load adjusting parameters are proportional adjusting parameters, the adjusting parameters of each load are consistent during each adjustment, and the actual adjusting quantity is in positive correlation with the adjustable range of each load;
setting the power regulating quantity of the ith load in the park at the time t as delta P relative to the initial powerload,i,t＝Pload,i,t-Pload,i,0The total power regulation quantity of the load in the park is delta P relative to the initial powerΣ,t＝P∑,t-P∑,0；
The aggregate characteristics of the respective loads can be obtained from the load regulation characteristics (30), (33), (38):
Figure FDA0002266180090000093
and also
Figure FDA0002266180090000101
The load polymerization characteristics of the industrial park were:
Figure FDA0002266180090000102
wherein the content of the first and second substances,
Figure FDA0002266180090000103
n is to betSubstituting 1 into the expressions (47) and (48) can respectively obtain the maximum upward adjustment capacity delta P of each load and industrial parkload,i,up max、ΔPΣ,up max(ii) a N is atThe maximum downward regulating capacity delta P of each load and industrial park can be respectively obtained by substituting the value of-1load,i,down max、ΔPΣ,down max；
Step 2.2, the power selling company responds to the aggregation algorithm of the power grid requirement;
step 2.2.1, an off-line polymerization stage;
step 2.2.1.1, solving the load power characteristic of the industrial park: respectively determining power characteristics of electrolytic aluminum, a submerged arc furnace and a polycrystalline silicon load according to formulas (1), (7), (8), (20) and (22) based on the measured data;
step 2.2.1.2. solving the total load power characteristic: obtaining the load power characteristic parameters of the expressions (31), (35), (36) and (39) and obtaining the total load power characteristic of the expression (41);
step 2.2.1.3, solving the voltage limit regulation range: the voltage limit variation range of each load is obtained by combining equations (4), (5), (19) and (29) with equations (32), (37) and (40)
Figure FDA0002266180090000104
Step 2.2.1.4, determining the actual voltage regulation range: recording the voltage regulation range and the regulation dead zone reported by each load independently, and solving the actual maximum upward regulation quantity of the voltage of each load by the formulas (42) and (43);
step 2.2.1.5. calculating the load polymerization characteristics and the maximum regulating capacity of the industrial park: the load polymerization characteristics DeltaP of the industrial park are obtained from the formulae (48) and (49)∑,t-nt(ii) a Let n betDetermining the maximum upward and downward regulating capacity DeltaP as 1 or-1Σ,up max、ΔPΣ,down max；
Step 2.2.2, an on-line calculation stage;
setting the power of the industrial park to be delta P at the previous moment of load controlΣ,t-1Initial adjustment time Δ PΣ,0When the maximum upward and downward regulation capacity of the industrial park is equal to 0, the maximum upward and downward regulation capacity of the industrial park at the current time t is delta PΣ,up,t max＝ΔPΣ,up max-ΔPΣ,t-1；
Step 2.2.3, reporting the adjustable capacity of the power grid at the current moment;
the electricity selling company adjusts the maximum upward and downward capacity delta P at the momentΣ,up,t max、ΔPΣ,down,t maxAnd uploading to a power grid dispatching center.
4. The electric utility company response power grid control method based on industrial park load aggregation as claimed in claim 3, wherein the implementation of step 3 comprises the following steps:
step 3.1, a load distribution principle;
when the power grid issues a control target power P to the power selling companyt netThen, the power selling company distributes power to each load in the industrial park;
setting the total power regulating quantity of the industrial park at t moment as delta P∑,tThe solution is as follows:
ΔP∑,t＝ΔPt net＝Pt net-Pt-1 net(50)
the power regulation quantity Delta P is obtained from the equation (48)∑,tCorresponding aggregate instruction nt：
Figure FDA0002266180090000111
From equations (44) and (46), the voltage control target of each load is:
Figure FDA0002266180090000112
the actual power of each load after participating in the regulation can be obtained according to the formulas (1), (7), (8) and the formulas (20) and (22) respectively;
step 3.2, a load control instruction is carried out;
step 3.2.1, controlling the load of the electrolytic aluminum;
the direct current bus voltage V of the electrolytic cellB,i,t＝VBN,i·U* i,tSubstituting into formula (2) can solve the saturated reactor equivalent value L of the electrolytic aluminum loadSR,iComprises the following steps:
Figure FDA0002266180090000121
direct current bus of electrolytic cellVoltage VB,i,t＝VBN,i·U* i,tSubstituting the formula (1) into the load to obtain the corresponding rectified current value I of the loadd,i,tComprises the following steps:
Figure FDA0002266180090000122
step 3.2.2, a control instruction of the submerged arc furnace load;
using electrode current command value IrefTo control the lifting of the electrode; the electrode hydraulic lifting model comprises the following steps:
Figure FDA0002266180090000123
wherein L is the distance between the electrode and the furnace charge, is equal to the arc length of the electric arc, and L0As the initial position of the electrode, IrefIs an electrode current command value; v. ofupAnd vdownMaximum values representing the electrode rising and falling speeds, both constant, cupAnd cdownIs a proportionality coefficient between the current difference and the speed;
when changing the arc impedance, correspondingly adjusting I in the lifting model according to the following formularefAnd the electrode lifting can be realized:
Figure FDA0002266180090000124
when the voltage of the low-voltage side of the submerged arc furnace is
Figure FDA0002266180090000125
The electrode control command is then obtained by:
Figure FDA0002266180090000131
step 3.2.3, controlling the polysilicon load;
from the equation (20), it can be seen that the polysilicon load single-phase voltage U is changedvalThe polysilicon load power can be realized(iii) adjustment of (c); controlling U by wave splicing technologyvalThe following are:
Figure FDA0002266180090000132
in the formula, the voltage angular frequency omega is constant, U1、U2To the splicing voltage, t1The wave splicing time;
when obtaining the polysilicon load single-phase voltage Uval,i,t＝UvalN,i·U* i,tThen, firstly, the wave splicing voltage U is determined1、U2: the value of the wave-splicing voltage is U1,U2E is (0,380,600,800,1500) V, and the fixed splicing wave voltage U is obtained1、U2Are respectively a target voltage Uval,i,tTwo adjacent voltage levels U1＜U2(ii) a Then put Uval,i,t、U1、U2Substitution formula (58) for calculating splicing wave voltage time t1(ii) a Will Uval,i,tSubstituting formula (26) and solving the surface temperature T of the silicon rodxCooling water flow rate α, wherein the surface temperature T of the silicon rod is maintained preferentially during regulationx＝1080℃。
5. The electric utility company response power grid control method based on industrial park load aggregation as claimed in claim 4, wherein the implementation of step 4 comprises the following steps:
step 4.1, controlling the influence of the strategy on the load production benefit;
according to different directions of the load participation of the industrial park in the upward/downward power regulation, the load control cost is divided into:
step 4.1.1, when load power regulating quantity delta PloadWhen less than 0, the loss of unit value FvCalculated from the following formula:
Fv＝(Fp-Fc)/CE(59)
wherein FvIs the loss of value per unit load, unit cell/kilowatt-hour, FpIs sold per unit load, unit/ton, FcProduction cost per unit load, unit/ton, CEIs a unit yieldPower consumption, unit kilowatt-hour/ton;
step 4.1.2 when load power regulating quantity delta PloadAnd when the load is more than 0, modeling the loss of the overload operation based on the life cycle model:
the maintenance cost of the load is set as lambdaiUnit cell, rated maintenance life cycle τi,NThe unit h; when the load operation is set at time t, the service life changes to tau due to the overload operationi,tUnit h, the amount of change ρ in the hourly converted maintenance cost of the loadi,tUnit cell/h, calculated by:
Figure FDA0002266180090000141
recording the load power as PN,iIn kW unit, the reduced maintenance cost of the unit power consumption of the load is the maintenance unit price FreUnit cell/kwh, is solved by:
Figure FDA0002266180090000142
in summary, the control cost F of the unit electric quantity of the load0Unit cell/kwh, is:
Figure FDA0002266180090000143
step 4.2, controlling a cost model in the industrial park;
setting the control cost of each high energy consumption load in an industrial park as F0,iWherein i is 1,2, … NAL+NSAF+NPCSThe electrolytic aluminum, the submerged arc furnace and the polycrystalline silicon are sequentially numbered according to the load sequence;
the control cost of each load is related to the adjustment power by the following formula:
Figure FDA0002266180090000144
wherein Fload,i,tThe control cost for the ith load;
power regulation of industrial park
Figure FDA0002266180090000145
Let the total control cost of the industrial park be FIP,t：
Figure FDA0002266180090000146
Solving the control cost F of the industrial parkIP,tAnd adjusting the electric quantity delta PΣ,tFunctional relationship F ofIP,t-ΔPΣ,tIn combination with the load distribution and characteristic Δ Pload,i,t-ntSolving the control cost of the industrial park-polymerization model FIP,tNt, and then the model F of the power regulation cost of the industrial park is solvedIP,t-ΔPΣ,tNamely:
Figure FDA0002266180090000151
the power adjustment for the ith load on the industrial park is determined by:
Figure FDA0002266180090000152
from the equations (65) and (66), the polymerization model F, which is the control cost of the industrial park, can be determinedIP,tNt is as follows:
Figure FDA0002266180090000153
wherein the content of the first and second substances,
Figure FDA0002266180090000154
when n istWhen 1, the cost F of providing the maximum upward adjustable power for the industrial park can be obtainedIP,t max(ii) a When n istWhen is equal to-1The cost F of providing the maximum downwardly adjustable power for the industrial park is obtainedIP,t min；
The power regulation cost model F of the industrial park can be obtained by the formulas (51) and (67)IP,t-ΔPΣ,tThe following formula:
Figure FDA0002266180090000155
wherein the content of the first and second substances,
Figure FDA0002266180090000161
6. the electric utility company response power grid control method based on industrial park load aggregation as claimed in claim 5, wherein the implementation of step 5 comprises the following steps:
step 5.1, a control cost model and a coordination control method of the power selling company;
the number of the management industrial parks of the electricity selling company is NIPThe industrial parks are sorted according to the minimum price of the times of power up regulation, and the serial numbers are recorded as u1 and u2 … uNIP(ii) a The times of price minimum sequencing of each industrial park are regulated downwards according to power, and the serial numbers are recorded as d1 and d2 … dNIPThus, the control cost model F of the electricity selling companyPSC-ΔPPSCComprises the following steps:
Figure FDA0002266180090000162
in the formula, FPSCIs the control cost, Δ P, of the electricity selling companyPSCIs the sum of the total regulated power of the power selling company, i.e. the regulated power of each industrial park, delta PΣ,up ui,max、ΔPΣ,down di,maxThe maximum upward regulating capacity of the industrial park numbered ui and the maximum downward regulating capacity of the industrial park numbered di are respectively;
according to a bidding preference strategy, the power of each industrial park is regulated in a sequence from low to high in control cost according to regulation requirements; the relationship between the total regulated power of the power selling company and the power of each park is as follows:
Figure FDA0002266180090000171
the total power of each industrial park can be determined by equation (70), and the aggregate demand n for each parktThe voltage target value of each load and the actual control command of each load can be obtained by the equation (51) and the equations (52), (54), (57), (58) and (26);
step 5.2, the coordination control algorithm of the power selling company to the multiple parks;
step 5.2.1, an off-line polymerization stage;
step 5.2.1.1. each park load reports the voltage regulation range independently
Figure FDA0002266180090000172
Calculating control cost by the formula (62);
step 5.2.1.2. evaluation of the formula (43) confirms the respective load control ranges, and calculation of the polymerization characteristics Δ P of the respective parks from the formula (48)Σ,t-ntCalculating a power regulation cost model F for each park from equation (67)IP-ΔPΣ；
Step 5.2.1.3, calculating the control cost model F of the power selling company by the formula (69)PSC-ΔPPSC；
Step 5.2.2, an on-line polymerization calculation stage;
solving the maximum upward/downward adjustable capacity delta P of each park at the current time t from the last control time t-1Σ,up,t ui,max＝ΔPΣ,up ui,max-ΔPΣ,t-1、ΔPΣ,down,t di,max＝ΔPΣ,down di,max-ΔPΣ,t-1And the control cost of the power selling company at the time t is FPSC,t＝FPSC(ΔPPSC-ΔPΣ,i)；
Step 5.2.2, reporting to the power grid; updating the adjustable capacity and cost model and reporting to the power grid;
step 5.2.3, receiving a power grid power regulation requirement;
step 5.2.3.1 receiving a grid power regulation demand Δ Pnet,tThen, the total power regulating quantity delta P of each industrial park is confirmed in sequence according to the formula (70) and the current adjustable capacity of each parkΣ,t；
Step 5.2.3.2, distributing the load power in each industrial park according to the step 3, and determining each control instruction;
and 5.2.3.3, recording and updating the current load state and waiting for the next adjustment.
CN201911088544.7A 2019-11-08 2019-11-08 Electricity selling company response power grid control method based on industrial park load aggregation Active CN110880772B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911088544.7A CN110880772B (en) 2019-11-08 2019-11-08 Electricity selling company response power grid control method based on industrial park load aggregation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911088544.7A CN110880772B (en) 2019-11-08 2019-11-08 Electricity selling company response power grid control method based on industrial park load aggregation

Publications (2)

Publication Number Publication Date
CN110880772A true CN110880772A (en) 2020-03-13
CN110880772B CN110880772B (en) 2023-04-18

Family

ID=69729242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911088544.7A Active CN110880772B (en) 2019-11-08 2019-11-08 Electricity selling company response power grid control method based on industrial park load aggregation

Country Status (1)

Country Link
CN (1) CN110880772B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112417792A (en) * 2020-11-23 2021-02-26 云南电网有限责任公司 Electrolytic aluminum external characteristic modeling method based on controllable boundary of saturable reactor
CN113162065A (en) * 2021-03-11 2021-07-23 国网山东省电力公司 Power grid frequency modulation method and system based on cooperative coordination between production loads
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
CN113852093A (en) * 2021-07-28 2021-12-28 国网山东省电力公司 Energy storage device and polysilicon load aggregation frequency modulation method and system
CN115187091A (en) * 2022-07-18 2022-10-14 东南大学溧阳研究院 Method for evaluating and controlling online response capability of air conditioner load
CN116131292A (en) * 2023-04-17 2023-05-16 武汉大学 Power grid frequency modulation method and system based on cooperative coordination of electrolytic aluminum and polysilicon
CN116191480A (en) * 2023-04-28 2023-05-30 武汉大学 Method and system for determining frequency modulation capacity of submerged arc furnace based on voltage regulation and electrode lifting

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106124826A (en) * 2016-06-17 2016-11-16 华中科技大学 A kind of mineral hot furnace working of a furnace digital management system, control method and monitoring device thereof
CN107425521A (en) * 2017-06-28 2017-12-01 国网江苏省电力公司电力科学研究院 A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method
CN109103912A (en) * 2018-07-18 2018-12-28 合肥工业大学 Consider the industrial park active distribution system method for optimizing scheduling of peaking demand of power grid
CN110112728A (en) * 2019-05-10 2019-08-09 四川大学 A kind of probabilistic more garden microgrid cooperative game methods of consideration wind-powered electricity generation robust

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106124826A (en) * 2016-06-17 2016-11-16 华中科技大学 A kind of mineral hot furnace working of a furnace digital management system, control method and monitoring device thereof
CN107425521A (en) * 2017-06-28 2017-12-01 国网江苏省电力公司电力科学研究院 A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method
CN109103912A (en) * 2018-07-18 2018-12-28 合肥工业大学 Consider the industrial park active distribution system method for optimizing scheduling of peaking demand of power grid
CN110112728A (en) * 2019-05-10 2019-08-09 四川大学 A kind of probabilistic more garden microgrid cooperative game methods of consideration wind-powered electricity generation robust

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIAN XU: "An Isolated Industrial Power System Driven by Wind-Coal Power for Aluminum Productions: A Case Study of Frequency Control" *
刘文颖: "考虑风电消纳的电力系统源荷协调多目标优化方法" *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112417792A (en) * 2020-11-23 2021-02-26 云南电网有限责任公司 Electrolytic aluminum external characteristic modeling method based on controllable boundary of saturable reactor
CN112417792B (en) * 2020-11-23 2023-05-26 云南电网有限责任公司 Electrolytic aluminum external characteristic modeling method based on controllable boundary of saturation reactor
CN113162065A (en) * 2021-03-11 2021-07-23 国网山东省电力公司 Power grid frequency modulation method and system based on cooperative coordination between production loads
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
CN113852093A (en) * 2021-07-28 2021-12-28 国网山东省电力公司 Energy storage device and polysilicon load aggregation frequency modulation method and system
CN115187091A (en) * 2022-07-18 2022-10-14 东南大学溧阳研究院 Method for evaluating and controlling online response capability of air conditioner load
CN116131292A (en) * 2023-04-17 2023-05-16 武汉大学 Power grid frequency modulation method and system based on cooperative coordination of electrolytic aluminum and polysilicon
CN116191480A (en) * 2023-04-28 2023-05-30 武汉大学 Method and system for determining frequency modulation capacity of submerged arc furnace based on voltage regulation and electrode lifting
CN116191480B (en) * 2023-04-28 2023-08-15 武汉大学 Method and system for determining frequency modulation capacity of submerged arc furnace based on voltage regulation and electrode lifting

Also Published As

Publication number Publication date
CN110880772B (en) 2023-04-18

Similar Documents

Publication Publication Date Title
CN110880772B (en) Electricity selling company response power grid control method based on industrial park load aggregation
Wang et al. Optimal sizing of distributed generations in DC microgrids with comprehensive consideration of system operation modes and operation targets
CN110400059B (en) Power control method for tie line of comprehensive energy system
CN103368173B (en) Active power flow optimizing distribution method for alternating-direct current parallel system containing soft direct current power transmission
CN111950808B (en) Comprehensive energy system random robust optimization operation method based on comprehensive demand response
CN111725841A (en) Photovoltaic inverter-based power quality optimization method for distribution transformer area
CN110729726B (en) Intelligent community energy optimization scheduling method and system
CN107863773B (en) Multi-energy-flow energy supply method and system
CN112990523B (en) Hierarchical optimization operation method for regional comprehensive energy system
CN107769234B (en) Frequency modulation capability evaluation method of multi-point energy storage coordination distributed power supply
CN110474370B (en) Cooperative control system and method for air conditioner controllable load and photovoltaic energy storage system
CN114169236A (en) Phase-change heat storage type electric heating negative control system control method based on LSTM algorithm
CN102751724A (en) Prediction-based three-phase load scheduling method and device responding to demand side
CN116365506A (en) Energy-saving and loss-reducing optimization method and system for active power distribution network containing distributed photovoltaic
CN117081143A (en) Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion
CN113765120A (en) Power grid frequency modulation method based on cooperative coordination of industrial park submerged arc furnace and polycrystalline silicon
CN116706943B (en) Multi-electrolytic aluminum load coordination control method based on segmented difference adjustment coefficient
Yu et al. Adaptive load control of electrolytic aluminum for power system frequency regulation based on the aluminum production operation state
Tziovani et al. Grid friendly operation of a PV-storage system with profit maximization and reliability enhancement
CN114825371B (en) Aggregation temperature control load multi-layer regulation and control method based on node voltage constraints before and after regulation
CN115879714A (en) Power system optimization scheduling method for electrolytic aluminum load participating in demand response
CN116316553A (en) Multi-time scale layered operation control method for hydrogen electric coupling system
CN115149540A (en) Multi-region virtual power plant collaborative optimization operation method based on partition autonomy
CN115483692A (en) Coordination control method and system based on electrolytic aluminum and polycrystalline silicon load polymerization
CN114399162A (en) Rolling optimization scheduling method based on energy scheduling time adaptive change

Legal Events

Date Code Title Description
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