廣東省綜合單位線在深圳茅洲河流域的應用與評價

張徐杰 程開宇 岳青華 張發鴻 富強

摘要：隨著城市建設的發展，不透水覆蓋面積不斷擴大，用于洪水計算的單位線方法已不能完全適用于城市匯流計算，現行單位線方法推算的基礎條件和城市化后的流域匯流條件已有較大區別。以深圳茅洲河流域為例，對廣東省綜合單位線在茅洲河流域的適用性進行了驗證和評價，提出了廣東省綜合單位線在城市化發展背景下可能存在的問題和需要改進的方向。

關鍵詞：廣東省綜合單位線;洪水計算;城市匯流;適用性評價;茅洲河流域;深圳

中圖法分類號：TV214文獻標志碼：A DOI：10.15974/j.enki.slsdkb.2019.05.004

目前，城市小流域水文計算方法有推理公式、單位線和水文模型等方法，其中以推理公式和單位線方法應用較多。但隨著城市建設的發展，不透水覆蓋面積和排水系統不斷擴大，單位線方法已不能完全適用于城市匯流的計算，這主要是因為現行單位線方法推算的基礎條件和城市化后的流域匯流條件已有較大區別。近年來，很多學者對暴雨推求洪水的方法進行了新的探討。黃啟有等建立了不同形狀小流域推理公式法設計流量表達式;黃國如等探討了城市中小流域設計洪峰流量計算方法的適應性;李整等對廣東省改進綜合單位線法在洪峰流量計算中的應用進行了研究。

本次研究以深圳茅洲河流域為例，基于中國電建集團華東勘測設計研究院有限公司（以下簡稱“華東院”）在茅洲河流域設立的實測水文雨量站點數據，對廣東省綜合單位線在茅洲河流域的適用性進行了驗證和評價，提出了廣東省綜合單位線在城市化發展背景下可能存在的問題和需要改進的方向。

1流域概況

茅洲河發源于深圳市境內的羊臺山北麓，地跨深圳、東莞兩市，流經石巖、光明、公明、松崗、沙井和東莞長安鎮等，在沙井民主村匯入伶仃洋。茅洲河總流域面積388.23km，干流全長約41.61km，其中深圳境內流域面積310.85km，東莞境內流域面積77.38km。由于茅洲河上游石巖河及其支流匯入石巖水庫后，石巖水庫下泄經溢洪道排入鐵崗水庫，直接排入下游珠江口，并不進入茅洲河流域，因此除石巖水庫流域面積（44km）外，茅洲河流域面積為344.23km，干流全長30.69km。感潮河段長約13km，下游河口段11.4km為深圳市與東莞市界河。茅洲河口水系見圖1。

2廣東省綜合單位線

廣東省綜合單位線是通過深入研究納希瞬時單位線方法，汲取國內外經驗，提出的一套具有廣東特色的綜合單位線。選用了廣東省內集水面積為2.3-950km的50個水文站共639場雨洪對應資料進行分析綜合，資料使用到1983年。其基本思路主要包括以下兩點。

（1）對于解決水文資料短缺地區中小流域的設計洪水計算問題，綜合單位線法是一種可行的辦法。該法通過綜合流域特征指標與單位線要素或參數，得出經驗公式和圖表，供相似地區移用。

（2）為消除不同面積和不同計算時段長度對綜合單位線的影響，采用無因次單位線ui-Xi，是一種較好的方法。具體做法是將各站納希單位線的峰頂對齊，求出一條平均單位線，再轉換成無因次單位線ui-xi，最后再按全省分區分類平均，得出分區分類的典型無因次單位線供設計使用。

廣東省綜合單位線示意見圖2。其中，時段單位線的定義如下：河流上某一斷面的時段單位線是指由t△內凈雨強度不變，集水區內空間分布均勻的單位凈雨（即1mm）在出口斷面形成的地表徑流過程線。

3廣東省綜合單位線適用性評價

3.1基礎資料

