黃 璐 張家年 吳昊雨 黃河清

(安徽工業(yè)大學(xué)環(huán)境流體研究所 安徽馬鞍山 243002)

0 引言

濁流是在諸如地震、海底滑坡等觸發(fā)機(jī)制下誘發(fā)的由水流湍流支撐的載有顆粒的、重力驅(qū)動(dòng)的底流，是將大陸架上的碎屑沉積物遷移到深海的主要作用力之一［1，2］。海底濁流具有極大的破壞作用，可造成海底電纜的折斷和一些海底設(shè)備及儀表的破壞［3］。濁流傳播過程中會(huì)引發(fā)海底地形演變，形成溝渠、堤壩和海底扇沉積系統(tǒng)。同時(shí)海底濁流在運(yùn)動(dòng)過程中也會(huì)發(fā)生沉積，可能會(huì)形成含有油藏的濁積砂體［2，4］。近20年來，國內(nèi)外和古代濁流沉積相關(guān)的深水油氣探測(cè)已經(jīng)取得了很大的進(jìn)展，在南美、墨西哥灣、北海以及澳大利亞西北大陸架等海域的濁積巖中相繼發(fā)現(xiàn)了許多大型海底油氣田，勘探領(lǐng)域已擴(kuò)展到以西非、巴西坎波斯盆地和墨西哥灣這3大熱點(diǎn)地區(qū)為代表的水深3 000 m 以上的深海區(qū)［5，6］。

由于對(duì)濁流及其沉積的研究具有重要的經(jīng)濟(jì)、環(huán)境意義，國內(nèi)外對(duì)其都十分重視。21世紀(jì)以來，濁流的三維研究逐漸增多。Peakall等［7］用一系列物理實(shí)驗(yàn)研究了具有恒定寬度的彎曲海底峽谷的沉積模式，發(fā)現(xiàn)在相同的地形條件下，海底彎曲峽谷中內(nèi)岸的點(diǎn)壩反而比河流中的點(diǎn)壩所處的位置更遠(yuǎn)。Huang等［8］分析了橫剖面為梯形的海底直峽谷中濁流的流動(dòng)形態(tài)及沉積，通過數(shù)值模擬再現(xiàn)了海底濁流的自我建堤過程，并發(fā)現(xiàn)濁流沉積物平均粒徑云圖和沉積厚度云圖在形狀上相似，但在峽谷方向上略有延長(zhǎng)。Amos等［9］通過在沉積物底床上釋放鹽水流的方式產(chǎn)生濁流，進(jìn)行了不同彎曲度下的濁流流動(dòng)實(shí)驗(yàn)，發(fā)現(xiàn)彎曲度足夠大時(shí)，上游的溢流會(huì)重新進(jìn)入下游的渠道中，產(chǎn)生異常的流動(dòng)和沉積模式。Straub等［10］通過一系列實(shí)驗(yàn)，模擬了彎曲海底峽谷中濁流的流動(dòng)和沉積，提出了一個(gè)控制濁流在峽谷內(nèi)流動(dòng)的速度最大值公式。盡管這一系列的實(shí)驗(yàn)研究對(duì)加深海底濁流的理解有幫助，但實(shí)驗(yàn)室實(shí)驗(yàn)由于規(guī)模和實(shí)驗(yàn)條件的限制，模型常常過于簡(jiǎn)化，一般模擬渠道的長(zhǎng)度都不超過10 m;同時(shí)實(shí)驗(yàn)?zāi)M耗資大，不易模擬大尺度的濁流。數(shù)值計(jì)算模擬海底濁流的三維流動(dòng)形態(tài)及其沉積，有利于縮小模型與自然界真實(shí)濁流的差距，加深我們對(duì)自然界濁流的理解，同時(shí)具有相對(duì)于實(shí)驗(yàn)?zāi)M的經(jīng)濟(jì)性和簡(jiǎn)便性。

本文采用一經(jīng)多項(xiàng)實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證的計(jì)算模型對(duì)深海中橫剖面為圓弧形的彎曲海底峽谷中的濁流進(jìn)行模擬并分析濁流在峽谷內(nèi)部和漫灘上的流動(dòng)、密度變化以及二次流和一些沉積特征。這些特征對(duì)于我們根據(jù)對(duì)濁流的沉積觀察推測(cè)其形成環(huán)境及油氣儲(chǔ)層的調(diào)查等方面有一定的幫助作用。

1 模型概述

如圖1所示，濁流流經(jīng)的海底峽谷地形多為彎曲狀，和地表的河流有許多相似之處(Wynne等，2007［11］)。

圖1 典型的海底峽谷示意圖(據(jù)Wynn等，2007［11］)Fig.1 Schematic diagram of typical submarine canyon(after Wynn，et al.，2007)

根據(jù)Straub等、Langbein和Leopold及Pirmez的對(duì)地表及水下渠道的觀察及實(shí)驗(yàn)研究［10，12，13］，本研究采用可以很好地模擬典型地面和水下渠道的理想正弦曲線［10，12，13］作為模型彎道的峽谷中心線，該曲線由兩個(gè)波長(zhǎng)為50 m的彎道組成，彎曲度為1.028，底坡坡度設(shè)置為0.001 5°。在彎道上游設(shè)有70 m的直道，以使?jié)崃髂芤越咏匀坏臓顟B(tài)進(jìn)入彎道。兩個(gè)彎道后的下游與100 m左右長(zhǎng)的直道相連。計(jì)算模型與原型的長(zhǎng)度比尺為1∶100，如圖2所示，整個(gè)模型模擬的區(qū)域長(zhǎng)約270 m，寬為170 m，峽谷的兩邊為漫灘;峽谷橫剖面為曲率半徑為10 m的圓弧，圓弧頂寬為10 m。計(jì)算模型可看作模擬一個(gè)在小坡度的彎曲峽谷內(nèi)流動(dòng)的海底濁流，計(jì)算區(qū)域沿x(流向)、y(垂向)、z(橫向)方向的網(wǎng)格數(shù)分別為 100，40，120，共含4.8×105個(gè)計(jì)算網(wǎng)格點(diǎn)。

深海中濁流流速小到不足1 m/s，大到十幾米/秒［1］。Prior等［14］對(duì)位于英國哥倫比亞的比特島入口一個(gè)細(xì)長(zhǎng)的渠道扇系統(tǒng)內(nèi)長(zhǎng)達(dá)1年的檢測(cè)得出的數(shù)據(jù)表明，濁流最大速度為 3.35 m/s，Xu等［15］于2004年觀測(cè)到濁流的速度最大值為1.9 m/s。從1929年紐芬蘭大淺灘濁流事件對(duì)濁流的破壞估計(jì)的頭部最大速度約為19 m/s［3］，并且觀察到的濁流中的顆粒物體積比濃度一般都較低，為0.1% ～7%［2］。在數(shù)值模擬計(jì)算中，應(yīng)用密度弗雷德數(shù)相似準(zhǔn)則，當(dāng)長(zhǎng)度比為1∶100時(shí)，設(shè)入口處的流速為0.26 m/s，入流含比重為2.65平均粒徑為25 μm沉積物的體積比濃度為2%，相當(dāng)于原型濁流的平均流速為2.6 m/s，在頂寬1 km，長(zhǎng)約27 km彎曲峽谷內(nèi)以同等的沉積物體積比流動(dòng)，和上面提到的自然界中觀察到的濁流速度接近。計(jì)算模擬的一些主要參數(shù)見表1。模擬濁流流動(dòng)時(shí)間為600 s，入流面面積為6.40 m2，占峽谷橫剖面面積的70.62%，位于峽谷橫剖面的底部。

2 計(jì)算模型簡(jiǎn)述

含低密度沉積物的海底濁流可將其視為不可壓縮的連續(xù)體，并應(yīng)用鮑辛尼斯克假設(shè)(Boussinesq assumption)對(duì)其進(jìn)行模擬［8］。在上述條件下模型的雷諾平均質(zhì)量和動(dòng)量方程如下:

圖2 模型設(shè)置圖(a).整個(gè)峽谷的俯視圖及橫剖面分布;(b).峽谷橫剖面圖;(c).峽谷處(b)圖中的方框部分的放大圖Fig.2 Schematic diagram of the model(a).Planform view of the canyon and the location of the cross-sections，(b).a transect of the canyon，(c).Close-up view of the rectangular part in(b)

表1 計(jì)算模擬主要參數(shù)Table 1 Main parameters of numerical simulation

其中，ui，uj為 x1，xj方向上的雷諾平均速度，P為雷諾平均壓力，ρ，ρs分別為環(huán)境流體和沉積物密度，c為懸浮沉積物的雷諾平均體積比濃度，gi為重力加速度。υi為渦黏度系數(shù)，通過求解經(jīng)浮力項(xiàng)修正的如下湍流k-ε模型，得

式中k為湍流動(dòng)能，ε為湍流動(dòng)能耗散率。沉積物質(zhì)量守恒方程為

式中，υs為沉積物顆粒的下降速度，δj2為克羅內(nèi)克符，Sct為施密特系數(shù)。在模型中假設(shè)懸浮沉積物隨流體一同流動(dòng)，但在垂直方向有一下降速度，此速度根據(jù)由實(shí)驗(yàn)確定的經(jīng)驗(yàn)方程來計(jì)算。

假定到達(dá)航道入口的船舶數(shù)量為一個(gè)計(jì)數(shù)過程，每艘船舶到達(dá)時(shí)間是相互獨(dú)立的，而且在充分短的時(shí)間區(qū)間上最多只到達(dá)一艘船舶。在數(shù)學(xué)上，這一過程一般采用泊松分布來描述，其中單位時(shí)間到達(dá)的船舶數(shù)量用參數(shù)λ表示。

本模型計(jì)算模擬海底濁流的幾乎所有重要的物理過程，在河床的底部邊界采用Exner方程模擬由于沉積物的沉積、夾帶和推移質(zhì)搬運(yùn)導(dǎo)致的底床平面的改變［16］。沉積物沉積速度根據(jù)由方程(4)計(jì)算的河床底部沉積物體積比濃度與沉積物沉降速度相乘得出，沉積物夾帶速度和推移質(zhì)搬運(yùn)分別根據(jù)Smith和Mclean［17］以及 van Rijn［18］實(shí)驗(yàn)得出的經(jīng)驗(yàn)公式計(jì)算，詳細(xì)請(qǐng)參見所給出的參考文獻(xiàn)。

應(yīng)用了基于非正交網(wǎng)格的有限體積法來對(duì)上述模型方程進(jìn)行數(shù)值求解，具體解法請(qǐng)參見參考文獻(xiàn)8及16。此模型經(jīng)實(shí)驗(yàn)的連續(xù)入流的水下鹽水密度流、含單粒徑及多粒徑濁流、突然釋放型的重力流以及三維直渠道內(nèi)的濁流流動(dòng)及沉積驗(yàn)證，均取得了理想的模擬結(jié)果［8，16，19，20］。

3 模擬結(jié)果及討論

應(yīng)用上述數(shù)值計(jì)算模型模擬了參數(shù)如表1所示的濁流在如圖2所示水下彎曲渠道內(nèi)的流動(dòng)，觀察到的濁流在海底峽谷水平面上的密度分布、峽谷橫剖面的直道及彎道頂點(diǎn)中心沿流向速度分布、直道及彎道頂點(diǎn)處的密度分布、彎道頂點(diǎn)處的二次流特征等如下。

3.1 海底峽谷水平面上的濁流密度分布

濁流的驅(qū)動(dòng)力來自于與環(huán)境流體的密度差，濁流在流動(dòng)過程中會(huì)逐漸形成密度分層，還會(huì)伴隨著顆粒沉積及夾帶，形成一些類似于河流點(diǎn)壩的沉積形態(tài)。圖3顯示了峽谷頂部上方0.005 m處的密度分布云圖，從圖中可以看出:

(1)在彎道前的直道尾部，流體已經(jīng)有明顯的外溢，因?yàn)闈崃髟谶\(yùn)動(dòng)過程中不斷夾帶上層周圍流體，使?jié)崃骱穸炔粩嘣黾?，?dāng)濁流厚度超過峽谷深度時(shí)就會(huì)外溢，到達(dá)漫灘區(qū)域，其中的懸浮沉積物將在漫灘上沉積;

(2)在彎道頂點(diǎn)外岸處，濁流發(fā)生了明顯的流動(dòng)剝離(Kassem和Irman［21］將流體由于離心力導(dǎo)致的在彎道處溢出峽谷外稱為“剝離(stripping)”，將峽谷任意處由于與周圍流體的夾帶帶來的流體厚度的增加所導(dǎo)致的流體溢出峽谷外稱為“溢流(spilling)”，在這里我們也采用此種說法)，使?jié)崃鳟a(chǎn)生大量的動(dòng)量和密度損失，從而使峽谷內(nèi)的沉積型濁流的速度和密度減小;

(3)彎道外岸和內(nèi)岸產(chǎn)生了明顯的密度差，彎道內(nèi)岸基本保持了與環(huán)境流體相近的密度，而在彎道外岸，濁流密度明顯增大。每個(gè)彎道頂點(diǎn)處，漫灘區(qū)域的流體密度在靠近彎道外岸處達(dá)到最大，其中無量綱過量密度最大可達(dá)入流的37.5%，然后隨著與外岸距離的增加而減小;

(4)隨著沿流向距離的增加，濁流剝離產(chǎn)生的密度損失越來越小，由區(qū)域B至E，最大無量綱過量密度從入流的37.5%逐漸減小為18.8%。區(qū)域A的特殊性是因?yàn)閯傔M(jìn)入彎道時(shí)流體所受的離心力較小，導(dǎo)致剝離較少。

這些模擬結(jié)果和美國的Straub等［10］以及英國的Kane等［22］的實(shí)驗(yàn)觀察基本是一致的。

圖3 峽谷頂部上方0.005 m處的密度(kg/m3)云圖Fig.3 Density(kg/m3)contour of the current at 0.005 m above the canyon top

3.2 峽谷橫剖面直道及彎道頂點(diǎn)中心沿流向速度分布

圖4顯示了峽谷直道末端x=70 m處和三個(gè)彎道頂點(diǎn)，即x=95，120，145 m處的橫剖面中心線的縱向速度結(jié)構(gòu)(對(duì)應(yīng)圖2a橫剖面2～5)，橫坐標(biāo)表示濁流水平速度，縱坐標(biāo)表示距峽谷底部的高度。直道上濁流的縱向速度結(jié)構(gòu)與實(shí)驗(yàn)室測(cè)得的直渠道內(nèi)濁流基本一致，即在距底部不遠(yuǎn)處，存在一速度極大值，其下為壁面邊界層流，其上為自由剪切層流［23］。但在三個(gè)彎道頂點(diǎn)處，最大速度上方部分，速度結(jié)構(gòu)曲線平滑，最大值下方，流向速度曲線出現(xiàn)一明顯的可能是由于二次流引起的波動(dòng)，與Kane等［22］的實(shí)驗(yàn)觀測(cè)值一致;并且有越往下游的彎道頂點(diǎn)處沿流向最大速度越小，其下波動(dòng)位置越往上移的特征，反映了沉積型濁流漸漸減弱的特性。第一個(gè)彎道頂點(diǎn)處中心線的最大流速較前面直道上的要大，是因?yàn)闈崃鞯倪\(yùn)動(dòng)受到離心力影響，流體整體朝外岸傾斜，外彎部分的流體速度密度都比內(nèi)彎大，因此峽谷中心線的速度相對(duì)于直道上有所增大。

3.3 彎道頂點(diǎn)處的濁流密度分布

濁流與環(huán)境流體的密度差相對(duì)較小，所以在彎道處濁流在外岸處的升高會(huì)比明渠流中的水要高得多。Keevil等［24］于2006年用實(shí)驗(yàn)說明了控制海底和河流峽谷的演化和遷移的流動(dòng)過程可能存在極大的差異，并指出海底峽谷的形成和演化也可能顯示出與具有相似地形的河流峽谷不同的特征。圖5a～d所示的圖2中1，3，4，5橫剖面的密度云圖(圖中I表示內(nèi)彎，O表示外彎)有如下特征:

圖4 三個(gè)彎道轉(zhuǎn)換處的橫剖面中心線處的縱向速度結(jié)構(gòu)Fig.4 Longitudinal velocity profile at the canyon centerline of each apex of the three bends

(1)直道內(nèi)的濁流基本上是對(duì)稱的(圖5a);

(2)濁流在彎道處峽谷上部顯示出向著外岸傾斜的特征，在和岸高相同的水平高度處彎道外彎的密度明顯比內(nèi)彎的大(圖5b，c，d)，這與圖3所顯示的密度分布特征是一致的;

(3)在峽谷的底部，濁流密度顯示了與上部不同的分布特征，靠近彎道內(nèi)彎處底部的密度要高一些(圖5b，d)，這是由于受到下節(jié)所分析的彎道中的二次流影響的原因;

(4)流體自上而下出現(xiàn)密度分層，密度層具有越往底部密度增加越快的特征，這是由于受到顆粒沉積的影響。密度層上層等高線厚度逐漸增加，這是因?yàn)闈崃髟谶\(yùn)動(dòng)過程中不斷夾帶周圍流體，與周圍流體進(jìn)行交換，使本身密度減小。

圖5 直道1(a)及彎道橫剖面3-5(b，c，d)的密度(kg/m3)云圖Fig.5 Density(kg/m3)contours of cross-section 1(a)and cross-sections 3 to 5(b，c，d)

3.4 彎道頂點(diǎn)處的二次流特征

濁流在彎曲峽谷中的二次流非常明顯，最近有一系列關(guān)于水下重力流在彎曲峽谷內(nèi)的二次流的實(shí)驗(yàn)、計(jì)算模擬及理論分析論文，討論彎道頂點(diǎn)處的濁流流動(dòng)方向究竟和河流一樣還是相反［25～28］，Abad 等［27］根據(jù)實(shí)驗(yàn)和理論分析得出彎曲矩形渠道內(nèi)的鹽水重力流二次流的方向會(huì)因?yàn)榱鲬B(tài)的不同而不同，Serchi等［28］則指出二次流的方向會(huì)因?yàn)槠露鹊淖兓l(fā)生改變。圖6a和圖6b表示第一個(gè)彎道頂點(diǎn)(x=95 m)處峽谷橫剖面以及彎道外岸峽谷和漫灘交界處的速度矢量圖(I表示內(nèi)岸，O表示外岸)，由圖中可見:

(1)對(duì)于我們所模擬的小坡度、低彎曲度的圓弧形峽谷來說，彎道頂點(diǎn)底部的二次流是和河流類似的順時(shí)針方向的(圖6a)，而在 Serchi等［28］的模擬結(jié)果中也曾提到在低坡度下，彎曲渠道中的二次流方向是與河流中的一致的，Kassem和Irman［21］所模擬的小坡度矩形水下彎曲渠道中的二次流方向也是如此，和我們的模擬結(jié)果相符合;Corney等［25，26］論及的水下濁流在彎道頂點(diǎn)處是二次流和河流相反的觀點(diǎn)可能并不具一般性;

(2)圖6b顯示:彎道外岸，濁流在峽谷與漫灘交界處發(fā)生分流，一部分沿著峽谷邊壁向下流動(dòng)，另一部分則流向漫灘;

(3)正是由于底部二次流動(dòng)的方向指向內(nèi)岸，從而在彎道頂點(diǎn)處使底部高密度濁流向內(nèi)彎傾斜，和頂部濁流的向外岸傾斜不同(參見圖5b，d)。

圖6 第一個(gè)彎道頂點(diǎn)(x=95 m)橫剖面(a)及其外岸局部(b)的速度矢量圖Fig.6 Velocity vector of the first bend apex cross-section(a)and a close-up of its outer bank region(b)

總的說來，根據(jù)我們模擬的結(jié)果及綜合其它的實(shí)驗(yàn)和模擬結(jié)果［25～28］來看，海底濁流在彎曲峽谷內(nèi)的二次流的方向及形態(tài)不僅和流態(tài)而且和峽谷的形狀及寬深比等都有一定的關(guān)系，不可一概而論。

3.5 基本沉積特征

圖7 t=600 s時(shí)的底床上顆粒物沉積厚度Fig.7 Sedimentary thickness over the bed at t=600 s

