CN110752605B - Optimal power flow calculation method for electric-thermal coupling comprehensive energy system - Google Patents

Optimal power flow calculation method for electric-thermal coupling comprehensive energy system Download PDF

Info

Publication number
CN110752605B
CN110752605B CN201911040179.2A CN201911040179A CN110752605B CN 110752605 B CN110752605 B CN 110752605B CN 201911040179 A CN201911040179 A CN 201911040179A CN 110752605 B CN110752605 B CN 110752605B
Authority
CN
China
Prior art keywords
power
node
representing
pipeline
power flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911040179.2A
Other languages
Chinese (zh)
Other versions
CN110752605A (en
Inventor
王罡
侯磊
杨超
张靖
苏琨
夏天
张正文
姚大伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xiongan New Area Power Supply Company State Grid Hebei Electric Power Co
Tsinghua University
State Grid Corp of China SGCC
Original Assignee
Xiongan New Area Power Supply Company State Grid Hebei Electric Power Co
Tsinghua University
State Grid Corp of China SGCC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xiongan New Area Power Supply Company State Grid Hebei Electric Power Co, Tsinghua University, State Grid Corp of China SGCC filed Critical Xiongan New Area Power Supply Company State Grid Hebei Electric Power Co
Priority to CN201911040179.2A priority Critical patent/CN110752605B/en
Publication of CN110752605A publication Critical patent/CN110752605A/en
Application granted granted Critical
Publication of CN110752605B publication Critical patent/CN110752605B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to an optimal power flow calculation method of an electrothermal coupling comprehensive energy system based on a generalized inverse method, and belongs to the technical field of operation regulation and control of electric power/comprehensive energy systems. The method comprises the steps of firstly establishing a steady-state model of the regional electric-thermal coupling comprehensive energy system, then establishing a multi-target optimal power flow model of the regional comprehensive energy system based on a generalized inverse matrix theory, wherein multiple targets comprise energy use cost minimization and power distribution system network loss minimization, solving is carried out iteratively based on a Newton-Raphson method, and information transmission between an electric power system and a thermodynamic system is realized through conversion of running power of coupling equipment. The method provided by the invention can accurately reflect the operation condition of the regional comprehensive energy system, has good numerical stability and convergence, the optimization process is a feasible solution optimization process, each step can obtain a more optimal feasible solution, the calculation can be stopped at any time according to the actual condition, and the method is suitable for the real-time system optimization scheduling application.

Description

Optimal power flow calculation method for electric-thermal coupling comprehensive energy system
Technical Field
The invention relates to an optimal power flow calculation method of an electrothermal coupling comprehensive energy system based on a generalized inverse method, and belongs to the technical field of operation regulation and control of electric energy systems.
Background
The energy internet is a deep integration and development of a new generation of energy system and internet technology, and is also the focus and innovation frontier of current academic and industrial circles at home and abroad. On the basis of the energy internet, the national grid company also provides a ubiquitous power internet of things, and defines that a power system is used as a core and a backbone, various distributed resources are gathered together to form a regional comprehensive energy system, and interaction is established with a superior power grid. The research and construction of regional comprehensive energy systems is a key factor for the implementation of the landing of energy Internet in China, and has important significance.
In the field of optimal power flow calculation of regional integrated energy systems, the optimal power flow calculation of a power system is mature at present, but the research on multi-energy flow optimal power flow, particularly the optimal power flow research of a regional electric heating coupling system, is shallow. The defects of the prior art are mainly that the equipment modeling of the thermodynamic system is simple, and the situation of flow distribution of a thermal pipeline network system is not considered.
Disclosure of Invention
The invention aims to provide an optimal power flow calculation method of an electrothermal coupling comprehensive energy system based on a generalized inverse method, and overcomes the defects of the prior art. The method comprises the steps of firstly establishing a steady-state model of the regional electric-thermal coupling comprehensive energy system, then establishing a multi-target optimal power flow model of the regional comprehensive energy system based on a generalized inverse matrix theory, wherein multiple targets comprise energy use cost minimization and power distribution system network loss minimization, solving is carried out iteratively based on a Newton-Raphson method, and information transmission between an electric power system and a thermodynamic system is realized through conversion of running power of coupling equipment.
Specifically, the invention provides an optimal power flow calculation method of an electrothermal coupling comprehensive energy system based on a generalized inverse method, which comprises the following steps:
(1) establishing a steady-state power flow model of the electrothermal coupling comprehensive energy system, comprising the following steps of:
(1-1) establishing a steady-state power flow model of the power system;
(1-2) establishing a steady-state power flow model of the thermodynamic system;
(1-3) establishing a steady-state power flow model of the electric-thermal coupling equipment according to the power conversion relation between the electricity and the heat energy flow;
(2) establishing an optimal power flow model of the thermodynamic system, comprising the following steps:
(2-1) determining control variables and variable ranges thereof of the thermodynamic system;
(2-2) determining the control variable of the thermodynamic system together according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable;
(2-3) establishing an optimal power flow objective function of the thermodynamic system, wherein the objective is that the energy cost of the system is minimized;
(2-4) establishing an augmented constraint power flow equation;
(3) establishing an optimal power flow model of the power system, comprising:
(3-1) determining control variables of the power system and variable ranges thereof;
(3-2) jointly determining the control variable of the power system according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable;
(3-3) establishing an optimal power flow objective function of the power system, wherein the objective is that the power loss of a system network is minimized;
(3-4) establishing an augmented constraint power flow equation;
(4) solving the optimal power flow model of the thermodynamic system in the step (2) based on a generalized inverse method;
(5) solving a steady-state load flow model of the coupling equipment according to the calculation result in the step (4);
(6) and (4) solving the optimal power flow model of the power system in the step (3) based on the generalized inverse method.
Further, the air conditioner is provided with a fan,
(1) establishing a steady-state power flow model of the electrothermal coupling comprehensive energy system, which specifically comprises the following steps:
(1-1) establishing a steady-state power flow model of the power system, wherein the expression is as follows:
Figure BDA0002252621580000021
in the above formula, Pi、QiRespectively representing active and reactive power injected at node i, Gij、BijRespectively representing the conductance and susceptance between node i and node j, thetaijRepresenting the phase angle difference, V, between node i and node ji、VjRespectively representing the voltage amplitudes of the nodes i and j;
(1-2) establishing a steady-state power flow model of the thermodynamic system, wherein the expression is as follows:
Figure BDA0002252621580000022
in the above formula, A represents the node-pipe association matrix of the thermodynamic system, Au、AdRespectively representing an upper and a lower correlation matrix, BfRepresenting a loop matrix, M representing a pipeline working medium flow vector, MlRepresents the injection flow of the node, deltaH represents the pressure difference vector of the head end and the tail end of the pipeline, Z represents the height difference vector of the head end and the tail end of the pipeline, and T represents the pressure difference vector of the head end and the tail end of the pipelineeIndicating the temperature vector, T, of the working medium at the end of the pipenRepresenting the node temperature vector, PhRepresenting the node thermal power vector, TaRepresenting the external ambient temperature vector, HpRepresenting the pressure difference vector of the inlet and the outlet of the pipeline pump, HvRepresenting the pressure difference vector of the inlet and the outlet of the pipeline valve, F representing the friction coefficient vector of the pipeline, D representing the pipe diameter vector of the pipeline, S representing the cross-sectional area vector of the pipeline, L representing the length vector of the pipeline, CpRepresents a constant pressure specific heat capacity vector of the node,
Figure BDA0002252621580000023
represents the medium density vector of the pipeline working medium, E represents the diagonal matrix of the pipeline temperature attenuation coefficient, and has:
Figure BDA0002252621580000024
wherein λ isiWhich represents the thermal conductivity of the pipe i,
Figure BDA0002252621580000025
represents the medium value constant pressure specific heat capacity of the working medium of the pipeline iiDenotes the length, m, of the pipe iiThe flow of the pipeline i is represented, and the lower corner mark B represents B pipelines in total;
(1-3) establishing a steady-state power flow model of the electric-thermal coupling equipment according to the power conversion relation between the electric energy flow and the heat energy flow, wherein the expression is as follows:
Figure BDA0002252621580000031
wherein P represents the electrical power of the node, PpRepresenting electric power of the water pump, etapIndicates the efficiency of the water pump, hpThe pump head of the water pump is represented, k represents a power conversion coefficient, f represents an electric-thermal power conversion relation function of the coupling equipment, and m represents the flow of the pipeline;
(2) establishing an optimal power flow model of the thermodynamic system, which specifically comprises the following steps:
(2-1) determining the control variable and the variable range of the thermodynamic system, wherein the expression is as follows:
Figure BDA0002252621580000032
Figure BDA0002252621580000033
Figure BDA0002252621580000034
Figure BDA0002252621580000035
Figure BDA0002252621580000036
Figure BDA0002252621580000037
Figure BDA0002252621580000038
in the formula, mjDenotes the flow rate, T, of the pipe je.jRepresenting the temperature vector, xi, of the working medium at the j end of the pipelinejIs the opening of the valve on the jth pipe, ωjIs the shaft speed of the pump on the jth pipe, pn.iIs the pressure at the ith node, Tn.iIndicating the temperature of the i-th node, Ph.iRepresenting the thermal power vector of the ith node, b is the total quantity of the pipeline, n is the total quantity of the nodes, v is the total quantity of the valves on the pipeline, p is the total quantity of the pumps on the pipeline, the lower limit value of the variable is represented by the variable symbol with lower lines, and the upper limit value of the variable is represented by the variable symbol with upper lines;
(2-2) according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable, rewriting the control variable of the thermodynamic system and the variable range expression thereof as follows:
Figure BDA0002252621580000039
Figure BDA00022526215800000310
Figure BDA00022526215800000311
Figure BDA00022526215800000312
Figure BDA00022526215800000313
Figure BDA00022526215800000314
Figure BDA00022526215800000315
wherein alpha isj、
Figure BDA00022526215800000317
γj、τi、σi、βi、δiAre respectively corresponding to mj、Tej、ξj、ωj、pni、Tni、Phi1/2 with the variable symbol of superscript (1) representing the sum of the upper and lower limits of the variable, 1/2 with the variable symbol of superscript (2) representing the difference between the upper and lower limits of the variable;
(2-3) establishing an optimal power flow objective function of the thermodynamic system, wherein the objective is that the cost of the energy consumption of the system is minimized, and the expression is as follows:
Figure BDA00022526215800000316
in the formula, b0i、b1i、b2iIs a thermodynamic system sectionZero, first, second fitting coefficients of a cost-power relationship quadratic polynomial function for point i, FhsIs the upper bound of the objective function, yhbIs an aid decision variable, Δ FhIs an objective function FhThe amount of deviation of (d);
(2-4) establishing an augmented constraint power flow equation, wherein the expression is as follows:
Figure BDA0002252621580000041
in the above formula,. DELTA.mliRepresents the node injection flow deviation, Δ h, of node ijRepresenting the pressure difference vector of the head and the tail of the pipeline j, delta Te.jIndicating the temperature deviation of the working medium at the j end of the pipeline, delta PhiIndicating the thermal power deviation, Z, of node ijDenotes the head and tail end height difference of the pipe j, djDenotes the pipe diameter vector, s, of the pipe jjDenotes the cross-sectional area, f, of the conduit jjRepresenting the friction coefficient vector of the pipe j,
Figure BDA0002252621580000045
represents the median density of j working medium in the pipeline, ljWhich represents the length of the pipe j,
Figure BDA0002252621580000046
expressing the medium value constant pressure specific heat capacity of the working medium of the pipeline j, ad.jiLower correlation matrix, a, representing pipe i and pipe ju.jiAn upper correlation matrix, c, representing pipe i and pipe jp.iRepresents the constant-pressure specific heat capacity of the node i, aijIndicating connection information of the ith node and the jth branch, ajiIndicating connection information of the jth branch to the ith node, prefRepresents the reference node pressure value, k0j、k1j、k2jFitting coefficients of 0, 1 and 2 in a quadratic polynomial representing the relation between pump pressure and flow ratej.sExpressing the working medium density, lambda, at the location of valve s installation in the jth pipejRepresenting the heat conductivity coefficient of the jth pipeline;
linearizing the formula to obtain a correction equation of an iterative process by using a Newton Raphson method, wherein the expression is as follows:
Figure BDA0002252621580000042
(3) establishing an optimal power flow model of the power system, which specifically comprises the following steps:
(3-1) determining control variables and variable ranges thereof of the power system, wherein the expression is as follows:
Figure BDA0002252621580000043
Figure BDA0002252621580000044
wherein n is the total number of nodes, the lower-limit value of the variable is represented by the underlined variable symbol, and the upper-limit value of the variable is represented by the underlined variable symbol;
(3-2) rewriting the above formula as:
Figure BDA0002252621580000051
phi and upsilon are decision variables corresponding to V, Q respectively, a variable symbol with an upper label (1) represents 1/2 of the sum of an upper limit value and a lower limit value of the variable, and a variable symbol with an upper label (2) represents 1/2 of the difference of the upper limit value and the lower limit value of the variable;
(3-3) establishing an optimal power flow objective function of the power system, wherein the objective is that the power loss of the system network is minimized, and the expression is as follows:
Figure BDA0002252621580000052
in which b denotes a balanced node, Ps,bIs the sum of the active power of the power supply of the balanced node b, Pd,bIs the sum of the active power of the load of the balancing nodes, VbAnd VjIdentifying the voltages at balanced node b and node j, respectively, FesIs the upper bound of the objective function, yebIs an aid decision variable, θbjRepresenting the phase angle difference, G, between node j and equilibrium node bbj、BbjDenotes the conductance and susceptance, Δ F, between the balanced node b and the node j, respectivelyeIs the deviation of the objective function Fe;
(3-4) establishing an augmented constraint power flow equation, wherein the expression is as follows:
Figure BDA0002252621580000053
in the above formula, Ps,iIs the sum of the active power of the power supply of node i, Pd,iIs the sum of the i load active power of the node, Qs,iIs the sum of the active power of the node i power supply, Qd,iIs the sum of the i load active power of the node, Gij、BijRespectively representing the conductance and susceptance between node i and node j, thetaijRepresenting the phase angle difference between node j and node i;
linearizing the formula to obtain a correction equation of an iterative process by using a Newton Raphson method, wherein the expression is as follows:
Figure BDA0002252621580000054
(4) based on the generalized inverse method, solving the optimal power flow model of the thermodynamic system in the step (2), specifically comprising:
(4-1) giving values of contraction step length delta s, contraction multiple ns and calculation error epsilon;
(4-2) solving the augmented constrained power flow equation in the step (2-4) according to the correction equation in the step (2-4);
(4-3) if the calculation is converged,
1) solving a group of solutions X by an augmented constrained flow equationhAnd Fhs，Xh＝[α,σ,γ,τ,
Figure BDA0002252621580000061
β,δ,yhb]T；
2) Storing the tidal current solution;
3) according to Fhs＝FhsΔ s shrinkage FhsAfter k successive convergence, let Δ s be Δ s × ns, Fhs＝Fhs-Δs；
4) Turning to the step (4-2), solving the augmented constraint power flow equation again;
(4-4) if the calculation does not converge,
1) if the initial load flow has no solution, the calculation is finished;
2) recovering the last power flow solution;
3) according to Fhs＝FhsΔ s shrinkage FhsWherein Δ s ═ Δ s/ns;
4) turning to the step (4-2), solving the augmented constraint power flow equation again;
(4-5) when the delta s is less than or equal to the epsilon, finishing the calculation;
(5) according to the heat power P of the coupling equipment in the calculation result of the step (4)hCalculating electric power P by using the steady-state load flow model equation of the coupling equipment in the step (1-3), and calculating the pump lift h of the pipeline water pump according to the calculation result in the step (4)pCalculating the electric power P of the water pump by using the steady-state load flow model equation of the coupling equipment in the step (1-3)p；
(6) Based on the generalized inverse method, solving the optimal power flow model of the power system in the step (3), specifically comprising:
(6-1) giving values of contraction step length delta s, contraction multiple ns and calculation error epsilon;
(6-2) solving the augmented constrained power flow equation in the step (3-4) according to the modified equation in the step (3-4);
(6-3) if the calculation is converged,
1) solving a group of solutions X by an augmented constrained flow equationeAnd Fes，Xe＝[θ,φ,υ,yeb]T；
2) Storing the tidal current solution;
3) according to Fes＝Fes- Δ s contracts Fes, and when converging for k consecutive times, let Δ s ═ Δ s × ns, Fes＝Fes-Δs；
4) Turning to the step (6-2), solving the augmented constraint power flow equation again;
(6-4) if the calculation does not converge,
1) if the initial load flow has no solution, the calculation is finished;
2) recovering the last power flow solution;
3) according to Fes＝FesΔ s shrinkage FesWherein Δ s ═ Δ s/ns;
4) turning to the step (6-2), solving the augmented constraint power flow equation again;
(6-5) when Δ s ≦ ε, the calculation is ended.
Further, the air conditioner is provided with a fan,
the step (4-2) specifically comprises the following steps:
1) given a maximum number of iterations kmaxAnd a set of initial values Xh＝[α,σ,γ,τ,
Figure BDA0002252621580000062
β,δ,yhb]TWherein y ishb0, and let k be 0;
2) according to XhThe solution of a and deltab is carried out,
Figure BDA0002252621580000071
Δb＝[Δml,Δh,ΔTe,ΔPh,ΔFh]T
3) solving the generalized inverse matrix A of A+，A+＝AT(AAT)-1；
4) Solving for Δ XhAnd update Xh，ΔXh＝A+Δb，Xh＝Xh-ΔXhAnd let k be k + 1;
5) if Δ XhEpsilon is less than or equal to or the number of iterations k is more than kmaxAnd finishing the calculation, otherwise, turning to the step 2).
Further, the step (6-2) specifically comprises:
1) given epsilon, maximum number of iterations kmaxAnd a set of initial values Xe＝[θ,φ,υ,yeb]TWherein y iseb0, and let k be 0;
2) according to XeThe solution of a and deltab is carried out,
Figure BDA0002252621580000072
Δb＝[ΔP,ΔQ,ΔFe]T
3) solving the generalized inverse matrix A of A+，A+＝AT(AAT)-1；
4) Solving for Δ XeAnd update Xe，ΔXe＝A+Δb，Xe＝Xe-ΔXeAnd let k be k + 1;
5) if Δ XeEpsilon is less than or equal to or the number of iterations k is more than kmaxAnd finishing the calculation, otherwise, turning to the step 2).
Further, the air conditioner is provided with a fan,
the heat transfer working medium of the pipe network in the thermodynamic system is liquid, and the pipeline i
Figure BDA0002252621580000073
And
Figure BDA0002252621580000074
is a constant.
Further, the air conditioner is provided with a fan,
the electric-thermal power conversion relation function f of the coupling equipment is in the form of a quadratic polynomial:
Ph＝khe2*Pe*Pe+khe1*Pe+khe0，khe2、khe1、khe0as fitting coefficient, PhIndicating the thermal power of the device, PeRepresenting the plant electrical power.
Further, the air conditioner is provided with a fan,
the electric-thermal power conversion relation function f of the coupling equipment is in the form of a direct proportional relation.
The method provided by the invention can accurately reflect the operation condition of the regional comprehensive energy system, and considers the situation of the trend distribution of the thermal power pipeline system. The method has good numerical stability and convergence, the optimization process is a feasible solution optimization process, a more optimal feasible solution can be obtained in each step, the calculation can be stopped at any time according to actual conditions, and the method is suitable for real-time system optimization scheduling application.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a flow chart of an optimal power flow calculation method for an electrothermal coupling integrated energy system according to an embodiment of the invention;
fig. 2 shows a flow chart of solving an augmented constrained power flow equation according to a step-modified equation according to an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides an optimal power flow calculation method of an electrothermal coupling comprehensive energy system based on a generalized inverse method, which comprises the following steps of:
(1) establishing a steady-state power flow model of the electrothermal coupling comprehensive energy system, comprising the following steps of:
(1-1) establishing a steady-state power flow model of the power system;
(1-2) establishing a steady-state power flow model of the thermodynamic system;
the load flow calculation method considers the load flow distribution of the thermal power pipe network system, more accords with the actual application scene of the comprehensive energy system, and has more accurate and comprehensive calculation results.
(1-3) establishing a steady-state power flow model of the electric-thermal coupling equipment according to the power conversion relation between the electricity and the heat energy flow;
(2) establishing an optimal power flow model of the thermodynamic system, comprising the following steps:
(2-1) determining control variables and variable ranges thereof of the thermodynamic system;
(2-2) determining the control variable of the thermodynamic system together according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable;
(2-3) establishing an optimal power flow objective function of the thermodynamic system, wherein the objective is that the energy cost of the system is minimized;
(2-4) establishing an augmented constraint power flow equation;
(3) establishing an optimal power flow model of the power system, comprising:
(3-1) determining control variables of the power system and variable ranges thereof;
(3-2) jointly determining the control variable of the power system according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable;
(3-3) establishing an optimal power flow objective function of the power system, wherein the objective is that the power loss of a system network is minimized;
(3-4) establishing an augmented constraint power flow equation;
(4) solving the optimal power flow model of the thermodynamic system in the step (2) based on a generalized inverse method;
(5) solving a steady-state load flow model of the coupling equipment according to the calculation result in the step (4);
(6) and (4) solving the optimal power flow model of the power system in the step (3) based on the generalized inverse method.
The invention does not limit the sequence of the steps, and the steps without parameter reference relation can be exchanged in sequence or executed at the same time. The above calculation steps are further detailed below:
(1) establishing a steady-state power flow model of the electrothermal coupling comprehensive energy system, which specifically comprises the following steps:
(1-1) establishing a steady-state power flow model of the power system, wherein the expression is as follows:
Figure BDA0002252621580000091
in the above formula, Pi、QiRespectively representing active and reactive power injected at node i, Gij、BijRespectively representing the conductance and susceptance between node i and node j, thetaijRepresenting the phase angle difference, V, between node i and node ji、VjRespectively representing the voltage amplitudes of the nodes i and j;
(1-2) establishing a steady-state power flow model of the thermodynamic system, wherein the expression is as follows:
Figure BDA0002252621580000092
in the above formula, A represents the node-pipe association matrix of the thermodynamic system, Au、AdRespectively representing an upper and a lower correlation matrix, BfRepresenting a loop matrix, M representing a pipeline working medium flow vector, MlRepresents the injection flow of the node, deltaH represents the pressure difference vector of the head end and the tail end of the pipeline, Z represents the height difference vector of the head end and the tail end of the pipeline, and T represents the pressure difference vector of the head end and the tail end of the pipelineeIndicating the temperature vector, T, of the working medium at the end of the pipenRepresenting the node temperature vector, PhRepresenting the node thermal power vector, TaRepresenting the external ambient temperature vector, HpRepresenting the pressure difference vector of the inlet and the outlet of the pipeline pump, HvExpressing the pressure difference vector of the inlet and the outlet of the pipeline valve, F expressing the friction coefficient vector of the pipeline, D expressing the pipe diameter vector of the pipeline, S expressing the cross-sectional area vector of the pipeline, and L expressing the length direction of the pipelineAmount, CpRepresents a constant pressure specific heat capacity vector of the node,
Figure BDA0002252621580000093
represents the medium density vector of the pipeline working medium, E represents the diagonal matrix of the pipeline temperature attenuation coefficient, and has:
Figure BDA0002252621580000094
wherein λ isiWhich represents the thermal conductivity of the pipe i,
Figure BDA0002252621580000095
represents the medium value constant pressure specific heat capacity of the working medium of the pipeline iiDenotes the length, m, of the pipe iiThe flow of the pipeline i is represented, and the lower corner mark B represents B pipelines in total;
(1-3) establishing a steady-state power flow model of the electric-thermal coupling equipment according to the power conversion relation between the electric energy flow and the heat energy flow, wherein the expression is as follows:
Figure BDA0002252621580000101
wherein P represents the electrical power of the node, PpRepresenting electric power of the water pump, etapIndicates the efficiency of the water pump, hpThe pump head of the water pump is represented, k represents a power conversion coefficient, f represents an electric-thermal power conversion relation function of the coupling equipment, and m represents the flow of the pipeline;
(2) establishing an optimal power flow model of the thermodynamic system, which specifically comprises the following steps:
(2-1) determining the control variable and the variable range of the thermodynamic system, wherein the expression is as follows:
Figure BDA0002252621580000102
Figure BDA0002252621580000103
Figure BDA0002252621580000104
Figure BDA0002252621580000105
Figure BDA0002252621580000106
Figure BDA0002252621580000107
Figure BDA0002252621580000108
in the formula, mjDenotes the flow rate, T, of the pipe je.jRepresenting the temperature vector, xi, of the working medium at the j end of the pipelinejIs the opening of the valve on the jth pipe, ωjIs the shaft speed of the pump on the jth pipe, pn.iIs the pressure at the ith node, Tn.iIndicating the temperature of the i-th node, Ph.iRepresenting the thermal power vector of the ith node, b is the total quantity of the pipeline, n is the total quantity of the nodes, v is the total quantity of the valves on the pipeline, p is the total quantity of the pumps on the pipeline, the lower limit value of the variable is represented by the variable symbol with lower lines, and the upper limit value of the variable is represented by the variable symbol with upper lines;
(2-2) according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable, rewriting the control variable of the thermodynamic system and the variable range expression thereof as follows:
Figure BDA0002252621580000109
Figure BDA00022526215800001010
Figure BDA00022526215800001011
Figure BDA00022526215800001012
Figure BDA00022526215800001013
Figure BDA00022526215800001014
Figure BDA00022526215800001015
wherein alpha isj、
Figure BDA00022526215800001016
γj、τi、σi、βi、δiAre respectively corresponding to mj、Tej、ξj、ωj、pni、Tni、Phi1/2 with the variable symbol of superscript (1) representing the sum of the upper and lower limits of the variable, 1/2 with the variable symbol of superscript (2) representing the difference between the upper and lower limits of the variable;
(2-3) establishing an optimal power flow objective function of the thermodynamic system, wherein the objective is that the cost of the energy consumption of the system is minimized, and the expression is as follows:
Figure BDA00022526215800001017
in the formula, b0i、b1i、b2iIs the zero, first and second fitting coefficient of the cost-power relation quadratic polynomial function of the thermodynamic system node i, FhsIs the upper bound of the objective function, yhbIs an aid decision variable, Δ FhIs an objective function FhThe amount of deviation of (d); Δ FhIndicating the deviation that occurs during the solution, in general if the calculation converges, Δ Fh should be 0.
(2-4) establishing an augmented constraint power flow equation, wherein the expression is as follows:
Figure BDA0002252621580000111
in the above formula,. DELTA.mliRepresents the node injection flow deviation, Δ h, of node ijRepresenting the pressure difference vector of the head and the tail of the pipeline j, delta Te.jIndicating the temperature deviation of the working medium at the j end of the pipeline, delta PhiIndicating the thermal power deviation, Z, of node ijDenotes the head and tail end height difference of the pipe j, djDenotes the pipe diameter vector, s, of the pipe jjDenotes the cross-sectional area, f, of the conduit jjRepresenting the friction coefficient vector of the pipe j,
Figure BDA0002252621580000115
represents the median density of j working medium in the pipeline, ljWhich represents the length of the pipe j,
Figure BDA0002252621580000116
expressing the medium value constant pressure specific heat capacity of the working medium of the pipeline j, ad.jiLower correlation matrix, a, representing pipe i and pipe ju.jiAn upper correlation matrix, c, representing pipe i and pipe jp.iRepresents the constant-pressure specific heat capacity of the node i, aijIndicating connection information of the ith node and the jth branch, ajiIndicating connection information of the jth branch to the ith node, prefRepresents the reference node pressure value, k0j、k1j、k2jFitting coefficients of 0, 1 and 2 in a quadratic polynomial representing the relation between pump pressure and flow ratej.sExpressing the working medium density, lambda, at the location of valve s installation in the jth pipejRepresenting the heat conductivity coefficient of the jth pipeline;
linearizing the formula to obtain a correction equation of an iterative process by using a Newton Raphson method, wherein the expression is as follows:
Figure BDA0002252621580000112
(3) establishing an optimal power flow model of the power system, which specifically comprises the following steps:
(3-1) determining control variables and variable ranges thereof of the power system, wherein the expression is as follows:
Figure BDA0002252621580000113
Figure BDA0002252621580000114
wherein n is the total number of nodes, the lower-limit value of the variable is represented by the underlined variable symbol, and the upper-limit value of the variable is represented by the underlined variable symbol;
(3-2) rewriting the above formula as:
Figure BDA0002252621580000121
phi and upsilon are decision variables corresponding to V, Q respectively, a variable symbol with an upper label (1) represents 1/2 of the sum of an upper limit value and a lower limit value of the variable, and a variable symbol with an upper label (2) represents 1/2 of the difference of the upper limit value and the lower limit value of the variable;
(3-3) establishing an optimal power flow objective function of the power system, wherein the objective is that the power loss of the system network is minimized, and the expression is as follows:
Figure BDA0002252621580000122
in which b denotes a balanced node, Ps,bIs the sum of the active power of the power supply of the balanced node b, Pd,bIs the sum of the active power of the load of the balancing nodes, VbAnd VjIdentifying the voltages at balanced node b and node j, respectively, FesIs the upper bound of the objective function, yebIs an aid decision variable, θbjRepresenting the phase angle difference, G, between node j and equilibrium node bbj、BbjDenotes the conductance and susceptance, Δ F, between the balanced node b and the node j, respectivelyeIs the deviation of the objective function Fe;
(3-4) establishing an augmented constraint power flow equation, wherein the expression is as follows:
Figure BDA0002252621580000123
in the above formula, Ps,iIs the sum of the active power of the power supply of node i, Pd,iIs the sum of the i load active power of the node, Qs,iIs the sum of the active power of the node i power supply, Qd,iIs the sum of the i load active power of the node, Gij、BijRespectively representing the conductance and susceptance between node i and node j, thetaijRepresenting the phase angle difference between node j and node i;
linearizing the formula to obtain a correction equation of an iterative process by using a Newton Raphson method, wherein the expression is as follows:
Figure BDA0002252621580000124
(4) based on the generalized inverse method, solving the optimal power flow model of the thermodynamic system in the step (2), specifically comprising:
(4-1) giving values of contraction step length delta s, contraction multiple ns and calculation error epsilon;
(4-2) solving the augmented constrained power flow equation in the step (2-4) according to the correction equation in the step (2-4);
(4-3) if the calculation is converged,
1) solving a group of solutions X by an augmented constrained flow equationhAnd Fhs，Xh＝[α,σ,γ,τ,
Figure BDA0002252621580000125
β,δ,yhb]T；
2) Storing the tidal current solution;
3) according to Fhs＝FhsΔ s shrinkage FhsAfter k successive convergence, let Δ s be Δ s × ns, Fhs＝Fhs-Δs；
4) Turning to the step (4-2), solving the augmented constraint power flow equation again;
(4-4) if the calculation does not converge,
1) if the initial load flow has no solution, the calculation is finished;
2) recovering the last power flow solution;
3) according to Fhs＝FhsΔ s shrinkage FhsWherein Δ s ═ Δ s/ns;
4) turning to the step (4-2), solving the augmented constraint power flow equation again;
(4-5) when the delta s is less than or equal to the epsilon, finishing the calculation;
and when the current solution is continuously converged, performing recalculation after expansion on the current solution, and if the current solution is not converged, performing contraction on the current solution so as to continuously approach the optimal solution meeting the convergence.
(5) According to the heat power P of the coupling equipment in the calculation result of the step (4)hCalculating electric power P by using the steady-state load flow model equation of the coupling equipment in the step (1-3), and calculating the pump lift h of the pipeline water pump according to the calculation result in the step (4)pCalculating the electric power P of the water pump by using the steady-state load flow model equation of the coupling equipment in the step (1-3)p；
(6) Based on the generalized inverse method, solving the optimal power flow model of the power system in the step (3), specifically comprising:
(6-1) giving values of contraction step length delta s, contraction multiple ns and calculation error epsilon;
(6-2) as shown in FIG. 2, solving the augmented constrained power flow equation in the step (3-4) according to the modified equation in the step (3-4);
(6-3) if the calculation is converged,
1) restraint tide by augmentationSolving a set of solutions X by a flow equationeAnd Fes，Xe＝[θ,φ,υ,yeb]T；
2) Storing the tidal current solution;
3) according to Fes＝Fes- Δ s contracts Fes, and when converging for k consecutive times, let Δ s ═ Δ s × ns, Fes＝Fes-Δs；
4) Turning to the step (6-2), solving the augmented constraint power flow equation again;
(6-4) if the calculation does not converge,
1) if the initial load flow has no solution, the calculation is finished;
2) recovering the last power flow solution;
3) according to Fes＝FesΔ s shrinkage FesWherein Δ s ═ Δ s/ns;
4) turning to the step (6-2), solving the augmented constraint power flow equation again;
(6-5) when Δ s ≦ ε, the calculation is ended.
Further, the step (4-2) specifically includes:
1) based on a given epsilon, a maximum number of iterations k is specifiedmaxAnd a set of initial values Xh＝[α,σ,γ,τ,
Figure BDA0002252621580000131
β,δ,yhb]TWherein y ishb0, and let k be 0;
2) according to XhThe solution of a and deltab is carried out,
Figure BDA0002252621580000141
Δb＝[Δml,Δh,ΔTe,ΔPh,ΔFh]T
3) solving the generalized inverse matrix A of A+，A+＝AT(AAT)-1；
4) Solving for Δ XhAnd update Xh，ΔXh＝A+Δb，Xh＝Xh-ΔXhAnd is combined withLet k be k + 1;
5) if Δ XhEpsilon is less than or equal to or the number of iterations k is more than kmaxAnd finishing the calculation, otherwise, turning to the step 2).
Further, the step (6-2) specifically comprises:
1) given epsilon, maximum number of iterations kmaxAnd a set of initial values Xe [ theta, phi, upsilon, yeb]TWherein y iseb0, and let k be 0;
2) according to XeThe solution of a and deltab is carried out,
Figure BDA0002252621580000142
Δb＝[ΔP,ΔQ,ΔFe]T
3) solving the generalized inverse matrix A of A+，A+＝AT(AAT)-1；
4) Solving for Δ XeAnd update Xe，ΔXe＝A+Δb，Xe＝Xe-ΔXeAnd let k be k + 1;
5) if Δ XeEpsilon is less than or equal to or the number of iterations k is more than kmaxAnd finishing the calculation, otherwise, turning to the step 2).
In the embodiment of the invention, the heat transfer working medium of the pipe network in the thermodynamic system is in a liquid state, and the pipeline i
Figure BDA0002252621580000144
And
Figure BDA0002252621580000143
is a constant.
In the embodiment of the invention, the electric-thermal power conversion relation function f of the coupling equipment is in the form of a quadratic polynomial:
Ph＝khe2*Pe*Pe+khe1*Pe+khe0，khe2、khe1、khe0as fitting coefficient, PhIndicating the thermal power of the device, PeRepresenting the plant electrical power.
In another embodiment, the coupling device electrical-thermal power conversion function f may also be in the form of a direct proportional relation.
The variables in the above model expressions are per unit values, and specifically, all variables are processed as dimensionless variables except for the radian system of the control variable related to the angle.
Although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (7)

1. An optimal power flow calculation method for an electrothermal coupling comprehensive energy system is characterized by comprising the following steps:
(1) establishing a steady-state power flow model of the electrothermal coupling comprehensive energy system, comprising the following steps of:
(1-1) establishing a steady-state power flow model of the power system;
(1-2) establishing a steady-state power flow model of the thermodynamic system;
(1-3) establishing a steady-state power flow model of the electric-thermal coupling equipment according to the power conversion relation between the electricity and the heat energy flow;
(2) establishing an optimal power flow model of the thermodynamic system, comprising the following steps:
(2-1) determining control variables and variable ranges thereof of the thermodynamic system;
(2-2) determining the control variable of the thermodynamic system together according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable;
(2-3) establishing an optimal power flow objective function of the thermodynamic system, wherein the objective is that the energy cost of the system is minimized;
(2-4) establishing an augmented constraint power flow equation;
(3) establishing an optimal power flow model of the power system, comprising:
(3-1) determining control variables of the power system and variable ranges thereof;
(3-2) jointly determining the control variable of the power system according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable;
(3-3) establishing an optimal power flow objective function of the power system, wherein the objective is that the power loss of a system network is minimized;
(3-4) establishing an augmented constraint power flow equation;
(4) solving the optimal power flow model of the thermodynamic system in the step (2) based on a generalized inverse method;
(5) solving a steady-state load flow model of the coupling equipment according to the calculation result in the step (4);
(6) and (4) solving the optimal power flow model of the power system in the step (3) based on the generalized inverse method.
2. The optimal power flow calculation method for the electrothermal coupling integrated energy system according to claim 1,
(1) establishing a steady-state power flow model of the electrothermal coupling comprehensive energy system, which specifically comprises the following steps:
(1-1) establishing a steady-state power flow model of the power system, wherein the expression is as follows:
Figure FDA0002883512920000011
in the above formula, Pi、QiRespectively representing active and reactive power injected at node i, Gij、BijRespectively representing the conductance and susceptance between node i and node j, thetaijRepresenting the phase angle difference, V, between node i and node ji、VjRespectively representing the voltage amplitudes of the nodes i and j;
(1-2) establishing a steady-state power flow model of the thermodynamic system, wherein the expression is as follows:
Figure FDA0002883512920000021
in the above formula, A representsNode-pipeline association matrix of thermodynamic system, Au、AdRespectively representing an upper and a lower correlation matrix, BfRepresenting a loop matrix, M representing a pipeline working medium flow vector, MlRepresents the injection flow of the node, deltaH represents the pressure difference vector of the head end and the tail end of the pipeline, Z represents the height difference vector of the head end and the tail end of the pipeline, and T represents the pressure difference vector of the head end and the tail end of the pipelineeIndicating the temperature vector, T, of the working medium at the end of the pipenRepresenting the node temperature vector, PhRepresenting the node thermal power vector, TaRepresenting the external ambient temperature vector, HpRepresenting the pressure difference vector of the inlet and the outlet of the pipeline pump, HvRepresenting the pressure difference vector of the inlet and the outlet of the pipeline valve, F representing the friction coefficient vector of the pipeline, D representing the pipe diameter vector of the pipeline, S representing the cross-sectional area vector of the pipeline, L representing the length vector of the pipeline, CpRepresents a constant pressure specific heat capacity vector of the node,
Figure FDA00028835129200000212
represents the medium density vector of the pipeline working medium, E represents the diagonal matrix of the pipeline temperature attenuation coefficient, and has:
Figure FDA0002883512920000022
wherein λ isiWhich represents the thermal conductivity of the pipe i,
Figure FDA0002883512920000023
represents the medium value constant pressure specific heat capacity of the working medium of the pipeline iiDenotes the length, m, of the pipe iiThe flow of the pipeline i is represented, and the lower corner mark B represents B pipelines in total;
(1-3) establishing a steady-state power flow model of the electric-thermal coupling equipment according to the power conversion relation between the electric energy flow and the heat energy flow, wherein the expression is as follows:
Figure FDA0002883512920000024
wherein P represents the electrical power of the node, PpRepresenting electric power of the water pump, etapIndicates the efficiency of the water pump, hpThe pump head of the water pump is represented, k represents a power conversion coefficient, f represents an electric-thermal power conversion relation function of the coupling equipment, and m represents the flow of the pipeline;
(2) establishing an optimal power flow model of the thermodynamic system, which specifically comprises the following steps:
(2-1) determining the control variable and the variable range of the thermodynamic system, wherein the expression is as follows:
Figure FDA0002883512920000025
Figure FDA0002883512920000026
Figure FDA0002883512920000027
Figure FDA0002883512920000028
Figure FDA0002883512920000029
Figure FDA00028835129200000210
Figure FDA00028835129200000211
in the formula, mjDenotes the flow rate, T, of the pipe je.jRepresenting the temperature vector, xi, of the working medium at the j end of the pipelinejIs the opening of the valve on the jth pipe,ωjis the shaft speed of the pump on the jth pipe, pn.iIs the pressure at the ith node, Tn.iIndicating the temperature of the i-th node, Ph.iRepresenting the thermal power vector of the ith node, b is the total quantity of the pipeline, n is the total quantity of the nodes, v is the total quantity of the valves on the pipeline, p is the total quantity of the pumps on the pipeline, the lower limit value of the variable is represented by the variable symbol with lower lines, and the upper limit value of the variable is represented by the variable symbol with upper lines;
(2-2) according to the upper limit value and the lower limit value of the control variable and the decision value of the control variable, rewriting the control variable of the thermodynamic system and the variable range expression thereof as follows:
mj＝mj (1)+mj (2)sinαj,(j＝1,...,b)
Te.j＝Te.j (1)+Te.j (2)sinφj,(j＝1,...,b)
ξj＝ξj (1)+ξj (2)sinγj,(j＝1,...,v)
ωj＝ωj (1)+ωj (2)sinτj,(j＝1,...,p)
pn.i＝pn.i (1)+pn.i (2)sinσi,(i＝1,...,n)
Tn.i＝Tn.i (1)+Tn.i (2)sinβi,(i＝1,...,n)
Ph.i＝Ph.i (1)+Ph.i (2)sinδi,(i＝1,...,n)
wherein alpha isj、
Figure FDA0002883512920000039
γj、τi、σi、βi、δiAre respectively corresponding to mj、Tej、ξj、ωj、pni、Tni、PhiWith the variable symbol of the superscript (1) representing the decision variable1/2 for the sum of the upper and lower limits of the variable, 1/2 for the difference between the upper and lower limits of the variable being indicated by the variable symbol with superscript (2);
(2-3) establishing an optimal power flow objective function of the thermodynamic system, wherein the objective is that the cost of the energy consumption of the system is minimized, and the expression is as follows:
Figure FDA0002883512920000031
in the formula, b0i、b1i、b2iIs the zero, first and second fitting coefficient of the cost-power relation quadratic polynomial function of the thermodynamic system node i, FhsIs the upper bound of the objective function, yhbIs an aid decision variable, Δ FhIs an objective function FhThe amount of deviation of (d);
(2-4) establishing an augmented constraint power flow equation, wherein the expression is as follows:
Figure FDA0002883512920000032
Figure FDA0002883512920000033
Figure FDA0002883512920000034
Figure FDA0002883512920000035
Figure FDA0002883512920000036
in the above formula,. DELTA.mliRepresents the node injection flow deviation, Δ h, of node ijRepresenting the pressure difference vector of the head and the tail of the pipeline j, delta Te.jIndicating the temperature deviation of the working medium at the j end of the pipeline, delta PhiIndicating the thermal power deviation, z, of node ijDenotes the head and tail end height difference of the pipe j, djDenotes the pipe diameter vector, s, of the pipe jjDenotes the cross-sectional area, f, of the conduit jjRepresenting the friction coefficient vector of the pipe j,
Figure FDA0002883512920000037
represents the median density of j working medium in the pipeline, ljWhich represents the length of the pipe j,
Figure FDA0002883512920000038
expressing the medium value constant pressure specific heat capacity of the working medium of the pipeline j, ad.jiLower correlation matrix, a, representing pipe i and pipe ju.jiAn upper correlation matrix, c, representing pipe i and pipe jp.iRepresents the constant-pressure specific heat capacity of the node i, aijIndicating connection information of the ith node and the jth branch, ajiIndicating connection information of the jth branch to the ith node, prefRepresents the reference node pressure value, k0j、k1j、k2jFitting coefficients of 0, 1 and 2 in a quadratic polynomial representing the relation between pump pressure and flow ratej.sExpressing the working medium density, lambda, at the location of valve s installation in the jth pipejRepresenting the heat conductivity coefficient of the jth pipeline;
linearizing the formula to obtain a correction equation of an iterative process by using a Newton Raphson method, wherein the expression is as follows:
Figure FDA0002883512920000041
(3) establishing an optimal power flow model of the power system, which specifically comprises the following steps:
(3-1) determining control variables and variable ranges thereof of the power system, wherein the expression is as follows:
Figure FDA0002883512920000042
Figure FDA0002883512920000043
wherein n is the total number of nodes, the lower-limit value of the variable is represented by the underlined variable symbol, and the upper-limit value of the variable is represented by the underlined variable symbol;
(3-2) rewriting the above formula as:
Figure FDA0002883512920000044
phi and upsilon are decision variables corresponding to V, Q respectively, a variable symbol with an upper label (1) represents 1/2 of the sum of an upper limit value and a lower limit value of the variable, and a variable symbol with an upper label (2) represents 1/2 of the difference of the upper limit value and the lower limit value of the variable;
(3-3) establishing an optimal power flow objective function of the power system, wherein the objective is that the power loss of the system network is minimized, and the expression is as follows:
Figure FDA0002883512920000045
in which b denotes a balanced node, Ps，bIs the sum of the active power of the power supply of the balanced node b, Pd，bIs the sum of the active power of the load of the balancing nodes, VbAnd VjIdentifying the voltages at balanced node b and node j, respectively, FesIs the upper bound of the objective function, yebIs an aid decision variable, θbjRepresenting the phase angle difference, G, between node j and equilibrium node bbj、BbjDenotes the conductance and susceptance, Δ F, between the balanced node b and the node j, respectivelyeIs the deviation of the objective function Fe;
(3-4) establishing an augmented constraint power flow equation, wherein the expression is as follows:
Figure FDA0002883512920000046
in the above formula, Ps，iIs the sum of the active power of the power supply of node i, Pd，iIs the sum of the i load active power of the node, Qs，iIs the sum of the active power of the node i power supply, Qd，iIs the sum of the i load active power of the node, Gij、BijRespectively representing the conductance and susceptance between node i and node j, thetaijRepresenting the phase angle difference between node j and node i;
linearizing the formula to obtain a correction equation of an iterative process by using a Newton Raphson method, wherein the expression is as follows:
Figure FDA0002883512920000051
(4) based on the generalized inverse method, solving the optimal power flow model of the thermodynamic system in the step (2), specifically comprising:
(4-1) giving values of contraction step length delta s, contraction multiple ns and calculation error epsilon;
(4-2) solving the augmented constrained power flow equation in the step (2-4) according to the correction equation in the step (2-4);
(4-3) if the calculation is converged,
1) solving a group of solutions X by an augmented constrained flow equationhAnd Fhs，
Figure FDA0002883512920000052
Figure FDA0002883512920000053
2) Storing the tidal current solution;
3) according to Fhs＝FhsΔ s shrinkage FhsAfter k successive convergence, let Δ s be Δ s × ns, Fhs＝Fhs-Δs；
4) Turning to the step (4-2), solving the augmented constraint power flow equation again;
(4-4) if the calculation does not converge,
1) if the initial load flow has no solution, the calculation is finished;
2) recovering the last power flow solution;
3) according to Fhs＝FhsΔ s shrinkage FhsWherein Δ s ═ Δ s/ns;
4) turning to the step (4-2), solving the augmented constraint power flow equation again;
(4-5) when the delta s is less than or equal to the epsilon, finishing the calculation;
(5) according to the heat power P of the coupling equipment in the calculation result of the step (4)hCalculating electric power P by using the steady-state load flow model equation of the coupling equipment in the step (1-3), and calculating the pump lift h of the pipeline water pump according to the calculation result in the step (4)pCalculating the electric power P of the water pump by using the steady-state load flow model equation of the coupling equipment in the step (1-3)p；
(6) Based on the generalized inverse method, solving the optimal power flow model of the power system in the step (3), specifically comprising:
(6-1) giving values of contraction step length delta s, contraction multiple ns and calculation error epsilon;
(6-2) solving the augmented constrained power flow equation in the step (3-4) according to the modified equation in the step (3-4);
(6-3) if the calculation is converged,
1) solving a group of solutions X by an augmented constrained flow equationeAnd Fes，Xe＝[θ,φ,υ,yeb]T；
2) Storing the tidal current solution;
3) according to Fes＝Fes- Δ s contracts Fes, and when converging for k consecutive times, let Δ s ═ Δ s × ns, Fes＝Fes-Δs；
4) Turning to the step (6-2), solving the augmented constraint power flow equation again;
(6-4) if the calculation does not converge,
1) if the initial load flow has no solution, the calculation is finished;
2) recovering the last power flow solution;
3) according to Fes＝FesΔ s shrinkage FesWherein Δ s ═ Δ s/ns;
4) turning to the step (6-2), solving the augmented constraint power flow equation again;
(6-5) when Δ s ≦ ε, the calculation is ended.
3. The optimal power flow calculation method of the electrothermal coupling comprehensive energy system according to claim 2, characterized by comprising the following steps: the step (4-2) specifically comprises the following steps:
1) given a maximum number of iterations kmaxAnd a set of initial values
Figure FDA0002883512920000063
Figure FDA0002883512920000064
Wherein y ishb0, and let k be 0;
2) according to XhThe solution of a and deltab is carried out,
Figure FDA0002883512920000061
Δb＝[Δml,Δh,ΔTe,ΔPh,ΔFh]T
3) solving the generalized inverse matrix A of A+，A+＝AT(AAT)-1；
4) Solving for Δ XhAnd update Xh，ΔXh＝A+Δb，Xh＝Xh-ΔXhAnd let k be k + 1;
5) if Δ XhEpsilon is less than or equal to or the number of iterations k is more than kmaxAnd finishing the calculation, otherwise, turning to the step 2).
4. The optimal power flow calculation method of the electrothermal coupling comprehensive energy system according to claim 2 or 3, characterized by comprising the following steps of: the step (6-2) specifically comprises the following steps:
1) given epsilon, maximum number of iterations kmaxAnd a set of initial values Xe＝[θ,φ,υ,yeb]TWherein y iseb0, and let k be 0;
2) according to XeThe solution of a and deltab is carried out,
Figure FDA0002883512920000062
Δb＝[ΔP,ΔQ,ΔFe]T
3) solving the generalized inverse matrix A of A+，A+＝AT(AAT)-1；
4) Solving for Δ XeAnd update Xe，ΔXe＝A+Δb，Xe＝Xe-ΔXeAnd let k be k + 1;
5) if Δ XeEpsilon is less than or equal to or the number of iterations k is more than kmaxAnd finishing the calculation, otherwise, turning to the step 2).
5. The optimal power flow calculation method of the electrothermal coupling comprehensive energy system according to claim 2, characterized by comprising the following steps: the heat transfer working medium of the pipe network in the thermodynamic system is liquid, and the pipeline i
Figure FDA0002883512920000072
And
Figure FDA0002883512920000071
is a constant.
6. The optimal power flow calculation method of the electrothermal coupling comprehensive energy system according to claim 2 or 3, characterized by comprising the following steps of: the electric-thermal power conversion relation function f of the coupling equipment is in the form of a quadratic polynomial:
Ph＝khe2*Pe*Pe+khe1*Pe+khe0，khe2、khe1、khe0as fitting coefficient, PhIndicating the thermal power of the device, PeRepresenting the plant electrical power.
7. The optimal power flow calculation method of the electrothermal coupling comprehensive energy system according to claim 2 or 3, characterized by comprising the following steps of: the electric-thermal power conversion relation function f of the coupling equipment is in the form of a direct proportional relation.
CN201911040179.2A 2019-10-29 2019-10-29 Optimal power flow calculation method for electric-thermal coupling comprehensive energy system Active CN110752605B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911040179.2A CN110752605B (en) 2019-10-29 2019-10-29 Optimal power flow calculation method for electric-thermal coupling comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911040179.2A CN110752605B (en) 2019-10-29 2019-10-29 Optimal power flow calculation method for electric-thermal coupling comprehensive energy system

