波流共同作用下河口泥沙運動及地形演變數值模擬研究進展

呂曉瑩+文藝+陳昌亮+林偉龍

【摘 要】河口區域動力環境復雜，河流動力、海洋動力在此交匯。波浪、潮流是泥沙運動的主要動力，河口區域泥沙活動劇烈，長期來看影響河口地形發育演變，因此對波流耦合作用下河口泥沙運動及地形演變進行數值模擬具有重要的實用價值及工程意義。本文介紹了在河口區域波流共同作用數值模擬、在波流共同作用下泥沙運動及地形演變模擬方面的研究進展。

【關鍵詞】波流共同作用；泥沙運動；數值模擬；地形演變

河口區域動力環境復雜，河流動力、海洋動力在此交匯，泥沙活動劇烈。波浪、潮流是泥沙運動的主要動力，“波浪掀沙、潮流輸沙”是波流共同作用下的河口泥沙運動的重要機制，在波流作用下泥沙起動、輸運、沉積、再懸浮，長期來看影響著河口地形的發育演變。目前國內外許多地區已開展關于波流共同作用對河口演變影響的研究。

一、國外研究進展

Wright（1977，1980）描述并分析了波浪作用下的河口地貌特征及河口沉積過程[1-2]。Howard N. Southgate（1995）研究了波浪對中到長周期（幾周到幾十年）地貌演變產生的影響，為分析波序列以及波浪對海床和海灘平面發育演變的影響提供了方法[3]。Bernabeuet al.（2003）將波浪與潮流耦合，建立了一個模擬海灘剖面地貌變化的模型[4]。Matthieuet al.（2013）以印度海的馬約特島（Mayotte Islands）為例，關注了無河流注入的海岸的泥沙運移過程和地貌演變特征。他們通過分析波浪在近岸處能量衰減的特征，研究了地貌對波浪作用的響應[5]。Saruwatariet al.（2013）在奧克尼群島（Orkney Islands）應用SWAN模型對波流共同作用下海洋能與電能轉化進行研究。通過改變流的強度，以及相應的波向，分析不同波氣候時潮能和波能的變化[6]。Bola?os et al.（2014）結合實測數據，使用先進的模型對潮控河口進行了對波流共同作用下波場和流場特征的模擬。水動力模型建立時考慮了潮流、徑流以及氣壓引起的正壓力與斜壓力，并與波譜模型和湍流模型耦合。模擬結果表明河口處潮流對波浪的改變明顯，有效波高和波周期主要由隨時間變化的水深控制，同時波周期還受到流產生的多普勒頻移影響[7]。

二、國內研究進展

國內對波流共同作用的研究主要集中在黃河、南渡江、長江等大河河口。朱大奎（1997）應用大量實測波浪資料，從波浪動力的角度揭示了波浪作用對黃河三角洲地形演變的影響[8]。虞志英等（2002）調查分析了廢黃河河口水下三角洲的水動力及泥沙特征，對水下三角洲的平坦海床、水下斜坡和近岸淺灘三個地貌單元分別計算波流和潮流作用下底部泥沙的橫向沖刷率，并預測了水下三角洲的地貌演變模式。龔文平、王寶燦（1998）；陳沈良（1998）；戴志軍、陳子粲、歐素英（2000）；羅憲林，李春初，羅章仁（2000）等從不同方面分析研究了南渡江三角洲海岸的波浪動力過程，及其引起的沿岸泥沙運移，探討了泥沙在波浪作用下的沿岸運動狀態與地貌發育演變規律[9-[16][17][18]13]。劉紅，何青，吉曉強等（2008）通過對海灘剖面、懸沙粒度與表層沉積物的實測資料分析，研究了長江口崇明東灘在波流共同作用下潮灘剖面沉積物以及地貌的分異規律，并分析了沉積物的運移機制。他們認為隨波浪的分選作用由破波帶向兩側逐漸變差，表層沉積物中值粒徑隨之變細。破波帶內泥沙以“波浪掀沙”機制引起的分選運移為主，破波帶兩側處泥沙以潮流的“平流輸沙”機制為主[14]。苗麗敏、楊世倫、朱琴等（2016）以長江口與杭州灣兩大水系交匯處的南匯潮灘為例，探討了波流聯合剪切應力作用下潮灘懸沙濃度和懸沙輸運對風暴事件的響應過程及其動力機制[15]。

