CN111401647B - Distributed optimal scheduling method for electric coupling system considering uncertainty transfer - Google Patents

Abstract

The invention relates to an electric coupling system distributed optimization scheduling method considering uncertainty transfer, and belongs to the technical field of operation control of comprehensive energy systems. The method fully considers the transmission of uncertainty in the power grid and the natural gas network in the optimization process of the electric coupling system, and establishes a distributed optimization model of the power grid and the natural gas network, so that more reasonable parameters for optimizing the operation of the electric coupling system are obtained. In the method, a power grid constraint condition and a natural gas grid constraint condition which take the uncertainty of the wind power active power into consideration are established; establishing a power grid and natural gas grid distributed optimization model; the method for optimizing information interaction of the power grid and the natural gas grid is provided, and by continuously interacting optimization information, a plurality of uncertainties caused by injection of high-proportion renewable energy sources in the electric coupling system are overcome, so that distributed optimization of the electric coupling system is finally realized, and safe, reliable and economic operation of the electric coupling system is realized.

Description

Distributed optimal scheduling method for electric coupling system considering uncertainty transfer

Technical Field

The invention relates to an electric coupling system distributed optimization scheduling method considering uncertainty transfer, and belongs to the technical field of operation control of comprehensive energy systems.

Technical Field

With the development of renewable energy sources and distributed power generation technologies, the existing power grid is more and more difficult to meet the requirements of people on high efficiency and greenization of energy sources. In the traditional energy system, different energy industries such as power supply, heat supply, cold supply, gas supply and the like are relatively closed, the interconnection degree is limited, and the improvement of energy efficiency and the consumption of renewable energy are not facilitated. Therefore, how to realize the comprehensive utilization of multiple energy sources such as electricity, heat, cold, gas, oil, traffic and the like to form an open interconnected multifunctional coupling system taking electricity as a core has become a new focus of attention in the international academic world and the industrial world at present.

The natural gas network becomes a main research object of the multi-energy coupling system due to the characteristics of wide distribution, huge volume, considerable optimization space, high coupling degree with the power grid and the like. The natural gas grid is coupled to the power grid primarily through a gas power plant. The gas power station has the advantages of low investment cost, high energy utilization efficiency, high flexibility, low price and the like, and the installed capacity of the gas power station is rapidly increased in the world. As consumers in a natural gas network and producers in a power grid, a gas power station creates possibility for coordinated operation of an electrical coupling system and improvement of overall benefits, but also brings new risks, such as sufficient natural gas supply, fluctuation of natural gas market price, pipeline accidents and the like, which directly affect the safety and economy of power grid operation, and changes in power load requirements also cause changes in natural gas flow in the natural gas network. Therefore, how to realize safe, reliable and economic operation of the electrical coupling system becomes a hot point of research. In addition, high proportion of renewable energy injection is the development trend of energy systems, and future electrical coupling systems will contain a lot of uncertainty. At present, the optimization research on the electrical coupling system does not consider the influence of uncertainty transmitted in a power grid and a natural gas grid.

Disclosure of Invention

The invention aims to provide an electric coupling system distributed optimization scheduling method considering uncertainty transfer, which improves the existing electric coupling system scheduling method, fully considers the uncertainty transfer in a power grid and a natural gas grid in the electric coupling system optimization process, and establishes a power grid and natural gas grid distributed optimization model so as to obtain more reasonable parameters for optimizing the operation of the electric coupling system.

The invention provides an electric coupling system distributed optimization scheduling method considering uncertainty transfer, which comprises the following steps:

(1) the electric coupling system is divided into a power grid and a natural gas grid, and the power grid and the natural gas grid are coupled through H gas power stations;

(2) recording the value of the gas consumption of the gas power station h in the power grid

The value of the gas consumption of the gas power station h in the natural gas network is recorded as

And

the following relation is satisfied:

(3) establishing a power grid constraint condition, comprising the following steps:

(3-1) establishing the electric quantity balance constraint of the power grid node as follows:

in the formula, M is the node serial number in the power grid, M is the total number of nodes in the power grid,

the active power of the gas power station h, the variable to be solved,

the sum of the active power of all the gas power stations connected with the node m in the power grid;

is the active power of the non-gas power station i, is a variable to be solved,

the sum of the active power of all the non-gas power stations connected with the node m in the power grid;

is a predicted value of the active power of the wind turbine generator j, is a known quantity and is given by the dispatching of a power grid, WS_jThe active power is used for abandoning the wind turbine generator j, the variable to be solved is obtained,

the sum of active power actually injected into the power grid by all wind turbine generators connected with the node m in the power grid is obtained; PD (photo diode)_kFor the predicted value of the active power of the electrical load k, given by the grid schedule, LS, for a known quantity_kThe power of the electric load k is the power-abandoning active power, which is a variable to be solved,

the sum of the actual active power of all the electric loads connected with the node m in the power grid; pf_lmThe active power of the branch between the node l and the node m is defined as a variable to be solved, the flow from the node l to the node m is positive, the flow from the node m to the node l is negative,

the sum of the active power flowing to the node m from all the nodes l connected with the node m in the power grid;

(3-2) establishing a power grid direct current power flow constraint as follows:

in the formula, theta_lAnd theta_mVoltage phase angles of a node l and a node m are respectively used as variables to be solved; x is the number of_lmRepresenting the reactance of a branch connected between the node l and the node m, wherein the reactance is a known quantity and is given by power grid dispatching;

(3-3) establishing a voltage phase angle constraint of a reference node in the power grid as follows:

θ_n＝0,n∈REF^p

in the formula, theta_nRepresenting the phase angle, REF, of node n in the grid^pIn the representation of the electric networkThe reference node set is given by power grid dispatching;

(3-4) establishing the upper limit constraint and the lower limit constraint of the active power of the gas power station in the power grid as follows:

in the formula (I), the compound is shown in the specification,

and

respectively setting the upper limit of active power and the lower limit of active power of the gas power station h by power grid dispatching;

(3-5) establishing the upper limit constraint and the lower limit constraint of the active power of the non-gas power station in the power grid as follows:

in the formula (I), the compound is shown in the specification,

and

respectively setting an upper active power limit and a lower active power limit of the non-gas power station i by power grid dispatching;

(3-6) establishing the upper limit constraint and the lower limit constraint of the abandoned power of the wind turbine generator as follows:

in the formula (I), the compound is shown in the specification,

the upper limit of the active power for abandoning the wind turbine generator j is given by the dispatching of a power grid;

(3-7) establishing the upper limit constraint and the lower limit constraint of the electric load electricity abandoning active power as follows:

0≤LS_k≤LS_k ^max

in the formula, LS_k ^maxAbandoning an upper limit of active power for the electric load k, and dispatching and giving the upper limit by a power grid;

(3-8) establishing upper limit constraint and lower limit constraint of branch active power in the power grid as follows:

-pf_lm ^max≤pf_lm≤pf_lm ^max

in the formula, pf_lm ^maxThe upper limit of the active power of the branch between the node l and the node m is given by the dispatching of a power grid;

(3-9) establishing the active power of the gas power plant h

And gas consumption

The constraints between are as follows:

in the formula, a_h ^ngu、b_h ^nguAnd c_h ^nguQuadratic term coefficients, primary term coefficients and constant terms which are quadratic relational expressions of the active power and the gas consumption of the gas power station h are respectively given by the gas power station;

(3-10) setting an active power variation interval of the wind turbine generator as follows:

wherein the content of the first and second substances,

the upper limit of the active power fluctuation of the wind turbine generator j is a known quantity and is given by the dispatching of the power grid,

the lower limit of the active power fluctuation of the wind turbine generator j is a known quantity and is given by the dispatching of the power grid,

the active power of the wind turbine generator j is obtained;

(3-11) setting quasi-steady-state output power transfer distribution factors of all power stations, namely gas power stations, non-gas power stations and wind power generation sets in power grid

In the formula (I), the compound is shown in the specification,

the vector is an Nx 1-dimensional column vector, N is the total number of power stations in a power grid, including a gas power station, a non-gas power station and a wind power generation unit, and the upper standard R is a quasi-steady-state identifier; h_lmFor each station n to the active power pf of the branch between node l and node m_lmN x 1-dimensional column vector, I, formed by the transfer distribution factors of_NIs an NxN dimensional identity matrix, alpha_NAn N x 1-dimensional column vector composed of the bearing coefficients bearing the unbalanced power for each station N,

is an N × 1 dimensional column vector with all values of 1;

above H_lmThe elements in (a) are represented as:

in the formula, n^pNumbering nodes of a power station n in a power grid, X being an impedance matrix of the power grid,

is the l-th row and the n-th row in the impedance matrix X^pThe elements of the column are,

is the m-th row and the n-th row in the impedance matrix X^pThe elements of the column are,

and

are all given by the power grid dispatching;

coefficient of bearing alpha_NIn the wind turbine, the bearing coefficient is 0, alpha_NThe bearing coefficient of the medium gas power station and the non-gas power station is more than 0, alpha_NGiven by the grid schedule and satisfying the following relationship:

(3-12) setting the active power adjustment vector of each power station in the power grid caused by the active power change of the wind turbine generator as

Is an N x 1-dimensional column vector,

the element values of the corresponding wind generating set are as follows:

j^psfor the numbering of the wind turbine j in each station in the grid,

taking the element values of corresponding gas power stations and non-gas power stations as 0;

(3-13) advantageCalculating the active power pf of the branch between the node l and the node m by using the following formula_lmAmount of change of

Is expressed as [ Delta pf ]_lm ^min,Δpf_lm ^max]，Δpf_lm ^minIs composed of

Lower limit value of value range, Δ pf_lm ^maxIs composed of

Taking the upper limit value of the value interval;

(3-14) establishing the active power upper limit and active power lower limit constraints of the branch in the power grid considering the active power change of the wind turbine generator as follows:

-pf_lm ^max≤pf_lm+Δpf_lm ^min≤pf_lm ^max

-pf_lm ^max≤pf_lm+Δpf_lm ^max≤pf_lm ^max

(3-15) calculating the active power regulating quantity of each power station in the power grid

In the formula (I), the compound is shown in the specification,

is an N x 1-dimensional column vector,

any of the elements of

Has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

Taking the upper limit value of the value interval;

(3-16) establishing the active power upper limit and the active power lower limit of a gas power station in a power grid considering the active power change of the wind turbine generator as follows:

in the formula, h^psThe serial numbers of the gas power stations h in each power station in the power grid,

is composed of

H in (1)^psThe lower limit value of the value interval of each element,

is composed of

H in (1)^psThe upper limit value of the value interval of each element;

(3-17) establishing the active power upper limit and the active power lower limit of a non-gas power station in the power grid considering the active power change of the wind turbine generator as follows:

in the formula i^psThe serial numbers of the non-gas power stations i in each power station in the power grid,

is composed of

I of (1)^psThe lower limit value of the value interval of each element,

is composed of

I of (1)^psThe upper limit value of the value interval of each element;

(3-18) establishing the following constraint between the active power variation and the air consumption variation of the gas power station h considering the active power variation of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

is the variation of the gas consumption caused by the variation of the active power of the gas power station h,

has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

The upper limit value of the value-taking interval,

and

the calculation formula is as follows:

(4) establishing natural gas network constraint conditions, comprising the following steps:

(4-1) establishing the natural gas flow balance constraint of the natural gas network node as follows:

in the formula, R is the node number in the natural gas network, and R is the node number in the natural gas network; s is the number of the gas well in the natural gas network, G_sThe outlet gas flow of the natural gas well s is used as a variable to be solved,

the sum of the gas outlet flow of all the natural gas wells connected with the node r in the natural gas network; t is the number of the gas load of residents in the natural gas network,

the gas consumption for the resident gas load t, which is a known quantity, is given by the natural gas network dispatch,

the gas consumption of the gas load of all residents connected with the node r in the natural gas network;

the sum of the gas consumption of all gas power stations connected with the node r in the natural gas network is a variable to be solved; gf_urFor the natural gas flow of the pipeline between the node u and the node r in the natural gas network, the gf when the natural gas flows from the node u to the gas point r is specified as a variable to be solved_urTake positive value, gf when flowing from node r to node u_urTaking the negative value of the reaction mixture,

the natural gas flow of all nodes connected with the node r in the natural gas network flows into the node r;

(4-2) establishing upper and lower constraints on node pressures in the natural gas network as follows:

in the formula (I), the compound is shown in the specification,

and

respectively representing the lower pressure limit and the upper pressure limit of a node r in the natural gas network, and being given by the scheduling of the natural gas network;

(4-3) establishing the relationship between the natural gas flow and pressure in the natural gas grid as follows:

in the formula, ω_uAnd ω_rPressure at node u and node r, sgn (ω), respectively, in the natural gas network_u,ω_r) Is about ω_u、ω_rFunction of when ω is_u＞ω_rThen, sgn (ω)_u,ω_r) Take 1 when ω_u≤ω_r，sgn(ω_u,ω_r) The value is 0; c_urIs the Welmos constant of the pipeline between node u and node r, a known quantity, given by the natural gas grid schedule, due to sgn (ω) in the constraint on the relationship between natural gas flow and pressure_u,ω_r) Is a binary variable, and an integer variable is introduced

The following relation is satisfied:

defining a mathematical variable F_ur，

And translating the constraint on the relationship between natural gas flow and pressure into the following expression:

in the formula (I), the compound is shown in the specification,

and

respectively setting a lower pressure limit and an upper pressure limit of the node u by scheduling of a natural gas network;

wherein the above-mentioned

The further relaxation is:

(4-4) establishing a pressure reference node constraint in the natural gas network as follows:

ω_v＝PR,v∈REF^g

in the formula, ω_vRepresenting the pressure at node v, PR is a constant given by the natural gas grid schedule, REF^gA set of reference nodes representing a natural gas network, given by a natural gas network schedule;

(4-5) establishing the upper limit constraint and the lower limit constraint of the gas well effluent flow in the natural gas network as follows:

wherein S is the number of natural gas wells,

the lower limit and the upper limit of the gas flow rate of the natural gas well are respectively given by the dispatching of a natural gas network;

(4-6) establishing a natural gas flow constraint for the pipes in the natural gas network as follows:

in the formula (I), the compound is shown in the specification,

the upper limit of the natural gas flow of the pipeline between the node u and the node r in the natural gas network is given by the scheduling of the natural gas network;

(4-7) defining nodes connected with the natural gas well and the gas power station in the natural gas network as gas injection quantity variable nodes, marking as w, and marking the number of all the gas injection quantity variable nodes in the natural gas network as Q;

(4-8) setting a transfer distribution factor of the quasi-steady-state gas outlet flow of each variable gas injection amount node in the natural gas network

In the formula (I), the compound is shown in the specification,

is a Q x 1-dimensional column vector, K_urIs a Q x 1-dimensional column vector, K_urThe natural gas flow gf of the pipeline between the node u and the node r in the natural gas network is jointed by each variable gas injection amount node_urComposition of transfer distribution factor of (I)_QIs a QxQ dimensional identity matrix, alpha_QIs a Q x 1-dimensional column vector, alpha_QThe device consists of bearing coefficients of unbalanced natural gas flow borne by each variable gas injection quantity node,

a Q × 1 dimensional column vector with all values being 1;

wherein K_urThe elements in (a) are represented as:

where k is the number of the pipe between node u and node r, T is the road-pipe correlation matrix, T is a known quantity, given by the natural gas network schedule,

for the qth in the road-pipe association matrix T^gRow, k column element, q^gNumbering nodes of the variable gas injection quantity node q in the natural gas network;

coefficient of bearing alpha_QThe value rule of the elements is as follows: bearing coefficient alpha of natural gas well_QGreater than 0, the gas power plant has a bearing coefficient alpha_QIs equal to 0, alpha_QGiven by the natural gas grid schedule and satisfying the following relation:

(4-9) recording the Q multiplied by 1 dimension row vector formed by the variable gas consumption interval of each gas injection variable node as

Wherein the row value of the gas injection quantity variable node corresponding to the gas power station h is

The corresponding row value of other natural gas wells is 0;

(4-10) UsingCalculating the natural gas flow gf of the pipeline between the node u and the node r in the natural gas network_urAmount of change of

In the formula (I), the compound is shown in the specification,

has a value interval of [ delta gf ]_ur ^min,Δgf_ur ^max]，Δgf_ur ^minIs composed of

Lower limit value of interval, delta gf_ur ^maxIs composed of

Taking the upper limit value of the value interval;

(4-11) establishing the constraint of the natural gas flow rate of the pipeline in the natural gas network considering the active power change of the wind turbine generator as follows:

(4-12) calculating the adjustment amount of the natural gas flow at each variable node of the gas injection amount by using the following formula

In the formula (I), the compound is shown in the specification,

is a Q x 1 dimensional column vector,

any one element of

Has a value interval of [ Delta G ]_q ^min,ΔG_q ^max]，ΔG_q ^minIs composed of

Lower limit of value range, Δ G_q ^maxIs composed of

Taking the upper limit value of the value interval;

(4-13) establishing the upper limit constraint and the lower limit constraint of the gas flow of the natural gas well in the natural gas network considering the active power change of the wind turbine generator as follows:

in the formula, s^gsNumbering all gas injection quantity variable nodes of the natural gas wells in a natural gas network;

is composed of

Zhongth s^gsThe lower limit value of the value interval of each element,

is composed of

S of (1)^gsThe upper limit value of the value interval of each element;

(4-14) establishing the change quantity of the gas consumption quantity of the gas power station in consideration of the active power change of the wind turbine generator

The safety constraint of the node pressure in the natural gas network when the boundary value is taken comprises the following steps:

(4-14-1) setting the gas consumption of the gas power station to rise

In the time, the Q multiplied by 1 dimension row vector formed by the variable gas consumption interval of each gas injection variable node is recorded as delta G⁽⁰),^upWherein the row value of the gas injection variable node corresponding to the gas power station h

The corresponding row value of the natural gas well is 0;

(4-14-2) calculating the natural gas flow rate gf of the pipeline between the node u and the node r in the natural gas network by using the following formula_urAmount of change of

(4-14-3) defining mathematical variables

The method comprises the following steps of establishing the following safety constraints of the pressure of the nodes in the natural gas network considering the active power change of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

and

gas consumption rise of node u and node r in gas power station

The pressure of time;

(4-14-4) setting the gas consumption variation of the gas power station

In the time, the Q multiplied by 1 dimension column vector formed by the variable gas consumption interval of each gas injection variable node is delta G⁽⁰),^downWherein the row value of the gas injection variable node corresponding to the gas power station h

The corresponding row value of the natural gas well is 0;

(4-14-5) calculating the natural gas flow gf of the pipeline between the node u and the node r in the natural gas network_urAmount of change of

(4-14-6) defining mathematical variables

The method comprises the following steps of establishing the following safety constraints of the pressure of the nodes in the natural gas network considering the active power change of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

and

change of gas consumption of gas power station respectively for node u and node r

The pressure of time;

(5) the method for establishing the distributed optimization scheduling model of the electric coupling system based on the consideration of uncertainty transfer comprises the following steps:

(5-1) setting the coordination vector sent by the power grid to the natural gas grid comprises the following steps: gas consumption vector GD of gas power station in power grid^ngu,pMinimum vector delta GD of variation vector of gas consumption of gas power station^ngu,minMaximum vector delta GD of variation vector of gas consumption of gas power station^ngu,max(ii) a Wherein, GD^ngu,pIs an H x 1 column vector, GD^ngu,pThe value of any element in is

h＝1,2…H；ΔGD^ngu,minIs an H1 column vector, Δ GD^ngu,minThe value of any element in is

h＝1,2…H；ΔGD^ngu,maxIs an H1 column vector, Δ GD^ngu,maxThe value of any element in is

H is 1,2 … H; defining a coordination vector sent by a natural gas network to a power grid as follows: gas consumption vector GD of gas power station in natural gas network^ngu,g，GD^ngu,gIs a column vector of dimension H x 1, GD^ngu,gThe value of any element in is

(5-2) defining a Lagrange multiplier column vector lambda, the dimension of the lambda is H multiplied by 1, and the value of any element in the lambda is marked as lambda_h，h＝1,2…H；

(5-3) defining a punishment factor column vector rho, the dimension of rho being H1, and the value of any element in rho is recorded as rho_h，h＝1,2…H，ρ_hThe value range is 0.1-10, and rho is given by power grid dispatching;

(5-4) defining a convergence threshold column vector ε₁And ε₂，ε₁Has dimension of H × 1, ε₁The value of any element is marked as epsilon_1,h，h＝1,2…H，ε_1,hThe value range is 0.001-0.1; epsilon₂The value of any element is marked as epsilon_2,h，h＝1,2…H，ε_2,hThe value range is 0.001-0.1; epsilon₁And ε₂Dispatching and giving by a power grid;

(5-5) initializing the power grid, setting the number of initialization iterations z to 0, and initializing

Initialization

Initialization

Initialization

And GD to be initialized^ngu,p,z、ΔGD^ngu,min,z、ΔGD^ngu,max,z、λ^zρ is sent to the natural gas network;

(5-6) GD sent by power grid received by natural gas network^ngu,p,z、ΔGD^ngu,min,z、ΔGD^ngu,max,z、λ^zAnd after rho, establishing a natural gas network optimization model, wherein the objective function of the natural gas network optimization model is as follows:

in the formula, GPC is the operation cost of the natural gas network, and the calculation formula is as follows:

in the formula, PRI_sThe unit gas production cost of the natural gas wells is given by the scheduling of a natural gas network; s is the number of natural gas wells;

the constraint condition of the natural gas network optimization model is the constraint condition established in the step (4);

(5-7) solving the natural gas network optimization model in the step (5-6) by using an interior point method to obtain the gas consumption vector GD of each gas power station^ngu,g，GD^ngu,gHas dimension of H × 1, GD^ngu,gThe numerical value of the element in (1) is

And GD is to be^{ngu ,g}Is marked as GD^ngu,g,z+1The natural gas network will consume the gas consumption GD of each gas power station^ngu,g,z+1Sending the data to a power grid;

(5-8) GD sent by natural gas network received by power network^ngu,g,z+1Then, establishing a power grid optimization model, wherein the objective function of the power grid optimization model is as follows:

in the formula, the PGC is the power grid operation cost, and the PGC calculation formula is as follows:

in the formula (I), the compound is shown in the specification,

the operating cost of the gas power station h;

the operation cost of the non-gas power station I, wherein I is the number of the non-gas power stations;

for wind-power generation unitsJ is a wind abandon penalty factor which is a known quantity and is given by power grid dispatching, and J is the number of the wind turbine generators;

the load abandoning factor of the electric load K is a known quantity and is given by the dispatching of the power grid, and K is the number of the electric loads;

operating cost of gas power plant h

The calculation method comprises the following steps:

in the formula (I), the compound is shown in the specification,

and

the secondary cost coefficient, the primary cost coefficient and the fixed cost coefficient of the gas power station h are provided by a gas power station manufacturer;

the calculation method comprises the following steps:

in the formula (I), the compound is shown in the specification,

and

the secondary cost coefficient, the primary cost coefficient and the fixed cost coefficient of the non-gas power station i are respectively provided by a non-gas power station manufacturer;

the constraint condition of the power grid optimization model is the constraint condition established in the step (3);

(5-9) solving the power grid optimization model obtained in the step (5-8) by using an interior point method to obtain the gas consumption vector GD of each gas power station^ngu,pVector GD of gas consumption^ngu,pIs marked as GD^ngu,p,z+1Obtaining the minimum value vector delta GD of the variation vector of the gas consumption of the gas power station^ngu,minVector of minimum value Δ GD^ngu,minIs noted as Δ GD^ngu,min,z+1Obtaining the maximum value vector Delta GD of the variation vector of the gas consumption of the gas power station^ngu,maxVector of maximum values Δ GD^ngu,maxIs noted as Δ GD^ngu,max,z+1；

(5-10) the power grid judges the operation result, if the operation result meets the following conditions:

and is

Performing the step (6);

if it satisfies

Or satisfy

Then the lagrange multiplier is updated as follows:

electric network to GD^ngu,p,z+1、ΔGD^ngu,min,z+1、ΔGD^ngu,max,z+1、λ^z+1Transferring to a natural gas network, and enabling z to be z +1, and returning to the step (5-6);

(6) solving the distributed optimized dispatching model of the electric coupling system based on the consideration of uncertainty transfer in the step (5) to obtain the values of the variables to be solved of the power grid and the natural gas grid, namely the active power of the gas power station h in the power grid

H gas consumption of gas power station

Active power of non-gas power station i

Abandon active power WS of wind turbine generator j_jElectric power LS of electric load k_kActive power pf of branch between node l and node m_lmPhase angle theta of voltage at node l_lGas outlet flow G of natural gas well s in natural gas network_sH gas consumption of gas power station

Natural gas flow gf of pipeline between node u and node r_urAnd pressure ω of node r_rAnd taking the value of the variable to be solved as a parameter for distributed optimized operation of the electrical coupling system, and realizing distributed optimized scheduling of the electrical coupling system considering uncertainty transfer.

The distributed optimal scheduling method of the electrical coupling system considering the uncertainty transfer, provided by the invention, has the advantages that:

the distributed optimal scheduling method of the electric coupling system considering uncertainty transfer improves the existing scheduling method of the electric coupling system, fully considers the uncertainty transfer in the electric network and the natural gas network in the optimization process of the electric coupling system, and establishes the distributed optimal model of the electric network and the natural gas network, thereby obtaining more reasonable parameters for optimal operation of the electric coupling system. In the method, a power grid constraint condition and a natural gas grid constraint condition which take the uncertainty of the wind power active power into consideration are established; establishing a power grid and natural gas grid distributed optimization model; the method for optimizing information interaction of the power grid and the natural gas grid is provided, and by continuously interacting optimization information, a plurality of uncertainties caused by injection of high-proportion renewable energy sources in the electric coupling system are overcome, so that distributed optimization of the electric coupling system is finally realized, and safe, reliable and economic operation of the electric coupling system is realized.

Detailed Description

The invention provides an electric coupling system distributed optimization scheduling method considering uncertainty transfer, which comprises the following steps:

(1) the electric coupling system is divided into a power grid and a natural gas grid, and the power grid and the natural gas grid are coupled through H gas power stations;

(2) recording the value of the gas consumption of the gas power station h in the power grid

The value of the gas consumption of the gas power station h in the natural gas network is recorded as

And

the following relation is satisfied:

(3) establishing a power grid constraint condition, comprising the following steps:

(3-1) establishing the electric quantity balance constraint of the power grid node as follows:

in the formula, M is the node serial number in the power grid, M is the total number of nodes in the power grid,

the active power of the gas power station h, the variable to be solved,

the sum of the active power of all the gas power stations connected with the node m in the power grid;

is the active power of the non-gas power station i, is a variable to be solved,

the sum of the active power of all the non-gas power stations connected with the node m in the power grid;

is a predicted value of the active power of the wind turbine generator j, is a known quantity and is given by the dispatching of a power grid, WS_jThe active power is used for abandoning the wind turbine generator j, the variable to be solved is obtained,

the sum of active power actually injected into the power grid by all wind turbine generators connected with the node m in the power grid is obtained; PD (photo diode)_kFor the predicted value of the active power of the electrical load k, given by the grid schedule, LS, for a known quantity_kThe power of the electric load k is the power-abandoning active power, which is a variable to be solved,

the sum of the actual active power of all the electric loads connected with the node m in the power grid; pf_lmThe active power of the branch between the node l and the node m is defined as a variable to be solved, the flow from the node l to the node m is positive, the flow from the node m to the node l is negative,

the sum of the active power flowing to the node m from all the nodes l connected with the node m in the power grid;

(3-2) establishing a power grid direct current power flow constraint as follows:

in the formula, theta_lAnd theta_mVoltage phase angles of a node l and a node m are respectively used as variables to be solved; x is the number of_lmRepresents the reactance of a branch connected between node l and node m, and is a known quantityThe power grid dispatching is given;

(3-3) establishing a voltage phase angle constraint of a reference node in the power grid as follows:

θ_n＝0,n∈REF^p

in the formula, theta_nRepresenting the phase angle, REF, of node n in the grid^pRepresenting a set of reference nodes in the grid, given by grid dispatch;

(3-4) establishing the upper limit constraint and the lower limit constraint of the active power of the gas power station in the power grid as follows:

in the formula (I), the compound is shown in the specification,

and

respectively setting the upper limit of active power and the lower limit of active power of the gas power station h by power grid dispatching;

(3-5) establishing the upper limit constraint and the lower limit constraint of the active power of the non-gas power station in the power grid as follows:

in the formula (I), the compound is shown in the specification,

and

respectively setting an upper active power limit and a lower active power limit of the non-gas power station i by power grid dispatching;

(3-6) establishing the upper limit constraint and the lower limit constraint of the abandoned power of the wind turbine generator as follows:

in the formula (I), the compound is shown in the specification,

the upper limit of the active power for abandoning the wind turbine generator j is given by the dispatching of a power grid;

(3-7) establishing the upper limit constraint and the lower limit constraint of the electric load electricity abandoning active power as follows:

0≤LS_k≤LS_k ^max

in the formula, LS_k ^maxAbandoning an upper limit of active power for the electric load k, and dispatching and giving the upper limit by a power grid;

(3-8) establishing upper limit constraint and lower limit constraint of branch active power in the power grid as follows:

-pf_lm ^max≤pf_lm≤pf_lm ^max

in the formula, pf_lm ^maxThe upper limit of the active power of the branch between the node l and the node m is given by the dispatching of a power grid;

(3-9) establishing the active power of the gas power plant h

And gas consumption

The constraints between are as follows:

in the formula, a_h ^ngu、b_h ^nguAnd c_h ^nguQuadratic term coefficients, primary term coefficients and constant terms which are quadratic relational expressions of the active power and the gas consumption of the gas power station h are respectively given by the gas power station;

(3-10) setting an active power variation interval of the wind turbine generator as follows:

wherein the content of the first and second substances,

the upper limit of the active power fluctuation of the wind turbine generator j is a known quantity and is given by the dispatching of the power grid,

the lower limit of the active power fluctuation of the wind turbine generator j is a known quantity and is given by the dispatching of the power grid,

the active power of the wind turbine generator j is obtained;

(3-11) setting quasi-steady-state output power transfer distribution factors of all power stations, namely gas power stations, non-gas power stations and wind power generation sets in power grid

In the formula (I), the compound is shown in the specification,

the vector is an Nx 1-dimensional column vector, N is the total number of power stations in a power grid, including a gas power station, a non-gas power station and a wind power generation unit, and the upper standard R is a quasi-steady-state identifier; h_lmFor each station n to the active power pf of the branch between node l and node m_lmN x 1-dimensional column vector, I, formed by the transfer distribution factors of_NIs an NxN dimensional identity matrix, alpha_NAn N x 1-dimensional column vector composed of the bearing coefficients bearing the unbalanced power for each station N,

is an N × 1 dimensional column vector with all values of 1;

above H_lmThe elements in (a) are represented as:

in the formula, n^pNumbering nodes of a power station n in a power grid, X being an impedance matrix of the power grid,

is the l-th row and the n-th row in the impedance matrix X^pThe elements of the column are,

is the m-th row and the n-th row in the impedance matrix X^pThe elements of the column are,

and

are all given by the power grid dispatching;

coefficient of bearing alpha_NIn the wind turbine, the bearing coefficient is 0, alpha_NThe bearing coefficient of the medium gas power station and the non-gas power station is more than 0, alpha_NGiven by the grid schedule and satisfying the following relationship:

(3-12) setting the active power adjustment vector of each power station in the power grid caused by the active power change of the wind turbine generator as

Is an N x 1-dimensional column vector,

the element values of the corresponding wind generating set are as follows:

j^psfor the numbering of the wind turbine j in each station in the grid,

taking the element values of corresponding gas power stations and non-gas power stations as 0;

(3-13) calculating the active power pf of the branch between the node l and the node m by using the following formula_lmAmount of change of

Is expressed as [ Delta pf ]_lm ^min,Δpf_lm ^max]，Δpf_lm ^minIs composed of

Lower limit value of value range, Δ pf_lm ^maxIs composed of

Taking the upper limit value of the value interval;

(3-14) establishing the active power upper limit and active power lower limit constraints of the branch in the power grid considering the active power change of the wind turbine generator as follows:

-pf_lm ^max≤pf_lm+Δpf_lm ^min≤pf_lm ^max

-pf_lm ^max≤pf_lm+Δpf_lm ^max≤pf_lm ^max

(3-15) calculating the active power regulating quantity of each power station in the power grid

In the formula (I), the compound is shown in the specification,

is an N x 1-dimensional column vector,

any of the elements of

Has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

Taking the upper limit value of the value interval;

(3-16) establishing the active power upper limit and the active power lower limit of a gas power station in a power grid considering the active power change of the wind turbine generator as follows:

in the formula, h^psFor gas power stations h, individual stations in the gridThe serial number in (1) is (d),

is composed of

H in (1)^psThe lower limit value of the value interval of each element,

is composed of

H in (1)^psThe upper limit value of the value interval of each element;

(3-17) establishing the active power upper limit and the active power lower limit of a non-gas power station in the power grid considering the active power change of the wind turbine generator as follows:

in the formula i^psThe serial numbers of the non-gas power stations i in each power station in the power grid,

is composed of

I of (1)^psThe lower limit value of the value interval of each element,

is composed of

I of (1)^psThe upper limit value of the value interval of each element;

(3-18) establishing the following constraint between the active power variation and the air consumption variation of the gas power station h considering the active power variation of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

is the variation of the gas consumption caused by the variation of the active power of the gas power station h,

has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

The upper limit value of the value-taking interval,

and

the calculation formula is as follows:

(4) establishing natural gas network constraint conditions, comprising the following steps:

(4-1) establishing the natural gas flow balance constraint of the natural gas network node as follows:

in the formula, R is the node number in the natural gas network, and R is the node number in the natural gas network; s is the number of the gas well in the natural gas network, G_sThe outlet gas flow of the natural gas well s is used as a variable to be solved,

the sum of the gas outlet flow of all the natural gas wells connected with the node r in the natural gas network; t is the number of the gas load of residents in the natural gas network,

the gas consumption for the resident gas load t, which is a known quantity, is given by the natural gas network dispatch,

the gas consumption of the gas load of all residents connected with the node r in the natural gas network;

the sum of the gas consumption of all gas power stations connected with the node r in the natural gas network is a variable to be solved; gf_urFor the natural gas flow of the pipeline between the node u and the node r in the natural gas network, the gf when the natural gas flows from the node u to the gas point r is specified as a variable to be solved_urTake positive value, gf when flowing from node r to node u_urTaking the negative value of the reaction mixture,

the natural gas flow of all nodes connected with the node r in the natural gas network flows into the node r;

(4-2) establishing upper and lower constraints on node pressures in the natural gas network as follows:

in the formula (I), the compound is shown in the specification,

and

respectively representing the lower pressure limit and the upper pressure limit of a node r in the natural gas network, and being given by the scheduling of the natural gas network;

(4-3) establishing the relationship between the natural gas flow and pressure in the natural gas grid as follows:

in the formula, ω_uAnd ω_rPressure at node u and node r, sgn (ω), respectively, in the natural gas network_u,ω_r) Is about ω_u、ω_rFunction of when ω is_u＞ω_rThen, sgn (ω)_u,ω_r) Take 1 when ω_u≤ω_r，sgn(ω_u,ω_r) The value is 0; c_urIs the Welmos constant of the pipeline between node u and node r, a known quantity, given by the natural gas grid schedule, due to sgn (ω) in the constraint on the relationship between natural gas flow and pressure_u,ω_r) Is a binary variable, and an integer variable is introduced

The following relation is satisfied:

defining a mathematical variable F_ur，

And translating the constraint on the relationship between natural gas flow and pressure into the following expression:

in the formula (I), the compound is shown in the specification,

and

respectively setting a lower pressure limit and an upper pressure limit of the node u by scheduling of a natural gas network;

wherein the above-mentioned

The further relaxation is:

(4-4) establishing a pressure reference node constraint in the natural gas network as follows:

ω_v＝PR,v∈REF^g

in the formula, ω_vRepresenting the pressure at node v, PR is a constant given by the natural gas grid schedule, REF^gA set of reference nodes representing a natural gas network, given by a natural gas network schedule;

(4-5) establishing the upper limit constraint and the lower limit constraint of the gas well effluent flow in the natural gas network as follows:

wherein S is the number of natural gas wells,

the lower limit and the upper limit of the gas flow rate of the natural gas well are respectively given by the dispatching of a natural gas network;

(4-6) establishing a natural gas flow constraint for the pipes in the natural gas network as follows:

in the formula (I), the compound is shown in the specification,

the upper limit of the natural gas flow of the pipeline between the node u and the node r in the natural gas network is given by the scheduling of the natural gas network;

(4-7) defining nodes connected with the natural gas well and the gas power station in the natural gas network as gas injection quantity variable nodes, marking as w, and marking the number of all the gas injection quantity variable nodes in the natural gas network as Q;

(4-8) setting a transfer distribution factor of the quasi-steady-state gas outlet flow of each variable gas injection amount node in the natural gas network

In the formula (I), the compound is shown in the specification,

is a Q x 1-dimensional column vector, K_urIs a Q x 1-dimensional column vector, K_urThe natural gas flow gf of the pipeline between the node u and the node r in the natural gas network is jointed by each variable gas injection amount node_urComposition of transfer distribution factor of (I)_QIs a QxQ dimensional identity matrix, alpha_QIs a Q x 1-dimensional column vector, alpha_QThe device consists of bearing coefficients of unbalanced natural gas flow borne by each variable gas injection quantity node,

a Q × 1 dimensional column vector with all values being 1;

wherein K_urThe elements in (a) are represented as:

where k is the number of the pipe between node u and node r, T is the road-pipe correlation matrix, T is a known quantity, given by the natural gas network schedule,

for the qth in the road-pipe association matrix T^gRow, k column element, q^gNumbering nodes of the variable gas injection quantity node q in the natural gas network;

coefficient of bearing alpha_QThe value rule of the elements is as follows: bearing coefficient alpha of natural gas well_QGreater than 0, the gas power plant has a bearing coefficient alpha_QIs equal to 0, alpha_QGiven by the natural gas grid schedule and satisfying the following relation:

(4-9) recording the Q multiplied by 1 dimension row vector formed by the variable gas consumption interval of each gas injection variable node as

Wherein the row value of the gas injection quantity variable node corresponding to the gas power station h is

The corresponding row value of other natural gas wells is 0;

(4-10) calculating the natural gas flow gf of the pipeline between the node u and the node r in the natural gas network using the following formula_urAmount of change of

In the formula (I), the compound is shown in the specification,

has a value interval of [ delta gf ]_ur ^min,Δgf_ur ^max]，Δgf_ur ^minIs composed of

Lower limit value of interval, delta gf_ur ^maxIs composed of

Taking the upper limit value of the value interval;

(4-11) establishing the constraint of the natural gas flow rate of the pipeline in the natural gas network considering the active power change of the wind turbine generator as follows:

(4-12) calculating the adjustment amount of the natural gas flow at each variable node of the gas injection amount by using the following formula

In the formula (I), the compound is shown in the specification,

is a Q x 1 dimensional column vector,

any one element of

Has a value interval of [ Delta G ]_q ^min,ΔG_q ^max]，ΔG_q ^minIs composed of

Lower limit of value range, Δ G_q ^maxIs composed of

Taking the upper limit value of the value interval;

(4-13) establishing the upper limit constraint and the lower limit constraint of the gas flow of the natural gas well in the natural gas network considering the active power change of the wind turbine generator as follows:

in the formula, s^gsNumbering all gas injection quantity variable nodes of the natural gas wells in a natural gas network;

is composed of

Zhongth s^gsThe lower limit value of the value interval of each element,

is composed of

S of (1)^gsThe upper limit value of the value interval of each element;

(4-14) establishing the change quantity of the gas consumption quantity of the gas power station in consideration of the active power change of the wind turbine generator

The safety constraint of the node pressure in the natural gas network when the boundary value is taken comprises the following steps:

(4-14-1) setting the gas consumption of the gas power station to rise

In the time, the Q multiplied by 1 dimension row vector formed by the variable gas consumption interval of each gas injection variable node is recorded as delta G^(0),upWherein the row value of the gas injection variable node corresponding to the gas power station h

The corresponding row value of the natural gas well is 0;

(4-14-2) calculating the natural gas flow rate gf of the pipeline between the node u and the node r in the natural gas network by using the following formula_urAmount of change of

(4-14-3) defining mathematical variables

The method comprises the following steps of establishing the following safety constraints of the pressure of the nodes in the natural gas network considering the active power change of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

and

gas consumption rise of node u and node r in gas power station

The pressure of time;

(4-14-4) setting the gas consumption variation of the gas power station

In the time, the Q multiplied by 1 dimension column vector formed by the variable gas consumption interval of each gas injection variable node is delta G^(0),downWherein the row value of the gas injection variable node corresponding to the gas power station h

The corresponding row value of the natural gas well is 0;

(4-14-5) calculating the natural gas flow gf of the pipeline between the node u and the node r in the natural gas network_urAmount of change of

(4-14-6) defining mathematical variables

The method comprises the following steps of establishing the following safety constraints of the pressure of the nodes in the natural gas network considering the active power change of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

and

change of gas consumption of gas power station respectively for node u and node r

The pressure of time;

(5) the method for establishing the distributed optimization scheduling model of the electric coupling system based on the consideration of uncertainty transfer comprises the following steps:

(5-1) setting the coordination vector sent by the power grid to the natural gas grid comprises the following steps: gas consumption vector GD of gas power station in power grid^ngu,pMinimum vector delta GD of variation vector of gas consumption of gas power station^ngu,minMaximum vector delta GD of variation vector of gas consumption of gas power station^ngu,max(ii) a Wherein, GD^ngu,pIs an H x 1 column vector, GD^ngu,pThe value of any element in is

h＝1,2…H；ΔGD^ngu,minIs an H1 column vector, Δ GD^ngu,minThe value of any element in is

h＝1,2…H；ΔGD^ngu,maxIs an H1 column vector, Δ GD^ngu,maxThe value of any element in is

H is 1,2 … H; defining a coordination vector sent by a natural gas network to a power grid as follows: gas consumption vector GD of gas power station in natural gas network^ngu,g，GD^ngu,gIs a column vector of dimension H x 1, GD^ngu,gThe value of any element in is

(5-2) defining a Lagrange multiplier column vector lambda, the dimension of the lambda is H multiplied by 1, and the value of any element in the lambda is marked as lambda_h，h＝1,2…H；

(5-3) defining a punishment factor column vector rho, wherein the dimension of rho is H multiplied by 1, and the value of any element in rho is marked as rho_h，h＝1,2…H，ρ_hThe value range is 0.1-10, and rho is given by power grid dispatching;

(5-4) defining a convergence threshold column vector ε₁And ε₂，ε₁Has dimension of H × 1, ε₁The value of any element is marked as epsilon_1,h，h＝1,2…H，ε_1,hThe value range is 0.001-0.1; epsilon₂The value of any element is marked as epsilon_2,h，h＝1,2…H，ε_2,hThe value range is 0.001-0.1; epsilon₁And ε₂Dispatching and giving by a power grid;

(5-5) initializing the power grid, setting the number of initialization iterations z to 0, and initializing

Initialization

Initialization

Initialization

And GD to be initialized^ngu,p,z、ΔGD^ngu,min,z、ΔGD^ngu,max,z、λ^zρ is sent to the natural gas network;

(5-6) GD sent by power grid received by natural gas network^ngu,p,z、ΔGD^ngu,min,z、ΔGD^ngu,max,z、λ^zAnd after rho, establishing a natural gas network optimization model, wherein the objective function of the natural gas network optimization model is as follows:

in the formula, GPC is the operation cost of the natural gas network, and the calculation formula is as follows:

in the formula, PRI_sThe unit gas production cost of the natural gas wells is given by the scheduling of a natural gas network; s is the number of natural gas wells;

the constraint condition of the natural gas network optimization model is the constraint condition established in the step (4);

(5-7) solving the natural gas network optimization model in the step (5-6) by using an interior point method to obtain the gas consumption vector GD of each gas power station^ngu,g，GD^ngu,gHas dimension of H × 1, GD^ngu,gThe numerical value of the element in (1) is

And GD is to be^{ngu ,g}Is marked as GD^ngu,g,z+1The natural gas network will consume the gas consumption GD of each gas power station^ngu,g,z+1Sending the data to a power grid;

(5-8) GD sent by natural gas network received by power network^ngu,g,z+1Then, establishing a power grid optimization model, wherein the objective function of the power grid optimization model is as follows:

in the formula, the PGC is the power grid operation cost, and the PGC calculation formula is as follows:

in the formula (I), the compound is shown in the specification,

the operating cost of the gas power station h;

the operation cost of the non-gas power station I, wherein I is the number of the non-gas power stations;

a wind abandon penalty factor of a wind turbine generator J is a known quantity and is given by power grid dispatching, and J is the number of the wind turbine generators;

the load abandoning factor of the electric load K is a known quantity and is given by the dispatching of the power grid, and K is the number of the electric loads;

operating cost of gas power plant h

The calculation method comprises the following steps:

in the formula (I), the compound is shown in the specification,

and

the secondary cost coefficient, the primary cost coefficient and the fixed cost coefficient of the gas power station h are provided by a gas power station manufacturer;

the calculation method comprises the following steps:

in the formula (I), the compound is shown in the specification,

and

the secondary cost coefficient, the primary cost coefficient and the fixed cost coefficient of the non-gas power station i are respectively provided by a non-gas power station manufacturer;

the constraint condition of the power grid optimization model is the constraint condition established in the step (3);

(5-9) solving the power grid optimization model obtained in the step (5-8) by using an interior point method to obtain the gas consumption vector GD of each gas power station^ngu,pVector GD of gas consumption^ngu,pIs marked as GD^ngu,p,z+1Obtaining the minimum value vector delta GD of the variation vector of the gas consumption of the gas power station^ngu,minVector of minimum value Δ GD^ngu,minIs noted as Δ GD^ngu,min,z+1Obtaining the maximum value vector Delta GD of the variation vector of the gas consumption of the gas power station^ngu,maxVector of maximum values Δ GD^ngu,maxIs noted as Δ GD^ngu,max,z+1；

(5-10) the power grid judges the operation result, if the operation result meets the following conditions:

and is

Performing the step (6);

if it satisfies

Or satisfy

Then the lagrange multiplier is updated as follows:

electric network to GD^ngu,p,z+1、ΔGD^ngu,min,z+1、ΔGD^ngu,max,z+1、λ^z+1Transferring to a natural gas network, and enabling z to be z +1, and returning to the step (5-6);

(6) solving forObtaining the values of variables to be solved of the power grid and the natural gas grid based on the distributed optimization scheduling model of the electric coupling system considering the uncertainty transfer in the step (5), namely the active power of the gas power station h in the power grid

H gas consumption of gas power station

Active power of non-gas power station i

Abandon active power WS of wind turbine generator j_jElectric power LS of electric load k_kActive power pf of branch between node l and node m_lmPhase angle theta of voltage at node l_lGas outlet flow G of natural gas well s in natural gas network_sH gas consumption of gas power station

Natural gas flow gf of pipeline between node u and node r_urAnd pressure ω of node r_rAnd taking the value of the variable to be solved as a parameter for distributed optimized operation of the electrical coupling system, and realizing distributed optimized scheduling of the electrical coupling system considering uncertainty transfer.

Claims (1)

1. An electric coupling system distributed optimization scheduling method considering uncertainty transfer is characterized by comprising the following steps:

(1) the electric coupling system is divided into a power grid and a natural gas grid, and the power grid and the natural gas grid are coupled through H gas power stations;

(2) recording the value of the gas consumption of the gas power station h in the power grid

The value of the gas consumption of the gas power station h in the natural gas network is recorded as

And

the following relation is satisfied:

(3) establishing a power grid constraint condition, comprising the following steps:

(3-1) establishing the electric quantity balance constraint of the power grid node as follows:

in the formula, M is the node serial number in the power grid, M is the total number of nodes in the power grid, and P_h ^nguThe active power of the gas power station h, the variable to be solved,

the sum of the active power of all the gas power stations connected with the node m in the power grid; p_i ^genIs the active power of the non-gas power station i, is a variable to be solved,

the sum of the active power of all the non-gas power stations connected with the node m in the power grid;

is a predicted value of the active power of the wind turbine generator j, is a known quantity and is given by the dispatching of a power grid, WS_jThe active power is used for abandoning the wind turbine generator j, the variable to be solved is obtained,

for all wind turbine generators connected with node m in power grid, actually injecting the wind turbine generators into the power gridThe sum of the work power; PD (photo diode)_kFor the predicted value of the active power of the electrical load k, given by the grid schedule, LS, for a known quantity_kThe power of the electric load k is the power-abandoning active power, which is a variable to be solved,

the sum of the actual active power of all the electric loads connected with the node m in the power grid; pf_lmThe active power of the branch between the node l and the node m is defined as a variable to be solved, the flow from the node l to the node m is positive, the flow from the node m to the node l is negative,

the sum of the active power flowing to the node m from all the nodes l connected with the node m in the power grid;

(3-2) establishing a power grid direct current power flow constraint as follows:

in the formula, theta_lAnd theta_mVoltage phase angles of a node l and a node m are respectively used as variables to be solved; x is the number of_lmRepresenting the reactance of a branch connected between the node l and the node m, wherein the reactance is a known quantity and is given by power grid dispatching;

(3-3) establishing a voltage phase angle constraint of a reference node in the power grid as follows:

θ_n＝0,n∈REF^p

in the formula, theta_nRepresenting the phase angle, REF, of node n in the grid^pRepresenting a set of reference nodes in the grid, given by grid dispatch;

(3-4) establishing the upper limit constraint and the lower limit constraint of the active power of the gas power station in the power grid as follows:

in the formula (I), the compound is shown in the specification,

and P_h ^ngu,maxThe lower limit of the active power and the upper limit of the active power of the gas power station h are respectively given by the dispatching of a power grid;

(3-5) establishing the upper limit constraint and the lower limit constraint of the active power of the non-gas power station in the power grid as follows:

P_i ^gen,min≤P_i ^gen≤P_i ^gen,max

in the formula, P_i ^gen,minAnd P_i ^gen,maxRespectively setting the lower limit of active power and the upper limit of active power of a non-gas power station i by power grid dispatching;

(3-6) establishing the upper limit constraint and the lower limit constraint of the abandoned power of the wind turbine generator as follows:

in the formula (I), the compound is shown in the specification,

the upper limit of the active power for abandoning the wind turbine generator j is given by the dispatching of a power grid;

(3-7) establishing the upper limit constraint and the lower limit constraint of the electric load electricity abandoning active power as follows:

0≤LS_k≤LS_k ^max

in the formula, LS_k ^maxAbandoning an upper limit of active power for the electric load k, and dispatching and giving the upper limit by a power grid;

(3-8) establishing upper limit constraint and lower limit constraint of branch active power in the power grid as follows:

-pf_lm ^max≤pf_lm≤pf_lm ^max

in the formula, pf_lm ^maxThe upper limit of the active power of the branch between the node l and the node m is given by the dispatching of a power grid;

(3-9) establishing the active power of the gas power plant h

And gas consumption

The constraints between are as follows:

in the formula, a_h ^ngu、b_h ^nguAnd c_h ^nguQuadratic term coefficients, primary term coefficients and constant terms which are quadratic relational expressions of the active power and the gas consumption of the gas power station h are respectively given by the gas power station;

(3-10) setting an active power variation interval of the wind turbine generator as follows:

wherein the content of the first and second substances,

the upper limit of the active power fluctuation of the wind turbine generator j is a known quantity and is given by the dispatching of the power grid,

the lower limit of the active power fluctuation of the wind turbine generator j is a known quantity and is given by the dispatching of the power grid,

the active power of the wind turbine generator j is obtained;

(3-11) setting quasi-steady-state output power transfer distribution factors of all power stations, namely gas power stations, non-gas power stations and wind power generation sets in power grid

In the formula (I), the compound is shown in the specification,

the vector is an Nx 1-dimensional column vector, N is the total number of power stations in a power grid, including a gas power station, a non-gas power station and a wind power generation unit, and the upper standard R is a quasi-steady-state identifier; h_lmFor each station n to the active power pf of the branch between node l and node m_lmN x 1-dimensional column vector, I, formed by the transfer distribution factors of_NIs an NxN dimensional identity matrix, alpha_NAn N x 1-dimensional column vector composed of the bearing coefficients bearing the unbalanced power for each station N,

is an N × 1 dimensional column vector with all values of 1;

above H_lmThe elements in (a) are represented as:

in the formula, n^pNumbering nodes of a power station n in a power grid, X being an impedance matrix of the power grid,

is the l-th row and the n-th row in the impedance matrix X^pThe elements of the column are,

is the m-th row and the n-th row in the impedance matrix X^pThe elements of the column are,

and

are all dispatched to by the power gridDetermining;

coefficient of bearing alpha_NIn the wind turbine, the bearing coefficient is 0, alpha_NThe bearing coefficient of the medium gas power station and the non-gas power station is more than 0, alpha_NGiven by the grid schedule and satisfying the following relationship:

(3-12) setting the active power adjustment vector of each power station in the power grid caused by the active power change of the wind turbine generator as

Is an N x 1-dimensional column vector,

the element values of the corresponding wind generating set are as follows:

j^psfor the numbering of the wind turbine j in each station in the grid,

taking the element values of corresponding gas power stations and non-gas power stations as 0;

(3-13) calculating the active power pf of the branch between the node l and the node m by using the following formula_lmAmount of change of

Is expressed as [ Delta pf ]_lm ^min,Δpf_lm ^max]，Δpf_lm ^minIs composed of

Lower limit value of value range, Δ pf_lm ^maxIs composed of

Taking the upper limit value of the value interval;

(3-14) establishing the active power upper limit and active power lower limit constraints of the branch in the power grid considering the active power change of the wind turbine generator as follows:

-pf_lm ^max≤pf_lm+Δpf_lm ^min≤pf_lm ^max

-pf_lm ^max≤pf_lm+Δpf_lm ^max≤pf_lm ^max

(3-15) calculating the active power regulating quantity of each power station in the power grid

In the formula (I), the compound is shown in the specification,

is an N x 1-dimensional column vector,

any of the elements of

Has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

Taking the upper limit value of the value interval;

(3-16) establishing the active power upper limit and the active power lower limit of a gas power station in a power grid considering the active power change of the wind turbine generator as follows:

in the formula, h^psThe serial numbers of the gas power stations h in each power station in the power grid,

is composed of

H in (1)^psThe lower limit value of the value interval of each element,

is composed of

H in (1)^psThe upper limit value of the value interval of each element;

(3-17) establishing the active power upper limit and the active power lower limit of a non-gas power station in the power grid considering the active power change of the wind turbine generator as follows:

in the formula i^psThe serial numbers of the non-gas power stations i in each power station in the power grid,

is composed of

I of (1)^psThe lower limit value of the value interval of each element,

is composed of

I of (1)^psThe upper limit value of the value interval of each element;

(3-18) establishing the following constraint between the active power variation and the air consumption variation of the gas power station h considering the active power variation of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

is the variation of the gas consumption caused by the variation of the active power of the gas power station h,

has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

The upper limit value of the value-taking interval,

and

the calculation formula is as follows:

(4) establishing natural gas network constraint conditions, comprising the following steps:

(4-1) establishing the natural gas flow balance constraint of the natural gas network node as follows:

in the formula, R is the node number in the natural gas network, and R is the node number in the natural gas network; s is the number of the gas well in the natural gas network, G_sThe outlet gas flow of the natural gas well s is used as a variable to be solved,

the sum of the gas outlet flow of all the natural gas wells connected with the node r in the natural gas network; t is the number of the gas load of residents in the natural gas network,

the gas consumption for the resident gas load t, which is a known quantity, is given by the natural gas network dispatch,

the gas consumption of the gas load of all residents connected with the node r in the natural gas network;

the sum of the gas consumption of all gas power stations connected with the node r in the natural gas network is a variable to be solved; gf_urFor the natural gas flow of the pipeline between the node u and the node r in the natural gas network, the gf when the natural gas flows from the node u to the gas point r is specified as a variable to be solved_urTake positive value, gf when flowing from node r to node u_urTaking the negative value of the reaction mixture,

the natural gas flow of all nodes connected with the node r in the natural gas network flows into the node r;

(4-2) establishing upper and lower constraints on node pressures in the natural gas network as follows:

in the formula (I), the compound is shown in the specification,

and

respectively representing the lower pressure limit and the upper pressure limit of a node r in the natural gas network, and being given by the scheduling of the natural gas network;

(4-3) establishing the relationship between the natural gas flow and pressure in the natural gas grid as follows:

in the formula, ω_uAnd ω_rPressure at node u and node r, sgn (ω), respectively, in the natural gas network_u,ω_r) Is about ω_u、ω_rFunction of when ω is_u＞ω_rThen, sgn (ω)_u,ω_r) Take 1 when ω_u≤ω_r，sgn(ω_u,ω_r) The value is 0; c_urIs the Welmos constant of the pipeline between node u and node r, a known quantity, given by the natural gas grid schedule, due to sgn (ω) in the constraint on the relationship between natural gas flow and pressure_u,ω_r) Is a binary variable, and an integer variable is introduced

The following relation is satisfied:

defining a mathematical variable F_ur，

And will beThe constraint on the relationship between natural gas flow and pressure translates into the following expression:

in the formula (I), the compound is shown in the specification,

and

respectively setting a lower pressure limit and an upper pressure limit of the node u by scheduling of a natural gas network;

wherein the above-mentioned

The further relaxation is:

(4-4) establishing a pressure reference node constraint in the natural gas network as follows:

ω_v＝PR,v∈REF^g

in the formula, ω_vRepresenting the pressure at node v, PR is a constant given by the natural gas grid schedule, REF^gA set of reference nodes representing a natural gas network, given by a natural gas network schedule;

(4-5) establishing the upper limit constraint and the lower limit constraint of the gas well effluent flow in the natural gas network as follows:

wherein S is the number of natural gas wells,

the lower limit and the upper limit of the gas flow rate of the natural gas well are respectively given by the dispatching of a natural gas network;

(4-6) establishing a natural gas flow constraint for the pipes in the natural gas network as follows:

in the formula (I), the compound is shown in the specification,

the upper limit of the natural gas flow of the pipeline between the node u and the node r in the natural gas network is given by the scheduling of the natural gas network;

(4-7) defining nodes connected with the natural gas well and the gas power station in the natural gas network as gas injection quantity variable nodes, marking as w, and marking the number of all the gas injection quantity variable nodes in the natural gas network as Q;

(4-8) setting a transfer distribution factor of the quasi-steady-state gas outlet flow of each variable gas injection amount node in the natural gas network

In the formula (I), the compound is shown in the specification,

is a Q x 1-dimensional column vector, K_urIs a Q x 1-dimensional column vector, K_urThe natural gas flow gf of the pipeline between the node u and the node r in the natural gas network is jointed by each variable gas injection amount node_urComposition of transfer distribution factor of (I)_QIs a QxQ dimensional identity matrix, alpha_QIs a Q x 1-dimensional column vector, alpha_QThe device consists of bearing coefficients of unbalanced natural gas flow borne by each variable gas injection quantity node,

a Q × 1 dimensional column vector with all values being 1;

wherein K_urThe elements in (a) are represented as:

where k is the number of the pipe between node u and node r, T is the road-pipe correlation matrix, T is a known quantity, given by the natural gas network schedule,

for the qth in the road-pipe association matrix T^gRow, k column element, q^gNumbering nodes of the variable gas injection quantity node q in the natural gas network;

coefficient of bearing alpha_QThe value rule of the elements is as follows: bearing coefficient alpha of natural gas well_QGreater than 0, the gas power plant has a bearing coefficient alpha_QIs equal to 0, alpha_QGiven by the natural gas grid schedule and satisfying the following relation:

(4-9) forming the gas consumption change interval of each variable gas injection nodeIs expressed as a Q × 1-dimensional column vector

Wherein the row value of the gas injection quantity variable node corresponding to the gas power station h is

The corresponding row value of other natural gas wells is 0;

(4-10) calculating the natural gas flow gf of the pipeline between the node u and the node r in the natural gas network using the following formula_urAmount of change of

In the formula (I), the compound is shown in the specification,

has a value interval of

Is composed of

The lower limit value of the value-taking interval,

is composed of

Taking the upper limit value of the value interval;

(4-11) establishing the constraint of the natural gas flow rate of the pipeline in the natural gas network considering the active power change of the wind turbine generator as follows:

(4-12) calculating the adjustment amount of the natural gas flow at each variable node of the gas injection amount by using the following formula

In the formula (I), the compound is shown in the specification,

is a Q x 1 dimensional column vector,

any one element of

Has a value interval of [ Delta G ]_q ^min,ΔG_q ^max]，ΔG_q ^minIs composed of

Lower limit of value range, Δ G_q ^maxIs composed of

Taking the upper limit value of the value interval;

(4-13) establishing the upper limit constraint and the lower limit constraint of the gas flow of the natural gas well in the natural gas network considering the active power change of the wind turbine generator as follows:

in the formula, s^gsNumbering all gas injection quantity variable nodes of the natural gas wells in a natural gas network;

is composed of

Zhongth s^gsThe lower limit value of the value interval of each element,

is composed of

S of (1)^gsThe upper limit value of the value interval of each element;

(4-14) establishing the change quantity of the gas consumption quantity of the gas power station in consideration of the active power change of the wind turbine generator

The safety constraint of the node pressure in the natural gas network when the boundary value is taken comprises the following steps:

(4-14-1) setting the gas consumption of the gas power station to rise

In the time, the Q multiplied by 1 dimension row vector formed by the variable gas consumption interval of each gas injection variable node is recorded as delta G^(0),upWherein the row value of the gas injection variable node corresponding to the gas power station h

The corresponding row value of the natural gas well is 0;

(4-14-2) calculating the natural gas flow rate gf of the pipeline between the node u and the node r in the natural gas network by using the following formula_urAmount of change of

(4-14-3) defining mathematical variables

The method comprises the following steps of establishing the following safety constraints of the pressure of the nodes in the natural gas network considering the active power change of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

and

gas consumption rise of node u and node r in gas power station

The pressure of time;

(4-14-4) setting the gas consumption variation of the gas power station

In the time, the Q multiplied by 1 dimension column vector formed by the variable gas consumption interval of each gas injection variable node is delta G^(0),downWherein the row value of the gas injection variable node corresponding to the gas power station h

The corresponding row value of the natural gas well is 0;

(4-14-5) calculating the natural gas flow gf of the pipeline between the node u and the node r in the natural gas network_urAmount of change of

(4-14-6) defining mathematical variables

The method comprises the following steps of establishing the following safety constraints of the pressure of the nodes in the natural gas network considering the active power change of the wind turbine generator:

in the formula (I), the compound is shown in the specification,

and

change of gas consumption of gas power station respectively for node u and node r

The pressure of time;

(5) the method for establishing the distributed optimization scheduling model of the electric coupling system based on the consideration of uncertainty transfer comprises the following steps:

(5-1) setting the coordination vector sent by the power grid to the natural gas grid comprises the following steps: gas consumption vector GD of gas power station in power grid^ngu,pMinimum vector delta of variation vector of gas consumption of gas power stationGD^ngu,minMaximum vector delta GD of variation vector of gas consumption of gas power station^ngu,max(ii) a Wherein, GD^ngu,pIs an H x 1 column vector, GD^ngu,pThe value of any element in is

h＝1,2…H；ΔGD^ngu,minIs an H1 column vector, Δ GD^ngu,minThe value of any element in is

h＝1,2…H；ΔGD^ngu,maxIs an H1 column vector, Δ GD^ngu,maxThe value of any element in is

H is 1,2 … H; defining a coordination vector sent by a natural gas network to a power grid as follows: gas consumption vector GD of gas power station in natural gas network^ngu,g，GD^ngu,gIs a column vector of dimension H x 1, GD^ngu,gThe value of any element in is

(5-2) defining a Lagrange multiplier column vector lambda, the dimension of the lambda is H multiplied by 1, and the value of any element in the lambda is marked as lambda_h，h＝1,2…H；

(5-3) defining a punishment factor column vector rho, wherein the dimension of rho is H multiplied by 1, and the value of any element in rho is marked as rho_h，h＝1,2…H，ρ_hThe value range is 0.1-10, and rho is given by power grid dispatching;

(5-4) defining a convergence threshold column vector ε₁And ε₂，ε₁Has dimension of H × 1, ε₁The value of any element is marked as epsilon_1,h，h＝1,2…H，ε_1,hThe value range is 0.001-0.1; epsilon₂The value of any element is marked as epsilon_2,h，h＝1,2…H，ε_2,hThe value range is 0.001-0.1; epsilon₁And ε₂Dispatching and giving by a power grid;

(5-5) initializing the power grid, setting the number of initialization iterations z to 0, and initializing

Initialization

Initialization

Initialization

And GD to be initialized^ngu,p,z、ΔGD^ngu,min,z、ΔGD^ngu,max,z、λ^zρ is sent to the natural gas network;

(5-6) GD sent by power grid received by natural gas network^ngu,p,z、ΔGD^ngu,min,z、ΔGD^ngu,max,z、λ^zAnd after rho, establishing a natural gas network optimization model, wherein the objective function of the natural gas network optimization model is as follows:

in the formula, GPC is the operation cost of the natural gas network, and the calculation formula is as follows:

in the formula, PRI_sThe unit gas production cost of the natural gas wells is given by the scheduling of a natural gas network; s is the number of natural gas wells;

the constraint condition of the natural gas network optimization model is the constraint condition established in the step (4);

(5-7) solving the natural gas network optimization model in the step (5-6) by using an interior point method to obtain the gas consumption vector GD of each gas power station^ngu,g，GD^ngu,gDimension of (2)Is H x 1, GD^ngu,gThe numerical value of the element in (1) is

And GD is to be^ngu,gIs marked as GD^ngu,g,z+1The natural gas network will consume the gas consumption GD of each gas power station^ngu,g,z+1Sending the data to a power grid;

(5-8) GD sent by natural gas network received by power network^ngu,g,z+1Then, establishing a power grid optimization model, wherein the objective function of the power grid optimization model is as follows:

in the formula, the PGC is the power grid operation cost, and the PGC calculation formula is as follows:

in the formula (I), the compound is shown in the specification,

the operating cost of the gas power station h;

the operation cost of the non-gas power station I, wherein I is the number of the non-gas power stations;

a wind abandon penalty factor of a wind turbine generator J is a known quantity and is given by power grid dispatching, and J is the number of the wind turbine generators;

the load abandoning factor of the electric load K is a known quantity and is given by the dispatching of the power grid, and K is the number of the electric loads;

operating cost of gas power plant h

The calculation method comprises the following steps:

in the formula (I), the compound is shown in the specification,

and f_h ^nguThe secondary cost coefficient, the primary cost coefficient and the fixed cost coefficient of the gas power station h are provided by a gas power station manufacturer;

the calculation method comprises the following steps:

in the formula (I), the compound is shown in the specification,

and f_i ^genThe secondary cost coefficient, the primary cost coefficient and the fixed cost coefficient of the non-gas power station i are respectively provided by a non-gas power station manufacturer;

the constraint condition of the power grid optimization model is the constraint condition established in the step (3);

(5-9) solving the power grid optimization model in the step (5-8) by using an interior point method to obtain a gas consumption vector GD of each gas power station^ngu,pVector GD of gas consumption^ngu,pIs marked as GD^ngu,p,z+1Obtaining the minimum value vector delta GD of the variation vector of the gas consumption of the gas power station^ngu,minVector of minimum value Δ GD^ngu,minIs noted as Δ GD^ngu,min,z+1Obtaining the maximum value vector Delta GD of the variation vector of the gas consumption of the gas power station^ngu,maxVector of maximum values Δ GD^ngu,maxIs noted as Δ GD^ngu,max,z+1；

(5-10) the power grid judges the operation result, if the operation result meets the following conditions:

and is

Performing the step (6);

if it satisfies

Or satisfy

Then the lagrange multiplier is updated as follows:

electric network to GD^ngu,p,z+1、ΔGD^ngu,min,z+1、ΔGD^ngu,max,z+1、λ^z+1Transferring to a natural gas network, and enabling z to be z +1, and returning to the step (5-6);

(6) solving the distributed optimized dispatching model of the electric coupling system based on the consideration of uncertainty transfer in the step (5) to obtain the values of the variables to be solved of the power grid and the natural gas grid, namely the active power of the gas power station h in the power grid

H gas consumption of gas power station

Active power P of non-gas power station i_i ^genActive power WS for abandoning wind turbine generator j_jElectric power LS of electric load k_kActive power pf of branch between node l and node m_lmPhase angle theta of voltage at node l_lGas outlet flow G of natural gas well s in natural gas network_sH gas consumption of gas power station

Natural gas flow gf of pipeline between node u and node r_urAnd pressure ω of node r_rAnd taking the value of the variable to be solved as a parameter for distributed optimized operation of the electrical coupling system, and realizing distributed optimized scheduling of the electrical coupling system considering uncertainty transfer.