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

    線(xiàn)性Radon變換噪音壓制法及其在古龍斷陷中的應(yīng)用

    2012-12-14 02:56:04王維紅陸加敏
    關(guān)鍵詞:同相軸古龍壓制

    石 穎,李 瑩,王維紅,陸加敏

    (1.東北石油大學(xué) 地球科學(xué)學(xué)院,黑龍江 大慶 163318; 2.大慶油田有限責(zé)任公司 勘探開(kāi)發(fā)研究院,黑龍江 大慶163712)

    線(xiàn)性Radon變換噪音壓制法及其在古龍斷陷中的應(yīng)用

    石 穎1,李 瑩1,王維紅2,陸加敏2

    (1.東北石油大學(xué) 地球科學(xué)學(xué)院,黑龍江 大慶 163318; 2.大慶油田有限責(zé)任公司 勘探開(kāi)發(fā)研究院,黑龍江 大慶163712)

    針對(duì)松遼盆地北部古龍斷陷地震資料信噪比低、線(xiàn)性干擾強(qiáng)特點(diǎn),提出應(yīng)用線(xiàn)性Radon變換進(jìn)行疊前線(xiàn)性噪音壓制的預(yù)處理方法,Radon變換可在炮集和CMP道集上進(jìn)行運(yùn)算,算法簡(jiǎn)單,易于編程實(shí)現(xiàn),其積分路徑的特點(diǎn)適合線(xiàn)性噪音壓制.模擬數(shù)據(jù)和實(shí)際地震資料應(yīng)用結(jié)果表明,線(xiàn)性Radon變換法能夠?qū)崿F(xiàn)保幅的線(xiàn)性噪音壓制,是疊前提高地震資料信噪比的實(shí)用方法,在地震資料預(yù)處理中具有應(yīng)用前景.

    Radon變換;線(xiàn)性噪音;壓制;信噪比;古龍斷陷;地震資料

    0 引言

    當(dāng)前油氣勘探以復(fù)雜構(gòu)造和巖性油氣藏為主,要求高保真和高信噪比的地震成果數(shù)據(jù)體,以利于進(jìn)行綜合地質(zhì)研究和勘探部署.地震資料噪音主要包括隨機(jī)噪音和規(guī)則噪音兩類(lèi),其中規(guī)則噪音中的線(xiàn)性噪音在地震資料中普遍存在,往往具有較強(qiáng)的能量,所以線(xiàn)性噪音壓制效果在一定程度上決定著成像數(shù)據(jù)體質(zhì)量[].

    根據(jù)有效信號(hào)與線(xiàn)性噪聲的特征差異,研究并提出了不同的線(xiàn)性噪音壓制方法,主要包括F-K濾波法、徑向?yàn)V波法、K-L變換法和Radon變換法等[2-6].其中F-K濾波法是頻率域方法,易于實(shí)現(xiàn);但是在有效波和線(xiàn)性干擾之間的視速度相近時(shí),難以有效實(shí)現(xiàn)有效波和噪音的分離.徑向?yàn)V波法根據(jù)選定的線(xiàn)性同相軸的視速度和濾波的徑向道數(shù),實(shí)現(xiàn)線(xiàn)性干擾的壓制;但是該方法需要過(guò)多的人為參與,對(duì)于海量地震數(shù)據(jù)處理,計(jì)算效率低,計(jì)算結(jié)果在很大程度上取決于濾波道數(shù)的選擇,且地震數(shù)據(jù)體中往往有較多的殘余噪音.K-L變換法是基于目標(biāo)統(tǒng)計(jì)特征的正交變換,能夠?qū)崿F(xiàn)隨機(jī)噪音、線(xiàn)性噪音和多次波等規(guī)則噪音壓制;但是該方法噪音壓制效果取決于本征值選擇,本征值選擇需要反復(fù)的數(shù)值試驗(yàn),若本征值選擇不當(dāng)將造成噪音的大量剩余,同時(shí)有效波的信號(hào)也受到一定的損失.

    線(xiàn)性Radon變換法算法簡(jiǎn)單,在地震資料處理中被廣泛應(yīng)用,如資料預(yù)處理中的線(xiàn)性噪音壓制、地震數(shù)據(jù)的插值和重建、VSP資料處理中的波場(chǎng)分離和疊前偏移領(lǐng)域等[6-9].線(xiàn)性Radon變換法能夠較好地實(shí)現(xiàn)有效波和相干線(xiàn)性噪音的分離,計(jì)算效率高,線(xiàn)性同相軸在Radon正變換域可變換為點(diǎn),濾波函數(shù)簡(jiǎn)單,同時(shí)經(jīng)噪音模型反變換后,從原始數(shù)據(jù)中將線(xiàn)性噪音減去,對(duì)有效波的振幅沒(méi)有任何傷害,是一種高效的相對(duì)保幅的線(xiàn)性噪音壓制技術(shù)[10-12].依據(jù)線(xiàn)性Radon變換原理,對(duì)大慶油田古龍斷陷區(qū)的地震資料進(jìn)行強(qiáng)線(xiàn)性干擾壓制,獲得良好成像效果,試算結(jié)果表明算法具有計(jì)算效率高、精度高和實(shí)用性強(qiáng)的特點(diǎn).

    1 Radon變換原理與實(shí)現(xiàn)

    對(duì)于二維地震資料,假設(shè)沿偏移距方向?yàn)榫鶆虿蓸?,則離散線(xiàn)性Radon正變換[9,13]可以寫(xiě)為

    式中:d(xi,t)為x-t域地震數(shù)據(jù)體,共有N 道數(shù)據(jù),xi為偏移距;m(τ,p)為Radon正變換域數(shù)據(jù)體,τ為地震數(shù)據(jù)中零偏移距道對(duì)應(yīng)的截距時(shí)間;p為射線(xiàn)參數(shù).

    對(duì)Radon變換的離散公式(1)兩端進(jìn)行Fourier變換,對(duì)任一頻率分量進(jìn)行計(jì)算[9,11]:

    式中:M(p,ω),D(xi,ω)分別為m(τ,p)和d(xi,t)相對(duì)應(yīng)的Fourier變換域的表示形式.將方程(2)寫(xiě)成矩陣形式:

    式中:d為時(shí)間空間域地震數(shù)據(jù)矩陣;m為Radon域矩陣;L為數(shù)據(jù)變換矩陣,L中的元素可表示為

    一般來(lái)說(shuō),方程(3)是欠定或超定的,采用最小二乘方法求其最優(yōu)解,其一般形式為

    式中:LH為矩陣L的共軛轉(zhuǎn)置;m,d分別為模型空間(Radon正變換域)和數(shù)據(jù)空間(時(shí)間空間域)的任一頻率成分的向量.

    式(5)和式(3)共同構(gòu)成線(xiàn)性Radon正反變換對(duì).在離散計(jì)算時(shí),Radon變換進(jìn)行參數(shù)的合理采樣.線(xiàn)性Radon變換的參數(shù)τ和時(shí)間域的參數(shù)t的采樣相等,其射線(xiàn)參數(shù)p的采樣間隔Δp為

    最大射線(xiàn)參數(shù)pmax的表達(dá)式為

    式(6-7)中:fmax為地震資料的最大有效頻率;Xmax為原始數(shù)據(jù)單炮或CMP道集的最大偏移距;Δx為道集中的道間距[14].

    線(xiàn)性Radon變換求解方法是矩陣求逆,式(5)可能是超定方程(未知數(shù)個(gè)數(shù)少于方程個(gè)數(shù)),也可能是欠定方程(未知數(shù)個(gè)數(shù)多于方程個(gè)數(shù)),對(duì)于文中線(xiàn)性噪音壓制問(wèn)題,單個(gè)頻率成分所形成的矩陣為T(mén)oeplitz矩陣,可以應(yīng)用共軛梯度法、Cholesky分解法或Levinson遞推法進(jìn)行求解,其中Levinson遞推算法的計(jì)算效率最高[13].

    2 模擬數(shù)據(jù)試驗(yàn)

    為驗(yàn)證Radon變換線(xiàn)性噪音壓制方法的有效性,應(yīng)用理論模擬的單炮數(shù)據(jù)進(jìn)行噪音壓制試算.數(shù)值模擬含一個(gè)有效波同相軸和一個(gè)線(xiàn)性同相軸的單炮數(shù)據(jù)(見(jiàn)圖1(a)).由圖1(a)可以看出,模擬線(xiàn)性噪音的能量很強(qiáng).模擬采用主頻為20 Hz的Ricker子波,共84道地震記錄,地震道采樣間隔為40 m,時(shí)間方向采樣間隔為4 ms,最大偏移距為3 320 m.采用式(5)計(jì)算最小二乘Radon域的正變換域見(jiàn)圖1(b).由圖1(b)可以看出,模擬炮記錄中線(xiàn)性同相軸在線(xiàn)性Radon正變換域表現(xiàn)為相對(duì)聚焦的剪刀狀特征,而有效波表現(xiàn)為橢圓形式.在Radon正變換的模型空間域設(shè)計(jì)濾波函數(shù),可以容易地濾除點(diǎn)狀形式的線(xiàn)性同相軸,將濾波的結(jié)果進(jìn)行反變換,即可得到線(xiàn)性同相軸(見(jiàn)圖1(c)),將圖1(c)從原始地震數(shù)據(jù)中減去,得到線(xiàn)性噪音壓制后的結(jié)果剖面(見(jiàn)圖1(d)).因此,線(xiàn)性Radon變換法可以有效壓制地震數(shù)據(jù)中的線(xiàn)性同相軸.

    圖1 Radon變換線(xiàn)性噪音壓制模擬數(shù)據(jù)試驗(yàn)

    線(xiàn)性Radon變換法的基本思想是對(duì)某個(gè)函數(shù)在給定的線(xiàn)性路徑上進(jìn)行積分運(yùn)算,基于平面波分解原理,也考慮地震波場(chǎng)的性質(zhì).根據(jù)線(xiàn)性Radon域的特點(diǎn),可以設(shè)計(jì)濾波器,在Radon正變換域?yàn)V波,以切除不需要的同相軸或者相干噪音,進(jìn)而大幅提高地震資料的信噪比,表明Radon變換法在線(xiàn)性噪音壓制計(jì)算中的有效性和實(shí)用性.

    3 現(xiàn)場(chǎng)資料應(yīng)用

    松遼盆地北部深層天然氣勘探在徐家圍子斷陷取得成功,獲得儲(chǔ)量的巨大發(fā)現(xiàn)[15],當(dāng)前大慶油田探區(qū)的古龍斷陷是重要的儲(chǔ)量接替區(qū).該區(qū)針對(duì)深層火山巖勘探的研究程度相對(duì)較低,其中影響因素主要包括火山巖目的層埋藏深、地震資料品質(zhì)差和信噪比低等.為提高古龍斷陷深層復(fù)雜構(gòu)造和復(fù)雜地質(zhì)體的成像精度,采用Kirchhoff積分的疊前深度偏移技術(shù)[16].該技術(shù)要求輸入信噪比高的疊前預(yù)處理數(shù)據(jù)體,所以疊前精細(xì)噪音壓制是提高成像精度的關(guān)鍵環(huán)節(jié).

    以古龍斷陷的葡南工區(qū)為例,原始地震資料分析表明,該地區(qū)地表?xiàng)l件復(fù)雜,折射干擾嚴(yán)重,發(fā)育兩組折射波,速度分別為1 750 m/s和2 200 m/s.線(xiàn)性噪音具有高能量、強(qiáng)振幅和分布范圍廣的特點(diǎn),且有較強(qiáng)的規(guī)律性.全區(qū)測(cè)線(xiàn)或多或少含有這種干擾,幾乎分布于整個(gè)單炮記錄,從而掩蓋有效波組.一般而言,根據(jù)線(xiàn)性干擾與有效波在速度、頻率、時(shí)空上的差別,進(jìn)行線(xiàn)性干擾的識(shí)別和壓制.在古龍斷陷區(qū),線(xiàn)性干擾波具有相對(duì)穩(wěn)定的視速度,所以可以應(yīng)用線(xiàn)性Radon變換法實(shí)現(xiàn)線(xiàn)性干擾的有效壓制.另外,在實(shí)現(xiàn)過(guò)程中同時(shí)采用噪音壓制的減去法,使得有效波的保幅性得到很大程度的增強(qiáng).

    Radon變換法線(xiàn)性干擾壓制前后的單炮效果見(jiàn)圖2.由圖2(a)可以看出,線(xiàn)性噪音呈排狀自上而下分布,幾乎掩蓋所有的有效波;由圖2(b)可以看出,壓制干擾后的炮記錄較為清晰地顯示拋物線(xiàn)型的有效波同相軸,單炮記錄的信噪比得到提高;由圖2(c)可以看出,該方法能夠有效壓制古龍斷陷的強(qiáng)能量線(xiàn)性噪音,同時(shí)有效波的振幅不受損失,即文中方法具有很好的保幅效果.

    圖2 線(xiàn)性干擾壓制前后炮記錄對(duì)比

    線(xiàn)性干擾壓制前后疊加剖面的效果見(jiàn)圖3.由圖3可以看出,經(jīng)過(guò)線(xiàn)性噪音壓制后,地震剖面信噪比明顯提高.通過(guò)圖2和圖3的對(duì)比,無(wú)論是從單炮還是從疊加剖面都可以看出,線(xiàn)性干擾得到有效壓制,突出有效波能量,地震資料的信噪比大幅提高,為疊前深度偏移成像提供資料保障.古龍斷陷的實(shí)際地震資料線(xiàn)性噪音壓制表明,文中Radon變換法具有很強(qiáng)的實(shí)用性,能夠在信噪比低的同類(lèi)地震資料處理中推廣應(yīng)用.

    圖3 線(xiàn)性干擾壓制前后疊加剖面對(duì)比

    4 結(jié)論

    (1)線(xiàn)性Radon變換法是對(duì)給定函數(shù)沿一定路徑的積分計(jì)算,實(shí)際上也考慮地震波場(chǎng)的性質(zhì),其變換形式可從平面波分解的原理進(jìn)行求解和推導(dǎo).線(xiàn)性Radon變換原理簡(jiǎn)單,在地震資料預(yù)處理中應(yīng)用廣泛.

    (2)線(xiàn)性Radon變換法能夠較好地實(shí)現(xiàn)有效波和相干線(xiàn)性噪音的分離.該方法計(jì)算效率高,可以用于炮集或CMP道集數(shù)據(jù),線(xiàn)性同相軸在Radon正變換域可變換為點(diǎn),濾波函數(shù)簡(jiǎn)單,同時(shí)經(jīng)噪音模型反變換后,從原始數(shù)據(jù)中減去,能夠?qū)崿F(xiàn)線(xiàn)性噪音的保幅壓制.

    (3)松遼盆地古龍斷陷地震資料的線(xiàn)性噪音較為發(fā)育,在部分炮記錄上具有很強(qiáng)的能量,應(yīng)用線(xiàn)性Radon變換法進(jìn)行線(xiàn)性噪音壓制,實(shí)際地震數(shù)據(jù)線(xiàn)性噪音壓制結(jié)果表明,文中給出的線(xiàn)性噪音壓制方法具有有效性和實(shí)用性,對(duì)后續(xù)地震波精確成像具有指導(dǎo)意義.

    [1]胡天躍.地震資料疊前去噪技術(shù)的現(xiàn)狀與未來(lái)[J].地球物理學(xué)進(jìn)展,2002,17(2):35-40.

    [2]閆立志,景新義,李剛.F-K濾波在噪音減去法中的應(yīng)用[J].海洋地質(zhì)動(dòng)態(tài),2006,22(10):28-32.

    [3]余波,黃中玉,談大龍,等.徑向道濾波法去線(xiàn)性干擾[J].石油物探,2005,44(2):109-112.

    [4]黃雪繼,劉來(lái)祥,王永勝.分頻徑向道中值濾波在地震資料處理中的應(yīng)用[J].物探與化探,2012,36(2):317-320.

    [5]劉洪林,張春堂,朱秋影,等.K-L變換在地震資料去噪中的應(yīng)用[J].大慶石油學(xué)院學(xué)報(bào),2007,31(4):19-21.

    [6]吳律,武克奮,孫立.變換方法及其在地震資料處理中的應(yīng)用[J].石油物探,1996,1(1):37-52.

    [7]Kabir M N,Verschuur D J.Restoration of missing offsets by parabolic Radon transform[J].Geophys.Prosp.,1995,43:347-368.

    [8]王維紅,劉洪.拋物 Radon變換法近偏移距波場(chǎng)外推[J].地球物理學(xué)進(jìn)展,2005,20(2):289-293.

    [9]Schultz P S,Claerbout J F.Velocity estimation and downward continuation by wavefront synthesis[J].Geophysics,1978,43(3):691-714.

    [10]林文,魏大力,王建民,等.Radon變換多次波壓制方法及應(yīng)用研究[J].物探化探計(jì)算技術(shù),2009,31(4):344-348.

    [11]王維紅,首皓,劉洪,等.線(xiàn)性同相軸波場(chǎng)分離的高分辨率τ-p變換法[J].地球物理學(xué)進(jìn)展,2006,21(1):74-78.

    [12]Cambois G.Preserved amplitude processing in the presence of noise[C].66th Annual Internat.Mtg.,Soc.Expl.Geophys.,SEG Expanded Abstracts,1996:1595-1598.

    [13]Kostov C.Toeplitz structure in slant-stack inversion[C].60th Annual Internat.Mtg.,Soc.Expl.Geophys.,SEG Expanded Abstracts,1990:1618-1621.

    [14]張旭東.幾種多次波壓制方法討論[J].內(nèi)蒙古石油化工,2011(1):84-85.

    [15]胡明,遲建功,呂延防,等.徐家圍子斷陷徐深21區(qū)塊蓋層封閉能力演化與天然氣成藏的關(guān)系[J].大慶石油學(xué)院學(xué)報(bào),2010,34(1):24-28.

    [16]王維紅,林春華,陳志德,等.古龍斷陷深層火山巖地震資料成像方法及應(yīng)用研究[J].地球物理學(xué)報(bào),2011,54(2):310-319.

    Prestack reverse time migration based on GPU parallel accelerating algorithm/2012,36(4):111-115

    SHI Ying1,LU Jia-min2,KE Xuan1,TIAN Dong-sheng1,WANG Fei3
    (1.School of Geosciences,Northeast Petroleum University,Daqing,Heilongjiang 163318,China;2.Exploration and Development Research Institute,Daqing Oilfield Co.Ltd.,Daqing,Heilongjiang 163712,China;3.School of Earth Science,China Petroleum University (Beijing),Beijing 102200,China)

    In order to improve the complex subsurface imaging accuracy and computational efficiency of the algorithm,this paper presents an algorithm of prestack reverse time migration based on GPU(Graphic Processing Unit)accelerating which can image the underground complex structure effectively and accurately.By two-way wave equation to calculate wave field extrapolation,prestack reverse-time migration can overcome the dip limit,and the imaging algorithm is performed by high order finite difference in the paper.Wave field extrapolation and imaging condition are calculated by GPU parallel accelerating technology,comparing to conventional algorithm,its computation efficiency has been greatly improved,and it meets large amount of computation requirement in prestack reverse-time migration.The random boundary condition approach is adopted to obtain wavefield information,which reduces the memory demand but sacrifices the computation cost,and it solves the massy memory problem in reverse time migration.The tests on model illustrate that this approach can imaging complicated geological body efficiently and precisely.

    reverse time migration;GPU;acceleration;high order finite difference;random boundary condition;complicated structure

    TE132.1

    A

    2095-4107(2012)04-0116-05

    DOI 10.3969/j.issn.2095-4107.2012.04.021

    2012-05-07;編輯:任志平

    國(guó)家自然科學(xué)基金青年基金項(xiàng)目(41104088,41004057);國(guó)家“863”高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2012AA061202);中國(guó)博士后科學(xué)基金項(xiàng)目(2011M501009);中國(guó)石油科技創(chuàng)新基金項(xiàng)目(2011D-5006-0304);黑龍江省教育廳科學(xué)技術(shù)研究項(xiàng)目(12511025);黑龍江省博士后科學(xué)基金項(xiàng)目(LBH-Z11272)

    石 穎(1976-),女,博士,副教授,主要從事地震資料處理方面的研究.

    猜你喜歡
    同相軸古龍壓制
    虛同相軸方法及其在陸上地震層間多次波壓制中的應(yīng)用
    一種新型無(wú)人機(jī)數(shù)據(jù)鏈抗壓制干擾技術(shù)的研究
    活到四十多歲,才知道古龍有多天才
    文苑(2018年20期)2018-11-09 01:36:02
    空射誘餌在防空壓制電子戰(zhàn)中的應(yīng)用
    一種改進(jìn)的相關(guān)法自動(dòng)拾取同相軸
    一種反射同相軸自動(dòng)拾取算法
    一種舊物品擠壓成型機(jī)
    科技資訊(2016年12期)2016-05-30 05:07:58
    朋友間的分寸
    當(dāng)代香帥古龍
    對(duì)GPS接收機(jī)帶限高斯噪聲壓制干擾的干擾帶寬選擇分析
    十分钟在线观看高清视频www| 欧美人与性动交α欧美精品济南到 | 两个人看的免费小视频| 亚洲人成电影观看| 欧美日韩视频精品一区| 麻豆精品久久久久久蜜桃| 国产精品99久久99久久久不卡 | 伊人久久大香线蕉亚洲五| 国产亚洲av片在线观看秒播厂| 国产精品蜜桃在线观看| 天天躁夜夜躁狠狠躁躁| 国产乱人偷精品视频| 一区福利在线观看| 美女午夜性视频免费| 日韩精品免费视频一区二区三区| 日韩,欧美,国产一区二区三区| 在线看a的网站| 免费人妻精品一区二区三区视频| 99香蕉大伊视频| 中文字幕色久视频| 国产精品亚洲av一区麻豆 | 日本-黄色视频高清免费观看| 久久女婷五月综合色啪小说| 女人久久www免费人成看片| videosex国产| 永久免费av网站大全| 熟女电影av网| 一本色道久久久久久精品综合| 国产精品久久久久久av不卡| 在线观看www视频免费| 国产一级毛片在线| 中文字幕人妻丝袜一区二区 | 精品国产乱码久久久久久小说| 99热网站在线观看| 啦啦啦在线免费观看视频4| 亚洲婷婷狠狠爱综合网| 国产野战对白在线观看| 欧美日韩综合久久久久久| 哪个播放器可以免费观看大片| 国产高清国产精品国产三级| 成人二区视频| 国产一区二区激情短视频 | av在线app专区| 国产精品一区二区在线观看99| 日本午夜av视频| 国产一区二区 视频在线| 国产精品久久久久久精品电影小说| 一个人免费看片子| 国产av精品麻豆| 久久久久精品久久久久真实原创| 热re99久久国产66热| 女人被躁到高潮嗷嗷叫费观| 亚洲,欧美,日韩| 这个男人来自地球电影免费观看 | 欧美日韩av久久| 免费少妇av软件| 国产免费福利视频在线观看| 精品午夜福利在线看| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产亚洲av高清一级| 十分钟在线观看高清视频www| 丰满少妇做爰视频| 日本欧美视频一区| 免费高清在线观看视频在线观看| 精品一区在线观看国产| 国产 精品1| 亚洲国产看品久久| 这个男人来自地球电影免费观看 | 麻豆乱淫一区二区| 亚洲精品美女久久久久99蜜臀 | 亚洲av日韩在线播放| 免费不卡的大黄色大毛片视频在线观看| 日日爽夜夜爽网站| videos熟女内射| 精品国产一区二区三区久久久樱花| 性少妇av在线| 日韩免费高清中文字幕av| 涩涩av久久男人的天堂| 十八禁高潮呻吟视频| 国产av一区二区精品久久| a级毛片在线看网站| 多毛熟女@视频| 97人妻天天添夜夜摸| 午夜福利乱码中文字幕| 免费人妻精品一区二区三区视频| 一二三四中文在线观看免费高清| 丁香六月天网| 中文天堂在线官网| 亚洲情色 制服丝袜| 久久久精品区二区三区| 性高湖久久久久久久久免费观看| 18禁裸乳无遮挡动漫免费视频| 一本久久精品| 80岁老熟妇乱子伦牲交| 中文字幕另类日韩欧美亚洲嫩草| 久久99蜜桃精品久久| 黄片无遮挡物在线观看| 免费高清在线观看视频在线观看| 一区二区三区精品91| 大香蕉久久成人网| 成人国产av品久久久| 亚洲成av片中文字幕在线观看 | 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 国产在线免费精品| 18禁裸乳无遮挡动漫免费视频| www日本在线高清视频| 亚洲第一av免费看| 国产成人精品久久二区二区91 | 日日爽夜夜爽网站| 久久国产精品男人的天堂亚洲| 大话2 男鬼变身卡| 免费久久久久久久精品成人欧美视频| 性色av一级| 久久精品国产亚洲av涩爱| 最黄视频免费看| 男男h啪啪无遮挡| 日日啪夜夜爽| 伊人亚洲综合成人网| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 纯流量卡能插随身wifi吗| 久久女婷五月综合色啪小说| 亚洲成国产人片在线观看| 色婷婷久久久亚洲欧美| 不卡视频在线观看欧美| 爱豆传媒免费全集在线观看| 秋霞在线观看毛片| 亚洲人成77777在线视频| 中文字幕色久视频| 天堂俺去俺来也www色官网| 国产一区二区三区综合在线观看| 色吧在线观看| 欧美成人午夜精品| 国产成人精品在线电影| www.精华液| 精品国产超薄肉色丝袜足j| 久久鲁丝午夜福利片| 搡女人真爽免费视频火全软件| 一本大道久久a久久精品| h视频一区二区三区| 婷婷色麻豆天堂久久| 青青草视频在线视频观看| 在线亚洲精品国产二区图片欧美| 精品国产一区二区三区久久久樱花| 又粗又硬又长又爽又黄的视频| 狂野欧美激情性bbbbbb| 人人妻人人爽人人添夜夜欢视频| 亚洲av.av天堂| www日本在线高清视频| 亚洲视频免费观看视频| 国产精品秋霞免费鲁丝片| 最近手机中文字幕大全| 国产精品无大码| 一区福利在线观看| 2022亚洲国产成人精品| 午夜福利视频精品| 两个人免费观看高清视频| 青草久久国产| 啦啦啦视频在线资源免费观看| 国产日韩欧美亚洲二区| av在线观看视频网站免费| 啦啦啦在线免费观看视频4| 午夜福利视频在线观看免费| 久久免费观看电影| 亚洲成人手机| 熟女av电影| 亚洲国产av新网站| 人妻系列 视频| 亚洲经典国产精华液单| 国产人伦9x9x在线观看 | 电影成人av| 国产成人av激情在线播放| 日韩制服丝袜自拍偷拍| av网站免费在线观看视频| 国产综合精华液| 黄色 视频免费看| 亚洲国产毛片av蜜桃av| 9色porny在线观看| 国产亚洲午夜精品一区二区久久| 街头女战士在线观看网站| 国产男女内射视频| 热99国产精品久久久久久7| 少妇猛男粗大的猛烈进出视频| 观看av在线不卡| 午夜激情久久久久久久| 午夜福利影视在线免费观看| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 成人国语在线视频| 久久久久国产一级毛片高清牌| 最近中文字幕高清免费大全6| 夜夜骑夜夜射夜夜干| 极品少妇高潮喷水抽搐| 国产福利在线免费观看视频| 精品一区二区免费观看| 超碰成人久久| 欧美成人午夜免费资源| 日韩,欧美,国产一区二区三区| 国产国语露脸激情在线看| 亚洲欧洲日产国产| 亚洲av日韩在线播放| 大香蕉久久成人网| 青春草亚洲视频在线观看| 三级国产精品片| 在线天堂中文资源库| 国产熟女欧美一区二区| 日韩av在线免费看完整版不卡| 各种免费的搞黄视频| 久久人人爽av亚洲精品天堂| 婷婷色av中文字幕| av视频免费观看在线观看| 精品国产乱码久久久久久男人| 熟女电影av网| 男女午夜视频在线观看| 中文字幕人妻丝袜一区二区 | 精品国产一区二区三区四区第35| 国产淫语在线视频| 在线亚洲精品国产二区图片欧美| 国产亚洲av片在线观看秒播厂| 久久亚洲国产成人精品v| 99久久精品国产国产毛片| av.在线天堂| 成人毛片a级毛片在线播放| 国产成人免费观看mmmm| 免费在线观看视频国产中文字幕亚洲 | 搡老乐熟女国产| 99热全是精品| 国产午夜精品一二区理论片| 久久免费观看电影| 精品久久蜜臀av无| 97人妻天天添夜夜摸| 精品国产一区二区三区四区第35| 永久网站在线| 午夜av观看不卡| 激情五月婷婷亚洲| 亚洲国产色片| 乱人伦中国视频| 色94色欧美一区二区| 国产精品无大码| 欧美成人午夜免费资源| 日韩在线高清观看一区二区三区| 久久久久久久久久久久大奶| 日本猛色少妇xxxxx猛交久久| 亚洲欧美中文字幕日韩二区| 精品一区二区三卡| 亚洲激情五月婷婷啪啪| 日本爱情动作片www.在线观看| 午夜福利视频在线观看免费| 在线看a的网站| 色视频在线一区二区三区| 激情视频va一区二区三区| 国产精品二区激情视频| 欧美精品人与动牲交sv欧美| 国产成人免费观看mmmm| 一本大道久久a久久精品| 天天影视国产精品| 高清黄色对白视频在线免费看| 精品国产一区二区久久| 性色av一级| 国产av码专区亚洲av| 国产男人的电影天堂91| 亚洲国产最新在线播放| 看非洲黑人一级黄片| 国产熟女午夜一区二区三区| 久久久久人妻精品一区果冻| 晚上一个人看的免费电影| 欧美黄色片欧美黄色片| 91午夜精品亚洲一区二区三区| 亚洲欧洲日产国产| 亚洲婷婷狠狠爱综合网| 国产又爽黄色视频| 香蕉丝袜av| 国产精品三级大全| 免费在线观看视频国产中文字幕亚洲 | 久久精品久久久久久久性| 精品国产一区二区三区四区第35| 人妻一区二区av| 99香蕉大伊视频| 欧美 亚洲 国产 日韩一| 国产免费又黄又爽又色| 欧美激情高清一区二区三区 | 有码 亚洲区| 最近手机中文字幕大全| 人妻 亚洲 视频| 在线观看免费日韩欧美大片| 午夜免费鲁丝| 老女人水多毛片| 免费观看在线日韩| 最近最新中文字幕免费大全7| 91精品国产国语对白视频| 晚上一个人看的免费电影| 亚洲人成77777在线视频| 日本wwww免费看| 高清欧美精品videossex| 男人添女人高潮全过程视频| 一级片免费观看大全| 免费播放大片免费观看视频在线观看| 2021少妇久久久久久久久久久| 亚洲伊人色综图| 欧美+日韩+精品| 天天躁夜夜躁狠狠久久av| 欧美国产精品一级二级三级| 亚洲精品国产av成人精品| 新久久久久国产一级毛片| 久久久亚洲精品成人影院| 观看美女的网站| 如何舔出高潮| 男女无遮挡免费网站观看| 这个男人来自地球电影免费观看 | 成人18禁高潮啪啪吃奶动态图| 丝袜在线中文字幕| 成人亚洲欧美一区二区av| 多毛熟女@视频| 下体分泌物呈黄色| av有码第一页| xxx大片免费视频| 成年人午夜在线观看视频| 日本欧美视频一区| 寂寞人妻少妇视频99o| 国产精品秋霞免费鲁丝片| 日韩 亚洲 欧美在线| 宅男免费午夜| 精品人妻偷拍中文字幕| 黄片播放在线免费| 自线自在国产av| 亚洲精品一二三| 女人精品久久久久毛片| 91精品国产国语对白视频| 青春草视频在线免费观看| 日韩一本色道免费dvd| 亚洲精品国产一区二区精华液| av女优亚洲男人天堂| 亚洲精品美女久久av网站| 日韩大片免费观看网站| 女性被躁到高潮视频| 人妻 亚洲 视频| 久久精品亚洲av国产电影网| 久久久精品94久久精品| 国产精品蜜桃在线观看| 国产白丝娇喘喷水9色精品| 日韩精品有码人妻一区| 曰老女人黄片| 国产亚洲av片在线观看秒播厂| 肉色欧美久久久久久久蜜桃| 成人国产麻豆网| 校园人妻丝袜中文字幕| 国产野战对白在线观看| 人人妻人人澡人人爽人人夜夜| 久久99蜜桃精品久久| 多毛熟女@视频| 欧美 亚洲 国产 日韩一| 蜜桃在线观看..| 麻豆av在线久日| 亚洲色图综合在线观看| 男女午夜视频在线观看| 91aial.com中文字幕在线观看| 中文字幕最新亚洲高清| 久久99热这里只频精品6学生| 最近2019中文字幕mv第一页| 人人妻人人澡人人爽人人夜夜| av网站免费在线观看视频| 男的添女的下面高潮视频| 免费久久久久久久精品成人欧美视频| a 毛片基地| 两个人免费观看高清视频| 中文欧美无线码| 王馨瑶露胸无遮挡在线观看| 日本欧美视频一区| 香蕉丝袜av| 亚洲国产欧美在线一区| 性色avwww在线观看| 自线自在国产av| 精品人妻偷拍中文字幕| 国产成人av激情在线播放| 精品国产一区二区三区四区第35| 波多野结衣一区麻豆| 亚洲精品久久午夜乱码| 欧美亚洲日本最大视频资源| 免费在线观看完整版高清| 免费高清在线观看视频在线观看| 在线观看国产h片| 一级黄片播放器| 一级片免费观看大全| 免费在线观看视频国产中文字幕亚洲 | 丰满少妇做爰视频| 国产午夜精品一二区理论片| 国产 一区精品| 亚洲综合精品二区| 18禁国产床啪视频网站| 亚洲伊人久久精品综合| 国产日韩欧美视频二区| 99精国产麻豆久久婷婷| 国产日韩一区二区三区精品不卡| 久久久久国产网址| av在线播放精品| 亚洲中文av在线| 嫩草影院入口| 免费观看av网站的网址| 婷婷色综合www| 色网站视频免费| 亚洲综合精品二区| 纵有疾风起免费观看全集完整版| 97人妻天天添夜夜摸| 亚洲激情五月婷婷啪啪| 亚洲精品在线美女| 精品国产乱码久久久久久小说| 成年av动漫网址| 少妇人妻久久综合中文| 巨乳人妻的诱惑在线观看| 欧美日韩亚洲国产一区二区在线观看 | 夫妻午夜视频| 1024香蕉在线观看| 亚洲色图综合在线观看| 欧美日韩综合久久久久久| 婷婷色综合www| 久久99一区二区三区| 黄色视频在线播放观看不卡| 亚洲精品,欧美精品| 美女视频免费永久观看网站| 又大又黄又爽视频免费| 极品少妇高潮喷水抽搐| 高清av免费在线| 国产午夜精品一二区理论片| 日韩成人av中文字幕在线观看| www.av在线官网国产| 在线观看免费高清a一片| 色网站视频免费| 99久久精品国产国产毛片| 纵有疾风起免费观看全集完整版| 国产乱来视频区| 日本色播在线视频| 精品亚洲乱码少妇综合久久| 成年女人毛片免费观看观看9 | 人妻少妇偷人精品九色| 国产极品天堂在线| 丝袜在线中文字幕| 亚洲综合色惰| 在线观看免费视频网站a站| 中文字幕最新亚洲高清| 黄色一级大片看看| 夜夜骑夜夜射夜夜干| 精品福利永久在线观看| 热re99久久精品国产66热6| 欧美bdsm另类| 国产一区有黄有色的免费视频| 国产男女内射视频| 中文字幕人妻丝袜一区二区 | 又粗又硬又长又爽又黄的视频| 国产色婷婷99| 国产日韩一区二区三区精品不卡| 黄色一级大片看看| 亚洲欧美清纯卡通| 亚洲精品自拍成人| 蜜桃国产av成人99| 亚洲国产av新网站| 久久精品熟女亚洲av麻豆精品| 欧美少妇被猛烈插入视频| 国产成人一区二区在线| 久久午夜福利片| 国产有黄有色有爽视频| 九色亚洲精品在线播放| 日本av免费视频播放| 九九爱精品视频在线观看| 女人精品久久久久毛片| 欧美变态另类bdsm刘玥| 精品少妇久久久久久888优播| 国产精品一区二区在线观看99| 免费在线观看视频国产中文字幕亚洲 | 久久韩国三级中文字幕| 久久人妻熟女aⅴ| 久久精品夜色国产| 免费观看a级毛片全部| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 永久免费av网站大全| 国产亚洲一区二区精品| 极品人妻少妇av视频| 欧美av亚洲av综合av国产av | 日日撸夜夜添| 少妇精品久久久久久久| 制服人妻中文乱码| 亚洲中文av在线| av不卡在线播放| 国产精品国产三级专区第一集| 丝袜在线中文字幕| 免费不卡的大黄色大毛片视频在线观看| 午夜福利一区二区在线看| 日韩成人av中文字幕在线观看| 色网站视频免费| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| 午夜福利视频在线观看免费| 国产高清国产精品国产三级| 人人妻人人澡人人爽人人夜夜| 亚洲一区中文字幕在线| 日本色播在线视频| videos熟女内射| 国产片内射在线| 国产不卡av网站在线观看| 男女边吃奶边做爰视频| 人人妻人人爽人人添夜夜欢视频| 一区二区三区乱码不卡18| 老鸭窝网址在线观看| 国产精品久久久久久久久免| 久久久久国产一级毛片高清牌| 欧美老熟妇乱子伦牲交| 色播在线永久视频| 91午夜精品亚洲一区二区三区| 在线亚洲精品国产二区图片欧美| av网站在线播放免费| 热99久久久久精品小说推荐| 国产人伦9x9x在线观看 | 最近最新中文字幕免费大全7| 亚洲欧洲精品一区二区精品久久久 | 成人黄色视频免费在线看| a级毛片黄视频| 精品国产乱码久久久久久男人| 丰满少妇做爰视频| 亚洲av电影在线进入| 在线观看免费高清a一片| 精品酒店卫生间| 午夜福利视频精品| 久久久国产精品麻豆| 91精品伊人久久大香线蕉| 亚洲av电影在线进入| 咕卡用的链子| 中文字幕另类日韩欧美亚洲嫩草| 夜夜骑夜夜射夜夜干| 成年美女黄网站色视频大全免费| 国产精品久久久久成人av| 精品国产乱码久久久久久男人| 不卡av一区二区三区| 最近中文字幕2019免费版| 性色avwww在线观看| 男女免费视频国产| 成年人午夜在线观看视频| 久久久久久久久久久久大奶| 一级a爱视频在线免费观看| 2022亚洲国产成人精品| 热re99久久精品国产66热6| 亚洲天堂av无毛| 熟妇人妻不卡中文字幕| 极品人妻少妇av视频| 亚洲欧美色中文字幕在线| 久久久久视频综合| 亚洲视频免费观看视频| 午夜老司机福利剧场| freevideosex欧美| 在线观看人妻少妇| 观看美女的网站| 最近最新中文字幕大全免费视频 | 国产精品99久久99久久久不卡 | 欧美日韩国产mv在线观看视频| 2022亚洲国产成人精品| 国产亚洲av片在线观看秒播厂| 国产精品女同一区二区软件| 亚洲第一区二区三区不卡| 国产精品久久久久久久久免| 亚洲精品国产一区二区精华液| 亚洲中文av在线| 免费黄色在线免费观看| 亚洲国产av影院在线观看| 18禁国产床啪视频网站| 99久久精品国产国产毛片| 一区在线观看完整版| 日韩精品免费视频一区二区三区| 国产精品久久久久久av不卡| 2022亚洲国产成人精品| av在线老鸭窝| 亚洲色图综合在线观看| 天天躁夜夜躁狠狠久久av| 日日啪夜夜爽| 天堂8中文在线网| 亚洲国产精品一区二区三区在线| 在线亚洲精品国产二区图片欧美| 欧美另类一区| 国产熟女欧美一区二区| 丝袜在线中文字幕| 男女边吃奶边做爰视频| 精品国产一区二区三区久久久樱花| 久久精品久久久久久久性| 可以免费在线观看a视频的电影网站 | www.av在线官网国产| 最近2019中文字幕mv第一页| 日日撸夜夜添| 91精品三级在线观看| 丝袜脚勾引网站| 国产成人91sexporn| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 中文字幕人妻丝袜一区二区 | 考比视频在线观看| 一本久久精品| 中文字幕人妻熟女乱码| 亚洲人成网站在线观看播放| 最近手机中文字幕大全| 日韩中字成人| 蜜桃在线观看..| 热99国产精品久久久久久7| 最新中文字幕久久久久| 日韩一区二区视频免费看| 亚洲国产欧美网| 如日韩欧美国产精品一区二区三区| 国产又爽黄色视频| 黄色怎么调成土黄色| 中文天堂在线官网| 男女午夜视频在线观看| 最近中文字幕高清免费大全6| 亚洲av日韩在线播放| 欧美成人午夜免费资源| 黑人欧美特级aaaaaa片| av网站免费在线观看视频| 国产野战对白在线观看|