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

    用曲波變換消除VSP數(shù)據(jù)處理中的空間假頻與采樣不規(guī)則

    2016-02-05 08:55:15彭軍韓立國周家雄王宇隋波李海鵬
    地質(zhì)與勘探 2016年3期
    關(guān)鍵詞:曲波同相軸插值

    彭軍,韓立國,周家雄,王宇,隋波,李海鵬

    (1.中海石油(中國)有限公司湛江分公司,廣東湛江 524057;2.吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院,吉林長春 130026)

    用曲波變換消除VSP數(shù)據(jù)處理中的空間假頻與采樣不規(guī)則

    彭 軍1,2,韓立國2,周家雄1,王 宇1,隋 波1,李海鵬1

    (1.中海石油(中國)有限公司湛江分公司,廣東湛江 524057;2.吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院,吉林長春 130026)

    由于條件與經(jīng)費的限制,在實際地震資料采集過程中,常常會遇到獲得的地震資料較稀疏或不規(guī)則,導(dǎo)致空間假頻、采樣不規(guī)則的問題,影響地震資料的成像效果??衫玫卣鸬揽臻g插值來解決這些問題,提升地震數(shù)據(jù)處理結(jié)果品質(zhì)。曲波變換的基函數(shù)是具有傾角、頻率、時間三方面同時局部化的函數(shù),它可以對空間信號做稀疏表達,適合地震資料的插值處理。本文采用該方法對理論模型數(shù)據(jù)和實際VSP數(shù)據(jù)作插值處理,解決不規(guī)則采樣、空間假頻問題。實際應(yīng)用表明該方法精度高,插值準(zhǔn)確,處理結(jié)果的振幅和同相軸連續(xù)性好。

    曲波變換 空間假頻 VSP 插值

    Peng Jun,Han Li-guo,Zhou Jia-xiong,Wang Yu,Sui Bo,Li Hai-peng. Removing spatial alias frequency and irregularity of sampling in VSP data processing by using curvelet transform[J].Geology and Exploration,2016,52(3):0537-0542.

    1 引言

    隨著油氣勘探開發(fā)不斷深入,油氣勘探、開發(fā)由構(gòu)造油藏轉(zhuǎn)向巖性油藏,由厚層油氣藏向薄層油氣藏邁進。在復(fù)雜的地質(zhì)條件下,提高鉆探成功率,鉆前地質(zhì)研究與認識對邊際油氣藏的勘探、開發(fā)顯得尤為重要。地震資料作為認識地下地質(zhì)情況的重要手段,地震資料的好壞,將極大的限制地質(zhì)人員對地下情況認知的質(zhì)量。良好的地震資料品質(zhì)將極大提高地質(zhì)人員對地下地質(zhì)情況認識的精度。無論海洋地震勘探、陸地地震勘探,還是VSP(垂直地震剖面)采集,由于資金和采集技術(shù)發(fā)展水平的限制,野外地震采集往往存在假頻和不規(guī)則采樣的問題,特別是海洋三維地震的聯(lián)絡(luò)線方向。海上油氣勘探、開發(fā)相對陸地油氣勘探、開發(fā)成本更高,對鉆前地下地質(zhì)條件認識要求更高,這就需要擁有質(zhì)量品質(zhì)優(yōu)良的地震資料為地質(zhì)研究人員提供較精細的地質(zhì)認知條件。VSP數(shù)據(jù)為時間與深度記錄格式,相比于傳統(tǒng)地震資料具有分辨率高、具有深度信息等優(yōu)勢,使VSP數(shù)據(jù)處理近年來受到極大關(guān)注。

    地震資料反射波勘探是對中深層地層識別的有效方法,廣泛應(yīng)用于油氣與金屬礦床成像(徐明才等,1999;梁光河等,2001;孫明等,2001;徐明才等,2005;李戰(zhàn)業(yè)等,2009),然而地震資料處理過程中的很多算法(如直達波分離處理、去多次波處理(SRME)、偏移成像處理等)都需要規(guī)則且密集的地震采樣,就需采用恰當(dāng)?shù)姆椒▽⒉灰?guī)則采樣處理為規(guī)則采樣。在VSP采集中由于成本控制,一般在目的層段采集樣點較為密集,在非目的層段采集樣點較為稀疏。由于采樣的不規(guī)則和空間采樣較稀疏將產(chǎn)生空間假頻??臻g假頻的存在會導(dǎo)致頻散現(xiàn)象,無法得到較好的地震資料成像,使解釋人員無法正確認識地下地質(zhì)情況。目前,地震道插值和面元均化等處理技術(shù),都在努力解決空間假頻和不規(guī)則采樣的問題,采用這些方法處理后的數(shù)據(jù)含有的地球物理信息更加真實地反映地下地質(zhì)體,更有利于進行構(gòu)造解釋和地震地層學(xué)的研究。所以地震資料去假頻和規(guī)則化是地球物理研究中的一個難點與熱點問題,近20多來一直受到專業(yè)人員的關(guān)注和重視。

    常見的地震資料插值技術(shù)主要有:T-X預(yù)測誤差濾波插值方法,抗假頻F-X域地震道插值方法、f-k域插值(李學(xué)聰?shù)龋?009)、Radon域插值等。這些傳統(tǒng)的插值方法都存在它們各自的問題:有的是計算量太大,在實際地震資料處理中可用性大大降低;有的是插值系數(shù)不能為任意數(shù);還有的是不能正確內(nèi)插空間假頻地震道。我們希望插值技術(shù)不僅可以解決不規(guī)則采樣問題,同時也解決空間假頻問題。針對以前算法的不足,最近也發(fā)展了很多新的算法(高建軍等,2011,劉喜武等,2011),希望尋找到一種既準(zhǔn)確又易于快速實現(xiàn)的計算方法。

    在1999年,Candesetal.(1999)提出了曲波變換(curvelet transform),構(gòu)建了曲波變換的數(shù)學(xué)理論基礎(chǔ)。2006年,Candes發(fā)布了曲波變換的快速算法(第二代曲波變換數(shù)值算法)(Emmanueletal.,2005),計算效率大大提高。2007年,Demanet和Ying發(fā)表了鏡像延伸的曲波變換算法。最近一些年,曲波變換不僅理論上得到了快速的發(fā)展,并在工業(yè)計算中發(fā)揮了重要作用,廣泛的應(yīng)用在數(shù)字信號處理方面。曲波變換具有時頻局部化、傾角局部化的性質(zhì);地震資料是分段連續(xù)、具有不同傾角的地震同相軸組成,地震資料也具有時傾角局部特性,因此將曲波變換能夠稀疏的表達地震資料,更加適合地震資料去噪、插值、多次波去除等方面的處理研究(Neelamani,Retal.,2008;彭才等,2008;Kumar Vetal.,2009;仝中飛等,2009;Mostafa Naghizadehetal.,2010;薛念,2010;周家雄等,2010;劉磊等,2011; Jing Wangetal.,2013)。

    2 基本原理

    第一代曲波變換需要進行子帶分解、平滑分塊、正則化和Ridgelet分解實現(xiàn)曲波變換,這就導(dǎo)致了巨大的數(shù)據(jù)冗余,計算效率低下。第二代曲波變換脫離了Ridgelet理論,計算速度和計算冗余都得到了改善。第二代曲波變換的計算復(fù)雜程度滿足O(n2logn),它的逆變換不僅準(zhǔn)確,并且逆變換的計算復(fù)雜程度也滿足O(n2logn)??焖匐x散曲波變換算法的實現(xiàn)有三種算法:一種是基于不規(guī)則采樣的快速傅立葉變換設(shè)計的算法,第二種是基于外推算法的快速曲波變換,還有一種是最新的基于鏡像外推的快速算法(周家雄等,2010)。

    (1)

    對于每一個j≥j0,可以用傅立葉變換定義頻率窗Uj,即

    (2)

    其中?j/2?是j/2的整數(shù)部分。所以,Uj是一個帶方向的楔狀體,它等于半徑窗函數(shù)和角度窗函數(shù)的乘積,且具有對稱性,Uj(r,θ)+Uj(r,θ+π),如圖1所示。

    圖1 頻域和時域曲波Fig.1 Curvelets in frequency domain and time domain a-頻域;b-時域a-frequency domain;b-time domain

    (3)

    其他2-j尺度的曲波函數(shù)都可以由該母函數(shù)φj(x)經(jīng)過旋轉(zhuǎn)、平移得到。

    (4)

    其中Rθ是旋轉(zhuǎn)函數(shù),θ是旋轉(zhuǎn)角,

    (5)

    把上式帶入到頻率域中,可以得到

    (6)

    從曲波變換的原理可知曲波變換具有多分辨率分析以及時頻局部性的特點,更重要的具有局部方向特性。在進行地震插值時,將缺失道標(biāo)記為0,非缺失道標(biāo)記為1,這樣在計算時只針對缺失道進行數(shù)據(jù)填充,已有數(shù)據(jù)不被改變。然后對地震數(shù)據(jù)進行曲波變換,得到曲波系數(shù)。因為地震數(shù)據(jù)在曲波域是稀疏的,在曲波局部方向特性的控制下,缺失的地震數(shù)據(jù)可以通過相近地震道曲波變換系數(shù)得到相應(yīng)的曲波系數(shù)。將得到的曲波系數(shù)經(jīng)過曲波逆變換,就可以將缺失的地震數(shù)據(jù)填補上,達到插值的目的。

    地震資料空間假頻是由于地震資料采樣不規(guī)則或采樣太稀疏導(dǎo)致的,因此當(dāng)?shù)卣饠?shù)據(jù)空間采樣規(guī)則、密度足夠時,就能有效去除空間假頻。

    3 模型實驗

    地震資料的空間假頻與地震同相軸的傾角關(guān)系密切。為了避免空間假頻的出現(xiàn),大傾角同相軸需要更小的空間采樣間隔。本文采用40赫茲的雷克子波所組成的不同傾角的同相軸,共150道,如圖 2a所示。模型數(shù)據(jù)的同相軸都是向右方向傾斜的,那么它的f-k譜的波數(shù)應(yīng)該是正值。但是由于空間采樣不足,大傾角的同相軸出現(xiàn)了假頻的現(xiàn)象,模型數(shù)據(jù)的f-k譜的部分能量被劃入了負象限。

    為了消除空間假頻現(xiàn)象,采用曲波進行均勻采樣數(shù)據(jù)內(nèi)插。通過1 ∶1的等間距插值處理,成功的解決了假頻問題,處理后的數(shù)據(jù)的f-k譜的能量都恢復(fù)到了正象限,如圖2b所示。

    圖2 去假頻處理實驗Fig.2 Experiment of removing alias a-存在假頻的地震數(shù)據(jù)及其頻譜;b-插值后的地震數(shù)據(jù)及其頻譜a-eismic data and its spectrum with alias;b-interpolated seismic data and its spectrum

    設(shè)計如圖3a所示模型,研究在地震道缺失的情況下,該算法同時處理空間假頻和不規(guī)則采樣的能力。該數(shù)據(jù)不僅缺失道,同時也存在假頻的影響,見圖 3a。圖 3b是插值處理后結(jié)果及其f-k譜。該結(jié)果表明,本文方法不僅可以消除假頻,而且可以把缺失的數(shù)據(jù)補齊,達到均勻采樣的目的。

    圖3 缺失數(shù)據(jù)的模型及其數(shù)據(jù)恢復(fù)處理結(jié)果Fig.3 Missing data model and data recovery processing results a-不規(guī)則地震數(shù)據(jù)及其頻譜;b-插值后的地震數(shù)據(jù)及其頻譜a-irregular seismic data and its spectrum;b-interpolated seismic data and its spectrum

    4 VSP數(shù)據(jù)處理

    為了節(jié)省成本與時間,在非目的層層段,VSP野外采集的空間采樣間隔都比較大。本次VSP采集用的是Sercel GeoWaves2級3分量檢波器,2級檢波點距離10m,在非目的層段炮檢距40m,在目的層段炮檢距20m。所以,非目的層的數(shù)據(jù)缺值達到了50%,同時這批數(shù)據(jù)也是不規(guī)則采樣。本次插值運算的目的是補齊非目的層段缺失的數(shù)據(jù)。如圖 4a所示,由于不規(guī)則采樣,地震數(shù)據(jù)的f-k譜出現(xiàn)了類似假頻的干擾。通過曲波變換插值處理,空間假頻消除了,f-k譜更加符合理論,同相軸和振幅連續(xù)性好,如圖4b所示。

    5 結(jié)論

    地震資料作為認知地下地質(zhì)情況的基礎(chǔ)資料,其品質(zhì)的好壞直接決定了地質(zhì)人員對地下地質(zhì)情況認知的精度。地震資料的規(guī)則化與去假頻是地震資料處理中的一個重要步驟。本文采用曲波變換進行VSP資料的規(guī)則化與去假頻處理,通過理論模型實驗表明,曲波變換能進行等間隔與非等間隔插值,采用該方法能有效去除假頻,保證插值后的地震數(shù)據(jù)相位的精度。插值后VSP數(shù)據(jù)同相軸、振幅連續(xù),頻譜圖清晰的顯示插值前后僅有假頻被去除,有效頻率未有變化。將曲波變換插值應(yīng)用于實際VSP數(shù)據(jù)表明,使用曲波變換的插值結(jié)果能有效去除空間假頻,VSP數(shù)據(jù)信噪比高,振幅的局部特性保持較好。采用曲波變換對模型與實際VSP數(shù)據(jù)插值結(jié)果表明該方法計算結(jié)果準(zhǔn)確,具有良好的實用性。曲波變換的高精度、振幅屬性的局部保持特性對于高保真地震資料的處理將帶來極大的幫助,有利于后續(xù)儲層預(yù)測的順利進行。

    圖4 插值處理前后VSP地震數(shù)據(jù)對比及頻譜的變化Fig.4 Comparison of VSP seismic and spectra before and after interpolation and their variations a-處理前VSP地震數(shù)據(jù)與頻譜;b-處理后VSP地震數(shù)據(jù)與頻譜a-before processing;b-after processing

    本文僅將該方法應(yīng)用于VSP數(shù)據(jù)處理,下一步我們研究把該方法應(yīng)用在海洋拖纜數(shù)據(jù)中。研究當(dāng)海纜采集地震資料連續(xù)缺失大塊數(shù)據(jù)時,如何使用曲波進行插值處理,重構(gòu)該區(qū)地震資料。

    Candes E J,Donoho D L. 1999. Curvelets: a surprisingly effective nonadaptive representation of objects with edges[C]. Curves and Surfaces,Vanderbilt University Press,Nashville,TN.: 1-16

    EmmanuelCandès,Laurent Demanet,David Donoho,Lexing Ying. 2005.Fast Discrete Curvelet Transforms[J]. Siam Journal on Multiscale Modeling & Simulation,5:861-899

    Gao Jian-jun,Chen Xiao-hong,Li Jing-ye. 2011. The reconstruction method for 3D irregular seismic data[J],Oil geophysical prospecting,46(1):40-47 (in Chinese with English abstract)

    Jing Wang,Xuekai Sun,Sam Zandong Sun. 2013. Seismic denoising in curvelet domain: an adaptive thresholding based on local variance analysis[C]. SEG Houston 2013 Annual Meeting: 4309-4313

    Kumar V.,Herrmann,F(xiàn). 2009. Incoherent noise suppression with curvelet-domain sparsity[C]. Technical Program Expanded Abstracts:3356-3360

    Li Xue-cong,Liu Yi-ke,Chang Xu,Jin De-gang,Zhang Heng-chao.2009. Interpolation for seismic traces based on F-K migration and demigration [J]. Progress in Geophysics,24(6):2020-2029 (in Chinese with English abstract)

    Li Zhan-ye,Yin Jun-jie,Wang Yun. 2009. The application of the numerical simulation and imaging of seismic scattered wave in exploration of metal mineral deposits [J]. Geology and Exploration,45(2):80-84 (in Chinese with English abstract)

    Liang Guang-he,Cai Xin-ping,Zhang Bao-lin,Xu Xing-wang. 2001. The application of seismic exploration method in deep prediction of gold deposits [J]. Geology and Exploration,37(6):29-33 (in Chinese with English abstract)

    Liu Lei,Liu Zhen,Zhang Jun-hua. 2011. Denoising and detecting seismic weak signal based on curvelete thresholding method. Progress in geophys,26(4):1415-1422 (in Chinese with English abstract)

    Liu Xi-wu,Liu Hong,Liu Bin. 2004. Astudy on algorithm for reconstruction of De-alias uneven seismic data[J]. Chinese Journal of geophysics,47(2) :299-305(in Chinese with English abstract)

    Mostafa Naghizadeh,Mauricio D. Sacchi. 2010. Beyond alias hierarchical scale curvelet interpolation of regularly and irregularly sampled seismic data. Geophysics,75(6):189-202

    Neelamani R,A. I. Baumstein,D. G. Gillard,M. T.Hadidi,W. L. Soroka. 2008. Coherent and random noise attenuation using the curvelet transform[J]. The Leading Edge,27(2):240-248

    Peng Cai,Chang Zhi,Zhu Shi-jun. 2008. Seismic data denoising by curvelet transformation[J]. Geophysical prospecting for petroleum,47(5):461-464(in Chinese with English abstract)

    Sun Ming,Lin Jun. 2001. Study of seismic scattering wave field numerical model for metallic ore exploration [J] . Geology and Exploration,37(4):68-70 (in Chinese with English abstract)

    Tong Zong-fei,Wang De-li,Liu Qi. 2009. The study on seismic data Interpolation with Curvelet thresholding iterative method[C].Chinese Geophysics:66(in Chinese with English abstract)

    Xu Ming-cai,Gao Jing-hua,Chai Ming-tao,Wang Guang-ke,Liu Jian-xun,Rong Li-xin. 1999. Seismic technique applied to active fault assessment [J]. Geology and Exploration,35(1):36-40 (in Chinese with English abstract)

    Xu Ming-cai,Gao Jing-hua,Rong Li-xin,Chai Ming-tao,Lv Qing-tian,Shi Da-nian. 2005. Combining exploration technique of the ground seismic tomography and the high-resolution seismic method[J]. Geology and Exploration,41(4): 83-87 (in Chinese with English abstract)

    Xue Nian.2010. Curvelet transformation for seismic data denoising and interpolating[D]. Xi’an Jiao Tong University master thesis:5-6(in Chinese with English abstract)

    Zhou Jia-xiog,Wang Yu,Guo Ai-hua,Sun Yue-cheng,Li Hai-peng.2010. A study on denoising poststack 3D seismic data by Curvelet transform [J]. China Offshore Oil And Gas,22(1):17-21(in Chinese with English abstract)

    [附中文參考文獻]

    高建軍,陳小宏,李景業(yè).2011.三維不規(guī)則地震數(shù)據(jù)重建方法[J]. 石油地球物理勘探,46(1):40-47

    李學(xué)聰,劉伊克,常 旭,金德剛,張恒超.2009.基于F-K偏移和反偏移的地震道插值方法研究[J]. 地球物理學(xué)進展,24(6):2020-2029

    李戰(zhàn)業(yè),尹軍杰,王 赟.2009.地震散射波模擬成像在金屬礦勘探中的應(yīng)用[J]. 地質(zhì)與勘探,45(2):80-84

    梁光河,蔡新平,張寶林,徐興旺.2001.淺層地震勘探方法在金礦深部預(yù)測中的應(yīng)用[J]. 地質(zhì)與勘探,37(6):29-33

    劉 磊,劉 振,張軍華.2011. 曲波閾值法地震弱信號識別及去噪方法研究[J].地球物理學(xué)進展,26(4):1415-1422

    劉喜武,劉 洪,劉 彬.2004.反假頻非均勻地震數(shù)據(jù)重建方法研究[J]. 地球物理學(xué)報,47(2):299-305

    彭 才,常 智,朱仕軍.2008.基于曲波變換的地震數(shù)據(jù)去噪方法[J]. 石油物探,47(5):461-464

    孫 明,林 君.2001.金屬礦地震散射波場的數(shù)值模擬研究[J]. 地質(zhì)與勘探,37(4):68-70

    仝中飛,王德利,劉 琦.2009.Curvelet閾值迭代法在地震數(shù)據(jù)插值中的方法研究[C]. 中國地球物理:66

    徐明才,高景華,柴明濤,王廣科,劉建勛,榮立新.1999.應(yīng)用于評價活動斷層的地震方法技術(shù)[J]. 地質(zhì)與勘探,35(1):36-40

    徐明才,高景華,榮立新,柴銘濤,呂慶田,史大年.2005.地面地震層析成像和高分辨率地震聯(lián)合勘探技術(shù)[J]. 地質(zhì)與勘探,41(4):83-87

    薛 念. 2010. 基于Curvelet變換的地震數(shù)據(jù)插值和去噪[D]. 西安交通大學(xué)碩士論文:5-6

    周家雄,王 宇,郭愛華,孫月成,李海鵬. 2010. 基于三維曲波變換的疊后三維地震資料去噪處理研究[J]. 中國海上油氣,22(1):17-21

    Removing Spatial Alias Frequency and Irregularity of Sampling in VSP Data Processing by Using Curvelet Transform

    PENG Jun1,2,HAN Li-guo2,ZHOU Jia-xiong1,WANG Yu1,SUI Bo1,LI Hai-peng1

    (1.ChinaNationalOffshoreOilCompany,ZhanjiangBranch,Zhanjiang,Guangdong5240571; 2.CollegeofGeoExplorationScienceandTechnology,JilinUniversity,Changchun,Jilin130026)

    Spatial alias frequency and irregularity of sampling are often encountered in seismic acquisition due to limited conditions and funds,which can affect many aspects of seismic processing. Seismic interpolation can be used to solve these issues to improve the quality of seismic data processing. To do so,the curvelet transform decomposes the seismic data into local wavelets which have local time-frequency and local dip attributes. Spatial signal can be expressed sparsely by curvelet transform,so is suitable for seismic data interpolation. In this paper,we employ the curvelet transform to solve the problems of spatial alias frequency and irregular spatial sampling. The tests on model data and real VSP field data show that this method can achieve high precision,accurate interpolation and good continuity in amplitude and co-phase axes.

    curvelet transform,spatial aliasing,VSP,interpolation

    2015-09-15;

    2016-04-15;[責(zé)任編輯]陳偉軍。

    中海石油(中國)有限公司南海西部海域典型低孔低滲油氣藏勘探開發(fā)關(guān)鍵技術(shù)研究及實踐(編碼:CNOOC-KJ125ZDXM07LTD02ZJ11)贊助。

    彭 軍(1984年-),在讀碩士,工程師,現(xiàn)從事地球物理解釋與處理技術(shù)的研究和應(yīng)用工作。E-mail:pengjun@cnooc.com.cn。

    P315

    A

    0495-5331(2016)03-0537-06

    猜你喜歡
    曲波同相軸插值
    林海雪原(五)
    林海雪原(三)
    虛同相軸方法及其在陸上地震層間多次波壓制中的應(yīng)用
    林海雪原(四)
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    曲波變換三維地震數(shù)據(jù)去噪技術(shù)
    一種改進的相關(guān)法自動拾取同相軸
    一種反射同相軸自動拾取算法
    一種改進FFT多譜線插值諧波分析方法
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    亚洲成国产人片在线观看| 香蕉国产在线看| avwww免费| 欧美激情高清一区二区三区| 久久99热这里只有精品18| 色播在线永久视频| 中文字幕精品免费在线观看视频| 激情在线观看视频在线高清| 久久精品国产综合久久久| 亚洲免费av在线视频| АⅤ资源中文在线天堂| 午夜日韩欧美国产| 好男人在线观看高清免费视频 | 亚洲国产精品sss在线观看| 波多野结衣高清无吗| 男人舔奶头视频| 97人妻精品一区二区三区麻豆 | 亚洲狠狠婷婷综合久久图片| 一个人观看的视频www高清免费观看 | 国产精品久久久人人做人人爽| 亚洲第一电影网av| 久久亚洲真实| 国产av在哪里看| 欧美日韩瑟瑟在线播放| 级片在线观看| 最好的美女福利视频网| 亚洲avbb在线观看| 亚洲精品色激情综合| 啦啦啦观看免费观看视频高清| 老熟妇仑乱视频hdxx| 亚洲成人久久性| 亚洲最大成人中文| 九色国产91popny在线| 香蕉丝袜av| 亚洲av中文字字幕乱码综合 | 亚洲欧美精品综合一区二区三区| 国产精品久久视频播放| 欧美成人性av电影在线观看| 香蕉丝袜av| 天堂动漫精品| 色综合亚洲欧美另类图片| 国产亚洲精品av在线| 女人被狂操c到高潮| a级毛片a级免费在线| 国产精品,欧美在线| 亚洲七黄色美女视频| 很黄的视频免费| 高潮久久久久久久久久久不卡| 久久久久久久久中文| 中文字幕精品免费在线观看视频| 亚洲成人免费电影在线观看| 成人三级黄色视频| 久久狼人影院| 男人舔女人的私密视频| 日韩欧美一区视频在线观看| 国产一区在线观看成人免费| 天堂√8在线中文| 黄色女人牲交| 他把我摸到了高潮在线观看| 国产午夜精品久久久久久| 国产蜜桃级精品一区二区三区| 美女高潮到喷水免费观看| 久久中文字幕人妻熟女| 夜夜爽天天搞| 国产极品粉嫩免费观看在线| 成人亚洲精品av一区二区| 国产又黄又爽又无遮挡在线| 欧美一级a爱片免费观看看 | 91在线观看av| 美女国产高潮福利片在线看| 国产免费男女视频| 老司机靠b影院| 此物有八面人人有两片| 国产真人三级小视频在线观看| 香蕉av资源在线| 午夜免费观看网址| 国产精品自产拍在线观看55亚洲| 国产成人影院久久av| 国产真实乱freesex| 午夜老司机福利片| 国产人伦9x9x在线观看| 国产伦一二天堂av在线观看| 国产亚洲精品久久久久久毛片| 欧美亚洲日本最大视频资源| 欧美日本亚洲视频在线播放| 中文字幕人妻丝袜一区二区| 成人午夜高清在线视频 | 欧美av亚洲av综合av国产av| 一区二区日韩欧美中文字幕| 久久香蕉国产精品| 亚洲av第一区精品v没综合| 亚洲欧美一区二区三区黑人| 国产一区二区三区视频了| 在线观看日韩欧美| 午夜免费成人在线视频| 国产精品 欧美亚洲| 成人国产一区最新在线观看| 91九色精品人成在线观看| 桃色一区二区三区在线观看| 精品久久久久久久久久久久久 | 国产成人系列免费观看| 免费高清视频大片| 99久久无色码亚洲精品果冻| 女警被强在线播放| 国产在线精品亚洲第一网站| 岛国在线观看网站| 亚洲精品中文字幕在线视频| 婷婷六月久久综合丁香| 精品免费久久久久久久清纯| 在线观看免费午夜福利视频| 久久国产精品人妻蜜桃| 色老头精品视频在线观看| 亚洲精品国产区一区二| 丰满的人妻完整版| av视频在线观看入口| 亚洲色图 男人天堂 中文字幕| 精品久久蜜臀av无| 成人精品一区二区免费| 身体一侧抽搐| 男女下面进入的视频免费午夜 | 亚洲第一青青草原| 听说在线观看完整版免费高清| 18禁国产床啪视频网站| 女生性感内裤真人,穿戴方法视频| 十八禁网站免费在线| av欧美777| 女性生殖器流出的白浆| 国产亚洲av嫩草精品影院| 99热6这里只有精品| 久久久久久亚洲精品国产蜜桃av| 可以在线观看毛片的网站| 91字幕亚洲| 制服丝袜大香蕉在线| 久久久久久久精品吃奶| 亚洲人成网站高清观看| 丝袜在线中文字幕| 成在线人永久免费视频| 欧美中文综合在线视频| 欧美国产精品va在线观看不卡| 精品久久久久久,| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清| 18禁观看日本| 久久人妻福利社区极品人妻图片| 久久伊人香网站| 国产精品99久久99久久久不卡| av在线播放免费不卡| av超薄肉色丝袜交足视频| 日本在线视频免费播放| 日韩欧美 国产精品| 色综合站精品国产| 老司机午夜福利在线观看视频| 亚洲午夜精品一区,二区,三区| 男女做爰动态图高潮gif福利片| 久久久久国内视频| 成人欧美大片| 天天躁狠狠躁夜夜躁狠狠躁| 18美女黄网站色大片免费观看| 国产av一区在线观看免费| 黑人巨大精品欧美一区二区mp4| 亚洲激情在线av| 国产色视频综合| 脱女人内裤的视频| 一级a爱片免费观看的视频| 亚洲一码二码三码区别大吗| 成人一区二区视频在线观看| 人人妻人人澡欧美一区二区| a级毛片在线看网站| 天堂影院成人在线观看| 老熟妇仑乱视频hdxx| 成在线人永久免费视频| 国产片内射在线| 亚洲 欧美 日韩 在线 免费| 国产av不卡久久| 精品一区二区三区四区五区乱码| 又黄又爽又免费观看的视频| АⅤ资源中文在线天堂| 欧美成人免费av一区二区三区| 亚洲专区国产一区二区| 久久香蕉国产精品| 国产精品久久电影中文字幕| 国产一区二区三区视频了| 天天添夜夜摸| 久久久久精品国产欧美久久久| 午夜两性在线视频| 女性生殖器流出的白浆| 校园春色视频在线观看| 不卡一级毛片| 日韩视频一区二区在线观看| 日韩欧美在线二视频| 日本黄色视频三级网站网址| 黄色视频,在线免费观看| 老司机福利观看| 日本成人三级电影网站| 淫秽高清视频在线观看| 国产高清videossex| 亚洲电影在线观看av| 国产精品二区激情视频| 性欧美人与动物交配| 午夜久久久久精精品| 日日爽夜夜爽网站| 久久午夜综合久久蜜桃| 可以在线观看毛片的网站| 欧美三级亚洲精品| 人人妻,人人澡人人爽秒播| 国产91精品成人一区二区三区| 国产一区二区三区在线臀色熟女| 日韩大码丰满熟妇| 欧美中文日本在线观看视频| 精品不卡国产一区二区三区| 亚洲狠狠婷婷综合久久图片| 在线播放国产精品三级| 国产三级在线视频| 一区二区三区精品91| 麻豆一二三区av精品| 啦啦啦免费观看视频1| 一区二区三区精品91| 国产熟女xx| 999精品在线视频| 叶爱在线成人免费视频播放| 大香蕉久久成人网| 亚洲人成网站高清观看| 亚洲午夜精品一区,二区,三区| 男女下面进入的视频免费午夜 | 十八禁网站免费在线| 视频在线观看一区二区三区| 午夜久久久在线观看| 色综合婷婷激情| 亚洲专区国产一区二区| 国产高清视频在线播放一区| av欧美777| 久久人妻福利社区极品人妻图片| 免费在线观看黄色视频的| 国产亚洲精品综合一区在线观看 | 性色av乱码一区二区三区2| 精品久久久久久久久久久久久 | 99riav亚洲国产免费| 中文字幕高清在线视频| 久久久国产欧美日韩av| 88av欧美| 啦啦啦免费观看视频1| 中文字幕精品免费在线观看视频| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av| 中出人妻视频一区二区| 免费在线观看视频国产中文字幕亚洲| 一进一出抽搐gif免费好疼| 久久精品国产清高在天天线| 黄色视频,在线免费观看| 黄色女人牲交| 国产精品日韩av在线免费观看| 母亲3免费完整高清在线观看| 国产精品 欧美亚洲| 嫩草影视91久久| 免费人成视频x8x8入口观看| 欧美日本视频| 久久久久久国产a免费观看| 国产成人啪精品午夜网站| 嫩草影视91久久| 99riav亚洲国产免费| 中文字幕精品免费在线观看视频| a级毛片a级免费在线| 观看免费一级毛片| 午夜福利欧美成人| 午夜老司机福利片| 悠悠久久av| 日日干狠狠操夜夜爽| 国产激情偷乱视频一区二区| 欧美激情 高清一区二区三区| 亚洲免费av在线视频| 久久久久久免费高清国产稀缺| 香蕉丝袜av| 精品国产乱子伦一区二区三区| 一级a爱片免费观看的视频| 真人一进一出gif抽搐免费| 国产精品乱码一区二三区的特点| 琪琪午夜伦伦电影理论片6080| 国产精品永久免费网站| 午夜福利18| 精品欧美国产一区二区三| 女性被躁到高潮视频| av有码第一页| 黄片播放在线免费| av在线播放免费不卡| 人人妻人人澡人人看| 亚洲专区国产一区二区| 免费搜索国产男女视频| 又黄又粗又硬又大视频| 久久久国产精品麻豆| 午夜成年电影在线免费观看| 午夜影院日韩av| 久久久国产成人免费| 老熟妇乱子伦视频在线观看| 日本a在线网址| 日日爽夜夜爽网站| 深夜精品福利| 高清在线国产一区| av在线播放免费不卡| 国产日本99.免费观看| 日本黄色视频三级网站网址| www.www免费av| 欧美日韩瑟瑟在线播放| 国产一区二区三区在线臀色熟女| 我的亚洲天堂| 18禁黄网站禁片免费观看直播| av视频在线观看入口| 久久精品aⅴ一区二区三区四区| 国产一区二区在线av高清观看| 最近在线观看免费完整版| 精品一区二区三区av网在线观看| 身体一侧抽搐| 美女高潮到喷水免费观看| 18禁黄网站禁片午夜丰满| 看片在线看免费视频| 午夜福利在线在线| 欧美性猛交黑人性爽| 久久香蕉激情| 久久久久精品国产欧美久久久| 国产亚洲精品久久久久久毛片| 欧美激情 高清一区二区三区| 亚洲人成伊人成综合网2020| 自线自在国产av| 亚洲五月婷婷丁香| 夜夜看夜夜爽夜夜摸| 日韩av在线大香蕉| 91在线观看av| 露出奶头的视频| 91老司机精品| 国产91精品成人一区二区三区| 亚洲av电影不卡..在线观看| 日韩 欧美 亚洲 中文字幕| 99久久国产精品久久久| 亚洲成人久久性| 一边摸一边抽搐一进一小说| 国产精品亚洲美女久久久| 啦啦啦 在线观看视频| 中文字幕最新亚洲高清| 亚洲第一av免费看| 十分钟在线观看高清视频www| 91字幕亚洲| 国产欧美日韩一区二区精品| 国产亚洲精品久久久久久毛片| 最近最新中文字幕大全免费视频| 亚洲专区字幕在线| 国产午夜精品久久久久久| 黄片播放在线免费| 免费在线观看完整版高清| 亚洲国产中文字幕在线视频| 亚洲精品一区av在线观看| 欧美色视频一区免费| 成年版毛片免费区| 搡老妇女老女人老熟妇| 久久精品国产综合久久久| 视频区欧美日本亚洲| 欧美一级a爱片免费观看看 | or卡值多少钱| 免费一级毛片在线播放高清视频| АⅤ资源中文在线天堂| 最新美女视频免费是黄的| 男人舔女人下体高潮全视频| 精品午夜福利视频在线观看一区| www日本在线高清视频| 国产真人三级小视频在线观看| 99riav亚洲国产免费| 久久 成人 亚洲| 国产男靠女视频免费网站| 国产高清有码在线观看视频 | 色尼玛亚洲综合影院| 搡老岳熟女国产| 99热6这里只有精品| 12—13女人毛片做爰片一| 又大又爽又粗| 亚洲人成网站在线播放欧美日韩| 午夜福利在线观看吧| 亚洲国产精品合色在线| 99在线人妻在线中文字幕| xxx96com| 老汉色av国产亚洲站长工具| 国产aⅴ精品一区二区三区波| 99re在线观看精品视频| 色综合欧美亚洲国产小说| 中文字幕人成人乱码亚洲影| 欧美中文综合在线视频| 日韩欧美一区视频在线观看| 夜夜夜夜夜久久久久| 亚洲中文字幕日韩| 亚洲av日韩精品久久久久久密| 欧美中文日本在线观看视频| 亚洲av五月六月丁香网| 日韩欧美一区视频在线观看| 给我免费播放毛片高清在线观看| 看免费av毛片| 亚洲欧美日韩无卡精品| 亚洲黑人精品在线| 91老司机精品| 免费av毛片视频| 成人18禁在线播放| 欧美午夜高清在线| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| 亚洲黑人精品在线| 美国免费a级毛片| 国产蜜桃级精品一区二区三区| 99热6这里只有精品| 麻豆成人av在线观看| 99国产综合亚洲精品| 国产成人一区二区三区免费视频网站| 欧美不卡视频在线免费观看 | 成人国产一区最新在线观看| 日本一区二区免费在线视频| 日本五十路高清| bbb黄色大片| 亚洲专区中文字幕在线| 国产日本99.免费观看| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 色在线成人网| 亚洲精品国产一区二区精华液| 丰满的人妻完整版| 一本精品99久久精品77| 欧美激情极品国产一区二区三区| 亚洲av片天天在线观看| 国产伦人伦偷精品视频| 曰老女人黄片| 99久久无色码亚洲精品果冻| 精品国产超薄肉色丝袜足j| tocl精华| 嫩草影视91久久| 这个男人来自地球电影免费观看| 国产伦一二天堂av在线观看| АⅤ资源中文在线天堂| 精品国产国语对白av| 久久人人精品亚洲av| 在线观看免费日韩欧美大片| 国产精品乱码一区二三区的特点| 亚洲性夜色夜夜综合| 久久久久免费精品人妻一区二区 | 成人欧美大片| 国产一区二区三区在线臀色熟女| 亚洲最大成人中文| 精品久久久久久久末码| 亚洲欧美精品综合一区二区三区| 丝袜人妻中文字幕| 亚洲 欧美 日韩 在线 免费| 中文字幕人妻熟女乱码| 亚洲精品在线观看二区| 中出人妻视频一区二区| 男人舔女人下体高潮全视频| 欧美中文综合在线视频| 成人三级黄色视频| 国产成人一区二区三区免费视频网站| 夜夜看夜夜爽夜夜摸| а√天堂www在线а√下载| 黄片播放在线免费| 国产免费av片在线观看野外av| 三级毛片av免费| 国产精品1区2区在线观看.| 国产精品精品国产色婷婷| 国产aⅴ精品一区二区三区波| 国产成人系列免费观看| 亚洲成人久久爱视频| 动漫黄色视频在线观看| 久久久久九九精品影院| 免费在线观看影片大全网站| 1024手机看黄色片| 国产熟女xx| 免费电影在线观看免费观看| 国产亚洲欧美98| 在线观看免费视频日本深夜| 国内毛片毛片毛片毛片毛片| 啦啦啦韩国在线观看视频| 久久狼人影院| 久久久精品国产亚洲av高清涩受| 在线观看66精品国产| 久久亚洲真实| 身体一侧抽搐| 日本撒尿小便嘘嘘汇集6| 欧美在线黄色| 一进一出抽搐动态| 国产亚洲精品久久久久久毛片| 女性生殖器流出的白浆| 法律面前人人平等表现在哪些方面| 精品久久久久久,| 国产精品 国内视频| 国产亚洲精品av在线| 欧美黄色片欧美黄色片| 国产区一区二久久| 日本成人三级电影网站| 国产三级黄色录像| 国产又黄又爽又无遮挡在线| 亚洲国产精品sss在线观看| 一进一出抽搐动态| 一区二区三区激情视频| 这个男人来自地球电影免费观看| 精品国产超薄肉色丝袜足j| 无人区码免费观看不卡| 免费看a级黄色片| 成人亚洲精品av一区二区| 久久天堂一区二区三区四区| 国产激情欧美一区二区| 51午夜福利影视在线观看| 国产伦一二天堂av在线观看| 91成年电影在线观看| 久久热在线av| 欧美又色又爽又黄视频| 精品熟女少妇八av免费久了| 亚洲精品中文字幕在线视频| 亚洲黑人精品在线| 中文字幕精品免费在线观看视频| 久久久久久久午夜电影| 人妻丰满熟妇av一区二区三区| 久热这里只有精品99| 人妻久久中文字幕网| 人人妻人人澡欧美一区二区| 久99久视频精品免费| 日韩欧美国产一区二区入口| 人妻久久中文字幕网| 日韩欧美一区二区三区在线观看| 看片在线看免费视频| 国产高清有码在线观看视频 | 亚洲中文字幕一区二区三区有码在线看 | 搡老岳熟女国产| 最新美女视频免费是黄的| 91成年电影在线观看| 黄色成人免费大全| 男人操女人黄网站| 村上凉子中文字幕在线| 成人18禁在线播放| 欧美性猛交黑人性爽| 欧美日韩精品网址| 中文字幕人妻熟女乱码| 欧美在线黄色| 国产亚洲欧美精品永久| 亚洲国产精品999在线| 精品久久蜜臀av无| 美女免费视频网站| 欧美丝袜亚洲另类 | 黄片小视频在线播放| 伦理电影免费视频| 婷婷六月久久综合丁香| 天堂影院成人在线观看| 国产一区二区三区在线臀色熟女| 国产免费男女视频| 51午夜福利影视在线观看| 757午夜福利合集在线观看| 国产成人精品无人区| 国产精品亚洲美女久久久| 成人亚洲精品一区在线观看| 日日夜夜操网爽| 天堂√8在线中文| 日韩大尺度精品在线看网址| 亚洲久久久国产精品| 一级毛片高清免费大全| 99国产精品一区二区蜜桃av| 999精品在线视频| 国产精品亚洲av一区麻豆| 久久天堂一区二区三区四区| 国产亚洲精品av在线| 亚洲五月婷婷丁香| 国产av在哪里看| 亚洲精品国产一区二区精华液| 一二三四社区在线视频社区8| 成人永久免费在线观看视频| 99riav亚洲国产免费| 亚洲精品一区av在线观看| 欧美大码av| 成在线人永久免费视频| 亚洲三区欧美一区| 91麻豆av在线| 国产黄a三级三级三级人| 777久久人妻少妇嫩草av网站| 国产麻豆成人av免费视频| 香蕉av资源在线| 女人爽到高潮嗷嗷叫在线视频| 桃红色精品国产亚洲av| 精华霜和精华液先用哪个| 久久久久久久久免费视频了| 成人午夜高清在线视频 | 亚洲黑人精品在线| 欧美成人午夜精品| 亚洲av第一区精品v没综合| 又黄又爽又免费观看的视频| 老鸭窝网址在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月天丁香| 欧美三级亚洲精品| e午夜精品久久久久久久| 欧美黄色片欧美黄色片| 一本综合久久免费| 香蕉国产在线看| 美女午夜性视频免费| 亚洲中文日韩欧美视频| 亚洲国产欧美一区二区综合| 久久久久亚洲av毛片大全| 最新美女视频免费是黄的| 国产黄色小视频在线观看| 国内少妇人妻偷人精品xxx网站 | 十分钟在线观看高清视频www| 深夜精品福利| 亚洲一区中文字幕在线| 亚洲人成电影免费在线| 1024香蕉在线观看| 窝窝影院91人妻| 国内毛片毛片毛片毛片毛片| 国产黄片美女视频| 一边摸一边抽搐一进一小说| 国产午夜福利久久久久久| 亚洲色图 男人天堂 中文字幕| 天天躁夜夜躁狠狠躁躁| 国产日本99.免费观看| 亚洲片人在线观看| 一级毛片精品| 亚洲五月天丁香| av天堂在线播放| 久久久久久亚洲精品国产蜜桃av|