CN108416507B - Static sensitivity analysis method for electricity-gas coupling comprehensive energy system - Google Patents
Static sensitivity analysis method for electricity-gas coupling comprehensive energy system Download PDFInfo
- Publication number
- CN108416507B CN108416507B CN201810123783.0A CN201810123783A CN108416507B CN 108416507 B CN108416507 B CN 108416507B CN 201810123783 A CN201810123783 A CN 201810123783A CN 108416507 B CN108416507 B CN 108416507B
- Authority
- CN
- China
- Prior art keywords
- gas
- node
- power
- sensitivity
- iegs
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000008878 coupling Effects 0.000 title claims abstract description 34
- 238000010168 coupling process Methods 0.000 title claims abstract description 34
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 34
- 230000003068 static effect Effects 0.000 title claims abstract description 22
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 19
- 230000035945 sensitivity Effects 0.000 claims abstract description 62
- 238000002347 injection Methods 0.000 claims abstract description 49
- 239000007924 injection Substances 0.000 claims abstract description 49
- 239000007789 gas Substances 0.000 claims description 215
- 239000011159 matrix material Substances 0.000 claims description 40
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 30
- 239000003345 natural gas Substances 0.000 claims description 16
- 150000001875 compounds Chemical class 0.000 claims description 4
- 239000000243 solution Substances 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 230000003993 interaction Effects 0.000 abstract description 14
- 230000007246 mechanism Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 17
- 238000004364 calculation method Methods 0.000 description 14
- 230000002411 adverse Effects 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 239000002737 fuel gas Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000013439 planning Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 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 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013097 stability assessment Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012384 transportation and delivery Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (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)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a static sensitivity analysis method for an electricity-gas coupling comprehensive energy system, which is used for analyzing an interaction mechanism between an electric power energy supply system and a gas energy supply system. Firstly, a unified power flow model of an electric-gas coupling comprehensive energy system is provided; on the basis, defining the gas pressure-node injection power sensitivity of the electric-gas coupling comprehensive energy system; and finally, analyzing the influence of the power grid node injection power on the gas pressure by combining the sensitivity index of the comprehensive energy system in a typical scene, and positioning weak links of the comprehensive energy system. The examples show that the method can provide auxiliary information for the safe and stable operation of the regional comprehensive energy system, and effectively improve the safety of the system.
Description
Technical Field
The invention relates to the field of an electric-gas coupling comprehensive energy system, in particular to a static sensitivity analysis method of the electric-gas coupling comprehensive energy system based on a unified power flow model.
Background
With the increasing scarcity of fossil energy and the continuous deterioration of the environment, energy transformation becomes a necessary way for realizing sustainable development of the economy and the society. The existing mode that each energy supply system is independently planned and independently operated is broken through, a brand new energy system with the energy system step towards multi-energy integration and integration complementation is realized, and the method is an important way for promoting the process. Therefore, the ' thirteen-five ' planning of energy development ' clearly puts forward ' implementation of multi-energy complementary integrated optimization engineering ' and ' overall planning of infrastructures such as electric power, gas, heat, cold supply and water supply pipe gallery ' and the like, and a first foundation is constructedThe integrated energy supply system promotes the development of comprehensive utilization of various energy sources from the national level[1]。
Integrated Energy Systems (IES) are known as one of the main ways to realize comprehensive utilization of multiple Energy sources, and can improve the efficient utilization efficiency of Energy sources and realize the full consumption of renewable Energy sources by scientific management and optimized scheduling of cold/heat/electricity/gas/hydrogen and other Energy sources of different grades in the system[2]. However, the IES also carries the integrity risk of safe operation of the system while achieving the above objectives. In the face of the tight coupling of multiple energy sources, different energy supply systems "pull one another to move the whole body", and the interactive shadow response among them is closely concerned. On one hand, the fault of the gas system is transmitted to the power system, and the power failure accident is directly caused, for example, the power failure accident of 8.15 times in the Taiwan area of China in 2017 is caused by the disconnection of 6 units caused by the interruption of the natural gas supply, the power failure of 668 ten thousand users is caused by the accident, and the influenced population is more than 85 percent of the total population of the Taiwan area of China[3]. Similar events also occur in the united states, where gas plant natural gas supply is inadequate due to aiso Canyon natural gas leaks 2015, which severely affect the proper operation of local power systems[4]. On the other hand, the adverse effects of the power system are also conducted to the natural gas system and jeopardize the safe operation of the natural gas system. Along with the improvement of the permeability of the American renewable energy, the frequent adjustment of a gas power plant as a main peak regulation resource causes the great fluctuation of the pressure of a gas pipe network, and the gas transmission safety of a natural gas system is directly influenced[5]. Meanwhile, the power link of the IES often includes a Distributed Generation (DG), and the intermittent fluctuation of the output will be transmitted to the whole IES through the coupling link (such as a gas turbine set) between the two, which will adversely affect the whole system[6]. In view of the above, the adverse factors of each energy supply system must cause great concern to the influence of the safe operation of the electric-Gas coupled Integrated Electric and Gas Systems (IEGS). Without loss of generality, the influence of adverse factors on the safe operation of the system is often determined by the safe operation margin of the weak link of the system, so that the system is fast and convenient to operateThe positioning of the weak link of the IEGS becomes a problem which needs to be solved urgently at present.
Document [7] establishes a steady-state analysis model of the electric-gas coupling system, and considers the interaction of the two systems from the energy flow perspective; further, a literature [8] provides an optimal energy flow model of the electric-gas coupling system on the basis of considering network operation constraints, and provides decision assistance for optimal scheduling of the system; document [9] analyzes the influence of the increase of the load of the heat supply network on the voltage of the nodes of the power grid and the heat supply temperature of the nodes of the heat supply network by adopting a point-by-point method based on load flow calculation; the research of the document [10] shows that the continuous change of the natural gas load can bring important influence on the pressure of a gas network in the IEGS system; document [11] analyzes the interaction characteristics of each electric-gas energy supply network in a steady state by using a hybrid power flow algorithm based on an energy hub model.
However, the above research focuses on the trend calculation of the IEGS, and does not fully consider the influence of different energy supply network states (such as network topology, pipeline structure, load level, etc.) on the interaction between different energy sources of the integrated energy system, and how to determine the weak link of the IEGS.
Disclosure of Invention
The invention provides a static sensitivity analysis method for an electricity-gas coupling comprehensive energy system, which defines a gas pressure-node injection power sensitivity index of the electricity-gas coupling comprehensive energy system; analyzing the influence of the power grid node injection power on the gas pressure by combining the sensitivity indexes of the comprehensive energy system in a typical scene; weak links of the comprehensive energy system are positioned to provide decision basis for operation control (especially site selection of a gas storage device) of the comprehensive energy system, and the detailed description is as follows:
a static sensitivity analysis method for an electric-gas coupling integrated energy system comprises the following steps:
establishing an IEGS unified power flow model consisting of a gas system model, a power distribution system model and an energy coupling link model;
acquiring a current stable operation point and a Jacobian matrix of the current operation point according to an IEGS unified power flow model;
calculating the gas pressure-node injection power sensitivity according to the Jacobian matrix; and (3) sequencing the gas pressure-node injection power sensitivity, wherein the node set with high sensitivity is the weak link of the IEGS.
The IEGS unified power flow model specifically comprises the following steps:
wherein x is [ theta, V, p ═ p]Representing the state variables of the IEGS, namely phase angles of nodes of the power distribution network except for a balance node, voltage of an electric load node and air pressure of a gas load node; u ═ PPDS，QPDS，LNGS]Representing disturbance variables of the IEGS, and respectively injecting active power, reactive power and gas load into the nodes; a represents an independent parameter; pPDSInjecting active power into the node;
is the grid node voltage; y is a node admittance matrix,*is complex conjugate; a. theNGSIs a node-branch incidence matrix of the gas network;
the output electric power of the light gas turbine; c. CgeIs the conversion factor.
The gas pressure-node injection power sensitivity calculated according to the Jacobian matrix is specifically as follows:
respectively calculating gas pressure-gas load sensitivity and output-node injection power sensitivity of the light gas turbine according to the Jacobi matrix;
and calculating the gas pressure-node injection power sensitivity through the gas pressure-gas load sensitivity and the output-node injection power sensitivity of the light gas turbine.
The gas pressure-node injection power sensitivity is as follows:
wherein, PPDSInjecting active power into the node; p is gas node pressure; l isNGSIs a gas load;
represents the output electrical power of the MT;
the gas consumption for MT.
Further, the gas pressure-gas load sensitivity SggEqual to Jacobian matrix J in natural gas system trend solutionggNegative number of inversion:
in the formula, ANGSIs a node-branch incidence matrix of the gas network; f. ofiThe flow rate of the pipeline i is measured; n shapeiIs the pressure difference of the pipeline i; n ispThe number of pipelines; t is matrix transposition; diag () represents the construction of a diagonal matrix.
In the specific implementation, if a certain node i injects power
Changes occur, for the distribution network:
in the formula (I), the compound is shown in the specification,
for the network loss of the distribution network, to
The derivation comprises:
in the formula, the second term on the right side of the equation is the net loss micro-gain rate, which represents the variation relationship between the net loss and the node injection power.
The technical scheme provided by the invention has the beneficial effects that:
1. the method defines the air pressure safety index of the comprehensive energy system and positions the weak link of the comprehensive energy system;
2. the method analyzes the action mechanism of unstable operation air pressure of the power system in the comprehensive energy system in a typical scene due to the fluctuation of load power or the fluctuation of the new energy generated output;
3. the method provides corresponding guidance for selecting the gas storage position of the comprehensive energy system.
Drawings
FIG. 1 is a flow chart of a static sensitivity analysis method for an electric-gas coupling integrated energy system;
FIG. 2 is a schematic diagram of an electric-to-gas coupled integrated energy system;
FIG. 3 is a schematic view of an IEGS algorithm;
FIG. 4 shows S of each node of the distribution networkeeA schematic diagram;
FIG. 5 shows S of each node of the gas networkggA schematic diagram;
FIG. 6 is a gas pressure-nodal injection power sensitivity diagram;
FIG. 7 is a schematic view of air pressure at different load growth levels;
fig. 8 is a schematic view of the gas pressure in different embodiments.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention are described in further detail below.
The sensitivity analysis technique has been widely used in electricityThe operation of the force system is optimized, wherein a static analysis method based on a power flow equation can analyze the change situation between different physical quantities through the differential relation between the two, so that the information of weak nodes, key branches and the like of the system is judged[12]。
Therefore, the embodiment of the invention is expected to provide a static sensitivity analysis method for an electricity-gas coupling integrated energy system aiming at IEGS (electronic-gas systems) by means of a mature sensitivity analysis method in a power system. The method considers the network state of the power-gas system, analyzes the interaction rule between the gas pipeline pressure and the electric power injection space, deeply excavates the interaction influence between different energy supply networks in the IEGS, and judges the weak area influencing the safe operation of the system.
Example 1
An analysis method of static sensitivity for an electric-gas coupling integrated energy system, referring to fig. 1, comprises the following steps:
101: establishing an IEGS unified power flow model consisting of a gas system model, a power distribution system model and an energy coupling link model;
102: acquiring a current stable operation point and a Jacobian matrix of the current operation point according to an IEGS unified power flow model;
103: respectively calculating gas pressure-gas load sensitivity and output-node injection power sensitivity of the light gas turbine according to the Jacobi matrix, so as to calculate the gas pressure-node injection power sensitivity;
104: and (3) sequencing the gas pressure-node injection power sensitivity, wherein the node set with higher sensitivity is the weak link of the IEGS.
In summary, in the embodiment of the present invention, the network state of the power-gas system is considered through the above steps 101 to 104, and the interaction rule between the gas pipeline pressure and the electric power injection space is analyzed, so that the method quickly and directly obtains the weak area of the IEGS, and provides a decision basis for the operation control of the IEGS.
Example 2
The scheme in example 1 is further described below with reference to a specific calculation formula, fig. 2 and a specific example, and is described in detail in the following description:
201: IEGS modeling;
the embodiment of the invention takes IEGS comprising a Natural Gas System (NGS), a Power Distribution System (PDS) and a coupling link as an example to verify the effectiveness of the method.
Wherein, the NGS comprises an air source, a gas pipeline, a gas load, a compressor, a valve and the like. The valves are used to control the flow or cut-off of the gas in the pipeline, assuming that the valves are only in either a fully open or fully closed state, and thus the network topology of the NGS is deterministic. The IEGS is connected to the large grid through distribution transformers and assumes that the IEGS contracts for power delivery with the upstream grid. The energy coupling link is a key link of interaction of all energy supply systems in the IEGS, and has the characteristic of various coupling equipment. For example, a light-duty gas turbine (Micro-turbines, MT) realizes the coupling of a gas system and a power grid; through the electric Gas conversion technology (Power to Gas, P2G), electric energy can be converted into natural Gas energy. Considering that the MT is widely used in the IEGS, the embodiment of the present invention sets the MT as a coupling link of the IEGS. An electric-gas coupling comprehensive energy system researched by the embodiment of the invention is shown in figure 2.
In the embodiment of the invention, the IEGS and the external large power grid are assumed to sign a power supply contract, and the contract specifies the interactive electric quantity of the power distribution network and the large power grid in different periods, so that unbalanced power brought by various adverse factors in the operation process is borne by the MT, namely the MT is used as a balancing unit of the power distribution network.
Consider a typical conducted scenario of IEGS penalties: after summer, the power shortage in the system is complemented by MT, which causes the change of gas load and further causes the change of gas pressure; or MT output fluctuation caused by DG short-time output fluctuation represented by wind power generation or photovoltaic, and further causes the change of gas pressure of a gas-gas network. In summary, changes in the injected power at any node of the distribution network affect the gas pressure level of the gas network through the coupling link. In some extreme scenarios, even the operation problems of the gas system, such as large fluctuation of gas pressure and gas pressure out-of-limit, may be caused.
Therefore, the gas pressure level of the gas network can be taken as one of important indexes for safe operation of the IEGS, and an area with fast gas pressure reduction in the IEGS is considered as a weak link of the IEGS when adverse factors act on the system. The weak link of the system is obtained, and the gas storage equipment is arranged at the weak link, so that the method is a preventive control measure for improving the IEGS gas pressure safety margin.
1) Gas system model
The model of the gas system is mainly a steady-state flow equation of the pipeline, wherein the model can be divided into a pipeline with a compressor and a pipeline without the compressor. The pressure drop exists at two ends of the pipeline, and particularly in a long-distance large-capacity medium-high pressure gas transmission network, a compressor is required to be configured to raise the pressure of the gas transmission pipeline. Considering that the IEGS studied by the embodiment of the invention belongs to the category of low-pressure distribution networks (0-75mbar), the embodiment of the invention does not consider the influence of the compressor for the moment[13]。
The node variables of the gas system comprise injected gas flow and node air pressure, the node classification of the power system is simulated, and the nodes can be divided into nodes with known pressure and nodes with known injected flow according to known variables. In a gas system, the gas supply is a balanced node with known gas pressure and unknown injection flow, and the gas load has unknown gas pressure and known gas demand. If the gas system has ngA node, npAnd the joint 1 of the strip pipeline is an air source joint, and the rest are gas load joints. For connecting nodes i, j (i, j ═ 1,2, ·, ng) K, a steady flow f in a low-pressure gas line (0-75mbar)k(k＝1,2,···,np) Can be described by the following formula[10,13]：
Π＝-(ANGS)Tp (3)
In the formula, λkIs the coefficient of friction of the pipe k; II typekRepresenting the pressure difference across the conduit k; dkAnd LkRepresents the diameter and length of the pipe k; g is the relative density of the fuel gas; p represents gas node pressure; a. theNGSIs a node-branch incidence matrix of the gas network. When pik>At 0, s k1 is ═ 1; when pik<At 0, sk＝-1。
By analogy with kirchhoff's first and second laws, the gas flow in the network should satisfy the following two conditions: the gas flow injected into a certain node is equal to the gas flow flowing out of the node; for any loop of the network, the sum of the pressure drops of the gas in the flowing process is zero. Thus, the flow behavior of the gas in the network can be described by the following equation:
ANGSf＝LNGS (4)
in the formula, f is a steady-state flow column vector of the pipeline; l isNGSIs a gas load.
2) Power distribution system model
The embodiment of the invention effectively solves the problem of three-phase imbalance of the power distribution network by assuming phase commutation[14]Therefore, the embodiment of the invention ignores the influence of the three-phase imbalance on the power flow calculation. If the distribution network has neA node, node 1 being a balanced node, node 2 to node 1+ npvThe PV node (the node given active power P and voltage V) and the remaining PQ node (the node given active power P and reactive power Q). The model of the power distribution network is a node power equation reflecting the relationship between node power and node voltage and phase angle[14]：
In the formula, PiInjecting active power for the node i; qiInjecting reactive power for node i;
active power and reactive power are generated by a generator on a node i;
the active power and the reactive power of the load on the node i; vi，VjIs the voltage at node i, j; gij，BijIs the admittance Y between nodes i, jijThe real part, the imaginary part; thetaijIs the phase angle difference between node i and node j (i, j ═ 1,2, ·, n)e)。
3) Energy coupling link model
The embodiment of the invention assumes that the output adjustment of MT is in the second level[15]When considering the interaction effect among different energy supply networks, the dynamic characteristics of the MT are ignored, and the interaction effect of the networks is considered from the energy interaction point of view. The input of the MT is fuel gas, and high-temperature gas generated by fuel gas combustion is expanded to do work so as to output electric energy. The model of the MT is therefore a static model characterizing its steady state input-output characteristics, as follows:
in the formula (I), the compound is shown in the specification,
represents the output electrical power of the MT;
represents the gas consumption of the MT; c. CgeThe conversion factor for MT is represented.
202: IEGS static sensitivity analysis method; the IEGS power flow solving model can be divided into a unified power flow model and a discrete solving model[10]The gas pressure-node injection power sensitivity index is further defined on the basis of a unified power flow model, and weak links of a positioning system are analyzed through static sensitivity.
On the basis of the IEGS unified power flow model, the gas pressure-node injection power sensitivity index is further defined, and weak links of a positioning system are analyzed through static sensitivity.
1) An IEGS unified power flow model;
for the IEGS, the essence of power flow calculation is the process of finding the stable operating point of the system given a series of conditions, and the unified power flow model can be described as:
wherein x is [ theta, V, p ═ p]Representing the state variables of the IEGS, namely phase angles of nodes of the power distribution network except for a balance node, voltage of an electric load node and air pressure of a gas load node; u ═ PPDS，QPDS，LNGS]Representing disturbance variables of the IEGS, and respectively injecting active power, reactive power, gas load and gas load into the nodes; and a represents independent parameters, such as power system node admittance, gas system network topology, pipeline parameters and the like. PPDSInjecting active power into the node;
is the grid node voltage; y is a node admittance matrix, which is illustrated in equations (5) - (8),*is a mathematical operation sign of complex conjugate.
The unified power flow model is solved through a Newton Raphson method, and the iteration form is as follows:
in the formula,. DELTA.x(k+1)The variation of the state variable of the IEGS in the iteration k + 1; x is the number of(k+1)State variables of the IEGS in this iteration are k + 1; delta F is a deviation value of the power flow equation set; j is the Jacobian matrix, which can be expressed as[7]：
And (3) solving a power flow equation set in the joint type (10) - (12) according to the initial condition given by the system to obtain the current stable operation point of the system.
The meanings of the elements in the jacobian matrix J are well known to those skilled in the art, and are not described in detail in the embodiments of the present invention.
2) Gas pressure-node injection power sensitivity matrix
The influence of adverse factors on different energy supply networks is equivalent to analyzing the change of the system operation state of the IEGS under the action of adverse factors (such as node injection power change), and the adverse factors can be generally characterized by disturbance variables of the system, so that the IEGS static sensitivity analysis is generally expressed as:
according to the variation relation of different physical quantities in different energy supply networks, the sensitivity indexes with various forms can be constructed. The interaction between the two energy supply networks is conducted via the MT, which is set as a balancing machine for the distribution network according to the scenario in section 1 and connected to node h in the gas network. Considering the typical conduction scene of adverse factors, the embodiment of the invention takes the fluctuation of the power grid node injection power as the adverse factor of the IEGS, and takes the gas pressure as the state variable of the IEGSDefining a gas pressure-node injection power sensitivity matrix S according to equation (13)geThe following were used:
wherein the content of the first and second substances,
is a matrix SgeRepresents the variation of the injection power of the distribution network node i (i ═ 1,2, ·, n)e) For gas node j (j ═ 1,2, ·, ng) The influence of the air pressure.
Definition of SggFor gas pressure-gas load sensitivity, SeeThe power injection sensitivity for MT force-node is given by equation (9)
Is cegThen S isgeIs the product of these three.
For the gas net, S can be found by analyzing the formulas (1) to (4)ggEqual to Jacobian matrix J in natural gas system trend solutionggNegative number of inversion:
in the formula, the function diag () represents constructing a diagonal matrix.
For the distribution network, if a certain node i injects power
Changes occur, for the distribution network:
in the formula (I), the compound is shown in the specification,
for distribution network loss, equation (19) injects power to the node considering other nodes injecting power is not changed
The derivation comprises:
in the formula, the second term on the right side of the equation is the net loss micro-gain rate, which represents the variation relationship between the net loss and the node injection power. Intuitively, the MT contribution for a lossless network is the same as the node injected power variation. However, for a real grid, the network topology, the load level will affect the net loss increase rate of the network and thus the energy interaction between the energy supply networks. The network loss is determined by the voltage and phase angle of each node[16]Further deducing that the net loss micro-increment rate is as follows:
note that the Jacobian matrix J of the power grid in equation (12)eeThe inversion can be expressed as:
comparing equations (20) and (22), the Jacobian matrix J can be seeneeCan be solved for SeeProviding corresponding information, the joint type (20) - (22) can solve See。
203: and (3) gas pressure-node injection power sensitivity analysis.
It should be noted that the difference of the MT access positions will affect the power flow distribution of the IEGS, so the deployment position of the MT needs to be given before calculating the sensitivity matrix, and the initial conditions of the power flow calculation are changed to form different stable operation points. The sensitivity analysis method comprises the following steps:
1) reading natural gas network status information, comprising: natural gas pipeline parameters, pipeline topology information, form ANGSReading the network information of the power distribution network to form Y;
2) solving a power flow equation set according to the formulas (10) to (12) to obtain a current stable operation point of the IEGS and a Jacobian matrix J of the current operation point;
3) from the Jacobian matrix J, S is calculated according to expressions (16), (17), (20) to (22)gg、See；
4) Gas pressure-node injection power sensitivity S is calculated according to equation (14)ge；
5) Sensitivity to gas pressure-node injected power SgeAnd (3) sorting the sensitivity values, wherein a set of nodes with larger sensitivity values (the specific value is set according to the requirement in practical application, the first 30% of the sensitivity values are generally selected, and the method is not limited in the embodiment of the invention) is a weak link of the IEGS.
In summary, in the embodiment of the present invention, the network state of the power-gas system is considered in the steps 201 to 203, and the interaction law between the gas pipeline pressure and the electric power injection space is analyzed, so that the method quickly and directly obtains the weak area of the IEGS, and provides a decision basis for the operation control of the IEGS (such as the site selection of the gas storage device).
Example 3
The feasibility of the schemes of examples 1 and 2 is demonstrated below in specific calculations, in conjunction with fig. 3-8, and described in detail below:
to verify the effectiveness of the method, a typical IEGS is taken as an example for simulation explanation. As shown in FIG. 3, the IEGS embodiment of the present invention is formed by coupling an IEEE-33 node power distribution system and a modified 11-node gas network through an MT[13,17]EBi and GBi denote grid nodes and gas nodes, respectively.
The power distribution network is connected with an external large power grid through EB1, the electric power obtained by the power distribution network and the external large power grid is assumed to be 3500kW, and EB1 is a PV node at the moment. The MT is connected to GB11 of the gas network and EB2 of the distribution grid as a balancing unit, so EB2 is the balancing node of the distribution grid. Assuming that the safety operation warning air pressure of load nodes (GB 2-GB 11) in the gas network is 20mbar, the air pressure of an air source node GB1 is 75mbar, cgeIs 6.1kW/m3. The gas network data are shown in tables A1 and A2.
TABLE A1 gas pipeline data
TABLE A2 gas load data
1) Sensitivity calculation and weak link analysis firstly, the operating point of the system and the Jacobian matrix J are calculated according to the initial conditions of the IEGS. Calculating S from equations (20) to (22)eeThe results are shown in fig. 4, where the size of the bar graph represents the sensitivity of MT output variation to grid node injected power variation. It can be seen from fig. 4 that the farther the electrical distance from the MT, the more pronounced the fluctuation of the MT output caused by the change in node injection power. As shown in the EB18 of fig. 3, the MT is far from the EB2 where the MT is located, so that when the load in the EB18 node increases by 1 unit of load demand, the MT should increase by 1.1404 units of active power. And the EB19 is closer to the MT in fig. 3, and the MT only needs to increase the active power by 1.0007 unit after the load requirement of 1 unit is increased. The power distribution network has different positions of the node injection power change and different contribution influences on MT, and the reason is that the network loss is different due to the fact that the distance between the electrical distance and the distance in the process of carrying out active power rebalancing on the power distribution network.
Calculating S from equations (16) and (17)ggAs shown in fig. 5, the size of the bar graph in the graph indicates the sensitivity of the gas network pressure change to the gas load change. Note that the calculation result in fig. 5 is a result when the MT accesses GB11 in the gas system.
When the MT gas load demand increases by 1 unit, the GB11 gas pressure decreases by 0.3175 units correspondingly, and the GB2 gas pressure close to the gas source decreases by 0.1155 units correspondingly. From equation (1), it can be seen that as the flow in the gas pipeline increases, the pressure difference across the pipeline will increase. When newly added 1 unit of fuel gas is conveyed to GB11 through the fuel gas network, the increase of the air pressure difference is caused by the increase of the flow in the pipeline in the network, and the increase of the air pressure difference is accumulated together, which is reflected that the air pressure drop of GB11 is most obvious and the air pressure drop of GB2 at the upstream of the network is not obvious.
Binding to See、SggFurther calculating S according to equation (14)geAs shown in fig. 6, for a certain point (x, y, z) in the graph, the sensitivity z of the air pressure of the air grid node y to the power fluctuation of the grid node x is shown. According to fig. 6, it can be known that the sensitivity values corresponding to the gas grid nodes 8, 9, and 11 are large, and the air pressure change speed caused by the power fluctuation of the grid nodes is higher than that of other nodes. And taking the node 30% before the sensitivity value as a weak node of the IEGS, wherein the weak node of the gas network is { GB8, GB9 and GB11 }.
To verify the effectiveness of the method, the changes in barometric pressure at different injection power increase levels are simulated by load flow calculation as shown in fig. 7. The active growth rate (compared with the original node load level) of each load node is 2%, and when the load growth rate of the node is 10%, the GB11 air pressure is out of limit (circle in the figure). As can be seen from the data in fig. 6, the MT output increases due to the increase of the load, the corresponding air pressure of the air network decreases, and the air pressure of the weak node reaches the boundary of safe operation first. On the other hand, from an actual physical system, GB11 is at the end of the line and is heavily loaded in fig. 3, and the gas load variation thereof has the greatest influence on the system gas pressure. The sensitivity analysis result can be seen to conform to the actual situation of the network.
2) Optimization of gas storage locations
The arrangement of the gas storage device at the weak node of the system is an effective means for improving the gas pressure level and ensuring the safe operation of the system. In the peak period of gas consumption in the short time of the system, the gas storage device plays a role in stabilizing the gas pressure level and ensuring continuous and stable gas flow by injecting gas into the system. However, it is not economical to provide gas storage devices for all nodes, and embodiments of the present invention provide gas storage at weak nodes of the system according to the sensitivity calculation result.
For example, a capacity of 20m is set for a load increase of 10%3The gas storage device. The selection of the gas storage position comprises the following 3 schemes:
1) without gas storage
2) Gas storage is arranged at a non-weak link such as GB 4;
3) and gas storage is arranged at the system weak node GB 11.
In the load flow calculation, the gas storage device can be processed by approximately reducing the corresponding gas load. The calculation of the air pressure level of the system at different access positions by tidal current is shown in fig. 8. Although the whole air pressure level is effectively raised by accessing the air storage equipment at the GB4 node, compared with the scheme 2 and the scheme 3, the air storage is arranged at the weak node, so that the air pressure lifting of the system is more effective, and the safe operation margin of the system is obviously improved. The static sensitivity analysis method based on the unified power flow model in the embodiment of the invention is based on the unified power flow model of the power system and the gas system, researches an IEGS gas pressure-node injection power sensitivity matrix, analyzes the influence of the power grid node injection power on the gas pressure, positions weak links of the regional comprehensive energy system, provides auxiliary information for the safe and stable operation of the comprehensive energy system, and effectively improves the safety of the system. The following conclusions can be drawn by example analysis:
1) the gas pressure-node injection power sensitivity of the node can effectively reflect the safe operation margin of the gas pressure of each load node of the gas network;
2) the sensitivity index can be used for quickly positioning weak links of the system, accurately determining the most easily out-of-range node of the air pressure in the system and avoiding mass calculation of a point-by-point method; 3) the gas pressure-node injection power sensitivity is real, the influence degree of fluctuation of the load node injection space on the gas network gas pressure level is effectively reflected, and the gas storage equipment is arranged at the weak link of the system, so that the gas pressure safe operation margin can be effectively improved, and the safety of the system is improved.
Reference to the literature
[1] Energy development of the national energy agency, thirteen five programs [ C ]2016:17-18.
[2] Gi Hongjie, Wang Dan, Xuxiong, et al, regional Integrated energy systems several problems research [ J ] electric power systems automation, 2015,39(7): 198-.
JIA H,Wang D,XU X,et al.Research on some key problems related to integrated energy systems[J].Automation of Electric Power Systems,2015,39(7):198-207.(in Chinese).
[3] Wangping, Taiwan with power failure, who should carry on the back of the pot [ N ], Japanese newspaper of people overseas edition, 2017-08-16(004).
[4]Fairley P.The troubled link between gas and electricity grids[News][J].IEEE Spectrum,2016,53(6):11-12.
[5]Chertkov M,Backhaus S,Lebedev V.Cascading of fluctuations ininterdependent energy infrastructures:Gas-grid coupling[J].Applied Energy,2015,160:541-551.
[6]Touretzky C R,McGuffin D L,ZIEGSmer J C,et al.The effect of distributed electricity generation using natural gas on the electric and natural gas grids[J].Applied Energy,2016,177:500-514.
[7]Martinez-Mares A,Fuerte-Esquivel C R.A unified gas and power flow analysis in natural gas and electricity coupled networks[J].IEEE Transactions on Power Systems,2012,27(4):2156-2166.
[8]Zlotnik A,Roald L,Backhaus S,et al.Coordinated scheduling for interdependent electric power and naturalgas infrastructures[J].IEEE Transactions onPower Systems,2017,32(1):600-610.
[9] Wang Ying Rui, Zengbo, Guo Jing, et al.
WANG Y,ZENG B,GUO J,et al.Multi-energy flow cal-culation method for integrated energy system containing electricity,heat and gas[J].Power System Technology,2016,40(10):2942-2951.
[10] Wangweiliang, Wangdan, Jiahongjie, etc. consider the steady state analysis [ J ] of the power-natural gas area comprehensive energy system of the natural gas network state, China Motor engineering report, 2017,37(5):1293 plus 1304.
Wang Weiliang,Wang Dan,Jia Hongjie,et al.Steady state analysis of electricity-gas regional integrated energy system with consideration of NGS network status[J]Proceedings of the CSEE,,2017,37(5):1293-1304.
[11] Xuxiong, Jia hongjie, long and so on, regional comprehensive energy system electricity/gas/heat hybrid trend algorithm research [ J ], Chinese electro-mechanical engineering report 2015,35(14): 3634-.
Xiandandong X, Hongjie J, Xiaolong J, et al, study on hybrid heat-gas-power flow algorithm for integrated communication energy system [ J ]. Proceedings of the CSEE,2015,35(14) 3634 and 3642 [12] Luo Yuan, Luo Ting, Suliuqing, Liubaiyang.
LUO Qingyue,LUO Ting,SUN Liuqing,et al.A multi-index integrated static voltage stability assessment method for weak buses identification[J].Power System Protection and Control,2017,V45(3):82-86.
[13]Abeysekera M,Wu J,Jenkins N,et al.Steady state analysis of gas networks with distributed injection of alternative gas[J].Applied Energy,2016,164:991-1002.
[14] The research on the real-time online treatment method of three-phase load imbalance in the power distribution station area [ J ]. the journal of the electrical engineering of China, 2015,35(9): 2185-plus 2193.
[15]Xu X,Jia H,Chiang H D,et al.Dynamic modeling and interaction of hybrid natural gas and electricity supply system in microgrid[J].IEEE Transactions on Power Systems,2015,30(3):1212-1221.
[16] The grid loss micro-augmentation rate algorithm [ J ] based on direct current flow is a grid technology, 2007,31(15):39-43.
KONG X,FANG D,HOU Y.An incremental transmission loss algorithm based on DC power flow[J].Power System Technology,2007,15:012.
In the embodiment of the present invention, except for the specific description of the model of each device, the model of other devices is not limited, as long as the device can perform the above functions.
Those skilled in the art will appreciate that the drawings are only schematic illustrations of preferred embodiments, and the above-described embodiments of the present invention are merely provided for description and do not represent the merits of the embodiments.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (5)
1. A static sensitivity analysis method for an electric-gas coupling comprehensive energy system is characterized by comprising the following steps:
establishing an IEGS unified power flow model consisting of a gas system model, a power distribution system model and an energy coupling link model;
acquiring a current stable operation point and a Jacobian matrix of the current operation point according to an IEGS unified power flow model;
calculating the gas pressure-node injection power sensitivity according to the Jacobian matrix; sequencing the gas pressure-node injection power sensitivity, wherein the node set with high sensitivity is the weak link of the IEGS;
the IEGS unified power flow model specifically comprises the following steps:
wherein x is [ theta, V, p ═ p]Representing the state variables of the IEGS, namely phase angles of nodes of the power distribution network except for a balance node, voltage of an electric load node and air pressure of a gas load node; u ═ PPDS，QPDS，LNGS]Representing disturbance variables of the IEGS, and respectively injecting active power, reactive power and gas load into the nodes; a represents an independent parameter; pPDSInjecting active power into the node;
is the grid node voltage; y is a node admittance matrix,*is complex conjugate; a. theNGSIs a node-branch incidence matrix of the gas network;
is lightElectric power output of the model gas turbine; c. CgeIs the conversion factor.
2. The static sensitivity analysis method for the electric-gas coupling integrated energy system as claimed in claim 1, wherein the calculating of the gas pressure-node injection power sensitivity according to the jacobian matrix specifically comprises:
respectively calculating gas pressure-gas load sensitivity and output-node injection power sensitivity of the light gas turbine according to the Jacobi matrix;
and calculating the gas pressure-node injection power sensitivity through the gas pressure-gas load sensitivity and the output-node injection power sensitivity of the light gas turbine.
3. The static sensitivity analysis method for the electric-gas coupling integrated energy system as claimed in claim 2, wherein the gas pressure-node injection power sensitivity is specifically:
wherein, PPDSInjecting active power into the node; p is gas node pressure; l isNGSIs a gas load;
represents the output electrical power of the MT;
the gas consumption for MT.
4. The method for analyzing the static sensitivity of the electric-gas coupling integrated energy system according to claim 2,
the gas pressure-gas load sensitivity SggEqual to Jacobian matrix J in natural gas system trend solutionggNegative number of inversion:
in the formula, ANGSIs a node-branch incidence matrix of the gas network; f. ofiIs the flow of the pipeline i; II typeiThe gas pressure difference of the pipeline i; n ispThe number of gas pipelines; t is matrix transposition; diag () represents the construction of a diagonal matrix.
5. The method for analyzing the static sensitivity of the electric-gas coupling integrated energy system according to any one of the claims 1 to 4,
if a certain node i injects power
Changes occur, for the distribution network:
in the formula (I), the compound is shown in the specification,
for the network loss of the distribution network, to
The derivation comprises:
in the formula, the second term on the right side of the equation is the net loss micro-gain rate, which represents the variation relationship between the net loss and the node injection power.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810123783.0A CN108416507B (en) | 2018-02-07 | 2018-02-07 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810123783.0A CN108416507B (en) | 2018-02-07 | 2018-02-07 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108416507A CN108416507A (en) | 2018-08-17 |
| CN108416507B true CN108416507B (en) | 2021-08-20 |
Family
ID=63127971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810123783.0A Active CN108416507B (en) | 2018-02-07 | 2018-02-07 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108416507B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110543695B (en) * | 2019-08-14 | 2023-08-01 | 天津大学 | Electric-gas coupling comprehensive energy system feasible region calculation method |
| CN110601185B (en) * | 2019-09-17 | 2021-07-06 | 武汉大学 | Unified power flow model-based comprehensive energy system weak point identification method |
| CN110647040B (en) * | 2019-10-10 | 2022-10-04 | 东北电力大学 | Safety control method and device of comprehensive energy system |
| CN111767660A (en) * | 2020-09-02 | 2020-10-13 | 国网江西省电力有限公司电力科学研究院 | Optimal arrangement method and device for gas storage devices of electric coupling energy system |
| CN112510680A (en) * | 2020-10-26 | 2021-03-16 | 华能定边新能源发电有限公司 | Collaborative planning method for distributed power supply and electric-to-gas equipment in comprehensive energy system |
| CN112964961B (en) * | 2021-02-04 | 2022-12-06 | 南方电网科学研究院有限责任公司 | Electric-gas coupling comprehensive energy system fault positioning method and system |
| CN113609651A (en) * | 2021-07-13 | 2021-11-05 | 华南理工大学 | Unified power flow model-based electric-gas interconnection system action mechanism analysis method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105186525A (en) * | 2015-10-29 | 2015-12-23 | 山东大学 | Reactive voltage control partitioning method under wind power integration |
| CN105356447A (en) * | 2015-10-27 | 2016-02-24 | 河海大学 | Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system |
| CN106777708A (en) * | 2016-12-21 | 2017-05-31 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2997387A4 (en) * | 2013-05-14 | 2017-01-18 | Rensselaer Polytechnic Institute | Methods of computing steady-state voltage stability margins of power systems |
-
2018
- 2018-02-07 CN CN201810123783.0A patent/CN108416507B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105356447A (en) * | 2015-10-27 | 2016-02-24 | 河海大学 | Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system |
| CN105186525A (en) * | 2015-10-29 | 2015-12-23 | 山东大学 | Reactive voltage control partitioning method under wind power integration |
| CN106777708A (en) * | 2016-12-21 | 2017-05-31 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
Non-Patent Citations (2)
| Title |
|---|
| Decentralized Optimization of Multi-Area Electricity-Natural Gas Flows Based on Cone Reformulation;Yubin He等;《IEEE》;20171229;全文 * |
| 一种考虑天然气系统动态过程的气电联合系统优化运行模型;艾小猛等;《电网技术》;20171103;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108416507A (en) | 2018-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108416507B (en) | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system | |
| Sun et al. | Day-ahead optimization schedule for gas-electric integrated energy system based on second-order cone programming | |
| Zhang et al. | Optimal allocation of CHP-based distributed generation on urban energy distribution networks | |
| CN106777708B (en) | Steady state analysis method of electric power-natural gas regional comprehensive energy system | |
| Majid et al. | Minimization of transmission loss in application of HVDC networks under load increase scenario | |
| Chen et al. | Fast decoupled multi-energy flow calculation for integrated energy system | |
| CN110647040B (en) | Safety control method and device of comprehensive energy system | |
| CN112035984B (en) | Collaborative planning method for comprehensive energy system of electricity-gas-storage area | |
| Maihemuti et al. | Dynamic security and stability region under different renewable energy permeability in IENGS system | |
| He et al. | Economic dispatch analysis of regional Electricity–Gas system integrated with distributed gas injection | |
| CN110543695B (en) | Electric-gas coupling comprehensive energy system feasible region calculation method | |
| Li et al. | Economic dispatch of an integrated heat-power energy distribution system with a concentrating solar power energy hub | |
| Arabkhaburi et al. | Optimal placement of UPFC in power systems using genetic algorithm | |
| Tang et al. | A composite voltage stability index for integrated energy systems based on L-index and the minimum eigenvalue of reduced Jacobian matrix | |
| Mu et al. | Research on sensitivity analysis of wind power consumption capability of integrated energy system based on unified optimal power flow model | |
| Tian et al. | Modeling and simulation for multi energy flow coupled network computing | |
| CN109742757A (en) | A kind of generation method of integrated energy system load security domain | |
| CN113517697A (en) | Gas-electricity combined network steady-state load flow calculation method based on all-pure embedding method | |
| Zhai et al. | Reliability analysis of power-gas integrated energy system based on dynamic simulation | |
| Zhao et al. | Static sensitivity analysis method of integrated energy system based on power flow model | |
| CN116739404A (en) | Construction method of networking feature module in comprehensive energy system entropy state model | |
| Liu et al. | Feasible Operation Regions of Electricity-gas Integrated Energy Distribution System | |
| Zhang et al. | Optimal energy flow of electricity-gas integrated energy system using second-order cone program | |
| Ayan et al. | An examination of the effects of distributed generation on distribution systems by load flow analysis | |
| Hashemi et al. | Impact of distributed generation on unbalanced distribution networks |
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 |