CN110492533B - Control method and device of combined cooling heating and power system, computer and storage medium - Google Patents
Control method and device of combined cooling heating and power system, computer and storage medium Download PDFInfo
- Publication number
- CN110492533B CN110492533B CN201910745485.XA CN201910745485A CN110492533B CN 110492533 B CN110492533 B CN 110492533B CN 201910745485 A CN201910745485 A CN 201910745485A CN 110492533 B CN110492533 B CN 110492533B
- Authority
- CN
- China
- Prior art keywords
- network
- variable
- equation
- heating
- cooling
- 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
Links
- 238000001816 cooling Methods 0.000 title claims abstract description 241
- 238000010438 heat treatment Methods 0.000 title claims abstract description 220
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000012937 correction Methods 0.000 claims abstract description 273
- 239000011159 matrix material Substances 0.000 claims abstract description 151
- 230000008878 coupling Effects 0.000 claims description 94
- 238000010168 coupling process Methods 0.000 claims description 94
- 238000005859 coupling reaction Methods 0.000 claims description 94
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 60
- 239000013598 vector Substances 0.000 claims description 51
- 238000004364 calculation method Methods 0.000 claims description 37
- 238000004590 computer program Methods 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 8
- 238000010248 power generation Methods 0.000 claims description 6
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005057 refrigeration Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 239000008400 supply water Substances 0.000 description 2
- 241000284466 Antarctothoa delta Species 0.000 description 1
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a control method, a device, a computer and a storage medium of a combined cooling heating and power system, which are used for converting a Jacobian matrix into a diagonal polygonal matrix equation, obtaining the correction value of each variable in the combined cooling and heating and power system by solving the diagonal polygonal matrix equation, calculating the energy flow of a cooling network, a heating network and a power supply network in the combined cooling and heating and power system according to the correction value and controlling the output power of the system. By implementing the embodiment of the invention, the energy flow of the combined cooling heating and power system can be rapidly and accurately calculated, and the output power of the combined cooling heating and power system can be controlled.
Description
Technical Field
The present invention relates to the field of power systems, and in particular, to a method, an apparatus, a computer, and a medium for controlling a combined cooling, heating, and power system.
Background
The combined cooling heating and power system is a total energy system which is established on the basis of an energy cascade utilization concept and integrates the processes of refrigeration, heating and power generation. A typical combined cooling heating and power system generally comprises three parts, namely a power generation device, a refrigeration system and a heating system. In order to calculate the energy flow of the combined cooling heating and power system, the following two methods are generally adopted in the related art: (1) the unified iteration method combines energy flow equations of a cooling network, a heating network and a power supply network to solve uniformly, a unified iterative Jacobian matrix of the whole comprehensive energy micro-grid needs to be formed, and inversion operation is carried out on the Jacobian matrix, so that the calculation scale is large, and the calculation time consumption is long. (2) And a decoupling iteration method, which is used for respectively solving the energy flows of a cooling network, a heating network, a power supply network and an energy source station based on the coupling relation among various networks in the combined cooling heating and power micro-grid. However, in the method, only the balance nodes of the cooling network, the heating network and the power supply network can be set as a combined cooling heating and power supply unit, and the unreasonable node setting results in larger calculation error.
Therefore, how to quickly and accurately calculate the energy flow of the combined cooling heating and power system is a problem to be solved at present.
Disclosure of Invention
Accordingly, there is a need for a method, an apparatus, a computer and a storage medium for controlling a combined cooling, heating and power system, which can quickly and accurately calculate the energy flow of the combined cooling, heating and power system and control the output power of the combined cooling, heating and power system.
In one embodiment, a method for controlling a combined cooling, heating and power system is provided, the method including: acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system; generating a correction equation comprising a jacobian matrix, a plurality of elements in the jacobian matrix being functions of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; performing matrix row-column transformation on the correction equation to obtain a diagonal and edge-shaped matrix equation; decoupling operation is carried out on the diagonal addition matrix equation, and a first variable correction value of a cooling supply network, a second variable correction value of a heating supply network and a third variable correction value of the power supply network are obtained respectively; respectively calculating the energy flows of a cooling network, a heating network and a power supply network according to the first variable correction value, the second variable correction value and the third variable correction value; and controlling the output power distribution of the combined cooling heating and power system according to the energy flow of the cooling network, the heating network and the power supply network.
According to the control method of the combined cooling heating and power system, the Jacobian matrix is converted into the diagonal adding-edge matrix equation, the correction values of all variables in the combined cooling and power system are obtained by solving the diagonal adding-edge matrix equation, the energy flows of the cooling network, the heating network and the power supply network in the combined cooling and power system are calculated according to the correction values, and the output power of the system is controlled. The method is based on the deformation of the Jacobian matrix, and the balance nodes in the system do not need to be set as a joint supply unit, so the node setting is reasonable, and the calculation accuracy is high.
In one embodiment, the decoupling operation is performed on the diagonal polygonal matrix equation to obtain the variable correction values of the cooling network, the heating network and the power supply network respectively, and the method includes: expanding a diagonal and polygonal matrix equation to obtain a plurality of linear equations with coupling relations; performing iterative solution on the linear equations to obtain a cooling supply network correction equation, a heating supply network correction equation and a power supply network correction equation; and solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation respectively to obtain a first variable correction value, a second variable correction value and a third variable correction value.
In one embodiment, the correction equation is a relational expression among a Jacobian matrix, a variable correction amount and an unbalance equation; calculating the energy flows of a cooling network, a heating network and a power supply network according to the first variable correction value, the second variable correction value and the third variable correction value respectively, wherein the energy flows comprise: respectively correcting a first variable of a cooling supply network, a second variable of a heating supply network and a third variable of a supply network according to the first variable correction value, the second variable correction value and the third variable correction value; calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable; when the maximum values of the imbalance vectors meet a preset convergence condition, calculating variable values of balance nodes in the combined cooling, heating and power system according to the corrected first variable, second variable and third variable, and calculating energy flows of a cooling network, a heating network and a power supply network according to the variable values of the balance nodes.
In one embodiment, the method further comprises: when the maximum value of the imbalance vectors does not satisfy the preset convergence condition, repeating the following steps until the maximum value of the imbalance vectors satisfies the preset convergence condition: generating a correction equation comprising the Jacobian matrix again according to the corrected first variable, the second variable and the third variable; recalculating the first variable correction value, the second variable correction value and the third variable correction value according to the correction equation; and correcting the first variable, the second variable and the third variable again according to the first variable correction value, the second variable correction value and the third variable correction value obtained by recalculation, and calculating the values of a plurality of unbalance vectors in the unbalance equation.
In one embodiment, acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system comprises: acquiring a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in the combined cooling heating and power system as a hydraulic model: acquiring a linear equation set of water supply temperature and return water temperature of a cooling network or a heating network in a combined cooling heating and power system as a thermal model; acquiring an active power flow equation and a reactive power flow equation of a power system in a combined cooling heating and power system as a power system model; obtaining a relation between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system as a coupling element model, wherein the plurality of energy coupling elements comprise a thermoelectric power combined generation CHP unit, a refrigerator, an electric water pump and a heat pump.
In one embodiment, a control device for a combined cooling heating and power system is provided, which includes: the acquisition module is used for acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system; the generating module is used for generating a correction equation comprising a Jacobian matrix, and a plurality of elements in the Jacobian matrix are functions of a plurality of variables in a hydraulic model, a thermal model, a power system model and a coupling element model; the transformation module is used for carrying out matrix row-column transformation on the correction equation to obtain a diagonal and edge-shaped matrix equation; the first calculation module is used for performing decoupling operation on the diagonal polygonal matrix equation to respectively obtain a first variable correction value of the cooling network, a second variable correction value of the heating network and a third variable correction value of the power supply network; the second calculation module is used for calculating the energy flows of a cooling network, a heating network and a power supply network according to the first variable correction value, the second variable correction value and the third variable correction value; and the control module is used for controlling the output power distribution of the combined cooling heating and power system according to the energy flow of the cooling network, the heating network and the power supply network.
In one embodiment, the first operation module includes: the expansion unit is used for expanding the diagonal and polygonal matrix equation to obtain a plurality of linear equations with coupling relations; the first calculation unit is used for carrying out iterative solution on the linear equations to obtain a cooling supply network correction equation, a heating supply network correction equation and a power supply network correction equation; and the second calculation unit is used for respectively solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation to obtain a first variable correction value, a second variable correction value and a third variable correction value.
In one embodiment, the correction equation is a relational expression among a Jacobian matrix, a variable correction amount and an unbalance equation; the second calculation module includes: the correction unit is used for respectively correcting the first variable of the cooling supply network, the second variable of the heating supply network and the third variable of the supply network according to the first variable correction value, the second variable correction value and the third variable correction value; the third calculating unit is used for calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable; and the fourth calculating unit is used for calculating the variable value of a balance node in the combined cooling, heating and power system according to the corrected first variable, second variable and third variable when the maximum value of the imbalance vectors meets the preset convergence condition, and calculating the energy flow of a cooling network, the energy flow of a heating network and the energy flow of a power supply network according to the variable value of the balance node.
In one embodiment, a computer is provided, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the computer program is executed by the processor, the steps of the method of any embodiment of the present application are implemented.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any of the embodiments of the application.
Drawings
Fig. 1 is a schematic flowchart of a control method of a combined cooling, heating and power system according to an embodiment of the present invention;
fig. 2 is a flowchart illustrating a step S140 in a control method of a combined cooling, heating and power system according to an embodiment of the present invention;
fig. 3 is a flowchart illustrating a step S150 in the control method of the combined cooling, heating and power system according to an embodiment of the invention;
fig. 4 is a schematic pipeline diagram of a cooling network and a heating network in the cogeneration system according to an embodiment of the invention;
fig. 5 is a schematic structural diagram of a control device of the combined cooling, heating and power system according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a computer according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
For example, the present invention provides a control method for a combined cooling heating and power system, including: acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system; generating a correction equation comprising a jacobian matrix, a plurality of elements in the jacobian matrix being functions of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; performing matrix row-column transformation on the correction equation to obtain a diagonal and edge-shaped matrix equation; decoupling operation is carried out on the diagonal addition matrix equation, and a first variable correction value of a cooling supply network, a second variable correction value of a heating supply network and a third variable correction value of the power supply network are obtained respectively; respectively calculating the energy flows of a cooling network, a heating network and a power supply network according to the first variable correction value, the second variable correction value and the third variable correction value; and controlling the output power distribution of the combined cooling heating and power system according to the energy flow of the cooling network, the heating network and the power supply network.
According to the control method of the combined cooling heating and power system, the Jacobian matrix is converted into the diagonal adding-edge matrix equation, the correction values of all variables in the combined cooling and power system are obtained by solving the diagonal adding-edge matrix equation, the energy flows of the cooling network, the heating network and the power supply network in the combined cooling and power system are calculated according to the correction values, and the output power of the system is controlled. The method is based on the deformation of the Jacobian matrix, and the balance nodes in the system do not need to be set as a joint supply unit, so the node setting is reasonable, and the calculation accuracy is high.
In one embodiment, a method for controlling a combined cooling, heating and power system is provided, as shown in fig. 1, the method includes the following steps:
s101, acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system.
The hydraulic model comprises a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in the combined cooling heating and power system; the thermodynamic model comprises a linear equation set of the water supply temperature and the return water temperature of a cooling network or a heating network in a combined cooling heating and power system; the power system model comprises an active power flow equation and a reactive power flow equation of a power system in a combined cooling heating and power system; the coupling element model comprises a relational expression between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system. The multiple energy coupling elements comprise a thermoelectric power combined generation CHP unit, a refrigerator, an electric water pump and a heat pump.
S102, generating a correction equation comprising a Jacobian matrix. The correction equation is a relational expression among a Jacobian matrix, a variable correction amount and an unbalance equation.
In the present embodiment, the variable correction amount refers to a correction amount of a plurality of variables in a hydraulic model, a thermal model, a power system model, and a coupling element model. The variable correction is represented by Δ X, the jacobian matrix is represented by J, and the imbalance equation is represented by Δ f (X), then the correction equation is: -J Δ X ═ Δ f (X).
Wherein, each element in the jacobian matrix J is a partial derivative of the corresponding unbalance equation to the variable.
And S103, performing matrix row-column transformation on the correction equation to obtain a diagonal and edge-shaped matrix equation.
Specifically, a node where the coupling element is located is selected as a network division point, and the formula correction equation is arranged into a matrix equation in the form of the following diagonal addition through matrix row-column transformation. The node where the coupling element is located is, for example, a node where an energy station is located.
The diagonal-adding-edge matrix equation comprises a plurality of diagonal block matrixes and a plurality of boundary block matrixes. The diagonal block matrixes are respectively Jacobian matrixes related to a power supply network, a heat supply network and a cold supply network after the coupling elements are removed, and the boundary block matrix comprises Jacobian matrix elements related to the coupling elements.
And S104, performing decoupling operation on the diagonal polygonal matrix equation to respectively obtain a first variable correction value of the cooling network, a second variable correction value of the heating network and a third variable correction value of the power supply network.
In the present embodiment, the variable correction value of the cooling network is referred to as a first variable correction value, the variable correction value of the heating network is referred to as a second variable correction value, and the variable correction value of the supply network is referred to as a third variable correction value.
Specifically, the diagonal edge-added matrix equation is expanded and subjected to iterative processing, so that three equations without coupling relation can be obtained, wherein the three equations are respectively the first variable correction quantity Δ XcAnd coupling variable correction quantity DeltaXBSecond variable correction amount DeltaXhAnd coupling variable correction quantity DeltaXBThe third variable correction amount DeltaXeAnd coupling variable correction quantity DeltaXBThe relational expression (c) of (c). Wherein the coupling variable correction quantity DeltaXBThe calculation can be carried out according to a node-pipeline incidence matrix of the cooling network or the heat network and a loop-pipeline incidence matrix of the cooling network or the heat network. The three equations without coupling relation are respectively solved in the step to obtain a first variable correction quantity delta XcSecond variable correction amount DeltaXhAnd a third variable correction amount DeltaXeThe values of (b) are referred to as a first variable correction value, a second variable correction value, and a third variable correction value, respectively, that is, the diagonal addition matrix is subjected to decoupling processing.
And S105, respectively calculating the energy flows of the cooling network, the heating network and the power supply network according to the first variable correction value, the second variable correction value and the third variable correction value.
Specifically, the plurality of variables in the cooling network, the heating network, and the power supply network may be respectively corrected according to the first variable correction value, the second variable correction value, and the third variable correction value, and the energy flows of the cooling network, the heating network, and the power supply network may be respectively calculated according to the corrected values.
And S106, controlling the output power distribution of the combined cooling heating and power system according to the energy flow of the cooling network, the heating network and the power supply network.
In this step, the output power distribution of the combined cooling heating and power system is controlled according to the energy flow of the cooling network, the heating network and the power supply network. For example, when the energy flow of the cooling network is less than the preset energy flow threshold, the power output of the cooling network is increased to make the energy flow of the cooling network reach the expectation. By controlling the output power distribution of the combined cooling, heating and power system, the output power distribution of the combined cooling, heating and power system is more reasonable, and the maximization of energy utilization is facilitated.
According to the control method of the combined cooling heating and power system, the Jacobian matrix is converted into the diagonal adding-edge matrix equation, the correction values of all variables in the combined cooling and power system are obtained by solving the diagonal adding-edge matrix equation, the energy flows of the cooling network, the heating network and the power supply network in the combined cooling and power system are calculated according to the correction values, and the output power of the system is controlled. The method is based on the deformation of the Jacobian matrix, and the balance nodes in the system do not need to be set as a joint supply unit, so the node setting is reasonable, and the calculation accuracy is high.
In one embodiment, step S101 includes:
1. acquiring a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in a combined cooling heating and power system as hydraulic models; for example, the cold/hot water flow and the pressure head in the hydraulic model should satisfy the following equations:
wherein, A is a node-pipeline incidence matrix of a cooling/heating network, m is a pipeline flow vector, mq is a cooling/heating network node flow vector, B is a loop-pipeline incidence matrix of the cooling/heating network, delta h is a pipeline head loss vector, and K is a pipeline resistance coefficient matrix.
2. And acquiring a linear equation set of the water supply temperature and the return water temperature of a cooling network or a heating network in the combined cooling heating and power system as a thermal model. For example, the cold/hot power of the cold/heat source and the load point can be calculated by the following formula:
Φ＝cwmqs(Ts-Tr)；
where Φ is the cold/hot power at the node, CwM is specific heat capacity of cooling/heating mediumqIs the traffic of the node, TsTemperature of water supply to the node, TrAnd s is the identifier of a cold/heat source or load. For ease of distinction, one embodiment is where S is +1 when a heating network is indicated and S is-1 when a cooling network is indicated.
For each pipeline, the temperature rise and temperature drop equation of the pipeline is as follows:
wherein, TinAnd ToutRespectively the inlet and outlet temperature of the pipeline, TaIs the ambient temperature, L is the length of the pipe, and λ is the heat transfer coefficient of the pipe.
For each node, there is a node temperature hybrid model as follows:
Σ(minTin)＝(Σmout)Tout；
according to the pipeline temperature rise and drop equation and the node temperature mixed model, a water supply temperature and return water temperature calculation equation can be listed for each node, and two linear equation sets related to the water supply temperature and the return water temperature can be obtained after arrangement.
Wherein, CsAnd CrFor a matrix of coefficients related to network structure and traffic, bsAnd brVectors associated with cold/heat source supply water temperature and load return water temperature.
3. And acquiring an active power flow equation and a reactive power flow equation of a power system in the combined cooling heating and power system as a power system model. For example, if the power flow of the power system is solved by the newton-pulling method, the active power flow P of the power system can be calculated according to the following formulaiAnd reactive power flow Qi：
4. Obtaining a relation between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system as a coupling element model, wherein the plurality of energy coupling elements comprise a thermoelectric power combined generation CHP unit, a refrigerator, an electric water pump and a heat pump. In a real-time scene, the CHP unit, the refrigerating machine and the heat pump are arranged in an energy station of a combined cooling, heating and power supply park, and the electric water pumps are respectively arranged at the energy station and a user side and are configured in two stages.
Wherein, the CHP unit adopts a back pressure unit with a thermoelectric ratio cmThe constant value is adopted, the constant heat power mode is adopted for working, so the relationship between the generated energy and the heat supply quantity is as follows: pCHP＝ΦCHP/cm。
Wherein, the refrigerator is an absorption refrigerator, and supplies cold by absorbing partial heat energy supplied by the CHP unit, so that the refrigerator consumes thermal power phihAnd output cold power phicThe relationship of (1) is: phic＝COPACΦhWherein, COPACIs the energy efficiency ratio of the absorption chiller.
Wherein, for the electric water pump, assume mpThe water flow through the water pump, H is the lift of the water pump, eta is the efficiency of the water pump,
for the power factor angle of the electric water pumps, the electric power consumed by each electric water pump is as follows:
wherein, for the heat pump, the relation between the consumed electric power and the produced thermal power is phiHP＝COPHPPHP，COPHPIs the heat-to-power ratio of the heat pump.
In one embodiment, at step S102: and (3) combining variables of a cooling network and a heating network, and performing unified iterative solution on the micro-grid energy flow of the combined cooling heating and power supply park by adopting a Newton Raphson method to obtain a correction equation-J delta X ═ delta F (X) comprising a Jacobian matrix. Wherein the unbalance amount equation Δ f (x) is as follows:
in one embodiment, in step S103, on the basis of a unified iterative algorithm, a large-scale unified iterative jacobian matrix is converted into a diagonal edge-added form (BBDF) matrix, so as to implement decoupling calculation of an energy flow, which is specifically as follows: selecting a node where a coupling element is located, selecting a node where an energy station is located as a network division point, and arranging the formula correction equation into the following diagonal edge-added form matrix equation through matrix row-column transformation:
wherein, the diagonal block matrix Jee、JhhAnd JccRespectively, Jacobian matrix and boundary block matrix J related to power supply network, heat supply network and cold supply network after removing coupling elementsBB、JiB、JBi(i ═ 1,2,3) contains the jacobian matrix elements associated with the coupling elements. Δ Xe、ΔXhAnd Δ XcCorrection quantities, DeltaX, of variables associated with the supply network, the heating network and the cooling network, respectively, after removal of the coupling variableBThe correction quantity of the coupling variable corresponding to the coupling element is obtained; Δ Fe、ΔFhAnd Δ FcRespectively, the imbalance equations, Δ F, associated with the supply network, the heating network, the cooling network after removal of the coupling variablesBThe imbalance equation for the coupling element.
In one embodiment, as shown in fig. 2, the step S104 may include the following steps:
and S1041, expanding a diagonal and polygonal matrix equation to obtain a plurality of linear equations with coupling relations. For example, the above diagonal form matrix equation is expanded to obtain the following four linear equations:
JeeΔXe+J1BΔXB＝ΔFe (1)
JhhΔXh+J2BΔXB＝ΔFh (2)
JccΔXc+J3BΔXB＝ΔFc (3)
JB1ΔXe+JB2ΔXh+JB3ΔXc+JBBΔXB＝ΔFB (4)
and S1042, performing iterative solution on the linear equations to obtain a cooling network correction equation, a heating network correction equation and a power supply network correction equation.
In this step, the formulae (1), (2) and (3) are substituted with the formula (4), and Δ X is eliminatede、ΔXh、ΔXcThe equation A.DELTA.X can be obtainedBB, wherein:
by solving the system of linear equations, the correction quantity DeltaX related to the coupling variable can be obtainedBThen Δ XBThe following equations can be obtained by substituting the equations (1) to (3) respectively:
the three equations in the equation set are respectively a cooling network correction equation, a heating network correction equation and a power supply network correction equation. Therefore, no coupling relation exists among the cooling supply network correction equation, the heating supply network correction equation and the power supply network correction equation.
And S1043, respectively solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation to obtain a first variable correction value, a second variable correction value and a third variable correction value.
In the step, a cooling network correction equation, a heating network correction equation and a power supply network correction equation are separately solved to obtain a correction value delta X of a cooling network variablecCorrected value delta X of heat supply network variablehAnd correction values DeltaX of supply network variablese。
In one embodiment, as shown in fig. 3, the step S105 may include the following steps:
and S1051, respectively correcting the first variable of the cooling network, the second variable of the heating network and the third variable of the power supply network according to the first variable correction value, the second variable correction value and the third variable correction value.
And S1052, calculating the values of the imbalance vectors in the imbalance equation according to the corrected first variable, second variable and third variable.
S1053, judging whether the maximum value of the imbalance vectors meets the preset convergence condition, if yes, executing the step S1054, otherwise, returning to the step S102.
That is, when the maximum value of the plurality of unbalance vectors does not satisfy the preset convergence condition, the following steps are repeated until the maximum value of the plurality of unbalance vectors satisfies the preset convergence condition: generating a correction equation comprising the Jacobian matrix again according to the corrected first variable, the second variable and the third variable; recalculating the first variable correction value, the second variable correction value and the third variable correction value according to the correction equation; and correcting the first variable, the second variable and the third variable again according to the first variable correction value, the second variable correction value and the third variable correction value obtained by recalculation, and calculating the values of a plurality of unbalance vectors in the unbalance equation.
And S1054, calculating the variable value of the balance node in the combined cooling, heating and power system according to the corrected first variable, second variable and third variable.
And S1055, calculating the energy flow of the cooling network, the energy flow of the heating network and the energy flow of the power supply network according to the variable values of the balance nodes.
In the present embodiment, it is preferred that,obtaining corrected values DeltaX for individual subnetwork variablese、ΔXhAnd Δ XcAnd then, correcting all variables, iterating again to obtain the unbalance amount, sequentially constructing a boundary block matrix and a diagonal block matrix in the BBDF according to the new unbalance amount, and performing decoupling operation again until a convergence output result.
According to the combined cooling heating and power control method, the operation principle of the unified iterative method is not changed, decoupling optimization is performed in the operation process, the node where the energy coupling element is located is taken as a matrix separation point in the calculation process, the diagonal block and the boundary block matrix with small scale are directly formed, the unified iterative Jacobian matrix with large scale does not need to be formed, inversion operation on the large-scale unified iterative Jacobian matrix is avoided in the calculation process, the operation time is obviously reduced, and the calculation efficiency is improved.
In addition, the method is not influenced by a network structure, can be suitable for cooling/heating and power supply networks in various forms such as a radiation type and a ring network type, is not influenced by the type and the operation mode of the coupling element of the combined cooling and heating micro-grid, and has wider application range. Theoretically, the larger the scale of the energy network is, the more obvious the advantages of the combined cooling heating and power supply control method provided by the invention are, and therefore, with the further development of a future comprehensive energy system, the algorithm has more prominent practical application value.
In one embodiment, the BBDF-based combined cooling heating and power park microgrid energy flow decoupling calculation method provided herein is verified by taking a certain combined cooling heating and power park microgrid as an example. The park is divided into two energy supply areas, which are two energy stations, energy coupling elements such as CHP units and hot water type absorption refrigerators are arranged in the energy stations, and electric circulating water pumps are respectively arranged at a cold/hot source end and a cold/hot load end. The cooling/heating pipelines in the microgrid are shown in fig. 4, and the cooling network and the heating network have the same structure and both comprise 46 nodes and 46 sections of pipelines. The water supply network and the water return network of the cold/heat supply pipeline have the same structure, and are not described in detail herein. Energy stations I and II are located at 46 and 45 nodes of the grid respectively, with energy station I being set as a balanced node of the cooling/heating grid. The microgrid power supply network comprises 90 nodes, wherein the 90 nodes are nodes connected with a public power distribution network, the nodes are assumed to be power supply network balance nodes, the energy stations I and II are assumed to be PV nodes, and the other nodes are assumed to be PQ nodes.
When the unified iteration method is adopted, a large-scale unified iterative jacobian matrix needs to be formed, and the size of the matrix is 444 by 444. The method comprises the steps of selecting energy stations I and II as network division points, converting an original large-scale unified iterative Jacobian matrix into a Jacobian matrix of a BBDF model through matrix row-column transformation, wherein a diagonal block matrix is the Jacobian matrix of a power supply network, a heat supply network and a cold supply network except coupling elements, a boundary block matrix comprises the coupling elements, and the rest parts except the diagonal block matrix and the boundary block matrix are zero matrices.
Given a convergence parameter of 10-4The energy flow calculation is carried out on the above calculation examples by respectively adopting a unified iteration method and a BBDF-based decoupling algorithm, the program is converged after 16 iterations, the calculation results and the iteration times of the unified iteration method and the BBDF-based decoupling algorithm are completely the same, and the obtained results are shown as follows.
(1) Calculation of energy flow in power supply network
The phase angle and amplitude of the power supply network node voltage are shown in table 1:
TABLE 1 phase angle and amplitude of partial node voltage of power supply network
| |
1 | 10 | 20 | 30 | 40 | 50 |
| Voltage angle (°) | 16.475 | 18.505 | 18.522 | 18.485 | 16.787 | 15.375 |
| Voltage amplitude (p.u.) | 1.002 | 1.031 | 1.031 | 1.029 | 1.002 | 1.008 |
| |
60 | 70 | 80 | 88 | 89 | 90 |
| Voltage angle (°) | 18.767 | 18.777 | 13.437 | 22.227 | 21.627 | 0.000 |
| Voltage amplitude (p.u.) | 1.031 | 1.032 | 1.001 | 1.050 | 1.050 | 1.050 |
(2) Heating network energy flow calculation results
The heat supply network pipeline flow and node supply and return water temperature meters 2,3 and 4 are as follows:
TABLE 2 flow rates of the pipes of the heating network
TABLE 3 supply water temperature of each node of the heating network
TABLE 4 return water temperature of each node of heat supply network
(3) Cooling network energy flow calculation results
The flow of each pipeline of the cooling supply network and the supply and return water temperature of each node are shown in the tables 5, 6 and 7:
TABLE 5 flow of each pipe of cooling network
TABLE 6 Water supply temperature at each node of cooling network
TABLE 7 Return water temperature of each node of cooling network
(4) Variable value of cold/hot side of energy station
TABLE 8 Cold/Heat output Power of energy station
(5) Variable values at generator set and distribution network nodes of energy station
TABLE 9 active power output of energy station and distribution network node
| Node name | Active output/MW |
| Energy station I | 1.778 |
| Energy station II | 2.000 |
| Power distribution network node | -1.979 |
(6) Variation of coupling element
TABLE 10 coupling element variable values
| Name of variable | Value of variable |
| Heat pump consuming electric power (energy station I) | 0.4445MW |
| Heat pump heating power (energy station I) | 1.3335MW |
| Heat pump consuming electric power (energy station II) | 0.5000MW |
| Heat pump heating power (energy station II) | 1.5000MW |
| Source water pump consumes electric power (energy station I) | 0.1404MW |
| Source water pump consumes electric power (energy station II) | 0.1456MW |
| Total electric power consumed by water pump at load end | 0.1541MW |
(7) Comparison of operation times
The total time consumption of the two algorithms in the iterative process is shown in table 11, and it can be seen that, in the example solution of the present invention, the total time consumption of the conventional unified iterative algorithm is about 0.3136s, while the total time consumption of the BBDF-based decoupling algorithm is about 0.1873 s. In contrast, the BBDF-based decoupling algorithm can reduce 37.41% of the computation time, and the two algorithms are completely consistent in the computation result. Therefore, the BBDF-based decoupling algorithm can effectively improve the calculation efficiency of the microgrid energy flow calculation in the combined cooling heating and power supply park, and has the same calculation accuracy as the unified iterative algorithm.
TABLE 11 comparison of operating times
| Method | Calculating time/s |
| Unified iteration | 0.313312 |
| BBDF | 0.196169 |
Particularly, as the network scale is further increased, if a unified iterative algorithm is adopted, the time consumption ratio of the formation of the large-scale unified iterative jacobian matrix and the inversion process of the large-scale unified iterative jacobian matrix is obviously increased, but the BBDF-based energy flow decoupling algorithm provided by the invention can directly form a diagonal block matrix and a boundary block matrix with smaller scale in the BBDF model, and can avoid inversion of the large-scale unified iterative jacobian matrix, so that the larger the scale of the combined cooling heating and power supply microgrid is, particularly the larger the scale of the cooling/heating pipe network is, the more the operation time consumption of the provided algorithm is reduced, the more the operation advantages are highlighted, and the calculation result is shown in table 12. Therefore, with the further development of the future comprehensive energy system, the algorithm has more prominent practical application value.
TABLE 12 BBDF energy flow decoupling algorithm reduction percentage consumed time under different energy grid scales
| Number of nodes of power grid | Number of nodes of cold/hot pipe network | The proposed algorithm reduces the percentage of time consumed |
| 90 | 46 | 37.41％ |
| 180 | 46 | 40.02％ |
| 180 | 121 | 67.91％ |
In one embodiment, a control device 60 of a combined cooling heating and power system is provided, as shown in fig. 5, the device includes: the acquisition module 601 is used for acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system; a generating module 602 for generating a correction equation comprising a jacobian matrix, a plurality of elements of the jacobian matrix being a function of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; a transformation module 603, configured to perform matrix row-column transformation on the modified equation to obtain a diagonal-added-edge matrix equation; the first calculating module 604 is configured to perform decoupling operation on the diagonal edge-added matrix equation, and obtain a first variable correction value of the cooling network, a second variable correction value of the heating network, and a third variable correction value of the power supply network, respectively; a second calculation module 605, configured to calculate energy flows of the cooling network, the heating network, and the power supply network according to the first variable correction value, the second variable correction value, and the third variable correction value; a control module 606, configured to control output power distribution of the combined cooling heating and power system according to energy flows of the cooling network, the heating network, and the power supply network.
The control device of the combined cooling heating and power system converts the Jacobian matrix into a diagonal edge-added matrix equation, obtains the correction value of each variable in the combined cooling and power system by solving the diagonal edge-added matrix equation, calculates the energy flow of a cooling network, a heating network and a power supply network in the combined cooling and power system according to the correction value and controls the output power of the system. The method is based on the deformation of the Jacobian matrix, and the balance nodes in the system do not need to be set as a joint supply unit, so the node setting is reasonable, and the calculation accuracy is high.
In one embodiment, the first operation module comprises: the expansion unit is used for expanding the diagonal and polygonal matrix equation to obtain a plurality of linear equations with coupling relations; the first calculation unit is used for carrying out iterative solution on the linear equations to obtain a cooling supply network correction equation, a heating supply network correction equation and a power supply network correction equation; and the second calculation unit is used for respectively solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation to obtain the first variable correction value, the second variable correction value and the third variable correction value.
In one embodiment, the correction equation is a relational expression between a jacobian matrix, a variable correction amount and an unbalance equation; the second calculation module includes: a correction unit, configured to correct the first variable of the cooling network, the second variable of the heating network, and the third variable of the supplying network according to the first variable correction value, the second variable correction value, and the third variable correction value, respectively; the third calculating unit is used for calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable; and the fourth calculating unit is used for calculating the variable value of a balance node in the combined cooling heating and power system according to the corrected first variable, the second variable and the third variable when the maximum value of the imbalance vectors meets a preset convergence condition, and calculating the energy flow of the cooling network, the energy flow of the heating network and the energy flow of the power supply network according to the variable value of the balance node.
In one embodiment, a control device of a combined cooling, heating and power system includes a functional module corresponding to the control method of the combined cooling, heating and power system according to any embodiment of the present invention.
In one embodiment, as shown in fig. 6, there is provided a computer, namely a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the computer program to implement the following steps: acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system; generating a modified equation comprising a jacobian matrix, a plurality of elements in the jacobian matrix being a function of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; performing matrix row-column transformation on the correction equation to obtain a diagonal and edge-shaped matrix equation; performing decoupling operation on the diagonal edge-added matrix equation to respectively obtain a first variable correction value of a cooling network, a second variable correction value of a heating network and a third variable correction value of the power supply network; calculating the energy flows of the cooling network, the heating network and the power supply network according to the first variable correction value, the second variable correction value and the third variable correction value; controlling the output power distribution of the combined cooling heating and power system according to the energy flow of the cooling network, the heating network and the power supply network.
In one embodiment, the processor, when executing the computer program, performs the steps of: the step of performing decoupling operation on the diagonal edge-added matrix equation to respectively obtain the variable correction values of the cooling network, the heating network and the power supply network comprises the following steps of: expanding the diagonal adding edge matrix equation to obtain a plurality of linear equations with coupling relations; performing iterative solution on the linear equations to obtain a cooling network correction equation, a heating network correction equation and a power supply network correction equation; and solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation respectively to obtain the first variable correction value, the second variable correction value and the third variable correction value.
In one embodiment, the processor, when executing the computer program, performs the steps of: the step of calculating the energy flows of the cooling network, the heating network and the supply network, respectively, from the first, second and third correction variables comprises: correcting a first variable of the cooling network, a second variable of the heating network and a third variable of the power supply network according to the first variable corrected value, the second variable corrected value and the third variable corrected value respectively; calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable; when the maximum value of the imbalance vectors meets a preset convergence condition, calculating the variable value of a balance node in the combined cooling heating and power system according to the corrected first variable, second variable and third variable, and calculating the energy flow of the cooling network, the energy flow of the heating network and the energy flow of the power supply network according to the variable value of the balance node.
In one embodiment, the processor, when executing the computer program, further performs the steps of: when the maximum value of the imbalance vectors does not satisfy a preset convergence condition, repeating the following steps until the maximum value of the imbalance vectors satisfies the preset convergence condition: generating a correction equation comprising a Jacobian matrix again according to the corrected first variable, the second variable and the third variable; recalculating the first variable correction value, the second variable correction value and the third variable correction value according to the correction equation; and correcting the first variable, the second variable and the third variable again according to the first variable correction value, the second variable correction value and the third variable correction value which are obtained through recalculation, and calculating the values of a plurality of unbalance vectors in the unbalance equation.
In one embodiment, the processor, when executing the computer program, performs the steps of: the method for acquiring the hydraulic model, the thermal model, the electric power system model and the coupling element model of the combined cooling heating and power system comprises the following steps: acquiring a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in the combined cooling heating and power system as a hydraulic model: acquiring a linear equation set of water supply temperature and return water temperature of a cooling network or a heating network in the combined cooling heating and power system as a thermodynamic model; acquiring an active power flow equation and a reactive power flow equation of a power system in the combined cooling heating and power system as a power system model; obtaining a relational expression between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system as a coupling element model, wherein the plurality of energy coupling elements comprise a combined thermal power generation (CHP) unit, a refrigerator, an electric water pump and a heat pump.
In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of: acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system; generating a modified equation comprising a jacobian matrix, a plurality of elements in the jacobian matrix being a function of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; performing matrix row-column transformation on the correction equation to obtain a diagonal and edge-shaped matrix equation; performing decoupling operation on the diagonal edge-added matrix equation to respectively obtain a first variable correction value of a cooling network, a second variable correction value of a heating network and a third variable correction value of the power supply network; calculating the energy flows of the cooling network, the heating network and the power supply network according to the first variable correction value, the second variable correction value and the third variable correction value; controlling the output power distribution of the combined cooling heating and power system according to the energy flow of the cooling network, the heating network and the power supply network.
In one embodiment, the computer program when executed by the processor implements the steps of: the step of performing decoupling operation on the diagonal edge-added matrix equation to respectively obtain the variable correction values of the cooling network, the heating network and the power supply network comprises the following steps of: expanding the diagonal adding edge matrix equation to obtain a plurality of linear equations with coupling relations; performing iterative solution on the linear equations to obtain a cooling network correction equation, a heating network correction equation and a power supply network correction equation; and solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation respectively to obtain the first variable correction value, the second variable correction value and the third variable correction value.
In one embodiment, the computer program when executed by the processor implements the steps of: the step of calculating the energy flows of the cooling network, the heating network and the supply network, respectively, from the first, second and third correction variables comprises: correcting a first variable of the cooling network, a second variable of the heating network and a third variable of the power supply network according to the first variable corrected value, the second variable corrected value and the third variable corrected value respectively; calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable; when the maximum value of the imbalance vectors meets a preset convergence condition, calculating the variable value of a balance node in the combined cooling heating and power system according to the corrected first variable, second variable and third variable, and calculating the energy flow of the cooling network, the energy flow of the heating network and the energy flow of the power supply network according to the variable value of the balance node.
In one embodiment, the computer program when executed by the processor further performs the steps of: when the maximum value of the imbalance vectors does not satisfy a preset convergence condition, repeating the following steps until the maximum value of the imbalance vectors satisfies the preset convergence condition: generating a correction equation comprising a Jacobian matrix again according to the corrected first variable, the second variable and the third variable; recalculating the first variable correction value, the second variable correction value and the third variable correction value according to the correction equation; and correcting the first variable, the second variable and the third variable again according to the first variable correction value, the second variable correction value and the third variable correction value which are obtained through recalculation, and calculating the values of a plurality of unbalance vectors in the unbalance equation.
In one embodiment, the computer program when executed by the processor implements the steps of: the method for acquiring the hydraulic model, the thermal model, the electric power system model and the coupling element model of the combined cooling heating and power system comprises the following steps: acquiring a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in the combined cooling heating and power system as a hydraulic model: acquiring a linear equation set of water supply temperature and return water temperature of a cooling network or a heating network in the combined cooling heating and power system as a thermodynamic model; acquiring an active power flow equation and a reactive power flow equation of a power system in the combined cooling heating and power system as a power system model; obtaining a relational expression between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system as a coupling element model, wherein the plurality of energy coupling elements comprise a combined thermal power generation (CHP) unit, a refrigerator, an electric water pump and a heat pump.
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A control method of a combined cooling heating and power system is characterized by comprising the following steps:
acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system;
generating a modified equation comprising a jacobian matrix, a plurality of elements in the jacobian matrix being a function of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; the correction equation is a relational expression among a Jacobian matrix J, a variable correction value delta X and an unbalance equation delta F (X): -J Δ X ═ Δ f (X); wherein the unbalance equation Δ f (x) has the following relation:
wherein, cwM is specific heat capacity of cooling/heating mediumhqAnd mhRespectively, flow value vectors of a heating network node and a pipeline, s is an identifier of a cold/heat source or a load, s is equal to +1 when representing a heating network, and is equal to-1 when representing a cooling network, ThsAnd ThrWater supply temperature and return temperature vectors, phi, of heat supply network nodes, respectivelyhspFor a given heat supply network node power vector, B is a loop-pipeline correlation matrix of a cold/heat supply network, K is a pipeline resistance coefficient matrix, ChsAnd ChrFor a coefficient matrix relating to heating network structure and flow, bhsAnd bhrVectors related to the water supply temperature and the load return water temperature of the heat supply network are respectively; m iscqAnd mcFlow value vectors, T, for cooling network nodes and pipes, respectivelycsAnd TcrRespectively the supply and return water temperature vectors, phi, of the cooling network nodescspFor a given cooling network node power vector, CcsAnd CcrFor a matrix of coefficients related to the cooling network structure and flow, bcsAnd bcrVectors related to the water supply temperature and the load return water temperature of the cooling network are respectively;
performing matrix row-column transformation on the correction equation to obtain a diagonal-angle and edge-shaped matrix equation:
wherein, the diagonal block matrix Jee、JhhAnd JccRespectively, Jacobian matrix, boundary block matrix J, associated with the power supply network, heat supply network and cold supply network after removal of coupling elementsBB、JiB、JBi(i ═ 1,2,3) comprising jacobian matrix elements, Δ X, associated with the coupling elementse、ΔXhAnd Δ XcCorrection quantities, DeltaX, of variables associated with the supply network, the heating network and the cooling network after removal of the coupling variableBThe correction quantity of the coupling variable corresponding to the coupling element is obtained; Δ Fe、ΔFhAnd Δ FcRespectively connected with a power supply network, a heat supply network and a cooling network after removing coupling variablesNetwork-related imbalance equation, Δ FBAn imbalance equation corresponding to the coupling element;
performing decoupling operation on the diagonal edge-added matrix equation to respectively obtain a first variable correction value delta X of the cooling networkcSecond variable correction value delta X of heat supply networkhAnd third variable correction Δ X of the supply networke；
According to said first variable correction value DeltaXcThe second variable correction value DeltaXhAnd said third variable correction value DeltaXeCorrecting a first variable of the cooling network, a second variable of the heating network and a third variable of the supply network, respectively;
calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable;
when the maximum value of the imbalance vectors meets a preset convergence condition, calculating the variable value of a balance node in the combined cooling heating and power system according to the corrected first variable, second variable and third variable, and calculating the energy flows of the cooling network, the heating network and the power supply network according to the variable value of the balance node;
controlling the output power distribution of the combined cooling heating and power system according to the energy flows of the cooling network, the heating network and the power supply network.
2. The method of claim 1, wherein performing a decoupling operation on the diagonal-polygonal matrix equation to obtain modified variable values for the cooling network, the heating network, and the power supply network, respectively, comprises:
expanding the diagonal adding edge matrix equation to obtain a plurality of linear equations with coupling relations;
performing iterative solution on the linear equations to obtain a cooling network correction equation, a heating network correction equation and a power supply network correction equation;
and respectively solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation to obtain the first variable correction value, the second variable correction value and the third variable correction value.
3. The method according to claim 2, wherein the diagonal-edged matrix equation is expanded to obtain a plurality of linear equations with coupling relations, respectively:
JeeΔXe+J1BΔXB＝ΔFe
JhhΔXh+J2BΔXB＝ΔFh
JccΔXc+J3BΔXB＝ΔFc
JB1ΔXe+JB2ΔXh+JB3ΔXc+JBBΔXB＝ΔFB
and the iterative solution is carried out on the linear equations to obtain a cooling network correction equation, a heating network correction equation and a power supply network correction equation which are respectively as follows:
4. the method of claim 1, further comprising:
when the maximum value of the imbalance vectors does not satisfy a preset convergence condition, repeating the following steps until the maximum value of the imbalance vectors satisfies the preset convergence condition:
generating a correction equation comprising a Jacobian matrix again according to the corrected first variable, the second variable and the third variable;
recalculating the first variable correction value, the second variable correction value and the third variable correction value according to the correction equation;
and correcting the first variable, the second variable and the third variable again according to the first variable correction value, the second variable correction value and the third variable correction value which are obtained through recalculation, and calculating the values of a plurality of unbalance vectors in the unbalance equation.
5. The method of claim 4, wherein the obtaining of the hydraulic model, the thermal model, the electric power system model and the coupling element model of the combined cooling heating and power system comprises:
acquiring a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in the combined cooling heating and power system as hydraulic models;
acquiring a linear equation set of water supply temperature and return water temperature of a cooling network or a heating network in the combined cooling heating and power system as a thermodynamic model;
acquiring an active power flow equation and a reactive power flow equation of a power system in the combined cooling heating and power system as a power system model;
obtaining a relational expression between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system as a coupling element model, wherein the plurality of energy coupling elements comprise a combined thermal power generation (CHP) unit, a refrigerator, an electric water pump and a heat pump.
6. A control device for a combined cooling heating and power system, the device comprising:
the acquisition module is used for acquiring a hydraulic model, a thermal model, an electric power system model and a coupling element model of the combined cooling heating and power system;
a generation module to generate a correction equation comprising a Jacobian matrix, a plurality of elements in the Jacobian matrix being a function of a plurality of variables in the hydraulic model, the thermal model, the power system model, and the coupling element model; the correction equation is a relational expression among a Jacobian matrix J, a variable correction value delta X and an unbalance equation delta F (X): -J Δ X ═ Δ f (X); wherein the unbalance equation Δ f (x) has the following relation:
wherein, cwM is specific heat capacity of cooling/heating mediumhqAnd mhRespectively, flow value vectors of a heating network node and a pipeline, s is an identifier of a cold/heat source or a load, s is equal to +1 when representing a heating network, and is equal to-1 when representing a cooling network, ThsAnd ThrWater supply temperature and return temperature vectors, phi, of heat supply network nodes, respectivelyhspFor a given heat supply network node power vector, B is a loop-pipeline correlation matrix of a cold/heat supply network, K is a pipeline resistance coefficient matrix, ChsAnd ChrFor a coefficient matrix relating to heating network structure and flow, bhsAnd bhrVectors related to the water supply temperature and the load return water temperature of the heat supply network are respectively; m iscqAnd mcFlow value vectors, T, for cooling network nodes and pipes, respectivelycsAnd TcrRespectively the supply and return water temperature vectors, phi, of the cooling network nodescspFor a given cooling network node power vector, CcsAnd CcrFor a matrix of coefficients related to the cooling network structure and flow, bcsAnd bcrVectors related to the water supply temperature and the load return water temperature of the cooling network are respectively;
the transformation module is used for performing matrix row-column transformation on the correction equation to obtain a diagonal-angle-added-edge matrix equation:
wherein, the diagonal block matrix Jee、JhhAnd JccRespectively, Jacobian matrix, boundary block matrix J, associated with the power supply network, heat supply network and cold supply network after removal of coupling elementsBB、JiB、JBi(i ═ 1,2,3) comprising jacobian matrix elements, Δ X, associated with the coupling elementse、ΔXhAnd Δ XcRespectively after removing coupling variablesCorrection of relevant variables, DeltaX, of power supply networks, heating networks, cooling networksBThe correction quantity of the coupling variable corresponding to the coupling element is obtained; Δ Fe、ΔFhAnd Δ FcRespectively, the imbalance equations, Δ F, associated with the supply network, the heating network, the cooling network after removal of the coupling variablesBAn imbalance equation corresponding to the coupling element;
a first calculation module, configured to perform decoupling operation on the diagonal edge-added matrix equation to obtain first variable correction values Δ X of the cooling network respectivelycSecond variable correction value delta X of heat supply networkhAnd third variable correction Δ X of the supply networke；
A second calculation module for calculating a correction value Δ X based on the first variablecThe second variable correction value DeltaXhAnd said third variable correction value DeltaXeCorrecting a first variable of the cooling network, a second variable of the heating network and a third variable of the supply network, respectively; calculating values of a plurality of unbalance vectors in the unbalance equation according to the corrected first variable, the corrected second variable and the corrected third variable; when the maximum value of the imbalance vectors meets a preset convergence condition, calculating the variable value of a balance node in the combined cooling heating and power system according to the corrected first variable, second variable and third variable, and calculating the energy flows of the cooling network, the heating network and the power supply network according to the variable value of the balance node;
and the control module is used for controlling the output power distribution of the combined cooling heating and power system according to the energy flows of the cooling network, the heating network and the power supply network.
7. The apparatus of claim 6, wherein the first computing module comprises:
the expansion unit is used for expanding the diagonal and polygonal matrix equation to obtain a plurality of linear equations with coupling relations;
the first calculation unit is used for carrying out iterative solution on the linear equations to obtain a cooling network correction equation, a heating network correction equation and a power supply network correction equation;
and the second calculation unit is used for respectively solving the power supply network correction equation, the heat supply network correction equation and the cold supply network correction equation to obtain the first variable correction value, the second variable correction value and the third variable correction value.
8. The device of claim 6, wherein the obtaining module is further configured to obtain a flow balance equation and a pipeline pressure drop balance equation of a cooling network or a heating network in the combined cooling, heating and power system as hydraulic models; acquiring a linear equation set of water supply temperature and return water temperature of a cooling network or a heating network in the combined cooling heating and power system as a thermodynamic model; acquiring an active power flow equation and a reactive power flow equation of a power system in the combined cooling heating and power system as a power system model; obtaining a relational expression between consumed power and output power of a plurality of energy coupling elements in the combined cooling heating and power system as a coupling element model, wherein the plurality of energy coupling elements comprise a combined thermal power generation (CHP) unit, a refrigerator, an electric water pump and a heat pump.
9. A computer comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any of claims 1 to 5 are implemented when the computer program is executed by the processor.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910745485.XA CN110492533B (en) | 2019-08-13 | 2019-08-13 | Control method and device of combined cooling heating and power system, computer and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910745485.XA CN110492533B (en) | 2019-08-13 | 2019-08-13 | Control method and device of combined cooling heating and power system, computer and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110492533A CN110492533A (en) | 2019-11-22 |
| CN110492533B true CN110492533B (en) | 2021-10-08 |
Family
ID=68550817
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910745485.XA Active CN110492533B (en) | 2019-08-13 | 2019-08-13 | Control method and device of combined cooling heating and power system, computer and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110492533B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112906220B (en) * | 2021-02-10 | 2023-04-07 | 海南省电力学校(海南省电力技工学校) | Method for estimating state of comprehensive energy microgrid park system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106844895A (en) * | 2016-12-30 | 2017-06-13 | 华南理工大学 | A kind of decoupling calculation method of supply of cooling, heating and electrical powers microgrid energy stream |
| CN108494021A (en) * | 2018-04-20 | 2018-09-04 | 东北大学 | The Stability Assessment of electric-thermal-gas integrated energy system and static cost control method |
-
2019
- 2019-08-13 CN CN201910745485.XA patent/CN110492533B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106844895A (en) * | 2016-12-30 | 2017-06-13 | 华南理工大学 | A kind of decoupling calculation method of supply of cooling, heating and electrical powers microgrid energy stream |
| CN108494021A (en) * | 2018-04-20 | 2018-09-04 | 东北大学 | The Stability Assessment of electric-thermal-gas integrated energy system and static cost control method |
Non-Patent Citations (2)
| Title |
|---|
| "Fine-Grained Fully Parallel Power Flow Calculation by Incorporating BBDF Method Into a Multistep NR Algorithm";Xueneng Su 等;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20181130;第33卷(第6期);全文 * |
| "冷热电联供微网能量流的解耦计算方法";杨智斌 等;《电网技术》;20171231;第41卷(第12期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110492533A (en) | 2019-11-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mazloum et al. | Exergy analysis and exergoeconomic optimization of a constant-pressure adiabatic compressed air energy storage system | |
| CN111191182B (en) | Static combined heat and power flow hybrid solving method based on linearized thermal model | |
| CN107947182B (en) | Dynamic power flow analysis method and dynamic power flow analysis system | |
| CN111082417A (en) | State estimation method based on comprehensive energy system electric and heat combined network | |
| CN110991092A (en) | Complementary energy integrated utilization calculation method and system | |
| CN110492533B (en) | Control method and device of combined cooling heating and power system, computer and storage medium | |
| CN112070358A (en) | Method and system for determining electric load adjustment interval of low-vacuum heat supply unit | |
| Li et al. | Graph theory-based heat current analysis method for supercritical CO2 power generation system | |
| CN107887912B (en) | Static power flow analysis method and static power flow analysis system | |
| CN111049144B (en) | Active power distribution network load flow calculation PV node processing method based on constant change rate matrix | |
| CN111463795B (en) | Electric-thermal system multi-energy-flow analysis method based on improved Suhoff cooling formula | |
| CN111310343B (en) | Heat supply network hot circuit modeling method for comprehensive energy system scheduling | |
| CN103235512A (en) | Power generator set operating method | |
| CN211044253U (en) | Complementary energy integrated utilization computing system | |
| CN115062555A (en) | Integrated energy system exergy flow direct calculation method based on unbalanced node exergy | |
| CN113324283B (en) | Central heating network equivalence method considering temperature semi-dynamic characteristics | |
| CN111310310B (en) | Thermodynamic system static state trend rapid decoupling calculation method for quantity adjustment | |
| CN114221346A (en) | Load flow calculation method of comprehensive energy system | |
| CN111950135B (en) | Electric-gas interconnection system probability energy flow calculation method based on network decoupling | |
| CN111507011B (en) | Method for correcting influence quantity of heat supply extraction steam on sliding pressure operation of steam turbine | |
| CN114629124A (en) | Comprehensive energy system load flow calculation method based on subareas | |
| CN114118818B (en) | Method and device for estimating running state of distributed comprehensive energy system | |
| CN113221428B (en) | Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system | |
| CN108280302A (en) | A kind of heat supply network hot-fluid computational methods | |
| Plis et al. | Mathematical modelling of single pressure heat recovery steam generator |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210222 Address after: 510620, No. two, No. 2, Tianhe South Road, Guangzhou, Guangdong, Tianhe District Applicant after: Guangzhou Power Supply Bureau of Guangdong Power Grid Co.,Ltd. Address before: 510620, No. two, No. 2, Tianhe South Road, Guangzhou, Guangdong, Tianhe District Applicant before: GUANGZHOU POWER SUPPLY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |