李長(zhǎng)征 李姝昱 王銳 楊勇

摘要：河底泥沙密度、粒徑等是河流動(dòng)力學(xué)研究的必備物性參數(shù)。機(jī)械取樣和淺地層剖面探測(cè)是獲取泥沙物性的重要途徑。將取樣測(cè)試的泥沙粒徑信息與淺地層剖面探測(cè)數(shù)據(jù)結(jié)合起來(lái)，研究了一種快速的河底泥沙物性參數(shù)（密度、孔隙度等）檢測(cè)方法。利用級(jí)配結(jié)果確定表征孔隙度和滲透率關(guān)系的Kozeny-Carman公式，將該式與孔隙介質(zhì)模型結(jié)合計(jì)算河底泥沙密度等參數(shù)。采用該方法獲取了三門峽庫(kù)區(qū)2個(gè)斷面的濕密度和波速等參數(shù)。將泥沙濕密度計(jì)算值與取樣測(cè)試值對(duì)比，發(fā)現(xiàn)2個(gè)測(cè)點(diǎn)的差異較小，證明了該方法的有效性。

關(guān)鍵詞：孔隙介質(zhì)；淺地層剖面；快速估算；物理特性；表層泥沙；河底

中圖分類號(hào)：P332.5 文獻(xiàn)標(biāo)志碼：A doi：10.3969/j.issn.1000-1379.2018.03.004

研究河底泥沙物理特性（密度、孔隙度等）對(duì)河流動(dòng)力學(xué)研究和庫(kù)區(qū)清淤等有重要意義。目前，河底泥沙物性參數(shù)的獲取方式主要有3種：簡(jiǎn)單機(jī)械取樣、保真采樣和聲學(xué)探測(cè)。簡(jiǎn)單機(jī)械取樣是在船上吊裝一定質(zhì)量的取樣器到水底，在水底取樣后將樣品運(yùn)送到實(shí)驗(yàn)室測(cè)試巖土力學(xué)參數(shù)。由于在取樣過(guò)程中不可避免地對(duì)樣品產(chǎn)生較大程度的擾動(dòng)，使樣品（尤其是稀軟的水底表層淤積物）的原始狀態(tài)發(fā)生變化，因此最終使測(cè)量結(jié)果產(chǎn)生較大的偏差。簡(jiǎn)單機(jī)械取樣為點(diǎn)測(cè)方式，速度較慢，取樣較淺（約為30cm），難以滿足底質(zhì)調(diào)查的需求。

以保真為目標(biāo)的取樣技術(shù)正在逐步完善和發(fā)展。為了避免沉積物樣品原始層序混雜或彎曲樣品形狀的擾動(dòng)，實(shí)現(xiàn)采集到的樣品保壓、保溫、無(wú)污染，保證樣品信息記錄的完整性，保真采樣設(shè)備研究經(jīng)歷了從國(guó)外研究到國(guó)內(nèi)自行設(shè)計(jì)和完善階段，浙江大學(xué)和國(guó)家海洋局第一海洋研究所設(shè)計(jì)的采樣器均可實(shí)現(xiàn)10 m以上的巖心采樣[1-2]。2013年，黃河水利科學(xué)研究院設(shè)計(jì)了針對(duì)河底淤積層的水下采樣器[3-4]，采用自動(dòng)切割的菊花瓣式的活塞刀口，保證了已取樣品的完整性，目前在細(xì)沙區(qū)（粒徑小于0.025mm）能夠獲取10m長(zhǎng)的泥沙樣品。

近年來(lái)，淺地層剖面儀逐步應(yīng)用在河流底質(zhì)調(diào)查中[5-6]。通過(guò)聲學(xué)換能器向水底發(fā)射寬帶調(diào)頻脈沖[7]，脈沖遇到聲阻抗突變界面（水底不同媒質(zhì)層的分界面）反射回到接收換能器，研究人員根據(jù)接收的聲吶信號(hào)可大致判斷河底淤積層的結(jié)構(gòu)。然而，河底泥沙淤積層的聲速仍是未知參數(shù)，僅依據(jù)聲學(xué)接收信號(hào)難以判斷淤積層的具體厚度，淤積層的其他巖土力學(xué)參數(shù)也不能精確檢測(cè)。

