CN106548416A - A kind of wind energy turbine set and electricity turn the collaboration Site planning method of gas plant stand - Google Patents

Abstract

The invention discloses a kind of wind energy turbine set and electricity turn the collaboration Site planning method of gas plant stand, comprise the following steps：(1) set up the Mathematical Modeling that electricity turns gas plant stand；(2) with historical wind speed time series as original scene, cluster analysis is carried out to numerous historical datas using a kind of density clustering algorithm, the output of wind electric field scene after being reduced；(3) the output of wind electric field scene obtained according to the Mathematical Modeling and step (2) of step (1), to maximize net investment profit as optimization aim, sets up the collaboration siteselecting planning Mathematical Modeling based on scene analysis；(4) according to step (1) step (3), the collaboration siteselecting planning model that constructed wind energy turbine set and electricity turn gas plant stand is solved using AMPL/BONMIN solvers, obtains program results.The economy of programme is assessed and analyzed to the collaboration Site planning method of the present invention, it is possible to which effectively reduction system abandons air quantity it can be considered that electricity turns the technical characteristic of gas plant stand operation.

Description

A kind of wind energy turbine set and electricity turn the collaboration Site planning method of gas plant stand

Technical field

The present invention relates to Power System Planning field, more particularly to a kind of wind energy turbine set and electricity turn the collaboration addressing of gas plant stand Planing method.

Background technology

In recent years, the environmental pressure that energy-consuming is brought is gradually increased, energy development and the center of gravity for utilizing progressively variable wind The regenerative resources such as energy.But, in the trans-utilization of regenerative resource, wind abandoned, abandon the phenomenons such as light and yet suffer from.To lift energy Source is changed and utilization ratio, and the coupling between different energy sources system gradually receives publicity.Traditional power system passes through gas engine Group is associated with natural gas system, and the electricity for occurring in recent years turns gas (Power to gas, PtG) technology and causes therebetween Closed loop coupling is possibly realized.PtG technologies can convert electric energy to hydrogen or natural gas so that the wind energy that power system cannot be dissolved Can be conveyed and be stored in being injected into natural gas network after converted etc. the generated output of regenerative resource, so as to substantially add The space of big optimized utilizing energy.

In recent years, PtG technologies and its application in power system progressively receive more and more extensive concern. Connection bridge of the PtG plant stands as two energy resource systems of electric power and natural gas, its siteselecting planning directly affect PtG plant stands at two The effect and contribution of synthetic operation in energy resource system.But up to the present, research in terms of PtG plant stand investment plannings report compared with It is few.In addition, individually carry out the electric economy for turning gas plant stand investment construction often poor, this explanation PtG plant stand is when addressing is built It is necessary to consider the new model of collaboration investment planning.

The content of the invention

The purpose of the present invention is for the deficiencies in the prior art, there is provided a kind of wind energy turbine set and electricity turn the collaboration addressing of gas plant stand Planing method, solves the problems, such as that wind energy turbine set and electricity turn the investment planning of gas plant stand, it is ensured that electricity turns the returns of investment of gas plant stand.

In order to achieve the above object, the technical solution adopted in the present invention is as follows：A kind of wind energy turbine set and electricity turn gas plant stand Collaboration Site planning method, comprises the following steps：

(1) set up the Mathematical Modeling that electricity turns gas plant stand；

(2) with historical wind speed time series as original scene, numerous historical datas are entered using a kind of density clustering algorithm Row cluster analysis, the output of wind electric field scene after being reduced；

(3) the output of wind electric field scene obtained according to the Mathematical Modeling and step (2) of step (1), it is net to maximize investment Income is optimization aim, sets up the collaboration siteselecting planning Mathematical Modeling based on scene analysis；

(4) according to step (1)-step (3), gas is turned to constructed wind energy turbine set and electricity using AMPL/BONMIN solvers The collaboration siteselecting planning model of plant stand is solved, and obtains program results.

Further, the step (1) is specially：

Set up the Mathematical Modeling that electricity turns gas plant stand：

The electric energy of input is converted into PtG plant stands the natural gas of output, and its energy input output relation is:

In formula:WithRepresent that PtG plant stands turn gas discharge and the consumption of gas generation in moment t by electricity respectively Electrical power；h^ptgTurn the energy conversion efficiency of gas for PtG plant stand electricity；a_gasRepresent gas discharge of the electric energy conversion for identical energy Units conversion factor.

To maintain the normal work of PtG plant stands, there is its electrical power for being consumed technology to exert oneself bound, and need to meet phase The start-stop time answered limits, i.e.,:

In formula:WithIt is expressed as maintaining the minimum electrical power and maximum electricity required for PtG plant stand normal works Power；u_tFor the working condition of t PtG unit, 1 and 0 represents start respectively and shuts down；WithPtG units are represented respectively Minimum allow the available machine time and minimum allow downtime.

Further, the step (2) is specially：

It is assumed that day wind speed time series isN_wFor wind speed historical data Sample size, wind speed sample period interval be 1h, i.e. T=24.Every wind speed time series can be considered a sample point, two The distance between sample point g and l d_glCalculated using Euclidean distance, i.e.,:

After the distance between wind speed time series is calculated, the local density s of sample point g_gFor:

In formula:d_cutTo block distance；c(d_gl-d_cut) value relation be:

The distance between sample point g to the sample point l bigger than its local density is:

After cluster, if the element number in g-th class is N_g, then the scene probability of happening is q_g=N_g/N_w。

The power output of blower fan with the functional relation between wind speed can approximate description be:

In formula:P^windFor wind power output power；For blower fan rated power；v、v_in、v_outAnd v_rateRespectively wind speed, cut Enter wind speed, cut-out wind speed and rated wind speed.After cluster analysis obtains typical wind speed time series, using formula (9) Wind power output sequence is obtained, so as to the output of wind electric field scene after being reduced.

Further, the step (3) is specially：

Build the collaboration siteselecting planning model of wind energy turbine set and PtG plant stands：

When investment siteselecting planning is carried out to PtG plant stands and wind energy turbine set, it is contemplated that with regenerative resource (wind-force) genco For investment subject, the collaboration investment planning of the two is carried out；

(3.1) object function

Wind energy turbine set and PtG plant stands collaboration siteselecting planning object function be:

In formula：

In formula:N_profitRepresent the net profit of renewable energy power generation company； WithRepresent wind energy turbine set respectively runs income, cost of investment, operating cost, the discount factor of cost of investment, fixed operating cost With coefficient, unit capacity investment cost and unit quantity of electricity operating cost；WithPoint Not Biao Shi PtG plant stands operation income, cost of investment, operating cost, the discount factor of cost of investment, fixed operating cost system Number, unit capacity investment cost and unit quantity of electricity operating cost.Convenient to discuss, unit quantity of electricity operating cost is set to constant.W^S And W^wScene set and wind energy turbine set node set to be selected are represented respectively；q_sRepresent the probability of s-th scene generation；T_aRepresent one Period sum in it；Dt represents period interval；WithRepresent wind-powered electricity generation online price, natural gas Selling price, SO₂、NO_xAnd CO₂Discharge transaction value；Represent moment t PtG plant stands as the standby of system reserve capacity Use price； WithThe SO of per unit electricity is represented respectively₂、NO_xAnd CO₂Reduce discharging conversion factor；With Represent respectively the moment t wind-powered electricity generation online under s-th scene of i-th node to be selected exert oneself, PtG plant stands dissolve abandon wind power and PtG plant stands turn the gas discharge of gas generation by electricity；J is time conversion coefficient, and value is 1/365；x_iWith y_iIt is illustrated respectively in wind energy turbine set rated capacity, PtG plant stand rated capacities, Construction of Wind Power 0-1 that i-th node to be selected can build The 0-1 decision variables that decision variable and PtG plant stands are built.

(3.2) power system constraints

In formula:WithRepresent respectively s-th Under scene wind energy turbine set surf the Net out force vector, conventional power unit go out force vector, Gas Generator Set go out force vector, load vector, Line Flow Vector, conventional power unit upper and lower limit vector sum Gas Generator Set of exerting oneself is exerted oneself upper and lower limit vector, and B represents electric power networks admittance matrix Imaginary part；θ represents node phase angle vector；Represent line power transmission limit vector；WithRespectively i-th routine The unit and Gas Generator Set moment t under s-th scene exerts oneself,Respectively i-th routine The climbing of unit and Gas Generator Set and landslide speed.

(3.3) natural gas system constraints

Natural gas line, pressurizing point and balance network traffic constraints need to be mainly considered in natural gas network；

For natural gas line, if the first and last node of the pipeline be respectively m and n, flow through pipeline steady state flow and Its constraints can use formula (23) and formula (24) to represent:

In formula:WithIt is that moment t flows to the pipeline flow of node n, saves from node m under s-th scene The air pressure of point m and node n；C_mnFor channel factor；The respectively upper and lower limit of pipeline flow；Take It is worth and is:

For pressurizing point, if the air inlet node of pressurizing point k is e, gas outlet node is o, then pressurizing point is consumed Equivalent flow and its constraints are:

In formula:t_k,s,tWithThe equivalent flow of respectively s-th scene lower moment t pressurizing points consumption and flow through pressurizing point Flow；B and h^cThe energy conversion factor and whole work efficiency of the equivalent flow that compressor is consumed are represented respectively；A represents many Varying index；g_k,maxAnd g_k,minUpper and lower limit is compared in the respectively pressurization of pressurizing point k.

For whole natural gas network, which needs the flow equilibrium relation for meeting be expressed as:

In formula:A, U and T be respectively the node-branch incident matrix of natural gas network, node-pressurizing point incidence matrix and The incidence matrix of node-pressurizing point consumed flow；And τ_sRespectively s-th scene In lower natural gas network, the flow vector of pipeline and pressurizing point, gas source and flow amount vector, load flow vector, PtG plant stands are generated The flow vector that the gas discharge vector sum pressurizing point that gas discharge vector, Gas Generator Set are consumed is consumed；WithPoint The air pressure of the air inflow and m-th node of m-th source of the gas node of moment t not under s-th scene；WithRespectively The bound of m source of the gas node air inflow；WithThe upper and lower limit of respectively m-th node air pressure；a_i、b_iAnd c_iRespectively For the gas consumption coefficient of discharge of Gas Generator Set

(3.4) the collaboration construction of wind energy turbine set and PtG plant stands and operation constraints

For the collaboration construction and operation of wind energy turbine set and PtG plant stands, should also there is following constraints:

In formula:Wind-powered electricity generation prediction for s-th scene lower moment t is exerted oneself.

Beneficial effects of the present invention are as follows：(1) collaboration Site planning method of the invention can turn gas plant stand for electricity Technical characteristic Rational Model；(2) collaboration Site planning method of the invention can obtain wind energy turbine set and electricity turns the planning of gas plant stand As a result, the economy of programme is calculated and is analyzed, solves the problems, such as that electricity turns the returns of investment of gas plant stand；(3) collaboration of the invention Site planning method can effectively reduce system and abandon air quantity, and lift system wind power output is dissolved level.

Accompanying drawing content

Fig. 1 is the electric-gas interacted system that 39 node power systems of IEEE and Belgian 20 node natural gas system are constituted Figure；

Fig. 2 is system daily load curve figure；

Fig. 3 is that wind speed is exerted oneself box diagram under 9 and 13 each scene of node；

Fig. 4 is wind energy turbine set income of the node 9 under each scene and PtG plant stands income figure；

Fig. 5 is wind energy turbine set income of the node 13 under each scene and PtG plant stands income figure；

Fig. 6 is the proportionate relationship figure of unit capacity PtG plant stands income and wind energy turbine set income；

Fig. 7 be under each scene the construction of PtG plant stands to abandoning the impact figure of air quantity；

Fig. 8 is that one each scene of scheme sells tolerance comparison diagram from extraneous purchase tolerance, Gas Generator Set air consumption and PtG plant stands；

Fig. 9 is the air pressure change situation map of scheme one and scheme three Scenes, 1 time natural gas system node 5.

Specific embodiment

The present invention is further illustrated with reference to the accompanying drawings and examples.

The present invention provides a kind of wind energy turbine set and electricity turns the collaboration Site planning method of gas plant stand, comprises the following steps：

(1) set up the Mathematical Modeling that electricity turns gas plant stand；

(2) with historical wind speed time series as original scene, numerous historical datas are entered using a kind of density clustering algorithm Row cluster analysis, the output of wind electric field scene after being reduced；

(3) the output of wind electric field scene obtained according to the Mathematical Modeling and step (2) of step (1), it is net to maximize investment Income is optimization aim, sets up the collaboration siteselecting planning Mathematical Modeling based on scene analysis；

(4) according to step (1)-step (3), gas is turned to constructed wind energy turbine set and electricity using AMPL/BONMIN solvers The collaboration siteselecting planning model of plant stand is solved, and obtains program results.

Further, the step (1) is specially：

Set up the Mathematical Modeling that electricity turns gas plant stand：

The electric energy of input is converted into PtG plant stands the natural gas of output, and its energy input output relation is:

In formula:WithRepresent that PtG plant stands turn gas discharge and the consumption of gas generation in moment t by electricity respectively Electrical power；h^ptgTurn the energy conversion efficiency of gas for PtG plant stand electricity；a_gasRepresent gas discharge of the electric energy conversion for identical energy Units conversion factor.

To maintain the normal work of PtG plant stands, there is its electrical power for being consumed technology to exert oneself bound, and need to meet phase The start-stop time answered limits, i.e.,:

In formula:WithIt is expressed as maintaining the minimum electrical power and maximum electricity required for PtG plant stand normal works Power；u_tFor the working condition of t PtG unit, 1 and 0 represents start respectively and shuts down；WithPtG units are represented respectively Minimum allow the available machine time and minimum allow downtime.

Further, the step (2) is specially：

Wind power output typical scene is clustered.Wind speed time series is clustered using fast search density clustering algorithm, Wind power output typical scene is generated on the basis of wind speed typical scene.

It is assumed that day wind speed time series isN_wFor wind speed historical data Sample size, wind speed sample period interval be 1h, i.e. T=24.Every wind speed time series can be considered a sample point, two The distance between sample point g and l d_glCalculated using Euclidean distance, i.e.,:

After the distance between wind speed time series is calculated, the local density s of sample point g_gFor:

In formula:d_cutTo block distance；c(d_gl-d_cut) value relation be:

The distance between sample point g to the sample point l bigger than its local density is:

After cluster, if the element number in g-th class is N_g, then the scene probability of happening is q_g=N_g/N_w。

The power output of blower fan with the functional relation between wind speed can approximate description be:

In formula:P^windFor wind power output power；For blower fan rated power；v、v_in、v_outAnd v_rateRespectively wind speed, cut Enter wind speed, cut-out wind speed and rated wind speed.After cluster analysis obtains typical wind speed time series, using formula (9) Wind power output sequence is obtained, so as to the output of wind electric field scene after being reduced.

Further, the step (3) is specially：

Build the collaboration siteselecting planning model of wind energy turbine set and PtG plant stands：

When investment siteselecting planning is carried out to PtG plant stands and wind energy turbine set, it is contemplated that with regenerative resource (wind-force) genco For investment subject, the collaboration investment planning of the two is carried out；

(3.1) object function

The object function of wind energy turbine set and PtG plant stands collaboration siteselecting planning is its day net investment profit of maximization.Wind energy turbine set and The net profit of PtG plant stands deducts investment and operating cost for its operation income.Wherein, the operation income of wind energy turbine set is by selling electric income And energy-saving and emission-reduction income composition；PtG plant stands operation income by natural gas system sell gas income, energy-saving and emission-reduction income and Standby income composition.Wind energy turbine set and PtG plant stands collaboration siteselecting planning object function be:

In formula：

In formula:N_profitThe net profit of renewable energy power generation company is represented, is final output amount； WithRepresent respectively the operation income of wind energy turbine set, cost of investment, operating cost, The discount factor of cost of investment, fixed operating cost coefficient, unit capacity investment cost and unit quantity of electricity operating cost；WithRepresent respectively the operation income of PtG plant stands, cost of investment, operating cost, The discount factor of cost of investment, fixed operating cost coefficient, unit capacity investment cost and unit quantity of electricity operating cost.To discuss Convenient, unit quantity of electricity operating cost is set to constant.W^SAnd W^wScene set and wind energy turbine set node set to be selected are represented respectively；q_s Represent the probability of s-th scene generation；T_aRepresent intraday period sum；Dt represents period interval；WithRepresent wind-powered electricity generation online price, the selling price of natural gas, SO₂、NO_xAnd CO₂Discharge hand over Easy price；Represent moment t PtG plant stands as the standby price of system reserve capacity；WithRepresent respectively every The SO of unit quantity of electricity₂、NO_xAnd CO₂Reduce discharging conversion factor；WithRepresent i-th node to be selected at s-th respectively The online of scene lower moment t wind-powered electricity generations is exerted oneself, PtG plant stands the are dissolved natural gas for abandoning that the power and PtG plant stands of wind turn gas generation by electricity Flow；J is time conversion coefficient, and value is 1/365；x_iAnd y_iIt is illustrated respectively in i-th node to be selected to build Wind energy turbine set rated capacity, the 0-1 decision-makings built of PtG plant stand rated capacities, Construction of Wind Power 0-1 decision variables and PtG plant stands Variable.

(3.2) power system constraints

In formula:WithRepresent respectively s-th Under scene wind energy turbine set surf the Net out force vector, conventional power unit go out force vector, Gas Generator Set go out force vector, load vector, Line Flow Vector, conventional power unit upper and lower limit vector sum Gas Generator Set of exerting oneself is exerted oneself upper and lower limit vector, and B represents electric power networks admittance matrix Imaginary part；θ represents node phase angle vector；Represent line power transmission limit vector；WithRespectively i-th routine The unit and Gas Generator Set moment t under s-th scene exerts oneself,Respectively i-th routine The climbing of unit and Gas Generator Set and landslide speed.

(3.3) natural gas system constraints

Natural gas line, pressurizing point and balance network traffic constraints need to be mainly considered in natural gas network；

For natural gas line, if the first and last node of the pipeline be respectively m and n, flow through pipeline steady state flow and Its constraints can use formula (23) and formula (24) to represent:

In formula:WithIt is that moment t flows to the pipeline flow of node n, saves from node m under s-th scene The air pressure of point m and node n；C_mnFor channel factor；The respectively upper and lower limit of pipeline flow；Take It is worth and is:

For pressurizing point, if the air inlet node of pressurizing point k is e, gas outlet node is o, then pressurizing point is consumed Equivalent flow and its constraints are:

In formula:t_k,s,tWithThe equivalent flow of respectively s-th scene lower moment t pressurizing points consumption and flow through pressurizing point Flow；B and h^cThe energy conversion factor and whole work efficiency of the equivalent flow that compressor is consumed are represented respectively；A represents many Varying index；g_k,maxAnd g_k,minUpper and lower limit is compared in the respectively pressurization of pressurizing point k.

For whole natural gas network, which needs the flow equilibrium relation for meeting be expressed as:

In formula:A, U and T be respectively the node-branch incident matrix of natural gas network, node-pressurizing point incidence matrix and The incidence matrix of node-pressurizing point consumed flow；And τ_sRespectively s-th scene In lower natural gas network, the flow vector of pipeline and pressurizing point, gas source and flow amount vector, load flow vector, PtG plant stands are generated The flow vector that the gas discharge vector sum pressurizing point that gas discharge vector, Gas Generator Set are consumed is consumed；WithPoint The air pressure of the air inflow and m-th node of m-th source of the gas node of moment t not under s-th scene；WithRespectively m The bound of individual source of the gas node air inflow；WithThe upper and lower limit of respectively m-th node air pressure；a_i、b_iAnd c_iRespectively The gas consumption coefficient of discharge of Gas Generator Set

(3.4) the collaboration construction of wind energy turbine set and PtG plant stands and operation constraints

For the collaboration construction and operation of wind energy turbine set and PtG plant stands, should also there is following constraints:

In formula:Wind-powered electricity generation prediction for s-th scene lower moment t is exerted oneself.

Embodiment：

Parameter setting:Be made up of 39 node power systems of IEEE and Belgian 20 node natural gas system, Fig. 1 Shown electric-gas interacted system is example, and the method to being proposed tests.In 39 node systems of IEEE, positioned at node 30th, 35,36,37,38 and 39 unit is conventional power unit, is Gas Generator Set positioned at the unit of node 31 to 34.To 39 node The unit and load parameter of system is modified, as shown in table 1, table 2 and Fig. 2.Assume the node 9 and 13 in power system Wind energy resources relatively enriches, and power system and natural gas system are coupled in the two near nodals.Jing analyses can the test System has following two siteselecting planning schemes:

