CN106339794A - Electric-thermal coupling multi-energy flow network node energy price calculation method - Google Patents

Abstract

The invention relates to an electric-thermal coupling multi-energy flow network node energy price calculation method, and belongs to the field of the energy market of a comprehensive energy system. According to the method, the coupling relationship of electricity price and heat price in an electric-thermal coupling network is considered, the present situation that joint cost pricing of a thermal-electric joint supply set can only be calculated in a decoupling way can be broken through, and the market interaction mechanism of the electric-thermal coupling system is restored so that the accuracy of marginal cost calculation can be enhanced in comparison with the existing method of respective pricing of the power supply and heat supply systems. Meanwhile, node energy price acts as one of real-time price and compensates the blank of real-time heat price in the heat supply network so that a certain thought is provided for congestion management and network loss allocation of the heat supply network.

Description

A kind of electric-thermal couples multipotency flow network node energy valency computational methods

Technical field

The present invention relates to a kind of electro thermal coupling multipotency flow network node energy valency computational methods, belong in integrated energy system Energy market field.

Background technology

Increasingly serious with energy supply, present in renewable energy utilization process, problem highlights, in order to optimize energy Source utilization ratio, building multipotency streaming system has become new means.In multipotency streaming system, the various forms of energy, such as electricity, Gas, hot and cold etc., producing, transmission, the links such as consumption realize coupling and interaction, improve the performance driving economy of whole system. In typical multipotency flow network, with the fastest developing speed with cogeneration system.Multipotency drift net currently for electric-thermal coupling The joint cost pricing problem of network has had a series of achievement in research, but research seeks to producing electricity price and caloric value Side decouples, then participates in bidding to respective energy market, and the process of this decoupling has simplification and error, deviates actual conditions.For Coupled relation between the accurate description multipotency stream market price, needs to explore a kind of pricing machine considering electric-thermal coupling effect System.

With the popularization of intelligent meter equipment, the user's excitation based on real time price and dsm possess enforcement base Plinth.Real time price improves the efficiency of load side demand management, improves the transparency of pricing mechanism, enhances use The property of participation at family and enthusiasm.Although the development level of heat metering is no more than electric-power metering, it is subject to energy saving, reduce temperature The continuous of the targets such as room gas, are increasingly received publicity based on the real-time caloric value of calorimeter.Heat supply market in world wide Develop leading country, Switzerland and Finland etc., when having occurred in that heating network runs, temperature and pressure reaches in service requirement The situation of limit, that is, there occurs choking phenomenon.As the node marginal pricing that can characterize Spatial Dimension in real time price, Ke Yiwei Heating network carries out congestion management provides certain thinking.

In order to consider electric-thermal coupling effect, accurately reflect user uses energy demand, proposes a kind of new pricing mechanism Node energy valency.Node energy valency is the popularization of node electricity price, its intension that is abundant and extending node electricity price, can preferably adapt to Needs in following comprehensive energy network Development.

Content of the invention

The purpose of the present invention is to propose to a kind of electric-thermal coupling multipotency flow network node energy valency computational methods are it is considered to electric-thermal system Electricity price in system, the influencing each other of caloric value, are that the foundation of electric-thermal System Market mechanism and improving lays the first stone.

Electro thermal coupling multipotency flow network node energy valency computational methods proposed by the present invention, comprise the following steps:

(1) set up the object function that an electric-thermal couples multipotency flow network optimal load flow, as follows:

$\min \underset{i\&element;n_{chp}}{\mathrm{\σ}}({\mathrm{\α}}_{5,i}{p}_{gi}{q}_{gi}+{\mathrm{\α}}_{4,i}{q}_{gi}^2+{\mathrm{\α}}_{3,i}{p}_{gi}^2+{\mathrm{\α}}_{2,i}{q}_{gi}+{\mathrm{\α}}_{1,i}{p}_{gi}+{\mathrm{\α}}_{0,i})$

