A dynamic configuration method for electric-pneumatic coupling system based on digital twin technology

Abstract

The present invention discloses a method for dynamic configuration of an electric-gas coupling system based on digital twin technology. The method comprises: establishing a digital twin system of the electric-gas coupling system; obtaining and updating device parameters; a configuration requirement assessment model outputting a comprehensive evaluation value of the configuration requirement, comparing until the value is less than an evaluation threshold, and placing the device model library into a problem library of a configuration subsystem; a dynamic configuration model of the electric-gas coupling system outputting construction data of the electric-gas coupling system, adjusting the current device model library based on the construction data; until the obtained comprehensive evaluation value of the configuration requirement is greater than the evaluation threshold, placing the current device model library into a project library of the configuration subsystem; and reconfiguring the electric-gas coupling system to achieve dynamic configuration of the electric-gas coupling system. Based on digital twin technology, the method of the present invention can improve the robustness of the configuration scheme in dealing with uncertainty, significantly reduce redundancy, and facilitate the absorption of new energy.

Description

Technical Field

[0001] The present invention relates to a dynamic configuration method for an electric-pneumatic coupling system, and in particular to a dynamic configuration method for an electric-pneumatic coupling system based on digital twin technology. Background Art

[0002] With the widespread deployment of natural gas generators, the coupling between power and gas systems is deepening. Increased uncertainty factors, such as renewable energy output and multi-energy loads, lead to increased redundancy in configuration to address these uncertainties. Therefore, a new configuration method for coupled power-gas systems is urgently needed to achieve optimal decision-making during the configuration phase.

[0003] Existing configuration methods often characterize uncertainty based on optimal power flow, using the "Nk check" to verify whether the system meets reliability requirements under typical scenarios. The uncertainty of renewable energy output and multi-energy loads significantly increases the uncertainty of the configuration boundaries, making existing configuration methods prone to suboptimal configuration solutions. Existing configurations are typically open-loop, long-term configurations or fixed-year rolling configurations. Due to a lack of data accumulation, long-term configurations often have high redundancy to account for uncertainty in the distant future. Rolling configurations often have rigid fixed-year configurations and fail to reflect the system's real-time state and actual configuration requirements. Existing configurations typically use factory values ​​for device parameters as configuration inputs. However, device parameters can change during operation due to factors such as aging and environmental degradation. Existing configurations cannot accurately capture these changes, resulting in suboptimal configuration results. Therefore, new configuration methods for coupled electrical and pneumatic systems are urgently needed to effectively address the increased redundancy caused by uncertainty and reflect the actual system state and configuration requirements. Summary of the Invention

[0004] To address the issues presented in the background art, the present invention provides a dynamic configuration method for an electric-gas coupling system based on digital twin technology. By modifying and dynamically evaluating device parameters, the method enables more timely and effective assessment of the urgency of electric-gas coupling system configuration, accurately determining configuration requirements. Through a problem library and project library, it effectively addresses the uncertainty of new energy sources in a vast array of scenarios, reduces the redundancy of configuration results, and promotes renewable energy consumption.

[0005] The technical solution adopted in the present invention is:

[0006] The method for dynamic configuration of an electric-pneumatic coupling system of the present invention comprises the following steps:

[0007] Step 1) Establish a digital twin system of the electro-pneumatic coupling system, including a device model subsystem, an evaluation subsystem, and a configuration subsystem; the device model subsystem includes a device model library, which includes device models of various devices in the electro-pneumatic coupling system; the evaluation subsystem includes a configuration requirement evaluation model of the electro-pneumatic coupling system; the configuration subsystem includes an electro-pneumatic coupling system dynamic configuration model, a problem library, and a project library.

[0008] A digital twin system is a digital space that can represent the true state of an electro-pneumatic coupling system, enabling the mapping of the electro-pneumatic coupling system to a digital twin system. This helps guide the configuration design and actual production of the electro-pneumatic coupling system from a full lifecycle perspective. Each device model in the device model library of the device model subsystem of the digital twin system is a mapping of each device in the electro-pneumatic coupling system in the digital space. The configuration requirement assessment model of the evaluation subsystem of the digital twin system can truly evaluate the configuration requirements of the electro-pneumatic coupling system. The dynamic configuration model of the configuration subsystem of the digital twin system can provide configuration solutions. Initially, no scenarios exist in the problem library and project library of the configuration subsystem, i.e., they are empty sets. During multiple dynamic configuration processes, scenarios that are inferior to the evaluation threshold are placed in the problem library, forming a problem library of the configuration subsystem that includes a set of scenarios that are inferior to the evaluation threshold. The dynamic configuration model then provides corresponding configuration solutions for scenarios that are inferior to the evaluation threshold. These configuration solutions are scenarios that are superior to the evaluation threshold. Scenarios that are superior to the evaluation threshold are placed in the project library, forming a project library of the configuration subsystem that includes a set of scenarios that are superior to the evaluation threshold.

[0009] Step 2) for each device in the electro-pneumatic coupling system and its device model in the device model library of the device model subsystem, the device model includes several device parameters of the device; after the electro-pneumatic coupling system runs for a preset time period, the observed value of each device parameter of each device at the current moment is obtained, and for each observed value of each device parameter at the current moment, the device parameter is corrected using the least squares method to obtain an estimated value of the device parameter and update the device parameter in the device model; state perception is achieved through real-time information interaction to obtain data, and a true mapping of the electro-pneumatic coupling system is obtained after correction.

[0010] At the initial moment, every device parameter of every device in the electro-pneumatic coupling system is an initial value. When establishing the device model subsystem of the digital twin system, every device parameter of every device in the device model subsystem is set to an initial value. When the electro-pneumatic coupling system starts running, the device parameters will deviate from the initial values ​​due to factors such as environmental changes and device aging. The electro-pneumatic coupling system can be dynamically configured. Through information exchange between the digital twin system and the electro-pneumatic coupling system, the operating data of the current electro-pneumatic coupling system is collected during multiple repeated configurations. The least squares method is used to correct the device parameters in the device model library during each repeated configuration, thereby realizing the digital twin system's state perception of the electro-pneumatic coupling system and ensuring that the digital twin system is a mapping of the real electro-pneumatic coupling system.

[0011] Information exchange between the electro-pneumatic coupling system and the digital twin system is achieved through sensors. Information exchange includes data acquisition and state perception. Data acquisition refers to collecting data from various devices in the electro-pneumatic coupling system and sending it to the digital twin system. State perception refers to the digital twin system updating the parameters and other information of the digital space based on the collected data to ensure the accuracy of the mapping.

[0012] Step 3) The estimated values ​​of the various device parameters in the device model of the device model library of the device model subsystem obtained in step 2) are used as the device parameters of the various devices in the electric-pneumatic coupling system in the current operating state. Several operating data of the electric-pneumatic coupling system at this time are obtained and input into the configuration requirement evaluation model of the electric-pneumatic coupling system of the evaluation subsystem. The configuration requirement evaluation model outputs a comprehensive evaluation value of the configuration requirements of the electric-pneumatic coupling system. The comprehensive evaluation value of the configuration requirements is compared with the preset evaluation threshold for triggering the configuration. When the comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, steps 2)-3) are repeated until the comprehensive evaluation value of the configuration requirements is less than the evaluation threshold. The device model library at this time and the comprehensive evaluation value of the configuration requirements calculated this time are then placed into the problem library of the configuration subsystem.

[0013] The preset evaluation threshold combines the development status and goals of the electric-gas coupling system. The role of the evaluation threshold is to compare the comprehensive evaluation value of the configuration requirements calculated by the evaluation subsystem of the digital twin system each time with the evaluation threshold to decide whether to trigger the configuration, that is, whether to continue with step 4).

[0014] Step 4) Obtain some operating data of the electric-gas coupling system at this time and input it into the dynamic configuration model of the electric-gas coupling system of the configuration subsystem of the digital twin system. The dynamic configuration model of the electric-gas coupling system outputs the construction data of the electric-gas coupling system at this time, and adjusts the equipment model library of the current equipment model subsystem according to the construction data.

[0015] Step 5) obtains each device model in the device model library of the device model subsystem adjusted in step 4) and performs the same operations on the device models in steps 2)-3). If the obtained comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, the current device model library is placed in the project library of the configuration subsystem; if the obtained comprehensive evaluation value of the configuration requirements is less than the evaluation threshold, repeat steps 4)-5) until the obtained comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, and the current device model library is placed in the project library of the configuration subsystem.

[0016] Each device model library and the comprehensive evaluation value of configuration requirements in the problem library correspond to each device model library in the project library. The problem library and the project library can be used as historical data sets. When the device model library of a certain electric-gas coupling system corresponds to one of the device model libraries in the problem library, the configuration of the electric-gas coupling system can be directly adjusted according to the device model library in the project library corresponding to the device model library in the problem library.

[0017] Step 6) reconfiguring the electric-pneumatic coupling system according to the current device model library put into the project library in step 5) to achieve dynamic configuration of the electric-pneumatic coupling system.

[0018] In the step 1), the electric-gas coupling system includes an electric power network and a natural gas network; the equipment in the electric power network includes several power stations, transmission lines, electric load equipment and energy storage equipment, the power station includes several gas generators and non-gas generators, each gas generator, non-gas generator, electric load equipment and energy storage equipment are interconnected through various transmission lines, the non-gas generator is specifically a coal-fired generator or a new energy generator; the electric load equipment is specifically an electric-to-gas device or an electric load equipment; the energy storage equipment is specifically an energy storage battery; the equipment in the natural gas network includes several gas production stations, gas transmission pipelines and gas load equipment, each gas production station and gas load equipment The equipment is connected through various gas pipelines; the gas production station includes several gas sources and power-to-gas equipment; the gas load equipment is specifically gas-to-electricity load equipment or gas load equipment, and the gas-to-electricity load equipment is specifically a gas generator; the power-to-gas equipment of the power network is connected to the nodes of the natural gas network, and the gas-to-electricity equipment of the natural gas network is connected to the nodes of the power network; the equipment in the power-gas coupling system includes gas generators, non-gas generators, electric load equipment, transmission lines and energy storage equipment in the power network, and gas sources, power-to-gas equipment, gas pipelines and gas load equipment in the natural gas network. Except for the transmission lines and gas pipelines, each equipment in the power-gas coupling system is located at its respective node.

[0019] In step 2), after the electric-gas coupling system runs for a preset period of time, the observed value of each device parameter of each device at the current moment is obtained. For each observed value of the device parameter at the current moment, the device parameter is corrected using the least squares method to obtain an estimated value of the device parameter, as follows:

[0020]

[0021] Among them, x i,k+1 represents the observed value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time. Represents the estimated value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time.

[0022] The estimated value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time is used Update the i-th device parameter in the device model library of the device model subsystem at the current k+1 time.

[0023] Device parameters in the electric-gas coupled system include the voltage amplitude / phase angle of power nodes, active / reactive power transmitted by transmission lines, rated power of generators, gas pressure at natural gas nodes, transmission capacity of natural gas pipelines, conductance and susceptance of electric lines, and transmission coefficients of natural gas pipelines. Least squares methods are used to process this large amount of observed device parameter data to restore their true values. These true values ​​are then used to correct the corresponding device parameters in the device library model. By dynamically updating device parameters, the digital twin system can perceive the state of the electric-gas coupled system, ensuring that the digital twin system is a reflection of the actual electric-gas coupled system.

[0024] In step 3), the configuration requirement evaluation model of the electric-pneumatic coupling system of the evaluation subsystem is specifically as follows:

[0025] Ω=f(EENS)+f(EGNS)+f(β car )+f(η re )+f(C inv )+f(C ope )

[0026] Where Ω represents the comprehensive evaluation value of configuration requirements, f() represents the normalization formula, EENS represents the expected value of electricity shortage in the electricity-gas coupling system, EGNS represents the expected value of natural gas shortage in the electricity-gas coupling system, and β car represents the carbon emissions of the electric-gas coupling system, η re represents the renewable energy consumption rate of the electricity-gas coupling system, C inv Indicates the construction measurement value of the electric-gas coupling system, C opeIndicates the operating measurement value of the electric-pneumatic coupling system.

[0027] The normalization processing formula is as follows:

[0028]

[0029] Where y represents the indicator to be processed in the normalized processing formula f(), which includes positive indicators and negative indicators. When y is a positive indicator, the positive indicator y includes the renewable energy consumption rate η of the electricity-gas coupling system. re , the larger the value of the positive index y, the higher the comprehensive evaluation value of the configuration demand Ω; when y is a negative index, the negative index y includes the expected value of electricity shortage EENS of the electric-gas coupling system, the expected value of natural gas shortage EGNS, and the carbon emission β car , Construction measurement value C inv and operating measurement value C ope , the smaller the value of the reverse index y, the higher the comprehensive evaluation value of the configuration demand and the higher Ω.

[0030] The expected electricity shortage value EENS, the expected natural gas shortage value EGNS, and the carbon emission value β of the electric-gas coupling system are car , renewable energy consumption rate η re , Construction measurement value C of electric-gas coupling system inv and the operating measurement value C of the electric-pneumatic coupling system ope , as follows:

[0031]

[0032] Wherein, T represents the time period between the current moment and the initial moment of operation of the electric-gas coupling system or the moment of the last repetition of step 2); the T period is divided into several time periods, and Ξs represents the set of the time periods divided by the T period; p s represents the probability of time period s within time period T, p s =s / T; It represents the power loss load of the electric-gas coupling system in the period s; represents the natural gas load loss of the electricity-gas coupling system in the period s; N represents the number of power stations in the power network of the electricity-gas coupling system, represents the carbon emission factor of the g-th power station among all power stations in the power network, represents the power generation of the gth power station in the power grid; M represents the number of gas production stations in the natural gas network of the power-gas coupling system, represents the carbon emission factor of the jth gas production station in the natural gas network, represents the gas production of the jth gas production station in the natural gas network; Q H Represents the total power generation of each power station in the power network; Q L Indicates the total gas production of each gas production station in the natural gas network; Q re represents the total renewable energy generation of each new energy generator in the power network; Q e Indicates the total amount of electricity transmitted by the electric-gas coupling system; t represents the net present value of the electric-gas coupling system at time t, o t =1 / (1+d) -t , d represents the discount rate of the electric-gas coupling system at time t; ψ (·) represents the residual value rate of fixed assets of the equipment to be built in the electric-gas coupling system; C (·) Indicates the construction measurement value of the equipment to be built in the electric-gas coupling system; (·),t and z (·),t-1 Respectively represent the status of the equipment to be built in the electric-gas coupling system at time t and time t-1. If the status of the equipment to be built is 1, it is 0 if it does not exist; z g,t Indicates the status of the g-th power station in the power grid at time t, which is 1 if it exists and 0 if it does not exist; C gen represents the unit capacity fuel measurement value of the g-th power station among the power stations in the power network; g,t represents the output power of the gth power station among the power stations in the power network at time t; z a,t Indicates the state of the ath gas source in the natural gas network, which is 1 if it exists and 0 if it does not exist; C well Indicates the unit capacity fuel measurement value of the ath gas source among all gas sources in the natural gas network; G a,t It represents the gas flow of the ath gas source among all gas sources in the natural gas network at time t.

[0033] The operating data of the electric-gas coupling system obtained in step 3) at this time includes the power loss load of the electric-gas coupling system in the time period s Natural gas load loss The power generation of the gth power station among all power stations in the power network The gas production of the jth gas production station in the natural gas network The total power generation Q of each power station in the power network H , the total gas production Q of each gas production station in the natural gas network L , the total renewable energy generation capacity Q of each new energy generator in the power network re And the total amount of electricity transmitted by the electric-gas coupling system Q e .

[0034] The probability of failure of the power-gas coupling system under uncertainty is characterized by the expected value of electricity shortage EENS and the expected value of natural gas shortage EGNS; the probability of failure of the power-gas coupling system under uncertainty is characterized by the expected value of carbon emissions β car and renewable energy consumption rate η re etc. to characterize the possible impact of the electric-gas coupling system on the environment; according to the construction measurement value C inv and operating measurement value C ope To characterize the cost of construction and operation of the electric-gas coupling system.

[0035] In step 4), the dynamic configuration model of the electric-pneumatic coupling system of the configuration subsystem of the digital twin system is specifically the electric-pneumatic coupling system objective function considering the constraints of the electric-pneumatic coupling system. The electric-pneumatic coupling system objective function is specifically as follows:

[0036]

[0037] in, The shedding measurement value representing the unit load shedding amount of the power network at time t; The EENS is the load shedding value of the natural gas network at time t. t represents the expected load shedding of the power network at time t; EGNS t represents the expected load shedding of the natural gas network at time t.

[0038] The several operating data of the electric-gas coupling system obtained in step 4) at this time include the unit capacity fuel metering value C of the g-th power station among the various power stations in the power network. gen , the unit capacity fuel measurement value C of the ath gas source among all gas sources in the natural gas network well , the cut-off value of the unit load shedding of the power network at time t And the removal measurement value of the unit load shedding of the natural gas network at time t The construction data of the electric-gas coupling system output by the dynamic configuration model of the electric-gas coupling system includes the state z of the equipment to be built in the electric-gas coupling system at time t. (·),t , the output power P of the gth power station among the power stations in the power network at time t g,t , the gas flow rate G of the ath gas source among the gas sources in the natural gas network at time t a,t , the expected load shedding EENS of the power network at time t t And the expected load shedding of the natural gas network at time t EGNS t .

[0039] The constraints of the electric-gas coupling system are as follows:

[0040] z(·),t-1 ≤z (·),t

[0041] EENS t ≤EENS set

[0042] EGNS t ≤EGNS set

[0043]

[0044] Among them, EENS set Indicates the expected value of the power network's preset power supply shortage; EGNS set represents the expected value of the preset natural gas supply shortage of the natural gas network; B r represents the power network constant term, B r is a very large number; z r,t Indicates the state of the r-th transmission line in the power network at time t, which is 1 if it exists and 0 if it does not exist. The r-th transmission line in the power network is located between the u-th node and the v-th node in the power network. uv,t P represents the active power transmitted from the uth node to the vth node on the rth transmission line in the power network at time t, vu,t represents the active power transmitted from the vth node to the uth node on the rth transmission line in the power network at time t; θ u,t represents the power angle of the u-th node in the power network at time t, θ v,t represents the power angle of the vth node in the power network at time t; w r represents the reactance of the rth transmission line in the power network; each gas generator, non-gas generator and energy storage device in the power network constitutes a first device set, the first device set includes a plurality of first devices, and They represent the minimum power and maximum power of the qth first device in the first device set, P q,t represents the output power of the qth first device in the first device set at time t, z q,t Indicates the status of the qth first device in the power network, which is 1 if it exists and 0 if it does not exist; ur K represents the element in row u and column r of the correlation matrix composed of each node and each transmission line in the power network. uq represents the element in the uth row and qth column of the association matrix formed by each node in the power network and each first device in the first device set, W uε Represents the element in the uth row and the εth column of the association matrix composed of each node and each electric load device in the power network, represents the normal load level of the εth load device in the power network at time t; V uδ Represents the element in the uth row and the δth column of the association matrix formed by each node in the power network and each power-to-gas device in the natural gas network; represents the power-to-gas load level of the δth power-to-gas device in the natural gas network at time t; z p,t Indicates the state of the pth gas pipeline in the electric-gas coupling system, which is 1 if it exists and 0 if it does not exist; C p represents the transmission coefficient of the pth gas pipeline in the electric-gas coupling system; and is a 0-1 variable. The pth gas transmission pipeline in the electric-gas coupling system connects the mth node and the nth node of the natural gas network. When the natural gas in the pth gas transmission pipeline flows out from the mth node, When the natural gas in the pth gas pipeline flows into the mth node, π m,t and π n,t They represent the square of the gas pressure at the mth node and the nth node in the natural gas network at time t; Z p,t represents the gas flow rate of the pth gas transmission pipeline in the electric-gas coupling system at time t; and They represent the maximum and minimum gas pressures of the mth node in the natural gas network respectively; and They represent the maximum transmission capacity and minimum transmission capacity of the pth gas pipeline in the electric-gas coupling system; z j,t Indicates the status of the j-th gas production station in the natural gas network, which is 1 if it exists and 0 if it does not exist. mj Represents the element in the mth row and jth column of the association matrix composed of each node and each gas production station in the natural gas network, G j,t represents the gas production of the jth gas production station in the natural gas network, z δ,t Indicates the status of the δth power-to-gas device in the natural gas network, which is 1 if it exists and 0 if it does not exist. mδ Represents the element in the mth row and the δth column of the association matrix composed of each node and each power-to-gas device in the natural gas network, represents the gas production of the δth power-to-gas device in the natural gas network at time t; Y mp represents the element in the mth row and pth column of the connection matrix composed of each node and each gas pipeline in the natural gas network; J mn Represents the element in the mth row and nth column of the connection matrix composed of each node and each gas load device in the natural gas network, It represents the gas load level of the nth gas load device in the natural gas network at time t, U mλRepresents the element in the mth row and λth column of the connection matrix composed of each node and each gas-to-electricity load device in the natural gas network, represents the gas-to-electricity load level of the λth gas-to-electricity load device in the natural gas network at time t; μ represents the conversion coefficient of electricity to natural gas; Γ represents the conversion coefficient of gas volume to electrical power; and η represents the energy conversion efficiency of the power-to-gas device in the natural gas network.

[0045] The several operating data of the electric-gas coupling system obtained in step 4) at this time also include the expected power shortage value EENS of the power network set , the expected gas supply shortage value EGNS of the natural gas network set , the conventional electric load level of the εth electric load device in the power network at time t The gas load level of the nth gas load device in the natural gas network at time t The construction data of the electric-gas coupling system output by the dynamic configuration model of the electric-gas coupling system at this time also includes the active power P transmitted from the uth node to the vth node on the rth transmission line in the power network at time t uv,t , the active power P transmitted from the vth node to the uth node on the rth transmission line in each transmission line in the power network at time t vu,t , the output power P of the qth first device in the first device set at time t q,t , the power-to-gas load level of the δth power-to-gas device in the natural gas network at time t The square of the gas pressure at the mth node and the nth node in the natural gas network at time t is π m,b,t and π n,b,t , the gas flow Z of the pth gas pipeline in the electric-gas coupling system at time t p,b,t , the gas production of the δth power-to-gas device in the natural gas network at time t And the gas-to-electricity load level of the λth gas-to-electricity load device in the natural gas network at time t

[0046] The device model library of the current device model subsystem is adjusted according to the construction data. Specifically, the device model of the device to be built with the existing status is added to the device model library of the device model subsystem according to the construction data, and the output power of each power station in the power network, the gas flow of each gas source in the natural gas network, the expected load shedding amount of the power network and the expected load shedding amount of the natural gas network, the active power on each transmission line in the power network, the output power of each first device in the first device set, the power-to-gas load level of each power-to-gas device in the natural gas network, the square of the gas pressure of each node in the natural gas network, the gas flow of each gas transmission pipeline in the electric-gas coupling system, the gas production of each power-to-gas device in the natural gas network, and the gas-to-power load level of each gas-to-power load device in the natural gas network are adjusted according to the output construction data.

[0047] In the step 5), after repeating steps 4)-5), the final construction data of the electric-gas coupling system is obtained, that is, the construction data of the current device model library, and the current device model library is placed in the project library of the configuration subsystem.

[0048] In step 6), the electric-gas coupling system is reconfigured according to the current device model library placed in the project library in step 5), that is, the electric-gas coupling system is constructed and adjusted according to the final construction data to achieve dynamic configuration of the electric-gas coupling system.

[0049] The beneficial effects of the present invention are:

[0050] The method of the present invention is a new method for realizing dynamic configuration of electric-gas coupling system under large-scale grid connection of new energy, overcoming the deficiency of previous configuration methods in dealing with uncertainty and high redundancy; the method of the present invention can accurately evaluate the configuration requirements of electric-gas coupling system based on digital twin technology, solve the problem that the configuration scheme easily loses optimality due to the uncertainty of new energy, ensure the optimal configuration of equipment, reduce the redundancy of configuration results, and improve the absorption rate of new energy. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a logic block diagram of the method of the present invention;

[0052] Figure 2 Schematic diagram of the electric-gas coupling system and the digital twin system in the present invention. DETAILED DESCRIPTION

