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

    基于神經(jīng)網(wǎng)絡(luò)與早至波的波方程反演

    2024-11-05 00:00:00李亞杰余瀚湯致文孫明皓
    科技創(chuàng)新與應(yīng)用 2024年30期

    摘 要:針對傳統(tǒng)波方程全波形反演步驟繁瑣、計(jì)算量大和難度高等問題,提出一種基于全連接神經(jīng)網(wǎng)絡(luò)的早至波反演方法。實(shí)驗(yàn)結(jié)果表明,所提方法在早至波反演中交并比IoU(Intersection over Union)和平均精度均值mAP(mean Average Precision)分別達(dá)到74.01%和73.72%,可對速度模型進(jìn)行有效重構(gòu)。

    關(guān)鍵詞:波方程;早至波;神經(jīng)網(wǎng)絡(luò);反演;交并比;平均精度均值

    中圖分類號:P631.4 文獻(xiàn)標(biāo)志碼:A 文章編號:2095-2945(2024)30-0040-04

    Abstract: Aiming at the problems of tedious steps, large amount of calculation and high difficulty in the full waveform inversion of traditional wave equation, a method for early-arrival wave inversion based on fully connected neural network is proposed. Experimental results show that the proposed method achieves 74.01% and 73.72% in the Intersection over Union(IoU) and mean Average Precision(mAP) of early-arrival wave inversion, which can effectively reconstruct the velocity model.

    Keywords: wave equation; early arrival wave; neural network; inversion; Intersection over Urion; mean Average Precision

    波方程主要描述自然界中的各種波動(dòng)現(xiàn)象,求解一維波動(dòng)方程的一般步驟是將其轉(zhuǎn)化為一個(gè)簡單的常微分方程或特殊的偏微分方程,然后通過求解該方程得到波動(dòng)解析表達(dá)式。我們通過有限差分法求解以下波動(dòng)方程[1]進(jìn)行數(shù)值正演得到仿真波形數(shù)據(jù)集

    ▽P=-ρ

    P=-ρc2▽u+Src(z,t) , (1)

    式中:z代表豎直方向的一維空間變量,t為時(shí)間項(xiàng),P=P(z,t)為壓強(qiáng),作為常規(guī)空間梯度算子,u(z,t)代表豎直方向上的位移,ρ=ρ(z)為壓強(qiáng)場密度分布,c=c(z)代表速度場,Src代表激發(fā)振動(dòng)源項(xiàng)。

    相較于正演,傳統(tǒng)波方程反演的具體公式較為復(fù)雜,其本質(zhì)框架是對以下殘差平方和目標(biāo)函數(shù)J的優(yōu)化[2]

    J=||dcalc-dobs||2 , (2)

    式中:dcalc代表通過反演計(jì)算得到的速度場數(shù)據(jù),dobs代表原始速度場數(shù)據(jù)。搭建神經(jīng)網(wǎng)絡(luò)實(shí)現(xiàn)反演,意義在于使用神經(jīng)網(wǎng)絡(luò)去替代反演的過程,減少計(jì)算難度和時(shí)間。

    1 方法

    訓(xùn)練神經(jīng)網(wǎng)絡(luò)反演波方程一般根據(jù)早至波或初至走時(shí),張利振等[3]提出的根據(jù)初至走時(shí)反演地層速度求解非線性地震走時(shí)問題以及Tarantola[4]提出的基于最小二乘法的全波形反演理論對波形反演起到了很大的推進(jìn)作用,通過早至波進(jìn)行反演也已經(jīng)擁有了一定研究基礎(chǔ)[5]。甚至也有相關(guān)研究[6]將早至波和初至走時(shí)聯(lián)合反演。本文將利用早至波信息,完成基于早至波的波形反演。

    近年來,深度學(xué)習(xí)在圖像識(shí)別、語音處理等領(lǐng)域得到成功應(yīng)用,為地震速度建模提供了新的方法[7-11],例如U-Net神經(jīng)網(wǎng)絡(luò)、生成對抗神經(jīng)網(wǎng)絡(luò)和全卷積神經(jīng)網(wǎng)絡(luò)等都對該領(lǐng)域的研究提供便利。然而,當(dāng)前可獲取地震數(shù)據(jù)集較少,構(gòu)建豐富且具有復(fù)雜特征的地震速度模型是當(dāng)前深度學(xué)習(xí)速度建模的一大挑戰(zhàn)。本文使用已搭建好的地震波仿真器依據(jù)一定的速度場計(jì)算得出地震早至波的早至壓強(qiáng),避免真實(shí)數(shù)據(jù)中噪聲對神經(jīng)網(wǎng)絡(luò)訓(xùn)練的干擾,使訓(xùn)練更加理想化。

    神經(jīng)網(wǎng)絡(luò)的具體優(yōu)化方式為反向傳播,是訓(xùn)練神經(jīng)網(wǎng)絡(luò)最有效且最常用的方法[12-14]:訓(xùn)練過程中神經(jīng)網(wǎng)絡(luò)的輸出結(jié)果與實(shí)際結(jié)果會(huì)產(chǎn)生誤差,計(jì)算過程輸出與實(shí)際輸出之間的誤差,并將該誤差從輸出層向隱藏層直至輸入層反向傳播[15]

    wl→wl-η∑x δl(al-1)T

    bl→bl-η∑x δl, (3)

    式中:wl表示第l層權(quán)重,al-1表示第l-1層輸出,δl表示第l層誤差,bl表示第l層偏置,η表示學(xué)習(xí)率,“→”表示對內(nèi)部參數(shù)的更新。

    本文計(jì)算不同速度場的波方程得出縱深-時(shí)間的時(shí)空壓強(qiáng)場數(shù)組,對于速度場而言:由長度為2 000 ft的縱深上的速度值構(gòu)成速度場;對于壓強(qiáng)場而言:每一列代表縱軸上某一個(gè)接收器在不同時(shí)間上的壓強(qiáng)大小,共有30列,每一行代表某個(gè)時(shí)間點(diǎn)上不同接收點(diǎn)處的壓強(qiáng)大小,根據(jù)波傳導(dǎo)到最深處的時(shí)間行數(shù)會(huì)有所變化,約在1 000行以內(nèi),且早至拾取完畢的時(shí)刻也會(huì)根據(jù)速度場的變化而變化,假定每0.01 s進(jìn)行一次壓強(qiáng)快照,接收點(diǎn)獲取一次壓強(qiáng)值。

    神經(jīng)網(wǎng)絡(luò)模型以地震早至波的早至壓強(qiáng)作為輸入標(biāo)簽,速度場作為輸出標(biāo)簽,通過兩層隱藏層將地震數(shù)據(jù)由(z,t)域映射到速度模型(z,c)域,并在訓(xùn)練過程中不斷更新隱藏層的權(quán)重w和偏置b,以達(dá)到使反演速度場貼近于真實(shí)速度場的目的。該過程可以用如下公式表示

    θ(w,b)=argmin[θ(d)-v] , (4)

    式中:θ表示神經(jīng)網(wǎng)絡(luò)計(jì)算,w表示神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程中的權(quán)重,b表示神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程中的偏置,d表示地震早至波數(shù)據(jù),θ(d)表示預(yù)測速度場,v表示真實(shí)速度場。

    使用網(wǎng)格搜索的方式進(jìn)行超參數(shù)調(diào)優(yōu),在訓(xùn)練過程中不斷優(yōu)化得到適合該網(wǎng)絡(luò)模型的超參數(shù):學(xué)習(xí)率為0.01,迭代次數(shù)為256次,每批次加載64組數(shù)據(jù)。輸入層神經(jīng)元數(shù)量為30,代表由30個(gè)接收器接收到的早至波壓強(qiáng),兩層隱藏神經(jīng)元數(shù)量分別為300和1 000,輸出層神經(jīng)元數(shù)量為2 000,代表波傳導(dǎo)的深度為2 000 ft。使用FCN神經(jīng)網(wǎng)絡(luò)處理回歸問題最常用的ReLU函數(shù)作為激活函數(shù),用于引入非線性因素,使得神經(jīng)網(wǎng)絡(luò)能夠擬合復(fù)雜的非線性關(guān)系。損失函數(shù)則使用普適的均方差損失函數(shù)。

    具體方法的流程如圖1所示。

    數(shù)據(jù)選?。河?jì)算不同速度場的波方程得出縱深-時(shí)間的時(shí)空數(shù)組,拾取每一列上的早至壓強(qiáng)數(shù)據(jù)。

    輸入層:輸入層神經(jīng)元數(shù)量為30,代表由30個(gè)接收器接收到的早至波壓強(qiáng),將64批次的早至波向量構(gòu)成的混淆矩陣作為神經(jīng)網(wǎng)絡(luò)的輸入。

    隱藏層:設(shè)置兩層隱藏層神經(jīng)元,數(shù)量分別為300和1 000,可以有效且平穩(wěn)地從輸入層30個(gè)神經(jīng)元傳播到輸出層2 000個(gè)神經(jīng)元。在訓(xùn)練過程中不斷計(jì)算梯度和更新神經(jīng)元層之間的偏置與權(quán)重,實(shí)現(xiàn)誤差的反向傳播。

    輸出層:輸出層神經(jīng)元數(shù)量為2 000,代表波傳導(dǎo)的深度為2 000 ft,以64批次的真實(shí)速度場作為輸出。

    反演對比:將神經(jīng)網(wǎng)絡(luò)反演出的速度場與真實(shí)速度場進(jìn)行對比,觀察神經(jīng)網(wǎng)絡(luò)是否能準(zhǔn)確反演波方程。

    2 數(shù)據(jù)來源和選取

    有限差分法有計(jì)算迅速、精度較高、易在計(jì)算機(jī)上實(shí)現(xiàn)的優(yōu)點(diǎn),故使用此方法求解波方程,對速度場進(jìn)行正演,得到其數(shù)值解,使用傳統(tǒng)STA/LTA方法拾取早至波壓強(qiáng),在理想的低信噪比環(huán)境下快速且簡易[16],最終將早至波壓強(qiáng)場與真實(shí)速度場構(gòu)成的數(shù)據(jù)對作為神經(jīng)網(wǎng)絡(luò)的訓(xùn)練標(biāo)簽,劃分80%的訓(xùn)練標(biāo)簽用于訓(xùn)練,20%的訓(xùn)練標(biāo)簽用于測試。

    神經(jīng)網(wǎng)絡(luò)的訓(xùn)練結(jié)果很大程度上依賴數(shù)據(jù)集規(guī)模,不合適的數(shù)據(jù)集往往會(huì)造成以下問題。

    當(dāng)數(shù)據(jù)量過多,模型可能會(huì)過度擬合訓(xùn)練數(shù)據(jù),過度學(xué)習(xí)訓(xùn)練數(shù)據(jù)中的噪聲和特定的樣本特征,而忽略了一般化的模式。這導(dǎo)致模型在未見過的數(shù)據(jù)上表現(xiàn)不佳。

    過少的數(shù)據(jù)集模型可能無法捕捉到數(shù)據(jù)的復(fù)雜結(jié)構(gòu)和泛化模式,導(dǎo)致模型在訓(xùn)練數(shù)據(jù)和測試數(shù)據(jù)上都表現(xiàn)不佳,從而無法在真實(shí)反演中推廣到新數(shù)據(jù)上。

    不合理的數(shù)據(jù)集需要更多的計(jì)算資源進(jìn)行訓(xùn)練,包括更多的內(nèi)存、更長的訓(xùn)練時(shí)間和更高的計(jì)算成本。本文通過多次實(shí)驗(yàn),最終選定數(shù)據(jù)集大小為1 280組,其次限制所有類型速度場的范圍為[1 219.2 m/s,2 743.2 m/s],速度場層數(shù)為4層。本文仿真了上述速度場類型,在此基礎(chǔ)上對速度場求解波方程的到壓強(qiáng)場。圖2展示了其中2種速度場類型對應(yīng)的壓強(qiáng)場的求取過程。

    3 結(jié)果

    3.1 神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果與分析

    如圖3所示,可以看到通過輸入早至波反演出速度場,將真實(shí)速度場與反演速度場進(jìn)行對比,顯示了該方法可以精細(xì)刻畫層狀速度場,整體而言較為接近真實(shí)速度場,反演效果較好,但加入小范圍凹陷后,其并不能有效刻畫凹陷處速度場,KnI/SPfLCXzvEwHxSknstA==對于凹陷處速度場反演并不夠充分,精確度有待提高。

    3.2 模型評估

    使用常見回歸模型評價(jià)指標(biāo),包括交并比IoU(Intersection over Union)、平均精度AP(Average Precision)和平均精度均值mAP(mean Average Precision)。這些指標(biāo)統(tǒng)稱為測試誤差,是測試集的預(yù)測值和真實(shí)值的統(tǒng)計(jì)量通過多方位標(biāo)準(zhǔn)評估該神經(jīng)網(wǎng)絡(luò)的泛化能力和預(yù)測高效性。將速度場與x軸形成的區(qū)域視為待預(yù)測物體,可得出本模型訓(xùn)練集平均IoU約為0.98,測試集平均IoU約為0.74,故將IoU閾值設(shè)定為0.8,IoU≥0.8的測試用例視為正確預(yù)測,將模型轉(zhuǎn)換為分類任務(wù),畫出精準(zhǔn)率-召回率PR(Precision-Recall)曲線、計(jì)算方波狀速度場和階梯狀速度場情況下的AP和總的mAP。

    PR曲線與坐標(biāo)軸圍成的面積是用于評估分類模型優(yōu)劣的常用方法,圖4表明本模型PR曲線較為理想,且表1中2種分類AP值波動(dòng)不大且得出平均值mAP符合預(yù)期,這說明本文神經(jīng)網(wǎng)絡(luò)具有一定泛化能力,在測試范圍內(nèi)預(yù)測精確性較高。

    但是本模型依然存在一定的局限性:在選取數(shù)據(jù)集時(shí),選定速度場為層狀速度場,且規(guī)定速度場的層數(shù)和范圍,這使本神經(jīng)網(wǎng)絡(luò)反演其他類型的速度場時(shí)效果不佳。其次,對于現(xiàn)實(shí)世界中的地震波而言,速度場不可避免會(huì)存在一定塌陷[17],即小范圍孤立異常體,在后續(xù)神經(jīng)網(wǎng)絡(luò)的優(yōu)化中,可針對這一問題進(jìn)行進(jìn)一步研究。

    4 結(jié)論

    本文通過神經(jīng)網(wǎng)絡(luò)研究了面向早至波的波方程反演方法。神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果顯示,早至波神經(jīng)網(wǎng)絡(luò)反演方法可以精確地刻畫層狀速度場,具有一定反演精度。隨機(jī)訓(xùn)練過程對相同層狀速度模型的建模結(jié)果顯示,模型預(yù)測精度不具有隨機(jī)性,是比較穩(wěn)定的,并且花費(fèi)較少的時(shí)間成本和運(yùn)算成本[18]。但在現(xiàn)實(shí)世界的地震波或者其他應(yīng)用場景中,需要獲取足夠數(shù)量且具有不同特征的數(shù)據(jù)集,以提高神經(jīng)網(wǎng)絡(luò)的的泛化能力和應(yīng)用范圍。

    參考文獻(xiàn):

    [1] 斯小琴,陳大偉.一維波動(dòng)方程的計(jì)算模擬[J].淮北師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2022,43(3):22-25.

    [2] 崔叢越,王彥賓.一維速度結(jié)構(gòu)的地震波全波形反演理論測試[J].北京大學(xué)學(xué)報(bào)自然科學(xué)版,2019,55(2):246-252.

    [3] 張利振,孫成禹,王志農(nóng),等.面波信息約束的初至波走時(shí)層析反演方法[J].物探與化探,2023,47(5):1198-1205.

    [4] TARANTOLA A. Inversion seismi reflection data in the acoustic approximation[J].Geophysics,1984,49(8):1259-1266.

    [5] 孫明皓,余瀚,陳雨青,等.基于U形多層感知機(jī)網(wǎng)絡(luò)的地震波初至拾取與反演[J/OL].計(jì)算機(jī)應(yīng)用:1-11[2024-04-15].http://kns.cnki.net/kcms/detail/51.1307.TP.20230828.1444.011.html.

    [6] 許祥,鄒志輝,韓明亮,等.聯(lián)合地震初至走時(shí)與早至波形的深度學(xué)習(xí)速度建模[J].地球物理學(xué)報(bào),2023,66(12):5107-5122.

    [7] 許祥,鄒志輝.地震早至波深度學(xué)習(xí)速度建模[C]//中國地球物理學(xué)會(huì)油氣地球物理專業(yè)委員會(huì).第五屆油氣地球物理學(xué)術(shù)年會(huì)論文集.中國海洋大學(xué)海洋地球科學(xué)學(xué)院海底科學(xué)與探測技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室深海圈層與地球系統(tǒng)教育部前沿科學(xué)中心;青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評價(jià)與探測技術(shù)功能實(shí)驗(yàn)室,2023:5.

    [8] ARAYA-POLO M, JENNINGS J, ADLER A, et al. Deep-learning tomography[J].The Leading Edge,2018,37(1):58-66.

    [9] MOSSER L, KIMMAN W, DRAMSCH J, et al. Rapid seismic domain transfer: Seismic velocity inversion and modeling using deep generative neural networks[C]//80th eage conference and exhibition. EAGE Publications BV,2018:1-5.

    [10] YANG F, MA J. Deep-learning inversion: a next generation seismic velocity-model building method[J]. Geophysics,2019,84(4):583-599.

    [11] 廖成旺,龐聰,江勇,等.基于多層全連接神經(jīng)網(wǎng)絡(luò)的6C地震波極化向量識(shí)別研究[J/OL].大地測量與地球動(dòng)力學(xué),1-11[2024-02-05].https://doi.org/10.14075/j.jgg.2023.06.158.

    [12] 馬雄,申天賜,李國發(fā),等.重構(gòu)震源振幅譜一致性的吸收參數(shù)反演方法[J].石油地球物理勘探,2024,59(1):122-132.

    [13] 胡豹,高永衛(wèi),昔華倩.基于誤差反向傳播神經(jīng)網(wǎng)絡(luò)的機(jī)彈分離軌跡預(yù)測研究[J].氣動(dòng)研究與試驗(yàn),2024,2(2):59-65.

    [14] 劉浩民,楊洪才,劉戰(zhàn),等.基于粒子群優(yōu)化算法的電弧增材制造焊道尺寸反向傳播神經(jīng)網(wǎng)絡(luò)預(yù)測模型[J].機(jī)械工程材料,2024,48(2):97-102.

    [15] 張馳,郭媛,黎明.人工神經(jīng)網(wǎng)絡(luò)模型發(fā)展及應(yīng)用綜述[J].計(jì)算機(jī)工程與應(yīng)用,2021,57(11):57-69.

    [16] 譚玉陽,于靜,馮剛,等.微地震事件初至拾取SLPEA算法[J].地球物理學(xué)報(bào),2016,59(1):185-196.

    [17] 趙俐紅,劉欣悅,支鵬遙,等.南海東北陸緣的地殼速度結(jié)構(gòu)及其構(gòu)造意義:來自廣角地震剖面的約束[J].地質(zhì)學(xué)報(bào),2023, 97(6):1766-1779.

    [18] 馬月寧,鄭晶. 基于深度神經(jīng)網(wǎng)絡(luò)的地震速度模型反演[C]//中國地球科學(xué)聯(lián)合學(xué)術(shù)年會(huì).2021年中國地球科學(xué)聯(lián)合學(xué)術(shù)年會(huì)論文集(三十六)-專題一百零六 地球科學(xué)大數(shù)據(jù)與人工智能、專題一百零七 地球生物學(xué)與天體生物學(xué)、專題一百零八 古脊椎動(dòng)物學(xué)最新進(jìn)展.北京伯通電子出版社,2021:4.

    基金項(xiàng)目:江蘇省大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目(202310293095Y,202210293098Y)

    *通信作者:余瀚(1984-),男,博士,副教授。研究方向?yàn)榉囱莩上瘛C(jī)器學(xué)習(xí)。

    亚洲自拍偷在线| 春色校园在线视频观看| 少妇熟女欧美另类| 久久午夜亚洲精品久久| 国模一区二区三区四区视频| av免费在线看不卡| 国产一区亚洲一区在线观看| 国内揄拍国产精品人妻在线| 亚洲av.av天堂| 国产精品亚洲一级av第二区| 欧美日韩一区二区视频在线观看视频在线 | 国产欧美日韩精品一区二区| 午夜福利高清视频| 最新中文字幕久久久久| 亚州av有码| 一区二区三区免费毛片| 真实男女啪啪啪动态图| 熟女电影av网| 亚洲熟妇中文字幕五十中出| 欧美一区二区亚洲| 久久午夜福利片| 男女那种视频在线观看| 最近中文字幕高清免费大全6| 日日撸夜夜添| 一级毛片久久久久久久久女| 精品午夜福利在线看| 搡老岳熟女国产| 久久久久久久久久黄片| 免费看美女性在线毛片视频| 亚洲国产欧洲综合997久久,| 俄罗斯特黄特色一大片| 99久久无色码亚洲精品果冻| 天美传媒精品一区二区| 欧美性猛交黑人性爽| 久久天躁狠狠躁夜夜2o2o| 春色校园在线视频观看| 99热6这里只有精品| 特大巨黑吊av在线直播| 美女 人体艺术 gogo| 免费看av在线观看网站| 日本一本二区三区精品| 免费看光身美女| 97超碰精品成人国产| 日韩成人伦理影院| 最近中文字幕高清免费大全6| 亚洲欧美清纯卡通| 男人和女人高潮做爰伦理| 中文字幕精品亚洲无线码一区| 国产v大片淫在线免费观看| 日本在线视频免费播放| 亚洲人成网站高清观看| 亚洲欧美成人精品一区二区| 不卡视频在线观看欧美| 免费av不卡在线播放| 国产高清视频在线播放一区| 97人妻精品一区二区三区麻豆| 国产午夜精品久久久久久一区二区三区 | 天堂√8在线中文| 亚洲四区av| 亚洲在线观看片| 日本黄大片高清| 亚洲熟妇中文字幕五十中出| 国产视频内射| 99久久九九国产精品国产免费| av国产免费在线观看| 天美传媒精品一区二区| 日韩欧美国产在线观看| 国产爱豆传媒在线观看| 日韩欧美免费精品| 成人亚洲精品av一区二区| 亚洲最大成人中文| 国产精品国产三级国产av玫瑰| 国产亚洲av嫩草精品影院| 精品久久久久久久久久免费视频| 大型黄色视频在线免费观看| 日本免费a在线| 久久精品国产99精品国产亚洲性色| 国产午夜精品论理片| 日本与韩国留学比较| 少妇裸体淫交视频免费看高清| 久久久国产成人精品二区| 日日摸夜夜添夜夜添av毛片| 国产精品福利在线免费观看| 亚洲国产精品成人综合色| 欧洲精品卡2卡3卡4卡5卡区| 日本黄大片高清| 亚洲国产欧美人成| 色视频www国产| 欧美人与善性xxx| 精华霜和精华液先用哪个| 亚洲国产精品成人综合色| 欧美国产日韩亚洲一区| 亚洲中文字幕日韩| 嫩草影院精品99| 亚洲一区高清亚洲精品| 国产老妇女一区| 久久中文看片网| 久久久色成人| 国产精品99久久久久久久久| av免费在线看不卡| 麻豆乱淫一区二区| 黄色一级大片看看| 国产欧美日韩精品亚洲av| 最近2019中文字幕mv第一页| 欧美激情久久久久久爽电影| 久久久久久伊人网av| 亚洲av免费高清在线观看| 99热精品在线国产| 亚洲精品日韩av片在线观看| 热99在线观看视频| 日日摸夜夜添夜夜爱| 亚洲第一电影网av| 亚洲激情五月婷婷啪啪| 夜夜爽天天搞| 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 嫩草影院新地址| 特大巨黑吊av在线直播| 1000部很黄的大片| 岛国在线免费视频观看| 亚洲欧美日韩高清在线视频| 国产精品无大码| 嫩草影院入口| 真人做人爱边吃奶动态| 日本 av在线| 插阴视频在线观看视频| 波多野结衣巨乳人妻| 国产 一区精品| 欧美成人a在线观看| 丝袜喷水一区| 国产一区二区三区av在线 | 国产av一区在线观看免费| 国内精品美女久久久久久| 色哟哟·www| 一级黄片播放器| 久久久久国产精品人妻aⅴ院| 日日撸夜夜添| 美女高潮的动态| 欧美不卡视频在线免费观看| 国产一区亚洲一区在线观看| 成年女人毛片免费观看观看9| 伦理电影大哥的女人| 久久久久久久久大av| 人妻少妇偷人精品九色| 老熟妇乱子伦视频在线观看| 长腿黑丝高跟| 99久国产av精品| 女同久久另类99精品国产91| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 亚洲第一区二区三区不卡| 麻豆国产97在线/欧美| 夜夜夜夜夜久久久久| 三级经典国产精品| 欧美性猛交黑人性爽| 亚洲av不卡在线观看| 全区人妻精品视频| 国产黄a三级三级三级人| 久久精品国产亚洲网站| 亚洲成人av在线免费| 国产黄a三级三级三级人| 欧美中文日本在线观看视频| 身体一侧抽搐| av视频在线观看入口| 国产aⅴ精品一区二区三区波| 不卡一级毛片| 亚洲av熟女| 在线播放国产精品三级| 欧美日韩国产亚洲二区| 在线观看av片永久免费下载| 看非洲黑人一级黄片| 日韩欧美三级三区| a级毛片免费高清观看在线播放| 一级黄片播放器| 免费看光身美女| 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 尤物成人国产欧美一区二区三区| 国产老妇女一区| 五月伊人婷婷丁香| 亚洲av一区综合| 三级毛片av免费| 老司机福利观看| 成人高潮视频无遮挡免费网站| 97人妻精品一区二区三区麻豆| 色哟哟·www| av女优亚洲男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 少妇熟女aⅴ在线视频| 美女免费视频网站| 国产欧美日韩精品亚洲av| 国模一区二区三区四区视频| 中出人妻视频一区二区| 免费av观看视频| 亚洲一区二区三区色噜噜| 男人舔女人下体高潮全视频| 亚洲最大成人av| 久久精品夜色国产| 国产av麻豆久久久久久久| 日韩欧美 国产精品| 亚洲自偷自拍三级| 亚洲av免费高清在线观看| 日韩大尺度精品在线看网址| 国内精品美女久久久久久| 长腿黑丝高跟| 亚洲最大成人手机在线| 国产69精品久久久久777片| 精品国内亚洲2022精品成人| 日韩人妻高清精品专区| 久久精品综合一区二区三区| 在线免费观看不下载黄p国产| 日韩大尺度精品在线看网址| 亚洲国产精品合色在线| 欧美又色又爽又黄视频| 久久人人爽人人爽人人片va| 亚洲四区av| 国产在视频线在精品| 人人妻人人看人人澡| 亚洲电影在线观看av| 国产伦精品一区二区三区视频9| 美女xxoo啪啪120秒动态图| 插阴视频在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 一个人看视频在线观看www免费| 此物有八面人人有两片| 久久中文看片网| 校园人妻丝袜中文字幕| 亚洲欧美中文字幕日韩二区| 中文字幕人妻熟人妻熟丝袜美| 白带黄色成豆腐渣| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲av香蕉五月| 免费观看人在逋| 最新在线观看一区二区三区| 日韩欧美国产在线观看| 国产精品电影一区二区三区| 午夜福利在线观看免费完整高清在 | 小说图片视频综合网站| av视频在线观看入口| av专区在线播放| 久久久成人免费电影| a级毛色黄片| 国产成人精品久久久久久| 日本a在线网址| 又爽又黄无遮挡网站| 国内精品美女久久久久久| 日韩av在线大香蕉| 欧美激情国产日韩精品一区| 欧美一区二区国产精品久久精品| 亚洲在线观看片| 午夜a级毛片| 少妇人妻精品综合一区二区 | 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看| 少妇熟女aⅴ在线视频| 精品久久国产蜜桃| 人妻丰满熟妇av一区二区三区| 色噜噜av男人的天堂激情| 一进一出抽搐动态| 精品一区二区三区人妻视频| 日日摸夜夜添夜夜爱| 国产亚洲精品综合一区在线观看| 国产一区二区在线av高清观看| 尤物成人国产欧美一区二区三区| 午夜a级毛片| 夜夜夜夜夜久久久久| 我的老师免费观看完整版| 99久久无色码亚洲精品果冻| 久久久久国产精品人妻aⅴ院| 国产黄a三级三级三级人| 精品无人区乱码1区二区| 波多野结衣巨乳人妻| 女同久久另类99精品国产91| 五月伊人婷婷丁香| 亚洲av熟女| 成年女人毛片免费观看观看9| 久久久久久久久久成人| 亚洲专区国产一区二区| 久久国产乱子免费精品| 毛片一级片免费看久久久久| 日本在线视频免费播放| 亚洲最大成人中文| 亚洲精品在线观看二区| 长腿黑丝高跟| 午夜老司机福利剧场| 亚洲在线自拍视频| 国产熟女欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 看免费成人av毛片| 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 日韩精品青青久久久久久| 午夜免费激情av| 国产白丝娇喘喷水9色精品| 俄罗斯特黄特色一大片| 91在线观看av| 91久久精品国产一区二区三区| 十八禁网站免费在线| 久久九九热精品免费| 黄色日韩在线| 亚洲精品影视一区二区三区av| 成人永久免费在线观看视频| 亚洲av成人av| 国产69精品久久久久777片| 一个人看视频在线观看www免费| 日韩av不卡免费在线播放| 国产黄a三级三级三级人| 免费观看的影片在线观看| 欧美国产日韩亚洲一区| 波野结衣二区三区在线| avwww免费| 国产黄色小视频在线观看| 99热全是精品| 中文资源天堂在线| 神马国产精品三级电影在线观看| 成人漫画全彩无遮挡| 女的被弄到高潮叫床怎么办| 久久久久国产网址| 男人舔女人下体高潮全视频| 午夜精品在线福利| 啦啦啦韩国在线观看视频| 欧美又色又爽又黄视频| 日日干狠狠操夜夜爽| 亚洲国产精品国产精品| 1000部很黄的大片| 国产黄色小视频在线观看| 成人一区二区视频在线观看| 有码 亚洲区| 精品久久久噜噜| 亚洲一级一片aⅴ在线观看| 久久精品国产亚洲av天美| 日本色播在线视频| 久久草成人影院| 国产三级在线视频| 欧美3d第一页| 亚洲成人中文字幕在线播放| 晚上一个人看的免费电影| 无遮挡黄片免费观看| 亚洲成a人片在线一区二区| 此物有八面人人有两片| 在线看三级毛片| 亚洲精品亚洲一区二区| 最近的中文字幕免费完整| 一进一出抽搐gif免费好疼| 麻豆国产av国片精品| 国产极品精品免费视频能看的| 亚洲精品国产av成人精品 | 美女cb高潮喷水在线观看| 欧美激情久久久久久爽电影| 国产一区二区亚洲精品在线观看| 十八禁国产超污无遮挡网站| 亚洲精品色激情综合| aaaaa片日本免费| 人人妻人人澡人人爽人人夜夜 | 久久久久免费精品人妻一区二区| 嫩草影院精品99| 亚洲成av人片在线播放无| 亚洲国产欧美人成| 日韩在线高清观看一区二区三区| 欧美一区二区亚洲| 少妇熟女aⅴ在线视频| 搡老妇女老女人老熟妇| 久久久久国产精品人妻aⅴ院| 国产高清有码在线观看视频| 国产毛片a区久久久久| 99国产精品一区二区蜜桃av| 黄色日韩在线| 日本一二三区视频观看| 免费在线观看影片大全网站| 一进一出抽搐动态| 亚洲精品一区av在线观看| 久久99热6这里只有精品| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| 最好的美女福利视频网| 欧美成人精品欧美一级黄| 搡老岳熟女国产| 色av中文字幕| 国产成人91sexporn| 国内久久婷婷六月综合欲色啪| 久久久成人免费电影| 亚洲国产精品国产精品| 国产精品人妻久久久影院| 天堂网av新在线| 国产av不卡久久| 亚洲三级黄色毛片| 欧美性猛交黑人性爽| 18+在线观看网站| 俄罗斯特黄特色一大片| 久久久久国产网址| 中国美女看黄片| 精品不卡国产一区二区三区| 亚洲欧美日韩高清专用| 男人舔奶头视频| a级毛色黄片| 欧美日韩在线观看h| 蜜桃亚洲精品一区二区三区| 久久婷婷人人爽人人干人人爱| 99久久成人亚洲精品观看| 国产精品久久久久久精品电影| 国产av一区在线观看免费| 美女黄网站色视频| 国产av麻豆久久久久久久| 九九热线精品视视频播放| 人人妻人人看人人澡| 99九九线精品视频在线观看视频| 黄色一级大片看看| 久久久久九九精品影院| 极品教师在线视频| 亚洲自拍偷在线| 亚洲欧美精品综合久久99| 少妇被粗大猛烈的视频| 乱码一卡2卡4卡精品| 亚洲av免费在线观看| 亚洲在线观看片| 国产精品三级大全| 国产在线男女| 亚洲av一区综合| 大香蕉久久网| 日本一二三区视频观看| 天天躁夜夜躁狠狠久久av| 在线看三级毛片| 精华霜和精华液先用哪个| 国产久久久一区二区三区| 国产私拍福利视频在线观看| 美女被艹到高潮喷水动态| 中国国产av一级| 特级一级黄色大片| 一个人观看的视频www高清免费观看| 在线观看66精品国产| 午夜福利视频1000在线观看| 91午夜精品亚洲一区二区三区| 成年免费大片在线观看| 搡女人真爽免费视频火全软件 | 少妇被粗大猛烈的视频| 最近最新中文字幕大全电影3| 免费高清视频大片| 亚洲无线观看免费| 天天躁日日操中文字幕| 寂寞人妻少妇视频99o| 欧美xxxx性猛交bbbb| 日本-黄色视频高清免费观看| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av| 在线观看av片永久免费下载| av在线天堂中文字幕| 国产精品一区二区免费欧美| av国产免费在线观看| 可以在线观看的亚洲视频| 三级经典国产精品| 亚洲国产精品成人综合色| 国产欧美日韩精品一区二区| 精品国产三级普通话版| 又爽又黄无遮挡网站| 日日啪夜夜撸| 亚洲欧美中文字幕日韩二区| 你懂的网址亚洲精品在线观看 | 日韩一本色道免费dvd| 日本与韩国留学比较| 午夜激情欧美在线| 一a级毛片在线观看| 国产 一区 欧美 日韩| 成年女人永久免费观看视频| 亚洲,欧美,日韩| 夜夜爽天天搞| 12—13女人毛片做爰片一| 日本黄色片子视频| 久久久久久伊人网av| 大又大粗又爽又黄少妇毛片口| 中文字幕熟女人妻在线| 亚洲色图av天堂| 久久午夜亚洲精品久久| 国产三级在线视频| av黄色大香蕉| 国产精品三级大全| 看免费成人av毛片| 老司机影院成人| 欧美日本亚洲视频在线播放| 日本在线视频免费播放| 国产黄a三级三级三级人| 欧美+日韩+精品| 午夜福利在线观看吧| 老熟妇仑乱视频hdxx| 男人狂女人下面高潮的视频| 国产精品国产三级国产av玫瑰| 小蜜桃在线观看免费完整版高清| 日韩av在线大香蕉| 国产高潮美女av| 国产单亲对白刺激| 美女黄网站色视频| 中文资源天堂在线| 午夜老司机福利剧场| 久久精品人妻少妇| 国产91av在线免费观看| 亚洲熟妇熟女久久| 国产探花极品一区二区| 美女cb高潮喷水在线观看| 亚洲精品在线观看二区| 久久精品国产鲁丝片午夜精品| 国产成人福利小说| 国产毛片a区久久久久| 男人和女人高潮做爰伦理| 亚洲乱码一区二区免费版| 久久精品国产鲁丝片午夜精品| 精品久久久久久久久av| 国内精品一区二区在线观看| 麻豆国产av国片精品| 深夜a级毛片| av在线老鸭窝| 中国国产av一级| 在线观看66精品国产| 国产麻豆成人av免费视频| 黄片wwwwww| 欧美一级a爱片免费观看看| 听说在线观看完整版免费高清| 黄色日韩在线| 在线免费观看不下载黄p国产| 成人国产麻豆网| av专区在线播放| 69av精品久久久久久| 身体一侧抽搐| 特大巨黑吊av在线直播| 哪里可以看免费的av片| 黄色一级大片看看| 久99久视频精品免费| 亚洲精品成人久久久久久| 又黄又爽又刺激的免费视频.| 亚洲性夜色夜夜综合| 国产高清不卡午夜福利| 久久久久性生活片| 免费av毛片视频| 赤兔流量卡办理| 亚洲精品影视一区二区三区av| 大型黄色视频在线免费观看| 久久精品国产亚洲av天美| 久久久精品94久久精品| 欧美最新免费一区二区三区| 精品无人区乱码1区二区| 久久久成人免费电影| 在线国产一区二区在线| 欧美成人a在线观看| 嫩草影视91久久| 亚洲第一电影网av| 中文字幕熟女人妻在线| 国产精品久久久久久久电影| 99热这里只有是精品在线观看| 国产精品,欧美在线| 亚洲成av人片在线播放无| 一级毛片我不卡| 嫩草影视91久久| 尤物成人国产欧美一区二区三区| 黄色视频,在线免费观看| 高清午夜精品一区二区三区 | 亚洲,欧美,日韩| 午夜精品在线福利| 看非洲黑人一级黄片| 国产激情偷乱视频一区二区| 六月丁香七月| 国产男靠女视频免费网站| 十八禁国产超污无遮挡网站| 最近中文字幕高清免费大全6| 两个人的视频大全免费| 内地一区二区视频在线| 国产精品爽爽va在线观看网站| 国产老妇女一区| 女人被狂操c到高潮| 日本精品一区二区三区蜜桃| 又黄又爽又免费观看的视频| 最近的中文字幕免费完整| 91麻豆精品激情在线观看国产| 91狼人影院| 1024手机看黄色片| 热99在线观看视频| 亚洲国产日韩欧美精品在线观看| 男插女下体视频免费在线播放| 美女大奶头视频| 国产av麻豆久久久久久久| 18禁在线无遮挡免费观看视频 | 国产在视频线在精品| 观看美女的网站| 日韩欧美 国产精品| 久久精品国产亚洲av天美| 听说在线观看完整版免费高清| 成年女人看的毛片在线观看| 色av中文字幕| or卡值多少钱| 国产一级毛片七仙女欲春2| 亚洲av电影不卡..在线观看| 淫秽高清视频在线观看| 精品熟女少妇av免费看| 午夜精品国产一区二区电影 | 午夜福利在线观看吧| 村上凉子中文字幕在线| 免费人成在线观看视频色| 久久午夜福利片| 少妇人妻一区二区三区视频| 日韩一区二区视频免费看| 日韩高清综合在线| 五月伊人婷婷丁香| 人人妻人人看人人澡| 日本欧美国产在线视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品人妻久久久久久| 亚洲av免费高清在线观看| 亚洲欧美成人综合另类久久久 | 亚洲av中文字字幕乱码综合| 亚洲va在线va天堂va国产| www日本黄色视频网| 18+在线观看网站| 精品欧美国产一区二区三| 男女做爰动态图高潮gif福利片| 搡老熟女国产l中国老女人| 在线免费十八禁| 国产精品一及| 大香蕉久久网|