龍建軍，周華建，李趕先，陳聰

(1.廣東工業(yè)大學(xué) 機(jī)電工程學(xué)院，廣東 廣州 510006；2.中國(guó)科學(xué)院 南海海洋研究所 中科院邊緣海地質(zhì)重點(diǎn)實(shí)驗(yàn)室，廣東 廣州，510301)

?

受控三軸應(yīng)力-應(yīng)變下沉積物聲速與物理力學(xué)性質(zhì)的關(guān)系

龍建軍1，周華建1，李趕先2，陳聰1

(1.廣東工業(yè)大學(xué) 機(jī)電工程學(xué)院，廣東 廣州 510006；2.中國(guó)科學(xué)院 南海海洋研究所 中科院邊緣海地質(zhì)重點(diǎn)實(shí)驗(yàn)室，廣東 廣州，510301)

對(duì)南海某海域深度100～400 m的海底淺層(約2 m埋深范圍)沉積物柱狀樣在接近海底水壓力下進(jìn)行三軸應(yīng)變-聲學(xué)同步測(cè)量，結(jié)果表明沉積物縱波聲速有兩個(gè)特征：(1)從應(yīng)變過程開始到結(jié)束，沉積物縱波聲速不斷變化；(2)平均聲速隨著平均靜彈性模量的增加，由大變小又由小變大，存在聲速最小值。這些結(jié)果與海底淺表層沉積物的物理力學(xué)性質(zhì)、圍壓、顆粒的結(jié)合狀態(tài)改變有關(guān)。此外，沉積物動(dòng)彈性模量和孔隙度呈良好的負(fù)相關(guān)性，這與孔隙度增大含水量增大有關(guān)；動(dòng)彈性模量是靜彈性模量的10～100倍，這主要與三軸應(yīng)變?cè)囼?yàn)的應(yīng)變數(shù)量級(jí)與聲波振動(dòng)產(chǎn)生的應(yīng)變數(shù)量級(jí)的差異大有關(guān)。采用本論文實(shí)驗(yàn)測(cè)量的數(shù)據(jù)分別建立了雙復(fù)合參數(shù)-聲速和孔隙度-聲速經(jīng)驗(yàn)公式，分析結(jié)果表明雙復(fù)合參數(shù)-聲速公式聲速預(yù)報(bào)誤差約是孔隙度-聲速公式的1/4，表明雙復(fù)合參數(shù)-聲速公式更加有效。

南海；海底沉積物；聲速；三軸應(yīng)變；經(jīng)驗(yàn)公式

1　引言

海底是實(shí)際載體的承受界面，在海洋工程中是前期勘探對(duì)象之一。對(duì)海底沉積物力學(xué)性質(zhì)的研究開始于Stoll等[1—2]。在垂直加載的受控應(yīng)力-應(yīng)變過程中，聲速隨應(yīng)變?cè)黾訒?huì)有不同的響應(yīng)，其物理機(jī)制已經(jīng)引起研究人員的注意[3—4]。海底沉積物物理性質(zhì)對(duì)聲速的影響，Hamilton和Bachman[5—6]、Anderson[7]和盧博等[8]學(xué)者都獨(dú)立得到了海底沉積物孔隙度和顆粒度等聲速的經(jīng)驗(yàn)公式。這些結(jié)果基本沒有考慮三軸應(yīng)力-應(yīng)變的影響，但實(shí)際海底沉積物使處于三軸應(yīng)力-應(yīng)變狀態(tài)下的，因此本文開展垂直加載的受控三軸應(yīng)力-應(yīng)變過程的聲學(xué)測(cè)量，探討三軸應(yīng)變下聲速與物理力學(xué)性質(zhì)之間關(guān)系的變化趨勢(shì)；另外三軸應(yīng)變-聲學(xué)同步測(cè)量可以獲得沉積物的體積壓縮模量等數(shù)據(jù)，滿足利用作者前期論文[9—10]提出的雙復(fù)合參數(shù)-聲速理論模型建立聲速經(jīng)驗(yàn)公式的需要，因此采用相同的實(shí)驗(yàn)數(shù)據(jù)分別建立雙復(fù)合參數(shù)-聲速和孔隙度-聲速經(jīng)驗(yàn)公式，將可探討它們的差別。

2　實(shí)驗(yàn)方法和數(shù)值計(jì)算方法

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

本文實(shí)驗(yàn)研究的受控軸向應(yīng)力-應(yīng)變過程是加圍壓、不排水狀態(tài)下的壓縮試驗(yàn)，附加一對(duì)100 kHz聲學(xué)換能器在被測(cè)量的沉積物試樣兩端，組成一個(gè)軸向應(yīng)力-應(yīng)變/聲學(xué)測(cè)量系統(tǒng)；上端聲學(xué)換能器增加浮力結(jié)構(gòu)使其在水中成懸浮狀態(tài)，避免聲學(xué)換能器自重影響測(cè)試結(jié)果，見圖1。測(cè)量與裝置是在文獻(xiàn)[3]的基礎(chǔ)上改進(jìn)，采用特別定制的土工三軸儀(江蘇省南京土壤儀器廠造)可以準(zhǔn)確自動(dòng)加載圍壓、軸向壓力與測(cè)量壓縮變形?？杉虞d的最大圍壓為2.5 MPa。聲學(xué)測(cè)量系統(tǒng)由DB-4型聲波參數(shù)測(cè)定儀和個(gè)人計(jì)算機(jī)組成，用來實(shí)現(xiàn)窄脈沖聲波的發(fā)射和接收，進(jìn)行信號(hào)放大、數(shù)據(jù)采集和數(shù)據(jù)儲(chǔ)存。

為模擬沉積物海底承載過程，本試驗(yàn)選擇不固結(jié)不排水剪切的UU試驗(yàn)方法，用三軸試驗(yàn)測(cè)定土的強(qiáng)度參數(shù)。樣品的安裝和剪切試驗(yàn)過程按國(guó)家土工試驗(yàn)方法標(biāo)準(zhǔn)[12]要求進(jìn)行。

根據(jù)沉積物取樣時(shí)站位水深及試驗(yàn)條件，試驗(yàn)所用三級(jí)圍壓分別為300 kPa、500 kPa、700 kPa。樣品初始直徑68 mm，經(jīng)過削土后的試樣直徑為38 mm，高度125 mm，高度與直徑之比為2.0～3.5。

圖1　三軸應(yīng)力-應(yīng)變/聲學(xué)測(cè)量系統(tǒng)Fig.1　Schematic diagram of acoustic measurement system in controlled three axis stress-strain

試驗(yàn)所用剪切應(yīng)變速率取0.9 mm/min。試驗(yàn)時(shí)每產(chǎn)生0.5 mm變形值，記錄一次測(cè)力計(jì)讀數(shù)和軸向變形值，同時(shí)測(cè)量一次聲波走時(shí)和聲波首波幅值，測(cè)量過程由計(jì)算機(jī)自動(dòng)記錄數(shù)據(jù)。當(dāng)測(cè)力計(jì)讀數(shù)達(dá)到穩(wěn)定或出現(xiàn)倒退時(shí)，本級(jí)圍壓下的剪切試驗(yàn)結(jié)束。換試樣按相同的方法進(jìn)行第二級(jí)、第三級(jí)圍壓下的剪切試驗(yàn)。

2.2數(shù)值計(jì)算方法

沉積物聲波速度為：

(1)

式中，l為不斷受壓的樣品的高度，單位：mm，該高度在應(yīng)變過程中是不斷變化的；t為聲波通過試樣的走時(shí)，單位：μs。

試樣的軸向變形、校正面積和軸向應(yīng)力計(jì)算采用土工三軸儀制造公司提供的軟件自動(dòng)完成。

由于部分試樣為水飽和的黏性土，在實(shí)測(cè)時(shí)，試樣的高度和直徑均有稍微變化，但在允許范圍之內(nèi)。

2.3實(shí)驗(yàn)樣品的制備

沉積物樣品來自南海某海域大陸架。由于海底沉積物的物理力學(xué)性質(zhì)與巖石和泥巖截然不同，特別是在抵抗單軸應(yīng)力的應(yīng)變過程中，抗剪強(qiáng)度比巖石和泥巖要低很多[13]。經(jīng)初步試驗(yàn)發(fā)現(xiàn)，海底沉積物樣品不適合做單一試樣的多級(jí)加荷試驗(yàn)。根據(jù)試驗(yàn)標(biāo)準(zhǔn)，把沉積物試驗(yàn)樣品分為6組，每組3個(gè)樣品，共18個(gè)樣品。每一組樣品中的3個(gè)沉積物試驗(yàn)樣品均取自同一海域，海底埋深相近，物理性質(zhì)相似。

表1是6組樣品沉積的物理參數(shù)、站位水深及試驗(yàn)圍壓。

表1　試驗(yàn)樣品分組編號(hào)、物理性質(zhì)、站位水深及試驗(yàn)圍壓

續(xù)表1

注：NONE表示沒有獲得該項(xiàng)實(shí)驗(yàn)數(shù)據(jù)的可靠值。

2.4聲速測(cè)量數(shù)據(jù)預(yù)處理

三軸試驗(yàn)過程中的聲速測(cè)量是在軸向壓縮條件下進(jìn)行的，樣品的長(zhǎng)度不斷地縮短，聲速也在不斷地變化。如樣品GD-10的聲速變化曲線見圖2。

圖2　樣品GD-10的聲速變化曲線Fig.2　Velocity curve of sample GD-10 sediments with deformation

把整個(gè)壓縮過程中的聲速取平均值，各組沉積物樣品在剪切壓縮過程中的平均聲速見表2。

表2　沉積物樣品平均聲速

續(xù)表2

3　實(shí)驗(yàn)結(jié)果和分析

3.1動(dòng)彈性模量與聲速關(guān)系分析

在聲波幅度較小時(shí)，平面縱波波速公式為：

(2)

式中，M為等效彈性模量。M=K+4/3 μ，K、μ和ρ分別為介質(zhì)的體積模量、切變模量和密度。

在三軸試驗(yàn)的壓縮過程中，通過聲學(xué)測(cè)量得到沉積物樣品在不同壓縮階段的聲速數(shù)據(jù)。然后根據(jù)式(2)得到沉積物在三軸壓縮變形過程中的動(dòng)彈性模量M。表3為各沉積物樣品的平均動(dòng)彈性模量。

表3　各沉積物樣品的平均動(dòng)彈性模量

以沉積物樣品的孔隙度和平均動(dòng)彈性模量進(jìn)行回歸分析，得到擬合曲線，見圖3。

圖3　沉積物孔隙度與動(dòng)彈性模量關(guān)系Fig.3　Relationship between the porosity and the dynamic modulus of the sediments

由圖3可知，沉積物的孔隙度和動(dòng)彈性模量呈良好的負(fù)相關(guān)關(guān)系，這與Hamilton[14]的結(jié)果一致。

3.2靜彈性模量與聲速關(guān)系分析

在三軸試驗(yàn)中，主應(yīng)力差σ1-σ3為已知測(cè)量值。試樣為圓柱體時(shí)泊松比ν表示為：

(3)

式中，ε3為側(cè)向應(yīng)變，ε1為軸向應(yīng)變。

沉積物樣品的體積應(yīng)變?chǔ)舦為：

(4)

根據(jù)三軸試驗(yàn)原理[11]，沉積物樣品在不同壓縮階段的靜彈性模量為:

(5)

應(yīng)用Hamilton的結(jié)果[15]，沉積物的站位水深對(duì)應(yīng)的泊松比如表4。

表4　沉積物水深及泊松比

由式(5)可知，沉積物靜彈性模量隨著壓縮過程中應(yīng)變的增大而減小，且不同圍壓σ3對(duì)沉積物靜彈性模量影響較大。聯(lián)立式(3)～(5)可以計(jì)算出沉積物樣品在三軸壓縮過程對(duì)應(yīng)的靜彈性模量E。

以下分別對(duì)不同圍壓、不同軸向變形的靜彈性模量E進(jìn)行對(duì)比分析。沉積物在軸向壓縮過程中，靜彈性模量是不斷地變化的，故取各沉積物樣品的平均靜彈性模量與孔隙度進(jìn)行數(shù)據(jù)擬合。圖4a～c為不同圍壓下沉積物靜彈性模量與孔隙度關(guān)系。

圖4　不同軸向變形下靜彈性模量與孔隙度關(guān)系Fig.4　Relationship between the static elastic modulus and the porosity under different axial deformation

由圖4a～c可知，在圍壓相同時(shí)，沉積物在不同的軸向變形下靜彈性模量值有明顯的差異，且靜彈性模量隨孔隙度的增加呈增大的趨勢(shì)。其主要原因是飽和沉積物剛性與含水量密切相關(guān)。

平均靜彈性模量和平均聲速的相關(guān)性曲線，見圖5。由圖5可知，沉積物的聲速先隨著靜彈性模量的增大而減小，在靜彈性模量4×107～5×107N/m2之間聲速取得最小值，爾后聲速隨著靜彈性模量的增大迅速增大。

圖5　沉積物靜彈性模量與聲速關(guān)系Fig.5　Relationship between the static elastic modulus and the acoustic velocity of the sediments

通過對(duì)比發(fā)現(xiàn)動(dòng)彈性模量是靜彈性模量的10～100倍。與Mavko等[16]的結(jié)果比對(duì)，本文的動(dòng)、靜彈性模量值在合理區(qū)間內(nèi)，它們數(shù)量級(jí)差異也與Mavko Gary的結(jié)果一致。

作者前期論文[9—10]已經(jīng)提出了雙復(fù)合參數(shù)-聲速理論模型，三軸應(yīng)變-聲學(xué)同步測(cè)量可以獲得沉積物的體壓縮模量等數(shù)據(jù)，滿足建立雙復(fù)合參數(shù)-聲速經(jīng)驗(yàn)公式的需要，第4節(jié)對(duì)此加以討論。

4　海底沉積物雙復(fù)合參數(shù)-聲速經(jīng)驗(yàn)公式

4.1固結(jié)等效彈性模量和固結(jié)密度計(jì)算

飽和海底沉積物濕密度[9]ρn為:

(6)

式中，n、ρw分別為孔隙度和海水的密度，ρs表示排水固結(jié)后砂、粉砂、黏土3種成分組成的沉積物的密度。

排水固結(jié)后以砂、粉砂、黏土組成的沉積物等效彈性模量[9]Ms為：

(7)

式中，M1、M2、M3分別表示沉積物的砂、粉砂、黏土的等效彈性模量，θ1、θ2、θ3分別表示砂、粉砂、黏土的體積分?jǐn)?shù)。

各成分的體積分?jǐn)?shù)和質(zhì)量分?jǐn)?shù)可以通過以下理論關(guān)系換算：

(8)

式中，參數(shù)R1、R2、R3分別表示表示砂、粉砂、黏土的質(zhì)量分?jǐn)?shù)，ρs1、ρs2、ρs3分別表示表示砂、粉砂、黏土的礦物顆粒密度。

砂、粉砂、黏土的等效彈性模量和礦物顆粒密度參考Hamilton[14]給出的大陸臺(tái)階(大陸架和大陸坡)沉積物對(duì)應(yīng)參數(shù)的計(jì)算值。各參數(shù)通過單位換算后的數(shù)值如表5。

表5　沉積物3種成分的等效彈性模量和顆粒密度

聯(lián)立式(6)～(8)可以計(jì)算出本實(shí)驗(yàn)所用的各海底沉積物的固結(jié)密度ρs和固結(jié)等效彈性模量Ms，其結(jié)果如表6(數(shù)據(jù)缺失或重復(fù)樣品未列出)。

表6　沉積物固結(jié)密度和等效固結(jié)彈性模量

等效固結(jié)模量Ms考慮的是沉積物固結(jié)后不含水無空隙的極限狀態(tài)。由式(9)可以換算出飽和沉積物的等效模量Mn。

(9)

式中，Mn為飽和沉積物的等效彈性模量，Mw為孔隙水的等效彈性模量，Ms為固結(jié)模量。

孔隙水的等效彈性模量[5]取Mw=2.374 0×109N/m2，算得飽和沉積物的等效彈性模量如表7。

表7　沉積物飽和等效彈性模量

本文計(jì)算得出的沉積物的等效模量值趨近于Hamilton和Bachman[5]的大陸臺(tái)階(大陸架和大陸坡)環(huán)境的模量值，在3.14×109～4.51×109N/m2之間。

4.2沉積物雙復(fù)合參數(shù)-聲速經(jīng)驗(yàn)公式的建立與分析

利用表7的數(shù)據(jù)、等效固結(jié)彈性模量Ms和等效固結(jié)體積模量Ks的換算關(guān)系，可得沉積物樣品的等效體積模量Ks如表8。

表8　沉積物固結(jié)體積模量

由理論關(guān)系Kn=nKw+(1-n)Ks，可求得飽和沉積物的體積模量Kn如表9。

表9　飽和沉積物體積模量

文獻(xiàn)[9]得到了海底沉積物聲速多項(xiàng)式模型，本文取其前4項(xiàng)得到：

CP=CPr(1+a1x+a2x2+a3x3)×

(1+b1y+b2y2+b3y3)，

(10)

式中，CPr=(Ms/ρs)1/2，x=n(1-Kw/Ks)，y=n(1-ρw/ρs)，a1、a2、a3和b1、b2、b3為待定多項(xiàng)式系數(shù)。

利用表6～表9數(shù)據(jù)，可以確定多項(xiàng)式(10)的待定系數(shù)，得到雙復(fù)合參數(shù)-聲速經(jīng)驗(yàn)公式為：

CP=2 860.9(1-0.667 8x+3.956 5x2-3.197 0x3)

×(1+1.513 7y+1.711 1y2-0.696 6y3).

(11)

根據(jù)實(shí)驗(yàn)數(shù)據(jù)，復(fù)合參數(shù)x=n(1-Kw/Ks)的取值范圍為0.147 9～0.419 1，復(fù)合參數(shù)y=n(1-ρw/ρs)的取值范圍為0.190 0～0.538 4，式(11)所刻畫的聲速曲面見圖6。圖中從藍(lán)色到黃色代表沉積物的聲速的增大過程。

圖6　沉積物聲速曲面Fig.6　Surface of velocity of sediments

利用表1數(shù)據(jù)，可以建立孔隙度-聲速回歸公式為:

CP=4 093.85-129.93n+2.231 2n2-0.013 29n3.

(12)

圖6中，沉積物聲速隨著復(fù)合參數(shù)x和復(fù)合參數(shù)y的增大而變小，并且在(x，y)≈(0.147 9，0.19)處聲速取得最大值。這和式(12)的沉積物聲速隨著孔隙度的增大而減小的結(jié)論是一致的。

4.3預(yù)報(bào)聲速與測(cè)量聲速誤差比較

把復(fù)合參數(shù)x=n(1-Kw/Ks)，y=n(1-ρw/ρs)數(shù)值代入曲面方程式(11)，及把孔隙度數(shù)值代入式(12)，可以分別計(jì)算出兩個(gè)公式的聲速預(yù)報(bào)值。與沉積物樣品測(cè)量聲速值比較，對(duì)應(yīng)的誤差如表10和表11。

表10　雙復(fù)合參數(shù)-聲速公式計(jì)算沉積物聲速預(yù)報(bào)值與誤差

通過對(duì)比可知，式(11)的雙復(fù)合參數(shù)-聲速公式

的絕對(duì)誤差和相對(duì)誤差最大值分別為：38.42 m/s、0.026；式(12)的孔隙度-聲速公式的絕對(duì)誤差和相對(duì)誤差最大值分別為：151.32 m/s、0.094 45?？梢娗罢呗曀兕A(yù)報(bào)誤差約是后者的1/4。

表11　孔隙度-聲速公式計(jì)算沉積物聲速預(yù)報(bào)值與誤差

5　結(jié)論