Scheme one:The wind energy turbine set of 400MW is built in the node 9 of power system；It is contemplated that the PtG factories of auxiliary construction 20MW Stand, be connected with the node 5 of natural gas system after building up.

Scheme two:The wind energy turbine set of 400MW is built in the node 13 of power system；It is contemplated that the PtG factories of auxiliary construction 20MW Stand, be connected with the node 10 of natural gas system after building up.

In order to probe into the impact of system Congestion, it is assumed that be possible to be transmitted electricity when at node 9 and 13, wind-powered electricity generation is sent The impact of obstruction, maximum delivery power are 300MW.The capacity that wind energy turbine set now can be built at such scheme interior joint 9 and 13 is constant, And the PtG plant stands capacity of auxiliary construction can expand as 50MW, both schemes are designated as into scheme three and scheme four respectively.

Using the wind speed historical data of American National Ministry of Energy regenerative resource laboratory (NREL) 2004-2006, if The incision wind speed of node 9 and 13 is 2m/s, and cut-out wind speed is 20m/s, and rated wind speed is respectively 12m/s and 10m/s, PtG The minimum operation power of plant stand is taken as the 5% of its rated power.Other price parameters, wind-powered electricity generation field parameters, PtG plant stand parameters etc. as Shown in table 3.

1 generator 's parameter of table

2 electric load parameter of table

3 wind energy turbine set of table, PtG plant stands and price parameter

* Annual Percentage Rate was set as 6%, by depreciation in 25 years.

The wind speed historical data of node 9 and 13 can respectively obtain 7 classes and 4 class wind speed seasonal effect in time series typical fields through cluster Scape.Fig. 3 a are that wind speed is exerted oneself box diagram under 9 each scene of node, and Fig. 3 b are that wind speed is exerted oneself box diagram under 13 each scene of node.Cluster The probability of happening of each scene marks the lower section of transverse axis correspondence scene in figure 3 respectively afterwards.As can be seen from Figure 3, to node 9 Speech, its typical wind power output scene are more, and the probability that each scene occurs is more or less the same.In the maximum scene 1 of probability of happening Wind power output level is higher and peak-valley difference is less.Compared with node 9, the wind power output scene of node 13 is less, and each scene is sent out Raw probability difference is very big；In the maximum scene 1 of probability of happening, wind power output in a few days changes acutely and fluctuating range is very big.

In the case of clog-free, the program results of scheme one and scheme two is as shown in table 4.Scheme one and scheme two are building wind After electric field, the PtG plant stands of auxiliary construction 20MW are considered.Scheme one day total income than more than 2 1,150,000 yuan of scheme, net profit ratio More than 2 1,100,000 yuan of scheme, scheme one are more preferable than the economy of scheme two, are considered as building wind energy turbine set auxiliary construction in node 9 Corresponding PtG plant stands.

4 each siteselecting planning scenario outcomes of table contrast (clog-free situation)

The wind energy turbine set income of two schemes and PtG plant stands are taken in as shown in Figure 4 and Figure 5.For ease of comparing, Fig. 4 and Fig. 5 are pressed Descending arrangement is carried out according to the height of wind energy turbine set income under each scene.Wind energy turbine set income is intraday power selling income and wind energy turbine set Income is reduced discharging, PtG plant stands income is taken in, reduces discharging income and standby income to sell gas.

For the composition of wind energy turbine set income and PtG plant stands income, it is observed that by the contrast of scheme one and scheme two Reduce discharging income and higher ratio is account in the two total income respectively, this is because SO₂、NO_xAnd CO₂The reduction of discharging price of three it About 0.8 yuan/kWh, the reduction of discharging price of per unit energy are higher than sale of electricity price and sell gas price lattice.The gas of selling of PtG plant stands is received Enter higher than its standby income, and this also with unit energy sell gas price lattice higher than PtG plant stands the standby price of unit energy it is relevant.

In the contrast of wind energy turbine set income and PtG plant stands income in scheme one and scheme two, wind energy turbine set income is long-range Take in PtG plant stands.For each node, PtG plant stands sell gas and PtG plant stands to be had as the income acquired in system reserve Limit, its major source of revenues remain wind energy turbine set income, and this is to sell gas price lattice and standby price less than single due to unit energy The sale of electricity price of potential energy amount, and the finite capacity of PtG plant stands causes the PtG plant stands income of each node to seem limited.Fig. 6 enters One step gives the proportionate relationship of each scheme unit capacity PtG plant stand income and wind energy turbine set income, and Fig. 6 a are scheme together scheme Three unit capacity income proportionate relationship comparison diagram, Fig. 6 b are unit capacity income proportionate relationship contrast of the scheme two with scheme four Figure, the proportionate relationship are defined as unit capacity PtG plant stand income and take in divided by unit capacity wind energy turbine set.It can be seen that, scheme one Income proportionate relationship with scheme two under some scenes is only 20%～30%, and unit capacity PtG plant stand income is not high.

