• <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
    国产欧美日韩一区二区三| 日本 欧美在线| 一区福利在线观看| 少妇高潮的动态图| 熟妇人妻久久中文字幕3abv| 手机成人av网站| 全区人妻精品视频| 亚洲av成人不卡在线观看播放网| 此物有八面人人有两片| 在线视频色国产色| 亚洲美女视频黄频| a级毛片a级免费在线| 十八禁网站免费在线| 男女床上黄色一级片免费看| 国产国拍精品亚洲av在线观看 | 国产精品日韩av在线免费观看| 一个人免费在线观看的高清视频| 国产精品三级大全| 国内毛片毛片毛片毛片毛片| 99riav亚洲国产免费| 国产精品美女特级片免费视频播放器| 首页视频小说图片口味搜索| 精品久久久久久久久久久久久| 可以在线观看毛片的网站| 香蕉久久夜色| 日韩国内少妇激情av| 熟女电影av网| 欧美+亚洲+日韩+国产| 3wmmmm亚洲av在线观看| 精华霜和精华液先用哪个| 国产主播在线观看一区二区| 高清日韩中文字幕在线| 在线免费观看的www视频| 久久久久久久久久黄片| 日韩欧美三级三区| 淫秽高清视频在线观看| 超碰av人人做人人爽久久 | 俺也久久电影网| 中文字幕精品亚洲无线码一区| 国产美女午夜福利| 99在线视频只有这里精品首页| 男女视频在线观看网站免费| 欧美日韩乱码在线| 超碰av人人做人人爽久久 | 午夜老司机福利剧场| 免费看十八禁软件| 中国美女看黄片| 午夜福利在线观看吧| 亚洲欧美日韩卡通动漫| 高清毛片免费观看视频网站| 99热6这里只有精品| 美女 人体艺术 gogo| 特大巨黑吊av在线直播| 成人18禁在线播放| 精品国内亚洲2022精品成人| 国产一区二区三区视频了| 午夜免费观看网址| 欧美成人一区二区免费高清观看| 高清在线国产一区| 小说图片视频综合网站| 久久精品国产亚洲av香蕉五月| 亚洲七黄色美女视频| 久久精品国产99精品国产亚洲性色| 国产熟女xx| 人人妻,人人澡人人爽秒播| 国产精品三级大全| 成年女人毛片免费观看观看9| 日日摸夜夜添夜夜添小说| 999久久久精品免费观看国产| 午夜免费观看网址| 中国美女看黄片| 亚洲av电影不卡..在线观看| 成人永久免费在线观看视频| 中文字幕av成人在线电影| 日本 av在线| 偷拍熟女少妇极品色| 成人性生交大片免费视频hd| 欧美最黄视频在线播放免费| 亚洲精品在线观看二区| 老汉色av国产亚洲站长工具| 久久性视频一级片| 99久久精品一区二区三区| 全区人妻精品视频| 美女高潮喷水抽搐中文字幕| 精品一区二区三区视频在线 | 最新中文字幕久久久久| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美激情综合另类| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 丰满人妻一区二区三区视频av | 黄色丝袜av网址大全| 亚洲欧美日韩东京热| 女人高潮潮喷娇喘18禁视频| 亚洲无线在线观看| 两人在一起打扑克的视频| 成人精品一区二区免费| 亚洲成人久久性| 国产精品免费一区二区三区在线| 母亲3免费完整高清在线观看| x7x7x7水蜜桃| 波多野结衣巨乳人妻| 悠悠久久av| 久久久久久久久大av| 欧美日韩乱码在线| 国产精品久久电影中文字幕| 欧美日韩综合久久久久久 | 一区二区三区免费毛片| 亚洲avbb在线观看| 一级黄色大片毛片| 亚洲精品粉嫩美女一区| 少妇的丰满在线观看| av女优亚洲男人天堂| 久久人妻av系列| 亚洲国产高清在线一区二区三| 黄片小视频在线播放| 最近在线观看免费完整版| 精品国内亚洲2022精品成人| 成人国产一区最新在线观看| 久久久久久久精品吃奶| 尤物成人国产欧美一区二区三区| 国内精品美女久久久久久| 有码 亚洲区| 三级毛片av免费| 看黄色毛片网站| 日本与韩国留学比较| 我要搜黄色片| 最新美女视频免费是黄的| 亚洲一区二区三区色噜噜| 搡老岳熟女国产| 中文字幕人妻丝袜一区二区| 露出奶头的视频| 午夜免费男女啪啪视频观看 | 中文字幕高清在线视频| 午夜视频国产福利| 老司机深夜福利视频在线观看| 天天一区二区日本电影三级| 亚洲精品在线美女| 此物有八面人人有两片| 午夜视频国产福利| 亚洲av五月六月丁香网| 亚洲一区二区三区色噜噜| 亚洲美女视频黄频| 超碰av人人做人人爽久久 | 欧美成人性av电影在线观看| 午夜福利视频1000在线观看| 免费看日本二区| 一区二区三区国产精品乱码| 在线a可以看的网站| 丁香欧美五月| 亚洲无线在线观看| 禁无遮挡网站| 亚洲国产中文字幕在线视频| 啦啦啦观看免费观看视频高清| 天天一区二区日本电影三级| av中文乱码字幕在线| 老汉色av国产亚洲站长工具| 最近最新免费中文字幕在线| 精品午夜福利视频在线观看一区| 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| 99精品久久久久人妻精品| 亚洲国产欧美网| 美女大奶头视频| 欧美成人一区二区免费高清观看| 人人妻人人看人人澡| 搡老岳熟女国产| 中国美女看黄片| 一级毛片女人18水好多| 在线观看一区二区三区| 男女视频在线观看网站免费| 色综合欧美亚洲国产小说| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 国产精品自产拍在线观看55亚洲| 色精品久久人妻99蜜桃| 日本三级黄在线观看| 内射极品少妇av片p| 免费看a级黄色片| 国产一区二区激情短视频| 在线观看日韩欧美| 精品一区二区三区av网在线观看| 免费av不卡在线播放| 国产高清激情床上av| 我要搜黄色片| 欧美丝袜亚洲另类 | 国产精品精品国产色婷婷| 亚洲人与动物交配视频| 欧美最黄视频在线播放免费| 12—13女人毛片做爰片一| 99精品欧美一区二区三区四区| 99久久无色码亚洲精品果冻| 婷婷丁香在线五月| 国产精品亚洲美女久久久| 男人舔奶头视频| 淫秽高清视频在线观看| bbb黄色大片| 人妻丰满熟妇av一区二区三区| 国产视频一区二区在线看| 久久亚洲精品不卡| 欧美黑人巨大hd| 午夜福利成人在线免费观看| 亚洲成av人片在线播放无| 91字幕亚洲| 国产成人福利小说| 99精品久久久久人妻精品| 成人一区二区视频在线观看| 身体一侧抽搐| 欧美黑人欧美精品刺激| 亚洲最大成人手机在线| 亚洲在线观看片| 国产精品98久久久久久宅男小说| 啦啦啦免费观看视频1| av视频在线观看入口| 久久精品国产99精品国产亚洲性色| 成人18禁在线播放| 国产成+人综合+亚洲专区| 久久人人精品亚洲av| 午夜福利免费观看在线| 在线观看日韩欧美| 在线观看日韩欧美| 最近最新中文字幕大全电影3| 欧美日韩国产亚洲二区| 人人妻,人人澡人人爽秒播| 久久久久久大精品| 99热6这里只有精品| 99久久99久久久精品蜜桃| 精品电影一区二区在线| 床上黄色一级片| 男人的好看免费观看在线视频| 看片在线看免费视频| 中文亚洲av片在线观看爽| 色吧在线观看| 性色avwww在线观看| 熟女少妇亚洲综合色aaa.| 午夜福利在线观看免费完整高清在 | 99久久99久久久精品蜜桃| 免费av毛片视频| 国产麻豆成人av免费视频| 在线观看66精品国产| 精品久久久久久久毛片微露脸| svipshipincom国产片| 91字幕亚洲| 亚洲精品久久国产高清桃花| 欧美日韩瑟瑟在线播放| 露出奶头的视频| 九色国产91popny在线| 精品久久久久久,| 午夜精品在线福利| 色综合婷婷激情| 人人妻人人澡欧美一区二区| 免费观看人在逋| 久久久久性生活片| a级一级毛片免费在线观看| 精品一区二区三区人妻视频| 少妇的逼好多水| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜添小说| 欧美+亚洲+日韩+国产| 最近视频中文字幕2019在线8| 国产视频内射| 久久久久免费精品人妻一区二区| 午夜影院日韩av| 男女之事视频高清在线观看| 熟女少妇亚洲综合色aaa.| 日本成人三级电影网站| 国产私拍福利视频在线观看| 好男人在线观看高清免费视频| 精品国内亚洲2022精品成人| 人人妻,人人澡人人爽秒播| а√天堂www在线а√下载| 校园春色视频在线观看| 欧美日韩国产亚洲二区| 最近最新免费中文字幕在线| 色尼玛亚洲综合影院| 999久久久精品免费观看国产| 亚洲av免费在线观看| 国产成人av教育| 少妇高潮的动态图| 国产精品久久久久久精品电影| 丝袜美腿在线中文| 老司机午夜福利在线观看视频| 国产熟女xx| 国产主播在线观看一区二区| 夜夜爽天天搞| 欧美日韩亚洲国产一区二区在线观看| 91九色精品人成在线观看| 国产真人三级小视频在线观看| 欧美在线一区亚洲| 日韩欧美在线二视频| 精品一区二区三区视频在线观看免费| 97超视频在线观看视频| 99精品久久久久人妻精品| 看黄色毛片网站| 亚洲美女黄片视频| 亚洲七黄色美女视频| 黄色丝袜av网址大全| 精品一区二区三区人妻视频| 国产精品爽爽va在线观看网站| 夜夜夜夜夜久久久久| 亚洲成人久久性| 免费看十八禁软件| 小说图片视频综合网站| 老鸭窝网址在线观看| 免费看光身美女| 国产精品爽爽va在线观看网站| 国产色婷婷99| 五月玫瑰六月丁香| 国产精品一区二区免费欧美| 国产精品,欧美在线| 亚洲人成网站在线播| 国产av在哪里看| 国产成人aa在线观看| 亚洲 欧美 日韩 在线 免费| 成人无遮挡网站| 国产精品国产高清国产av| 乱人视频在线观看| 9191精品国产免费久久| 99在线人妻在线中文字幕| 中文字幕人妻熟人妻熟丝袜美 | 在线免费观看不下载黄p国产 | 国产精品国产高清国产av| av在线天堂中文字幕| 一区二区三区国产精品乱码| 国产一区二区亚洲精品在线观看| 午夜久久久久精精品| 男人舔女人下体高潮全视频| 成熟少妇高潮喷水视频| 成人三级黄色视频| 成年女人毛片免费观看观看9| av在线蜜桃| 男人的好看免费观看在线视频| 亚洲av二区三区四区| 亚洲七黄色美女视频| 国产亚洲精品av在线| 久久久国产精品麻豆| 精品国内亚洲2022精品成人| 久久精品影院6| 精品一区二区三区av网在线观看| 日韩欧美三级三区| 久久久久精品国产欧美久久久| 日韩欧美精品v在线| 久久久久性生活片| 最好的美女福利视频网| 美女高潮的动态| 神马国产精品三级电影在线观看| 国产日本99.免费观看| 精品电影一区二区在线| 18禁黄网站禁片午夜丰满| 欧美一区二区精品小视频在线| 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 最近最新免费中文字幕在线| 国产精品98久久久久久宅男小说| 午夜福利在线观看吧| 白带黄色成豆腐渣| 中亚洲国语对白在线视频| 99久久九九国产精品国产免费| 欧美成人性av电影在线观看| 老司机深夜福利视频在线观看| 久久久精品欧美日韩精品| 欧美高清成人免费视频www| 国产精品三级大全| 少妇的逼水好多| 国产成人福利小说| a级一级毛片免费在线观看| 国产亚洲欧美在线一区二区| 中文字幕av在线有码专区| 国产一区二区在线av高清观看| 亚洲av成人精品一区久久| 国产精品免费一区二区三区在线| 国产成人a区在线观看| 黄片小视频在线播放| 九色国产91popny在线| 久久精品国产清高在天天线| 亚洲在线自拍视频| 天堂√8在线中文| 国产在线精品亚洲第一网站| www.999成人在线观看| 99热这里只有是精品50| АⅤ资源中文在线天堂| 99精品在免费线老司机午夜| 一个人免费在线观看的高清视频| 欧美日本视频| 国产高清有码在线观看视频| 国产精品一区二区三区四区免费观看 | 国产探花在线观看一区二区| 有码 亚洲区| 一级毛片女人18水好多| 露出奶头的视频| 国产亚洲精品综合一区在线观看| 制服丝袜大香蕉在线| 国产精品99久久久久久久久| 国产黄色小视频在线观看| 国产麻豆成人av免费视频| 大型黄色视频在线免费观看| 在线看三级毛片| 97碰自拍视频| 首页视频小说图片口味搜索| 3wmmmm亚洲av在线观看| svipshipincom国产片| 国产精品一及| 精品一区二区三区av网在线观看| 床上黄色一级片| 免费一级毛片在线播放高清视频| 亚洲av电影在线进入| 国产伦一二天堂av在线观看| av视频在线观看入口| h日本视频在线播放| 中亚洲国语对白在线视频| 亚洲精品一卡2卡三卡4卡5卡| 在线观看午夜福利视频| 精品国产亚洲在线| 久9热在线精品视频| 精品电影一区二区在线| 国产视频内射| 一区福利在线观看| 制服人妻中文乱码| 国产三级黄色录像| 91在线观看av| 女人十人毛片免费观看3o分钟| 亚洲欧美日韩东京热| 日日干狠狠操夜夜爽| 一个人看的www免费观看视频| 色精品久久人妻99蜜桃| 国产99白浆流出| 亚洲成人久久爱视频| 欧美又色又爽又黄视频| 老司机福利观看| 国产精品av视频在线免费观看| 一本精品99久久精品77| 国产精品一区二区三区四区久久| 99精品在免费线老司机午夜| 欧美一区二区亚洲| 757午夜福利合集在线观看| 老鸭窝网址在线观看| av福利片在线观看| 午夜福利视频1000在线观看| 97碰自拍视频| 中出人妻视频一区二区| 免费电影在线观看免费观看| 可以在线观看毛片的网站| 亚洲性夜色夜夜综合| 久久精品国产亚洲av涩爱 | 久久草成人影院| 国内揄拍国产精品人妻在线| 男女那种视频在线观看| 国产精品日韩av在线免费观看| 国产麻豆成人av免费视频| 欧美中文日本在线观看视频| 麻豆久久精品国产亚洲av| 老鸭窝网址在线观看| 国产亚洲精品一区二区www| 久久精品影院6| 国产精品永久免费网站| 村上凉子中文字幕在线| tocl精华| 亚洲专区国产一区二区| 国产高清有码在线观看视频| 午夜免费男女啪啪视频观看 | 18禁美女被吸乳视频| 久久久久久久久中文| 成年女人永久免费观看视频| 内射极品少妇av片p| 老鸭窝网址在线观看| 日本黄大片高清| 午夜激情福利司机影院| 3wmmmm亚洲av在线观看| 精品一区二区三区av网在线观看| 琪琪午夜伦伦电影理论片6080| 国产精品99久久久久久久久| 国产精品久久视频播放| 国产高潮美女av| 女人被狂操c到高潮| avwww免费| 夜夜躁狠狠躁天天躁| 国产精品爽爽va在线观看网站| 亚洲最大成人中文| 精品久久久久久久末码| 一区二区三区高清视频在线| 蜜桃亚洲精品一区二区三区| 宅男免费午夜| 国产亚洲精品综合一区在线观看| 成人午夜高清在线视频| 国产一级毛片七仙女欲春2| bbb黄色大片| 91av网一区二区| 日韩欧美国产在线观看| 欧美一区二区国产精品久久精品| 亚洲国产精品合色在线| 国产亚洲av嫩草精品影院| 亚洲精品亚洲一区二区| 欧洲精品卡2卡3卡4卡5卡区| 色吧在线观看| 欧美日本视频| 中亚洲国语对白在线视频| 91久久精品电影网| 国产精品亚洲av一区麻豆| 国产蜜桃级精品一区二区三区| 脱女人内裤的视频| 性色av乱码一区二区三区2| 中文亚洲av片在线观看爽| 一进一出抽搐动态| 热99在线观看视频| 国产精品电影一区二区三区| 精品免费久久久久久久清纯| 精品一区二区三区视频在线观看免费| 无人区码免费观看不卡| av欧美777| 亚洲在线自拍视频| 一夜夜www| 国产精品亚洲一级av第二区| 午夜激情福利司机影院| 在线观看一区二区三区| 人妻丰满熟妇av一区二区三区| 成人特级av手机在线观看| 久久久久国产精品人妻aⅴ院| 欧美乱妇无乱码| 久久久久亚洲av毛片大全| 一卡2卡三卡四卡精品乱码亚洲| 天堂动漫精品| 国产午夜精品论理片| a在线观看视频网站| 日韩精品青青久久久久久| 欧美极品一区二区三区四区| 五月玫瑰六月丁香| 搞女人的毛片| 欧美色视频一区免费| 国产野战对白在线观看| 国产aⅴ精品一区二区三区波| 丰满乱子伦码专区| 欧美最新免费一区二区三区 | www国产在线视频色| 国产精品免费一区二区三区在线| 免费av观看视频| 亚洲精品久久国产高清桃花| 久久欧美精品欧美久久欧美| 国内久久婷婷六月综合欲色啪| 亚洲内射少妇av| 手机成人av网站| 99在线视频只有这里精品首页| 亚洲自拍偷在线| av黄色大香蕉| 最近视频中文字幕2019在线8| 亚洲av五月六月丁香网| 无人区码免费观看不卡| 久久婷婷人人爽人人干人人爱| 人妻丰满熟妇av一区二区三区| 在线观看美女被高潮喷水网站 | 91字幕亚洲| 内地一区二区视频在线| 成年女人永久免费观看视频| 精品福利观看| 老司机午夜十八禁免费视频| 黄色女人牲交| 18+在线观看网站| 欧美黑人巨大hd| 欧美又色又爽又黄视频| 国产中年淑女户外野战色| 午夜福利18| 成年免费大片在线观看| 日韩欧美国产一区二区入口| 别揉我奶头~嗯~啊~动态视频| 夜夜看夜夜爽夜夜摸| 日本一二三区视频观看| 我的老师免费观看完整版| 亚洲中文字幕日韩| 亚洲精华国产精华精| 精品国产亚洲在线| 欧美成狂野欧美在线观看| 三级国产精品欧美在线观看| 国产欧美日韩一区二区精品| 精品一区二区三区人妻视频| 热99re8久久精品国产| 此物有八面人人有两片| 色视频www国产| 日韩 欧美 亚洲 中文字幕| 99热这里只有精品一区| 国产高清激情床上av| 午夜久久久久精精品| 少妇的逼好多水| 欧美高清成人免费视频www| 免费观看人在逋| 欧美最新免费一区二区三区 | 成年女人永久免费观看视频| 99久久99久久久精品蜜桃| 欧美成人性av电影在线观看| av天堂中文字幕网| 桃红色精品国产亚洲av| 12—13女人毛片做爰片一| 成年免费大片在线观看| 久久精品国产自在天天线| 国内久久婷婷六月综合欲色啪| 深夜精品福利| 国产欧美日韩精品一区二区| 亚洲精品影视一区二区三区av| 国产乱人视频| 观看免费一级毛片| 亚洲欧美一区二区三区黑人| 国内少妇人妻偷人精品xxx网站| 日韩欧美在线二视频| 偷拍熟女少妇极品色| 亚洲一区二区三区色噜噜| 蜜桃久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 国产一区二区三区视频了| 亚洲avbb在线观看| 欧美性猛交╳xxx乱大交人| 夜夜看夜夜爽夜夜摸| 校园春色视频在线观看| a级毛片a级免费在线| 国产一级毛片七仙女欲春2| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av电影在线进入| 精品国产超薄肉色丝袜足j|