基于現(xiàn)場實(shí)驗(yàn)結(jié)果對傾倒泥沙擴(kuò)散數(shù)值模型參數(shù)率定研究

郭玉臣，石志洲，王慶業(yè)，夏永健，高煒，宋微，吳毓儒（1.國家海洋局寧德海洋環(huán)境監(jiān)測中心站，福建　寧德　352100；2.海洋赤潮災(zāi)害立體監(jiān)測技術(shù)與應(yīng)用國家海洋局重點(diǎn)實(shí)驗(yàn)室，上海　200080）

?

基于現(xiàn)場實(shí)驗(yàn)結(jié)果對傾倒泥沙擴(kuò)散數(shù)值模型參數(shù)率定研究

郭玉臣1，2，石志洲1，2，王慶業(yè)1，2，夏永健1，2，高煒1，2，宋微1，2，吳毓儒1，2
（1.國家海洋局寧德海洋環(huán)境監(jiān)測中心站，福建寧德352100；2.海洋赤潮災(zāi)害立體監(jiān)測技術(shù)與應(yīng)用國家海洋局重點(diǎn)實(shí)驗(yàn)室，上海200080）

摘要：傾倒區(qū)傾倒過程泥沙輸運(yùn)擴(kuò)散數(shù)值模型理論較為成熟，但模型的一些相關(guān)參數(shù)給定范圍較大，且一般由實(shí)驗(yàn)室靜水或河口淡水實(shí)驗(yàn)獲取，用于海洋傾倒區(qū)的泥沙模型有一定的局限性，使得模型預(yù)測結(jié)果波動范圍也很大，具有較強(qiáng)的不確定性，模型參數(shù)亟需物理實(shí)驗(yàn)數(shù)據(jù)的率定和驗(yàn)證。2011-2014年在沙埕港臨時性海洋傾倒區(qū)開展了3次傾倒物海洋傾倒實(shí)驗(yàn)，對現(xiàn)場實(shí)驗(yàn)結(jié)果數(shù)據(jù)分析論證后，建立傾倒過程二維泥沙輸運(yùn)擴(kuò)散數(shù)值模型，用實(shí)驗(yàn)數(shù)據(jù)對模型中一些主要參數(shù)進(jìn)行計(jì)算和率定。計(jì)算率定結(jié)果表明：沙埕港傾倒區(qū)傾倒過程初期疏浚泥絮團(tuán)沉降速度可達(dá)3.0 cm/s，細(xì)沙顆粒充分混合后沉降速度減小，一般小于0.5 mm/s；懸沙紊動擴(kuò)散系數(shù)Dx、Dy取45 m2/s時計(jì)算懸浮物濃度增量曲線與實(shí)測結(jié)果吻合效果最佳，建議調(diào)試值域?yàn)?0～50 m2/s。模型參數(shù)率定后計(jì)算結(jié)果的客觀性和準(zhǔn)確性進(jìn)一步提高，可為海洋傾倒區(qū)選劃泥沙擴(kuò)散模型參數(shù)的選取提供一定參考。

關(guān)鍵詞：海洋傾倒區(qū)；傾倒實(shí)驗(yàn)；數(shù)值模型

Keywords：marine dumping site；dumping experiment；numerical model

海洋傾倒處理疏浚物過程中，疏浚物的運(yùn)動通常概化為沉降、擴(kuò)散和輸運(yùn)3個主要過程，許多學(xué)者針對疏浚泥沙運(yùn)動過程的不同階段開展了大量的研究。Koh等（1973）、王錫侯等（1989）、韓康等（1990）、程偉平等（2006）分別采用實(shí)驗(yàn)室水槽實(shí)驗(yàn)和數(shù)值模型方法，重點(diǎn)對疏浚物沉降過程進(jìn)行分析和模擬研究；Cole等（1983）、徐宏明（2000）、陳祖華等（2003）、吳修廣等（2008）、孫毛明等（2009）、李佳等（2011）、張世民等（2014）主要采用數(shù)值模型計(jì)算的方法，研究了傾倒后疏浚物水平運(yùn)動和濃度分布情況，并將模型應(yīng)用于河口、近岸海域及外海等不同水域的傾倒區(qū)進(jìn)行模擬預(yù)測，其研究則側(cè)重于傾倒后的擴(kuò)散和輸移過程。

