CN114498751A - Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid - Google Patents
Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid Download PDFInfo
- Publication number
- CN114498751A CN114498751A CN202210141478.0A CN202210141478A CN114498751A CN 114498751 A CN114498751 A CN 114498751A CN 202210141478 A CN202210141478 A CN 202210141478A CN 114498751 A CN114498751 A CN 114498751A
- Authority
- CN
- China
- Prior art keywords
- calculating
- power
- per unit
- energy
- phi
- 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.)
- Pending
Links
- 238000010248 power generation Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 29
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 78
- 239000003345 natural gas Substances 0.000 claims abstract description 39
- 239000007789 gas Substances 0.000 claims description 11
- 230000005855 radiation Effects 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 8
- 239000013256 coordination polymer Substances 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 7
- 238000004146 energy storage Methods 0.000 description 6
- 230000033228 biological regulation Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000007599 discharging Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/388—Islanding, i.e. disconnection of local power supply from the network
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Power Engineering (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a method for determining multi-energy stability of an outgoing power grid based on renewable energy power generation, which relates to the field of power systems based on renewable energy power generation, and comprises the steps of firstly analyzing wind-light output characteristics, obtaining influence factors influencing wind-light output according to the wind-light output characteristics, then calculating energy balance influence factors of a power system, a thermodynamic system and a natural gas system of a sending end power grid, finally integrating the influence factors of wind power generation and photovoltaic power generation and the energy balance influence factors of the power system, the thermodynamic system and the natural gas system of the sending end power grid, calculating stability indexes of the sending end system, and judging the stability of the sending end system; the stability of the system is judged according to the value of the stability index, the stability of the system can be conveniently and effectively determined, and the effect of adjusting the balance influence factor by adjusting energy to achieve the stability of the system is achieved.
Description
Technical Field
The invention relates to the field of power systems based on renewable energy power generation, in particular to a method for determining multi-energy stability of an outgoing power grid based on renewable energy power generation.
Background
With the development of economic society, the permeability of renewable energy sources is higher and higher, and for areas with rich power generation of renewable energy sources, if the local load demand cannot completely consume the electric energy generated by a high proportion of renewable energy sources, the electric energy needs to be sent out, so that the electric energy loss is reduced. And the output of renewable energy sources such as wind, light and the like has uncertainty, indirectness and instability, and various energy flows of electricity, heat and gas have more and more complex influence factors, so that the multi-energy of a power grid at a transmitting end is unstable.
Much research is being conducted at home and abroad on energy flow calculation. An extended Newton-Lavenson electric/heat/gas multi-energy flow calculation method is provided at present, a Jacobian matrix reflecting a multi-energy coupling relation is derived by establishing a mathematical model of an electric power network, a thermal power network and a natural gas network and considering the operation mode of an electric power system in multi-energy flow, a multi-energy flow calculation process is provided, and the method has effectiveness and practicability under two scenes of grid connection and island of the electric power system; an energy flow model of a combined cooling heating and power system is established; the energy flow calculation method is used as a coupling node of a power distribution network and a natural gas network, a finite element node method suitable for calculating the energy flow of the natural gas network is deduced, and the energy flow calculation method of the electro-pneumatic coupling micro energy network is provided, so that the rapid solving of the multi-time discontinuous steady-state energy flow is realized; however, for the outgoing power grid, the output of renewable energy has uncertainty, indirectness and instability, and the influence factors of the stability of the energy at the sending end are complex, so that the method cannot accurately determine the stability of the outgoing power grid.
Disclosure of Invention
In order to solve the problems, the stability of various energy flows of the power grid at the sending end is calculated based on renewable energy power generation, and when the various energy flows of the power grid at the sending end are unstable, the stability of the power grid at the sending end is adjusted by utilizing charging and discharging of stored energy.
In order to achieve the technical effects, the invention provides a method for determining the multi-energy stability of an outgoing power grid based on renewable energy power generation, which comprises the following steps:
step 1: analyzing the wind-solar output characteristics to obtain influence factors influencing the wind-solar output;
step 2: calculating energy balance influence factors of a power system, a thermodynamic system and a natural gas system of a sending end power grid;
and step 3: and calculating the stability index of the sending end system and judging the stability of the sending end system.
The step 1 comprises the following steps:
step 1.1: calculating wind power output characteristic influence factor phiw;
Step 1.2: calculating a photovoltaic output impact factor phil。
The step 2 comprises the following steps:
step 2.1: calculating an influence factor phi of energy balance of a power system of a transmission-end power gride;
Step 2.2: calculating an influence factor phi of energy balance of a power grid thermodynamic system at a sending endh;
Step 2.3: calculating an influence factor phi of energy balance of a natural gas system of a transmission-end power gridg。
The step 3 comprises the following steps:
step 3.1: calculating stability index phi of an outgoing end systemz:
Step 3.2: according to the stability index phizRegulating and controlling the energy balance of a sending end system: when the stability index satisfies delta1≤|φz|≤δ2When the system is in an unstable state, the system at the sending end is indicated; when the stability index is | phiz|>δ1Or | phiz|<δ2Time, indicates that the sending end system is in an unstable state, where, delta1、δ2To set the threshold.
The step 1.1 comprises the following steps:
step 1.1.1: collecting an influence parameter phi influencing wind power generation in a wind power generator set1,κ(ii) a The influence parameter Φ1,κComprising a wind turbine at twPower generated at any momentRadius R of effective area of fan blade and angular velocity omega of impellerlLocal air density rho, wind field wind speed v and wind energy utilization coefficient CPThermodynamic temperature T of airw;
Step 1.1.2: calculating an influencing parameter phi1,κPer unit value of the relevant parameters:
In the formula, pwFor the wind power output power reference value, rhowIs a reference value of air density, vwIs a wind field wind speed reference value, TwIs a thermodynamic temperature reference value of air, RwThe radius reference value is the effective area of the fan blade;
step 1.1.3: calculating the wind power output characteristic influence factor phiw:
The step 1.2 comprises the following steps:
step 1.2.1: collecting influence parameter phi of photovoltaic power generation2,κ(ii) a The influence parameter Φ2,κIncluding photovoltaic standard conditions at tlPower generated at any momenttlIntensity of solar radiation at time of dayPower temperature coefficient alpha of photovoltaic panelTAmbient temperature TtlA parameter epsilon related to the atmospheric quality, an altitude angle alpha of the sun in the area, a latitude phi of the area, a declination angle delta of the sun and a time angle omega of the sun;
step 1.2.2: calculating an influencing parameter phi2,κPer unit value of the relevant parameter;
In the formula, plIs a photovoltaic output power reference value; i islThe radiation intensity reference value of the sun; t islIs ambient temperatureA reference value; alpha is alphaBTThe power temperature coefficient reference value of the photovoltaic is obtained;
step 1.2.3: calculating a photovoltaic output impact factor phil:
The step 2.1 comprises the following steps:
step 2.1.1: determining an impact factor omega affecting energy balance of an electric power system1,κ(ii) a The influence factor omega1,κInjected active power P comprising a node iiAnd reactive power Qi(ii) a Voltage modulus V of node iiAnd the phase difference delta between the voltages of nodes i and kik(ii) a Grid synchronization coefficient lambdagActive control coefficient lambdapReactive power control coefficient lambdaq(ii) a Node number n of power system and wind power output characteristic influence factor phiw(ii) a Photovoltaic output impact factor phil;
Step 2.1.2: calculating the influence factor omega1,κPer unit value of the relevant parameters:
In the formula, peThe reference value of active power is obtained; qe,iIs a reactive power reference value; ve,iThe reference value is the voltage modulus; deltae,ikThe reference value is the voltage phase difference of the nodes i and k;
step 2.1.3: calculating an impact factor phi of an energy balance of an electric power systeme:
The step 2.2 comprises:
step 2.2.1: determining an influencing factor omega influencing thermal power transmission in a thermodynamic system2,κ(ii) a The influence factor omega2,κIncluding the specific heat capacity c of waterpTime thMass flow injected through internal pipeTime thHeat of internal thermodynamic systemCoefficient of thermal conductivity lambdahTemperature T of the beginning and end of the pipelines,hAnd Te,dLength l of the heat pipe;
step 2.2.2: calculating the influence factor omega2,κPer unit value of the relevant parameters:
In the formula, WhThe heat reference value of the thermodynamic system is obtained;a mass flow standard value for the pipeline injection; lambda [ alpha ]HThe reference value of the heat conduction coefficient is obtained;
step 2.2.3: calculating thermodynamic system energy balance influence factor phih:
The step 2.3 comprises:
step 2.3.1: determining an impact factor omega affecting the energy balance of a natural gas system3,κ(ii) a The influence factor omega3,κIncluding the time tgEnergy produced by internal natural gas systemsSlope correction value of pipeline of natural gas systemTime t of system operationgHeight difference H between any two nodes of pipelinegGas pressure P at any two nodesgAnd PkAverage pressure of pipelinePipe constant C between any two nodesgkTemperature T of pipelinegAverage gas flow temperature TaCoefficient of friction between pipes f, relative density of pipes SG;
Step 2.3.2: calculating the influence factor omega3,κPer unit value of the relevant parameters:
In the formula, Delta HGA reference value for the height of the pipeline; t isGThe reference value is the temperature of the pipeline; pAThe reference value is the average pressure of the pipeline; wG,tAs a natural gas systemSystem tgThe energy reference value generated in time;
step 2.3.3: calculating an influence factor phi of the energy balance of the natural gas systemg:
The invention has the beneficial effects that:
the invention provides a method for determining the multi-energy stability of an outgoing power grid based on renewable energy power generation, aiming at the problem of instability of fluctuation of power, heat and gas energy of a transmission-end power grid, firstly calculating influence factors of wind and light output of renewable energy, then calculating influence factors influencing energy balance, obtaining stability indexes of the outgoing system according to the influence factors, and judging the stability of the system according to values of the stability indexes.
Drawings
Fig. 1 is a flowchart of a method for determining multi-energy stability of an outgoing power grid based on renewable energy power generation according to the present invention.
Detailed Description
The invention is further described with reference to the following figures and specific examples. Because the renewable energy power generation refers to wind power and photovoltaic power generation, the influence factors of wind and light output are analyzed, and the influence factors of the wind power generation and the photovoltaic power generation are calculated according to the characteristics of the wind and light output; then analyzing the influence factors of the electricity, gas and heat of the power grid at the transmitting end, and respectively calculating the influence factors of the energy flows; and finally, integrating the influence factors of wind power generation and photovoltaic power generation and the energy balance influence factors of a power system, a thermodynamic system and a natural gas system of a sending end power grid to obtain a stability index of an external sending end system, judging the stability of the system according to the result of the index, and if the system is unstable, charging and discharging are carried out by using an energy storage system to adjust the balance influence factors so as to achieve the stability of the system.
As shown in fig. 1, a method for determining multi-energy stability of an outgoing power grid based on renewable energy power generation includes:
step 1: analyzing the wind-solar output characteristics to obtain influence factors influencing the wind-solar output; the method comprises the following steps:
step 1.1: calculating wind power output characteristic influence factors; the method comprises the following steps:
step 1.1.1: because wind and solar power generation uncertainty is large, in order to maintain the stability of the system, influence factors influencing wind power generation need to be calculated firstly, and influence parameters phi influencing wind power generation in a wind power generator set are collected1,κ(ii) a The influence parameter phi1,κComprising a wind turbine at twPower generated at any momentRadius R (rad/s) of effective area of fan blade and angular speed omega of impellerlLocal air density ρ (kg/m)3) Wind field wind speed v (m/s) and wind energy utilization coefficient CP(relating manufacturing parameters of the fan to local wind parameters, etc.), thermodynamic temperature T of the airw;
Step 1.1.2: calculating an influencing parameter phi1,κPer unit value of the relevant parameters:
In the formula, pwFor the wind power output power reference value, rhowIs a reference value of air density, vwIs a wind field wind speed reference value, TwIs a thermodynamic temperature reference value of air, RwThe radius reference value is the effective area of the fan blade;
step 1.1.3: calculating the wind power output characteristic influence factor phiw:
Firstly, aiming at a certain region of Mongolian, collecting influence parameters, wherein a certain wind turbine generator set is at twPower generated at any momentWind power twTime of dayThe reference value of the output power is Pw30000 kw; the local air density rho is 1.29kg/m3The reference value of the air density is 1.205kg/m3(ii) a Wind field wind speed is 10m/s, and the reference value of the wind field wind speed is 12 m/s; the thermodynamic temperature of local air is 293.15 ℃, and the thermodynamic temperature reference value of air is 298.15 ℃; the radius of the effective area of the fan blade is 52m, and the radius reference value of the effective area of the fan blade is 55 m; the angular speed of the fan impeller is 3405r/min, and the reference value of the angular speed of the fan impeller is 3500 r/min; efficiency of wind energy utilization CP=53.6%。
Calculating the per unit value of the wind power output characteristic influence parameter, as shown in table 2:
TABLE 2 wind power output characteristic influence parameter per unit value
Calculating the influence factor phi of the wind power output characteristic according to the collected parameterswIs composed of
Step 1.2: calculating a photovoltaic output influence factor; the method comprises the following steps:
step 1.2.1: determining influence factors of photovoltaic power generation, and collecting influence parameters phi of photovoltaic power generation2,κ(ii) a The influence parameter phi2,κIncluding photovoltaic standard conditions at tlPower generated at any momenttlIntensity of solar radiation at time of dayPower temperature coefficient alpha of photovoltaic panelTAmbient temperature TtlA parameter epsilon related to the atmospheric quality, an altitude angle alpha of the sun in the area, a latitude phi of the area, a declination angle delta of the sun and a time angle omega of the sun;
step 1.2.2: calculating an influencing parameter phi2,κPer unit value of the relevant parameter;
In the formula, plIs a photovoltaic output power reference value; i islThe radiation intensity reference value of the sun; t islIs an ambient temperature reference value; alpha is alphaBTThe power temperature coefficient reference value of the photovoltaic is obtained;
step 1.2.3: calculating a photovoltaic output impact factor phil:
Determining the influence factor of photovoltaic power generation, a certain light in the areaThe voltage generator set is at tlOutput power over timetl60s, photovoltaic generator set tlReference value p of the instantaneous output powerl30000kw, actual solar radiation intensity I of a photovoltaic generator setl,t=957W/m2Reference value I of intensity of solar radiationl=1000W/m2(ii) a Power temperature coefficient alpha of photovoltaic panelT-0.408%/DEG C, power temperature coefficient reference value alpha of photovoltaic panelBT-0.500%/° c; ambient temperature Ttl20 ℃, reference value T of ambient temperaturelThe altitude angle α of the sun in the area is 32 °, the latitude of the area is 38 °, the declination angle δ of the sun is 19 °, the time angle ω of the sun is 36 °, and the per unit values are as shown in table 3:
TABLE 3 photovoltaic output characteristic influence parameter per unit value
Step 2: calculating energy balance influence factors of a power system, a thermodynamic system and a natural gas system of a sending end power grid; the method comprises the following steps:
in order to maintain stable energy flow of the power grid at the sending end, influence factors of the power, heat and gas energy flow of the power grid at the sending end need to be analyzed to obtain influence factors of energy balance of a power system, a thermodynamic system and a natural gas system of the power grid at the sending end.
Step 2.1: calculating an influence factor of energy balance of a power system of a transmission-end power grid; the method comprises the following steps:
step 2.1.1: in a delivery system, in order to ensure energy balance between a transmitting-end power grid and a receiving-end power grid, an influence factor omega influencing energy balance of a power system is firstly determined1,κ(ii) a The influence factor omega1,κThe injection of node i is first determined as followsActive power PiAnd reactive power Qi(ii) a Then determining the voltage modulus V of the node iiAnd the phase difference delta between the voltages of nodes i and kik(ii) a Grid synchronization coefficient lambdagActive control coefficient lambdapReactive power control coefficient lambdaq(ii) a Node number n of power system and wind power output characteristic influence factor phiw(ii) a Photovoltaic output impact factor phil;
Step 2.1.2: calculating the influence factor omega1,κPer unit value of the relevant parameters:
In the formula, peThe reference value of active power is obtained; qe,iIs a reactive power reference value; ve,iThe reference value is the voltage modulus; deltae,ikIs a voltage of node i and kA reference value of the potential difference;
step 2.1.3: calculating an impact factor phi of an energy balance of an electric power systeme:
In the outgoing system, in order to ensure energy balance between a transmitting-end power grid and a receiving-end power grid, an influence factor influencing power balance is determined firstly. Firstly determining the injection active power P of the node ii22453kw, the reference value of the injected active power of the node i is 30000kw, and the injected reactive power Q of the node ii12253var, reference value Q of injected reactive power of node ii10000 var; then determining the voltage modulus V of the node ii218V, reference value V of voltage modulus of node ii220V and the phase difference delta between the voltages of the nodes i and kik23.7 DEG, reference value delta for the phase difference of the voltagee,ik30 °; grid synchronization coefficient lambdag1.034, active control coefficient λp0.896, reactive power control coefficient lambdaq0.824; the number of nodes n in the power system is 30, phil=1.067,φwThe influence factor of the energy balance of the power system can be obtained as 1.042.
The per unit values of the influence parameters for calculating the energy balance of the power system are shown in table 4:
TABLE 4 influence parameters per unit value of energy balance of electric power system
Influence factor of power system energy balance:
step 2.2: calculating an influence factor of energy balance of a power grid thermodynamic system at a sending end; the method comprises the following steps:
step 2.2.1: due to heat in the thermodynamic systemThe temperature of the force transmission medium is greatly influenced by the environment, and an influence factor omega influencing thermal transmission in a thermal system is determined2,κ(ii) a The influence factor omega2,κIncluding the specific heat capacity c of waterpTime thMass flow injected internally through a pipeTime thHeat of internal thermodynamic systemCoefficient of thermal conductivity lambdahTemperature T of the beginning and end of the pipelines,hAnd Te,dLength l of the heat pipe;
step 2.2.2: calculating the influence factor omega2,κPer unit value of the relevant parameters:
In the formula, WhThe heat reference value of the thermodynamic system is obtained;mass flow for pipe injectionA standard value; lambda [ alpha ]HThe reference value of the heat conduction coefficient is obtained;
step 2.2.3: calculating thermodynamic system energy balance influence factor phih:
Since the temperature of the heat transfer medium in the thermodynamic system is greatly influenced by the environment, the following parameters, namely the specific heat capacity c of water, need to be determinedp=4.2×103J/(kg. DEG C.) and time thMass flow injected internally through a pipeth60s, reference value of mass flowIs a time thHeat of internal thermodynamic systemTime thHeat reference value W of internal thermodynamic systemh350000J; coefficient of thermal conductivity lambdah45W/m.K, heat conduction coefficient reference value lambdaH50W/m.K, pipeline starting point Ts,h25 DEG end point temperature Te,dThe heat pipe length l is 10m at 135 °, per unit value set as shown in table 5:
TABLE 5 influence parameter per unit value of thermodynamic system energy balance
And (3) obtaining an energy balance influence factor of the thermodynamic system:
step 2.3: calculating an influence factor of energy balance of a natural gas system of a transmission-end power grid; the method comprises the following steps:
step 2.3.1: in the outward transmission system, in order to ensure the energy balance between the transmission end power grid and the receiving end power grid, finally determining an influence factor influencing the natural gas balance, and determining an influence factor omega influencing the energy balance of the natural gas system3,κ(ii) a The influence factor omega3,κIncluding the time tgEnergy produced by internal natural gas systemsSlope correction value of pipeline of natural gas systemTime t of system operationgHeight difference H between any two nodes of pipelinegGas pressure P at any two nodesgAnd PkAverage pressure of pipelinePipe constant C between any two nodesgkTemperature T of pipelinegAverage gas flow temperature TaCoefficient of friction between pipes f, relative density of pipes SG;
Step 2.3.2: calculating the influence factor omega3,κPer unit value of the relevant parameters:
In the formula, Delta HGA reference value for the height of the pipeline; t isGThe reference value is the temperature of the pipeline; pAThe reference value is the average pressure of the pipeline; wG,tFor natural gas system tgThe energy reference value generated in time;
step 2.3.3: calculating an impact factor phi of the energy balance of a natural gas systemg:
In the outward transmission system, in order to ensure the energy balance between the transmission end power grid and the receiving end power grid and finally determine the influence factors influencing the natural gas balance, the influence factors influencing the natural gas system energy balance are firstly determined, including tgEnergy produced by natural gas system over timetg=60s,tgReference value of energy generated by natural gas system in timeTime t of system operationg(second/S unit), height difference Δ H of two pointsg8.6m, the reference value delta H of the height difference of the pipeline between two pointsG10 m; mean pressure of pipelineReference value P of average pressure of pipelineA3.0MPa, pipeline constant C between nodesgk=1.086×10-5Temperature T of pipeg125 deg. reference value T of pipe temperatureG135 °; coefficient of friction between pipes f is 0.234, pipe relative density SG=1.293kg/m3Reference value S of relative density of pipelineG=1.293kg/m3(ii) a The per unit value settings are shown in table 6:
TABLE 6 sending end electric network natural gas system energy balance influence parameter per unit value
Thus, the natural gas system balance influence factor phi can be calculatedg:
Synthesizing the influence factors of wind power generation and photovoltaic power generation and the energy balance influence factors of a power system, a thermodynamic system and a natural gas system of a sending end power grid to obtain a stability index phi of an outgoing end systemzAnd according to the calculation result of the index, as shown in table 1, judging the stability of the system, and if the system is unstable, utilizing the energy storage system to charge and discharge so as to adjust the stability index of the system to achieve the stability of the system.
TABLE 1 stability index φ for delivery end systemz
φzValue range | System stability determination | System regulation |
0.25≤|φz|≤0.63 | System stabilization | Without energy storage regulation |
|φzLess than or equal to 0.25 or phiz|≥0.63 | Instability of the system | Energy storage charge-discharge regulation |
And step 3: calculating the stability index of the sending end system and judging the stability of the sending end system; the method comprises the following steps:
step 3.1: calculating stability index phi of an outgoing end systemz:
Step 3.2: according to the stability index phizRegulating and controlling the energy balance of a sending end system: when the stability index is more than or equal to 0.25 ≦ phizWhen the | is less than or equal to 0.63, the sending end system is in a stable state; when the stability index is | phiz| less than 0.25 or | φz|>When 0.63, the sending end system is in an unstable state; when phiz|>0.63, the stability index of the external sending end system is larger than the stable value, and the external sending end system needs to be charged through an energy storage device, so that the value of an influence factor is reduced; when phi iszAnd | < 0.25, the stability index of the sending end system is less than a stable value, the discharging is carried out through an energy storage device, and the value of the stability index of the sending end system is increased to maintain the balance of the sending end system.
Stability index phi of delivery end systemzThe calculation is as follows:
due to phiz≈0.341,0.25≤|φzAnd if the | is less than or equal to 0.63, obtaining that the system of the sending end outside the region is stable.
Claims (9)
1. A method for determining multi-energy stability of an outgoing power grid based on renewable energy power generation is characterized by comprising the following steps:
step 1: analyzing the wind-solar output characteristics to obtain influence factors influencing the wind-solar output;
step 2: calculating energy balance influence factors of a power system, a thermodynamic system and a natural gas system of a sending end power grid;
and step 3: and calculating the stability index of the sending end system and judging the stability of the sending end system.
2. The method for determining multi-energy stability of renewable energy power generation based delivery grid according to claim 1, wherein the step 1 comprises:
step 1.1: calculating wind power output characteristic influence factor phiw;
Step 1.2: calculating a photovoltaic output impact factor phil。
3. The method for determining multi-energy stability of renewable energy power generation based delivery grid according to claim 1, wherein the step 2 comprises:
step 2.1: calculating an influence factor phi of energy balance of a power system of a transmission-end power gride;
Step 2.2: calculating an influence factor phi of energy balance of a power grid thermodynamic system at a sending endh;
Step 2.3: calculating an influence factor phi of energy balance of a natural gas system of a transmission-end power gridg。
4. The method for determining multi-energy stability of renewable energy power generation based delivery grid according to claim 1, wherein the step 3 comprises:
step 3.1: computing distribution systemOverall stability index phiz:
Step 3.2: according to the stability index phizRegulating and controlling the energy balance of a sending end system: when the stability index satisfies delta1≤|φz|≤δ2When the system is in an unstable state, the system at the sending end is indicated; when the stability index is | phiz|>δ1Or | phiz|<δ2Time, indicates that the sending end system is in an unstable state, where, delta1、δ2To set the threshold.
5. The method for determining multi-energy stability of an outgoing power grid based on renewable energy power generation as claimed in claim 2, wherein said step 1.1 comprises:
step 1.1.1: collecting an influence parameter phi influencing wind power generation in a wind power generator set1,κ(ii) a The influence parameter phi1,κComprising a wind turbine at twPower generated at any momentRadius R of effective area of fan blade and angular velocity omega of impellerlLocal air density rho, wind field wind speed v and wind energy utilization coefficient CPThermodynamic temperature T of airw;
Step 1.1.2: calculating an influencing parameter phi1,κPer unit value of the relevant parameters:
In the formula, pwFor a wind power output power reference value, ρwIs a reference value of air density, vwIs a wind field wind speed reference value, TwIs a thermodynamic temperature reference value of air, RwThe radius reference value is the effective area of the fan blade;
step 1.1.3: calculating wind force outputCharacteristic influence factor phiw:
6. The method for determining multi-energy stability of renewable energy power generation based delivery grid according to claim 2, wherein the step 1.2 comprises:
step 1.2.1: collecting influence parameter phi of photovoltaic power generation2,κ(ii) a The influence parameter phi2,κIncluding photovoltaic standard conditions at tlPower generated at any momenttlIntensity of solar radiation at time of dayPower temperature coefficient alpha of photovoltaic panelTAmbient temperatureA parameter epsilon related to the atmospheric quality, the altitude angle alpha of the sun in the area, the latitude phi of the area, the declination angle delta of the sun and the time angle omega of the sun;
step 1.2.2: calculating an influencing parameter phi2,κPer unit value of the relevant parameter;
In the formula, plIs a photovoltaic output power reference value; i islThe radiation intensity reference value of the sun; t islIs an ambient temperature reference value; alpha is alphaBTThe power temperature coefficient reference value of the photovoltaic is obtained;
step 1.2.3: calculating a photovoltaic output impact factor phil:
7. The method according to claim 3, wherein the step 2.1 comprises:
step 2.1.1: determining an impact factor omega affecting energy balance of an electric power system1,κ(ii) a The influence factor omega1,κInjected active power P comprising a node iiAnd reactive power Qi(ii) a Voltage modulus V of node iiAnd the phase difference delta between the voltages of nodes i and kik(ii) a Grid synchronization coefficient lambdagActive control coefficient lambdapReactive power control coefficient lambdaq(ii) a Node number n of power system and wind power output characteristic influence factor phiw(ii) a Photovoltaic output impact factor phil;
Step 2.1.2: calculating the influence factor omega1,κPer unit value of the relevant parameters:
In the formula, peThe reference value of active power is obtained; qe,iIs a reactive power reference value; ve,iThe reference value is the voltage modulus; deltae,ikThe reference value is the voltage phase difference of the nodes i and k;
step 2.1.3: calculating an impact factor phi of an energy balance of an electric power systeme:
8. The method according to claim 3, wherein the step 2.2 comprises:
step 2.2.1: determining an influencing factor omega influencing thermal power transmission in a thermodynamic system2,κ(ii) a The influence factor omega2,κIncluding the specific heat capacity c of waterpTime thMass flow injected internally through a pipeTime thHeat of internal thermodynamic systemCoefficient of thermal conductivity lambdahTemperature T of the beginning and end of the pipelines,hAnd Te,dLength l of the heat pipe;
step 2.2.2: calculating the influence factor omega2,κPer unit value of the relevant parameters:
In the formula, WhThe heat reference value of the thermodynamic system is obtained;a mass flow standard value for the pipeline injection; lambda [ alpha ]HThe reference value of the heat conduction coefficient is obtained;
step 2.2.3: calculating thermodynamic system energy balance influence factor phih:
9. The method according to claim 3, wherein the step 2.3 comprises:
step 2.3.1: determining an impact factor omega affecting the energy balance of a natural gas system3,κ(ii) a The influence factor omega3,κIncluding the time tgEnergy produced by internal natural gas systemsSlope correction value of pipeline of natural gas systemTime t of system operationgHeight difference H between any two nodes of pipelinegGas pressure P at any two nodesgAnd PkAverage pressure of pipelinePipe between any two nodesConstant CgkTemperature T of pipelinegAverage gas flow temperature TaCoefficient of friction between pipes f, relative density of pipes SG;
Step 2.3.2: calculating the influence factor omega3,κPer unit value of the relevant parameters:
In the formula, Delta HGA reference value for the height of the pipeline; t isGThe reference value is the temperature of the pipeline; pAThe reference value is the average pressure of the pipeline; wG,tFor natural gas system tgThe energy reference value generated in time;
step 2.3.3: calculating natural gas system energyInfluence factor phi of the quantity balanceg:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210141478.0A CN114498751A (en) | 2022-02-16 | 2022-02-16 | Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210141478.0A CN114498751A (en) | 2022-02-16 | 2022-02-16 | Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114498751A true CN114498751A (en) | 2022-05-13 |
Family
ID=81480725
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210141478.0A Pending CN114498751A (en) | 2022-02-16 | 2022-02-16 | Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114498751A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116937576A (en) * | 2023-09-15 | 2023-10-24 | 国网山西省电力公司运城供电公司 | Large-scale active power distribution network stability evaluation method, system and equipment |
-
2022
- 2022-02-16 CN CN202210141478.0A patent/CN114498751A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116937576A (en) * | 2023-09-15 | 2023-10-24 | 国网山西省电力公司运城供电公司 | Large-scale active power distribution network stability evaluation method, system and equipment |
CN116937576B (en) * | 2023-09-15 | 2023-11-24 | 国网山西省电力公司运城供电公司 | Large-scale active power distribution network stability evaluation method, system and equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107808216B (en) | Method for constructing comprehensive minimum optimization model for wind curtailment and electrical thermal load reduction of electricity-gas-heat interconnection system | |
WO2019205561A1 (en) | Cchp micro-grid structure including compressed air energy storage and operation method therefor | |
CN105656064B (en) | Realize that thermoelectricity unit participates in the method that peak-load regulating is dispatched using heat supply time lag | |
CN106593771B (en) | Electric heating decoupled system and its control method in a kind of co-generation unit | |
CN105162167A (en) | Adaptive-droop-control-based wind-photovoltaic-energy-storage micro-grid frequency modulation method | |
CN108988356B (en) | Electric heating microgrid tie line power fluctuation stabilizing method based on virtual energy storage | |
CN111371104B (en) | Power grid frequency stability control method based on wind-storage combined power generation system | |
CN106953367B (en) | Wind power consumption method for thermoelectric system containing gas-steam combined cycle unit | |
CN113255198B (en) | Multi-objective optimization method for combined cooling heating and power supply micro-grid with virtual energy storage | |
CN114498751A (en) | Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid | |
Lin et al. | Analysis of capacity and control strategy for distributed energy system with hybrid energy storage system | |
CN113159983B (en) | Coordinated control method for off-grid phosgene combined heat and power generation comprehensive energy system | |
CN108964039A (en) | A kind of electric-gas coupled system steady stability of wind power plant containing DFIG domain tangent plane method for solving | |
CN111193261A (en) | Day-ahead optimization method of multi-energy flow system based on building equivalent heat energy storage | |
CN114362173B (en) | Digital twinning-based power grid load flow stability evaluation method for transmitting end | |
CN112528210B (en) | Combined cooling heating and power system and control method | |
CN113446656B (en) | Power-load matched photovoltaic photo-thermal PV/T combined cooling heating and power system regulation and control method | |
CN113761727B (en) | Method for constructing optimal scheduling model of combined heat and power system containing distributed electric heat pump | |
CN114218765A (en) | Comprehensive energy inertial support method for coping with power shortage of power grid | |
CN113595105A (en) | Method and system for realizing power grid frequency modulation by utilizing adiabatic compressed air energy storage | |
Chen et al. | Air source heat pump energy storage heating system for smart building | |
CN111030140A (en) | Battery energy storage power grid frequency stability control method based on big data | |
CN113483479B (en) | Auxiliary service method and system combining variable frequency air conditioner and energy storage battery | |
Lin et al. | MPPT and Energy Storage Control of a Hybrid PV Array and Heat Pump System: An Experimental Research | |
CN114221364B (en) | Photovoltaic-energy storage-load coordination operation method for low-voltage distribution network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |