• <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ù)的寬帶測量方法
    欧美精品啪啪一区二区三区 | 日韩视频一区二区在线观看| 欧美日韩黄片免| 国产视频一区二区在线看| 亚洲午夜精品一区,二区,三区| 欧美性长视频在线观看| 热99国产精品久久久久久7| 亚洲国产精品成人久久小说| 777久久人妻少妇嫩草av网站| 九色亚洲精品在线播放| 亚洲av电影在线观看一区二区三区| 成人三级做爰电影| 男女高潮啪啪啪动态图| 高清av免费在线| 免费久久久久久久精品成人欧美视频| 成人国语在线视频| av又黄又爽大尺度在线免费看| 大片电影免费在线观看免费| 国产97色在线日韩免费| av网站在线播放免费| 亚洲av男天堂| 首页视频小说图片口味搜索| 亚洲第一欧美日韩一区二区三区 | 国产成人精品无人区| 1024视频免费在线观看| 麻豆国产av国片精品| 视频区图区小说| 99热网站在线观看| 婷婷丁香在线五月| 国产精品免费大片| 中国国产av一级| 久久国产精品人妻蜜桃| 欧美精品av麻豆av| 久久青草综合色| 女人爽到高潮嗷嗷叫在线视频| 热99国产精品久久久久久7| 久久人人爽人人片av| 韩国高清视频一区二区三区| 老熟妇仑乱视频hdxx| 国产淫语在线视频| 777久久人妻少妇嫩草av网站| 国产国语露脸激情在线看| 极品人妻少妇av视频| 婷婷丁香在线五月| √禁漫天堂资源中文www| 亚洲情色 制服丝袜| 色播在线永久视频| 丰满少妇做爰视频| 老司机在亚洲福利影院| 国产日韩欧美视频二区| 成人av一区二区三区在线看 | 这个男人来自地球电影免费观看| 久9热在线精品视频| 亚洲精品日韩在线中文字幕| 久久久久国内视频| 麻豆国产av国片精品| 国产精品亚洲av一区麻豆| 大片电影免费在线观看免费| 久久人人爽av亚洲精品天堂| 啦啦啦视频在线资源免费观看| 国产成+人综合+亚洲专区| 亚洲一区中文字幕在线| 少妇 在线观看| 三上悠亚av全集在线观看| videos熟女内射| 男人添女人高潮全过程视频| 欧美激情 高清一区二区三区| 国产免费av片在线观看野外av| 宅男免费午夜| 欧美激情极品国产一区二区三区| 亚洲av日韩精品久久久久久密| 狠狠狠狠99中文字幕| 亚洲av国产av综合av卡| 亚洲伊人久久精品综合| 亚洲伊人久久精品综合| 成人国产av品久久久| 桃红色精品国产亚洲av| 热99re8久久精品国产| 777米奇影视久久| 日本av手机在线免费观看| 99国产精品一区二区三区| 国产精品一区二区精品视频观看| 18禁观看日本| 亚洲欧美色中文字幕在线| 久久亚洲国产成人精品v| 亚洲人成77777在线视频| 天天操日日干夜夜撸| 免费观看a级毛片全部| 性少妇av在线| 狂野欧美激情性bbbbbb| 亚洲va日本ⅴa欧美va伊人久久 | 精品人妻熟女毛片av久久网站| 叶爱在线成人免费视频播放| a级毛片在线看网站| 18禁黄网站禁片午夜丰满| 国产亚洲午夜精品一区二区久久| 一级毛片精品| 99国产精品一区二区三区| 丰满饥渴人妻一区二区三| 亚洲成人免费av在线播放| 考比视频在线观看| 日韩欧美免费精品| 伊人亚洲综合成人网| bbb黄色大片| 亚洲精品国产色婷婷电影| 午夜福利视频精品| 亚洲一区中文字幕在线| 99久久人妻综合| 久久天堂一区二区三区四区| 亚洲精品国产精品久久久不卡| 国产精品久久久久久人妻精品电影 | 91精品国产国语对白视频| 九色亚洲精品在线播放| 精品人妻1区二区| 久久精品国产亚洲av香蕉五月 | 岛国在线观看网站| 欧美成人午夜精品| 女人高潮潮喷娇喘18禁视频| 热99久久久久精品小说推荐| 国产一区二区三区在线臀色熟女 | 国产欧美日韩综合在线一区二区| 2018国产大陆天天弄谢| 啦啦啦中文免费视频观看日本| 9热在线视频观看99| 人人妻,人人澡人人爽秒播| 亚洲成人国产一区在线观看| 狠狠狠狠99中文字幕| 免费人妻精品一区二区三区视频| 国产97色在线日韩免费| 色视频在线一区二区三区| 五月天丁香电影| 伊人久久大香线蕉亚洲五| 国产91精品成人一区二区三区 | 男女之事视频高清在线观看| 午夜老司机福利片| 亚洲av电影在线观看一区二区三区| 国产成+人综合+亚洲专区| 老司机靠b影院| 考比视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩视频精品一区| 国产成人啪精品午夜网站| 免费久久久久久久精品成人欧美视频| 中文字幕色久视频| 爱豆传媒免费全集在线观看| 大片电影免费在线观看免费| 精品少妇内射三级| 老熟妇乱子伦视频在线观看 | 女性生殖器流出的白浆| 嫁个100分男人电影在线观看| 少妇的丰满在线观看| 久久国产亚洲av麻豆专区| 久热爱精品视频在线9| 十八禁高潮呻吟视频| 人人妻人人爽人人添夜夜欢视频| 欧美在线一区亚洲| 午夜老司机福利片| 亚洲欧美日韩高清在线视频 | 国产视频一区二区在线看| 在线观看一区二区三区激情| 久久久久精品人妻al黑| 亚洲av国产av综合av卡| 久久久久视频综合| 国产1区2区3区精品| 亚洲精品一卡2卡三卡4卡5卡 | 久久人妻福利社区极品人妻图片| a在线观看视频网站| 亚洲国产中文字幕在线视频| 王馨瑶露胸无遮挡在线观看| 亚洲国产欧美日韩在线播放| 日韩视频在线欧美| 一个人免费在线观看的高清视频 | 99久久精品国产亚洲精品| 国产高清国产精品国产三级| 老司机亚洲免费影院| 好男人电影高清在线观看| 国产精品二区激情视频| 亚洲成av片中文字幕在线观看| 满18在线观看网站| 日本vs欧美在线观看视频| 亚洲欧洲精品一区二区精品久久久| 又黄又粗又硬又大视频| 欧美日韩黄片免| 少妇被粗大的猛进出69影院| 满18在线观看网站| 国产一区二区在线观看av| 久久久久久久精品精品| 国产精品一区二区精品视频观看| 波多野结衣av一区二区av| 国产一区二区三区在线臀色熟女 | 搡老乐熟女国产| 精品一品国产午夜福利视频| 考比视频在线观看| 男女下面插进去视频免费观看| 国产成人免费观看mmmm| 久久国产亚洲av麻豆专区| 九色亚洲精品在线播放| 久久狼人影院| 日韩视频一区二区在线观看| 精品熟女少妇八av免费久了| 国产日韩一区二区三区精品不卡| 老司机影院成人| av线在线观看网站| 捣出白浆h1v1| 黄片大片在线免费观看| 777米奇影视久久| 我要看黄色一级片免费的| 在线观看人妻少妇| 性高湖久久久久久久久免费观看| 桃红色精品国产亚洲av| 国产高清videossex| 成年人午夜在线观看视频| 精品人妻1区二区| 国产免费现黄频在线看| 国产精品偷伦视频观看了| 法律面前人人平等表现在哪些方面 | 久久久久久久国产电影| 国产欧美日韩精品亚洲av| 亚洲精品日韩在线中文字幕| 亚洲欧洲日产国产| 欧美 亚洲 国产 日韩一| 黑人巨大精品欧美一区二区蜜桃| 一二三四社区在线视频社区8| 国产成人影院久久av| 国产片内射在线| 又紧又爽又黄一区二区| 久久亚洲精品不卡| cao死你这个sao货| 夫妻午夜视频| 手机成人av网站| 精品视频人人做人人爽| 亚洲精品中文字幕一二三四区 | 国产亚洲欧美精品永久| 大片免费播放器 马上看| 纯流量卡能插随身wifi吗| 国产91精品成人一区二区三区 | 久9热在线精品视频| 久久青草综合色| 自线自在国产av| 欧美日韩中文字幕国产精品一区二区三区 | 夫妻午夜视频| 国产精品免费视频内射| 考比视频在线观看| 国产精品久久久久久人妻精品电影 | 女性生殖器流出的白浆| 99热国产这里只有精品6| 50天的宝宝边吃奶边哭怎么回事| 丰满饥渴人妻一区二区三| 国产欧美日韩精品亚洲av| 69精品国产乱码久久久| 免费在线观看日本一区| 国产亚洲一区二区精品| 天天躁狠狠躁夜夜躁狠狠躁| 91成年电影在线观看| 日韩一卡2卡3卡4卡2021年| 一本久久精品| 日韩视频在线欧美| 日本a在线网址| 婷婷丁香在线五月| 黄色视频在线播放观看不卡| 久热这里只有精品99| av网站免费在线观看视频| 美女扒开内裤让男人捅视频| 精品福利观看| 另类亚洲欧美激情| 国产黄色免费在线视频| 国产人伦9x9x在线观看| 亚洲欧美激情在线| 精品少妇黑人巨大在线播放| 一区在线观看完整版| 国产99久久九九免费精品| a 毛片基地| 国产免费av片在线观看野外av| 又大又爽又粗| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产成人一精品久久久| 成人免费观看视频高清| 少妇人妻久久综合中文| 99国产精品一区二区三区| 国产精品麻豆人妻色哟哟久久| 亚洲精品中文字幕一二三四区 | 美女扒开内裤让男人捅视频| 一级a爱视频在线免费观看| 午夜福利一区二区在线看| 亚洲美女黄色视频免费看| av国产精品久久久久影院| 成年av动漫网址| 久久久久精品人妻al黑| 桃红色精品国产亚洲av| 精品少妇久久久久久888优播| 成人av一区二区三区在线看 | 丝袜美足系列| 国产日韩欧美在线精品| 色精品久久人妻99蜜桃| 精品久久蜜臀av无| 精品少妇一区二区三区视频日本电影| 国产精品九九99| 欧美日韩视频精品一区| 久久久久精品国产欧美久久久 | 国产xxxxx性猛交| 亚洲黑人精品在线| 中国国产av一级| 一边摸一边做爽爽视频免费| 久久亚洲国产成人精品v| 国产亚洲精品久久久久5区| 亚洲精品国产一区二区精华液| 午夜免费鲁丝| 午夜久久久在线观看| 久久香蕉激情| 啦啦啦 在线观看视频| 午夜视频精品福利| 建设人人有责人人尽责人人享有的| 男女床上黄色一级片免费看| 最新的欧美精品一区二区| 9色porny在线观看| 国产高清视频在线播放一区 | 黄色a级毛片大全视频| 伊人久久大香线蕉亚洲五| 天天添夜夜摸| av片东京热男人的天堂| 亚洲一区中文字幕在线| 在线永久观看黄色视频| 国产色视频综合| 日日夜夜操网爽| 日本五十路高清| 满18在线观看网站| 精品一区在线观看国产| 动漫黄色视频在线观看| 久久av网站| 免费观看a级毛片全部| 欧美乱码精品一区二区三区| 国产xxxxx性猛交| 亚洲av国产av综合av卡| 50天的宝宝边吃奶边哭怎么回事| 成人国产一区最新在线观看| 亚洲人成77777在线视频| av片东京热男人的天堂| 无限看片的www在线观看| 另类亚洲欧美激情| 69精品国产乱码久久久| 王馨瑶露胸无遮挡在线观看| 亚洲国产毛片av蜜桃av| 丝袜人妻中文字幕| 9色porny在线观看| 老熟女久久久| 少妇裸体淫交视频免费看高清 | av线在线观看网站| 精品免费久久久久久久清纯 | 欧美av亚洲av综合av国产av| 69av精品久久久久久 | 欧美另类亚洲清纯唯美| 一级黄色大片毛片| 啦啦啦 在线观看视频| 首页视频小说图片口味搜索| 国产精品一区二区精品视频观看| 亚洲精品一二三| 免费观看a级毛片全部| 丰满少妇做爰视频| 午夜福利一区二区在线看| 成人黄色视频免费在线看| 精品高清国产在线一区| 亚洲国产毛片av蜜桃av| 欧美在线一区亚洲| 久久综合国产亚洲精品| 午夜精品国产一区二区电影| 中文字幕另类日韩欧美亚洲嫩草| 首页视频小说图片口味搜索| 狂野欧美激情性bbbbbb| 十八禁网站网址无遮挡| 大码成人一级视频| 秋霞在线观看毛片| 亚洲国产日韩一区二区| 可以免费在线观看a视频的电影网站| 777久久人妻少妇嫩草av网站| 亚洲第一青青草原| 99热网站在线观看| 亚洲一区中文字幕在线| 亚洲,欧美精品.| 一本色道久久久久久精品综合| 精品福利观看| 欧美黄色片欧美黄色片| 日韩制服骚丝袜av| 中文字幕精品免费在线观看视频| 久久亚洲精品不卡| 99精品久久久久人妻精品| 久久久久网色| 日韩三级视频一区二区三区| 桃花免费在线播放| 久热爱精品视频在线9| 久久女婷五月综合色啪小说| 一本—道久久a久久精品蜜桃钙片| av免费在线观看网站| 成人亚洲精品一区在线观看| 国产精品欧美亚洲77777| 国产亚洲av高清不卡| 一个人免费在线观看的高清视频 | 亚洲专区字幕在线| 在线永久观看黄色视频| 免费观看a级毛片全部| 免费日韩欧美在线观看| 精品一品国产午夜福利视频| 男人操女人黄网站| 国产在视频线精品| 一二三四社区在线视频社区8| 欧美一级毛片孕妇| 欧美精品亚洲一区二区| 岛国毛片在线播放| 久久久精品94久久精品| 亚洲专区中文字幕在线| 最近最新中文字幕大全免费视频| 91大片在线观看| 免费黄频网站在线观看国产| 男女边摸边吃奶| 国产成人a∨麻豆精品| 精品人妻在线不人妻| 18禁黄网站禁片午夜丰满| 亚洲精品成人av观看孕妇| 五月开心婷婷网| 中文字幕最新亚洲高清| 免费黄频网站在线观看国产| 69av精品久久久久久 | 欧美精品av麻豆av| 午夜精品国产一区二区电影| 精品福利永久在线观看| 久久 成人 亚洲| 嫩草影视91久久| 麻豆av在线久日| 欧美乱码精品一区二区三区| 国产不卡av网站在线观看| 亚洲伊人久久精品综合| 久久人人爽人人片av| 视频区欧美日本亚洲| 成人影院久久| 天天躁狠狠躁夜夜躁狠狠躁| tocl精华| 一边摸一边抽搐一进一出视频| 精品国产超薄肉色丝袜足j| 嫁个100分男人电影在线观看| 又黄又粗又硬又大视频| 久久av网站| 亚洲免费av在线视频| 悠悠久久av| 最新的欧美精品一区二区| 久久久久网色| videosex国产| 欧美另类亚洲清纯唯美| 久久久久久人人人人人| 满18在线观看网站| 中文字幕高清在线视频| 91麻豆精品激情在线观看国产 | 中文精品一卡2卡3卡4更新| 久久久久久久大尺度免费视频| 久久免费观看电影| 国产在视频线精品| 精品少妇黑人巨大在线播放| 国产男人的电影天堂91| 男人舔女人的私密视频| 黄片大片在线免费观看| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 国产高清videossex| 欧美精品高潮呻吟av久久| 国产熟女午夜一区二区三区| 国产精品99久久99久久久不卡| 香蕉丝袜av| 91精品伊人久久大香线蕉| 午夜福利,免费看| www.精华液| 91成人精品电影| 中文字幕人妻熟女乱码| 精品少妇一区二区三区视频日本电影| 亚洲七黄色美女视频| 久久国产精品人妻蜜桃| 久久国产精品男人的天堂亚洲| 国产成人av教育| 久久久精品区二区三区| 美女大奶头黄色视频| 老司机深夜福利视频在线观看 | 亚洲精品久久午夜乱码| 午夜福利乱码中文字幕| av片东京热男人的天堂| 成人av一区二区三区在线看 | 亚洲精品一卡2卡三卡4卡5卡 | 在线av久久热| 成年av动漫网址| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久成人av| 19禁男女啪啪无遮挡网站| 精品一区在线观看国产| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆| 亚洲国产成人一精品久久久| 国产黄色免费在线视频| 黄色视频在线播放观看不卡| 中文字幕人妻丝袜制服| 精品卡一卡二卡四卡免费| 黄色视频不卡| 青青草视频在线视频观看| 国产精品国产av在线观看| 久久久久国产一级毛片高清牌| 18禁观看日本| 久久香蕉激情| 午夜两性在线视频| 亚洲一区二区三区欧美精品| 不卡av一区二区三区| 亚洲性夜色夜夜综合| 亚洲成人手机| 亚洲精品国产精品久久久不卡| 亚洲精品久久午夜乱码| 国产一卡二卡三卡精品| 欧美中文综合在线视频| 欧美97在线视频| av又黄又爽大尺度在线免费看| 男女午夜视频在线观看| 国产精品秋霞免费鲁丝片| 午夜免费成人在线视频| 国产成人一区二区三区免费视频网站| 在线亚洲精品国产二区图片欧美| 在线观看免费日韩欧美大片| 99热全是精品| av超薄肉色丝袜交足视频| 丝袜脚勾引网站| 51午夜福利影视在线观看| av电影中文网址| 久久中文字幕一级| 欧美乱码精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线免费观看视频4| 精品少妇黑人巨大在线播放| 人人澡人人妻人| 80岁老熟妇乱子伦牲交| svipshipincom国产片| 精品视频人人做人人爽| 日韩中文字幕视频在线看片| 在线观看一区二区三区激情| 少妇人妻久久综合中文| 99国产极品粉嫩在线观看| 在线观看免费午夜福利视频| 国产高清视频在线播放一区 | 午夜两性在线视频| 免费高清在线观看视频在线观看| 夜夜骑夜夜射夜夜干| 香蕉丝袜av| 欧美日韩国产mv在线观看视频| 日韩欧美免费精品| 热re99久久国产66热| 亚洲九九香蕉| 中文字幕精品免费在线观看视频| 国产免费现黄频在线看| 99久久综合免费| 免费高清在线观看日韩| 亚洲av成人一区二区三| 国产免费现黄频在线看| 窝窝影院91人妻| 国产欧美亚洲国产| 久久午夜综合久久蜜桃| 久久亚洲国产成人精品v| 精品第一国产精品| 99re6热这里在线精品视频| 黑人欧美特级aaaaaa片| 午夜福利,免费看| 久久国产精品人妻蜜桃| 黄片小视频在线播放| 一级片'在线观看视频| 国产在线观看jvid| 91麻豆av在线| 啦啦啦免费观看视频1| 亚洲av男天堂| 亚洲成人免费av在线播放| 少妇的丰满在线观看| 少妇精品久久久久久久| 色综合欧美亚洲国产小说| 久久亚洲精品不卡| 一级毛片电影观看| 黄片小视频在线播放| 我的亚洲天堂| 欧美日韩中文字幕国产精品一区二区三区 | 欧美性长视频在线观看| 亚洲欧美精品自产自拍| 亚洲自偷自拍图片 自拍| 日韩欧美免费精品| 欧美精品啪啪一区二区三区 | 日本猛色少妇xxxxx猛交久久| 精品视频人人做人人爽| 美女扒开内裤让男人捅视频| 捣出白浆h1v1| 亚洲精品久久午夜乱码| 91精品三级在线观看| 国产免费视频播放在线视频| 岛国毛片在线播放| 国产亚洲精品一区二区www | 97人妻天天添夜夜摸| 狠狠婷婷综合久久久久久88av| netflix在线观看网站| 满18在线观看网站| 一边摸一边抽搐一进一出视频| 日韩人妻精品一区2区三区| 老司机在亚洲福利影院| 欧美亚洲 丝袜 人妻 在线| 午夜两性在线视频| 欧美日韩一级在线毛片| 丝袜在线中文字幕| 精品国产一区二区久久| 久久久国产欧美日韩av| 午夜久久久在线观看| 国产不卡av网站在线观看| 一级毛片女人18水好多| 久久精品人人爽人人爽视色| 一本久久精品| 中亚洲国语对白在线视频| 建设人人有责人人尽责人人享有的| 国产av又大| www.自偷自拍.com| 18禁裸乳无遮挡动漫免费视频| 男人添女人高潮全过程视频|