Partition optimization method and system for electricity-gas integrated energy system

Abstract

The invention discloses an electricity-gas integrated energy system partition optimization method and system, and the method comprises the steps: constructing a multi-level evaluation index system according to the internal links and external influences of an electricity-gas integrated energy system; and calculating an evaluation index, and constructing a target function and a constraint condition. And establishing an internal and external operation risk-based power-gas integrated energy system partitioning optimization model, and solving the model to obtain an optimal post-disaster partitioning result of the power-gas integrated energy system. According to the electricity-gas integrated energy system partition optimization method and system, the topological structure and the resource abundance of a gas transmission system are considered, internal and external risks possibly faced by the system in the post-disaster recovery period can also be considered, and a more reasonable and effective post-disaster partition strategy can be formulated; and parallel recovery of the electricity-gas comprehensive energy system in different areas can be realized, and stable energy supply can be maintained while the whole post-disaster recovery process is accelerated.

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 have become more frequent in recent years, including typhoons and earthquakes. The operational safety and stability of the electricity-gas integrated energy system are facing huge threats, and energy outage crises are prone to occur. Therefore, accelerating the post-disaster recovery process of the electricity-gas integrated energy system and minimizing losses are of great economic and social significance.

[0003] By constructing a reasonable post-disaster zoning strategy, it will be possible to achieve parallel recovery of multiple regions and accelerate the post-disaster recovery process of the electricity-gas integrated energy system. However, in reality, existing studies have only focused on how to construct a post-disaster zoning strategy for the power system, but lack sufficient attention to the optimization of the post-disaster zoning of the electricity-gas integrated energy system. In addition, with the increasingly close coupling between the power transmission and gas transmission systems, their scale and topological structure are larger and more complex. Therefore, how to balance the supply of electricity and gas resources and accelerate the recovery process after a disaster will be the key issue in constructing the optimal post-disaster zoning strategy for the electricity-gas integrated energy system. Summary of the invention

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

[0005] Therefore, the technical problem solved by the present invention is that existing studies have only focused on how to construct a post-disaster zoning strategy for the power system, but have not paid enough attention to the optimization of the post-disaster zoning of the electricity-gas integrated energy. In addition, as the coupling between the power transmission and gas transmission systems becomes increasingly close, their scale and topological structure 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 according to the internal links and external influences of the electricity-gas integrated energy system.

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

[0008] A zoning optimization model for the electricity-gas integrated energy system 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 solution of the method for optimizing the zoning of the electric-gas integrated energy system described in the present invention, 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 major main indicators for evaluation. After classification, each main indicator expands downward to a secondary indicator layer, and each secondary indicator layer contains two secondary indicators.

[0010] As a preferred solution of the method for optimizing the zoning 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 zoning of the electric-gas integrated energy system described in the present invention, the calculation evaluation index includes: the resource abundance degree after the power disaster It is 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 node m is divided into area z. represents the total power load of node m in partition z, represents the total 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, It represents the post-disaster irreparable rate of the power load of node m in partition z.

[0014] Natural gas resource abundance after a disaster It is 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. Represents 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 It is 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 index It is expressed as:

[0021]

[0022] in, Represents the set of pipelines in 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 It is expressed as:

[0024]

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

[0026] Post-disaster pipeline damage rate index It is 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 partition of the electric-gas integrated energy system of the present invention, wherein: the objective function of constructing the objective function includes that the objective function of the partition optimization model of the electric-gas integrated energy system is the minimum comprehensive evaluation index R:

[0030]

[0031] In the formula, β 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 post-disaster tower failure rate index, β 6 It represents the weight of pipeline damage rate after disaster. Z represents the total number of partitions of the electricity-gas integrated energy system.

[0032] As a preferred scheme of the electricity-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 electricity-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 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.

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

[0037]

[0038] in, It 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, It represents the electric-driven gas components in the zone z of the electric-gas integrated energy system. Represents the set of gas units in the electricity-gas integrated energy system zone z. 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 electric-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 pipeline under earthquake disaster is expressed as:

[0043]

[0044] Among them, P k,z represents the total number of post-disaster damaged points of pipeline k in the zone z of the electricity-gas integrated energy system, κ 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 scheme of the electricity-gas integrated energy system zoning optimization method described in the present invention, 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, and calculating the optimal post-disaster zoning result of the electricity-gas integrated energy system.

[0046] An electric-gas integrated energy system partition optimization system, characterized in that it includes:

[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 comprises 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 view of 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 operating 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, which 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] In order 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. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. Among them:

[0054] Figure 1 An overall flow chart of an electric-gas integrated energy system zoning optimization method and system provided for the first embodiment of the present invention.

[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, which is a method and system for optimizing the zoning of an electric-gas integrated energy system provided in a second embodiment of the present invention. DETAILED DESCRIPTION

[0056] In order to make the above-mentioned purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are described in detail below in conjunction with the drawings of the specification. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in the art without creative work should fall within the scope of protection of the present invention.

[0057] Example 1, reference Figure 1 , which is an 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 the hierarchical risk indicator system for integrated energy systems is established in accordance with the principles of scientificity, pertinence, independence, practicality and comparability. Based on the internal links and external influences of the electricity-gas integrated energy system, it is judged whether the internal energy supply links and energy transmission network structure of the integrated energy system have the risk of affecting the stable energy supply recovery of the system, the efficiency of post-disaster recovery, or even another shutdown during the post-disaster recovery period, and whether the continuous changes in the post-disaster environment will bring additional recovery risks to the key links of the integrated energy system.

[0060] The multi-level evaluation index system uses energy supply risk index, pipeline structure risk index, and post-disaster environmental risk index as three main indicators for evaluation. After classification, each main indicator expands downward to a secondary indicator layer, where each secondary indicator layer contains two secondary indicators. Table 1 shows the overall framework of the post-disaster risk index system for the electricity-gas integrated 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] Power resource adequacy after a disaster It is 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 node m is divided into area z. represents the total power load of node m in partition z, represents the total 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, It represents the post-disaster irreparable rate of the power load of node m in partition z.

[0067] Natural gas resource abundance after a disaster It is 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. Represents 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 It is 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 partition z.

[0073] System connectivity index It is expressed as:

[0074]

[0075] in, Represents the set of pipelines in 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 It is expressed as:

[0077]

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

[0079] Post-disaster pipeline damage rate index It is 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 partition of an electricity-gas integrated energy system according to claim 4, characterized in that: the objective function of constructing the objective function includes that the objective function of the partition optimization model of the electricity-gas integrated energy system is the minimum comprehensive evaluation index R:

[0083]

[0084] In the formula, β 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 post-disaster tower failure rate index, β 6 It represents the weight of pipeline damage rate after disaster. Z represents the total number of partitions of 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] In the formula, β 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 post-disaster tower failure rate index, β 6 It represents the weight of pipeline damage rate after disaster. Z represents the total number of partitions of the electricity-gas integrated energy system.

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

[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 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.

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

[0093]

[0094] in, It 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, It represents the electric-driven gas components in the zone z of the electric-gas integrated energy system. Represents the set of gas units in the electricity-gas integrated energy system zone z. 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 electric-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 pipeline under earthquake disaster is expressed as:

[0099]

[0100] Among them, P k,z represents the total number of post-disaster damaged points of pipeline k in the zone z of the electricity-gas integrated energy system, κ 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 operation 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 minimum comprehensive weighted value of the aforementioned power post-disaster resource abundance index, natural gas post-disaster resource abundance index, system pipeline length index, system connectivity index, post-disaster tower failure rate index and post-disaster pipeline damage rate index is used as the objective function, and 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 operation risks, and the model is solved by the Gurobi 10.0.0 solver of MATLAB software to obtain the optimal post-disaster partition result of 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 operation risks are as follows:

[0103] Step 1: Based on the internal links and external influences of the electricity-gas integrated energy system, a multi-level evaluation index 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 resource abundance index for electricity, the post-disaster resource abundance index for natural gas, the system pipeline length index, the system connectivity index, the post-disaster tower failure rate index, and the post-disaster pipeline damage rate index, and assign weights to these indicators.

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

[0106] Step 4: Take the indicators and their weights constructed in step 1 as the objective function, and then take 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 topological structure and parameters of the actual electric-gas integrated energy system into the MATLAB software.

[0108] Step 6: Solve the model using Gurobi 10.0.0 solver in MATLAB software, and calculate the optimal post-disaster partition 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 for short) 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 a 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] A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment method can be completed 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 the memory, database or other medium used in the embodiments provided in the present application can 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. As an 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 each embodiment provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include distributed databases based on blockchains, etc., but are not limited to this. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but are not limited to this.