華東院在茅洲河流域設立了2個雨量站和2個水文站（見圖1），其中沙井雨量站位于流域下游位置，樓村雨量站位于流域中上游位置。2個水文站均位于茅洲河流域中上游，因此采用樓村雨量站和2個水文站的雨洪資料，進行雨洪相關分析。

利用2016年6-12月的實測降雨、流量資料，提取其中的所有洪水過程與對應降雨過程，組成72 h歷時的雨洪分析資料，干流水文站和新陂頭水文站分別有8個場次的雨洪資料，見圖3-4。由圖可見，干流水文站和新陂頭水文站第8個場次的洪水過程均存在較大的基流，與降雨過程不完全對應，因此不作為樣本進行研究。

3.2單位線法初始模擬成果

采用廣東省綜合單位線法，根據實測降雨過程來模擬流量過程。干流水文站和新陂頭水文站的7個場次流量模擬結果見圖3-4。由于每場洪水過程僅有72個流量值，計算單場次洪水過程的評價指標如納什效率系數（NsE）較不合適。因此，對于干流水文站和新陂頭水文站，分別將7個場次的流量過程合并起來后作為一個模擬和實測流量過程計算NSE，干流水文站和新陂頭水文站的NSE值分別為0.57和0.36。由此可見，廣東省綜合單位線法在干流水文站和新陂頭水文站模擬的流量成果仍有待提高。

3.3單位線法參數率定

對于廣東省綜合單位線法中的兩個重要參數（單位線滯時m1和24h穩滲率f24），本次研究采用拉丁超立方抽樣法（抽樣10000次）進行參數率定，采用NSE為目標函數。干流水文站和新陂頭水文站的NSE分布情況見圖5。由圖可見，干流水文站NSE大于0.5的模擬次數超過70%，最大NSE為0.63;新陂頭水文站NSE大于0.5的模擬次數不足20%，最大NSE為0.55。兩個站點的最大NSE雖沒有達到較理想值，但相對于廣東省綜合單位線法初始參數值模擬的成果，均有一定程度的提高。

干流水文站和新陂頭水文站達到最大目標函數值的前10組模擬成果分別見表l～2。由表可見：對于干流水文站，24h穩滲率f24率定結果大概在3.7-3.9范圍內，單位線滯時m1率定結果大概在4.13-4.14范圍內;對于新陂頭水文站，24h穩滲率f24率定結果大概在4.3～4.5范圍內，單位線滯時m1率定結果大概在3.56-3.58范圍內。

4結論與討論

根據廣東省綜合單位線法，茅洲河流域24h穩滲率f24=4.5mm，干流水文站單位線滯時m1初始值為3.56，新陂頭水文站單位線滯時m1初始值為2.58。因此，對于24h穩滲率f24，根據綜合單位線法，f24=4.5mm。從模型率定成果看，新陂頭水文站的f24值按綜合單位線法較合理，而干流水文站的f24值按綜合單位線法可能會偏大。對于單位線滯時m1，干流水文站和新陂頭水文站按綜合單位線法計算的m1值均可能偏小。

以上研究的初步結論是基于茅洲河流域實測水文資料而得出的，因此，水文資料的可靠性將是本次研究的關鍵。目前，華東院設立的干流水文站和支流水文站主要是根據2016年統計的水位一流量關系曲線，通過測量水位來推求流量的。而2016年的水位一流量關系，則是根據幾次測流成果而擬定的。根據測量成果，干流水文站的實測水位一流量點據相對較合理，擬合關系相對較好，而新陂頭水文站的實測水位一流量點據并沒有較理想的相關關系。因此，新陂頭水文站的流量數據可靠性存在一定問題，有待進一步分析和驗證。干流水文站的水位一流量關系相對較穩定，但實測數據仍存在偏少問題，同樣需要更多的實測資料進一步驗證其合理性和可靠性。

另外，本研究中采用的所有洪水場次中，可能會存在上游水庫泄洪的情況，但受于限制，沒有掌握上游水庫的泄洪情況和相關資料。如洪水過程中存在水庫泄洪的情況，將對研究造成一定影響。