疏浚物傾倒后運(yùn)動情況的研究除水槽實(shí)驗(yàn)和數(shù)值模擬兩種手段以外，海上現(xiàn)場傾倒實(shí)驗(yàn)方法應(yīng)該更為客觀和準(zhǔn)確。目前，疏浚物傾倒試驗(yàn)觀測方法主要有3類：船舶直接采樣分析法、聲學(xué)或光電拖曳儀器測量法和放射性示蹤法。受儀器設(shè)備和海上作業(yè)難度大等條件影響，關(guān)于海洋疏浚物傾倒實(shí)驗(yàn)的研究仍然較少。Proni等（1975）利用聲學(xué)測控系統(tǒng)進(jìn)行了海洋傾廢懸浮物的研究。張玉芬等（1992）在黃白咀預(yù)選疏浚物傾倒區(qū)進(jìn)行了3次疏浚物傾倒聲學(xué)追蹤試驗(yàn)。張效龍等（2011）采用示蹤沙方法研究長江口拋泥區(qū)泥沙運(yùn)移趨勢。何東海等（2013）在蒼南海域開展了疏浚物傾倒試驗(yàn)，直接采取水樣來觀測懸浮泥沙擴(kuò)散濃度。

本研究采用船舶直接采樣分析法，在福建北部沙埕港臨時性海洋傾倒區(qū)于2011-2014年共開展3次傾倒實(shí)驗(yàn)，對現(xiàn)場實(shí)驗(yàn)觀測數(shù)據(jù)進(jìn)行分析后，建立海洋傾倒泥沙擴(kuò)散數(shù)學(xué)模型進(jìn)行數(shù)值模擬試驗(yàn)，對模型中主要參數(shù)進(jìn)行計(jì)算和率定，提高了海洋傾倒泥沙擴(kuò)散數(shù)學(xué)模型對懸浮物擴(kuò)散規(guī)律和包絡(luò)范圍預(yù)測的準(zhǔn)確性和客觀性。

1　海洋傾倒現(xiàn)場實(shí)驗(yàn)

1.1實(shí)驗(yàn)概述

傾倒實(shí)驗(yàn)共開展3次，時間分別為2011年6月、2012年5月和2014年10月。實(shí)驗(yàn)海域在福建沙埕港臨時性海洋傾倒區(qū)（見圖1），傾倒區(qū)位于東海南部近岸海域，沙埕港東南約20 km處，海底地形變化平緩，水深約為25 m。據(jù)海流監(jiān)測分析結(jié)果，附近海域?yàn)橐?guī)則半日潮流，潮流運(yùn)動形式為旋轉(zhuǎn)流，垂線平均潮流流速一般在30～40 cm/s。

實(shí)驗(yàn)傾倒點(diǎn)設(shè)在臨時性海洋傾倒區(qū)中心位置，實(shí)驗(yàn)選在低平潮后漲潮過程拋泥，然后進(jìn)行跟蹤監(jiān)測。工程淤泥疏浚中采用的是抓斗式挖泥，疏浚泥傾倒方式為開底式泥駁船拋放，單船傾倒量約為600 m3。

圖1　傾倒區(qū)及周邊海域位置示意圖

1.2實(shí)驗(yàn)方法

傾倒實(shí)驗(yàn)中對傾倒形成的懸浮泥沙采用船舶直接采樣分析法。在傾倒區(qū)中心點(diǎn)布設(shè)兩條調(diào)查船，其中一條拋錨固定于傾倒區(qū)中心附近位置，在傾倒船投放后進(jìn)行定點(diǎn)采樣、觀測；另一條調(diào)查船則跟蹤海面漂浮標(biāo)識物進(jìn)行懸浮物的擴(kuò)散與沉降的動態(tài)觀測。

固定調(diào)查船觀測歷時2個小時。監(jiān)測內(nèi)容為懸浮物、海流、風(fēng)向、風(fēng)速同步觀測，監(jiān)測頻率先緊后松，在傾倒后第1個半小時內(nèi)每5 min 1次，第2個半小時內(nèi)觀測頻率變?yōu)?0 min 1次，第2個小時及之后為20 min 1次。懸浮物采樣及海流觀測按表、中、底3層采樣，其中，表層距海面0.2H，中層取現(xiàn)場水深0.6H，底層離底0.8H，H為測點(diǎn)水深。

