A comprehensive energy system optimization scheduling method and system considering exergy loss

Abstract

The present invention relates to a method and system for optimizing and scheduling an integrated energy system taking losses into account, comprising: determining the optimal voltage of each energy supply node of the power supply network in the integrated energy system according to the losses of the integrated energy system; determining the optimal transformer ratio of each energy supply node of the power supply network in the integrated energy system according to the optimal voltage of each energy supply node of the power supply network in the integrated energy system; and adjusting the transformer ratio of each energy supply node of the power supply network in the integrated energy system to the optimal transformer ratio. The technical solution provided by the present invention, as a standard for measuring the energy supply capacity and quality of the system, proposes a method for calculating the losses of the power network and the heat supply network of the integrated energy system, and optimizes and schedules the integrated energy system with the goal of minimizing the losses of the integrated energy system, so as to maximize the external work of the integrated energy system as much as possible, realize the "energy-quality matching" of the energy used by the integrated energy system, and improve the energy utilization efficiency of the integrated energy system.

Description

Technical Field

[0001] The present invention relates to the field of energy and thermodynamics, and in particular to a method for calculating An optimized scheduling method and system for an integrated energy system with low energy consumption. Background Art

[0002] In the integrated energy system, energy of different forms and properties will inevitably be dissipated during transmission and conversion, thereby reducing the system's overall energy conversion efficiency. The losses incurred during energy transmission for users of different energy levels, such as those in cities, parks, districts, and towns, are also different. Therefore, it is necessary to analyze the energy losses in the process and formulate relevant network optimization plans to reduce the losses.

[0003] On the one hand, when a system reacts spontaneously, the direction of the reaction must change in the direction of reducing the total energy quality of the system, which is the principle of energy degradation. Traditional optimization scheduling methods are based on the premise that energy utilization efficiency and energy loss satisfy the first law of thermodynamics (because no matter what kind of energy conversion occurs, the amount of energy entering and leaving the multi-energy flow network must be conserved). However, this optimization scheduling method can only determine the trend of energy quantity changes during the transmission process, but cannot determine the trend of energy quality changes during the transmission process (the gradual decline in the system's energy quality is equivalent to the gradual decline in the system's ability to perform external work). The accuracy of this type of optimization scheduling method needs to be improved;

[0004] On the other hand, the integrated energy system covers energy flow sub-networks with different energy flow properties, such as gas networks, power networks, and heating networks. However, different systems have different evaluation systems for energy efficiency. In traditional optimization and scheduling methods, the optimization results of a single subsystem are mostly achieved at the expense of the benefits of other subsystems, and its applicability needs to be improved. Summary of the Invention

[0005] In view of the deficiencies in the prior art, the present invention aims to provide a The optimized dispatching method of integrated energy system with low energy consumption is based on As a standard for measuring the energy supply capacity and quality of the system, it is proposed to calculate the power network and heating network of the integrated energy system loss method and integrated energy system The integrated energy system is optimized and dispatched with the goal of minimizing losses, maximizing the external work of the integrated energy system as much as possible, achieving "energy-quality matching" of the energy consumption of the integrated energy system, and improving the energy utilization efficiency of the integrated energy system.

[0006] The purpose of the present invention is achieved by adopting the following technical solutions:

[0007] The present invention provides a The improvement of the method for optimizing and dispatching a comprehensive energy system with low energy consumption is that the method comprises:

[0008] According to the integrated energy system Determine the optimal voltage of each energy supply node in the power supply network of the integrated energy system based on the loss;

[0009] Determining an optimal transformer ratio for each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system;

[0010] The transformer ratio of each energy supply node of the power supply network in the integrated energy system is adjusted to the optimal transformer ratio.

[0011] Preferably, the integrated energy system The optimal voltage of each energy supply node in the power supply network of the integrated energy system is determined by the loss, including:

[0012] Integrated energy system The objective function of the integrated energy system optimization scheduling model is established with the goal of minimizing the loss.

[0013] The objective function of the integrated energy system optimization scheduling model is solved based on the constraint conditions corresponding to the objective function of the integrated energy system optimization scheduling model to obtain the optimal voltage of each energy supply node in the power supply network in the integrated energy system.

[0014] Furthermore, the objective function of the integrated energy system optimization scheduling model is determined as follows:

[0015]

[0016] Where C loss For integrated energy systems loss, c1 is the power of Quality coefficient, c2 is heat of Quality coefficient, w is pressure With electricity The conversion factor, is the power of the i-th transmission line in the power supply network in the integrated energy system. Loss, is the pressure of the kth heating pipe in the heating network in the integrated energy system Loss, The heat of the kth heating pipe in the heating network in the integrated energy system Loss, k∈[1~S hl ], S hlis the number of heating pipes in the heating network of the integrated energy system, i∈[1~S el ], S el The number of transmission lines supplying the electricity network in an integrated energy system.

[0017] Furthermore, the constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraints, operating output inequality constraints of generator sets, gas units, and circulating water pumps in the power supply network, tolerance constraints, transmission capacity constraints, transformer tap position constraints, and pipeline temperature constraints.

[0018] The present invention provides a The improvement of the integrated energy system optimization and dispatching system is that the system includes:

[0019] The first determination module is used to determine the Determine the optimal voltage of each energy supply node in the power supply network of the integrated energy system based on the loss;

[0020] A second determination module is used to determine the optimal transformer ratio of each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system;

[0021] The adjustment module is used to adjust the transformer ratio of each energy supply node in the power supply network in the integrated energy system to the optimal transformer ratio.

[0022] Preferably, the first determining module includes:

[0023] Building blocks for integrated energy systems The objective function of the integrated energy system optimization scheduling model is established with the goal of minimizing the loss.

[0024] An acquisition unit is used to solve the objective function of the integrated energy system optimization scheduling model based on the constraint conditions corresponding to the objective function of the integrated energy system optimization scheduling model, and obtain the optimal voltage of each energy supply node in the power supply network in the integrated energy system.

[0025] Furthermore, the objective function of the integrated energy system optimization scheduling model is determined as follows:

[0026]

[0027] Where C loss For integrated energy systems loss, c1 is the power of Quality coefficient, c2 is heat of Quality coefficient, w is pressure With electricity The conversion factor, is the power of the i-th transmission line in the power supply network in the integrated energy system. Loss, is the pressure of the kth heating pipe in the heating network in the integrated energy system Loss, The heat of the kth heating pipe in the heating network in the integrated energy system Loss, k∈[1~S hl ], S hl is the number of heating pipes in the heating network of the integrated energy system, i∈[1~S el ], S el The number of transmission lines supplying the electricity network in an integrated energy system.

[0028] Furthermore, the constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraints, operating output inequality constraints of generator sets, gas units, and circulating water pumps in the power supply network, tolerance constraints, transmission capacity constraints, transformer tap position constraints, and pipeline temperature constraints.

[0029] Compared with the closest prior art, the present invention has the following beneficial effects:

[0030] The technical solution provided by the present invention is based on the integrated energy system The optimal voltage of each energy supply node in the power supply network of the integrated energy system is determined based on the loss; the optimal transformer ratio of each energy supply node in the power supply network of the integrated energy system is determined based on the optimal voltage of each energy supply node in the power supply network of the integrated energy system; the transformer ratio of each energy supply node in the power supply network of the integrated energy system is adjusted to the optimal transformer ratio; the scheme will As a standard for measuring the system's energy supply capacity and quality, the integrated energy system The integrated energy system is optimized and dispatched with the goal of minimizing losses, maximizing the external work of the integrated energy system as much as possible, achieving "energy-quality matching" of the energy consumption of the integrated energy system, and improving the energy utilization efficiency of the integrated energy system. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Figure 1 It is a consideration Flowchart of the optimized dispatching method for integrated energy system with low energy consumption;

[0032] Figure 2 This is a topological diagram of a comprehensive energy system on an island according to an embodiment of the present invention;

[0033] Figure 3 It is an algorithm iteration curve diagram in an embodiment of the present invention;

[0034] Figure 4 2 is a graph showing voltage amplitude changes in cooling and heating modes of a comprehensive energy system on an island according to an embodiment of the present invention;

[0035] Figure 5 1 is a graph showing changes in pipeline mass flow in cooling and heating modes of a certain island integrated energy system according to an embodiment of the present invention;

[0036] FIG6( a ) is a diagram showing the changes in heating temperature of each node of a comprehensive energy system on an island before and after optimization in an embodiment of the present invention;

[0037] FIG6( b ) is a diagram showing the changes in cooling temperature of each node in a comprehensive energy system on an island before and after optimization in an embodiment of the present invention;

[0038] Figure 7 is a comparison chart of the cooling / heating capacity of a certain island integrated energy system according to an embodiment of the present invention;

[0039] Figure 8 This is the comprehensive energy system of an island before and after optimization in the embodiment of the present invention. Loss comparison chart;

[0040] Figure 9 It is a consideration Structure diagram of the integrated energy system optimization and dispatching system with low energy consumption. DETAILED DESCRIPTION

[0041] The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0042] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0043] The present invention provides a The optimal dispatching method of integrated energy system with low energy consumption, such as Figure 1 As shown, the method includes:

[0044] Step 101. Based on the integrated energy system Determine the optimal voltage of each energy supply node in the power supply network of the integrated energy system based on the loss;

[0045] Step 102: Determine the optimal transformer ratio of each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system;

[0046] Step 103: Adjust the transformer ratio of each energy supply node in the power supply network in the integrated energy system to the optimal transformer ratio.

[0047] Specifically, step 101 includes:

[0048] Step a: Based on the integrated energy system The objective function of the integrated energy system optimization scheduling model is established with the goal of minimizing the loss.

[0049] Step b: Solve the objective function of the integrated energy system optimization scheduling model based on the constraint conditions corresponding to the objective function of the integrated energy system optimization scheduling model to obtain the optimal voltage of each energy supply node in the power supply network in the integrated energy system.

[0050] Furthermore, the objective function of the integrated energy system optimization scheduling model is determined as follows:

[0051]

[0052] Where C loss For integrated energy systems loss, c1 is the power of Quality coefficient, c2 is heat of Quality coefficient, w is pressure With electricity The conversion factor, is the power of the i-th transmission line in the power supply network in the integrated energy system. Loss, is the pressure of the kth heating pipe in the heating network in the integrated energy system Loss, The heat of the kth heating pipe in the heating network in the integrated energy system Loss, k∈[1~S hl ], S hl is the number of heating pipes in the heating network of the integrated energy system, i∈[1~S el ], S el the number of transmission lines supplying the electricity network in the integrated energy system;

[0053] The power of the i-th transmission line in the power supply network of the integrated energy system is determined as follows: Loss

[0054]

[0055] Where, The power flowing from the xth transmission line connected to the head end node of the i-th transmission line in the power supply network in the integrated energy system into the i-th transmission line in the power supply network in the integrated energy system is is the power flowing from the i-th transmission line in the power supply network of the integrated energy system to the c-th transmission line connected to its terminal node. R i is the resistance of the i-th transmission line in the power supply network in the integrated energy system, P i,s is the electric power flowing into the head end node of the i-th transmission line in the power supply network of the integrated energy system, U i,s is the voltage of the node at the head end of the ith transmission line in the power supply network in the integrated energy system, x∈[1~S x ], S x is the sum of the transmission lines connected to the head end node of the i-th transmission line in the power supply network of the integrated energy system, c∈[1~S c ], S c is the sum of the transmission lines connected to the terminal node of the ith transmission line in the power supply network of the integrated energy system;

[0056] In the best embodiment of the present invention, assuming that the transmission pipeline is uniform in texture, when the first and last ends of the transmission pipeline are only connected

[0057] When there is a transmission pipeline, the voltage at the head and end nodes of the transmission line of the power supply network in the integrated energy system satisfies the following formula:

[0058]

[0059] Where U B is the voltage at the head end of the transmission line of the power supply network in the integrated energy system, U C is the terminal node voltage of the transmission line of the power supply network in the integrated energy system, R is the resistance of the transmission line of the power supply network in the integrated energy system, P B Input power to the head-end node of the transmission line of the power supply network in the integrated energy system;

[0060] Therefore, the electric power of the head and terminal nodes of the transmission line of the power supply network in the integrated energy system satisfies the following formula:

[0061]

[0062] Where, P C Output power to the terminal nodes of the transmission lines of the power supply network in the integrated energy system;

[0063] Therefore, the head and end nodes of the transmission lines of the power supply network in the integrated energy system are Satisfy the following formula:

[0064]

[0065] Where, Input power to the head end node of the transmission line of the power supply network in the integrated energy system Output electricity to the terminal node of the transmission line of the power supply network in the integrated energy system U0 is the quiescent voltage of the node in the power supply network of the integrated energy system, which is 0;

[0066] When the head and the end of the transmission pipeline are only connected to multiple transmission pipelines, the head and end node electrical Satisfy the following formula:

[0067]

[0068] Where, P Bz Input power to the head end node of the transmission line of the power supply network in the integrated energy system The sum, that is For the transmission lines of the power supply network in the integrated energy system Change.

[0069] The pressure of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Loss

[0070]

[0071] Where R k is the fluid flow resistance of the kth heating pipe in the heating network in the integrated energy system, x k is the volume flow of the kth heating pipeline in the heating network of the integrated energy system;

[0072] In the preferred embodiment of the present invention, assuming that the transmission pipeline has a uniform texture, when the head end and the tail end of the transmission pipeline are each connected to only one transmission pipeline, the head-end node voltage of the heating pipeline of the heating network in the integrated energy system satisfies the following formula:

[0073] p C =p B -R y x Q

[0074] Where p Bis the node pressure at the head end of the heating pipe in the heating network of the integrated energy system, p C is the terminal node pressure of the heating pipeline in the integrated energy system, R y is the fluid flow resistance of the heating pipe in the integrated energy system, x Q is the volume flow of the heating pipeline in the integrated energy system;

[0075] Therefore, the pressure power of the head and end nodes of the heating pipes in the heating network of the integrated energy system satisfies the following formula:

[0076] P Cy =P By -R y x Q 2

[0077] Where, P By P is the node pressure power at the head end of the heating pipeline in the heating network of the integrated energy system; Cy The power output of the terminal node of the heating pipeline of the heating network in the integrated energy system;

[0078] Therefore, the pressure at the head and end nodes of the heating pipes in the heating network of the integrated energy system is Satisfy the following formula:

[0079]

[0080] Where, The node pressure of the heating pipe at the head end of the heating network in the integrated energy system P Cy The terminal node pressure of the heating pipeline in the heating network of the integrated energy system

[0081] When the head end and the tail end of the transmission pipeline are only connected to multiple transmission pipelines, the node pressure at the head end of the heating pipeline in the heating network of the integrated energy system is Satisfy the following formula:

[0082]

[0083] Where x Qz is the total volume flow input of the heating pipes of the heating network in the integrated energy system, i.e. The pressure of heating pipes in the heating network of the integrated energy system Change.

[0084] The heat of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Loss

[0085]

[0086] Where, The heat from the zth heating pipeline connected to the head end node of the kth heating pipeline in the heating network in the integrated energy system flows into the kth heating pipeline in the heating network in the integrated energy system. is the heat flow from the kth heating pipe in the heating network of the integrated energy system to the bth heating pipe connected to its terminal node. ρ is the density of the working medium transmitted by the heating pipeline, c a is the specific heat capacity of the working medium transmitted by the heating pipeline, is the temperature of the head end node of the kth heating pipe in the heating network of the integrated energy system, X T0 is the static value of the heating network temperature in the integrated energy system, ε k is the composite heat transfer coefficient of the kth heating pipe in the heating network of the integrated energy system, z∈[1~S z ], S z is the number of heating pipes connected to the head end node of the kth heating pipe in the heating network of the integrated energy system, b∈[1~S b ], S b is the number of heating pipelines connected to the end node of the kth heating pipeline in the heating network of the integrated energy system.

[0087] In the preferred embodiment of the present invention, assuming that the transmission pipeline has a uniform texture, when the head end and the tail end of the transmission pipeline are each connected to only one transmission pipeline, the head and tail node temperatures of the heating pipeline in the heating network of the integrated energy system satisfy the following formula:

[0088]

[0089] is the heat transfer coefficient of the insulated pipe with n layers of insulation material, τ i ,i∈n is the heat transfer coefficient of each layer of insulation material from the inside to the outside, d i , i∈n is the radius of each layer of material; x Q is the volume flow of the heating pipe of the heating network in the integrated energy system; ρ, c are the density and specific heat capacity of the hot / cold working medium respectively, X T0 is the static intensity value of thermal energy.

[0090] Therefore, the heat supply network in the integrated energy system is the first and last node of the heating pipeline. Satisfy the following formula:

[0091]

[0092] Where, The first end node heat supply of the heating pipe of the heating network in the integrated energy system The terminal node heat supply of the heating pipeline in the heating network of the integrated energy system

[0093] When the head and the end of the transmission pipeline are only connected to multiple transmission pipelines, the heating network in the integrated energy system

[0094] The node pressure at the head end of the heating pipe Satisfy the following formula:

[0095]

[0096] Where x Qz is the total volume flow input of the heating pipes of the heating network in the integrated energy system, i.e. The pressure of heating pipes in the heating network of the integrated energy system Change.

[0097] Further, the electric of Quality coefficient c1:

[0098]

[0099] Where, φ1 is the energy level evaluation factor of electric energy, μ1 is the energy level factor of electric energy;

[0100] Determine the heat by pressing the formula of Quality coefficient c2:

[0101]

[0102] Where, φ2 is the energy level evaluation factor of thermal energy, and μ2 is the energy level factor of thermal energy;

[0103] The resistance R of the i-th transmission line in the power supply network of the integrated energy system is determined by the following formula: i :

[0104]

[0105] Where, L i is the length of the i-th transmission line in the power supply network in the integrated energy system, A i is the transmission channel area of ​​the i-th transmission line in the power supply network of the integrated energy system, K i is the resistivity of the conductor material of the i-th transmission line of the power supply network in the integrated energy system;

[0106] The fluid flow resistance R of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k :

[0107]

[0108] Where μ k is the fluid viscosity of the kth heating pipe in the heating network of the integrated energy system, L k is the length of the kth heating pipeline in the heating network of the integrated energy system, r k The inner radius of the kth heating pipe in the heating network of the integrated energy system;

[0109] The volume flow rate x of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k :

[0110]

[0111] Where p k,s is the pressure at the head node of the kth heating pipeline in the heating network of the integrated energy system, p k,m is the pressure at the end node of the kth heating pipeline in the heating network of the integrated energy system;

[0112] The composite heat transfer coefficient ε of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k :

[0113]

[0114] Where, d k(σ+1) is the radius of the σ+1th layer of material of the kth heating pipe in the heating network of the integrated energy system, d kσ is the radius of the σ-th layer of material of the k-th heating pipe in the heating network of the integrated energy system, γ kσ is the heat transfer coefficient of the σth layer of the material of the kth heating pipe in the heating network of the integrated energy system, σ∈(1~Ψ-1), Ψ is the number of layers of the heat transfer material of the kth heating pipe in the heating network of the integrated energy system;

[0115] Furthermore, the constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraints, operating output inequality constraints of generator sets, gas units, and circulating water pumps in the power supply network, tolerance constraints, transmission capacity constraints, transformer tap position constraints, and pipeline temperature constraints.

[0116] Furthermore, the equality constraint condition is determined as follows:

[0117]

[0118] Where, is the input power of the qth power supply node in the power supply network of the integrated energy system, is the electricity demand of the fth load node in the power supply network in the integrated energy system, is the current of the i-th transmission line in the power supply network in the integrated energy system, is the voltage variation of the i-th transmission line in the power supply network of the integrated energy system, is the input thermal power of the dth heating node of the heating network in the integrated energy system, is the heat demand of the rth load node in the heating network in the integrated energy system, is the entropy of the kth heating pipe in the heating network in the integrated energy system, is the temperature variation of the kth heating pipe in the heating network in the integrated energy system, q∈[1~S q ], S q is the total number of power supply nodes in the power supply network of the integrated energy system, f∈[1~S f ], S f is the total number of load nodes in the power supply network of the integrated energy system, d∈[1~S d ], S d is the total number of heating nodes in the heating network of the integrated energy system, r∈[1~S r ] is the total number of load nodes in the heating network of the integrated energy system;

[0119] The operating output inequality constraints of the generator set, gas unit, and circulating water pump in the power supply network are determined as follows:

[0120]

[0121] Where, P FD P is the operating output of the generator set in the power supply network of the integrated energy system. FD,max is the upper limit of the operating output of the generator set in the power supply network of the integrated energy system, P FD,min is the lower limit of the operating output of the generator set in the power supply network of the integrated energy system, P RQ The operating output of the gas-fired unit in the heating network of the integrated energy system, P RQ,max is the upper limit of the operating output of the gas-fired units in the heating network of the integrated energy system, P RQ,min is the lower limit of the operating output of the gas-fired units in the heating network of the integrated energy system, P XH is the operating output of the circulating water pump in the heating network of the integrated energy system, P XH,max is the upper limit of the operating output of the circulating water pump in the heating network of the integrated energy system, P XH,minIt is the lower limit of the operating output of the circulating water pump in the heating network of the integrated energy system;

[0122] The tolerance constraint condition is determined as follows:

[0123]

[0124] Where, is the voltage of the hth node in the power supply network of the integrated energy system, is the maximum tolerance value of the hth node in the power supply network of the integrated energy system, is the minimum tolerance value of the hth node in the power supply network of the integrated energy system, is the temperature of the Jth node in the heating network of the integrated energy system, is the maximum temperature tolerance of the J-th node in the heating network of the integrated energy system, is the minimum temperature tolerance value of the J-th node in the heating network of the integrated energy system, is the working fluid pressure of the kth heating pipe in the heating network of the integrated energy system, is the maximum tolerable pressure of the working fluid in the kth heating pipe in the heating network of the integrated energy system, is the minimum tolerable pressure of the working fluid in the kth heating pipe in the heating network of the integrated energy system, is the gas pressure of the t-th gas unit in the heating network of the integrated energy system, is the maximum tolerable gas pressure of the t-th gas unit in the heating network of the integrated energy system, is the minimum pressure tolerance value of the t-th gas unit in the heating network of the integrated energy system, J∈(1~S J ), S J is the total number of nodes in the heating network of the integrated energy system, h∈(1~S h ), S h is the total number of nodes in the power supply network of the integrated energy system, t∈(1~S t ), S t is the total number of gas units in the heating network of the integrated energy system;

[0125] The transmission capacity constraint condition is determined as follows:

[0126]

[0127] Where, is the current transmission capacity of the i-th transmission line of the power supply network in the integrated energy system, is the minimum current transmission capacity of the i-th transmission line in the power supply network of the integrated energy system, is the maximum current transmission capacity of the i-th transmission line in the power supply network in the integrated energy system, is the volume flow transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the maximum volume flow transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the minimum volume flow transmission capacity of the kth heating pipeline in the heating network of the integrated energy system, is the heat energy transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the maximum heat energy transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the minimum heat energy transmission capacity of the kth heating pipeline in the heating network of the integrated energy system;

[0128] The transformer tap position constraint condition is determined as follows:

[0129]

[0130] Where, Tap δ is the δth transformer tap position of the power supply network in the integrated energy system, The adjustable lower limit value of the δth transformer tap position in the power supply network in the integrated energy system. is the adjustable upper limit value of the δth transformer tap position in the power supply network in the integrated energy system, δ∈[1~S Tap ], S Tap Number of transformer taps for the power supply network in the integrated energy system;

[0131] The pipeline temperature constraint condition is determined as follows:

[0132]

[0133] Where C s is the first coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, C r is the first coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system, X Ts is the temperature matrix of the load heat flow inlet, X Tr is the temperature matrix of the load heat flow outlet, X T0 is the temperature static value matrix, b s is the second coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, b r It is the second coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system.

[0134] Furthermore, the first coefficient matrix C of the node temperature between the heating pipes of the heating network in the integrated energy system is determined as follows: s :

[0135]

[0136] Where, is the node υ in the heating network of the integrated energy system a Temperature and Node γ a The first coefficient between temperatures, υ a , γ a ∈(γ N ), γ N It is the set of nodes between heating pipes in the heating network of the integrated energy system;

[0137] The second coefficient matrix b of the node temperature between the heating pipes of the heating network in the integrated energy system is determined by the following formula: s :

[0138]

[0139] Where, is the node υ in the heating network of the integrated energy system a The second coefficient of temperature, T is the transposed sign;

[0140] Among them, when the heat medium is supplied from node γ a Through pipe M a Flow to node υ a hour:

[0141] If the node υ a The node γ a The supplied heat medium and node γ a is a load node, then:

[0142]

[0143]

[0144] If the node υ a The node γ a The supplied heat medium and node γ a is not a load node, then:

[0145]

[0146]

[0147]

[0148] If the node υ aThe node γ a The supplied heat medium is

[0149]

[0150]

[0151] In the above formula, The heating pipe M of the heating network in the integrated energy system a The volume flow rate, The heating pipe M of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe M of the heating network in the integrated energy system a length, is the node υ in the heating network of the integrated energy system a The total number of heating pipes supplying hot working medium, X T0 is the quiescent value of temperature, Node γ in the heating network of the integrated energy system a temperature, is the node υ in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature;

[0152] The first coefficient matrix C of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r :

[0153]

[0154] Where, is the node η in the heating network of the integrated energy system a Temperature and node φ a The first coefficient between temperatures, φ a ,η a ∈(1~χ N ), χ N is the total number of nodes between return pipes in the heating network of the integrated energy system;

[0155] The second coefficient matrix b of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r

[0156]

[0157] Where, is the node η in the heating network of the integrated energy system a The second coefficient of temperature;

[0158] Among them, when the heat medium is supplied from node φ a Through pipe H a Flow to node η a hour:

[0159] If the node η a Not only by node φ a When the hot working medium is supplied, then:

[0160]

[0161]

[0162]

[0163] If the node η a Only by node φ a When hot working medium is supplied, then:

[0164]

[0165]

[0166] In the above formula, is the node φ in the heating network of the integrated energy system a The outlet temperature, is the node φ in the heating network of the integrated energy system a The volume flow rate, The heating pipe H of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe H of the heating network in the integrated energy system a length, is the node η in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature, is the node η in the heating network of the integrated energy system a The total number of return pipes supplying hot working medium.

[0167] Specifically, step 102 includes:

[0168] The optimal transformer ratio of the hth energy supply node in the power supply network of the integrated energy system is determined by the following formula:

[0169]

[0170] Where U eh,b is the voltage of the h-th energy supply node connected to the bus in the power supply network of the integrated energy system, is the optimal voltage of the hth energy supply node in the power supply network of the integrated energy system, h∈(1~S h ), S h is the total number of nodes in the power supply network of the integrated energy system.

[0171] In a specific embodiment of the present invention, an integrated energy system for providing electricity and heat to an island is analyzed. Figure 2 The topology diagram of an island's integrated energy system is shown in Table 1, where the node loads for the power network and the heating network are shown in Table 2. The power network and the heating network are interconnected through three CHP units and circulating water pumps to achieve energy coupling.

[0172] Table 1

[0173] Power network node number 1 2 3 4 5 6 Load / (MW) 0.2 0 0.5 0.5 0.2 0.2

[0174] Table 2

[0175]

[0176] exist Figure 2 In this example, units G1, G2, and G3 are gas turbines, extraction steam turbines, and reciprocating internal combustion engines, respectively. The power network consists of nine buses and five loads, and the heating network consists of 32 nodes and 32 pipelines. The state variables and control variables of the example are shown in Tables 3 and 4, respectively.

[0177] Table 3

[0178]

[0179]

[0180] Table 4

[0181]

[0182] In this example, the dimension of the particle swarm algorithm is 7, the learning factors c1 and c2 of the particle swarm algorithm are set to 2, the inertia factor w is set to 1, the number of iterations of the algorithm is 60, and the number of particles is 40.

[0183] Algorithm iteration: In the example, the fitness function is selected as the system line transmission caused by The inverse of the loss can be used to obtain the convergence characteristics by depicting the best fitness value of each generation in the population. Figure 3 Iteration curves of the particle swarm algorithm for 60 iterations in cooling and heating modes.

[0184] Depend on Figure 3It can be seen that when the particle swarm algorithm is used for calculation, the fitness value of each generation of particles is shown as the small solid points. The small solid points are connected by lines to form the algorithm convergence characteristic curve, which reflects the convergence trend of the algorithm. The fitness curve of this example basically reaches the optimal value after more than 30 iterations, so it has strong convergence ability and fast convergence speed.

[0185] Results comparison: The particle swarm algorithm can obtain the optimal value of the objective function and the optimized value of the control variable. Under the optimal particle result, the voltage amplitude of the integrated energy system in the cooling and heating modes is as follows: Figure 4 As shown in the figure, the pipeline mass flow rate in the cooling and heating modes of the integrated energy system is as follows Figure 5 As shown in Figure 6(a), the changes in the heating temperature of each node in the integrated energy system before and after optimization are shown in Figure 6(b);

[0186] Figure 4 It can be seen that the voltage amplitude of each node in the power system has been improved after the optimization calculated by the particle swarm algorithm. The voltage amplitude of each power node has decreased significantly after optimization compared with that before optimization, which is helpful to reduce the system power. From the data analysis, the optimization results of the power network for the regional integrated energy system in the heating and cooling modes are very similar.

[0187] Figure 5 It can be seen that in the heating mode, the mass flow rate in the heating system pipes after optimization by the particle swarm algorithm has decreased to a certain extent, but the decrease is small. The mass flow rate in the cooling system pipes after optimization by the particle swarm algorithm has decreased significantly. Therefore, the mass flow rate of the regional integrated energy system in the cooling / heating mode has decreased, thereby reducing the system pressure. It reduces losses and also helps improve system security.

[0188] As can be seen from Figure 6(a) and Figure 6(b), the node temperature curves of the two modes are basically consistent before and after optimization. However, since the mass flow rate of each node has decreased after optimization, in order to ensure that the system's cooling / heating load requirements are met, the quality of the system's energy supply must be improved. Therefore, after optimization, the temperature of most of the system's thermal nodes in the heating mode has increased slightly compared to before optimization, while the temperature of each thermal node in the cooling mode has decreased by an average of 2°C compared to before optimization. In different modes, by adjusting the temperature of the thermal network nodes, the system is guaranteed to be The total amount is stable.

[0189] From the optimization results of the cooling and heating modes of the above regional integrated energy system, it can be seen that under the same power conditions, the mass flow rate required for the transmission of cold energy in the cooling mode is much greater than the mass flow rate required for the transmission of heat energy in the heating mode. This is because under certain ambient temperature conditions, the cooling temperature of the cold source node is generally a few degrees Celsius to more than ten degrees Celsius, while the heating temperature of the heat source node is generally more than 60 degrees Celsius or even more than 90 degrees Celsius, which results in a higher cooling rate per unit mass flow rate. Generally less than hot That is, the work capacity of cold water is generally less than that of hot water. Figure 7 As shown in the figure, if the source node temperature remains constant, when the cooling / heating demand gradually increases, the increase in mass flow rate caused by the increased cooling demand in the cooling mode is much greater than the increase in mass flow rate caused by the increased heating demand in the heating mode. In this case, when the pipeline capacity is fixed, the cooling capacity of the cooling system is very limited, which restricts the cooling range of the cooling system.

[0190] It can be seen that the use of After the loss-based integrated energy system operation optimization method is applied, the mass flow rate of the system in different modes has decreased, which effectively increases the safety margin of the system and enhances the safety and reliability of the system. When the system parameters do not change much, such as Figure 8 As shown in the figure, under the heating mode of the original island integrated energy system, The loss is 0.0306MW. After optimization The loss is 0.0157MW, and the total system The loss reduction ratio reached 48.91%; in cooling mode, The loss is 0.0534MW. After optimization The loss is 0.0349MW, and the total system The loss reduction ratio reached 34.59%, which effectively improved the energy utilization efficiency of the integrated energy system and realized the "energy-quality matching" principle of the system, proving the effectiveness of the proposed method.

[0191] The present invention provides a Loss-based integrated energy system optimization dispatching system, such as Figure 9 As shown, the system includes:

[0192] The first determination module is used to determine the Determine the optimal voltage of each energy supply node in the power supply network of the integrated energy system based on the loss;

[0193] A second determination module is used to determine the optimal transformer ratio of each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system;

[0194] The adjustment module is used to adjust the transformer ratio of each energy supply node in the power supply network in the integrated energy system to the optimal transformer ratio.

[0195] Specifically, the first determining module includes:

[0196] Establishment of units for integrated energy systems The objective function of the integrated energy system optimization scheduling model is established with the goal of minimizing the loss.

[0197] A solving unit is used to solve the objective function of the integrated energy system optimization scheduling model based on the constraint conditions corresponding to the objective function of the integrated energy system optimization scheduling model, and obtain the optimal voltage of each energy supply node in the power supply network in the integrated energy system.

[0198] Furthermore, the objective function of the integrated energy system optimization scheduling model is determined as follows:

[0199]

[0200] Where C loss For integrated energy systems loss, c1 is the power of Quality coefficient, c2 is heat of Quality coefficient, w is pressure With electricity The conversion factor, is the power of the i-th transmission line in the power supply network in the integrated energy system. Loss, is the pressure of the kth heating pipe in the heating network in the integrated energy system Loss, The heat of the kth heating pipe in the heating network in the integrated energy system Loss, k∈[1~S hl ], S hl is the number of heating pipes in the heating network of the integrated energy system, i∈[1~S el ], S el the number of transmission lines supplying the electricity network in the integrated energy system;

[0201] The power of the i-th transmission line in the power supply network of the integrated energy system is determined as follows: Loss

[0202]

[0203] Where, The power flowing from the xth transmission line connected to the head end node of the i-th transmission line in the power supply network in the integrated energy system into the i-th transmission line in the power supply network in the integrated energy system is is the power flowing from the i-th transmission line in the power supply network of the integrated energy system to the c-th transmission line connected to its terminal node. R i is the resistance of the i-th transmission line in the power supply network in the integrated energy system, P i,s is the electric power flowing into the head end node of the i-th transmission line in the power supply network of the integrated energy system, U i,s is the voltage of the node at the head end of the ith transmission line in the power supply network in the integrated energy system, x∈[1~S x ], S x is the sum of the transmission lines connected to the head end node of the i-th transmission line in the power supply network of the integrated energy system, c∈[1~S c ], S c is the sum of the transmission lines connected to the terminal node of the ith transmission line in the power supply network of the integrated energy system;

[0204] The pressure of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Loss

[0205]

[0206] Where R k is the fluid flow resistance of the kth heating pipe in the heating network in the integrated energy system, x k is the volume flow of the kth heating pipeline in the heating network of the integrated energy system;

[0207] The heat of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Loss

[0208]

[0209] Where, The heat from the zth heating pipeline connected to the head end node of the kth heating pipeline in the heating network in the integrated energy system flows into the kth heating pipeline in the heating network in the integrated energy system. is the heat flow from the kth heating pipe in the heating network of the integrated energy system to the bth heating pipe connected to its terminal node. ρ is the density of the working medium transmitted by the heating pipeline, c a is the specific heat capacity of the working medium transmitted by the heating pipeline, is the temperature of the head end node of the kth heating pipe in the heating network of the integrated energy system, X T0 is the static value of the heating network temperature in the integrated energy system, ε k is the composite heat transfer coefficient of the kth heating pipe in the heating network of the integrated energy system, z∈[1~S z ], S z is the number of heating pipes connected to the head end node of the kth heating pipe in the heating network of the integrated energy system, b∈[1~S b ], S b is the number of heating pipelines connected to the end node of the kth heating pipeline in the heating network of the integrated energy system.

[0210] Further, the electric of Quality coefficient c1:

[0211]

[0212] Where, φ1 is the energy level evaluation factor of electric energy, μ1 is the energy level factor of electric energy;

[0213] Determine the heat by pressing the formula of Quality coefficient c2:

[0214]

[0215] Where, φ2 is the energy level evaluation factor of thermal energy, and μ2 is the energy level factor of thermal energy;

[0216] The resistance R of the i-th transmission line in the power supply network of the integrated energy system is determined by the following formula: i :

[0217]

[0218] Where, L i is the length of the i-th transmission line in the power supply network in the integrated energy system, A i is the transmission channel area of ​​the i-th transmission line in the power supply network of the integrated energy system, K i is the resistivity of the conductor material of the i-th transmission line of the power supply network in the integrated energy system;

[0219] The fluid flow resistance R of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k :

[0220]

[0221] Where μ k is the fluid viscosity of the kth heating pipe in the heating network of the integrated energy system, Lk is the length of the kth heating pipeline in the heating network of the integrated energy system, r k The inner radius of the kth heating pipe in the heating network of the integrated energy system;

[0222] The volume flow rate x of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k :

[0223]

[0224] Where p k,s is the pressure at the head node of the kth heating pipeline in the heating network of the integrated energy system, p k,m is the pressure at the end node of the kth heating pipeline in the heating network of the integrated energy system;

[0225] The composite heat transfer coefficient ε of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k :

[0226]

[0227] Where, d k(σ+1) is the radius of the σ+1th layer of material of the kth heating pipe in the heating network of the integrated energy system, d kσ is the radius of the σ-th layer of material of the k-th heating pipe in the heating network of the integrated energy system, γ kσ is the heat transfer coefficient of the σth layer material of the kth heating pipe in the heating network of the integrated energy system, σ∈(1~Ψ-1), Ψ is the number of layers of heat transfer material of the kth heating pipe in the heating network of the integrated energy system.

[0228] Furthermore, the constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraints, operating output inequality constraints of generator sets, gas units, and circulating water pumps in the power supply network, tolerance constraints, transmission capacity constraints, transformer tap position constraints, and pipeline temperature constraints.

[0229] Furthermore, the equality constraint condition is determined as follows:

[0230]

[0231] Where, is the input power of the qth power supply node in the power supply network of the integrated energy system, is the electricity demand of the fth load node in the power supply network in the integrated energy system, is the current of the i-th transmission line in the power supply network in the integrated energy system, is the voltage variation of the i-th transmission line in the power supply network of the integrated energy system, is the input thermal power of the dth heating node of the heating network in the integrated energy system, is the heat demand of the rth load node in the heating network in the integrated energy system, is the entropy of the kth heating pipe in the heating network in the integrated energy system, is the temperature variation of the kth heating pipe in the heating network in the integrated energy system, q∈[1~S q ], S q is the total number of power supply nodes in the power supply network of the integrated energy system, f∈[1~S f ], S f is the total number of load nodes in the power supply network of the integrated energy system, d∈[1~S d ], S d is the total number of heating nodes in the heating network of the integrated energy system, r∈[1~S r ] is the total number of load nodes in the heating network of the integrated energy system;

[0232] The operating output inequality constraints of the generator set, gas unit, and circulating water pump in the power supply network are determined as follows:

[0233]

[0234] Where, P FD P is the operating output of the generator set in the power supply network of the integrated energy system. FD,max is the upper limit of the operating output of the generator set in the power supply network of the integrated energy system, P FD,min is the lower limit of the operating output of the generator set in the power supply network of the integrated energy system, P RQ The operating output of the gas-fired unit in the heating network of the integrated energy system, P RQ,max is the upper limit of the operating output of the gas-fired units in the heating network of the integrated energy system, P RQ,min is the lower limit of the operating output of the gas-fired units in the heating network of the integrated energy system, P XH is the operating output of the circulating water pump in the heating network of the integrated energy system, P XH,max is the upper limit of the operating output of the circulating water pump in the heating network of the integrated energy system, P XH,min It is the lower limit of the operating output of the circulating water pump in the heating network of the integrated energy system;

[0235] The tolerance constraint condition is determined as follows:

[0236]

[0237] Where, is the voltage of the hth node in the power supply network of the integrated energy system, is the maximum tolerance value of the hth node in the power supply network of the integrated energy system, is the minimum tolerance value of the hth node in the power supply network of the integrated energy system, is the temperature of the Jth node in the heating network of the integrated energy system, is the maximum temperature tolerance of the J-th node in the heating network of the integrated energy system, is the minimum temperature tolerance value of the J-th node in the heating network of the integrated energy system, is the working fluid pressure of the kth heating pipe in the heating network of the integrated energy system, is the maximum tolerable pressure of the working fluid in the kth heating pipe in the heating network of the integrated energy system, is the minimum tolerable pressure of the working fluid in the kth heating pipe in the heating network of the integrated energy system, is the gas pressure of the t-th gas unit in the heating network of the integrated energy system, is the maximum tolerable gas pressure of the t-th gas unit in the heating network of the integrated energy system, is the minimum pressure tolerance value of the t-th gas unit in the heating network of the integrated energy system, J∈(1~S J ), S J is the total number of nodes in the heating network of the integrated energy system, h∈(1~S h ), S h is the total number of nodes in the power supply network of the integrated energy system, t∈(1~S t ), S t is the total number of gas units in the heating network of the integrated energy system;

[0238] The transmission capacity constraint condition is determined as follows:

[0239]

[0240] Where, is the current transmission capacity of the i-th transmission line of the power supply network in the integrated energy system, is the minimum current transmission capacity of the i-th transmission line in the power supply network of the integrated energy system, is the maximum current transmission capacity of the i-th transmission line in the power supply network in the integrated energy system, is the volume flow transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the maximum volume flow transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the minimum volume flow transmission capacity of the kth heating pipeline in the heating network of the integrated energy system, is the heat energy transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the maximum heat energy transmission capacity of the kth heating pipeline in the heating network in the integrated energy system, is the minimum heat energy transmission capacity of the kth heating pipeline in the heating network of the integrated energy system;

[0241] The transformer tap position constraint condition is determined as follows:

[0242]

[0243] Where, Tap δ is the δth transformer tap position of the power supply network in the integrated energy system, The adjustable lower limit value of the δth transformer tap position in the power supply network in the integrated energy system. is the adjustable upper limit value of the δth transformer tap position in the power supply network in the integrated energy system, δ∈[1~S Tap ], S Tap Number of transformer taps for the power supply network in the integrated energy system;

[0244] The pipeline temperature constraint condition is determined as follows:

[0245]

[0246] Where C s is the first coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, C r is the first coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system, X Ts is the temperature matrix of the load heat flow inlet, X Tr is the temperature matrix of the load heat flow outlet, X T0 is the temperature static value matrix, b s is the second coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, b r It is the second coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system.

[0247] Furthermore, the first coefficient matrix C of the node temperature between the heating pipes of the heating network in the integrated energy system is determined as follows: s :

[0248]

[0249] Where, is the node υ in the heating network of the integrated energy system a Temperature and Node γ aThe first coefficient between temperatures, υ a , γ a ∈(γ N ), γ N It is the set of nodes between heating pipes in the heating network of the integrated energy system;

[0250] The second coefficient matrix b of the node temperature between the heating pipes of the heating network in the integrated energy system is determined by the following formula: s :

[0251]

[0252] Where, is the node υ in the heating network of the integrated energy system a The second coefficient of temperature, T is the transposed sign;

[0253] Among them, when the heat medium is supplied from node γ a Through pipe M a Flow to node υ a hour:

[0254] If the node υ a The node γ a The supplied heat medium and node γ a is a load node, then:

[0255]

[0256]

[0257] If the node υ a The node γ a The supplied heat medium and node γ a is not a load node, then:

[0258]

[0259]

[0260]

[0261] If the node υ a The node γ a The supplied heat medium is

[0262]

[0263]

[0264] In the above formula, The heating pipe M of the heating network in the integrated energy system aThe volume flow rate, The heating pipe M of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe M of the heating network in the integrated energy system a length, is the node υ in the heating network of the integrated energy system a The total number of heating pipes supplying hot working medium, X T0 is the quiescent value of temperature, Node γ in the heating network of the integrated energy system a temperature, is the node υ in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature;

[0265] The first coefficient matrix C of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r :

[0266]

[0267] Where, is the node η in the heating network of the integrated energy system a Temperature and node φ a The first coefficient between temperatures, φ a ,η a ∈(1~χ N ), χ N is the total number of nodes between return pipes in the heating network of the integrated energy system;

[0268] The second coefficient matrix b of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r

[0269]

[0270] Where, is the node η in the heating network of the integrated energy system a The second coefficient of temperature;

[0271] Among them, when the heat medium is supplied from node φ a Through pipe H a Flow to node η a hour:

[0272] If the node η a Not only by node φ a When the hot working medium is supplied, then:

[0273]

[0274]

[0275]

[0276] If the node η a Only by node φ a When hot working medium is supplied, then:

[0277]

[0278]

[0279] In the above formula, is the node φ in the heating network of the integrated energy system a The outlet temperature, is the node φ in the heating network of the integrated energy system a The volume flow rate, The heating pipe H of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe H of the heating network in the integrated energy system a length, is the node η in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature, is the node η in the heating network of the integrated energy system a The total number of return pipes supplying hot working medium.

[0280] Specifically, the second determining module is used to:

[0281] The optimal transformer ratio of the hth energy supply node in the power supply network of the integrated energy system is determined by the following formula:

[0282]

[0283] Where U eh,b is the voltage of the h-th energy supply node in the power supply network of the integrated energy system connected to the bus, #imgpt338# is the optimal voltage of the h-th energy supply node in the power supply network of the integrated energy system, h∈(1~S h ), S h is the total number of nodes in the power supply network of the integrated energy system.

[0284] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0285] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the steps in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0286] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0287] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0288] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered by the scope of protection of the claims of the present invention.

Claims

1. A method for optimizing and dispatching an integrated energy system taking into account exergy loss, characterized in that: The method comprises: Determine the optimal voltage of each energy supply node in the power supply network of the integrated energy system according to the exergy loss of the integrated energy system; Determining an optimal transformer ratio for each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system; Adjusting the transformer ratio of each energy supply node in the power supply network of the integrated energy system to the optimal transformer ratio; The determining of the optimal voltage of each energy supply node of the power supply network in the integrated energy system according to the exergy loss of the integrated energy system includes: The objective function of the integrated energy system optimization scheduling model is established with the goal of minimizing the exergy loss of the integrated energy system; Solving the objective function of the integrated energy system optimization scheduling model based on the constraint conditions corresponding to the objective function of the integrated energy system optimization scheduling model to obtain the optimal voltage of each energy supply node of the power supply network in the integrated energy system; The objective function of the integrated energy system optimization scheduling model is determined as follows: Where C loss is the exergy loss of the integrated energy system, c1 is the exergy mass coefficient of electrical exergy, c2 is the exergy mass coefficient of thermal exergy, w is the conversion coefficient between pressure exergy and electrical exergy, is the exergy loss of the i-th transmission line in the power supply network of the integrated energy system, is the pressure exergy loss of the kth heating pipeline in the heating network of the integrated energy system, is the heat exergy loss of the kth heating pipe in the heating network of the integrated energy system, k∈[1~S hl ], S hl is the number of heating pipes in the heating network of the integrated energy system, i∈[1~S el ], S el The number of transmission lines supplying the electricity network in an integrated energy system.

2. The method according to claim 1, wherein The exergy loss of the i-th transmission line in the power supply network of the integrated energy system is determined by the following formula: Where, is the exergy flowing from the xth transmission line connected to the head end node of the i-th transmission line in the power supply network of the integrated energy system into the i-th transmission line in the power supply network of the integrated energy system, is the exergy flowing from the i-th transmission line in the power supply network of the integrated energy system to the c-th transmission line connected to its terminal node, R i is the resistance of the i-th transmission line in the power supply network in the integrated energy system, P i,s is the electric power flowing into the head end node of the i-th transmission line in the power supply network of the integrated energy system, U i,s is the voltage of the node at the head end of the ith transmission line in the power supply network in the integrated energy system, x∈[1~S x ], S x is the sum of the transmission lines connected to the head end node of the i-th transmission line in the power supply network of the integrated energy system, c∈[1~S c ], S c is the sum of the transmission lines connected to the terminal node of the ith transmission line in the power supply network of the integrated energy system; The pressure exergy loss of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Where R k is the fluid flow resistance of the kth heating pipe in the heating network in the integrated energy system, x k is the volume flow of the kth heating pipeline in the heating network of the integrated energy system; The heat exergy loss of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Where, is the heat exergy flowing from the zth heating pipeline connected to the head end node of the kth heating pipeline in the heating network in the integrated energy system into the kth heating pipeline in the heating network, is the heat exergy flowing from the kth heating pipeline in the heating network of the integrated energy system to the bth heating pipeline connected to its terminal node, ρ is the density of the working fluid transmitted by the heating pipeline, c a is the specific heat capacity of the working medium transmitted by the heating pipeline, is the temperature of the head end node of the kth heating pipe in the heating network of the integrated energy system, X T0 is the static value of the heating network temperature in the integrated energy system, ε k is the composite heat transfer coefficient of the kth heating pipe in the heating network of the integrated energy system, z∈[1~S z ], S z is the number of heating pipes connected to the head end node of the kth heating pipe in the heating network of the integrated energy system, b∈[1~S b ], S b is the number of heating pipelines connected to the end node of the kth heating pipeline in the heating network of the integrated energy system.

3. The method according to claim 2, wherein The exergy coefficient c1 of the electric exergy is determined by the following formula: Where, φ1 is the energy level evaluation factor of electric energy, μ1 is the energy level factor of electric energy; The exergy coefficient c2 of thermal exergy is determined by the following formula: Where, φ2 is the energy level evaluation factor of thermal energy, and μ2 is the energy level factor of thermal energy; The resistance R of the i-th transmission line in the power supply network of the integrated energy system is determined by the following formula: i : Where, L i is the length of the i-th transmission line in the power supply network in the integrated energy system, A i is the transmission channel area of ​​the i-th transmission line in the power supply network of the integrated energy system, K i is the resistivity of the conductor material of the i-th transmission line of the power supply network in the integrated energy system; The fluid flow resistance R of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k : Where μ k is the fluid viscosity of the kth heating pipe in the heating network of the integrated energy system, L k is the length of the kth heating pipeline in the heating network of the integrated energy system, r k The inner radius of the kth heating pipe in the heating network of the integrated energy system; The volume flow rate x of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k : Where p k,s is the pressure at the head node of the kth heating pipeline in the heating network of the integrated energy system, p k,m is the pressure at the end node of the kth heating pipeline in the heating network of the integrated energy system; The composite heat transfer coefficient ε of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k : Where, d k(σ+1) is the radius of the σ+1th layer of material of the kth heating pipe in the heating network of the integrated energy system, d kσ is the radius of the σ-th layer of material of the k-th heating pipe in the heating network of the integrated energy system, γ kσ is the heat transfer coefficient of the σth layer material of the kth heating pipe in the heating network of the integrated energy system, σ∈(1~Ψ-1), Ψ is the number of layers of heat transfer material of the kth heating pipe in the heating network of the integrated energy system.

4. The method according to claim 1, wherein The constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraints, operating output inequality constraints of generator sets, gas units, and circulating water pumps in the power supply network, tolerance constraints, transmission capacity constraints, transformer tap position constraints, and pipeline temperature constraints.

5. The method according to claim 4, wherein The equality constraints are determined as follows: Where, is the input power of the qth power supply node in the power supply network of the integrated energy system, is the electricity demand of the fth load node in the power supply network in the integrated energy system, is the current of the i-th transmission line in the power supply network in the integrated energy system, is the voltage variation of the i-th transmission line in the power supply network of the integrated energy system, is the input thermal power of the dth heating node of the heating network in the integrated energy system, is the heat demand of the rth load node in the heating network in the integrated energy system, is the entropy of the kth heating pipe in the heating network in the integrated energy system, is the temperature variation of the kth heating pipe in the heating network in the integrated energy system, q∈[1~S q ], S q is the total number of power supply nodes in the power supply network of the integrated energy system, f∈[1~S f ], S f is the total number of load nodes in the power supply network of the integrated energy system, d∈[1~S d ], S d is the total number of heating nodes in the heating network of the integrated energy system, r∈[1~S r ] is the total number of load nodes in the heating network of the integrated energy system; The pipeline temperature constraint condition is determined as follows: Where C s is the first coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, C r is the first coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system, X Ts is the temperature matrix of the load heat flow inlet, X Tr is the temperature matrix of the load heat flow outlet, X T0 is the temperature static value matrix, b s is the second coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, b r It is the second coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system.

6. The method according to claim 5, wherein The first coefficient matrix C of the node temperature between the heating pipes of the heating network in the integrated energy system is determined by the following formula: s : Where, is the node υ in the heating network of the integrated energy system a Temperature and Node γ a The first coefficient between temperatures, υ a , γ a ∈(γ N ), γ N It is the set of nodes between heating pipes in the heating network of the integrated energy system; The second coefficient matrix b of the node temperature between the heating pipes of the heating network in the integrated energy system is determined by the following formula: s : Where, is the node υ in the heating network of the integrated energy system a The second coefficient of temperature, T is the transposed sign; Among them, when the heat medium is supplied from node γ a Through pipe M a Flow to node υ a hour: If the node υ a The node γ a The supplied heat medium and node γ a is a load node, then: If the node υ a The node γ a The supplied heat medium and node γ a is not a load node, then: If the node υ a The node γ a The supplied heat medium is In the above formula, The heating pipe M of the heating network in the integrated energy system a The volume flow rate, The heating pipe M of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe M of the heating network in the integrated energy system a length, is the node υ in the heating network of the integrated energy system a The total number of heating pipes supplying hot working medium, X T0 is the quiescent value of temperature, Node γ in the heating network of the integrated energy system a temperature, is the node υ in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature; The first coefficient matrix C of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r : Where, is the node η in the heating network of the integrated energy system a Temperature and node φ a The first coefficient between temperatures, φ a ,η a ∈(1~χ N ), χ N is the total number of nodes between return pipes in the heating network of the integrated energy system; The second coefficient matrix b of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r Where, is the node η in the heating network of the integrated energy system a The second coefficient of temperature; Among them, when the heat medium is supplied from node φ a Through pipe H a Flow to node η a hour: If the node η a Not only by node φ a When the hot working medium is supplied, then: If the node η a Only by node φ a When hot working medium is supplied, then: In the above formula, is the node φ in the heating network of the integrated energy system a The outlet temperature, is the node φ in the heating network of the integrated energy system a The volume flow rate, The heating pipe H of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe H of the heating network in the integrated energy system a length, is the node η in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature, is the node η in the heating network of the integrated energy system a The total number of return pipes supplying hot working medium.

7. The method according to claim 1, wherein The determining of the optimal transformer ratio of each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system includes: The optimal transformer ratio of the hth energy supply node in the power supply network of the integrated energy system is determined by the following formula: Where U eh,b is the voltage of the h-th energy supply node connected to the bus in the power supply network of the integrated energy system, is the optimal voltage of the hth energy supply node in the power supply network of the integrated energy system, h∈(1~S h ), S h is the total number of nodes in the power supply network of the integrated energy system.

8. An integrated energy system optimization and dispatching system taking into account exergy loss, characterized in that: The system comprises: A first determination module is used to determine the optimal voltage of each energy supply node of the power supply network in the integrated energy system according to the exergy loss of the integrated energy system; A second determination module is used to determine the optimal transformer ratio of each energy supply node in the power supply network of the integrated energy system according to the optimal voltage of each energy supply node in the power supply network of the integrated energy system; an adjustment module, configured to adjust the transformer ratio of each energy supply node in the power supply network of the integrated energy system to the optimal transformer ratio; The first determining module includes: A construction unit for establishing an objective function of an integrated energy system optimization scheduling model with the goal of minimizing the exergy loss of the integrated energy system; An acquisition unit, configured to solve the objective function of the integrated energy system optimization scheduling model based on the constraint conditions corresponding to the objective function of the integrated energy system optimization scheduling model, and obtain the optimal voltage of each energy supply node of the power supply network in the integrated energy system; The objective function of the integrated energy system optimization scheduling model is determined as follows: Where C loss is the exergy loss of the integrated energy system, c1 is the exergy mass coefficient of electrical exergy, c2 is the exergy mass coefficient of thermal exergy, w is the conversion coefficient between pressure exergy and electrical exergy, is the exergy loss of the i-th transmission line in the power supply network of the integrated energy system, is the pressure exergy loss of the kth heating pipeline in the heating network of the integrated energy system, is the heat exergy loss of the kth heating pipe in the heating network of the integrated energy system, k∈[1~S hl ], S hl is the number of heating pipes in the heating network of the integrated energy system, i∈[1~S el ], S el The number of transmission lines supplying the electricity network in an integrated energy system.

