• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于Buckingham模型和Biot-Stoll模型的南沙海域沉積物聲速分布特征

    2016-07-04 03:22:05王景強郭常升劉保華陳文景侯正瑜韓國忠
    地球?qū)W報 2016年3期
    關(guān)鍵詞:聲速聲學沉積物

    王景強, 郭常升, 劉保華, 陳文景, 侯正瑜, 韓國忠

    1)國家海洋局第一海洋研究所, 海洋沉積與環(huán)境地質(zhì)國家海洋局重點實驗室, 山東青島 266061; 2)中國科學院海洋研究所, 中國科學院海洋地質(zhì)與環(huán)境重點實驗室, 山東青島 266071; 3)國家深?;毓芾碇行? 山東青島 266061; 4)青島海洋科學與技術(shù)國家實驗室, 海洋地質(zhì)過程與環(huán)境功能實驗室, 山東青島 266061; 5)中國科學院大學, 北京 100049

    ?

    基于Buckingham模型和Biot-Stoll模型的南沙海域沉積物聲速分布特征

    王景強1, 4), 郭常升2, 4)*, 劉保華3, 4), 陳文景2, 5), 侯正瑜2, 5), 韓國忠1, 4)

    1)國家海洋局第一海洋研究所, 海洋沉積與環(huán)境地質(zhì)國家海洋局重點實驗室, 山東青島 266061; 2)中國科學院海洋研究所, 中國科學院海洋地質(zhì)與環(huán)境重點實驗室, 山東青島 266071; 3)國家深海基地管理中心, 山東青島 266061; 4)青島海洋科學與技術(shù)國家實驗室, 海洋地質(zhì)過程與環(huán)境功能實驗室, 山東青島 266061; 5)中國科學院大學, 北京 100049

    摘 要:為了選擇合理的南海南部海底沉積物的聲速反演模型, 分別利用Wood方程、Buckingham模型、Biot-Stoll模型等聲波傳播理論模型計算了沉積物的縱波聲速, 結(jié)合沉積物實測聲速和物理參數(shù)等, 對比了模型計算聲速與實測聲速的變化特征。聲速實測值與模型計算值的對比表明, Buckingham模型計算值與實測值的平均偏差最小。研究區(qū)實測聲速比和模型計算聲速比的分布表明, 聲速的變化特征與沉積環(huán)境具有相關(guān)性, 聲速模型計算值與實測值具有相似的變化特征; 模型計算聲速與實測聲速的垂向變化規(guī)律不相同,模型計算聲速隨深度基本不變, 而實測聲速在陸架、陸坡和海槽等沉積環(huán)境中呈現(xiàn)出多種垂向分布特征。研究結(jié)果認為, 側(cè)重于沉積物孔隙性的Biot-Stoll理論僅適用于陸架區(qū)粗顆粒沉積物, 兼具流體和彈性固體雙重性質(zhì)的Buckingham理論具有更廣的適用范圍。該研究對獲取南海南部海域沉積物聲速分布特征有重要意義, 可為多波束測深或淺地層精細解釋等提供基礎(chǔ)資料。

    關(guān)鍵詞:海底沉積物; Buckingham模型; Biot-Stoll模型; 聲速分布特征; 聲速比

    本文由海洋公益性行業(yè)科研專項經(jīng)費(編號: 201405032; 200905025)和國家自然科學基金重點項目(編號: 41330965)聯(lián)合資助。

    海底松散沉積物通常被作為一種液固兩相介質(zhì), 主要由砂、粉砂、黏土和孔隙流體(海水)等物質(zhì)組成。海底表層沉積物構(gòu)成了海水的下邊界和沉積物的上邊界, 從海底表層到深層, 沉積物由未固結(jié)松散介質(zhì)變?yōu)榘牍探Y(jié), 乃至固結(jié)沉積物介質(zhì)(鄒大鵬等, 2012a)。相對于結(jié)構(gòu)密實的深層沉積物, 作為海水與深層沉積物分界面的表層松散沉積物(1 m埋深范圍內(nèi)的沉積層), 具有更復雜的聲學特性和獨特的聲學指示意義, 其聲速的分布特征是海域聲場環(huán)境的重要部分, 也是進行海洋多波束或淺地層精細解釋的基礎(chǔ)資料。目前獲取海底沉積物聲學參數(shù)的方法, 主要包括理論模型預測、經(jīng)驗方程預報和聲學參數(shù)測量等。許多研究學者利用取樣測量方法(Hamilton, 1985; Lu et al., 2010)和原位測量技術(shù)(Briggs and Richardson, 1997; Liu et al., 2013; Hou et al., 2014)進行了大量的海底沉積物聲學測量工作,建立了聲速與物理參數(shù)之間的經(jīng)驗公式, 探討了聲學性質(zhì)與物理性質(zhì)之間的關(guān)系, 這些聲學實測工作為驗證聲波傳播理論模型的有效性提供了基礎(chǔ)資料。目前建立的聲波傳播理論模型主要包括Biot-Stoll模型(Biot, 1956a, b; Stoll, 1989; Stoll and Bautista, 1998; Chotiros and Isakson, 2004; 朱祖揚等, 2013)、Wood方程(Jackson and Richardson, 2007) 和Buckingham模型(Buckingham, 1997, 1998)等。聲波傳播理論模型的有效性得到了聲學實測的分析和驗證, Williams等(2002)基于SAX99聲學試驗, 比較了理論模型預測值與海底沉積物聲學實測值, 結(jié)果表明Biot理論和等效流體近似理論能夠匹配聲速的頻散特征, 但不能準確表述衰減隨頻率的變化趨勢,而衰減和頻率的關(guān)系與Buckingham理論預測結(jié)果能夠?qū)崿F(xiàn)最佳擬合。Buckingham(2005)系統(tǒng)地比較了聲波速度和衰減的模型預測值與實測值, 認為高頻段時衰減與頻率之間具有近似的線性關(guān)系, 而低頻段的聲衰減與頻率呈近似平方關(guān)系。這些聲學測量和模型預測的對比基于實驗室沉積物模型, 對聲速和聲衰減頻散特征進行比較, 缺乏與沉積物所處海域的沉積環(huán)境相聯(lián)系。因此, 需要分析聲波傳播模型在復雜沉積環(huán)境聲學參數(shù)預測中的適用性, 從而為海底沉積物選擇更有效的聲學反演模型, 有利于實現(xiàn)快速反演復雜環(huán)境沉積物的聲學特性。

    南海南部海域覆蓋大陸架、大陸坡及南海海槽等多種地貌單元, 其沉積物具有復雜的分層結(jié)構(gòu),從而影響了物理性質(zhì)和聲學特性的分布特征。因此,需要結(jié)合實際聲學測量分析復雜沉積環(huán)境對縱波聲速分布特征的影響, 從而選擇合理的聲速計算預測模型。張福生等(2001)根據(jù)沉積物聲速與底層海水聲速的比值(聲速比R)將該海域的表層沉積物分為高聲速和低聲速兩類, 但沒有分析聲速的空間分布特征。針對該研究區(qū)Biot-Stoll模型參數(shù)的選擇, 陳靜等(2012)分析了Stoll參數(shù)和Schock參數(shù)的模擬匹配效果, 但缺乏與其他理論模型的對比。本文將基于研究區(qū)沉積物的聲學測量資料, 分別利用Wood方程、Biot-Stoll模型和Buckingham模型等聲波傳播理論模型進行縱波聲速計算, 并與實際測量獲取的沉積物聲速值進行對比, 分析不同地貌單元中模型預測聲速的有效性。根據(jù)模型計算聲速比和實測聲速比的分布特征的對比, 探討理論模型對復雜沉積環(huán)境聲速分布特征的適用性。

    1 海底沉積物聲波傳播理論

    海底沉積物聲波傳播理論主要包括流體理論、彈性理論及多孔彈性理論三種類型, 分別把海底沉積物看成流體、彈性固體和多孔彈性介質(zhì)。上述三種類型的聲波傳播理論均可歸入“等效介質(zhì)”理論的范疇。由于聲波波長遠大于沉積物顆粒的尺寸,可以把海底沉積物看作連續(xù)介質(zhì), 該連續(xù)介質(zhì)的物理性質(zhì)參數(shù)可以通過求取一定范圍內(nèi)的不同的微觀測量值的平均值而獲得。如果聲波頻率很高, 聲波波長與顆粒的尺寸及孔隙的尺寸差不多時, 此時等效介質(zhì)理論不能正確地描述沉積物中的波動現(xiàn)象。

    1.1流體理論

    流體理論假定海底沉積物中聲波引起的應力可以用壓力場和相應的波動方程來描述, 其所需要的主要參數(shù)包括聲速、聲衰減及沉積物和海水的密度等。該理論將海底沉積物看成連續(xù)介質(zhì), 把一定體積內(nèi)包含許多沉積物顆粒的物理參數(shù)平均值作為沉積物物理參數(shù)的數(shù)值, 對一定體積的沉積物單元體的等效密度和等效體積模量建立方程, 并得到沉積物聲速的計算公式, 即Wood方程(Jackson and Richardson, 2007)。該方程忽略了沉積物顆粒間相互接觸產(chǎn)生的力, 也沒有考慮孔隙流體與介質(zhì)骨架間的相對運動, 即忽略了慣性力和黏滯力, 因此其具有很大的聲速預測誤差和應用的局限性, 其對不太注重沉積物顆粒接觸力和密度差異影響的實驗室模型測量可以做出合理的解釋。沉積物單元體的等效密度ρ0和體積模量Kb的表達式分別為:

    式中, N為孔隙度, ρw為水密度, ρg為沉積物顆粒密度; Kw為水的體積模量, Kg為沉積物顆粒的體積模量。根據(jù)式(1)、(2)和(3)可得出沉積物聲速c0計算公式:

    1.2彈性理論

    聲波傳播彈性理論考慮了固體顆粒的彈性性質(zhì)和流體理論所未考慮的剪切力的效應。對于黏土、粉砂和砂質(zhì)沉積物, 剪切波傳播速度非常低, 遠低于海水的聲速。Gassmann方程(Berryman, 1999)基于彈性理論, 考慮了流體飽和多孔介質(zhì)的固體顆粒骨架的彈性性質(zhì), 但忽略了孔隙介質(zhì)的孔隙流體和介質(zhì)骨架間的相對運動, 不適于討論滲透系數(shù)較大的沉積物的聲波傳播特性。基于彈性理論和流體理論, Buckingham模型假定沉積物顆粒相互接觸但沒有膠結(jié), 顆粒間存在一種類似于黏滯作用的力, 因此該理論兼具流體和彈性固體的雙重性質(zhì)(Buckingham, 1997, 1998)。該模型以顆粒粒徑為主要參數(shù), 建立了單一頻率下近似流體沉積物的聲速cp預測方程:

    式中, μ0為縱波摩擦剛性常數(shù), 取值為2×109Pa; ug為沉積物顆粒粒徑; u0為參照粒徑, 取值為1 000 μm; ρ0和c0與Wood方程相同。

    1.3多孔彈性理論

    Biot理論是一種經(jīng)典的多孔彈性理論(Biot, 1956a, b), 同時考慮介質(zhì)的孔隙性與彈性, 用來描述彈性波在固體和流體應力相互作用下的傳播特征。Stoll(1989)將Biot理論應用于海底沉積物介質(zhì)中的聲速和聲衰減計算, 建立了Biot-Stoll理論。該理論認為固體顆粒構(gòu)成了與孔隙流體相耦合的彈性骨架, 流體中的聲波與骨架中的縱波和剪切波混合在一起, 構(gòu)成了三種波: 快縱波、慢縱波和剪切波。該理論結(jié)合沉積物顆粒間膠結(jié)差、骨架模量低的特點, 認為骨架耗散是聲波能量衰減的重要原因, 骨架耗散的大小與頻率無關(guān), 而流體黏滯耗散則隨頻率的變化而變化。理論本構(gòu)方程涉及的13個參數(shù)通過測試或者曲線擬合的方法進行確定。Biot理論彈性模量、附加彈性模量、和復彈性模量可利用孔隙度N, 骨架體積模量, 骨架剪切模量孔隙流體體積模量Kw和顆粒體積模量Kg表示為:

    式中, Kw, Kg為實數(shù)。和為復數(shù), = K0(1 + iδk) , = μ0(1 + i δμ) , δk為體積耗散系數(shù), δμ為剪切耗散系數(shù)。

    根據(jù)Biot-Stoll模型理論, 頻率域的縱波傳播方程可以寫成公式:

    式中, ρ為沉積物介質(zhì)密度, ρw為孔隙流體密度;t為復波數(shù), ?為黏滯系數(shù), κ為滲透系數(shù), ω為角頻率, m為有效密度, m=αρw/N, α為孔隙彎曲因子, F為高頻校正因子。其中, 復波數(shù)k=ω/ cb+ i α , 可以得到聲速cb的計算公式:

    其中, Re(x)為復數(shù)x的實部。

    2 試驗測量與理論計算

    2.1區(qū)域概況及數(shù)據(jù)獲取

    研究區(qū)域位于南海南部海域, 覆蓋了大陸架、南沙臺階和南沙海槽等多種海底地貌單元, 水深100~3 000 m(圖1)。陸架區(qū)沉積物一部分來源于現(xiàn)代淺海相沉積, 一部分屬于早期殘留沉積, 其沉積物組成既受到物源區(qū)沉積物自身特征的影響, 還受到后期改造作用的影響, 含量較多的粗顆粒物質(zhì)主要是高能環(huán)境中長期沉積的產(chǎn)物, 在冰后期海平面上升的過程中, 可能被細顆粒沉積物覆蓋(楊濤等, 2003), 呈現(xiàn)出粗顆粒與細顆粒沉積物相互交錯沉積的特征。水深大于1 000 m的陸坡主要為深水區(qū),地形和物源等環(huán)境因素更為復雜, 沉積物主要是細顆粒的黏土質(zhì)粉砂和粉砂質(zhì)黏土, 物源以生物源碎屑和陸源風成黏土等為主(李傲仙等, 2005; Liu et al., 2013), 由于陸坡水深變化較大且沉積速率較快(付淑清等, 2008), 沉積搬運作用強烈, 加上東北島礁區(qū)覆蓋粗顆粒珊瑚砂, 使得整個陸坡沉積物物理性質(zhì)分布的局部變化特征更為突出(李趕先和龍建軍, 2014)。南沙海槽的槽谷較為平坦, 水動力條件較弱, 沉積物埋藏后期剝蝕和搬運作用較弱, 主要受壓實作用影響, 覆蓋有較厚的沉積物, 細顆粒物質(zhì)含量高, 以粉砂質(zhì)黏土和黏土質(zhì)粉砂為主。

    圖1顯示了研究海域的位置, 站位深度136~2 843 m。沉積物柱狀樣品由科考船“科學一號”于2012年5月采用重力柱狀取樣器進行采集, 每根樣品長100 cm。樣品取回甲板后, 以15 cm為間距,采用WSD-3數(shù)字聲波儀和同軸差距法進行了聲學測量, 測量頻率為100 kHz。其中樣品長度測量精度為0.05 cm, 聲時測量精度為±0.1 μs。聲學測量完成后封存沉積物樣品, 后期在實驗室進行樣品的濕密度、孔隙度和粒度分析測量測試。其中, 沉積物的顆粒組成和粒徑利用Cilas 940L型激光粒度儀進行測試。圖2為聲學測量和物理參數(shù)測量的沉積物層位圖。通過聲學參數(shù)測量和物理參數(shù)測量, 獲得15 cm、30 cm、60 cm、75 cm、90 cm等6個埋深層位共計126組樣品的聲速和孔隙度、濕密度、粒徑、顆粒組分含量、顆粒密度等物理參數(shù)。

    圖1 研究區(qū)域示意圖Fig. 1 The study area

    由于甲板實驗室聲學測量時, 沉積物的溫度和壓力等環(huán)境參數(shù)會發(fā)生改變, 需要獲取沉積物原位的溫度壓力參數(shù)(Hamilton, 1971; 鄒大鵬等, 2012a),校正計算得到沉積物的原位聲速和近底上覆海水的聲速。根據(jù)CTD實測資料確定21個站位的底層海水溫度在1 . 5 ~ 1 6℃, 沉積物的測量溫度在12~18℃。通過海水聲速剖面、底層海水的溫度、深度和鹽度等參數(shù), 分別利用Hamilton校正法(Hamilton, 1971)和Wilson海水聲速公式(周豐年等, 2001), 精確計算校正了沉積物的聲速和近底上覆海水的聲速, 并根據(jù)校正的沉積物聲速和近底上覆海水聲速計算得到沉積物的聲速比, 其中聲速的范圍為1 402~1 627 m/s, 聲速比的范圍為0.94~1.07,該結(jié)果與張福生等(2001)在該海域測量結(jié)果的變化范圍相一致。

    2.2海底沉積物的物理參數(shù)

    由于理論模型均涉及較多的參數(shù), 參數(shù)的選取對計算結(jié)果有不同程度的影響, 所以針對研究區(qū)域選取合理的參數(shù)是必要的。Wood方程和Buckingham模型的主要參數(shù)為孔隙水密度、顆粒密度、孔隙度、粒徑、顆粒體積模量、孔隙水體積模量等。除上述參數(shù)外, Biot-Stoll模型還需要輸入更多的參數(shù), 包括滲透率、流體黏滯系數(shù)、骨架體積模量、骨架剪切模量、孔隙因子、高頻校正因子等。通過物理性質(zhì)測試獲取了顆粒密度、孔隙度、濕密度、粒徑等參數(shù), 其余參數(shù)均需要利用經(jīng)驗公式計算或從文獻資料中獲取。表1為本文詳細的參數(shù)取值, 其中滲透率、沉積物孔隙因子、顆粒體積模量和剪切模量、孔隙水體積模量等參數(shù)根據(jù)Stoll和Bautista(1998)、Buckingham(1998, 2005)模型和Jackson和Richardson(2007)的專著“High-Frequency Seafloor Acoustics”(高頻海底聲學)等文獻資料的參數(shù)變化范圍進行搜索優(yōu)化擬合確定或基于Schock(2004)公式利用實測物理參數(shù)計算得到。

    3 結(jié)果與討論

    3.1聲速實測值與聲速預測計算值比較

    圖2 沉積物柱狀樣聲學測量層位Fig. 2 The acoustic measured horizons of sediment samples

    表1 海底沉積物的物理參數(shù)Table 1 Physical parameters of marine sediments

    圖3 聲速實測值與聲速預測計算值比較Fig. 3 The comparison between the measured and the calculated sound speeds

    為比較聲速預測值與實測值, 分析聲速預測模型的適用性, 分別采用Wood方程、Buckingham模型和Biot-Stoll模型聲速預測方程進行了聲速計算,計算所采用的參數(shù)見表1。圖3顯示了聲速實測值與計算值的比較結(jié)果: 三種模型計算聲速的變化曲線與實測聲速具有大體相同的變化趨勢, 陸架區(qū)沉積物的聲速明顯高于海槽和陸坡沉積物的聲速。

    經(jīng)過差值計算, 實測值與Wood方程計算值的平均偏差為3.17%, 實測值與Biot-Stoll模型計算值的平均偏差為2.15%, 實測值與Buckingham模型計算值的平均偏差為0.87%。從圖3中可以看出, 相對于聲速實測值, Wood方程計算的聲速值偏小, Biot-Stoll模型計算的聲速值偏大。在陸坡和海槽區(qū), Biot-Stoll模型計算值為實測值值域的上限, Wood方程計算值為實測值值域的下限, Buckingham模型計算值介于二者之間。在陸架區(qū), Wood方程計算值低于實測值的值域范圍, Biot-Stoll和Buckingham模型的預測值與實測值更為接近。相對于低聲速段, Biot-Stoll模型在高聲速段的反演效果較好, 這與陳靜等(2012)在該海域的計算結(jié)果相一致。

    相對于另外兩種模型, Wood方程對沉積物聲速的預測誤差較大, 這與過往研究的結(jié)果相一致(Hamilton and Bachman, 1982), 主要原因在于該模型忽略了海底沉積物中顆粒間相互接觸的力, 還未考慮孔隙流體和顆粒骨架相對運動的黏滯力, 而這兩種力是沉積物的高頻聲學測量需要著重考慮的內(nèi)容。海底沉積物作為一種彈性介質(zhì), 具有一定的剛性, 并不是純粹的流體, 其剛性與沉積物顆粒間的接觸作用有關(guān), 因而Wood方程適合于解釋不注重顆粒接觸力和密度差異影響的沉積物的聲學特性(Jackson and Richardson, 2007)。不同于Wood方程, Buckingham模型和Biot-Stoll模型既考慮了固體顆粒骨架的彈性性質(zhì), 還考慮了孔隙流體和骨架相對運動的黏滯力對聲波傳播產(chǎn)生的影響, 因而這兩種模型更適用于表達沉積物的聲波傳播特征。陸架區(qū)的砂質(zhì)沉積物作為多孔彈性介質(zhì), 陸坡和海槽的松軟泥質(zhì)沉積物作為流體飽和多孔介質(zhì), 單一模型還不能完全解釋和準確表達聲波在多類型沉積物中的傳播特征, 基于介質(zhì)的多孔彈性效應且側(cè)重孔隙性的Biot-Stoll理論僅適用于陸架區(qū)沉積物, 而兼具流體和彈性固體雙重性質(zhì)的Buckingham理論, 具有相對較廣的適用范圍。

    3.2沉積物聲速變化特征

    沉積物聲速的空間變化特征包括橫向和垂向兩部分, 是海底聲場結(jié)構(gòu)、海域地聲模型的重要組成部分。聲速的橫向變化特征指示聲速與沉積物物理性質(zhì)和沉積環(huán)境的關(guān)系, 而聲速垂向的變化規(guī)律是建立聲速結(jié)構(gòu)模型的重要依據(jù)。圖4、圖5和圖6分別顯示了實測聲速比、Buckingham模型計算聲速比和Biot-Stoll模型計算聲速比的變化規(guī)律, 可以分析出研究區(qū)沉積物聲速的橫向和垂向變化特征, 識別研究海域的高聲速區(qū)和低聲速區(qū); 還可以通過對比實測聲速與計算聲速的變化規(guī)律, 探討利用Buckingham模型和Biot-Stoll模型預測分析復雜環(huán)境沉積物聲速變化特征的適用性。

    3.2.1沉積物聲速平面的變化規(guī)律

    根據(jù)該海域沉積物實測聲速比變化特征圖(圖4)可以看出, 西南陸架區(qū)在15~45 cm深度范圍的聲速比小于1, 屬于低聲速區(qū); 在60~90 cm深度范圍的聲速比大于1, 屬于高聲速區(qū)。陸坡和海槽區(qū)的大部分聲速比低于1, 局部區(qū)域甚至低于0.95, 均屬于低聲速區(qū)。由此可知, 研究海域的高聲速沉積物主要分布于西南陸架和東北島礁區(qū), 而低聲速沉積物廣泛分布于陸坡和南沙海槽。

    圖4 沉積物的實測聲速比的變化特征Fig. 4 The distribution of the measured sound speed ratios

    物理參數(shù)測量結(jié)果顯示, 陸坡和海槽沉積物類型以粉砂質(zhì)黏土和黏土質(zhì)粉砂為主, 其中粉砂含量為40%~60%, 粒徑范圍為3~5 μm, 濕密度范圍為1.1~1.4 g/cm3; 西南陸架的沉積物類型則以砂質(zhì)粉砂為主, 砂含量為20%~50%, 粒徑范圍為30~70 μm, 濕密度范圍為1.46~1.83 g/cm3。根據(jù)大陸架、陸坡和海槽沉積物實測物理參數(shù), 主要的物理性質(zhì)呈現(xiàn)有規(guī)律的變化: 大陸架沉積物表現(xiàn)為粗顆粒、低孔隙度、高密度; 陸坡和海槽沉積物則呈現(xiàn)細顆粒、高孔隙度、低密度。這表明聲速比的分布特征與物理性質(zhì)、沉積環(huán)境等因素具有相關(guān)性。究其原因, 沉積環(huán)境決定了沉積物的孔隙度、粒徑、密度等物理性質(zhì), 而物理性質(zhì)又是影響沉積物聲速的重要因素, 因此聲速的分布特征宏觀上受到沉積環(huán)境的控制。

    圖5 Buckingham模型計算聲速比變化特征Fig. 5 The distribution of the calculated sound speed ratios based on Buckingham model

    圖6 Biot-Stoll模型計算聲速比分布Fig. 6 The distribution of the calculated sound speed ratios based on Biot-Stoll model

    3.2.2沉積物聲速垂向的變化規(guī)律

    通常認為, 由于上覆沉積物載荷壓力的增加和固結(jié)作用導致沉積物結(jié)構(gòu)發(fā)生變化, 其密度、粒徑等具有增大的變化特征, 而孔隙度、含水量呈現(xiàn)減少的趨勢, 從而使得沉積物聲速呈現(xiàn)隨深度增大的變化特征。然而這種認識與實際聲學測量的研究結(jié)果不一致: Hamilton(1980)建立了淺海大陸架和深海丘陵、深海平原的地聲模型, Lu等(2003)提出了南海北部大陸架低聲速和高聲速海底的聲速結(jié)構(gòu)模式, 這些研究均認為沉積物聲速及物理性質(zhì)呈現(xiàn)多種垂向的變化模式。然而這些地聲模型基于多個海域數(shù)米乃至數(shù)十米埋深的沉積物, 采用的測量手段也不相同, 具有較大的聲學測量尺度, 是否能闡述表層松散沉積物的聲速分布特征有待驗證。

    從研究區(qū)實測聲速比分布圖(圖4)可以看到,研究區(qū)沉積物聲速整體上隨深度呈現(xiàn)增大的趨勢,但在不同的沉積區(qū)域呈現(xiàn)出多種變化規(guī)律。西南部陸架區(qū)和東北島礁區(qū)沉積物的聲速在15~90 cm深度呈現(xiàn)逐漸增大的變化趨勢, 這是由于陸架區(qū)主要為陸源碎屑沉積物, 粗顆粒成分在該區(qū)沉降下來,而細顆粒成分被搬運至深水區(qū), 先沉降沉積物受到后期沉積物載荷壓力的作用, 結(jié)構(gòu)密實且孔隙度相對較小, 而表層沉積物的結(jié)構(gòu)較為松散且具有較大的孔隙度。東北區(qū)島礁周圍覆蓋粗顆粒的淺海珊瑚砂, 與鈣質(zhì)成分的深海泥質(zhì)沉積物相比, 島礁區(qū)沉積物的孔隙度和抗壓強度相對較低(李趕先和龍建軍, 2014), 表層沉積物固結(jié)時間較短其結(jié)構(gòu)較為松散, 而隨著埋深增加, 沉積物結(jié)構(gòu)逐漸被壓實。由孔隙度和聲速的相關(guān)關(guān)系可知, 這兩個區(qū)域沉積物的聲速呈現(xiàn)出增大的垂向分布規(guī)律。陸坡和海槽沉積物的聲速隨深度的變化特征較為復雜: 部分區(qū)域的沉積物聲速呈基本不變的趨勢, 這是由于水深變化相對較小且沉積物來源單一, 后期改造作用較弱,沉積物的物理性質(zhì)隨深度的變化幅度較小; 而部分區(qū)域沉積物的聲速隨深度呈現(xiàn)增大的趨勢, 主要原因與陸架區(qū)相似, 受上覆沉積物載荷壓力的增加的影響。上述變化趨勢表明, 研究海區(qū)海底的聲速結(jié)構(gòu)模式既包括垂向增大模式, 還包括垂向不變和垂向復雜變化等多種模式, 這與過往的研究結(jié)果相一致(Lu et al., 2003; 鄒大鵬等, 2012b; Wang et al., 2014)。究其原因, 沉積物的聲速與物理性質(zhì)具有相關(guān)性, 而物理性質(zhì)受沉積環(huán)境和沉積結(jié)構(gòu)的分層性所影響, 物理性質(zhì)尤其是孔隙度易受壓實作用的影響, 隨深度呈現(xiàn)出多種變化趨勢, 從而決定了聲速隨深度變化的復雜特征。因此, 較之深層沉積物受到壓實和固結(jié)作用的影響從而具有較為密實的結(jié)構(gòu),表層沉積物的結(jié)構(gòu)相對較為松散, 其聲學性質(zhì)介于流體(海水等)和固體(密實的沉積物、巖石等)之間,不能利用單一聲波傳播理論描述聲波在表層沉積物中傳播的特征, 也不能用單一的聲速垂向增大模式描述所有海域沉積物的聲速結(jié)構(gòu)。

    3.2.3Buckingham模型計算聲速與實測聲速的變化特征對比

    根據(jù)Buckingham模型計算聲速比變化特征圖(圖5), 西南陸架區(qū)沉積物聲速比大于1, 屬于高聲速區(qū)。陸坡和南沙海槽沉積物的聲速比普遍小于1,屬于低聲速區(qū)。與實測聲速比的分布特征相類似,計算聲速的平面分布特征與物理性質(zhì)、沉積環(huán)境具有相關(guān)性, 但較之實測聲速分布特征的復雜性和局部性, 計算聲速比的平面變化特征較為單一且缺乏突變的特點。從聲速比垂向的分布角度來看, 陸架區(qū)沉積物的聲速比垂向上呈現(xiàn)小幅度增大的趨勢,陸坡和海槽區(qū)沉積物的聲速比垂向上基本不變, 這與實測聲速比的垂向分布特征不一致。這表明, 利用Buckingham模型計算的聲速比可以匹配實測聲速比的橫向分布特征, 但不能體現(xiàn)實測聲速比分布的復雜性及變化幅度大等特點, 也不具備實測聲速比垂向分布規(guī)律的多樣性。

    3.2.4Biot-Stoll模型計算聲速與實測聲速的變化特征對比

    根據(jù)Biot-Stoll模型計算聲速比變化特征圖(圖6)分析, 與實測聲速比和Buckingham模型計算聲速比分布相似, 陸架區(qū)沉積物具有較高的聲速比, 陸坡和海槽沉積物的聲速比相對較低。大部分區(qū)域沉積物的聲速比大于1, 屬于高聲速區(qū), 僅陸坡和海槽局部區(qū)域的沉積物聲速比小于1, 屬于低聲速區(qū)。這表明Biot-Stoll模型計算聲速符合陸架區(qū)聲速的變化特征, 但不能匹配陸坡和海槽區(qū)聲速比的平面變化規(guī)律。計算聲速比垂向上呈現(xiàn)基本不變的特征,這與實測聲速比垂向的變化規(guī)律不一致。由上述分析可知, 與Buckingham模型相類似, Biot-Stoll模型計算聲速比能匹配實測聲速比的平面分布特征, 但不能體現(xiàn)實測聲速比垂向變化的復雜性及局部突變等特征。

    綜合圖4至圖6進行分析, 雖然兩種模型計算聲速的分布均與實測聲速的分布特征具有相似性,但實測聲速的分布具有復雜性和局部性。這表明理論模型預測不能完全體現(xiàn)復雜沉積環(huán)境對聲速分布的影響。究其原因, 由于砂質(zhì)和泥質(zhì)等類型沉積物的物理性質(zhì)之間存在著較大的差異, 不同類型沉積物的聲速分布特征也不一致。即使沉積物的類型相同, 但沉積結(jié)構(gòu)受外部環(huán)境的影響, 也是造成沉積物聲速分布的復雜特征的重要因素, 因此不可能存在單一模型能夠廣泛適用于描述各種類型沉積物的聲速分布規(guī)律。另一方面, 模型參數(shù)的選擇多基于沉積物的孔隙性或者顆粒大小等物理因素, 而這些參數(shù)往往難以直接測量, 或無法測量, 部分參數(shù)甚至沒有明確的物理意義, 大部分參數(shù)由均勻各向同性介質(zhì)得到而又被運用到非均勻沉積物中。模型參數(shù)的選擇一般采用曲線擬合和混合優(yōu)化的方法; Stoll參數(shù)將模型參數(shù)分為砂質(zhì)和泥質(zhì)沉積物兩類,均采用單數(shù)值形式; Schock參數(shù)把滲透率、孔隙因子和骨架體積模量等分別看作是孔隙度、粒徑和有效壓力的函數(shù), 建立出這些未知參數(shù)的計算公式,這些模型參數(shù)選擇和確定方法均得到了驗證和應用。但對于復雜沉積環(huán)境中的沉積物, 其模型參數(shù)的選擇仍需要基于大量的沉積物物理性質(zhì)和聲速測試, 按照沉積物類型和沉積環(huán)境等進行更為精細的分類, 并且考慮沉積物的孔隙形狀、堆壘方式和礦物成分, 總結(jié)實測物理參數(shù)與無法測量的物理參數(shù)之間的函數(shù)關(guān)系, 代替利用單一數(shù)值表示整個區(qū)域沉積物的簡化方法, 建立更精確的參數(shù)選擇方法。

    4 結(jié)論

    通過對南海南部表層沉積物進行聲學測量試驗和利用聲波傳播理論進行聲速計算, 對比了測量試驗與三種理論模型的結(jié)果, 分析了該海域的聲速空間變化特征, 得到以下認識:

    (1)實測聲速值與模型聲速計算值的對比分析表明, Wood方程計算值與實測值的相對偏差較大, Biot-Stoll模型適合于陸架高聲速區(qū)沉積物的聲速預測, Buckingham模型適合于更廣范圍海域沉積物聲速的計算。模型計算聲速不能體現(xiàn)實測聲速分布的復雜性和局部性。相比側(cè)重于沉積物孔隙性的Biot-Stoll模型, 兼具流體和彈性固體雙重性質(zhì)的Buckingham模型更適用于解釋該海域沉積物的聲波傳播特征。

    (2)通過分析聲速平面和垂向的變化特征, 表明沉積物實測聲速的空間變化特征與沉積環(huán)境和物理性質(zhì)具有相關(guān)性。由于壓實作用而造成的聲速垂向增大模式, 不能代表所有海域淺表層沉積物的聲速結(jié)構(gòu)。

    (3)在確定模型的參數(shù)時, 由于模型涉及參數(shù)較多, 一些參數(shù)還僅依靠文獻資料和值域范圍進行優(yōu)化擬合確定, 因此需要針對不同類型的海洋沉積物進行大量的物理性質(zhì)和聲速測試, 建立實測參數(shù)和未知參數(shù)的函數(shù)關(guān)系, 從而更精確的選擇模型參數(shù)。Acknowledgements:

    This study was supported by the Public Science and Technology Research Funds Projects of Ocean (Nos. 201405032 and 200905025), and the National Natural Science Foundation of China (No. 41330965).

    參考文獻:

    陳靜, 閻貧, 王彥林, 金丹, 林秋金. 2012. 基于Biot-Stoll模型聲速反演中的參數(shù)選擇——以南海南部沉積物為例[J]. 熱帶海洋學報, (1): 50-54.

    付淑清, 朱照宇, 邱燕, 歐陽婷萍. 2008. 南海南部陸坡晚第四紀沉積環(huán)境變化[J]. 熱帶地理, 28(2): 99-102.

    李傲仙, 李延河, 樂國良. 2005. 深海沉積物中碲異常的成因[J].地球?qū)W報, 26(S1): 186 - 189.

    李趕先, 龍建軍. 2014. 南海南部海域島礁區(qū)海底珊瑚砂聲速影響因素的初步研究[J]. 海洋學報, 36(5): 152-160.

    楊濤, 薛紫晨, 楊競紅, 蔣少涌. 2003. 南海北部地區(qū)海洋沉積物中孔隙水的氫、氧同位素組成特征[J]. 地球?qū)W報, 24(6): 511-514.

    張福生, 盧博, 李趕先, 黃韶健. 2001. 南沙群島附近海域海底聲學性質(zhì)的分析與研究[J]. 湛江海洋大學學報, 21(S1): 116-122.

    周豐年, 趙建虎, 周才揚. 2001. 多波束測深系統(tǒng)最優(yōu)聲速公式的確定[J]. 臺灣海峽, 20(4): 411-419.

    朱祖揚, 王東, 周建平, 王秀明. 2012. 基于非飽和Biot-Stoll模型的海底沉積物介質(zhì)聲頻散特性研究[J]. 地球物理學報, 55(1): 180-188.

    鄒大鵬, 盧博, 閻貧, 吳百海. 2012a. 南海北部海底沉積物在溫度變化下的三種聲速類型[J]. 地球物理學報, 55(3): 1017-1024.

    鄒大鵬, 閻貧, 盧博. 2012b. 基于海底表層沉積物聲速特征的南海地聲模型[J]. 海洋學報: 中文版, 34(3): 80-86.

    References:

    BERRMAN J G. 1999. Origin of Gassmann’s equations[J]. Geophysics, 64: 1627-1629.

    BIOT M A. 1956a. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low frequency range[J]. J. Acoust. Soc. Am., 28: 168-178.

    BIOT M A. 1956b. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range[J]. J. Acoust. Soc. Am., 28: 179-191.

    BRIGGS K B, RICHARDSON M D. 1997. Small-scale fluctuations in acoustic and physical properties in surficial carbonate sediments and their relationship to bioturbation[J]. Geo-Marine Letters, 17(4): 306-315.

    BUCKINGHAM M J. 1997. Theory of acoustic attenuation, dispersion, and pulse propagation in unconsolidated granular materials including marine sediments[J]. J. Acoust. Soc. Am., 102(5): 2579-2596.

    BUCKINGHAM M J. 1998. Theory of compressional and shear waves in fluidlike marine sediments[J]. J. Acoust. Soc. Am., 103(1): 288-299.

    BUCKINGHAM M J. 2005. Compressional and shear wave properties of marine sediments: comparisons between theory and data[J]. J. Acoust. Soc. Am., 117(1): 137-152.

    CHEN Jing, YAN Pin, WANG Yan-lin, LIN Qiu-jin. 2012. Choice of Parameters for Biot-Stoll Model- based Inversion of Sound Velocity of Seafloor Sediments in the Southern South China Sea[J]. Journal of Tropical Oceanography, (1): 50-54(in Chinese with English abstract).

    CHOTIROS N P, ISAKSON M J. 2004. A broadband model of sandy ocean sediments: Biot-Stoll with contact squirtflow and shear drag[J]. J. Acoust. Soc. Am., 116: 2011-2022.

    FU Shu-qing, ZHU Zhao-yu, QIU Yan, OUYANG Ting-ping. 2008. The Changes of the Depositional Environment of the Southern South China Sea Continental Slope during Late Quaternary[J]. Tropical Geography, 28(2): 99-102(in Chinese with English abstract).

    HAMILTON E L. 1971. Prediction of in-situ acoustic and elastic properties of marine sediments[J]. Geophysics, 36(2): 266-284.

    HAMILTON E L. 1980. Geoacoustic modeling of the sea floor[J]. J. Acoust. Soc. Am., 68: 1313-1340.

    HAMILTON E L, BACHMAN R T. 1982. Sound velocity and related properties of marine sediments[J]. J. Acoust. Soc. Am., 72: 1891-1904.

    HAMILTON E L. 1985. Sound velocity as a function of depth in marine sediments[J]. J. Acoust. Soc. Am., 78: 1348-1355.

    HOU Zheng-yu, GUO Chang-sheng, WANG Jing-qing, LI Hui-yin, LI Tie-gang. 2014. Tests of new in-situ seabed acoustic measurement system in Qingdao[J]. Chinese Journal of Oceanology and Limnology, 32(5): 1172-1178.

    JACKSON D R, RICHARDSON M D. 2007. High-frequency seafloor acoustics[M]. New York, NY, USA: Springer Science + Business Media.

    LI Ao-xian, LI Yan-he, YUE Guo-liang. 2005. The Cause for Tellurium Enrichment in Deep-Sea Sediments[J]. Acta Geoscientica Sinica, 26(S1): 186-189(in Chinese with English abstract).

    LI Gan-xian, LONG Jian-jun. 2014. A preliminary study of the sound velocity influence factors of submarine coral sand of islands sea area in southern South China Sea[J]. Acta Oceanologica Sinica, 36(5): 152-160(in Chinese with English abstract).

    LIU Bao-hua, HAN Tong-cheng, KAN Guang-ming, LI Guan-bao. 2013. Correlations between the in situ acoustic properties and geotechnical parameters of sediments in the Yellow Sea[J]. China. J. Asian Earth Sci, 77: 83-90.

    LIU Jian-guo, XIANG Rong, CHEN Zhong, CHEN Mu-hong, YAN Wen, ZHANG Lan-lan, CHEN H. 2013. Sources, transport and deposition of surface sediments from the South China Sea[J]. Deep Sea Research Part I: Oceanographic Research Papers, 71: 92-102.

    LU Bo, HUANG Shao-jian, LI Gan-xian, ZHANG Fu-sheng. 2003. Vertical variations of core sound velocity: Evidence of paleooceanographic history since the Pleistocene epoch[J]. Marine Georesources and Geotechnology, 21(2): 63-71.

    LU Bo, LIU Qiang, LI Gan-xian. 2010. Grain and pore factors in acoustic response to seafloor sediments[J]. Marine Georesources and Geotechnology, 28(2): 115-129.

    SCHOCK S G. 2004. A method for estimating the physical and acoustic properties of the sea bed using chirp sonar data[J]. IEEE Journal of Oceanic Engineering, 29(4): 1200-1217.

    STOLL R D. 1989. Sediment Acoustics[M]. New York: Lecture Notes in Earth Sciences, Springer-Verlag.

    STOLL R D, BAUTISTA E O. 1998. Using the Biot theory to establish a baseline geoacoustic model for seafloor sediments[J]. Continental shelf research, 18(14): 1839-1857.

    WANG Jing-qiang, GUO Chang-sheng, HOU Zheng-yu, FU Yong-tao, YAN Jun. 2014. Distributions and vertical variation patterns of sound speed of surface sediments in South China[J]. Sea. J. Asian Earth Sci., 89: 46-53.

    WILLIAMS K L, JACKSON D R, THORSOS E I, TANG D J, SCHOCK S S. 2002. Comparison of sound speed and attenuation measured in a sandy sediment to predictions based on the Biot theory of porous media[J]. IEEE Journal of Oceanic Engineering, 27(3): 413-428.

    YANG Tao, XUE Zi-chen, YANG Jing-hong, JIANG Shao-yong. 2003. Oxygen and Hydrogen Isotopic Compositions of Pore Water from Marine Sediments in the Northern South China Sea[J]. Acta Geoscientica Sinica, 24(6): 511-514(in Chinese with English abstract).

    ZHANG Fu-sheng, LU Bo, LI Gan-xian, HUANG Shao-jian. 2001. A Study on the Seafloor Sediments Acoustic Characteristics of Nansha Islands Sea Area[J]. Journal of Zhanjiang Ocean University, 21(Sl): 116-122(in Chinese with English abstract).

    ZHOU Feng-nian, ZHAO Jian-hu, ZHOU Cai-yang. 2001. Determination of Classic Experiential Sound Speed Formulae in Multibeam Echo Sounding System[J]. Journal of Oceanography in Taiwan Strait, 20(4): 411-419(in Chinese with English abstract).

    ZHOU Zu-yang, WANG Dong, ZHOU Jian-ping, WANG Xiu-ming. 2012. Acoustic Wave Dispersion and Attenuation in Marine Sediment Based on Partially Gas-saturated Biot-Stoll Model[J]. Chinese J. Geophys., 55(1): 180-188(in Chinese with English abstract).

    ZOU Da-peng, LU Bo, YAN Pin, WU Bo-hai. 2012. Three Kinds of Acoustic Speeds of Seafloor Sediments in the Northern South China Sea with Temperature Variation[J]. Chinese J. Geophys, 55(3): 1017-1024(in Chinese with English abstract).

    ZOU Da-peng, YAN Pin, LU Bo. 2012. A Geoacoustic Model Based on Sound Speed Characteristic of Seafloor Surface Sediments of the South China Sea[J]. Acta Oceanologica Sinaica, 34(3): 80-86(in Chinese with English abstract).

    Sound Speed Distribution of Seafloor Sediments in Nansha Islands Sea Based on Buckingham Model and Biot-Stoll Model

    WANG Jing-qiang1, 4), GUO Chang-sheng2, 4)*, LIU Bao-hua3, 4), CHEN Wen-jing2, 5), HOU Zheng-yu2, 5), HAN Guo-zhong1, 4)
    1) Key Laboratory of Marine Sedimentology and Environment Geology, First Institute of Oceanography,
    State Oceanic Administration, Qingdao, Shandong 266061; 2) Key Laboratory of Marine Geology and Environment, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, Shandong 266071; 3) National Deep Sea Center, State Oceanic Administration, Qingdao, Shandong 266061; 4) Laboratory for Marine Geology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266061; 5) University of Chinese Academy of Sciences, Beijing 100049

    Abstract:In order to analyze the application of the speed theoretical models to predicting the sound speed of sediments in the southern South China Sea, the authors calculated the sound speed by using Wood equation, Buckingham model, and Biot-Stoll model respectively. The comparison between the calculated and the measured sound speeds shows that the average deviation between the calculated sound speeds with Buckingham model and the measured sound speeds is smaller than the calculated values with Wood equation and Biot-Stoll model. The comparison of the sound speed ratio distributions indicate that the calculated sound speed ratios with Buckingham model and Biot-Stoll model have the similar horizontal distributions with the measured sound speed ratios. Thesound speeds have different horizontal distributions in different sedimentary environments and are closely related to the physical properties and the deposition of sediments. The calculated sound speed ratios are invariant with the depth, while the measured values have different vertical distributions in several different geomorphic units. According to the comparison results, the Biot-Stoll model is only suitable for the coarse sediments in the continental shelf, and the Buckingham model is more applicable to predicating the sound speed of sediments in the study area. The research is significant in obtaining the sound speed distribution of the sediments in the southern South China Sea, and can provide reference for the multi-beam echo sounding and the sub-bottom survey.

    Key words:seafloor sediment; Buckingham model; Biot-Stoll model; sound speed distribution; sound speed ratio

    中圖分類號:O422.1; O424

    文獻標志碼:A

    doi:10.3975/cagsb.2016.03.13

    收稿日期:2016-02-19; 改回日期: 2016-03-25。責任編輯: 魏樂軍。

    第一作者簡介:王景強, 男, 1989年生。博士后。主要從事海洋沉積聲學研究。E-mail: wang_jingqiang@aliyun.com。

    *通訊作者:郭常升, 男, 1964年生。研究員。主要從事海洋地質(zhì)聲學和海洋地球物理研究。E-mail: guochine@qdio.ac.cn。

    猜你喜歡
    聲速聲學沉積物
    晚更新世以來南黃海陸架沉積物源分析
    海洋通報(2022年2期)2022-06-30 06:07:04
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對磷的吸收及釋放研究進展
    愛的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學)Balthus 70
    Acoustical Treatment Primer:Diffusion談談聲學處理中的“擴散”
    Acoustical Treatment Primer:Absorption談談聲學處理中的“吸聲”(二)
    Acoustical Treatment Primer:Absorption 談談聲學處理中的“吸聲”
    聲速是如何測定的
    討論用ICP-AES測定土壤和沉積物時鈦對鈷的干擾
    跨聲速風洞全模顫振試驗技術(shù)
    航空學報(2015年4期)2015-05-07 06:43:32
    人人妻人人添人人爽欧美一区卜| 在线观看免费午夜福利视频| 国产亚洲欧美在线一区二区| www.999成人在线观看| 国产成人系列免费观看| 国产野战对白在线观看| 久久婷婷成人综合色麻豆| 男人舔女人的私密视频| 亚洲国产看品久久| 亚洲,欧美精品.| 久久中文看片网| 日本精品一区二区三区蜜桃| 日韩熟女老妇一区二区性免费视频| 国产成人精品在线电影| 51午夜福利影视在线观看| 国产成人影院久久av| 交换朋友夫妻互换小说| 久久免费观看电影| 欧美精品一区二区免费开放| 日本a在线网址| av超薄肉色丝袜交足视频| 一区二区三区精品91| tocl精华| 亚洲精华国产精华精| 亚洲性夜色夜夜综合| 午夜福利欧美成人| 精品一区二区三区视频在线观看免费 | 欧美日韩福利视频一区二区| 97在线人人人人妻| 久久亚洲真实| 午夜91福利影院| 国产欧美日韩一区二区三| 久久ye,这里只有精品| 日本av免费视频播放| 大香蕉久久网| 18禁观看日本| 在线观看免费日韩欧美大片| 一本综合久久免费| 久久精品国产亚洲av高清一级| 国产1区2区3区精品| 国产av又大| 欧美黑人欧美精品刺激| 国产高清激情床上av| 在线永久观看黄色视频| 国产伦理片在线播放av一区| 免费看a级黄色片| 黄色a级毛片大全视频| 国产视频一区二区在线看| 天堂8中文在线网| 成人18禁在线播放| www.熟女人妻精品国产| av天堂久久9| 91老司机精品| 国产成人av激情在线播放| 91老司机精品| 一区二区三区国产精品乱码| 国产成人精品久久二区二区免费| 亚洲专区字幕在线| 美女国产高潮福利片在线看| 黑人巨大精品欧美一区二区mp4| 妹子高潮喷水视频| 黑人巨大精品欧美一区二区蜜桃| 精品国产一区二区三区四区第35| 啪啪无遮挡十八禁网站| 国产精品 欧美亚洲| 亚洲一区中文字幕在线| 搡老岳熟女国产| 成在线人永久免费视频| 日韩大片免费观看网站| 老熟妇乱子伦视频在线观看| 亚洲精品一二三| 18在线观看网站| 一级毛片电影观看| 超碰97精品在线观看| 国产精品免费一区二区三区在线 | cao死你这个sao货| 正在播放国产对白刺激| 午夜两性在线视频| 国产欧美亚洲国产| 桃花免费在线播放| 色老头精品视频在线观看| 国产在线观看jvid| av网站免费在线观看视频| 大码成人一级视频| 侵犯人妻中文字幕一二三四区| 日韩 欧美 亚洲 中文字幕| 亚洲av日韩在线播放| 精品欧美一区二区三区在线| 国产精品一区二区精品视频观看| 欧美黑人欧美精品刺激| 视频区图区小说| 免费一级毛片在线播放高清视频 | 人妻 亚洲 视频| 1024视频免费在线观看| 精品久久蜜臀av无| 国产精品九九99| 亚洲自偷自拍图片 自拍| 欧美乱妇无乱码| 在线av久久热| 午夜激情久久久久久久| 黄网站色视频无遮挡免费观看| 久久中文字幕一级| 最近最新中文字幕大全电影3 | 岛国毛片在线播放| 午夜福利免费观看在线| 国产高清videossex| 高清在线国产一区| 天天躁日日躁夜夜躁夜夜| 日韩精品免费视频一区二区三区| 久久性视频一级片| 国产不卡一卡二| 飞空精品影院首页| 精品国内亚洲2022精品成人 | 激情在线观看视频在线高清 | 国产亚洲午夜精品一区二区久久| av天堂久久9| 脱女人内裤的视频| 99精品在免费线老司机午夜| 欧美成狂野欧美在线观看| 欧美日韩亚洲综合一区二区三区_| 在线观看一区二区三区激情| 日韩成人在线观看一区二区三区| 日本五十路高清| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品国产av在线观看| 久久精品aⅴ一区二区三区四区| 亚洲中文字幕日韩| 一本大道久久a久久精品| 最近最新中文字幕大全电影3 | 久久99热这里只频精品6学生| 日本撒尿小便嘘嘘汇集6| 欧美黑人精品巨大| 久久久久久久大尺度免费视频| videos熟女内射| 亚洲av成人不卡在线观看播放网| 亚洲av成人不卡在线观看播放网| 国产精品.久久久| 亚洲国产av影院在线观看| 91精品国产国语对白视频| 亚洲专区中文字幕在线| 日韩欧美一区视频在线观看| 老汉色av国产亚洲站长工具| 他把我摸到了高潮在线观看 | 美女高潮到喷水免费观看| 1024香蕉在线观看| videos熟女内射| av线在线观看网站| 日本a在线网址| 午夜老司机福利片| 法律面前人人平等表现在哪些方面| 亚洲成a人片在线一区二区| 国产一区二区三区综合在线观看| 三上悠亚av全集在线观看| 国产欧美日韩一区二区三区在线| 大型黄色视频在线免费观看| 999精品在线视频| 亚洲成人国产一区在线观看| 欧美乱妇无乱码| 黄色毛片三级朝国网站| 乱人伦中国视频| 亚洲精品久久成人aⅴ小说| 另类亚洲欧美激情| 久久精品亚洲熟妇少妇任你| tocl精华| 国产xxxxx性猛交| 午夜激情av网站| 国产伦理片在线播放av一区| 大香蕉久久网| 在线天堂中文资源库| 日韩视频在线欧美| 久久国产亚洲av麻豆专区| 亚洲性夜色夜夜综合| a级毛片黄视频| 亚洲avbb在线观看| 亚洲免费av在线视频| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人一区二区三| 色视频在线一区二区三区| 亚洲国产欧美在线一区| 狂野欧美激情性xxxx| 国产精品欧美亚洲77777| 多毛熟女@视频| 在线 av 中文字幕| 日韩三级视频一区二区三区| 久久九九热精品免费| 日韩大片免费观看网站| 最新的欧美精品一区二区| 90打野战视频偷拍视频| 黄色视频不卡| 大陆偷拍与自拍| 丰满迷人的少妇在线观看| 后天国语完整版免费观看| 精品福利永久在线观看| 国产精品偷伦视频观看了| 国产区一区二久久| 国产精品一区二区在线不卡| 国产精品98久久久久久宅男小说| 男女午夜视频在线观看| 99re6热这里在线精品视频| 亚洲国产欧美一区二区综合| 老司机午夜福利在线观看视频 | 又紧又爽又黄一区二区| 欧美精品高潮呻吟av久久| 欧美日韩中文字幕国产精品一区二区三区 | 黄网站色视频无遮挡免费观看| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 久久精品国产综合久久久| 亚洲国产看品久久| 成年人黄色毛片网站| 亚洲精品久久午夜乱码| 在线观看人妻少妇| 久久天躁狠狠躁夜夜2o2o| 我的亚洲天堂| 精品一区二区三卡| 国产国语露脸激情在线看| 久久精品亚洲精品国产色婷小说| a级片在线免费高清观看视频| av福利片在线| 曰老女人黄片| 久久久欧美国产精品| 99精品久久久久人妻精品| 成人永久免费在线观看视频 | 精品少妇久久久久久888优播| 日韩制服丝袜自拍偷拍| a级毛片在线看网站| 国精品久久久久久国模美| 亚洲第一av免费看| 99香蕉大伊视频| 99国产精品一区二区三区| 91麻豆精品激情在线观看国产 | 成年人黄色毛片网站| 亚洲专区字幕在线| 亚洲精品美女久久av网站| 日本撒尿小便嘘嘘汇集6| 色老头精品视频在线观看| 欧美激情 高清一区二区三区| 精品一区二区三区四区五区乱码| 免费在线观看日本一区| 飞空精品影院首页| 热re99久久精品国产66热6| 50天的宝宝边吃奶边哭怎么回事| 美女国产高潮福利片在线看| 国产精品久久久av美女十八| 久久久欧美国产精品| 久久中文字幕一级| 国产国语露脸激情在线看| 精品高清国产在线一区| 精品亚洲乱码少妇综合久久| 亚洲熟妇熟女久久| 精品久久久久久电影网| av天堂久久9| 美女视频免费永久观看网站| 国产成人精品在线电影| 亚洲人成伊人成综合网2020| 免费在线观看视频国产中文字幕亚洲| 国产高清视频在线播放一区| 久久久欧美国产精品| 免费在线观看黄色视频的| 熟女少妇亚洲综合色aaa.| 久久国产精品影院| 精品一品国产午夜福利视频| 美女高潮到喷水免费观看| 亚洲人成伊人成综合网2020| 欧美国产精品va在线观看不卡| 一区二区日韩欧美中文字幕| cao死你这个sao货| 操美女的视频在线观看| 老司机亚洲免费影院| 国产主播在线观看一区二区| 黄片大片在线免费观看| 国产一区二区 视频在线| 欧美 亚洲 国产 日韩一| 国产不卡一卡二| 中文欧美无线码| 看免费av毛片| 精品国产超薄肉色丝袜足j| 考比视频在线观看| 欧美日韩视频精品一区| 精品国产国语对白av| 丁香欧美五月| cao死你这个sao货| 亚洲国产av新网站| 欧美日韩亚洲高清精品| h视频一区二区三区| 两性夫妻黄色片| 国产精品一区二区免费欧美| 午夜精品国产一区二区电影| 18禁美女被吸乳视频| 日本av手机在线免费观看| 亚洲欧美激情在线| 久久久欧美国产精品| 中亚洲国语对白在线视频| 亚洲欧洲日产国产| 少妇 在线观看| 国产精品国产av在线观看| 久久精品亚洲av国产电影网| 免费在线观看视频国产中文字幕亚洲| 国产精品av久久久久免费| 免费在线观看日本一区| 久久久精品免费免费高清| 成人手机av| 精品欧美一区二区三区在线| 99久久99久久久精品蜜桃| 老司机靠b影院| 99精品欧美一区二区三区四区| 老司机福利观看| 极品少妇高潮喷水抽搐| 日本vs欧美在线观看视频| 在线亚洲精品国产二区图片欧美| 如日韩欧美国产精品一区二区三区| 久热爱精品视频在线9| 在线看a的网站| 国产在视频线精品| 女性生殖器流出的白浆| 国产1区2区3区精品| 一本一本久久a久久精品综合妖精| 国产精品亚洲av一区麻豆| 亚洲精品成人av观看孕妇| 国产精品久久久久久精品古装| 最新美女视频免费是黄的| 色在线成人网| 波多野结衣一区麻豆| 黄色毛片三级朝国网站| 国产精品香港三级国产av潘金莲| 久久午夜亚洲精品久久| 超碰成人久久| 啦啦啦 在线观看视频| 一本—道久久a久久精品蜜桃钙片| 亚洲三区欧美一区| 日韩欧美免费精品| 正在播放国产对白刺激| 久久久久久人人人人人| 狠狠婷婷综合久久久久久88av| 丰满人妻熟妇乱又伦精品不卡| 好男人电影高清在线观看| 久久久久久久精品吃奶| 色视频在线一区二区三区| 男人舔女人的私密视频| 免费看十八禁软件| 黑人操中国人逼视频| 久久毛片免费看一区二区三区| 久久香蕉激情| 亚洲va日本ⅴa欧美va伊人久久| 国产精品免费一区二区三区在线 | 国产主播在线观看一区二区| 国产精品久久电影中文字幕 | 99re6热这里在线精品视频| 久久ye,这里只有精品| 成年动漫av网址| 精品久久蜜臀av无| 男女无遮挡免费网站观看| 日韩中文字幕视频在线看片| 国产不卡一卡二| 一级片'在线观看视频| tube8黄色片| 国产成人欧美在线观看 | 亚洲欧美一区二区三区黑人| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人巨大精品欧美一区二区蜜桃| 极品人妻少妇av视频| 天天躁夜夜躁狠狠躁躁| 99热国产这里只有精品6| 欧美 亚洲 国产 日韩一| 91字幕亚洲| av一本久久久久| 国产深夜福利视频在线观看| 男女下面插进去视频免费观看| 高清在线国产一区| 亚洲全国av大片| 久久久久国内视频| 国产亚洲精品一区二区www | 色在线成人网| 亚洲一区二区三区欧美精品| 中文字幕av电影在线播放| 极品少妇高潮喷水抽搐| 免费人妻精品一区二区三区视频| 男男h啪啪无遮挡| 国产精品香港三级国产av潘金莲| 成人手机av| 欧美精品啪啪一区二区三区| www.999成人在线观看| 精品第一国产精品| 中文字幕人妻丝袜一区二区| 啦啦啦视频在线资源免费观看| 国产欧美日韩精品亚洲av| 我的亚洲天堂| 久久久国产精品麻豆| 王馨瑶露胸无遮挡在线观看| 国产有黄有色有爽视频| a级毛片黄视频| 久久亚洲真实| 国产亚洲一区二区精品| 麻豆国产av国片精品| 老司机靠b影院| 国产欧美日韩综合在线一区二区| 国产精品国产高清国产av | 国产精品久久电影中文字幕 | 欧美大码av| 建设人人有责人人尽责人人享有的| 精品国产乱码久久久久久男人| 久久热在线av| 亚洲少妇的诱惑av| 制服诱惑二区| 国产精品亚洲一级av第二区| 亚洲九九香蕉| 久久国产亚洲av麻豆专区| 久久久欧美国产精品| 亚洲成a人片在线一区二区| 亚洲精品一卡2卡三卡4卡5卡| 99re6热这里在线精品视频| tocl精华| 9热在线视频观看99| 国产男女超爽视频在线观看| 欧美另类亚洲清纯唯美| 国产一卡二卡三卡精品| 天堂8中文在线网| 亚洲熟妇熟女久久| 香蕉久久夜色| 久久久久久久久久久久大奶| 亚洲精品在线美女| 国产精品秋霞免费鲁丝片| 老熟女久久久| 大码成人一级视频| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 欧美激情极品国产一区二区三区| 99久久99久久久精品蜜桃| 黄频高清免费视频| 侵犯人妻中文字幕一二三四区| a在线观看视频网站| 久久久久久免费高清国产稀缺| 国产欧美日韩一区二区精品| 十八禁高潮呻吟视频| 免费观看av网站的网址| 老熟妇乱子伦视频在线观看| 高清毛片免费观看视频网站 | 两个人看的免费小视频| 女性生殖器流出的白浆| 超碰成人久久| 亚洲av国产av综合av卡| 纯流量卡能插随身wifi吗| 午夜老司机福利片| 少妇被粗大的猛进出69影院| av一本久久久久| 亚洲成av片中文字幕在线观看| 国产免费视频播放在线视频| 18禁美女被吸乳视频| 国产精品免费视频内射| 91国产中文字幕| 好男人电影高清在线观看| 一本—道久久a久久精品蜜桃钙片| 男女无遮挡免费网站观看| 成人18禁高潮啪啪吃奶动态图| 啦啦啦 在线观看视频| 亚洲国产欧美网| 男女之事视频高清在线观看| 亚洲专区国产一区二区| 如日韩欧美国产精品一区二区三区| 欧美+亚洲+日韩+国产| 亚洲人成伊人成综合网2020| 在线观看免费视频网站a站| 成人影院久久| 国产在线免费精品| 五月天丁香电影| 黄片大片在线免费观看| av线在线观看网站| 国产亚洲午夜精品一区二区久久| 高清视频免费观看一区二区| 在线观看舔阴道视频| 亚洲黑人精品在线| 不卡一级毛片| 国产精品美女特级片免费视频播放器 | 亚洲国产欧美一区二区综合| 亚洲精品国产区一区二| 免费不卡黄色视频| 久久久精品免费免费高清| 亚洲精品久久午夜乱码| 少妇裸体淫交视频免费看高清 | 精品久久蜜臀av无| 欧美精品高潮呻吟av久久| 在线av久久热| 在线亚洲精品国产二区图片欧美| 成年人黄色毛片网站| 手机成人av网站| 人成视频在线观看免费观看| 亚洲一区中文字幕在线| 狠狠婷婷综合久久久久久88av| 91九色精品人成在线观看| 精品少妇黑人巨大在线播放| 色尼玛亚洲综合影院| 国产又爽黄色视频| 日韩大片免费观看网站| 国产欧美日韩一区二区三区在线| tocl精华| 欧美日韩黄片免| 免费人妻精品一区二区三区视频| 国产视频一区二区在线看| 丝袜人妻中文字幕| 国产精品1区2区在线观看. | 国产在线视频一区二区| 国产精品成人在线| 高潮久久久久久久久久久不卡| 中文字幕av电影在线播放| videosex国产| 精品久久久精品久久久| 大型av网站在线播放| 欧美日韩精品网址| 亚洲精品自拍成人| 超碰97精品在线观看| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 精品国产一区二区久久| 一个人免费看片子| 美女扒开内裤让男人捅视频| 日韩欧美一区视频在线观看| 亚洲精品国产一区二区精华液| 国产成人系列免费观看| 亚洲成a人片在线一区二区| 女人高潮潮喷娇喘18禁视频| 一区二区av电影网| 黑人巨大精品欧美一区二区蜜桃| 国产成人一区二区三区免费视频网站| 丝瓜视频免费看黄片| 精品一区二区三区四区五区乱码| 在线观看一区二区三区激情| 亚洲情色 制服丝袜| 国产有黄有色有爽视频| 在线观看www视频免费| av片东京热男人的天堂| 国产一卡二卡三卡精品| av电影中文网址| 这个男人来自地球电影免费观看| 妹子高潮喷水视频| 18禁黄网站禁片午夜丰满| 亚洲成人免费av在线播放| 久久热在线av| 免费观看av网站的网址| 日本黄色视频三级网站网址 | 午夜激情久久久久久久| 国产av一区二区精品久久| 久久国产精品人妻蜜桃| 国产精品自产拍在线观看55亚洲 | 999久久久国产精品视频| 桃花免费在线播放| 久久中文看片网| 香蕉久久夜色| 如日韩欧美国产精品一区二区三区| 亚洲欧美色中文字幕在线| 精品一区二区三区av网在线观看 | 王馨瑶露胸无遮挡在线观看| a在线观看视频网站| 久久久国产欧美日韩av| 亚洲欧美日韩高清在线视频 | 超碰97精品在线观看| 可以免费在线观看a视频的电影网站| 中文字幕人妻丝袜一区二区| 五月天丁香电影| 国产精品美女特级片免费视频播放器 | 亚洲色图 男人天堂 中文字幕| 丁香欧美五月| 一级毛片精品| 久久亚洲真实| 97在线人人人人妻| 久久影院123| 亚洲九九香蕉| 久久婷婷成人综合色麻豆| 老司机午夜福利在线观看视频 | 搡老乐熟女国产| 少妇裸体淫交视频免费看高清 | 在线观看免费高清a一片| 久久精品91无色码中文字幕| 国产日韩欧美亚洲二区| 黄色a级毛片大全视频| 99国产精品一区二区蜜桃av | 黄色视频不卡| 在线观看66精品国产| 两人在一起打扑克的视频| 99re6热这里在线精品视频| 一边摸一边抽搐一进一出视频| 女警被强在线播放| 18禁国产床啪视频网站| 久久精品亚洲熟妇少妇任你| 在线观看免费午夜福利视频| 他把我摸到了高潮在线观看 | 一区二区三区精品91| 久久久国产欧美日韩av| 日日夜夜操网爽| 中国美女看黄片| 少妇精品久久久久久久| 丁香六月欧美| 欧美日韩成人在线一区二区| 巨乳人妻的诱惑在线观看| 午夜福利视频精品| 大香蕉久久成人网| 国产精品久久久av美女十八| 免费一级毛片在线播放高清视频 | 韩国精品一区二区三区| 两个人看的免费小视频| 人成视频在线观看免费观看| 精品人妻1区二区| 久久影院123| 在线看a的网站| 国产精品影院久久| 亚洲欧美日韩高清在线视频 | 999久久久国产精品视频| 中文字幕最新亚洲高清| 国产高清激情床上av| 王馨瑶露胸无遮挡在线观看| 久久午夜亚洲精品久久| 精品人妻在线不人妻| 国产成+人综合+亚洲专区| 99热国产这里只有精品6| 久久中文看片网| 老熟妇乱子伦视频在线观看|