動態(tài)調(diào)查船主要對疏浚泥擴(kuò)散團(tuán)進(jìn)行跟蹤采樣，在疏浚泥傾倒后，沿現(xiàn)場潮流跟蹤懸浮泥團(tuán)采集水樣，并同步進(jìn)行DGPS定位，定位精度為1 m；采樣頻次與定點(diǎn)調(diào)查船同步，采樣方法與定點(diǎn)船一致。

1.3水樣采集與處理

水體懸浮物濃度使用顛倒采水器采水樣1000mL，懸浮物濃度自下而上分層采樣，抽濾時用蒸餾水洗三次，水樣用0.45 μm微孔濾膜抽濾后烘干測定。量積誤差不超過千分之五，稱量用萬分之一電子天平。

1.4現(xiàn)場實(shí)驗(yàn)結(jié)果分析

沙埕港臨時性海洋傾倒區(qū)3次現(xiàn)場傾倒試驗(yàn)定點(diǎn)調(diào)查船和跟蹤調(diào)查船監(jiān)測結(jié)果分別見圖2-圖4。

圖2　2011年傾倒區(qū)中心HQ07站懸浮物濃度隨時間變化圖

圖3　2012年傾倒區(qū)中心HQ07站懸浮物濃度隨時間變化圖

圖4　2014年傾倒區(qū)中心HQ07站懸浮物濃度隨時間變化圖

從定點(diǎn)觀測數(shù)據(jù)來看，實(shí)驗(yàn)開始第1組采樣數(shù)據(jù)為未傾倒時懸浮物濃度，與實(shí)驗(yàn)前在附近海域取樣調(diào)查的懸浮物濃度比較一致，傾倒完成后的10min內(nèi)測點(diǎn)海水懸浮物濃度逐漸增大，0.8H層濃度增大最為明顯，其中2011年實(shí)驗(yàn)過程中0.8H層濃度最大值為157.7 mg/L，較本底值增量為69.6 mg/L，0.2H層和0.6H層濃度變化相對較小。10 min后各層懸浮物濃度均迅速降低，在20 min后測點(diǎn)懸浮物濃度基本接近本底值狀態(tài)。

從跟蹤觀測數(shù)據(jù)來看，2011年傾倒實(shí)驗(yàn)的跟蹤船舶追蹤路徑為N'向，2012年和2014年兩次實(shí)驗(yàn)中追蹤路徑均為WSW向，3次實(shí)驗(yàn)追蹤路徑見圖5。傾倒后的20 min內(nèi)，采樣海水懸浮物濃度均未出現(xiàn)明顯的特征性變化，20 min后0.6H層和0.8H層海水懸浮物濃度有一個波動增大過程，0.2H層濃度變化不明顯。

圖5　傾倒試驗(yàn)跟蹤觀測船路徑圖

對比分析3次實(shí)驗(yàn)過程和結(jié)果見表1，可以看出：2011年沙埕港臨時傾倒區(qū)傾倒試驗(yàn)采樣本底懸浮物濃度最大，各站垂線平均濃度均值為57.0 mg/L，2014年實(shí)驗(yàn)過程中本底采樣濃度最小。傾倒點(diǎn)附近定點(diǎn)采樣最大濃度出現(xiàn)在2011年實(shí)驗(yàn)過程的0.8H層，最大濃度增量值為69.6 mg/L，另外兩次實(shí)驗(yàn)最大濃度增量分別為39.8 mg/L和44.8 mg/L。根據(jù)傾倒過程潮流分析，傾倒泥沙應(yīng)隨潮流向W'向輸運(yùn)擴(kuò)散，2012年和2014年兩次實(shí)驗(yàn)最總路徑相對合理為WSW向，而2011年跟蹤路徑為N'向，與潮流方向并不吻合，分析原因可能是表層漂浮標(biāo)識物受到風(fēng)的影響發(fā)生偏移。

表1　三次實(shí)驗(yàn)監(jiān)測結(jié)果數(shù)據(jù)表

1.5實(shí)驗(yàn)結(jié)果合理性分析

葉慧明等（2008）對洋山臨時海洋傾倒區(qū)疏浚物傾倒懸浮物動態(tài)跟蹤監(jiān)測結(jié)果顯示，傾倒引起的懸浮物增量范圍為37.5～126.5 mg/L，傾倒1 h后各觀測站點(diǎn)的懸浮物增量效應(yīng)已明顯減弱。本研究3次實(shí)驗(yàn)定點(diǎn)觀測懸浮物濃度增量變化規(guī)律具有較好的重現(xiàn)性，測值較其他研究成果來看基本合理。

2　泥沙擴(kuò)散數(shù)值模型建立

2.1二維泥沙輸移擴(kuò)散模型

2.1.1基本方程

述二維懸沙輸移擴(kuò)散采用以下方程：

式中：s為懸浮物濃度（kg/ m3）；u、v分別為x、y向的水流速度（m/ s）；Dx、Dy分別為x、y向的泥沙紊動擴(kuò)散系數(shù)（m2/s）；H為總水深（m）；Fs為源匯函數(shù)（kg/m2·s），采用挾沙力法其表達(dá)式為：

當(dāng)模擬拋泥時，源匯函數(shù)可表示為如下形式：

式（2）、（3）中：α為泥沙沉降機(jī)率，ω為泥沙沉降速度（m/s），s*為挾沙力（kg/m3），F(xiàn)s′為拋泥源函數(shù)（kg/m2·s）。

2.1.2邊界條件和初始條件

泥沙固邊界法向泥沙通量為0：

泥沙水邊界懸沙按入流和出流情況分別控制，入流時：

出流時：

式中：Γ為水邊界，s*為已知懸浮物濃度（kg/m3），un為法向流速（m/s）。

初始條件由下式確定：

模型泥沙初始條件采用零初始條件，即懸浮物濃度初始為0，計(jì)算懸浮物濃度增量。

2.1.3方程離散求解

泥沙輸運(yùn)擴(kuò)散模型采用均勻正交網(wǎng)格上的有限差分ADI法進(jìn)行離散，并用追趕法對離散方程求解，網(wǎng)格步長為50 m，網(wǎng)格點(diǎn)u，v及ζ由潮流動力模擬結(jié)果插值獲取。

2.2模型參數(shù)選取和率定

模型計(jì)算過程中，根據(jù)不同參數(shù)特點(diǎn)其確定方法也不相同，泥沙干容重、沉降速度等參數(shù)依據(jù)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行相關(guān)計(jì)算確定；疏浚泥懸移質(zhì)比例、沉降幾率等參數(shù)相關(guān)文獻(xiàn)結(jié)論較為明確，取值相對容易確定，故參考相關(guān)文獻(xiàn)經(jīng)驗(yàn)取值；泥沙紊動擴(kuò)散系數(shù)相關(guān)規(guī)范給定值域較大，相關(guān)研究也很少給出明確結(jié)論，故采用數(shù)值試驗(yàn)進(jìn)行率定，模型主要參數(shù)及其確定方法見表2。

圖2　傾倒過程二維泥沙擴(kuò)散數(shù)值模型主要參數(shù)

2.2.1模型源強(qiáng)

拋泥進(jìn)入水體后，大部分將以比單顆粒泥沙沉速大得多的速率沉入海底，形成懸移質(zhì)參與輸運(yùn)擴(kuò)散的比例一般只有5%～10%。單次拋泥源強(qiáng)計(jì)算按如下公式：

式中：sc為拋泥源強(qiáng)（kg/s）；Q為單次拋泥量（m3）；T為傾倒時間（s）；P為疏浚泥產(chǎn)生的懸移質(zhì)比例（%）；γ0為泥沙干容重（kg/m3）；s0為單位源強(qiáng)（kg/ m3·s）；Ai為網(wǎng)格點(diǎn)代表的面積（m2）；Hi為網(wǎng)格點(diǎn)水深（m）。

本研究中參考已有研究成果和實(shí)驗(yàn)監(jiān)測數(shù)據(jù)確定，懸疑質(zhì)比例P取8%，傾倒時間T取5 min，泥沙干容重γ0取1 300 kg/m3。

2.2.2懸沙紊動擴(kuò)散系數(shù)

海岸與河口潮流泥沙模擬技術(shù)規(guī)程（2010）提出懸沙紊動擴(kuò)散系數(shù)可取與相應(yīng)的水流紊動粘性系數(shù)相同數(shù)值，水流紊動粘性系數(shù)宜由試驗(yàn)確定或通過驗(yàn)證計(jì)算確定，可取0～100 m2/s；海洋傾倒區(qū)選劃技術(shù)導(dǎo)則（2009）提出泥沙紊動擴(kuò)散系數(shù)Dx、Dy由試驗(yàn)確定，或采用Elder公式：

式中：C為謝才系數(shù)，C = n-1H1/6，n為曼寧系數(shù)。利用公式計(jì)算的該海域泥沙紊動擴(kuò)散系數(shù)Dx、Dy為35～78 m2/s，數(shù)值試驗(yàn)中利用實(shí)驗(yàn)數(shù)據(jù)對Dx、Dy進(jìn)行調(diào)試率定，試驗(yàn)過程不同取值擬合曲線見圖6，最終模型調(diào)試率定Dx、Dy取45 m2/s時擬合效果最佳，建議調(diào)試值域30～50 m2/s。

圖6　懸沙紊動擴(kuò)散系數(shù)不同取值擬合曲線

2.2.3泥沙沉降速度

泥沙沉降速度通常根據(jù)武水公式計(jì)算，本研究疏浚物樣品分析結(jié)果為粘土質(zhì)粉沙，粒徑0.004～0.032 mm占比最大，按公式計(jì)算沉降速度變化范圍很大8×10-3～0.5 mm/s。而根據(jù)現(xiàn)場試驗(yàn)監(jiān)測結(jié)果顯示，拋泥后海水中絮團(tuán)在10 min前后即沉降到0.8H層或更深的海底，據(jù)此推算，傾倒后絮凝疏浚泥沙沉速應(yīng)≥3.0 cm/s。分析原因主要是傾倒開始的一段時間內(nèi)細(xì)顆粒泥沙在海水環(huán)境會發(fā)生絮凝，絮凝泥沙塊沉速增大。據(jù)相關(guān)研究：受潮流動力影響絮凝泥沙沉速可達(dá)絮凝前分散泥沙沉速的10余倍，乃至100余倍（嚴(yán)鏡海，1988）。通常未發(fā)生絮凝的懸浮泥沙沉降速度ω取值在0.3～0.5mm/s，本研究泥沙模型主要考慮懸沙輸運(yùn)和擴(kuò)散過程，故計(jì)算調(diào)試后取泥沙沉速0.35 mm/s。

2.2.4沉降幾率

沉降系數(shù)是反映懸沙運(yùn)動的一個綜合系數(shù)，表示懸沙沉降的幾率。關(guān)于的取值目前沒有定論，丁興平（2003）在長江河口波流共同作用下的全沙數(shù)值模擬中取α=0.75，張華慶等（2002）在海河口潮流泥沙運(yùn)動數(shù)值模擬及清淤積方案研究中取α= 0.5。本研究考慮疏浚泥沙在傾倒時顆粒間相互碰撞比較劇烈，并易受上層泥沙的重力影響，增大了泥沙沉降幾率，故取α=0.9。

3　數(shù)值模型率定驗(yàn)證結(jié)果分析

3.1潮位、潮流驗(yàn)證

準(zhǔn)確、合理的潮流動力數(shù)值模擬結(jié)果是泥沙輸運(yùn)擴(kuò)散模擬的基本條件，采用2014年8月大、小潮期附近海域2個長期潮位站和4個周日連續(xù)海流測站的實(shí)測資料對模型進(jìn)行了驗(yàn)證，結(jié)果表明：各潮位驗(yàn)證點(diǎn)潮位計(jì)算值與實(shí)測值的誤差很小，相對誤差均小于10%；計(jì)算潮流流速、流向曲線與實(shí)測值吻合較好，能夠很好的反映各站潮流變化規(guī)律。三沙站（ss）潮位驗(yàn)證結(jié)果見圖7，HQ07、HQ11測站潮流驗(yàn)證結(jié)果見圖8。

