(國家海洋局 南海工程勘察中心,廣東 廣州 510300)

沉積物的粒度分析對(duì)于研究沉積物的來源、運(yùn)移機(jī)制,了解沉積環(huán)境具有十分重要的意義。沉積物顆粒大小受流水營力的控制,粒度分布特征受物源和沉積水動(dòng)力環(huán)境控制。從沉積物的粒度特征中能提取沉積環(huán)境的有關(guān)信息,或進(jìn)行沉積環(huán)境的辨識(shí)[1-4]。如石英顆粒平均粒徑與重礦物顆粒平均粒徑比率可用來區(qū)分海灘砂和沙丘砂,并作為一種區(qū)分風(fēng)成成因還是水成成因的途徑。利用C-M圖區(qū)分出砂的推移和懸移搬運(yùn)方式[5]。隨著激光粒度儀的廣泛應(yīng)用,使得沉積物粒級(jí)頻率分布的細(xì)微變化也能反映出來。海洋沉積物粒度與沉積環(huán)境的關(guān)系的研究也逐漸多樣化。包括利用沉積物參數(shù)特征探討水動(dòng)力環(huán)境[6-7]和物源[8];利用粒徑趨勢(shì)分析方法來研究沉積物輸運(yùn)途徑[9-12]。而多元統(tǒng)計(jì)方法和地統(tǒng)計(jì)分析方法引入到沉積物粒度分析中[13-15],使得定量研究沉積物粒度與水動(dòng)力的關(guān)系成為可能。與此同時(shí),對(duì)潮間帶(潮灘)沉積物粒度也開展了大量的研究工作,如潮灘沉積物粒度的時(shí)空分布特征[16-19]、潮間帶水動(dòng)力特征[20-22]等。

本文利用粒度參數(shù)、粒徑分布和多元統(tǒng)計(jì)分析方法來探討湛江海區(qū)海島潮間帶表層沉積物的運(yùn)移方式和沉積環(huán)境。

1 研究區(qū)概況及分析方法

1.1 研究區(qū)概況

湛江海區(qū)地處109°30′～110°40′E,21°10′～21°30′N,該海區(qū)內(nèi)分布大小島嶼(不含礁石)共 31個(gè)(據(jù)廣東省908海島統(tǒng)計(jì)數(shù)據(jù)),絕大多數(shù)海島屬于堆積成因島嶼。湛江海區(qū)屬不規(guī)則半日潮,平均潮差2 m左右,潮流類型為往復(fù)流,流向受岸線和海底地形控制。漲潮流歷時(shí)長于落潮流歷時(shí),落潮流速大于漲潮流速[23]。沿岸表層海流,冬季以SW向漂流為主;夏季以NE向漂流為主,近岸范圍有一股狹窄的SW向沿岸密度流[23]。

現(xiàn)場(chǎng)調(diào)查表明,湛江灣及附近海域,除特呈島出露的潮間帶主要為紅樹林灘外,其他海島如東海島、南三島、硇洲島等出露的潮間帶類型以沙灘為主;雷州半島東部海域,北莉島、東松島、六極島等諸島出露潮間帶以粉砂-淤泥質(zhì)潮灘為主,局部岸段出露紅樹林潮灘。

1.2 分析方法