[0053] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0054] like Figure 1 As shown, the method for dynamic configuration of the electric-pneumatic coupling system of the present invention includes the following steps:

[0055] Step 1) Figure 2 As shown, a digital twin system of the electric-pneumatic coupling system is established, including a device model subsystem, an evaluation subsystem and a configuration subsystem; the device model subsystem includes a device model library, which includes device models of various devices in the electric-pneumatic coupling system; the evaluation subsystem includes a configuration requirement evaluation model of the electric-pneumatic coupling system; the configuration subsystem includes an electric-pneumatic coupling system dynamic configuration model, a problem library and a project library of the electric-pneumatic coupling system.

[0056] In step 1), in the step 1), the electric-gas coupling system includes an electric power network and a natural gas network; the equipment in the electric power network includes several power stations, transmission lines, electric load equipment and energy storage equipment, the power station includes several gas generators and non-gas generators, each gas generator, non-gas generator, electric load equipment and energy storage equipment are interconnected through various transmission lines, the non-gas generator is specifically a coal-fired generator or a new energy generator; the electric load equipment is specifically an electric-to-gas equipment or an electric load equipment; the energy storage equipment is specifically an energy storage battery; the equipment in the natural gas network includes several gas production stations, gas transmission pipelines and gas load equipment, each gas production station and gas load equipment They are connected through various gas pipelines; the gas production station includes several gas sources and power-to-gas equipment; the gas load equipment is specifically gas-to-electricity load equipment or gas load equipment, and the gas-to-electricity load equipment is specifically a gas generator; the power-to-gas equipment of the power network is connected to the nodes of the natural gas network, and the gas-to-electricity equipment of the natural gas network is connected to the nodes of the power network; the equipment in the power-gas coupling system includes gas generators, non-gas generators, electric load equipment, transmission lines and energy storage equipment in the power network, and gas sources, power-to-gas equipment, gas pipelines and gas load equipment in the natural gas network. Except for the transmission lines and gas pipelines, each equipment in the power-gas coupling system is located at its respective node.

[0057] A digital twin system is a digital space that can represent the true state of an electro-pneumatic coupling system, enabling the mapping of the electro-pneumatic coupling system to a digital twin system. This helps guide the configuration design and actual production of the electro-pneumatic coupling system from a full lifecycle perspective. Each device model in the device model library of the device model subsystem of the digital twin system is a mapping of each device in the electro-pneumatic coupling system in the digital space. The configuration requirement assessment model of the evaluation subsystem of the digital twin system can truly evaluate the configuration requirements of the electro-pneumatic coupling system. The dynamic configuration model of the configuration subsystem of the digital twin system can provide configuration solutions. Initially, no scenarios exist in the problem library and project library of the configuration subsystem, i.e., they are empty sets. During multiple dynamic configuration processes, scenarios that are inferior to the evaluation threshold are placed in the problem library, forming a problem library of the configuration subsystem that includes a set of scenarios that are inferior to the evaluation threshold. The dynamic configuration model then provides corresponding configuration solutions for scenarios that are inferior to the evaluation threshold. These configuration solutions are scenarios that are superior to the evaluation threshold. Scenarios that are superior to the evaluation threshold are placed in the project library, forming a project library of the configuration subsystem that includes a set of scenarios that are superior to the evaluation threshold.

[0058] Step 2) for each device in the electro-pneumatic coupling system and its device model in the device model library of the device model subsystem, the device model includes several device parameters of the device; after the electro-pneumatic coupling system runs for a preset time period, the observed value of each device parameter of each device at the current moment is obtained, and for each observed value of each device parameter at the current moment, the device parameter is corrected using the least squares method to obtain an estimated value of the device parameter and update the device parameter in the device model; state perception is achieved through real-time information interaction to collect data, and a true mapping of the electro-pneumatic coupling system is obtained after correction.

[0059] In step 2), after the electric-gas coupling system runs for a preset period of time, the observed value of each device parameter of each device at the current moment is obtained. For the observed value of each device parameter at the current moment, the device parameter is corrected using the least squares method to obtain an estimated value of the device parameter, as follows:

[0060]

[0061] Among them, x i,k+1 represents the observed value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time. Represents the estimated value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time.

[0062] The estimated value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time is used Update the i-th device parameter in the device model library of the device model subsystem at the current k+1 time.

[0063] Device parameters in the electric-gas coupled system include the voltage amplitude / phase angle of power nodes, active / reactive power transmitted by transmission lines, rated power of generators, gas pressure at natural gas nodes, transmission capacity of natural gas pipelines, conductance and susceptance of electric lines, and transmission coefficients of natural gas pipelines. Least squares methods are used to process this large amount of observed device parameter data to restore their true values. These true values ​​are then used to correct the corresponding device parameters in the device library model. By dynamically updating device parameters, the digital twin system can perceive the state of the electric-gas coupled system, ensuring that the digital twin system is a reflection of the actual electric-gas coupled system.

[0064] At the initial moment, every device parameter of every device in the electro-pneumatic coupling system is an initial value. When establishing the device model subsystem of the digital twin system, every device parameter of every device in the device model subsystem is set to an initial value. When the electro-pneumatic coupling system starts running, the device parameters will deviate from the initial values ​​due to factors such as environmental changes and device aging. The electro-pneumatic coupling system can be dynamically configured. Through information exchange between the digital twin system and the electro-pneumatic coupling system, the operating data of the current electro-pneumatic coupling system is collected during multiple repeated configurations. The least squares method is used to correct the device parameters in the device model library during each repeated configuration, thereby realizing the digital twin system's state perception of the electro-pneumatic coupling system and ensuring that the digital twin system is a mapping of the real electro-pneumatic coupling system.

[0065] Information exchange between the electro-pneumatic coupling system and the digital twin system is achieved through sensors. Information exchange includes data acquisition and state perception. Data acquisition refers to collecting data from various devices in the electro-pneumatic coupling system and sending it to the digital twin system. State perception refers to the digital twin system updating the parameters and other information of the digital space based on the collected data to ensure the accuracy of the mapping.

[0066] Step 3) The estimated values ​​of the various device parameters in the device model of the device model library of the device model subsystem obtained in step 2) are used as the device parameters of the various devices in the electric-pneumatic coupling system in the current operating state. Several operating data of the electric-pneumatic coupling system at this time are obtained and input into the configuration requirement evaluation model of the electric-pneumatic coupling system of the evaluation subsystem. The configuration requirement evaluation model outputs a comprehensive evaluation value of the configuration requirements of the electric-pneumatic coupling system. The comprehensive evaluation value of the configuration requirements is compared with the preset evaluation threshold for triggering the configuration. When the comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, steps 2)-3) are repeated until the comprehensive evaluation value of the configuration requirements is less than the evaluation threshold. The device model library at this time and the comprehensive evaluation value of the configuration requirements calculated this time are then placed into the problem library of the configuration subsystem.

[0067] In step 3), the configuration requirement evaluation model of the electric-pneumatic coupling system of the evaluation subsystem is as follows:

[0068] Ω=f(EENS)+f(EGNS)+f(β car )+f(η re )+f(C inv )+f(C ope )

[0069] Where Ω represents the comprehensive evaluation value of configuration requirements, f() represents the normalization formula, EENS represents the expected value of electricity shortage in the electricity-gas coupling system, EGNS represents the expected value of natural gas shortage in the electricity-gas coupling system, and β car represents the carbon emissions of the electric-gas coupling system, η re represents the renewable energy consumption rate of the electricity-gas coupling system, C inv Indicates the construction measurement value of the electric-gas coupling system, C ope Indicates the operating measurement value of the electric-pneumatic coupling system.

[0070] The normalization processing formula is as follows:

[0071]

[0072] Where y represents the indicator to be processed in the normalized processing formula f(), which includes positive indicators and negative indicators. When y is a positive indicator, the positive indicator y includes the renewable energy consumption rate η of the electricity-gas coupling system. re , the larger the value of the positive index y, the higher the comprehensive evaluation value of the configuration demand Ω; when y is a negative index, the negative index y includes the expected value of electricity shortage EENS of the electric-gas coupling system, the expected value of natural gas shortage EGNS, and the carbon emission β car , Construction measurement value C inv and operating measurement value C ope , the smaller the value of the reverse index y, the higher the comprehensive evaluation value of the configuration demand and the higher Ω.

