二維泥沙數學模型在唐河倒虹吸工程設計中的應用

張衛東

(河北省水利水電勘測設計研究院 天津 300250)

1 工程概況

南水北調中線工程唐河倒虹吸位于定州市辛莊村附近,在其西北約 150m處穿越唐河,渠線走向NE57°,與唐河基本正交。工程以上控制流域面積4578km2,距西大洋水庫 16.4km,距京廣鐵路14.1km。工程等別為 I等。主體長度 1155m,設計流量為 135m3/s;加大流量為 160m3/s。

唐河倒虹吸附近河床為寬淺游蕩性沙質河床,河寬 3～5km,河段比降為 0.13%,灘地寬闊,目前主槽位置是 1963年洪水以后形成的。河道除兩側上部有部分砂壤土外,表層均為砂層,厚度 4～6m,平均粒徑為 0.24～0.44m;下部為礫石層,平均粒徑為 6.0～8.7m,河床覆蓋層厚度大于30m。

2 1963年洪水概況

1963年洪水是唐河有實測資料以來的最大洪水,西大洋水庫至京廣鐵路橋區間暴雨歷時長,強度大,出現 “平地漫流”,鐵路橋附近 6km路基路軌被沖毀沖走。1992年河北省水利水電勘測設計研究院對唐河倒虹吸工程至鐵路橋河段進行了“63.8”歷史洪水調查,共調查到洪痕 3處;搜集到 1967年調繪的 1∶10000地形圖上標測的 “63.8”洪痕 10處。

采用該模型進行泥沙模擬計算,確定 “63.8”洪峰流量是關鍵。由于西大洋水庫至鐵路橋河段寬淺游蕩,沖淤變化劇烈,雖有一些洪痕仍難以直接推算 “63.8”洪峰流量,在此采用了鐵路橋調查流量 4100m3/s按河段長度比率折算出工程斷面“63.8”洪峰流量為 2980m3/s,后用多種方法復核后確定工程斷面 “63.8”洪峰流量為 3000m3/s。

此次建立唐河交叉斷面泥沙數學模型采用“63.8”洪峰流量為 3000m3/s,選用了模型范圍內“63.8”洪痕 5處。

3 泥沙數學模型理論基礎

采用天津大學王尚義教授提出的泥沙數學模型,參與河床演變過程的物理參數共計 9個,系統要求建立 9個方程式,包括水流連續方程式、河床演變方程式、水流阻力方程式、水流挾沙能力方程式、懸沙連續性方程式、床面穩定性條件式等。

其中糙率公式為

式中:A0為水流計算調試系數,是模型計算的主要調試系數。

4 模型建立

4.1 模型計算范圍及模型剖分

參照南水北調中線總干渠的線路布置和跨河建筑物型式,模型上邊界為釣魚臺村北,下邊界為大流村南。模擬范圍最大寬度 5.5km,最小寬度1.6km,長度 10km,總干渠距上邊界5.5km,距下邊界 4.5km,范圍約 40km2。

為適應模型計算范圍內的地形、鐵路橋、堤防等的布置,提高計算精度以適應有限元計算方法的要求,模型中采用三角形網格剖分,最小步長80m,最大步長 200m,天然情況共劃分節點 1203個,三角形單元 2271個,工程布置后節點數和單元數略有增加。

4.2 模型的調試與驗證

模型的調試主要依據 “63.8”洪峰流量(3000m3/s)和5處洪痕點。在 1967年地形圖和其他有關資料的基礎上,調試模型上下邊界水位及糙率系數 A0。根據網格和水流條件的要求,水流計算步長取 0.5s。調試結果以符合多數洪痕為原則,并對模型中各節點的流速、水位及床面穩定性等參數進行合理性分析。

為使模型進一步符合 “63.8”洪水的實際情況,在對地形條件進行調整后,重新調試各種模型參數,以便將模型成果應用于工程方案。

通過調試,糙率范圍為 0.02～0.10, “63.8”洪水驗證值和洪痕值比較接近,多數差值在 0.10m左右?！?3.8”洪水計算淹沒范圍與實際洪水基本一致,京廣鐵路橋處河床變形與實際情況基本符合。洪痕點驗證結果見表 1。

表1 “63.8”洪水洪痕點驗證成果 單位:m

在此采用泥沙數學模型計算的結果與 “63.8”洪水實際情況基本吻合,表明采用該模型進行數值模擬是可靠的,可用于南水北調中線工程總干渠與唐河交叉工程的分析研究。

5 模型應用

將該模型應用到唐河工程前后,對無工程情況、有工程情況 (工程斷面寬度 1000m)下 50年一遇、100年一遇和 300年一遇洪水分別進行模型計算。計算內容包括有工程、無工程兩個方案三種洪水標準的水位及沖刷等情況的分析。

模型計算是在現有河床地形地物條件下進行的。通過對比 1967年地形圖和 1995年地形圖,可以看出河道地形無明顯變化,且 “96.8”洪水未出河槽,因此現狀河道行洪條件基本上與 1963年洪水所造成的地形相同。

在模型計算中,通過直接移植模型驗證的糙率系數 A0和邊界條件,并且進行水流和沖淤變化的交替計算,保證水流和沖淤的雙重穩定。

有工程情況 (工程斷面寬度 1000m)方案,100年一遇、300年一遇洪水最大壅水高度分別為1.6m和 1.9m,均位于軸線上游 300m處。工程斷面下游水位明顯跌落,并很快與天然水位相銜接。由于唐河河床組成為中細紗,抗沖能力弱,工程斷面附近形成明顯的沖坑。300年一遇洪水工程斷面附近最大沖深為 5.2m,沖刷情況比較嚴重。計算成果見表 2。

表2 唐河倒虹吸有工程、無工程方案軸線水位、沖刷深度成果

6 結語

通過模型的建立、調試與驗證和應用,可得出以下結論。

(1)根據南水北調中線總干渠唐河工程提供的水文、地質、地形條件和工程方案的具體要求,采用二維泥沙數學模型進行分析是適用的。

(2)結合該模型過去的實踐成果與經驗,利用1963年實際洪水資料和行洪后1967年的地形資料建立模型,并據此分析研究的工程問題和得出的成果反映了唐河交叉斷面天然和工程條件下的水力特性。

(3)該模型經過驗證后,應用于計算河道沖淤情況、沿程最大水位變化情況及壅水等,計算結果合理、可靠,為唐河倒虹吸工程設計提供了依據。

1 楊國錄 .河流數學模型 [M].北京:海洋出版社,1993.

2 白玉川.水流泥沙水質數學模型理論及應用 [M].天津:天津大學出版社,2005.

3 張紅武.沖積河床糙率問題的探討 [J].武漢水利電力學院學報,1986,(3).

4 周建軍,林秉南,王連祥.平面二維泥沙數學模型研究及其應用 [J].水利學報,1993,(11).

5 韋直林 .黃河下游水流泥沙平面二維數學模型研究[J].武漢水利電力學院學報,1990,(6).