三、波流耦合數學模型研究進展

自20世紀60年代起，潮流的研究開始以數值模擬為主[16]，在Longuet-Higgins等（1960，1962）提出輻射應力理論后[17-18]，波浪對潮流驅動作用的研究取得突破性的進展，為波流共同作用的研究提供了理論基礎。近年來，波流耦合數學模型蓬勃發展，在河口研究方面得到了廣泛的應用。河口海岸泥沙運動的機理為“波浪掀沙、潮流輸沙”。由于波浪的周期較潮汐周期小，數學模型對波流共同作用通常使用兩種不同的方法對波或流加以概化處理。其一是將周期變化的潮流概化為具有某一特征的恒定流，與波浪運動方程疊加來模擬波流運動結構在短時間內發生的變化；另一種方法是把波浪過程概化為一個潮周期中具有平均意義的波浪流要素，與潮流運動方程疊加，計算長時段的水流及泥沙運動的變化。后者能更好地體現懸沙濃度變化和地形沖淤過程的周期性規律。王厚杰，李瑞杰（1999）在前人研究的基礎上考慮了波生流場對波浪場的影響，建立了一個討論近岸波流耦合作用的二維數學模型，并采用Dingemans的地形對模型進行驗證分析，結果表明近岸處波流隨時耦合，在研究近岸泥沙運移、岸灘演變和近岸工程設計時， 應予以充分的考慮[19]。Prandle et al（2000）以Holderness海岸為例，建立了波流耦合的二維懸沙模型，分析了近岸處潮流和波浪的相互作用及二者的基本特性[20]。王彪等（2012）考慮了波浪和潮流對泥沙運動的共同作用，概化波浪運動為潮周期中時均的波浪流分布場，耦合潮流運動方程、泥沙運動方程與波浪摩阻力和波浪挾沙力等要素，建立了一套波流共同作用下的泥沙運動模型[21]。孫麗等（2014）將波浪的輻射應力和底切應力同潮流運動方程和懸沙輸運方程耦合，在大連灣建立波流共同作用的二維懸沙數學模型，并在純潮流、常浪和極端天氣三種狀況下分析對比波浪對流場、懸沙濃度和底床沖淤演變的影響[22]。何用等（2014）對珠江口泥沙的起動、絮凝以及波浪作用下的挾沙力等問題進行探討，建立了珠江口波浪、潮流、鹽度、泥沙的耦合模型，除了增加波浪作用下泥沙的起動判別外，還考慮了鹽度對泥沙絮凝的影響[23]。李大鳴等（2014）考慮了潮流、波浪和泥沙的相互作用機制，將浪場產生的“波浪輻射應力”和“波流挾沙能力”作為對潮流場和泥沙運動的重要影響因子，構建了一套考慮波浪動力的潮流泥沙數學模型，并應用于對錦州港港區潮流場和地形沖淤變化的模擬，取得了良好的結果[16]。許婷等（2015）采用波流共同作用下二維水沙模型模擬了韓江河口附近大型離岸人工島群工程建設后潮流和泥沙的運動及海床沖淤變化規律[24]。張心鳳（2017）以為例，利用波流耦合作用下三維懸沙數學模型模擬珠海高欄港區深水航道在大風天的泥沙運動和驟淤規律，并預測了50年一遇情況下航道的驟淤情況[25]。endprint