[0073] Expected electricity shortage EENS, expected natural gas shortage EGNS, and carbon emissions β of the electricity-gas coupling system car , renewable energy consumption rate η re , Construction measurement value C of electric-gas coupling system inv and the operating measurement value C of the electric-pneumatic coupling system ope , as follows:

[0074]

[0075]

[0076] Wherein, T represents the time period between the current moment and the initial moment of operation of the electric-gas coupling system or the moment of the last repetition of step 2); the T period is divided into several time periods, and Ξs represents the set of the time periods divided by the T period; p s represents the probability of time period s within time period T, p s =s / T; It represents the power loss load of the electric-gas coupling system in the period s; represents the natural gas load loss of the electricity-gas coupling system in the period s; N represents the number of power stations in the power network of the electricity-gas coupling system, represents the carbon emission factor of the g-th power station among all power stations in the power network, represents the power generation of the gth power station in the power grid; M represents the number of gas production stations in the natural gas network of the power-gas coupling system, represents the carbon emission factor of the jth gas production station in the natural gas network, represents the gas production of the jth gas production station in the natural gas network; Q H Represents the total power generation of each power station in the power network; Q L Indicates the total gas production of each gas production station in the natural gas network; Q re represents the total renewable energy generation of each new energy generator in the power network; Q e Indicates the total amount of electricity transmitted by the electric-gas coupling system; t represents the net present value of the electric-gas coupling system at time t, o t =1 / (1+d) -t , d represents the discount rate of the electric-gas coupling system at time t; ψ (·) represents the residual value rate of fixed assets of the equipment to be built in the electric-gas coupling system; C (·) Indicates the construction measurement value of the equipment to be built in the electric-gas coupling system; (·),t and z (·),t-1 Respectively represent the status of the equipment to be built in the electric-gas coupling system at time t and time t-1. If the status of the equipment to be built is 1, it is 0 if it does not exist; z g,t Indicates the status of the g-th power station in the power grid at time t, which is 1 if it exists and 0 if it does not exist; C gen represents the unit capacity fuel measurement value of the g-th power station among the power stations in the power network; g,t represents the output power of the gth power station among the power stations in the power network at time t; z a,t Indicates the state of the ath gas source in the natural gas network, which is 1 if it exists and 0 if it does not exist; C wellIndicates the unit capacity fuel measurement value of the ath gas source among all gas sources in the natural gas network; G a,t It represents the gas flow of the ath gas source among all gas sources in the natural gas network at time t.

[0077] The probability of failure of the power-gas coupling system under uncertainty is characterized by the expected value of electricity shortage EENS and the expected value of natural gas shortage EGNS; the probability of failure of the power-gas coupling system under uncertainty is characterized by the expected value of carbon emissions β car and renewable energy consumption rate η re etc. to characterize the possible impact of the electric-gas coupling system on the environment; according to the construction measurement value C inv and operating measurement value C ope To characterize the cost of construction and operation of the electric-gas coupling system.

[0078] The preset evaluation threshold combines the development status and goals of the electric-gas coupling system. The role of the evaluation threshold is to compare the comprehensive evaluation value of the configuration requirements calculated by the evaluation subsystem of the digital twin system each time with the evaluation threshold to decide whether to trigger the configuration, that is, whether to continue with step 4).

[0079] Step 4) Obtain some operating data of the electric-gas coupling system at this time and input it into the dynamic configuration model of the electric-gas coupling system of the configuration subsystem of the digital twin system. The dynamic configuration model of the electric-gas coupling system outputs the construction data of the electric-gas coupling system at this time, and adjusts the equipment model library of the current equipment model subsystem according to the construction data.

[0080] In step 4), the dynamic configuration model of the electric-pneumatic coupling system of the configuration subsystem of the digital twin system is specifically the electric-pneumatic coupling system objective function considering the constraints of the electric-pneumatic coupling system. The electric-pneumatic coupling system objective function is as follows:

[0081]

[0082] in, The shedding measurement value representing the unit load shedding amount of the power network at time t; The EENS is the load shedding value of the natural gas network at time t. t It represents the expected load shedding of the power network at time t; EGNS t represents the expected load shedding of the natural gas network at time t.

[0083] The several operating data of the electric-gas coupling system obtained in step 4) at this time include the unit capacity fuel metering value C of the g-th power station among the various power stations in the power network. gen , the unit capacity fuel measurement value C of the ath gas source among all gas sources in the natural gas network well , the cut-off value of the unit load shedding of the power network at time t And the removal measurement value of the unit load shedding of the natural gas network at time t The construction data of the electric-gas coupling system output by the dynamic configuration model of the electric-gas coupling system includes the state z of the equipment to be built in the electric-gas coupling system at time t. (·),t , the output power P of the gth power station among the power stations in the power network at time t g,t , the gas flow rate G of the ath gas source among the gas sources in the natural gas network at time t a,t , the expected load shedding EENS of the power network at time t t And the expected load shedding of the natural gas network at time t EGNS t .

[0084] The constraints of the electric-gas coupling system are as follows:

[0085] z (·),t-1 ≤z (·),t

[0086] EENS t ≤EENS set

[0087] EGNS t ≤EGNS set

[0088]

[0089] Among them, EENS set Indicates the expected value of the power network's preset power supply shortage; EGNS set represents the expected value of the preset natural gas supply shortage of the natural gas network; B r represents the power network constant term, B r is a very large number; z r,t Indicates the state of the r-th transmission line in the power network at time t, which is 1 if it exists and 0 if it does not exist. The r-th transmission line in the power network is located between the u-th node and the v-th node in the power network. uv,t P represents the active power transmitted from the uth node to the vth node on the rth transmission line in the power network at time t, vu,t represents the active power transmitted from the vth node to the uth node on the rth transmission line in the power network at time t; θ u,t represents the power angle of the u-th node in the power network at time t, θ v,t represents the power angle of the vth node in the power network at time t; w rrepresents the reactance of the rth transmission line in the power network; each gas generator, non-gas generator and energy storage device in the power network constitutes a first device set, the first device set includes a plurality of first devices, and They represent the minimum power and maximum power of the qth first device in the first device set, P q,t represents the output power of the qth first device in the first device set at time t, z q,t Indicates the status of the qth first device in the power network, which is 1 if it exists and 0 if it does not exist; ur K represents the element in row u and column r of the correlation matrix composed of each node and each transmission line in the power network. uq represents the element in the uth row and qth column of the association matrix formed by each node in the power network and each first device in the first device set, W uε Represents the element in the uth row and the εth column of the association matrix composed of each node and each electric load device in the power network, represents the normal load level of the εth load device in the power network at time t; V uδ Represents the element in the uth row and the δth column of the association matrix formed by each node in the power network and each power-to-gas device in the natural gas network; represents the power-to-gas load level of the δth power-to-gas device in the natural gas network at time t; z p,t Indicates the state of the pth gas pipeline in the electric-gas coupling system, which is 1 if it exists and 0 if it does not exist; C p represents the transmission coefficient of the pth gas pipeline in the electric-gas coupling system; and is a 0-1 variable. The pth gas transmission pipeline in the electric-gas coupling system connects the mth node and the nth node of the natural gas network. When the natural gas in the pth gas transmission pipeline flows out from the mth node, When the natural gas in the pth gas pipeline flows into the mth node, π m,t and π n,t They represent the square of the gas pressure at the mth node and the nth node in the natural gas network at time t; Z p,t represents the gas flow rate of the pth gas transmission pipeline in the electric-gas coupling system at time t; and They represent the maximum and minimum gas pressures of the mth node in the natural gas network respectively; and They represent the maximum transmission capacity and minimum transmission capacity of the pth gas pipeline in the electric-gas coupling system; z j,tIndicates the status of the j-th gas production station in the natural gas network, which is 1 if it exists and 0 if it does not exist. mj Represents the element in the mth row and jth column of the association matrix composed of each node and each gas production station in the natural gas network, G j,t represents the gas production of the jth gas production station in the natural gas network, z δ,t Indicates the status of the δth power-to-gas device in the natural gas network, which is 1 if it exists and 0 if it does not exist. mδ Represents the element in the mth row and the δth column of the association matrix composed of each node and each power-to-gas device in the natural gas network, represents the gas production of the δth power-to-gas device in the natural gas network at time t; Y mp represents the element in the mth row and pth column of the connection matrix composed of each node and each gas pipeline in the natural gas network; J mn Represents the element in the mth row and nth column of the connection matrix composed of each node and each gas load device in the natural gas network, It represents the gas load level of the nth gas load device in the natural gas network at time t, U mλ Represents the element in the mth row and λth column of the connection matrix composed of each node and each gas-to-electricity load device in the natural gas network, represents the gas-to-electricity load level of the λth gas-to-electricity load device in the natural gas network at time t; μ represents the conversion coefficient of electricity to natural gas; Γ represents the conversion coefficient of gas volume to electrical power; and η represents the energy conversion efficiency of the power-to-gas device in the natural gas network.

[0090] The several operating data of the electric-gas coupling system obtained in step 4) at this time also include the expected power shortage value EENS of the power network set , the expected gas supply shortage value EGNS of the natural gas network set , the conventional electric load level of the εth electric load device in the power network at time t The gas load level of the nth gas load device in the natural gas network at time t The construction data of the electric-gas coupling system output by the dynamic configuration model of the electric-gas coupling system at this time also includes the active power P transmitted from the uth node to the vth node on the rth transmission line in the power network at time t uv,t , the active power P transmitted from the vth node to the uth node on the rth transmission line in each transmission line in the power network at time t vu,t , the output power P of the qth first device in the first device set at time t q,t , the power-to-gas load level of the δth power-to-gas device in the natural gas network at time t The square of the gas pressure at the mth node and the nth node in the natural gas network at time t is πm,b,t and π n,b,t , the gas flow Z of the pth gas pipeline in the electric-gas coupling system at time t p,b,t , the gas production of the δth power-to-gas device in the natural gas network at time t And the gas-to-electricity load level of the λth gas-to-electricity load device in the natural gas network at time t

[0091] The device model library of the current device model subsystem is adjusted according to the construction data. Specifically, the device model of the device to be built with the existing status is added to the device model library of the device model subsystem according to the construction data, and the output power of each power station in the power network, the gas flow of each gas source in the natural gas network, the expected load shedding amount of the power network and the expected load shedding amount of the natural gas network, the active power on each transmission line in the power network, the output power of each first device in the first device set, the power-to-gas load level of each power-to-gas device in the natural gas network, the square of the gas pressure of each node in the natural gas network, the gas flow of each gas transmission pipeline in the electric-gas coupling system, the gas production of each power-to-gas device in the natural gas network, and the gas-to-power load level of each gas-to-power load device in the natural gas network are adjusted according to the output construction data.

[0092] Step 5) obtains each device model in the device model library of the device model subsystem adjusted in step 4) and performs the same operations on the device models in steps 2)-3). If the obtained comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, the current device model library is placed in the project library of the configuration subsystem; if the obtained comprehensive evaluation value of the configuration requirements is less than the evaluation threshold, repeat steps 4)-5) until the obtained comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, and the current device model library is placed in the project library of the configuration subsystem.

[0093] In step 5), after repeating steps 4)-5), the final construction data of the electric-gas coupling system is obtained, that is, the construction data of the current device model library, and the current device model library is placed in the project library of the configuration subsystem.

[0094] Each device model library and the comprehensive evaluation value of configuration requirements in the problem library correspond to each device model library in the project library. The problem library and the project library can be used as historical data sets. When the device model library of a certain electric-gas coupling system corresponds to one of the device model libraries in the problem library, the configuration of the electric-gas coupling system can be directly adjusted according to the device model library in the project library corresponding to the device model library in the problem library.

[0095] Step 6) reconfiguring the electric-pneumatic coupling system according to the current device model library put into the project library in step 5) to achieve dynamic configuration of the electric-pneumatic coupling system.

[0096] In step 6), the electric-gas coupling system is reconfigured according to the current device model library put into the project library in step 5), that is, the electric-gas coupling system is constructed and adjusted according to the final construction data to realize the dynamic configuration of the electric-gas coupling system.

[0097] The specific embodiments of the present invention are as follows:

[0098] Taking the electric-gas coupling system composed of IEEE 30 nodes and Belgium 20 nodes as an example, the specific implementation of the present invention is described in detail in combination with the technical solution and the drawings.

[0099] The coupled power-gas system consists of an electric power system and a natural gas system. The electric power system includes 30 electric nodes and 41 electric power lines. The natural gas system includes 20 gas nodes and 19 natural gas pipelines. Six generators are connected to electric nodes 1, 2, 5, 8, 11, and 13. The generators at electric nodes 5, 8, and 13 are gas-fired generators, supplied by natural gas nodes 1, 3, and 5, respectively. The remaining generators are non-gas-fired generators. The gas source for natural gas system nodes 13 and 14 is a power-to-gas converter, supplied by electric nodes 2 and 5.

[0100] The equipment to be selected includes 10 power lines, 10 natural gas pipelines, 3 coal-fired generators, 2 gas-fired generators, and 2 energy storage devices. The equipment parameters of the equipment to be selected are shown in Table 1-4:

[0101] Table 1 Parameters of the selected power line equipment

[0102] Selected power lines First Node End Node Resistance / Ω Reactance / Ω Transmission capacity / MW Length / km Selected power line 1 1 2 0.02 0.06 160 62 Selected power line 2 2 4 0.03 0.07 120 86 Selected power line 3 2 5 0.07 0.15 120 74 Selected power line 4 2 6 0.01 0.02 120 84 Selected power lines 5 4 6 0.06 0.2 160 62 Selected power lines 6 5 7 0.0497 0.192 175 50 Power lines to be selected 7 6 7 0.0057 0.0447 100 62 Power lines to be selected 8 10 19 0.0057 0.0447 100 62 Power lines to be selected 9 10 24 0.0057 0.0447 100 62 Selected power lines 10 14 18 0.0078 0.0606 100 84

[0103] Table 2 Parameters of selected natural gas pipeline equipment

[0104]

[0105]

[0106] Table 3 Parameters of the selected generators and storage devices

[0107] dynamo Electrical Node Capacity (MW) Coal-fired generator to be selected 1 30 80 Coal-fired generator 2 26 60 Coal-fired generator 3 17 50 Gas generator to be selected 1 15 90 Gas generator 2 17 50 Selected energy storage equipment 1 4 30 Selected energy storage equipment 2 10 20

[0108] Table 4 Parameters of the selected gas production station equipment

[0109] Selected gas production stations Gas Node <![CDATA[Capacity (10 3 m 3 )]]> Candidate gas production station 1 7 35 Candidate gas production station 2 17 28 Candidate gas production station 3 16 22 Gas production station 4 to be selected 20 22 Candidate gas production station 5 4 27

[0110] In step 2), the device parameters before and after correction are shown in Table 5-6 below:

[0111] Table 5 Equipment parameters of original power lines 1-10 before and after modification

[0112]

[0113] Table 6 Parameters of original natural gas pipelines 1-10 before and after correction

[0114]

[0115]

[0116] In step 3), using a preset one-year time period as the evaluation interval, the comprehensive evaluation value of the configuration requirements for the first year is 0.918. This indicates that the comprehensive evaluation value for the configuration requirements for the first year exceeds the evaluation threshold, so no equipment is invested in the first year. If the comprehensive evaluation value for the configuration requirements falls below the evaluation threshold, the current scenario is added to the problem library of the digital twin system, and dynamic configuration of the electrical-gas coupling system is performed to obtain a configuration solution for the current scenario. Since the results for years 2, 4, 5, 7, 9, and 10 do not meet the requirements, these years are added to the problem library.

[0117] The configuration scheme is used as the boundary, that is, it is invested and constructed in the device model library of the digital twin system. If the comprehensive evaluation value of the configuration requirements of the post-configuration scenario is better than the evaluation threshold, the configuration scheme is placed in the project library of the digital twin system, as shown in Table 7 below:

[0118] Table 7 Dynamic configuration scheme

[0119]

[0120] Finally, construction is carried out in the electric-gas coupling system according to the configuration plan in the project library, and the investment and construction of new equipment are completed within the specified period.

Claims