濁流中的懸浮沉積物在運(yùn)動(dòng)遷移過程中會(huì)發(fā)生沉積，同時(shí)快速流動(dòng)的濁流會(huì)侵蝕底床，重塑底床形態(tài)。圖7顯示了所模擬濁流在t=600 s時(shí)的沉積厚度分布。峽谷底部中心線附近兩邊有較多的沉積，由于邊壁較陡，峽谷邊壁上則出現(xiàn)小規(guī)模的侵蝕，這將造成峽谷橫剖面形狀的輕微改變;進(jìn)入彎道前，直道上的沉積厚度分布基本是對(duì)稱的;峽谷外漫灘部分由于濁流有外溢，也產(chǎn)生了少量的沉積，并隨著與峽谷中心線距離的增加而遞減;進(jìn)入彎道后，在彎道內(nèi)部，每個(gè)彎道頂點(diǎn)外岸上游，即內(nèi)岸下游處有明顯的由于迴流導(dǎo)致的沉積，而在外岸下游處則發(fā)生侵蝕，與A-mos等［9］實(shí)驗(yàn)觀測(cè)到的一致，長(zhǎng)期的內(nèi)岸下游沉積和外岸下游侵蝕會(huì)造成峽谷整體形狀的變化，即彎曲度和峽谷寬度的變化，以及峽谷軸線的遷移。在峽谷彎道頂點(diǎn)的外岸漫灘及其下游方向上，由于離心力的作用產(chǎn)生較多的沉積，并呈由上游向下游逐漸遞減的趨勢(shì)。本研究給出了海底彎曲峽谷中一次濁流事件的基本沉積特征，可見其趨勢(shì)與實(shí)驗(yàn)觀察及海底實(shí)測(cè)的一些研究相符合［9～11］，今后我們將進(jìn)一步研究多次濁流事件所形成的沉積序列。

4 總結(jié)

本研究對(duì)圓弧形彎曲海底峽谷中的濁流的流動(dòng)及其沉積進(jìn)行了三維數(shù)值計(jì)算模擬和分析，得到如下主要結(jié)論:濁流在運(yùn)動(dòng)過程中不斷夾帶周圍流體使?jié)崃骱穸仍黾?，隨著濁流厚度超過峽谷深度，產(chǎn)生溢出;濁流到達(dá)彎道部后，由于離心力的作用會(huì)產(chǎn)生更多的溢出，其最大過量密度可達(dá)入流的37.5%;彎道頂點(diǎn)橫剖面處的二次流在峽谷底部形成了一個(gè)類似于河流二次流的循環(huán)，靠近底床處的方向從外岸指向內(nèi)岸;彎道外彎處岸邊觀察到分流。觀察到的濁流在水下彎曲峽谷的沉積特征為:每個(gè)彎道頂點(diǎn)內(nèi)岸下游處有由于回流產(chǎn)生的沉積，而在外岸下游處則發(fā)生侵蝕，在峽谷彎道頂點(diǎn)的外岸漫灘及其下游方向上，由于離心力的作用產(chǎn)生較多的沉積，并呈現(xiàn)由上游向下游逐漸遞減的趨勢(shì)。本研究的模擬結(jié)果不僅確認(rèn)了一些已知的實(shí)驗(yàn)事實(shí)，并且提供了濁流在海底峽谷內(nèi)流動(dòng)及沉積相關(guān)聯(lián)的一些新的信息，有助于進(jìn)一步加深對(duì)濁流運(yùn)動(dòng)和沉積關(guān)系的理解和認(rèn)識(shí)。

References)

1 Middleton G V.Sediment deposition from turbidity currents［J］.Annual Review of Earth and Planetary Sciences，1993，21:89-114

2 Meiberg E，Kneller B.Turbidity currents and their deposits［J］.Annual Review of Fluid Mechanisms，2010，42:135-156

3 Heezen C B，Ewing M.Turbidity currents and submarine slumps，and the 1929 Grand Banks earthquake［J］.American Journal of Science，1952，250:849-873

4 姜輝.濁流沉積的動(dòng)力學(xué)機(jī)制與響應(yīng)［J］.石油與天然氣地質(zhì)，2010，31(4):428-435［Jiang Hui.Dynamical mechanism and depositional responses of turbidity current sedimentation［J］.Oil＆ Gas Geology，2010，31(4):428-435］

5 楊川恒，杜栩，潘和順，等.國外深水領(lǐng)域油氣勘探新進(jìn)展及我國南海北部陸坡深水區(qū)油氣勘探潛力［J］.地學(xué)前緣，2000，7(3):247-256［Yang Chuanheng，Du Xu，Pan Heshun，et al.Advances in worldwide deep water hydrocarbon exploration and oil and gas exploration potential in the northen continental slope in South China Sea［J］.Earth Science Frontiers，2000，7(3):247-256］

