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

    三維地震資料疊前隨機(jī)噪音壓制

    2015-04-24 05:57:30王維紅宋元東林春華劉瑞有姜艷影
    關(guān)鍵詞:同相軸壓制噪音

    王維紅, 宋元東, 林春華, 劉瑞有, 姜艷影

    ( 1. 東北石油大學(xué) 地球科學(xué)學(xué)院,黑龍江 大慶 163318; 2. 大慶油田有限責(zé)任公司 勘探開發(fā)研究院,黑龍江 大慶 163712; 3. 梧州中國(guó)石油昆侖燃?xì)庥邢薰荆瑥V西 梧州 543002; 4. 大慶油田有限責(zé)任公司 第二采油廠,黑龍江 大慶 163414 )

    ?

    三維地震資料疊前隨機(jī)噪音壓制

    王維紅1, 宋元東1, 林春華2, 劉瑞有3, 姜艷影4

    ( 1. 東北石油大學(xué) 地球科學(xué)學(xué)院,黑龍江 大慶 163318; 2. 大慶油田有限責(zé)任公司 勘探開發(fā)研究院,黑龍江 大慶 163712; 3. 梧州中國(guó)石油昆侖燃?xì)庥邢薰荆瑥V西 梧州 543002; 4. 大慶油田有限責(zé)任公司 第二采油廠,黑龍江 大慶 163414 )

    為提高疊前道集數(shù)據(jù)的信噪比,將頻率域基于雙向預(yù)測(cè)誤差濾波的三維疊后隨機(jī)噪音衰減方法(3D RNA),經(jīng)過地震資料的道頭修改和重置后,應(yīng)用于疊前地震資料,實(shí)現(xiàn)三維疊前地震資料隨機(jī)噪音的壓制處理.某工區(qū)CMP道集和CRP道集隨機(jī)噪音壓制的試處理結(jié)果表明:該方法能夠提高三維地震數(shù)據(jù)信噪比,得到的地震道集同相軸更清晰,可在不同的處理流程和步驟中串行應(yīng)用,并且?guī)缀蹩蛇m用于所有三維地震資料.為提高速度分析精度、改善疊加效果、保持地震數(shù)據(jù)的AVO振幅特性和提高成果數(shù)據(jù)體的成像精度提供指導(dǎo).

    疊前; 隨機(jī)噪音; 壓制; 道頭修改; 地震資料

    0 引言

    隨機(jī)噪音在地震資料中普遍存在,主要來源于采集過程和處理過程.在資料采集過程中,隨機(jī)噪音主要來源于地面微震、炮點(diǎn)激發(fā)和儀器接收等因素造成的不規(guī)則干擾[1-2];在資料處理過程中,隨機(jī)噪音主要來源于偏移成像和正演模擬過程,如逆時(shí)偏移中采用隨機(jī)邊界條件引起的隨機(jī)噪音[3-6].在地震記錄上,隨機(jī)噪音沒有一定的視速度和傳播方向,是一種具有各態(tài)歷經(jīng)性質(zhì)的平穩(wěn)隨機(jī)過程[7].隨機(jī)噪音的存在不但引起地震資料信噪比的降低,而且嚴(yán)重制約地震資料后續(xù)的高分辨率處理、目標(biāo)準(zhǔn)確成像、屬性提取和反演等[8],對(duì)當(dāng)前高分辨率和巖性油氣藏勘探有較大影響.在地震資料的“三高”(高信噪比、高分辨率和高保真度)處理過程中,信噪比是基礎(chǔ),是保證地震資料分辨率和保幅的關(guān)鍵;三維地震勘探是地震勘探的主要方法,隨機(jī)噪音的存在嚴(yán)重制約后續(xù)的地震資料處理步驟,如速度分析、AVO分析、疊加效果及偏移成像的效果等.

    疊加是壓制隨機(jī)噪音的有效方法,當(dāng)前大多數(shù)地震處理技術(shù)應(yīng)用于疊前,因此需要在疊前進(jìn)行隨機(jī)噪音壓制,以提高疊前處理的精度,進(jìn)而提高勘探成功率.疊前隨機(jī)噪音壓制的方法主要是濾波法,如中值濾波、f-x域預(yù)測(cè)濾波、多項(xiàng)式擬合、K-L變換和τ-p變換去噪等[9-16].其中中值濾波法計(jì)算簡(jiǎn)單,易于實(shí)現(xiàn),但是平滑處理相當(dāng)于低通濾波,對(duì)有效信號(hào)的高頻成分造成損失,難以滿足巖性勘探的保幅高分辨率處理的要求[9];K-L變換法的基本假設(shè)是水平地層,對(duì)實(shí)際存在的傾斜地層的噪音壓制處理,常導(dǎo)致有效波的振幅受到較大影響[11-12].非局部平均濾波(NLM)方法是近年來發(fā)展起來的隨機(jī)噪音壓制的有效方法[17],對(duì)二維數(shù)據(jù)和疊后地震取得理想的效果;處理三維疊前地震數(shù)據(jù)時(shí),由于三維權(quán)系數(shù)計(jì)算較為困難,并且計(jì)算效率低,所以在三維疊前地震數(shù)據(jù)的隨機(jī)噪音壓制中未得到廣泛應(yīng)用.

    目前廣泛應(yīng)用的是基于疊后地震數(shù)據(jù)的三維隨機(jī)噪音衰減方法(3D RNA).它根據(jù)疊后資料有效波同相軸的相干性[18],應(yīng)用二維預(yù)測(cè)算子實(shí)現(xiàn)隨機(jī)噪音的壓制,對(duì)提高疊后資料信噪比、壓制高頻隨機(jī)噪音非常有效.筆者將三維疊后隨機(jī)噪音衰減技術(shù)應(yīng)用于三維疊前處理,即對(duì)疊前地震數(shù)據(jù)體,通過地震記錄中道頭字修改與重排,形成三維疊后數(shù)據(jù)體的結(jié)構(gòu)形式,應(yīng)用3D RNA方法完成三維疊前隨機(jī)噪音壓制,能夠有效提高疊前資料的信噪比,為三維地震數(shù)據(jù)高精度成像奠定技術(shù)基礎(chǔ).

    1 疊后3D RNA原理

    疊后三維隨機(jī)噪音壓制方法(3D RNA)可以有效壓制隨機(jī)噪音等非相干干擾,提高疊后地震資料的信噪比.在實(shí)際計(jì)算時(shí),需要將三維地震數(shù)據(jù)體進(jìn)行Fourier變換,得到頻率—空間(f-x,y)域的地震數(shù)據(jù).

    三維地震數(shù)據(jù)對(duì)時(shí)間進(jìn)行Fourier正變換的計(jì)算公式為

    (1)

    式中:D(x,y,ω)為三維時(shí)間域地震數(shù)據(jù)d(x,y,t)的Fourier正變換;ω為圓頻率.

    在頻率空間域,對(duì)于給定的任一頻率成分ωi,式(1)所給出的地震數(shù)據(jù)D(x,y,ωi)只是空間位置坐標(biāo)x和y的函數(shù),將其記為f(x,y).設(shè)計(jì)一個(gè)矩形窗,即二維預(yù)測(cè)算子,設(shè)在該矩形窗內(nèi)有m組不同視速度的反射波同相軸Ri(t),各同相軸在兩個(gè)空間方向的反射時(shí)差表示為Δtx(i)和Δty(i),其中i=1,2,…m,則特定頻率成分的二維地震信號(hào)f(x,y)可以表示為二維Z變換的形式:

    (2)

    式中:F(Z1,Z2)為二維信號(hào)f(x,y)的Z變換;R(f)為反射同相軸R(t)的頻率域表示形式.

    將式(2)右端展開,整理可得

    E(Z1,Z2)=F(Z1,Z2)P(Z1,Z2),

    (3)

    式中:E(Z1,Z2)為預(yù)測(cè)誤差的Z變換;P(Z1,Z2)為二維預(yù)測(cè)誤差濾波器;F(Z1,Z2)為輸入的離散信號(hào).式(3)將信號(hào)F(Z1,Z2)應(yīng)用于二維預(yù)測(cè)誤差濾波器P(Z1,Z2),進(jìn)行濾波得到E(Z1,Z2)結(jié)果.

    實(shí)際上,預(yù)測(cè)誤差濾波算子P(Z1,Z2)是未知的,設(shè)三維地震數(shù)據(jù)經(jīng)過矩形預(yù)測(cè)算子預(yù)測(cè)后,輸出的預(yù)測(cè)誤差能量最小,則構(gòu)建目標(biāo)函數(shù),并令目標(biāo)函數(shù)對(duì)預(yù)測(cè)算子的偏導(dǎo)數(shù)為零,得到矩陣方程為

    RH·P=R,

    (4)

    式中:RH為頻率—空間域地震記錄的多道自相關(guān)的Hermite矩陣;R為頻率—空間域地震記錄的多道自相關(guān)列陣;P為二維預(yù)測(cè)算子各分量的列陣.

    根據(jù)式(4)求出對(duì)應(yīng)于不同頻率成分的二維預(yù)測(cè)算子后,對(duì)地震數(shù)據(jù)進(jìn)行維納濾波處理,即可得到頻率域地震信號(hào)的最小平方近似[18-19];再進(jìn)行Fourier反變換,得到時(shí)間空間域隨機(jī)噪音壓制后的三維疊后地震數(shù)據(jù)體.該方法的本質(zhì)是徑向預(yù)測(cè)濾波,利用相鄰地震道有效波的相似性,實(shí)現(xiàn)三維疊后地震數(shù)據(jù)的隨機(jī)噪音壓制處理.

    圖1 三維疊前隨機(jī)噪音壓制道頭修改示意Fig.1 Trace header modification schematic diagram of 3D prestack random noise suppression

    2 疊前3D RNA方法

    一般疊前地震數(shù)據(jù)是以共炮點(diǎn)道集或共中心點(diǎn)道集(CMP)依次存放的,其有效波同相軸的時(shí)距曲線近似為雙曲線.對(duì)疊前地震數(shù)據(jù)隨機(jī)噪音壓制,有效波同相軸具有雙曲性質(zhì),所以不能直接應(yīng)用疊后3D RNA方法.為有效實(shí)現(xiàn)疊前地震數(shù)據(jù)隨機(jī)噪音的壓制,對(duì)疊前地震數(shù)據(jù)的道頭字進(jìn)行修改(見圖1),使它與疊后三維地震數(shù)據(jù)的道頭形式基本一致;修改后的一個(gè)空間坐標(biāo)軸表示炮號(hào)或CMP號(hào)等,另一個(gè)空間坐標(biāo)軸表示道號(hào),修改前后時(shí)間坐標(biāo)軸保持不變(見圖1).由圖1可知,在疊前應(yīng)用3D RNA方法進(jìn)行隨機(jī)噪音衰減,在共炮點(diǎn)道集和共中心點(diǎn)道集可以實(shí)現(xiàn).

    以共炮點(diǎn)道集數(shù)據(jù)為例,說明道頭修改和隨機(jī)噪音壓制的實(shí)現(xiàn)過程.首先將疊前共炮點(diǎn)道集數(shù)據(jù)中單炮記錄作為縱軸,將每炮的地震道作為橫軸,形成一個(gè)平面,每炮中的地震道與按順序排列的炮記錄組成的數(shù)據(jù)體類似于三維疊后數(shù)據(jù)體;然后應(yīng)用疊后3D RNA方法進(jìn)行隨機(jī)噪音壓制處理,再將地震數(shù)據(jù)道頭修改為原來疊前共炮點(diǎn)道集的形式,即實(shí)現(xiàn)疊前地震數(shù)據(jù)的隨機(jī)噪音壓制處理.在疊前隨機(jī)噪音壓制計(jì)算中,只要正確完成道頭字的修改,就可以應(yīng)用疊后3D RNA方法實(shí)現(xiàn)噪音壓制處理;因此提出的疊前隨機(jī)噪音壓制的實(shí)現(xiàn)方法,對(duì)于任何三維疊前地震數(shù)據(jù)體都是適用的.另外,該方法是相對(duì)保幅的隨機(jī)噪音壓制技術(shù),可為后續(xù)的疊前成像算法提供信噪比高的數(shù)據(jù)體.

    3 應(yīng)用效果

    在地震資料處理過程中,疊前在共炮點(diǎn)道集和共檢波點(diǎn)道集處理的步驟較多,一般CMP道集中隨機(jī)噪音壓制對(duì)于提高資料信噪比和速度分析精度等具有重要意義;偏移后CRP道集上,由于存在偏移算法和數(shù)據(jù)不規(guī)則等,也產(chǎn)生隨機(jī)噪音,為提高后續(xù)儲(chǔ)層預(yù)測(cè)、優(yōu)化疊加等處理精度也需要噪音壓制處理.

    大慶油田某工區(qū)三維地震資料于2011年采集,滿覆蓋面積為300 km2,覆蓋次數(shù)為96,面元大小為25 m×25 m,最大偏移距為4 500 m.該工區(qū)勘探目的是斷陷區(qū)深層火山巖和致密砂礫巖儲(chǔ)層,目的層埋藏深、波場(chǎng)復(fù)雜、地震資料信噪比低.為提高地震資料成像精度,對(duì)疊前的CMP道集應(yīng)用提出的疊前3D RNA方法,實(shí)現(xiàn)隨機(jī)噪音壓制處理,提高疊前地震數(shù)據(jù)的信噪比,進(jìn)而使靜校正、速度分析等精度得到提高,為后續(xù)偏移成像精度的提高提供技術(shù)基礎(chǔ).

    對(duì)研究區(qū)的某一CMP道集,應(yīng)用提出的三維疊前隨機(jī)噪音壓制方法得到的噪音壓制前后結(jié)果見圖2,其中動(dòng)校正后含噪音的CMP道集見圖2(a),噪音壓制后的CMP道集見圖2(b).由圖2可知,隨機(jī)噪音壓制后,在1.8~3.0 s之間目的層信噪比得到有效提高,有效波同相軸連續(xù)性增強(qiáng).為更好地說明隨機(jī)噪音壓制后的道集對(duì)后續(xù)處理的影響和數(shù)據(jù)的應(yīng)用,給出該CMP道集位置處疊前隨機(jī)噪音壓制前后的速度譜(見圖3).其中原始含隨機(jī)噪音的CMP道集速度譜見圖3(a),隨機(jī)噪音壓制后的CMP道集速度譜見圖3(b).由圖3可知,經(jīng)過疊前隨機(jī)噪音壓制后,速度譜中能量和聚焦效果得到大幅改善,可以提高速度分析的精度,進(jìn)而提高成像精度.

    圖2 CMP道集疊前隨機(jī)噪音壓制前后結(jié)果Fig.2 The comparison of CMP gathers before and after 3D prestack random noise suppression

    圖3 CMP道集三維疊前隨機(jī)噪音壓制前后速度譜Fig.3 The comparison of velocity spectrum before and after 3D prestack random noise suppression

    CRP道集是疊前時(shí)間或深度偏移后的成果道集,對(duì)CRP道集進(jìn)行疊加,得到偏移的最終成果數(shù)據(jù)體.同時(shí),CRP道集也可用于AVO分析、疊前反演和儲(chǔ)層定量預(yù)測(cè),所以道集質(zhì)量對(duì)成果數(shù)據(jù)體的解釋、應(yīng)用和儲(chǔ)層預(yù)測(cè)等具有決定作用.由于存在偏移計(jì)算等原因,CRP道集含有一定的隨機(jī)噪音.應(yīng)用三維疊前隨機(jī)噪音壓制方法可以提高CRP道集的信噪比.研究區(qū)某一CRP道集隨機(jī)噪音壓制前后結(jié)果見圖4.由圖4可知,噪音壓制后CRP道集同相軸連續(xù)性好、信噪比高,可以有效提高AVO分析和疊前反演的計(jì)算精度.

    研究區(qū)一條主測(cè)線CMP道集疊前噪音壓制前后的疊加結(jié)果見圖5,為便于顯示,抽取其中的一段進(jìn)行對(duì)比.其中常規(guī)處理得到的疊加剖面見圖5(a),對(duì)CMP道集進(jìn)行三維疊前隨機(jī)噪音壓制后的疊加剖面見圖5(b).由圖5可知,噪音壓制后地震數(shù)據(jù)體信噪比有一定程度提高,同相軸連續(xù)性增強(qiáng),能夠?yàn)楹罄m(xù)處理提供較好的數(shù)據(jù)基礎(chǔ).

    圖4 CRP道集疊前隨機(jī)噪音壓制前后結(jié)果Fig.4 The comparison of CRP gathers before and after prestack random noise suppression

    圖5 CMP道集疊前噪音壓制前后疊加數(shù)據(jù)結(jié)果Fig.5 The comparison of stack sections before and after prestack random noise suppression in CMP domain

    4 結(jié)論

    (1)基于疊后三維地震數(shù)據(jù)隨機(jī)噪音壓制方法,修改疊前數(shù)據(jù)的道頭字,形成三維地震數(shù)據(jù)的疊前噪音壓制方法.疊前隨機(jī)噪音壓制方法可以在不同數(shù)據(jù)道集中實(shí)現(xiàn)隨機(jī)噪音壓制,對(duì)于CMP道集,除提高道集的信噪比外,還可提高速度分析精度,為后續(xù)基于CMP的資料處理提供基礎(chǔ)數(shù)據(jù).

    (2)在CRP道集上進(jìn)行疊前隨機(jī)噪音壓制,可以有效提高地震資料的信噪比,為最終成像疊加及成像數(shù)據(jù)體的解釋、反演和應(yīng)用提供基礎(chǔ),同時(shí)也可為AVO分析提供高精度道集數(shù)據(jù).

    (3)該方法思路簡(jiǎn)單,計(jì)算效率較高,易于實(shí)現(xiàn),可以有效提高地震資料的信噪比和成像精度,適合在低信噪比資料的地區(qū)推廣應(yīng)用.

    [1] Ristau J P, Wooil M M. Adaptive filtering of random noise in 2-D geophysical data [J]. Geophysics, 2001,66:342-349.

    [2] 胡天躍.地震資料疊前去噪技術(shù)的現(xiàn)狀與未來[J].地球物理學(xué)進(jìn)展,2002,17(2):218-223. Hu Tianyue. The current situation and future of Seismic data prestack noise attenuation techniques [J]. Progress in Geophysics, 2002,17(2):218-223.

    [3] 王維紅,郭雪豹,石穎.非局部平均濾波噪聲壓制方法及其在VSP資料逆時(shí)偏移中的應(yīng)用[J].石油物探,2015,54(2):165-171. Wang Weihong , Guo Xuebao, Shi Ying. Nonlocal means filtering denoising approach and its application in VSP reverse time migration seismic data [J]. Geophysical Prospecting for Petroleum, 2015,54(2):165-171.

    [4] 郭雪豹,王建民,王維紅,等.基于GPU并行加速的VSP數(shù)據(jù)逆時(shí)偏移[J].東北石油大學(xué)學(xué)報(bào),2014,38(2):58-62. Guo Xuebao, Wang Jianmin, Wang Weihong, et al. VSP data reverse-time migration based on GPU parallel acceleration [J]. Journal of Northeast Petroleum University, 2014,38(2):58-62.

    [5] 石穎,陸加敏,柯璇,等.基于GPU并行加速的疊前逆時(shí)偏移方法研究[J].東北石油大學(xué)學(xué)報(bào),2012,36(4):111-115. Shi Ying, Lu Jiamin, Ke Xuan, et al. Prestack reverse time migration based on GPU parallel accelerating algorithm [J]. Journal of Northeast Petroleum University, 2012,36(4):111-115.

    [6] 田東升,王云專,李義鵬,等.單程和雙程波動(dòng)方程疊前深度偏移方法[J].東北石油大學(xué)學(xué)報(bào),2014,38(4):39-44. Tian Dongsheng, Wang Yunzhuan, Li Yipeng, et al. One-way and two-way wave equation pre-stack depth migration approaches [J]. Journal of Northeast Petroleum University, 2014,38(4):39-44.

    [7] 陸基孟,王永剛.地震勘探原理[M].東營(yíng):中國(guó)石油大學(xué)出版社,1993:420-428. Lu Jimeng, Wang Yonggang. The principle of seismic exploration [M]. Dongying: China University of Petroleum Press, 1993:420-428.

    [8] 熊翥.地震數(shù)據(jù)處理應(yīng)用技術(shù)[M].北京:石油工業(yè)出版社,1993:141-224. Xiong Zhu. Seismic data processing technique application [M]. Beijing: Petroleum Industry Press, 1993:141-224.

    [9] Luis L C. Random noise reduction [C]. SEG Technical Program Expanded Abstracts, 1984:525-527.

    [10] Yu S, Cai X, Su Y. Seismic signal enhancement by polynomial fitting [J]. Applied Geophysics, 1989,1:57-65.

    [11] Jones I F, Levy S. Signal-to-noise ratio enhancement in multichannel seismic data via the Karhunen-Loeve transform [J]. Geophysical Prospecting, 1987,35:12-32.

    [12] 劉洪林,張春堂,朱秋影,等.K-L變換在地震資料去噪中的應(yīng)用[J].大慶石油學(xué)院學(xué)報(bào),2007,31(4):19-21. Liu Honglin, Zhang Chuntang, Zhu Qiuying, et al. Application of K-L transformation to the noise attenuation in geological data [J]. Journal of Daqing Petroleum Institute, 2007,31(4):19-21.

    [13] 石穎,李瑩,王維紅,等.線性Radon變換噪音壓制法及其在古龍斷陷中的應(yīng)用[J].東北石油大學(xué)學(xué)報(bào),2012,36(4):116-120. Shi Ying, Li Ying, Wang Weihong, et al. Approach of linear noise suppression using radon transform and its application in Gulong fault depression [J]. Journal of Northeast Petroleum University, 2012,36(4):116-120.

    [14] Al-Yahya K M. Application of the partial Karhunen-Loeve transform to suppress random noise in seismic sections [J]. Geophysical Prospecting, 1993,39:77-93.

    [15] 陳遵德,段天友,朱廣生,等.SVD濾波方法的改進(jìn)及應(yīng)用[J].石油地球物理勘探,1994,29(6):783-792. Chen Zunde, Duan Tianyou, Zhu Guangsheng, et al. Improvement of singular-value decomposition filtering and the application [J]. OGP, 1994,29(6):783-792.

    [16] Alan R M, Panos G K. Efficient tau-p hyperbolic velocity filtering [J]. Geophysics, 1990,55(5):619-625.

    [17] Bonar D, Sacchi M. Denoising seismic data using the nonlocal means algorithm [J]. Geophysics, 2012,77(1):A5-A8.

    [18] 李國(guó)發(fā).全三維隨機(jī)噪聲衰減技術(shù)[J].石油地球物理勘探,1995,30(3):310-318. Li Guofa. Full 3-D random-noise attenuation [J]. OGP, 1995,30(3):310-318.

    [19] 國(guó)九英,周興元,楊慧珠.三維f-x,y域隨機(jī)噪音衰減[J].石油地球物理勘探,1995,30(2):207-215. Guo Jiuying, Zhou Xingyuan, Yang Huizhu. Attenuation of random noise in (f-x,y) domain [J]. OGP, 1995,30(2):207-215.

    2015-04-20;編輯:任志平

    國(guó)家自然科學(xué)基金項(xiàng)目(41474118);國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2012AA061202)

    王維紅(1975-),男,博士,副教授,主要從事地震資料數(shù)字處理方面的研究.

    TE132.1,P631

    A

    2095-4107(2015)03-0025-07

    DOI 10.3969/j.issn.2095-4107.2015.03.004

    猜你喜歡
    同相軸壓制噪音
    噪音,總是有噪音!
    無法逃避的噪音
    虛同相軸方法及其在陸上地震層間多次波壓制中的應(yīng)用
    一種新型無人機(jī)數(shù)據(jù)鏈抗壓制干擾技術(shù)的研究
    空射誘餌在防空壓制電子戰(zhàn)中的應(yīng)用
    無人機(jī)(2018年1期)2018-07-05 09:51:02
    一種改進(jìn)的相關(guān)法自動(dòng)拾取同相軸
    噪音的小把戲
    白噪音的三種用法
    Coco薇(2017年9期)2017-09-07 22:09:28
    一種反射同相軸自動(dòng)拾取算法
    一種舊物品擠壓成型機(jī)
    科技資訊(2016年12期)2016-05-30 05:07:58
    久久久久国产精品人妻一区二区| 国产精品av久久久久免费| 电影成人av| 操美女的视频在线观看| 亚洲国产欧美网| 在线观看www视频免费| 在线观看一区二区三区激情| 免费观看av网站的网址| 日韩中文字幕视频在线看片| 美国免费a级毛片| 在线亚洲精品国产二区图片欧美| 老汉色av国产亚洲站长工具| 别揉我奶头~嗯~啊~动态视频 | 欧美变态另类bdsm刘玥| 一区福利在线观看| 另类亚洲欧美激情| 久久亚洲国产成人精品v| 国产一卡二卡三卡精品 | 国产精品一国产av| bbb黄色大片| √禁漫天堂资源中文www| 在线观看国产h片| 看免费av毛片| 男男h啪啪无遮挡| 不卡av一区二区三区| 一级a爱视频在线免费观看| 99热全是精品| 悠悠久久av| 免费少妇av软件| 午夜福利,免费看| 久久人妻熟女aⅴ| 操美女的视频在线观看| 2021少妇久久久久久久久久久| 久久精品国产a三级三级三级| 国产欧美亚洲国产| 日本wwww免费看| 黄色毛片三级朝国网站| 亚洲少妇的诱惑av| 亚洲欧美一区二区三区黑人| www.精华液| 美女国产高潮福利片在线看| 80岁老熟妇乱子伦牲交| 十八禁人妻一区二区| 伦理电影免费视频| 亚洲av综合色区一区| 亚洲国产精品999| 日韩大码丰满熟妇| xxxhd国产人妻xxx| 悠悠久久av| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 久久韩国三级中文字幕| 免费黄网站久久成人精品| 国产一区二区激情短视频 | 亚洲av国产av综合av卡| 我的亚洲天堂| 自线自在国产av| 日韩熟女老妇一区二区性免费视频| 看非洲黑人一级黄片| 久久97久久精品| 女的被弄到高潮叫床怎么办| 国产免费现黄频在线看| 国产一区亚洲一区在线观看| 久久精品亚洲熟妇少妇任你| 日韩欧美一区视频在线观看| 亚洲精品,欧美精品| 欧美日韩视频高清一区二区三区二| 天堂俺去俺来也www色官网| 欧美激情 高清一区二区三区| 国产精品国产三级专区第一集| 亚洲av电影在线观看一区二区三区| 久久久久久久国产电影| 女人久久www免费人成看片| 成人国产麻豆网| 精品少妇内射三级| 巨乳人妻的诱惑在线观看| 大话2 男鬼变身卡| 久久午夜综合久久蜜桃| 亚洲精品中文字幕在线视频| 午夜福利在线免费观看网站| 成人三级做爰电影| 日韩中文字幕欧美一区二区 | 日本色播在线视频| 一区二区三区四区激情视频| 午夜福利视频在线观看免费| 一级片免费观看大全| 亚洲精品国产色婷婷电影| 大片免费播放器 马上看| 欧美日韩亚洲综合一区二区三区_| av有码第一页| 亚洲一码二码三码区别大吗| 亚洲精品自拍成人| 欧美日韩一级在线毛片| 亚洲国产欧美在线一区| 亚洲精品久久午夜乱码| www日本在线高清视频| 亚洲国产欧美日韩在线播放| 美女中出高潮动态图| 国产成人精品福利久久| 中文天堂在线官网| 亚洲美女视频黄频| 精品少妇久久久久久888优播| 亚洲成人手机| 国产淫语在线视频| 哪个播放器可以免费观看大片| 啦啦啦视频在线资源免费观看| 99久久人妻综合| 香蕉丝袜av| 精品午夜福利在线看| 黄色视频在线播放观看不卡| 亚洲成av片中文字幕在线观看| 欧美日韩一级在线毛片| 国产精品熟女久久久久浪| avwww免费| 亚洲精品自拍成人| 午夜日韩欧美国产| 免费少妇av软件| 亚洲激情五月婷婷啪啪| 少妇猛男粗大的猛烈进出视频| 新久久久久国产一级毛片| 男女免费视频国产| 国产日韩一区二区三区精品不卡| 国产日韩欧美亚洲二区| 日韩大码丰满熟妇| 精品久久久精品久久久| 黄片小视频在线播放| 日本色播在线视频| 亚洲一区中文字幕在线| 日韩,欧美,国产一区二区三区| 精品久久久久久电影网| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 中文字幕色久视频| 国产无遮挡羞羞视频在线观看| 大香蕉久久网| 夫妻午夜视频| 国产亚洲一区二区精品| 香蕉丝袜av| 免费av中文字幕在线| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载| www.熟女人妻精品国产| 老司机靠b影院| 国产精品一区二区在线观看99| 中文字幕色久视频| 亚洲国产欧美日韩在线播放| 欧美精品一区二区大全| av一本久久久久| 国产成人欧美| 欧美在线一区亚洲| 国产成人av激情在线播放| 亚洲国产欧美日韩在线播放| 尾随美女入室| 美国免费a级毛片| 国产精品免费视频内射| 亚洲专区中文字幕在线 | 超色免费av| 可以免费在线观看a视频的电影网站 | 老司机亚洲免费影院| 国产在线一区二区三区精| 在线精品无人区一区二区三| 日韩av在线免费看完整版不卡| 国产精品人妻久久久影院| 精品国产国语对白av| 欧美亚洲日本最大视频资源| 国产成人一区二区在线| 又大又爽又粗| 国产亚洲最大av| 日韩 亚洲 欧美在线| 国产一区有黄有色的免费视频| 国产精品 国内视频| 免费观看性生交大片5| 欧美日韩av久久| 高清av免费在线| 国产av国产精品国产| 国产精品熟女久久久久浪| 男女高潮啪啪啪动态图| 免费观看性生交大片5| 狂野欧美激情性bbbbbb| 青春草视频在线免费观看| 你懂的网址亚洲精品在线观看| 午夜激情久久久久久久| 亚洲国产最新在线播放| 在线观看一区二区三区激情| 麻豆av在线久日| 国产视频首页在线观看| 国产亚洲av高清不卡| 日韩精品免费视频一区二区三区| 男女之事视频高清在线观看 | 在线精品无人区一区二区三| 国产国语露脸激情在线看| 毛片一级片免费看久久久久| 日本91视频免费播放| 叶爱在线成人免费视频播放| 亚洲国产av新网站| 国产免费视频播放在线视频| 婷婷成人精品国产| a级片在线免费高清观看视频| 在线观看免费日韩欧美大片| 久久综合国产亚洲精品| 热re99久久精品国产66热6| 精品人妻熟女毛片av久久网站| 午夜av观看不卡| 91老司机精品| 999精品在线视频| 国产免费现黄频在线看| 精品一区二区免费观看| 一级毛片我不卡| 黑人巨大精品欧美一区二区蜜桃| 99精国产麻豆久久婷婷| 在线观看免费日韩欧美大片| 嫩草影院入口| 免费少妇av软件| 这个男人来自地球电影免费观看 | 精品视频人人做人人爽| 黑丝袜美女国产一区| av免费观看日本| 久久久久精品国产欧美久久久 | 久久99一区二区三区| av网站在线播放免费| 成年人午夜在线观看视频| 国产免费又黄又爽又色| 天天添夜夜摸| 日韩精品有码人妻一区| 国产av国产精品国产| 女人高潮潮喷娇喘18禁视频| 久久精品国产亚洲av涩爱| 亚洲精品日韩在线中文字幕| 久久久国产一区二区| 久久精品国产综合久久久| 极品少妇高潮喷水抽搐| 精品少妇黑人巨大在线播放| 免费观看a级毛片全部| 久久午夜综合久久蜜桃| 最近的中文字幕免费完整| 国产不卡av网站在线观看| 18在线观看网站| 两个人看的免费小视频| 性色av一级| 亚洲精品美女久久久久99蜜臀 | 国产一区二区在线观看av| 国产免费一区二区三区四区乱码| 乱人伦中国视频| 精品国产乱码久久久久久男人| 一区二区av电影网| 久久精品国产综合久久久| 国产精品麻豆人妻色哟哟久久| 狠狠婷婷综合久久久久久88av| 成人三级做爰电影| 久久精品久久久久久久性| 国产探花极品一区二区| 丝袜在线中文字幕| 午夜福利免费观看在线| 我的亚洲天堂| 熟女av电影| 免费看不卡的av| 国产精品偷伦视频观看了| 国产精品熟女久久久久浪| 五月开心婷婷网| 久久久久精品国产欧美久久久 | 我要看黄色一级片免费的| 老汉色av国产亚洲站长工具| 国产一区二区 视频在线| 国产日韩欧美在线精品| 亚洲免费av在线视频| 国产无遮挡羞羞视频在线观看| 人成视频在线观看免费观看| 夜夜骑夜夜射夜夜干| 日韩,欧美,国产一区二区三区| 黑人猛操日本美女一级片| 亚洲av在线观看美女高潮| 51午夜福利影视在线观看| 久久99精品国语久久久| 肉色欧美久久久久久久蜜桃| 久久久精品国产亚洲av高清涩受| 亚洲国产欧美网| 在线观看免费视频网站a站| 大香蕉久久成人网| 午夜日韩欧美国产| 在线观看人妻少妇| 国产免费福利视频在线观看| 高清av免费在线| 黄网站色视频无遮挡免费观看| 欧美激情 高清一区二区三区| 亚洲天堂av无毛| 中文字幕高清在线视频| 波多野结衣av一区二区av| 国产福利在线免费观看视频| 国产精品三级大全| 婷婷成人精品国产| 久久国产精品大桥未久av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品成人av观看孕妇| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| a 毛片基地| 激情视频va一区二区三区| 黄色视频不卡| 亚洲成国产人片在线观看| 亚洲欧美日韩另类电影网站| 欧美日本中文国产一区发布| 日韩制服骚丝袜av| 亚洲四区av| 亚洲国产精品999| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久精品电影小说| 夫妻性生交免费视频一级片| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| av在线播放精品| 国产伦理片在线播放av一区| 国产精品三级大全| 日本欧美视频一区| 成人三级做爰电影| 美女中出高潮动态图| 一本—道久久a久久精品蜜桃钙片| 自线自在国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久久久精品古装| 亚洲精品视频女| 五月天丁香电影| 美女脱内裤让男人舔精品视频| 热re99久久国产66热| 久久鲁丝午夜福利片| 久久午夜综合久久蜜桃| 国产精品久久久av美女十八| 亚洲精品,欧美精品| 9色porny在线观看| 精品午夜福利在线看| 曰老女人黄片| 国产一区二区激情短视频 | 在线观看一区二区三区激情| 香蕉国产在线看| 亚洲四区av| 久久精品国产亚洲av涩爱| 午夜激情久久久久久久| 天天躁夜夜躁狠狠躁躁| 日韩av不卡免费在线播放| 超碰成人久久| 久久久国产精品麻豆| 国产片内射在线| videos熟女内射| 男人操女人黄网站| 老司机影院成人| 亚洲久久久国产精品| 啦啦啦在线观看免费高清www| 亚洲专区中文字幕在线 | 国产乱来视频区| 国产一区二区激情短视频 | 国产福利在线免费观看视频| 国产一区二区 视频在线| 伊人久久国产一区二区| 91aial.com中文字幕在线观看| 国产1区2区3区精品| 亚洲国产成人一精品久久久| 51午夜福利影视在线观看| 久久久久精品久久久久真实原创| 黄色怎么调成土黄色| 99久久人妻综合| 国产女主播在线喷水免费视频网站| 精品少妇久久久久久888优播| av在线老鸭窝| 香蕉国产在线看| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av涩爱| 777米奇影视久久| 在线亚洲精品国产二区图片欧美| 纵有疾风起免费观看全集完整版| 啦啦啦中文免费视频观看日本| 国产精品香港三级国产av潘金莲 | 亚洲国产av新网站| 日韩中文字幕欧美一区二区 | 国产精品国产三级专区第一集| 亚洲欧美一区二区三区久久| 免费观看人在逋| 婷婷色综合www| 一级毛片黄色毛片免费观看视频| 2021少妇久久久久久久久久久| 精品免费久久久久久久清纯 | 青春草视频在线免费观看| 在线观看三级黄色| 久热这里只有精品99| 熟女少妇亚洲综合色aaa.| 日韩 欧美 亚洲 中文字幕| av网站在线播放免费| 欧美日韩亚洲高清精品| 在线精品无人区一区二区三| 亚洲人成网站在线观看播放| 青草久久国产| 又黄又粗又硬又大视频| 深夜精品福利| 日本一区二区免费在线视频| 亚洲欧美精品综合一区二区三区| 国产成人精品福利久久| svipshipincom国产片| 精品亚洲乱码少妇综合久久| 国产精品成人在线| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 高清av免费在线| 在线亚洲精品国产二区图片欧美| 男女之事视频高清在线观看 | 80岁老熟妇乱子伦牲交| 日本色播在线视频| 中文乱码字字幕精品一区二区三区| 精品少妇黑人巨大在线播放| 国产一区二区在线观看av| 一级黄片播放器| 亚洲av日韩精品久久久久久密 | 亚洲一码二码三码区别大吗| 国产成人av激情在线播放| 国产又爽黄色视频| 又大又爽又粗| 一本大道久久a久久精品| 亚洲精品美女久久av网站| 亚洲国产中文字幕在线视频| 亚洲少妇的诱惑av| av免费观看日本| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 99久久99久久久精品蜜桃| 久久久久精品人妻al黑| 亚洲精品国产色婷婷电影| 精品国产一区二区三区四区第35| 精品一区二区免费观看| 老司机靠b影院| 最新的欧美精品一区二区| 老司机影院成人| 国产男女超爽视频在线观看| 久久久久久久精品精品| 人妻 亚洲 视频| 中国三级夫妇交换| 亚洲av成人不卡在线观看播放网 | 欧美在线黄色| 在线观看三级黄色| 欧美老熟妇乱子伦牲交| 美女福利国产在线| 满18在线观看网站| 色婷婷久久久亚洲欧美| 波野结衣二区三区在线| 欧美精品一区二区大全| svipshipincom国产片| 午夜福利视频在线观看免费| 一区在线观看完整版| 免费黄频网站在线观看国产| 国产精品成人在线| 精品亚洲成国产av| 国产福利在线免费观看视频| av国产精品久久久久影院| 亚洲成人国产一区在线观看 | 99久久精品国产亚洲精品| av在线老鸭窝| 亚洲精品国产色婷婷电影| 国产探花极品一区二区| 女人爽到高潮嗷嗷叫在线视频| 男女下面插进去视频免费观看| 亚洲精品成人av观看孕妇| 69精品国产乱码久久久| 天天躁狠狠躁夜夜躁狠狠躁| √禁漫天堂资源中文www| 亚洲久久久国产精品| 午夜免费鲁丝| 亚洲精品一二三| 亚洲一区二区三区欧美精品| 男女边吃奶边做爰视频| 亚洲第一av免费看| 亚洲成国产人片在线观看| 精品人妻一区二区三区麻豆| 成人黄色视频免费在线看| 性色av一级| 一区二区三区精品91| 麻豆精品久久久久久蜜桃| 亚洲精品视频女| 国产黄色视频一区二区在线观看| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说| 国产一区二区激情短视频 | 午夜福利在线免费观看网站| 黄色怎么调成土黄色| 欧美变态另类bdsm刘玥| 成年女人毛片免费观看观看9 | 国产熟女午夜一区二区三区| 韩国高清视频一区二区三区| 又大又爽又粗| 一级a爱视频在线免费观看| 人成视频在线观看免费观看| 中文乱码字字幕精品一区二区三区| 久久人妻熟女aⅴ| 中文字幕精品免费在线观看视频| 欧美日韩综合久久久久久| 国产99久久九九免费精品| 亚洲精品美女久久久久99蜜臀 | 精品卡一卡二卡四卡免费| 性高湖久久久久久久久免费观看| 亚洲精品av麻豆狂野| 捣出白浆h1v1| 精品少妇久久久久久888优播| 亚洲精品美女久久久久99蜜臀 | 亚洲国产最新在线播放| 人体艺术视频欧美日本| 欧美人与性动交α欧美软件| 18在线观看网站| 韩国精品一区二区三区| 国产精品 欧美亚洲| 国产在线免费精品| 亚洲国产精品成人久久小说| 99香蕉大伊视频| 国产日韩欧美在线精品| 黄色怎么调成土黄色| 肉色欧美久久久久久久蜜桃| 在线观看www视频免费| 桃花免费在线播放| 欧美激情 高清一区二区三区| 色播在线永久视频| 午夜日本视频在线| 亚洲av成人不卡在线观看播放网 | 可以免费在线观看a视频的电影网站 | 亚洲熟女精品中文字幕| 一级爰片在线观看| 久久久久久久精品精品| 久久99一区二区三区| 亚洲国产欧美日韩在线播放| 午夜激情久久久久久久| 国产在线免费精品| 午夜激情av网站| 久久久久精品久久久久真实原创| 欧美黑人精品巨大| 丝袜美足系列| 男女边摸边吃奶| 丝瓜视频免费看黄片| 黑人猛操日本美女一级片| videosex国产| 精品少妇一区二区三区视频日本电影 | 免费观看性生交大片5| 久久精品国产亚洲av高清一级| 超碰97精品在线观看| 国产精品人妻久久久影院| 国产探花极品一区二区| 亚洲一区中文字幕在线| 国产极品粉嫩免费观看在线| avwww免费| 久久国产亚洲av麻豆专区| 电影成人av| 国产深夜福利视频在线观看| 亚洲精品国产色婷婷电影| 大陆偷拍与自拍| 日韩中文字幕欧美一区二区 | 欧美国产精品va在线观看不卡| 久久性视频一级片| 一级爰片在线观看| 亚洲,一卡二卡三卡| 亚洲精品aⅴ在线观看| 三上悠亚av全集在线观看| 欧美中文综合在线视频| 国产精品一区二区在线观看99| 亚洲五月色婷婷综合| 极品人妻少妇av视频| 久久亚洲国产成人精品v| 在线观看免费高清a一片| 国产高清不卡午夜福利| 欧美国产精品va在线观看不卡| 黄色 视频免费看| 国产精品成人在线| 欧美成人精品欧美一级黄| 亚洲美女视频黄频| 国产野战对白在线观看| 久久久国产一区二区| 亚洲精华国产精华液的使用体验| 国精品久久久久久国模美| a级毛片在线看网站| 日日啪夜夜爽| 国产成人av激情在线播放| 满18在线观看网站| 自线自在国产av| 在线观看国产h片| av.在线天堂| 18禁动态无遮挡网站| 国产精品三级大全| 99国产综合亚洲精品| 国产精品久久久久久精品古装| 国产亚洲欧美精品永久| 女人久久www免费人成看片| 一级a爱视频在线免费观看| 一级毛片 在线播放| 精品人妻熟女毛片av久久网站| 9热在线视频观看99| 精品一区二区三区四区五区乱码 | 赤兔流量卡办理| 婷婷色综合www| 看非洲黑人一级黄片| 免费看不卡的av| 免费在线观看视频国产中文字幕亚洲 | 欧美黑人欧美精品刺激| 中文乱码字字幕精品一区二区三区| 亚洲av成人精品一二三区| 少妇被粗大的猛进出69影院| 夜夜骑夜夜射夜夜干| 操出白浆在线播放| 国产不卡av网站在线观看| 久久精品亚洲熟妇少妇任你| 日韩精品有码人妻一区| 青春草视频在线免费观看| 天堂8中文在线网| 国产熟女欧美一区二区| 国产精品99久久99久久久不卡 | 国产亚洲av高清不卡| av免费观看日本| 色视频在线一区二区三区| 人人妻人人澡人人爽人人夜夜| 老司机影院成人| 黄片无遮挡物在线观看| 精品亚洲乱码少妇综合久久| 国产精品免费大片| 男女床上黄色一级片免费看| 亚洲国产最新在线播放|