隨著聲學(xué)理論的發(fā)展，通過(guò)分析聲學(xué)信號(hào)獲取表層泥沙物性參數(shù)已具備可行性。本文利用保真取樣獲取的泥沙參數(shù)作為先驗(yàn)信息，結(jié)合孔隙介質(zhì)聲學(xué)理論快速得到水底淤積泥沙密度等參數(shù)，對(duì)提高河底淤積表層探測(cè)效率有重要意義。

1 聲學(xué)模型參數(shù)

1.1 模型參數(shù)

Biot M.A.[8]的孔隙介質(zhì)理論較好地解釋了由泥沙顆粒和孔隙水組成的雙相介質(zhì)系統(tǒng)中的聲學(xué)現(xiàn)象。因此，本文基于Biot孔隙介質(zhì)聲學(xué)模型（簡(jiǎn)稱Biot模型）計(jì)算河底表層泥沙密度等參數(shù)。Biot模型包含的參數(shù)主要分兩類：第一類是各組成相的參數(shù)，該類參數(shù)是物理常數(shù)，可以通過(guò)試驗(yàn)測(cè)定或參考相關(guān)文獻(xiàn)，見(jiàn)表1；第二類是與孔隙空間中流體流動(dòng)相關(guān)的參數(shù)，其中流體黏滯系數(shù)較容易測(cè)定，而滲透率、孔隙半徑和有效物質(zhì)常數(shù)較難精確測(cè)定，但可以通過(guò)經(jīng)驗(yàn)公式與其他參數(shù)建立聯(lián)系。

現(xiàn)場(chǎng)探測(cè)時(shí)需要接收信號(hào)有足夠的時(shí)長(zhǎng)，以便有效接收二次波信號(hào)。

2.2 泥沙物性參數(shù)計(jì)算方法

由表1確定Biot模型第一類參數(shù)，并且第二類參數(shù)相互關(guān)系已知（式（1）～式（6）），因此Biot模型里僅孔隙度和滲透率為未知參數(shù)。理論聲波反射系數(shù)計(jì)算公式涉及的輔助公式較多，本文不詳細(xì)列出，具體可參見(jiàn)文獻(xiàn)[10]。

計(jì)算過(guò)程和方法：①利用機(jī)械取樣器對(duì)表層淤積泥沙取樣后，通過(guò)實(shí)驗(yàn)室測(cè)試得到樣品的粒徑d；②將粒徑d代入式（7）；③從淺地層剖面數(shù)據(jù)中提取各測(cè)點(diǎn)的聲波反射系數(shù)；④聯(lián)合式（7）和式（12），以及參考文獻(xiàn)[10]中的水-沙界面反射系數(shù)，得到孔隙度和滲透率；⑤由式（9）計(jì)算泥沙的濕密度等參數(shù)。

3 探測(cè)實(shí)例

3.1 探測(cè)儀器

利用美國(guó)Edgetech公司生產(chǎn)的3200-XS型淺地層剖面儀進(jìn)行現(xiàn)場(chǎng)探測(cè)，儀器包括拖魚、拖纜和甲板單元。拖魚型號(hào)為SB-216S，由1個(gè)發(fā)射換能器、2個(gè)接收換能器組成。發(fā)射信號(hào)的波束角（-3dB）為24°。換能器之間用隔聲板隔開(kāi)，其作用是阻止向下傳播的聲波對(duì)接收信號(hào)產(chǎn)生干擾，且盡量減小表面反射。換能器向其正下方水底發(fā)射調(diào)頻信號(hào)，該信號(hào)具有較強(qiáng)的地層穿透能力，在遇到阻抗不同的界面時(shí)將發(fā)生反射。接收換能器收到水底泥沙的反射信號(hào)，信號(hào)經(jīng)拖纜傳送到甲板單元。利用匹配濾波技術(shù)對(duì)回波信號(hào)脈沖進(jìn)行壓縮，并在顯示器上顯示出來(lái)。

試驗(yàn)中用船體牽引拖魚，為了得到較穩(wěn)定的信號(hào)，選擇風(fēng)浪較小且水面平靜的情況下探測(cè)。為了降低隨機(jī)噪聲的影響，GPS與甲板單元的主機(jī)相連，拖魚在行進(jìn)過(guò)程中，GPS數(shù)據(jù)與反射信號(hào)被記錄在同一文件中，因此可提取每個(gè)測(cè)點(diǎn)對(duì)應(yīng)的GPS坐標(biāo)。

3.2 探測(cè)結(jié)果

2015年11月，在黃河三門峽庫(kù)區(qū)進(jìn)行了機(jī)械取樣和淺地層剖面儀探測(cè)試驗(yàn)。根據(jù)現(xiàn)場(chǎng)淤積泥沙特征，選取的線性調(diào)頻信號(hào)頻帶為700～12000Hz，脈沖重復(fù)頻率為5Hz，拖魚水平行進(jìn)速度約為1.2m/s。測(cè)量時(shí)根據(jù)水流情況船速作相應(yīng)改變。垂直水流方向，連接左右岸的直線段為1個(gè)斷面（兩個(gè)斷面之間的距離約為2km），取樣和探測(cè)以斷面為單元進(jìn)行研究。

對(duì)三門峽庫(kù)區(qū)斷面2和斷面11進(jìn)行了探測(cè)，測(cè)線方向?yàn)樽蟀吨劣野叮瑴y(cè)線長(zhǎng)度為197m，斷面2水深約為25m。圖2為斷面2主河槽的聲納探測(cè)圖像，由圖像可見(jiàn)清晰的河底反射信號(hào)和表層泥沙底界面的反射信號(hào)，其中表層泥沙厚度為3.8～6.0m。斷面11的測(cè)線方向?yàn)橛野吨磷蟀叮瑴y(cè)線長(zhǎng)度約為201m，聲吶圖像見(jiàn)圖30主河槽水深約為15m，可見(jiàn)泥沙層界面的反射和多次波。表層泥沙厚0.5～2.0m，第2層泥沙厚4.5～5.0m。

在斷面2和斷面11分別選取6個(gè)、8個(gè)測(cè)點(diǎn)，提取到聲波反射系數(shù)。泥沙參數(shù)計(jì)算結(jié)果見(jiàn)表2和表3。斷面2的反射系數(shù)較斷面11的小，對(duì)應(yīng)較大的孔隙度和較低的泥沙密度。斷面2的泥沙干密度大部分低于1000kg/m3，斷面11的泥沙干密度為1100～1649kg/m3。在斷面2和斷面11進(jìn)行了機(jī)械取樣，取樣點(diǎn)的位置分別為A和B，測(cè)試到的泥沙中值粒徑分別為0.0110、0.0237mm。A、B點(diǎn)樣品的測(cè)試密度分別為1270、2150kg/m3，與附近測(cè)點(diǎn)密度計(jì)算結(jié)果差異較小。

4 結(jié)語(yǔ)

