A method and system for optimizing partitioning of an electric-gas integrated energy system

Abstract

The present invention discloses a method and system for optimizing the zoning of an integrated electricity-gas energy system, comprising: constructing a multi-level evaluation index system based on the internal links and external influences of the integrated electricity-gas energy system. Calculating evaluation indicators, and constructing objective functions and constraints. Establishing an optimization model for the zoning of an integrated electricity-gas energy system based on internal and external operational risks, and solving the model to obtain the optimal post-disaster zoning results for the integrated electricity-gas energy system. The method and system for optimizing the zoning of an integrated electricity-gas energy system provided by the present invention take into account the topological structure and resource abundance of the gas transmission system, as well as the internal and external risks that the system may face during post-disaster recovery. This helps to formulate a more reasonable and effective post-disaster zoning strategy, and can achieve parallel recovery of integrated electricity-gas energy systems in different regions, accelerating the entire post-disaster recovery process while maintaining stable energy supply.

Description

Technical Field

[0001] The present invention relates to the technical field of electricity-gas integrated energy, and in particular to a partition optimization method and system for an electricity-gas integrated energy system. Background Art

[0002] Extreme natural disasters, including typhoons and earthquakes, have become increasingly frequent in recent years. These events pose significant threats to the operational safety and stability of integrated power and gas energy systems, making them prone to energy outages. Therefore, accelerating the post-disaster recovery of integrated power and gas energy systems and minimizing losses are of great economic and social significance.

[0003] By developing a rational post-disaster zoning strategy, it is possible to achieve parallel recovery in multiple regions, accelerating the post-disaster recovery of the integrated power-gas energy system. However, existing research has focused solely on developing post-disaster zoning strategies for the power system, while insufficiently addressing the optimization of post-disaster zoning for the integrated power-gas energy system. Furthermore, as the coupling between power and gas transmission systems becomes increasingly tight, their scale and topology are becoming larger and more complex. Therefore, balancing the supply of power and gas resources and accelerating recovery after a disaster is a key issue in developing an optimal post-disaster zoning strategy for the integrated power-gas energy system. Summary of the Invention

[0004] In view of the above-mentioned problems, the present invention is proposed.

[0005] Therefore, the technical problem addressed by this invention is that existing research has focused solely on developing post-disaster zoning strategies for power systems, but has insufficiently addressed the optimization of post-disaster zoning for integrated power and gas energy systems. Furthermore, as the coupling between power and gas transmission systems becomes increasingly close, their scale and topology are becoming larger and more complex.

[0006] To solve the above technical problems, the present invention provides the following technical solutions: a method for optimizing the zoning of an electricity-gas integrated energy system, comprising: constructing a multi-level evaluation index system based on the internal links and external influences of the electricity-gas integrated energy system.

[0007] Calculate evaluation indicators and construct objective functions and constraints.

[0008] An electricity-gas integrated energy system zoning optimization model based on internal and external operation risks is established, and the optimal post-disaster zoning results of the electricity-gas integrated energy system are obtained by solving the model.

[0009] As a preferred embodiment of the zoning optimization method for the integrated electricity-gas energy system described in the present invention, the multi-level evaluation index system includes three main categories of indicators: energy supply risk indicator, pipeline network structure risk indicator, and post-disaster environmental risk indicator. After classification, each main indicator is further expanded into a sub-indicator layer, each of which contains two sub-indicators.

[0010] As a preferred solution of the method for optimizing the partitioning of the electric-gas integrated energy system described in the present invention, the energy supply risk index includes the power disaster resource abundance R 1 and natural gas post-disaster resource abundance R 2 The pipeline network structure risk indicators include the system pipeline length R 3 and system connectivity R 4 , post-disaster environmental risk indicators include post-disaster tower failure rate R 5 and post-disaster pipeline damage rate R 6 .

[0011] As a preferred solution of the method for optimizing the partitioning of the electric-gas integrated energy system described in the present invention, the calculation evaluation index includes the resource abundance after the power disaster. Expressed as:

[0012]

[0013] in, Represents the set of all nodes in the integrated energy system partition z. Represents the 0-1 partition state variable of node m, when When node m is divided into area z. represents the total power load of node m in partition z, represents the total power generation capacity of node m in partition z. It represents the power consumption of the electric drive gas component of node m in partition z at rated output. represents the post-disaster failure rate of the unit at node m in partition z. When node m has no power generation resources, represents the post-disaster unrepairability rate of the power load of node m in partition z.

