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 PDF

Info

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
Application number
CN202210141478.0A
Other languages
Chinese (zh)
Inventor
项颂
陈璐
吴坚
刘鑫
马继涛
李昱材
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Shenyang Institute of Engineering
State Grid Eastern Inner Mongolia Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Shenyang Institute of Engineering
State Grid Eastern Inner Mongolia Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Shenyang Institute of Engineering, State Grid Eastern Inner Mongolia Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202210141478.0A priority Critical patent/CN114498751A/en
Publication of CN114498751A publication Critical patent/CN114498751A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/388Islanding, i.e. disconnection of local power supply from the network
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/22The renewable source being solar energy
    • H02J2300/24The renewable source being solar energy of photovoltaic origin
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/40Systems 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

Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid
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
Figure BDA0003506569190000021
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 moment
Figure BDA0003506569190000022
Radius 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:
calculating per unit value of wind power output power
Figure BDA0003506569190000023
Figure BDA0003506569190000024
Calculating per unit value of air density
Figure BDA0003506569190000025
Figure BDA0003506569190000026
Calculating the per unit value of wind speed of wind field
Figure BDA0003506569190000027
Figure BDA0003506569190000028
Calculating the per unit value of the thermodynamic temperature of air
Figure BDA0003506569190000029
Figure BDA00035065691900000210
Calculating per unit value of radius of effective area of fan blade
Figure BDA00035065691900000211
Figure BDA00035065691900000212
Calculating per unit value of angular velocity of fan impeller
Figure BDA0003506569190000031
Figure BDA0003506569190000032
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
Figure BDA0003506569190000033
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 moment
Figure BDA0003506569190000034
tlIntensity of solar radiation at time of day
Figure BDA0003506569190000035
Power 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;
calculating the per unit value of the photovoltaic output power
Figure BDA0003506569190000036
Figure BDA0003506569190000037
Calculating the per unit value of the radiation intensity of the sun
Figure BDA0003506569190000038
Figure BDA0003506569190000039
Calculating per unit value of ambient temperature
Figure BDA00035065691900000310
Figure BDA00035065691900000311
Calculating the per unit value of the power temperature coefficient of the photovoltaic panel
Figure BDA00035065691900000312
Figure BDA00035065691900000313
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
Figure BDA0003506569190000041
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:
calculating the per unit value of the wind power output active power
Figure BDA0003506569190000042
Figure BDA0003506569190000043
Calculating wind power output reactive power per unit value
Figure BDA0003506569190000044
Figure BDA0003506569190000045
Calculating per unit value of voltage modulus value
Figure BDA0003506569190000046
Figure BDA0003506569190000047
Calculating per unit value of i and k voltage phase difference of node
Figure BDA0003506569190000048
Figure BDA0003506569190000049
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
Figure BDA00035065691900000410
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 pipe
Figure BDA0003506569190000051
Time thHeat of internal thermodynamic system
Figure BDA0003506569190000052
Coefficient 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:
calculating the per unit value of heat of thermodynamic system
Figure BDA0003506569190000053
Figure BDA0003506569190000054
Calculating the per unit value of mass flow injected into the pipeline
Figure BDA0003506569190000055
Figure BDA0003506569190000056
Calculating the per unit value of the heat transfer coefficient
Figure BDA0003506569190000057
Figure BDA0003506569190000058
In the formula, WhThe heat reference value of the thermodynamic system is obtained;
Figure BDA0003506569190000059
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
Figure BDA00035065691900000510
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 systems
Figure BDA00035065691900000511
Slope correction value of pipeline of natural gas system
Figure BDA00035065691900000512
Time t of system operationgHeight difference H between any two nodes of pipelinegGas pressure P at any two nodesgAnd PkAverage pressure of pipeline
Figure BDA00035065691900000513
Pipe 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:
calculating the per unit value of the height difference of the pipeline
Figure BDA00035065691900000514
Figure BDA00035065691900000515
Calculating the per unit value of the pipeline temperature
Figure BDA00035065691900000516
Figure BDA00035065691900000517
Calculating per unit value of average pressure of pipeline
Figure BDA0003506569190000061
Figure BDA0003506569190000062
Calculating per unit value of energy generated by natural gas system
Figure BDA0003506569190000063
Figure BDA0003506569190000064
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
Figure BDA0003506569190000065
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 moment
Figure BDA0003506569190000071
Radius 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:
calculating per unit value of wind power output power
Figure BDA0003506569190000072
Figure BDA0003506569190000073
Calculating per unit value of air density
Figure BDA0003506569190000074
Figure BDA0003506569190000075
Calculating the per unit value of wind speed of wind field
Figure BDA0003506569190000076
Figure BDA0003506569190000077
Calculating the per unit value of the thermodynamic temperature of air
Figure BDA0003506569190000078
Figure BDA0003506569190000079
Calculating per unit value of radius of effective area of fan blade
Figure BDA00035065691900000710
Figure BDA00035065691900000711
Calculating per unit value of angular velocity of fan impeller
Figure BDA00035065691900000712
Figure BDA00035065691900000713
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
Figure BDA00035065691900000714
Firstly, aiming at a certain region of Mongolian, collecting influence parameters, wherein a certain wind turbine generator set is at twPower generated at any moment
Figure BDA0003506569190000081
Wind 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
Figure BDA0003506569190000082
Calculating the influence factor phi of the wind power output characteristic according to the collected parameterswIs composed of
Figure BDA0003506569190000083
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 moment
Figure BDA0003506569190000084
tlIntensity of solar radiation at time of day
Figure BDA0003506569190000085
Power 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;
calculating the per unit value of the photovoltaic output power
Figure BDA0003506569190000091
Figure BDA0003506569190000092
Calculating the per unit value of the radiation intensity of the sun
Figure BDA0003506569190000093
Figure BDA0003506569190000094
Calculating per unit value of ambient temperature
Figure BDA0003506569190000095
Figure BDA0003506569190000096
Calculating the power temperature coefficient per unit value of the photovoltaic panel
Figure BDA0003506569190000097
Figure BDA0003506569190000098
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
Figure BDA0003506569190000099
Determining the influence factor of photovoltaic power generation, a certain light in the areaThe voltage generator set is at tlOutput power over time
Figure BDA00035065691900000910
tl60s, 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
Figure BDA0003506569190000101
Figure BDA0003506569190000102
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:
calculating the per unit value of the wind power output active power
Figure BDA0003506569190000103
Figure BDA0003506569190000104
Calculating wind power output reactive power per unit value
Figure BDA0003506569190000105
Figure BDA0003506569190000106
Calculating per unit value of voltage modulus value
Figure BDA0003506569190000107
Figure BDA0003506569190000108
Calculating per unit value of i and k voltage phase difference of node
Figure BDA0003506569190000111
Figure BDA0003506569190000112
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
Figure BDA0003506569190000113
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
Figure BDA0003506569190000114
Influence factor of power system energy balance:
Figure BDA0003506569190000115
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 pipe
Figure BDA0003506569190000121
Time thHeat of internal thermodynamic system
Figure BDA0003506569190000122
Coefficient 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:
calculating the per unit value of heat of thermodynamic system
Figure BDA0003506569190000123
Figure BDA0003506569190000124
Calculating the per unit value of mass flow injected into the pipeline
Figure BDA0003506569190000125
Figure BDA0003506569190000126
Calculating the per unit value of the heat conduction coefficient
Figure BDA0003506569190000127
Figure BDA0003506569190000128
In the formula, WhThe heat reference value of the thermodynamic system is obtained;
Figure BDA0003506569190000129
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
Figure BDA00035065691900001210
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 pipe
Figure BDA00035065691900001211
th60s, reference value of mass flow
Figure BDA00035065691900001212
Is a time thHeat of internal thermodynamic system
Figure BDA00035065691900001213
Time 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
Figure BDA0003506569190000131
And (3) obtaining an energy balance influence factor of the thermodynamic system:
Figure BDA0003506569190000132
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 systems
Figure BDA0003506569190000133
Slope correction value of pipeline of natural gas system
Figure BDA0003506569190000134
Time t of system operationgHeight difference H between any two nodes of pipelinegGas pressure P at any two nodesgAnd PkAverage pressure of pipeline
Figure BDA0003506569190000135
Pipe 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:
calculating the per unit value of the height difference of the pipeline
Figure BDA0003506569190000136
Figure BDA0003506569190000137
Calculating per unit value of pipeline temperature
Figure BDA0003506569190000138
Figure BDA00035065691900001313
Calculating per unit value of average pressure of pipeline
Figure BDA0003506569190000139
Figure BDA00035065691900001310
Calculating the per unit value of energy generated by the natural gas system
Figure BDA00035065691900001311
Figure BDA00035065691900001312
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
Figure BDA0003506569190000141
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 time
Figure BDA0003506569190000142
tg=60s,tgReference value of energy generated by natural gas system in time
Figure BDA0003506569190000143
Time 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 pipeline
Figure BDA0003506569190000144
Reference 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
Figure BDA0003506569190000145
Thus, the natural gas system balance influence factor phi can be calculatedg
Figure BDA0003506569190000151
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
Figure BDA0003506569190000152
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:
Figure BDA0003506569190000161
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
Figure FDA0003506569180000011
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 moment
Figure FDA0003506569180000012
Radius 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:
calculating per unit value of wind power output power
Figure FDA0003506569180000021
Figure FDA0003506569180000022
Calculating per unit value of air density
Figure FDA0003506569180000023
Figure FDA0003506569180000024
Calculating the per unit value of wind speed of wind field
Figure FDA0003506569180000025
Figure FDA0003506569180000026
Calculating the per unit value of the thermodynamic temperature of air
Figure FDA0003506569180000027
Figure FDA0003506569180000028
Calculating per unit value of radius of effective area of fan blade
Figure FDA0003506569180000029
Figure FDA00035065691800000210
Calculating per unit value of angular velocity of fan impeller
Figure FDA00035065691800000211
Figure FDA00035065691800000212
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
Figure FDA00035065691800000213
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 moment
Figure FDA00035065691800000214
tlIntensity of solar radiation at time of day
Figure FDA00035065691800000215
Power temperature coefficient alpha of photovoltaic panelTAmbient temperature
Figure FDA00035065691800000216
A 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;
calculating the per unit value of the photovoltaic output power
Figure FDA0003506569180000031
Figure FDA0003506569180000032
Calculating the per unit value of the radiation intensity of the sun
Figure FDA0003506569180000033
Figure FDA0003506569180000034
Calculating per unit value of ambient temperature
Figure FDA0003506569180000035
Figure FDA0003506569180000036
Calculating the power temperature coefficient per unit value of the photovoltaic panel
Figure FDA0003506569180000037
Figure FDA0003506569180000038
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
Figure FDA0003506569180000039
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:
calculating the per unit value of the wind power output active power
Figure FDA00035065691800000310
Figure FDA00035065691800000311
Calculating wind power output reactive power per unit value
Figure FDA00035065691800000312
Figure FDA00035065691800000313
Calculating per unit value of voltage modulus value
Figure FDA0003506569180000041
Figure FDA0003506569180000042
Calculating per unit value of i and k voltage phase difference of node
Figure FDA0003506569180000043
Figure FDA0003506569180000044
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
Figure FDA0003506569180000045
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 pipe
Figure FDA0003506569180000046
Time thHeat of internal thermodynamic system
Figure FDA0003506569180000047
Coefficient 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:
calculating the per unit value of heat of thermodynamic system
Figure FDA0003506569180000048
Figure FDA0003506569180000049
Calculating the per unit value of mass flow injected into the pipeline
Figure FDA00035065691800000410
Figure FDA00035065691800000411
Calculating the per unit value of the heat conduction coefficient
Figure FDA00035065691800000412
Figure FDA00035065691800000413
In the formula, WhThe heat reference value of the thermodynamic system is obtained;
Figure FDA00035065691800000414
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
Figure FDA0003506569180000051
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 systems
Figure FDA0003506569180000052
Slope correction value of pipeline of natural gas system
Figure FDA0003506569180000053
Time t of system operationgHeight difference H between any two nodes of pipelinegGas pressure P at any two nodesgAnd PkAverage pressure of pipeline
Figure FDA0003506569180000054
Pipe 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:
calculating the per unit value of the height difference of the pipeline
Figure FDA0003506569180000055
Figure FDA0003506569180000056
Calculating per unit value of pipeline temperature
Figure FDA0003506569180000057
Figure FDA0003506569180000058
Calculating the per unit value of the average pressure of the pipeline
Figure FDA0003506569180000059
Figure FDA00035065691800000510
Calculating per unit value of energy generated by natural gas system
Figure FDA00035065691800000511
Figure FDA00035065691800000512
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
Figure FDA0003506569180000061
CN202210141478.0A 2022-02-16 2022-02-16 Renewable energy power generation-based method for determining multi-energy stability of outgoing power grid Pending CN114498751A (en)

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)

* Cited by examiner, † Cited by third party
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

Cited By (2)

* Cited by examiner, † Cited by third party
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