9. The system according to claim 8, wherein The exergy loss of the i-th transmission line in the power supply network of the integrated energy system is determined by the following formula: Where, is the exergy flowing from the xth transmission line connected to the head end node of the i-th transmission line in the power supply network of the integrated energy system into the i-th transmission line in the power supply network of the integrated energy system, is the exergy flowing from the i-th transmission line in the power supply network of the integrated energy system to the c-th transmission line connected to its terminal node, R i is the resistance of the i-th transmission line in the power supply network in the integrated energy system, P i,s is the electric power flowing into the head end node of the i-th transmission line in the power supply network of the integrated energy system, U i,s is the voltage of the node at the head end of the ith transmission line in the power supply network in the integrated energy system, x∈[1~S x ], S x is the sum of the transmission lines connected to the head end node of the i-th transmission line in the power supply network of the integrated energy system, c∈[1~S c ], S c is the sum of the transmission lines connected to the terminal node of the ith transmission line in the power supply network of the integrated energy system; The pressure exergy loss of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Where R k is the fluid flow resistance of the kth heating pipe in the heating network in the integrated energy system, x k is the volume flow of the kth heating pipeline in the heating network of the integrated energy system; The heat exergy loss of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: Where, is the heat exergy flowing from the zth heating pipeline connected to the head end node of the kth heating pipeline in the heating network in the integrated energy system into the kth heating pipeline in the heating network, is the heat exergy flowing from the kth heating pipeline in the heating network of the integrated energy system to the bth heating pipeline connected to its terminal node, ρ is the density of the working medium transmitted by the heating pipeline, c a is the specific heat capacity of the working medium transmitted by the heating pipeline, is the temperature of the head end node of the kth heating pipe in the heating network of the integrated energy system, X T0 is the static value of the heating network temperature in the integrated energy system, ε k is the composite heat transfer coefficient of the kth heating pipe in the heating network of the integrated energy system, z∈[1~S z ], S z is the number of heating pipes connected to the head end node of the kth heating pipe in the heating network of the integrated energy system, b∈[1~S b ], S b is the number of heating pipelines connected to the end node of the kth heating pipeline in the heating network of the integrated energy system.

10. The system according to claim 9, wherein: The exergy coefficient c1 of the electric exergy is determined by the following formula: Where, φ1 is the energy level evaluation factor of electric energy, μ1 is the energy level factor of electric energy; The exergy coefficient c2 of thermal exergy is determined by the following formula: Where, φ2 is the energy level evaluation factor of thermal energy, and μ2 is the energy level factor of thermal energy; The resistance R of the i-th transmission line in the power supply network of the integrated energy system is determined by the following formula: i : Where, L i is the length of the i-th transmission line in the power supply network in the integrated energy system, A i is the transmission channel area of ​​the i-th transmission line in the power supply network of the integrated energy system, K i is the resistivity of the conductor material of the i-th transmission line of the power supply network in the integrated energy system; The fluid flow resistance R of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k : Where μ k is the fluid viscosity of the kth heating pipe in the heating network of the integrated energy system, L k is the length of the kth heating pipeline in the heating network of the integrated energy system, r k The inner radius of the kth heating pipe in the heating network of the integrated energy system; The volume flow rate x of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k : Where p k,s is the pressure at the head node of the kth heating pipeline in the heating network of the integrated energy system, p k,m is the pressure at the end node of the kth heating pipeline in the heating network of the integrated energy system; The composite heat transfer coefficient ε of the kth heating pipe in the heating network of the integrated energy system is determined by the following formula: k : Where, d k(σ+1) is the radius of the σ+1th layer of material of the kth heating pipe in the heating network of the integrated energy system, d kσ is the radius of the σ-th layer of material of the k-th heating pipe in the heating network of the integrated energy system, γ kσ is the heat transfer coefficient of the σth layer material of the kth heating pipe in the heating network of the integrated energy system, σ∈(1~Ψ-1), Ψ is the number of layers of heat transfer material of the kth heating pipe in the heating network of the integrated energy system.

11. The system according to claim 8, wherein The constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraints, operating output inequality constraints of generator sets, gas units, and circulating water pumps in the power supply network, tolerance constraints, transmission capacity constraints, transformer tap position constraints, and pipeline temperature constraints.

12. The system according to claim 11, wherein The equality constraints are determined as follows: Where, is the input power of the qth power supply node in the power supply network of the integrated energy system, is the electricity demand of the fth load node in the power supply network in the integrated energy system, is the current of the i-th transmission line in the power supply network in the integrated energy system, is the voltage variation of the i-th transmission line in the power supply network of the integrated energy system, is the input thermal power of the dth heating node of the heating network in the integrated energy system, is the heat demand of the rth load node in the heating network in the integrated energy system, is the entropy of the kth heating pipe in the heating network in the integrated energy system, is the temperature variation of the kth heating pipe in the heating network in the integrated energy system, q∈[1~S q ], S q is the total number of power supply nodes in the power supply network of the integrated energy system, f∈[1~S f ], S f is the total number of load nodes in the power supply network of the integrated energy system, d∈[1~S d ], S d is the total number of heating nodes in the heating network of the integrated energy system, r∈[1~S r ] is the total number of load nodes in the heating network of the integrated energy system; The pipeline temperature constraint condition is determined as follows: Where C s is the first coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, C r is the first coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system, X Ts is the temperature matrix of the load heat flow inlet, X Tr is the temperature matrix of the load heat flow outlet, X T0 is the temperature static value matrix, b s is the second coefficient matrix of the node temperature between the heating pipes of the heating network in the integrated energy system, b r It is the second coefficient matrix of the node temperature between the return pipes of the heating network in the integrated energy system.

13. The system according to claim 12, wherein: The first coefficient matrix C of the node temperature between the heating pipes of the heating network in the integrated energy system is determined by the following formula: s : Where, is the node υ in the heating network of the integrated energy system a Temperature and Node γ a The first coefficient between temperatures, υ a , γ a ∈(γ N ), γ N is the set of nodes between heating pipes in the heating network of the integrated energy system, C s Initially, all elements in the matrix are 0; The second coefficient matrix b of the node temperature between the heating pipes of the heating network in the integrated energy system is determined by the following formula: s : Where, is the node υ in the heating network of the integrated energy system a The second coefficient of temperature, T is the transposed sign, b s Initially, all elements in the matrix are 0; Among them, when the heat medium is supplied from node γ a Through pipe M a Flow to node υ a hour: If the node υ a The node γ a The supplied heat medium and node γ a is a load node, then: If the node υ a The node γ a The supplied heat medium and node γ a is not a load node, then: If the node υ a The node γ a The supplied heat medium is In the above formula, The heating pipe M of the heating network in the integrated energy system a The volume flow rate, The heating pipe M of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe M of the heating network in the integrated energy system a length, is the node υ in the heating network of the integrated energy system a The total number of heating pipes supplying hot working medium, X T0 is the quiescent value of temperature, Node γ in the heating network of the integrated energy system a temperature, is the node υ in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature; The first coefficient matrix C of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r : Where, is the node η in the heating network of the integrated energy system a Temperature and node φ a The first coefficient between temperatures, φ a ,η a ∈(1~χ N ), χ N is the total number of nodes between return pipes in the heating network of the integrated energy system, C r Initially, all elements in the matrix are 0; The second coefficient matrix b of the node temperature between the return pipes of the heating network in the integrated energy system is determined by the following formula: r Where, is the node η in the heating network of the integrated energy system a The second coefficient of temperature, b r Initially, all elements in the matrix are 0; Among them, when the heat medium is supplied from node φ a Through pipe H a Flow to node η a hour: If the node η a Not only by node φ a When the hot working medium is supplied, then: If the node η a Only by node φ a When hot working medium is supplied, then: In the above formula, is the node φ in the heating network of the integrated energy system a The outlet temperature, is the node φ in the heating network of the integrated energy system a The volume flow rate, The heating pipe H of the heating network in the integrated energy system a The composite heat transfer coefficient, The heating pipe H of the heating network in the integrated energy system a length, is the node η in the heating network of the integrated energy system a The first coefficient between the temperature and its own node temperature, is the node η in the heating network of the integrated energy system a The total number of return pipes supplying hot working medium.

14. The system according to claim 8, wherein The second determining module is configured to: The optimal transformer ratio of the hth energy supply node in the power supply network of the integrated energy system is determined by the following formula: Where U eh,b is the voltage of the h-th energy supply node connected to the bus in the power supply network of the integrated energy system, is the optimal voltage of the hth energy supply node in the power supply network of the integrated energy system, h∈(1~S h ), S h is the total number of nodes in the power supply network of the integrated energy system.

Images