CN115276051A - Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics - Google Patents

Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics Download PDF

Info

Publication number
CN115276051A
CN115276051A CN202210870906.3A CN202210870906A CN115276051A CN 115276051 A CN115276051 A CN 115276051A CN 202210870906 A CN202210870906 A CN 202210870906A CN 115276051 A CN115276051 A CN 115276051A
Authority
CN
China
Prior art keywords
power grid
typhoon
energy storage
model
elasticity
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
CN202210870906.3A
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
State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202210870906.3A priority Critical patent/CN115276051A/en
Publication of CN115276051A publication Critical patent/CN115276051A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28—Arrangements for balancing of the load in a network by storage of energy
    • 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/067—Enterprise or organisation modelling
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • H02J3/144—Demand-response operation of the power transmission or distribution 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/381—Dispersed generators
    • 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
    • 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/48—Controlling the sharing of the in-phase component
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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/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
    • 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
    • H02J2310/00—The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/50—The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads
    • H02J2310/56—The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads characterised by the condition upon which the selective controlling is based
    • H02J2310/58—The condition being electrical
    • H02J2310/60—Limiting power consumption in the network or in one section of the network, e.g. load shedding or peak shaving

Abstract

The invention relates to a receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics, which comprises the following steps of: determining an extreme climate disaster model; determining a power transmission channel fault model, and acquiring a vulnerability curve of a power grid element; determining an autonomous shutdown model of the new energy unit under the influence of the extreme climate disasters; determining an energy storage rate discharge model; determining a resistance response phase model and a recovery phase model; acquiring an elasticity evaluation index, and evaluating the elasticity of the power grid according to the size of the evaluation index: the larger the evaluation index value is, the smaller the power grid elasticity is; according to the method, typhoon is taken as an extreme event, on the basis of the existing elastic trapezoid evaluation model, the scene of full generation or autonomous shutdown of the wind turbine generator caused by typhoon crossing is considered, the models of a resistance response stage and a recovery stage are established, the receiving-end power grid elasticity evaluation is carried out, and the power grid elasticity related condition can be accurately obtained.

Description

Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics
Technical Field
The invention belongs to the technical field of power grid evaluation methods, and particularly relates to a receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics.
Background
In an electrical power system, resiliency may be defined as the ability of the grid to resist response and recover as much of the original operating state as possible in extreme events. Along with the increase of the frequency of extreme events, the influence on the power grid is more and more obvious, meanwhile, the new energy grid-connected scale in China is continuously increased, and a large-scale new energy unit is fed into the receiving end power grid, so that the valley-cutting and peak-shaving capacity of the power grid is reduced, and the safe and stable operation of the system is influenced. Although the installed total capacity of the new energy is increased, the phenomena of 'wind abandoning' and 'light abandoning' are serious, and particularly, extreme climates such as ice disasters, high temperature, strong wind and the like have great influence on the new energy machine set. Therefore, it is important to accurately define the grid elasticity and build a suitable evaluation model.
Disclosure of Invention
The technical problem to be solved by the invention is realized by adopting the following technical scheme:
the receiving-end urban power grid elasticity evaluation method considering the response characteristics of new energy and energy storage comprises the following steps:
determining an extreme climate disaster model;
determining a power transmission channel fault model, and acquiring a power grid element vulnerability curve;
determining an autonomous shutdown model of the new energy unit under the influence of the extreme climate disasters;
determining an energy storage rate discharge model;
determining a resistance response phase model and a recovery phase model;
acquiring an elasticity evaluation index, and evaluating the elasticity of the power grid according to the size of the evaluation index: the larger the evaluation index value is, the smaller the grid elasticity is.
Further, the extreme climate disaster is a typhoon disaster, and the typhoon disaster model is as follows:
Figure BDA0003760606850000021
V10＝Kv·Vr
in the formula: t is a unit oftypIs a parameter for adjusting the wind speed distribution; r is the distance from the center of the typhoon; vrIs the wind speed of the typhoon in the area; vmaxIs the maximum wind speed of a typhoon; r ismaxIs the typhoon maximum wind speed radius; kvIs a correction factor;
Figure BDA0003760606850000022
V＝Vd+V10
in the formula: vdIs an axisymmetric wind velocity V10The correction vector of (2); vcIs the moving velocity vector of the typhoon, the correction process is equal to zero at the center of the typhoon, at RmaxTo a maximum of 0.5VcThen decreases radially outward to zero; and V is a calculated wind speed vector.
Further, the grid element vulnerability curve in typhoon can be expressed as:
Figure BDA0003760606850000023
in the formula: v is the maximum wind speed of the line; α is the vulnerability curve coefficient.
Further, the active power output of the wind power station at different wind speeds is as follows:
Figure BDA0003760606850000024
in the formula: pwiThe active power output of the ith wind power station; cwiThe output coefficient of the ith wind power station is obtained; pwmaxiThe maximum output is the maximum output when the ith wind power generation station is full; vin、Vout、VeRespectively the cut-in wind speed, the cut-out wind speed and the rated wind speed of the wind turbine generator.
Further, the energy storage rate charge-discharge model is as follows:
Figure BDA0003760606850000031
Figure BDA0003760606850000032
Figure BDA0003760606850000033
Figure BDA0003760606850000034
in the formula:
Figure BDA0003760606850000035
Pcha max、Pcha mincharging power and upper and lower limits for the energy storage device at the moment t;
Figure BDA0003760606850000036
Prel max、Prel minthe discharge power and the upper and lower limits of the energy storage device at the moment t;
Figure BDA0003760606850000037
Estomaxthe energy and the upper limit of the energy storage device at the time t are stored; etacha、ηrelThe charge-discharge efficiency of the energy storage device; Δ tstoThe charging and discharging time of the energy storage device is long.
Further, a resistance response model of the power grid in the typhoon is established, as shown in the formula:
Figure BDA0003760606850000038
in the formula:
Figure BDA0003760606850000039
the running cost of the unit m at the moment t is calculated;
Figure BDA00037606068500000310
the actual active power output of the unit m at the moment t is obtained; t is the total number of time periods during the typhoon crossing; n is a radical of hydrogenDAnd NGIs the node and generator total;
Figure BDA00037606068500000311
the load shedding amount of the node i at the moment t is shown; μ is a penalty factor.
Further, the recovery phase model is as follows:
Tr＝kwTn
in the formula: t isnThe restoration time of the power grid in a normal state is set; k is a radical ofwIs a weather influence factor, kwTaking a random number of the maximum actual wind speed suffered by the line in the typhoon passing process:
Figure BDA00037606068500000312
in the formula: v is the typhoon real-time velocity; u (x, y) is a random number in the interval (x, y).
Further, the evaluation indexes of the grid elasticity are as follows:
Figure BDA00037606068500000313
in the formula, whereinIIs the ideal curve of the system, lRIs a curve of the actual operation of the system.
The invention has the advantages and positive effects that:
the method takes typhoon as an extreme event, considers the scene of full power generation or autonomous shutdown of the wind turbine generator set caused by typhoon crossing, establishes the models of a resistance response stage and a recovery stage on the basis of the existing elastic trapezoid evaluation model, performs receiving-end power grid elasticity evaluation, and can accurately obtain the related conditions of the power grid elasticity.
Drawings
The technical solutions of the present invention will be described in further detail below with reference to the accompanying drawings and examples, but it should be understood that these drawings are designed for illustrative purposes only and thus do not limit the scope of the present invention. Furthermore, unless otherwise indicated, the drawings are intended to be illustrative of the structural configurations described herein only, and are not necessarily drawn to scale.
FIG. 1 is a schematic diagram of an IEEE14 system for typhoon safety ratio transit according to an embodiment of the present invention;
FIG. 2 is an IEEE14 system elasticity ladder according to an embodiment of the present invention.
Detailed Description
First, it should be noted that the specific structures, features, advantages, etc. of the present invention will be specifically described below by way of example, but all the descriptions are for illustrative purposes only and should not be construed as limiting the invention in any way. Furthermore, any single feature described or implicit in any embodiment or any single feature shown or implicit in any drawing may still be combined or subtracted between any of the features (or equivalents thereof) to obtain still further embodiments of the invention that may not be directly mentioned herein.
It should be noted that, in the present application, the embodiments and features of the embodiments may be combined with each other without conflict.
Example 1
The receiving-end urban power grid elasticity evaluation method considering the response characteristics of new energy and energy storage comprises the following steps:
determining an extreme climate disaster model;
determining a power transmission channel fault model, and acquiring a vulnerability curve of a power grid element;
determining an autonomous shutdown model of the new energy unit under the influence of the extreme climate disasters;
determining an energy storage rate discharge model;
determining a resistance response phase model and a recovery phase model;
and acquiring an elasticity evaluation index, and evaluating the elasticity of the power grid according to the size of the evaluation index. The larger the evaluation index value is, the smaller the elasticity of the power grid is proved under the condition.
Further, the extreme climate disaster is a typhoon disaster, and the typhoon disaster model is as follows:
Figure BDA0003760606850000051
V10＝Kv·Vr
in the formula: t istypThe parameter is a parameter for adjusting the wind speed distribution, and the numerical value is set to be 0.5; r is the distance from the center of the typhoon; vrIs the wind speed of the typhoon in the area; vmaxIs the maximum wind speed of a typhoon; rmaxIs the typhoon maximum wind speed radius; kvThe value is 0.8 for the correction coefficient;
Figure BDA0003760606850000052
V＝Vd+V10
in the formula: vdIs an axisymmetric wind velocity V10The correction vector of (2); vcIs the moving velocity vector of the typhoon, the correction process is equal to zero at the center of the typhoon, at RmaxTo a maximum of 0.5VcThen decreases radially outward to zero; and V is a calculated wind speed vector.
Further, the grid element vulnerability curve in typhoon can be expressed as:
Figure BDA0003760606850000061
in the formula: v is the maximum wind speed of the line, and is generally 30m/s; it can be seen that when V ismaxLess than v, the line will not occurA failure; when V ismaxBetween v and 2v, the fault rate of the line shows exponential multiple increase; when the voltage is more than 2v, the line can be in failure certainly;
wherein α is the vulnerability curve coefficient; compared with the transmission lines in the region, the extra-high voltage transmission lines outside the region are huge and complex, strong in specialization and wide in coverage, and the high reliability of the extra-high voltage transmission lines is the premise of ensuring the safety and stability of a receiving end power grid, so that the transmission lines in the region adopt alpha =2/3, and the transmission lines outside the region adopt alpha =1/3.
Further, the active power output of the wind power station at different wind speeds is as follows:
Figure BDA0003760606850000062
in the formula: pwiThe active power output of the ith wind power station; cwiThe output coefficient of the ith wind power station is obtained; pwmaxiThe maximum output of the ith wind power generation station when the station is full; vin、Vout、VeThe cut-in wind speed, the cut-out wind speed and the rated wind speed of the wind turbine generator are respectively 3 m/s, 14 m/s and 25m/s.
Further, in the above-mentioned case,
the energy storage multiplying power charge-discharge model is as follows:
Figure BDA0003760606850000063
Figure BDA0003760606850000064
Figure BDA0003760606850000065
Figure BDA0003760606850000066
in the formula:
Figure BDA0003760606850000067
Pcha max、Pcha mincharging power and upper and lower limits for the energy storage device at the moment t;
Figure BDA0003760606850000068
Prel max、Prel minthe discharge power and the upper and lower limits of the energy storage device at the moment t;
Figure BDA0003760606850000069
Estomaxthe energy and the upper limit of the energy storage device at the time t are stored; etacha、ηrelThe charge-discharge efficiency of the energy storage device; Δ tstoThe charging and discharging time of the energy storage device is long.
In a further aspect of the present invention,
establishing a resistance response model of the power grid in the typhoon, as shown in the formula:
Figure BDA0003760606850000071
in the formula:
Figure BDA0003760606850000072
the running cost of the unit m at the moment t is calculated;
Figure BDA0003760606850000073
the actual active power output of the unit m at the moment t is obtained; t is the total number of time periods during the typhoon crossing; n is a radical ofDAnd NGIs the node and generator total;
Figure BDA0003760606850000074
the load shedding amount of the node i at the moment t is shown; μ is a penalty factor.
Wherein the constraint conditions are as follows:
Figure BDA0003760606850000075
Figure BDA0003760606850000076
Figure BDA0003760606850000077
Figure BDA0003760606850000078
Figure BDA0003760606850000079
Figure BDA00037606068500000710
Figure BDA00037606068500000711
Figure BDA00037606068500000712
Figure BDA00037606068500000713
Figure BDA00037606068500000714
Figure BDA00037606068500000715
Figure BDA00037606068500000716
in the formula:
Figure BDA00037606068500000717
the actual active power output of the unit m at the moment t is obtained; t is the total number of the time periods during the typhoon crossing; n is a radical of hydrogenDAnd NGIs the node and generator total;
Figure BDA0003760606850000081
the load shedding amount of the node i at the moment t is obtained; μ is a penalty coefficient; cg，CftConnecting a generator connection matrix and a line connection matrix;
Figure BDA0003760606850000082
Figure BDA0003760606850000083
the generator output power matrix, the line tide matrix, the load shedding matrix, the load matrix and the energy storage output matrix at the moment t are respectively; pgmmax，PgmminThe maximum output power and the minimum output power of the unit m are respectively; delta Pgmmax，ΔPgmminRespectively limiting the maximum power change and the minimum power change of the unit m; mGTmThe minimum startup time of the unit m; p isftmaxThe maximum power flow transmitted for the transmission line; omegaG，ΩW，ΩBThe system is a thermal power generator, a wind power generator and a node set.
Further, the recovery phase model is as follows:
Tr＝kwTn
in the formula: t isnThe restoration time of the power grid in a normal state is set; k is a radical ofwAs a weather influence factor, kwTaking a random number of the maximum actual wind speed suffered by the line in the typhoon passing process:
Figure BDA0003760606850000084
in the formula: v is the typhoon real-time velocity; u (x, y) is a random number in the interval (x, y).
Further, the evaluation indexes of the grid elasticity are as follows:
Figure BDA0003760606850000085
in the formula, whereinIIs the ideal curve of the system,/RIs a curve of the actual operation of the system.
Example 2
In the embodiment, an IEEE14 standard node model is improved, and 1-2 generator nodes are selected as the characteristics of a wind turbine generator set for simulating new energy access of a receiving-end power grid; adding 1 power generator node simulation area external incoming call, and connecting a circuit to simulate an extra-high voltage direct current transmission line; and 1-2 generator nodes are selected and set as nodes for installing the energy storage device, so that the output adjustment range of the generator is expanded. The parameters of the IEEE14 system generator node are shown in Table 1.
TABLE 1 IEEE14 System Generator node parameters
Figure BDA0003760606850000091
Typhoon data adopts typhoon safety ratio (Ampil) data in a tropical cyclone data center of the China Meteorological office, wherein a typhoon center coordinate adopts simplified coordinates, the moving distance proportion of typhoon in each time period is guaranteed to be unchanged, the step length of the typhoon data provided by the data is 3-6 h, the simulation step length is 1h, intermediate data is obtained by adopting a linear interpolation method, and the maximum wind speed and the radius of the typhoon in the period are assumed to be unchanged. The coordinates of each node in the system and a simulation graph of a typhoon transit path are shown in fig. 1, a dotted line represents a moving path of a typhoon center, and a circle represents the radius of the maximum wind speed of the typhoon at two moments:
and simulating the expected failure rate value of each line at each moment according to the data analysis of the typhoon passing through, wherein the maximum wind resistance speed of the line is set to be 30m/s. According to the fault rate of each line, simulating the fault condition of an IEEE14 system in a typhoon climate by adopting a Monte Carlo simulation method, obtaining the load change of the system by utilizing a resistance response model and a recovery model, drawing an elastic trapezoidal curve, wherein the curve is the elastic trapezoidal model for evaluating the elasticity of the power grid, and referring to fig. 2, the typhoon influence time, the time spent on the distance of a maintenance team and the repair time, and the repair time of a single element is assumed to be 48h.
The maximum wind-resistant speed of the line is 30m/s, and the restoring time is 48h, the elasticity of the IEEE14 system, namely the loss load in the process of simulating typhoon passing, is 14134MW.
The present invention has been described in detail with reference to the above examples, but the description is only for the preferred examples of the present invention and should not be construed as limiting the scope of the present invention. All equivalent changes and modifications made within the scope of the present invention shall fall within the scope of the present invention.

Claims (8)

1. The receiving-end urban power grid elasticity evaluation method considering the response characteristics of new energy and energy storage is characterized by comprising the following steps of:
determining an extreme climate disaster model;
determining a power transmission channel fault model, and acquiring a power grid element vulnerability curve;
determining an autonomous shutdown model of the new energy unit under the influence of the extreme climate disasters;
determining an energy storage multiplying power discharge model;
determining a resistance response phase model and a recovery phase model;
acquiring an elasticity evaluation index, and evaluating the elasticity of the power grid according to the size of the evaluation index: the larger the evaluation index value is, the smaller the grid elasticity is.
2. The method for assessing the elasticity of the receiving-end urban power grid considering the response characteristics of the new energy and the stored energy according to claim 1, wherein the method comprises the following steps: the extreme climate disaster is a typhoon disaster, and the typhoon disaster model comprises the following steps:
Figure FDA0003760606840000011
V10＝Kv·Vr
in the formula: t is a unit oftypIs a parameter for adjusting the wind speed distribution; r is the distance from the center of the typhoon; vrIs the wind speed of the typhoon in the area; vmaxIs the maximum wind speed of a typhoon; rmaxIs the typhoon maximum wind speed radius; kvIs a correction factor;
Figure FDA0003760606840000012
V＝Vd+V10
in the formula: vdIs an axisymmetric wind velocity V10The correction vector of (2); vcIs the moving velocity vector of the typhoon, the correction process is equal to zero at the center of the typhoon, at RmaxIncreased to maximum 0.5VcThen decreases radially outward to zero; and V is a calculated wind speed vector.
3. The method for assessing the elasticity of the receiving-end urban power grid considering the response characteristics of the new energy and the stored energy according to claim 2, wherein the method comprises the following steps: the grid element vulnerability curve in typhoon can be expressed as:
Figure FDA0003760606840000021
in the formula: v is the maximum wind speed of the line; α is the vulnerability curve coefficient.
4. The method for assessing the flexibility of the receiving-end urban power grid considering the new energy and the energy storage response characteristics according to claim 3, wherein the active power output of the wind power generation station at different wind speeds is as follows:
Figure FDA0003760606840000022
in the formula: pwiActive power output of the ith wind power station is obtained; cwiThe output coefficient of the ith wind power station is the output coefficient of the ith wind power station; pwmaxiThe maximum output of the ith wind power generation station when the station is full; vin、Vout、VeThe cut-in wind speed, the cut-out wind speed and the rated wind speed of the wind turbine generator are respectively.
5. The method for assessing the flexibility of a receiving-end urban power grid considering new energy and energy storage response characteristics according to claim 4,
the energy storage multiplying power charge-discharge model is as follows:
Figure FDA0003760606840000023
Figure FDA0003760606840000024
Figure FDA0003760606840000025
Figure FDA0003760606840000026
in the formula:
Figure FDA0003760606840000027
Pchamax、Pchamincharging power and upper and lower limits for the energy storage device at the moment t;
Figure FDA0003760606840000028
Prelmax、Prelminthe discharge power and the upper and lower limits of the energy storage device at the moment t;
Figure FDA0003760606840000029
Estomaxthe energy and the upper limit of the energy storage device at the time t are stored; etacha、ηrelThe charge-discharge efficiency of the energy storage device is improved; Δ tstoThe charging and discharging time of the energy storage device is long.
6. The method for evaluating the flexibility of the receiving-end urban power grid according to claim 5, wherein the characteristics of response to new energy and stored energy are considered,
establishing a resistance response model of the power grid in the typhoon, as shown in the formula:
Figure FDA0003760606840000031
in the formula:
Figure FDA0003760606840000032
the running cost of the unit m at the moment t is calculated;
Figure FDA0003760606840000033
the actual active power output of the unit m at the moment t is obtained; t is the total number of time periods during the typhoon crossing; n is a radical ofDAnd NGIs the node and generator total;
Figure FDA0003760606840000034
the load shedding amount of the node i at the moment t is shown; μ is a penalty factor.
7. The method for assessing the elasticity of the receiving-end urban power grid considering the new energy and the energy storage response characteristics according to claim 6, wherein the recovery phase model is as follows:
Tr＝kwTn
in the formula: t isnThe restoration time of the power grid in a normal state is set; k is a radical ofwIs a weather influence factor, kwTaking a random number of the maximum actual wind speed suffered by the line in the typhoon passing process:
Figure FDA0003760606840000035
in the formula: v is the typhoon real-time velocity; u (x, y) is a random number in the interval (x, y).
8. The method for evaluating the elasticity of the receiving-end urban power grid considering the new energy and the energy storage response characteristics according to claim 7, wherein the evaluation indexes of the elasticity of the power grid are as follows:
Figure FDA0003760606840000036
in the formula, wherein lIIs the ideal curve of the system,/RIs a curve of the actual operation of the system.
CN202210870906.3A 2022-07-22 2022-07-22 Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics Pending CN115276051A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210870906.3A CN115276051A (en) 2022-07-22 2022-07-22 Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210870906.3A CN115276051A (en) 2022-07-22 2022-07-22 Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics

Publications (1)

Publication Number Publication Date
CN115276051A true CN115276051A (en) 2022-11-01

Family

ID=83769916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210870906.3A Pending CN115276051A (en) 2022-07-22 2022-07-22 Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics

Country Status (1)

Country Link
CN (1) CN115276051A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117688475A (en) * 2024-02-04 2024-03-12 山东电工时代能源科技有限公司 Disaster prediction-based energy network assessment method, system, terminal and storage medium

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117688475A (en) * 2024-02-04 2024-03-12 山东电工时代能源科技有限公司 Disaster prediction-based energy network assessment method, system, terminal and storage medium

Similar Documents

Publication Publication Date Title
US9618914B2 (en) Energy resource-grid-load automatic control system of smart microgrid and control methods thereof
CN107944757A (en) Electric power interacted system regenerative resource digestion capability analysis and assessment method
CN102682222B (en) Continuous tide calculation method based on wind power fluctuation rule
Shourangiz-Haghighi et al. Developing more efficient wind turbines: A survey of control challenges and opportunities
CN104463697A (en) Risk assessment method for power system including large-scale wind power
CN109787297A (en) A kind of extra-high voltage direct-current sending end reactive power optimization method considering transient characterisitics
CN109919398A (en) The zonal reserve Optimal Configuration Method of electric system containing wind-powered electricity generation based on figure partitioning algorithm
CN112287531B (en) Offshore wind farm multi-state reliability obtaining method considering typhoon influence
CN113659620A (en) Water-wind complementary power generation system day-ahead scheduling method based on dynamic frequency constraint
CN115276051A (en) Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics
CN109818347A (en) A kind of appraisal procedure of electric system wind electricity digestion capability
CN109921459B (en) Active frequency control method after direct current blocking of photovoltaic high-ratio transmitting-end power grid
CN108335004A (en) A kind of wind generator system method for evaluating reliability equal based on the electric energy that is obstructed
CN114117730A (en) Elasticity evaluation method for power distribution network under typhoon disaster
Enslin et al. Integrated approach to network stability and wind energy technology for on-shore and offshore applications
CN109995084B (en) Cascade hydropower station-thermal power plant combined optimization scheduling method and system
CN116562563A (en) Power system optimization operation method based on minimum inertia demand evaluation
CN107886227B (en) Method for evaluating disaster resistance improvement degree of power distribution network
CN113346541B (en) Wind power prevention control optimization method under typhoon disaster
CN115882486A (en) Power distribution network comprehensive benefit optimal energy storage configuration method based on double-layer model
CN109787217A (en) Based on the modified spare clearing method out of wind-powered electricity generation multistate model and opportunity cost
CN115313508A (en) Microgrid energy storage optimal configuration method, device and storage medium
Li et al. Reliability evaluation of distribution network considering distributed intermittent renewable energy access
CN109274112B (en) Energy storage power station scheduling method considering regional stability
CN104993526A (en) Multi-type power supply joint operation system assembly ratio calculation method

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