[0014] Natural gas resource abundance after a disaster Expressed as:

[0015]

[0016] in, represents the total gas load of node m in the integrated energy system partition z, Represents the total gas supply capacity of node m in zone z of the integrated energy system. Indicates the gas consumption of the gas-driven electrical component at node m in partition z at rated output. represents the post-disaster failure rate of the gas supply device of node m in partition z. When node m has no gas supply resources, It represents the post-disaster irreparable rate of the gas load of node m in partition z.

[0017] System pipeline length indicator Expressed as:

[0018]

[0019] in, It is represented as the shortest geographical distance between node m and the black start unit in the integrated energy system zone z.

[0020] System connectivity indicators Expressed as:

[0021]

[0022] in, Represents the pipeline set of the integrated energy system partition z. It is represented as the 0-1 partition state variable of pipeline k, when When pipeline k is divided into area z.

[0023] Post-disaster tower failure rate index Expressed as:

[0024]

[0025] in, It is expressed as the post-disaster tower failure rate of line k in the integrated energy system zone z.

[0026] Post-disaster pipeline damage rate indicator Expressed as:

[0027]

[0028] in, It is expressed as the post-disaster pipeline damage rate of pipeline k in zone z of the integrated energy system.

[0029] As a preferred solution of the method for optimizing the partitioning of the electric-gas integrated energy system of the present invention, wherein: the constructing of the objective function includes: the objective function of the optimization model for the partitioning of the electric-gas integrated energy system is to minimize the comprehensive evaluation index R:

[0030]

[0031] Where, β 1represents the weight of the power resource abundance index after a disaster, β 2 represents the weight of natural gas post-disaster resource abundance index, β 3 represents the weight of the system pipeline length indicator, β 4 System connectivity index weight, β 5 represents the weight of the tower failure rate index after the disaster, β 6 represents the weight of the pipeline damage rate indicator after the disaster. Z represents the total number of zones in the electricity-gas integrated energy system.

[0032] As a preferred solution of the electric-gas integrated energy system zoning optimization method described in the present invention, the constraints include: zoning connectivity constraints for constructing post-disaster zoning of the electric-gas integrated energy system, zoning classification constraints, transmission line tower failure constraints under typhoon weather, and gas pipeline damage constraints under earthquake disasters.

[0033] The partition connectivity constraint constructed based on the state constraint virtual network flow method is expressed as:

[0034]

[0035] in, Indicates the black start unit in the zone z of the electric-gas integrated energy system, Represents the set of pipelines connected to node m in partition z of the electricity-gas integrated energy system. It is represented as the set of nodes in partition z that are connected to pipeline k. a 0-1 partition state variable representing the end of pipeline k in partition z, and Represents the 0-1 partition state variable of the head-end node of pipeline k in partition z. represents the virtual network flow in pipeline k with the black start unit as the head node in partition z represents the virtual network flow in pipeline k with the black start unit as the terminal node in partition z. M is a sufficiently large constant.

[0036] The partition classification constraint based on the state constraint virtual network flow method is expressed as:

[0037]

[0038] in, Indicates the black start unit in the zone z of the electric-gas integrated energy system. Represents the gas supply node in the electricity-gas integrated energy system zone z, Represents the electric-driven gas components in the electric-gas integrated energy system zone z, Represents the set of gas generator sets in the electricity-gas integrated energy system zone z. N min Indicates the minimum number of nodes allowed for each partition.

[0039] The failure constraint of transmission line towers under typhoon weather is expressed as:

[0040]

[0041] in, It represents the actual wind speed that line k in zone z of the electricity-gas integrated energy system bears. It represents the design wind speed that line k in zone z of the electric-gas integrated energy system can withstand. k,z Represents the tower failure rate model coefficient of line k in zone z.

[0042] The damage constraint of gas pipelines under earthquake disasters is expressed as:

[0043]

[0044] Among them, P k,z represents the total number of post-disaster damaged points of pipeline k in the electricity-gas integrated energy system zone z, κ k,z It represents the number and total length of pipeline k in zone z of the electricity-gas integrated energy system.

[0045] As a preferred solution of the electricity-gas integrated energy system zoning optimization method described in the present invention, the solution model obtains the optimal post-disaster zoning result of the electricity-gas integrated energy system, including inputting the optimization model and the topological structure and parameters of the actual electricity-gas integrated energy system into the solver to calculate the optimal post-disaster zoning result of the electricity-gas integrated energy system.