分析樣品來自 2008年6～12月在廣東省湛江市管轄包括東海島、南三島、硇洲島、北莉島、東松島等在內(nèi)的29個(gè)海島上布設(shè)60條調(diào)查剖面,分別編號(hào)1～60,其中1～43號(hào)剖面采集高、中、低潮表層樣品(每條剖面3個(gè)樣品),44～60號(hào)剖面僅采集中潮表層樣品(每條剖面1個(gè)樣品),共計(jì)136個(gè)樣品,采樣剖面編號(hào)示意圖見圖1。沉積物樣品均用采樣勺采集,為保證粒度分析的有效性,所有沉積物樣品都在低潮剖面出露時(shí)采集,并控制在表層下 0～10cm 深度范圍內(nèi)。室內(nèi)分析流程如下：首先利用篩析法分析大于 0.063 mm的粗顆粒物不同粒級(jí)的含量,將小于0.063 mm的細(xì)顆粒物置于燒杯中,加入15 mL 3%的雙氧水浸泡 24 h,去除有機(jī)質(zhì),然后加入 5 mL 3 mol/L的稀鹽酸浸泡24 h,去除沉積物中的鈣質(zhì)膠結(jié)物及生物貝殼,其后將樣品進(jìn)行反復(fù)離心、洗鹽直至溶液呈中性為止。處理后的樣品經(jīng)超聲波振蕩分散后上機(jī)測(cè)試。粒度分析使用的儀器為英國Malvern公司生產(chǎn)的Mastersizer 2000 型激光粒度儀,測(cè)量范圍為0.02～2 000 μm。粒級(jí)統(tǒng)一使用尤登—溫德華氏等比值粒級(jí)標(biāo)準(zhǔn)表示,平均粒徑(Mz)、標(biāo)準(zhǔn)偏差(σi)、偏態(tài)(Ski)、峰態(tài)(Kg)等粒度參數(shù)采用Folk—Ward圖解法公式進(jìn)行計(jì)算[24-25]。

圖1 湛江海區(qū)海島潮間帶表層沉積物剖面編號(hào)Fig.1 The map showing the coding of island intertidal surface sediment sampling section,Zhanjiang offshore area

2 結(jié)果與討論

2.1 沉積物命名

研究區(qū)沉積物的分類和命名采用?？说热说某练e物粒度三角圖解法對(duì)不含礫和含礫沉積物分別命名[26]。根據(jù) Folk沉積物分類結(jié)果,研究區(qū)潮間帶不含礫沉積物主要為砂和粉砂質(zhì)砂;含礫沉積物主要為礫質(zhì)砂,其次為含礫泥質(zhì)砂和含礫砂。粉砂質(zhì)砂主要發(fā)現(xiàn)于中潮帶和低潮帶,高潮帶沉積物類型基本為砂。含礫泥質(zhì)砂和含礫砂主要發(fā)現(xiàn)于中潮帶和低潮帶,高潮帶主要沉積物類型為礫質(zhì)砂(圖2)。

2.2 表層沉積物粒度參數(shù)特征

沉積物粒度參數(shù)的計(jì)算,Folk-Ward 公式物理意義明確,精確度很高,應(yīng)用最為廣泛[27]。計(jì)算的結(jié)果表明研究區(qū)潮間帶表層沉積物平均粒徑為-0.36φ～5.98φ,平均值為2.55φ;標(biāo)準(zhǔn)偏差為0.43～4.76,平均值為 1.38;偏態(tài)為-4.39～0.72,平均值為0.10;峰態(tài)為0.10～6.63,平均值為1.36。

Folk和 Ward建立了一套使用標(biāo)準(zhǔn)偏差、偏態(tài)和峰態(tài)來評(píng)價(jià)沉積物分選性、粒度分布對(duì)稱性和粒度頻率曲線尖銳程度標(biāo)準(zhǔn)[28],本文使用該標(biāo)準(zhǔn)評(píng)價(jià)研究區(qū)潮間帶表層沉積物粒度特征。從圖3可以看出沉積物平均粒徑主要分布在 1φ～3φ之間;多數(shù)沉積物分選較差,其次為分選較好和分選差,少數(shù)為分選好和分選中等;粒度分布曲線有一半以上呈近對(duì)稱狀,其余曲線基本呈正偏態(tài),尖銳度多數(shù)呈中等尖銳狀,其次呈很尖銳狀和尖銳狀。

2.3 表層沉積物粒度C-M圖解

以累積 1%最粗顆粒粒徑值和中值粒徑值投影在對(duì)數(shù)坐標(biāo)平面上,得到研究區(qū)海島潮間帶表層沉積物粒度C-M圖[5,25]。

