CN103605837A - Multi-heat-source looped pipe network accessibility analysis and determination method, and multi-heat-source looped pipe network accessibility improvement method - Google Patents

Multi-heat-source looped pipe network accessibility analysis and determination method, and multi-heat-source looped pipe network accessibility improvement method Download PDF

Info

Publication number
CN103605837A
CN103605837A CN201310537328.2A CN201310537328A CN103605837A CN 103605837 A CN103605837 A CN 103605837A CN 201310537328 A CN201310537328 A CN 201310537328A CN 103605837 A CN103605837 A CN 103605837A
Authority
CN
China
Prior art keywords
network
matrix
accessibility
pipeline section
node
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.)
Granted
Application number
CN201310537328.2A
Other languages
Chinese (zh)
Other versions
CN103605837B (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.)
Henan University of Technology
Original Assignee
Henan University of Technology
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 Henan University of Technology filed Critical Henan University of Technology
Priority to CN201310537328.2A priority Critical patent/CN103605837B/en
Publication of CN103605837A publication Critical patent/CN103605837A/en
Application granted granted Critical
Publication of CN103605837B publication Critical patent/CN103605837B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Pipeline Systems (AREA)

Abstract

The invention discloses a multi-heat-source looped pipe network accessibility analysis and determination method, and a multi-heat-source looped pipe network accessibility improvement method. According to the analysis and determination method, on the premise that parameters and operating condition of a looped pipe network structure is determined, accessibility of the looped pipe network is determined under this condition by quantitatively analyzing whether or not circulating heads provided by heat sources can meet available heads required by a heat supply network to distribute predetermined flows; according to results of accessibility analysis and determination, the accessibility of the looped pipe network can simultaneously meet the needs of various operating conditions by adjusting structure parameters of the looped pipe network, such as pipe segment diameter and circulating water pump lift, in a certain sequence. Accessibility of the looped pipe network and main factors influencing the accessibility are visually judged by the established accessibility analysis method for the multi-heat-source looped pipe network; after final pipe diameter adjustment and circulating water pump adjustment, the accessibility of the looped pipe network can be satisfactory under any condition.

Description

Many thermals source looping network accessibility is analyzed decision method and improvement method thereof
Technical field
The present invention relates to judgement and the improvement method of central heating system accessibility, the accessibility that relates in particular to many thermals source looping network central heating system is judged and improvement method.
Background technology
Many thermals source looping network is mainly by the operating mechanism of " according to the variation of heating season thermal load; thermal source drops into and be out of service in certain sequence by stages, in batches ", realize the preferential oepration at full load of base-load heat source that efficiency of energy utilization is high, operating cost is low, and the thermal source that efficiency of energy utilization is low, operating cost is high is only used as adjusting thermal source when thermal load is larger, to make up the deficiency of base-load heat source heat capacity, thereby realize the economy of many thermals source looping network heating system.According to the thermal source operation conditions of many thermals source looping network, its operating condition can be divided into base regime, regulate operating mode and design conditions.Under base regime, only have base-load heat source to put into operation; Regulating under operating mode has one or more to regulate thermal source respectively or participation operation simultaneously; Under design conditions, on looping network, all thermals source all put into operation.Therefore, under different operating conditions, the hydrodynamic parameters of many thermals source looping network, as the flow of each pipeline section, pipeline section flow velocity,, all can there is larger variation, thereby cause the often also larger variation of corresponding generation of the accessibility of certain customers under different operating conditions in each node pressure, node flow etc., and under different operating conditions, the accessibility of many thermals source looping network is often difficult to meet simultaneously.The accessibility analysis of so-called pipe network, refers to that pipe network is under given customer flow, analyzes pipe network and whether has the ability that realizes this assignment of traffic.The pipe network structure of single thermal source branched network is simple, and its accessibility analysis ratio is easier to realize, and many thermals source ring-type is due to the complicacy of pipe network structure and network hydraulic operation state, still there is no at present rational accessibility analysis and improvement method.
With the immediate prior art of the application be < < Heating,Ventilating and Air Conditioning (HVAC) > > the 29th the 1st phase of volume disclosed " the accessibility analysis of a central heating net " literary composition in 1999.Following technology is disclosed in this article, the technical measures that adopt " on some pipeline section of looping network, valve being set; also by optimizing these valve openings of calculative determination ", change the hydraulic condition of some the unreachable user on looped network, thereby improve its accessibility.According to the concept of pipe network accessibility, under the accessibility of the many thermals source looping network principle that should to be looping network distribute by rated flow each user, pipe network can be realized ability and the state of this distribution, be pipe network be build-in attribute.Only have can and pipe network, could by the mode regulating, solve the imbalance of pipe network, realize all users can and; And untouchable pipe network regulates its valve opening in any case, all can not guarantee all use realize per family can and.Therefore, the accessibility improvement method of the looping network proposing in document, is actually a kind of improvement method of pipe network hydraulic misadjustment situation, and does not belong to the category that improves pipe network accessibility.
Summary of the invention
In order to overcome existing central heating net accessibility analytical technology, can not carry out to the accessibility of many thermals source looping network the technological deficiency of Correct Analysis judgement, the accessibility that the invention provides a kind of rational many thermals source looping network is analyzed decision method, this analysis decision method is based upon pipe network structure parameter and (comprises pipeline section caliber, water circulating pump performance parameter) and in the definite situation of operating condition, can the circulating head that can provide by each thermal source of quantitative test meet heat supply network and under this operating mode, realize the required available pressure head of predetermined assignment of traffic, judge the accessibility of pipe network under this operating mode.
The present invention provides the improvement method of many thermals source looping network accessibility simultaneously, the method is analyzed result of determination according to many thermals source looping network in the accessibility of various operating conditions, by the structural parameters of looping network (pipeline section caliber, circulating pump head etc.) are adjusted in certain sequence, make its accessibility can meet the requirement of various operating conditions simultaneously.
For achieving the above object, the technical scheme that accessibility analysis decision method of the present invention is taked realizes on computers, and detailed step is as follows:
The first step: determine the operating mode that need to carry out accessibility assay, and many thermals source looping network is decomposed into for water ring net and backwater looped network, draw respectively the network calculations figure that supplies water ring net and backwater looped network under this operating mode; Select the confession backwater looped network node of the base-load heat source that is provided with constant pressure arrangement as the pressure reference point for water ring net and backwater looped network, and respectively the node of water supply network and return pipe net and pipeline section are carried out to corresponding numbering, wherein:
For the node serial number of water ring net G1 respectively, G2 ... represent, the respective nodes numbering of backwater looped network is used respectively H1, H2 ... represent; Pipeline section for water ring net is numbered 1G respectively, 2G ..., the respective tube segment number of backwater looped network is used respectively 1H, 2H ... represent; During numbering, pressure reference point is numbered maximum numbering; On this basis, write out respectively the connection matrix A for water ring net g, transposition connection matrix A ' g, resistance matrix number
Figure BDA0000407828650000031
node rated flow matrix
Figure BDA0000407828650000032
pipe network circulation area square formation F g, and the connection matrix A of backwater looped network h, transposition connection matrix A ' h, resistance matrix number
Figure BDA0000407828650000033
node rated flow matrix
Figure BDA0000407828650000034
pipe network circulation area square formation F h;
Each element of the above-mentioned traffic matrix for water ring net or backwater looped network, flows into as negative value, flow out on the occasion of;
The above-mentioned pipe network circulation area square formation F for water ring net gor the pipe network circulation area square formation F of backwater looped network hthe symmetry square matrix being formed by the circulation area for water ring net or each pipeline section of backwater looped network, the circulation area that the elements in a main diagonal is each pipeline section
Figure BDA0000407828650000035
d in formula irepresent pipeline section internal diameter, non-the elements in a main diagonal is zero;
Second step: set up many thermals source looping network with respect to the node pressure drop solving model of reference point;
Y &CenterDot; p &OverBar; = q &OverBar;
In formula,
Y---the matrix of coefficients of looping network node pressure drop equation,
Figure BDA0000407828650000037
The connection matrix of A---looping network ((n-1) * b rank matrix, n is nodes, b is pipeline section number);
The transposed matrix of A '---looping network connection matrix A (b * (n-1) rank matrix);
Figure BDA0000407828650000038
---looping network pipeline section resistance consumption square formation,
Figure BDA0000407828650000039
Figure BDA00004078286500000310
---the symmetry square matrix that each pipeline section resistance number of looping network forms, on principal diagonal, element is respectively each pipeline section resistance number, and non-the elements in a main diagonal is 0;
Figure BDA00004078286500000311
---the symmetry square matrix that each pipeline flow absolute value of looping network forms, on principal diagonal, element is respectively each pipeline flow absolute value, and non-the elements in a main diagonal is 0;
Figure BDA0000407828650000041
---the node pressure drop matrix of looping network;
Figure BDA0000407828650000042
Figure BDA0000407828650000043
---each node rated flow matrix of looping network,
Figure BDA0000407828650000044
The 3rd step: solve the node pressure drop matrix for water ring net
Figure BDA0000407828650000045
The node pressure drop matrix of the solving model of application second step to water supply looping network
Figure BDA0000407828650000046
while solving, due to the node rated flow matrix of water supply network
Figure BDA0000407828650000047
known quantity, and its matrix of coefficients Y gbe by
Figure BDA0000407828650000048
determine as calculated, because pipeline flow is unknown quantity, its matrix of coefficients Y galso be unknown quantity; Therefore, utilize above-mentioned equation direct solution node pressure drop matrix
Figure BDA0000407828650000049
impossible; In this case, utilize following iterative algorithm can obtain node pressure drop matrix
Figure BDA00004078286500000410
detailed algorithm is as follows:
The 3.1st step: the initial flow rate matrix of setting arbitrarily pipe network
Figure BDA00004078286500000411
each element of matrix represents each pipeline section initial flow rate, according to try to achieve the initial flow matrix of pipeline section and the symmetry square matrix of corresponding pipeline section initial flow formation
Above-mentioned
Figure BDA00004078286500000415
represent pipe network traffic matrix, element is each pipeline flow;
Above-mentioned F represents pipe network circulation area square formation, for water ring net, is F g, return pipe net is F h;
Above-mentioned represent pipe network flow velocity matrix, element is each pipeline section flow velocity;
The 3.2nd step: according to formula s j = 8 ( &lambda; j l j d j + &Sigma; i &zeta; j ) &pi; 2 &CenterDot; d j 4 &CenterDot; g , &lambda; i = 0.11 ( K d j + 68 Re j ) 0.25 , Re j = v j d j &upsi; , Determine matrix
Figure BDA00004078286500000418
and then according to
Figure BDA00004078286500000419
determine pipe network resistance consumption matrix
Figure BDA00004078286500000420
s in formula j, l j, d j,
Figure BDA00004078286500000421
λ i, Re j, v jrepresent respectively pipeline section resistance number, length of pipe section, pipeline section internal diameter, the total coefficient of shock resistance of pipeline section, pipeline section coefficient of frictional resistance, pipeline section Reynolds number and the pipeline section flow velocity of j pipeline section; K represents the absolute coefficient of roughness of pipeline; υ represents the kinematic viscosity of fluid;
The 3.3rd step: according to
Figure BDA00004078286500000422
determine the matrix of coefficients Y of node pressure drop equation 0;
The 3.4th step: utilize formula determine the node pressure drop matrix with respect to pressure reference point
Figure BDA00004078286500000424
The 3.5th step: utilize formula
Figure BDA00004078286500000425
determine pipeline section pressure drop matrix
The 3.6th step: utilize formula
Figure BDA00004078286500000519
determine pipeline flow matrix
Figure BDA0000407828650000051
The 3.7th step: according to
Figure BDA0000407828650000052
determine pipeline section flow velocity matrix
Figure BDA0000407828650000053
The 3.8th step: relatively
Figure BDA0000407828650000054
check whether its precision meets the demands;
Preset accuracy requirement as follows:
Figure BDA0000407828650000055
the corresponding element of two matrixes subtracts each other in the matrix of differences of gained, if the absolute value of the element of absolute value maximum is less than or equal to 0.01, thinks that precision meets the demands; As precision meets the demands, iteration finishes; If do not met, will
Figure BDA0000407828650000056
each element value is replaced accordingly each element value, repeats 3.1-3.7 step process, until meet predefined accuracy requirement;
The 3.9th step: output pipeline section flow velocity matrix
Figure BDA0000407828650000058
pipeline flow matrix
Figure BDA0000407828650000059
pipeline section pressure drop matrix
Figure BDA00004078286500000510
and node pressure drop matrix
Figure BDA00004078286500000511
Above-mentioned iterative process finishes the node pressure drop matrix obtaining be the node pressure drop matrix for water ring net
Figure BDA00004078286500000513
The 4th step: the node pressure drop matrix that solves backwater looped network
Adopt and the node pressure drop matrix solving for water ring net
Figure BDA00004078286500000515
identical alternative manner, can obtain the node pressure drop matrix of backwater looped network
Figure BDA00004078286500000516
the matrix relating in iterative process all adopts the corresponding matrix of return pipe net;
The 5th step: the least favorable node of determining many thermals source looping network;
According to the water supply looped network node pressure drop matrix of many thermals source looping network
Figure BDA00004078286500000517
and backwater looped network node pressure drop matrix
Figure BDA00004078286500000518
utilize following formula can determine the least favorable node J of many thermals source looping network;
p JG-p JH=max(p iG-p iH)i=1,2,…,n
P in formula jG---the water supply looped network node pressure drop of least favorable node J, mH 2o;
P jH---the backwater looped network node pressure drop of least favorable node J, mH 2o;
The node number of n---many thermals source looping network;
The 6th step: the required available pressure head H providing of each thermal source of many thermals source looping network is provided rk;
Due to the accessibility analysis of many thermals source looping network mainly with the confession of base-load heat source, backwater looped network node as reference point, therefore, the user one that reference point is corresponding is decided to be base-load heat source, the required available pressure head providing of base-load heat source can be determined by following formula:
H r1=Δp 0+(p JG-p JH)+H 0
H in formula r1---the required available pressure head providing of base-load heat source, mH 2o;
Δ p 0---the available pressure head of minimum permission (pressing regulation and stipulation value) of least favorable node J, mH 2o;
H 0---the drag losses of thermal source inside under rated flow (pressing regulation and stipulation value), mH 2o;
For other thermal source, its required available pressure head providing can be determined by following formula:
H rk=Δp 0+(p JG-p JH)-(p iG-p iH)+H 0k=2,…,M
H in formula rk---the available pressure head that k thermal source need provide, mH 2o;
P iG---the node pressure drop of the water supply looped network node that k thermal source is corresponding, mH 2o;
P iH---the node pressure drop of the backwater looped network node that k thermal source is corresponding, mH 2o;
M---the number of many thermals source looping network operation thermal source;
The 7th step: the accessibility of judging heating network;
If the circulating head that many thermals source looping network thermal source can provide is H ' r1or H ' rk, wherein the former represents the circulating head that base-load heat source can provide, the latter represents the circulating head that other thermal source can provide, the available pressure head H that this thermal source of determining according to the 6th step need to provide r1or H rkif, H ' r1>=H r1or H ' rk>=H rk, this thermal source meets the accessibility of pipe network under this operating mode; Otherwise this thermal source does not meet the accessibility of pipe network under this operating mode.Can judge thus whether all thermals source meet the accessibility of pipe network under this operating mode.
Above-mentioned analysis judges that the step of many thermals source looping network accessibility can analyze judgement in the accessibility of any single operating mode to many thermals source looping network.But many thermals source looping network has multiple operating condition, and from the Investigation and analysis to existing many thermals source looping network, pipe network is often difficult to meet in the accessibility of different operating conditions simultaneously.And utilize following improvement method, can realize the accessibility of pipe network under different operating modes and all meet the demands.
When the accessibility of many thermals source looping network is improved, first judge that pipe network meets the complexity of accessibility under different operating conditions, then according to this complexity, determine that accessibility is judged and the operating mode order of steam pipeline adjustment.Generally, the thermal source putting into operation is more, and under this operating mode, the accessibility of pipe network more easily meets.Therefore, should first the accessibility of pipe network under design conditions be judged and be adjusted, secondly, to regulating the accessibility of operating mode judge and adjust, finally the accessibility of pipe network under base regime is judged and adjusted.Given this, the step of accessibility improvement method of the present invention is as follows:
The first step: the accessibility to many thermals source looping network under design conditions is judged and adjusts;
After pipe network accessibility to many thermals source looping network under design conditions is analyzed and judged, according to the flow speed value of each pipeline section, the pipeline section that first flow velocity is greater than to 0.8m/s is adjusted, and by increasing the caliber of these pipeline sections, makes its flow velocity be not more than 0.8m/s; And then the accessibility of the pipe network after adjusting is judged, if still have the flow velocity of part pipeline section to be greater than 0.8m/s, need again to corresponding pipeline section, adjust, so repeatedly, until the flow velocity of all pipeline sections is not more than 0.8m/s under design conditions;
Second step: to regulating the accessibility of many thermals source looping network under operating mode judge and adjust;
On the basis of the pipe network after above-mentioned adjustment, according to the method for the first step, to many thermals source looping network, regulating the accessibility under operating mode analyze judgement and the caliber not meeting the demands is adjusted.By repeatedly judging and adjusting, make each pipeline section of pipe network regulate the flow velocity under operating mode to be not more than 0.8m/s;
The 3rd step: the accessibility to many thermals source looping network under base regime is judged and adjusts;
On the basis that regulates the adjustment of operating mode pipeline section to complete, still by the method for the first step, many thermals source looping network is analyzed and judged and the caliber not meeting the demands is adjusted in the accessibility under base regime.By repeatedly judging and adjusting, make each pipeline section of pipe network regulate the flow velocity under operating mode to be not more than 0.8m/s;
The 4th step: the adjustment to pipe network water circulating pump under base regime;
To having passed through the whole pipe network of above-mentioned three steps, the accessibility analytical approach of many thermals source looping network that recycling is introduced above, many thermals source looping network is analyzed to judgement in the accessibility under base regime, according to analyzing result of determination, many thermals source looping network base-load heat source required available pressure head providing under base regime can be provided, if the circulating head that base-load heat source can provide under this operating mode is less than the required available pressure head of base regime, should change the corresponding water circulating pump of base-load heat source, make under base regime, the circulating head that base-load heat source can provide is not less than its available pressure head requirement required under base regime,
The 5th step: to regulating the adjustment of pipe network water circulating pump under operating mode;
On the basis completing in the adjustment of base regime water circulating pump, to many thermals source looping network, regulating the accessibility under operating mode to analyze judgement; Then according to the method for the 4th step, judge and adjust the water circulating pump of each thermal source putting into operation under this operating mode, make regulating under operating mode, the circulating head that pipe network respectively moves the water circulating pump of thermal source can be provided is not less than this thermal source and is regulating required available pressure head requirement under operating mode;
The 6th step: the adjustment to pipe network water circulating pump under design conditions;
On the basis that regulates the adjustment of Operation mode cycle water pump to complete, many thermals source looping network is analyzed to judgement in the accessibility under design conditions; Then according to the method for the 4th step, judge and adjust the water circulating pump of each thermal source putting into operation under this operating mode, make under design conditions, the circulating head that pipe network respectively moves the water circulating pump of thermal source can be provided is not less than this thermal source in the required available pressure head requirement of design conditions.
Advantage of the present invention is:
(1) accessibility of many thermals source looping network of setting up analyze decision method with each node pressure head of pipe network for solving variable, compare as solving the method for solving of variable hydraulic parameters with traditional pipeline flow of take, solving equation number is relatively less, do not need complicated adjustment process, solve speed, only need the iteration of a small amount of number of times, can obtain satisfied convergence result.
(2) the accessibility analytical approach of many thermals source looping network of utilize setting up, can obtain the results such as available pressure that flow, flow velocity, drag losses and each thermal source of each pipeline section under various operating modes must provide.According to these results, can judge intuitively the accessibility of pipe network and affect the principal element of accessibility.
(3) adjustment to many thermals source looping network caliber, can be from accessibility the most satisfied operating mode, adopt the strategy of " limit is analyzed; limit is adjusted ", successively the accessibility under various operating modes is carried out to analysis and adjustment (the accessibility analysis of rear a kind of operating mode judges it is all to carry out on the basis of having carried out adjusting at the pipe network caliber to front a kind of operating mode).Thereby, often once adjust (caliber of part pipeline section increases), the accessibility situation of pipe network will inevitably make moderate progress.Therefore,, after last caliber adjustment and water circulating pump adjustment, the accessibility of pipe network under various operating modes all can meet the demands.
Accompanying drawing explanation
Fig. 1 is the water supply network network diagram of a certain many thermals source looping network system;
Fig. 2 is the return pipe net network diagram of embodiment in Fig. 1;
Fig. 3 is that accessibility of the present invention is analyzed determination methods flow chart;
Fig. 4 is that accessibility of the present invention is improved method program block diagram.
In Fig. 1: 1G-19G represents the pipeline section numbering of water supply network;
G1-G17 represents the node serial number of water supply network, the numbering that wherein G17 is pressure reference point.
In Fig. 2: 1H-19H represents the pipeline section numbering of return pipe net;
H1-H17 represents the node serial number of return pipe net, the numbering that wherein 17H is pressure reference point.
Embodiment
Accessibility of the present invention is analyzed the flow chart of determination methods and accessibility improvement method and is seen respectively Fig. 3 and Fig. 4, and its detailed implementation step has been done detailed record in summary of the invention, is not described in detail in this.
The a certain many thermals source looping network system of take is below example, illustrates according to the inventive method and carries out accessibility judgement and improve the process of adjusting.
A certain many thermals source looping network heating system, pipe network comprises three thermals source and ten hot users (heat exchange station), under base regime, only has base-load heat source 1 to put into operation; Regulating under operating mode base-load heat source 1 and adjusting thermal source 2 cooperations; Under design conditions, base-load heat source 1, adjusting thermal source 2 and adjusting thermal source 3 all put into operation.Base-load heat source 1 arranges a water circulating pump, and its rated flow is 294m 3/ h, rated head is 65mH 2o; Regulate thermal source 2 that a water circulating pump is set, its rated flow is 214m 3/ h, rated head is 77mH 2o; Regulate thermal source 3 that a water circulating pump is set, its rated flow is 198m 3/ h, rated head is 74mH 2o.The network structure of this many thermals source looping network as shown in Figures 1 and 2.As can be seen from Figure 2, base-load heat source 1 is from node G17(H17) network, regulate thermal source 2 from node G9(H9) network, regulate thermal source 3 from node G11(H11) network.The caliber of the pipeline section 1G-19G of this many thermals source looping network and pipe range are in Table 1; The rated flow of pipe network node G1-G17 under three kinds of operating modes is in Table 2.The pipeline section 1H-19H caliber of return pipe net and pipe range are identical with caliber and the pipe range of corresponding pipeline section 1G-19G, and the rated flow of return pipe net node G1-G17 under three kinds of operating modes equate with the rated flow numerical value of corresponding node G1-G17, and symbol is contrary.
The caliber of each pipeline section of table 1 water supply network
Figure BDA0000407828650000101
Each node rated flow of water supply network under tri-kinds of operating modes of table 2()
Figure BDA0000407828650000111
The decision method of the pipe network accessibility that summary of the invention is recorded according to the present invention, can obtain following result of determination:
Under design conditions, all use of many thermals source looping network per family can and; Regulating under operating mode, the certain customers in base-load heat source 1 service range are unreachable, regulate use in thermal source 2 service ranges per family can and; Under base regime, only have base-load heat source 1 operation, unreachable compared with multi-user.
Accessibility method of adjustment by introduction of the present invention is adjusted as follows to above-mentioned pipe network, can meet the accessibility of various users under different operating modes:
(1) pipeline section adjustment: the caliber of pipeline section 3G and 3H is adjusted into Ф 219 * 4.5 by Ф 159 * 4.5; The caliber of pipeline section 4G and 4H is adjusted into Ф 219 * 4.5 by Ф 73 * 3.5; The caliber of pipeline section 6G and 6H is adjusted into Ф 219 * 4.5 by Ф 73 * 3.5; The caliber of pipeline section 4G and 4H is adjusted into Ф 159 * 4.5 by Ф 73 * 3.5.
(2) water circulating pump adjustment: under design conditions, the water circulating pump of base-load heat source 1, adjusting thermal source 2, adjusting thermal source 3 all meets accessibility requirement; Regulate under operating mode, base-load heat source 1 need be set up a water circulating pump moving under this operating mode, and its type selecting flow is constant, and lift is 68mH 2o, regulates thermal source 2 need set up a water circulating pump moving under this operating mode, and its type selecting flow is constant, and lift is 58mH 2o; Under base regime, base-load heat source 1 is still needed and set up a water circulating pump moving under this operating mode, and its type selecting flow is constant, and lift is 85mH 2o.

Claims (2)

1. the looping network of thermal source more than accessibility is analyzed a decision method, it is characterized in that, this analytical approach realizes on computers, and detailed step is as follows:
The first step: determine the operating mode that need to carry out accessibility assay, and many thermals source looping network is decomposed into for water ring net and backwater looped network, draw respectively the network calculations figure that supplies water ring net and backwater looped network under this operating mode; Select the confession backwater looped network node of the base-load heat source that is provided with constant pressure arrangement as the pressure reference point for water ring net and backwater looped network, and respectively the node of water supply network and return pipe net and pipeline section are carried out to corresponding numbering, wherein:
For the node serial number of water ring net G1 respectively, G2 ... represent, the respective nodes numbering of backwater looped network is used respectively H1, H2 ... represent; Pipeline section for water ring net is numbered 1G respectively, 2G ..., the respective tube segment number of backwater looped network is used respectively 1H, 2H ... represent; During numbering, pressure reference point is numbered maximum numbering; On this basis, write out respectively the connection matrix A for water ring net g, transposition connection matrix A ' g, resistance matrix number node rated flow matrix
Figure FDA0000407828640000016
pipe network circulation area square formation F g, and the connection matrix A of backwater looped network h, transposition connection matrix A ' h, resistance matrix number
Figure FDA0000407828640000017
node rated flow matrix
Figure FDA0000407828640000018
pipe network circulation area square formation F h;
Each element of the above-mentioned traffic matrix for water ring net or backwater looped network, flows into as negative value, flow out on the occasion of;
The above-mentioned pipe network circulation area square formation F for water ring net gor the pipe network circulation area square formation F of backwater looped network hthe symmetry square matrix being formed by the circulation area for water ring net or each pipeline section of backwater looped network, the circulation area that the elements in a main diagonal is each pipeline section
Figure FDA0000407828640000011
d in formula irepresent pipeline section internal diameter, non-the elements in a main diagonal is zero;
Second step: set up many thermals source looping network with respect to the node pressure drop solving model of reference point;
Y &CenterDot; p &OverBar; = q &OverBar;
In formula,
Y---the matrix of coefficients of looping network node pressure drop equation,
Figure FDA0000407828640000013
The connection matrix of A---looping network;
A '---the transposed matrix of looping network connection matrix A;
---looping network pipeline section resistance consumption square formation,
Figure FDA0000407828640000014
Figure FDA0000407828640000021
---the symmetry square matrix that each pipeline section resistance number of looping network forms, on principal diagonal, element is respectively each pipeline section resistance number, and non-the elements in a main diagonal is 0;
Figure FDA0000407828640000022
---the symmetry square matrix that each pipeline flow absolute value of looping network forms, on principal diagonal, element is respectively each pipeline flow absolute value, and non-the elements in a main diagonal is 0;
Figure FDA0000407828640000023
---the node pressure drop matrix of looping network;
Figure FDA0000407828640000024
Figure FDA0000407828640000025
---each node rated flow matrix of looping network,
The 3rd step: solve the node pressure drop matrix for water ring net
Figure FDA0000407828640000027
Utilize following iterative algorithm to obtain node pressure drop matrix
Figure FDA0000407828640000028
detailed algorithm is as follows:
The 3.1st step: the initial flow rate matrix of setting arbitrarily pipe network
Figure FDA0000407828640000029
each element of matrix represents each pipeline section initial flow rate, according to
Figure FDA00004078286400000210
try to achieve the initial flow matrix of pipeline section
Figure FDA00004078286400000211
and the symmetry square matrix of corresponding pipeline section initial flow formation
Figure FDA00004078286400000212
Above-mentioned
Figure FDA00004078286400000213
represent pipe network traffic matrix, element is each pipeline flow;
Above-mentioned F represents pipe network circulation area square formation, for water ring net, is F g, return pipe net is F h;
Above-mentioned
Figure FDA00004078286400000214
represent pipe network flow velocity matrix, element is each pipeline section flow velocity;
The 3.2nd step: according to formula s j = 8 ( &lambda; j l j d j + &Sigma; i &zeta; j ) &pi; 2 &CenterDot; d j 4 &CenterDot; g , &lambda; i = 0.11 ( K d j + 68 Re j ) 0.25 , Re j = v j d j &upsi; , Determine matrix
Figure FDA00004078286400000216
and then according to
Figure FDA00004078286400000217
determine pipe network resistance consumption matrix
Figure FDA00004078286400000218
s in formula j, l j, d j,
Figure FDA00004078286400000219
λ i, Re j, v jrepresent respectively pipeline section resistance number, length of pipe section, pipeline section internal diameter, the total coefficient of shock resistance of pipeline section, pipeline section coefficient of frictional resistance, pipeline section Reynolds number and the pipeline section flow velocity of j pipeline section; K represents the absolute coefficient of roughness of pipeline; υ represents the kinematic viscosity of fluid;
The 3.3rd step: according to determine the matrix of coefficients Y of node pressure drop equation 0;
The 3.4th step: utilize formula
Figure FDA00004078286400000221
determine the node pressure drop matrix with respect to pressure reference point
Figure FDA00004078286400000222
The 3.5th step: utilize formula determine pipeline section pressure drop matrix
Figure FDA00004078286400000224
The 3.6th step: utilize formula
Figure FDA0000407828640000031
determine pipeline flow matrix
Figure FDA0000407828640000032
The 3.7th step: according to determine pipeline section flow velocity matrix
Figure FDA0000407828640000034
The 3.8th step: relatively
Figure FDA0000407828640000035
check whether its precision meets the demands;
Preset accuracy requirement as follows: the corresponding element of two matrixes subtracts each other in the matrix of differences of gained, if the absolute value of the element of absolute value maximum is less than or equal to 0.01, thinks that precision meets the demands; As precision meets the demands, iteration finishes; If do not met, will each element value is replaced accordingly
Figure FDA0000407828640000038
each element value, repeats 3.1-3.7 step process, until meet predefined accuracy requirement;
The 3.9th step: output pipeline section flow velocity matrix
Figure FDA0000407828640000039
pipeline flow matrix
Figure FDA00004078286400000310
pipeline section pressure drop matrix
Figure FDA00004078286400000311
and node pressure drop matrix
Figure FDA00004078286400000312
Above-mentioned iterative process finishes the node pressure drop matrix obtaining
Figure FDA00004078286400000313
be the node pressure drop matrix for water ring net
Figure FDA00004078286400000314
The 4th step: the node pressure drop matrix that solves backwater looped network
Figure FDA00004078286400000315
Adopt and the node pressure drop matrix solving for water ring net
Figure FDA00004078286400000316
identical alternative manner, can obtain the node pressure drop matrix of backwater looped network the matrix relating in iterative process all adopts the corresponding matrix of return pipe net;
The 5th step: the least favorable node of determining many thermals source looping network;
According to the water supply looped network node pressure drop matrix of many thermals source looping network
Figure FDA00004078286400000318
and backwater looped network node pressure drop matrix
Figure FDA00004078286400000319
utilize following formula can determine the least favorable node J of many thermals source looping network;
p JG-p JH=max(p iG-p iH)i=1,2,…,n
P in formula jG---the water supply looped network node pressure drop of least favorable node J, mH 2o;
P jH---the backwater looped network node pressure drop of least favorable node J, mH 2o;
The node number of n---many thermals source looping network;
The 6th step: the required available pressure head H providing of each thermal source of many thermals source looping network is provided rk;
Due to the accessibility analysis of many thermals source looping network mainly with the confession of base-load heat source, backwater looped network node as reference point, therefore, the user one that reference point is corresponding is decided to be base-load heat source, the required available pressure head providing of base-load heat source is determined by following formula:
H r1=Δp 0+(p JG-p JH)+H 0
H in formula r1---the required available pressure head providing of base-load heat source, mH 2o;
Δ p 0---the available pressure head of minimum permission of least favorable node J, mH 2o;
H 0---the drag losses of thermal source inside under rated flow, mH 2o;
For other thermal source, its required available pressure head providing can be determined by following formula:
H rk=Δp 0+(p JG-p JH)-(p iG-p iH)+H 0k=2,…,M
H in formula rk---the available pressure head that k thermal source need provide, mH 2o;
P iG---the node pressure drop of the water supply looped network node that k thermal source is corresponding, mH 2o;
P iH---the node pressure drop of the backwater looped network node that k thermal source is corresponding, mH 2o;
M---the number of many thermals source looping network operation thermal source;
The 7th step: the accessibility of judging heating network;
If the circulating head that many thermals source looping network thermal source can provide is H ' r1or H ' rk, wherein the former represents the circulating head that base-load heat source can provide, the latter represents the circulating head that other thermal source can provide, the available pressure head H that this thermal source of determining according to the 6th step need to provide r1or H rkif, H ' r1>=H r1or H ' rk>=H rk, this thermal source meets the accessibility of pipe network under this operating mode; Otherwise this thermal source does not meet the accessibility of pipe network under this operating mode; Judge thus all thermals source and under this operating mode, whether meet the accessibility of pipe network.
2. utilize the accessibility described in claim 1 to analyze the method that decision method improves many thermals source looping network accessibility, it is characterized in that,
The first step: the accessibility to many thermals source looping network under design conditions is judged and adjusts;
After pipe network accessibility to many thermals source looping network under design conditions is analyzed and judged, according to the flow speed value of each pipeline section, the pipeline section that first flow velocity is greater than to 0.8m/s is adjusted, and by increasing the caliber of these pipeline sections, makes its flow velocity be not more than 0.8m/s; And then the accessibility of the pipe network after adjusting is judged, if still have the flow velocity of part pipeline section to be greater than 0.8m/s, need again to corresponding pipeline section, adjust, so repeatedly, until the flow velocity of all pipeline sections is not more than 0.8m/s under design conditions;
Second step: to regulating the accessibility of many thermals source looping network under operating mode judge and adjust;
On the basis of the pipe network after above-mentioned adjustment, according to the method for the first step, to many thermals source looping network, regulating the accessibility under operating mode analyze judgement and the caliber not meeting the demands is adjusted.By repeatedly judging and adjusting, make each pipeline section of pipe network regulate the flow velocity under operating mode to be not more than 0.8m/s;
The 3rd step: the accessibility to many thermals source looping network under base regime is judged and adjusts;
On the basis that regulates the adjustment of operating mode pipeline section to complete, still by the method for the first step, many thermals source looping network is analyzed and judged and the caliber not meeting the demands is adjusted in the accessibility under base regime.By repeatedly judging and adjusting, make each pipeline section of pipe network regulate the flow velocity under operating mode to be not more than 0.8m/s;
The 4th step: the adjustment to pipe network water circulating pump under base regime;
To having passed through the whole pipe network of above-mentioned three steps, the accessibility analytical approach of many thermals source looping network that recycling is introduced above, many thermals source looping network is analyzed to judgement in the accessibility under base regime, according to analyzing result of determination, many thermals source looping network base-load heat source required available pressure head providing under base regime can be provided, if the circulating head that base-load heat source can provide under this operating mode is less than the required available pressure head of base regime, should change the corresponding water circulating pump of base-load heat source, make under base regime, the circulating head that base-load heat source can provide is not less than its available pressure head requirement required under base regime,
The 5th step: to regulating the adjustment of pipe network water circulating pump under operating mode;
On the basis completing in the adjustment of base regime water circulating pump, to many thermals source looping network, regulating the accessibility under operating mode to analyze judgement; Then according to the method for the 4th step, judge and adjust the water circulating pump of each thermal source putting into operation under this operating mode, make regulating under operating mode, the circulating head that pipe network respectively moves the water circulating pump of thermal source can be provided is not less than this thermal source and is regulating required available pressure head requirement under operating mode;
The 6th step: the adjustment to pipe network water circulating pump under design conditions;
On the basis that regulates the adjustment of Operation mode cycle water pump to complete, many thermals source looping network is analyzed to judgement in the accessibility under design conditions; Then according to the method for the 4th step, judge and adjust the water circulating pump of each thermal source putting into operation under this operating mode, make under design conditions, the circulating head that pipe network respectively moves the water circulating pump of thermal source can be provided is not less than this thermal source in the required available pressure head requirement of design conditions.
CN201310537328.2A 2013-11-04 2013-11-04 Multi-heat source looping network accessibility analyzes decision method and ameliorative way thereof Expired - Fee Related CN103605837B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310537328.2A CN103605837B (en) 2013-11-04 2013-11-04 Multi-heat source looping network accessibility analyzes decision method and ameliorative way thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310537328.2A CN103605837B (en) 2013-11-04 2013-11-04 Multi-heat source looping network accessibility analyzes decision method and ameliorative way thereof

Publications (2)

Publication Number Publication Date
CN103605837A true CN103605837A (en) 2014-02-26
CN103605837B CN103605837B (en) 2016-08-17

Family

ID=50124058

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310537328.2A Expired - Fee Related CN103605837B (en) 2013-11-04 2013-11-04 Multi-heat source looping network accessibility analyzes decision method and ameliorative way thereof

Country Status (1)

Country Link
CN (1) CN103605837B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103939968A (en) * 2014-04-15 2014-07-23 河南理工大学 Method for arranging circulating water pumps of multi-heat-source looped network heating supply system and adjusting pipe network parameters
CN104462777A (en) * 2014-11-13 2015-03-25 浙江工业大学 Fluid pipe minimum feasible energy requirement measuring method
CN105160062A (en) * 2015-07-24 2015-12-16 北京市建筑工程研究院有限责任公司 Hydraulic check method for reverse return pipe network
CN105160063A (en) * 2015-07-24 2015-12-16 北京市建筑工程研究院有限责任公司 Hydraulic check method for direct return pipe network
CN107885085A (en) * 2017-11-14 2018-04-06 宁波大学 Complicated conduit running control method based on deep learning
CN109785185A (en) * 2019-01-09 2019-05-21 清华大学 The level of factory complex energy management system of more grade heat supplies is dispatched in steam power plant
CN110473118A (en) * 2019-08-14 2019-11-19 华北电力大学(保定) A kind of generation method of multi-heat source tree-shaped heating network incidence matrix
CN111379966A (en) * 2018-12-27 2020-07-07 中国石油化工股份有限公司 Memory, hydrogen pipe network working condition data generation method, device and equipment
CN111396986A (en) * 2020-03-26 2020-07-10 河南理工大学 Impedance-based manual static hydraulic balance adjusting method for central heating pipe network
CN112417632A (en) * 2020-12-01 2021-02-26 南京益美沃电子科技有限公司 Simulation calculation method and system for town gas pipe network

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101211169A (en) * 2007-12-21 2008-07-02 哈尔滨工业大学 Minimum conveying energy consumption heat supply regulating method
CN101846355A (en) * 2010-02-26 2010-09-29 新疆孚德节能工程有限公司 Arrangement method based on dynamic balance unit technology in hot water heating system
CN102261691A (en) * 2011-05-23 2011-11-30 河南理工大学 Power transmission system of multi-heat source loop pipe network, system configuration method and operation mode

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101211169A (en) * 2007-12-21 2008-07-02 哈尔滨工业大学 Minimum conveying energy consumption heat supply regulating method
CN101846355A (en) * 2010-02-26 2010-09-29 新疆孚德节能工程有限公司 Arrangement method based on dynamic balance unit technology in hot water heating system
CN102261691A (en) * 2011-05-23 2011-11-30 河南理工大学 Power transmission system of multi-heat source loop pipe network, system configuration method and operation mode

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
XU WENZHONG,ET AL.: "Genetic Algorithm’s Optimization for Parameters of Ring-shaped Heat-supply Network with Multi-heat Sources", 《2009 INTERNATIONAL CONFERENCE ON ENERGY AND ENVIRONMENT TECHNOLOGY》 *
徐宝萍: "大型多热源环状热水网水力计算方法", 《区域供热》 *
徐文忠,等。: "遗传算法在多热源环状管网优化中的应用研究", 《2011 INTERNATIONAL CONFERENCE ON FUTURE COMPUTER SCIENCE AND APPLICATION》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103939968B (en) * 2014-04-15 2016-06-22 河南理工大学 Multi-heat source looped network heating system water circulating pump is arranged and pipeline parameter method of adjustment
CN103939968A (en) * 2014-04-15 2014-07-23 河南理工大学 Method for arranging circulating water pumps of multi-heat-source looped network heating supply system and adjusting pipe network parameters
CN104462777A (en) * 2014-11-13 2015-03-25 浙江工业大学 Fluid pipe minimum feasible energy requirement measuring method
CN105160062A (en) * 2015-07-24 2015-12-16 北京市建筑工程研究院有限责任公司 Hydraulic check method for reverse return pipe network
CN105160063A (en) * 2015-07-24 2015-12-16 北京市建筑工程研究院有限责任公司 Hydraulic check method for direct return pipe network
CN105160063B (en) * 2015-07-24 2018-02-02 北京市建筑工程研究院有限责任公司 A kind of DRS hydraulic pipeline check method
CN105160062B (en) * 2015-07-24 2018-02-02 北京市建筑工程研究院有限责任公司 A kind of same journey hydraulic pipeline check method
CN107885085A (en) * 2017-11-14 2018-04-06 宁波大学 Complicated conduit running control method based on deep learning
CN111379966A (en) * 2018-12-27 2020-07-07 中国石油化工股份有限公司 Memory, hydrogen pipe network working condition data generation method, device and equipment
CN111379966B (en) * 2018-12-27 2021-07-09 中国石油化工股份有限公司 Memory, hydrogen pipe network working condition data generation method, device and equipment
CN109785185A (en) * 2019-01-09 2019-05-21 清华大学 The level of factory complex energy management system of more grade heat supplies is dispatched in steam power plant
CN110473118A (en) * 2019-08-14 2019-11-19 华北电力大学(保定) A kind of generation method of multi-heat source tree-shaped heating network incidence matrix
CN110473118B (en) * 2019-08-14 2023-04-14 华北电力大学(保定) Generation method of multi-heat-source dendritic heat supply pipe network incidence matrix
CN111396986A (en) * 2020-03-26 2020-07-10 河南理工大学 Impedance-based manual static hydraulic balance adjusting method for central heating pipe network
CN112417632A (en) * 2020-12-01 2021-02-26 南京益美沃电子科技有限公司 Simulation calculation method and system for town gas pipe network
CN112417632B (en) * 2020-12-01 2023-03-10 南京益美沃电子科技有限公司 Simulation calculation method and system for town gas pipe network

Also Published As

Publication number Publication date
CN103605837B (en) 2016-08-17

Similar Documents

Publication Publication Date Title
CN103605837A (en) Multi-heat-source looped pipe network accessibility analysis and determination method, and multi-heat-source looped pipe network accessibility improvement method
CN105243438B (en) A kind of probabilistic carry-over storage Optimization Scheduling of consideration runoff
CN107391807B (en) Plate-fin heat exchanger heat transfer flow performance value analogy method based on transient technology
Chen et al. An entransy dissipation-based optimization principle for building central chilled water systems
CN108592165A (en) A kind of heat exchange station system optimal control method
CN102261691B (en) Power transmission system of multi-heat source loop pipe network, system configuration method and operation mode
Wang et al. Experimental study on the dynamics, quality and impacts of using variable-speed pumps in buildings for frequency regulation of smart power grids
Kiss et al. Modeling of an electric vehicle thermal management system in MATLAB/Simulink
CN102095447B (en) Multi-medium combined heat exchanger test bed
Wemhoff et al. Predictions of energy savings in HVAC systems by lumped models
CN111222213A (en) Dynamic simulation method and device for thermodynamic network
CN110555264A (en) Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method
CN103558046A (en) Heat exchanger energy efficiency evaluation system
Kotcioglu et al. Heat transfer analysis of a rectangular channel having tubular router in different winglet configurations with Taguchi method
CN113190999A (en) Electric heating coordination method and device for regulating flow of heat supply pipe network to improve wind power consumption
CN104101105A (en) Method and device used for compound control of crude oil heater temperature
CN105258396B (en) A kind of surface water source heat pump system energy consumption optimization method and device based on MATLAB
Aleksiejuk et al. Analog model of dynamics of a flat-plate solar collector
CN103939968B (en) Multi-heat source looped network heating system water circulating pump is arranged and pipeline parameter method of adjustment
Viswanathan et al. Dynamic model for a small-capacity ammonia–water absorption chiller
CN109670273A (en) A kind of coach air conditioner Performance Match method
CN106066937B (en) It is a kind of using heat-conduction oil heat as the forming machine hot blast temperature evaluation method of heat source
CN111625906B (en) Method for selecting fan type in power battery air cooling system by combining test and simulation
CN202033612U (en) Optimized guide system for cold-end equipment of thermal power generating unit
CN113050450B (en) Parallel variable frequency pump distribution system simulation module compiling method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20191104