[0046] An electric-gas integrated energy system partition optimization system, characterized by comprising:

[0047] Construct an evaluation index system module and build a multi-level evaluation index system based on the internal links and external influences of the electricity-gas integrated energy system.

[0048] The calculation module calculates the evaluation indicators and constructs the objective function and constraints.

[0049] The solution module establishes an optimization model for the zoning of the electricity-gas integrated energy system based on internal and external operation risks, and solves the model to obtain the optimal post-disaster zoning result of the electricity-gas integrated energy system.

[0050] A computer device includes a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above method when executing the computer program.

[0051] A computer-readable storage medium stores a computer program, which implements the steps of the method described above when executed by a processor.

[0052] Beneficial effects of the present invention: In response to the serious consequences caused by large-scale energy outages in electricity-gas integrated energy systems under various extreme natural disasters, the present invention proposes a zoning optimization method for electricity-gas integrated energy systems based on internal and external operational risks. Compared with the post-disaster zoning optimization method that only considers the topological structure of the power system, the post-disaster zoning optimization method for the electricity-gas integrated energy system proposed by the present invention can not only take into account the topological structure and resource abundance of the gas transmission system, but also take into account the internal and external risks that the system may face during the post-disaster recovery period. It helps to formulate a more reasonable and effective post-disaster zoning strategy, and can realize the parallel recovery of electricity-gas integrated energy systems in different regions, accelerating the entire post-disaster recovery process while maintaining stable energy supply. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] To more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for describing the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. Those skilled in the art can also derive other drawings based on these drawings without inventive effort. Among them:

[0054] Figure 1 The first embodiment of the present invention provides an overall flow chart of an electric-gas integrated energy system partition optimization method and system.

[0055] Figure 2 A topological diagram of a 125-node electric system and a 49-node natural gas system in a certain region of China is provided for a method and system for optimizing the zoning of an electric-gas integrated energy system according to a second embodiment of the present invention. DETAILED DESCRIPTION

[0056] To make the above-mentioned objects, features, and advantages of the present invention more clearly understood, the following detailed description of the specific embodiments of the present invention is given in conjunction with the accompanying drawings. It is obvious that the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in this field without creative work should fall within the scope of protection of the present invention.

[0057] Example 1, reference Figure 1 , as one embodiment of the present invention, provides a method for optimizing the partitioning of an electric-gas integrated energy system, comprising:

[0058] S1: Construct a multi-level evaluation index system based on the internal links and external influences of the electricity-gas integrated energy system.

[0059] It should be noted that a hierarchical integrated energy system risk indicator system has been established, adhering to the principles of scientificity, pertinence, independence, practicality, and comparability. Based on both the internal links and external influences of the electricity-gas integrated energy system, this system assesses whether its internal energy supply links and energy transmission network structure pose risks that could affect the system's stable energy supply recovery, post-disaster recovery efficiency, or even lead to further outages during post-disaster recovery. Furthermore, it assesses whether continued environmental changes after a disaster will create additional recovery risks for key links in the integrated energy system.

[0060] The multi-level evaluation index system assesses energy supply risk, pipeline network structure risk, and post-disaster environmental risk as three main categories. Each of these categories is further categorized into sub-indicator layers, each containing two sub-indicators. Table 1 shows the overall framework of the post-disaster risk index system for the integrated electricity and gas energy system.

[0061] Table 1 Overall framework of the post-disaster risk indicator system for the electricity-gas integrated energy system

[0062]

[0063] S2: Calculate evaluation indicators and construct objective functions and constraints.

[0064] Post-disaster power resource adequacy Expressed as:

[0065]

[0066] in, Represents the set of all nodes in the integrated energy system partition z. Represents the 0-1 partition state variable of node m, when When node m is divided into area z. represents the total power load of node m in partition z, represents the total power generation capacity of node m in partition z. It represents the power consumption of the electric drive gas component of node m in partition z at rated output. represents the post-disaster failure rate of the unit at node m in partition z. When node m has no power generation resources, represents the post-disaster unrepairability rate of the power load of node m in partition z.

[0067] Natural gas resource abundance after a disaster Expressed as:

[0068]

[0069] in, represents the total gas load of node m in the integrated energy system partition z, Represents the total gas supply capacity of node m in zone z of the integrated energy system. Indicates the gas consumption of the gas-driven electrical component at node m in partition z at rated output. represents the post-disaster failure rate of the gas supply device of node m in partition z. When node m has no gas supply resources, It represents the post-disaster irreparable rate of the gas load of node m in partition z.

[0070] System pipeline length indicator Expressed as:

[0071]

[0072] in, It is represented as the shortest geographical distance between node m and the black start unit in the integrated energy system zone z.

[0073] System connectivity indicators Expressed as:

[0074]

[0075] in, Represents the pipeline set of the integrated energy system partition z. It is represented as the 0-1 partition state variable of pipeline k, when When pipeline k is divided into area z.

[0076] Post-disaster tower failure rate index Expressed as:

[0077]

[0078] in, It is expressed as the post-disaster tower failure rate of line k in the integrated energy system zone z.

[0079] Post-disaster pipeline damage rate indicator Expressed as:

[0080]

[0081] in, It is expressed as the post-disaster pipeline damage rate of pipeline k in zone z of the integrated energy system.

[0082] 5. The method for optimizing the zoning of an integrated electricity-gas energy system according to claim 4, wherein the objective function is constructed such that the objective function of the zoning optimization model of the integrated electricity-gas energy system is to minimize the comprehensive evaluation index R:

[0083]

[0084] Where, β 1 represents the weight of the power resource abundance index after a disaster, β 2 represents the weight of natural gas post-disaster resource abundance index, β 3 represents the weight of the system pipeline length indicator, β 4 System connectivity index weight, β 5 represents the weight of the tower failure rate index after the disaster, β 6 represents the weight of the pipeline damage rate indicator after the disaster. Z represents the total number of zones in the electricity-gas integrated energy system.

[0085] The objective function of the zoning optimization model of the electricity-gas integrated energy system is to minimize the comprehensive evaluation index R:

[0086]

[0087] Where, β 1 represents the weight of the power resource abundance index after a disaster, β 2 represents the weight of natural gas post-disaster resource abundance index, β 3 represents the weight of the system pipeline length indicator, β 4 System connectivity index weight, β 5 represents the weight of the tower failure rate index after the disaster, β 6 represents the weight of the pipeline damage rate indicator after the disaster. Z represents the total number of zones in the electricity-gas integrated energy system.

[0088] Construct post-disaster zoning of the electricity-gas integrated energy system based on zoning connectivity constraints, zoning classification constraints, transmission line tower failure constraints under typhoon weather, and gas pipeline damage constraints under earthquake disasters.

[0089] The partition connectivity constraint constructed based on the state constraint virtual network flow method is expressed as:

[0090]

[0091] in, Indicates the black start unit in the zone z of the electric-gas integrated energy system, Represents the set of pipelines connected to node m in partition z of the electricity-gas integrated energy system. It is represented as the set of nodes in partition z that are connected to pipeline k. a 0-1 partition state variable representing the end of pipeline k in partition z, and Represents the 0-1 partition state variable of the head-end node of pipeline k in partition z. represents the virtual network flow in pipeline k with the black start unit as the head node in partition z represents the virtual network flow in pipeline k with the black start unit as the terminal node in partition z. M is a sufficiently large constant.

[0092] The partition classification constraint based on the state constraint virtual network flow method is expressed as:

[0093]

[0094] in, Indicates the black start unit in the zone z of the electric-gas integrated energy system. Represents the gas supply node in the electricity-gas integrated energy system zone z, Represents the electric-driven gas components in the electric-gas integrated energy system zone z, Represents the set of gas generator sets in the electricity-gas integrated energy system zone z. N min Indicates the minimum number of nodes allowed for each partition.

[0095] The failure constraint of transmission line towers under typhoon weather is expressed as:

[0096]

[0097] in, It represents the actual wind speed that line k in zone z of the electricity-gas integrated energy system bears. It represents the design wind speed that line k in zone z of the electric-gas integrated energy system can withstand. k,z Represents the tower failure rate model coefficient of line k in zone z.

[0098] The damage constraint of gas pipelines under earthquake disasters is expressed as:

[0099]

[0100] Among them, P k,z represents the total number of post-disaster damaged points of pipeline k in the electricity-gas integrated energy system zone z, κ k,z It represents the number and total length of pipeline k in zone z of the electricity-gas integrated energy system.

