CN113704977A - Decomposition method based integrated energy system event driving type simulation method - Google Patents
Decomposition method based integrated energy system event driving type simulation method Download PDFInfo
- Publication number
- CN113704977A CN113704977A CN202110893014.0A CN202110893014A CN113704977A CN 113704977 A CN113704977 A CN 113704977A CN 202110893014 A CN202110893014 A CN 202110893014A CN 113704977 A CN113704977 A CN 113704977A
- Authority
- CN
- China
- Prior art keywords
- network
- disturbance
- simulation
- solving
- energy system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 26
- 230000008878 coupling Effects 0.000 claims abstract description 52
- 238000010168 coupling process Methods 0.000 claims abstract description 52
- 238000005859 coupling reaction Methods 0.000 claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 238000005094 computer simulation Methods 0.000 abstract description 10
- 230000003993 interaction Effects 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 16
- 230000008859 change Effects 0.000 description 13
- 238000005485 electric heating Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 239000008400 supply water Substances 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000005338 heat storage Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a comprehensive energy system event driving type simulation method based on a decomposition method, belonging to the field of energy system simulation; a method for event-driven simulation of an integrated energy system based on decomposition comprises the following steps: s1, an electric, gas and thermal energy network model and an equipment model are established, S2 a comprehensive energy system model is established, the model is divided into different modules according to the type of the energy network and the characteristics of the equipment, S3 system parameters, system normal state load and system fault conditions are obtained, the solving sequence of each module is determined according to different fault types, and the parameters are transmitted to the next module after the modules are solved independently until the system is solved completely. The method fully considers the characteristics of strong coupling of the comprehensive energy system and frequent parameter interaction between equipment and networks in dynamic simulation, divides the whole system into different modules by adopting a decomposition method, and determines the solving sequence of the modules according to different fault positions to realize an event-driven simulation strategy, thereby effectively improving the simulation efficiency.
Description
Technical Field
The disclosure belongs to the technical field of energy system simulation, and particularly relates to a comprehensive energy system event driving type simulation method based on a decomposition method.
Background
The traditional energy system is mainly operated in a discrete mode, and most mainstream energy systems such as coal, petroleum, natural gas, electric power and cooling/heating systems are operated respectively, so that a closed industry and interest barrier are formed. Due to the barriers, the mutual complementary interaction characteristic and the cascade utilization potential of various energy sources cannot be exerted, so that the comprehensive utilization efficiency of the energy sources is reduced, the carbon emission is increased to a certain extent, and the problem of environmental pollution is aggravated.
Under the background that the world Energy System is in urgent need of transformation and upgrading, the concept of Integrated Energy System (IES) is widely supported by academia and industry. The comprehensive energy system can break the discrete operation barriers of the traditional energy system, fully exerts complementary interaction characteristics and cascade utilization potential among various forms of energy (including coal, petroleum, natural gas, electric power, heat/cold, wind, light, water, hydrogen and the like), thereby improving the comprehensive energy efficiency of the energy system, reducing carbon emission and pollutant emission, improving the stability and reliability of the energy supply system, and promoting the absorption of renewable energy, and has bright development prospect.
The simulation analysis of the comprehensive energy system is an experimental verification of a system construction scheme, and the real physical process occurring in the actual system is reproduced by establishing a device and pipe network working condition model according with the actual condition, the state of the system under each scene is simulated, so that safety check and operation analysis are provided for IES, and the simulation analysis plays a key role in the development and maintenance of the comprehensive energy system. The deep research of the simulation algorithm can improve the accuracy and the reliability of system simulation, provide powerful guarantee for the effectiveness of the construction scheme and have profound influence on the development of the comprehensive energy system.
The simulation of the comprehensive energy system is divided into a steady state and a dynamic state, the steady state emphasizes the working condition of the system after stabilization in a certain scene, the change process of the state along with time is ignored in the simulation process, the change condition of the depicting working condition along with time is emphasized in the dynamic simulation, the refinement degree of a model considered in the simulation process is higher, the problems of pipeline section crossing, overload and the like can be judged by analyzing the change of the working condition under disturbance, the safety and the stability of the system are further improved, meanwhile, the fineness of multi-energy cooperative scheduling can be improved through the dynamic simulation, the energy interaction between energy subsystems under a strong coupling scene is enhanced, the advantages of IES energy complementation and cascade utilization are fully played, and the new energy consumption capability and the user experience degree of the system are improved. In the comprehensive energy system, a strong coupling relation exists in the system, and the significance of searching for a proper method to process the coupling in the system so as to realize dynamic simulation is great. Meanwhile, the dynamic simulation mainly focuses on solving system parameters under disturbance, and the high-precision and rapid solving of the comprehensive energy system under disturbance is also a core problem of the dynamic simulation.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a comprehensive energy system event driving type simulation method based on a decomposition method, which effectively improves the simulation efficiency, ensures the solving speed of the system in various scenes and provides powerful guarantee for the safety analysis of the comprehensive energy system.
The purpose of the disclosure can be realized by the following technical scheme:
a method for event-driven simulation of an integrated energy system based on decomposition comprises the following steps: the method comprises the steps of synthesizing an energy system decomposition method, an event driving type simulation strategy, a conventional state solving strategy, a solving strategy under heat network disturbance, a solving strategy under air network disturbance and a solving strategy under power network disturbance;
the comprehensive energy system decomposition method divides a system into different modules according to different network types and equipment characteristics, then obtains corresponding module solving sequences according to different disturbance types by adopting a time-driven strategy, can divide the disturbance into heat supply network disturbance, air network disturbance and power grid disturbance according to the disturbance occurrence position, and can cause the known quantity of system load and different changes of a pipe network topological structure due to the disturbance at different positions, so that the corresponding solving sequences are obtained.
Further, the method for decomposing the comprehensive energy system comprises the following specific steps of:
step 1: determining the composition of a comprehensive energy system;
step 2: inputting network parameters, initial values and equipment operation modes;
and step 3: based on step 1 and step 2, the system is divided into different modules according to type, characteristics and operation mode, such as: heat supply network, gas network, power network module, etc.;
and 4, step 4: after the single module is solved, the parameters are transmitted to the next module.
Further, the conventional state solving strategy refers to: the system does not encounter disturbance and runs stably, at the moment, the system decomposition and solution has no hard sequence requirement, the coupling relation of the system needs to be deeply analyzed, the module with the weakest coupling and the lowest relevance with other modules is searched, and the solution is gradually carried out from the position.
Further, the event-driven type refers to: determining a module solving sequence according to a disturbance condition of a system in a simulation period and following a certain mechanism;
when the system is always in a conventional state, the simulation precision and speed can be ensured by solving the strategy according to the conventional state; if disturbance exists in the system, the system can be divided into different types of disturbance according to different disturbance positions, when the different types of disturbance occur, the disturbance network module needs to be solved preferentially as much as possible, then the coupling relation between the known quantity and the system module after the module solution is considered, a proper solving mode is selected according to the coupling relation, and if the disturbance network still has unknown quantity under the fault, the network where the unknown quantity exists is solved first, so that the solving priority of the fault network is ensured.
Further, the solving strategy under the disturbance of the heat supply network refers to: when the system generates heat supply network disturbance, whether the heat supply network can be directly solved or not is judged firstly based on the event driving type simulation strategy, the heat supply network disturbance can be directly simulated and solved generally, the coupling strength of other networks is analyzed after the solution, for an energy system with electricity, gas and heat networks, the coupling strength of the gas network is generally strong, the power grid can be solved firstly, and then the gas network is calculated
Further, the solving strategy under the disturbance of the air network refers to: when the air network is disturbed, whether the air network can be directly solved or not is judged firstly based on the event driving type simulation strategy, the air network can be directly calculated under common faults, and the coupling of the electric network and the heat network is equivalent after the working condition of the air network is obtained, so that the decomposition iteration solution through the thermoelectric network can be considered, and the high-precision calculation of the system is realized, and the specific steps are as follows:
further, the solving strategy under the power grid disturbance means that: when the system is disturbed by the power grid, whether the power grid can be directly solved or not is judged firstly based on the event-driven simulation strategy, and the power grid generally bears the task of supplying energy to other networks, so that direct calculation cannot be carried out, the network serving as the electric load is searched for to be solved, the power grid calculation is carried out after the solution is completed, and the remaining networks are calculated if the network which is not solved still exists.
The beneficial effect of this disclosure: aiming at the problems that parameters among subsystems are frequently interacted, the coupling of working conditions is strong, the solving sequence is difficult to determine in a disturbance scene and the simulation is difficult to implement in the dynamic simulation of the comprehensive energy system, the invention divides the system into a plurality of modules by adopting a decomposition method according to the characteristics of the system, provides an event-driven strategy for determining the simulation sequence according to the disturbance event type, determines the solving sequence of the modules and realizes the high-precision rapid solving of the whole system.
Drawings
In order to more clearly illustrate the embodiments or technical solutions in the prior art of the present disclosure, the drawings used in the description of the embodiments or prior art will be briefly described below, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a flow chart of a decomposition-based event-driven simulation strategy for an integrated energy system;
FIG. 2 is a block diagram of an embodiment of the present invention;
FIG. 3 is a flow chart of a simulation of final determination in a conventional state according to an embodiment of the present invention;
FIG. 4 is a flow chart of a simulation of a final determination under a disturbance of an air network according to an embodiment of the present invention;
FIG. 5 is a diagram illustrating simulation results of a conventional state in an embodiment of the present invention;
FIG. 6 is a diagram of simulation results under thermal network disturbance in an embodiment of the present invention;
FIG. 7 is a diagram of simulation results under grid disturbance in the embodiment of the present invention;
FIG. 8 is a diagram of simulation results under disturbance of the air network in the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present disclosure will be clearly and completely described below with reference to the drawings in the embodiments of the present disclosure, and it is obvious that the described embodiments are only a part of the embodiments of the present disclosure, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments disclosed herein without making any creative effort, shall fall within the protection scope of the present disclosure.
As shown in fig. 2, in this embodiment, there are three energy networks of a heat supply network, an air network and a power network, and a plurality of energy coupling devices such as a gas boiler, a steam turbine and an electric boiler, and parameters of each energy network and device do not substantially affect the method of the present invention, and therefore, the details are not described herein again. In the embodiment, the co-simulation duration is 6000s, and four scenes of conventional state, heat supply network, power grid and air network disturbance are simulated respectively. The conventional state simulation flow obtained by the event-driven policy flow shown in fig. 1 is shown in fig. 3, and includes the following specific steps:
step 1: building a comprehensive energy system;
step 2: inputting heat supply network, gas network, power network topology, pipeline, equipment parameters and load;
and step 3: the system is divided into a power grid, a heat supply network, an air network module, a steam-water heat exchanger, a water-water heat exchanger, an electric boiler and other equipment modules.
And 4, step 4: the system is found to have no fault, conventional state simulation is carried out, the heat supply network is found to appear in the system only as a load, the coupling is the weakest, the air supply network is used as electricity, the energy supply coupling of the heat supply network is strong, the coupling is centered even if the energy supply end of the heat supply network is the load of the air supply network, and therefore the heat supply network simulation is carried out firstly.
And 5: inputting an initial value of a heat supply network, and carrying out iteration in a heat supply network module;
step 6: transmitting the parameters of the heat supply network to the electric heating coupling equipment module and the gas-heat coupling module;
and 7: simulating the electric heating coupling equipment module, and transmitting the output parameters to the power grid module;
and 8: carrying out iterative calculation in the power grid module;
and step 9: transmitting the power grid parameters to an electrical coupling module;
step 10: simulating the electric coupling module and the gas-heat coupling module, and transmitting the output parameters to the gas network module;
step 11: simulating the air network, and finishing the simulation of the system after the simulation is finished;
the total simulation duration of the conventional state in this embodiment is 6000s, the space-time step length of the conventional state of the heat supply network is 4s/5m, the simulation time step length of the conventional time of the power network is 750s, the space-time step length of the conventional time of the air network is 40s/1000m, and the simulation result is shown in fig. 5. Graph (a) shows the supply and return water temperature of the heat source node 27, graph (b) shows the supply and return water temperature of the heat source node 28, graph (c) shows the generated power of the generator, and graph (d) shows the pressure change of each node of the air network. The simulation actual time of the reference example is 500s, and because the heat load is constant, the water supply temperature is kept stable after the influence of the initial temperature of the heat supply network is eliminated, the corresponding power grid load flow calculation result does not change too much, and the dynamic change process of the air pressure of each node in the air network simulation result is obvious and accords with the change condition of the node load.
The embodiment of the setting is that the heat supply network has faults: no. 11 node has working fault in 2000-2600s, and cannot supply heat, and the simulation steps are as follows:
step 1: building a comprehensive energy system;
step 2: inputting heat supply network, gas network, power network topology, pipeline, equipment parameters and load;
and step 3: the system is divided into a power grid, a heat supply network, an air network module, a steam-water heat exchanger, a water-water heat exchanger, an electric boiler and other equipment modules.
And 4, step 4: and (5) discovering the system heat supply network fault and carrying out fault state simulation.
And 5: according to the event driving type simulation strategy, the heat supply network can be directly simulated, the initial value of the heat supply network is input, and iteration in a heat supply network module is carried out;
step 6: analyzing other modules, wherein a heat supply network in the system is the electric load of a power grid, and the power grid and the heat supply network are both the gas loads of a gas network, so that the coupling of the power grid is weak, and the parameters of the heat supply network are transmitted to the electric heating coupling equipment module and the gas-heat coupling module;
and 7: simulating the electric heating coupling equipment module, and transmitting the output parameters to the power grid module;
and 8: carrying out iterative calculation in the power grid module;
and step 9: transmitting the power grid parameters to an electrical coupling module;
step 10: simulating the electric coupling module and the gas-heat coupling module, and transmitting the output parameters to the gas network module;
step 11: simulating the air network, and finishing the simulation of the system after the simulation is finished;
the dynamic simulation step length of each module under the heat supply network fault is consistent with the conventional state, and the simulation result is shown in fig. 6. Due to the particularity of the heat supply network fault, the fault simulation sequence is consistent with the conventional state. In the graphs (a) and (b), the supply and return water temperatures of the two heat source nodes 27 and 28 are changed respectively, the return water temperature is not influenced because the node 28 is far away from the fault node, but the return water temperature of the node 27 is increased due to the return water temperature of the node 11. The graph (c) is the generated power change condition of the generator in the fault, the graph (d) is the air pressure change of the air network node, the obvious pressure fluctuation of the air network can be observed, the depiction is accurate, and the air pressure is slightly reduced at the moment because the power of the generator is reduced and the air load is reduced, and the fault is ended and is increased again. And (e) the temperature of the supply water and the return water of the fault node is changed, and the temperature of the supply water of the fault node is stable after the load is lost and the temperature of the return water is greatly increased because the temperature of the heat supplied by the heat source is unchanged.
The power grid fault of the embodiment is as follows: the other generator of the power grid loses the working capacity in 2000-2600s, and the simulation steps are as follows:
step 1: building a comprehensive energy system;
step 2: inputting heat supply network, gas network, power network topology, pipeline, equipment parameters and load;
and step 3: the system is divided into a power grid, a heat supply network, an air network module, a steam-water heat exchanger, a water-water heat exchanger, an electric boiler and other equipment modules.
And 4, step 4: and discovering a system power grid, and performing fault state simulation.
And 5: although the system fault is located in the power grid, the heat supply network is used as the electric load of the power grid, if the heat supply network cannot be solved, the power grid calculation cannot be carried out, so that the heat supply network simulation is firstly carried out, parameters are transmitted to the power grid after the simulation is finished, initial values of the heat supply network are input, and the iteration in a heat supply network module is carried out;
step 6: transmitting the parameters of the heat supply network to the electric heating coupling equipment module and the gas-heat coupling module;
and 7: simulating the electric heating coupling equipment module, and transmitting the output parameters to the power grid module;
and 8: carrying out iterative calculation in the power grid module;
and step 9: transmitting the power grid parameters to an electrical coupling module;
step 10: simulating the electric coupling module and the gas-heat coupling module, and transmitting the output parameters to the gas network module;
step 11: simulating the air network, and completing system simulation;
the dynamic simulation step length of each module under the power grid fault is consistent with the conventional state, and the simulation result is shown in fig. 7. The graphs (a) and (b) respectively show the temperature changes of the water supply and return of the two heat source nodes 27 and 28, the heat source temperature does not change greatly because the heat supply network is only the electric load of the power grid and the coupling relation is not strong, and the graph (c) shows the power of the generator, and the power of the generator is increased and the load of the air grid is increased because the other generator cannot supply power during the fault period so as to meet the power balance of the system. The graph (d) shows the change of the air pressure of each node of the air network, and the air pressure changes correspondingly under the influence of the rise of the air load.
Example air network failure: the gas network energy coupling node is influenced by the limitation of gas supply amount, the gas mass flow is limited to a fixed value smaller than the conventional requirement in 2000-2600s, meanwhile, in order to meet the power balance of the system, a power grid cuts off certain electric load, a simulation flow obtained by event driving is shown in fig. 4, and simulation steps are as follows:
step 1: building a comprehensive energy system;
step 2: inputting heat supply network, gas network, power network topology, pipeline, equipment parameters and load;
and step 3: the system is divided into a power grid, a heat supply network, an air network module, a steam-water heat exchanger, a water-water heat exchanger, an electric boiler and other equipment modules.
And 4, step 4: and (5) finding that the system has the air network fault, and performing fault state simulation.
And 5: the air network can still be simulated under the condition that the output of a certain node is limited, so that the parameter change condition of an air network module is solved firstly;
step 6: transmitting the parameters of the gas network to a gas-heat and gas-electricity coupling equipment module;
and 7: simulating the gas-heat coupling equipment module, and transmitting the output parameters to the heat supply network module;
and 8: carrying out iterative calculation in a heat supply network module;
and step 9: transmitting the parameters of the heat supply network to an electric heating coupling module;
step 10: simulating the electric heating coupling module, and transmitting the output parameters to the power grid module;
step 11: simulating the power grid, and transmitting the output parameters to the gas-electric coupling module;
step 12: comparing the parameter with the gas-electric coupling module parameter in the step 6, observing whether a convergence condition is met, if so, ending the system simulation, otherwise, iterating the gas-electric coupling parameter to obtain a new gas-heat coupling parameter, and performing a step 8;
the dynamic simulation step length of each module under the air network fault is consistent with the conventional state, and the simulation result is shown in fig. 8. 1510s are actually used in the simulation, the whole iteration is carried out for 8 times, and under the condition that the air supply capacity of the 2000-2600s air network is limited, although a certain load is cut off by the power network and the power generation pressure of the generator is reduced, the heat storage characteristic of the heat supply network pipeline needs to be utilized because the air consumption of the air network is reduced more, so that the heat supply power of the heat supply network is reduced and the power balance of the system is ensured. The graphs (a) and (b) are the temperature graphs of the supply water and the return water of the heat source nodes 27 and 28 respectively, and in order to meet the power balance of the system, the heat supply power of the two heat sources is reduced during the heat supply period. Graph (c) shows that the generator generated power is reduced due to the generator loss load and the electric boiler load drop. And (d) is the pressure change situation of each node of the air network.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited by the foregoing examples, which are provided to illustrate the principles of the invention, and that various changes and modifications may be made without departing from the spirit and scope of the invention, which is also intended to be covered by the appended claims. The scope of the invention is defined by the appended claims and equivalents thereof.
In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed.
Claims (7)
1. A method for event-driven simulation of an integrated energy system based on decomposition, the method comprising: the method comprises the steps of synthesizing an energy system decomposition method, an event driving type simulation strategy, a conventional state solving strategy, a solving strategy under heat network disturbance, a solving strategy under air network disturbance and a solving strategy under power network disturbance;
the comprehensive energy system decomposition method divides a system into different modules according to different network types and equipment characteristics, then obtains corresponding module solving sequences according to different disturbance types by adopting a time-driven strategy, can divide the disturbance into heat supply network disturbance, air network disturbance and power grid disturbance according to the disturbance occurrence position, and can cause the known quantity of system load and different changes of a pipe network topological structure due to the disturbance at different positions, so that the corresponding solving sequences are obtained.
2. The method for event-driven simulation of an integrated energy system based on decomposition method according to claim 1, wherein the integrated energy system decomposition method comprises the following specific steps:
step 1: determining the composition of a comprehensive energy system;
step 2: inputting network parameters, initial values and equipment operation modes;
and step 3: based on step 1 and step 2, the system is divided into different modules according to type, characteristics and operation mode, such as: heat supply network, gas network, power network module, etc.;
and 4, step 4: after the single module is solved, the parameters are transmitted to the next module.
3. The method for event-driven simulation of an integrated energy system based on decomposition method according to claim 1, wherein the normal state solution strategy is: the system does not encounter disturbance and runs stably, at the moment, the system decomposition and solution has no hard sequence requirement, the coupling relation of the system needs to be deeply analyzed, the module with the weakest coupling and the lowest relevance with other modules is searched, and the solution is gradually carried out from the position.
4. The method for integrated energy system event-driven simulation based on decomposition method according to claim 1, wherein the event-driven simulation refers to: determining a module solving sequence according to a disturbance condition of a system in a simulation period and following a certain mechanism;
when the system is always in a conventional state, the simulation precision and speed can be ensured by solving the strategy according to the conventional state; if disturbance exists in the system, the system can be divided into different types of disturbance according to different disturbance positions, when the different types of disturbance occur, the disturbance network module needs to be solved preferentially as much as possible, then the coupling relation between the known quantity and the system module after the module solution is considered, a proper solving mode is selected according to the coupling relation, and if the disturbance network still has unknown quantity under the fault, the network where the unknown quantity exists is solved first, so that the solving priority of the fault network is ensured.
5. The method for event-driven simulation of an integrated energy system based on decomposition method according to claim 4, wherein solving strategy under heat supply network disturbance refers to: when the system generates heat supply network disturbance, whether the heat supply network can be directly solved or not is judged firstly based on the event driving type simulation strategy, the heat supply network disturbance can be directly simulated and solved generally, the coupling strength of other networks is analyzed after the solution, for an energy system with electricity, gas and heat networks, the coupling strength of the gas network is generally strong, the power grid can be solved firstly, and then the gas network is calculated
6. The method for event-driven simulation of an integrated energy system based on decomposition method according to claim 4, wherein solving strategy under air network disturbance refers to: when the air network is disturbed, whether the air network can be directly solved or not is judged firstly based on the event driving type simulation strategy, the air network can be directly calculated under common faults, and the coupling of the electric network and the heat network is equivalent after the working condition of the air network is obtained, so that the decomposition iteration solution through the thermoelectric network can be considered, and the high-precision calculation of the system is realized, and the specific steps are as follows:
step 1, determining the structure and parameters of a comprehensive energy system;
step 2, determining the specific disturbance condition of the air network;
step 3, solving the dynamic working condition of the air network, and transmitting parameters to other modules;
step 4, after the known quantity is determined, assuming an initial value of an iterative heat supply network;
step 5, solving the working condition of the heat supply network and transmitting the parameters to the power grid;
step 6, solving the working condition of the power grid, and transmitting parameters to the heat supply network;
step 7, comparing whether the heat supply network meets the convergence condition, if so, ending the simulation, otherwise, calculating a new iteration value, and performing the step 4;
7. the method for event-driven simulation of an integrated energy system based on decomposition method according to claim 4, wherein solving strategy under power grid disturbance comprises: when the system is disturbed by the power grid, whether the power grid can be directly solved or not is judged firstly based on the event-driven simulation strategy, and the power grid generally bears the task of supplying energy to other networks, so that direct calculation cannot be carried out, the network serving as the electric load is searched for to be solved, the power grid calculation is carried out after the solution is completed, and the remaining networks are calculated if the network which is not solved still exists.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110893014.0A CN113704977A (en) | 2021-08-04 | 2021-08-04 | Decomposition method based integrated energy system event driving type simulation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110893014.0A CN113704977A (en) | 2021-08-04 | 2021-08-04 | Decomposition method based integrated energy system event driving type simulation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113704977A true CN113704977A (en) | 2021-11-26 |
Family
ID=78651589
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110893014.0A Pending CN113704977A (en) | 2021-08-04 | 2021-08-04 | Decomposition method based integrated energy system event driving type simulation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113704977A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114818274A (en) * | 2022-03-31 | 2022-07-29 | 东南大学 | Variable-step-size dynamic simulation method for comprehensive energy system |
| CN115659680A (en) * | 2022-11-03 | 2023-01-31 | 国网江苏省电力有限公司电力科学研究院 | Method for decomposition and variable step size dynamic simulation of large thermoelectric coupling system |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107291990A (en) * | 2017-05-24 | 2017-10-24 | 河海大学 | Energy stream emulation mode based on electrical interconnection integrated energy system transient Model |
| CN108596453A (en) * | 2018-04-10 | 2018-09-28 | 山东大学 | Consider integrated energy system Optimization Scheduling and the system a few days ago of network dynamics |
| CN108717598A (en) * | 2018-04-25 | 2018-10-30 | 南京工程学院 | Index calculating method temporarily drops in the network voltage of electric-gas interconnection integrated energy system |
| CN109190785A (en) * | 2018-07-06 | 2019-01-11 | 东南大学 | A kind of electro thermal coupling integrated energy system running optimizatin method |
| CN110263435A (en) * | 2019-06-20 | 2019-09-20 | 燕山大学 | Dual-layer optimization fault recovery method based on electric-gas coupling integrated energy system |
| CN110705066A (en) * | 2019-09-20 | 2020-01-17 | 天津大学 | Projection integral-based dynamic simulation method for integrated energy system of gas-electricity coupling park |
| CN111082417A (en) * | 2019-12-01 | 2020-04-28 | 国网辽宁省电力有限公司经济技术研究院 | State estimation method based on comprehensive energy system electric and heat combined network |
| CN111832161A (en) * | 2020-06-29 | 2020-10-27 | 山东大学 | Real-time simulation method and system for comprehensive energy system |
| CN112182905A (en) * | 2020-10-16 | 2021-01-05 | 北京科东电力控制系统有限责任公司 | Heat supply pipe network simulation method and device for comprehensive energy system |
| CN112364493A (en) * | 2020-10-30 | 2021-02-12 | 国网福建省电力有限公司厦门供电公司 | Intelligent scheduling method and device for urban comprehensive energy system |
-
2021
- 2021-08-04 CN CN202110893014.0A patent/CN113704977A/en active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107291990A (en) * | 2017-05-24 | 2017-10-24 | 河海大学 | Energy stream emulation mode based on electrical interconnection integrated energy system transient Model |
| CN108596453A (en) * | 2018-04-10 | 2018-09-28 | 山东大学 | Consider integrated energy system Optimization Scheduling and the system a few days ago of network dynamics |
| CN108717598A (en) * | 2018-04-25 | 2018-10-30 | 南京工程学院 | Index calculating method temporarily drops in the network voltage of electric-gas interconnection integrated energy system |
| CN109190785A (en) * | 2018-07-06 | 2019-01-11 | 东南大学 | A kind of electro thermal coupling integrated energy system running optimizatin method |
| CN110263435A (en) * | 2019-06-20 | 2019-09-20 | 燕山大学 | Dual-layer optimization fault recovery method based on electric-gas coupling integrated energy system |
| CN110705066A (en) * | 2019-09-20 | 2020-01-17 | 天津大学 | Projection integral-based dynamic simulation method for integrated energy system of gas-electricity coupling park |
| CN111082417A (en) * | 2019-12-01 | 2020-04-28 | 国网辽宁省电力有限公司经济技术研究院 | State estimation method based on comprehensive energy system electric and heat combined network |
| CN111832161A (en) * | 2020-06-29 | 2020-10-27 | 山东大学 | Real-time simulation method and system for comprehensive energy system |
| CN112182905A (en) * | 2020-10-16 | 2021-01-05 | 北京科东电力控制系统有限责任公司 | Heat supply pipe network simulation method and device for comprehensive energy system |
| CN112364493A (en) * | 2020-10-30 | 2021-02-12 | 国网福建省电力有限公司厦门供电公司 | Intelligent scheduling method and device for urban comprehensive energy system |
Non-Patent Citations (1)
| Title |
|---|
| AOBO GUAN，等: "An optimal step-size simulation framework for large-scale heat-electric integrated energy system considering fault states", ELECTRIC POWERSYSTEMSRESEARCH, vol. 223, 31 December 2023 (2023-12-31), pages 1 - 10 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114818274A (en) * | 2022-03-31 | 2022-07-29 | 东南大学 | Variable-step-size dynamic simulation method for comprehensive energy system |
| CN115659680A (en) * | 2022-11-03 | 2023-01-31 | 国网江苏省电力有限公司电力科学研究院 | Method for decomposition and variable step size dynamic simulation of large thermoelectric coupling system |
| CN115659680B (en) * | 2022-11-03 | 2023-11-17 | 国网江苏省电力有限公司电力科学研究院 | Large-scale thermoelectric coupling system decomposition and variable step dynamic simulation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108921727B (en) | Regional comprehensive energy system reliability assessment method considering thermal load dynamic characteristics | |
| CN111008468B (en) | Test method and test system of comprehensive energy management system | |
| CN106230028B (en) | A kind of Multipurpose Optimal Method of wind-powered electricity generation-water-storage association system | |
| CN113704977A (en) | Decomposition method based integrated energy system event driving type simulation method | |
| CN107516895B (en) | Power distribution network rapid simulation method, device, storage medium and its computer equipment | |
| CN104734147A (en) | Probability energy flow analysis method for integrated energy system (IES) | |
| CN111313429A (en) | Reliability assessment method and system for comprehensive energy system | |
| CN109524959B (en) | Power generation and transmission system abundance assessment method considering natural gas network fault | |
| Maihemuti et al. | Dynamic security and stability region under different renewable energy permeability in IENGS system | |
| CN111401476B (en) | Transient state safety evaluation method based on boundary region importance sampling and kernel vector machine | |
| CN103149840B (en) | Semanteme service combination method based on dynamic planning | |
| CN110266046B (en) | Electric heating micro-grid topology comprehensive diagnosis method and system based on complex network | |
| CN113536591A (en) | Variable-step-size dynamic simulation method for comprehensive energy system | |
| Huang et al. | A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale | |
| CN110705066B (en) | Dynamic simulation method for gas-electricity coupling park comprehensive energy system based on projection integration | |
| CN110429591B (en) | Power transmission network utilization rate evaluation method based on power system time sequence coupling | |
| Zhang et al. | Risk assessment of offshore micro integrated energy system based on fluid mosaic model | |
| Guo et al. | A thermal response time ahead energy demand prediction strategy for building heating system using machine learning methods | |
| Groumpos et al. | An overview of fuzzy cognitive maps for energy efficiency in intelligent buildings | |
| CN112670997B (en) | Electric heating energy system time sequence probability power flow calculation method considering photovoltaic uncertainty | |
| Tang et al. | Multi rate dynamic hybrid simulation of integrated energy system | |
| Aryan et al. | Simulation support for explainable cyber-physical energy systems | |
| CN106027301A (en) | Method for searching key failure nodes in power heterogeneous communication network system | |
| CN115935711B (en) | Multi-energy complementary active distribution network system reliability assessment method based on graph theory | |
| CN109740976A (en) | A kind of integrated energy system cascading failure source tracing method based on Markov process |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |