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

?

基于現(xiàn)場實驗結(jié)果對傾倒泥沙擴散數(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)用國家海洋局重點實驗室，上海200080）

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

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

Keywords：marine dumping site；dumping experiment；numerical model

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

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

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

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

1.1實驗概述

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

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

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

1.2實驗方法

傾倒實驗中對傾倒形成的懸浮泥沙采用船舶直接采樣分析法。在傾倒區(qū)中心點布設(shè)兩條調(diào)查船，其中一條拋錨固定于傾倒區(qū)中心附近位置，在傾倒船投放后進行定點采樣、觀測；另一條調(diào)查船則跟蹤海面漂浮標(biāo)識物進行懸浮物的擴散與沉降的動態(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為測點水深。

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

1.3水樣采集與處理

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

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

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

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

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

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

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

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

圖5　傾倒試驗跟蹤觀測船路徑圖

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

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

1.5實驗結(jié)果合理性分析

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

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

2.1二維泥沙輸移擴散模型

2.1.1基本方程

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

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

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

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

2.1.2邊界條件和初始條件

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

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

出流時：

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

初始條件由下式確定：

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

2.1.3方程離散求解

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

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

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

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

2.2.1模型源強

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

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

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

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

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

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

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

2.2.3泥沙沉降速度

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

2.2.4沉降幾率

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

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

3.1潮位、潮流驗證

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

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

圖8　大潮HQ07、HQ11潮流驗證

3.2懸浮物濃度過程曲線率定驗證

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

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

圖9　3次實驗定點船測點實測懸浮物濃度與計算懸浮物濃度過程對比

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

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

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

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

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

4　結(jié)論與討論

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

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

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

參考文獻

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)用.水利水運工程學(xué)報，（1）：44-48.

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

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

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

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

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

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

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

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

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

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

徐宏明，2000.疏浚土擴散數(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.海河口潮流泥沙運動數(shù)值模擬及清淤方案研究.水動力學(xué)研究與進展A輯，17（3）：318-326.

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

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

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

張玉芬，蔡思忠，1992.疏浚物在海洋環(huá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

文獻標(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

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

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