6 何家雄，夏斌，施小斌，等.世界深水油氣勘探進(jìn)展與南海深水油氣勘探前景［J］.天然氣地球科學(xué)，2006，17(6):747-752［He Jiaxiong，Xia Bin，Shi Xiaobin，et al.Prospect and progress for oil and gas in deep waters of the world，and the potential and prospect foreground for oil and gas in deep waters of the South China Sea［J］.Natural Gas Geoscience，2006，17(6):747-752］

7 Peakall J，Amos K J，Keevil G M，et al.Flow processes and sedimentation in submarine channel bends［J］.Marine and Petroleum Geology，2007，24:470-486

8 Huang H Q，Imran J，Pirmez C.Numerical model of turbidity currents with a deforming bottom boundary［J］.Journal of Hydraulic Engineer-ing，2005，131(4):283-293

9 Amos K J，Peakall J，Bradbury P W，et al.The influence of bend amplitude and planform morphology on flow and sedimentation in submarine channels［J］.Marine and Petroleum Geology，2010，27:1431-1447

10 Straub K M，Mohrig D，Buttles J，et al.Quantifying the influence of Channel sinuosity on the depositional mechanics of channelized turbidity currents:A laboratory study［J］.Marine and Petroleum Geology，2011，28:744-760

11 Wynn R B，Cronin B T，Peakall J.Sinuous deep-water channels:Genesis，geometry and architecture［J］.Marine and Petroleum Geology，2007，24:341-387

12 Leopold L B，Langbein W B.River meanders［J］.Geological Society of America Bulletin，1966，71:769-793

13 Pirmez C.Growth of a submarine meandering channel-Levee system on the Amazon Fan［D］.New York:Columbia University，1994

14 Prior D B，Bornhold B D，Wiseman W J，et al.Turbidity current activity in a British Columbia Fjord［J］.Science，1987，237:1330-1333

15 XU J P，Noble A M，Rosenfeld K L.In-situ measurements of velocity structure within turbidity currents［J］.Geophysical Research Letters，2004，31:1-4

16 Huang H Q，Imran J，Pirmez C.Numerical modeling of poorly sorted depositional turbidity currents［J］.Journal of Geophysical Research，2007，112:1-15

17 Smith J D，McLean S R.Spatially averaged flow over a wavy surface［J］.Journal of Geophysical Research，1977，82(12):1735-1746

18 Rijn，L C van.Sediment transport，part II:suspended load transport［J］.Journal of Hydraulic Engineering，1984，110(11):1613-1641

19 Huang H Q，Imran J，Pirmez C.Numerical study of turbidity currents with sudden-release and sustained-inflow mechanisms［J］.Journal of Hydraulic Engineering，2008，134(9):1199-1209

20 Huang H Q，Imran J，Pirmez C.The depositional characteristics of turbidity currents in submarine sinuous channels［J/OL］.Marine Geology，2012，329-331:93-102

21 Kassem A，Irman J.Three-dimensional modeling of density current.II.Flow in sinuous confined and unconfined channels［J］.Journal of Hydraulic Research，2004，42(6):591-602

22 Kane A I，McCaffrey W D，Peakall J，et al.Submarine channel levee shape and sediment waves from physical experiments［J］.Sedimentary Geology，2010，223:75-85

23 Garcia M H.Depositional turbidity currents laden with poorly sorted sediment［J］.Journal of Hydraulic Engineering，1994，120(11):1240-1263

24 Keevil G M，Peakall J，Best J L，et al.Flow structure in sinuous submarine channels:velocity and turbulence structure of an experimental submarine channel turbulence structure of an experimental submarine channel［J］.Marine Geology，2006，229:241-257

25 Corney K P，Peakall J，Parsons D R，et al.Reply to discussion of Imran et al.on"the orientation of helical flow in curved channels"by Corney et al.，Sedimentology，53，249-257［J］.Sedimentology，2008，55:241-247

26 Corney K P，Peakall J，Parsons D R，et al.The orientation of helical flow in curved channels［J］.Sedimentology，2006，53:249-257

27 Abad J D，Sequeiros O E，Spinewine B，et al.Secondary current of saline underflow in a highly meandering channel:experiments and theory［J］.Journal of Sedimentary Research，2011，81:787-813

28 Serchi G F，Peakall J，Ingham B D，et al.A unifying computational fluid dynamics investigation on the rive-like secondary circulation in submarine channel bends［J］.Journal of Geophysical Research，2011，116:1-19