CN111832138B - Regional pipe network topology optimization method - Google Patents

Regional pipe network topology optimization method Download PDF

Info

Publication number
CN111832138B
CN111832138B CN202010425292.9A CN202010425292A CN111832138B CN 111832138 B CN111832138 B CN 111832138B CN 202010425292 A CN202010425292 A CN 202010425292A CN 111832138 B CN111832138 B CN 111832138B
Authority
CN
China
Prior art keywords
pipe network
steam
unit
cost
yuan
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010425292.9A
Other languages
Chinese (zh)
Other versions
CN111832138A (en
Inventor
高为
王世朋
张海珍
阮慧锋
周宇昊
王彧斐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huadian Electric Power Research Institute Co Ltd
Original Assignee
Huadian Electric Power Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huadian Electric Power Research Institute Co Ltd filed Critical Huadian Electric Power Research Institute Co Ltd
Priority to CN202010425292.9A priority Critical patent/CN111832138B/en
Publication of CN111832138A publication Critical patent/CN111832138A/en
Application granted granted Critical
Publication of CN111832138B publication Critical patent/CN111832138B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes

Abstract

The invention discloses a regional pipe network topology optimization method, which is used for obtaining a better pipe network topology type and providing theoretical guidance for regional pipe network design optimization. The invention takes the total annual cost as an objective function, provides a new pipe network connection type (Steiner tree type), increases energy loss, namely operation cost, adopts different graph theory algorithms to respectively solve and compare and analyze different pipe network type algorithms so as to obtain a better pipe network topology type, and provides theoretical guidance for the pipe network design optimization of the regional system. Therefore, the method not only can minimize the annual total operation cost of pipe network projects, but also can achieve the purpose of energy conservation. The invention provides an effective method proposal for planning and optimizing the regional pipe network.

Description

Regional pipe network topology optimization method
Technical Field
The invention relates to a regional pipe network topology optimization method, which is used for obtaining a better pipe network topology type and providing theoretical guidance for regional pipe network design optimization.
Background
The renewable clean energy sources such as natural gas, biomass energy, solar energy, hydrogen energy and the like are taken as primary energy sources, and a small or miniature independent electric and thermal (cold) energy output system at a user side is used for realizing energy cascade utilization, and a support and supplement comprehensive energy production system is provided through a central energy supply system, wherein a transmission pipe network is an important component part in a distributed energy system, and the arrangement of pipes, pipe diameter optimization, selection of thickness of a heat preservation layer, determination of a transmission distance, hydraulic balance calculation and the like of the system have great influence on efficiency and economy of the system.
Aiming at pipe diameter optimization, heat preservation layer thickness selection, conveying distance determination and the like, students consider the changes of cold and heat supply loads and the operation adjustment of different time periods throughout the year, and adopt a genetic algorithm to solve the optimal combination of discrete diameters. Also, the learner can determine the pipe diameter by using integer coding genetic algorithm through improving genetic operators such as recombination and the like. Xia Bo an optimization model of the regional cooling area and the furthest conveying distance is established based on a life cycle cost analysis method, and an optimization model limiting condition and a solving method are provided. And the scholars also research the influence relationship of the base station action radius, the secondary pump position and the like on the system, construct the energy loss and investment function of the pipe network, and solve the problems by adopting an interior point algorithm. The method also has the advantages that a learner establishes a pipe network optimization model of the regional cooling system, and calculates and analyzes design parameters of the pipe network by applying an orthogonal optimization method.
And for the optimization of pipe network arrangement, a learner combines the load characteristics of the urban comprehensive cold and heat source system and the influence factors of the characteristics thereof to build a model for optimizing the pipeline of the urban comprehensive cold and heat source system, and adopts a genetic algorithm to solve the dendritic pipe network and the annular pipe network. Chen Juan and other scholars provide a decision-making and system-based Regional Distributed Energy System (RDES) double-layer optimization structure for the distributed energy system in the planning stage, and perform the pareto optimal solution of the multi-objective problem by adopting an NSGA-II multi-objective genetic algorithm for system optimization; for the layout and optimization of the station network, an energy distance concept of the annual value of the pipe network construction cost between supply and demand energy nodes is provided, a layout model for solving the bits in the RDES station network diagram P is constructed based on the energy distance concept, an RDES addressing network distribution algorithm based on the energy distance concept is adopted, and a determined optimal solution is finally obtained on the basis of traversing all feasible solutions in a solution set space. Feng Xiao the scholars of equal degree establish the overall optimization design model of the regional cooling system, select genetic algorithm to solve, put forward integer coding to genetic algorithm design and operation according to the research problem, adopt the optimal individual preservation strategy, the improved strategies such as the cross rate and the variation rate self-adaptation, use the heuristic back-pass method or the double-pass method to select loops at the same time so as to improve the pipe network resolving speed, apply the improved single-parent genetic algorithm to carry out the optimization arrangement of the dendritic pipe network, and design the corresponding fitness function, single-parent transposition operator and inversion operator. The Wang Chinese scholars combine the pipe network layout genetic optimization algorithm and the pipe diameter successive approximation optimization algorithm to construct a new algorithm to solve the problem of optimizing a huge amount of dimensional solution space in the overall optimization design of the complex dendritic pipe network, so that the pipe network optimization design efficiency is improved.
For optimization of a distributed energy system pipe network, most of students aim at lowest cost, constraint conditions such as user demands, flow rates and pipe diameters are adopted, a layout optimization problem is generally converted into a minimum tree (MST) problem in graph theory, the pipe network is regarded as a dendritic pipe network by default, and a genetic algorithm or an improved genetic algorithm is used for solving. The invention provides a new pipe network connection type (Steiner tree type) based on the method, takes annual total cost as an objective function, increases and considers energy loss, namely operation cost, and adopts different graph theory algorithms to respectively solve and compare and analyze different pipe network type algorithms so as to obtain a better pipe network topology type, thereby providing theoretical guidance for pipe network design optimization of a regional pipe network system.
The Kruskal algorithm is a classical algorithm for solving a minimum spanning tree of a weighted connected graph, based on greedy thought, the calculation process of the Kruskal algorithm takes edges as a main factor, and the Prim algorithm is also corresponding to the Kruskal algorithm, and the calculation process of the Kruskal algorithm takes vertexes as a main factor. Since the former can avoid multiple ordering and is more efficient than the latter, the Kruskal algorithm is chosen for calculation here, and the flow is as follows:
(1) Knowing the weighted connectivity graph g= (V, E, W), where V, E and W are the set of vertices, edges and weights of edges in the connectivity graph, respectively, the edges in the initial minimum spanning tree are empty sets;
(2) Sequencing each edge in the connected graph according to the weight from small to large;
(3) Sequentially adding edges into the minimum spanning tree according to the sequence and judging whether the minimum spanning tree is looped, if yes, discarding the added edges and entering the step (4), otherwise, directly entering the step (4);
(4) Judging whether the number of edges in the minimum spanning tree is Np-1 (Np is the number of vertices of the connected graph), if so, outputting the result of the minimum spanning tree after the calculation is finished, otherwise, turning to the step (3) until the calculation is finished.
The geoSteiner algorithm is a relatively efficient method for accurately solving the Steiner tree at present. The generation of ESMT and RSMT involves two very important steps, namely the generation and concatenation of the Man Sitan nanotree. A preprocessing stage is performed before Man Sitan nanotree generation to reduce the complexity of the computation process. Then, the Mantanin tree is gradually generated by a recursion algorithm, and only the results after the distance test, the long-side test, the branch tree test and the Mantanin tree type test are passed are reserved in the final connection stage. The Man Sitan nanotree connection can be regarded as a problem of finding the minimum spanning tree in the hypergraph, the process can be solved through integer programming modeling, and also can be solved through dynamic programming.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides a regional pipe network topology optimization method which comprehensively considers all limiting factors by using the pipe network connection type of the Steiner tree.
The invention solves the problems by adopting the following technical scheme: the regional pipe network topology optimization method is characterized by comprising the following steps of:
the objective function is shown in formula (1):
min C=C pipe +C pressureloss +C heatloss (27)
Figure BDA0002498485830000031
Figure BDA0002498485830000032
Figure BDA0002498485830000033
wherein C (yuan-year) -1 ) Is the total annual cost of the device,
Figure BDA0002498485830000034
and->
Figure BDA0002498485830000035
The investment cost, the pressure loss cost and the heat loss cost of the annual steam pipe network are respectively, and n is the number of branch pipe sections in the pipe network system.
Investment cost of steam pipe network
Figure BDA0002498485830000036
Obtained by the following formula: />
Figure BDA0002498485830000037
Wherein t is year (year) is the period of use of the pipe network system, I is the annual interest rate, a k (yuan·m -1 ) Is the unit price of a steam pipeline, L k (m) is the length of each section of branch line;
a k (yuan·m -1 ) The inner diameter of the pipeline required in the calculation process of (1) is obtained by the formula (6):
Figure BDA0002498485830000038
in the method, in the process of the invention,
Figure BDA0002498485830000039
is the inner diameter of the steam pipeline, W k (kg·s -1 ) Is the mass flow rate ρ (kg.m) -3 ) For the steam density, u (m.s) -1 ) Is the flow rate; wherein each section of branch line W k Is different in the whole pipe network system; therefore, W is calculated by solving the following linear programming model k Is a value of (2);
Figure BDA00024984858300000310
constraint conditions:
Figure BDA00024984858300000311
Figure BDA00024984858300000312
Figure BDA0002498485830000041
w α,β,γ ≥0 (37)
w β,α,γ ≥0 (38)
wherein w is α,β,γ (kg·s -1 ) Is the mass flow rate in the branch line connecting the vertices α and β, b binary Is a binary variable, which represents whether two vertexes are connected by a pipeline or not, and is obtained according to a pipe network topology result by a formula (9); w (W) γ Is the mass flow rate of the apex γ; if the user produces steam, W γ Is a negative value, otherwise, it is a positive value; the steam flow in each branch line has two directions from alpha to beta and from beta to alpha (-1) in equation (8) r For controlling the flow direction of the steam;
for medium-low pressure steam, the steam is heated,
Figure BDA0002498485830000042
and Wt k Calculated from the following formula:
Figure BDA0002498485830000044
Figure BDA0002498485830000045
in the case of high-pressure steam,
Figure BDA0002498485830000046
and Wt k Calculated from the following formula:
Figure BDA0002498485830000047
Figure BDA0002498485830000048
the system steam belongs to medium-low pressure steam and is calculated according to formulas (13) - (14);
a k obtained from equation (17):
Figure BDA0002498485830000049
in the method, in the process of the invention,
Figure BDA00024984858300000410
is the outer diameter of the branch line, wt k (kg·m -1 ) Is the weight of the branch line per unit length, A 1 (yuan·kg -1 ) Representing the cost per unit weight of the pipeline, A 2 (yuan·m -0.48 ) Representing the installation cost, A 3 (yuan·m -1 ) Representing the road use cost A 4 (yuan·m -1 ) Representing the cost of maintaining the temperature.
The cost of the pressure loss is as follows:
the steam has pressure, so that the auxiliary conveying of a compressor is not needed, but certain pressure loss is caused in the conveying process, the steam enters the steam production device after being condensed and pressurized by a water pump, and therefore, the pressure loss cost of a pipe network is low
Figure BDA00024984858300000413
Calculating according to the power consumption of the water pump;
Figure BDA00024984858300000411
the method is obtained by the following formula:
Figure BDA00024984858300000412
wherein a is E (yuan·kW -1 ·h -1 ) Is the electricity charge, t time (h) Is the number of years of operation of the device, N k (W) is the shaft power against resistance at the time of steam delivery, and is obtained by the following formula:
Figure BDA0002498485830000051
where η is the efficiency of the delivery device, ne k And (W) is the effective power, calculated by the following formula:
Figure BDA00024984858300000513
Figure BDA00024984858300000514
representing the head loss due to the on-way resistance, calculated according to equation (21):
Figure BDA0002498485830000052
Figure BDA0002498485830000053
Figure BDA0002498485830000054
Figure BDA0002498485830000055
in zeta k Zeta is the local resistance coefficient when the section suddenly increases or decreases E Is the local resistance coefficient at the bend,
Figure BDA0002498485830000056
the sectional area of the branch line, σ is the coefficient of friction of the pipeline.
The heat loss costs are as follows:
Figure BDA0002498485830000057
calculated according to equation (25):
Figure BDA0002498485830000058
wherein a (yuan. Kg) -1 ) And q (kJ.kg) -1 ) Steam unit price and latent heat of steam, Q k (kJ·m -1 ·s -1 ) Is heat loss;
in the presence of insulation, the heat loss of the pipe network is calculated by the following formula:
Figure BDA0002498485830000059
/>
wherein T (DEG C) is the outer surface temperature of the pipe, T a The temperature (DEG C) is the ambient temperature,
Figure BDA00024984858300000510
and->
Figure BDA00024984858300000511
The outer diameter and the inner diameter of the heat insulation layer are respectively;
λ(W·m -1 ·k -1 ) Is the heat conductivity coefficient of the thermal insulation material product at the average temperature, epsilon (W.m) -2 ·k -1 ) Is the heat exchange coefficient between the outer surface of the heat-insulating layer and the atmosphere, and the inner diameter of the heat-insulating layer
Figure BDA00024984858300000512
Taking the steam pipe outside diameter, the insulation layer thickness is related to the pipe diameter and the temperature of the fluid to be delivered, the invention is selected according to table 1.
Table 1 table for selecting thickness of insulation layer
Figure BDA0002498485830000061
Compared with the prior art, the invention has the following advantages and effects: the annual total cost is used as an objective function, a new pipe network connection type (Steiner tree type) is provided, energy loss, namely operation cost, is increased, different pipe network type algorithms are respectively solved and compared and analyzed by adopting different graph theory algorithms, so that a better pipe network topology type is obtained, and theoretical guidance is provided for regional system pipe network design optimization. The invention can save energy and simultaneously minimize the annual total operation cost of pipe network projects.
Drawings
FIG. 1 is a diagram of Star-type pipe network in an embodiment of the present invention.
FIG. 2 is a diagram of an MST-type pipe network in accordance with an embodiment of the present invention.
FIG. 3 is a diagram of an RMST-type pipe network in an embodiment of the invention.
FIG. 4 is a diagram of an EMST-type pipe network according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail by way of examples with reference to the accompanying drawings, which are illustrative of the present invention and not limited to the following examples.
Examples
In this embodiment, a topology optimization method for a regional pipe network specifically includes:
the objective function is shown in formula (1):
min C=C pipe +C pressureloss +C heatloss (53)
Figure BDA0002498485830000071
Figure BDA0002498485830000072
Figure BDA0002498485830000073
wherein C (yuan-year) -1 ) Is the total annual cost of the device,
Figure BDA0002498485830000074
and->
Figure BDA0002498485830000075
The investment cost, the pressure loss cost and the heat loss cost of the annual steam pipe network are respectively, and n is the number of branch pipe sections in the pipe network system.
Investment cost of steam pipe network
Figure BDA0002498485830000076
Obtained by the following formula:
Figure BDA0002498485830000077
wherein t is year (year) is the period of use of the pipe network system, I is the annual interest rate, a k (yuan·m -1 ) Is the unit price of a steam pipeline, L k (m) is the length of each section of branch line;
a k (yuan·m -1 ) The inner diameter of the pipeline required in the calculation process of (1) is obtained by the formula (6):
Figure BDA0002498485830000078
in the method, in the process of the invention,
Figure BDA0002498485830000079
is the inner diameter of the steam pipeline, W k (kg·s -1 ) Is the mass flow rate ρ (kg.m) -3 ) For the steam density, u (m.s) -1 ) Is the flow rate; wherein each section of branch line W k Is different in the whole pipe network system; therefore, W is calculated by solving the following linear programming model k Is a value of (2);
Figure BDA00024984858300000710
constraint conditions:
Figure BDA00024984858300000711
Figure BDA00024984858300000712
Figure BDA00024984858300000713
w α,β,γ ≥0 (63)
w β,α,γ ≥0 (64)
wherein w is α,β,γ (kg·s -1 ) Is the mass flow rate in the branch line connecting the vertices α and β, b binary Is a binary variable representing two verticesWhether pipelines are connected or not is obtained according to a pipe network topology result by a formula (9); w (W) γ Is the mass flow rate of the apex γ; if the user produces steam, W γ Is a negative value, otherwise, it is a positive value; the steam flow in each branch line has two directions from alpha to beta and from beta to alpha (-1) in equation (8) r For controlling the flow direction of the steam;
for medium-low pressure steam, the steam is heated,
Figure BDA0002498485830000081
and Wt k Calculated from the following formula:
Figure BDA0002498485830000082
Figure BDA0002498485830000083
in the case of high-pressure steam,
Figure BDA0002498485830000084
and Wt k Calculated from the following formula:
Figure BDA0002498485830000085
Figure BDA0002498485830000086
the system steam belongs to medium-low pressure steam and is calculated according to formulas (13) - (14);
a k obtained from equation (17):
Figure BDA0002498485830000087
in the method, in the process of the invention,
Figure BDA0002498485830000088
is the outer diameter of the branch line, wt k (kg·m -1 ) Is the weight of the branch line per unit length, A 1 (yuan·kg -1 ) Representing the cost per unit weight of the pipeline, A 2 (yuan·m -0.48 ) Representing the installation cost, A 3 (yuan·m -1 ) Representing the road use cost A 4 (yuan·m -1 ) Representing the cost of maintaining the temperature.
The cost of the pressure loss is as follows:
the steam has pressure, so that the auxiliary conveying of a compressor is not needed, but certain pressure loss is caused in the conveying process, the steam enters the steam production device after being condensed and pressurized by a water pump, and therefore, the pressure loss cost of a pipe network is low
Figure BDA00024984858300000811
Calculating according to the power consumption of the water pump;
Figure BDA0002498485830000089
the method is obtained by the following formula:
Figure BDA00024984858300000810
wherein a is E (yuan·kW -1 ·h -1 ) Is the electricity charge, t time (h) Is the number of years of operation of the device, N k (W) is the shaft power against resistance at the time of steam delivery, and is obtained by the following formula:
Figure BDA0002498485830000091
where η is the efficiency of the delivery device, ne k And (W) is the effective power, calculated by the following formula:
Figure BDA0002498485830000092
Figure BDA00024984858300000914
representing the head loss due to the on-way resistance, calculated according to equation (21):
Figure BDA0002498485830000093
Figure BDA0002498485830000094
Figure BDA0002498485830000095
Figure BDA0002498485830000096
/>
in zeta k Zeta is the local resistance coefficient when the section suddenly increases or decreases E Is the local resistance coefficient at the bend,
Figure BDA0002498485830000097
the sectional area of the branch line, σ is the coefficient of friction of the pipeline.
The heat loss costs are as follows:
Figure BDA0002498485830000098
calculated according to equation (25):
Figure BDA0002498485830000099
wherein a (yuan. Kg) -1 ) And q (kJ.kg) -1 ) Steam unit price and latent heat of steam, Q k (kJ·m -1 ·s -1 ) Is heat loss;
in the presence of insulation, the heat loss of the pipe network is calculated by the following formula:
Figure BDA00024984858300000910
wherein T (DEG C) is the outer surface temperature of the pipe, T a The temperature (DEG C) is the ambient temperature,
Figure BDA00024984858300000911
and->
Figure BDA00024984858300000912
The outer diameter and the inner diameter of the heat insulation layer are respectively;
λ(W·m -1 ·k -1 ) Is the heat conductivity coefficient of the thermal insulation material product at the average temperature, epsilon (W.m) -2 ·k -1 ) Is the heat exchange coefficient between the outer surface of the heat-insulating layer and the atmosphere, and the inner diameter of the heat-insulating layer
Figure BDA00024984858300000913
Taking the steam pipe outside diameter, the insulation layer thickness is related to the pipe diameter and the temperature of the fluid to be delivered, the invention is selected according to table 1.
Table 1 table for selecting thickness of insulation layer
Figure BDA0002498485830000101
In this embodiment, there are four types of pipe network connections between the heat source and each user, namely Star, MST, RSMT and ESMT. Wherein MST is solved by Kruskal algorithm, RSMT and ESMT are solved by GeoSteiner algorithm, and the whole optimization process is realized by C++ programming.
The case data is from a heat supply pipe network project, and the geodetic coordinates (longitude and latitude are converted into geodetic coordinates if longitude and latitude are adopted) and the heat load of the heat source and each user are extracted. The heat load of the user is specified to be positive and the heat load of the heat source point is specified to be negative. The data are summarized in table 2.
Table 2 heat source raw data
Figure BDA0002498485830000102
Figure BDA0002498485830000111
The required parameters in the model are shown in table 3 below.
TABLE 3 parameter Table
Figure BDA0002498485830000112
The case calculations are summarized below. FIGS. 1-4 are respectively topology schemes of four pipe network types, in which the diameter of the branch line varies with the thickness of the line segment.
Table 4 calculation results
Figure BDA0002498485830000113
As shown in Table 4, in this example, the EMST-type pipe network has the lowest total cost, while the RMST-type pipe network has the highest annual total cost. Star and MST type networks are more similar in cost, but still 13.6% and 10.1% higher than EMST type networks, respectively. The star-shaped pipe network connection mode is the simplest as can be seen from the pipe network topology scheme, but the layout is the least reasonable in engineering practice, and especially when a heat source is out of a user area, the pipeline length can be greatly increased, so that the pipe network conveying efficiency is reduced. In contrast, the schemes of MST, RSMT and EMST are reasonable. The MST pipe network is relatively simple in connection mode, but the length of the pipeline can be increased to a certain extent, so that energy conservation and consumption reduction are not facilitated. The RSMT pipeline has regular path, is beneficial to the maintenance and management of a pipe network to a certain extent, but has the highest total cost. In the comprehensive view, the EMST type pipe network has the shortest length and the lowest cost, and is the optimal pipe network topology type.
What is not described in detail in this specification is all that is known to those skilled in the art.
Although the present invention has been described with reference to the above embodiments, it should be understood that the invention is not limited to the embodiments described above, but is capable of modification and variation without departing from the spirit and scope of the present invention.

Claims (1)

1. The regional pipe network topology optimization method is characterized by comprising the following steps of:
the objective function is shown in formula (1):
min C=C pipe +C pressureloss +C heatloss (1)
Figure FDA0004144486040000011
Figure FDA0004144486040000012
Figure FDA0004144486040000013
wherein C is the annual total cost in units of yuan-year -1
Figure FDA0004144486040000014
And->
Figure FDA0004144486040000015
The investment cost, the pressure loss cost and the heat loss cost of the annual steam pipe network are respectively expressed in the units of yuan-year -1 The method comprises the steps of carrying out a first treatment on the surface of the n is the number of branch pipe sections in the pipe network system;
investment cost of steam pipe network
Figure FDA0004144486040000016
Obtained by the following formula:
Figure FDA0004144486040000017
wherein t is year The service cycle of the pipe network system is given by year; i is annual interest rate; a, a k Is the unit price of the steam pipeline, and the unit is yuan.m -1 ;L k Is the length of each section of branch line, and the unit is m;
a k the inner diameter of the pipeline is required in the calculation process of (a) and is expressed in terms of yuan.m -1 Obtained by the formula (6):
Figure FDA0004144486040000018
in the method, in the process of the invention,
Figure FDA0004144486040000019
is the inner diameter of the steam pipeline, and the unit is m; w (W) k The mass flow rate is in kg.s -1 The method comprises the steps of carrying out a first treatment on the surface of the ρ is the steam density in kg.m -3 The method comprises the steps of carrying out a first treatment on the surface of the u is the flow rate in m.s -1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein each section of branch line W k Is different in the whole pipe network system; therefore, W is calculated by solving the following linear programming model k Is a value of (2);
Figure FDA00041444860400000110
constraint conditions:
Figure FDA00041444860400000111
Figure FDA00041444860400000112
Figure FDA0004144486040000021
w α,β,γ ≥0 (11)
w β,α,γ ≥0 (12)
wherein w is α,β,γ Is the mass flow rate in kg.s in the branch line connecting the vertices α and β -1 ;b binary Is a binary variable, which represents whether two vertexes are connected by a pipeline or not, and is obtained according to a pipe network topology result by a formula (9); w (W) γ Is the mass flow rate of the apex γ; if the user produces steam, W γ Is a negative value, otherwise, it is a positive value; the steam flow in each branch line has two directions from alpha to beta and from beta to alpha (-1) in equation (8) r For controlling the flow direction of the steam;
for medium-low pressure steam, the steam is heated,
Figure FDA0004144486040000022
and Wt k Calculated from the following formula:
Figure FDA0004144486040000023
Figure FDA0004144486040000024
in the case of high-pressure steam,
Figure FDA0004144486040000025
and Wt k Calculated from the following formula:
Figure FDA0004144486040000026
Figure FDA0004144486040000027
the system steam belongs to medium-low pressure steam and is calculated according to formulas (13) - (14);
a k obtained from equation (17):
Figure FDA0004144486040000028
in the method, in the process of the invention,
Figure FDA0004144486040000029
is the outer diameter of the branch line in m; wt (Wt) k Is the weight of the branch line per unit length, and the unit is kg.m -1 :A 1 Represents the cost per unit weight of pipeline in terms of yuan kg -1 ;A 2 Representing the installation cost in units of yuan.m -0.48 ;A 3 Represents road use cost, and the unit is yuan.m -1 ;A 4 Represents the heat-preserving cost, and the unit is yuan.m -1
The cost of the pressure loss is as follows:
the steam has pressure, so that the auxiliary conveying of a compressor is not needed, but certain pressure loss is caused in the conveying process, the steam enters the steam production device after being condensed and pressurized by a water pump, and therefore, the pressure loss cost of a pipe network is low
Figure FDA00041444860400000210
Calculating according to the power consumption of the water pump;
Figure FDA00041444860400000211
the method is obtained by the following formula:
Figure FDA0004144486040000031
wherein a is E The unit is yuan kW for electricity charge -1 ·h -1 ;t time The unit is h, which is the number of annual operation hours of the device; n (N) k Is the shaft power which overcomes the resistance in the steam delivery, and the unit is W; the method is obtained by the following formula:
Figure FDA0004144486040000032
where η is the efficiency of the delivery device, ne k Is effective power, the unit is W, calculated by the following formula:
Figure FDA0004144486040000033
Figure FDA0004144486040000034
representing the head loss caused by the on-way resistance, the unit is m, calculated according to formula (21):
Figure FDA0004144486040000035
Figure FDA0004144486040000036
Figure FDA0004144486040000037
Figure FDA0004144486040000038
/>
in zeta k A local resistance coefficient when the section suddenly increases or decreases; zeta type E Is the local resistance coefficient at the elbow;
Figure FDA0004144486040000039
the unit of the cross-sectional area of the branch line is m 2 The method comprises the steps of carrying out a first treatment on the surface of the Sigma is the friction coefficient of the pipeline;
the heat loss costs are as follows:
Figure FDA00041444860400000310
calculated according to equation (25):
Figure FDA00041444860400000311
wherein a is a steam unit price and the unit is yuan.kg -1 The method comprises the steps of carrying out a first treatment on the surface of the q is the latent heat of steam in kJ.kg -1 ;Q k Is the heat loss in kJ.m -1 ·s -1
In the presence of insulation, the heat loss of the pipe network is calculated by the following formula:
Figure FDA00041444860400000312
wherein T is the temperature of the outer surface of the pipeline and the unit is the temperature; t (T) a Is ambient temperature in degrees celsius;
Figure FDA00041444860400000313
and->
Figure FDA00041444860400000314
The outside diameter and the inside diameter of the heat insulation layer are respectively m; lambda is the coefficient of thermal conductivity of the thermal insulation material product at average temperature, and the unit is W.m -1 ·k -1 The method comprises the steps of carrying out a first treatment on the surface of the Epsilon is the heat exchange coefficient between the outer surface of the heat-insulating layer and the atmosphere, and the unit is W.m -2 ·k -1 The method comprises the steps of carrying out a first treatment on the surface of the Inner diameter of heat insulation layer>
Figure FDA0004144486040000041
Taking the outer diameter of the steam pipeline, the thickness of the heat preservation layer is related to the diameter of the pipeline and the temperature of the conveyed fluid. />
CN202010425292.9A 2020-05-19 2020-05-19 Regional pipe network topology optimization method Active CN111832138B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010425292.9A CN111832138B (en) 2020-05-19 2020-05-19 Regional pipe network topology optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010425292.9A CN111832138B (en) 2020-05-19 2020-05-19 Regional pipe network topology optimization method

Publications (2)

Publication Number Publication Date
CN111832138A CN111832138A (en) 2020-10-27
CN111832138B true CN111832138B (en) 2023-05-09

Family

ID=72913713

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010425292.9A Active CN111832138B (en) 2020-05-19 2020-05-19 Regional pipe network topology optimization method

Country Status (1)

Country Link
CN (1) CN111832138B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112395718B (en) * 2020-11-16 2023-12-19 新奥新智科技有限公司 Pipe network layout generation method and device of energy system
CN112565272B (en) * 2020-12-09 2022-05-17 中国人民解放军国防科技大学 Method and device for blocking minimum Steiner tree of double-layer network and computer equipment
FR3133648A1 (en) * 2022-03-16 2023-09-22 Renault Sas Method for generating functional cooling circuit architectures
CN116910951B (en) * 2023-07-04 2024-01-23 天津大学 Regional heating pipe network optimization design method and device

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7593839B1 (en) * 2002-03-07 2009-09-22 Bentley Systems, Incorporated Method for optimizing design and rehabilitation of water distribution systems
CN110717226A (en) * 2019-09-30 2020-01-21 国网浙江杭州市萧山区供电有限公司 Regional comprehensive energy system layout planning method considering energy network topological characteristics

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6826737B2 (en) * 2000-12-06 2004-11-30 Cadence Design Systems, Inc. Recursive partitioning placement method and apparatus

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7593839B1 (en) * 2002-03-07 2009-09-22 Bentley Systems, Incorporated Method for optimizing design and rehabilitation of water distribution systems
CN110717226A (en) * 2019-09-30 2020-01-21 国网浙江杭州市萧山区供电有限公司 Regional comprehensive energy system layout planning method considering energy network topological characteristics

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Qiang Liu 等.An intelligent optimization method for oil-gas gathering and transportation pipeline network layout.《2016 Chinese Control and Decision Conference (CCDC) IEEE, 2016》.2016,4621-4626. *
俞铁铭 等.天然气分布式能源项目现状分析及建议.《电力与能源》.2016,第37卷(第5期),587-589. *
林健良 等.求解欧氏距离斯坦纳最小树的逐步调优法.《2011年青年通信国际会议论文集》.2011,284-288. *
赵太飞 等.树状注水管网中高可靠度的拓扑优化方法研究.《水资源与水工程学报》.2014,第25卷(第4期),44-49. *

Also Published As

Publication number Publication date
CN111832138A (en) 2020-10-27

Similar Documents

Publication Publication Date Title
CN111832138B (en) Regional pipe network topology optimization method
CN103426032B (en) A kind of economic optimization dispatching method of cogeneration units
CN105576657B (en) Random fuzzy continuous power flow modeling method for active power distribution network and power transmission network integrated system
CN109255471A (en) A kind of hot integrated energy system Expansion Planning optimization method of electric-gas-containing wind-powered electricity generation
CN110135631A (en) Electrical integrated energy system dispatching method based on information gap decision theory
CN107609680B (en) Hydraulic working condition optimization scheduling method for multi-heat-source annular centralized heat supply pipe network
CN110620397B (en) Peak regulation balance evaluation method for high-proportion renewable energy power system
CN113642184B (en) Natural gas pipeline operation optimization method under carbon neutralization target
CN109102104A (en) A kind of multi-energy system synergism low-carbon planing method based on carbon emission flow
CN112650888B (en) Regional comprehensive energy system site selection planning method and system based on graph theory
CN116402210A (en) Multi-objective optimization method, system, equipment and medium for comprehensive energy system
CN113379565A (en) Comprehensive energy system optimization scheduling method based on distributed robust optimization method
CN117035202B (en) Double-layer collaborative expansion planning method for electric heating comprehensive energy system considering demand response
Yang et al. Collaborative planning of integrated natural gas and power supply system considering P2G technique
CN114254890A (en) Regional energy network resource reallocation method based on supply and demand node adjustment
CN114048908A (en) Multi-time scale strengthening interval optimization method for unified power grid-distributed heat supply network system
CN116341881B (en) Robust advanced scheduling method and system for electric-thermal system considering flexibility of heat supply network
CN114266382A (en) Two-stage optimal scheduling method for cogeneration system considering thermal inertia
CN113240204A (en) Energy station capacity optimal configuration method and system considering renewable energy consumption area
CN111931977A (en) Virtual power plant extension planning model construction method considering electric-heat energy transmission value input
CN111768036A (en) Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid
CN109472716A (en) A kind of production capacity end power distribution method multifactor based on the considerations of genetic algorithm
CN112116131B (en) Multi-level optimization method for comprehensive energy system considering carbon emission
CN109948827A (en) It is a kind of based on the flow-optimized method of industrial user's multipotency provided multiple forms of energy to complement each other
Feng et al. Optimal design of pipe network of district cooling system based on genetic algorithm

Legal Events

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