為滿足河口研究的需求，許多成熟的數學模型被開發為軟件的形式，已被應用到實際的研究和工程中去，如DELFT3D（劉曙光，2010[26]），MIKE21（許婷，2010[27]），SWAN（王殿志，2004[28]；郝宇馳，陶建華，2007[29]；王道龍，華鋒，江志輝，2010[30]）等。TELEMAC數值模擬系統基于非結構化網格，采用有限元方法，對地形復雜的河口海岸區域十分適用。Andreas Malcherek（2000）將TELEMAC-2D模型用于Weser河口潮能耗散的模擬[31]；R.Kopmann， M. Markofsky（2000）利用TELEMAC-3D建立了三維水質模型[32]；J. Eric Jones， Alan M. Davies（2005）通過對英國西海岸潮汐的模擬，對比了有限差和有限元（TELEMAC）方法模擬結果的不同[33]。J. Eric Jones， Alan M. Davies（2006）利用TELEMAC模擬了愛爾蘭海的風生環流[34]； Samaras et al.（2013）以意大利南部的布林迪西口岸為研究區域，對比TELEMAC和MIKE-21兩種數學模型對波流共同作用和波浪傳播的模擬，發現在破波帶處對波高模擬的差異最為明顯[35]。范平易（2010）利用TELEMAC-2D模型建立了長江河口的二維水動力模型，并以蘇通大橋為研究對象，模擬建橋后對橋區河段通航的影響[36]。

波流共同作用下泥沙運動及地形發育演變數值模擬是前沿課題之一，由于河口區域人類活動劇烈、動力環境復雜，且泥沙運動受波流耦合作用顯著，因此更要加強對河口區域波流及泥沙運動的數值模擬研究，為河口海岸工程提供重要指導。

參考文獻：

[1]Wright L. D.Thom G.Higgins R.J.Wave Influences on River Mouth Depositional Process： Examples from Australia and Papua New Guinen[J].Estuarine and coastal marine science ，1980，03：263-277.

[2]Wright L. D.Morphodynamics of a Wave Dominated River Mouth[J].River mouth morphodynamics，1977，04：1721-1737.

[3]Howard N.Southgate.The effects of wave chronology on medium and long term coastal morphology[J].Coastal Engineering，1995，26：251-270.

[4]Bernabeu A. M.， Medina R.， Vidal C. A morphological model of the beach pro?le integrating wave and tidal influences[J].Marine Geology，2003，197：95-116.

[5]Matthieu Jeanson， Edward J.，Anthony.，Francket Dolique al.Wave characteristics and morphological variations of pocket beaches in a coral reef-lagoon setting， Mayotte Island， Indian Ocean[J].Geomorphology，2013，182：190-209 .

[6]Ayumi Saruwatari，David M. Ingram， Lucy Cradden.Wave-current interaction effects on marine energy converters[J].Ocean Engineering，2013，73：106-118.

[7]Bola?os Rodolfo， Brown Jennifer M.， Souza Alejandro J.Wave–current interactions in a tide dominated estuary[J].Continental Shelf Research，2014，87：109-123.

[8]朱大奎，李平.波浪在黃河三角洲形成中的作用[J].海洋地質與第四紀地質，1997（02）：40-47.

[9]虞志英，張國安，金镠等.波流共同作用下廢黃河河口水下三角洲地形演變預測模式[J].海洋與湖沼，2002，33（06）：583-590.

[10]龔文平，王寶燦.南渡江三角洲北岸的海岸演變及其機制分析[J].海洋學報，1998（03）：140-148.

[11]陳沈良，龔文平，王寶燦.南渡江三角洲海岸泥沙縱向運移與岸灘演變的響應[J].海洋湖沼通報，1998（01）：23-32.

[12]戴志軍，陳子燊，歐素英.海南島南渡江三角洲海岸演變的波浪作用分析[J].臺灣海峽，2000，19（04）：413-418.

[13]羅憲林，李春初，羅章仁.海南島南渡江三角洲的廢棄與侵蝕[J].海洋學報，2000，22（03）：55-60.