利用粒度象分析沉積作用水動(dòng)力狀況,可以反映一定粒級(jí)顆粒主要搬運(yùn)機(jī)制。根據(jù)常見的搬運(yùn)和沉積類型,利用C=1000、M=200,100和15 μm將沉積物分為Ⅰ～Ⅸ區(qū)[25]。從C-M圖上看(圖4),研究區(qū)海島潮間帶沉積類型主要分為 3類：第一類分布在Ⅰ區(qū)、Ⅱ區(qū)和Ⅲ區(qū),含滾動(dòng)顆粒,為靠近源區(qū)處沉積;第二類分布在Ⅳ區(qū)和Ⅴ區(qū),以中等—高擾動(dòng)遞變懸浮沉積為主;第三類分布在Ⅵ區(qū)和Ⅶ區(qū),為遞變懸浮和均勻懸浮沉積。

圖2 湛江海島潮間帶表層沉積物?？朔诸惾菆DFig.2 The Folk classification map of island intertidal surface sediments,Zhanjiang offshore area

圖3 湛江海島潮間帶表層沉積物平均粒徑、分選性、分布曲線對(duì)稱性和尖銳度頻率分布圖Fig.3 The frequency distribution maps of mean grain-size,grading and symmetry and acutance of distribution curve of the island intertidal surface sediment of Zhanjiang offshore area

圖4 湛江海島潮間帶表層沉積物粒度C-M圖解Fig.4 The grain-size C-M diagram of island intertidal surface sediments of Zhanjiang offshore area

2.4 表層沉積物粒度頻率曲線與概率累積曲線

研究區(qū)海島潮間帶表層沉積區(qū)粒度頻率曲線主要為單峰態(tài)和雙峰態(tài),兩者所占比例大致相等,三峰態(tài)則相對(duì)較少。單峰態(tài)、雙峰態(tài)對(duì)應(yīng)的概率累積曲線多數(shù)為三段式,三峰態(tài)對(duì)應(yīng)的概率累積曲線多為二段式。典型粒度頻率曲線與概率累積曲線見圖5。

三段式曲線由躍移組分和懸移組分組成,普遍缺乏滾動(dòng)組分,從躍移組分延伸的情況看,即使存在滾動(dòng)組分,其含量一般在 1%～2%以下;躍移組分具有明顯的兩段式特征;粗細(xì)兩段的截點(diǎn)位于4φ～6φ之間;粗粒段的斜率明顯大于細(xì)粒段,即粗粒段分選比細(xì)粒段好。這反映了潮間帶雙向流沉積作用;粗粒段代表動(dòng)力較強(qiáng)的落潮流沉積,細(xì)粒段代表動(dòng)力較弱的漲潮流沉積。懸移組分斜率明顯地小于躍移組分,其含量一般在 10%以下,躍—懸移組分分界點(diǎn)一般位于 6φ～8φ之間。

圖5 湛江海島潮間帶表層沉積物典型粒度頻率曲線與概率累積曲線Fig.5 Typical grain-size frequency curves and probability cumulative curves of island intertidal surface sediments of the Zhanjiang offshore area

二段式曲線由躍移組分和懸移組分組成,其截點(diǎn)變化區(qū)間較寬(一般介于 2φ～4φ)。截點(diǎn)的大小可以反映搬運(yùn)介質(zhì)的擾動(dòng)強(qiáng)度,強(qiáng)度高的在較粗粒度上發(fā)生截?cái)?。從截點(diǎn)區(qū)間來看,調(diào)查區(qū)二段式曲線所對(duì)應(yīng)沉積物搬運(yùn)強(qiáng)度為中等～高。

一段式由單一的直線段組成,為單一方式搬運(yùn),這種直線的斜度都小于 50°,在寬的區(qū)間內(nèi)延伸,本類曲線可能代表水速及密度均較大,全部沉積物均呈懸浮負(fù)載。

