張偉超, 宋 策, 郭夢京, 任 雷

(西安理工大學 陜西省西北旱區(qū)生態(tài)水利工程重點實驗室, 陜西 西安 710048)

基于MIKE建模的城市生態(tài)公園行洪能力分析

張偉超, 宋 策, 郭夢京, 任 雷

(西安理工大學 陜西省西北旱區(qū)生態(tài)水利工程重點實驗室, 陜西 西安 710048)

[目的] 預測并分析工程對河道行洪能力變化的影響，為城市生態(tài)公園工程防洪影響評價提供研究技術方法。[方法] 基于MIKE 21建立城市生態(tài)公園工程河段的平面二維水流數學模型，以渭河寶雞段為例，選取1954型洪水過程驗證模型合理性，模擬3種典型流量下河道流場和水位的變化特征，討論城市生態(tài)公園修建后對河道行洪能力的安全影響。[結果] 3種典型流量洪水下沿線水位均低于堤頂高程，對城市行洪安全不構成威脅，因此行洪安全可以保證。但由于河道行洪寬度變窄左岸流速增大，加劇了對岸堤的沖刷，采取岸堤加固措施是極為必要的。[結論] 河道形狀的改變會影響水的流態(tài)，水流速度過大會影響河堤結構穩(wěn)定性，最終影響河道的行洪能力。

MIKE模型; 行洪能力; 流場分析; 城市生態(tài)公園

文獻參數： 張偉超, 宋策, 郭夢京, 等.基于MIKE建模的城市生態(tài)公園行洪能力分析[J].水土保持通報，2017,37(1)：128-131.DOI:10.13961/j.cnki.stbctb.2017.01.023； Zhang Weichao, Song Ce, Guo Mengjing, et al. Flood discharging capacity of urban ecological park based on MIKE modeling[J]. Bulletin of Soil and Water Conservation, 2017,37(1)：128-131.DOI:10.13961/j.cnki.stbctb.2017.01.023

城市生態(tài)公園是近年來備受關注的一個研究領域[1-3]。城市生態(tài)公園經過合理的開發(fā)利用，形成了集生態(tài)保護和休閑于一體的綜合性公園[4]，不但可以改善現在城市普遍面臨的水體污染、水環(huán)境質量變差等問題，而且可以充分體現其生態(tài)環(huán)境保護和社會生態(tài)經濟效益。但是城市生態(tài)公園的建設必然會使河道受到影響，導致河道行洪斷面減小，水位抬高，對河道行洪安全帶來一定的風險。城市生態(tài)公園與其他天然濕地最大不同在于其緊鄰城市[5-6]，考慮到其特殊的地理位置因素和安全因素，在城市生態(tài)公園建設時，對其安全性有著非常高的要求。由于公園的修建將對河道行洪能力和流場產生影響，因此其是否會對河堤抗洪能力產生影響是工程設計中最關心的問題。河道行洪能力分析常用的數學模擬軟件主要有MIKE，FLUENT，SMS等，其中MIKE軟件應用較為廣泛，成功進行了河流水源地的水質風險評估、水庫下游淹沒區(qū)模，電廠擴建溫排水等方面的模擬研究[7-9]。

基于此，本文擬采用MIKE平面二維水流數學模型，研究渭河寶雞段城市生態(tài)公園修建后，河道洪水過程中的流場狀況，河道沖刷變形情況，以期為河道防洪影響評價提供理論依據。

1 材料與方法

天然水流是三維的，需要用三維的Navier-Stokes方程來描述。對于具有600 m堤距的渭河寶雞段河道來說，由于流動在垂向的尺度比平面尺度小得多，物理量沿水深分布不均勻的影響是次要的，動水壓力分布與靜水壓力十分接近，因此可以用靜壓假定對流動控制方程進行簡化，使用沿水深方向平均的淺水方程及平面二維水流數學模型來計算河道內的流場。這里采用以雷諾平均化N-S方程為基礎的平面二維水流數學模型。

1.1 基本方程

求解水流運動的基本方程為：

(1) 水流連續(xù)性方程:

(1)

式中：ξ(x,y,t)——水流表面波動(m)；p——壓力；q——流量；SS——單位水平面積的源量值(m/s)。

(2) 水流動量方程。

①x方向動量方程:

(2)

②y方向動量方程：

(3)

式中：u——速度分量；h——水深(m)；Fx,Fy——波浪應力所產生的分量；Six,Siy——在X軸和Y軸方向的源匯量值(m2/s2)；C——Chezy阻力常數(m1/2/s)；E(x,y)——旋渦或動量擴散系數(m2/s)；f——風的摩阻力常數(m1/2/s)；Vx,Vy(x,y,t) ——在X軸和Y軸方向上的風速分量(m2/s)； Ω——旋轉角度(rad/s);g——重力加速度(m/s2);V——流速(m/s)。

1.2 地形條件

地形由渭河寶雞段城市生態(tài)公園工程沿線實測河道地形資料獲得，并依據此地形建立模型計算網格，采用三角網格，確定計算網格大小為100 m2，節(jié)點數為10 884，網格數為21 143。

1.3 初始及邊界條件

模型初始條件的變量包括流量、流速、水位，計算區(qū)域內流量和水位為計算相應的初值，流速為零。即：u(ξ,η,t0)=0，v(ξ,η,t0)=0，w(ξ,η,t0)=0。

模型計算時上游進口邊界采用流量條件，下游出口邊界采用水位條件，河床糙率取值根據河床結構取0.033。

1.4 模型驗證

根據河道上下游斷面資料，采用曼寧公式計算得出水位流量關系，代入模型進行計算，其計算結果與實測資料基本一致，說明本模型的水動力學理論符合實際，可用來模擬不同工況下河道的流場[10-14]。根據河道實測地形資料和河道平面圖將其劃分為9個斷面，對本模型計算的水面線成果進行合理性分析。

由表1可以看出，在D4(550m)斷面與D6(1 750m)斷面之間，河道主槽縮窄，河道過流斷面變小，且縮窄段長約1.2km，河道水面比降因壅水在D5以上變化平緩；經過D6斷面，河道在突然收縮之后又快速拓寬，拓寬段長約400m，河道水位在D5(1 300m)與D6(1 650m)區(qū)間內出現跌水，D7(2 200m)斷面與D9(2 850m)斷面之間平面地形變化不大，河道過流斷面順直，水位變化相對比較平穩(wěn)，說明計算水面線形態(tài)符合水面線與河槽平面形態(tài)的對應關系，本模型計算結果合理。此外，通過模型計算的百年一遇流量洪水水面線和水文局百年流量實測基本一致，證明了模型的準確性，可以進行城市生態(tài)公園的行洪能力模擬與分析。

表1 百年一遇流量條件下水面線對比

注: ① 水文局百年流量觀測水面線(1972年)； ② 模型計算百年一遇流量洪水水面線。

2 結果與分析

為了研究渭河寶雞段城市生態(tài)公園工程沿線河道在不同流量量級條件下的水位、流場，分析渭河寶雞段城市生態(tài)公園工程的防洪影響，擬定如下工況進行模擬計算。 ① 設計流量(3 000 m3/s)條件下模擬河道的流場； ② 1954型特征洪水(5 030 m3/s)條件下模擬河道的流場； ③ 百年一遇流量(7 260 m3/s)條件下模擬河道的流場，分別模擬3種工況的河道水位。

城市生態(tài)公園子堤設計流量為3 000 m3/s,對此工況下行洪能力進行模擬，模擬結果如圖1所示。由圖1可知，設計流量洪水水位沿線均低于子堤堤頂高程，表明設計流量洪水發(fā)生時不會漫過子堤，對城市行洪安全不構成威脅，因此設計條件下子堤行洪是安全的。若發(fā)生1954年型的大洪水，在距離子堤首段850～1 300 m范圍內子堤堤頂的高程不足，屆時洪水會經由此處出現漫溢，威脅到城市行洪的安全，更危險時可能會在此處出現決口，造成更大程度的損失。因此，應該對該段河堤進行加固與增高，以應對特征流量洪水帶來的威脅，保證河道行洪能力。

圖1 河道設計水位與堤防高程的對比分析

雖然河道縮窄使河道行洪能力在一定程度上受到影響，但研究表明百年一遇洪水水位并沒有超過渭河寶雞段左右大堤的實測堤頂高程(圖2)，大堤行洪安全沒有受到威脅，說明城市生態(tài)公園工程不對大堤行洪安全產生影響，城市生態(tài)公園在安全性方面具有可靠性。

圖2 河道100年一遇流量水位與大堤高程比分析

根據上述模型對3種工況下河道流場變化進行模擬計算，河道流場計算成果如圖3所示。

圖3 3種工況下流速及流場變化

由圖3可以看出，在工況1(設計流量)條件下，受子堤的影響河道過流面積變小，在D1—D3河段寬度逐漸減小但變化幅度不大，河道流速也受河寬影響有所波動；在D4—D6段河道主槽變窄，河道過流面積變小，流速逐漸增大，最大流速出現在2 250 m里程處約為4.7 m/s；D6-D8河段，河道主槽寬度增大，河道流速減小為2.1 m/s；D9斷面處河道主槽變窄河道流速增大到3.5 m/s。河道左岸流場密集，較大的流速會對左岸護坡造成影響，子堤附近的水位較低，大部分流場的流速約為2.0 m/s不會對子堤構成威脅。在工況2(1954型洪水)條件下，隨著上游來水增加，河道流速明顯加大，尤其在D4—D6河段，河道左岸流速最大約為5.4 m/s，與工況1相比，河道左岸流場密集度有所增加，高流速帶加長，對左岸護坡影響變大。對于工況3(100年一遇流量)情況，來水條件大于工況1和工況2，河道流速明顯增強，D4—D6河段左岸流速提高到約5.7 m/s，對左岸堤腳構成威脅，應對河岸堤腳進行加固措施，D9斷面處流速同樣很高，但不靠近堤岸，不會對堤岸造成影響。

3 討論與結論

建立MIKE平面二維水流數學模型，以渭河寶雞段城市生態(tài)公園為實例，通過對不同工況下的水位、流場，及行洪能力的模擬分析，評價工程對河道行洪能力的影響。根據模擬結果，在1954型洪水(Q=5 030 m3/s)和100年一遇流量下(Q=7 260 m3/s)，河道流速分別達到5.4和5.7 m/s,從流場分布圖可以明顯看出,河道左岸流場密集度有所增加，高流速帶加長，對左岸護坡影響變大。在D4—D6段由于河道寬度變窄流速到達峰值，因此在該河段應設置相應的護坡和堤腳保護措施，預防由于高流速洪水的沖刷，河堤被掏空。結果表明河道形狀的改變會影響水的流態(tài)，水流速度過大會影響河堤結構穩(wěn)定性最終影響河道的行洪能力，本研究成果可以為流域內類似工程的修建提供借鑒意義。

[1] Duan N, Liu X D, Dai J, et al. Evaluating the environmental impacts of an urban wetland park based on emergy accounting and life cycle assessment: A case study in Beijing[J]. Ecological Modelling, 2011,222(2):351-359.[2] Li Yufeng, Liu Hongyu, Zheng Nan, et al. Analysis of trophic status and its influence factors of different water body types in Xixi National Wetland Park, China[J]. Procedia Environmental Sciences, 2010(2):768-780.

[3] Li Xiaoping, Chen Manman, Anderson B C. Design and performance of a water quality treatment wetland in a public park in Shanghai, China[J]. Ecological Engineering, 2009,35(1):18-24.

[4] Johnston J. Nature Areas For City People[M]. London: Ecology Unit, 1990.

[5] 鄧毅.城市生態(tài)公園的發(fā)展及其概念之探討[J].中國園林,2003,19(12):51-53.

[6] 唐銘.西北地區(qū)城市生態(tài)公園評價體系研究:以蘭州銀灘生態(tài)公園為例[J].山東農業(yè)大學學報:自然科學版,2010,41(1):80-86.

[7] 宋策,譚奇林.水電工程干擾下飲用水水源地的水質風險評估[J].長江流域資源與環(huán)境,2013,22(1):59-65.

[8] 劉冀,李偉,張弛,等.碧流河水庫下游河道行洪能力及洪水淹沒模擬[J].中國農村水利水電,2008(2):22-25.

[9] 王慶改,戴文楠,趙曉宏,等.基于Mike21 FM的來賓電廠擴建工程溫排水數值模擬研究[J].環(huán)境科學研究,2009,22(3):332-336.

[10] 李大鳴,管永寬,李玲玲,等.蓄滯洪區(qū)洪水演進數學模型研究及應用[J].水利水運工程學報,2011(3):27-35.

[11] 李大鳴,林毅,徐亞男,等.河道、滯洪區(qū)洪水演進數學模型[J].天津大學學報:自然科學與工程技術版,2009,42(1):47-55.

[12] Savage B M, Johnson M C. Flow over ogee spillway: Physical and numerical model case study[J]. Journal of Hydraulic Engineering, 2001,127(8): 640-649.

[13] 何杰,辛文杰.含有紊動黏性項淺水方程的數值求解[J].水利水運工程學報,2010(3):95-100.

[14] 高慶先,胡銘,楊新興,等,湖泊流體動力學模型及其應用[J].環(huán)境科學研究,2001,14(6):29-32.

Flood Discharging Capacity of Urban Ecological Park Based on MIKE Modeling

ZHANG Weichao, SONG Ce, GUO Mengjing, REN Lei

(StateKeyLaboratoryofEco-hydraulicEngineeringofNorthwestAridAreainShaanxi,Xi’anUniversityofTechnology,Xi’an,Shaanxi710048,China)

[Objective] The study aims to predict and analyze the project impacts on the river flood discharging capacity, and provide methods for assessing flood controlling impacts of urban ecological park projects.[Methods] Based on MIKE 21, we constructed a 2-dimension water flow model for urban ecological park section in Weihe River in Baoji City, selected 1954 flooding process to verify the rationality of the model, and simulated three typical flow characteristics of flow field and changes in water level and discuss the urban ecological park impacts on river flood discharging capacity. [Results] Under the three typical flow flooding events, the water level along the river was below embankment and not a threat to urban flood discharging safety, and the flooding safety can be guaranteed. But as a result of flooding passage narrowing and increasing velocity, the erosion of the embankment was intensified, and reinforcement measures was urgently needed. [Conclusion] The river shape will affect the flow regime. The higher water flow speed will affect the stability of embankment structure thus influencing the river flood discharging capacity.

MIKE modeling; flood discharging capacity; flow flood analysis; urban ecology garden

2016-08-19

2016-09-09

國家重點基礎科學研究發(fā)展計劃(973)項目“西北旱區(qū)山洪致災機理及防治研究”(2012CB723201); 西安理工大學水利水電學院開放基金 (2016ZZKT-37)

張偉超(1992—),男(漢族),陜西省西安市人,碩士研究生,研究方向為水環(huán)境模擬。E-mail:lan_yan_lei@126.com。

郭夢京(1986—),男(漢族),山西省運城市人,博士,講師，研究方向為水環(huán)境模擬。E-mail:916874373@qq.com。

B

1000-288X(2017)01-0128-04

TV131.4