張明志　仇銀豪　何一飛等

摘要應用MIKE21模型，模擬膠州灣紅島濕地在陸域淡水和潮汐流交匯作用下水體鹽分梯度的分布規(guī)律，并預測濕地水體及土壤的環(huán)境狀況，為紅島濕地整體規(guī)劃設計中的植栽選型提供必要的理論依據和技術手法。

關鍵詞 紅島濕地；MIKE21；潮汐；鹽度；植栽

中圖分類號S11；Q948文獻標識碼A文章編號0517-6611（2015）16-179-04

Application of MIKE 21 Model in the Wetland Planting of Hongdao，Jiaozhou Bay

ZHANG Mingzhi， QIU Yinhao， HE Yifei et al

（Orient United Design Group， Shanghai 200233）

AbstractThis article tried to use MIKE21 model to imitate the distribution of salinity gradient， with the effect of river freshwater and ocean tidal in Hongdao wetland， Jiaozhou Bay. By imitated the salinity gradient distribution， the conditions of water and soil environment in Hongdao wetland could be predicted， and theories and technologies were able be provided for landscape plant selection in future integrated planning of Hongdao wetland.

Key wordsHongdao wetland； MIKE21； Tide； Salinity； Plants

濱海濕地作為陸地生態(tài)系統(tǒng)和海洋生態(tài)系統(tǒng)的交錯過渡帶，目前正受到日益增強的圍墾、污染、生物資源過度利用等因素威脅，濕地面積正快速減少和逐步破碎化。如何保護濱海濕地，恢復健康的濕地生態(tài)系統(tǒng)，已成為解決濕地、海洋問題的關鍵所在。

紅島濕地（圖1）位于山東省青島市紅島經濟區(qū)，西側為

青島市“母親河”大沽河，南側為膠州灣，面積約25 km2，是目前膠州灣最大最重要的濱海濕地。受常年人類生產活動的影響，鹽田、蝦池等人工濕地大量存在，場地已形成以大量埤塘型人工濕地為主，少量鹽沼濕地和河口濕地構成的體系。紅島濕地的建設旨在解決以下四大主要問題：淡水嚴重不足；土壤鹽漬化程度加??；埤塘濕地比重過大且污染嚴重；物

圖1紅島濕地區(qū)位

種結構單一，生物多樣性受到威脅。

1 MIKE21模型簡介

丹麥水力研究所（DHI）開發(fā)的平面二維數學模型MIKE21是一款應用廣泛的商業(yè)模型。目前該軟件在國內的應用發(fā)展很快，并在一些大型工程中廣泛應用。MIKE21是DHI開發(fā)的系列水動力學軟件之一，屬于平面二維自由表面流模型，可用于模擬河流、湖泊、河口、海灣、海岸及海洋的水流、波浪、泥沙及環(huán)境場，為水利工程、海岸規(guī)劃等提供完備、有效的設計條件和參數[1]。

2 MIKE21HD及Transport計算原理

MIKE21軟件中的水動力學模塊（HD模塊）是其最核心的基礎模塊，可以模擬水位和水流變化，以及任何忽略分層的二維自由表面流[1]?；贖D模塊的水動力條件及運行結果，應用Transport模塊模擬物質組分在時間和空間上的輸移特性。

2.1控制方程組

對于水平尺度遠大于垂直尺度的情況，由于水深、流速等水力參數沿垂直方向的變化比沿水平方向的變化要小得多，因此，將三維流動的控制方程沿水深積分，并取水深平均，可得到沿水深平均的二維淺水流動質量和動量守恒控制方程組。其連續(xù)性方程、X和Y方向動量方程可分別表示為[2]：

ζt+px+qy=ht

pt+x（p2H）+y（pqH）+gHζx+gpp2+q2C2H2-1ρ[x（Hτxx）+y（Hτxy）]-fq-fw|W|Wx=0

qt+y（q2H）+x（pqH）+gHζy+gpp2+q2C2H2-1ρ[y（Hτyy）+x（Hτxy）]-fq-fw|W|Wy=0

式中，H為水深，H=h+ ζ，其中h、ζ分別為凈水深和水位；t為時間；p、q為X、Y方向上的流通通量；C為謝才系數；g為重力加速度，f為科氏力系數；ρ為水的密度；W、Wx、Wy為風速及在X、Y方向上的分量；fw為風阻力系數；τxx、τxy、τyy為有效剪切力分量。

模型在二維X、Y方向求解的對流擴散方程式為：

ct+pcx+qcy=Dx2cx2+Dy2cy2

式中，c為組分質量濃度；t為時間；p、q為X、Y方向上的流通通量；Dx、Dy為組分在X、Y方向上的分子擴散系數。

2.2定解條件

2.2.1

邊界條件。 開邊界：按潮位過程和流速或流量過程計算； 閉邊界：利用岸壁法，取法向不可入條件，即法向流速為零[1]。

2.2.2

初始條件。 模型計算的初始條件有3種：①恒定初始值在全模擬區(qū)域從靜止狀態(tài)開始，取流速和水位為定值；②對模擬區(qū)域內各網格點或區(qū)域指定不同初始水位和流速；③“熱啟動”（hot start）方式，即利用之前計算結果作為本次計算的初始條件[2]。

2.2.3

計算時間步長。 根據Courant數（CR）選取時間步長：

Δt=CRΔxgh

式中，CR為Courant數；h為凈水深；△t為時間步長；△x為空間步長；g為重力加速度。

3 MIKE21模型應用

3.1水文氣象

膠州灣潮汐屬規(guī)則半日潮類型，設計高潮位1.921 m，設計低潮位-2.049 m（以1985年國家高程基準面起算）。膠州灣全年以風生浪為主，灣內常見波的波高不足0.5 m。秋、冬兩季全海域盛行偏北向浪，以NW向出現頻率最多（17%～31%）；春、夏兩季盛行偏南向浪，以SE向出現頻率為最多（30%～33%）。

3.2參數設置

采用MIKE21 FM非結構三角形網格構建模型，網格數量約3萬個，數值模擬計算采用MIKE21模型的水動力模塊（HD）及Transport模塊。

3.2.1

邊界條件。膠州灣開邊界采用天文潮實測潮位（圖2），大沽河入流邊界流量數據采用該河段的多年平均流量，內陸其他入流量主要來自市政雨水排放，依據紅島經濟區(qū)防洪排澇規(guī)劃取值。

3.2.2

初始條件。濕地底部高程及水位采用實測數據（圖2）。

3.2.3

參數處理。風向采用不利風向東南風，風速取模擬月份平均風速；干濕邊界、阻力系數等參數取模型自定義經驗數值。

3.2.4

鹽度設置。內河入流初始鹽度取0.1 g/kg，膠州灣開邊界初始鹽度取35 g/kg。

圖22014年膠州灣潮汐表及邊界條件

3.3模擬結果

3.3.1

流場模擬。根據膠州灣長系列潮汐表進行分析，8月份膠州灣潮位偏高，因此模型計算膠州灣開邊界選取2014年8月上旬實測潮位值，計算時長取15 d。紅島濕地水動力模擬結果（圖3）表明：低潮位及平均潮位時，紅島濕地受潮汐流影響較小，當潮位高于1.05 m時，海水漫入紅島濕地，濕地受潮汐漫灘影響時間約占總模擬時長的1/3。

3.3.2

鹽度模擬。選取2014年8月15日低潮位（-2.02 m）及高潮位（2.03 m）時段。相應潮位對應水體鹽分布見圖4。

根據濱海地區(qū)海水入侵及土壤鹽漬化研究[3]、濱海鹽漬土水鹽運動模擬研究[4]等相關研究內容以及濕地現狀用地分析，按《巖土工程勘察規(guī)范》（GB50021-2001）（2009年版）鹽漬土含鹽量分類標準[5]，將紅島濕地土壤鹽漬化程度分為以下4類（圖5）。

3.4結果驗證

根據膠州灣紅島濕地土樣含鹽量檢測報告，結合濕地現狀植被，復核MIKE21模擬預測結果。

3.4.1

土樣分析檢測。中國海洋大學土工試驗室出具的紅島濕地中心區(qū)兩處土樣的檢測報告[6]（表1）。

由表可知，濕地中部取樣點（1#）不同深度的土壤含鹽量在0.78%～2.04%，臨海區(qū)取樣點（2#）在1.6%～3.77%，按《巖土工程勘察規(guī)范》（GB50021-2001），1#、2#點均為中鹽漬土，與模擬預測結果基本一致，誤差控制在合理范圍內。

3.4.2

濕地現狀植被。通過現場勘察，遠離海灣的濕地北部區(qū)域水生植物較為豐富，景觀性較好，且有少量耕地，以玉米種植為主。濕地中部以鹽田為主（東風鹽場），地表裸露，植被稀少。濕地南側臨海，以蝦池、海參池為主，四周較高臺地生長有少量蘆葦及堿生植物（圖6）。

總體而言，紅島濕地現狀植被分布與MIKE21模擬得到的鹽度分布規(guī)律較為一致，部分區(qū)域植被分布有所出入，或與地形數據缺失及MIKE21模型網格劃分精度不足有關。總體而言，模擬預測結果對植物栽植選型具有參考價值。

圖3濕地流場

圖4水體鹽度分布

圖5濕地鹽漬區(qū)分類

43卷16期

張明志等MIKE21模型在膠州灣紅島濕地植栽中的應用研究

表1土樣腐蝕性分析成果

序號土樣編號取樣深度∥m電導率ms/cm易溶鹽含量mg/kg含鹽量%

11#-10.3～0.42.227 763.000.78

21#-20.7～0.86.1920 403.002.04

31#-31.7～1.85.7917 910.001.79

42#-10.3～0.45.9715 955.001.60

52#-20.8～1.012.5237 689.003.77