Publications (2)

Publication Number Publication Date
CN110752605A CN110752605A (en) 2020-02-04
CN110752605B true CN110752605B (en) 2021-05-14

Family

ID=69280957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911040179.2A Active CN110752605B (en) 2019-10-29 2019-10-29 Optimal power flow calculation method for electric-thermal coupling comprehensive energy system

Country Status (1)

Country Link
CN (1) CN110752605B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112000012B (en) * 2020-08-14 2022-04-22 浙江浙能绍兴滨海热电有限责任公司 Unit sliding pressure control optimization method and system based on thermoelectric load condition
CN112329217B (en) * 2020-10-23 2022-07-15 重庆大学 Quasi-steady state energy flow analysis method of electricity-water comprehensive energy system
CN112421635B (en) * 2020-10-29 2022-12-27 上海电力大学 Method for calculating multi-energy flow load flow inside and outside multi-energy micro-grid cluster network
CN113407896B (en) * 2021-06-15 2022-10-14 电子科技大学 Optimal stress and sample distribution determining method in nuclear power valve accelerated life test

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104636829A (en) * 2015-02-06 2015-05-20 河海大学 Decoupling algorithm for increasing temperature optimal power flow (OPF) calculation efficiency of electric power system
CN107887912A (en) * 2017-12-14 2018-04-06 清华大学 A kind of Static Power Flow analysis method and Static Power Flow analysis system
CN107947182A (en) * 2017-12-14 2018-04-20 清华大学 A kind of dynamic power flow analysis method and dynamic power flow analysis system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241662A (en) * 2018-09-29 2019-01-18 河海大学 A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies
CN109242365B (en) * 2018-11-06 2022-01-14 国电南瑞科技股份有限公司 Interval power flow calculation method of electricity-heat interconnection comprehensive energy system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104636829A (en) * 2015-02-06 2015-05-20 河海大学 Decoupling algorithm for increasing temperature optimal power flow (OPF) calculation efficiency of electric power system
CN107887912A (en) * 2017-12-14 2018-04-06 清华大学 A kind of Static Power Flow analysis method and Static Power Flow analysis system
CN107947182A (en) * 2017-12-14 2018-04-20 清华大学 A kind of dynamic power flow analysis method and dynamic power flow analysis system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Dynamic power flow calculation and open simulation of electro-thermal coupled system;jie hao 等;《2017 IEEE Conference on Energy Internet and Energy System Integration (EI2)》;20180104;第1-6页 *
考虑供热系统建模的综合能源系统最优能流计算方法;张义志 等;《电工技术学报》;20190228;第34卷(第3期);第562-570页 *