[14]劉紅，何青，吉曉強等.波流共同作用下潮灘剖面沉積物和地貌分異規律——以長江口崇明東灘為例[J].沉積學報，2008，26（05）：833-843.

[15]苗麗敏，楊世倫，朱琴等.風暴過程中潮灘懸沙濃度和懸沙輸運的變化及其動力機制——以長江三角洲南匯潮灘為例[J].海洋學報，2016，38（05）：158-165.endprint

[16]李大鳴，歐陽錫鈺，潘番等.考慮波浪輻射應力的潮流泥沙數學模型[J].海洋通報，2014（06）：703-711.

[17]Stewart R. W. Longuet-Higgins M S.Changes in the form of short gravity waves on long waves and tidal currents[J]. Journal of Fluid Mechanics，1960，8（04）565-583.

[18]Stewart R. W. Longuet Higgins M S. Radiation stress and mass transport in gravity waves on long waves ，with application to ‘surf beats[J].Journal of Fluid Mechanics，1962，13（04）：481-504.

[19]王厚杰，李瑞杰.近岸區域波流耦合作用的數學模型[J].海洋湖沼通報，1999，3：1-9.

[20]David Prandle Julia C. Wolf.Tide，wave and suspended sediment modeling on an open coastal—Holderness[J].Coastal Engineering，2000，41（01）：237-267.

[21]王彪，沈永明，王亮.波浪潮流共同作用下長興島海區二維懸沙輸運數值模擬研究[J].水動力學研究與進展A輯，2012，27（03）：321-330.

[22]孫麗，王永學，王國玉.波流共同作用下大連灣海床沖淤變化的數值研究[J].中國水運（下半月），2014（04）：91-95.

[23]何用，高時友，徐峰俊等.珠江口潮流泥沙數學模型關鍵問題探討[J].泥沙研究，2014（04）：60-66.

[24]許婷，劉國亭，溫春鵬.波流共同作用下大型離岸人工島群工程海床沖淤預測[J].海岸工程，2015，34（04）：11-22.

[25]張心鳳.大風浪作用下河口區深水航道驟淤預測模式研究[J].水動力學研究與進展，2017，32（03）：365-373.

[26]劉曙光，郁微微，匡翠萍.三峽工程對長江口南匯邊灘近期演變影響初步預測[J].同濟大學學報，2010，38（05）：679-684.

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

[28]王殿志，張慶河，時鐘.渤海灣風浪場的數值模擬[J].海洋通報，2004，23（05）：10-17.

[29]郝宇馳，陶建華.用SWAN模型模擬近岸波破碎[J].港工技術，2007（01）：4-6.

[30]王道龍，華鋒，江志輝.SWAN近岸海浪模式在遼東灣的應用[J].海洋科學進展，2010，28（03）：25-291.

[31]Malcherek Andreas.Application of TELEMAC-2D in a Narrow Estuarine Tributary[J].H ydrological processes， 2000， 14（13）： 2293-2300.

[32]Kopmann R.， Markofsky M.Three-dimensional water quality modelling with TELEMAC-3D[J]. Hydrological processes，2000，14（13）：2279-2292.

[33]J. Eric Jones， Davies Alan M.An intercomparison between finite difference and finite element （TELEMAC） approaches to modelling west coast of Britain tides[J].Ocean Dynamics，2005，55（3-4）：178-198.

[34]J. Eric Jones， Davies Alan M.Application of a finite element model （TELEMAC） to computing the wind induced response of the Irish Sea[J].Continental Shelf Research，2006，26（12）：1519-1541.

[35]A. G. Samaras， M Vacchi， R. Archettiet al.Wave and hydrodynamics modelling in coastal areas with TELEMAC and MIKE21，2013.

[36]范平易.橋區流場數值模擬及建橋對通航的影響研究.碩士論文.上海交通大學，2010.endprint