62#-32.0～2.18.9727 918.002.79

4 植物栽植設計

植物作為生產者，在濕地的恢復、重建過程中發(fā)揮著重大作用。在濱海濕地鹽漬化的背景下，為保障植物前期成活率及后期生長效果，科學分析濕地鹽分梯度分布對植栽具有重要的影響。通過MIKE21模型對濕地鹽分梯度的模擬、分析和預測，該研究將從陸生、水生植物兩個方面對濕地植物進行合理優(yōu)化配置。

4.1陸生植物

樹木的耐鹽能力是指造林后1～3年內，幼樹對土壤鹽堿的適應性，是鹽堿地上樹木忍受鹽漬化程度并產生產量的能力。把樹木生長受到鹽堿抑制，但不明顯降低樹木成活率及生長量時的土壤含鹽量作為該樹種的耐鹽能力。

不同物種具有不同的耐鹽能力，同時耐鹽能力也與樹齡大小、樹木強弱、土壤質地、含水率及鹽分種類有關。一般將樹木的耐鹽能力劃分為強、中、弱3級[7]，并將相應樹木[8-9]應用到紅島濕地中（表2）。

圖6現狀植被分布

表2樹種耐鹽能力分級

耐鹽力等級含鹽量%出苗率%植物

弱0.2～0.3≥70毛白楊、垂柳、小葉楊等

中0.3～0.550～70刺槐、苦楝、臭椿、紫穗槐、合歡、杞柳等

強0.5～0.7≤50柳、枸杞、檉柳等

重≥0.7堿蓬、獐毛、馬藺、結縷草、珊瑚菜等

4.2水生植物

水生植物作為水生態(tài)系統(tǒng)的重要部分，可通過發(fā)達的枝葉和根系形成天然過濾層，吸附、分解和轉化水中污染物質，促進水域養(yǎng)分平衡，兼具生態(tài)、經濟、景觀、養(yǎng)殖及氣候效益等多種功能。

紅島濕地在陸域淡水及潮汐流的交匯下，形成多梯度鹽分的水體結構。通過MIKE21模型的前期分析，并結合水生植物生長習性，在濕地水分鹽度低于0.2%的淡水及微咸區(qū)配置蘆葦、茭草、慈姑、菖蒲、美人蕉、芡實等挺、浮水植物，以及菹草、黑藻、苦草、眼子菜等沉水植物；濕地水分鹽度不高于1.0%的微咸水種植蘆葦、千屈菜、蓖齒眼子菜、金魚藻、水鱉等水生植物；高于1.0%的咸水種植堿蓬、鹽角草、川蔓藻等水生植物。

同時通過以植物蛋白（輪葉黑藻、苦草、伊樂藻等沉水植物莖葉）代替動物蛋白作為餌料及多漁業(yè)品種混養(yǎng)、輪養(yǎng)、輪捕輪放的生態(tài)養(yǎng)殖模式等途徑，豐富埤塘種群結構，減輕傳統(tǒng)人工埤塘的污染負荷。

5 結語

該研究通過MIKE21模型的水動力模塊及Transport模塊，科學地模擬、分析和預測膠州灣紅島濕地水體及土壤的鹽分梯度分布，為紅島濕地植被恢復提供必要的理論依據。通過多層次、多品種的濕地植物，提高淡水利用率，削減污染負荷和豐富生物種群結構，以期達到構建健康、完整、自我持續(xù)的濕地生態(tài)系統(tǒng)的目的，為解決濱海濕地及海洋問題找到一條行之有效的途徑。

參考文獻

[1]

安永寧，楊琨，王瑩，等.MIKE21模型在海洋工程研究中的應用[J].海岸工程，2013，32（3）：1-10.

[2] 許婷.丹麥MIKE21模型概述及應用實例[J].水利科技與經濟， 2010，16（8）：867-869.

[3] 張國明，顧衛(wèi)，吳之正，等.渤海灣風暴潮倒灌對沿岸農田土壤鹽分的影響[J].地球科學進展，2006，21（2）：157-160.

[4] 趙耕毛，劉兆普，陳銘達，等.海水倒灌濱海鹽漬土的水鹽運動模擬研究[J].中國農業(yè)科學， 2003， 36（6）：686-680.

[5] 中華人民共和國建設部.巖土工程勘察規(guī)范（GB50021-2001）[S].北京：中國建筑工業(yè)出版社.2009：82-83.

[6] 中國海洋大學土工實驗室.濕地公園工程土壤腐蝕性分析成果表[R].青島，2013.

[7] 王治國，張云龍，劉徐師，等.林業(yè)生態(tài)工程學[M].北京：中國林業(yè)出版社，2000：397-399.

[8] 趙可夫，李法曾，樊守金，等.中國的鹽生植物[J].植物學通報，1999，16（3）：201-207.

[9] 李紅麗，丁國棟，董智，等.中捷農場濱海鹽堿地立地類型劃分及其植被恢復技術[J].中國水土保持科學， 2010，8（5）：86-91.