发明名称 多热源环状管网可及性分析判定方法及其改善方法
摘要 本发明公开了一种多热源环状管网可及性分析判定方法及其改善方法,这种分析判定方法建立在管网结构参数和运行工况确定的情况下,通过定量分析各热源所能提供的循环压头能否满足热网在该工况下实现预定的流量分配所需的资用压头,来判定管网在这种工况下的可及性;并依据可及性分析判定结果,通过对环状管网的结构参数包括管段管径、循环水泵扬程等,按一定顺序进行调整,使其可及性能够同时满足各种运行工况的要求。本发明利用建立的多热源环状管网的可及性分析方法可以直观地判断出管网的可及性及影响可及性的主要因素;经过最后的管径调整和循环水泵调整后,管网在各种工况下的可及性则均能满足要求。
申请公布号 CN103605837A 申请公布日期 2014.02.26
申请号 CN201310537328.2 申请日期 2013.11.04
申请人 河南理工大学 发明人 徐文忠;吕鲲;张娅玲;张永胜;张丹
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 济南金迪知识产权代理有限公司 37219 代理人 段毅凡
主权项 1.一种多热源环状管网可及性分析判定方法,其特征在于,该分析方法是在计算机上实现的,详细步骤如下:第一步:确定需要进行可及性分析评价的工况,并将多热源环状管网分解为供水环网和回水环网,分别绘出该工况下供水环网和回水环网的网络计算图;选定设有定压装置的基本热源的供回水环网节点作为供水环网和回水环网的压力基准点,并分别对供水管网和回水管网的节点和管段进行相应编号,其中:供水环网的节点编号分别用G1,G2,...表示,回水环网的相应节点编号分别用H1,H2,...表示;供水环网的管段编号分别用1G,2G,...,回水环网的相应管段编号分别用1H,2H,...表示;编号时,压力基准点编号为最大编号;在此基础上,分别写出供水环网的连接矩阵A<sub>G</sub>,转置连接矩阵A′<sub>G</sub>、阻力数矩阵<img file="FDA0000407828640000015.GIF" wi="83" he="62" />节点额定流量矩阵<img file="FDA0000407828640000016.GIF" wi="80" he="61" />管网流通面积方阵F<sub>G</sub>,以及回水环网的连接矩阵A<sub>H</sub>、转置连接矩阵A′<sub>H</sub>、阻力数矩阵<img file="FDA0000407828640000017.GIF" wi="85" he="65" />节点额定流量矩阵<img file="FDA0000407828640000018.GIF" wi="82" he="57" />管网流通面积方阵F<sub>H</sub>;上述供水环网或回水环网的流量矩阵的各元素,流入为负值,流出为正值;上述供水环网的管网流通面积方阵F<sub>G</sub>或回水环网的管网流通面积方阵F<sub>H</sub>是由供水环网或回水环网各管段的流通面积构成的对称方阵,主对角线元素为各管段的流通面积<img file="FDA0000407828640000011.GIF" wi="219" he="104" />式中d<sub>i</sub>表示管段内径,非主对角线元素为零;第二步:建立多热源环状管网相对于基准点的节点压降求解模型;<maths num="0001"><![CDATA[<math><mrow><mi>Y</mi><mo>&CenterDot;</mo><mover><mi>p</mi><mo>&OverBar;</mo></mover><mo>=</mo><mover><mi>q</mi><mo>&OverBar;</mo></mover></mrow></math>]]></maths>式中,Y——环状管网节点压降方程的系数矩阵,<img file="FDA0000407828640000013.GIF" wi="349" he="77" />A——环状管网的连接矩阵;A′——环状管网连接矩阵A的转置矩阵;<img file="FDA0000407828640000019.GIF" wi="44" he="61" />——环状管网管段阻耗方阵,<img file="FDA0000407828640000014.GIF" wi="238" he="99" /><img file="FDA0000407828640000021.GIF" wi="46" he="68" />——环状管网各管段阻力数构成的对称方阵,主对角线上元素分别为各管段阻力数,非主对角线元素为0;<img file="FDA0000407828640000022.GIF" wi="67" he="99" />——环状管网各管段流量绝对值构成的对称方阵,主对角线上元素分别为各管段流量绝对值,非主对角线元素为0;<img file="FDA0000407828640000023.GIF" wi="49" he="65" />——环状管网的节点压降矩阵;<img file="FDA0000407828640000024.GIF" wi="405" he="84" /><img file="FDA0000407828640000025.GIF" wi="44" he="65" />——环状管网各节点额定流量矩阵,<img file="FDA0000407828640000026.GIF" wi="393" he="84" />第三步:求解供水环网的节点压降矩阵<img file="FDA0000407828640000027.GIF" wi="90" he="67" />利用如下的迭代算法获得节点压降矩阵<img file="FDA0000407828640000028.GIF" wi="86" he="67" />详细算法如下:第3.1步:任意设定管网的初始流速矩阵<img file="FDA0000407828640000029.GIF" wi="436" he="86" />矩阵各元素代表各管段初始流速,根据<img file="FDA00004078286400000210.GIF" wi="208" he="77" />求得管段的初始流量矩阵<img file="FDA00004078286400000211.GIF" wi="458" he="86" />及相应的管段初始流量构成的对称方阵<img file="FDA00004078286400000212.GIF" wi="112" he="96" />上述<img file="FDA00004078286400000213.GIF" wi="48" he="77" />代表管网流量矩阵,元素为各管段流量;上述F代表管网流通面积方阵,供水环网为F<sub>G</sub>,回水管网为F<sub>H</sub>;上述<img file="FDA00004078286400000214.GIF" wi="44" he="55" />代表管网流速矩阵,元素为各管段流速;第3.2步:根据公式<maths num="0002"><![CDATA[<math><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><mn>8</mn><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>j</mi></msub><mfrac><msub><mi>l</mi><mi>j</mi></msub><msub><mi>d</mi><mi>j</mi></msub></mfrac><mo>+</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&zeta;</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><msup><mi>&pi;</mi><mn>2</mn></msup><mo>&CenterDot;</mo><msup><msub><mi>d</mi><mi>j</mi></msub><mn>4</mn></msup><mo>&CenterDot;</mo><mi>g</mi></mrow></mfrac><mo>,</mo><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>=</mo><mn>0.11</mn><msup><mrow><mo>(</mo><mfrac><mi>K</mi><msub><mi>d</mi><mi>j</mi></msub></mfrac><mo>+</mo><mfrac><mn>68</mn><msub><mi>Re</mi><mi>j</mi></msub></mfrac><mo>)</mo></mrow><mn>0.25</mn></msup><mo>,</mo><msub><mi>Re</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>v</mi><mi>j</mi></msub><msub><mi>d</mi><mi>j</mi></msub></mrow><mi>&upsi;</mi></mfrac><mo>,</mo></mrow></math>]]></maths>确定矩阵<img file="FDA00004078286400000216.GIF" wi="82" he="78" />进而根据<img file="FDA00004078286400000217.GIF" wi="258" he="99" />确定管网阻耗矩阵<img file="FDA00004078286400000218.GIF" wi="84" he="77" />公式中s<sub>j</sub>、l<sub>j</sub>、d<sub>j</sub>、<img file="FDA00004078286400000219.GIF" wi="146" he="99" />λ<sub>i</sub>、Re<sub>j</sub>、v<sub>j</sub>分别代表j管段的管段阻力数、管段长度、管段内径、管段总局部阻力系数、管段摩擦阻力系数、管段雷诺数和管段流速;K代表管道的绝对粗糙系数;υ代表流体的运动粘度;第3.3步:根据<img file="FDA00004078286400000220.GIF" wi="372" he="77" />确定节点压降方程的系数矩阵Y<sub>0</sub>;第3.4步:利用公式<img file="FDA00004078286400000221.GIF" wi="225" he="67" />确定相对于压力基准点的节点压降矩阵<img file="FDA00004078286400000222.GIF" wi="112" he="65" />第3.5步:利用公式<img file="FDA00004078286400000223.GIF" wi="238" he="76" />确定管段压降矩阵<img file="FDA00004078286400000224.GIF" wi="84" he="76" />第3.6步:利用公式<img file="FDA0000407828640000031.GIF" wi="306" he="77" />确定管段流量矩阵<img file="FDA0000407828640000032.GIF" wi="84" he="77" />第3.7步:根据<img file="FDA0000407828640000033.GIF" wi="272" he="77" />确定管段流速矩阵<img file="FDA0000407828640000034.GIF" wi="97" he="66" />第3.8步:比较<img file="FDA0000407828640000035.GIF" wi="176" he="67" />检查其精度是否满足要求;预先设定精度要求如下:<img file="FDA0000407828640000036.GIF" wi="159" he="66" />两矩阵的对应元素相减所得的差值矩阵中,如果绝对值最大的元素的绝对值小于或等于0.01,即认为精度满足要求;如精度满足要求,则迭代结束;如不满足,则将<img file="FDA0000407828640000037.GIF" wi="48" he="67" />各元素值对应地替换<img file="FDA0000407828640000038.GIF" wi="49" he="67" />各元素值,重复第3.1-3.7步过程,直至满足预先设定的精度要求为止;第3.9步:输出管段流速矩阵<img file="FDA0000407828640000039.GIF" wi="75" he="66" />管段流量矩阵<img file="FDA00004078286400000310.GIF" wi="87" he="77" />管段压降矩阵<img file="FDA00004078286400000311.GIF" wi="46" he="70" />及节点压降矩阵<img file="FDA00004078286400000312.GIF" wi="81" he="65" />上述迭代过程结束获得的节点压降矩阵<img file="FDA00004078286400000313.GIF" wi="89" he="65" />即为供水环网的节点压降矩阵<img file="FDA00004078286400000314.GIF" wi="95" he="66" />第四步:求解回水环网的节点压降矩阵<img file="FDA00004078286400000315.GIF" wi="94" he="66" />采用与求解供水环网的节点压降矩阵<img file="FDA00004078286400000316.GIF" wi="71" he="67" />相同的迭代方法,可获得回水环网的节点压降矩阵<img file="FDA00004078286400000317.GIF" wi="84" he="66" />迭代过程中涉及的矩阵均采用回水管网的相应矩阵;第五步:确定多热源环状管网的最不利节点;根据多热源环状管网的供水环网节点压降矩阵<img file="FDA00004078286400000318.GIF" wi="66" he="67" />及回水环网节点压降矩阵<img file="FDA00004078286400000319.GIF" wi="104" he="76" />利用下式可判定出多热源环状管网的最不利节点J;p<sub>JG</sub>-p<sub>JH</sub>=max(p<sub>iG</sub>-p<sub>iH</sub>)i=1,2,…,n式中p<sub>JG</sub>——最不利节点J的供水环网节点压降,mH<sub>2</sub>O;p<sub>JH</sub>——最不利节点J的回水环网节点压降,mH<sub>2</sub>O;n——多热源环状管网的节点个数;第六步:确定多热源环状管网各热源所需提供的资用压头H<sub>rk</sub>;由于多热源环状管网的可及性分析多以基本热源的供、回水环网节点作为基准点,因此,基准点对应的用户一定为基本热源,则基本热源所需提供的资用压头按下式确定:H<sub>r1</sub>=Δp<sub>0</sub>+(p<sub>JG</sub>-p<sub>JH</sub>)+H<sub>0</sub>式中H<sub>r1</sub>——基本热源所需提供的资用压头,mH<sub>2</sub>O;Δp<sub>0</sub>——最不利节点J的最低允许资用压头,mH<sub>2</sub>O;H<sub>0</sub>——在额定流量下热源内部的阻力损失,mH<sub>2</sub>O;对于其它热源,其所需提供的资用压头可按下式确定:H<sub>rk</sub>=Δp<sub>0</sub>+(p<sub>JG</sub>-p<sub>JH</sub>)-(p<sub>iG</sub>-p<sub>iH</sub>)+H<sub>0</sub>k=2,…,M式中H<sub>rk</sub>——k热源需提供的资用压头,mH<sub>2</sub>O;p<sub>iG</sub>——k热源对应的供水环网节点的节点压降,mH<sub>2</sub>O;p<sub>iH</sub>——k热源对应的回水环网节点的节点压降,mH<sub>2</sub>O;M——多热源环状管网运行热源的个数;第七步:判定供热管网的可及性;设多热源环状管网某热源能够提供的循环压头为H′<sub>r1</sub>或H′<sub>rk</sub>,其中前者代表基本热源能够提供的循环压头,后者代表其它热源能够提供的循环压头,根据第六步确定的该热源需要提供的资用压头H<sub>r1</sub>或H<sub>rk</sub>,若H′<sub>r1</sub>≥H<sub>r1</sub>或H′<sub>rk</sub>≥H<sub>rk</sub>,则该热源满足管网在该工况下的可及性;反之,该热源不满足管网在该工况下的可及性;由此判断出所有热源在该工况下是否满足管网的可及性。
地址 454000 河南省焦作市高新区世纪大道2001号河南理工大学