• <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
    国产久久久一区二区三区| 激情在线观看视频在线高清| 日本黄大片高清| 日韩成人在线观看一区二区三区| 99热只有精品国产| 夜夜爽天天搞| 首页视频小说图片口味搜索| 国产探花在线观看一区二区| 黄色成人免费大全| 伊人久久大香线蕉亚洲五| 亚洲片人在线观看| 日本免费a在线| 可以在线观看的亚洲视频| 成熟少妇高潮喷水视频| 亚洲国产日韩欧美精品在线观看 | xxx96com| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩精品青青久久久久久| 午夜免费观看网址| 大型av网站在线播放| 久久久久性生活片| 一本精品99久久精品77| 好看av亚洲va欧美ⅴa在| 色哟哟哟哟哟哟| 国产主播在线观看一区二区| 91av网站免费观看| 午夜精品久久久久久毛片777| 国内久久婷婷六月综合欲色啪| 特级一级黄色大片| 亚洲精品一区av在线观看| 久久精品国产亚洲av香蕉五月| 欧美性猛交╳xxx乱大交人| 男人的好看免费观看在线视频 | 国产成年人精品一区二区| 一个人观看的视频www高清免费观看 | 老鸭窝网址在线观看| e午夜精品久久久久久久| 日日爽夜夜爽网站| 色哟哟哟哟哟哟| 国产精品一及| 麻豆av在线久日| 成人三级做爰电影| 国产伦在线观看视频一区| 91字幕亚洲| www.999成人在线观看| 午夜福利在线观看吧| 一卡2卡三卡四卡精品乱码亚洲| 无限看片的www在线观看| 97碰自拍视频| 国产亚洲欧美在线一区二区| 久久香蕉激情| 啦啦啦观看免费观看视频高清| 一边摸一边做爽爽视频免费| 久久这里只有精品中国| 岛国在线免费视频观看| 国产v大片淫在线免费观看| 搡老熟女国产l中国老女人| 精品电影一区二区在线| 久久精品影院6| 成人手机av| 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区免费观看 | 脱女人内裤的视频| 不卡一级毛片| 中文字幕人妻丝袜一区二区| 亚洲成av人片免费观看| 淫妇啪啪啪对白视频| 亚洲中文字幕日韩| 18禁国产床啪视频网站| 嫁个100分男人电影在线观看| 国产久久久一区二区三区| 手机成人av网站| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 精品国产亚洲在线| 免费在线观看完整版高清| 这个男人来自地球电影免费观看| 在线播放国产精品三级| 变态另类成人亚洲欧美熟女| 他把我摸到了高潮在线观看| 日韩大码丰满熟妇| 国产又黄又爽又无遮挡在线| 国产亚洲精品久久久久久毛片| 男女做爰动态图高潮gif福利片| 最近最新中文字幕大全免费视频| 成熟少妇高潮喷水视频| 婷婷亚洲欧美| 草草在线视频免费看| 国产成人精品无人区| 精品日产1卡2卡| 国产免费男女视频| 国产真人三级小视频在线观看| 搡老妇女老女人老熟妇| 精品欧美一区二区三区在线| 色老头精品视频在线观看| av国产免费在线观看| 亚洲av第一区精品v没综合| 又爽又黄无遮挡网站| 老司机深夜福利视频在线观看| 亚洲av成人不卡在线观看播放网| 午夜福利免费观看在线| 又爽又黄无遮挡网站| 欧美绝顶高潮抽搐喷水| 色在线成人网| 国产高清有码在线观看视频 | 此物有八面人人有两片| 成在线人永久免费视频| 中出人妻视频一区二区| 成人国产综合亚洲| 亚洲片人在线观看| 日韩高清综合在线| 日本一本二区三区精品| 老汉色∧v一级毛片| АⅤ资源中文在线天堂| 三级毛片av免费| 9191精品国产免费久久| 美女大奶头视频| 久久精品国产亚洲av高清一级| 成人亚洲精品av一区二区| 亚洲午夜精品一区,二区,三区| 男人舔女人下体高潮全视频| 日韩高清综合在线| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看 | 国产精品一区二区三区四区免费观看 | 亚洲欧美日韩东京热| 一级a爱片免费观看的视频| 成人午夜高清在线视频| 老司机福利观看| 亚洲午夜精品一区,二区,三区| 两个人免费观看高清视频| 中亚洲国语对白在线视频| 精品第一国产精品| 日本免费一区二区三区高清不卡| 国产精品av久久久久免费| 法律面前人人平等表现在哪些方面| 久热爱精品视频在线9| 色综合亚洲欧美另类图片| 又大又爽又粗| 日韩大尺度精品在线看网址| 嫩草影视91久久| 久久久久久国产a免费观看| 精品一区二区三区四区五区乱码| 欧美 亚洲 国产 日韩一| 五月伊人婷婷丁香| 国产成人精品久久二区二区91| 欧洲精品卡2卡3卡4卡5卡区| 妹子高潮喷水视频| 久久久久久国产a免费观看| 国产视频内射| 激情在线观看视频在线高清| 美女高潮喷水抽搐中文字幕| 无人区码免费观看不卡| 在线十欧美十亚洲十日本专区| 久久精品国产亚洲av香蕉五月| 桃红色精品国产亚洲av| 天堂av国产一区二区熟女人妻 | 久久九九热精品免费| 窝窝影院91人妻| 国产精品影院久久| 大型av网站在线播放| 国产又黄又爽又无遮挡在线| www.自偷自拍.com| 岛国在线免费视频观看| 黑人欧美特级aaaaaa片| 99久久久亚洲精品蜜臀av| 国产精品爽爽va在线观看网站| 性欧美人与动物交配| 免费在线观看日本一区| 麻豆久久精品国产亚洲av| 露出奶头的视频| 99久久精品国产亚洲精品| 亚洲国产精品合色在线| 18禁国产床啪视频网站| 精品无人区乱码1区二区| 欧美 亚洲 国产 日韩一| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 在线国产一区二区在线| 一区福利在线观看| 欧美日韩福利视频一区二区| 少妇熟女aⅴ在线视频| 国产不卡一卡二| 男人的好看免费观看在线视频 | 国产精品久久久久久人妻精品电影| 成人18禁在线播放| 日本免费一区二区三区高清不卡| 大型av网站在线播放| 亚洲av成人精品一区久久| 欧美+亚洲+日韩+国产| 亚洲av熟女| 国产在线精品亚洲第一网站| 国产av在哪里看| 啪啪无遮挡十八禁网站| 精品国内亚洲2022精品成人| 天天躁狠狠躁夜夜躁狠狠躁| 俺也久久电影网| 一区二区三区国产精品乱码| 久久亚洲精品不卡| 日韩有码中文字幕| 亚洲av五月六月丁香网| 国产黄a三级三级三级人| 成人国语在线视频| 精品国产超薄肉色丝袜足j| 中文字幕熟女人妻在线| 丰满人妻熟妇乱又伦精品不卡| 老司机靠b影院| 国产伦人伦偷精品视频| 亚洲欧美日韩高清在线视频| 欧美在线黄色| 99精品久久久久人妻精品| 久久精品亚洲精品国产色婷小说| 欧美午夜高清在线| 欧美色欧美亚洲另类二区| av免费在线观看网站| 麻豆国产97在线/欧美 | 男女午夜视频在线观看| 高清毛片免费观看视频网站| 好看av亚洲va欧美ⅴa在| 国产精品一及| 国产成人欧美在线观看| 精品一区二区三区视频在线观看免费| 国产精品精品国产色婷婷| 欧美国产日韩亚洲一区| 国产男靠女视频免费网站| 这个男人来自地球电影免费观看| 国产午夜精品论理片| 18禁国产床啪视频网站| 精品久久蜜臀av无| 亚洲成人久久性| а√天堂www在线а√下载| 最新美女视频免费是黄的| www.自偷自拍.com| 国产成人影院久久av| 99久久99久久久精品蜜桃| 黄色女人牲交| 亚洲成a人片在线一区二区| 女同久久另类99精品国产91| 久久中文字幕人妻熟女| 精品国产乱子伦一区二区三区| 国产亚洲精品综合一区在线观看 | 香蕉国产在线看| 日日爽夜夜爽网站| av欧美777| 免费在线观看日本一区| 亚洲色图 男人天堂 中文字幕| 在线播放国产精品三级| 国产成人系列免费观看| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 69av精品久久久久久| 精品久久久久久久人妻蜜臀av| av天堂在线播放| 久久性视频一级片| 国语自产精品视频在线第100页| 久久久久久大精品| 国产精品一区二区免费欧美| 国产精品香港三级国产av潘金莲| aaaaa片日本免费| 国产午夜福利久久久久久| 在线观看美女被高潮喷水网站 | 大型黄色视频在线免费观看| 一级片免费观看大全| 真人做人爱边吃奶动态| 精品乱码久久久久久99久播| 校园春色视频在线观看| 欧美日本亚洲视频在线播放| 精品熟女少妇八av免费久了| 一级黄色大片毛片| 制服丝袜大香蕉在线| 一级毛片女人18水好多| 黄色女人牲交| 成人三级做爰电影| 亚洲美女黄片视频| 亚洲人与动物交配视频| 日韩中文字幕欧美一区二区| 欧美成人性av电影在线观看| 亚洲国产日韩欧美精品在线观看 | 国产精品香港三级国产av潘金莲| 欧美色欧美亚洲另类二区| 久久婷婷成人综合色麻豆| 最近最新免费中文字幕在线| 岛国在线观看网站| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 午夜福利在线在线| 日本一区二区免费在线视频| 久久婷婷成人综合色麻豆| 久久精品91无色码中文字幕| 在线十欧美十亚洲十日本专区| a级毛片a级免费在线| 欧美一级a爱片免费观看看 | 波多野结衣巨乳人妻| 99热只有精品国产| 日韩欧美国产在线观看| 欧美黑人精品巨大| 欧美极品一区二区三区四区| 午夜福利18| 在线观看日韩欧美| 动漫黄色视频在线观看| 这个男人来自地球电影免费观看| 国产黄片美女视频| 性色av乱码一区二区三区2| 午夜日韩欧美国产| 亚洲中文av在线| 精品一区二区三区av网在线观看| 久久天堂一区二区三区四区| 亚洲国产精品久久男人天堂| 午夜a级毛片| 欧美3d第一页| 婷婷精品国产亚洲av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 中文亚洲av片在线观看爽| 欧美zozozo另类| 国产精品亚洲美女久久久| 久久久国产精品麻豆| 脱女人内裤的视频| 长腿黑丝高跟| 精品国内亚洲2022精品成人| 天堂√8在线中文| 国产aⅴ精品一区二区三区波| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 国产野战对白在线观看| 欧美色欧美亚洲另类二区| 天天躁夜夜躁狠狠躁躁| 国产成人aa在线观看| 亚洲精品中文字幕在线视频| 成年人黄色毛片网站| 在线a可以看的网站| 91av网站免费观看| 亚洲国产欧美网| 19禁男女啪啪无遮挡网站| 久久久精品大字幕| 国产精品一区二区三区四区久久| 亚洲av成人精品一区久久| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 欧美在线黄色| 欧美极品一区二区三区四区| 久久婷婷成人综合色麻豆| 国产黄色小视频在线观看| 搡老熟女国产l中国老女人| 91成年电影在线观看| 一级毛片精品| av有码第一页| 日韩欧美在线乱码| 国产精品自产拍在线观看55亚洲| 香蕉久久夜色| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 母亲3免费完整高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久国产精品人妻蜜桃| 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3| 久久这里只有精品中国| 国产高清videossex| 婷婷六月久久综合丁香| 成人三级做爰电影| 男女做爰动态图高潮gif福利片| 99热这里只有精品一区 | 亚洲国产欧洲综合997久久,| 欧美最黄视频在线播放免费| 国产亚洲欧美在线一区二区| 亚洲成av人片在线播放无| 91成年电影在线观看| 色综合婷婷激情| 露出奶头的视频| 国产成人精品久久二区二区免费| 91成年电影在线观看| 中文字幕av在线有码专区| videosex国产| 好看av亚洲va欧美ⅴa在| 国内毛片毛片毛片毛片毛片| 老司机深夜福利视频在线观看| 日本免费a在线| 不卡av一区二区三区| 51午夜福利影视在线观看| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 国产伦在线观看视频一区| 少妇裸体淫交视频免费看高清 | 国产亚洲欧美在线一区二区| 国产亚洲精品av在线| www国产在线视频色| 国产成人欧美在线观看| av天堂在线播放| 99久久精品国产亚洲精品| 亚洲真实伦在线观看| 欧美3d第一页| 亚洲av成人不卡在线观看播放网| 性欧美人与动物交配| 亚洲成人精品中文字幕电影| 亚洲av电影不卡..在线观看| 无遮挡黄片免费观看| 无人区码免费观看不卡| 国产99白浆流出| 久久久国产成人免费| 91字幕亚洲| 啦啦啦免费观看视频1| 精品午夜福利视频在线观看一区| 成人三级黄色视频| 国产久久久一区二区三区| 色综合亚洲欧美另类图片| 观看免费一级毛片| 亚洲精品国产精品久久久不卡| 麻豆一二三区av精品| 亚洲人成77777在线视频| 亚洲国产日韩欧美精品在线观看 | 日韩欧美一区二区三区在线观看| 99热这里只有是精品50| 视频区欧美日本亚洲| 亚洲成av人片在线播放无| 少妇裸体淫交视频免费看高清 | 亚洲熟妇熟女久久| 91在线观看av| 大型av网站在线播放| av福利片在线观看| 色精品久久人妻99蜜桃| 精品国产超薄肉色丝袜足j| 成人av在线播放网站| 91字幕亚洲| 亚洲成a人片在线一区二区| 丁香欧美五月| 欧美另类亚洲清纯唯美| 亚洲欧美日韩无卡精品| 欧美中文综合在线视频| 99精品欧美一区二区三区四区| 在线播放国产精品三级| 天堂√8在线中文| 色综合婷婷激情| 亚洲av成人av| 色综合婷婷激情| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 制服丝袜大香蕉在线| 成人国产综合亚洲| 日本a在线网址| av免费在线观看网站| 久久九九热精品免费| 日日干狠狠操夜夜爽| 欧美日本视频| 亚洲av电影在线进入| 欧美绝顶高潮抽搐喷水| 亚洲七黄色美女视频| 国产亚洲精品一区二区www| 欧美成人性av电影在线观看| 日韩大尺度精品在线看网址| a级毛片a级免费在线| 国产精品乱码一区二三区的特点| 久久久久国内视频| 麻豆成人av在线观看| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 可以在线观看的亚洲视频| 日本五十路高清| 色综合欧美亚洲国产小说| 国产三级在线视频| √禁漫天堂资源中文www| 亚洲av美国av| 国内久久婷婷六月综合欲色啪| 天堂动漫精品| 欧美日本视频| 亚洲九九香蕉| 亚洲一区二区三区不卡视频| 国产精品av久久久久免费| 757午夜福利合集在线观看| 十八禁人妻一区二区| 国产乱人伦免费视频| 久久亚洲真实| 国产精品综合久久久久久久免费| 桃色一区二区三区在线观看| 精品人妻1区二区| 国产激情久久老熟女| 男女床上黄色一级片免费看| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 国产精品国产高清国产av| 欧美zozozo另类| 亚洲精品在线美女| 精品一区二区三区四区五区乱码| 人妻夜夜爽99麻豆av| 99久久精品国产亚洲精品| avwww免费| 欧美久久黑人一区二区| 亚洲成人国产一区在线观看| 国产亚洲精品综合一区在线观看 | 亚洲一区二区三区不卡视频| 免费看十八禁软件| 成人国产一区最新在线观看| 99精品久久久久人妻精品| 亚洲国产精品sss在线观看| 国产高清videossex| 久久精品国产清高在天天线| 欧美高清成人免费视频www| 久久久久久久精品吃奶| 欧美性长视频在线观看| 亚洲黑人精品在线| 午夜视频精品福利| 成年女人毛片免费观看观看9| 日韩免费av在线播放| 精品日产1卡2卡| 五月玫瑰六月丁香| 国产成+人综合+亚洲专区| 一夜夜www| 亚洲精品在线美女| 成人国语在线视频| 亚洲国产精品sss在线观看| 国产aⅴ精品一区二区三区波| 淫秽高清视频在线观看| 色综合婷婷激情| a级毛片在线看网站| 18禁观看日本| 亚洲第一欧美日韩一区二区三区| 欧美黑人精品巨大| 麻豆一二三区av精品| 人成视频在线观看免费观看| 一二三四在线观看免费中文在| 国产成人aa在线观看| 日韩精品免费视频一区二区三区| 国产精品久久电影中文字幕| 久久99热这里只有精品18| 波多野结衣巨乳人妻| 母亲3免费完整高清在线观看| 韩国av一区二区三区四区| 一级片免费观看大全| 男人舔女人的私密视频| 99热这里只有精品一区 | 精品一区二区三区视频在线观看免费| 国产区一区二久久| 国产精品一及| 亚洲av片天天在线观看| 午夜亚洲福利在线播放| 日本黄大片高清| 久久久久国产一级毛片高清牌| 午夜a级毛片| 看片在线看免费视频| a级毛片在线看网站| 在线十欧美十亚洲十日本专区| 俺也久久电影网| 亚洲乱码一区二区免费版| 中文字幕最新亚洲高清| 免费电影在线观看免费观看| 国产伦在线观看视频一区| 禁无遮挡网站| 久久婷婷成人综合色麻豆| 色综合欧美亚洲国产小说| 日韩三级视频一区二区三区| 老司机在亚洲福利影院| www.www免费av| 亚洲乱码一区二区免费版| 久久久久久久久免费视频了| 免费观看人在逋| 岛国在线观看网站| 日韩有码中文字幕| 久久99热这里只有精品18| 亚洲av片天天在线观看| 国产成人精品久久二区二区91| 午夜视频精品福利| svipshipincom国产片| 美女大奶头视频| 亚洲欧美日韩高清在线视频| 99国产精品一区二区蜜桃av| 日韩免费av在线播放| 国内精品久久久久久久电影| 久久精品国产清高在天天线| 一本大道久久a久久精品| 香蕉国产在线看| 日韩精品青青久久久久久| 欧美大码av| 男女之事视频高清在线观看| 人妻夜夜爽99麻豆av| e午夜精品久久久久久久| 男女床上黄色一级片免费看| 亚洲一区中文字幕在线| 可以在线观看的亚洲视频| 无遮挡黄片免费观看| 看免费av毛片| 亚洲精华国产精华精| 99国产极品粉嫩在线观看| 人妻久久中文字幕网| 亚洲中文av在线| 国产视频一区二区在线看| 国产精品亚洲av一区麻豆| 草草在线视频免费看| 久久久久免费精品人妻一区二区| 青草久久国产| 成人av一区二区三区在线看| 国产精品永久免费网站| 俄罗斯特黄特色一大片| 老司机靠b影院| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧洲综合997久久,| 精品福利观看| xxxwww97欧美| 欧美+亚洲+日韩+国产| 日韩欧美国产在线观看| 一个人免费在线观看电影 | 欧美午夜高清在线| 久久精品亚洲精品国产色婷小说| 最近最新中文字幕大全免费视频| www.自偷自拍.com| 国产精品一区二区免费欧美| 亚洲黑人精品在线| 老汉色∧v一级毛片| 成人18禁在线播放| 国产精品爽爽va在线观看网站| 搡老熟女国产l中国老女人| 久久久国产成人免费| 免费av毛片视频| 天堂av国产一区二区熟女人妻 | 国产亚洲av嫩草精品影院|