圖7　三沙站大潮潮位驗(yàn)證結(jié)果

圖8　大潮HQ07、HQ11潮流驗(yàn)證

3.2懸浮物濃度過程曲線率定驗(yàn)證

基于對沙埕港臨時性海洋傾倒區(qū)3次傾倒實(shí)驗(yàn)結(jié)果的分析，對泥沙輸運(yùn)擴(kuò)散模型的率定和驗(yàn)證主要采用定點(diǎn)船的觀測數(shù)據(jù)。

二維泥沙擴(kuò)散模型計(jì)算結(jié)果給出的泥沙變化是經(jīng)過垂向積分的垂線平均懸浮物濃度，實(shí)驗(yàn)過程監(jiān)測的0.2H、0.6H、0.8H層懸浮物濃度經(jīng)加權(quán)平均法計(jì)算得垂線平均懸浮物濃度后，繪制計(jì)算值與實(shí)測值比較結(jié)果見圖9。其中，2014年實(shí)驗(yàn)數(shù)據(jù)用于泥沙模型主要參數(shù)的率定，率定對比見圖9（a）；模型調(diào)試后對2011年和2012年兩次實(shí)驗(yàn)過程進(jìn)行后報驗(yàn)證，驗(yàn)證結(jié)果見圖9（b、c）分析率定后數(shù)學(xué)模擬驗(yàn)證曲線對比結(jié)果可知，模型結(jié)果能較好的反映傾倒區(qū)拋泥后懸浮物濃度變化過程，除個別時刻量值差異較大外，曲線吻合較好，尤其在泥沙濃度時間相位擬合較為準(zhǔn)確。

圖9　3次實(shí)驗(yàn)定點(diǎn)船測點(diǎn)實(shí)測懸浮物濃度與計(jì)算懸浮物濃度過程對比

3.32014年傾倒試驗(yàn)懸沙擴(kuò)散過程模擬結(jié)果分析

對模型模擬的2014年沙埕港臨時海洋傾倒區(qū)傾倒實(shí)驗(yàn)過程進(jìn)行分析：拋泥后形成的懸浮物高濃度區(qū)僅出現(xiàn)在拋泥點(diǎn)附近50 m半徑范圍內(nèi)，且在5 min后消失，懸浮物云團(tuán)隨潮流向偏西方向移動，濃度降低較快；30 min后大于20 mg/L的濃度增量區(qū)已經(jīng)消失，10mg/L濃度增量區(qū)域半徑約為200m，區(qū)域中心位于傾倒點(diǎn)偏西向500 m附近；60 min后5 mg/L懸浮物增量區(qū)域也逐漸消失，后期較低濃度的云團(tuán)消散速度變慢。傾倒后30 min和60 min時刻垂線平均懸浮物濃度增量等值線見圖10。

圖10　傾倒過程垂線平均懸浮物濃度增量等值線分布圖（單位mg/L）

根據(jù)傾倒過程數(shù)值模擬計(jì)算結(jié)果，繪制各類水質(zhì)標(biāo)準(zhǔn)濃度值包絡(luò)范圍見圖11，可以看出，100mg/L（超三類水質(zhì)標(biāo)準(zhǔn)）及更高濃度增量出現(xiàn)在傾倒點(diǎn)很小范圍內(nèi)，50 mg/L濃度增量包絡(luò)范圍最大直徑約為80 m，10 mg/L（超一、二類水質(zhì)標(biāo)準(zhǔn)）濃度增量包絡(luò)范圍相對較大，沿潮流方向在傾倒點(diǎn)西側(cè)延伸分布，最遠(yuǎn)離傾倒點(diǎn)約800 m。

圖11　傾倒過程懸浮物濃度增量包絡(luò)范圍圖

4　結(jié)論與討論

（1）沙埕港臨時性海洋傾倒實(shí)驗(yàn)結(jié)果顯示，傾倒完成后的10min內(nèi)測點(diǎn)海水懸浮物濃度逐漸增大，0.8H層濃度增大最為明顯，其中2011年實(shí)驗(yàn)過程中0.8H層濃度最大值為157.7 mg/L，較本底值增量為69.6 mg/L，2012年和2014年2次實(shí)驗(yàn)最大濃度增量分別為39.8 mg/L和44.8 mg/L。傾倒點(diǎn)附近海水懸浮物濃度10 min后開始迅速降低，在20 min后測點(diǎn)懸浮物濃度基本接近本底值狀態(tài)。