[0114] Example 2, reference Figure 2 , which is an embodiment of the present invention, provides a method and system for optimizing the zoning 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] The present invention uses an actual 125-node power system (represented as black lines) and a 49-node natural gas system (represented as red lines) in a certain region of China for example analysis. The topological structure is shown in Figure 2 In addition, in this implementation example, 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 the present 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 the zones, and another 4 lines cannot be included in the recovery zone due to the continuous typhoon disaster after the shutdown (shown as black dotted lines); node 121 cannot be included in any recovery zone due to the typhoon disaster; no gas pipeline is not included in the recovery zone due to the earthquake disaster; the average distance between each node and the black start unit is 25.47km; the absolute difference between the power generation capacity and the power load is 0.38pu; the absolute difference between the natural gas supply and the gas demand is 7840m3; the wind speed borne by each transmission line is 19.63m / s.

[0117] Likewise Figure 2 As shown, according to the number of black start units, the electric-gas integrated energy system is divided into 5 partitions in total. More specifically, partition 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 partition 1, respectively, while 2 gas sources and 1 compressor are located at nodes 132, 139 and 133 of the gas supply system, respectively; partition 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 wind turbine unit 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 batteries are located at nodes 76 and 84 respectively; in zone 4, one black start 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 black start unit and three non-black start units are located at nodes 30, 34, 42 and 110 respectively, while one gas source is located at node 153. It should be noted that the three black start units in zones 1, 4 and 5 are all gas units, and the three non-black start units on nodes 42, 44 and 110 are also gas units. 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 of these 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 stages of the electric-gas integrated energy system. Therefore, most natural gas pipelines and nodes should be allocated to partitions with gas black start units to provide a stable and sufficient supply of natural gas. Some electrically driven gas components, such as P2G devices and electrically 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 partition strategy that takes into account the respective topological structures and resource abundance 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 rather than to limit it. 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 electric-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; Calculate evaluation indicators, construct objective functions and constraints; A zoning optimization model for the electricity-gas integrated energy system 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.

2. The method for optimizing the zoning of an electric-gas integrated energy system according to claim 1, characterized in that: 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 indicators for evaluation; After classification, each main indicator expands downward to a secondary indicator layer, and each secondary indicator layer contains two secondary indicators.

3. The method for optimizing the zoning of an electric-gas integrated energy system according to claim 2, characterized in that: 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 .

4. The method for optimizing the zoning of an electric-gas integrated energy system according to claim 3, characterized in that: The calculation evaluation index includes: the power resource sufficiency after the disaster It is 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 node m is divided into region 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 It is 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 element 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 It is 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 index It is expressed as: in, represents the set of pipelines in 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 It is expressed as: in, It is expressed as the post-disaster tower failure rate of line k in zone z of the integrated energy system; Post-disaster pipeline damage rate index It is expressed as: in, It is expressed as the post-disaster pipeline damage rate of pipeline k in zone z of the integrated energy system.

5. The method for optimizing the zoning of an electric-gas integrated energy system according to claim 4, characterized in that: The objective function of constructing the objective function includes that the objective function of the electricity-gas integrated energy system partition optimization model is to minimize the comprehensive evaluation index R: In the formula, β 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 post-disaster tower failure rate index, β 6 represents the weight of the pipeline damage rate indicator after the disaster; Z represents the total number of partitions in the electricity-gas integrated energy system.

6. The method for optimizing the zoning of an electric-gas integrated energy system according to claim 5, characterized in that: The constraints include: constituting post-disaster zoning zoning of the electricity-gas integrated energy system, zoning classification constraints, power 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 zone z of the electricity-gas integrated energy system; 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, It 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, It represents the electric-driven gas components in the zone z of the electric-gas integrated energy system. represents the set of gas units in the electricity-gas integrated energy system zone z; N 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 electric-gas integrated energy system bears. represents the design wind speed that line k in zone z of the electricity-gas integrated energy system bears; α k,z represents the tower failure rate model coefficient of line k in zone z; The damage constraint of gas pipeline under earthquake disaster is expressed as: Among them, P k,z represents the total number of post-disaster damaged points of pipeline k in the zone z of the electricity-gas integrated energy system, κ k,z It represents the number and total length of pipeline k in zone z of the electricity-gas integrated energy system.

7. The method for optimizing the zoning of an electric-gas integrated energy system according to claim 6, characterized in that: 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, and calculating the optimal post-disaster zoning result of the electricity-gas integrated energy system.

8. A system for optimizing the zoning of an electric-gas integrated energy system using the method according to any one of claims 1 to 7, 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.

9. 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 7 are implemented.

10. 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 7 are implemented.

Images