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 PDFInfo
- 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
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
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 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:
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:
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:
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:
the maximum adjustable power range of the electrolytic aluminum load obtained by the formula (1) is as follows:
step 1.2, the power characteristics of the submerged arc furnace load are as follows:
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:
arc static resistance RarcThe limiting range of (2):
Rarc min≤Rarc≤Rarc max(12)
upper limit of arc partial power factorThere is a lower limit to the arc segment power factorArc partial power factorThe limitations of (2) are:
namely, it is
Wherein the content of the first and second substances,
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:
the power characteristic parameter of the ore furnace during rated operation isThe power regulation range at this time is converted into an equivalent low-voltage side voltage variation range, namely
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:
step 1.3, the polysilicon load power meets the following electrical relationship:
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:
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 radiationCorresponds 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:
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 α
αmin≤α≤αmax(25)
Wherein, αmin=90%,αmaxWhen the operation is rated, α is 100%;
the following equations (20) and (21) can be obtained:
from equations (23) and (26), the adjustment range of the polysilicon load power can be determined as follows:
effective value U of single-phase voltage of polysilicon loadvalThe adjusting range of (A) is as follows:
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
Step 1.4.1. for electrolytic aluminum load i ═ {1,2, …, NALHas the following components:
PAL,i=m1,i·U*2+m2,iU*(30)
in the formula, PAL,iThe power of the ith electrolytic aluminum load,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:
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)
In the formula, PSAF,i、QSAF,iRespectively the load active power and the load reactive power of the ith submerged arc furnace,is the low side voltage rating of the ith furnace,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:
step 1.4.3. load i ═ N to polysiliconAL+NSAF+1,NAL+NSAF+2,…,NAL+NSAF+NPCS}:
PPCS,i=m1,i·U*2(38)
In the formula, PPCS,iThe power is loaded for the ith poly-silicon,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:
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,
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 asVoltage regulation dead zone ofConsidering 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:
the voltage of each load can be adjusted to a practically adjustable minimum/maximum value of
The actual maximum allowable up/down adjustment amount of each load voltage is respectively as follows:
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
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):
and also
The load polymerization characteristics of the industrial park were:
wherein the content of the first and second substances,
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)
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:
From equations (44) and (46), the voltage control target of each load is:
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:
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:
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:
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:
when the voltage of the low-voltage side of the submerged arc furnace isThe electrode control command is then obtained by:
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:
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:
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:
in summary, the control cost F of the unit electric quantity of the load0Unit cell/kwh, is:
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:
wherein Fload,i,tThe control cost for the ith load;
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:
the power adjustment for the ith load on the industrial park is determined by:
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:
wherein the content of the first and second substances,
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:
wherein the content of the first and second substances,
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:
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:
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 independentlyCalculating 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 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
For electrolytic aluminum load (i ═ {1,2, …, N)AL}) there are:
PAL,i=m1,i·U*2+m2,iU*(30)
in the formula, PAL,iThe power of the ith electrolytic aluminum load,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:
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)
In the formula, PSAF,i、QSAF,iRespectively the load active power and the load reactive power of the ith submerged arc furnace,is the low side voltage rating of the ith furnace,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:
for polysilicon loading (i ═ N)AL+NSAF+1,NAL+NSAF+2,…,NAL+NSAF+NPCS}):
PPCS,i=m1,i·U*2(38)
In the formula, PPCS,iThe power is loaded for the ith poly-silicon,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:
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
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
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)
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,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 regulationVoltage regulation dead bandAs shown in table 4 below.
TABLE 4 load Voltage Limit Change constraints and autonomous reporting parameters
In the table, the number of the first and second,is the minimum value of the limit of the load voltage,the maximum value of the limit voltage of the load,reporting the minimum value of the voltage for the load,reporting the maximum value of the voltage for the load,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
The actual maximum allowable up/down adjustment amount of each load voltage is respectively as follows:
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
In the table, the number of the first and second,the maximum upward adjustment is actually allowed for the load,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
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):
and also
The load polymerization characteristics of the industrial park were:
wherein the content of the first and second substances,
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:
From equations (44) and (46), the voltage control target of each load is:
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:
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:
in summary, the control cost F of the unit electric quantity of the load0(yuan/kwh) is:
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
Power regulation cost model F of industrial parkIP,t-ΔPΣ,tThe following formula:
wherein the content of the first and second substances,
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:
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:
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:
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:
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:
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:
the maximum adjustable power range of the electrolytic aluminum load obtained by the formula (1) is as follows:
step 1.2, the power characteristics of the submerged arc furnace load are as follows:
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:
arc static resistance RarcThe limiting range of (2):
Rarc min≤Rarc≤Rarc max(12)
upper limit of arc partial power factorThere is a lower limit to the arc segment power factorArc partial power factorThe limitations of (2) are:
namely, it is
Wherein the content of the first and second substances,
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:
the power characteristic parameter of the ore furnace during rated operation isAdjusting the power at the momentInto an equivalent low-side voltage variation range, i.e.
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:
step 1.3, the polysilicon load power meets the following electrical relationship:
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:
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 radiationCorresponds 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:
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 α
αmin≤α≤αmax(25)
Wherein, αmin=90%,αmaxWhen the operation is rated, α is 100%;
the following equations (20) and (21) can be obtained:
from equations (23) and (26), the adjustment range of the polysilicon load power can be determined as follows:
effective value U of single-phase voltage of polysilicon loadvalThe adjusting range of (A) is as follows:
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
Step 1.4.1. for electrolytic aluminum load i ═ {1,2, …, NALHas the following components:
PAL,i=m1,i·U*2+m2,iU*(30)
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:
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)
In the formula, PSAF,i、QSAF,iRespectively the load active power and the load reactive power of the ith submerged arc furnace,is the low side voltage rating of the ith furnace,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:
step 1.4.3. load i ═ N to polysiliconAL+NSAF+1,NAL+NSAF+2,…,NAL+NSAF+NPCS}:
PPCS,i=m1,i·U*2(38)
In the formula, PPCS,iThe power is loaded for the ith poly-silicon,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:
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,
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 asVoltage regulation dead zone ofConsidering 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:
the voltage of each load can be adjusted to a practically adjustable minimum/maximum value of
The actual maximum allowable up/down adjustment amount of each load voltage is respectively as follows:
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
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):
and also
The load polymerization characteristics of the industrial park were:
wherein the content of the first and second substances,
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)
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:
From equations (44) and (46), the voltage control target of each load is:
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:
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:
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:
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:
when the voltage of the low-voltage side of the submerged arc furnace isThe electrode control command is then obtained by:
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:
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:
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:
in summary, the control cost F of the unit electric quantity of the load0Unit cell/kwh, is:
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:
wherein Fload,i,tThe control cost for the ith load;
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:
the power adjustment for the ith load on the industrial park is determined by:
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:
wherein the content of the first and second substances,
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:
wherein the content of the first and second substances,
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:
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:
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 independentlyCalculating 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.
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)
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)
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 |
-
2019
- 2019-11-08 CN CN201911088544.7A patent/CN110880772B/en active Active
Patent Citations (4)
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)
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)
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 | |
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 | |
CN112990523B (en) | Hierarchical optimization operation method for regional comprehensive energy system | |
CN110729726B (en) | Intelligent community energy optimization scheduling method and system | |
CN107863773B (en) | Multi-energy-flow energy supply method and system | |
CN110474370B (en) | Cooperative control system and method for air conditioner controllable load and photovoltaic energy storage system | |
CN107769234B (en) | Frequency modulation capability evaluation method of multi-point energy storage coordination distributed power supply | |
CN114169236A (en) | Phase-change heat storage type electric heating negative control system control method based on LSTM algorithm | |
CN116365506A (en) | Energy-saving and loss-reducing optimization method and system for active power distribution network containing distributed photovoltaic | |
CN112200348A (en) | Regional comprehensive energy system multi-target operation decision method considering comprehensive demand response | |
CN112993987A (en) | Active adjustable capacity cooperative electrolytic aluminum load coordination power grid control method | |
CN114825371B (en) | Aggregation temperature control load multi-layer regulation and control method based on node voltage constraints before and after regulation | |
CN102751724A (en) | Prediction-based three-phase load scheduling method and device responding to demand side | |
CN115149540A (en) | Multi-region virtual power plant collaborative optimization operation method based on partition autonomy | |
CN116316553A (en) | Multi-time scale layered operation control method for hydrogen electric coupling system | |
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 | |
CN115879714A (en) | Power system optimization scheduling method for electrolytic aluminum load participating in demand response | |
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 |