Additionally, in the contrast of scheme one and scheme two, 1 time wind energy turbine set of Fig. 4 Scenes and PtG plant stands take in highest, this It is that, as under the scene, wind power output is relatively stable, amplitude is big, and the wind power output that power system and PtG plant stands are dissolved is more；And The wind energy turbine set income of scene 3 is greatly reduced on the contrary higher than scene 7, PtG plant stands income.This is due to load valley in scene 3 Overlap with wind power output low ebb, the obvious situation for crossing Sheng of wind-powered electricity generation generated output do not occur.Similar, 2 times wind energy turbine sets of Fig. 5 Scenes It is most with PtG plant stands income, and 4 times PtG plant stands of scene 3 and scene do not obtain income.It can be seen that, the income of PtG plant stands is not only received The impact of related output of wind electric field, also the fluctuation situation with load is closely related.

The impact of air quantity is abandoned in the construction that PtG plant stands are analyzed by taking scheme one as an example to system.Fig. 7 enters according to the size for abandoning air quantity Row sequence, gives the impact whether PtG plant stands construction in scheme one abandons air quantity to system.In scene 1, PtG plant stands are built After operation compared with the situation of PtG plant stands is not built, that can reduce 129MWh abandons air quantity；In scene 4, under abandoning wind ratio even more Drop 73.71%.Therefore, the construction of PtG plant stands can effectively reduce system with operation and abandon air quantity, improve efficiency of energy utilization.

Fig. 8 give in scheme one under each scene electric-gas interacted system from extraneous purchase tolerance, Gas Generator Set air consumption and PtG plant stands sell the schematic diagram of tolerance.Whole system is still relied primarily on from the external world and buys natural gas as can be seen from Figure 8, in PtG factories Station is sold in the more scene of tolerance, and system is slightly reduced from the natural gas that the external world buys.

In the case of system congestion, the siteselecting planning result of scheme three and scheme four is as shown in table 5.As can be seen from Table 5, The corresponding PtG plant stands of auxiliary construction are still considered in scheme three and scheme four.Contrast table 5 and table 4 it can be found that, three Tongfang of scheme Case one is compared net profit and have dropped 1.22%, and scheme four improves 0.74% than net profit with scheme two-phase.It can be seen that, system congestion And the lifting of PtG plant stand capacity generates certain impact to the income of programme.

5 each siteselecting planning scenario outcomes of table contrast (having congestion situations)

As can be seen from Figure 4 and Figure 5, Congestion is taken in and is not made significant difference to the wind energy turbine set under most of scene；And For PtG plant stands income, the appearance of Congestion causes system to abandon the increase of wind period.Air quantity increase is abandoned, the lifting of capacity makes The generated output of renewable energy source increase that PtG plant stands are dissolved is obtained, therefore the income of PtG plant stands is obviously improved.PtG factories in scheme four Increases of income of standing has been even more than the wind energy turbine set income reduction caused because of system congestion, and this also explains why scheme four Lifted than the net profit of scheme two.From Fig. 6 also it can be seen that when system occur block when, it is jumbo under most of scene PtG plant stands income proportionate relationship is improve on the contrary.When system occurs to block, unit capacity PtG plant stand is taken in significantly for this explanation Degree is improved.

Fig. 7 gives scheme one and the landscape condition of abandoning of scheme three is contrasted.Under different scenes, the air quantity of abandoning of system has Reduced.The air quantity of abandoning of some scenes significantly declines, and abandons air quantity and have dropped 47.53%, abandon under air quantity in scene 1 in scene 7 Drop 18.92%, part scene (scene 4 and scene 5) or even occurred in that zero abandons wind.This means when system occurs to block, The lifting of PtG capacity can be significantly improved abandons landscape condition, reduces the wind power output that each moment cannot dissolve because of system congestion.

Fig. 9 give scheme one and scheme three in scene 1 the air pressure change situation of natural gas system node 5 and with this The pipeline flow situation of change that node is connected.As can be seen from Figure 9,9 up to 14 when, the PtG plant stands in scheme three are in Peak power running status, therefore the air pressure and pipeline flow change of scheme three is more gentle.In addition, system congestion and PtG The lifting of plant stand capacity causes more synthetic natural gas injection natural gas networks, so that scheme three is compared with scheme one, Node air pressure amplitude rises, and pipeline flow increases.

Claims (4)

1. a kind of wind energy turbine set and electricity turn the collaboration Site planning method of gas plant stand, it is characterised in that comprise the following steps：

(1) set up the Mathematical Modeling that electricity turns gas plant stand；

(2) with historical wind speed time series as original scene, numerous historical datas are gathered using a kind of density clustering algorithm Alanysis, the output of wind electric field scene after being reduced；

(3) the output of wind electric field scene obtained according to the Mathematical Modeling and step (2) of step (1), to maximize net investment profit For optimization aim, the collaboration siteselecting planning Mathematical Modeling based on scene analysis is set up；

(4) according to step (1)-step (3), gas plant stand is turned to constructed wind energy turbine set and electricity using AMPL/BONMIN solvers Collaboration siteselecting planning model solved, obtain program results.

2. wind energy turbine set according to claim 1 and electricity turn the collaboration Site planning method of gas plant stand, it is characterised in that described Step (1) is specially：

Set up the Mathematical Modeling that electricity turns gas plant stand：

The electric energy of input is converted into PtG plant stands the natural gas of output, and its energy input output relation is:

$Q_t^{ptg}=h^{ptg}{a}_{gas}{P}_t^{ptg}---\left(1\right)$

In formula:WithRepresent that PtG plant stands turn the gas discharge and the electric work for consuming of gas generation in moment t by electricity respectively Rate；h^ptgTurn the energy conversion efficiency of gas for PtG plant stand electricity；a_gasRepresent list of the electric energy conversion for the gas discharge of identical energy Position conversion coefficient.

To maintain the normal work of PtG plant stands, its electrical power for being consumed that there is technology to exert oneself bound, and need to meet corresponding Start-stop time limits, i.e.,:

In formula:WithIt is expressed as maintaining the minimum electrical power and maximum electric power required for PtG plant stand normal works； u_tFor the working condition of t PtG unit, 1 and 0 represents start respectively and shuts down；WithRepresent PtG units most respectively Little permission available machine time and minimum permission downtime.

3. wind energy turbine set according to claim 1 and electricity turn the collaboration Site planning method of gas plant stand, it is characterised in that described Step (2) is specially：

It is assumed that day wind speed time series isG=1,2 ..., N_w。N_wFor the sample number of wind speed historical data Amount, wind speed sample period interval are 1h, i.e. T=24.Every wind speed time series can be considered a sample point, two sample point g With the distance between l d_glCalculated using Euclidean distance, i.e.,:

$d_{gl}=\sqrt{\underset{t=1}{\overset{T}{\Σ}}{(v_{g,t}-v_{l,t})}^2}---\left(5\right)$

After the distance between wind speed time series is calculated, the local density s of sample point g_gFor:

$s_g=\underset{l,,g}{\Σ}c(d_{gl}-d_{cut})---\left(6\right)$

In formula:d_cutTo block distance；c(d_gl-d_cut) value relation be:

The distance between sample point g to the sample point l bigger than its local density is:

After cluster, if the element number in g-th class is N_g, then the scene probability of happening is q_g=N_g/N_w。

The power output of blower fan with the functional relation between wind speed can approximate description be:

In formula:P^windFor wind power output power；For blower fan rated power；v、v_in、v_outAnd v_rateRespectively wind speed, incision wind Speed, cut-out wind speed and rated wind speed.After cluster analysis obtains typical wind speed time series, it is obtained using formula (9) Wind power output sequence, so as to the output of wind electric field scene after being reduced.

4. wind energy turbine set according to claim 1 and electricity turn the collaboration Site planning method of gas plant stand, it is characterised in that described Step (3) is specially：

Build the collaboration siteselecting planning model of wind energy turbine set and PtG plant stands：

When investment siteselecting planning is carried out to PtG plant stands and wind energy turbine set, it is contemplated that with regenerative resource (wind-force) genco as throwing Money main body, carries out the collaboration investment planning of the two；

(3.1) object function

Wind energy turbine set and PtG plant stands collaboration siteselecting planning object function be:

$\max N_{profit}=R_{rev}^{wind}+R_{rev}^{ptg}-(C_{inv}^{wind}+C_{inv}^{ptg})-(C_{opr}^{wind}+C_{opr}^{ptg})---\left(10\right)$

In formula：

$R_{rev}^{wind}=\underset{i.W^w}{\Σ}\underset{s.W^s}{\Σ}{q}_s\underset{t=1}{\overset{T_a}{\Σ}}\{p_{price}^{wind}{P}_{i,s,t}^{wind}Dt+(p_{price}^{{\mathrm{SO}}_2}{a}_{{\mathrm{SO}}_2}+p_{price}^{{\mathrm{NO}}_x}{a}_{{\mathrm{NO}}_x}+p_{price}^{{\mathrm{CO}}_2}{a}_{{\mathrm{CO}}_2})P_{i,s,t}^{wind}Dt\}---\left(11\right)$

$R_{rev}^{ptg}=\underset{i.W^w}{\Σ}\underset{s.W^s}{\Σ}{q}_s\underset{t=1}{\overset{T_a}{\Σ}}\{p_{price}^{gas}{Q}_{i,s,t}^{ptg}Dt+(p_{price}^{{\mathrm{SO}}_2}{a}_{{\mathrm{SO}}_2}+p_{price}^{{\mathrm{NO}}_x}{a}_{{\mathrm{NO}}_x}+p_{price}^{{\mathrm{CO}}_2}{a}_{{\mathrm{CO}}_2})P_{i,s,t}^{ptg}Dt+p_t^{spr}{P}_{i,s,t}^{ptg}\}---\left(12\right)$

$C_{inv}^{wind}=j(e_{dis}^{wind}+e_{fix}^{wind})\underset{i.W^w}{\Σ}{x}_i{c}_{inv}^{wind}{P}_{i,rate}^{wind}---\left(13\right)$

$C_{inv}^{ptg}=(e_{dis}^{ptg}+e_{fix}^{ptg})\underset{i.W^w}{\Σ}{y}_i{c}_{inv}^{ptg}{P}_{i,rate}^{ptg}---\left(14\right)$

$C_{opr}^{wind}=\underset{i.W^w}{\Σ}\underset{s.W^s}{\Σ}{q}_s\underset{t=1}{\overset{T_a}{\Σ}}{c}_{opr}^{wind}{P}_{i,s,t}^{wind}Dt---\left(15\right)$

$C_{opr}^{ptg}=\underset{i.W^w}{\Σ}\underset{s.W^s}{\Σ}{q}_s\underset{t=1}{\overset{T_a}{\Σ}}{c}_{opr}^{ptg}{P}_{i,s,t}^{ptg}Dt---\left(16\right)$

