CN111400937A - Load flow calculation method of comprehensive energy system - Google Patents
Load flow calculation method of comprehensive energy system Download PDFInfo
- Publication number
- CN111400937A CN111400937A CN202010354139.1A CN202010354139A CN111400937A CN 111400937 A CN111400937 A CN 111400937A CN 202010354139 A CN202010354139 A CN 202010354139A CN 111400937 A CN111400937 A CN 111400937A
- Authority
- CN
- China
- Prior art keywords
- power
- node
- energy
- power flow
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 19
- 230000008878 coupling Effects 0.000 claims abstract description 18
- 238000010168 coupling process Methods 0.000 claims abstract description 18
- 238000005859 coupling reaction Methods 0.000 claims abstract description 18
- 238000005485 electric heating Methods 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 238000010438 heat treatment Methods 0.000 claims description 7
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims 1
- 239000002994 raw material Substances 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 239000007789 gas Substances 0.000 description 39
- 238000004146 energy storage Methods 0.000 description 11
- 238000003860 storage Methods 0.000 description 8
- 238000010248 power generation Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 239000002737 fuel gas Substances 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Images
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a load flow calculation method of an integrated energy system, which comprises the following steps that a processing module obtains load flow models of all energy subsystems; the acquisition module acquires an energy system structure and establishes an electric heating and gas coupling equipment model through the processing module; constructing a comprehensive energy system power flow model covering each energy subsystem by using the processing module; the processing module obtains a parameter result of the energy system node through an expansibility Newton iteration method algorithm; and judging whether the parameter result is converged or not, and judging the state of the comprehensive energy system according to the output result. The invention has the beneficial effects that: the unified solution of the comprehensive energy system power flow is realized, the power flows of different energy systems are mutually associated and linked under the connection of the coupling equipment model, the relation between various energy supply and requirements can be correctly reflected, the support is provided for effectively arranging the output of each unit, maintaining the integral supply and demand balance of the system, and the safe and efficient operation of the comprehensive energy system is ensured.
Description
Technical Field
The invention relates to the technical field of comprehensive energy, in particular to a load flow calculation method of a comprehensive energy system.
Background
In recent years, under the double pressure of energy safety and resource shortage, it has become a consensus in the energy community by promoting the joint supply between various energy sources and improving the flexibility between systems. The demands of Integrated Energy Systems (IES) research are proposing and developing profound background appeals reflecting many driving forces such as environment, economy, society, technology and policy. The unified modeling of the comprehensive energy system is used as the unified description of different energy systems, is the research basis of the planning, scheduling, control and interaction of the multi-energy system, and describes the operation and complementary conversion characteristics of each energy system. For the power flow calculation of an energy system, deep research is carried out at home and abroad on the power flow calculation of a single system, for example, an expansibility Newton algorithm is mainly adopted for calculation in a power system, or combined operation is carried out by combining with a Gaussian method, and research related to combined modeling and unified solution of a plurality of systems is rarely carried out, or only aiming at an electric heating coupling system, a gas system is not considered, and the coverage is not wide enough. Some existing researches are aimed at modeling of a specific system and are not provided, when new elements are added into the system, the system cannot be directly accessed into an original load flow calculation model, and the universality is not achieved.
Disclosure of Invention
This section is for the purpose of summarizing some aspects of embodiments of the invention and to briefly introduce some preferred embodiments. In this section, as well as in the abstract and the title of the invention of this application, simplifications or omissions may be made to avoid obscuring the purpose of the section, the abstract and the title, and such simplifications or omissions are not intended to limit the scope of the invention.
The present invention has been made in view of the above-mentioned conventional problems.
Therefore, the technical problem solved by the invention is as follows: the comprehensive energy system load flow calculation method can deal with the addition of new nodes or elements of an energy system in calculation, and has universality and more accurate calculation results.
In order to solve the technical problems, the invention provides the following technical scheme: a load flow calculation method of an integrated energy system comprises the steps that a processing module obtains load flow models of all energy subsystems; the acquisition module acquires an energy system structure and establishes an electric heating and gas coupling equipment model through the processing module; constructing a comprehensive energy system power flow model covering each energy subsystem by using the processing module; the processing module obtains a parameter result of the energy system node through an expansibility Newton iteration method algorithm; and judging whether the parameter result is converged or not, and judging the state of the comprehensive energy system according to the output result.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the power flow model of the energy subsystem comprises a power system power flow model, a gas system power flow model and a thermal system power flow model.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the obtaining of the power system load flow model further comprises the steps of determining the positions and node types of all nodes in the power system, including a newly added unit or equipment; the power flow model of the power subsystem is represented by an alternating current power flow model, the node power of the power flow model is as follows,
where Real represents the Real part, Imag represents the imaginary part, P, Q are the active power and reactive power vectors of the node, Y is the node admittance matrix, and U is the node voltage phasor.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the obtaining of the gas system power flow model further comprises obtaining a gas system network topological structure; calculating a node-pipeline incidence matrix and a loop-pipeline incidence matrix according to the network structure; and obtaining a gas system power flow model according to the incidence matrix.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the obtaining of the thermodynamic system power flow model further comprises obtaining a thermodynamic system pipe network topological structure; calculating a node-pipeline incidence matrix and a loop-pipeline incidence matrix according to a pipe network topological structure; and obtaining a thermodynamic system power flow model according to the incidence matrix.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the electric heating and gas coupling equipment is used as coupling link equipment and comprises a CHP combined supply unit, a gas boiler, an electric boiler and a heat pump, and the established electric heating and gas coupling equipment model is an equipment model based on energy equivalent balance.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: and if the elements in the objective function matrix accord with the set accuracy, outputting a load flow calculation result, otherwise, performing iteration by using an iterative formula until the elements in the matrix accord with the set accuracy.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the calculation formula of the power system flow model is as follows,
wherein, Δ F represents the difference between each energy supply parameter and the known parameter, P is the active power of the power system, Q is the reactive power of the power system, Φ is the node pressure, P is the thermal power, T is the powersFor heating temperature, TrThe return water temperature is set, and M is the pipeline flow.
As a preferable scheme of the integrated energy system load flow calculation method of the present invention, wherein: the iterative formula is calculated as follows,
and the delta F is a difference result of each energy supply parameter and a known parameter obtained by the power flow calculation model.
The invention has the beneficial effects that: the load flow calculation method provided by the invention realizes the unified solution of the load flow of the comprehensive energy system, and the load flows of different energy systems are mutually associated and linked under the connection of the coupling equipment model, so that the relation between various energy supply and demands can be correctly reflected, the support is provided for effectively arranging the output of each unit, maintaining the integral supply and demand balance of the system, and the safe and efficient operation of the comprehensive energy system is ensured.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise. Wherein:
fig. 1 is a schematic overall flow chart of a method for calculating a power flow of an integrated energy system according to a first embodiment of the present invention.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, specific embodiments accompanied with figures are described in detail below, and it is apparent that the described embodiments are a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making creative efforts based on the embodiments of the present invention, shall fall within the protection scope of the present invention.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways than those specifically described and will be readily apparent to those of ordinary skill in the art without departing from the spirit of the present invention, and therefore the present invention is not limited to the specific embodiments disclosed below.
Furthermore, reference herein to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one implementation of the invention. The appearances of the phrase "in one embodiment" in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments.
The present invention will be described in detail with reference to the drawings, wherein the cross-sectional views illustrating the structure of the device are not enlarged partially in general scale for convenience of illustration, and the drawings are only exemplary and should not be construed as limiting the scope of the present invention. In addition, the three-dimensional dimensions of length, width and depth should be included in the actual fabrication.
Meanwhile, in the description of the present invention, it should be noted that the terms "upper, lower, inner and outer" and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of describing the present invention and simplifying the description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation and operate, and thus, cannot be construed as limiting the present invention. Furthermore, the terms first, second, or third are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
The terms "mounted, connected and connected" in the present invention are to be understood broadly, unless otherwise explicitly specified or limited, for example: can be fixedly connected, detachably connected or integrally connected; they may be mechanically, electrically, or directly connected, or indirectly connected through intervening media, or may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Example 1
Referring to the schematic diagram of fig. 1, a flow chart diagram illustrating a method for calculating a power flow of an integrated energy system according to the embodiment specifically includes the following steps,
s1: the processing module acquires a power flow model of each energy subsystem; the power flow model of the energy subsystem comprises a power system power flow model, a gas system power flow model and a thermal system power flow model;
specifically, the obtaining of the power flow model of the power system further comprises determining the position and the type of each node in the power system, including a newly added unit or equipment; because the equipment types of the power system in the integrated energy system are richer, and uncertainty or mobility exists in units and equipment accessed by the system, before load flow calculation, the node types of each node need to be determined, and for power load flow, the node types mainly include 3 node types: and the balance nodes, the PV nodes and the PQ nodes are used for selecting typical equipment and units on a source-load-storage side of an energy system respectively and analyzing the node types of the typical equipment and the typical units in load flow calculation.
The source side is an available energy generating set and comprises a wind power generation node and a photovoltaic power generation node. The node type of the wind power generation node in the power system load flow calculation is related to the type of the generator. The wind driven generator mainly comprises four types of constant-speed type wind driven generators, slip type wind driven generators, synchronous direct-drive type wind driven generators and double-fed wind driven generators, wherein the constant-speed type wind driven generators and the slip type wind driven generators are regarded as P-Q (V) nodes in power flow calculation, and the synchronous direct-drive type wind driven generators and the double-fed wind driven generators are regarded as PQ nodes. The P-Q (V) node means P deterministic, V deterministic, Q is limited to P and V. When the wind farm has the capability of providing real-time reactive compensation and the reactive demand exceeds the compensation capability during the operation of the node, the reactive power is out of limit, at the moment, the wind power generation node is treated as a PV node, and at the moment, the PV node is converted into a PQ node.
The photovoltaic power generation node is connected to a power grid through an inverter, in the embodiment, the node type is divided according to the type of the adopted inverter, and the inverter includes a current control type inverter and a voltage control type inverter. The current control type is a main current and outputs stable current, so that the current control type can be used as a PI node of known active power and current; the voltage controlled inverter serves as the PV node. Wherein, the PI node type can be converted into PQ node through calculation, the calculation formula is as follows,
wherein Q isk+1The reactive power values of the PI nodes for the (k + 1) th iteration I, P are respectivelyConstant current amplitude and active power value, V, of PI nodekThe node voltage amplitude obtained for the kth iteration.
The storage refers to an energy storage node, and the embodiment comprises an energy storage battery and an electric automobile. The control mode of the energy storage battery connected to the power grid is basically the same as that of the photovoltaic power generation system, the energy storage battery is connected to the power grid through the inverter, the energy storage battery is used as a PI node for processing when the current is adopted to control the inverter, and the energy storage battery is used as a PV node for processing when the voltage is adopted to control the inverter. The difference is that the energy storage battery has bidirectional property, and can be used as a source side to supply power to a power grid and can also be used as a charge side to store electric energy. When the energy storage battery is used as a power supply, the battery discharges to the power grid and feeds back the electric energy to the power grid; when the energy storage battery is used as a charge side to store energy, the energy flows from the grid side to the battery.
In the V1G mode, the electric vehicle handles as a load. In the V2G mode, the electric vehicle, like the energy storage system, can be used as both a load of the power grid to consume electric energy and a backup power supply to supply electric energy to the power grid. When the electric vehicle is used as a load, the electric vehicle is treated as a PQ node; when the electric automobile is used as a standby power supply, the electric automobile has the same processing method as an energy storage system and is processed as a PI or PV node.
The load refers to an adjustable load node, and in order to improve flexibility and economy of system operation, demand response on a load side and application of emerging technologies such as a virtual power plant, the load needs to be classified so as to judge which loads have adjusting capacity and adjusting space. In the embodiment, the load with adjustable capacity is classified as adjustable load, and such load is treated as PQ node in load flow calculation as the common load, but the node power is increased or decreased accordingly.
S2: the acquisition module acquires an energy system structure and establishes an electric heating and gas coupling equipment model through the processing module;
specifically, the established electric heating and gas coupling equipment model further comprises a power subsystem power flow model, a gas system power flow model and a gas system power flow model. The method comprises the steps that a power system load flow model is obtained, wherein the power system load flow model also comprises the steps of determining the positions and node types of all nodes in the power system, including a newly added unit or equipment; the power flow model of the power subsystem is represented by an alternating current power flow model, the node power of the power flow model is as follows,
where Real represents the Real part, Imag represents the imaginary part, P, Q are the active power and reactive power vectors of the node, Y is the node admittance matrix, and U is the node voltage phasor.
The obtaining of the gas system flow model further comprises,
acquiring a network topological structure of a gas system; the network topology structure comprises a ring pipe network and a branch pipe network.
Computing a node-pipeline incidence matrix A according to a network structure1And a loop-pipe association matrix B1(ii) a Node-pipe incidence matrix A1The formula (2) is calculated as follows,
wherein A is1When +1 is taken, the node i is the end point of the pipeline j, when 0 is taken, the node i is irrelevant to the pipeline j, and when-1 is taken, the node i is the starting point of the pipeline j.
Loop-pipe association matrix B1The formula (2) is calculated as follows,
wherein, B1When +1 is taken, the direction of the branch pipeline is the same as that of the loop, when 0 is taken, the branch pipeline is not in the loop, and when-1 is taken, the direction of the branch pipeline is opposite to that of the loop.
Specifically, the gas network needs to satisfy a hydraulic model, and the hydraulic model needs to satisfy a flow continuity equation, a pressure drop equation and an energy equation, wherein the steady-state flow M of the natural gas pipeline LLCan be expressed as a number of times,
wherein, KrIs the constant of the pipeline and is,
is the pressure drop, S, of the line LijIndicates the flow direction of the gas, when pi>pjIf yes, taking +1, otherwise, taking-1.
Specifically, the flow of each pipeline at each node satisfies a node flow continuity equation, i.e., the flow of the node inflow is equal to the flow of the node outflow, and satisfies the following formula,
AgM＝mg
wherein A isgRepresenting a node-pipe correlation matrix in the gas supply network, M representing the flow rate of each gas pipe, MgIndicating the traffic flowing out of each node, i.e. the traffic required by the user.
The pressure drop equation is satisfied in the branched pipe network, namely the calculated total pressure drop is within the allowable total pressure drop range, the following formula is satisfied,
wherein the content of the first and second substances,
is the voltage drop between the user nodes i, j,
and
representing the pressure between nodes i, j, respectively.
In the annular pipe network, the sum of the pressure losses of the fuel gas flowing in the pipeline is 0 because of the requirement of an energy equation, and the following formula is satisfied,
wherein, BgFor a loop-pipe correlation matrix in a gas supply network,
is the voltage drop between the user nodes i, j.
The obtaining of the thermodynamic system power flow model further comprises the following steps,
acquiring a topological structure of a thermodynamic system pipe network; the network topology structure comprises a ring pipe network and a branch pipe network.
Calculating a node-pipeline incidence matrix and a loop-pipeline incidence matrix according to a pipe network topological structure; calculating a node-pipeline incidence matrix A according to a pipe network structure of a thermodynamic system2And a loop-pipe association matrix B2The calculation method is the same as the calculation method of the gas system node-pipeline incidence matrix A1And a loop-pipe association matrix B1And (4) calculating.
And obtaining a thermodynamic system power flow model according to the incidence matrix. The thermodynamic network needs to satisfy a hydraulic model and a thermodynamic model, the hydraulic model is the same as that in the gas network, and the thermodynamic model needs to satisfy a temperature heat calculation equation and a path heat loss equation. Wherein the temperature parameter of the thermodynamic system comprises a heating temperature TsOutput temperature ToAnd return water temperature TrThree parameters, heating temperature TsIndicating the temperature of the hot water injected into the user node, the output temperature ToIndicating the temperature of hot water flowing out of the user node, the return water temperature TrIndicating the temperature before flowing out of the customer node into other pipes.
The relationship between the thermal power of each node in the thermal power subsystem and the temperature and the flow is as follows,
Φ＝Cpmq(Ts-T0)
where Φ is the thermal power of the node, mqFor traffic flowing into the node, CpIs the specific heat capacity of water, TsFor heating temperature, ToIs the output temperature.
The temperature relationship between the beginning and the end of the pipeline is as follows,
wherein, TendIs the temperature at the end of the pipe, TstartIs the temperature at the beginning of the pipe, TaLambda is the heat transfer coefficient of the pipe, which is the length of the pipe, for ambient temperature. Line of T'end＝Tend-Ta，T′start＝Tstart-Ta，
The temperature relationship between the beginning and the end of the pipe can be simplified to,
T′end＝T′startψ
for a node with multiple pipe injections, the relationship between the input temperature and the output temperature of the hot water is as follows,
(∑mout)Tout＝∑minTin
wherein m isoutFlow of water in a pipe to an outflow node, ToutTemperature of the outflowing water, minFor flow into the pipes of the node, TinIs the temperature of the incoming water.
The electric heating and gas coupling equipment is used as coupling link equipment and comprises a CHP combined supply unit, a gas boiler, an electric boiler and a heat pump, and the established electric heating and gas coupling equipment model is an equipment model based on energy equivalent balance. Wherein, the relationship among the gas consumption, the thermal power and the electric power of the CHP combined supply unit is as follows,
wherein, FinFor gas consumption of gas-fired units, ηeTo the generating efficiency of the unit, PCHPAnd phiCHPElectric and thermal power, C, of the unit, respectivelymIs the thermoelectric ratio of the unit.
The exhaust-heat boiler in the CHP unit represented by the above formula is not afterburning-free, and if the exhaust-heat boiler with afterburning is adopted, a model of an afterburning-type exhaust-heat boiler should be added, as shown below,
QGB＝kfGBηGB+QCHPηy
wherein Q isGBFor total heating power of the boiler, ηGBFor the thermal conversion efficiency of the gas, fGBThe unit is the gas consumption of the natural gas, k represents whether the afterburning mode is started or not, k is 1 when the afterburning is started, k is 0 when the afterburning is started, ηyIndicating the thermal conversion efficiency of the waste heat.
S3: constructing a comprehensive energy system power flow model covering each energy subsystem by using the processing module;
and if the elements in the objective function matrix accord with the set accuracy, outputting a load flow calculation result, otherwise, performing iteration by using an iterative formula until the elements in the matrix accord with the set accuracy.
Specifically, the calculation formula of the energy system flow model is as follows,
wherein, Δ F represents the difference between each calculated energy supply parameter and a known parameter, the superscript sp represents the established, P is the active power of the power system, Q is the reactive power of the power system, Φ is the node pressure, P is the thermal power, TsFor heating temperature, TrThe return water temperature is set, and M is the pipeline flow.
For power flow analysis, vectors in active power mismatch are specified. Whereas for integrated electro-hydro-thermal calculations, the elements of the error deltaf vector are determined by the thermal power generated at the thermal relaxation node, and are expressed as a function of the thermal network,
is the state quantity of the system.
S4: the processing module obtains a parameter result of the energy system node through an expansibility Newton iteration method algorithm; in the embodiment, the expansibility Newton iteration algorithm is expanded to solve the energy source flow of the comprehensive energy system, the calculation formula of the iteration formula is,
wherein, Δ F is the difference result between each energy supply parameter and the known parameter obtained by the power flow calculation model, so that Δ Fe＝[ΔP,ΔQ]T、
And Δ FeΔ f represents the deviation amount related to electricity and heat, respectively, and xe＝[θ、|V|]T、xh＝[m,(T′s,load,T′r,load)]TAnd
respectively, representing the state quantities related to electricity, heat, gas, where the jacobian matrix J can be expressed as,
wherein, Jee、JhhAnd JggRespectively representing the relation between the self tidal current and the self state quantity of the independent electric, thermal and gas systems, wherein the expression is the partial derivative of the subsystem energy flow model relative to the state parameter of the subsystem, JehAnd JegRespectively representing the influence of the energy flows of the thermodynamic system and the gas system on the power flow of the power system, JheAnd JhgRespectively representing the influence of the energy flows of the power system and the gas system on the state of the thermodynamic system, JgeAnd JghRespectively shows the influence of the energy flow of the electric power system and the thermodynamic system on the state of the gas system, and for the comprehensive energy system based on CHP combined supply, when the supply and demand balance of the internal nodes of the gas system changes, the balance nodes, namely the gas supply sources, can stabilize, so JegAnd JhgIs zero. For the systemFor the system operation, 2 cases of power-on-heat and power-on-heat need to be distinguished to determine the mutual influence between subsystems.
S5: and judging whether the parameter result is converged or not, and judging the state of the comprehensive energy system according to the output result. Specifically, the difference precision of the data is set according to the requirement of the system on the energy supply quality, when the parameter result meets the precision requirement, the parameter result is judged to be convergence and a result is output, and the state of the comprehensive energy system is normal; otherwise, the state of the comprehensive energy system is considered to be abnormal, and technicians are required to adjust the system and update the system for processing and calculation until the system is normal.
Scene one:
the comprehensive energy system load flow calculation method provided by the embodiment can be used for integrating a plurality of energy systems to perform combined modeling and unified solution in practical application, and when new elements are added into the energy systems, the new elements can be accessed into the original load flow calculation model, so that the comprehensive energy system load flow calculation method has higher universality in practical application.
In the traditional method, when the energy system is solved, only an electric heating coupling system is usually used, other energy systems such as a gas system and the like are not considered, load flow calculation in a comprehensive energy system is usually carried out at a cost, or modeling is only carried out for a specific system, and a general method is not provided, so that in practical application, the method is difficult to adapt to change of nodes in the energy system, and the flexibility is poor.
In order to verify the beneficial effect of the comprehensive energy system load flow calculation method provided by the embodiment in the practical application compared with the traditional energy system load flow calculation method, the following experiments are carried out: selecting a growing intelligent town in a certain planning of Guizhou as an example, setting that only small enterprises exist at the initial stage and the electric load is small; with the planning according to the small town, new enterprises with larger power loads can live in successively, the power loads are larger and larger, the cold and hot load requirements of the air conditioner are increased, and energy is supplied by adopting a mode of increasing distributed energy sources. The small towns depend on a local power distribution network, the small towns are divided into 3 different scenes according to the near-term, medium-term and long-term energy planning of the small towns, and a roof photovoltaic system is installed in the scene 1; scene 2 is used for increasing the power of a roof photovoltaic power generation system, and meanwhile, a cold, heat and electricity triple power supply system formed by distributed gas turbines is built for carrying out combined power supply; and in the scene 3, the capacity of the energy station unit is expanded to meet the subsequent increase of load demand. The results obtained by performing calculations using the method provided in this example are shown in table 1 below,
table 1: energy parameter table of system under different periods
It can be seen that the method provided by the embodiment is used for calculating the comprehensive energy system under different conditions, the self equipment change of the energy system can be simply adjusted, the complexity of the calculation process is reduced, and the method has strong universality for the phenomenon that the comprehensive energy system is changeable and complex in condition in modern practical application.
It should be recognized that embodiments of the present invention can be realized and implemented by computer hardware, a combination of hardware and software, or by computer instructions stored in a non-transitory computer readable memory. The methods may be implemented in a computer program using standard programming techniques, including a non-transitory computer-readable storage medium configured with the computer program, where the storage medium so configured causes a computer to operate in a specific and predefined manner, according to the methods and figures described in the detailed description. Each program may be implemented in a high level procedural or object oriented programming language to communicate with a computer system. However, the program(s) can be implemented in assembly or machine language, if desired. In any case, the language may be a compiled or interpreted language. Furthermore, the program can be run on a programmed application specific integrated circuit for this purpose.
Further, the operations of processes described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The processes described herein (or variations and/or combinations thereof) may be performed under the control of one or more computer systems configured with executable instructions, and may be implemented as code (e.g., executable instructions, one or more computer programs, or one or more applications) collectively executed on one or more processors, by hardware, or combinations thereof. The computer program includes a plurality of instructions executable by one or more processors.
Further, the method may be implemented in any type of computing platform operatively connected to a suitable interface, including but not limited to a personal computer, mini computer, mainframe, workstation, networked or distributed computing environment, separate or integrated computer platform, or in communication with a charged particle tool or other imaging device, and the like. Aspects of the invention may be embodied in machine-readable code stored on a non-transitory storage medium or device, whether removable or integrated into a computing platform, such as a hard disk, optically read and/or write storage medium, RAM, ROM, or the like, such that it may be read by a programmable computer, which when read by the storage medium or device, is operative to configure and operate the computer to perform the procedures described herein. Further, the machine-readable code, or portions thereof, may be transmitted over a wired or wireless network. The invention described herein includes these and other different types of non-transitory computer-readable storage media when such media include instructions or programs that implement the steps described above in conjunction with a microprocessor or other data processor. The invention also includes the computer itself when programmed according to the methods and techniques described herein. A computer program can be applied to input data to perform the functions described herein to transform the input data to generate output data that is stored to non-volatile memory. The output information may also be applied to one or more output devices, such as a display. In a preferred embodiment of the invention, the transformed data represents physical and tangible objects, including particular visual depictions of physical and tangible objects produced on a display.
As used in this application, the terms "component," "module," "system," and the like are intended to refer to a computer-related entity, either hardware, firmware, a combination of hardware and software, or software in execution. For example, a component may be, but is not limited to being: a process running on a processor, an object, an executable, a thread of execution, a program, and/or a computer. By way of example, both an application running on a computing device and the computing device can be a component. One or more components can reside within a process and/or thread of execution and a component can be localized on one computer and/or distributed between two or more computers. In addition, these components can execute from various computer readable media having various data structures thereon. The components may communicate by way of local and/or remote processes such as in accordance with a signal having one or more data packets (e.g., data from one component interacting with another component in a local system, distributed system, and/or across a network such as the internet with other systems by way of the signal).
It should be noted that the above-mentioned embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, which should be covered by the claims of the present invention.
Claims (9)
1. A load flow calculation method of an integrated energy system is characterized by comprising the following steps: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
the processing module acquires a power flow model of each energy subsystem;
the acquisition module acquires an energy system structure and establishes an electric heating and gas coupling equipment model through the processing module;
constructing a comprehensive energy system power flow model covering each energy subsystem by using the processing module;
the processing module obtains a parameter result of the energy system node through an expansibility Newton iteration method algorithm;
and judging whether the parameter result is converged or not, and judging the state of the comprehensive energy system according to the output result.
2. The integrated energy system power flow calculation method of claim 1, wherein: the power flow model of the energy subsystem comprises a power system power flow model, a gas system power flow model and a thermal system power flow model.
3. The integrated energy system power flow calculation method of claim 2, wherein: the obtaining of the power system load flow model further comprises the steps of determining the positions and node types of all nodes in the power system, including a newly added unit or equipment; the power flow model of the power subsystem is represented by an alternating current power flow model, the node power of the power flow model is as follows,
where Real represents the Real part, Imag represents the imaginary part, P, Q are the active power and reactive power vectors of the node, Y is the node admittance matrix, and U is the node voltage phasor.
4. The integrated energy system power flow calculation method of claim 2, wherein: the obtaining of the gas system flow model further comprises,
acquiring a network topological structure of a gas system;
calculating a node-pipeline incidence matrix and a loop-pipeline incidence matrix according to the network structure;
and obtaining a gas system power flow model according to the incidence matrix.
5. The integrated energy system power flow calculation method of claim 2, wherein: the obtaining of the thermodynamic system power flow model further comprises,
acquiring a topological structure of a thermodynamic system pipe network;
calculating a node-pipeline incidence matrix and a loop-pipeline incidence matrix according to a pipe network topological structure;
and obtaining a thermodynamic system power flow model according to the incidence matrix.
6. The integrated energy system power flow calculation method of claim 2, wherein: the electric heating and gas coupling equipment is used as coupling link equipment and comprises a CHP combined supply unit, a gas boiler, an electric boiler and a heat pump, and the established electric heating and gas coupling equipment model is an equipment model based on energy equivalent balance.
7. The integrated energy system power flow calculation method of claim 2, wherein: and if the elements in the objective function matrix accord with the set accuracy, outputting a load flow calculation result, otherwise, performing iteration by using an iterative formula until the elements in the matrix accord with the set accuracy.
8. The integrated energy system power flow calculation method of claim 2, wherein: the calculation formula of the power system flow model is as follows,
wherein, Δ F represents the difference between each energy supply parameter and the known parameter, P is the active power of the power system, Q is the reactive power of the power system, Φ is the node pressure, P is the thermal power, T is the powersFor heating temperature, TrThe return water temperature is set, and M is the pipeline flow.
9. The integrated energy system power flow calculation method of claim 2, wherein: the iterative formula is calculated as follows,
and the delta F is a difference result of each energy supply parameter and a known parameter obtained by the power flow calculation model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010354139.1A CN111400937B (en) | 2020-04-29 | 2020-04-29 | Comprehensive energy system tide calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010354139.1A CN111400937B (en) | 2020-04-29 | 2020-04-29 | Comprehensive energy system tide calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111400937A true CN111400937A (en) | 2020-07-10 |
| CN111400937B CN111400937B (en) | 2023-12-01 |
Family
ID=71433778
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010354139.1A Active CN111400937B (en) | 2020-04-29 | 2020-04-29 | Comprehensive energy system tide calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111400937B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112488535A (en) * | 2020-12-01 | 2021-03-12 | 中国科学院地理科学与资源研究所 | Energy network security measure method, system and data platform |
| CN112994020A (en) * | 2021-03-31 | 2021-06-18 | 南京信息工程大学 | Multi-energy system load flow decomposition calculation method |
| CN113255105A (en) * | 2021-04-26 | 2021-08-13 | 上海电力大学 | Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure |
| CN115562029A (en) * | 2022-10-17 | 2023-01-03 | 杭州天然气有限公司 | Intelligent control method and system for natural gas turbine expansion generator set |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108875269A (en) * | 2018-07-09 | 2018-11-23 | 重庆大学 | A kind of multi-period energy flux computation method of electric-gas interacted system considering electric system multi-balancing machine and natural gas system slow motion step response |
| WO2019200662A1 (en) * | 2018-04-20 | 2019-10-24 | 东北大学 | Stability evaluation and static control method for electricity-heat-gas integrated energy system |
| CN110543661A (en) * | 2019-07-18 | 2019-12-06 | 国网江苏省电力有限公司 | correlation-considered probability energy flow calculation method for electricity-heat interconnection comprehensive energy system |
| CN110601185A (en) * | 2019-09-17 | 2019-12-20 | 武汉大学 | Unified power flow model and random matrix-based comprehensive energy system weak point identification method |
| CN110955954A (en) * | 2019-08-19 | 2020-04-03 | 天津大学 | Optimal load reduction method for layered decoupling electric and thermal comprehensive energy system |
-
2020
- 2020-04-29 CN CN202010354139.1A patent/CN111400937B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019200662A1 (en) * | 2018-04-20 | 2019-10-24 | 东北大学 | Stability evaluation and static control method for electricity-heat-gas integrated energy system |
| CN108875269A (en) * | 2018-07-09 | 2018-11-23 | 重庆大学 | A kind of multi-period energy flux computation method of electric-gas interacted system considering electric system multi-balancing machine and natural gas system slow motion step response |
| CN110543661A (en) * | 2019-07-18 | 2019-12-06 | 国网江苏省电力有限公司 | correlation-considered probability energy flow calculation method for electricity-heat interconnection comprehensive energy system |
| CN110955954A (en) * | 2019-08-19 | 2020-04-03 | 天津大学 | Optimal load reduction method for layered decoupling electric and thermal comprehensive energy system |
| CN110601185A (en) * | 2019-09-17 | 2019-12-20 | 武汉大学 | Unified power flow model and random matrix-based comprehensive energy system weak point identification method |
Non-Patent Citations (3)
| Title |
|---|
| 刘聪等: "《电/ 热/ 气综合能源系统混合潮流计算方法》", 《山东工业技术》 * |
| 王英瑞等: "电–热–气综合能源系统多能流计算方法" * |
| 罗雯清: "计及风光储发电和电动汽车的配电网规划" * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112488535A (en) * | 2020-12-01 | 2021-03-12 | 中国科学院地理科学与资源研究所 | Energy network security measure method, system and data platform |
| CN112994020A (en) * | 2021-03-31 | 2021-06-18 | 南京信息工程大学 | Multi-energy system load flow decomposition calculation method |
| CN113255105A (en) * | 2021-04-26 | 2021-08-13 | 上海电力大学 | Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure |
| CN115562029A (en) * | 2022-10-17 | 2023-01-03 | 杭州天然气有限公司 | Intelligent control method and system for natural gas turbine expansion generator set |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111400937B (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111400937A (en) | Load flow calculation method of comprehensive energy system | |
| Jiang et al. | Modelling and operation optimization of an integrated energy based direct district water-heating system | |
| CN111428351B (en) | Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method | |
| CN107808216B (en) | Method for constructing comprehensive minimum optimization model for wind curtailment and electrical thermal load reduction of electricity-gas-heat interconnection system | |
| CN103439941B (en) | A kind of gas internal-combustion engine cold, heat and power triple supply system optimizing operation method | |
| CN108241782B (en) | Mixed transient analysis method and mixed transient analysis system | |
| CN108628176B (en) | Thermoelectric combined optimization scheduling method considering heat supply network flexibility | |
| CN106844895B (en) | Decoupling calculation method for combined cooling heating and power micro-grid energy flow | |
| CN107947182B (en) | Dynamic power flow analysis method and dynamic power flow analysis system | |
| CN110232217B (en) | Comprehensive energy power distribution system operation domain modeling method | |
| CN103248048B (en) | Abandoned wind electric quantity estimation method | |
| CN113159983B (en) | Coordinated control method for off-grid phosgene combined heat and power generation comprehensive energy system | |
| CN110955954A (en) | Optimal load reduction method for layered decoupling electric and thermal comprehensive energy system | |
| WO2024109327A1 (en) | Integrated energy operation control method and integrated energy system based on multi-energy complementation | |
| Li et al. | Classification and location scheme selection of coupling components in integrated electrical and heating systems with renewable energy | |
| CN113806972B (en) | Comprehensive energy system reliability analysis method considering supply and demand bilateral flexibility | |
| Li et al. | Dynamic modeling and operations of a heat-power station system based on renewable energy | |
| Man et al. | State estimation for integrated energy system containing electricity, heat and gas | |
| CN114444850A (en) | Park energy storage planning method, device, equipment and medium based on digital twin | |
| CN113156835A (en) | Modeling simulation method for operation control of electricity-heat comprehensive energy system | |
| CN114549232A (en) | Hybrid energy flow calculation method for electricity-heat comprehensive energy system | |
| CN109345030B (en) | Multi-microgrid comprehensive energy system thermoelectric energy flow distribution type optimization method and device | |
| CN113191638A (en) | Electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network | |
| CN114971000B (en) | Comprehensive energy system configuration optimization method considering unbalance rate of supply and demand | |
| Yin et al. | A data analysis system for thermal energy loss in electrical automation management |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |