利用積分原理計算非恒定流倒灌水量的探討

周天喜

摘 要：本文圍繞實際工作中遇到的河涌水通過截污管倒灌形成非恒定流的水量計算問題，采用積分原理尋求解決問題的方法，擬定了水量計算的積分公式，并對公式計算的可行性和精度進行了驗算和探討。

關鍵詞：積分；非恒定流；水量；流量；潮汐

0 前言

在通常的水量計算中，當忽略流速水頭損失和局部水頭損失時，對于恒定流常采用基本的流量計算公式“流量=流速×過水面積”；但當過水面積隨時間變化時，則屬于非恒定流的流量計算問題，需要對時間進行積分運算。本文就是針對這類非恒定流的計算進行探討。

1 問題提出的背景

位于廣州市珠江水系的馬涌流域截污管已投入使用多年，因防洪的需要，部分截污管網與河涌和珠江通過拍門連通；但經過長期的運行后，部分拍門已損壞和丟失，造成截污管與河涌水處于直接連通狀態，形成非正常的水體交換。

馬涌截污工程完成并運行一段時間后，馬涌截污系統內的實際水量是設計污水量的2倍以上，導致泵站在正常的工作條件下其水位的變化更多與河涌的水位變化一致，因此懷疑有河水倒灌入污水管網，增加了污水管網和泵站的工作負荷。為了查清馬涌截污管網與珠江水系連通點的情況，以及江水倒灌入截污管網的水量，需要對每個連通點進行監測和流量計算，然后累計得到總的倒灌水量。

然而截污管的倒灌流量并不是一個恒定值：當珠江水位在低潮位時，截污管管底高于外部水面，不存在倒灌現象；在一個潮汐周期，隨著潮位的上漲，當潮位高于截污管管底時才開始出現倒灌現象，且整個倒灌過程是一個動態變化的過程，灌入的水量隨著外部水位的上漲而增加（如圖1所示）。所以采用通常的水量計算公式“倒灌水量V=流量Q×時間t=過水面積S×水流速度v1×時間t”是不能精確的反映實際情況的，為了解決這個水量測定問題，需要采用積分原理進行計算。

該公式取水位上漲到達管底倒灌開始時為0時刻，V為k時刻累計倒灌水量，k的最大取值為水位從管底上升到管頂所經歷的時間。由于該積分公式解析過于復雜，將各參數具體數據代入公式簡化，可得其數值解。

3 實例分析與對比

傳統的倒灌水量計算方式通常采用“定時測量→分段估算→累加”的方式進行，本文以珠江某連通點實際監測數據為研究對象，通過傳統計算方式與積分公式分別計算某時段內累計倒灌水量，并進行分析對比。

表1通過計算各300N秒時刻（N=0、1、2、3……）的流量，用該流量乘以時間300s得出300N秒～300（N+1）秒時刻內倒灌的水量，再進行累加，得k時刻內總倒灌的水量，易知此結果比真實值略大。

當采用積分計算時，由于v1、v2隨時間的變化幅度較小，可認為v1、v2為恒值，故取算術平均值v1=0.15744m/s、

將傳統計算方式與積分公式計算方式所得的流量與時間關系繪圖（見圖3）：

如圖3所示，倒灌水量值等于流量曲線Q（t）與t軸所圍成的不規則圖形的面積。傳統的計算方法將該不規則圖形隨時間t進行粗略的等分，用各分段矩形面積代替梯形面積，再進行累加得出總面積；而采用積分方法的計算過程是將該不規則圖形隨時間t進行無限量細分，再累加面積的過程。顯而易見，采用積分公式所得結果比傳統計算方式所得結果更加接近于真實情況。

以2小時為例，通過傳統計算方法所得累積倒灌水量V傳統=464.370m3，通過積分公式所得累積倒灌水量V積分=391.915m3，兩者差值高達15.6%.由此可見，傳統的水量計算方法過于簡略，易嚴重造成統計值的不準確。

4 結語

通過以上分析與對比，利用積分計算倒灌水量是可行的，且相較于傳統計算方式能大幅減少計算工作量，提高工作效率，精度更接近于真實情況，能滿足工程要求，達到了預期的目標。我國沿海地區幅員廣闊，城市的截污系統大多存在類似的海水、河水倒灌的工程計算問題，開展這類方法的研究，具有很強的現實意義，能廣泛應用于沿海城市的環境治理、截污改造等領域。如何更好地解決倒灌流量計算問題，有待我們去探討和研究。本文旨在起到拋磚引玉的作用.由于時間倉促，應用實踐較少，不免存在不足之處，懇請各位同仁指出謬誤并提出寶貴意見，以期在今后類似的項目應用中取得更理想的效果。

參考文獻：

[1] 張志軍，熊德之.微積分及其應用[M]. 北京：科學出版社， 2007.

[2] 馬崢，陳紅勛.曲面擬合和曲面積分在流量計算中的應用[J]. 排灌機械， 2004， （3）.

[3] 劉曉波.管道流體流量測量實驗系統構建與分析[D]. 太原：中北大學， 2014.

[4] 李亮，方習勇.用曲線積分法計算洪水流量的分析[J]， 安徽地質， 2002， 12（4）.