Wherein: i is the numbering of thermoelectricity unit, n_chpFor the quantity of thermoelectricity unit, α_0,i～α_5,iFor i-th thermoelectricity unit Fuel cost coefficient, α_0,i～α_5,iObtain from the product description of i-th thermoelectricity unit, p_gi、q_giIt is respectively i-th thermoelectricity The electricity of unit is exerted oneself and heat is exerted oneself；

(2) set up the constraints of electro thermal coupling multipotency flow network optimal load flow, comprising:

The feasible region constraint of cogeneration unit in (2 1) electro thermal coupling multipotency flow network, as follows:

$a_l{p}_{gi}+b_l{q}_{gi}\&greaterequal;c_l,\∀l\&element;n_l$

Wherein: l is the line segment numbering of composition cogeneration unit feasible zone, n_lFor line segment sum, a_l,b_l,c_lFor each bar line Section parameter, a_l,b_l,c_lObtain from the product description of cogeneration unit,Represent arbitrarily,Represent for arbitrary l all There is above-mentioned linear restriction；

In (2 2) electric-thermal coupling multipotency flow network, electrical network adopts DC flow model, and its active power balance constrains such as Under:

$p_i=\underset{j\&element;i}{\mathrm{\σ}}{b}_{ij}{\mathrm{\θ}}_j,(i=1,2,3,...,n)$

Wherein: p_iFor the injection active power of electrical network interior joint i, θ_jFor the voltage phase angle of node j, b_ijFor node admittance square Battle array y i-th row, the imaginary part of jth column element, grid nodes admittance matrix y couples the EMS of multipotency streaming system from electric-thermal Middle acquisition；

The through-put power that (2 3) electric-thermal couples electrical network in multipotency flow network constrains:

p_ij=-b_ij(θ_i-θ_j)≤p_ijmax

Wherein: p_ijFor the line transmission active power of electrical network interior joint i to node j, p_ijmaxFor electrical network interior joint i to section The line transmission active power upper limit of point j；

(2-4) in electric-thermal coupling multipotency flow network, the thermodynamic equilibrium equation of heat supply network constrains:

$q=a_{a1}{g}_c{a}_{a1}^{t}t-a_{a2}{g}_c{t}_c$

g₀=a_a·g_c

$t_c=(a_{a1}^{t}t-t_a)\×(1-\mathrm{\λ}l/c_p\stackrel{\·}{m})+t_a$

Wherein: q is the injection thermal power column vector of heat supply network node, g_cFor heat supply network pipeline flow heat equivalent diagonal matrix, g₀For heat Net node flow heat equivalent column vector, t flows out temperature column vector, t for heat supply network node_cFor heat supply network pipeline section outlet temperature column vector, t_a The environment temperature being located for heat supply network pipeline section, c_pFor the specific heat capacity of heat supply network WATER AS FLOW MEDIUM, value is 4182 joules/(kilogram Celsius Degree), l is heat supply network length of pipe section, and λ is the thermal conductivity factor column vector of heat supply network pipeline section unit length, l and λ is from electro thermal coupling multipotency stream Obtain in the EMS of system,For the WATER AS FLOW MEDIUM mass flow of heat supply network pipeline section, a_a1Starting point augmentation for heat supply network pipeline section Incidence matrix, a_a2For the terminal augmented incidence matrix of heat supply network pipeline section, a_aFor the augmented incidence matrix of heat supply network pipeline section, a_a1、a_a2And a_a Obtain from the EMS of electro thermal coupling multipotency streaming system, t is matrix transposition；

(2-5) in electric-thermal coupling multipotency flow network, the operation of heat supply network constrains:

$\left\{\begin{array}{c}{g}_{cmin}\≤g_c\≤g_{cmax}\\ t_{cmin}\≤t_c\≤t_{cmax}\\ t_{min}\≤t\≤t_{\max }\end{array}\right.$

Wherein: g_cminAnd g_cmaxFor the operation bound of heat supply network pipeline flow heat equivalent, t_cminAnd t_cmaxGo out for heat supply network pipeline section The bound of mouth temperature, t_minAnd t_maxFlow out the bound of temperature, t for heat supply network node_cmin、t_cmax、t_minAnd t_maxFrom electric heating coupling Close in the EMS of multipotency streaming system and obtain.

(3) utilize method of Lagrange multipliers that the constraints in the object function in step (1) and step (2) is constituted one Individual Lagrangian, solves this Lagrangian using prim al- dual interior point m ethod, then power system node power equilibrium equation It is respectively node electricity price and the node of each node with the Lagrange multiplier corresponding to therrmodynamic system node power equilibrium equation Caloric value.

Electro thermal coupling multipotency flow network node energy valency computational methods proposed by the present invention, its feature and advantage are: this method Consider the coupled relation of electricity price and caloric value in electro thermal coupling network, having broken cogeneration unit joint cost price can only solve The present situation that coupling calculates, reduces the market interaction mechanism in electro thermal coupling system, compared to existing for power supply, heating system The method fixed a price respectively, improves the accuracy of marginal cost calculating.Meanwhile, node energy valency is as in real time price Kind, compensate for the blank of real-time caloric value in heat supply network, be that heating network carries out congestion management and Loss Allocation provides certain think of Road.

Specific embodiment

Electro thermal coupling multipotency flow network node energy valency computational methods proposed by the present invention, comprise the following steps:

(1) set up the object function that an electric-thermal couples multipotency flow network optimal load flow, as follows:

$\min \underset{i\&element;n_{chp}}{\mathrm{\σ}}({\mathrm{\α}}_{5,i}{p}_{gi}{q}_{gi}+{\mathrm{\α}}_{4,i}{q}_{gi}^2+{\mathrm{\α}}_{3,i}{p}_{gi}^2+{\mathrm{\α}}_{2,i}{q}_{gi}+{\mathrm{\α}}_{1,i}{p}_{gi}+{\mathrm{\α}}_{0,i})$

Wherein: i is the numbering of thermoelectricity unit, n_chpFor the quantity of thermoelectricity unit, α_0,i～α_5,iFor i-th thermoelectricity unit Fuel cost coefficient, α_0,i～α_5,iObtain from the product description of i-th thermoelectricity unit, p_gi、q_giIt is respectively i-th thermoelectricity The electricity of unit is exerted oneself and heat is exerted oneself；

(2) set up the constraints of electro thermal coupling multipotency flow network optimal load flow, comprising:

The feasible region constraint of cogeneration unit in (2 1) electro thermal coupling multipotency flow network, as follows:

$a_l{p}_{gi}+b_l{q}_{gi}\&greaterequal;c_l,\∀l\&element;n_l$

Wherein: l is the line segment numbering of composition cogeneration unit feasible zone, n_lFor line segment sum, a_l,b_l,c_lFor each bar line Section parameter, a_l,b_l,c_lObtain from the product description of cogeneration unit,Represent arbitrarily,Represent for arbitrary l all There is above-mentioned linear restriction；

In (2 2) electric-thermal coupling multipotency flow network, electrical network adopts DC flow model, and its active power balance constrains such as Under:

$p_i=\underset{j\&element;i}{\mathrm{\σ}}{b}_{ij}{\mathrm{\θ}}_j,(i=1,2,3,...,n)$

Wherein: p_iFor the injection active power of electrical network interior joint i, θ_jFor the voltage phase angle of node j, b_ijFor node admittance square Battle array y i-th row, the imaginary part of jth column element, grid nodes admittance matrix y couples the EMS of multipotency streaming system from electric-thermal Middle acquisition；

The through-put power that (2 3) electric-thermal couples electrical network in multipotency flow network constrains:

p_ij=-b_ij(θ_i-θ_j)≤p_ijmax

Wherein: p_ijFor the line transmission active power of electrical network interior joint i to node j, p_ijmaxFor electrical network interior joint i to section The line transmission active power upper limit of point j；

(2-4) in electric-thermal coupling multipotency flow network, the thermodynamic equilibrium equation of heat supply network constrains:

$q=a_{a1}{g}_c{a}_{a1}^{t}t-a_{a2}{g}_c{t}_c$

g₀=a_a·g_c

$t_c=(a_{a1}^{t}t-t_a)\×(1-\mathrm{\λ}l/c_p\stackrel{\·}{m})+t_a$

Wherein: q is the injection thermal power column vector of heat supply network node, g_cFor heat supply network pipeline flow heat equivalent diagonal matrix, g₀For heat Net node flow heat equivalent column vector, t flows out temperature column vector, t for heat supply network node_cFor heat supply network pipeline section outlet temperature column vector, t_a The environment temperature being located for heat supply network pipeline section, c_pFor the specific heat capacity of heat supply network WATER AS FLOW MEDIUM, value is 4182 joules/(kilogram Celsius Degree), l is heat supply network length of pipe section, and λ is the thermal conductivity factor column vector of heat supply network pipeline section unit length, l and λ is from electro thermal coupling multipotency stream Obtain in the EMS of system,For the WATER AS FLOW MEDIUM mass flow of heat supply network pipeline section, a_a1Starting point augmentation for heat supply network pipeline section Incidence matrix, a_a2For the terminal augmented incidence matrix of heat supply network pipeline section, a_aFor the augmented incidence matrix of heat supply network pipeline section, a_a1、a_a2And a_a Obtain from the EMS of electro thermal coupling multipotency streaming system, t is matrix transposition；

(2-5) in electric-thermal coupling multipotency flow network, the operation of heat supply network constrains:

$\left\{\begin{array}{c}{g}_{cmin}\≤g_c\≤g_{cmax}\\ t_{cmin}\≤t_c\≤t_{cmax}\\ t_{min}\≤t\≤t_{\max }\end{array}\right.$

Wherein: g_cminAnd g_cmaxFor the operation bound of heat supply network pipeline flow heat equivalent, t_cminAnd t_cmaxGo out for heat supply network pipeline section The bound of mouth temperature, t_minAnd t_maxFlow out the bound of temperature, t for heat supply network node_cmin、t_cmax、t_minAnd t_maxFrom electric heating coupling Close in the EMS of multipotency streaming system and obtain.

(3) utilize method of Lagrange multipliers that the constraints in the object function in step (1) and step (2) is constituted one Individual Lagrangian, solves this Lagrangian using prim al- dual interior point m ethod, then power system node power equilibrium equation It is respectively node electricity price and the node of each node with the Lagrange multiplier corresponding to therrmodynamic system node power equilibrium equation Caloric value.

Claims (1)

1. a kind of electro thermal coupling multipotency flow network node energy valency computational methods are it is characterised in that the method comprises the following steps:

(1) set up the object function that an electric-thermal couples multipotency flow network optimal load flow, as follows:

$\min \underset{i\&element;n_{chp}}{\σ}({\mathrm{\α}}_{5,i}{p}_{gi}{q}_{gi}+{\mathrm{\α}}_{4,i}{q}_{gi}^2+{\mathrm{\α}}_{3,i}{p}_{gi}^2+{\mathrm{\α}}_{2,i}{q}_{gi}+{\mathrm{\α}}_{1,i}{p}_{gi}+{\mathrm{\α}}_{0,i})$

Wherein: i couples the numbering of thermoelectricity unit in multipotency flow network, n for electric-thermal_chpFor the quantity of thermoelectricity unit, α_0,i～α_5,i For the fuel cost coefficient of i-th thermoelectricity unit, α_0,i～α_5,iObtain from the product description of i-th thermoelectricity unit, p_gi、 q_giThe electricity of respectively i-th thermoelectricity unit is exerted oneself and heat is exerted oneself；

(2) set up the constraints of an electro thermal coupling multipotency flow network optimal load flow, comprising:

The feasible region constraint of cogeneration unit in (2 1) electro thermal coupling multipotency flow network, as follows:

$\begin{array}{cc}{a}_l{p}_{gi}+b_l{q}_{gi}\&greaterequal;c_l& \∀l\&element;n_l\end{array}$

Wherein: l is the line segment numbering of composition cogeneration unit feasible zone, n_lFor line segment sum, a_l,b_l,c_lFor each bar line segment ginseng Number, a_l,b_l,c_lObtain from the product description of cogeneration unit,Represent arbitrarily,Represent on arbitrary l is had The linear restriction stated；

In (2 2) electric-thermal coupling multipotency flow network, electrical network adopts DC flow model, and the active power in DC flow model is put down Weighing apparatus constraint is as follows:

$p_i=\underset{j\&element;i}{\mathrm{\σ}}{b}_{ij}{\mathrm{\θ}}_j,(i=1,2,3,...,n)$

Wherein: p_iFor the injection active power of electrical network interior joint i, θ_jFor the voltage phase angle of electrical network interior joint j, b_ijFor grid nodes I-th row of admittance matrix y, the imaginary part of jth column element, grid nodes admittance matrix y couples the energy of multipotency streaming system from electric-thermal Obtain in management system；

The through-put power that (2 3) electric-thermal couples electrical network in multipotency flow network constrains:

p_ij=-b_ij(θ_i-θ_j)≤p_ijmax

Wherein: p_ijFor the line transmission active power of electrical network interior joint i to node j, p_ijmaxFor electrical network interior joint i to node j's The line transmission active power upper limit；

(2-4) in electric-thermal coupling multipotency flow network, the thermodynamic equilibrium equation of heat supply network constrains:

$q=a_{a1}{g}_c{a}_{a1}^{t}t-a_{a2}{g}_c{t}_c$

g₀=a_a·g_c

$t_c=(a_{a1}^{t}t-t_a)\×(1-\mathrm{\λ}l/c_p\stackrel{\·}{m})+t_a$

Wherein: q is the injection thermal power column vector of heat supply network node, g_cFor heat supply network pipeline flow heat equivalent diagonal matrix, g₀For heat supply network section Point flow heat equivalent column vector, t flows out temperature column vector, t for heat supply network node_cFor heat supply network pipeline section outlet temperature column vector, t_aFor heat The environment temperature that webmaster section is located, c_pFor the specific heat capacity of heat supply network WATER AS FLOW MEDIUM, value is 4182 joules/(kilogram degree Celsius), l is Heat supply network length of pipe section, λ is the thermal conductivity factor column vector of heat supply network pipeline section unit length, l and λ is from electro thermal coupling multipotency streaming system Obtain in EMS,For the WATER AS FLOW MEDIUM mass flow of heat supply network pipeline section, a_a1Starting point augmentation for heat supply network pipeline section associates square Battle array, a_a2For the terminal augmented incidence matrix of heat supply network pipeline section, a_aFor the augmented incidence matrix of heat supply network pipeline section, a_a1、a_a2And a_aFrom electric heating Obtain in the EMS of coupling multipotency streaming system, t is matrix transposition；

(2-5) in electric-thermal coupling multipotency flow network, the operation of heat supply network constrains:

$\left\{\begin{array}{c}{g}_{c\min }\≤g_c\≤g_{c\max }\\ t_{c\min }\≤t_c\≤t_{c\max }\\ t_{\min }\≤t\≤t_{\max }\end{array}\right.$

Wherein: g_cminAnd g_cmaxIt is respectively the operation upper and lower bound of heat supply network pipeline flow heat equivalent, t_cminAnd t_cmaxFor heat respectively The upper and lower bound of webmaster section outlet temperature, t_minAnd t_maxIt is respectively the upper and lower bound that heat supply network node flows out temperature, t_cmin、 t_cmax、t_minAnd t_maxObtain from the EMS of electro thermal coupling multipotency streaming system respectively；

(3) utilize method of Lagrange multipliers, the constraints in the object function in step (1) and step (2) is constituted one Lagrangian, solves this Lagrangian using prim al- dual interior point m ethod, obtains power system node power equilibrium equation With the Lagrange multiplier corresponding to therrmodynamic system node power equilibrium equation, that is, respectively electro thermal coupling multipotency flow network is each The node electricity price of node and node caloric value.