（2）利用現(xiàn)場傾倒實(shí)驗(yàn)結(jié)果對傾倒泥沙擴(kuò)散數(shù)值模型進(jìn)行計(jì)算和率定，該海域泥沙紊動擴(kuò)散系數(shù)Dx、Dy取值45m2/s時模擬計(jì)算效果最佳，建議值域?yàn)?0～50 m2/s，傾倒后絮凝疏浚泥沙沉速可達(dá)3.0 cm/s以上，參與懸浮的泥沙沉降速度依據(jù)疏浚物粒徑特征計(jì)算并調(diào)試后取0.35 mm/s。率定調(diào)試后數(shù)值模型對拋泥過程模擬結(jié)果更具客觀性和準(zhǔn)確性，可為傾倒區(qū)選劃泥沙擴(kuò)散數(shù)值模型參數(shù)選取提供一定參考價值。

（3）現(xiàn)場傾倒實(shí)驗(yàn)過程中懸浮物濃度監(jiān)測方法為現(xiàn)場采樣法，存在對實(shí)驗(yàn)過程數(shù)據(jù)的獲取不夠全面完整，無法非?？陀^、精確刻畫反映泥沙運(yùn)動擴(kuò)散情況的局限性，在今后的相關(guān)研究中引進(jìn)更先進(jìn)儀器，實(shí)現(xiàn)對懸浮泥沙運(yùn)動擴(kuò)散場全方位動態(tài)監(jiān)測，可進(jìn)一步提高數(shù)值模型參數(shù)率定研究成果的可靠性和科學(xué)性。

參考文獻(xiàn)

Cole P，Miles G V，1983. Two-dimensional model of mud transport. J Hydr Eng，109（1）：1-12.

Proni J R，Rona D C，Lauter C A，et al，1975. Acoustic observations of suspended particulate matter in the ocean. Nature，254（5499），413-415.

Koh R C Y，Chang Y C，1973. Mathematical model for barged ocean disposal of wastes. EPA Grant No.16070 FBY.

陳祖華，王光謙，王志石，2003.河口拋泥數(shù)學(xué)模型及應(yīng)用.水利水運(yùn)工程學(xué)報，（1）：44-48.

程偉平，毛根海，章軍軍，2006.拋泥過程中泥沙云團(tuán)運(yùn)動的雙流體大渦模擬與試驗(yàn)研究.水力發(fā)電學(xué)報，25（3）：110-115.

丁平興，胡克林，孔亞珍，等，2003.長江河口波-流共同作用下的全沙數(shù)值模擬.海洋學(xué)報，（5）：113-124.

韓康，黃河寧，王錫侯，1990.疏浚物海洋傾倒沉降過程的實(shí)驗(yàn)和數(shù)值模擬研究.海洋環(huán)境科學(xué)，9（3）：1-5.

何東海，何琴燕，吳光榮，等，2013.蒼南海域疏浚物傾倒懸浮物擴(kuò)散特征現(xiàn)場試驗(yàn)分析.海洋工程，31（3）：101-106.

交通部天津水運(yùn)工程科學(xué)研究所，2010.海岸與河口潮流泥沙模擬技術(shù)規(guī)程.人民交通出版社.

李佳，姚炎明，孫志林，等，2011.大型海洋傾倒區(qū)懸浮物遷移擴(kuò)散的數(shù)值模擬.浙江大學(xué)學(xué)報（工學(xué)版），45（7）：1319-1328.

全國海洋標(biāo)準(zhǔn)化技術(shù)委員會，2009.海洋傾倒區(qū)選劃技術(shù)導(dǎo)則.

孫毛明，吳修廣，倪勇強(qiáng)，等，2009.疏浚土傾倒后懸浮泥沙擴(kuò)散輸移的數(shù)值模擬.海洋學(xué)研究，27（3）：22-30.

王錫侯，黃河寧，韓康，1989.傾廢物云團(tuán)沉降過程的數(shù)學(xué)模型.海洋學(xué)報，11（1）：63-69.

吳修廣，楊潤偉，史英標(biāo)，等，2008.外海拋泥過程中懸浮泥沙擴(kuò)散輸移的數(shù)值模擬.海洋環(huán)境科學(xué)，27（3）：236-241.

徐宏明，2000.疏浚土擴(kuò)散數(shù)學(xué)模型及其應(yīng)用.海洋環(huán)境科學(xué)，19 （2）：34-37.

嚴(yán)鏡海，1988.潮汐水流中細(xì)顆粒泥沙絮凝沉降的初步探討.泥沙研究，（4）：10-22.

葉慧明，2008.洋山臨時海洋傾倒區(qū)使用對所在海域及其周邊環(huán)境的影響分析.海洋環(huán)境科學(xué)，27（S1）：29-33.

張華慶，李華國，岳翠平，2002.海河口潮流泥沙運(yùn)動數(shù)值模擬及清淤方案研究.水動力學(xué)研究與進(jìn)展A輯，17（3）：318-326.

張世民，陳德文，鄧兆青，2014.東碇傾倒區(qū)疏浚泥懸沙輸移和海床沖淤數(shù)值模擬.海洋通報，33（5）：541-551.

張效龍，劉敦武，王慧艷，2010.采用示蹤砂方法研究長江口北槽下航道南側(cè)拋泥區(qū)泥沙運(yùn)移趨勢.海洋通報，29（3）：277-282.

張效龍，劉敦武，王慧艷，2010.采用示蹤砂方法研究長江口北槽下航道南側(cè)拋泥區(qū)泥沙運(yùn)移趨勢.海洋通報，29（3）：277-282.

張玉芬，蔡思忠，1992.疏浚物在海洋環(huán)境中的物理沉降、擴(kuò)散和輸運(yùn).海洋通報，11（6）：40-76.

（本文編輯：袁澤軼）

Parameter determination of the numerical model for dumping sediment transport based on the field experimental results

GUO Yu-chen1，2，SHI Zhi-zhou1，2，WANG Qing-ye1，2，XIA Yong-jian1，2，GAO Wei1，2，SONG Wei1，2，WU Yu-ru1，2
（1. Ningde Marine Environmental Monitoring Center，SOA，Ningde 352100，China；2. Key Laboratory of Integrated Marine Monitoring and Applied Technologies for Harmful Algal Blooms，SOA，Shanghai 200080，China）

Abstract：Sediment transport numerical model theory for ocean dumping site is relatively mature，but the range of some parameters in the model is given too wide，and some parameter values are obtained by sediment experiments in the laboratory hydrostatic or estuarine fresh water in general. Using them for the marine dumping area has some limitations，which makes the model prediction results not accurate and vary considerably，so some model parameters urgently need the determination and verification by field experiment data. We did three ocean dumping experiments at the temporary ocean dumping site of Shacheng Harbour during 2011-2014，and based on the field experiment data，we established an two-dimensional sediment transport and diffusion numerical model for the ocean dumping，using observation data of experiments to determine and calculate parameters for the model. The results show that: in the temporary ocean dumping site of Shacheng Harbour，the sedimentation rate of floccule dredged mud may be bigger than 3.0 cm/s at the beginning of the process of dumping，and the velocity of mix sand particles decreases，and generally is less than 0.5 mm/s. The calculate sediment concentration increment curve agreeswith the measured results best while the sediment turbulent diffusion coefficient DxDy values 45 m2/s，so we suggest that the Dx Dy's debugging domain should be of 30～50 m2/s. After the parameter determination，the model prediction results are more objective and accurate，which can provide some reference for the parameter determination of sediment diffusion model in selecting marine dumping sites.

中圖分類號：P731.1

文獻(xiàn)標(biāo)識碼：A

文章編號：1001-6932（2016）02-0170-08

Doi：10.11840/j.issn.1001-6392.2016.02.007

收稿日期：2015-03-12；

修訂日期：2015-07-02

基金項(xiàng)目：國家海洋局東海分局青年科技基金（201206）。

作者簡介：郭玉臣（1981-），男，碩士，工程師，主要從事海洋環(huán)境預(yù)報和數(shù)值模擬研究。電子郵箱：gyc@eastsea.gov.cn。