海洋泥沙濃度及輸移的研究方法綜述

王文柳 孫蓓蓓 丁圓圓 焦嬌

【摘要】本文在綜述了海洋泥沙濃度及輸移的研究方法的基礎上，分析了各種研究方法的優劣，并詳細介紹了MIKE-21數值計算與分析軟件，以期為相關科研工作者提供幫助。

【關鍵字】泥沙濃度；研究方法；MIKE-21

海洋泥沙濃度及輸運過程的變化，直接影響著岸灘地貌的演變過程，尤其是近岸工程的實施，或改變水動力環境，進而影響泥沙輸運過程。李建超等[1]研究發現，1855年黃河北徙人渤海后，黃海泥沙補給驟減，改變了水沙之間原有的平衡狀態，使蘇北海岸沖淤分布格局發生改變。對于海洋泥沙濃度及輸運研究，應用較多的是建立模型來反演輸運與沉積關系[2]。而建立泥沙模型的方法有很多。其中比較常見的有VB程序模式和遙感。但是，這些方法都有些不足之處。VB程序過于簡單，無法應對海洋環境比較惡劣的海域沉積研究，而遙感主要是偏向于宏觀觀測，對于具體的海洋沉積情況不能做出更加準確的判斷。本文介紹了各種泥沙濃度及輸運的研究方法，重點敘述MIKE-21泥沙輸移模型的應用。

1. GIS與遙感

GIS與遙感作為20世紀快速盛行的地理系軟件，應用十分廣泛。遙感技術具有實時、連續、準確地獲取大范圍地表信息的能力，而GIS 因為其有力的分析能力被廣泛應用。特別在海岸帶地區，海水和陸地具有明顯的光譜特性差異，因而在遙感圖像上可以非常清楚地反映出來[3]；而區域性地表測量資料往往由于水陸性質不同及其他限制因素取得較為困難，因此綜合GIS與遙感，使二者相互補充、相互利用，所取得的結果具有較高時間與空間分辨率，并具有直觀、準確、易于應用的特點[4]。遙感技術憑借著獨特的宏觀優勢和實時精確的特點，在近岸泥沙輸運中運用廣泛。余佳 [5]應用MODIS L1B數據完成了對于黃海泥沙懸浮體的數值反演，建立了單波段懸浮體反演模型。通過建立這種懸浮體的反演模型，可以清晰地發現遙感波段與實測數據的各種相關性，進一步地建立懸浮體泥沙濃度反演模型，定量地對懸浮體泥沙與其他因素的變化關系進行研究。

2. 現場激光粒度儀

近年來，隨著觀測技術的不斷進步，現場激光粒度分析儀（LISST）被廣泛應用，用來實時地完成對泥沙顆粒級配與濃度的測量，其具有測驗效率高、垂向分辨率高、實現分粒級測量泥沙懸浮體濃度的能力強等優點[1]。LISST在測定水體懸浮體體積分數時會將所有顆粒都計入最終結果中，如含有較多礦物聚合體或生物碎屑時會有較大的LISST數值，在測量粒徑的上界，產生“上升尾”現象。因此，先對數據進行處理，將明顯異常值（如具有較高抽濾數據，但較低LISST數據點）剔除掉。為了使結果更直觀且能利用于傳統的基于質量濃度的懸浮體分析，且削弱大顆粒、低密度懸浮體對LISST觀測結果的影響，將LISST數據與抽濾數據對比，將體積濃度轉化為質量濃度。

3. MIKE-21

3.1 泥沙輸移模塊簡介

泥沙輸移模塊隸屬于MIKE-21水動力模型里面的ST 模塊。MIKE -21的ST模塊根據水流作用或波流的共同作用來計算輸沙率。在波流共同作用時，模型利用DHI模型STP的波周期公式計算輸沙率。如果只有水流作用，用戶可選擇多種常用的輸沙計算公式。初始的淤積/沖刷利用泥沙的質量守衡確定。

3.2 泥沙模塊步驟

