由世俊，米雷洋，王雅然，張?歡，鄭雪晶，鄭萬冬

?

集中供熱管網(wǎng)的非穩(wěn)態(tài)水力建模與動態(tài)響應(yīng)分析

由世俊，米雷洋，王雅然，張?歡，鄭雪晶，鄭萬冬

(天津大學(xué)環(huán)境科學(xué)與工程學(xué)院，天津 300350)

建立了單個供熱管段的非穩(wěn)態(tài)水力模型，導(dǎo)出了該非穩(wěn)態(tài)水力模型的解析解，分析了管道長度、管徑、進(jìn)出口壓差變化時管道流體的動態(tài)響應(yīng)特性．當(dāng)單個管道的進(jìn)出口壓差變化相同，管道的長度越長、管徑越大時，管道的流量響應(yīng)時間越慢；對于相同規(guī)格的管道，當(dāng)管道的進(jìn)出口壓差變化越大時，管道流體的動態(tài)響應(yīng)速度越快．在單個管段的非穩(wěn)態(tài)水力模型的基礎(chǔ)上，結(jié)合圖論的方法建立了集中供熱管網(wǎng)的非穩(wěn)態(tài)水力模型，并構(gòu)建了該模型的數(shù)值求解方法．該模型考慮管網(wǎng)運(yùn)行調(diào)節(jié)時閥門開度或水泵轉(zhuǎn)速改變等影響動態(tài)響應(yīng)時間的因素，提出的非穩(wěn)態(tài)水力模型和求解方法為熱網(wǎng)進(jìn)行頻繁、大范圍水力工況調(diào)節(jié)過程中各熱力站流量的動態(tài)響應(yīng)分析提供了支持．對一個實際熱網(wǎng)的動態(tài)水力進(jìn)行了分析，結(jié)果表明：當(dāng)熱網(wǎng)中的閥門開度不變、熱源泵轉(zhuǎn)速發(fā)生變化時，距離熱源越遠(yuǎn)的熱力站的流量動態(tài)響應(yīng)時間越長；同樣地，當(dāng)熱網(wǎng)中熱源泵的轉(zhuǎn)速不變、改變某一熱力站的閥門開度時，距離熱源越遠(yuǎn)的熱力站的流量動態(tài)響應(yīng)時間越長．

非穩(wěn)態(tài)水力模型；圖論；集中供熱管網(wǎng)；動態(tài)響應(yīng)時間；運(yùn)行調(diào)節(jié)

集中供熱管網(wǎng)的動態(tài)水力和熱力特性對系統(tǒng)的運(yùn)行調(diào)節(jié)至關(guān)重要．熱網(wǎng)的動態(tài)熱力特性主要體現(xiàn)在供熱管道中熱媒溫度的傳輸延遲[1]，動態(tài)水力特性主要表現(xiàn)為閥門和水泵調(diào)節(jié)后，各個支路流量的動態(tài)響應(yīng)．熱網(wǎng)的穩(wěn)態(tài)水力特性可以由節(jié)點流量平衡和回路壓力平衡結(jié)合圖論的方法來描述[2]，并且基于穩(wěn)態(tài)水力模型數(shù)值計算的仿真和優(yōu)化調(diào)度技術(shù)已經(jīng)在集中供熱系統(tǒng)的運(yùn)行分析和調(diào)節(jié)中得到廣泛而有效的應(yīng)用[2]．

