CN107491849A - Meter and the plant gas gas storage capacity calculation methods of electrical couplings system restriction - Google Patents
Meter and the plant gas gas storage capacity calculation methods of electrical couplings system restriction Download PDFInfo
- Publication number
- CN107491849A CN107491849A CN201710927821.3A CN201710927821A CN107491849A CN 107491849 A CN107491849 A CN 107491849A CN 201710927821 A CN201710927821 A CN 201710927821A CN 107491849 A CN107491849 A CN 107491849A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- gas
- natural gas
- node
- 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
- 238000003860 storage Methods 0.000 title claims abstract description 71
- 238000004364 calculation method Methods 0.000 title claims abstract description 9
- 230000008878 coupling Effects 0.000 title claims abstract description 6
- 238000010168 coupling process Methods 0.000 title claims abstract description 6
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 6
- 239000007789 gas Substances 0.000 claims abstract description 187
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 174
- 239000003345 natural gas Substances 0.000 claims abstract description 87
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 239000003245 coal Substances 0.000 claims description 31
- 239000011159 matrix material Substances 0.000 claims description 12
- 239000000567 combustion gas Substances 0.000 claims description 7
- 239000003795 chemical substances by application Substances 0.000 claims description 6
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000002485 combustion reaction Methods 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 239000002699 waste material Substances 0.000 abstract description 3
- 230000005611 electricity Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
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/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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
Counted and the plant gas gas storage capacity calculation methods of electrical couplings system restriction the present invention relates to a kind of, belong to the plan optimization technical field of multipotency stream coupled system.For this method using the method optimized, the capacity that gas storage equipment is configured for plant gas provides preferred plan.On the one hand power system and the operation of natural gas grid constraint, the potential safety hazard for avoiding the gas storage scheme of low capacity from being brought to power system are considered simultaneously;On the other hand the Cost Problems of gas storage equipment are considered, avoid economic waste caused by the gas storage scheme of excess capacity.This method can apply in the planning of plant gas gas storage facility, to manage the configuration foundation that gas storage equipment capacity is provided with designer, be favorably improved the security and economy of plant gas power supply.
Description
Technical field
The present invention relates to a kind of meter and the plant gas gas storage capacity calculation methods of electric-gas coupled system constraint, belong to more
The plan optimization technical field of coupled system can be flowed.
Background technology
Because the cost of gas electricity generator is relatively low, fast response time smaller to the destructive power of environment, the construction week of gas station
The huge advantages such as the phase is shorter, worldwide turn into the important component of power supply energy.And as plant gas is in electricity
Ratio in Force system energy resource supply increases sharply, and the degree of coupling of power system and natural gas system is gradually deepened, power train
The dependence that the safety and stability of system supplies safety to natural gas system is also deepened therewith.
On the other hand, fluctuation of the natural gas load between year, month, day is very violent, and in many national regulations, its
He business, civil natural gas load priority be higher than plant gas Gas Load.If there is spike in other natural gas loads
When, due to the security constraint of natural gas system, the Gas Load of plant gas may be caused insufficient, therefore in plant gas
Configuration certain capacity natural gas gas storage facility be a kind of common and effectively to improve Operation of Electric Systems safety side
Method.However, gas storage facility cost be usually proportionate with its capacity, therefore how reasonably to configure gas storage equipment capacity into
For urgent problem to be solved.
The content of the invention
The purpose of the present invention is to propose to count and electric-gas coupled system constraint plant gas gas storage capacity calculation methods, with
Avoid the power system energy supply deficiency caused by gas storage capacity configuration is too small, or the economy caused by gas storage capacity configuration is excessive
Waste.
The plant gas gas storage capacity calculation methods of meter and electric-gas coupled system constraint proposed by the present invention, including it is following
Step:
(1) optimization object function for the plant gas cost of electricity-generating established in an electric-gas coupled system is as follows:
Ccoal(Gcoal)=Pcoal×Gcoal
Cgas(Lgas)=Pgas×Lgas
Wherein, C be Operation of Electric Systems totle drilling cost, Ccoal(Gcoal) it is coal-fired cost, GcoalFor coal consumption, PcoalFor
The price of coal, Cgas(Lgas) it is combustion gas cost, LgasFor air consumption, PgasFor the purchasing price of natural gas, Cstorage(V) it is storage
Gas cost, V be gas storage equipment volume, TlifeFor the service life of gas storage equipment;
(2) the equality constraint equation of electric-gas coupled system steady state Safe Operation is established, including:
The power flow equation of power system is as follows in (2-1) electric-gas coupled system:
Wherein, PiFor the injection active power of i-th of node in power system, QiFor the note of i-th of node in power system
Enter reactive power, GijFor the corresponding conductance of the i-th row, jth row in the bus admittance matrix Y with power system, BijFor with electric power
I-th row in the bus admittance matrix Y of system, jth arrange corresponding susceptance, and power system bus admittance matrix Y is from power system
Control centre obtains;
The waterpower equation of pipeline is as follows in natural gas grid in (2-2) electric-gas coupled system:
Wherein, fkmFor natural gas volume flow in the pipeline in natural gas grid between k-th of node and m-th of node, pk,
pmThe respectively pressure of k-th of node and m-th of node, DkmAnd LkmPipeline km between respectively k-th of node and m-th of node
Pipe diameter and duct length, F be inner-walls of duct coefficient of friction, F is by formulaCalculate
Arrive, EfFor the efficiency factor of pipeline, EfValue is that 0.92, Re is Reynolds number, by formulaIt is calculated, ρ
For natural gas density, μmFor the coefficient of kinetic viscosity (relevant with natural gas density, can table look-up acquisition) of natural gas, γGFor natural gas
Proportion, 0 < γG< 1, TaFor natural gas mean temperature, TnAnd pnThe temperature and pressure of natural gas, Z respectively under standard stategFor
The average compressible coefficient of natural gas, 0.9 < Zg< 1.5, in above-mentioned natural gas grid in the waterpower equation of pipeline, when
When, the sgn in above formulap(pk,pm)=1, whenWhen, sgnp(pk,pm)=- 1;
The energy expenditure equation of compressor is as follows in the natural gas grid of (2-3) electric-gas coupled system:
Wherein, pk, pmThe pressure of k-th of node and m-th of node respectively in natural gas grid, BHPkmFor k-th of node
The energy expenditure of compressor between m-th of node,For the inlet volumetric flow of the compressor, ηcFor the total of the compressor
Efficiency, ckFor the polytropic coefficient of the compressor, ηcAnd ckObtained from the shop instructions of compressor;
Coupling between the power system coupled in (2-4) electric-gas coupled system by gas turbine and natural gas grid
It is as follows to close equation:
μPTur=Hgas×[fTur+(Vt-1-Vt)],
Wherein, fTurThe volume flow of the combustion gas inputted for gas turbine by pipeline, PTurFor the active power of gas turbine
Output, VtIt is gas storage equipment in the gas reserves of t, HgasFor the combustion heat value of natural gas, value 37.59MJ/m3, μ
For the efficiency factor of gas turbine, obtained by the shop instructions of gas turbine;
The node air balance equation of natural gas grid is as follows in (2-5) electric-gas coupled system:
AGF=L,
Wherein, AGFor node-branch road matrix of natural gas grid, f is the volume flow of natural gas grid branch road, and L is natural gas
The gas load of net node, L obtain according to natural gas grid history data;
(3) the inequality constraints condition of electric-gas coupled system steady state Safe Operation is set, including:
The power output P of (3-1) generators in power systems groupi genMore than or equal to 0, less than or equal to the generating set
Dispatch from the factory the peak power provided on nameplateI.e.:
The voltage magnitude U of i-th of node of (3-2) power systemiIn the upper limit of the safe operation of power system voltage of setting
ValueWith lower limit UiBetween run, UiFor 0.95 times of i-th of node rated voltage,For i-th node rated voltage
1.05 times, i.e.,:
The transmission capacity S of the l articles circuit in (3-3) power systemlLess than or equal to the safe operation of power system of setting
The maximum of transmission capacityI.e.:
The pressure p of k-th of node in (3-4) natural gas gridkThe higher limit of air pressure, lower limit are run in the pipe safety of setting
Valuep k、It is interior, i.e.,:
The flow f of b-th of pipeline in (3-5) natural gas gridbHigher limit, lower limit in the pipe safety operating flux of setting
Valuef b、It is interior, i.e.,:
Gas source feed amount f in (3-6) natural gas gridsLess than or equal to the maximum that the source of the gas can provide natural gas flow
fs,max, i.e.,:
fs≤fs,max;
The constraints of compressor safe operation in (3-7) natural gas grid:
Wherein:S is the step-up ratio of the compressor, SmaxIt is the maximum step-up ratio of the compressor, SmaxBy dispatching from the factory for compressor
Nameplate obtains,For the volume flow of the entrance of the compressor,For the maximum allowable volume flow of entrance of the compressor,Obtained by the nameplate that dispatches from the factory of compressor, poutFor the outlet pressure of compressor, pc,maxIt is maximum allowable for the outlet of compressor
Pressure, pc,maxObtained by the nameplate that dispatches from the factory of compressor;
Gas storage equipment should be greater than being equal to 0 in the reserves of t in (3-8) plant gas, less than or equal to gas storage equipment most
Big volume:
0≤Vt≤Vmax;
(4) optimized algorithm is utilized, it is that object function, above-mentioned steps (2) and step (3) are about to solve by above-mentioned steps (1)
The Optimized model of beam condition, obtain the moon optimal gas storage capacity curve of the gas storage equipment under prediction year load, optimal gas storage capacity
Abscissa in curve is 12 months, and the ordinate in optimal gas storage capacity curve is moon gas storage capacity;
(5) according to the optimal gas storage capacity curve of above-mentioned steps (4), the desired value of capacity is asked for, will be closest to desired value
Optimal capacity of the capacity as gas storage equipment in plant gas.
It is proposed by the present invention meter and electric-gas coupled system constraint plant gas gas storage capacity calculation methods, its feature and
Effect is:For the present invention using the method optimized, the capacity that gas storage equipment is configured for plant gas provides preferred plan.On the one hand
Power system and the operation of natural gas grid constraint are considered simultaneously, avoid what the gas storage scheme of low capacity was brought to power system
Potential safety hazard;On the other hand the Cost Problems of gas storage equipment are considered, avoid economy caused by the gas storage scheme of excess capacity
Waste.This method be can apply in the planning of plant gas gas storage facility, and gas storage equipment is provided for management and designer
The configuration foundation of capacity, it is favorably improved the security and economy of plant gas power supply.
Brief description of the drawings
Fig. 1 is the structure rough schematic view for the electric-gas coupled system that the inventive method is related to.
Embodiment
The plant gas gas storage capacity calculation methods of meter and electric-gas coupled system constraint proposed by the present invention, are directed to
Electric-gas coupled system structural representation as shown in figure 1, this method comprises the following steps:
(1) optimization object function for the plant gas cost of electricity-generating established in an electric-gas coupled system is as follows:
Ccoal(Gcoal)=Pcoal×Gcoal
Cgas(Lgas)=Pgas×Lgas
Wherein, C be Operation of Electric Systems totle drilling cost, Ccoal(Gcoal) it is coal-fired cost, GcoalFor coal consumption, (unit is
Ton), PcoalFor the price (unit for yuan/ton, coal price decision is sold by local then) of coal, Cgas(Lgas) for combustion gas into
This, LgasFor air consumption (unit is mark cubic meter), PgasFor the purchasing price of natural gas, (unit is member/mark cubic meter, by then
Local sells gas price decision), Cstorage(V) (provided for gas storage cost by production and installation producer), V is the appearance of gas storage equipment
Product, can be by obtaining in the shop instructions of gas storage equipment, Tlife, can be by gas storage equipment for the service life of gas storage equipment
Product nameplate obtains;
(2) the equality constraint equation of electric-gas coupled system steady state Safe Operation is established, including:
The power flow equation of power system is as follows in (2-1) electric-gas coupled system:
Wherein, PiFor the injection active power of i-th of node in power system, QiFor the note of i-th of node in power system
Enter reactive power, GijFor the corresponding conductance of the i-th row, jth row in the bus admittance matrix Y with power system, BijFor with electric power
I-th row in the bus admittance matrix Y of system, jth arrange corresponding susceptance, and power system bus admittance matrix Y is from power system
Control centre obtains;
The waterpower equation of pipeline is as follows in natural gas grid in (2-2) electric-gas coupled system:
Wherein, fkmFor natural gas volume flow in the pipeline in natural gas grid between k-th of node and m-th of node, pk,
pmThe respectively pressure of k-th of node and m-th of node, DkmAnd LkmPipeline km between respectively k-th of node and m-th of node
Pipe diameter and duct length, F be inner-walls of duct coefficient of friction, F is by formulaCalculate
Arrive, EfFor the efficiency factor of pipeline, EfValue is that 0.92, Re is Reynolds number, by formulaIt is calculated, ρ
For natural gas density, μmFor the coefficient of kinetic viscosity (relevant with natural gas density, can table look-up acquisition) of natural gas, γGFor natural gas
Proportion, 0 < γG< 1, TaFor natural gas mean temperature, TnAnd pnThe temperature and pressure of natural gas respectively under standard state,
In one embodiment of the present of invention, TnAnd pnValue be respectively 288K and 0.1Mpa, ZgFor the average compressible coefficient of natural gas,
0.9 < Zg< 1.5, in above-mentioned natural gas grid in the waterpower equation of pipeline, whenWhen, the sgn in above formulap(pk,
pm)=1, whenWhen, sgnp(pk,pm)=- 1;
The energy expenditure equation of compressor is as follows in the natural gas grid of (2-3) electric-gas coupled system:
Wherein, pk, pmThe pressure of k-th of node and m-th of node respectively in natural gas grid, BHPkmFor k-th of node
The energy expenditure of compressor between m-th of node,For the inlet volumetric flow of the compressor, ηcFor the total of the compressor
Efficiency, ckFor the polytropic coefficient of the compressor, ηcAnd ckObtained from the shop instructions of compressor;
Coupling between the power system coupled in (2-4) electric-gas coupled system by gas turbine and natural gas grid
It is as follows to close equation:
μPTur=Hgas×[fTur+(Vt-1-Vt)],
Wherein, fTurThe volume flow of the combustion gas inputted for gas turbine by pipeline, PTurFor the active power of gas turbine
Output, VtIt is gas storage equipment in the gas reserves of t, HgasFor the combustion heat value of natural gas, value 37.59MJ/m3, μ
For the efficiency factor of gas turbine, obtained by the shop instructions of gas turbine;
The node air balance equation of natural gas grid is as follows in (2-5) electric-gas coupled system:
AGF=L,
Wherein, AGFor node-branch road matrix of natural gas grid, f is the volume flow of natural gas grid branch road, and L is natural gas
The gas load of net node, L obtain according to natural gas grid history data;
(3) the inequality constraints condition of electric-gas coupled system steady state Safe Operation is set, including:
The power output P of (3-1) generators in power systems groupi genMore than or equal to 0, less than or equal to the generating set
Dispatch from the factory the peak power provided on nameplateI.e.:
The voltage magnitude U of i-th of node of (3-2) power systemiIn the upper limit of the safe operation of power system voltage of setting
ValueAnd lower limitU iBetween run,U iFor 0.95 times of i-th of node rated voltage,For i-th node rated voltage
1.05 times, i.e.,:
The transmission capacity S of the l articles circuit in (3-3) power systemlLess than or equal to the safe operation of power system of setting
The maximum of transmission capacityI.e.:
The pressure p of k-th of node in (3-4) natural gas gridkThe higher limit of air pressure, lower limit are run in the pipe safety of setting
Valuep k、It is interior, i.e.,:
The flow f of b-th of pipeline in (3-5) natural gas gridbHigher limit, lower limit in the pipe safety operating flux of setting
Valuef b、It is interior, i.e.,:
Gas source feed amount f in (3-6) natural gas gridsLess than or equal to the maximum that the source of the gas can provide natural gas flow
fs,max, i.e.,:
fs≤fs,max;
The constraints of compressor safe operation in (3-7) natural gas grid:
Wherein:S is the step-up ratio of the compressor, SmaxIt is the maximum step-up ratio of the compressor, SmaxBy dispatching from the factory for compressor
Nameplate obtains,For the volume flow of the entrance of the compressor,For the maximum allowable volume flow of entrance of the compressor,Obtained by the nameplate that dispatches from the factory of compressor, poutFor the outlet pressure of compressor, pc,maxFor the maximum allowable pressure in outlet of compressor
Power, pc,maxObtained by the nameplate that dispatches from the factory of compressor;
Gas storage equipment should be greater than being equal to 0 in the reserves of t in (3-8) plant gas, less than or equal to gas storage equipment most
Big volume:
0≤Vt≤Vmax;
(4) optimized algorithm is utilized, interior point method is used in one embodiment of the present of invention, it is mesh to solve by above-mentioned steps (1)
Scalar functions, above-mentioned steps (2) and step (3) are the Optimized model of constraints, obtain gas storage equipment under prediction year load
Month optimal gas storage capacity curve, the abscissa in optimal gas storage capacity curve is 12 months, vertical in optimal gas storage capacity curve
Coordinate is moon gas storage capacity;
(5) according to the optimal gas storage capacity curve of above-mentioned steps (4), the desired value of capacity is asked for, it is contemplated that gas storage equipment
Typically there is certain design specification (such as 1000m3,2000m3,5000m3Deng), using closest to the capacity of desired value as combustion gas
The optimal capacity of gas storage equipment in power plant.
Claims (1)
1. a kind of meter and the plant gas gas storage capacity calculation methods of electric-gas coupled system constraint, it is characterised in that this method bag
Include following steps:
(2) optimization object function for the plant gas cost of electricity-generating established in an electric-gas coupled system is as follows:
<mrow>
<mi>C</mi>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>y</mi>
<mi>e</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</munder>
<msub>
<mi>C</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>g</mi>
<mi>a</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>g</mi>
<mi>a</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>o</mi>
<mi>r</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>f</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
Ccoal(Gcoal)=Pcoal×Gcoal
Cgas(Lgas)=Pgas×Lgas
Wherein, C be Operation of Electric Systems totle drilling cost, Ccoal(Gcoal) it is coal-fired cost, GcoalFor coal consumption, PcoalFor coal
Price, Cgas(Lgas) it is combustion gas cost, LgasFor air consumption, PgasFor the purchasing price of natural gas, Cstorage(V) for gas storage into
This, V be gas storage equipment volume, TlifeFor the service life of gas storage equipment;
(2) the equality constraint equation of electric-gas coupled system steady state Safe Operation is established, including:
The power flow equation of power system is as follows in (2-1) electric-gas coupled system:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>P</mi>
<mi>i</mi>
</msup>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>cos</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>sin</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Q</mi>
<mi>i</mi>
</msup>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>sin</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>cos</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mi>j</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
Wherein, PiFor the injection active power of i-th of node in power system, QiFor the injection nothing of i-th of node in power system
Work(power, GijFor the corresponding conductance of the i-th row, jth row in the bus admittance matrix Y with power system, BijFor with power system
Bus admittance matrix Y in the i-th row, jth arrange corresponding susceptance, power system bus admittance matrix Y is from electric power system dispatching
Center obtains;
The waterpower equation of pipeline is as follows in natural gas grid in (2-2) electric-gas coupled system:
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>sgn</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>77.54</mn>
<msub>
<mi>T</mi>
<mi>n</mi>
</msub>
</mrow>
<msub>
<mi>p</mi>
<mi>n</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>D</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>2.5</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<mi>L</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
</msup>
<msub>
<mi>&gamma;</mi>
<mi>G</mi>
</msub>
<msubsup>
<mi>T</mi>
<mi>a</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
</msubsup>
<msub>
<mi>FZ</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msqrt>
<mrow>
<mo>(</mo>
<msubsup>
<mi>p</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</msqrt>
<mo>,</mo>
</mrow>
Wherein, fkmFor natural gas volume flow in the pipeline in natural gas grid between k-th of node and m-th of node, pk, pmPoint
Not Wei k-th of node and m-th of node pressure, DkmAnd LkmPipeline km between respectively k-th of node and m-th of node
Pipe diameter and duct length, F are the coefficient of friction of inner-walls of duct, and F is by formulaCalculate
Arrive, EfFor the efficiency factor of pipeline, EfValue is that 0.92, Re is Reynolds number, by formulaIt is calculated, ρ
For natural gas density, μmFor the coefficient of kinetic viscosity (relevant with natural gas density, can table look-up acquisition) of natural gas, γGFor natural gas
Proportion, 0 < γG< 1, TaFor natural gas mean temperature, TnAnd pnThe temperature and pressure of natural gas, Z respectively under standard stategFor
The average compressible coefficient of natural gas, 0.9 < Zg< 1.5, in above-mentioned natural gas grid in the waterpower equation of pipeline, when
When, the sgn in above formulap(pk,pm)=1, whenWhen, sgnp(pk,pm)=- 1;
The energy expenditure equation of compressor is as follows in the natural gas grid of (2-3) electric-gas coupled system:
<mrow>
<msup>
<mi>BHP</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>f</mi>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msubsup>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>m</mi>
</msub>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
Wherein, pk, pmThe pressure of k-th of node and m-th of node respectively in natural gas grid, BHPkmFor k-th of node and m
The energy expenditure of compressor between individual node,For the inlet volumetric flow of the compressor, ηcFor the gross efficiency of the compressor,
ckFor the polytropic coefficient of the compressor, ηcAnd ckObtained from the shop instructions of compressor;
Coupling side between the power system coupled in (2-4) electric-gas coupled system by gas turbine and natural gas grid
Journey is as follows:
μPTur=Hgas×[fTur+(Vt-1-Vt)],
Wherein, fTurThe volume flow of the combustion gas inputted for gas turbine by pipeline, PTurExported for the active power of gas turbine,
VtIt is gas storage equipment in the gas reserves of t, HgasFor the combustion heat value of natural gas, value 37.59MJ/m3, μ are combustion
The efficiency factor of gas-turbine, obtained by the shop instructions of gas turbine;
The node air balance equation of natural gas grid is as follows in (2-5) electric-gas coupled system:
AGF=L,
Wherein, AGFor node-branch road matrix of natural gas grid, f is the volume flow of natural gas grid branch road, and L is natural gas grid node
Gas load, L obtains according to natural gas grid history data;
(3) the inequality constraints condition of electric-gas coupled system steady state Safe Operation is set, including:
The power output P of (3-1) generators in power systems groupi genMore than or equal to 0, dispatched from the factory less than or equal to the generating set
The peak power provided on nameplateI.e.:
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>;</mo>
</mrow>
The voltage magnitude U of i-th of node of (3-2) power systemiIn the higher limit of the safe operation of power system voltage of settingWith
Lower limitU iBetween run,U iFor 0.95 times of i-th of node rated voltage,For 1.05 times of i-th of node rated voltage,
I.e.:
<mrow>
<msub>
<munder>
<mi>U</mi>
<mo>&OverBar;</mo>
</munder>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>;</mo>
</mrow>
The transmission capacity S of the l articles circuit in (3-3) power systemlSafe operation of power system less than or equal to setting transmits appearance
The maximum of amountI.e.:
<mrow>
<msub>
<mi>S</mi>
<mi>l</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&OverBar;</mo>
</mover>
<mi>l</mi>
</msub>
<mo>;</mo>
</mrow>
The pressure p of k-th of node in (3-4) natural gas gridkHigher limit, the lower limit of air pressure are run in the pipe safety of settingp k、It is interior, i.e.,:
<mrow>
<msub>
<munder>
<mi>p</mi>
<mo>&OverBar;</mo>
</munder>
<mi>k</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>;</mo>
</mrow>
The flow f of b-th of pipeline in (3-5) natural gas gridbIn higher limit, the lower limit of the pipe safety operating flux of settingf b、It is interior, i.e.,:
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>f</mi>
<mi>b</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mi>b</mi>
</msub>
<mo>;</mo>
</mrow>
Gas source feed amount f in (3-6) natural gas gridsLess than or equal to the maximum f that the source of the gas can provide natural gas flows,max, i.e.,:
fs≤fs,max;
The constraints of compressor safe operation in (3-7) natural gas grid:
<mrow>
<mn>1</mn>
<mo><</mo>
<mi>S</mi>
<mo><</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>=</mo>
<mfrac>
<msup>
<mi>p</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msup>
<msup>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>m</mi>
<mi>p</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>,</mo>
<msup>
<mi>p</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
Wherein:S is the step-up ratio of the compressor, SmaxIt is the maximum step-up ratio of the compressor, SmaxBy the nameplate that dispatches from the factory of compressor
Obtain,For the volume flow of the entrance of the compressor,For the maximum allowable volume flow of entrance of the compressor,
Obtained by the nameplate that dispatches from the factory of compressor, poutFor the outlet pressure of compressor, pc,maxFor the outlet maximum allowble pressure of compressor,
pc,maxObtained by the nameplate that dispatches from the factory of compressor;
Gas storage equipment should be greater than being equal to 0 in the reserves of t in (3-8) plant gas, and the maximum less than or equal to gas storage equipment is held
Product:
0≤Vt≤Vmax;
(4) optimized algorithm is utilized, it is that object function, above-mentioned steps (2) and step (3) are constraint bar to solve by above-mentioned steps (1)
The Optimized model of part, obtain the moon optimal gas storage capacity curve of the gas storage equipment under prediction year load, optimal gas storage capacity curve
In abscissa be 12 months, the ordinate in optimal gas storage capacity curve is moon gas storage capacity;
(5) according to the optimal gas storage capacity curve of above-mentioned steps (4), the desired value of capacity is asked for, by closest to the appearance of desired value
Measure the optimal capacity as gas storage equipment in plant gas.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710927821.3A CN107491849B (en) | 2017-10-09 | 2017-10-09 | Gas power plant gas storage capacity calculation method considering constraint of electricity-gas coupling system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710927821.3A CN107491849B (en) | 2017-10-09 | 2017-10-09 | Gas power plant gas storage capacity calculation method considering constraint of electricity-gas coupling system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107491849A true CN107491849A (en) | 2017-12-19 |
CN107491849B CN107491849B (en) | 2020-10-20 |
Family
ID=60653773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710927821.3A Active CN107491849B (en) | 2017-10-09 | 2017-10-09 | Gas power plant gas storage capacity calculation method considering constraint of electricity-gas coupling system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107491849B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109740242A (en) * | 2018-12-29 | 2019-05-10 | 重庆大学 | Consider that the electric-gas integrated energy system of natural gas thermal procession unifies energy flux computation method |
CN113570115A (en) * | 2021-07-01 | 2021-10-29 | 东方电气集团东方电机有限公司 | Comprehensive energy system P2G plant station planning method suitable for bidirectional energy flow |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958480A (en) * | 2016-05-27 | 2016-09-21 | 清华大学 | Combined static safety analysis method for electrical-gas coupling multiple energy flow system |
CN106096777A (en) * | 2016-06-12 | 2016-11-09 | 清华大学 | A kind of Optimization Scheduling of electrical couplings multipotency streaming system |
CN106096269A (en) * | 2016-06-12 | 2016-11-09 | 清华大学 | The Interval Power Flow computational methods of natural gas grid in a kind of electrical couplings system |
-
2017
- 2017-10-09 CN CN201710927821.3A patent/CN107491849B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958480A (en) * | 2016-05-27 | 2016-09-21 | 清华大学 | Combined static safety analysis method for electrical-gas coupling multiple energy flow system |
CN106096777A (en) * | 2016-06-12 | 2016-11-09 | 清华大学 | A kind of Optimization Scheduling of electrical couplings multipotency streaming system |
CN106096269A (en) * | 2016-06-12 | 2016-11-09 | 清华大学 | The Interval Power Flow computational methods of natural gas grid in a kind of electrical couplings system |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109740242A (en) * | 2018-12-29 | 2019-05-10 | 重庆大学 | Consider that the electric-gas integrated energy system of natural gas thermal procession unifies energy flux computation method |
CN113570115A (en) * | 2021-07-01 | 2021-10-29 | 东方电气集团东方电机有限公司 | Comprehensive energy system P2G plant station planning method suitable for bidirectional energy flow |
CN113570115B (en) * | 2021-07-01 | 2023-04-28 | 东方电气集团东方电机有限公司 | Comprehensive energy system P2G station planning method applicable to bidirectional energy flow |
Also Published As
Publication number | Publication date |
---|---|
CN107491849B (en) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106096269B (en) | The Interval Power Flow calculation method of natural gas grid in a kind of electric-gas coupled system | |
CN108667007A (en) | Meter and the voltage stability margin computational methods of electric-gas coupled system constraint | |
CN103490410B (en) | Micro-grid planning and capacity allocation method based on multi-objective optimization | |
CN104734155B (en) | A kind of method obtaining electrical interconnection energy resource system available transmission capacity | |
CN105005940A (en) | Correlation-considered GEPOPF calculation method | |
CN103034961B (en) | A kind of decision method of electrical network energy-saving horizontal | |
CN107451680A (en) | A kind of electrical interconnection system random optimum trend heuristic calculation methods | |
CN106096777A (en) | A kind of Optimization Scheduling of electrical couplings multipotency streaming system | |
CN102103720A (en) | Risk-based micro power grid distributed power generation standby optimized configuration method | |
Abdali et al. | Techno-Econo-Environmental study on the use of domestic-scale wind turbines in Iran | |
CN107491849A (en) | Meter and the plant gas gas storage capacity calculation methods of electrical couplings system restriction | |
CN107947245A (en) | Consider the equivalent optimal load flow model building method of natural gas system constraint | |
Liu et al. | Energy system optimization under uncertainties: A comprehensive review | |
Marneni et al. | Loss reduction and voltage profile improvement in a rural distribution feeder using solar photovoltaic generation and rural distribution feeder optimization using HOMER | |
Jirdehi et al. | A low‐carbon strategy using INVELOX turbines in the presence of real‐time energy price uncertainty | |
CN103761680B (en) | Grid and provincial dispatching method and system for AC/DC interconnected large power grid with wind farm | |
Ajami et al. | Optimal power flow for power system interconnection considering wind power plants intermittency | |
CN104732302A (en) | Multistage electromagnetic loop circuit optimized loop-opening method based on immunity algorithm | |
Qing et al. | Integrated optimal planning of distribution network with geographical-zone-restricted renewable energy sources | |
Attardo et al. | Urban energy hubs economic optimization and environmental comparison in Italy and Vietnam | |
CN105529747B (en) | A kind of modeling method rationalized uniform reallocation sub-area division wind-powered electricity generation and receive power | |
Silva et al. | Simulating an ocean wave power plant with Homer | |
Zhou et al. | Co-planning and feasibility assessment of an integrated energy system embedded with power-to-gas plants | |
Yu et al. | Stochastic optimal dispatch of virtual power plant considering correlation of distributed generations | |
Lukutin et al. | Stand-alone power supply system with DC photo-diesel source |
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 |