CN110263387A - The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established - Google Patents
The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established Download PDFInfo
- Publication number
- CN110263387A CN110263387A CN201910456160.XA CN201910456160A CN110263387A CN 110263387 A CN110263387 A CN 110263387A CN 201910456160 A CN201910456160 A CN 201910456160A CN 110263387 A CN110263387 A CN 110263387A
- Authority
- CN
- China
- Prior art keywords
- constraint
- natural gas
- power
- node
- energy
- 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
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 162
- 239000003345 natural gas Substances 0.000 title claims abstract description 81
- 238000009826 distribution Methods 0.000 title claims abstract description 51
- 238000003672 processing method Methods 0.000 title claims abstract description 10
- 239000007789 gas Substances 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 9
- 230000005611 electricity Effects 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 3
- 238000012986 modification Methods 0.000 claims description 3
- 238000000746 purification Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000012384 transportation and delivery Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 11
- 238000010438 heat treatment Methods 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000003786 synthesis reaction Methods 0.000 description 4
- 238000005094 computer simulation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000002040 relaxant effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Quality & Reliability (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The integrated energy system optimal operation model and its linear processing methods established the present invention relates to a kind of power distribution network based on energy hub, natural gas grid, it can be by energy hub model that power distribution network, natural gas grid couple, the lowest cost is purchased as optimization aim using the energy, the safe operation constraint of power distribution network, natural gas grid and energy hub, including the constraint of power distribution network Branch Power Flow, node voltage constraint, branch power constraint are considered simultaneously;The constraint of natural gas line air balance, the constraint of gas source gas output, the constraint of node air pressure size constraint, natural gas line Weymouth steady-state load flow, pressurizer boosting relation constraint;The heating power balance of energy hub constrains and the operation power constraint of each energy conversion apparatus, to realize the optimization operation of integrated energy system.Before solving model, linearization process is carried out to initial model, finally to solve integrated energy system optimal operation model by CPLEX.
Description
Technical field
The present invention relates to field of power electronics, especially a kind of power distribution network based on energy hub, natural gas grid are established
Integrated energy system optimal operation model and its linear processing methods.
Background technique
Environment worsening and the energy are petered out the problem of having become global concern, how to make full use of clean energy resource and
Guarantee that the mankind's uses energy demand simultaneously again, it has also become the common focus of attention of countries in the world today.For the association for realizing various energy resources
It adjusts management and makes full use of the complementarity of different energy sources, integrated energy system comes into being.Integrated energy system refers to certain
Coordinated planning, the optimization operation between the various energy resources subsystems such as coal, petroleum, natural gas, electric energy, thermal energy are realized in region.
Therefore, while meeting diversification energy demand in system, it is necessary to effectively promote efficiency of energy utilization, promote the energy can
Sustainable development.
It by power grid, the combined operating of gas net, runs while firepower unit can be reduced, and then alleviates the further of fossil energy
The further deterioration of consumption and environmental pollution.Power grid can not only realize that the energy is complementary with the joint of gas net, and high efficiency may be used also
Sufficiently meeting using for user side can demand.Accordingly, it is considered to the integrated energy system optimization of new generation in combination of power distribution network, natural gas grid
Operation, the operating mechanism of comprehensive study energy resource system also help it is following set up a unified market value measurement standard,
So that the conversion and complementation of the energy fully demonstrate economy and society value, enable electric power promote to more cleaning, it is more green
The energy transformation, and will for future bring a more efficient, economic, renewable, reliable, sustainable and stable synthesis
Energy resource system brings long-term and positive height interests.
Summary of the invention
In view of this, the purpose of the present invention is to propose to a kind of power distribution network based on energy hub, natural gas grids to establish
Integrated energy system optimal operation model and its linear processing methods can consider the mutual shadow of power network, natural gas grid simultaneously
It rings, under the premise of meeting electricity, air and heat workload demand, integrated energy system is more economically and safely run.
The present invention is realized using following scheme: the synthesis that a kind of power distribution network based on energy hub, natural gas grid are established
Energy resource system optimal operation model, comprising the following steps:
Step S1: providing one can be by energy hub model that power distribution network, natural gas grid couple;
Step S2: objective function is established, and about according to the safe operation of power distribution network, natural gas grid and energy hub
Beam completes the modeling of integrated energy system optimal operation model.
Further, objective function described in step S2 is shown below for energy acquisition cost, including purchases strategiesWith purchase gas costTo make integrated energy system in dry run
Total cost during scheduling is minimum;
In formula, TR indicates the set of substation's injecting power node in power distribution network;GS indicates gas source node in natural gas grid
Set;Unit price of power, gas price in distribution expression period t;Function injects in substation in respectively period t
The purchase of electricity of rate node j;For the purchase tolerance of the gas source node j in period t.
Further, the safe operation of power distribution network described in step S2 constraint includes: node power Constraints of Equilibrium, node voltage
Filters with Magnitude Constraints, the constraint of distribution line maximum carrying capacity and substation's power constraint;
The node power Constraints of Equilibrium are as follows:
In formula, rij、xijThe respectively resistance of branch ij, reactance;δ (j) is indicated using j as the end of the branch of headend node
The set of node, π (j) are indicated using j as the set of the headend node of the branch of endpoint node;Pij,t、Qij,tRespectively branch ij's
Active power and reactive power;Respectively indicate major network power output and electric load active power;Respectively indicate use
The active power output of P2G device electric conversion in Hub;Indicate Hub in CHP device injection power distribution network etc.
Imitate active power;Respectively indicate the reactive power value of major network power output, electric load;
The node voltage amplitude constraint are as follows:
Umin≤Uj,t≤Umax
In formula, Uj,t、UminAnd UmaxRespectively node voltage amplitude and its lower and upper limit;
The distribution line maximum carrying capacity constraint are as follows:
0≤Iij,t≤Imax
In formula, Iij,t、ImaxRespectively branch current magnitudes and its upper limit;
Substation's power constraint are as follows:
Further, the particular content natural gas network constrained in step S2 are as follows:
The constraint of node air balance:
In formula, σ (jk) is indicated using j as the pipeline set of headend node;μ (ij) respectively indicates the pipe using j as endpoint node
Road set;fp,t/fl,tFor the flow of pipeline p/l;The respectively gentle source node gas output of gas load;Respectively indicate the amount of natural gas for GF, CHP energy conversion in (comprehensive energy unit) Hub;P2G device electricity is converted into the equivalent natural gas injection rate of gas in expression Hub;
The constraint of gas source gas output:
The natural gas being produced from gas well needs the purification by refinery, due to by gas at place capacity and gas well
The limitation of pressure, gas source gas output should meet following constraint in the unit time:
In formula,For the output amount of natural gas of gas source j,Respectively gas source j exports amount of natural gas
Bound;
Node air pressure constraint:
The air pressure of each node of natural gas network must be maintained in safe, reasonable range of operation, and mathematic(al) representation is such as
Under:
In formula,Respectively represent air pressure upper lower limit value at node i;
The constraint of natural gas line Weymouth steady-state load flow:
The relationship of gas flow and both ends air pressure, expression are portrayed using Weymouth steady-state load flow model
Are as follows:
Wherein,
In addition pipeline flow also should ensure that within safe and reasonable range of operation:
In formula, fp,tRepresent the natural gas flow that pipeline p flows through;πi,tAnd πj,tRespectively indicate pipeline p both ends air pressure size;φp
Indicate the air stream transportation parameter of pipeline p;sgnpIndicate the air flow direction of pipeline p;Indicate maximum pipeline transmission capacity.
Pressurizer boosting relation constraint:
πj,t=Γcπi,t
In formula, fc,tIndicate the delivery air of pressurizer c;πi,tAnd πj,tRespectively represent pressurizer c inlet end and outlet side
Air pressure;Γc、The respectively boosting ratio and maximum transfer capacity of pressurizer c.
Further, the operation constraint of energy hub described in step S2 are as follows:
The constraint of the energy hub thermal power equilibrium of supply and demand:
The unidirectional constraint of energy hub energy conversion:
Further, the present invention also provides a kind of power distribution network based on energy hub, natural gas grid establish it is comprehensive
Close the linear processing methods of energy resource system optimal operation model, comprising the following steps:
Step S3: linearization process is carried out to the model established in step S2;
Step S4: using CPLEX solver, to linearisation in step S3, treated that model solves.
Further, the particular content of the step S3 are as follows: will using second order cone relaxation and increment piecewise-linear techniques
Model conversation carries out linearization process;Due to containing nonlinear terms in node power Constraints of Equilibrium, using second order cone relaxation method
Abbreviation is carried out to constraint, new variables is introduced and eliminates voltage and current quadratic term, be shown below:
Node power Constraints of Equilibrium can eliminate the nonlinear quadratic item of voltage and current, as follows:
Voltage, current amplitude constraint modification as a result, are as follows:
The third formula of node power constraint is further converted to second order tapered by relaxation;Concrete form is as follows:
Line is carried out to the pipeline Weymouth steady-state load flow constraint in natural gas network using increment piecewise-linear techniques
Property, it is non-linear by the quadratic term bring for being introduced into air pressure in new variable cancellation pipeline Weymouth steady-state load flow constraint,
Shown in specific as follows:
1、
Therefore, pipeline Weymouth steady-state load flow constraint can tentatively rewrite are as follows:
Due to sgnp(πi,t,πj,t) it is a sign function, work as πi,tGreater than πj,tWhen take 1, otherwise take -1;Therefore, formula
The left side write as absolute value expression, even Yp,t=fp,t|fp,t|;That is Yp,t=fp,t|fp,t| it can also be expressed as f (x)=x | x
| non-linearization formula, the step of increment piece-wise linearization is as follows:
Step SA: number of segment n is averagely segmented according to model built setting;
Step SB: carrying out n equal part for the value interval of x, and this makes it possible to obtain n+1 discrete points, i.e. x0,x1...xn;
Step SC: discrete point x is calculated0,x1...xnCorresponding f (x) functional value;
Step SD: introducing new auxiliary variable, and guarantees to meet following constraint:
In conjunction with the above method, the constraint of pipeline Weymouth steady-state load flow is finally described as:
Compared with prior art, the invention has the following beneficial effects:
The present invention considers that power network, natural gas grid interact simultaneously, in the premise for meeting electricity, air and heat workload demand
Under, integrated energy system is more economically and safely run.Plan model has fully considered the characteristic flexibly converted between electricity, gas,
On the basis of meeting safe operation constraint, totle drilling cost is purchased as optimization aim using the energy of integrated energy system, promotes synthesis
Economy in energy system operation level.
Detailed description of the invention
Fig. 1 is the energy hub schematic diagram of the embodiment of the present invention.
Fig. 2 is the natural gas network rough schematic of the embodiment of the present invention.
Fig. 3 is the temporal characteristics figure of the electricity of the embodiment of the present invention, air and heat load.
Fig. 4 is the IEEE-14 node power distribution net topology of the embodiment of the present invention.
Fig. 5 is 20 node natural gas network topologies of the embodiment of the present invention.
Fig. 6 is purchase of electricity after the electricity price of the embodiment of the present invention, gas price curve and optimization, purchase tolerance curve.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and embodiments.
Present embodiments provide a kind of power distribution network based on energy hub, integrated energy system that natural gas grid is established it is excellent
Change moving model, comprising the following steps:
Step S1: establish it is a kind of can by energy hub that power distribution network, natural gas grid couple (energy hub,
EH) model completes integrated energy system modeling;
Step S2: considering the safe operation constraint of power distribution network, natural gas grid and energy hub, establish optimization aim,
Complete integrated energy system modeling.
The following are the specific implementation processes of the present embodiment.
(1) as shown in Fig. 1, establish it is a kind of can be by energy hub model that power distribution network, natural gas grid couple.
(2) using energy the lowest cost as target, the objective function of Optimized model is established.
(3) comprehensively consider power distribution network safe operation constraint.
(4) natural gas network reduction models and considers that natural gas Cybersecurity Operation constrains.
(5) comprehensively consider the safe operation constraint of energy hub.
In the present embodiment, terminal comprehensive energy model of element is established, step S1's is implemented as follows:
The energy demand of the forms such as electricity, heat, gas load is existed simultaneously in integrated energy system.To meet load need as far as possible
Ask, establish it is a kind of can be by energy hub (energy hub, EH) model that power distribution network and natural gas grid couple.The energy
Hub mainly by cogeneration of heat and power (combined heat and power, CHP), gas fired-boiler (gas furnace, GF) and
P2G (power to gas, P2G) device composition, based on the energy hub that this is established, not only between achievable electricity, gas
Mutually conversion, also can be supplied to thermic load.
In the present embodiment, step S2 is implemented as follows:
2, the objective function of Optimized model is established
The main target of the Optimized model of proposition is the total cost for making integrated energy system during dry run scheduling
It is minimum, therefore objective function mainly considers energy acquisition cost, including purchases strategies and purchase gas cost.
In formula, TR indicates the set of substation's injecting power node in power distribution network;GS indicates gas source node in natural gas grid
Set;Unit price of power, gas price in distribution expression period t;Function injects in substation in respectively period t
The purchase of electricity of rate node j;For the purchase tolerance of the gas source node j in period t.
3, consider power distribution network operation constraint:
Node power Constraints of Equilibrium
In formula, rij、xijThe respectively resistance of branch ij, reactance;δ (j) is indicated using j as the end of the branch of headend node
The set of node, π (j) are indicated using j as the set of the headend node of the branch of endpoint node;Pij,t、Qij,tRespectively branch ij's
Active power and reactive power;Major network power output and electric load active power are respectively indicated,It respectively indicates
Active power output for P2G device electric conversion in Hub;Indicate CHP device injection power distribution network in Hub
Equivalent active power,Respectively indicate the reactive power value of major network power output, electric load.
Node voltage amplitude constraint
Umin≤Uj,t≤Umax
In formula, Uj,t、UminAnd UmaxRespectively node voltage amplitude and its lower and upper limit.
The constraint of distribution line maximum carrying capacity
0≤Iij,t≤Imax
In formula, Iij,t、ImaxRespectively branch current magnitudes and its upper limit
Substation's power constraint
4, natural gas network model is established, and considers natural gas grid operation constraint
As shown in Fig. 2, natural gas network is mainly made of gas source, natural gas line, pressurizer and natural gas load.
Always mainly frictional dissipation inevitably occurs for gas flow in the flow process from exploitation to transport, by certain
After conveying distance (generally 80~160km), air pressure can be gradually decreased.To avoid excessive frictional dissipation, guarantee natural gas energy
It is transported to enough high-quality and high-efficiencies load side, the installation pressurizer on certain pipelines is generally required in natural gas system and is boosted,
To maintain node air pressure safety operation level.It can thus be seen that the effect of pressurizer is analogous in fact in natural gas network
Transformer in electric system.Natural gas grid network is constrained below:
The constraint of node air balance
In formula, σ (jk) indicates to respectively indicate the pipe using j as endpoint node by the pipeline set of headend node, μ (ij) of j
Road set;fp,t/fl,tFor the flow of pipeline p/l;The respectively gentle source node gas output of gas load,Respectively indicate the amount of natural gas for GF, CHP energy conversion in comprehensive energy unit (Hub);P2G device electricity is converted into the equivalent natural gas injection rate of gas in expression Hub.
The constraint of gas source gas output
The natural gas being produced from gas well needs the purification by refinery.Due to by gas at place capacity and gas well
The limitation of pressure, gas source gas output should meet following constraint in the unit time:
In formula,For the output amount of natural gas of gas source j,Respectively gas source j exports amount of natural gas
Bound.
Node air pressure constraint
The air pressure of each node of natural gas network must be maintained in safe, reasonable range of operation, and mathematic(al) representation is such as
Under:
In formula,Respectively represent air pressure upper lower limit value at node i.
The constraint of natural gas line Weymouth steady-state load flow
This project portrays the relationship of gas flow and both ends air pressure, specific table using Weymouth steady-state load flow model
Up to formula are as follows:
Wherein,
In addition pipeline flow also should ensure that within safe and reasonable range of operation:
In formula, fp,tRepresent the natural gas flow that pipeline p flows through;πi,tAnd πj,tRespectively indicate pipeline p both ends air pressure size;φp
Indicate the air stream transportation parameter of pipeline p;sgnpIndicate the air flow direction of pipeline p;Indicate maximum pipeline transmission capacity.
Pressurizer boosting relation constraint
Pressurizer model is simplified herein, has ignored the energy consumed when its operation, and only retains the air inlet of pressurizer
The transmission capacity constraint of boosting ratio relationship and pressurizer between mouth and gas outlet.
πj,t=Γcπi,t
In formula, fc,tIndicate the delivery air of pressurizer c;πi,tAnd πj,tRespectively represent pressurizer c inlet end and outlet side
Air pressure;Γc、The respectively boosting ratio and maximum transfer capacity of pressurizer c.
5, energy hub operation constraint
The constraint of the energy hub thermal power equilibrium of supply and demand
The unidirectional constraint of energy hub energy conversion
Preferably, the synthesis energy that the present embodiment also provides a kind of power distribution network based on energy hub, natural gas grid is established
The linear processing methods of source system optimized operation model, comprising the following steps:
Step S3: carrying out linearization process to built Optimized model, to solve integrated energy system optimization by CPLEX
Moving model.
Step S4: integrated energy system optimal operation model is solved by CPLEX.
I.e. are as follows: (6) linearization process is carried out to built Optimized model.
(7) integrated energy system optimal operation model is solved by CPLEX.
Step S3's is implemented as follows:
The model that step S2 is proposed is a comprehensive non-convex problem.In order to solve this problem, before solving model,
First linear problem is converted by master mould with increment piecewise-linear techniques using second order cone is loose.
6, model linearization
1) second order cone relaxing techniques
Due to containing nonlinear terms in node power Constraints of Equilibrium, in order to solve this problem, using second order cone relaxing techniques
Abbreviation is carried out to constraint.It introduces new variables and eliminates voltage and current quadratic term, be shown below:
Node power Constraints of Equilibrium can eliminate the nonlinear quadratic item of voltage and current, as follows:
Voltage, current amplitude constraint can also be modified as a result, are as follows:
The third formula of node power constraint can be further converted to second order tapered by relaxation.The following institute of concrete form
Show:
2) increment piecewise-linear techniques
Pipeline Weymouth steady-state load flow constraint in natural gas network is a more intractable Nonlinear Nonconvex expression
Formula, it is difficult to which therefore direct solution also needs to linearize it.
It is possible, firstly, to pass through the quadratic term band for being introduced into air pressure in new variable cancellation pipeline Weymouth steady-state load flow constraint
That comes is non-linear, specific as follows shown:
Therefore, pipeline Weymouth steady-state load flow constraint can tentatively rewrite are as follows:
Due to sgnp(πi,t,πj,t) it is a sign function, work as πi,tGreater than πj,tWhen take 1, otherwise take -1.It therefore, can be
Absolute value expression is write as in the left side of formula, even Yp,t=fp,t|fp,t|。
Such as f (x)=x | x | this kind of non-linearization formula, the step of increment piece-wise linearization approximately as:
Step1: the characteristics of according to model built, suitable segmentation number of segment n is set;
Step2: carrying out n equal part for the value interval of x, and this makes it possible to obtain n+1 discrete points, i.e. x0,x1...xn;
Step3: discrete point x is calculated0,x1...xnCorresponding f (x) functional value;
Step4: introducing new auxiliary variable, and guarantees to meet following constraint:
In conjunction with the above method, the constraint of pipeline Weymouth steady-state load flow can finally be described as:
So far, all nonlinear terms in model are converted into linear term, can directly be asked by CPLEX solver
Solution.
CPlEX is that one kind can be to the tool that Optimized model effectively solves.First according to required majorization of solutions model,
Program is write under MATLAB environment, must define decision variable, objective function and linear constraint condition, i.e. root respectively in program
The model content solved according to foregoing description;The CPLEX being then based in the tool box this routine call YALMIP is had
Effect solves.
First by matlab environment, program is write based on the tool box yalmip, solution procedure passes through yalmip work
Having case calls CPLEX optimization solver to model solution.
It is illustrated below by specific example:
The present embodiment uses a comprehensive energy system being coupled to form by 14 node power distribution nets and 20 node natural gas networks
System carries out simulation analysis.Power distribution network has 1 power transformation tiny node, 7 electric load nodes, 8 nodes and 16 lines coupled
Road.Natural gas network has 2 gas source nodes, 10 gas load bus and 8 switching nodes and 25 pipelines.It is worth noting
, it can be also its supply thermal power that 8 switching nodes, which can not only meet electricity, gas workload demand,.The present embodiment passes through simulation one
The operation of 24 hours in it, with the minimum target of day operation expense, to obtain economic scheduling scheme.
Electric (Fig. 3 a), gas (Fig. 3 b), thermic load (Fig. 3 c) timing curve as shown in figure 3, related example topological diagram such as Fig. 4
With shown in Fig. 5, Hub node coupling information is as shown in table 1.
Each Hub cell node corresponding relationship of table 1 and load parameter
In conjunction with the corresponding example of the present embodiment, in power distribution network, node voltage amplitude (pu) bound is respectively 1.05 Hes
0.95.In natural gas grid, pipeline air pressure bound is respectively 60 and 50bar.Simultaneously, it is assumed that each each energy conversion apparatus of Hub
Efficiency it is identical, P2G 75%, cogeneration efficiency of fuel cell generation is 37%, heat production 37%, and boiler heat production efficiency is 96%, gas
The conversion coefficient for turning electricity is 0.0096MW/m3。
Purchase of electricity, purchase tolerance curve after the electricity price of the embodiment of the present invention, gas price curve and optimization is as shown in Figure 6.It can
See, the integrated energy system optimization of the power distribution network based on energy hub, natural gas grid foundation that are proposed according to the present embodiment
Moving model and linear processing methods, the electricity price shown in Fig. 6 (a), gas price change under influence, which will
It rationally responds, with decision reasonably purchase of electricity and purchase tolerance (such as Fig. 6 (b) shown in), minimizes energy purchase to realize
Cost.
The foregoing is merely presently preferred embodiments of the present invention, all equivalent changes done according to scope of the present invention patent with
Modification, is all covered by the present invention.
Claims (7)
1. the integrated energy system optimal operation model that a kind of power distribution network based on energy hub, natural gas grid are established, special
Sign is: the following steps are included:
Step S1: providing one can be by energy hub model that power distribution network, natural gas grid couple;
Step S2: establishing objective function, and constrained according to the safe operation of power distribution network, natural gas grid and energy hub, complete
At the modeling of integrated energy system optimal operation model.
2. the comprehensive energy system that a kind of power distribution network based on energy hub according to claim 1, natural gas grid are established
System optimal operation model, it is characterised in that: objective function described in step S2 is that energy acquisition cost is shown below, including purchases
Electric costWith purchase gas costSimulating integrated energy system
Total cost during traffic control is minimum;
In formula, TR indicates the set of substation's injecting power node in power distribution network;GS indicates the collection of gas source node in natural gas grid
It closes;Unit price of power, gas price in distribution expression period t;Substation's injecting power section in respectively period t
The purchase of electricity of point j;For the purchase tolerance of the gas source node j in period t.
3. a kind of power distribution network based on energy hub according to claim 1, the integrated energy system of natural gas grid are excellent
Change moving model, it is characterised in that: the safe operation constraint of power distribution network described in step S2 includes: node power Constraints of Equilibrium, section
The constraint of point voltage magnitude, the constraint of distribution line maximum carrying capacity and substation's power constraint;
The node power Constraints of Equilibrium are as follows:
In formula, rij、xijThe respectively resistance of branch ij, reactance;δ (j) is indicated using j as the endpoint node of the branch of headend node
Set, π (j) are indicated using j as the set of the headend node of the branch of endpoint node;Pij,t、Qij,tThe respectively wattful power of branch ij
Rate and reactive power;Respectively indicate major network power output and electric load active power;It respectively indicates in Hub
The active power of P2G device electric conversion exports;Indicate the equivalent wattful power of CHP device injection power distribution network in Hub
Rate;Respectively indicate the reactive power value of major network power output, electric load;
The node voltage amplitude constraint are as follows:
Umin≤Uj,t≤Umax
In formula, Uj,t、UminAnd UmaxRespectively node voltage amplitude and its lower and upper limit;
The distribution line maximum carrying capacity constraint are as follows:
0≤Iij,t≤Imax
In formula, Iij,t、ImaxRespectively branch current magnitudes and its upper limit;
Substation's power constraint are as follows:
4. a kind of power distribution network based on energy hub according to claim 1, the integrated energy system of natural gas grid are excellent
Change moving model, it is characterised in that: the particular content constrained in step S2 the natural gas network are as follows:
The constraint of node air balance:
In formula, σ (jk) is indicated using j as the pipeline set of headend node;μ (ij) respectively indicates the pipeline collection using j as endpoint node
It closes;fp,t/fl,tFor the flow of pipeline p/l;The respectively gentle source node gas output of gas load;
Respectively indicate the amount of natural gas converted for GF, CHP energy in Hub;Indicate P2G device electricity conversion in Hub
At the equivalent natural gas injection rate of gas;
The constraint of gas source gas output:
The natural gas being produced from gas well needs the purification by refinery, due to by air pressure at place capacity and gas well
Limitation, gas source gas output should meet following constraint in the unit time:
In formula,For the output amount of natural gas of gas source j,Above and below respectively gas source j output amount of natural gas
Limit;
Node air pressure constraint:
The air pressure of each node of natural gas network must be maintained in safe, reasonable range of operation, and mathematic(al) representation is as follows:
In formula,Respectively represent air pressure upper lower limit value at node i;
The constraint of natural gas line Weymouth steady-state load flow:
The relationship of gas flow and both ends air pressure, expression are portrayed using Weymouth steady-state load flow model are as follows:
Wherein,
In addition pipeline flow also should ensure that within safe and reasonable range of operation:
In formula, fp,tRepresent the natural gas flow that pipeline p flows through;πi,tAnd πj,tRespectively indicate pipeline p both ends air pressure size;φpIt indicates
The air stream transportation parameter of pipeline p;sgnpIndicate the air flow direction of pipeline p;Indicate maximum pipeline transmission capacity.
Pressurizer boosting relation constraint:
πj,t=Γcπi,t
In formula, fc,tIndicate the delivery air of pressurizer c;πi,tAnd πj,tRespectively represent the gas of pressurizer c inlet end and outlet side
Pressure;Γc、The respectively boosting ratio and maximum transfer capacity of pressurizer c.
5. a kind of power distribution network based on energy hub according to claim 1, the integrated energy system of natural gas grid are excellent
Change moving model, it is characterised in that: the operation of energy hub described in step S2 constrains are as follows:
The constraint of the energy hub thermal power equilibrium of supply and demand:
The unidirectional constraint of energy hub energy conversion:
6. a kind of a kind of power distribution network based on energy hub according to any one of claims 1 to 5, natural gas grid are established
Integrated energy system optimal operation model linear processing methods, which comprises the following steps:
Step S3: linearization process is carried out to the model established in step S2;
Step S4: using CPLEX solver, to linearisation in step S3, treated that model solves.
7. the comprehensive energy system that a kind of power distribution network based on energy hub according to claim 6, natural gas grid are established
The linear processing methods of system optimal operation model, it is characterised in that: the particular content of the step S3 are as follows: using second order cone relaxation
Model conversation is subjected to linearization process with increment piecewise-linear techniques;
Due to containing nonlinear terms in node power Constraints of Equilibrium, abbreviation is carried out to constraint using second order cone relaxation method, is introduced
New variables eliminates voltage and current quadratic term, is shown below:
Node power Constraints of Equilibrium can eliminate the nonlinear quadratic item of voltage and current, as follows:
Voltage, current amplitude constraint modification as a result, are as follows:
The third formula of node power constraint is further converted to second order tapered by relaxation;Concrete form is as follows:
The pipeline Weymouth steady-state load flow constraint in natural gas network is linearized using increment piecewise-linear techniques,
It is non-linear by the quadratic term bring for being introduced into air pressure in new variable cancellation pipeline Weymouth steady-state load flow constraint, specifically such as
Shown in lower:
Therefore, pipeline Weymouth steady-state load flow constraint can tentatively rewrite are as follows:
Due to sgnp(πi,t,πj,t) it is a sign function, work as πi,tGreater than πj,tWhen take 1, otherwise take -1;Therefore, the left side of formula
Write as absolute value expression, even Yp,t=fp,t|fp,t|;That is Yp,t=fp,t|fp,t| can also be expressed as f (x)=x | x | it is non-
The step of linearisation formula, increment piece-wise linearization, is as follows:
Step SA: number of segment n is averagely segmented according to model built setting;
Step SB: carrying out n equal part for the value interval of x, and this makes it possible to obtain n+1 discrete points, i.e. x0,x1...xn;
Step SC: discrete point x is calculated0,x1...xnCorresponding f (x) functional value;
Step SD: introducing new auxiliary variable, and guarantees to meet following constraint:
In conjunction with the above method, the constraint of pipeline Weymouth steady-state load flow is finally described as:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910456160.XA CN110263387B (en) | 2019-05-29 | 2019-05-29 | Energy system operation device based on power distribution network and natural gas network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910456160.XA CN110263387B (en) | 2019-05-29 | 2019-05-29 | Energy system operation device based on power distribution network and natural gas network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110263387A true CN110263387A (en) | 2019-09-20 |
CN110263387B CN110263387B (en) | 2023-04-07 |
Family
ID=67915690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910456160.XA Active CN110263387B (en) | 2019-05-29 | 2019-05-29 | Energy system operation device based on power distribution network and natural gas network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110263387B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110689206A (en) * | 2019-10-09 | 2020-01-14 | 国电南瑞南京控制系统有限公司 | Energy internet multi-energy main body energy purchasing and conversion operation method |
CN111401647A (en) * | 2020-03-23 | 2020-07-10 | 清华大学 | Distributed optimal scheduling method for electric coupling system considering uncertainty transfer |
CN111414721A (en) * | 2020-02-22 | 2020-07-14 | 清华大学 | Heat supply network waterway modeling method for comprehensive energy system scheduling |
CN111681130A (en) * | 2020-06-15 | 2020-09-18 | 西安交通大学 | Comprehensive energy system optimization scheduling method considering condition risk value |
CN111768036A (en) * | 2020-06-29 | 2020-10-13 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid |
CN111950122A (en) * | 2020-07-08 | 2020-11-17 | 国网(苏州)城市能源研究院有限责任公司 | Operation optimization method for park comprehensive energy system |
CN112016195A (en) * | 2020-08-10 | 2020-12-01 | 浙江大学 | Flexible planning method of electrical coupling system considering energy supply reliability |
CN112069634A (en) * | 2020-08-14 | 2020-12-11 | 广东工业大学 | Gas network capacity expansion planning system and method based on relaxation constraint |
CN112861292A (en) * | 2021-01-12 | 2021-05-28 | 浙江大学 | Recovery improvement method for electricity-gas comprehensive energy system |
CN113537618A (en) * | 2021-07-29 | 2021-10-22 | 中国电建集团河南省电力勘测设计院有限公司 | Comprehensive energy system optimization scheduling method considering resident user demand response |
CN113570117A (en) * | 2021-07-02 | 2021-10-29 | 浙江华云电力工程设计咨询有限公司 | Electricity-gas comprehensive energy microgrid optimal scheduling method based on second-order cone relaxation conversion method |
CN113592149A (en) * | 2021-07-01 | 2021-11-02 | 四川大学 | Optimization and reconstruction method for coupled new energy comprehensive energy power distribution network |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105703368A (en) * | 2016-02-04 | 2016-06-22 | 马瑞 | Multiple uncertain energy flow modeling and calculation method for integrated system of active distribution network (ADN) and power transmission network under energy interconnection |
CN107506851A (en) * | 2017-07-26 | 2017-12-22 | 河海大学 | A kind of multizone virtual plant comprehensive energy coordinated scheduling Optimized model |
CN107769215A (en) * | 2018-01-19 | 2018-03-06 | 国网天津市电力公司 | Garden energy mix system optimization dispatching method based on energy hub |
WO2018059096A1 (en) * | 2016-09-30 | 2018-04-05 | 国电南瑞科技股份有限公司 | Combined decision method for power generation plans of multiple power sources, and storage medium |
-
2019
- 2019-05-29 CN CN201910456160.XA patent/CN110263387B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105703368A (en) * | 2016-02-04 | 2016-06-22 | 马瑞 | Multiple uncertain energy flow modeling and calculation method for integrated system of active distribution network (ADN) and power transmission network under energy interconnection |
WO2018059096A1 (en) * | 2016-09-30 | 2018-04-05 | 国电南瑞科技股份有限公司 | Combined decision method for power generation plans of multiple power sources, and storage medium |
CN107506851A (en) * | 2017-07-26 | 2017-12-22 | 河海大学 | A kind of multizone virtual plant comprehensive energy coordinated scheduling Optimized model |
CN107769215A (en) * | 2018-01-19 | 2018-03-06 | 国网天津市电力公司 | Garden energy mix system optimization dispatching method based on energy hub |
Non-Patent Citations (1)
Title |
---|
倪伟 等: "基于机会约束规划的能源集线器系统气电购置优化建模", 《电网技术》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110689206A (en) * | 2019-10-09 | 2020-01-14 | 国电南瑞南京控制系统有限公司 | Energy internet multi-energy main body energy purchasing and conversion operation method |
CN110689206B (en) * | 2019-10-09 | 2022-11-01 | 国电南瑞南京控制系统有限公司 | Energy Internet multi-energy main body energy purchasing and converting operation method |
CN111414721A (en) * | 2020-02-22 | 2020-07-14 | 清华大学 | Heat supply network waterway modeling method for comprehensive energy system scheduling |
CN111401647B (en) * | 2020-03-23 | 2022-04-08 | 清华大学 | Distributed optimal scheduling method for electric coupling system considering uncertainty transfer |
CN111401647A (en) * | 2020-03-23 | 2020-07-10 | 清华大学 | Distributed optimal scheduling method for electric coupling system considering uncertainty transfer |
CN111681130A (en) * | 2020-06-15 | 2020-09-18 | 西安交通大学 | Comprehensive energy system optimization scheduling method considering condition risk value |
CN111681130B (en) * | 2020-06-15 | 2024-04-16 | 西安交通大学 | Comprehensive energy system optimal scheduling method considering conditional risk value |
CN111768036A (en) * | 2020-06-29 | 2020-10-13 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid |
CN111768036B (en) * | 2020-06-29 | 2023-11-03 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid |
CN111950122A (en) * | 2020-07-08 | 2020-11-17 | 国网(苏州)城市能源研究院有限责任公司 | Operation optimization method for park comprehensive energy system |
CN112016195A (en) * | 2020-08-10 | 2020-12-01 | 浙江大学 | Flexible planning method of electrical coupling system considering energy supply reliability |
CN112069634A (en) * | 2020-08-14 | 2020-12-11 | 广东工业大学 | Gas network capacity expansion planning system and method based on relaxation constraint |
CN112861292A (en) * | 2021-01-12 | 2021-05-28 | 浙江大学 | Recovery improvement method for electricity-gas comprehensive energy system |
CN113592149A (en) * | 2021-07-01 | 2021-11-02 | 四川大学 | Optimization and reconstruction method for coupled new energy comprehensive energy power distribution network |
CN113570117A (en) * | 2021-07-02 | 2021-10-29 | 浙江华云电力工程设计咨询有限公司 | Electricity-gas comprehensive energy microgrid optimal scheduling method based on second-order cone relaxation conversion method |
CN113570117B (en) * | 2021-07-02 | 2024-02-09 | 浙江华云电力工程设计咨询有限公司 | Electric-gas comprehensive energy microgrid optimal scheduling method based on second order cone relaxation conversion method |
CN113537618A (en) * | 2021-07-29 | 2021-10-22 | 中国电建集团河南省电力勘测设计院有限公司 | Comprehensive energy system optimization scheduling method considering resident user demand response |
Also Published As
Publication number | Publication date |
---|---|
CN110263387B (en) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110263387A (en) | The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established | |
CN108173282B (en) | A kind of consideration electricity turns gas operating cost integrated energy system Optimization Scheduling | |
CN110110913B (en) | Optimal configuration method for energy station of large-scale park comprehensive energy system | |
CN110163443B (en) | Natural gas pressure regulating station micro-energy network optimization scheduling method considering electricity-gas comprehensive demand response | |
CN110288152B (en) | Regional comprehensive energy system energy storage configuration method considering electric/thermal flexible load | |
CN105576710B (en) | A kind of configuration method of distributed generation resource in integrated energy system | |
CN109494809B (en) | Electric-gas interconnection system comprising electric gas conversion device and operation method | |
CN109978362A (en) | A kind of modeling of gas power grid joint multizone integrated energy system and systems organization method | |
CN111815068B (en) | Optimization method for solving urban comprehensive energy network by two-stage constraint boundary tightening algorithm | |
CN109740955A (en) | A kind of electric-gas integrated energy system planing method counted and improve staged carbon transaction | |
CN111799777A (en) | Comprehensive energy planning method considering coupling of natural gas and electric power | |
CN113315242A (en) | Virtual wind abandoning-hydrogen production combination for promoting wind abandoning consumption based on hydrogen energy economy | |
CN107067108A (en) | Meter and the garden integrated energy system energy efficiency evaluating method of regenerative resource access | |
CN109685291A (en) | Based on the regional complex energy resource system provided multiple forms of energy to complement each other method for optimization analysis a few days ago | |
Wei et al. | Energy-efficient distribution in smart grid | |
CN113128868B (en) | Regional comprehensive energy system scheduling optimization method and device | |
CN113806952B (en) | Cold-hot electricity comprehensive energy system considering source-charge-storage and optimal operation method thereof | |
Liu et al. | Low-carbon economic dispatch of integrated electricity-gas energy system considering carbon capture, utilization and storage | |
CN111476394B (en) | Robust operation optimization method suitable for multi-energy systems such as electric heating gas system | |
CN110377973B (en) | Construction method of standard linear comprehensive energy system model | |
CN109066695A (en) | A kind of electrical optimal energy flux computation method of two stages linearisation | |
CN107248032A (en) | A kind of comprehensive assessment index system construction method of regional complex energy resource system | |
Liu et al. | Influence Evaluation of Integrated Energy System on the Unit Commitment in Power System | |
CN110516868B (en) | Comprehensive energy system optimization operation model considering network constraints | |
CN110070216A (en) | A kind of industrial park integrated energy system economic operation optimization method |
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 |