熱網(wǎng)在運(yùn)行過程中，熱源需要通過調(diào)節(jié)熱源泵的轉(zhuǎn)速實現(xiàn)熱源流量的調(diào)度，各個熱力站則利用本地閉環(huán)控制系統(tǒng)通過改變一次側(cè)閥門的開度實現(xiàn)流量的自動調(diào)節(jié)[3-4]．而閥門和水泵的動作會使熱網(wǎng)中各個管段的流量經(jīng)歷一個動態(tài)過程，一方面，熱網(wǎng)流量的動態(tài)響應(yīng)時間將對熱網(wǎng)的水力調(diào)節(jié)產(chǎn)生影響；另一方面，隨著通信技術(shù)和自動化技術(shù)的發(fā)展，一些先進(jìn)的通信技術(shù)和控制技術(shù)在不斷地應(yīng)用到熱網(wǎng)當(dāng)中，熱網(wǎng)智能化趨勢不可避免，為了在滿足用戶熱舒適的前提下，最大限度地節(jié)省熱網(wǎng)的輸配能耗，依賴于精細(xì)化的供熱調(diào)節(jié)，而熱網(wǎng)流量的動態(tài)響應(yīng)時間也是精細(xì)化供熱調(diào)節(jié)時間間隔的依據(jù)．因此，建立管道和熱網(wǎng)的非穩(wěn)態(tài)水力模型，并分析管道和熱網(wǎng)的動態(tài)響應(yīng)特性具有重要的意義．

目前關(guān)于供熱管網(wǎng)的非穩(wěn)態(tài)水力特性的研究主要集中在供熱管網(wǎng)的水錘問題及管網(wǎng)的慢變流問題．蔡啟林等[5]建立了供熱管網(wǎng)的水錘計算的數(shù)學(xué)模型，通過實例驗證了該模型的正確性；還研究了供熱系統(tǒng)水錘防護(hù)的措施，為供熱管網(wǎng)的設(shè)計和運(yùn)行提供了參考．曹慧哲等[6]基于圖論建立了環(huán)狀管網(wǎng)慢變流的水力計算模型，給出了管網(wǎng)慢變流模型的數(shù)值模擬方法；并基于慢變流理論提出了一種管網(wǎng)泄漏診斷與定位的模型，對管網(wǎng)阻抗的辨識和管網(wǎng)的故障診斷具有重要的參考價值．

本文從最基本的質(zhì)量守恒和能量守恒原理出發(fā)，建立了單個管段和集中供熱管網(wǎng)的動態(tài)水力模型，給出了一種新的集中供熱管網(wǎng)慢變流建模和數(shù)值求解方法，為集中供熱系統(tǒng)的動態(tài)水力特性分析提供補(bǔ)充．所建立的集中熱網(wǎng)的動態(tài)水力模型既適用于規(guī)模不同的枝狀管網(wǎng)，也適用于環(huán)狀管網(wǎng)的水力特性分析．同時，本文推導(dǎo)出了管段流量的瞬態(tài)響應(yīng)時間的計算公式，分析了不同管徑、管長的管道的瞬態(tài)響應(yīng)時間．最后，根據(jù)所建立的集中供熱管網(wǎng)的動態(tài)水力模型，通過一個實例分析了熱網(wǎng)運(yùn)行過程中熱源泵的轉(zhuǎn)速變化和熱力站閥門開度變化時熱網(wǎng)的非穩(wěn)態(tài)水力響應(yīng)特性．

1?單個管段的非穩(wěn)態(tài)水力特性

1.1 單個管段的非穩(wěn)態(tài)水力建模

圖1是單個管段流體非穩(wěn)態(tài)流動的動量方程推導(dǎo)示意．坐標(biāo)軸沿管道中心線方向，控制體的標(biāo)高為，根據(jù)牛頓第二定律得到以下方程：

式中：為流體沿坐標(biāo)軸的壓力分布；為斷面面積，m2；δ為控制體的厚度，m；0為管道壁面對控制體的切應(yīng)力，Pa；為控制體的密度，kg/m3；為重力加速度，取9.8m/s2；為管道中心線與水平面的夾角；為控制體的加速度，m/s2．

圖1 ?單個管段流體非穩(wěn)態(tài)流動動量方程推導(dǎo)示意

消去δ，式(1)可以化簡為[7]

在上述兩式中，切應(yīng)力0可以按照穩(wěn)態(tài)工況下的達(dá)西公式計算：

式中：為沿程阻力損失系數(shù)；為斷面平均流速，m/s；為管長，m；為管道直徑，m．

在穩(wěn)態(tài)工況下，對于長度為、管徑為的水平管道，壁面切應(yīng)力引起的阻力與按照達(dá)西公式計算得到的阻力相等且滿足下式：

由此可以解得切應(yīng)力0為

加速度可以分解成如下形式：

式中：v為速度對位置的偏導(dǎo)數(shù)；v為速度對時間的偏導(dǎo)數(shù)．

在上述分析當(dāng)中，忽略壓力波在管道中的傳播和衰減過程，而且管道中的壓力波的傳播速度要遠(yuǎn)快于流量的動態(tài)變化速度；又由于供熱管道中的熱水近似為不可壓縮流體，因此可將v近似為0．對于水平管段，＝0°；由于供熱管段的坡度[8]一般為0.002～0.005，故對于有坡度的供熱管段，sin≈0．根據(jù)以上條件，將式(5)和式(6)代入到式(2)中得

將上式從＝0到＝沿管長進(jìn)行積分[9]，得

式中Δ＝0－L為管段的進(jìn)出口壓差，0為進(jìn)口壓力，L為出口壓力．

式(8)的物理意義是：管段進(jìn)出口壓差Δ是使流體加速的動力，而管壁的黏性切應(yīng)力則對流體起減速作用．當(dāng)流體的運(yùn)動狀態(tài)達(dá)到穩(wěn)態(tài)時，v＝0，方程(8)退化為描述管道穩(wěn)態(tài)流動沿程阻力的達(dá)西公式．將式(8)中的流速改寫成流量的形式，得到單個管段的非穩(wěn)態(tài)水力模型表達(dá)式為

式(9)表示管段中流體流動的非穩(wěn)態(tài)特性與管段的管徑、長度、管段的阻力系數(shù)等有關(guān)．

1.2 單個管段非穩(wěn)態(tài)水力建模的解析解

方程(9)描述了管道流體在非穩(wěn)態(tài)流動過程中壓力和流量變化之間的關(guān)系，它是1階非線性常微分方程．設(shè)某管段在＝0時刻之前流體處于穩(wěn)態(tài)，管段流體的流量為0，管段流體的進(jìn)出口壓差為Δ0．當(dāng)＝0時刻時管段的進(jìn)出口壓差產(chǎn)生了一個階躍，由Δ0變化為Δ1，此時刻之后管段的流量將經(jīng)過一個瞬態(tài)的變化過程．

對式(10)兩端進(jìn)行積分，得

式(11)中，的表達(dá)式為

1.3 單個管段的非穩(wěn)態(tài)水力特性分析

基于式(11)的關(guān)系對管道的非穩(wěn)態(tài)水力特性進(jìn)行分析，同時根據(jù)式(11)對管段內(nèi)流體流動的響應(yīng)速度進(jìn)行研究，并探索管段內(nèi)流體從靜止過渡到另一個穩(wěn)態(tài)時，管段流體流量的響應(yīng)時間與管徑、管長等管段物理參數(shù)之間的變化規(guī)律．

若管段初始狀態(tài)的流量0＝0，則管段初始狀態(tài)的進(jìn)出口壓差Δ0＝0．在＝0時刻，管段兩端的壓差從0階躍為Δ1，由式(11)導(dǎo)出管段內(nèi)流量隨時間的變化規(guī)律如下：

達(dá)到穩(wěn)態(tài)時的流量為

定義管段流量的瞬態(tài)響應(yīng)時間r：管段初始狀態(tài)的流量為0，當(dāng)管段進(jìn)出口壓差從0階躍為Δ1時，管段流量達(dá)到新的穩(wěn)態(tài)流量(∞)的95%所需的時間r為管段流量的瞬態(tài)響應(yīng)時間．根據(jù)式(14)和式(15)以及管段流量的瞬態(tài)響應(yīng)時間的定義，可得

求解式(16)，得出管段流量的瞬態(tài)響應(yīng)時間

由式(17)可以看出：管段的管徑越大、長度越長，則管段流量的瞬態(tài)響應(yīng)時間越長，流體從靜止到新的穩(wěn)態(tài)的過渡過程越緩慢；管段進(jìn)出口壓差越大，管段內(nèi)流體達(dá)到穩(wěn)態(tài)的時間越短．

圖2所示為管徑DN800、管長100m，管徑DN800、管長200m和管徑DN600、管長100m的3種不同規(guī)格的供熱管段的瞬態(tài)響應(yīng)時間．初始時刻，3根管道中的熱水均處于靜止?fàn)顟B(tài)，當(dāng)3根管道的進(jìn)出口壓力差在＝0時刻階躍為10kPa并保持恒定，3根管段的流量均經(jīng)過一個暫態(tài)過程后達(dá)到穩(wěn)定．由圖2可以看出，管徑為DN800、管長為100m的管段的瞬態(tài)響應(yīng)時間為r＝55.6s，管徑為DN800、管長為200m的管段的瞬態(tài)響應(yīng)時間為r＝78.6s，管徑為DN600、管長為100m的管段的瞬態(tài)響應(yīng)時間為r＝46.4s．仿真結(jié)果表明：管段的長度越長，管徑越大，管段的流量響應(yīng)時間越慢．

圖2? 不同管徑、管長的管道流體動態(tài)響應(yīng)過程

圖3所示是相同管段在不同進(jìn)出口壓差變化情況下所引起的管段流體動態(tài)變化過程中的流量響應(yīng)．可以看出，管徑為DN800、長度為100m的供熱管段，進(jìn)出口壓差從0Pa變?yōu)?0kPa引起的瞬態(tài)流量響應(yīng)時間為r＝55.6s，進(jìn)出口壓差從初態(tài)的0Pa變?yōu)?0kPa引起的瞬態(tài)流量響應(yīng)時間為r＝39.3s.表明進(jìn)出口壓差變化越大，動態(tài)響應(yīng)速度越快．

圖3 不同進(jìn)出口壓差階躍變化下的管道流體動態(tài)響應(yīng)過程

圖4所示是不同公稱直徑、不同長度管段中的熱水由靜止?fàn)顟B(tài)開始變化，在進(jìn)出口壓差從0Pa瞬間增加到10kPa時，各個管段中流體的瞬態(tài)響應(yīng)時間．可以看出，不同規(guī)格管段流體的瞬態(tài)響應(yīng)時間在幾分鐘到十幾分鐘之間變化．

圖4 ?不同管徑、不同管長的管道流體響應(yīng)時間

2 ?集中供熱管網(wǎng)的非穩(wěn)態(tài)水力建模與分析

本文第1節(jié)給出了單個管段的非穩(wěn)態(tài)水力特性描述，并得出了單個供熱管段由一個穩(wěn)態(tài)向另一個穩(wěn)態(tài)過渡時的瞬態(tài)過程的解析解，即式(11)．但是當(dāng)管網(wǎng)中有閥門的開度變化或者水泵轉(zhuǎn)速的變化時，式(11)就不適用于分析管網(wǎng)的非穩(wěn)態(tài)水力響應(yīng)了．在考慮上述兩個因素對管網(wǎng)的非穩(wěn)態(tài)水力響應(yīng)的影響時，必須聯(lián)立各個管段的非穩(wěn)態(tài)水力方程，即式(9).本節(jié)將建立集中供熱管網(wǎng)的非穩(wěn)態(tài)水力模型，并給出數(shù)值求解方法，進(jìn)而分析管網(wǎng)的非穩(wěn)態(tài)水力響應(yīng)．

2.1?集中供熱管網(wǎng)的非穩(wěn)態(tài)水力模型的建立

與供熱管網(wǎng)的穩(wěn)態(tài)水力模型類似，集中供熱管網(wǎng)的非穩(wěn)態(tài)水力模型也可以采用圖論的方法描述．由于熱水在熱網(wǎng)中的非穩(wěn)態(tài)水力過程中是不可壓縮的，故節(jié)點流量平衡方程和回路壓力平衡方程[10]逐時地成立，滿足以下兩個方程式：

式中：k為管網(wǎng)的基本關(guān)聯(lián)矩陣；f為管網(wǎng)的獨立回路矩陣[10]；()為各個分支在時刻的流量向量；Δ()為各個分支在時刻的壓差向量．

根據(jù)熱網(wǎng)的節(jié)點流量平衡方程和回路壓力平衡方程以及水泵揚(yáng)程的特性方程p,i(n，Q)，管道的非穩(wěn)態(tài)水力特性方程(9)可以寫成如下形式：

式中：Δp為各個管段的進(jìn)出口壓差；S為與局部阻力損失和沿程阻力損失相關(guān)的管路阻抗，即

f,i的表達(dá)式為

閥門開度變化與閥門阻抗(x)之間的關(guān)系可以通過下式表示：

式中：(x)為等百分比閥門的阻抗，x為等百分比閥門的開度；R為閥門的可調(diào)比；v,i為閥門的流通能力.

若集中供熱管網(wǎng)中有個管段，則個管段的動態(tài)水力特性方程可寫成向量形式，即

聯(lián)立式(19)和式(24)，消掉Δ，得

整個管網(wǎng)各個分支的流量向量也可以由余枝流量表示為如下形式：

其中

式中r是與余枝數(shù)階數(shù)相等的單位矩陣．

將式(32)與式(33)代入式(29)中，得到整個管網(wǎng)的非穩(wěn)態(tài)水力模型

上述各式中未注釋項代表的意義見上節(jié)．

式(34)是非線性常微分方程，該式中方程的數(shù)目與管網(wǎng)中的余枝數(shù)相等，即在節(jié)點數(shù)為、分支數(shù)為的管網(wǎng)中，余枝數(shù)為＋1．通過求解式(34)可以得到余枝流量隨時間的變化，再根據(jù)式(31)可以得到樹枝流量隨時間的變化．

2.2?管網(wǎng)非穩(wěn)態(tài)水力模型求解方法

由于集中供熱管網(wǎng)的非穩(wěn)態(tài)水力方程是非線性常微分方程組，直接獲得其解析解很困難．為方便進(jìn)行供熱管網(wǎng)非穩(wěn)態(tài)水力響應(yīng)的分析，本文利用4階Runge-Kutta法構(gòu)造式(34)的數(shù)值求解格式[11]．設(shè)供熱管網(wǎng)的初始狀態(tài)下，余枝流量為r,0．設(shè)時間步長為Δ，第步的余枝流量為r,k，則第＋1步的余枝流量r,k＋1可由下式給出：

并且，

3?集中供熱管網(wǎng)的非穩(wěn)態(tài)水力分析

3.1?水泵轉(zhuǎn)速變化對熱網(wǎng)流量水力響應(yīng)時間的影響

圖5所示是某城鎮(zhèn)的集中供熱管網(wǎng)，該熱網(wǎng)有一個環(huán)、兩個熱源和28個熱力站，是典型的多熱源環(huán)狀熱網(wǎng)，總供熱負(fù)荷為190MW．熱源1是熱電廠，承擔(dān)基本負(fù)荷，其供熱能力為150MW．熱源2是燃?xì)庹{(diào)峰鍋爐房，承擔(dān)峰值負(fù)荷，其供熱能力為55MW．在低負(fù)荷時，熱源1承擔(dān)全部熱力站的負(fù)荷，熱源1循環(huán)泵的調(diào)節(jié)會對各個熱力站的流量產(chǎn)生影響，在分析熱網(wǎng)各部分流量隨著熱源泵轉(zhuǎn)速變化的響應(yīng)特性時，考慮如下典型的情形：當(dāng)熱源1的水泵電機(jī)頻率在1min內(nèi)從50Hz線性地減為25Hz，而各個熱力站的閥門開度保持不變時，分析各熱力站的流量隨時間的變化情況．本節(jié)將基于熱網(wǎng)非穩(wěn)態(tài)水力建模及求解方法，分析上述情形下各熱力站的流量隨時間的變化情況．

在MATLAB中編制上節(jié)所提出的基于4階Runge-Kutta法的熱網(wǎng)非穩(wěn)態(tài)水力特性計算程序，并導(dǎo)入圖5所示熱網(wǎng)對應(yīng)式(34)形式的非穩(wěn)態(tài)水力模型，可以實現(xiàn)該熱網(wǎng)動態(tài)水力工況的仿真計算．

圖5 ?某城鎮(zhèn)集中供熱管網(wǎng)

在進(jìn)行熱網(wǎng)的非穩(wěn)態(tài)水力模型求解時，為保證Runge-Kutta算法計算的數(shù)值穩(wěn)定性，時間步長Δ不能選得過大，本例選取Δ＝0.1s，數(shù)值計算的總時長為30min．熱源1處的循環(huán)水泵電機(jī)頻率在＝0～5min內(nèi)維持在50Hz不變，從＝5min開始，到＝6min時，熱源1處循環(huán)水泵電機(jī)的頻率從50Hz線性地減為25Hz，并從＝6min開始、到＝30min時，維持熱源1處循環(huán)水泵電機(jī)的頻率在25Hz不變，在整個過程中，各個熱力站的閥門開度均維持在100%不變．部分熱力站的流量隨時間變化的計算結(jié)果見圖6．

圖6 ?熱源泵調(diào)節(jié)時部分熱力站的流量響應(yīng)

從圖6可以看出，當(dāng)熱源1的循環(huán)泵轉(zhuǎn)速降低以后，熱力站的流量經(jīng)過一個瞬態(tài)過程，達(dá)到新的穩(wěn)態(tài)值，且新的穩(wěn)態(tài)值較熱源1循環(huán)泵轉(zhuǎn)速變化之前有所減少．從圖6還可以得出，按管路距離計算，離熱源1近的熱力站流量的響應(yīng)速度都比較快，例如1～4號熱力站，其流量都在1min左右達(dá)到新的穩(wěn)態(tài)值，見圖6(a)～(d)．但隨著與熱源1管路距離的增大，熱力站的流量響應(yīng)速度逐漸變慢，例如10～13號熱力站，其流量的響應(yīng)時間在3～4min，見圖6(e)～(h)．對于距離熱源1較遠(yuǎn)的熱力站，如26～28號熱力站，其流量的響應(yīng)時間在5～8min，見圖6(i)～(k)．

3.2?閥門開度變化對熱網(wǎng)流量水力響應(yīng)時間的影響

對于熱力站閥門開度變化下的各熱力站流量響應(yīng)特性，考慮令熱源1的水泵電機(jī)頻率保持50Hz不變，熱力站1的閥門開度在1min內(nèi)由100%線性地變化為50%．?dāng)?shù)值計算的總時長也為30min．熱源1的循環(huán)水泵電機(jī)頻率在＝0～30min內(nèi)維持在50Hz不變．熱力站1的閥門開度在＝0～5min內(nèi)保持在100%不變，從＝5min開始、到＝6min時，熱力站1的閥門開度從100%線性地變?yōu)?0%，并從＝6min時開始、到＝30min時，維持熱力站1的閥門開度在50%不變．在整個過程中，其余熱力站的閥門開度均維持在100%不變．部分熱力站的流量隨時間變化的計算結(jié)果見圖7．

從圖7可以看出，當(dāng)熱力站1的閥門開度減小以后，其余熱力站的流量都經(jīng)過一個瞬態(tài)過程，達(dá)到新的穩(wěn)態(tài)值，且新的穩(wěn)態(tài)值較熱力站1閥門動作之前都有所增加．同樣地，從圖7中也可以得出，按管路距離計算，離熱力站1近的熱力站流量響應(yīng)速度都比較快，例如2～5號熱力站，其流量都在1min內(nèi)達(dá)到新的穩(wěn)態(tài)值，見圖7(a)～(d)．但隨著與熱源1管路距離的增大，熱力站的流量響應(yīng)速度逐漸變慢，例如9～13號熱力站，其流量的響應(yīng)時間在2min左右，見圖7(e)～(i)．對于距離熱源1較遠(yuǎn)的熱力站，如26～28號熱力站，其流量的響應(yīng)時間在4～6min左右，見圖7(j)～圖7(l)．

圖7 ?閥門調(diào)節(jié)時部分熱力站的流量響應(yīng)

圖7給出的結(jié)果表明：若熱網(wǎng)中的閥門不動，熱源泵轉(zhuǎn)速發(fā)生變化以后，按管路距離計算，距離該熱源越遠(yuǎn)的熱力站，其流量的動態(tài)響應(yīng)時間越長，反之則越短．同樣地，若熱網(wǎng)中的水泵不調(diào)節(jié)，而某一熱力站的閥門開度發(fā)生變化后，按管路距離計算，距離閥門動作位置越遠(yuǎn)的熱力站，其流量的動態(tài)響應(yīng)時間越長，反之則越短．當(dāng)熱源循環(huán)泵的轉(zhuǎn)速進(jìn)行大范圍調(diào)整時，各熱力站的流量響應(yīng)時間比熱力站進(jìn)行閥門調(diào)節(jié)所引起的流量響應(yīng)時間長．

這些規(guī)律可以為大型熱網(wǎng)的水力工況調(diào)節(jié)提供指導(dǎo)，并且，本文給出的建模與數(shù)值求解方法也為熱網(wǎng)動態(tài)水力特性的量化分析和精細(xì)化調(diào)節(jié)提供支持.

4?結(jié)?論

(1) 本文建立的單個管段的非穩(wěn)態(tài)水力模型可以用來計算分析單個不同規(guī)格管段的瞬態(tài)響應(yīng)時間．利用該模型還可以得到單個管段的瞬態(tài)水力特性：管道的管徑越大、長度越長，則管道兩端壓差發(fā)生變化后，流量的動態(tài)響應(yīng)時間越長，流體從靜止到新的穩(wěn)態(tài)的過渡過程越緩慢；而若管道進(jìn)出口壓差變化越大，則管道流量的動態(tài)響應(yīng)時間越短，流量達(dá)到新穩(wěn)態(tài)的過程越快．

(2) 本文建立的集中供熱管網(wǎng)的非穩(wěn)態(tài)水力模型可以用來計算分析不同規(guī)模熱網(wǎng)在水泵轉(zhuǎn)速變化和閥門開度變化時各個熱力站的動態(tài)響應(yīng)時間及管網(wǎng)的動態(tài)水力特性．

(3) 算例中熱源循環(huán)泵轉(zhuǎn)速的變化和熱力站閥門開度的變化過程所引起的動態(tài)水力響應(yīng)的時間尺度在10min以內(nèi)．對于更大規(guī)模的城市級的集中供熱管網(wǎng)，動態(tài)水力響應(yīng)時間會更長．為了保證供熱管網(wǎng)的水力穩(wěn)定性同時最大限度節(jié)省供熱管網(wǎng)的輸配能耗，當(dāng)進(jìn)行管網(wǎng)的運(yùn)行調(diào)節(jié)時，可將本文所建立的管網(wǎng)非穩(wěn)態(tài)水力模型和分析方法應(yīng)用到熱網(wǎng)中，根據(jù)計算出的動態(tài)響應(yīng)時間作為制定水力工況調(diào)節(jié)時間間隔的依據(jù)．

［1］ Wang Yaran，You Shijun，Zhang Huan，et al. Thermal transient prediction of district heating pipeline：Optimal selection of the time and spatial steps for fast and accurate calculation[J]. Applied Energy，2017，206：900-910.

［2］ Wang Yaran，You Shijun，Zhang Huan，et al. Hydraulic performance optimization of meshed district heating network with multiple heat sources[J]. Energy，2017，126：603-621.

［3］ Wang Yaran，You Shijun，Zheng Wandong，et al. State space model and robust control of plate heat exchanger for dynamic performance improvement[J]. Applied Thermal Engineering，2018，158：1588-1604.

［4］ Wang Yaran，You Shijun，Zhang Huan，et al. Operation stability analysis of district heating substation from the control perspective[J]. Energy and Buildings，2017，154：373-390.

［5］ 蔡啟林，李?銳. 熱水供熱網(wǎng)水錘計算數(shù)學(xué)模型[J]. 區(qū)域供熱，1995(5)：28-32.

