李樹平

（同濟大學 環(huán)境科學與工程學院，上海 200092）

給水管網建模中，管網簡化是不可避免的。通過簡化，可以減少模型開發(fā)和維護工作量，提高運行效率。相同直徑、材料和埋設年代的串聯(lián)管道合并是常用簡化方法之一。當中間節(jié)點有需水量的多段串聯(lián)管道簡化為單一管段時，中間節(jié)點需水量稱作管段的沿線流量。這些沿線流量在管網建模中，常集中分配到管段的上端節(jié)點和下端節(jié)點，采用的方法是將沿線流量均分到兩端節(jié)點［1-2］。這種沿線流量均分方法優(yōu)點是簡單、易于操作，尤其在管段方向不定的情況下；缺點是在管段兩端節(jié)點分配流量后，引起的管段水頭損失與合并前在沿線流量下引起的管段水頭損失具有差異性。

因此，在管線沿線流量簡化處理中，需要解決管段沿線流量按怎樣的分配比例，才可滿足水力等效性的問題。針對該問題，1988年Hamberg等［3］利用謝才-曼寧公式，討論了管線中間單點出流量、等間距多點出流量下流量簡化處理問題。2013年Izquierdo等［4］利用達西-魏斯巴赫公式討論了管線內等間距等出流的流量簡化處理及水力坡度線特征。

已有文獻很少結合給水管網管段沿線流量簡化處理討論管段內水流時間的變化。通常管段內水流時間越長，管段內水質越可能變差。此外，隨著給水管網水力模擬的技術進步，希望計算的節(jié)點水頭和管段流速應更接近實際情況。因此需要從理論上分析給水管網管段沿線流量簡化處理后帶來的影響。本研究首先分析管段內一般沿線流量簡化處理前后的水頭與流行時間，然后分析管段內均勻出流的流量簡化處理前后的水頭與流行時間。

1 管段內一般沿線流量簡化

1.1 管道沿線水頭與流行時間

給水管網管道沿程水頭損失公式和水流時間計算公式分別為［5］

式（1）、（2）中：hf為管道沿程水頭損失，m；t為管道內水流時間，s；k、n、m為參數(shù)，不同計算公式對應的參數(shù)值見表1；Q為管道流量，m3·s－1；D為管道直徑，m；L為管道長度，m。

表1 管道沿程水頭損失計算公式參數(shù)［6］Tab.1 Pipe friction headloss formula parameters[6]

由式（1）和式（2）可以看出，當管道內流量Q和直徑D一定時，沿程水頭隨管道長度線性下降，水流時間隨管道長度線性增加。

設管段總長度為L，直徑為D，上端節(jié)點水頭為H0，進流量為Qin。管段內N個出流點距上端節(jié)點的距離依次為L1，L2，…，LN（0＜L1＜L2＜…＜LN＜L）。各出流點相應需水量分別為Qi（Qi＞0，i=1，2，…，N）。設QL為各出流點需水量之和，即QL=為沿線流量占比，0≤γ≤1。

忽略局部水頭損失，則距管段上端節(jié)點距離為0，L1，L2，…，LN，L處的水頭，以及從上端節(jié)點到各處的水流時間見表2。流量差異使組成管段的各管道內水頭和水流時間的變化斜率各不相同，整個管段上的水頭和水流時間將以多段折線形式變化。

表2 管段兩端和各出流點處的水頭，以及從管段上端節(jié)點至各處的水流時間Tab.2 Water heads of end nodes and outlets along the link,and traveling times from upstream nodes to each node or outlet

1.2 沿線流量分配到管端兩端

設管段沿線流量為QL，經簡化處理，將其分配到管段兩端。令α為管段沿線流量QL的上端節(jié)點流量分配比，則通過管段的流量為Qin-αQL=(1-αγ)Qin。于是管段下端節(jié)點的水頭H′(L)和水流時間T′(L)分別為

當將沿線流量平均分配到管段上下兩端節(jié)點時，上端節(jié)點流量分配比α=0.5，則由式（3）和式（4），得到管段下端節(jié)點的水頭H"(L)和水流時間T"(L)分別為

由式（3）～（6）可知，當沿線流量分配到管段兩端節(jié)點后，沿管段的水頭和水流時間將以線性形式變化，不再是多段折線形式。

1.3 示例1

如圖1所示，某管段進流量Qin=400 L·s－1，內部節(jié)點（1）距離上端節(jié)點（0）為250 m，需水量為60 L·s－1；內部節(jié)點（2）距離上端節(jié)點（0）為400 m，需水量 為180 L·s－1；于 是 可 知 沿 線 流 量 占 比γ=節(jié)點（0）處水頭H0=80 m，管段總長為500 m，海曾-威廉系數(shù)C取110，管段直徑D為500 mm。試分析將內部節(jié)點（1）和（2）的需水量分配至上端節(jié)點（0）和下端節(jié)點（3）情況下的管段水頭和流行時間。

圖1 內部兩個出流節(jié)點的管段示意圖Fig.1 Schematic diagram of two outlets along the link

1.3.1 簡化處理前水頭和水流時間

1.3.2 沿線流量等分后水頭和水流時間

當將內部節(jié)點（1）和節(jié)點（2）的流量等分至節(jié)點（0）和節(jié)點（3）后，由式（5）和式（6）得節(jié)點（3）水頭H"(500)=77.56 m，該值大于表3中所列水頭76.42 m，說明采用流量等分方法將會樂觀估計管段下端節(jié)點的水頭。

表3 各節(jié)點水頭和從節(jié)點（0）算起的水流時間Tab.3 Water heads of nodes and traveling times from node(0)

從節(jié)點（0）至節(jié)點（3）的水流時間T"(500)=350 s，該值大于表3中所列的流行時間332 s，說明采用流量等分方法時水在管段內的停留時間將增長，水質分析時可能得出較差的水質結果。

1.3.3 水力等效條件下水頭和水流時間

為滿足水力等效性，使沿線流量分配到管段兩端節(jié)點后下端節(jié)點（3）處的水頭與簡化處理前相等，則式（3）中H′(500)=76.42 m。解得管段上端節(jié)點流量分配比α=0.23。說明水力等效性條件下，節(jié)點（0）分配23%的沿線流量（56 L·s－1），節(jié)點（3）分配77%的沿線流量（184 L·s－1）。

由式（4），管段內水流時間計算為T′(500)=286 s，該值小于表3中所列的流行時間332 s，說明采用水力等效方式簡化處理后，水在管段內的停留時間縮短，會樂觀估計管段水質。

由示例1可知，采用沿線流量平均分配到管段兩端節(jié)點時，計算出的下端節(jié)點水頭和水流時間均與簡化處理前不同。水力等效條件盡管是在滿足管段下端節(jié)點水頭相同條件下的簡化處理，但計算出的水流時間與簡化處理前存在差異。此外水力等效性盡管滿足簡化處理前后管段上端和下端節(jié)點水頭相同，但由于管段沿程水頭在簡化處理前為多段折線下降，而在簡化處理后呈直線下降，所以造成管段內部的水頭在簡化處理前后產生差異。

2 管段內均勻出流的沿線流量簡化

2.1 水頭和流行時間

管段內均勻出流為沿線出流的特例。設管段內有N個等間距出流點，則各點需水量相同，為QL/N；管段兩端和各出流點到管段上端節(jié)點的距離依次為0，L/（N+1），2L/（N+1），…，NL/（N+1），L。令沿線流量占比γ=QL/Qin（0≤γ≤1），則管段兩端節(jié)點及各出流點水頭分別為

從上端節(jié)點至各出流點處的水流時間分別為

當N→∞時，將沿管段連續(xù)均勻出流。令x為到管段上端節(jié)點的距離，則沿線dx長度內水頭變化為

于是距離管段上端節(jié)點x（0≤x≤L）處的水頭為

同理可得當N→∞時，沿管段連續(xù)出流下距離管段上端節(jié)點x（0≤x≤L）處的水流時間為

由式（13）和式（15）可知，在沿管段連續(xù)均勻出流情況下，整個管段內水頭的下降和水流時間的增加均將以連續(xù)曲線形式變化。

2.2 管段上端節(jié)點流量分配比

為滿足水力等效性，在管段均勻出流情況下，使沿線流量分配到管段兩端節(jié)點后，下游節(jié)點處的水頭相等。由式（3）和式（9），令H′(L)=H(L)，得管段上端節(jié)點流量分配比為

式（17）說明，管段均勻出流水力等效簡化處理后，上端節(jié)點流量分配比α僅與管段內部出流點個數(shù)N和沿線流量占比γ相關，而與管道直徑D、管段上端節(jié)點水頭H0、管段進流量Qin、管段長度L和系數(shù)k無關。這樣，將管段沿線均勻出流的流量分配到兩端節(jié)點時，它不僅可用于穩(wěn)態(tài)模擬，也可用于H0、Qin隨時間變化的延時模擬。

當N→∞時，管段連續(xù)均勻出流，為滿足水力等效性，由式（3）和式（14），令H′(L)=H(L)，得

式中：γ∈(0，1]。當γ=1時，管段進流量Qin全部沿本管段配水，沒有向下游管段轉輸?shù)亩嘤嗔髁俊．敠谩?時，管段進流量Qin在本管段分配流量很少，絕大部分流量轉輸至下游管段。當n=2和n=1.852時，γ=1和γ→0時管段均勻出流的極限α值見表4。

表4 管段均勻出流的極限α值Tab.4 Limitsαof uniform flow along the link

式（17）中取不同N值和γ值，當n=2和n=1.852時，管段上游節(jié)點流量分配比α的變化情況見圖2。

由圖2看出：①當管段中出流點數(shù)量N不變時，上端節(jié)點流量分配比α隨沿線流量占比γ的增加而減小。α值總是小于極值0.5，即管段內出流量總是按下端節(jié)點流量較大，上端節(jié)點流量較小的方式進行水力等效簡化處理。γ值很小情況下，即管道內流量主要是轉輸?shù)较掠喂芏味皇窃诒竟芏窝鼐€配水時，α接近0.5。②當沿線流量占比γ不變時，上端節(jié)點流量分配比α隨管段中出流點數(shù)量N的增加而增加。③當上端節(jié)點流量分配比α相同時，沿線流量占比γ隨管段中出流點數(shù)量N的增加而增加。

2.3 示例2

設管段進流量Qin=100 L·s－1，長度L=1 000 m，直徑D=300 mm，上端節(jié)點處水頭H0=30 m；采用海曾-威廉公式計算沿程水頭損失，m=4.87，0.001 503。設沿線流量占比從0.01～0.90，N值分別取1，2，5和∞，試分析沿線均勻出流簡化處理后管段下端節(jié)點水頭和從管段起點到該點的水流時間。

由式（9）和式（12）計算沿線流量簡化處理前管段下端節(jié)點的水頭和水流時間。由式（17）獲得管段上端節(jié)點流量分配比后，再由式（3）和式（4）得到水力等效簡化處理后管段下端節(jié)點的水頭和流行時間。由式（5）和式（6）得到流量等分（即α=0.5）時下端節(jié)點水頭和流行時間。不同出流點個數(shù)N和不同沿線流量占比γ的計算結果如圖3、圖4所示。

由圖3和圖4看出，不同N值均勻出流情況下，簡化處理前后管段下端節(jié)點水頭和流行時間隨不同沿線流量比γ的變化趨勢相似。

圖3 不同N值管段下端節(jié)點水頭Fig.3 Downstream node water heads of the link at different N values

圖4 不同N值管段下端流行時間Fig.4 Downstream node traveling times of the link at different N values

（1）因為等分流量處理總是比水力等效簡化處理在管段下端節(jié)點分配的流量要少，管段下端節(jié)點處簡化前的水頭（水力等效簡化處理后的水頭）總是小于等分流量處理后的水頭。

（2）等分流量處理后的管段流行時間高于水力等效簡化處理的但低于簡化處理前的。

（3）當沿線流量占比γ＜0.2時，簡化處理前后管段下端節(jié)點水頭和管段流行時間相差不大。說明當管段內沿線分配流量遠小于它的轉輸流量時（在供水管網距離水廠較近的管段內比較常見），采用將沿線流量以水力等效方式或等分方式分配到管段的兩端節(jié)點具有合理性。當沿線流量占比γ值較大時（距離水廠較遠的供水管網末梢或多水廠供水分界線附近管段處比較常見），在簡化處理前后管段下端節(jié)點水頭和管段流行時間出現(xiàn)較大差異。

3 結論

相同直徑、材料和埋設年代的串聯(lián)管道合并是給水管網常用簡化方式之一。根據能量守恒原理，給出了單管段沿線流量利用水力等效方式或流量等分方式向兩端節(jié)點分配，簡化處理前后水頭損失和水流時間的計算方法。水力等效簡化處理前后僅保證了管段上端節(jié)點和下端節(jié)點的水頭是相同的，而管段內部各處的水頭通常存在差異。

管段沿線均勻出流流量經水力等效簡化處理后，上端節(jié)點流量分配比與管段內部出流點個數(shù)和沿線流量占比相關，而與管段直徑、上端節(jié)點水頭、管段進流量、管段長度和管段摩阻系數(shù)無關，使得該簡化處理方法適用于上端節(jié)點水頭或管段進流量具有變化性的延時模擬模型。

管段沿線均勻出流情況下，距離水廠較近的管段沿線流量（當管段內部流量占比較小時）適合采用水力等效方式（或流量等分方式）簡化處理，沿線流量簡化處理前后管段水頭損失和水流時間差距均較小。而對于距離水廠較遠的管網末梢或多水廠供水分界線附近（當管段內部沿線流量占比較大時），經流量等分方式簡化處理后的管段水頭損失，以及經水力等效方式（或流量等分方式）簡化處理后的水流時間，與簡化前相比差距均較大。