In formula:N_profitRepresent the net profit of renewable energy power generation company；WithRepresent that run income, cost of investment, operating cost, the discount factor of cost of investment, the fixed operating cost of wind energy turbine set are used respectively Coefficient, unit capacity investment cost and unit quantity of electricity operating cost；WithRespectively Represent PtG plant stands operation income, cost of investment, operating cost, the discount factor of cost of investment, fixed operating cost coefficient, Unit capacity investment cost and unit quantity of electricity operating cost.Convenient to discuss, unit quantity of electricity operating cost is set to constant.W^SAnd W^w Scene set and wind energy turbine set node set to be selected are represented respectively；q_sRepresent the probability of s-th scene generation；T_aRepresent one day in Period sum；Dt represents period interval；WithRepresent wind-powered electricity generation online price, the pin of natural gas Price lattice, SO₂、NO_xAnd CO₂Discharge transaction value；Represent moment t PtG plant stands as the standby valency of system reserve capacity Lattice； WithThe SO of per unit electricity is represented respectively₂、NO_xAnd CO₂Reduce discharging conversion factor；WithRespectively Represent that exert oneself under s-th scene, PtG plant stands are dissolved and abandon power and the PtG factories of wind by the online of moment t wind-powered electricity generation for i-th node to be selected Stand and turn the gas discharge of gas generation by electricity；J is time conversion coefficient, and value is 1/365；x_iAnd y_iRespectively Represent that the wind energy turbine set rated capacity that can be built in i-th node to be selected, PtG plant stand rated capacities, Construction of Wind Power 0-1 decision-makings become The 0-1 decision variables that amount and PtG plant stands are built.

(3.2) power system constraints

$P_s^{wind}+P_s^{fossil}+P_s^{gas}-P_s^{load}=B\mathrm{\θ}---\left(17\right)$

In formula:WithS-th scene is represented respectively Lower wind energy turbine set surf the Net out force vector, conventional power unit go out force vector, Gas Generator Set go out force vector, load vector, Line Flow vector, Conventional power unit upper and lower limit vector sum Gas Generator Set of exerting oneself is exerted oneself upper and lower limit vector, and B represents the imaginary part of electric power networks admittance matrix； θ represents node phase angle vector；Represent line power transmission limit vector；WithRespectively i-th conventional power unit and Gas Generator Set moment t under s-th scene exerts oneself,Respectively i-th conventional power unit and The climbing of Gas Generator Set and landslide speed.

(3.3) natural gas system constraints

Natural gas line, pressurizing point and balance network traffic constraints need to be mainly considered in natural gas network；

For natural gas line, if the first and last node of the pipeline is respectively m and n, the steady state flow and its about of pipeline is flowed through Beam condition can use formula (23) and formula (24) to represent:

$f_{m,n,s,t}^{gas}=S(r_{m,s,t}^{gas},r_{n,s,t}^{gas})C_{mn}\sqrt{S(r_{m,s,t}^{gas},r_{n,s,t}^{gas})|{\left(r_{m,s,t}^{gas}\right)}^2-{\left(r_{n,s,t}^{gas}\right)}^2)}---\left(23\right)$

In formula:WithBe under s-th scene moment t from node m flow to the pipeline flow of node n, node m and The air pressure of node n；C_mnFor channel factor；The respectively upper and lower limit of pipeline flow；Value is:

For pressurizing point, if the air inlet node of pressurizing point k is e, gas outlet node is o, then pressurizing point consumption is equivalent Flow and its constraints are:

$t_{k,s,t}=b\frac{f_{k,s,t}^{gas}a}{h^c(a-1)}\[{\left(\frac{r_{o,s,t}^{gas}}{r_{e,s,t}^{gas}}\right)}^{(a-1)/a}-1\]---\left(26\right)$

In formula:t_k,s,tWithThe equivalent flow of respectively s-th scene lower moment t pressurizing points consumption and the stream for flowing through pressurizing point Amount；B and h^cThe energy conversion factor and whole work efficiency of the equivalent flow that compressor is consumed are represented respectively；A represents changeable finger Number；g_k,maxAnd g_k,minUpper and lower limit is compared in the respectively pressurization of pressurizing point k.

For whole natural gas network, which needs the flow equilibrium relation for meeting be expressed as:

In formula:A, U and T respectively node-branch incident matrix of natural gas network, node-pressurizing point incidence matrix and node- The incidence matrix of pressurizing point consumed flow；And τ_sIt is natural under respectively s-th scene The natural gas that the flow vector of pipeline and pressurizing point, gas source and flow amount vector, load flow vector, PtG plant stands are generated in gas network The flow vector that the gas discharge vector sum pressurizing point that flow vector, Gas Generator Set are consumed is consumed；WithRespectively in s The air pressure of the air inflow and m-th node of individual m-th source of the gas node of scene lower moment t；WithRespectively m-th source of the gas The bound of node air inflow；WithThe upper and lower limit of respectively m-th node air pressure；a_i、b_iAnd c_iRespectively gas engine The gas consumption coefficient of discharge of group

(3.4) the collaboration construction of wind energy turbine set and PtG plant stands and operation constraints

For the collaboration construction and operation of wind energy turbine set and PtG plant stands, should also there is following constraints:

In formula:Wind-powered electricity generation prediction for s-th scene lower moment t is exerted oneself.