在運用泥沙模塊前，要明確要解決的問題；選用適當的模型；選擇維度（一維、二維或三維）；選用適當的模塊、水動力、泥沙或波浪；模擬區域范圍、區域大小、空間步長；坐標投影（模型轉角）。

3.2.1 數據收集與前處理

地形數據（C-map，XYZ，GIS，CAD，紙制海圖等）；邊界實測數據和驗證數據；風場等參數數據。數據驗證利用MIKE.21軟件驗證其合理性與一致性。

1. 搭建模型

輸入模型參數有地形文件、邊界文件、糙率值、渦粘系數時間步長、風 場、干濕水深、初始水位等。

2. 運行模型

3.3.2 泥沙模型實例分析

MIKE.21軟件對于泥沙的分析，學者們在遼寧黃海在花園口工程海域實測水文測驗資料和地形資料的基礎上，從平面二維水沙模型的基本方程和基本理論出發，用MIKE21建立大尺度黃渤海潮流模型，采用嵌套的方法建立了計算區域內的小區域潮流耦合波浪的泥沙數學模型，三角單元有限體積法并對潮流運動和波浪潮流作用下泥沙的輸移進行了數值模擬。

在運用軟件時要注意對模型中關鍵參數（如剪切應力，水流挾沙能力等）進行重點研究，使地區模型能更好地適應實地情況，并根據計算的實測資料驗證計算模式和資料處理的合理性。為工程設計、日后規劃發展提供必要的參考依據。對實地地區進行了流場分析，解明其海域漲、落潮流特性。在此基礎上建立了波浪和潮流作用下的泥沙數學模型，對懸沙濃度和地形演變兩方面進行驗證，若懸沙濃度計算值和實測值吻合較為良好，則該泥沙計算模型沖淤是合理的。除此之外，還要關注地形演變反映出的海域的海床原始等深線是否處于動態平衡，這也是影響計算模型合理性的重要因素之一。完成后就可以在此基礎上得出年淤積厚度，海岸線短期內向近岸移動不是很大。通過分析海床長期沖淤處于動態平衡狀態，在處于動態平衡的海岸上修建建筑物不會引起大范圍的地形演變，只需注意建筑物造成的局部沖刷和淤積問題。這樣在只需要相關的實測資料就可以通過建立數字計算模型，快速有效地了解測定海域的泥沙輸運與沉積狀況，從而制定更加有效的措施。

4. 總結

泥沙沉積作為一種在近海地區起巨大作用的沉積分類，對自然環境與人類生活和建設都有著非凡的影響。我國沿海泥沙沉積問題一直是人們關注的話題，也有眾多學者采取各種方法進行研究。隨著快捷、高效的MIKE.21軟件技術的發展，可以通過運用泥沙模型分析方法對泥沙進行分析與模擬，對中長期的變化趨勢進行預測，從而擬定解決措施。同時，在泥沙沉積研究方面GIS、遙感與激光粒度儀測泥沙粒徑法也有著不俗的作用。

參考文獻：

[1]李建超，喬璐璐，李廣雪，等. 基于LISST數據的冬季南黃海懸浮體分布[J]. 海洋地質與第四紀地質，2013，05：13-25.

[2]王嶸，張永戰，夏非，等. 南黃海輻射沙脊群海域底質粒度特征及其輸運趨勢[J]. 海洋地質與第四紀地質，2012，06：1-8.

[3]樊輝，黃海軍. 黃、東海二類水體春季表觀光譜特性與表層懸浮體濃度反演模式[J]. 海洋與湖沼，2010，02：161-166.

[4]李海宇，王穎. GIS與遙感支持下的南黃海輻射沙脊群現代演變趨勢分析[J]. 海洋科學，2002，09：61-65

[5]余佳，王厚杰，畢乃雙，等. MODIS L1B數據的黃海懸浮體季節性分布的反演[J]. 海洋地質與第四紀地質，2014，01：1-9.

基金項目：江蘇省高等學校大學生實踐創新創業訓練計劃項目（201410324023Y）