粒度累積曲線與C-M圖存在較好的對(duì)應(yīng)關(guān)系,如一段式全部分布在C-M圖中的Ⅰ區(qū)、Ⅱ區(qū)和Ⅲ區(qū)內(nèi);二段式主要分布在C-M圖中的Ⅳ區(qū)和Ⅴ區(qū)內(nèi);三段式較復(fù)雜,一部分分布在C-M圖中Ⅳ區(qū)和Ⅴ區(qū)內(nèi),另一部分分布在Ⅵ區(qū)和Ⅶ區(qū)內(nèi)。

2.5 表層沉積物粒度聚類分析和因子分析

2.5.1 聚類分析

為分析表層沉積物的粒級(jí)組成及其相互關(guān)系,對(duì)研究區(qū)海島潮間帶表層沉積物粒度數(shù)據(jù)進(jìn)行了 R型聚類分析和 R 型因子分析。沉積物φ粒度聚類分析結(jié)果見圖6,以距離15為界,14個(gè)粒級(jí)基本可劃分為 4 組：＜-1φ～2φ、3φ、4φ和 5φ～＞11φ。沉積物相關(guān)性分析結(jié)果見表1,根據(jù)相關(guān)性可以將粒徑分為兩組：一組為細(xì)粒組分(5φ～＞11φ),另一組為粗粒組分(＜-1φ～4φ)。細(xì)粒組分之間的相關(guān)性較大,該粒徑范圍代表粉砂和黏土;粗粒組分之間的相關(guān)性較小,該粒徑范圍代表礫石和砂;粗、細(xì)組分之間的相關(guān)性也較小。

粉砂和黏土為沉積物中的懸移(包括遞變懸移和均勻懸移)組分,系弱水動(dòng)力沉積的產(chǎn)物,含量較穩(wěn)定,因此相關(guān)性也很大。礫石和砂為沉積物中的推移組分,在機(jī)械搬運(yùn)和沉積作用過程中,一般呈側(cè)向加積,系較強(qiáng)水動(dòng)力沉積的產(chǎn)物,可能由于海灘受海水作用而被擾動(dòng),再沉積過程中,粗粒物質(zhì)來不及分異,而混雜在一起,含量變化較大,因此相關(guān)性也較小。

圖6 湛江海島潮間帶沉積物粒度R型聚類分析Fig.6 The grain-size R-type clustering analysis of island intertidal surface sediment of the Zhanjiang offshore area

表1 湛江海區(qū)海島表層沉積物各粒級(jí)相關(guān)系數(shù)Tab.1 The correlation coefficients among different grain-size grades of island intertidal surface sediments of the Zhanjiang offshore area

2.5.2 因子分析

研究區(qū)海島潮間帶表層沉積物粒度資料的 R型因子分析主要結(jié)果見表2。表2中方差極大因子載荷表明,因子1特征值比例占49.89%,因子2特征值比例占16.88%,因子3特征值比例占14.397%,因子4特征值比例占7.84%,前4個(gè)因子特征值比例累計(jì)占88.58%。將粒級(jí)變量與 4個(gè)主因子載荷值投影到坐標(biāo)平面上,得到因子載荷圖(圖7)。

表2 湛江海島潮間帶表層沉積物粒度方差極大旋轉(zhuǎn)因子載荷Tab.2 The grain-size varimax rotation factor loading values of island intertidal surface sediments,the Zhanjiang offshore area

圖7 湛江海島潮間帶表層沉積物粒度因子載荷圖Fig.7 The grain-size factor loading map of island intertidal surface sediments of the Zhanjiang offshore area

因子1主要由5φ～＞11φ的正載荷組成,該粒徑范圍代表黏土和粉砂組分;因此,因子1在黏土和粉砂組分上有較大的正載荷。從因子載荷圖上看,該粒徑范圍以懸移沉積為主,為極細(xì)組分的均勻懸移沉積。

