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

    基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法*

    2010-09-08 05:58:42張世鑫印興耀孔國英張繁昌
    中國海上油氣 2010年5期
    關(guān)鍵詞:反射系數(shù)反演巖石

    張世鑫印興耀孔國英張繁昌

    (1.中國石油大學(xué)(華東)地球資源與信息學(xué)院; 2.中海油研究總院)

    基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法*

    張世鑫1印興耀1孔國英2張繁昌1

    (1.中國石油大學(xué)(華東)地球資源與信息學(xué)院; 2.中海油研究總院)

    由于反射系數(shù)線性近似公式受到小角度入射和界面彈性參數(shù)差異小的基本假設(shè)的限制,常規(guī)AVO反演在大角度入射和界面彈性參數(shù)差異大的情況下存在嚴重不確定性。提出了基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法,該方法將研究工區(qū)的特定巖石物理條件作為先驗軟約束融合到反射系數(shù)近似公式中,以最優(yōu)反射系數(shù)近似為基礎(chǔ)實現(xiàn)了最終的三參數(shù)同步反演。從應(yīng)用效果看,基于巖石物理模型的最優(yōu)化三參數(shù)同步反演具有更強的工區(qū)適用性,減小了常規(guī)反射系數(shù)近似公式的基本假設(shè)對反演結(jié)果的影響,而且正交基函數(shù)提高了AVO反演的可靠性。

    最優(yōu)化近似公式 巖石物理模型 AVO 三參數(shù)同步反演 反射系數(shù)

    AVO分析的理論基礎(chǔ)是Zoeppritz方程,但由于Zoeppritz方程的復(fù)雜性與物理意義的不直觀性,在實際應(yīng)用中普遍使用更為簡便的反射系數(shù)線性近似公式[1](如Aki和Richard近似、Shuey近似等)。但是,常規(guī)線性近似公式都是建立在一系列假設(shè)基礎(chǔ)之上的,主要假設(shè)有2點:一是相鄰地層介質(zhì)的彈性參數(shù)差異較小;二是入射角小于臨界角。隨著油氣勘探尋找的目的儲層越來越復(fù)雜,很多砂巖、碳酸鹽巖等高速層界面兩側(cè)介質(zhì)的彈性參數(shù)差異較大,以及對大角度地震資料的重視,只適用于小角度入射的常規(guī)線性近似公式必然會影響參數(shù)估計的精度和有效性。為解決這些問題,Causse等人[2-3]提出了基于巖石物理模型的最優(yōu)化近似公式,可以較好地克服常規(guī)反射系數(shù)近似的局限性,進行儲層類型劃分與常規(guī)AVO屬性提取;Varela等人[4]提出了基于最優(yōu)化近似公式的最優(yōu)系數(shù),可以反演裂縫密度。筆者則提出以巖石物理模型的最優(yōu)化反射系數(shù)近似公式為基礎(chǔ)進行常規(guī)彈性參數(shù)同步反演的方法,經(jīng)過理論模型分析與實際應(yīng)用,表明該方法是有效、可靠的。

    1 反演方法原理

    1.1 最優(yōu)化反射系數(shù)近似公式

    最優(yōu)化近似公式將反射系數(shù)用一個線性通式表示[2-3],即

    式(1)中:fi(θ)和Ci分別是最優(yōu)基函數(shù)和系數(shù)。

    若已知地層的巖性信息,假設(shè)該巖性地層的物性參數(shù)服從某一合理的概率分布,根據(jù)該分布構(gòu)建N個屬于該巖性地層的物性參數(shù)模型(主要包括縱橫波速度和密度);若假設(shè)有M個離散的入射角度,利用Zoeppritz方程可以得到這 N個參考模型的反射系數(shù)曲線分布(如圖1,在此假設(shè)只有含水砂巖與泥巖以及含氣砂巖與泥巖這2種類型的巖性分界面,為方便說明方法的優(yōu)越性,2種界面均屬于第一類AVO情況,實際情況要比這復(fù)雜得多),構(gòu)建一個M×N維的反射系數(shù)矩陣R(注意建立模型時M >N)。

    對反射系數(shù)矩陣R進行奇異值分解運算,即

    式(2)中:矩陣V,D和F的維數(shù)分別為N×N、N× N和M×N。根據(jù)奇異值分解的性質(zhì)可知,V為正交矩陣,D為對角陣,F的列是標準正交的。對于給定的模型,相關(guān)AVO曲線的主要特征信息隨基函數(shù)次序的增加而降低,并且考慮到奇異值遞減速度很快,忽略高次項基函數(shù)對近似精度影響很小,本文研究主要針對三項基函數(shù)近似。

    經(jīng)奇異值分解得到基函數(shù)矩陣 F后,選擇合適項數(shù)的最優(yōu)基函數(shù),結(jié)合對應(yīng)界面的準確反射系數(shù),采用最小二乘法計算得到最優(yōu)系數(shù) Ci,即

    1.2 巖石物理模型的建立

    基于最優(yōu)化反射系數(shù)進行疊前AVO反演提取屬性參數(shù)時,首先建立工區(qū)內(nèi)各種巖相的縱橫波速度與密度隨深度的變化趨勢,以此作為巖相模型的平均值;然后假設(shè)取值服從某一合理的概率分布,構(gòu)建各類巖相的巖石物理模型;最后結(jié)合研究工區(qū)的地質(zhì)規(guī)律,對巖相模型進行相應(yīng)組合后計算模型界面的準確反射系數(shù)。

    本次研究考慮到實際工區(qū)埋深較淺,采用Avseth等人提出的趨勢模型的建立方法[5-6],主要步驟如下:

    (1)建立孔隙度隨深度的變化趨勢

    Ramm和Bj?rlykke提出了一種砂巖孔隙度(φ)與深度(z)并依賴粘土的回歸模型,即

    式(4)中:φc表示臨界(即沉淀)孔隙度;α和β是回歸系數(shù),分別代表純砂巖骨架顆粒穩(wěn)定因子(Cl=0)和描述隨粘土含量增加的敏感性因子;Cl表示泥質(zhì)因子,定義為相對于穩(wěn)定骨架顆??偤康目傉惩梁俊J褂煤线m的地面臨界孔隙度值(砂巖一般為0.40~0.45,泥巖為0.6~0.8)和任意埋藏深度反演的密度測井值校正局部孔隙度-深度趨勢。

    (2)根據(jù)孔隙度隨深度的變化趨勢直接得到密度變化趨勢,即

    式(5)中:ρsat是飽和巖石密度;ρfl是流體密度;ρm是礦物密度;φ是孔隙度。對于干巖石來說,流體密度等于零。

    (3)根據(jù)Hertz-Mindlin接觸理論與Gassmann理論計算得到縱橫波速度隨深度的變化趨勢

    首先根據(jù) Hertz-Mindlin理論計算未固結(jié)沉積物的彈性模量(根據(jù)該理論可知,彈性模量是孔隙度和壓力的函數(shù)),得到彈性模量隨深度的變化趨勢,進而可以計算縱橫波速度隨深度的變化趨勢。

    Hertz-Mindlin接觸理論中定義干巖樣的彈性模量通過下式計算

    式(6)、(7)中:KHM和μHM分別表示孔隙度為φ的體積模量和剪切模量;P表示壓差,對于不同深度的有效壓力,P=g∫z0(ρb-ρl)dz,g是重力常數(shù),ρb和ρl分別是給定深度時的巖石固體密度和流體密度;K、μ和σ分別表示固體相的體積模量、剪切模量和泊松比;n表示配位數(shù)(每個顆粒的接觸顆粒數(shù)目),配位數(shù)與孔隙度之間的關(guān)系可近似為 n=20-34φ+ 14φ2。飽和巖樣的彈性模量 Ksat和μsat可以用 Gassmann方程計算得到。

    得到 Ksat,μsat和ρsat之后,便可以計算飽和巖樣的縱波速度Vp和橫波速度Vs,即

    1.3 彈性參數(shù)優(yōu)化

    最優(yōu)系數(shù)和常規(guī)近似系數(shù)之間存在一定的關(guān)系,Causse等人[2-3]研究了Aki-Richard近似公式中彈性參數(shù)的優(yōu)化問題,給出了如下公式

    式(11)中:L為所選基函數(shù)個數(shù);N為所建模型的個數(shù);hjk表示矩陣 H=W-1=VTD-1的第 j行第k列元素。

    1.4 反演方程的建立

    以矩陣形式表示最優(yōu)化反射系數(shù)為

    式(12)中:R是反射系數(shù)矩陣;F是最優(yōu)基函數(shù)矩陣;C是最優(yōu)系數(shù)矩陣。

    如果考慮褶積模型,地震記錄表示為子波與反射系數(shù)的褶積,式(12)用矩陣形式表示為

    式(13)中:Ti表示地震道的第i個采樣點;wj表示提取地震子波的第j個采樣點;Ri表示第i個采樣點的反射系數(shù)。綜合式(13)可以得到如下公式:

    由于最優(yōu)系數(shù)Ci(i=1,2,3)之間是正交的,所以在求解反演系數(shù)時不必計算協(xié)方差矩陣去相關(guān),直接采用最小二乘法計算 Ci(i=1,2,3)。得到 Ci(i=1, 2,3)之后,利用公式(11)可以進一步優(yōu)化得到彈性參數(shù)相對變化量(Rp=ΔVp/Vp,Rs=ΔVs/Vs,Rd= Δ ρ/ρ)。

    根據(jù)關(guān)系

    式(15)中:Ri表示第i層和第i+1層之間界面的反射系數(shù)。如果我們考慮 N個采樣點,則式(15)以矩陣形式表示為

    式(16)中:Li=lnZi。

    分別針對 Rp、Rs和 Rd求解矩陣(16),可以最終得到參數(shù)Vp、Vs和ρ。

    2 模型試算

    圖2 模型數(shù)據(jù)提取的彈性參數(shù)相對變化量的三維交會顯示

    以圖1所示含水砂巖與含氣砂巖(第一類AVO)為模型進行相應(yīng)彈性參數(shù)提取的試算。分別采用解析法、基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法(即本文方法)和基于常規(guī)Aki-Richard近似公式的L2模反演方法對該模型進行彈性參數(shù)相對變化量(ΔVp/Vp,ΔVs/Vs,Δ ρ/ρ)的提取并計算常規(guī)屬性;分別對計算的3個屬性進行三維交會顯示,如圖2所示。

    在圖2中,解析解表示模型彈性參數(shù)變化量的準確分布(圖2a);將圖2b、圖2c與圖2a比較可知,在滿足一定誤差條件下,基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法提取的3個彈性參數(shù)相對變化量與準確值基本吻合(彈性參數(shù)相對變化量是基于常規(guī)AVO屬性與彈性參數(shù)之間的近似關(guān)系間接計算得到的,不可避免的會引入一定誤差),但基于常規(guī)Aki-Richard近似公式的L2模反演方法計算的彈性參數(shù)相對變化量與解析解之間的誤差非常大。這說明,基于常規(guī)的AVO近似公式的反演方法無法滿足彈性參數(shù)差異大的界面以及需要大角度信息參與情況下的彈性參數(shù)的可靠反演,而本文提出的基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法可以較好地彌補這一不足,應(yīng)用該反演方法所得到的儲層的彈性參數(shù)能夠有效地表征儲層信息。

    3 實際應(yīng)用

    將該方法應(yīng)用到某勘探區(qū)塊。反演之前首先要對地震數(shù)據(jù)和測井數(shù)據(jù)進行預(yù)處理[7],對地震記錄采用去多次波、幾何擴散校正、地表一致性反褶積、振幅補償、保幅疊前偏移等技術(shù)進行處理之后得到疊前道集。

    反演中首先根據(jù)測井和巖性解釋資料將研究工區(qū)巖相劃分為4類模型(圖3),建立各種巖相隨深度變化的巖石物理模型;然后兩兩組合建立相應(yīng)的界面,對應(yīng)得到準確的反射系數(shù)矩陣;對反射系數(shù)矩陣進行奇異值分解,得到不同深度最優(yōu)反射系數(shù)的基函數(shù),最終反演得到縱波速度、橫波速度和密度等3個彈性參數(shù)體。

    圖3 根據(jù)測井資料劃分的研究工區(qū)4種巖相類型

    由于受到地球物理模型精度的限制,只能在較小的時窗內(nèi)進行彈性參數(shù)反演,本次實際應(yīng)用選擇的時間段是1.6~2.0s,該時間段內(nèi)有2個含油目的層段。圖4是其中一條測線的反演結(jié)果,圖中測井曲線是電阻率曲線,可以看到,油層處縱波速度降低,密度降低,相應(yīng)的橫波速度并沒有較大變化。通過巖石物理分析可知,上述特征符合研究工區(qū)的儲層特征,并且反演的井旁道與測井曲線在誤差允許的范圍內(nèi)匹配較好,表明了本文提出的彈性參數(shù)反演方法的有效性。

    圖4 研究工區(qū)彈性參數(shù)反演剖面

    4 結(jié)論

    與常規(guī)AVO彈性參數(shù)反演方法相比,本文提出的基于巖石物理模型的最優(yōu)化AVO三參數(shù)同步反演方法克服了常規(guī)反射系數(shù)近似公式存在的缺陷,適用于近臨界角入射以及界面物性差異大的情況。由于本文方法中基函數(shù)是正交的,所以不同系數(shù)所包含的噪聲不相關(guān),而且有著相同的幅度和方差(在近似公式中選擇的系數(shù)少,相當于利用低通濾波器進行去噪處理),因此應(yīng)用本文方法估算的參數(shù)的抗噪能力更強。另外,由于巖石物理模型作為先驗信息在公式推導(dǎo)中起到了軟約束的作用,因此所得到的最優(yōu)反射系數(shù)近似公式具有工區(qū)適定性,提高了AVO定性和定量解釋的精確性,最終得到了穩(wěn)定可靠的反演結(jié)果。

    需要指出的是,巖石物理模型建立的正確與否對本文方法起到至關(guān)重要的作用,而且影響因素的多樣性和復(fù)雜性制約了巖石物理模型建立的可靠性,因此在研究中必須結(jié)合研究工區(qū)特征靈活地選擇最為合適的模型建立方法。

    [1] TARANTALA A.Inversion problem theory and methods for model parameter estimation[M].Philadelphia:Society for Industry and Applied Mathematics,2004.

    [2] CAUSSE E,RIEDE M,WIJ INGAARDEN A J,et al.Amplitude analysis with an optimal model-based linear AVO approximation:PartⅠ—Theory[J].Geophysics,2007,72:59-69.

    [3] CAUSSE E,RIEDE M,WIJ INGAARDEN A J,et al.Amplitude analysis with an optimal model-based linear AVO approximation:PartⅡ—Field data example[J].Geophysics,2007, 72:71-79.

    [4] VARELA I,SONJA M,LI Xiangyang.Fracture density inversion from a physical geological model using azimuthal AVO with optimal basis function:79th annual international meeting [C].SEG,Expanded Abstracts,2009:2075-2079.

    [5] AVSETH P,FLESCHE H,WIJ INGAARDEN A J.AVO classification of lithology and pore fluids constrained by rock physics depth trends[J].The Leading Edge,2003,22:1004-1011.

    [6] AVSETH P,MU KERJ I T,MAVKO G,Quantitative seismic interpretation[M].Cambridge:Cambridge Univ.Press,2005.

    [7] SHERIFF R E.Factors affecting seismic amplitudes[J].Geophys Prospect,1975,23:125-138.收稿日期:2010-02-03 改回日期:2010-05-24

    (編輯:周雯雯)

    Abstract:Because the reflection coefficient approximation is limited by the basic assumptions of a low-angle incidence and small differentiation in interfacialelasticity coefficient,any conventional AVO inversion will be seriously uncertain with a high-angle incidence and large differentiation in interfacial elasticity coefficient.Therefore,an optimal simultaneous inversion method for AVO with three parameters has been developeds on a petrophysics-model basis,in which the special petrophysical conditions in a study area can be integrated into the reflection coefficient approximation as a prior soft constraint,and a simultaneous inversion with three parameters is implemented on a basis of the optimal reflection coefficient approximation. The application results have shown that this AVO inversion is higher in adaptability for study areas, and can decreasse the impacts of the basic assumptions for the conventional reflection coefficient approximation on the inversion result,and that the orthogonal basic functions can make the AVO inversion more reliable.

    Key words:optimal approximation;petrophysics model;AVO;simultaneous inversion with three parameters;reflection coefficient

    An optimal simultaneous inversion method for AVO with three parameters on a petrophysics-model basis

    Zhang Shixin1Yin Xingyao1Kong Guoying2Zhang Fanchang1

    (1.College of Geo-Resources and Inf ormation,China University ofPetroleum(East China),Shandong,266555;
    2.CNOOC Research Institute,Beijing,100027)

    *國家重大專項(2008ZX05030-04-001-001)和中國石油大學(xué)優(yōu)秀博士學(xué)位論文培育資助項目(z10-01)部分研究成果。

    張世鑫,男,中國石油大學(xué)(華東)在讀博士研究生,主要從事地震反演、儲層預(yù)測方面的研究。E-mail:zsxvegata@126.com。

    猜你喜歡
    反射系數(shù)反演巖石
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    第五章 巖石小專家
    3深源巖石
    多道隨機稀疏反射系數(shù)反演
    石油物探(2020年6期)2020-11-25 02:38:46
    一種叫做煤炭的巖石
    海藻與巖石之間
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    球面波PP反射系數(shù)的頻變特征研究
    沙質(zhì)沉積物反射系數(shù)的寬帶測量方法
    成人三级黄色视频| 成人国产综合亚洲| 午夜福利在线观看免费完整高清在 | 国产极品精品免费视频能看的| av国产免费在线观看| 成人特级av手机在线观看| 97热精品久久久久久| 老司机福利观看| 成人一区二区视频在线观看| av天堂中文字幕网| 久久亚洲精品不卡| 国产精品久久视频播放| av黄色大香蕉| 婷婷精品国产亚洲av| 久久国产乱子免费精品| 欧美成人性av电影在线观看| 中文字幕免费在线视频6| 可以在线观看毛片的网站| 小说图片视频综合网站| 99热网站在线观看| 桃红色精品国产亚洲av| 内地一区二区视频在线| 亚洲国产精品合色在线| 精品久久久久久,| 国产一区二区在线av高清观看| 99久久精品国产国产毛片| 欧美国产日韩亚洲一区| 全区人妻精品视频| 色吧在线观看| 成年人黄色毛片网站| 中文在线观看免费www的网站| 又黄又爽又免费观看的视频| 亚洲成人中文字幕在线播放| 亚洲黑人精品在线| 国产精品永久免费网站| 一区二区三区激情视频| 黄色女人牲交| 欧美成人a在线观看| 极品教师在线视频| 欧美绝顶高潮抽搐喷水| 成人三级黄色视频| 亚洲五月天丁香| 一卡2卡三卡四卡精品乱码亚洲| 伦精品一区二区三区| 波多野结衣高清无吗| 有码 亚洲区| 禁无遮挡网站| 人妻丰满熟妇av一区二区三区| av福利片在线观看| 国语自产精品视频在线第100页| 韩国av在线不卡| 夜夜夜夜夜久久久久| 亚洲精品乱码久久久v下载方式| 亚洲精品乱码久久久v下载方式| 桃色一区二区三区在线观看| 88av欧美| 老熟妇乱子伦视频在线观看| 伊人久久精品亚洲午夜| 精品不卡国产一区二区三区| 大型黄色视频在线免费观看| 男女啪啪激烈高潮av片| 99热网站在线观看| 久久久久久国产a免费观看| 亚洲一级一片aⅴ在线观看| 国产精品爽爽va在线观看网站| 两个人视频免费观看高清| 真人一进一出gif抽搐免费| 免费看av在线观看网站| 男女那种视频在线观看| 日韩强制内射视频| 99热只有精品国产| 美女xxoo啪啪120秒动态图| 嫁个100分男人电影在线观看| 99热这里只有是精品在线观看| 国产在视频线在精品| 中文字幕久久专区| 欧美+日韩+精品| 日韩在线高清观看一区二区三区 | 精品国内亚洲2022精品成人| 亚洲国产精品sss在线观看| 国产欧美日韩精品亚洲av| 国产精品精品国产色婷婷| 丝袜美腿在线中文| 51国产日韩欧美| 乱码一卡2卡4卡精品| 91在线观看av| 国产熟女欧美一区二区| 久久国产乱子免费精品| 在线免费观看的www视频| or卡值多少钱| 欧美在线一区亚洲| 亚洲av免费高清在线观看| 麻豆精品久久久久久蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 日韩大尺度精品在线看网址| 国产精品99久久久久久久久| 在线a可以看的网站| 一区二区三区高清视频在线| 精品久久久久久久久亚洲 | 亚洲精品久久国产高清桃花| 亚洲人与动物交配视频| 亚洲久久久久久中文字幕| 人妻丰满熟妇av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 一级黄色大片毛片| 日日啪夜夜撸| 国产精品一及| av在线天堂中文字幕| 高清毛片免费观看视频网站| 欧美丝袜亚洲另类 | 久久99热这里只有精品18| 舔av片在线| 国内精品美女久久久久久| 91麻豆精品激情在线观看国产| 色播亚洲综合网| 亚洲男人的天堂狠狠| 黄色一级大片看看| 国产成年人精品一区二区| 亚洲欧美清纯卡通| 国产伦精品一区二区三区四那| 99久久久亚洲精品蜜臀av| 亚洲av日韩精品久久久久久密| 国产中年淑女户外野战色| 特级一级黄色大片| 欧美bdsm另类| 亚洲av免费高清在线观看| 中文字幕av在线有码专区| 日本与韩国留学比较| 欧美成人性av电影在线观看| 大又大粗又爽又黄少妇毛片口| 国产高潮美女av| 亚洲av日韩精品久久久久久密| 在线看三级毛片| netflix在线观看网站| 国产男靠女视频免费网站| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 国语自产精品视频在线第100页| 国产一区二区在线观看日韩| 久久九九热精品免费| 国产伦精品一区二区三区四那| 日韩精品青青久久久久久| 免费看美女性在线毛片视频| 特大巨黑吊av在线直播| 国产综合懂色| 国语自产精品视频在线第100页| 91久久精品电影网| 亚洲精品日韩av片在线观看| 欧美黑人欧美精品刺激| 舔av片在线| 中亚洲国语对白在线视频| 国产成人aa在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产综合懂色| 不卡视频在线观看欧美| 亚洲真实伦在线观看| 色综合亚洲欧美另类图片| 乱人视频在线观看| 别揉我奶头~嗯~啊~动态视频| 中出人妻视频一区二区| 国产午夜精品论理片| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩无卡精品| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器| 精品一区二区三区av网在线观看| 欧美日韩综合久久久久久 | 狂野欧美激情性xxxx在线观看| 欧美成人a在线观看| 97超级碰碰碰精品色视频在线观看| 女的被弄到高潮叫床怎么办 | 久久国内精品自在自线图片| 91狼人影院| 精华霜和精华液先用哪个| 最后的刺客免费高清国语| 99精品在免费线老司机午夜| 搡女人真爽免费视频火全软件 | 午夜激情福利司机影院| 国产视频一区二区在线看| 精品久久久久久久久av| 国产单亲对白刺激| 亚洲成人久久性| 乱人视频在线观看| 久久久久精品国产欧美久久久| 一区二区三区高清视频在线| 国产精品久久视频播放| 欧美日韩国产亚洲二区| 久久99热这里只有精品18| 国产真实乱freesex| 久久精品国产亚洲av涩爱 | 色5月婷婷丁香| 国产精品久久久久久久电影| 九九在线视频观看精品| 精品久久久久久久久亚洲 | 长腿黑丝高跟| 欧美zozozo另类| 久久久国产成人精品二区| 国产精品不卡视频一区二区| 欧美又色又爽又黄视频| 国产私拍福利视频在线观看| 国产一区二区在线av高清观看| 狂野欧美激情性xxxx在线观看| 在线观看午夜福利视频| 午夜激情福利司机影院| 午夜久久久久精精品| 久久久久国产精品人妻aⅴ院| 国产男靠女视频免费网站| or卡值多少钱| 亚洲人成伊人成综合网2020| 久久久久久久久中文| 国产女主播在线喷水免费视频网站 | 少妇人妻精品综合一区二区 | 男女做爰动态图高潮gif福利片| videossex国产| 天堂av国产一区二区熟女人妻| 少妇人妻一区二区三区视频| 午夜福利视频1000在线观看| 亚洲国产精品合色在线| 亚洲最大成人av| 日韩欧美在线二视频| 欧美一区二区精品小视频在线| 亚洲欧美日韩卡通动漫| 久久久精品大字幕| 搡女人真爽免费视频火全软件 | 亚洲av电影不卡..在线观看| 久久欧美精品欧美久久欧美| 国产亚洲精品久久久com| 亚洲,欧美,日韩| 国产精品一区二区三区四区久久| 男人狂女人下面高潮的视频| 午夜a级毛片| 一本一本综合久久| 亚洲av二区三区四区| 亚洲精品国产成人久久av| 亚洲第一电影网av| 成人三级黄色视频| 国产中年淑女户外野战色| 在线a可以看的网站| 舔av片在线| 观看免费一级毛片| 无遮挡黄片免费观看| 日韩亚洲欧美综合| 欧美日韩国产亚洲二区| 国产高清视频在线观看网站| 高清在线国产一区| 亚洲成人精品中文字幕电影| 1024手机看黄色片| 黄色一级大片看看| 俄罗斯特黄特色一大片| 99久久精品一区二区三区| 一a级毛片在线观看| 亚洲人与动物交配视频| 精品一区二区三区人妻视频| 精品久久久久久久末码| 少妇熟女aⅴ在线视频| 久久中文看片网| 欧美丝袜亚洲另类 | 一夜夜www| 国产大屁股一区二区在线视频| 国产真实伦视频高清在线观看 | 婷婷精品国产亚洲av在线| www.www免费av| 国产伦在线观看视频一区| 欧美精品国产亚洲| 在线播放国产精品三级| 亚洲精品日韩av片在线观看| 桃色一区二区三区在线观看| 日韩人妻高清精品专区| 午夜精品在线福利| 99视频精品全部免费 在线| 欧美性感艳星| 中文字幕高清在线视频| 国产午夜精品久久久久久一区二区三区 | 亚洲人成网站在线播| 亚洲av五月六月丁香网| 日韩欧美三级三区| 国产精华一区二区三区| 久久精品国产亚洲av天美| 中国美白少妇内射xxxbb| 国产精品久久久久久亚洲av鲁大| 国产国拍精品亚洲av在线观看| 日本一二三区视频观看| 国内精品久久久久精免费| 一进一出抽搐gif免费好疼| 国产精品爽爽va在线观看网站| 久久久国产成人精品二区| 国产探花极品一区二区| 搡老岳熟女国产| 天天躁日日操中文字幕| 日韩欧美国产在线观看| 99久久久亚洲精品蜜臀av| 极品教师在线免费播放| 免费黄网站久久成人精品| 亚洲成人免费电影在线观看| 99在线人妻在线中文字幕| 一进一出抽搐动态| 搞女人的毛片| 亚洲av.av天堂| 能在线免费观看的黄片| 黄色一级大片看看| 男女边吃奶边做爰视频| 久久久国产成人精品二区| 老熟妇仑乱视频hdxx| 18禁裸乳无遮挡免费网站照片| 国内精品一区二区在线观看| 国产一区二区三区在线臀色熟女| 欧美最黄视频在线播放免费| 天堂√8在线中文| 级片在线观看| 香蕉av资源在线| 一卡2卡三卡四卡精品乱码亚洲| 国产三级在线视频| 久久精品国产亚洲网站| 亚洲 国产 在线| 久久久色成人| 国产精品一区二区性色av| 一边摸一边抽搐一进一小说| 亚洲专区国产一区二区| 中国美女看黄片| 人妻制服诱惑在线中文字幕| 一夜夜www| 亚洲国产日韩欧美精品在线观看| 我要看日韩黄色一级片| 亚洲午夜理论影院| 舔av片在线| 国产精品野战在线观看| 九色成人免费人妻av| 十八禁国产超污无遮挡网站| 久久久久久久久久黄片| 亚洲在线自拍视频| 亚洲av成人av| 国产精品电影一区二区三区| 九色成人免费人妻av| 亚洲av成人精品一区久久| 久久久久久伊人网av| 亚洲美女搞黄在线观看 | 国产伦人伦偷精品视频| 久久精品国产亚洲av涩爱 | 九色成人免费人妻av| 亚洲av成人精品一区久久| 91麻豆精品激情在线观看国产| 国产精品久久电影中文字幕| 亚洲美女视频黄频| 直男gayav资源| 成人国产麻豆网| 国产v大片淫在线免费观看| 村上凉子中文字幕在线| 日本黄色片子视频| 亚洲成人精品中文字幕电影| 久久草成人影院| 国产精品无大码| 最新在线观看一区二区三区| 国产一区二区三区av在线 | 日韩亚洲欧美综合| 亚洲成人久久性| 久久久色成人| 久久久久久九九精品二区国产| 真人一进一出gif抽搐免费| 亚洲电影在线观看av| 日韩高清综合在线| 亚洲第一区二区三区不卡| 高清毛片免费观看视频网站| 999久久久精品免费观看国产| 亚洲va在线va天堂va国产| 一级黄色大片毛片| 男插女下体视频免费在线播放| 中国美白少妇内射xxxbb| 午夜福利18| 国产熟女欧美一区二区| 亚洲精品粉嫩美女一区| 久久99热6这里只有精品| 99视频精品全部免费 在线| 国产精品人妻久久久影院| 亚洲国产色片| 国产精品一区二区三区四区免费观看 | 亚洲三级黄色毛片| 一区福利在线观看| av女优亚洲男人天堂| 麻豆一二三区av精品| 亚洲成人精品中文字幕电影| 嫁个100分男人电影在线观看| 亚洲av成人精品一区久久| 国产免费av片在线观看野外av| 91精品国产九色| 国产一区二区三区视频了| 啪啪无遮挡十八禁网站| 又爽又黄a免费视频| 免费观看在线日韩| 久久久久久久久大av| 久久草成人影院| av中文乱码字幕在线| 久久久色成人| 色综合色国产| 噜噜噜噜噜久久久久久91| 欧美最黄视频在线播放免费| 国产精品一及| 欧美中文日本在线观看视频| 久久久久久九九精品二区国产| 嫁个100分男人电影在线观看| 精品福利观看| 免费看美女性在线毛片视频| 大又大粗又爽又黄少妇毛片口| 久久精品久久久久久噜噜老黄 | av在线老鸭窝| 无遮挡黄片免费观看| 成人性生交大片免费视频hd| 网址你懂的国产日韩在线| 日韩强制内射视频| 99热这里只有是精品50| 成人国产麻豆网| 亚洲成a人片在线一区二区| 国产高潮美女av| a在线观看视频网站| 天天一区二区日本电影三级| 中文字幕av在线有码专区| 麻豆成人av在线观看| 久久6这里有精品| 亚洲,欧美,日韩| 国产一区二区亚洲精品在线观看| 黄色丝袜av网址大全| av国产免费在线观看| 国产亚洲欧美98| 日韩欧美国产一区二区入口| 三级毛片av免费| 成年女人看的毛片在线观看| 男女视频在线观看网站免费| 人妻制服诱惑在线中文字幕| 久久久久久国产a免费观看| 亚洲人成网站高清观看| 久久久成人免费电影| 亚洲男人的天堂狠狠| 91av网一区二区| 国内精品宾馆在线| 变态另类成人亚洲欧美熟女| 丰满乱子伦码专区| 淫妇啪啪啪对白视频| 丰满的人妻完整版| 亚洲精品456在线播放app | 日韩欧美免费精品| 可以在线观看毛片的网站| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩无卡精品| av女优亚洲男人天堂| 亚洲经典国产精华液单| 亚洲av一区综合| 99热这里只有是精品50| 热99在线观看视频| 美女 人体艺术 gogo| 国产精品98久久久久久宅男小说| 免费观看人在逋| 午夜视频国产福利| 日韩欧美三级三区| 成人二区视频| 亚洲精品一区av在线观看| 久久精品久久久久久噜噜老黄 | 国产 一区精品| 亚洲欧美日韩高清专用| 亚洲无线观看免费| 搡女人真爽免费视频火全软件 | 亚洲va在线va天堂va国产| 伊人久久精品亚洲午夜| 九色成人免费人妻av| av黄色大香蕉| 在线观看一区二区三区| 一级av片app| 成人精品一区二区免费| 日韩中文字幕欧美一区二区| 搡老妇女老女人老熟妇| 亚洲午夜理论影院| 国产一区二区激情短视频| 亚洲国产精品久久男人天堂| 精品久久久久久,| x7x7x7水蜜桃| 成人毛片a级毛片在线播放| 亚洲av熟女| 欧美高清成人免费视频www| 久久中文看片网| 国产单亲对白刺激| 人人妻人人澡欧美一区二区| 禁无遮挡网站| 国产视频一区二区在线看| 欧美激情在线99| 成人高潮视频无遮挡免费网站| 淫秽高清视频在线观看| 国产精品野战在线观看| 天堂√8在线中文| 亚洲在线自拍视频| 日韩欧美三级三区| 观看免费一级毛片| 日日撸夜夜添| 级片在线观看| 一进一出好大好爽视频| 内地一区二区视频在线| 久久精品夜夜夜夜夜久久蜜豆| 99久国产av精品| 少妇人妻一区二区三区视频| 久久久国产成人精品二区| 亚洲美女黄片视频| 久久精品国产亚洲av天美| 国产国拍精品亚洲av在线观看| 国产伦一二天堂av在线观看| 中亚洲国语对白在线视频| 国产精品综合久久久久久久免费| 蜜桃亚洲精品一区二区三区| 日本精品一区二区三区蜜桃| 日日摸夜夜添夜夜添av毛片 | 天天一区二区日本电影三级| 毛片女人毛片| 深夜a级毛片| 国产aⅴ精品一区二区三区波| 国产精品三级大全| 白带黄色成豆腐渣| 国产真实乱freesex| 美女免费视频网站| 亚洲在线观看片| 18禁在线播放成人免费| 熟女电影av网| 麻豆av噜噜一区二区三区| 亚洲av.av天堂| a级一级毛片免费在线观看| 亚洲av第一区精品v没综合| 日韩中文字幕欧美一区二区| 美女高潮的动态| 美女黄网站色视频| 尾随美女入室| 嫁个100分男人电影在线观看| 亚洲精品在线观看二区| 中文字幕人妻熟人妻熟丝袜美| a级毛片a级免费在线| 国内精品一区二区在线观看| 在线免费十八禁| 欧美日韩精品成人综合77777| 天堂网av新在线| 精品免费久久久久久久清纯| 不卡一级毛片| 特大巨黑吊av在线直播| 九色成人免费人妻av| 又黄又爽又免费观看的视频| 成人av在线播放网站| 男女下面进入的视频免费午夜| 亚洲欧美日韩高清在线视频| 两个人视频免费观看高清| av在线蜜桃| 男女下面进入的视频免费午夜| 精品国内亚洲2022精品成人| 亚洲av熟女| 欧美一区二区精品小视频在线| 欧美另类亚洲清纯唯美| 亚洲内射少妇av| 伦理电影大哥的女人| 一进一出抽搐动态| 色综合亚洲欧美另类图片| 在线播放无遮挡| 春色校园在线视频观看| 亚洲精品日韩av片在线观看| 国产探花在线观看一区二区| 亚洲欧美日韩卡通动漫| 22中文网久久字幕| 直男gayav资源| 最近视频中文字幕2019在线8| 欧美区成人在线视频| 欧美色欧美亚洲另类二区| 久久精品国产清高在天天线| 国产精品亚洲美女久久久| 一本一本综合久久| 男女下面进入的视频免费午夜| 麻豆一二三区av精品| 欧美成人一区二区免费高清观看| 精品人妻1区二区| 又黄又爽又刺激的免费视频.| 一级毛片久久久久久久久女| 麻豆成人av在线观看| 国产亚洲欧美98| 五月玫瑰六月丁香| 国产精品一区二区三区四区免费观看 | 99久久无色码亚洲精品果冻| 久久热精品热| 日韩欧美在线二视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美最新免费一区二区三区| 亚洲欧美激情综合另类| 深夜精品福利| 老司机深夜福利视频在线观看| av女优亚洲男人天堂| 国产精品av视频在线免费观看| 91在线观看av| 在现免费观看毛片| 久久精品夜夜夜夜夜久久蜜豆| 在线国产一区二区在线| 国产精华一区二区三区| 欧美bdsm另类| 人妻丰满熟妇av一区二区三区| 国产精品精品国产色婷婷| 久久久久久久久久成人| 成人无遮挡网站| 国产精品综合久久久久久久免费| 啦啦啦啦在线视频资源| 日韩欧美一区二区三区在线观看| 国内精品久久久久精免费| 嫩草影院精品99| h日本视频在线播放| 久久久国产成人免费| 日韩欧美精品v在线| 精品人妻一区二区三区麻豆 | eeuss影院久久| 中文字幕人妻熟人妻熟丝袜美| 有码 亚洲区| 不卡视频在线观看欧美| a在线观看视频网站| 91午夜精品亚洲一区二区三区 | 亚洲精品日韩av片在线观看| 色尼玛亚洲综合影院| 国产欧美日韩精品亚洲av| 色5月婷婷丁香| 一边摸一边抽搐一进一小说|