Cai Qilin，Li Rui. Mathematical model of water hammer in hot water heating network[J]. District Heating，1995(5)：28-32(in Chinese).

［6］ 曹慧哲，賀志宏，何鐘怡. 基于圖論的環(huán)狀管網(wǎng)慢變流的計算研究[J]. 哈爾濱工業(yè)大學(xué)學(xué)報，2007，39(10)：1559-1563.

Cao Huizhe，He Zhihong，He Zhongyi. Analysis and calculation of the slow transient flow in the multi-loops network based on graph theory[J]. Journal of Harbin Institute of Technology，2007，39(10)：1559-1563(in Chinese).

［7］ Wylie E B，Streeter V L. Fluid Transients[M]. Osborne：McGraw-Hill，1978.

［8］ 賀?平，孫?剛，王?飛，等. 供熱工程[M]. 4版. 北京：中國建筑工業(yè)出版社，2009.

He Ping，Sun Gang，Wang Fei，et al. Heating Engineering[M]. 4th ed. Beijing：China Architecture and Building Press，2009(in Chinese).

［9］ American Society of Heating，Refrigerating and Air Conditioning Engineers. ASHRAE Handbook of Fundamentals[M]. American：ASHRAE Handbook Committee，2013.

［10］ 付祥釗，肖益民. 流體輸配管網(wǎng)[M]. 3版. 北京：中國建筑工業(yè)出版社，2010.

Fu Xiangzhao，Xiao Yimin. Fluid Supply Pipework[M]. 3rd ed. Beijing：China Architecture and Building Press，2010(in Chinese).

［11］ 李榮華，劉?播. 微分方程數(shù)值解法[M]. 3版. 北京：高等教育出版社，2009.

Li Ronghua，Liu Bo. Numerical Methods for Solution of Differential Equations[M]. 3rd ed. Beijing：Higher Education Press，2009(in Chinese).

Unsteady Hydraulic Modeling and Dynamic Response Analysis of a District Heating Network

You Shijun，Mi Leiyang，Wang Yaran，Zhang Huan，Zheng Xuejing，Zheng Wandong

(School of Environmental Science and Engineering，Tianjin University，Tianjin 300350，China)

In this paper，the unsteady hydraulic model of a single pipeline is established，and the analytical solution of the unsteady hydraulic model is derived．The dynamic response characteristics of the pipeline fluid are analyzed when the pipeline length，pipeline diameter and differential pressure of the inlet and outlet vary．When the differential pressure between the inlet and outlet of a single pipeline changes in the same way，the longer the length of the pipeline and the larger the diameter of the pipeline，the slower the flow response time of the pipeline；moreover，for the same size of the pipeline，the greater the differential pressure between the inlet and outlet of the pipeline，the higher the dynamic response speed of the fluid in the pipeline．Based on the dynamic model of a single pipe and graph theory，the dynamic hydraulic model of the whole district heating network is established，and the numerical solution method of the model is constructed．The model considers the factors affecting the dynamic response time，such as the opening of the valve or the rotational speed change of the pump，when the pipe network operation is adjusted．The proposed unsteady hydraulic model and solution method provide support for dynamic response analysis of the flow of each thermal station in the process of frequent and large-scale hydraulic regulation of the heat network．Through the dynamic hydraulic analysis of a practical heat network，it is determined that when the opening of the valve in the district heating network is constant and the speed of the heat source pump changes，the dynamic response time of the flow rate in the heating substation increases，if the distance of the heating substation away from the heat source increases．Similarly，when the speed of the heat source pump is constant and the valve opening of a heating substation is changed，the dynamic response time of the flow rate in the heating substation increases，if the distance of the heating substation away from the heat source increases．

unsteady hydraulic model；graph theory；district heating network；dynamic response time；operating regulation

TU995.3

A

0493-2137(2019)08-0849-08

10.11784/tdxbz201808091

2018-08-31；

2019-01-09.

由世?。?955—??），男，博士，教授，yousj@tju.edu.cn.

王雅然，wangyaran19900101@126.com.

(責(zé)任編輯：田?軍)