[0101] S3: Establish an optimization model for the zoning of the electricity-gas integrated energy system based on internal and external operational risks, and solve the model to obtain the optimal post-disaster zoning results of the electricity-gas integrated energy system.

[0102] It should be noted that the objective function is to minimize the combined weighted value of the aforementioned post-disaster power resource abundance index, post-disaster natural gas resource abundance index, system pipeline length index, system connectivity index, post-disaster tower failure rate index, and post-disaster pipeline damage rate index. Furthermore, the aforementioned partition connectivity constraints, partition classification constraints, transmission line tower failure constraints under typhoon weather, and gas pipeline damage constraints under earthquake disasters are used as constraints to construct an electricity-gas integrated energy system partition optimization model based on internal and external operational risks. The model is solved using the Gurobi 10.0.0 solver in MATLAB software to obtain the optimal post-disaster partition results for the electricity-gas integrated energy system. The construction and solution process of the electricity-gas integrated energy system partition optimization model based on internal and external operational risks is as follows:

[0103] Step 1: Based on the internal links and external impacts of the electricity-gas integrated energy system, a multi-level evaluation indicator system is constructed with evaluation indicators such as energy supply, pipeline network structure and post-disaster environmental risks as the framework.

[0104] Step 2: Construct the post-disaster power resource abundance index, post-disaster natural gas resource abundance index, system pipeline length index, system connectivity index, post-disaster tower failure rate index, and post-disaster pipeline damage rate index, and assign weights to these indicators.

[0105] Step 3: Construct regional connectivity constraints, regional classification constraints, transmission line tower failure constraints under typhoon weather, and gas pipeline damage constraints under earthquake disasters.

[0106] Step 4: Use the indicators and their weights constructed in step 1 as the objective function, and then use the constraints obtained in step 2 as constraints to construct an electricity-gas integrated energy system zoning optimization model based on internal and external operation risks.

[0107] Step 5: Input the optimization model obtained in step 3 and the topology and parameters of the actual electric-gas integrated energy system into MATLAB software.

[0108] Step 6: Solve the model using the Gurobi 10.0.0 solver in MATLAB software, and calculate the optimal post-disaster partitioning result of the electricity-gas integrated energy system according to formula (7). The above embodiment also includes an electricity-gas integrated energy system partition optimization system, specifically:

[0109] Construct an evaluation index system module and build a multi-level evaluation index system based on the internal links and external influences of the electricity-gas integrated energy system.

[0110] The calculation module calculates the evaluation indicators and constructs the objective function and constraints.

[0111] The solution module establishes an optimization model for the zoning of the electricity-gas integrated energy system based on internal and external operation risks, and solves the model to obtain the optimal post-disaster zoning result of the electricity-gas integrated energy system.

[0112] The computer device may be a server. The computer device includes a processor, a memory, an input / output interface (I / O) and a communication interface. The processor, the memory and the input / output interface are connected via a system bus, and the communication interface is connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store data cluster data of the power monitoring system. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, a method for optimizing the partitioning of an electric-gas integrated energy system is implemented.

[0113] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, database or other media used in the embodiments provided in this application may include at least one of non-volatile and volatile memory. Non-volatile memory may include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory may include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in the various embodiments provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, distributed databases based on blockchains. The processor involved in the various embodiments provided in this application may be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic unit, a data processing logic unit based on quantum computing, etc., but are not limited to these.

[0114] Example 2, reference Figure 2 , which is an embodiment of the present invention, provides a method and system for optimizing the partitioning of an electric-gas integrated energy system. In order to verify the beneficial effects of the present invention, scientific demonstration is carried out through simulation experiments.

[0115] This paper uses an actual 125-node power system (represented by black lines) and a 49-node natural gas system (represented by red lines) in a certain region of China for example analysis. The topology is shown in the figure below. Figure 2 In addition, in this embodiment, only typhoon disasters with a high frequency occurred, and no earthquake disasters occurred (the pipeline damage rate after the disaster was zero).

[0116] The optimal post-disaster zoning strategy for the electricity-gas integrated energy system solved by this invention is as follows: Figure 2 As shown in the figure, under the optimal post-disaster zoning strategy, there are a total of 20 connecting pipelines between zones, and another four lines cannot be included in the recovery zone due to the continued typhoon disaster after the outage (represented by the black dashed line). Node 121 is also affected by the typhoon and cannot be included in any recovery zone. No gas pipelines are excluded from the recovery zone due to the earthquake. The average distance between each node and the black start unit is 25.47 km. The absolute difference between power generation capacity and power load is 0.38 pu. The absolute difference between natural gas supply and gas demand is 7840 m3. The wind speed experienced by each transmission line is 19.63 m / s.

[0117] Likewise Figure 2 As shown in the figure, according to the number of black start units, the electric-gas integrated energy system is divided into 5 zones in total. More specifically, zone 1 contains 1 black start unit and 1 non-black start unit, which are located at nodes 1 and 13 of the transmission network in zone 1, respectively, while 2 gas sources and 1 compressor are located at nodes 132, 139 and 133 of the gas supply system, respectively; zone 2 has only 1 black start unit and 1 electric drive compressor, which are located at nodes 15 and 60, respectively; 1 black start unit, 2 non-black start units, 1 electric drive compressor, 1 wind turbine unit, 1 photovoltaic unit and 1 gas source are located at nodes 132, 139 and 133 of the gas supply system, respectively. In zone 3, nodes 19, 28, 89, 107, 78, and 74 are located, while two battery storage units are located at nodes 76 and 84, respectively. In zone 4, one blackstart unit is located at node 36, one P2G unit is located at node 102, and one compressor and one gas storage tank are located at nodes 143 and 146, respectively. Finally, in zone 5, one blackstart unit and three non-blackstart units are located at nodes 30, 34, 42, and 110, respectively, and one gas source is located at node 153. It should be noted that the three blackstart units in zones 1, 4, and 5 are all gas-fired units, as are the three non-blackstart units at nodes 42, 44, and 110. Among all the five partitions, the gas supply nodes and gas pipelines of the gas supply system are mainly allocated to three of them, because all three partitions contain at least one gas unit as a black start unit, and the gas black start unit requires natural gas supply to increase power generation to support all recovery phases of the electric-gas integrated energy system. Therefore, most natural gas pipelines and nodes should be allocated to the partitions with gas black start units to provide a stable and sufficient supply of natural gas. Some electric-driven gas components, such as P2G devices and electric-driven compressors, are allocated to other partitions to avoid a lack of sufficient power to maintain the normal operation of the gas supply system during the recovery period. In summary, the electric-gas integrated energy system partition optimization method based on connectivity and classification constraints proposed in the present invention can calculate a reasonably balanced partitioning strategy that takes into account the respective topological structures and resource abundances of the power system and the natural gas system.