Also Published As

Publication number Publication date
CN110752605A (en) 2020-02-04

Similar Documents

Publication Publication Date Title
CN110752605B (en) Optimal power flow calculation method for electric-thermal coupling comprehensive energy system
CN110135631B (en) Electric comprehensive energy system scheduling method based on information gap decision theory
CN109242365A (en) A kind of Interval Power Flow calculation method of electric-thermal interconnection integrated energy system
CN111222213B (en) Dynamic simulation method and device for thermodynamic network
CN109492325B (en) Flow analysis method of multi-energy coupling system based on extended energy concentrator
CN108241782B (en) Mixed transient analysis method and mixed transient analysis system
CN110955954B (en) Method for reducing optimal load of layered decoupling electric heat comprehensive energy system
CN112417662B (en) Method for realizing dynamic hydraulic optimization of centralized heat supply pipe network system
CN111428351B (en) Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method
CN107947182B (en) Dynamic power flow analysis method and dynamic power flow analysis system
CN112016033B (en) Electric-thermal-gas comprehensive energy system tide calculation method based on forward-push back substitution method
CN111898816A (en) Dynamic state estimation method for comprehensive energy system
CN107887912B (en) Static power flow analysis method and static power flow analysis system
CN111463795B (en) Electric-thermal system multi-energy-flow analysis method based on improved Suhoff cooling formula
CN113379565A (en) Comprehensive energy system optimization scheduling method based on distributed robust optimization method
CN113255105B (en) Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure
CN110647040B (en) Safety control method and device of comprehensive energy system
CN108521125A (en) A kind of modeling of electric-gas coupled system probability continuous tide and its method for solving
Tian et al. Modeling and simulation for multi energy flow coupled network computing
CN108964039A (en) A kind of electric-gas coupled system steady stability of wind power plant containing DFIG domain tangent plane method for solving
CN111783309A (en) Dynamic simulation method of steam heating network based on internal conservation
CN111310310A (en) Thermodynamic system static power flow fast decoupling calculation method for quantity adjustment
CN110737993A (en) Multi-energy complementary system operation boundary analysis method considering load uncertainty
CN116542447A (en) Optimal scheduling method for electric heating system
CN107783941B (en) Solar energy compensation type electric boiler heat supply control method based on valley electricity price change

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant