常世杰 袁銘潤

(佳木斯大學信息電子技術(shù)學院，黑龍江 佳木斯 154007)

0 引言

伴隨著我國經(jīng)濟的飛速發(fā)展與科技水平的提高，人民生活水平得到顯著改善，針對于改善民生問題，熱網(wǎng)可以有效解決人民冬季冷暖問題，使人民不再飽受冬季寒冷之苦，在改善人民生活中的作用不可小覷。伴隨著熱網(wǎng)需求的不斷增多與可持續(xù)發(fā)展戰(zhàn)略的不斷推進，加快熱網(wǎng)基礎(chǔ)建設(shè)的同時，更要做到清潔供暖，節(jié)能減排。既要溫暖過冬，又要藍天白云，將黨和政府對人民的關(guān)懷“干凈”地送到千家萬戶，是供熱企業(yè)所共同奮斗的目標。

要將熱網(wǎng)做到“清潔、安全、高效、節(jié)能”四個方面，供暖方式的選擇，能源結(jié)構(gòu)的調(diào)整不可忽略，但是熱網(wǎng)分布更為重要。優(yōu)化熱網(wǎng)分布路徑，降低不必要的冗余熱網(wǎng)分布，縮短熱網(wǎng)管網(wǎng)輸送距離，可有效改善熱網(wǎng)因輸熱管距離較長導(dǎo)致熱量散失過多的問題，從本質(zhì)上提高熱網(wǎng)運行效率，同時降低熱網(wǎng)管網(wǎng)鋪設(shè)中能源設(shè)施消耗問題，減少熱網(wǎng)鋪設(shè)工作量，大幅度降低熱網(wǎng)投運成本，減少不必要的經(jīng)濟損失。

1 熱網(wǎng)換熱站分布模型建立

在實際生活中，熱網(wǎng)換熱站的分布需經(jīng)現(xiàn)場實地考察，測量以及地形分析，繪制出熱網(wǎng)換熱站擬分布圖形與方案，根據(jù)不同熱網(wǎng)換熱站節(jié)點位置與傳輸距離，篩選出較短的熱網(wǎng)路線，盡可能減小熱網(wǎng)在管網(wǎng)輸送載熱體中因傳輸距離導(dǎo)致的熱量散失。假設(shè)通過現(xiàn)場測量與數(shù)據(jù)分析，得到如圖1所示的供熱區(qū)域圖。其中熱源為火力發(fā)電廠，位于居民住宅區(qū)的邊緣，每個居民區(qū)都應(yīng)分布至少1個換熱站，居民區(qū)1的換熱站預(yù)設(shè)方案為2,3,5，居民區(qū)2的換熱站預(yù)設(shè)方案為6,4,7，居民區(qū)3的換熱站出于地形分布問題只能設(shè)置在8的位置。

熱源與換熱站間，換熱站與換熱站間的距離如表1所示。其中換熱站間距離單位為km，換熱站間無直接距離時為“—”。例如：熱源1到熱源1的距離為0，換熱站2到換熱站3間的距離為3 km，換熱站3到換熱站4間的距離為—，即沒有直接通路。

表1 熱網(wǎng)換熱站分布距離(預(yù)設(shè)方案)表

通過實際模型簡化，應(yīng)用圖論簡化實際熱網(wǎng)分布，將換熱站簡化成每一個帶標號的節(jié)點，其中節(jié)點1為熱源，節(jié)點2,3,4,5,6,7,8代表換熱站2,3,4,5,6,7,8，簡化后熱網(wǎng)換熱站分布模型圖如圖2所示。現(xiàn)需以熱源1為起始點，流經(jīng)居民區(qū)1中的換熱站(2/3/5)，再過居民區(qū)2中的換熱站(6/4/7)，最后到達居民區(qū)3，換熱站8為終止點，要通過應(yīng)用計算機輔助設(shè)計，規(guī)劃出熱網(wǎng)換熱站分布最短路線。

2 Matlab輔助設(shè)計

通過應(yīng)用計算機輔助設(shè)計軟件Matlab對熱網(wǎng)換熱站分布進行輔助設(shè)計，該問題可簡化成G(V(G),E(G))，其中V(G)代表圖論節(jié)點集;E(G)代表圖論邊長集，或權(quán)重集。通過應(yīng)用graph函數(shù)繪制G的圖像，graph(s,t,w)表示從節(jié)點s到極點t以w的權(quán)重創(chuàng)建邊，并在Matlab中繪制該圖像。Matlab程序如下：

>>s=[1 1 1 2 2 3 3 3 3 4 4 4 5 6 7];

t=[5 3 2 3 6 5 7 4 6 7 8 6 7 8 8];

w=[6 4 8 3 5 4 6 7 3 8 3 2 2 4 3];

G=graph(s,t,w);

plot(G,′EdgeLabel′,G.Edges.Weight,′linewidth′,2)

Matlab程序所繪制的熱網(wǎng)分布圖像如圖3所示。通過將該圖與圖2對比可得，兩者間差異不大，可全面的表示出各個換熱站節(jié)點與換熱站間距離，保證了下一步計算機輔助分析的不失真。

現(xiàn)要在節(jié)點1～節(jié)點8之間尋求一最短路徑，用來建立換熱站以及熱網(wǎng)的鋪設(shè)。通過應(yīng)用[P,d]=shortestpath(G,start,end[,′Method′,algorithm])函數(shù)，其中G為graph函數(shù)，start為起始熱源1，end為終止換熱站8，[,′Method′,algorithm]為最短路徑的求取方法，此處采用迪杰斯特拉Dijkstra算法，P為熱網(wǎng)分布最短路徑經(jīng)過的換熱站站點，d為最短路徑換熱站間的距離。Matlab程序如下：

>> clc,clear

s=[1 1 1 2 2 3 3 3 3 4 4 4 5 6 7];

t=[5 3 2 3 6 5 7 4 6 7 8 6 7 8 8];

w=[6 4 8 3 5 4 6 7 3 8 3 2 2 4 3];

G=graph(s,t,w);

[P,d]=shortestpath(G, 1, 8);

plot=plot(G, ′EdgeLabel′, G.Edges.Weight, ′linewidth′, 2);

highlight(plot, P, ′EdgeColor′, ′r′)

Matlab程序運行結(jié)果如圖4所示，求取最短路為節(jié)點1到節(jié)點3，再到節(jié)點6最后終止于節(jié)點8，此時熱網(wǎng)一次管網(wǎng)鋪設(shè)總距離為8 km，換熱站分別為3號換熱站，6號換熱站和8號換熱站，分別覆蓋三個居民區(qū)，在預(yù)設(shè)方案路徑中，該路徑所需鋪設(shè)管網(wǎng)最短，熱量在傳輸中的損失最少，大幅提升能源利用率，同時降低供熱工程所需材料設(shè)備，滿足設(shè)計要求。

3 設(shè)計校驗分析

在城市熱網(wǎng)鋪設(shè)時，往往會涉及到城鎮(zhèn)與城鎮(zhèn)之間，換熱站與熱源之間，換熱站與換熱站之間的路徑分布問題，從多條支路方案中篩選出距離相對較近的路線，是整個熱網(wǎng)換熱站分布的核心。通過應(yīng)用迪杰斯特拉算法確定熱網(wǎng)最短分布路線，迪杰斯特拉算法是目前在生活中針對分布規(guī)劃，最優(yōu)路徑的選取等問題應(yīng)用較多的一種算法，例如：鐵路物資調(diào)配問題，城市軌道交通路線鋪設(shè)問題，公交車站分布問題等，該算法的核心思想是以起始點為中心向外廣度優(yōu)先搜索至終點，從而尋找到從指定的起始點到指定的終止點間的最短路徑。應(yīng)用該算法將換熱站與熱源作為節(jié)點，預(yù)選方案路線距離作為參考，實現(xiàn)最短路徑選取。

通過上述Matlab程序設(shè)計及結(jié)果分析，在該熱網(wǎng)分布預(yù)設(shè)方案中，在諸多組合方案中選出最短路徑方案，該方案為熱源1到終點換熱站8的最短路徑，同時針對于路徑上所經(jīng)過的換熱站3和6，該路徑仍是熱源1到換熱站3和6的最短路徑，即該路徑上所經(jīng)過的換熱站到熱源1的距離也是最短的。例：從熱源1到換熱站3的方案有：1—3(4 km)，1—2—3(11 km)，1—5—3(10 km)，從熱源1到換熱站6之間的最短路徑如圖5所示。相比較而言仍是該路徑上的方案距離最短。同時，應(yīng)用計算機輔助設(shè)計熱網(wǎng)分布，工作效率有較大提升，在面臨較多數(shù)據(jù)與多種預(yù)設(shè)方案時，可通過計算機代替人工進行繁雜計算，從而挑選出較優(yōu)的預(yù)設(shè)方案，實際熱網(wǎng)鋪設(shè)工程中仍需考慮較多外界因素。

4 結(jié)語

通過應(yīng)用Matlab建立熱網(wǎng)分布優(yōu)化數(shù)學模型，借助計算機強大的計算能力，應(yīng)用迪杰斯特拉Dijkstra算法在諸多預(yù)設(shè)熱網(wǎng)方案中選出最優(yōu)路徑，同時可推廣至較多預(yù)設(shè)方案的分析縮短熱網(wǎng)鋪設(shè)距離，可以有效降低熱網(wǎng)架設(shè)中的成本，降低熱網(wǎng)維護工作人員的任務(wù)量，同時減少熱網(wǎng)傳輸過程中的熱量損失，從根本上解決熱量散失與傳輸距離的矛盾問題，提升能源利用率，可為供熱帶來較高的經(jīng)濟效益與社會效益，具有較高的實用性。