CN111444593B - Method for improving vulnerability of elements of electricity-gas comprehensive energy system - Google Patents

Method for improving vulnerability of elements of electricity-gas comprehensive energy system Download PDF

Info

Publication number
CN111444593B
CN111444593B CN202010135812.2A CN202010135812A CN111444593B CN 111444593 B CN111444593 B CN 111444593B CN 202010135812 A CN202010135812 A CN 202010135812A CN 111444593 B CN111444593 B CN 111444593B
Authority
CN
China
Prior art keywords
natural gas
power
gas pipeline
node
fault
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.)
Active
Application number
CN202010135812.2A
Other languages
Chinese (zh)
Other versions
CN111444593A (en
Inventor
丁一
桑茂盛
包铭磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202010135812.2A priority Critical patent/CN111444593B/en
Publication of CN111444593A publication Critical patent/CN111444593A/en
Application granted granted Critical
Publication of CN111444593B publication Critical patent/CN111444593B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for improving the vulnerability of elements of an electric-gas comprehensive energy system. Inputting original data, establishing an operation state model of the electricity-gas integrated energy system, and obtaining an operation result in an initial normal state; extracting a network topological structure of the comprehensive energy system, and processing to obtain structural vulnerability parameters of the elements; simulating element faults to obtain an operation result of the system in a fault state; then obtaining the airflow distribution entropy of the natural gas system and the power flow distribution entropy of the power system to obtain the operation vulnerability parameters of the elements; and calculating the weight of each parameter by using an optimal combination weight method to obtain the comprehensive vulnerability of the element. The invention can quickly perform key protection on elements with higher vulnerability, improves the reliability of the elements by thickening the power line and reinforcing the natural gas pipeline, effectively improves the overall operation reliability and robustness of the system, and reduces the probability of large-area faults of the system.

Description

Method for improving vulnerability of elements of electric-gas integrated energy system
Technical Field
The invention belongs to the technical field of power systems and comprehensive energy systems, and relates to a method for improving the vulnerability of elements of an electricity-gas comprehensive energy system under the condition of deep coupling of the power system and a natural gas system.
Background
In recent years, with the gradual increase of the specific gravity of natural gas power generation and the application of gas source equipment depending on power supply of a power grid, the coupling operation of a natural gas system and a power system is increasingly tight. Obviously, the interdependence and interplay of the natural gas system and the electrical power system can increase the vulnerability of the electric-gas integrated energy system. On one hand, the stoppage of a gas transmission natural gas pipeline in a natural gas system may cause the interruption of the supply of natural gas fuel of a gas turbine unit, further cause the insufficient power generation capacity of a power system and cause power failure accidents; on the other hand, the power line is broken, so that the working power supply of the gas source equipment which depends on the power supply of the power grid is lost, the gas source equipment cannot normally work, and the gas supply capacity of the natural gas system is insufficient. Therefore, it is necessary to provide a method for improving vulnerability of components of an electricity-gas integrated energy system by comprehensively considering the coupling operation characteristics of an electric power system and a natural gas system.
At present, a vulnerability improving method, especially a power system vulnerability improving method, mainly aims at improving the vulnerability of the whole system, for example, researching the whole vulnerability of the power system based on a complex network theory or a power system operation analysis theory, and then providing a measure for improving the vulnerability of the system. However, these methods do not analyze the components in the system, resulting in poor vulnerability improvement of the system components. In addition, the current research method only considers the scene of independent operation of the power system, does not consider the influence of other energy networks closely coupled with the power system, enables the two power systems and the natural gas system to be closely connected through the coupled operation of the two power systems and the natural gas system, and ensures that the vulnerability of elements depends on not only the self system but also other coupled systems.
Therefore, it is necessary to comprehensively consider the coupling operation characteristics of the power system and the natural gas system, and to provide a method for improving the vulnerability of elements of the electric-gas integrated energy system, find out the vulnerable elements, and perform important protection on the elements to improve the operation reliability and robustness of the whole system.
Disclosure of Invention
In order to solve the problems of the background art, the invention aims to provide a method for improving the vulnerability of elements of an electric-gas integrated energy system. The invention can quickly find the elements with higher vulnerability, improve the reliability of the elements by thickening the power line and reinforcing the key protection of the natural gas pipeline, effectively improve the overall operation reliability and robustness of the system and reduce the probability of large-area faults of the system.
In order to achieve the purpose, the invention adopts the specific technical scheme that the method comprises the following steps:
step 1, inputting original data, establishing an operation state model of the electricity-gas integrated energy system, and calculating to obtain an operation result of the system in an initial normal state;
the original data refers to the power, load and other data of the electricity-gas comprehensive energy system in the current operation state.
Step 2, extracting a network topological structure of the electricity-gas integrated energy system according to an operation result in an initial normal state, and processing to obtain structural vulnerability parameters of elements;
step 3, simulating element faults including power line outage and natural gas pipeline outage, and obtaining an operation result of the electricity-gas integrated energy system in a fault state by using an electricity-gas integrated energy system operation state model;
step 4, processing and obtaining a natural gas system airflow distribution entropy and a power system power flow distribution entropy based on an operation result in a fault state, so as to obtain an operation vulnerability parameter of the element;
step 5, calculating the weight of each parameter by using an optimal combination weight method based on each vulnerability parameter obtained in the step 2 and the step 4, thereby obtaining the comprehensive vulnerability of the element; and selecting elements with higher comprehensive weakness for processing, and improving the reliability of the elements by thickening power lines and reinforcing natural gas pipelines.
The electricity-gas comprehensive energy system consists of an electric power system and a natural gas system. The natural gas node is present in a natural gas system and refers to gas source equipment, compressor equipment and load equipment in the natural gas system, wherein the gas source equipment comprises an independent power supply gas source and a power supply gas source depending on a power grid. The power grid nodes exist in a power system and refer to power stations, substations and load equipment in the power system, wherein the power stations comprise two types of conventional units and gas units. The natural gas nodes are connected through natural gas pipelines, the power grid nodes are connected through power lines, and the two end nodes of each natural gas pipeline and each power line are respectively called as a head end node and a tail end node. The natural gas system and the electric power system are coupled through a gas unit and a gas source depending on a power grid for power supply, natural gas fuel consumed by the gas unit is sourced from gas source equipment of the natural gas system, and a working power source depending on the gas source supplied by the power grid is sourced from a power station of the electric power system.
The step 1 is specifically as follows:
step 1.1, establishing a minimum Minf of the sum of the operating energy consumption of the natural gas system and the power system as an objective function of an operating state model of the electricity-gas comprehensive energy system:
Figure GDA0002522606730000021
in the formula phig、ΦeRespectively a natural gas node set and a power grid node set;
Figure GDA0002522606730000022
respectively indicating the gas supply amount and the gas load removal amount of a gas source of a natural gas node i in a state l;
Figure GDA0002522606730000023
the energy consumption for the operation of the natural gas node i in the state l is generally
Figure GDA0002522606730000024
And
Figure GDA0002522606730000025
a quadratic or linear function of;
Figure GDA0002522606730000026
and
Figure GDA0002522606730000027
the active output power of the conventional unit and the active output power of the gas unit of the power grid node j in the state l are respectivelyAnd electrical load power shedding amount;
Figure GDA0002522606730000028
the energy consumption of the grid node j in the state l is generally
Figure GDA0002522606730000029
And
Figure GDA00025226067300000210
a quadratic or linear function of; i. j is the ordinal number of the natural gas node and the power grid node respectively; l is an ordinal number of the system running state, and the initial normal running state is represented when l is 0;
step 1.2, establishing equivalent constraint of the following electric-gas integrated energy system operation state model:
Figure GDA0002522606730000031
Figure GDA0002522606730000032
Figure GDA0002522606730000033
Figure GDA0002522606730000034
Figure GDA0002522606730000035
Figure GDA0002522606730000036
in the formula, Fload,i、
Figure GDA0002522606730000037
Respectively the gas load of a natural gas node i and the gas consumption of the gas unit in a state l; (i, t) denotes a natural gas pipeline connecting the head end node i and the tail end node t; phipipeRepresents a collection of natural gas pipelines;
Figure GDA0002522606730000038
representing the amount of gas flowing through the natural gas pipeline (i, t) in state l; m is a group ofitRepresenting a transmission parameter of the natural gas pipeline (i, t);
Figure GDA0002522606730000039
respectively representing the pressure square values of natural gas nodes i and t in the state l; pload,jRepresenting the active power load of the grid node j; (j, s) represents a power line connecting the head end node j and the tail end node s; philineRepresenting a set of power lines;
Figure GDA00025226067300000310
represents the active power flowing through the power line (j, s) in state l; x is the number ofjsRepresents the reactance of the power line (j, s);
Figure GDA00025226067300000311
and
Figure GDA00025226067300000312
respectively representing voltage phase angles of the grid nodes j and s in the state l; gHV、ηG2PRespectively representing the heat value of natural gas and the power generation efficiency of the gas engine set; j belongs to i and represents that the natural gas fuel of the gas turbine set on the power grid node j is provided by a natural gas node i; p issour,jThe active power which is consumed by the gas source connected to the grid node j to maintain normal work is represented; etaP2GRepresenting the conversion coefficient of electricity to gas; i belongs to j and represents that the working power supply of the gas source on the natural gas node i is provided by the power grid node j;
step 1.3, establishing an inequality constraint of an electric-gas integrated energy system operation state model:
Figure GDA00025226067300000313
Figure GDA00025226067300000314
Figure GDA00025226067300000315
Figure GDA00025226067300000316
Figure GDA00025226067300000317
Figure GDA00025226067300000318
Figure GDA00025226067300000319
Figure GDA00025226067300000320
Figure GDA00025226067300000321
Figure GDA00025226067300000322
Figure GDA0002522606730000041
in the formula (I), the compound is shown in the specification,
Figure GDA0002522606730000042
the working state of an air source on a natural gas node i under the state l is shown, when the normal working value of the air source is 1, otherwise, the working state is 0; (8) - (9) when the electrical load shedding of the grid node j is greater than a certain threshold, i.e.
Figure GDA0002522606730000043
Then
Figure GDA0002522606730000044
Indicating that the supplied gas source is shut down; the superscripts of the letters carry "max" and "min" respectively, which represent the upper and lower limits of the variable;
and finally, inputting the original data into the model, and solving by adopting a linear programming method to obtain an operation result of the system in an initial normal state.
The operation result comprises the air supply quantity of an air source of a natural gas node i in a state l
Figure GDA0002522606730000045
Gas load removal amount of natural gas node i in state l
Figure GDA0002522606730000046
Active power output power of conventional unit of power grid node j in state l
Figure GDA0002522606730000047
Gas turbine unit active output power of power grid node j in state l
Figure GDA0002522606730000048
And the electric load power cut-off quantity of the grid node j under the state l
Figure GDA0002522606730000049
The state l is a normal state or a fault state. The conventional units are referred to as coal-fired units.
The step 2 is specifically as follows: the structural vulnerability comprises a capability number centrality parameter and a capability number centrality parameter of the natural gas pipeline and an electrical number centrality parameter of the power line.
Step 2.1, according to the upper limit of the output of the gas source at the natural gas node i
Figure GDA00025226067300000410
Gas consumption of gas unit under initial normal operation state at natural gas node t
Figure GDA00025226067300000411
Gas load F at natural gas node tload,tProcessing and obtaining the capacity betweenness centrality parameter CBC of the natural gas pipeline according to the following formulap：
Figure GDA00025226067300000412
Figure GDA00025226067300000413
In the formula (I), the compound is shown in the specification,
Figure GDA00025226067300000414
representing a set of shortest paths between natural gas nodes i and t;
Figure GDA00025226067300000415
the weight of the mth shortest path between the nodes i and t is represented, and the calculation formula is shown in (20);
Figure GDA00025226067300000416
representing the gas consumption of the gas unit of the natural gas node i in an initial normal operation state;
Figure GDA00025226067300000417
the transmission capacity of the mth shortest path between the nodes i and t is represented and is the minimum value of the transmission capacities of all the natural gas pipelines contained in the shortest path;
Figure GDA00025226067300000418
a decision variable indicating whether the natural gas pipeline p is in the shortest path, when the natural gas pipeline p is in the mth shortest path between the nodes i and t,
Figure GDA00025226067300000419
otherwise
Figure GDA00025226067300000420
m is the ordinal number of the shortest path between the natural gas nodes i and t; p is the ordinal number of the natural gas pipeline;
step 2.2, according to the upper limit of the output of the gas source at the natural gas node i
Figure GDA00025226067300000421
Gas consumption of gas unit in initial normal operation state at natural gas node i
Figure GDA00025226067300000422
Gas load F at natural gas node iload,iProcessing and obtaining the capability degree centrality parameter CDC of the natural gas pipeline according to the following formulap：
Figure GDA00025226067300000423
CDCp＝(dcap,p1+dcap,p2)/2(22)
In the formula, c(i,t)Representing the transmission capacity of the natural gas pipeline (i, t); the t epsilon i represents that the natural gas node t is connected with the node i; p1 and p2 are respectively a head end node and a tail end node of the natural gas pipeline p;
step 2.3, according to the upper limit of the output of the conventional unit at the j position of the power grid node
Figure GDA0002522606730000051
Upper limit of gas turbine set output
Figure GDA0002522606730000052
At grid node s there isWork power load Pload,sThe active power P consumed by the gas source connected with the power grid node s to maintain normal worksour,sThe following formula is adopted to process and obtain the electrical betweenness centrality parameter EBC of the power linek：
Figure GDA0002522606730000053
Figure GDA0002522606730000054
In the formula (I), the compound is shown in the specification,
Figure GDA0002522606730000055
representing a set of shortest paths between grid nodes j, s, one path possibly comprising a plurality of power lines;
Figure GDA0002522606730000056
the weight of the nth shortest path between the nodes j and s is represented, and a calculation formula is shown in (24);
Figure GDA0002522606730000057
the transmission capacity of the nth shortest path between the grid nodes j and s is represented and is the minimum value of the transmission capacities of all the power lines contained in the shortest path;
Figure GDA0002522606730000058
as a decision variable whether or not power line k is in the shortest path, when line k is in the nth shortest path between nodes j, s,
Figure GDA0002522606730000059
otherwise
Figure GDA00025226067300000510
n is the ordinal number of the shortest path between the grid nodes j and s; k is the ordinal number of the power line;
step 2.4, according to the output of the conventional unit at the j position of the power grid nodeLimit of
Figure GDA00025226067300000511
Upper limit of gas turbine set output
Figure GDA00025226067300000512
Active power load P at grid node jload,jAnd the active power P consumed by the gas source connected to the power grid node j to maintain normal worksour,jProcessing and obtaining the electrical degree centrality parameter EDC of the power line according to the following formulak：
Figure GDA00025226067300000513
EDCk＝(dcapa,k1+dcapa,k2)/2(26)
In the formula, c(j,s)Representing the transmission capacity of the power line (j, s); s belongs to j and represents that the power grid node s is connected with the node j; k1 and k2 are respectively the head end node and the tail end node of the power line k.
The step 3 is specifically as follows: running state model of electricity-gas comprehensive energy system
Step 3.1, according to the simulated natural gas pipeline outage fault, updating the network topology structure of the electricity-gas integrated energy system in the simulated natural gas pipeline outage fault process, inputting the updated network topology structure of the electricity-gas integrated energy system and original data into the operation state model of the electricity-gas integrated energy system to calculate and solve to obtain operation data of the system in the state, wherein the operation data comprises gas flow rate in the natural gas pipeline p fault state
Figure GDA00025226067300000514
And active power flow of natural gas pipeline p in fault state
Figure GDA00025226067300000515
And further calculating the airflow increment of the natural gas pipeline and the active power flow increment of the power line under the outage fault of the natural gas pipeline according to the following formulas according to the operation data:
Figure GDA00025226067300000516
Figure GDA0002522606730000061
in the formula (I), the compound is shown in the specification,
Figure GDA0002522606730000062
respectively representing the airflow increment of the natural gas pipeline q and the active power flow increment of the power line b caused by the shutdown fault of the natural gas pipeline p;
Figure GDA0002522606730000063
and
Figure GDA0002522606730000064
respectively representing the transmission capacity of the natural gas pipeline q, the gas flow in the initial normal state and the gas flow in the failure state of the natural gas pipeline p;
Figure GDA0002522606730000065
and
Figure GDA0002522606730000066
respectively representing the transmission capacity of the power line b, the active power flow in the initial normal state and the active power flow in the fault state of the natural gas pipeline p;
step 3.2, simulating the outage fault of the power line, updating the network topology structure of the electric-gas integrated energy system in the process of simulating the outage fault of the power line, inputting the updated network topology structure of the electric-gas integrated energy system and original data into the operation state model of the electric-gas integrated energy system, calculating and solving to obtain the operation data of the system in the state, wherein the operation data comprises the gas flow rate of the natural gas pipeline q in the k fault state of the power line
Figure GDA0002522606730000067
And power line b in power line k fault conditionActive power flow of
Figure GDA0002522606730000068
And further calculating the airflow increment of the natural gas pipeline and the active power flow increment of the power line under the outage fault of the power line according to the following formulas according to the operation data:
Figure GDA0002522606730000069
Figure GDA00025226067300000610
in the formula (I), the compound is shown in the specification,
Figure GDA00025226067300000611
respectively representing the airflow increment of the natural gas pipeline q and the active power flow increment of the power line b caused by the disconnection fault of the power line k;
Figure GDA00025226067300000612
representing the gas flow of the natural gas pipeline q in the fault state of the power line k;
Figure GDA00025226067300000613
representing the active power flow of the power line b in a fault condition of the power line k.
The natural gas pipeline is shut down, which means that the natural gas pipeline is attacked by people or disconnected due to natural disasters. The power line outage is caused by tripping, disconnection and the like of the power line due to human attack or natural disaster.
The step 4 is specifically as follows:
step 4.1, processing according to the airflow increment of the natural gas pipeline under the outage fault of the natural gas pipeline and the power line and the following formula to obtain the vulnerability parameter PFI of the outage fault of the natural gas pipelinep：
Figure GDA00025226067300000614
In the formula (I), the compound is shown in the specification,
Figure GDA00025226067300000615
indicating the gas flow rate of the natural gas pipeline p in the initial normal state;
Figure GDA00025226067300000616
respectively representing natural gas system airflow distribution entropy and power system load flow distribution entropy caused by natural gas pipeline p outage fault; u. of1、u2Respectively representing the influence weight of the natural gas pipeline outage fault on a natural gas system and an electric power system;
the airflow distribution entropy and the load flow distribution entropy caused by the natural gas pipeline p fault are calculated as follows:
Figure GDA00025226067300000617
Figure GDA0002522606730000071
Figure GDA0002522606730000072
Figure GDA0002522606730000073
in the formula etaqp、ηbpRespectively representing the airflow impact ratio caused by the shutdown fault of the natural gas pipeline p to the natural gas pipeline q and the power flow impact ratio caused by the shutdown fault of the natural gas pipeline p to the line b;
step 4.2, processing according to the active power flow increment of the power line under the outage faults of the natural gas pipeline and the power line and the following formula to obtain the fault vulnerability parameter LFI of the power linek：
Figure GDA0002522606730000074
In the formula (I), the compound is shown in the specification,
Figure GDA0002522606730000075
representing the active power flow power of the power line k in an initial normal state;
Figure GDA0002522606730000076
respectively representing the power flow distribution entropy of the power system and the airflow distribution entropy of the natural gas system caused by the k fault of the line; v. of1、v2Respectively representing the influence weight of the line fault on the power system and the natural gas system;
the power flow distribution entropy and the airflow distribution entropy caused by the fault of the line k are calculated as follows:
Figure GDA0002522606730000077
Figure GDA0002522606730000078
Figure GDA0002522606730000079
Figure GDA00025226067300000710
in the formula etabk、ηqkAnd respectively representing the power flow impact ratio caused by the disconnection fault of the line k to the line b and the airflow impact ratio caused to the natural gas pipeline q.
The step 5 is specifically as follows:
step 5.1, establishing the following optimal combination weight model, and calculating the optimal weight of each vulnerability parameter:
Figure GDA00025226067300000711
ρs,β＝τs,β/(τs,β+τo,β)(42)
ρo,β＝τo,β/(τs,β+τo,β)(43)
Figure GDA00025226067300000712
wβ≥0(45)
in the formula, wβAn optimal weight representing the parameter β; tau iss,β、τo,βRespectively representing subjective weight obtained by calculating the parameter beta according to an analytic hierarchy process and objective weight obtained by calculating according to an entropy weight method; rhos,β、ρo,βThe bias weight coefficients respectively representing the parameter beta, the subjective weight and the objective weight are obtained by calculation through (42) and (43); m is the total number of all vulnerability parameters; beta is the ordinal number of the vulnerability parameter;
and 5.2, integrating the structural vulnerability parameters and the operation vulnerability parameters obtained in the step 2 and the step 4, and obtaining the integrated vulnerability of each natural gas pipeline and each power line by adopting the following formula:
VG(p)＝ωg1CBCp+ωg2CDCp+ωg3PFIp(46)
VE(k)＝ωe1EBCk+ωe2EDCk+ωe3LFIk(47)
in the formula, ωg1、ωg2And ωg3Are respectively CBCp、CDCpAnd PFIpThe parameters are calculated according to the step 5.1 to obtain the optimal weight; omegae1、ωe2And ωe3Are respectively EBCk、EDCkAnd LFIkThe parameters are calculated according to the step 5.1 to obtain the optimal weight;
step 5.3, calculating comprehensive fragility values of elements of the natural gas pipelines/power lines according to the step 5.2, improving the fragility of the electricity-gas comprehensive energy system according to the comprehensive fragility values, selecting the elements with high comprehensive fragility for key processing, and improving the reliability of the elements by thickening the power lines and reinforcing the natural gas pipelines, wherein the method specifically comprises the following steps of:
for elements with an integrated vulnerability above a preset vulnerability threshold,
if the element is a power line, measures of thickening the power line, increasing the split number of the line and replacing the line with a double-circuit line are taken to improve the transmission capacity of the power line and reduce the vulnerability;
if the element is a natural gas pipeline, the reliability of the natural gas pipeline is improved and the vulnerability is reduced by adopting measures of thickening the wall of the natural gas pipeline and reducing the gas transmission pressure of the pipeline.
And for the elements with the comprehensive vulnerability not higher than the preset vulnerability threshold value, not processing.
In the step 2.1: the natural gas pipeline capacity betweenness centrality parameter CBCpThe capacity and distribution condition of gas source equipment and load equipment, the coupling condition of a gas unit and a natural gas node, the transmission capacity of a natural gas pipeline and the like are considered, the role of the natural gas pipeline in the network transmission function is revealed from the perspective of the global network, the larger the centrality of the natural gas pipeline capacity betweenness is, the more the number of times of the shortest paths passing through the natural gas pipeline is, and the larger the transmission capacity of each shortest path is, so that the larger the transmission function borne by the natural gas pipeline in the global structure is, and the more the key is in the network.
In the step 2.2: the natural gas pipeline capability degree centrality parameter CDCpThe capacity and distribution of gas source equipment and load equipment, the coupling condition of a gas unit and a natural gas node, the transmission capacity of a natural gas pipeline and the like are considered, the role of the natural gas pipeline in the network transmission function is revealed from the local view of a network, the larger the centrality of the natural gas pipeline capacity is, the more the number of the natural gas pipelines connected with the natural gas pipeline is, and the larger the capacity of the connected natural gas pipeline is, so that the larger the transmission function of the natural gas pipeline in a local structure is, the more the key is in the network.
In said step 2.3: the electric power circuit gas dielectric constant centrality parameter EBCkConsidering the capacities and distribution conditions of power generation equipment and load equipment, the coupling condition of air source equipment and a power grid node, the transmission capacity of a power line and the like, the role of the power line in the network transmission function is revealed from the perspective of the global network, the larger the centrality of the electrical medium number of the power line is, the more the number of times of the shortest path passing through the line is, and the larger the transmission capacity of each shortest path is, so that the larger the transmission role of the line in the global structure is, and the more critical the transmission role of the line in the network is.
In the step 2.4: the electric power line electrical degree centrality parameter EDCkConsidering the capacities and distribution of power generation equipment and load equipment, the coupling condition of air source equipment and a power grid node, the transmission capacity of a power line and the like, the role of the power line in the network transmission function is revealed from the local view of the network, the larger the electrical degree centrality of the power line is, the larger the number of lines connected with the line is, and the larger the capacity of the connected line is, so that the larger the transmission role of the line in a local structure is, and the more critical the transmission role of the line in the network is.
In the step 4.1: the vulnerability parameter PFI of the natural gas pipeline outage faultpThe disturbance degree of the natural gas pipeline fault on the natural gas system and the disturbance degree of the natural gas pipeline fault on the power system are comprehensively reflected, the greater the vulnerability of the natural gas pipeline outage fault is, the greater the natural gas airflow impact and line tide impact caused by the natural gas pipeline fault are, and the more uneven the impact distribution is, the more easily the full load of other natural gas pipelines and lines is caused, so that the cascade fault is caused, and therefore, the more fragile the natural gas pipeline is.
In the step 4.2: the power line disconnection fault vulnerability parameter LFIkThe disturbance degree of a line fault on the power system and the disturbance degree of a natural gas system are comprehensively reflected, the greater the vulnerability of the line disconnection fault is, the greater the tidal current impact and the greater the airflow impact caused by the line fault are, and the more uneven the impact distribution is, the more easily the full load of other lines and natural gas pipelines is caused to cause cascade faults, so the more fragile the line is.
In the step 5.2: the comprehensive weakness V of the natural gas pipelineG(p) and the combined vulnerability V of the power lineE(k) The vulnerability of natural gas pipelines and power lines in the electric-gas integrated energy system is revealed in two aspects of the integrated network topology structure and the system operation state. The greater the comprehensive vulnerability of the element is, the greater the transmission role the element bears in the global and local structures of the network is, the greater the disturbance impact of the element fault on the self system and the coupled system is, and the higher the vulnerability of the element is.
The invention has the following beneficial effects:
the method of the invention is a method for improving the vulnerability of elements of the electricity-gas comprehensive energy system under the condition of deep coupling of the power system and the natural gas system, and overcomes the defect that the traditional vulnerability improvement method only analyzes and improves the whole vulnerability of a single energy system; the method can evaluate the vulnerability of the elements from two angles of network topology and system running state, more accurately find out the vulnerable elements in the system, and can effectively improve the running reliability and robustness of the whole system by improving the vulnerability of the vulnerable elements.
The invention can meet the engineering application requirements under the condition of deep coupling of a power system and a natural gas system in the future, quickly improve the vulnerability of each fragile element in the electricity-gas comprehensive energy system, perform key protection on the fragile element, improve the reliability of the element by thickening a power line and reinforcing a natural gas pipeline, and effectively reduce the probability of large-area faults of the electricity-gas comprehensive energy system.
Drawings
FIG. 1 is a schematic structural view of an embodiment of an electric-gas integrated energy system of the present invention;
FIG. 2 is a flow chart of the method of the present invention.
Detailed Description
The invention is further illustrated by the following figures and examples.
The specific embodiment of the complete method according to the invention is as follows:
the invention is described in detail by taking an electricity-gas comprehensive energy system composed of an IEEE30 node power system and a Belgian 20 node natural gas system as an example and combining a technical scheme and a drawing.
As shown in fig. 1, the IEEE30 node power system includes 7 generator sets and 41 power lines, and the belgium 20 node natural gas system includes 6 gas source devices and 19 natural gas pipelines. The generator sets on the power grid nodes 5, 7, 8, 11 and 13 are gas generator sets, natural gas fuel is supplied by the natural gas nodes 3, 12, 6, 10 and 15 respectively, and the rest generator sets are conventional generator sets. The gas sources on the natural gas nodes 1, 2, 5 and 8 are power supply sources depending on a power grid, the power supply sources are respectively supplied by the power grid nodes 14, 21, 30 and 26 to maintain normal operation, and the rest gas sources are independent gas sources. In addition, for the convenience of calculation, the belgium 20 node system is modified as follows: combining the two-loop natural gas pipelines into a single loop; taking 1.3 times of natural gas pipeline transmission airflow in the initial operation state as the natural gas pipeline capacity, wherein the natural gas pipeline capacity is less than 5 multiplied by 103m3·h-1By 5X 103m3·h-1And (4) calculating. The IEEE30 node system is modified as follows: simplifying all lines into non-directional weighted edges and not counting parallel capacitor branches; the line capacity is 1.5 times of the transmission power of the line in the initial operation state, and the less than 15MW is calculated by 15 MW.
The implementation flow of the invention is shown in fig. 2, and the specific steps are as follows:
step 1, inputting original data, establishing an operation state model of the electricity-gas integrated energy system, and calculating to obtain an operation result of the system in an initial normal state;
step 1.1, establishing a minimum Minf of the sum of the operating energy consumption of the natural gas system and the power system as an objective function of an operating state model of the electricity-gas comprehensive energy system:
step 1.2, establishing equivalent constraint of the following electric-gas integrated energy system operation state model:
step 1.3, establishing an inequality constraint of an electric-gas integrated energy system operation state model:
the initial normal operating state of the electricity-gas comprehensive energy system obtained by the solution in the step 1 is shown in table 1.
TABLE 1 initial Normal running State of the Electricity-gas Integrated energy System
Figure GDA0002522606730000101
Figure GDA0002522606730000111
Step 2, extracting a network topological structure of the system based on an operation result in an initial normal state, and processing to obtain a structural vulnerability parameter of the element;
step 2.1, processing and obtaining the capability betweenness centrality parameter CBC of the natural gas pipelinep：
Step 2.2, processing and obtaining capability degree centrality parameter CDC of the natural gas pipelinep：
Step 2.3, processing and obtaining the electrical betweenness centrality parameter EBC of the power linek：
Step 2.4, processing and obtaining electrical degree centrality parameter EDC of power linek：
Simulating element faults including power line tripping and natural gas pipeline outage, and calculating the operation result of the system in the fault state;
step 3.1, simulating the outage fault of the natural gas pipeline, and calculating the airflow increment of the natural gas pipeline and the active power flow increment of the power line:
step 3.2, simulating the power line disconnection fault (or the power line tripping caused by the misoperation of the relay protection device), and calculating the airflow increment of the natural gas pipeline and the active power flow increment of the power line:
step 4, processing and obtaining a natural gas system airflow distribution entropy and a power system power flow distribution entropy based on an operation result in a fault state, so as to obtain an operation vulnerability parameter of the element;
step 4.1, processing and obtaining the vulnerability parameter PFI of the outage fault of the natural gas pipelinep：
Step 4.2,Processing and obtaining a fault vulnerability parameter LFI of the power line according to the following formulak：
And 5, calculating the weight of each parameter by using an optimal combination weight method based on each vulnerability parameter obtained in the step 2 and the step 4, thereby obtaining the comprehensive vulnerability of the element.
Step 5.1, establishing the following optimal combination weight model, and calculating the optimal weight of each vulnerability parameter:
and 5.2, integrating the structural vulnerability parameters and the operation vulnerability parameters obtained in the step 2 and the step 4 to obtain the integrated vulnerability of the natural gas pipeline and the power line:
and 5.3, calculating comprehensive weakness values of the natural gas pipelines and the lines according to the step 5.2, selecting elements with high comprehensive weakness for key protection, and improving the reliability of the elements by thickening the power lines and reinforcing the natural gas pipelines.
And (3) calculating the comprehensive fragility of each element in the electricity-gas comprehensive energy system, sorting the elements from large to small, and screening the natural gas pipelines and the power lines which are ranked in the top five as the fragile elements, wherein the comprehensive fragility values of the fragile elements are shown in a table 2.
TABLE 2 comprehensive vulnerability value of vulnerable component of electro-pneumatic comprehensive energy system
Figure GDA0002522606730000121
It can be seen that the screened vulnerable elements are all key elements playing an important role in the network transmission function, and at the same time, once the elements fail, the elements are very easy to have a great influence on the own system or the coupled system, so that a large-scale failure is caused, and the safe operation of the system is endangered, for example, the natural gas pipeline P9-10、P8-9And an electric power line L1-2. Therefore, in the context of deep coupling of an electric power system and a natural gas system, system operators need to pay attention not only to the influence of the self system on the vulnerability of the elements, but also to the influence of the coupling system on the vulnerability of the elements.
According to the method, the comprehensive weakness of each element in the system is calculated quantitatively, the weak elements are found out and are subjected to key protection, so that the overall operation reliability and robustness of the system are improved, and the probability of large-area faults of the system is reduced.

Claims (5)

1. A method for improving the vulnerability of elements of an electricity-gas integrated energy system is characterized by comprising the following steps: the method comprises the following steps:
step 1, inputting original data, establishing an operation state model of the electricity-gas integrated energy system, and calculating to obtain an operation result of the system in an initial normal state;
step 2, extracting a network topological structure of the electricity-gas integrated energy system according to an operation result in an initial normal state, and processing to obtain structural vulnerability parameters of elements;
step 3, simulating element faults including power line outage and natural gas pipeline outage, and obtaining an operation result of the electricity-gas integrated energy system in a fault state by using an electricity-gas integrated energy system operation state model;
step 4, processing and obtaining a natural gas system airflow distribution entropy and a power system power flow distribution entropy based on an operation result in a fault state, so as to obtain an operation vulnerability parameter of the element;
step 5, calculating the weight of each parameter by using an optimal combination weight method based on each vulnerability parameter obtained in the step 2 and the step 4, thereby obtaining the comprehensive vulnerability of the element; selecting an element with higher comprehensive weakness for processing, and improving the reliability of the element by thickening a power line and reinforcing a natural gas pipeline;
the step 1 is specifically as follows:
step 1.1, establishing a minimum Minf of the sum of the operating energy consumption of the natural gas system and the power system as an objective function of an operating state model of the electricity-gas comprehensive energy system:
Figure FDA0003468477850000011
in the formula phig、ΦeRespectively a natural gas node set and a power grid node set;
Figure FDA0003468477850000012
respectively indicating the gas supply amount and the gas load removal amount of a gas source of a natural gas node i in a state l;
Figure FDA0003468477850000013
the energy consumption of the natural gas node i in the state l is
Figure FDA0003468477850000014
And
Figure FDA0003468477850000015
a quadratic or linear function of;
Figure FDA0003468477850000016
and
Figure FDA0003468477850000017
respectively obtaining the active output power of a conventional unit, the active output power of a gas unit and the cut-off quantity of the power load power of a power grid node j in a state l;
Figure FDA0003468477850000018
for the operating energy consumption of the grid node j in the state l
Figure FDA0003468477850000019
And
Figure FDA00034684778500000110
a quadratic or linear function of (a); i. j is the ordinal number of the natural gas node and the power grid node respectively; l is an ordinal number of the system running state, and the initial normal running state is represented when l is 0;
step 1.2, establishing equivalent constraint of the following electric-gas integrated energy system operation state model:
Figure FDA00034684778500000111
Figure FDA00034684778500000112
Figure FDA00034684778500000113
Figure FDA00034684778500000114
Figure FDA0003468477850000021
Figure FDA0003468477850000022
in the formula, Fload,i、
Figure FDA0003468477850000023
Respectively the gas load of a natural gas node i and the gas consumption of the gas unit under the state l; (i, t) denotes a natural gas pipeline connecting the head end node i and the tail end node t; phipipeRepresents a collection of natural gas pipelines;
Figure FDA0003468477850000024
representing the amount of gas flowing through the natural gas pipeline (i, t) in state l; mitRepresenting a transmission parameter of the natural gas pipeline (i, t);
Figure FDA0003468477850000025
respectively indicating that natural gas nodes i and t are in a state lThe air pressure square value of (a); pload,jRepresenting the active power load of the grid node j; (j, s) represents a power line connecting the head end node j and the tail end node s; philineRepresenting a set of power lines;
Figure FDA0003468477850000026
represents the active power flowing through the power line (j, s) in state l; x is the number ofjsRepresents the reactance of the power line (j, s);
Figure FDA0003468477850000027
and
Figure FDA0003468477850000028
respectively representing voltage phase angles of the grid nodes j and s in the state l; gHV、ηG2PRespectively representing the heat value of natural gas and the power generation efficiency of the gas engine set; j belongs to i and represents that the natural gas fuel of the gas turbine set on the power grid node j is provided by a natural gas node i; psour,jThe active power which is consumed by the gas source connected to the grid node j to maintain normal work is represented; etaP2GRepresenting the conversion coefficient of electricity to gas; i belongs to j and represents that the working power supply of the gas source on the natural gas node i is provided by the power grid node j;
step 1.3, establishing an inequality constraint of an electric-gas integrated energy system operation state model:
Figure FDA0003468477850000029
Figure FDA00034684778500000210
Figure FDA00034684778500000211
Figure FDA00034684778500000212
Figure FDA00034684778500000213
Figure FDA00034684778500000214
Figure FDA00034684778500000215
Figure FDA00034684778500000216
Figure FDA00034684778500000217
Figure FDA00034684778500000218
Figure FDA00034684778500000219
in the formula (I), the compound is shown in the specification,
Figure FDA00034684778500000220
the working state of an air source on a natural gas node i under the state l is shown, when the normal working value of the air source is 1, otherwise, the working state is 0; formula (II)
Figure FDA00034684778500000221
And
Figure FDA00034684778500000222
indicating when the electrical load shedding of grid node j is greater than a certain threshold, i.e.
Figure FDA00034684778500000223
Then
Figure FDA00034684778500000224
Figure FDA00034684778500000225
Indicating that the supplied gas source is shut down; the superscripts of the letters carry "max" and "min" respectively, which represent the upper and lower limits of the variable;
and finally, inputting the original data into the model to solve to obtain the running result of the system in the initial normal state.
2. The method of claim 1, wherein the method further comprises: the step 2 is specifically as follows:
step 2.1, according to the upper limit of the air supply quantity of the air source at the natural gas node i
Figure FDA0003468477850000031
Gas consumption of gas unit under initial normal operation state at natural gas node t
Figure FDA0003468477850000032
Gas load F at natural gas node tload,tProcessing and obtaining the capacity betweenness centrality parameter CBC of the natural gas pipeline according to the following formulap：
Figure FDA0003468477850000033
Figure FDA0003468477850000034
In the formula (I), the compound is shown in the specification,
Figure FDA0003468477850000035
representing a set of shortest paths between natural gas nodes i and t;
Figure FDA0003468477850000036
representing the weight of the mth shortest path between the nodes i and t;
Figure FDA0003468477850000037
representing the gas consumption of the gas unit of the natural gas node i in an initial normal operation state;
Figure FDA0003468477850000038
the transmission capacity of the mth shortest path between the nodes i and t is represented and is the minimum value of the transmission capacities of all the natural gas pipelines contained in the shortest path;
Figure FDA0003468477850000039
a decision variable indicating whether the natural gas pipeline p is in the shortest path, when the natural gas pipeline p is in the mth shortest path between the nodes i and t,
Figure FDA00034684778500000310
otherwise
Figure FDA00034684778500000311
m is the ordinal number of the shortest path between the natural gas nodes i and t; p is the ordinal number of the natural gas pipeline;
step 2.2, according to the upper limit of the air supply quantity of the air source at the natural gas node i
Figure FDA00034684778500000312
Gas consumption of gas unit in initial normal operation state at natural gas node i
Figure FDA00034684778500000313
Gas load F at natural gas node iload,iProcessing and obtaining the capability degree centrality parameter CDC of the natural gas pipeline according to the following formulap：
Figure FDA00034684778500000314
CDCp＝(dcap,p1+dcap,p2)/2
In the formula, c(i,t)Representing the transmission capacity of the natural gas pipeline (i, t); the t epsilon i represents that the natural gas node t is connected with the node i; p1 and p2 are respectively a head end node and a tail end node of the natural gas pipeline p;
step 2.3, according to the upper limit of the active output power of the conventional unit at the power grid node j
Figure FDA00034684778500000315
Gas unit active output power upper limit
Figure FDA00034684778500000316
Active power load P at grid node sload,sAnd the active power P consumed by the air source connected to the power grid node s for maintaining normal worksour,sThe following formula is adopted to process and obtain the electrical betweenness centrality parameter EBC of the power linek：
Figure FDA00034684778500000317
Figure FDA00034684778500000318
In the formula (I), the compound is shown in the specification,
Figure FDA0003468477850000041
representing the shortest path between grid nodes j, sA path set;
Figure FDA0003468477850000042
representing the weight of the nth shortest path between the nodes j and s;
Figure FDA0003468477850000043
the transmission capacity of the nth shortest path between the grid nodes j and s is represented and is the minimum value of the transmission capacities of all the power lines contained in the shortest path;
Figure FDA0003468477850000044
as a decision variable whether or not power line k is in the shortest path, when line k is in the nth shortest path between nodes j, s,
Figure FDA0003468477850000045
otherwise
Figure FDA0003468477850000046
n is the ordinal number of the shortest path between the grid nodes j and s; k is the ordinal number of the power line;
step 2.4, according to the upper limit of the active output power of the conventional unit at the power grid node j
Figure FDA0003468477850000047
Gas unit active output power upper limit
Figure FDA0003468477850000048
Active power load P at grid node jload,jAnd the active power P consumed by the air source connected to the power grid node j to maintain normal worksour,jProcessing and obtaining the electrical degree centrality parameter EDC of the power line according to the following formulak：
Figure FDA0003468477850000049
EDCk＝(dcapa,k1+dcapa,k2)/2
In the formula, c(j,s)Representing the transmission capacity of the power line (j, s); s belongs to j and represents that the power grid node s is connected with the node j; k1 and k2 are respectively the head end node and the tail end node of the power line k.
3. The method of claim 1, wherein the vulnerability of the element of the electric-gas integrated energy system is: the step 3 is specifically as follows:
step 3.1, according to the simulated natural gas pipeline outage fault, updating the network topology structure of the electricity-gas integrated energy system in the simulated natural gas pipeline outage fault process, inputting the updated network topology structure of the electricity-gas integrated energy system and original data into the operation state model of the electricity-gas integrated energy system to calculate and solve to obtain operation data of the system in the state, wherein the operation data comprises gas flow rate in the natural gas pipeline p fault state
Figure FDA00034684778500000410
And active power flow of natural gas pipeline p in fault state
Figure FDA00034684778500000411
And further calculating the airflow increment of the natural gas pipeline and the active power flow increment of the power line under the outage fault of the natural gas pipeline according to the following formulas according to the operation data:
Figure FDA00034684778500000412
Figure FDA00034684778500000413
in the formula (I), the compound is shown in the specification,
Figure FDA00034684778500000414
respectively due to the natural gas pipeline p being out of serviceThe fault causes the airflow increment of the natural gas pipeline q and the active power flow increment of the power line b;
Figure FDA00034684778500000415
and
Figure FDA00034684778500000416
respectively representing the transmission capacity of the natural gas pipeline q, the gas flow in the initial normal state and the gas flow in the failure state of the natural gas pipeline p;
Figure FDA00034684778500000417
and
Figure FDA00034684778500000418
respectively representing the transmission capacity of the power line b, the active power flow in the initial normal state and the active power flow in the fault state of the natural gas pipeline p;
step 3.2, simulating the outage fault of the power line, updating the network topology structure of the electric-gas integrated energy system in the process of simulating the outage fault of the power line, inputting the updated network topology structure of the electric-gas integrated energy system and original data into the operation state model of the electric-gas integrated energy system, calculating and solving to obtain the operation data of the system in the state, wherein the operation data comprises the gas flow rate of the natural gas pipeline q in the k fault state of the power line
Figure FDA0003468477850000051
And the active power flow of the power line b in the fault state of the power line k
Figure FDA0003468477850000052
And further calculating the airflow increment of the natural gas pipeline and the active power flow increment of the power line under the outage fault of the power line according to the following formulas according to the operation data:
Figure FDA0003468477850000053
Figure FDA0003468477850000054
in the formula (I), the compound is shown in the specification,
Figure FDA0003468477850000055
respectively representing the airflow increment of the natural gas pipeline q and the active power flow increment of the power line b caused by the disconnection fault of the power line k;
Figure FDA0003468477850000056
representing the gas flow of the natural gas pipeline q in the fault state of the power line k;
Figure FDA0003468477850000057
representing the active power flow of the power line b in a fault condition of the power line k.
4. The method of claim 1, wherein the vulnerability of the element of the electric-gas integrated energy system is: the step 4 is specifically as follows:
step 4.1, processing according to the airflow increment of the natural gas pipeline under the outage fault of the natural gas pipeline and the power line and the following formula to obtain the vulnerability parameter PFI of the outage fault of the natural gas pipelinep：
Figure FDA0003468477850000058
In the formula (I), the compound is shown in the specification,
Figure FDA0003468477850000059
indicating the gas flow rate of the natural gas pipeline p in an initial normal state;
Figure FDA00034684778500000510
respectively representing natural gas systems caused by shutdown faults of natural gas pipelines pThe system airflow distribution entropy and the power flow distribution entropy of the power system are calculated; u. of1、u2Respectively representing the influence weight of the natural gas pipeline outage fault on a natural gas system and an electric power system;
the airflow distribution entropy and the load flow distribution entropy caused by the natural gas pipeline p fault are calculated as follows:
Figure FDA00034684778500000511
Figure FDA00034684778500000512
Figure FDA00034684778500000513
Figure FDA00034684778500000514
in the formula etaqp、ηbpRespectively representing the airflow impact ratio caused by the shutdown fault of the natural gas pipeline p to the natural gas pipeline q and the power flow impact ratio caused by the shutdown fault of the natural gas pipeline p to the line b;
Figure FDA00034684778500000515
respectively representing the airflow increment of the natural gas pipeline q and the active power flow increment of the power line b caused by the shutdown fault of the natural gas pipeline p;
step 4.2, processing according to the active power flow increment of the power line under the outage faults of the natural gas pipeline and the power line and the following formula to obtain the fault vulnerability parameter LFI of the power linek：
Figure FDA0003468477850000061
In the formula，
Figure FDA0003468477850000062
Representing the active power flow power of the power line k in an initial normal state;
Figure FDA0003468477850000063
respectively representing the power flow distribution entropy of the power system and the airflow distribution entropy of the natural gas system caused by the k fault of the line; v. of1、v2Respectively representing the influence weight of the line fault on the power system and the natural gas system;
the power flow distribution entropy and the airflow distribution entropy caused by the fault of the line k are calculated as follows:
Figure FDA0003468477850000064
Figure FDA0003468477850000065
Figure FDA0003468477850000066
Figure FDA0003468477850000067
in the formula etabk、ηqkRespectively representing the power flow impact ratio caused by the disconnection fault of the line k to the line b and the airflow impact ratio caused by the disconnection fault of the line k to the natural gas pipeline q;
Figure FDA0003468477850000068
respectively representing the increment of the airflow of the natural gas pipeline q and the increment of the active power flow of the power line b caused by the disconnection fault of the power line k.
5. The method of claim 1, wherein the vulnerability of the element of the electric-gas integrated energy system is: the step 5 is specifically as follows:
step 5.1, establishing the following optimal combination weight model, and calculating the optimal weight of each vulnerability parameter:
Figure FDA0003468477850000069
ρs,β＝τs,β/(τs,β+τo,β)
ρo,β＝τo,β/(τs,β+τo,β)
Figure FDA00034684778500000610
wβ≥0
in the formula, wβAn optimal weight representing the parameter β; tau iss,β、τo,βRespectively representing subjective weight obtained by calculating the parameter beta according to an analytic hierarchy process and objective weight obtained by calculating according to an entropy weight method; rhos,β、ρo,βBias coefficients representing a parameter beta supervisor weight and an objective weight, respectively; m is the total number of all vulnerability parameters; beta is the ordinal number of the vulnerability parameter;
and 5.2, integrating the structural vulnerability parameters and the operation vulnerability parameters obtained in the step 2 and the step 4, and obtaining the integrated vulnerability of each natural gas pipeline and each power line by adopting the following formula:
VG(p)＝ωg1CBCp+ωg2CDCp+ωg3PFIp
VE(k)＝ωe1EBCk+ωe2EDCk+ωe3LFIk
in the formula, ωg1、ωg2And ωg3Are respectively CBCp、CDCpAnd PFIpThe parameters are calculated according to the step 5.1 to obtain the optimal weight; omegae1、we2And ωe3Are respectively EBCk、EDCkAnd LFIkThe parameters are calculated according to the step 5.1 to obtain the optimal weight; CBCpRepresenting the capability mesomeric centrality parameter, CDC, of a natural gas pipelinepRepresenting the power degree centrality parameter, PFI, of a natural gas pipelinepRepresenting the vulnerability parameter of the natural gas pipeline to out-of-service faults, EBCkIndicating the electrical neutral center parameter, EDC, of the power linekIndicating the electrical degree centrality parameter, LFI, of the power linekA fault vulnerability parameter representing a power line fault;
step 5.3, calculating the comprehensive vulnerability value of each natural gas pipeline/power line element according to the step 5.2, and improving the vulnerability of the electricity-gas comprehensive energy system according to the comprehensive vulnerability value specifically comprises the following steps: for the elements with the comprehensive weakness degree higher than the preset weakness degree threshold value, if the elements are power lines, measures of thickening the power lines, increasing the splitting number of the lines and replacing the elements with double-circuit lines are adopted; if the element is a natural gas pipeline, measures of thickening the wall of the natural gas pipeline and reducing the gas transmission pressure of the pipeline are taken.
CN202010135812.2A 2020-03-02 2020-03-02 Method for improving vulnerability of elements of electricity-gas comprehensive energy system Active CN111444593B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010135812.2A CN111444593B (en) 2020-03-02 2020-03-02 Method for improving vulnerability of elements of electricity-gas comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010135812.2A CN111444593B (en) 2020-03-02 2020-03-02 Method for improving vulnerability of elements of electricity-gas comprehensive energy system

Publications (2)

Publication Number Publication Date
CN111444593A CN111444593A (en) 2020-07-24
CN111444593B true CN111444593B (en) 2022-05-03

Family

ID=71652649

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010135812.2A Active CN111444593B (en) 2020-03-02 2020-03-02 Method for improving vulnerability of elements of electricity-gas comprehensive energy system

Country Status (1)

Country Link
CN (1) CN111444593B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112260261B (en) * 2020-09-08 2022-05-27 浙江大学 Method for evaluating and improving vulnerability of electricity-gas comprehensive energy system
CN112434442B (en) * 2020-12-08 2022-06-03 湘潭大学 Elasticity evaluation method for electricity-gas region comprehensive energy system
CN112861292B (en) * 2021-01-12 2022-08-05 浙江大学 Recovery improvement method for electricity-gas comprehensive energy system
CN113410842B (en) * 2021-06-18 2022-07-05 天津大学 Method and device for calculating probability energy flow of comprehensive energy system considering linkage fault
CN113420937B (en) * 2021-07-05 2022-12-30 山西大学 Comprehensive energy switching method and system considering vulnerability of conveying path

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934246A (en) * 2017-03-21 2017-07-07 广东电网有限责任公司惠州供电局 The computational methods and device of power network line fragility
CN107622360A (en) * 2017-10-20 2018-01-23 广东电网有限责任公司电力调度控制中心 A kind of critical circuits recognition methods for considering subjective and objective factor
CN107623319A (en) * 2017-08-17 2018-01-23 广东电网有限责任公司惠州供电局 A kind of power network critical circuits discrimination method based on more evaluation indexes
CN110096764A (en) * 2019-04-12 2019-08-06 浙江大学 A kind of identification of electric-gas coupled system vulnerable line and optimization method
CN110705879A (en) * 2019-09-30 2020-01-17 国网山东省电力公司滨州供电公司 Power grid vulnerability assessment method under high-proportion renewable energy access

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170288253A1 (en) * 2016-03-29 2017-10-05 The Board Of Trustees Of The Leland Stanford Junior University Thermoelectrochemical Heat Converter

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934246A (en) * 2017-03-21 2017-07-07 广东电网有限责任公司惠州供电局 The computational methods and device of power network line fragility
CN107623319A (en) * 2017-08-17 2018-01-23 广东电网有限责任公司惠州供电局 A kind of power network critical circuits discrimination method based on more evaluation indexes
CN107622360A (en) * 2017-10-20 2018-01-23 广东电网有限责任公司电力调度控制中心 A kind of critical circuits recognition methods for considering subjective and objective factor
CN110096764A (en) * 2019-04-12 2019-08-06 浙江大学 A kind of identification of electric-gas coupled system vulnerable line and optimization method
CN110705879A (en) * 2019-09-30 2020-01-17 国网山东省电力公司滨州供电公司 Power grid vulnerability assessment method under high-proportion renewable energy access

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A Multi-State Model for Reliability Assessment of Integrated Gas and Power Systems Utilizing Universal Generating Function Techniques;Minglei Bao et al;《 IEEE Transactions on Smart Grid》;20191130;第10卷(第6期);全文 *
Simplified operation models of integrated power and gas systems for vulnerability analysis;Hui Zhang et al;《Physica A: Statistical Mechanics and its Applications》;20191031;第531卷;全文 *
基于受冲击与断开后果脆弱度的电网关键线路识别;倪良华 等;《电力系统保护与控制》;20200101;第48卷(第1期);全文 *
复杂电力系统脆弱性评估方法研究;谭玉东;《中国博士学位论文全文数据库电子期刊 工程科技II辑》;20140915;第2014年卷(第9期);全文 *
考虑天然气网影响的电网脆弱线路辨识;桑茂盛 等;《电力系统自动化》;20191110;第43卷(第21期);全文 *

Also Published As

Publication number Publication date
CN111444593A (en) 2020-07-24

Similar Documents

Publication Publication Date Title
CN111444593B (en) Method for improving vulnerability of elements of electricity-gas comprehensive energy system
CN110210659B (en) Power distribution network planning method considering reliability constraint
CN108631306B (en) Method for evaluating recovery capability of power system after disaster
CN112260261B (en) Method for evaluating and improving vulnerability of electricity-gas comprehensive energy system
CN110096764B (en) Method for identifying and optimizing fragile line of electric-gas coupling system
CN112736923B (en) Natural gas network cascading failure assessment control method considering power grid failure influence
CN112271726B (en) Power distribution system fault recovery method considering electricity-water-gas coupling relation
CN113742917B (en) Comprehensive energy system toughness improvement method considering multi-stage recovery process
CN113346484B (en) Power distribution network elastic lifting method and system considering transient uncertainty
CN107622360A (en) A kind of critical circuits recognition methods for considering subjective and objective factor
Jin et al. Planning for distribution resilience under variable generation: prevention, surviving and recovery
CN115345391B (en) Post-disaster recovery method and device for electric-gas energy system and storage medium
CN111126841A (en) Power distribution network load to energy conversion power evaluation method based on random power flow
Gibson et al. MID-LIFE REFURBISHMENT OF AN HVDC INTERCONNECTOR CONTROL SYSTEM
CN113094854B (en) Backbone net rack reconstruction method for electric-gas coupling system
Tian et al. Resilience-based optimal placement method for integrated electricity and gas energy system
Tuinema et al. Network redundancy versus generation reserve in combined onshore-offshore transmission networks
CN115511171A (en) Two-stage fault recovery and maintenance method for regional electricity-gas comprehensive energy system
Niu et al. Multistage Condition-Based Resilient Maintenance Scheduling for Power Systems during Ice Storms
Gui et al. Transmission Adaptive Capacity-Based Resilience Metrics for Power Grid Contingency Analysis
CN117251971A (en) Distributed load recovery method for multi-energy coupling power distribution network
CN112861292B (en) Recovery improvement method for electricity-gas comprehensive energy system
CN116894342B (en) Toughness improving method for electric-gas comprehensive energy system based on natural gas network pipe storage
Bai et al. A Multi-Stage Resilient Optimization Method for Coupled Electric-Gas System Based on Island Division
CN115510618A (en) Two-stage toughness recovery method for electrical coupling system considering component importance degree identification

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
GR01 Patent grant
GR01 Patent grant