1. A dynamic configuration method for an electric-pneumatic coupling system based on digital twin technology, characterized by: The steps include: Step 1) Establishing a digital twin system of the electro-pneumatic coupling system, including a device model subsystem, an evaluation subsystem, and a configuration subsystem; the device model subsystem includes a device model library, which includes device models of various devices in the electro-pneumatic coupling system; the evaluation subsystem includes a configuration requirement evaluation model for the electro-pneumatic coupling system; and the configuration subsystem includes an electro-pneumatic coupling system dynamic configuration model, a problem library, and a project library; Step 2) for each device in the electric-pneumatic coupling system and its device model in the device model library of the device model subsystem, the device model including several device parameters of the device; after the electric-pneumatic coupling system has been running for a preset period of time, an observed value of each device parameter of each device at the current moment is obtained, and for each observed value of each device parameter at the current moment, the device parameter is corrected using a least squares method to obtain an estimated value of the device parameter and the device parameter in the device model is updated; Step 3) The estimated values ​​of the various device parameters in the device model of the device model library of the device model subsystem obtained in step 2) are used as the device parameters of the various devices in the electric-pneumatic coupling system in the current operating state. Several operating data of the electric-pneumatic coupling system at this time are obtained and input into the configuration requirement evaluation model of the electric-pneumatic coupling system of the evaluation subsystem. The configuration requirement evaluation model outputs a comprehensive evaluation value of the configuration requirements of the electric-pneumatic coupling system. The comprehensive evaluation value of the configuration requirements is compared with a preset evaluation threshold. When the comprehensive evaluation value of the configuration requirements is greater than the evaluation threshold, steps 2)-3) are repeated until the comprehensive evaluation value of the configuration requirements is less than the evaluation threshold. The device model library at this time and the comprehensive evaluation value of the configuration requirements calculated this time are then placed into the question library of the configuration subsystem. Step 4) Acquire some operating data of the electric-pneumatic coupling system at this time and input it into the electric-pneumatic coupling system dynamic configuration model of the configuration subsystem of the digital twin system. The electric-pneumatic coupling system dynamic configuration model outputs the construction data of the electric-pneumatic coupling system at this time, and adjusts the device model library of the current device model subsystem according to the construction data; Step 5) Obtain each device model in the device model library of the device model subsystem adjusted in step 4) and perform the same operations as for the device models in steps 2)-3). If the obtained configuration requirement comprehensive evaluation value is greater than the evaluation threshold, the current device model library is placed in the project library of the configuration subsystem; if the obtained configuration requirement comprehensive evaluation value is less than the evaluation threshold, repeat steps 4)-5) until the obtained configuration requirement comprehensive evaluation value is greater than the evaluation threshold, and the current device model library is placed in the project library of the configuration subsystem; Step 6) reconfiguring the electric-pneumatic coupling system according to the current device model library put into the project library in step 5) to achieve dynamic configuration of the electric-pneumatic coupling system.

2. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 1, characterized in that: In the above step 1), the electric-gas coupling system includes an electric power network and a natural gas network; the equipment in the electric power network includes a number of power stations, transmission lines, electric load equipment and energy storage equipment; the power station includes a number of gas generators and non-gas generators, each of which is interconnected through each transmission line, and the non-gas generator is specifically a coal-fired generator or a new energy generator; the electric load equipment is specifically a power-to-gas device or an electric load equipment; the energy storage equipment is specifically an energy storage battery; the equipment in the natural gas network includes a number of gas production stations, gas transmission pipelines and gas load equipment, and each gas production station and gas load equipment is connected through each gas transmission pipeline; A gas production station includes several gas sources and power-to-gas equipment; gas load equipment is specifically gas-to-electricity load equipment or gas load equipment, and gas-to-electricity load equipment is specifically a gas generator; the power-to-gas equipment of the power network is connected to the nodes of the natural gas network, and the gas-to-electricity equipment of the natural gas network is connected to the nodes of the power network; the equipment in the power-gas coupling system includes gas generators, non-gas generators, electric load equipment, transmission lines and energy storage equipment in the power network, and gas sources, power-to-gas equipment, gas pipelines and gas load equipment in the natural gas network. Except for the transmission lines and gas pipelines, each equipment in the power-gas coupling system is located at its respective node.

3. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 1, characterized in that: In step 2), after the electric-gas coupling system runs for a preset period of time, the observed value of each device parameter of each device at the current moment is obtained. For each observed value of the device parameter at the current moment, the device parameter is corrected using the least squares method to obtain an estimated value of the device parameter, as follows: Among them, x i,k+1 represents the observed value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time. Represents the estimated value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time; The estimated value of the i-th device parameter in the device model library of the device model subsystem at the current k+1 time is used Update the i-th device parameter in the device model library of the device model subsystem at the current k+1 time.

4. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 2, characterized in that: In step 3), the configuration requirement evaluation model of the electric-pneumatic coupling system of the evaluation subsystem is specifically as follows: Ω=f(EENS)+f(EGNS)+f(β car )+f(η re )+f(C inv )+f(C ope ) Where Ω represents the comprehensive evaluation value of configuration requirements, f() represents the normalization formula, EENS represents the expected value of electricity shortage in the electricity-gas coupling system, EGNS represents the expected value of natural gas shortage in the electricity-gas coupling system, and β car represents the carbon emissions of the electric-gas coupling system, η re represents the renewable energy consumption rate of the electricity-gas coupling system, C inv Indicates the construction measurement value of the electric-gas coupling system, C ope Indicates the operating measurement value of the electric-pneumatic coupling system; The normalization processing formula is as follows: Where y represents the indicator to be processed in the normalized processing formula f(), which includes positive indicators and negative indicators. When y is a positive indicator, the positive indicator y includes the renewable energy consumption rate η of the electricity-gas coupling system. re , the larger the value of the positive index y, the higher the comprehensive evaluation value of the configuration demand Ω; when y is a negative index, the negative index y includes the expected value of electricity shortage EENS of the electric-gas coupling system, the expected value of natural gas shortage EGNS, and the carbon emission β car , Construction measurement value C inv and operating measurement value C ope , the smaller the value of the reverse index y, the higher the comprehensive evaluation value of the configuration demand and the higher Ω.

5. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 4, characterized in that: The expected electricity shortage value EENS, the expected natural gas shortage value EGNS, and the carbon emission value β of the electric-gas coupling system are car , renewable energy consumption rate η re , Construction measurement value C of electric-gas coupling system inv and the operating measurement value C of the electric-pneumatic coupling system ope , as follows: Wherein, T represents the time period between the current moment and the initial moment of operation of the electric-gas coupling system or the moment of the last repetition of step 2); the T period is divided into several time periods, and Ξs represents the set of the time periods divided by the T period; p s represents the probability of time period s within time period T, p s =s / T; It represents the power loss load of the electric-gas coupling system in the period s; represents the natural gas load loss of the electricity-gas coupling system in the period s; N represents the number of power stations in the power network of the electricity-gas coupling system, represents the carbon emission factor of the g-th power station among all power stations in the power network, represents the power generation of the gth power station in the power grid; M represents the number of gas production stations in the natural gas network of the power-gas coupling system, represents the carbon emission factor of the jth gas production station in the natural gas network, represents the gas production of the jth gas production station in the natural gas network; Q H Represents the total power generation of each power station in the power network; Q L Indicates the total gas production of each gas production station in the natural gas network; Q re represents the total renewable energy generation of each new energy generator in the power network; Q e Indicates the total amount of electricity transmitted by the electric-gas coupling system; t represents the net present value of the electric-gas coupling system at time t; ψ (·) Indicates the residual value rate of the equipment to be built in the electric-gas coupling system; C (·) Indicates the construction measurement value of the equipment to be built in the electric-gas coupling system; (·),t and z (·),t-1 Respectively represent the status of the equipment to be built in the electric-gas coupling system at time t and time t-1. If the status of the equipment to be built is 1, it is 0 if it does not exist; z g,t Indicates the status of the g-th power station in the power grid at time t, which is 1 if it exists and 0 if it does not exist; C gen represents the unit capacity fuel measurement value of the g-th power station among the power stations in the power network; g,t represents the output power of the gth power station among the power stations in the power network at time t; z a,t Indicates the state of the ath gas source in the natural gas network, which is 1 if it exists and 0 if it does not exist; C well Indicates the unit capacity fuel measurement value of the ath gas source among all gas sources in the natural gas network; G a,t represents the gas flow rate of the ath gas source among all gas sources in the natural gas network at time t; The operating data of the electric-gas coupling system obtained in step 3) at this time includes the power loss load of the electric-gas coupling system in the time period s Natural gas load loss The power generation of the gth power station among all power stations in the power network The gas production of the jth gas production station in the natural gas network The total power generation Q of each power station in the power network H , the total gas production Q of each gas production station in the natural gas network L , the total renewable energy generation capacity Q of each new energy generator in the power network re And the total amount of electricity transmitted by the electric-gas coupling system Q e .

6. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 2, characterized in that: In step 4), the dynamic configuration model of the electric-pneumatic coupling system of the configuration subsystem of the digital twin system is specifically the electric-pneumatic coupling system objective function considering the constraints of the electric-pneumatic coupling system. The electric-pneumatic coupling system objective function is specifically as follows: in, The shedding measurement value representing the unit load shedding amount of the power network at time t; The EENS is the load shedding value of the natural gas network at time t. t represents the expected load shedding of the power network at time t; EGNS t represents the expected load shedding of the natural gas network at time t; The several operating data of the electric-gas coupling system obtained in step 4) at this time include the unit capacity fuel metering value C of the g-th power station among the various power stations in the power network. gen , the unit capacity fuel measurement value C of the ath gas source among all gas sources in the natural gas network well , the cut-off value of the unit load shedding of the power network at time t And the removal measurement value of the unit load shedding of the natural gas network at time t The construction data of the electric-gas coupling system output by the dynamic configuration model of the electric-gas coupling system includes the state z of the equipment to be built in the electric-gas coupling system at time t. (·),t , the output power P of the gth power station among the power stations in the power network at time t g,t , the gas flow rate G of the ath gas source among the gas sources in the natural gas network at time t a,t , the expected load shedding EENS of the power network at time t t And the expected load shedding of the natural gas network at time t EGNS t .

7. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 6, characterized in that: The constraints of the electric-gas coupling system are as follows: With (·),t-1 ≤of (·),t AT SOMETIME t ≤AGREE set COUNTRY t ≤EGNS set Among them, EENS set Indicates the expected value of the power network's preset power supply shortage; EGNS set represents the expected value of the preset natural gas supply shortage of the natural gas network; B r represents the power network constant term; z r,t Indicates the state of the r-th transmission line in the power network at time t, which is 1 if it exists and 0 if it does not exist. The r-th transmission line in the power network is located between the u-th node and the v-th node in the power network. uv,t P represents the active power transmitted from the uth node to the vth node on the rth transmission line in the power network at time t, vu,t represents the active power transmitted from the vth node to the uth node on the rth transmission line in the power network at time t; θ u,t represents the power angle of the u-th node in the power network at time t, θ v,t represents the power angle of the vth node in the power network at time t; w r represents the reactance of the rth transmission line in the power network; each gas generator, non-gas generator and energy storage device in the power network constitutes a first device set, the first device set includes a plurality of first devices, and They represent the minimum power and maximum power of the qth first device in the first device set, P q,t represents the output power of the qth first device in the first device set at time t, z q,t Indicates the status of the qth first device in the power network, which is 1 if it exists and 0 if it does not exist; ur K represents the element in row u and column r of the correlation matrix composed of each node and each transmission line in the power network. uq represents the element in the uth row and qth column of the association matrix formed by each node in the power network and each first device in the first device set, W uε Represents the element in the uth row and the εth column of the association matrix composed of each node and each electric load device in the power network, V represents the load level of the εth load device in the power network at time t; uδ Represents the element in the uth row and the δth column of the association matrix formed by each node in the power network and each power-to-gas device in the natural gas network; represents the power-to-gas load level of the δth power-to-gas device in the natural gas network at time t; z p,t Indicates the state of the pth gas pipeline in the electric-gas coupling system, which is 1 if it exists and 0 if it does not exist; C p represents the transmission coefficient of the pth gas pipeline in the electric-gas coupling system; and is a 0-1 variable. The pth gas transmission pipeline in the electric-gas coupling system connects the mth node and the nth node of the natural gas network. When the natural gas in the pth gas transmission pipeline flows out from the mth node, When the natural gas in the pth gas pipeline flows into the mth node, π m,t and π n,t They represent the square of the gas pressure at the mth node and the nth node in the natural gas network at time t; Z p,t represents the gas flow rate of the pth gas transmission pipeline in the electric-gas coupling system at time t; and They represent the maximum and minimum gas pressures of the mth node in the natural gas network respectively; and They represent the maximum transmission capacity and minimum transmission capacity of the pth gas pipeline in the electric-gas coupling system; z j,t Indicates the status of the j-th gas production station in the natural gas network, which is 1 if it exists and 0 if it does not exist. mj Represents the element in the mth row and jth column of the association matrix composed of each node and each gas production station in the natural gas network, G j,t represents the gas production of the jth gas production station in the natural gas network, z δ,t Indicates the status of the δth power-to-gas device in the natural gas network, which is 1 if it exists and 0 if it does not exist. mδ Represents the element in the mth row and the δth column of the association matrix composed of each node and each power-to-gas device in the natural gas network, represents the gas production of the δth power-to-gas device in the natural gas network at time t; Y mp represents the element in the mth row and pth column of the connection matrix composed of each node and each gas pipeline in the natural gas network; J mn Represents the element in the mth row and nth column of the connection matrix composed of each node and each gas load device in the natural gas network, It represents the gas load level of the nth gas load device in the natural gas network at time t, U mλ Represents the element in the mth row and λth column of the connection matrix composed of each node and each gas-to-electricity load device in the natural gas network, represents the gas-to-electricity load level of the λth gas-to-electricity load device in the natural gas network at time t; μ represents the conversion coefficient of electricity to natural gas; Γ represents the conversion coefficient of gas volume to electrical power; η represents the energy conversion efficiency of the power-to-gas device in the natural gas network; The several operating data of the electric-gas coupling system obtained in step 4) at this time also include the expected power shortage value EENS of the power network set , the expected gas supply shortage value EGNS of the natural gas network set , the load level of the εth load device in the power network at time t The gas load level of the nth gas load device in the natural gas network at time t The construction data of the electric-gas coupling system output by the dynamic configuration model of the electric-gas coupling system at this time also includes the active power P transmitted from the uth node to the vth node on the rth transmission line in the power network at time t uv,t , the active power P transmitted from the vth node to the uth node on the rth transmission line in each transmission line in the power network at time t vu,t , the output power P of the qth first device in the first device set at time t q,t , the power-to-gas load level of the δth power-to-gas device in the natural gas network at time t The square of the gas pressure at the mth node and the nth node in the natural gas network at time t is π m,b,t and π n,b,t , the gas flow Z of the pth gas pipeline in the electric-gas coupling system at time t p,b,t , the gas production of the δth power-to-gas device in the natural gas network at time t And the gas-to-electricity load level of the λth gas-to-electricity load device in the natural gas network at time t The device model library of the current device model subsystem is adjusted according to the construction data. Specifically, the device model of the device to be built with the existing status is added to the device model library of the device model subsystem according to the construction data, and the output power of each power station in the power network, the gas flow of each gas source in the natural gas network, the expected load shedding amount of the power network and the expected load shedding amount of the natural gas network, the active power on each transmission line in the power network, the output power of each first device in the first device set, the power-to-gas load level of each power-to-gas device in the natural gas network, the square of the gas pressure of each node in the natural gas network, the gas flow of each gas transmission pipeline in the electric-gas coupling system, the gas production of each power-to-gas device in the natural gas network, and the gas-to-power load level of each gas-to-power load device in the natural gas network are adjusted according to the output construction data.

8. The method for dynamic configuration of an electric-pneumatic coupling system based on digital twin technology according to claim 7, characterized in that: In step 5), after repeating steps 4)-5), the final construction data of the electric-gas coupling system is obtained, that is, the construction data of the current device model library, and the current device model library is placed in the project library of the configuration subsystem; In step 6), the electric-gas coupling system is reconfigured according to the current device model library placed in the project library in step 5), that is, the electric-gas coupling system is constructed and adjusted according to the final construction data to achieve dynamic configuration of the electric-gas coupling system.

Images