基于Flood Area 模型的模袋混凝土襯砌渠道糙率系數(shù)優(yōu)化方案初探

張海玲　霍軼珍　郭彥芬　王文達

摘 要：本文以河套灌區(qū)建設一分干模袋混凝土襯砌渠道為研究對象，引用Flood Area模型，將糙率系數(shù)的設置做區(qū)分化處理;對渠道區(qū)域進行劃分和賦值，在經(jīng)驗參數(shù)的基礎上，結(jié)合建設一分干模袋混凝土襯砌渠道實際情況，設置柵格糙率系數(shù)，增加糙率系數(shù)的精度。結(jié)果表明：考慮模擬中糙率系數(shù)的動態(tài)性以及渠道糙率系數(shù)較坡面糙率系數(shù)更為復雜的特點，將糙率系數(shù)的設置做區(qū)分化處理，改進后的方法簡化了率定流程，以解決模擬流速普遍偏低的問題。利用高程與糙率系數(shù)的函數(shù)關系，以較全面、真實地仿真模擬模袋混凝土襯砌渠道的流場。

關鍵詞：Flood Area模型;模袋混凝土;糙率;TOPMODEL;區(qū)分化處理

中圖分類號：TV135.3? ? ? ? ? ?文獻標識碼：A? ? ? ? ? ? 文章編號：1006—7973（2019）01-0078-02

模袋混凝土是我國20世紀80年代初從國外引進的一項現(xiàn)澆混凝土新技術，它采用織物模袋做軟模具，通過混凝土泵將砂漿或混凝土充灌進模袋成型，起到護坡、護底、防滲等作用[1]。目前專家學者采用室內(nèi)試驗、現(xiàn)場原型測試與數(shù)值模擬相結(jié)合的方法，對糙率的取值及影響因素進行了深入研究。但是，目前河套灌區(qū)模袋混凝土襯砌渠道數(shù)值模擬還沒有引用Flood Area模型的相關研究。

劉鋒[2]利用Fluent軟件對植物壩覆蓋區(qū)域流場進行數(shù)值模擬，胡靜[3]對微厚透水物糙率處理方法進行了數(shù)值模擬。但是，由于渠道邊界層的內(nèi)部流動很復雜 以及現(xiàn)有計算機硬件限制， 目前數(shù)值模擬采用的計算方法還比較單一，絕大多數(shù)數(shù)值仿真均采用簡化的物理模型，沒有考慮真實的物理情況。如果單純采用FLUENT軟件模擬模袋混凝土襯砌渠道的流場時，存在較大的難度，而且計算時程序不易收斂。

本文以河套灌區(qū)建設一分干模袋混凝土襯砌渠道為研究對象，引用Flood Area模型，將糙率系數(shù)的設置做區(qū)分化處理以期解決模擬流速普遍偏低的問題，便于全面、真實地仿真模擬模袋混凝土襯砌渠道的流場。

1選擇渠道

建設一分干模袋混凝土襯砌渠道的測流渠段較平直、水流均勻，縱橫斷面比較規(guī)則、穩(wěn)定 [4]。故選取建設一分干模袋混凝土襯砌渠道作為全模袋混凝土襯砌斷面典型渠道。

建設一分干模袋混凝土襯砌渠道，樁號9+715～13+000為鋪設12cm厚土工全模袋混凝土防滲襯砌，混凝土強度等級為C25，其他參數(shù)見表1。

2 “Flood Area”計算原理

“Flood Area”為德國 Geomer公司開發(fā)的模型，采用水動力方法計算，充分利用1個柵格附件的8個單元，對周圍柵格單元的瀉入量采用manning-stricker 公式計算：

V=kst·r2/3·I1/2? ? ? ?（1）

式中：kst是反映明渠渠道粗糙度對水流影響的系數(shù);r是水力半徑;I為明渠渠道的坡度。

Flood Area針對各個時段水流的運行全過程，對于二維非恒定流，采用DEM水動力模型進行建模，考慮了坡度以及粗糙度等要素對水流數(shù)值模擬的影響，以柵格方式導出，全面地、動態(tài)地反映流速、水深、流量等參數(shù)和流場[5-6]。

“Flood Area”在每個時相的運行過程中，即運行時間與相應水流范圍和水深，以柵格形式表示，簡潔直觀，便于查看[7-8]。

3引用Flood Area模型改進方法的運用

針對傳統(tǒng)糙率系數(shù)設置方法引起模擬流速偏低的問題，優(yōu)化方法主要考慮模擬中糙率系數(shù)的動態(tài)變化以及渠道糙率系數(shù)較為復雜的特點，將糙率系數(shù)的設置做區(qū)分化處理[9]。

結(jié)合建設一分干模袋混凝土襯砌渠道，輸入柵格糙率系數(shù)，壁面處采用局部加密處理，將糙率系數(shù)假設為水深的函數(shù)，在明渠流動中糙率系數(shù)會隨著水深的增加而呈現(xiàn)減小的趨勢[10]，因此利用其與水深的變化關系進行賦值。在模擬過程中水深是隨著時間動態(tài)變化的，且 Flood Area 模型限制糙率系數(shù)只能作為常量輸入。因此，需要建立高程與糙率系數(shù)的變化關系，在忽略水力坡度影響的情況下，水流區(qū)域高程與水深呈負相關性關系，建立如下函數(shù)方程：

Kst= Kmin+ [（ Kmax- Kmin） /（ Emax- Emin） ]* （ Emax- E） ）[11]

其中，Kst是Strickler 糙率系數(shù)，Kmax代表糙率系數(shù)最大值，Kmin代表糙率系數(shù)最小值，E高程值，Emax是最大高程值，Emin是最小高程值。不同的模擬情景參數(shù)不同，Kmax、Kmin在實際模擬中作為率定參數(shù)進行糙率系數(shù)的調(diào)節(jié)。

4研究技術路線（如圖2）

5結(jié)論

以建設一分干模袋混凝土襯砌渠道為研究對象，引用Flood Area模型，將糙率系數(shù)的設置做區(qū)分化處理，改進后的方法簡化了率定流程，增加糙率系數(shù)的精度，以解決模擬流速普遍偏低的問題。

Flood Area模型中糙率系數(shù)作為常量輸入，對渠道區(qū)域與坡面區(qū)域進行劃分和賦值，結(jié)合建設一分干模袋混凝土襯砌渠道實際情況，獲取有效的下斷面信息，利用高程與糙率系數(shù)的函數(shù)關系，以全面、真實的仿真模擬模袋混凝土襯砌渠道的流場。

參考文獻：

[1]鄭子成，何淑勤，吳發(fā)啟，等.地表糙度與水力糙率系數(shù)的關系[J].山地學報，2004，22（2） ： 236-239.

[2]劉鋒，邱秀云，周著. 利用Fluent軟件對植物壩覆蓋區(qū)域流場進行數(shù)值模擬[J].水利水電科技進展，2010，30 （2）： 32-35.

[3]胡靜，魏清福，曾婧揚. 數(shù)值模擬中微厚透水物糙率處理方法探討[J].人民長江，2015，46 （17）： 40-42.

[4]何建京，王惠民.流動型態(tài)對曼寧糙率系數(shù)的影響研究[J]. 水文，2002，22（6） ： 22-24.

[5]張磊，王文，文明章.基于“Flood Area”模型的山洪災害精細化預警方法研究[J].復旦學報（自然科學版），2015，54（3）： 282-287.

[6]張科利，唐克麗.黃土坡面細溝侵蝕能力的水動力學試驗研究[J].土壤學報，2000，37（1）： 9-15.

[7]蘇布達，姜彤，郭業(yè)友，等.基于 GIS 柵格數(shù)據(jù)的洪水風險動態(tài)模擬模型及其應用[J].河海大學學報（ 自然科學版） ，2005，33（4）： 370-374.

[8]張明達，李蒙，戴叢蕊，等.基于 Flood Area 模型的云南山洪淹沒模擬研究[J].災害學，2016，31（1） ： 78-82.

[9]張洪江，北原曜，解明曙，等.晉西幾種不同土地利用狀況下糙率系數(shù)的研究[J].北京林業(yè)大學學報，1994，16（4） ： 86-91.

[10]孫麗.復式渠道水流特性數(shù)值模擬研究[D].江蘇：河海大學.2007.

[11]羅尊驊，林廣發(fā)，張容焱，等. Flood Area 模型糙率系數(shù)改進方法[J].海南師范大學學報（自然科學版），2018，31 （1）： 98-104.