波流共同作用下河口泥沙運(yùn)動(dòng)及地形演變數(shù)值模擬研究進(jìn)展

呂曉瑩 文 藝 陳昌亮 林偉龍/.廣東省海洋發(fā)展規(guī)劃研究中心 .廣東省珠海航道局

波流共同作用下河口泥沙運(yùn)動(dòng)及地形演變數(shù)值模擬研究進(jìn)展

呂曉瑩1文 藝2陳昌亮1林偉龍1/1.廣東省海洋發(fā)展規(guī)劃研究中心 2.廣東省珠海航道局

河口區(qū)域動(dòng)力環(huán)境復(fù)雜，河流動(dòng)力、海洋動(dòng)力在此交匯。波浪、潮流是泥沙運(yùn)動(dòng)的主要?jiǎng)恿?，河口區(qū)域泥沙活動(dòng)劇烈，長(zhǎng)期來(lái)看影響河口地形發(fā)育演變，因此對(duì)波流耦合作用下河口泥沙運(yùn)動(dòng)及地形演變進(jìn)行數(shù)值模擬具有重要的實(shí)用價(jià)值及工程意義。本文介紹了在河口區(qū)域波流共同作用數(shù)值模擬、在波流共同作用下泥沙運(yùn)動(dòng)及地形演變模擬方面的研究進(jìn)展。

波流共同作用；泥沙運(yùn)動(dòng)；數(shù)值模擬；地形演變

河口區(qū)域動(dòng)力環(huán)境復(fù)雜，河流動(dòng)力、海洋動(dòng)力在此交匯，泥沙活動(dòng)劇烈。波浪、潮流是泥沙運(yùn)動(dòng)的主要?jiǎng)恿?，“波浪掀沙、潮流輸沙”是波流共同作用下的河口泥沙運(yùn)動(dòng)的重要機(jī)制，在波流作用下泥沙起動(dòng)、輸運(yùn)、沉積、再懸浮，長(zhǎng)期來(lái)看影響著河口地形的發(fā)育演變。目前國(guó)內(nèi)外許多地區(qū)已開(kāi)展關(guān)于波流共同作用對(duì)河口演變影響的研究。

一、國(guó)外研究進(jìn)展

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

二、國(guó)內(nèi)研究進(jìn)展

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

三、波流耦合數(shù)學(xué)模型研究進(jìn)展

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

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

波流共同作用下泥沙運(yùn)動(dòng)及地形發(fā)育演變數(shù)值模擬是前沿課題之一，由于河口區(qū)域人類活動(dòng)劇烈、動(dòng)力環(huán)境復(fù)雜，且泥沙運(yùn)動(dòng)受波流耦合作用顯著，因此更要加強(qiáng)對(duì)河口區(qū)域波流及泥沙運(yùn)動(dòng)的數(shù)值模擬研究，為河口海岸工程提供重要指導(dǎo)。

[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.Wavecurrent 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].海洋地質(zhì)與第四紀(jì)地質(zhì),1997（02）：40-47.

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

[10]龔文平,王寶燦.南渡江三角洲北岸的海岸演變及其機(jī)制分析[J].海洋學(xué)報(bào),1998（03）：140-148.

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

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

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

[14]劉紅,何青,吉曉強(qiáng)等.波流共同作用下潮灘剖面沉積物和地貌分異規(guī)律——以長(zhǎng)江口崇明東灘為例[J].沉積學(xué)報(bào),2008,26（05）：833-843.

[15]苗麗敏，楊世倫，朱琴等.風(fēng)暴過(guò)程中潮灘懸沙濃度和懸沙輸運(yùn)的變化及其動(dòng)力機(jī)制——以長(zhǎng)江三角洲南匯潮灘為例[J].海洋學(xué)報(bào),2016，38（05）：158-165.

[16]李大鳴,歐陽(yáng)錫鈺,潘番等.考慮波浪輻射應(yīng)力的潮流泥沙數(shù)學(xué)模型[J].海洋通報(bào),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]王厚杰,李瑞杰.近岸區(qū)域波流耦合作用的數(shù)學(xué)模型[J].海洋湖沼通報(bào),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]王彪,沈永明,王亮.波浪潮流共同作用下長(zhǎng)興島海區(qū)二維懸沙輸運(yùn)數(shù)值模擬研究[J].水動(dòng)力學(xué)研究與進(jìn)展A輯,2012,27（03）：321-330.

[22]孫麗,王永學(xué),王國(guó)玉.波流共同作用下大連灣海床沖淤變化的數(shù)值研究[J].中國(guó)水運(yùn)(下半月),2014（04）：91-95.

[23]何用,高時(shí)友,徐峰俊等.珠江口潮流泥沙數(shù)學(xué)模型關(guān)鍵問(wèn)題探討[J].泥沙研究,2014（04）：60-66.

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

[25]張心鳳.大風(fēng)浪作用下河口區(qū)深水航道驟淤預(yù)測(cè)模式研究[J].水動(dòng)力學(xué)研究與進(jìn)展,2017,32（03）：365-373.

[26]劉曙光,郁微微,匡翠萍.三峽工程對(duì)長(zhǎng)江口南匯邊灘近期演變影響初步預(yù)測(cè)[J].同濟(jì)大學(xué)學(xué)報(bào),2010,38（05）：679-684.

[27]許婷.丹麥MIKE21 模型概述及應(yīng)用實(shí)例[J].水利科技與經(jīng)濟(jì),2010,16（08）:867-869.

[28]王殿志,張慶河,時(shí)鐘.渤海灣風(fēng)浪場(chǎng)的數(shù)值模擬[J].海洋通報(bào),2004,23（05）：10-17.

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

[30]王道龍,華鋒,江志輝.SWAN近岸海浪模式在遼東灣的應(yīng)用[J].海洋科學(xué)進(jìn)展,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]范平易.橋區(qū)流場(chǎng)數(shù)值模擬及建橋?qū)νê降挠绊懷芯?碩士論文.上海交通大學(xué),2010.