圍壓相同時(shí)軸向應(yīng)變?cè)龃箪o彈性模量也增大；沉積物動(dòng)彈性模量與孔隙度的相關(guān)性很高，相關(guān)系數(shù)達(dá)0.97，動(dòng)彈性模量是靜彈性模量的10～100倍。在應(yīng)變加載到結(jié)束過程聲速不斷變化；本次試驗(yàn)樣品的平均聲速隨著平均靜彈性模量的增加，由大變小又由小變大，存在聲速最小值。分別建立了雙復(fù)合參數(shù)-聲速和孔隙度-聲速經(jīng)驗(yàn)公式，分析結(jié)果表明雙復(fù)合參數(shù)-聲速經(jīng)驗(yàn)公式預(yù)報(bào)聲速的誤差遠(yuǎn)小于孔隙度-聲速經(jīng)驗(yàn)公式預(yù)報(bào)聲速的誤差。

[1]Stoll R D. Measuring parameters that control propagation in granular sediments near the sea floor[C]//Proceedings of the 135th Meeting of the Acoustical. Seattle: Society of America, 1998: 1611-1612.

[2]Stoll R D, Sun Yuefeng, Bitte I. Seafloor properties from penetrometer tests[J]. IEEE Journal of Oceanic Engineering, 2007, 32(1): 57-63.

[3]Lu Bo, Li Ganxian, Huang Shaojian, et al. The three results of sound velocity variance in sediments during controlled stress-strain testing from the South China Sea[J]. Marine Georesources and Geotechnology, 2002, 20(3): 177-186.

[4]Chen Minpen, Shieh Y T, Chyan J M. Acoustic and physical properties of surface sediments in northern Taiwan strait[J]. Acta Oceanographica Taiwanica, 1988(21): 92-118.

[5]Hamilton E L, Bachman R T. Sound velocity and related properties of marine sediments[J]. Journal of Acoustical Society of America, 1982, 72(6): 1891-1904.

[6]Hamilton E L. Geoacoustic modeling of the sea floor[J]. Journal of Acoustical Society of America, 1980, 68(5): 1313-1340.

[7]Anderson R S. Statistical correlation of physical properties and sound velocity in sediments[M]//Hampton L. Physics of Sound in Marine Sediments. US: Springer, 1974: 48l-5l8.

[8]盧博, 李趕先, 孫東懷, 等. 中國(guó)東南近海海底沉積物聲學(xué)物理性質(zhì)及其相關(guān)關(guān)系[J]. 熱帶海洋學(xué)報(bào), 2006, 25(2): 12-17.

Lu Bo, Li Ganxian, Sun Donghuai, et al. Acoustic-physical properties of seafloor sediments from nearshore southeast China and their correlations[J]. Journal of Tropical Oceanography, 2006, 25(2): 12-17.

[9]龍建軍, 李趕先. 海底沉積物聲速與物理性質(zhì)的理論關(guān)系[J]. 聲學(xué)學(xué)報(bào), 2015, 40(3): 462-468.

Long Jianjun, Li Ganxian. Theoretical relations between sound velocity and physical-mechanical properties for seafloor sediments[J]. Acta Acustica, 2015, 40(3): 462-468.

[10]龍建軍, 李趕先, 鄒大鵬. 海底沉積物縱波波速與物理-力學(xué)性質(zhì)的理論關(guān)系和比較[J]. 海洋學(xué)報(bào), 2014, 36(7): 111-117.

Long Jianjun, Li Ganxian, Zou Dapeng. Theoretical relations of longitudinal wave velocity and physical-mechanical properties for seafloor sediments and comparison[J]. Haiyang Xuebao, 2014, 36(7): 111-117.

[11]朱思哲, 劉虔, 包承綱, 等. 三軸試驗(yàn)原理與應(yīng)用技術(shù)[M]. 北京: 中國(guó)電力出版社, 2003.

Zhu Sizhe, Liu Wen, Bao Chenggang, et al. Three Axis Test Principle and Application Technology[M]. Beijing: China Power Press, 2003.

[12]國(guó)家質(zhì)量技術(shù)監(jiān)督局, 中華人民共和國(guó)建設(shè)部. GB/T 50123-1999, 土工試驗(yàn)方法標(biāo)準(zhǔn)[S]. 北京: 中國(guó)計(jì)劃出版社, 2007.

The State Bureau of Quality and Technical Supervision,Ministry of Construction of the People's Republic of China. Soil Test Method Standard (GB/T50123-1999)[S]. Beijing: China Planning Publishing House,2007.

[13]劉寄能, 梁元博, 盧博, 等. 南海兩海區(qū)海底沉積物的聲速和應(yīng)力-應(yīng)變特征[J]. 南海海洋科學(xué)刊, 1989(9): 87-94.

Liu Jineng, Liang Yuanbo, Lu Bo, et al. Acoustic velocity and stress-strain characteristics of seabed sediments in two areas of South China Sea[J]. The Volume of South China Sea Marine Science, 1989(9): 87-94.

[14]Hamilton E L. Elastic properties of marine sediments[J]. Journal of Geophysical Research, 1971, 76(2): 579-604.

[15]Hamilton E L. Vp/Vs and Poisson's ratios in marine sediments and rocks[J]. Journal of Acoustical Society of America, 1979, 66(4): 1093-1101.

[16]Mavko G, Mukerji T, Dvorkin J. The rock physics handbook: tools for seismic analysis in porous media[M]. London: Cambridge University Press, 2003.

索安寧，王鵬，袁道偉，等. 基于高空間分辨率衛(wèi)星遙感影像的圍填海存量資源監(jiān)測(cè)與評(píng)估研究——以營(yíng)口市南部海岸為例[J]. 海洋學(xué)報(bào)， 2016， 38(9):54-63，doi:10.3969/j.issn.0253-4193.2016.09.006

Suo Anning, Wang Peng, Yuan Daowei, et al. Study on monitoring and analysis of existing sea reclamation resource based on high resolution satellite remote sensing imagery: A case in south coast of Yingkou[J]. Haiyang Xuebao， 2016， 38(9):54-63， doi：10.3969/j.issn.0253-4193.2016.09.006

Relationship between sound velocity and physical mechanical properties of seafloor sediments in controlled three axis stress-strain

Long Jianjun1， Zhou Huajian1， Li Ganxian2， Chen Cong1

(1.SchoolofElectromechanicalEngineering，GuangdongUniversityofTechnology，Guangzhou510006，China;2.KeyLaboratoryofMarginalSeaGeology，SouthChinaSeaInstituteofOceanology，ChineseAcademyofSciences，Guangzhou510301，China)

The sediment cores， from shallow seafloor in depth of 100-400 meters (about 2 m depth buried) in an area of the South China Sea， have been performed a simultaneous triaxial strain-acoustic measurement near seabed water pressure. The results show that the sediment longitudinal wave velocity has two characteristics: (1) the longitudinal wave velocity in sediments keeps changing in the process of strain applied from start to end; (2) the average of the speed of sound changes from large to small and then from small to large with the average of the static elastic modulus increasing， and there is a minimum during the course of the average of the speed of sound change. These results are related to physical and mechanical properties and confining pressure of the shallow seafloor sediments， and the changes of the state in which particles of these sediments combine together. In addition， sediment dynamic elastic modulus and porosity has a good negative correlation， and that results from increases of the water content with increases of porosity. Dynamic elastic modulus is 10 to 100 times the static elastic modulus， since the large difference between the strain magnitude resulted from triaxial strain-test and the one generated by wave vibrations is primarily related to that. Two empirical formulas of dual complex parameters-velocity and porosity-velocity were established with the experimentally measured data in this paper. The analysis of results shows that the error in sound speed predicted by the formula of composite double parameters-velocity is approximately a quarter of that predicted by the formula of the porosity-velocity， and it shows that the formula of composite double parameters-velocity is more effective．

South China Sea; seafloor sediment; sound velocity; three axis stress-strain; empirical formula

2015-12-10；

2016-03-26。

國(guó)家自然科學(xué)基金(41176034, 41476028)。

龍建軍(1960—)，男，廣西壯族自治區(qū)玉林市人，博士，副教授，主要從事海洋沉積聲學(xué)、海洋聲學(xué)測(cè)量技術(shù)研究。E-mail:longjj@gdut.edu.cn

P733.23

A

0253-4193(2016)09-0046-08