因子 2主要由＜-1φ～1φ的正載荷組成,該粒徑范圍為沉積物中的極粗砂到粗砂組分;因此,因子 2在極粗砂到粗砂組分上有較大的正載荷。從因子載荷圖上看,該粒徑范圍以推移沉積為主。

因子 3主要由 4φ的正載荷和 2φ的負(fù)載荷組成,4φ是極細(xì)砂與粗粉砂的分界,2φ是中砂與細(xì)砂的分界;因此,因子3在極細(xì)砂與粗粉砂的分界處有較大的正載荷,而在中砂與細(xì)砂的分界處有較大的負(fù)載荷。從因子載荷圖上看,該粒徑范圍代表以躍移沉積為主。

因子 4主要由 3φ的負(fù)載荷組成,3φ是細(xì)砂與極細(xì)砂的分界;因此因子 4在細(xì)砂與極細(xì)砂的分界上有較大的負(fù)載荷。從因子載荷圖上看,該粒徑范圍以躍移沉積為主。

3 結(jié)論

(1)粒度概率累積曲線反映出 3類水動(dòng)力環(huán)境：第一類曲線為三段式,搬運(yùn)方式包括懸移和躍移,躍—懸移分界點(diǎn)為 6φ～8φ,由于潮間帶雙向流沉積,根據(jù)水動(dòng)力強(qiáng)弱,分為兩段,截點(diǎn)為 4φ～6φ;第二類曲線為二段式,搬運(yùn)方式包括懸移和躍移,躍—懸移分界點(diǎn)2φ～4φ,反映的沉積物搬運(yùn)強(qiáng)度中等;第三類曲線為單一直線段,水動(dòng)力強(qiáng),僅懸浮一種搬運(yùn)方式。第一、二類曲線代表強(qiáng)、弱水動(dòng)力共同作用的沉積;第三類曲線為單一強(qiáng)水動(dòng)力沉積。(2)粒度R型聚類分析結(jié)果表明弱水動(dòng)力沉積的細(xì)顆粒物含量穩(wěn)定,不同粒級(jí)間相關(guān)性大;強(qiáng)動(dòng)力沉積的粗顆粒物由于海水?dāng)_動(dòng),含量變化大,不同粒級(jí)間相關(guān)性小。(3)R型因子分析與粒度概率曲線有較好的一致性。從粒度R型因子分析來看,因子1由5φ～＞11φ的正載荷組成,為均勻懸移沉積;因子2由＜-1φ～1φ的正載荷組成,以推移沉積為主;因子 3由 4φ的正載荷和2φ的負(fù)載荷組成;因子4由3φ的負(fù)載荷組成,以躍移沉積為主。

致謝：兩位審稿專家提出了寶貴的修改意見,在此表示衷心的感謝！

[1]Visher G S.Grain seze distributions and depositional processes[J].Journal of Sedimentary Petrology,1969,39：1074-1106.

[2]Friedman G M.Address of the retiring president of the international association of sedimentology：difference in size distributions of populations of particles among sands from various origins[J].Sedimentology,1979,26：3-32.

[3]MeLaren P,Bowles D.The effects of sediment transport on grain-size distributions[J].Journal of Sedimentary Petrology,1985,55：457-470.

[4]張存勇,馮秀麗.連云港近岸海域沉積物粒度空間分布特征及其分析[J].海洋學(xué)報(bào),2009,4：120-127.

[5]Passage R.Grain size representation by CM pattern as a geologic tool[J].Journal of Sedimentary Petrology,1964,34：830-847.

[6]王偉,李安春,徐方建,等.北黃海表層沉積物粒度分布特征及其沉積環(huán)境分析[J].海洋與湖沼,2009,40(5)：525-531.

[7]楊世倫.長江口沉積物粒度參數(shù)的統(tǒng)計(jì)規(guī)律及其沉積動(dòng)力學(xué)解釋[J].泥沙研究,1994,3：23-31.

