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

    基于深度學(xué)習(xí)的地震速度譜自動(dòng)拾取研究

    2022-12-03 09:36:18崔家豪楊平王洪強(qiáng)邊策胡揚(yáng)潘海俠
    地球物理學(xué)報(bào) 2022年12期
    關(guān)鍵詞:人工能量自動(dòng)

    崔家豪,楊平,王洪強(qiáng),邊策,胡揚(yáng),潘海俠*

    1 北京航空航天大學(xué)軟件學(xué)院,北京 100191 2 中國(guó)石油天然氣集團(tuán)東方地球物理勘探有限責(zé)任公司,河北涿州 072750

    0 引言

    疊加速度分析是常規(guī)地震數(shù)據(jù)處理中的重要環(huán)節(jié)之一,也是層析反演、疊前偏移、阻抗反演和全波形反演等初始速度模型建立的基礎(chǔ)(潘宏勛等,2008;Taner and Koehler,2012).由于拾取疊加速度需要考慮地震噪聲、多次波、繞射波和構(gòu)造復(fù)雜性等諸多因素的影響(戴曉峰等,2020),因此拾取工作需要具有豐富地震處理經(jīng)驗(yàn)和地質(zhì)學(xué)術(shù)背景的人員來(lái)完成.隨著石油勘探技術(shù)的不斷進(jìn)步,地震數(shù)據(jù)的規(guī)模和復(fù)雜性在近些年來(lái)一直在不斷增長(zhǎng),而傳統(tǒng)人工能量團(tuán)識(shí)別與疊加速度拾取的工作模式容易受拾取標(biāo)準(zhǔn)不一致及人為因素干擾而產(chǎn)生誤差,難以滿(mǎn)足實(shí)際生產(chǎn)過(guò)程中對(duì)拾取精度和效率的要求.如何從速度譜中自動(dòng)拾取疊加速度成為該行業(yè)面臨的挑戰(zhàn)之一.當(dāng)前速度譜自動(dòng)拾取的方法主要包括優(yōu)化搜索方法和機(jī)器學(xué)習(xí)方法.

    基于優(yōu)化搜索的方法是在層速度約束條件下,采用迭代或搜索算法尋找地質(zhì)上合理的速度模型.算法以層速度作為初始模型,通過(guò)反演的方式來(lái)得到最優(yōu)的速度模型.Toldi(1989)最先提出基于速度譜中疊加能量的自動(dòng)速度拾取方法.Zhang(1991)將Viterbi算子引入到地震資料處理中并解決了關(guān)于拾取的一些問(wèn)題,林年添等(2013)對(duì)該方法進(jìn)行了改進(jìn),將用于圖像處理的光順技術(shù)引用到速度譜能量團(tuán)曲面的處理環(huán)節(jié),借以提高速度譜的分辨率,然后再通過(guò)Viterbi算法自動(dòng)搜索及獲取最優(yōu)解.Lumley(1997)和陳志德等(2002)提出了一種利用蒙特卡洛方法對(duì)初始速度加以隨機(jī)擾動(dòng)來(lái)自動(dòng)拾取層速度.Almarzoug和Ahmed(2012)將速度拾取視為變分問(wèn)題,通過(guò)有限差分法求解最優(yōu)拾取軌跡.Leite和Vieira(2019)將非線(xiàn)性?xún)?yōu)化方法應(yīng)用于速度拾取中,在滿(mǎn)足邊界條件的前提下設(shè)置合理的搜索帶寬,然后運(yùn)行非線(xiàn)性函數(shù)智能拾取.Li和Biondi(2009)以及Oliveira等(2013)分別將模擬退火算法和遺傳算法應(yīng)用于速度拾取中.這些方法在獲取相對(duì)簡(jiǎn)單地質(zhì)體的速度時(shí)已經(jīng)取得了一定的效果,但是對(duì)于橫向速度突變的地區(qū)拾取準(zhǔn)確性相對(duì)不高,且需要過(guò)多設(shè)置算法超參數(shù)、初始化模型,這并不能實(shí)現(xiàn)完全的自動(dòng)取速過(guò)程.

    基于機(jī)器學(xué)習(xí)的方法是通過(guò)學(xué)習(xí)速度譜中能量團(tuán)的特征來(lái)獲取疊加速度.在近期的研究中,Huang和Yang(2015)采用Hopfield神經(jīng)網(wǎng)絡(luò)(Hopfield Neural Network,HNN)對(duì)速度譜圖像中的速度進(jìn)行自動(dòng)拾取.Smith(2017)和Song等(2018)將基于無(wú)監(jiān)督學(xué)習(xí)的聚類(lèi)分析技術(shù)應(yīng)用于速度譜自動(dòng)拾取過(guò)程中.該方法依據(jù)不同的屬性將速度譜劃分為若干個(gè)類(lèi)別,每個(gè)類(lèi)別中的數(shù)據(jù)具有一些共同的特征,通常根據(jù)某種距離度量進(jìn)行分組,然后將聚類(lèi)中心映射為正確的“時(shí)間-速度”對(duì).K-means聚類(lèi)是最常用的聚類(lèi)方法之一(Hung et al.,2013;Ahmad and Hashmi,2016),該算法使用歐幾里得距離度量將每個(gè)數(shù)據(jù)點(diǎn)分配到最近的聚類(lèi)中心.對(duì)于每個(gè)聚類(lèi)中心,計(jì)算分配給它的所有數(shù)據(jù)點(diǎn)的平均值,這個(gè)平均值成為聚類(lèi)中心的新位置,每個(gè)聚類(lèi)的聚類(lèi)中心即為我們所要拾取的疊加速度點(diǎn).然而傳統(tǒng)的K-means聚類(lèi)方法容易受到離群點(diǎn)的影響,使得聚類(lèi)中心偏離主要數(shù)據(jù)區(qū)域.Liu等(2018)提出了一種相似加權(quán)K-means聚類(lèi)的新方法,該方法對(duì)不同信噪比的數(shù)據(jù)點(diǎn)賦予不同的權(quán)值,使聚類(lèi)中心更接近權(quán)重大的數(shù)據(jù)點(diǎn).由于不同速度譜所需拾取的能量團(tuán)個(gè)數(shù)不一致,Chen和Schuster(2018)以及王迪等(2021)新提出了一種自底向上的策略來(lái)實(shí)現(xiàn)K-means迭代過(guò)程,該方法可以逐步迭代找到聚類(lèi)類(lèi)別最少的聚類(lèi)分布.除此之外,加速聚類(lèi)算法(Zhang and Lu,2016)以及DBSCAN聚類(lèi)(Schubert et al.,2017;Bin Waheed et al.,2019)也被應(yīng)用于智能速度譜解釋過(guò)程中.然而傳統(tǒng)機(jī)器學(xué)習(xí)方法也存在一些問(wèn)題,如易受到噪聲和離群點(diǎn)影響,部分拾取點(diǎn)偏離主要數(shù)據(jù)區(qū)域,進(jìn)而影響到速度拾取精度.近年來(lái),隨著計(jì)算機(jī)性能的提高,應(yīng)用深度學(xué)習(xí)方法拾取疊加速度的研究越來(lái)越多.Ma等(2018)以及Biswas(2018)和Fabien-Ouellet等(2020)分別采用卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Network,CNN)和循環(huán)神經(jīng)網(wǎng)絡(luò)(Recurrent Neural Network,RNN)來(lái)拾取速度譜中的疊加速度.在信噪比較高的工區(qū)中,以上深度學(xué)習(xí)方法可以提高疊加速度的拾取精度,但對(duì)于信噪比低的工區(qū),速度拾取的準(zhǔn)確性不足且受到噪聲的干擾比較大.

    為了有效解決速度拾取準(zhǔn)確率不高的問(wèn)題,本文提出了一種深度學(xué)習(xí)的方法.該方法將速度譜視為圖像進(jìn)行處理,結(jié)合目標(biāo)任務(wù)特點(diǎn)改進(jìn)了基于卷積神經(jīng)網(wǎng)絡(luò)的FCOS目標(biāo)檢測(cè)模型(Tian et al.,2019)用于速度譜中疊加速度的自動(dòng)拾取.與傳統(tǒng)的目標(biāo)檢測(cè)網(wǎng)絡(luò)不同,改進(jìn)后的FCOS檢測(cè)器具有以下優(yōu)點(diǎn):①顯著減少了設(shè)計(jì)參數(shù)的數(shù)目,一方面提升了對(duì)能量團(tuán)檢測(cè)的識(shí)別效果,同時(shí)構(gòu)建的混合損失函數(shù)增強(qiáng)了模型對(duì)能量團(tuán)檢測(cè)的敏感度,使得能量團(tuán)識(shí)別和疊加速度拾取僅通過(guò)深度學(xué)習(xí)模型即可完成,保證了識(shí)別和拾取處理結(jié)果的一致性.②由于預(yù)測(cè)目標(biāo)為疊加速度點(diǎn),每個(gè)點(diǎn)在速度譜上的邊界框具有固定的大小,使用改進(jìn)后的FCOS可直接對(duì)邊界框進(jìn)行回歸,從而避免了復(fù)雜的計(jì)算,極大地減小了內(nèi)存占用并使得訓(xùn)練速度更快.

    另外,當(dāng)?shù)卣鹳Y料受多次波干擾較強(qiáng)時(shí),除去一次波反射能量團(tuán)以外,多次波也能形成能量團(tuán),尤其當(dāng)多次波與一次波速度相差較大時(shí),速度譜能量團(tuán)聚焦特征變差,在同一時(shí)刻出現(xiàn)多個(gè)能量團(tuán).由于能量團(tuán)較為分散,改進(jìn)后的FCOS模型只能夠拾取到近地表淺層的疊加速度,對(duì)于地下深層區(qū)間的疊加速度僅通過(guò)單一的FCOS模型未能有效識(shí)別,因此在處理低信噪比探區(qū)時(shí)本文在FCOS后又加入基于深度神經(jīng)網(wǎng)絡(luò)(DNN)的線(xiàn)性回歸模型以擬合出全局速度曲線(xiàn).我們從方法原理及模型訓(xùn)練與測(cè)試等角度闡述了本文方法的操作流程,并應(yīng)用Ma等(2018)提出的傳統(tǒng)的CNN方法與本文方法處理相同的數(shù)據(jù)集來(lái)對(duì)比驗(yàn)證本文方法的有效性.模型和實(shí)際地震數(shù)據(jù)的測(cè)試結(jié)果表明該網(wǎng)絡(luò)模型準(zhǔn)確度較高,模型魯棒性強(qiáng),能夠較好的實(shí)現(xiàn)對(duì)復(fù)雜工區(qū)的速度譜中疊加速度的自動(dòng)拾取.

    1 方法原理

    1.1 模型整體結(jié)構(gòu)設(shè)計(jì)

    結(jié)合目標(biāo)任務(wù)特點(diǎn)設(shè)計(jì)的滿(mǎn)足上述需求的網(wǎng)絡(luò)模型的整體設(shè)計(jì)流程如圖1所示.該模型分為卷積神經(jīng)網(wǎng)絡(luò)和線(xiàn)性回歸網(wǎng)絡(luò)兩部分.首先對(duì)包含速度譜圖像及其標(biāo)簽數(shù)據(jù)的數(shù)據(jù)集進(jìn)行處理使其滿(mǎn)足網(wǎng)絡(luò)的輸入要求,然后將標(biāo)簽數(shù)據(jù)集輸入到基于CNN的FCOS模型中提取能量團(tuán)的特征并計(jì)算能量團(tuán)的位置.該模型能夠?qū)μ崛〉降亩喑叨忍卣鬟M(jìn)行融合,從而提升能量團(tuán)檢測(cè)精度.卷積神經(jīng)網(wǎng)絡(luò)階段輸出的中間結(jié)果為預(yù)測(cè)的“時(shí)間-速度”對(duì)序列.受多次波及噪聲等影響,速度譜圖像中存在信噪比較低且能量團(tuán)不夠聚焦的區(qū)域,在該區(qū)域FCOS模型只能完成淺層區(qū)域內(nèi)疊加速度的自動(dòng)拾取.針對(duì)這一情況我們引入人工判別機(jī)制,當(dāng)預(yù)測(cè)結(jié)果滿(mǎn)足拾取要求時(shí),可直接輸出中間結(jié)果為最終的智能拾取結(jié)果;當(dāng)預(yù)測(cè)結(jié)果不能滿(mǎn)足拾取要求時(shí),需要加入基于DNN的線(xiàn)性回歸模型來(lái)擬合出能量團(tuán)聚焦特性較差區(qū)域的速度曲線(xiàn).本章2、3小節(jié)分別對(duì)FCOS和線(xiàn)性回歸網(wǎng)絡(luò)展開(kāi)詳細(xì)介紹.

    圖1 模型整體設(shè)計(jì)流程圖Fig.1 Overall design flow chart of the model

    1.2 FCOS網(wǎng)絡(luò)結(jié)構(gòu)與工作原理

    本文針對(duì)能量團(tuán)結(jié)構(gòu)復(fù)雜、體積小,在速度譜中難以檢測(cè)等難點(diǎn)對(duì)FCOS網(wǎng)絡(luò)進(jìn)行了改進(jìn),改進(jìn)后的模型更適用于識(shí)別和定位速度譜中能量團(tuán)的位置.如圖2所示,首先通過(guò)一個(gè)卷積層和一個(gè)最大池化層將原始圖片轉(zhuǎn)化為特征圖.卷積操作(Convolution)采用卷積核提取能量團(tuán)的特征,該操作能突破傳統(tǒng)機(jī)器學(xué)習(xí)方法的限制,通過(guò)共享多個(gè)卷積核的權(quán)值,來(lái)學(xué)習(xí)每張速度譜圖片的局部信息,并用于構(gòu)建抽象特征圖譜.卷積核共享特性大大降低了訓(xùn)練網(wǎng)絡(luò)所需的參數(shù)量.最大池化操作(Maxpooling)將卷積處理后生成的特征圖劃分為若干個(gè)矩形區(qū)域,對(duì)每個(gè)子區(qū)域輸出最大值,從而有效的壓縮數(shù)據(jù)空間的大小,避免在模型訓(xùn)練過(guò)程中出現(xiàn)過(guò)擬合.然后將最大池化操作產(chǎn)生的特征圖輸入到4級(jí)降采樣(Down-sample)模塊中,每一級(jí)包含有若干個(gè)降采樣過(guò)程,用于逐級(jí)提取數(shù)據(jù)特征.降采樣模塊的設(shè)計(jì)可以降低特征的維度并保留有效信息,一定程度上減少了計(jì)算量,有效解決了傳統(tǒng)CNN網(wǎng)絡(luò)隨層數(shù)的加深容易出現(xiàn)的梯度爆炸或梯度消失等問(wèn)題.選用1×1大小的卷積核對(duì)第4級(jí)降采樣模塊輸出的特征圖進(jìn)行卷積操作,得到P5特征層.然后使用步長(zhǎng)為2大小為3×3的卷積核對(duì)P5特征層進(jìn)行兩次降采樣得到P6、P7特征層.同時(shí)對(duì)P5特征層進(jìn)行兩次上采樣(Up-sample),上采樣階段通過(guò)轉(zhuǎn)置卷積操作(Deconvolution)來(lái)提升特征圖的大小.分別將上采樣后得到的高層特征與對(duì)應(yīng)的殘差網(wǎng)絡(luò)模塊相同放大倍數(shù)的低層特征進(jìn)行融合(Concatenate)得到P3、P4特征層,融合操作可將低層次的細(xì)節(jié)信息與高層次的語(yǔ)義信息結(jié)合到一起,從而使特征提取更完備.通過(guò)上述操作可在5個(gè)不同尺寸的特征層上進(jìn)行目標(biāo)檢測(cè),從而實(shí)現(xiàn)對(duì)較小目標(biāo)與多尺度目標(biāo)的準(zhǔn)確識(shí)別.本文中卷積層和上采樣層均采用ReLU非線(xiàn)性激活函數(shù),并在卷積階段對(duì)圖像的邊緣進(jìn)行填充,使邊緣信息提取更加充分.

    圖2 FCOS結(jié)構(gòu)圖左側(cè)降采樣模塊用于逐層提取特征;中間部分進(jìn)行多尺度特征融合;右側(cè)是多個(gè)共享權(quán)重的FCOS檢測(cè)頭.參數(shù)含義:通道數(shù)Channel和降采樣率Down-sample Rate(標(biāo)注于每層左側(cè)),H×W×C表示圖像的高×寬×通道數(shù).Fig.2 FCOS structure diagramThe left Up-sample module is used to extract features layer by layer;Multi-scale feature fusion was performed in the middle part;On the right are multiple FCOS detection heads with shared weights.Parameter meaning:Channel number and Down-sample Rate (marked on the left of each layer).H×W×C represents the height × width × number of channels of the image.

    FCOS的檢測(cè)模塊總體分為類(lèi)別判斷分支(Classification)與位置回歸分支(Regression),它作用在特征融合后的多個(gè)不同層級(jí)的特征圖上,同時(shí)這些檢測(cè)模塊的權(quán)值共享.模型的損失函數(shù)如(1)式展示:

    (1)

    (2)

    (3)

    整體損失函數(shù)由類(lèi)別損失Lcls與位置損失Lreg組成.對(duì)于類(lèi)別損失,采用交叉熵?fù)p失(focal loss)的設(shè)計(jì),其含義為目標(biāo)與預(yù)測(cè)值之間概率分布的距離,交叉熵越小,概率分布越接近.它能夠較好地解決樣本不均衡問(wèn)題.對(duì)于位置損失,采用交并比損失(IOU loss)的設(shè)計(jì),其含義為預(yù)測(cè)框與標(biāo)注框的交集與并集的比值,它能夠編碼空間相關(guān)性信息輔助更加準(zhǔn)確地定位.

    1.3 線(xiàn)性回歸網(wǎng)絡(luò)結(jié)構(gòu)與工作原理

    在不考慮速度倒轉(zhuǎn)的情況下疊加速度會(huì)隨著時(shí)間的增大而增大,在處理能量團(tuán)聚焦特性較差的速度譜時(shí),經(jīng)驗(yàn)較為豐富的專(zhuān)家一般會(huì)根據(jù)拾取到的淺層疊加速度的線(xiàn)性趨勢(shì)來(lái)完成對(duì)較深地層中疊加速度的拾取.而僅僅依靠單個(gè)目標(biāo)檢測(cè)網(wǎng)絡(luò)很難處理復(fù)雜地質(zhì)體條件下的能量團(tuán)識(shí)別問(wèn)題.為了完善本文的算法方案,我們?cè)贔COS模型的基礎(chǔ)上設(shè)計(jì)了一個(gè)線(xiàn)性回歸網(wǎng)絡(luò),其本質(zhì)為深度神經(jīng)網(wǎng)絡(luò)(DNN)模型.該網(wǎng)絡(luò)主要由輸入層、隱藏層和輸出層等層級(jí)結(jié)構(gòu)組成(圖3),可通過(guò)學(xué)習(xí)淺層區(qū)域內(nèi)疊加速度的線(xiàn)性趨勢(shì)來(lái)完成對(duì)較深地層中疊加速度的預(yù)測(cè).網(wǎng)絡(luò)的輸入為FCOS模型自動(dòng)拾取的“時(shí)間-速度”對(duì),隱藏層由若干個(gè)全連接層構(gòu)成,全連接層的作用為擬合經(jīng)插值處理后形成的全局速度曲線(xiàn),其前向計(jì)算公式為

    (4)

    (5)

    使用ReLU激活函數(shù)可以克服訓(xùn)練過(guò)程中出現(xiàn)的梯度消失等問(wèn)題,加快訓(xùn)練速度.損失函數(shù)采用l2Loss,其形式如下:

    Loss(z[l],y)=0.5×(z[l]-y)2,

    (6)

    式中y表示網(wǎng)絡(luò)的輸入值,z[l]表示第l層輸出的數(shù)值矩陣.全連接層的反向傳播過(guò)程如下:

    dx=?Loss/?u,

    (7)

    其中u表示任意變量,包括z,x,W,b等.通過(guò)一系列前向和反向傳播過(guò)程來(lái)不斷更新網(wǎng)絡(luò)模型的參數(shù),直到Loss層的梯度保持平穩(wěn)不再下降.網(wǎng)絡(luò)最后一層采用sigmoid激活函數(shù)進(jìn)行二分類(lèi)任務(wù),使輸出轉(zhuǎn)換為0~1之間的概率分布,公式表示如下:

    (8)

    圖3 線(xiàn)性回歸網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.3 Linear fitting network structure diagram

    2 實(shí)驗(yàn)與結(jié)果分析

    為了驗(yàn)證提出方法的可行性,采用Marmousi模型數(shù)據(jù)和兩個(gè)實(shí)際探區(qū)數(shù)據(jù)對(duì)本文方法進(jìn)行測(cè)試.在訓(xùn)練階段采用小批量梯度下降法,我們?cè)O(shè)置batchsize為2,其含義為進(jìn)行一次訓(xùn)練所抓取的數(shù)據(jù)集樣本數(shù)量.訓(xùn)練過(guò)程中以一個(gè)batchsize為最小結(jié)構(gòu)單元進(jìn)行損失函數(shù)計(jì)算和參數(shù)梯度更新,并采用雙線(xiàn)程計(jì)算以提升程序的運(yùn)行效率.實(shí)驗(yàn)選取Adam算法來(lái)解決參數(shù)更新過(guò)程中的優(yōu)化問(wèn)題,基礎(chǔ)學(xué)習(xí)率設(shè)置為0.01,共計(jì)訓(xùn)練24個(gè)輪次,分別在第8輪和第11輪進(jìn)行學(xué)習(xí)率衰減,權(quán)重衰減系數(shù)設(shè)置為0.0001.實(shí)驗(yàn)應(yīng)用Pytorch深度學(xué)習(xí)框架,并使用GPU來(lái)加速訓(xùn)練過(guò)程.

    本文使用平均準(zhǔn)確度均值(mAP)來(lái)評(píng)估目標(biāo)檢測(cè)器的有效性,mAP的計(jì)算方法如公式(9)—(10)所示:

    (9)

    (10)

    其中P表示精準(zhǔn)度,其含義為所有預(yù)測(cè)為正樣本的樣本中正樣本所占比例;R表示召回率,其含義為正樣本被預(yù)測(cè)正確的比例.TP表示既是正樣本又被預(yù)測(cè)為正樣本的個(gè)數(shù);FP表示負(fù)樣本被預(yù)測(cè)為正樣本的個(gè)數(shù);FN表示正樣本被預(yù)測(cè)為負(fù)樣本的個(gè)數(shù).精準(zhǔn)度和召回率始終處于矛盾的狀態(tài),當(dāng)其中一個(gè)增高,必然意味著另一個(gè)降低.以PR為橫縱坐標(biāo)建立坐標(biāo)系,AP(Average Precision)計(jì)算方法為PR曲線(xiàn)下覆蓋的面積,所有類(lèi)別的AP進(jìn)行平均,得到的即為mAP (mean Average Precision).對(duì)于目標(biāo)檢測(cè)器的評(píng)估方法,其mAP越高,檢測(cè)方法越有效.

    2.1 模型數(shù)據(jù)

    首先應(yīng)用二維的Marmousi模型對(duì)本文方法進(jìn)行測(cè)試,從該模型中等間隔抽取了900道作為本次測(cè)試的真實(shí)層速度場(chǎng)(圖4).選用主頻為35 Hz的雷克子波通過(guò)聲波波動(dòng)方程有限差分法對(duì)該速度模型正演得到一系列CMP道集,然后對(duì)CMP道集進(jìn)行速度掃描,共計(jì)算得到900張速度譜圖像.再經(jīng)由專(zhuān)家來(lái)人工拾取速度譜中的疊加速度,并從某商業(yè)軟件中導(dǎo)出統(tǒng)一大小的速度譜圖像及其對(duì)應(yīng)的疊加速度做為標(biāo)簽數(shù)據(jù).將數(shù)據(jù)集按8∶1∶1分為訓(xùn)練集、驗(yàn)證集和測(cè)試集三部分,用于訓(xùn)練速度譜識(shí)別模型并對(duì)該方法進(jìn)行檢驗(yàn).在訓(xùn)練階段使用梯度反向傳播的方法來(lái)不斷地更新優(yōu)化Loss函數(shù),并使用mAP來(lái)衡量檢測(cè)模型在驗(yàn)證集上精度表現(xiàn)的優(yōu)劣,實(shí)驗(yàn)結(jié)果如圖5所示.

    圖4 層速度場(chǎng)模型Fig.4 Interval velocity model

    從圖5中可以看出通過(guò)梯度反向傳播的方式來(lái)不斷地更新和優(yōu)化Loss后,本文方法相較傳統(tǒng)CNN方法得到的Loss曲線(xiàn)更為收斂,且應(yīng)用本文方法訓(xùn)練的目標(biāo)檢測(cè)器在Marmousi數(shù)據(jù)集上的mAP精度最高為74.1,而傳統(tǒng)CNN方法在Marmousi數(shù)據(jù)集上的mAP精度最高僅為54.5.上述結(jié)果驗(yàn)證了本文方法在速度譜能量團(tuán)檢測(cè)上的有效性.

    選擇訓(xùn)練好的FCOS模型對(duì)測(cè)試集中的速度譜圖片進(jìn)行預(yù)測(cè).圖6(a,b)分別展示了應(yīng)用傳統(tǒng)CNN方法及本文方法拾取的疊加速度和真實(shí)均方根速度的對(duì)比.從圖中可以看出在淺層區(qū)間不存在多次波的干擾,能量團(tuán)集中在一個(gè)較窄的速度范圍內(nèi),聚焦較好.而在2.5 s以后多次波干擾較大,多次波速度與一次波速度疊加在一起導(dǎo)致能量團(tuán)逐漸發(fā)散,在同一時(shí)刻出現(xiàn)多個(gè)能量團(tuán).受深層多次波干擾的影響,傳統(tǒng)CNN方法(圖6a)僅能拾取到0.4~3.0 s時(shí)間范圍內(nèi)的速度,且該方法在2.3~3.0 s區(qū)間內(nèi)拾取結(jié)果遠(yuǎn)偏離真實(shí)值;而FCOS方法(圖6b)得到的速度曲線(xiàn)與真實(shí)值非常近似,證明該方法能夠有效避免多次波的影響,使拾取結(jié)果具有更高的精度.圖6(c—f)分別展示了Marmousi模型中第400道動(dòng)校正前后的CMP道集,從動(dòng)校正后的道集中可以看出,F(xiàn)COS方法(圖6e)拾取的速度拉平了淺層的同相軸,而深層區(qū)間受多次波剩余時(shí)差較大影響,導(dǎo)致同相軸隨偏移距增大具有下拉現(xiàn)象;而傳統(tǒng)CNN方法(圖6d)在2.3~3.0 s區(qū)間選取比疊加速度小的多次波速度進(jìn)行動(dòng)校正,導(dǎo)致同相軸被錯(cuò)誤的拉平.將FCOS方法得到的動(dòng)校正后道集與真實(shí)速度動(dòng)校正后道集(圖6f)進(jìn)行對(duì)比后發(fā)現(xiàn)二者主要反射同相軸拉平程度相一致,也驗(yàn)證了本文算法的有效性.

    圖5 Marmousi模型Loss曲線(xiàn)圖(a)及mAP曲線(xiàn)圖(b)Fig.5 Marmousi model loss curve (a) and mAP curve (b)

    圖6 Marmousi模型預(yù)測(cè)結(jié)果(a) 傳統(tǒng)CNN方法拾取的速度曲線(xiàn);(b) 本文方法拾取的速度曲線(xiàn);(c) CMP400道集;(d) 傳統(tǒng)CNN方法拾取速度動(dòng)校正;(e) 本文方法拾取速度動(dòng)校正;(f) 人工拾取速度動(dòng)校正.Fig.6 Prediction results of Marmousi model(a) Velocity curve picked up by traditional CNN method;(b) Velocity curves picked up by the method in this paper;(c) CMP400 gather;(d) NMO gather corrected of the stacked velocity picked up by traditional CNN method;(e) NMO gather corrected of the stacked velocity picked up by the method in this paper;(f) NMO gather corrected of the stacked velocity picked up by manual method.

    圖7 Marmousi模型疊加速度場(chǎng)(a) 人工方法拾取的速度場(chǎng);(b) 傳統(tǒng)CNN方法自動(dòng)拾取的速度場(chǎng);(c) 本文方法自動(dòng)拾取的速度場(chǎng).Fig.7 Stacking velocity field of Marmousi model(a) The velocity field picked up by manual method;(b) Velocity field automatically picked up by traditional CNN method;(c) The velocity field automatically picked by the method in this paper.

    圖8 Marmousi模型疊加剖面(a) 人工方法拾取的速度場(chǎng)疊加;(b) 傳統(tǒng)CNN方法自動(dòng)拾取的速度場(chǎng)疊加;(c) 本文方法自動(dòng)拾取的速度場(chǎng)疊加.Fig.8 Stacked sections of Marmousi model(a) The Stacked section of velocity field picked up by manual method;(b) The Stacked section of velocity field picked up by traditional CNN method;(c) The Stacked section of velocity field picked up by the method in this paper.

    圖9 實(shí)際探區(qū)數(shù)據(jù)Loss曲線(xiàn)圖(a)及mAP曲線(xiàn)圖(b)Fig.9 Actual exploration area data loss curve (a) and mAP curve (b)

    圖10 實(shí)際數(shù)據(jù)CMP道集(a) 傳統(tǒng)CNN方法拾取的速度曲線(xiàn);(b) 本文方法拾取的速度曲線(xiàn);(c) 2560測(cè)線(xiàn)CMP2240道集;(d) 傳統(tǒng)CNN方法拾取速度動(dòng)校正;(e) 本文方法拾取速度動(dòng)校正;(f) 人工拾取速度動(dòng)校正.Fig.10 CMP gather of actual data(a) Velocity curve picked up by traditional CNN method;(b) Velocity curves picked up by the method in this paper;(c) CMP2240 gather in the2560Line;(d) NMO gather corrected of the stacked velocity picked up by traditional CNN method ;(e) NMO gather corrected of the stacked velocity picked up by the method in this paper;(f) NMO gather corrected of the stacked velocity picked up by manual method.

    對(duì)整個(gè)Marmousi模型數(shù)據(jù)集進(jìn)行預(yù)測(cè),并將本文方法自動(dòng)拾取的疊加速度場(chǎng)(圖7c)與傳統(tǒng)CNN方法得到的疊加速度場(chǎng)(圖7b)和人工方法得到的速度場(chǎng)(圖7a)進(jìn)行對(duì)比.總體上看,本文方法拾取的速度場(chǎng)在構(gòu)造刻畫(huà)與連續(xù)性上與人工方法得到的速度場(chǎng)基本相同;而傳統(tǒng)CNN方法拾取的速度場(chǎng)在深層區(qū)間(2.5~4.0 s)內(nèi)速度值偏小,且橫向連續(xù)性較差.以CMP間隔為10將動(dòng)校正后的道集進(jìn)行同相疊加形成疊加剖面.與傳統(tǒng)CNN方法得到的疊加剖面(圖8b)相比,使用本文方法得到的疊加剖面(圖8c)在淺層成像更加清晰,對(duì)構(gòu)造細(xì)節(jié)刻畫(huà)地更為準(zhǔn)確,而傳統(tǒng)CNN方法得到的疊加剖面受多次波等因素的影響,在淺層部分區(qū)間同相軸橫向波動(dòng)較大,在3.0 s以后未見(jiàn)明顯地震反射(如紅色框所示);與人工拾取拾取速度場(chǎng)得到的剖面相比(圖8a),本文方法得到的疊加剖面在構(gòu)造連續(xù)性與清晰度上與人工方法基本一致,這進(jìn)一步證明了本文方法的有效性.

    2.2 實(shí)際數(shù)據(jù)

    基于模型數(shù)據(jù)的實(shí)驗(yàn)結(jié)果,分別選用不同信噪比的實(shí)際探區(qū)數(shù)據(jù)對(duì)傳統(tǒng)CNN方法和本文方法進(jìn)行了測(cè)試.首先選用中國(guó)西部某陸上三維探區(qū)真實(shí)數(shù)據(jù)集,該探區(qū)信噪比較高,偏移距展布范圍為0~8 km,數(shù)據(jù)的時(shí)間總長(zhǎng)度7 s,采樣間隔為2 ms.從中等間隔抽取若干條測(cè)線(xiàn),并以相同間隔從每條測(cè)線(xiàn)上抽取CMP道集,共獲得包含1000張速度譜圖像的數(shù)據(jù)集.將生成的數(shù)據(jù)集按 8∶1∶1分為訓(xùn)練集、驗(yàn)證集和測(cè)試集三部分,用于訓(xùn)練傳統(tǒng) CNN 方法和本文方法設(shè)計(jì)的速度譜識(shí)別模型并對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行測(cè)試.訓(xùn)練過(guò)程中梯度下降損失和推理過(guò)程中mAP精度表現(xiàn)如圖9所示.從圖9中可以看出本文方法計(jì)算梯度下降得到的Loss曲線(xiàn)更為收斂,且應(yīng)用本文方法訓(xùn)練的目標(biāo)檢測(cè)器在Marmousi數(shù)據(jù)集上的mAP精度最高為75.7,而傳統(tǒng)CNN方法在Marmousi數(shù)據(jù)集上的mAP精度最高僅為57.4.

    圖10展示了采用兩種方法得到的2560測(cè)線(xiàn)上CMP2240道集的速度拾取結(jié)果及動(dòng)校正前后的道集.為了增加對(duì)比效果,將人工拾取的速度曲線(xiàn)與網(wǎng)絡(luò)模型自動(dòng)拾取的速度曲線(xiàn)可視化到同一個(gè)速度譜中.圖10(a,b)分別展示了兩種方法的拾取結(jié)果,從速度譜中可以看出在3.5 s后雖然受到多次波的干擾,且在5 s以后能量團(tuán)橫向較為分散,但FCOS方法自動(dòng)拾取的疊加速度與人工拾取的疊加速度經(jīng)插值處理后的曲線(xiàn)基本重合.而傳統(tǒng)的CNN方法得到的速度曲線(xiàn)在3.5~4.5 s內(nèi)波動(dòng)較大,拾取結(jié)果明顯偏離人工拾取速度值.圖10c展示了動(dòng)校正前的CMP2240道集圖像,從中可以看出該道集淺層反射波時(shí)距曲線(xiàn)成像較為清晰,在深層區(qū)間(5 s以后)未見(jiàn)明顯的反射雙曲線(xiàn).分別對(duì)兩種方法的拾取結(jié)果進(jìn)行動(dòng)校正,采用FCOS網(wǎng)絡(luò)模型自動(dòng)拾取的疊加速度的動(dòng)校正的結(jié)果如圖10e所示,從圖中可知FCOS方法拾取的疊加速度動(dòng)校正結(jié)果的主要反射同相軸拉平程度與人工拾取速度動(dòng)校正結(jié)果(圖10f)相一致,這進(jìn)一步展現(xiàn)了本文方法的穩(wěn)健性.圖10d展示了傳統(tǒng)CNN方法得到的疊加速度的動(dòng)校正結(jié)果,在3.7 s處該方法受低速值的影響錯(cuò)誤的拉平了道集中的時(shí)距曲線(xiàn),而在4.9 s處該方法又受低速值的影響,導(dǎo)致同相軸隨偏移距增大出現(xiàn)上翹現(xiàn)象.

    對(duì)第2560測(cè)線(xiàn)上的CMP道集進(jìn)行測(cè)試,然后將應(yīng)用兩種方法得到的疊加速度場(chǎng)與人工拾取速度場(chǎng)進(jìn)行對(duì)比.總體上看基于FCOS方法拾取的速度場(chǎng)(圖11c)與人工拾取速度場(chǎng)(圖11a)在構(gòu)造上基本一致,且通過(guò)本文方法得到的速度場(chǎng)在橫向上速度走勢(shì)更加平穩(wěn);而傳統(tǒng)CNN方法得到的速度場(chǎng)(圖11b)在4.5 s以后速度值偏小,且在橫向分布上較為散亂,未出現(xiàn)明顯的速度分界面.為了進(jìn)一步展現(xiàn)結(jié)果的準(zhǔn)確性,本文展示了采用FCOS方法所得到的第2560測(cè)線(xiàn)上局部疊加剖面(圖12c)與和人工拾取疊加剖面(圖12a)的對(duì)比,從圖中可以看出FCOS方法得到的疊加剖面與人工拾取結(jié)果非常近似;而傳統(tǒng)CNN方法得到的局部疊加剖面(圖12b)與本文方法相比,其在3.5 s后由于受多次波能量團(tuán)的影響導(dǎo)致拾取的速度值偏低,在地震剖面上出現(xiàn)多次波能量聚焦效應(yīng),同相軸更為清晰(如紅色框所示).

    然后,采用中國(guó)西部某三維探區(qū)的信噪比較低的CMP道集數(shù)據(jù)進(jìn)一步測(cè)試本文方法.該探區(qū)偏移距展布范圍為0~7 km,數(shù)據(jù)的時(shí)間總長(zhǎng)度8 s,采樣間隔為2 ms.從中等間隔抽取若干條測(cè)線(xiàn),并從每條測(cè)線(xiàn)上等間距抽取CMP道集,共獲得包含1000張速度譜圖像的數(shù)據(jù)集,進(jìn)而采用在上述實(shí)際探區(qū)中所訓(xùn)練模型來(lái)測(cè)試本文方法的泛化性能.由于該探區(qū)信噪比低,在疊前道集中未見(jiàn)明顯的反射波時(shí)距曲線(xiàn),本文只展示出傳統(tǒng)CNN方法與本文方法得到的預(yù)測(cè)結(jié)果與人工拾取結(jié)果的對(duì)比(圖13).由圖中可以看出該探區(qū)數(shù)據(jù)質(zhì)量不高,速度譜上的能量團(tuán)聚焦特性較差.傳統(tǒng)CNN方法僅拾取到了0~3 s區(qū)間內(nèi)的疊加速度(圖13a),且在1.1 s處和2.0 s處拾取到部分高速異常值,在3~8 s區(qū)間范圍內(nèi)由于速度譜能量團(tuán)較為分散,傳統(tǒng)CNN方法未能自動(dòng)拾取.圖13b展示了僅使用FCOS方法(圖中黃線(xiàn))的拾取結(jié)果,該方法僅能拾取到0~4 s區(qū)間內(nèi)的疊加速度,由于其不滿(mǎn)足時(shí)間域全局的拾取要求,本文在FCOS網(wǎng)絡(luò)的基礎(chǔ)上加入一個(gè)基于DNN的線(xiàn)性回歸模型,即通過(guò)學(xué)習(xí)0~4 s區(qū)間內(nèi)的疊加速度的線(xiàn)性趨勢(shì)來(lái)完成對(duì)4~8 s區(qū)間內(nèi)的疊加速度的預(yù)測(cè),最終的預(yù)測(cè)結(jié)果如圖13b(綠線(xiàn)+黃線(xiàn))所示.由圖可知在加入線(xiàn)性歸回模型后本文算法的輸出結(jié)果與人工拾取結(jié)果基本相符.

    圖14展示了人工拾取速度場(chǎng)(圖14a)、傳統(tǒng)CNN方法拾取的速度場(chǎng)(圖14b)和本文方法拾取的速度場(chǎng)(圖14c)的比較.受信噪比較低的影響,傳統(tǒng)CNN方法只能拾取到0~3 s區(qū)間的速度場(chǎng),且在部分道存在低速或高速速度異常導(dǎo)致橫向連續(xù)性較差,在3.0 s以后的速度相較人工拾取速度場(chǎng)中的速度值普遍偏低.而本文方法拾取的速度場(chǎng)在速度變化趨勢(shì)上與人工拾取速度場(chǎng)相一致,但由于應(yīng)用本文方法拾取速度場(chǎng)并未充分考慮道間相似性,拾取結(jié)果在橫向連續(xù)性上相較人工拾取速度場(chǎng)較差,在速度分界面處的速度值較為發(fā)散.對(duì)上述三種方法拾取的道集進(jìn)行動(dòng)校正和疊加處理后得到疊加剖面.總體來(lái)看,本文方法的疊加剖面(圖15c)與人工拾取結(jié)果(圖15a)相近;在紅色框所示位置使用本文方法得到的同相軸相較傳統(tǒng)CNN方法(圖15b)表現(xiàn)更為連續(xù),這進(jìn)一步檢驗(yàn)了本文拾取結(jié)果的準(zhǔn)確性.

    圖11 實(shí)際探區(qū)速度場(chǎng)(a) 人工方法拾取的速度場(chǎng);(b) 傳統(tǒng)CNN方法自動(dòng)拾取的速度場(chǎng);(c) 本文方法自動(dòng)拾取的速度場(chǎng).Fig.11 The actual velocity field in the exploration area(a) The velocity field picked up by manual method;(b) Velocity field automatically picked up by traditional CNN method;(c) The velocity field automatically picked by the method in this paper.

    圖12 實(shí)際地震資料疊加剖面(a) 人工方法拾取的速度場(chǎng)疊加;(b) 傳統(tǒng)CNN方法自動(dòng)拾取的速度場(chǎng)疊加;(c) 本文方法自動(dòng)拾取的速度場(chǎng)疊加.Fig.12 Stacked sections ofactual seismic data(a) The Stacked section of velocity field picked up by manual method;(b) The Stacked section of velocity field picked up by traditional CNN method;(c) The Stacked section of velocity field picked up by the method in this paper.

    圖13 實(shí)際數(shù)據(jù)速度曲線(xiàn) (a) 傳統(tǒng)CNN方法拾取的速度曲線(xiàn);(b) 本文方法拾取的速度曲線(xiàn).Fig.13 Velocity curve of actual data(a) Velocity curve picked up by traditional CNN method;(b) Velocity curves picked up by the method in this paper.

    圖14 拾取的速度場(chǎng)與真實(shí)速度場(chǎng)比較(a) 人工方法拾取的速度場(chǎng);(b) 傳統(tǒng)CNN方法自動(dòng)拾取的速度場(chǎng);(c) 本文方法自動(dòng)拾取的速度場(chǎng).Fig.14 The actual velocity field in the exploration area(a) The velocity field picked up by manual method;(b) Velocity field automatically picked up by traditional CNN method;(c) The velocity field automatically picked by the method in this paper.

    圖15 實(shí)際地震資料疊加剖面(a) 人工方法拾取的速度場(chǎng)疊加;(b) 傳統(tǒng)CNN方法自動(dòng)拾取的速度場(chǎng)疊加;(c) 本文方法自動(dòng)拾取的速度場(chǎng)疊加.Fig.15 Stacked sections of actual seismic data(a) The Stacked section of velocity field picked up by manual method;(b) The Stacked section of velocity field picked up by traditional CNN method;(c) The Stacked section of velocity field picked up by the method in this paper.

    3 結(jié)論

    本文對(duì)FCOS目標(biāo)檢測(cè)模型進(jìn)行了優(yōu)化,分別采用模型數(shù)據(jù)集和實(shí)際工區(qū)數(shù)據(jù)集對(duì)該網(wǎng)絡(luò)模型進(jìn)行測(cè)試,并與傳統(tǒng)CNN方法與人工拾取方法得到的結(jié)果做對(duì)比來(lái)展示本文算法的應(yīng)用效果.通過(guò)對(duì)實(shí)驗(yàn)結(jié)果做分析得到以下結(jié)論:

    (1) 優(yōu)化后的FCOS模型能夠較為精準(zhǔn)的實(shí)現(xiàn)速度譜中能量團(tuán)的定位與識(shí)別,從而使速度拾取結(jié)果具有更高的準(zhǔn)確率.針對(duì)低信噪比探區(qū)中存在能量團(tuán)聚焦特征較差的情況又加入線(xiàn)性回歸模型來(lái)完成對(duì)復(fù)雜地質(zhì)體的速度拾取.

    (2) 本文設(shè)計(jì)的網(wǎng)絡(luò)模型能夠較為快速的拾取疊加速度,預(yù)測(cè)10000張速度譜圖片僅需耗時(shí)30 min.預(yù)測(cè)過(guò)程不需要任何人工干預(yù),其自動(dòng)化程度及準(zhǔn)確度較高,為解決地震資料速度譜自動(dòng)拾取問(wèn)題提供了一種新的方案.

    由于本文在進(jìn)行研究時(shí)所用數(shù)據(jù)集數(shù)據(jù)量比較少,一般情況下目標(biāo)檢測(cè)數(shù)據(jù)集圖片的數(shù)據(jù)量都在幾萬(wàn)張以上.在實(shí)際應(yīng)用中,可以通過(guò)增加數(shù)據(jù)量重新訓(xùn)練網(wǎng)絡(luò)模型,使其應(yīng)用于疊加速度拾取時(shí)具有更高的精度.

    猜你喜歡
    人工能量自動(dòng)
    人工3D脊髓能幫助癱瘓者重新行走?
    軍事文摘(2022年8期)2022-11-03 14:22:01
    人工,天然,合成
    人工“美顏”
    自動(dòng)捕盜機(jī)
    能量之源
    基于STM32的自動(dòng)喂養(yǎng)機(jī)控制系統(tǒng)
    詩(shī)無(wú)邪傳遞正能量
    新型多孔鉭人工種植牙
    關(guān)于自動(dòng)駕駛
    Stefan Greiner:我們?yōu)槭裁葱枰詣?dòng)駕駛?
    五月天丁香电影| 国语对白做爰xxxⅹ性视频网站| 我要看黄色一级片免费的| 一级av片app| 国产淫语在线视频| .国产精品久久| 寂寞人妻少妇视频99o| 一区二区三区精品91| 亚洲av.av天堂| www.色视频.com| 美女主播在线视频| 亚洲精品日韩在线中文字幕| 亚洲精品亚洲一区二区| 国产一区二区在线观看日韩| 日韩中文字幕视频在线看片 | 国产欧美日韩一区二区三区在线 | 国产成人aa在线观看| 久久久久久久亚洲中文字幕| 熟女电影av网| 美女国产视频在线观看| 国产精品久久久久成人av| 女的被弄到高潮叫床怎么办| 女人久久www免费人成看片| 五月伊人婷婷丁香| 91在线精品国自产拍蜜月| 国产69精品久久久久777片| 欧美变态另类bdsm刘玥| 免费人妻精品一区二区三区视频| 久久久久人妻精品一区果冻| 久久久久性生活片| 麻豆成人av视频| 亚洲精品日韩在线中文字幕| 视频中文字幕在线观看| 狠狠精品人妻久久久久久综合| 久久久久久久久久成人| 亚洲精品久久久久久婷婷小说| 麻豆国产97在线/欧美| 女性被躁到高潮视频| 国产精品麻豆人妻色哟哟久久| 成人亚洲欧美一区二区av| 看非洲黑人一级黄片| 国产黄片视频在线免费观看| 香蕉精品网在线| 在线精品无人区一区二区三 | 久久精品国产鲁丝片午夜精品| 性高湖久久久久久久久免费观看| 色婷婷av一区二区三区视频| 一本一本综合久久| 国产精品麻豆人妻色哟哟久久| 国产av码专区亚洲av| 亚洲国产精品专区欧美| 丰满少妇做爰视频| 久久久久久久亚洲中文字幕| 亚洲激情五月婷婷啪啪| 狂野欧美激情性xxxx在线观看| 97精品久久久久久久久久精品| h视频一区二区三区| 三级国产精品欧美在线观看| 少妇的逼好多水| 国产精品嫩草影院av在线观看| 久久久a久久爽久久v久久| 欧美日韩一区二区视频在线观看视频在线| 九色成人免费人妻av| 在线 av 中文字幕| 日韩在线高清观看一区二区三区| 哪个播放器可以免费观看大片| 校园人妻丝袜中文字幕| 久久精品国产亚洲av天美| 亚洲色图av天堂| 日韩不卡一区二区三区视频在线| 美女cb高潮喷水在线观看| 乱码一卡2卡4卡精品| 精品国产露脸久久av麻豆| 精华霜和精华液先用哪个| 日韩三级伦理在线观看| 色视频在线一区二区三区| 久久 成人 亚洲| 国产综合精华液| 国产精品一区二区在线不卡| 丰满乱子伦码专区| 人人妻人人添人人爽欧美一区卜 | 涩涩av久久男人的天堂| 国产成人午夜福利电影在线观看| 波野结衣二区三区在线| 中文字幕av成人在线电影| 精品久久久久久电影网| 伦理电影免费视频| 青春草亚洲视频在线观看| 国产精品久久久久久av不卡| 99久国产av精品国产电影| 三级国产精品片| 老女人水多毛片| 91精品一卡2卡3卡4卡| 人人妻人人澡人人爽人人夜夜| 日韩av在线免费看完整版不卡| 亚洲欧美成人综合另类久久久| 欧美少妇被猛烈插入视频| av专区在线播放| 九草在线视频观看| 91午夜精品亚洲一区二区三区| 亚洲伊人久久精品综合| 亚洲国产精品一区三区| 国产老妇伦熟女老妇高清| 国产爽快片一区二区三区| 另类亚洲欧美激情| 免费观看在线日韩| 久久精品国产自在天天线| 多毛熟女@视频| 啦啦啦中文免费视频观看日本| 黄色视频在线播放观看不卡| 老熟女久久久| 熟女av电影| 永久网站在线| 久久国产乱子免费精品| 一本色道久久久久久精品综合| 免费久久久久久久精品成人欧美视频 | 久久久久久久久久成人| 精品酒店卫生间| 婷婷色麻豆天堂久久| 欧美日韩在线观看h| 亚洲成人手机| 最近的中文字幕免费完整| 国产一区有黄有色的免费视频| 久久国产精品大桥未久av | 国产精品偷伦视频观看了| 久久久国产一区二区| 好男人视频免费观看在线| 嘟嘟电影网在线观看| 欧美人与善性xxx| 日本黄大片高清| 在线观看三级黄色| 亚洲精品国产av成人精品| 亚洲精品中文字幕在线视频 | 乱系列少妇在线播放| av线在线观看网站| 尾随美女入室| 成人毛片a级毛片在线播放| 国内揄拍国产精品人妻在线| 欧美一区二区亚洲| av女优亚洲男人天堂| av线在线观看网站| 精品少妇黑人巨大在线播放| 在线观看av片永久免费下载| 亚州av有码| 男男h啪啪无遮挡| 精品亚洲乱码少妇综合久久| 2022亚洲国产成人精品| 97热精品久久久久久| 黄色欧美视频在线观看| 精品国产乱码久久久久久小说| 18禁裸乳无遮挡免费网站照片| 免费人成在线观看视频色| 日韩人妻高清精品专区| 免费不卡的大黄色大毛片视频在线观看| 国产精品人妻久久久影院| 黄片无遮挡物在线观看| 国产一区亚洲一区在线观看| 欧美zozozo另类| 国产亚洲一区二区精品| kizo精华| 久久久亚洲精品成人影院| 黑人高潮一二区| 一区二区三区乱码不卡18| 老师上课跳d突然被开到最大视频| 国产男女超爽视频在线观看| 男人舔奶头视频| 国产一区亚洲一区在线观看| 91精品伊人久久大香线蕉| 久久精品国产亚洲网站| 在线观看免费视频网站a站| 美女脱内裤让男人舔精品视频| 日韩中字成人| 欧美日韩视频高清一区二区三区二| 亚洲人成网站在线播| 在线看a的网站| 一级毛片久久久久久久久女| 精品一区二区三卡| 草草在线视频免费看| 国产精品熟女久久久久浪| 日韩在线高清观看一区二区三区| 丰满人妻一区二区三区视频av| 日韩 亚洲 欧美在线| 亚洲国产日韩一区二区| 多毛熟女@视频| 女的被弄到高潮叫床怎么办| 日韩av在线免费看完整版不卡| 91久久精品国产一区二区成人| 成人综合一区亚洲| 亚洲美女黄色视频免费看| videossex国产| 人人妻人人看人人澡| 亚洲美女视频黄频| 三级国产精品片| 黑丝袜美女国产一区| 精品久久久久久久久av| 日日啪夜夜爽| 精品人妻视频免费看| 久久国内精品自在自线图片| 色综合色国产| 久久ye,这里只有精品| 99热这里只有精品一区| 国内少妇人妻偷人精品xxx网站| 日本av手机在线免费观看| 亚洲第一av免费看| 一区二区av电影网| 天堂俺去俺来也www色官网| 老司机影院毛片| 国国产精品蜜臀av免费| 亚洲av在线观看美女高潮| 男女边摸边吃奶| 国产成人a∨麻豆精品| 国产欧美日韩精品一区二区| 国产一区二区在线观看日韩| 建设人人有责人人尽责人人享有的 | 高清午夜精品一区二区三区| 久久6这里有精品| 成人美女网站在线观看视频| 蜜臀久久99精品久久宅男| videossex国产| kizo精华| 一级a做视频免费观看| 国产欧美日韩一区二区三区在线 | 亚洲国产色片| 大香蕉久久网| 啦啦啦在线观看免费高清www| 干丝袜人妻中文字幕| 亚洲国产高清在线一区二区三| 免费大片18禁| 亚洲av中文字字幕乱码综合| 日韩大片免费观看网站| 欧美xxxx性猛交bbbb| 日韩国内少妇激情av| 成人国产麻豆网| 国产一区二区三区综合在线观看 | 肉色欧美久久久久久久蜜桃| 久久人人爽人人爽人人片va| 亚洲欧美中文字幕日韩二区| h视频一区二区三区| 黑丝袜美女国产一区| 久久99热这里只有精品18| 国产高清有码在线观看视频| 人妻系列 视频| av视频免费观看在线观看| 午夜福利在线在线| 深爱激情五月婷婷| 人妻系列 视频| 久久久久久久久久成人| 观看av在线不卡| 啦啦啦视频在线资源免费观看| 18禁动态无遮挡网站| 亚洲色图av天堂| 七月丁香在线播放| 在线 av 中文字幕| 欧美日韩一区二区视频在线观看视频在线| 毛片一级片免费看久久久久| 亚洲,欧美,日韩| 高清日韩中文字幕在线| 欧美高清性xxxxhd video| 自拍偷自拍亚洲精品老妇| 色婷婷久久久亚洲欧美| 欧美区成人在线视频| 在线观看美女被高潮喷水网站| 夜夜爽夜夜爽视频| 高清日韩中文字幕在线| 在线亚洲精品国产二区图片欧美 | 身体一侧抽搐| 亚洲欧美清纯卡通| 黄色日韩在线| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验| 多毛熟女@视频| videos熟女内射| av免费观看日本| 少妇 在线观看| 午夜福利在线观看免费完整高清在| 啦啦啦视频在线资源免费观看| 国产又色又爽无遮挡免| 国产成人精品婷婷| 97在线人人人人妻| 日本wwww免费看| 婷婷色av中文字幕| 午夜福利高清视频| 热99国产精品久久久久久7| 亚洲成色77777| 欧美一区二区亚洲| 国产精品人妻久久久久久| 亚洲国产欧美在线一区| 国产色爽女视频免费观看| 日韩精品有码人妻一区| 国产v大片淫在线免费观看| 亚洲精品自拍成人| 最近的中文字幕免费完整| 亚洲欧美精品自产自拍| 中文字幕精品免费在线观看视频 | 国产片特级美女逼逼视频| .国产精品久久| 最近中文字幕高清免费大全6| 午夜视频国产福利| 国产亚洲午夜精品一区二区久久| 免费少妇av软件| 久久99热这里只有精品18| 晚上一个人看的免费电影| 亚洲精品日韩在线中文字幕| 久久久久国产精品人妻一区二区| 国产乱人视频| 国产人妻一区二区三区在| 最近中文字幕2019免费版| 国产av国产精品国产| 久久久精品免费免费高清| 亚洲精品色激情综合| 国产深夜福利视频在线观看| 美女主播在线视频| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 干丝袜人妻中文字幕| 成人毛片60女人毛片免费| 免费大片18禁| 在线亚洲精品国产二区图片欧美 | 丝袜脚勾引网站| 亚洲在久久综合| 各种免费的搞黄视频| 韩国av在线不卡| 欧美区成人在线视频| 大片电影免费在线观看免费| 午夜福利影视在线免费观看| 亚洲国产毛片av蜜桃av| 久久精品人妻少妇| 久久久久久久精品精品| 涩涩av久久男人的天堂| 成年美女黄网站色视频大全免费 | 欧美最新免费一区二区三区| 九色成人免费人妻av| 亚洲精品456在线播放app| 久久国产精品大桥未久av | 黄色怎么调成土黄色| 亚洲图色成人| 男人舔奶头视频| 亚洲欧美中文字幕日韩二区| 久久热精品热| 秋霞在线观看毛片| 日本与韩国留学比较| 韩国av在线不卡| 男女边摸边吃奶| 国产 精品1| av在线蜜桃| 日本一二三区视频观看| 永久网站在线| 国产淫片久久久久久久久| 制服丝袜香蕉在线| 男人爽女人下面视频在线观看| 午夜福利在线观看免费完整高清在| 精品久久久精品久久久| 全区人妻精品视频| 一区二区三区乱码不卡18| 黑丝袜美女国产一区| 王馨瑶露胸无遮挡在线观看| 久久女婷五月综合色啪小说| 一区二区三区免费毛片| 夜夜看夜夜爽夜夜摸| 最近手机中文字幕大全| 身体一侧抽搐| 日韩 亚洲 欧美在线| 精品人妻偷拍中文字幕| 久久精品国产亚洲网站| 国产真实伦视频高清在线观看| 性色av一级| 黄片无遮挡物在线观看| 中文字幕亚洲精品专区| 一区二区av电影网| 精品亚洲成国产av| 国产成人精品婷婷| 亚洲人成网站在线观看播放| 18禁在线播放成人免费| 一个人看视频在线观看www免费| 人人妻人人添人人爽欧美一区卜 | 丰满乱子伦码专区| 亚洲无线观看免费| 久久影院123| 国产午夜精品一二区理论片| 美女xxoo啪啪120秒动态图| 人妻一区二区av| 久久综合国产亚洲精品| 国产在线免费精品| 国产亚洲午夜精品一区二区久久| 少妇人妻精品综合一区二区| 久久久色成人| 亚洲无线观看免费| 久久精品国产自在天天线| 成年人午夜在线观看视频| 内射极品少妇av片p| 亚洲精品色激情综合| 香蕉精品网在线| 人妻夜夜爽99麻豆av| 亚洲aⅴ乱码一区二区在线播放| 国产成人a区在线观看| 自拍偷自拍亚洲精品老妇| 亚洲精品国产av成人精品| 国产精品99久久久久久久久| 国产精品国产三级国产av玫瑰| 热re99久久精品国产66热6| 精品国产露脸久久av麻豆| 一级二级三级毛片免费看| 国产成人精品一,二区| 亚洲伊人久久精品综合| 亚洲av成人精品一区久久| 中文精品一卡2卡3卡4更新| 国产成人a区在线观看| 亚洲内射少妇av| 久久毛片免费看一区二区三区| 我要看日韩黄色一级片| 日韩国内少妇激情av| 国产免费视频播放在线视频| 一级毛片黄色毛片免费观看视频| 亚洲国产色片| 自拍偷自拍亚洲精品老妇| 91精品国产九色| 日韩不卡一区二区三区视频在线| 亚洲国产毛片av蜜桃av| 黑人猛操日本美女一级片| 91精品一卡2卡3卡4卡| 色视频在线一区二区三区| 国产成人精品福利久久| 国产在视频线精品| 男的添女的下面高潮视频| 99热全是精品| 久久精品久久久久久噜噜老黄| 尤物成人国产欧美一区二区三区| 大片电影免费在线观看免费| 久久亚洲国产成人精品v| 久久青草综合色| 亚洲经典国产精华液单| 亚洲精品日本国产第一区| 国产乱人视频| 欧美成人精品欧美一级黄| 亚洲欧洲国产日韩| 少妇精品久久久久久久| 国产精品三级大全| 一个人看视频在线观看www免费| 在线观看免费高清a一片| 欧美区成人在线视频| 国产片特级美女逼逼视频| 久久久午夜欧美精品| 久热久热在线精品观看| 一级毛片 在线播放| 免费人成在线观看视频色| 亚洲va在线va天堂va国产| 香蕉精品网在线| 在线亚洲精品国产二区图片欧美 | 国产成人免费无遮挡视频| 亚洲av中文字字幕乱码综合| 中文精品一卡2卡3卡4更新| 精品久久久噜噜| 国产一区亚洲一区在线观看| 丝袜喷水一区| 99久久中文字幕三级久久日本| 干丝袜人妻中文字幕| 欧美精品一区二区免费开放| 妹子高潮喷水视频| 91午夜精品亚洲一区二区三区| 黄色视频在线播放观看不卡| 三级国产精品欧美在线观看| 日韩免费高清中文字幕av| 伦理电影大哥的女人| 久久久久久久精品精品| 亚洲无线观看免费| 18禁在线播放成人免费| 夜夜爽夜夜爽视频| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 亚洲精品日韩在线中文字幕| 男人爽女人下面视频在线观看| 久久人妻熟女aⅴ| 亚洲精品国产av成人精品| 国产男女超爽视频在线观看| 久久韩国三级中文字幕| 内射极品少妇av片p| 一级毛片久久久久久久久女| 蜜桃在线观看..| 精品一区在线观看国产| 国产淫语在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产日韩欧美在线精品| 最近最新中文字幕免费大全7| 久久久久久久久久久丰满| 亚洲欧美成人综合另类久久久| 香蕉精品网在线| 亚洲美女黄色视频免费看| 亚洲图色成人| 97在线视频观看| 十八禁网站网址无遮挡 | 亚洲人成网站高清观看| 99九九线精品视频在线观看视频| 内射极品少妇av片p| 日本av手机在线免费观看| 亚洲精品乱码久久久v下载方式| 成人午夜精彩视频在线观看| 各种免费的搞黄视频| 如何舔出高潮| 国产一区二区三区av在线| 干丝袜人妻中文字幕| 少妇被粗大猛烈的视频| 麻豆乱淫一区二区| 国产高潮美女av| 亚洲四区av| 九九在线视频观看精品| 国产精品99久久99久久久不卡 | 色网站视频免费| 男人和女人高潮做爰伦理| 亚洲av中文av极速乱| 欧美日韩亚洲高清精品| 免费不卡的大黄色大毛片视频在线观看| 国产午夜精品一二区理论片| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区性色av| 欧美3d第一页| 97超视频在线观看视频| 中文资源天堂在线| 能在线免费看毛片的网站| 日韩一本色道免费dvd| 中文字幕亚洲精品专区| 久久人人爽人人爽人人片va| 一个人看视频在线观看www免费| 性色avwww在线观看| 婷婷色麻豆天堂久久| 亚洲美女搞黄在线观看| 夫妻午夜视频| 欧美激情极品国产一区二区三区 | 一区在线观看完整版| 国产 一区 欧美 日韩| 在线观看国产h片| av福利片在线观看| 最近2019中文字幕mv第一页| 日韩不卡一区二区三区视频在线| 精品一区二区三卡| av国产精品久久久久影院| 久久人妻熟女aⅴ| 美女高潮的动态| 日韩大片免费观看网站| 色5月婷婷丁香| 免费观看的影片在线观看| 在线观看免费日韩欧美大片 | 国产久久久一区二区三区| 97在线视频观看| 久久女婷五月综合色啪小说| 免费在线观看成人毛片| 美女国产视频在线观看| 人妻一区二区av| 毛片女人毛片| 五月伊人婷婷丁香| 联通29元200g的流量卡| 久久韩国三级中文字幕| 在线播放无遮挡| 最近手机中文字幕大全| 久久国产乱子免费精品| 国产淫语在线视频| 天天躁日日操中文字幕| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看| 日韩欧美 国产精品| 精品国产乱码久久久久久小说| 18禁裸乳无遮挡免费网站照片| 美女脱内裤让男人舔精品视频| 国产在视频线精品| 两个人的视频大全免费| 成年av动漫网址| 亚洲精品国产成人久久av| 美女脱内裤让男人舔精品视频| 欧美丝袜亚洲另类| 亚州av有码| 国产欧美亚洲国产| 欧美xxxx性猛交bbbb| 精品一区二区免费观看| 人妻一区二区av| 国产老妇伦熟女老妇高清| 国内揄拍国产精品人妻在线| 蜜桃在线观看..| 久久国内精品自在自线图片| 国产一区二区在线观看日韩| 国产午夜精品一二区理论片| 亚洲欧美清纯卡通| 精华霜和精华液先用哪个| 久热这里只有精品99| 最后的刺客免费高清国语| 国产 一区精品| 国产淫语在线视频| 王馨瑶露胸无遮挡在线观看| 欧美精品一区二区免费开放| 青春草视频在线免费观看| 亚洲成人手机| 亚洲精品亚洲一区二区| 一个人免费看片子| 欧美zozozo另类| 日日啪夜夜撸| 国产成人a区在线观看| 国产免费视频播放在线视频| 免费高清在线观看视频在线观看| 乱码一卡2卡4卡精品| 观看美女的网站| 天天躁夜夜躁狠狠久久av| 97在线人人人人妻| 亚洲精品久久午夜乱码| 亚洲激情五月婷婷啪啪| 少妇精品久久久久久久| 少妇猛男粗大的猛烈进出视频| 欧美亚洲 丝袜 人妻 在线| 日韩视频在线欧美| 王馨瑶露胸无遮挡在线观看| 少妇人妻一区二区三区视频| 国产av一区二区精品久久 | 免费人成在线观看视频色| 成人无遮挡网站| 日本av免费视频播放| 人妻 亚洲 视频| 精品久久久精品久久久| 3wmmmm亚洲av在线观看|