[0118] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. A method for optimizing the partitioning of an electricity-gas integrated energy system, characterized in that: include: Construct a multi-level evaluation index system based on the internal links and external influences of the electricity-gas integrated energy system; The multi-level evaluation index system includes: the multi-level evaluation index system uses energy supply risk index, pipeline network structure risk index, and post-disaster environmental risk index as three main categories of indicators for evaluation; After classification, each main indicator is expanded downward to a sub-indicator layer, and each sub-indicator layer contains two sub-indicators; Calculate evaluation indicators, construct objective functions and constraints; The calculation evaluation index includes: the power resource adequacy after the disaster Expressed as: in, represents the set of all nodes in the integrated energy system partition z; Represents the 0-1 partition state variable of node m, when When node m is divided into area z; represents the total power load of node m in partition z, represents the total power generation capacity of node m in partition z; represents the power consumption of the electric drive gas component of node m in partition z at rated output; represents the post-disaster failure rate of the unit at node m in partition z. When node m has no power generation resources, represents the post-disaster irreparable rate of the power load of node m in partition z; Natural gas resource abundance after a disaster Expressed as: in, represents the total gas load of node m in the integrated energy system partition z, represents the total gas supply capacity of node m in zone z of the integrated energy system; It represents the gas consumption of the gas-driven electric component of node m in partition z at rated output; represents the post-disaster failure rate of the gas supply device of node m in partition z. When node m has no gas supply resources, represents the post-disaster irreparable rate of the gas load of node m in partition z; System pipeline length indicator Expressed as: in, It is represented as the shortest geographical distance between node m and the black start unit in the integrated energy system zone z; System connectivity indicators Expressed as: in, represents the pipeline set of the integrated energy system partition z; It is represented as the 0-1 partition state variable of pipeline k, when When pipeline k is divided into area z; Post-disaster tower failure rate index Expressed as: in, It is expressed as the post-disaster tower failure rate of line k in the integrated energy system zone z; Post-disaster pipeline damage rate indicator Expressed as: in, It is expressed as the post-disaster pipeline damage rate of pipeline k in zone z of the integrated energy system; The objective function of the constructed objective function includes: the objective function of the electricity-gas integrated energy system partition optimization model is to minimize the comprehensive evaluation index R: Where, β 1 represents the weight of the power resource abundance index after a disaster, β 2 represents the weight of natural gas post-disaster resource abundance index, β 3 represents the weight of the system pipeline length indicator, β 4 System connectivity index weight, β 5 represents the weight of the tower failure rate index after the disaster, β 6 represents the weight of the pipeline damage rate indicator after the disaster; Z represents the total number of partitions of the electricity-gas integrated energy system; The constraints include: post-disaster zoning connectivity constraints for the electricity-gas integrated energy system, zoning classification constraints, transmission line tower failure constraints under typhoon weather, and gas pipeline damage constraints under earthquake disasters; The partition connectivity constraint constructed based on the state constraint virtual network flow method is expressed as: in, Indicates the black start unit in the zone z of the electric-gas integrated energy system, represents the set of pipelines connected to node m in the electricity-gas integrated energy system partition z; It is represented as the set of nodes in partition z connected to pipeline k; a 0-1 partition state variable representing the end of pipeline k in partition z, and represents the 0-1 partition state variable of the head-end node of pipeline k in partition z; represents the virtual network flow in pipeline k with the black start unit as the head node in partition z represents the virtual network flow in pipeline k with the black start unit as the terminal node in partition z; M is represented by a sufficiently large constant; The partition classification constraint based on the state constraint virtual network flow method is expressed as: in, Indicates the black start unit in the zone z of the electric-gas integrated energy system. Represents the gas supply node in the electricity-gas integrated energy system zone z, Represents the electric-driven gas components in the electric-gas integrated energy system zone z, N represents the set of gas generator sets in the electricity-gas integrated energy system partition z; min Indicates the minimum number of nodes allowed for each partition; The failure constraint of transmission line towers under typhoon weather is expressed as: in, It represents the actual wind speed that line k in zone z of the electricity-gas integrated energy system bears. represents the design wind speed that line k in zone z of the electricity-gas integrated energy system can withstand; α k,z represents the tower failure rate model coefficient of line k in zone z; The damage constraint of gas pipelines under earthquake disasters is expressed as: Among them, P k,z represents the total number of post-disaster damaged points of pipeline k in the electricity-gas integrated energy system zone z, κ k,z represents the total length of pipeline k in zone z of the electricity-gas integrated energy system; Establish an optimization model for the zoning of an electricity-gas integrated energy system based on internal and external operational risks, and solve the model to obtain the optimal post-disaster zoning results of the electricity-gas integrated energy system; According to the constructed objective function and constraints, a zoning optimization model of the electric-gas integrated energy system based on internal and external operation risks is constructed.

2. The method for optimizing the zoning of an integrated electricity-gas energy system according to claim 1, wherein: The energy supply risk indicators include the power resource abundance after the disaster R 1 and natural gas post-disaster resource abundance R 2 The pipeline network structure risk indicators include the system pipeline length R 3 and system connectivity R 4 , post-disaster environmental risk indicators include post-disaster tower failure rate R 5 and post-disaster pipeline damage rate R 6 .

3. The method for optimizing the zoning of an integrated electricity-gas energy system according to claim 2, wherein: The solution model obtains the optimal post-disaster zoning result of the electricity-gas integrated energy system, including inputting the optimization model and the topological structure and parameters of the actual electricity-gas integrated energy system into the solver to calculate the optimal post-disaster zoning result of the electricity-gas integrated energy system.

4. A system for optimizing the zoning of an electric-gas integrated energy system using the method according to any one of claims 1 to 3, characterized in that: Construct an evaluation index system module and build a multi-level evaluation index system based on the internal links and external influences of the electricity-gas integrated energy system; Calculation module, calculates evaluation indicators, and constructs objective functions and constraints; The solution module establishes an optimization model for the zoning of the electricity-gas integrated energy system based on internal and external operation risks, and solves the model to obtain the optimal post-disaster zoning result of the electricity-gas integrated energy system.

5. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 3 are implemented.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 3 are implemented.

Images