[8]王偉偉,趙全民,王權(quán)明,等.長興島鄰近 海域的沉積物粒度參數(shù)特征及物源分析[J].海洋環(huán)境科學(xué),2010,29(2)：272-275.

[9]Gao S(高抒),Michael C.Net sediment transport patterns inferred from grain size trends,based upon definition of transport vectors[J].Sedimentary Geology,1992,81：47-60.

[10]梁娟,李春初,王世俊.珠江磨刀門河口底質(zhì)沉積特征及其泥沙運(yùn)移趨勢(shì)[J].海洋通報(bào),2006,25(5)：57-63.

[11]王國慶,石學(xué)法,劉焱光,等.粒徑趨勢(shì)分析對(duì)長江南支口外沉積物輸運(yùn)的指示意義[J].海洋學(xué)報(bào),2007,29(6)：161-166.

[12]程鵬,高抒.北黃海西部海底沉積物的粒度特征和凈輸運(yùn)趨勢(shì).[J].海洋與湖沼,2000,31：604-615.

[13]高為利,張富元,章偉艷,等.海南島周邊海域表層沉積物粒度分布特征[J].海洋通報(bào),2009,28(2)：71-80.

[14]李玉中,陳沈良.系統(tǒng)聚類分析在現(xiàn)代沉積環(huán)境劃分中的應(yīng)用——以崎嶇列島海區(qū)為例[J].沉積學(xué)報(bào).2003,21(3)：487-494.

[15]劉付程,張存勇,俊 彭.海州灣表層沉積物粒度的空間變異特征[J].海洋科學(xué),2010,34(7)：54-58.

[16]李平,陳沈良,谷國傳,等.長江口九段沙上沙潮灘沉積物短期時(shí)空變化特征[J].海洋通報(bào),2008,27(6)：60-67.

[17]劉旭英,高建華,白鳳龍,等.蘇北新洋港潮灘柱狀沉積物粒度分布特征 [J].海洋地質(zhì)與第四紀(jì)地質(zhì),2008,28(4)：27-35.

[18]曹惠美,蔡鋒,蘇賢澤.華南沿海若干砂質(zhì)海灘沉積物粒度特征的分析[J].海洋通報(bào),2005,24(4)：36-45.

[19]薛允傳,尹延鴻.黃河三角洲北部潮間帶沉積物的粒度特征[J].海洋科學(xué),2001,25(5)：50-54.

[20]黃廣,陳沈良,胡靜.南匯東灘沉積物粒度特征及其與水動(dòng)力的關(guān)系[J].海洋湖沼通報(bào),2008,1：32-38.

[21]蔣國俊.潮灘懸沙粒度參數(shù)的動(dòng)力沉積學(xué)意義[J].海洋與湖沼,1995,26(1)：90-97.

[22]王柯道,王建,何加武.茅家港灘面沉積物粒度變化及其與水動(dòng)力的關(guān)系[J].海洋科學(xué),2005,29(5)：73-79.

[23]趙煥庭,張喬民,宋朝景,等.華南海岸和南海諸島地貌與環(huán)境[M].北京：科學(xué)出版社,1999.

[24]Folk R L.A review of grain-size parameters[J].Sedimentology,1966,6：73-97.

[25]成都地質(zhì)學(xué)院陜北隊(duì).沉積巖(物)粒度分析及其應(yīng)用[M].北京：地質(zhì)出版社,1978.

[26]Folk R L,Andrews P B,Lewis D W.Detrital sedimentary rock classification and nomenclature for use in New Zealan[J].New Zealand Journal of Geology and Geophysics.1970,13：937-968.

[27]賈建軍,高抒,薛允傳.圖解法與矩法沉積物粒度參數(shù)的對(duì)比[J].海洋與湖沼,2002,33(6)：577-582.

[28]Folk R L,Ward W C.Brazos river bar：a study in the signification of grain size parameters[J].Journal of Sedimentary Petrology,1957,27：3-27.