在黃河三門峽庫(kù)區(qū)進(jìn)行了水底泥沙的機(jī)械取樣和淺地層剖面探測(cè)，基于Biot模型將取樣測(cè)試的級(jí)配結(jié)果和聲吶數(shù)據(jù)結(jié)合起來(lái)，計(jì)算了表層泥沙的密度和波速等參數(shù)，并與取樣測(cè)試的結(jié)果進(jìn)行了對(duì)比，結(jié)果證明了本文計(jì)算方法的正確性。研究發(fā)現(xiàn)：①機(jī)械取樣對(duì)泥沙的級(jí)配測(cè)試結(jié)果影響較小，因此級(jí)配測(cè)試結(jié)果可作為聲學(xué)反演的先驗(yàn)信息；②本文提出的計(jì)算方法能夠計(jì)算得到河底泥沙的密度、波速和表層淤積厚度等參數(shù)，這對(duì)取樣測(cè)試結(jié)果起驗(yàn)證作用；③本文的探測(cè)方法能夠減少取樣工作量，并且得到多個(gè)泥沙物性參數(shù)。

本文采用的計(jì)算方法基于以下前提條件：①計(jì)算時(shí)考慮了聲學(xué)反射信號(hào)的弱頻率獨(dú)立性；②在探測(cè)過(guò)程中忽略了拖魚在水下姿態(tài)變化的影響；③利用的Kozzy-Carman方程為經(jīng)驗(yàn)公式?？紤]以上因素，計(jì)算中將產(chǎn)生一定誤差。

獲取較深層泥沙物性參數(shù)還存在一定的挑戰(zhàn)性，期望通過(guò)進(jìn)一步研究結(jié)合機(jī)械取樣得到更多的泥沙（聲吶圖像識(shí)別的第二層和更深層）物性參數(shù)，以滿足河流研究和治理的需求。

參考文獻(xiàn)：

[1]朱亮，顧臨怡，秦華偉.深海沉積物保真采樣技術(shù)及應(yīng)用[J].浙江大學(xué)學(xué)報(bào)（工學(xué)版），2005，39（7）：1060-1063.

[2]秦華偉，陳鷹，顧臨怡，等.海底沉積物保真采樣技術(shù)研究進(jìn)展[J].熱帶海洋學(xué)報(bào)，2009，28（4）：42-48.

[3]鄭軍，唐華，郭維克，等.小浪底庫(kù)區(qū)深層淤積泥沙物理特性分析[J].人民黃河，2014，36（10）：23-25.

[4]張雷，余孝志，鄭軍，等.小浪底庫(kù)區(qū)深層泥沙基本特性分析[C]//中國(guó)水利學(xué)會(huì).2014學(xué)術(shù)年會(huì)論文集（下冊(cè)）.南京：河海大學(xué)出版社，2014；919-923.

[5]孫永福，董立峰，宋玉鵬.黃河水下三角洲粉質(zhì)土擾動(dòng)土層特征及成因探析[J].巖土力學(xué)，2儀粥，四（6）：1494-1499.

[6]許國(guó)輝，衛(wèi)聰聰，孫永福，等.黃河水下三角洲淺表局部擾動(dòng)地層工程特性與成因[J].海洋地質(zhì)與第四紀(jì)地質(zhì)，2009，28（6）：19-25.

[7]謝向文，郭玉松，張曉予，等.水聲勘探中的淺地層剖面技術(shù)應(yīng)用[J].港工技術(shù)，2009，46（5）：53-55.

[8]RIOT M A.Theory of Propagation of Elastic Waves in a Fluid-Saturated Porous Solid.I.Low-Frequency Range[J].Journal ofthe Acoustical Society of America，1956，28（2）：179-191.

[9]WILLIAMS K L，JACKSON D R，THORSOS E I，et al.Comparison of Sound Speed and Attenuation Measured in aSandy Sediment to Predictions Based on the Biot Theory ofPorous Media[J].IEEE Journal of Oceanic Engineering，2002，27（3）：413-428.

[10]SCHOCK S G.A Method for Estimating the Physical andAcoustic Properties of the Sea Bed Using Chirp Sonar Data[J].IEEE Journal of Oceanic Engineering，2004，29（4）：1200-1217.

[11]Hovem J M，INGRAM G D.Viscous Attenuation of Soundin Saturated Sand[J].Journal of the Acoustical Society ofAmerica，1979，66（6）：1807-1812.