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

    Kirchhoff疊前時(shí)間偏移關(guān)鍵參數(shù)分析與研究

    2011-10-17 08:40:28王鐸翰閻建國(guó)
    物探化探計(jì)算技術(shù) 2011年2期
    關(guān)鍵詞:射線孔徑傾角

    顧 雯,王鐸翰,閻建國(guó)

    (1.成都理工大學(xué) 信息工程學(xué)院,四川成都610059;2.中石油川慶鉆探工程有限公司 地球物理勘探公司,四川成都610213;3.成都理工大學(xué) 地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,四川成都610059)

    Kirchhoff疊前時(shí)間偏移關(guān)鍵參數(shù)分析與研究

    顧 雯1,王鐸翰2,閻建國(guó)3

    (1.成都理工大學(xué) 信息工程學(xué)院,四川成都610059;2.中石油川慶鉆探工程有限公司 地球物理勘探公司,四川成都610213;3.成都理工大學(xué) 地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,四川成都610059)

    疊前時(shí)間偏移已經(jīng)成為地震資料處理的常規(guī)手段之一,在構(gòu)造復(fù)雜但速度橫向變化不大時(shí),應(yīng)用克?;舴虔B前時(shí)間偏移方法,一般可以得到較好的成像效果。但是疊前時(shí)間偏移是一項(xiàng)系統(tǒng)工程,一些配套技術(shù)的選取和一些關(guān)鍵參數(shù)的選取,直接影響到成像結(jié)果。這里根據(jù)OMEGA2處理系統(tǒng)制定的疊前時(shí)間偏移的處理方法和流程,對(duì)疊前時(shí)間偏移處理的配套技術(shù)和關(guān)鍵參數(shù)(疊前道集偏移距分組、偏移孔徑和傾角選取、均方根速度建模、旅行時(shí)算法)進(jìn)行了分析總結(jié),結(jié)合試驗(yàn)程序進(jìn)行了研究和討論,并提出了關(guān)鍵參數(shù)選取的基本原則,對(duì)實(shí)際處理工作具有一定的應(yīng)用價(jià)值。

    克?;舴虔B前時(shí)間偏移;偏移成像;關(guān)鍵參數(shù)

    0 前言

    隨著勘探與開發(fā)要求的不斷提高,地震勘探越來越關(guān)注復(fù)雜構(gòu)造的精確成像。影響復(fù)雜介質(zhì)成像的主要原因在于:地層傾角大,橫向速度變化劇烈,斷層發(fā)育,埋藏較深等方面。地震偏移技術(shù)正是為了滿足這一需求而不斷發(fā)展和完善起來的。地下介質(zhì)的偏移歸位處理,既可以在疊加之前進(jìn)行,也可在疊加之后進(jìn)行。前者叫疊前偏移,后者叫疊后偏移。按其實(shí)施的領(lǐng)域,還可分為時(shí)間偏移和深度偏移。

    在四川地區(qū),工區(qū)地貌為典型的丘陵和高山地貌。由于地面海拔相對(duì)高差大,地勢(shì)險(xiǎn)惡,溝壑縱橫,地形切割厲害,縱向、橫向速度變化大,所以使用疊后偏移難以得到理想的成像效果,影響地質(zhì)任務(wù)的完成。作者在本文中,結(jié)合四川工區(qū)實(shí)際地震資料,主要圍繞影響Kirchhoff疊前時(shí)間偏移效果的關(guān)鍵參數(shù)進(jìn)行分析探討,總結(jié)了影響其成像效果的原因,實(shí)現(xiàn)了試射法旅行時(shí)計(jì)算程序,研究了不同旅行時(shí)計(jì)算方法的應(yīng)用效果,最后提出了關(guān)鍵參數(shù)選取的基本原則,這對(duì)于實(shí)際資料處理具有一定的參考價(jià)值。

    1 Kirchhoff疊前時(shí)間偏移原理和方法

    1.1 偏移的概念

    地震波偏移成像即通過波場(chǎng)外推,將炮檢點(diǎn)逐漸向繞射點(diǎn)延拓。其中,繞射曲線會(huì)逐漸收斂,最終波場(chǎng)外推到繞射點(diǎn)的深度,繞射曲線收斂到繞射點(diǎn)上,其本質(zhì)就是把數(shù)據(jù)域中的繞射波,收斂到模型空間中,并產(chǎn)生繞射的空間點(diǎn)。簡(jiǎn)而言之,地震偏移成像就是一個(gè)收斂繞射波的過程[2]。

    1.2 Kirchhoff繞射積分中三個(gè)因子的含義

    作者基于二維處理的Kirchhoff的積分,在二維情況下,可將積分公式簡(jiǎn)化為式(1)[1]。

    式(1)中結(jié)合了傾斜因子cosφ,該因子表示為傳播方向與垂直軸z之間夾角的余弦[4];球面擴(kuò)散因子在二維波動(dòng)空間中用表示[4];子波整形因子利用計(jì)算波場(chǎng)的半導(dǎo)數(shù),如式(1)中項(xiàng),這同45°相移和振幅譜用頻率的校正是等價(jià)的(此時(shí)的速度V為均方根速度)[4]。

    (1)傾斜因子。點(diǎn)孔(反射或繞射)得到的波前振幅與角度有關(guān),所以在求和前要考慮到傾斜因素構(gòu)成的振幅依賴角度的關(guān)系。

    (3)子波整形因子是關(guān)于惠更斯二次震源波形的相干性?;莞苟握鹪醋鳛樽硬憫?yīng),它在雙曲線軌跡各點(diǎn)必然具有一定的相位和頻率特性,否則,當(dāng)它們互相靠近時(shí),就不會(huì)有振幅相消的現(xiàn)象。對(duì)2D偏移,設(shè)計(jì)一個(gè)45°常相位譜,振幅譜正比于頻率平方根;對(duì)3D偏移,這個(gè)因子的相移為90°,振幅譜與頻率成正比[4]。

    1.3 Kirchhoff繞射偏移中的孔徑概念

    Kirchhoff偏移是沿著繞射雙曲線作振幅疊加。在理論上,一條繞射雙曲線在時(shí)間和距離上可以無(wú)限延伸,而實(shí)際中,我們只能截取有限段作為求和路徑,截取空間范圍,即實(shí)際求和路徑寬度,稱為偏移孔徑[4],如圖1所示。

    圖1 虛線范圍內(nèi)為計(jì)算的偏移孔徑Fig.1 The diagram showing the migration aperture selection:dotted lines for the calculation of the offset within the aperture

    繞射雙曲線的曲率由速度函數(shù)決定。低速雙曲線的孔徑與高速雙曲線的孔徑相比窄一些,這一點(diǎn)符合我們的直覺。即速度高意味著偏移量大,而一般高速度層伴隨著的是地層深度向下縱伸。

    2 Kirchhoff偏移處理

    2.1 資料處理流程

    此次地震數(shù)據(jù)來自四川某工區(qū)2008年AAA三維地震資料,我們選取第20線(主測(cè)線)為例,結(jié)合工區(qū)地貌為典型的丘陵和高山地貌,地面海拔相對(duì)高差大,溝壑縱橫,地形切割厲害,縱向、橫向速度變化大的特點(diǎn),采用疊前偏移處理,取得了較好的成像。

    從圖2可以看出,在疊前時(shí)間偏移處理時(shí),速度模型需要反復(fù)修改,它是影響資料處理進(jìn)度的關(guān)鍵。

    圖2 疊前時(shí)間偏移流程圖Fig.2 Flow chart of pre-stack timemigration

    將08AAA20地震測(cè)線數(shù)據(jù)進(jìn)行預(yù)處理,可以得到水平疊加剖面(見圖3)。由圖3中可見,剖面的橫向分辨率十分模糊。

    圖3 08AAA20線的水平疊加剖面Fig.3 Stack profile of08AAA20 line

    在CMP點(diǎn)從1 700和2 100附近可以看見層位有明顯的交叉,這是背斜構(gòu)造的明顯標(biāo)志。

    圖3的疊加剖面中背斜構(gòu)造變寬,但同向軸在結(jié)合部發(fā)生了交叉,不能準(zhǔn)確反映地下構(gòu)造情況。

    2.2 Kirchhoff偏移處理成像結(jié)果

    在Kirchhoff偏移模塊選取偏移距分組、偏移孔徑和傾角選取、旅行時(shí)算法的參數(shù)進(jìn)行計(jì)算,可得到如圖4所示的偏移結(jié)果。

    圖4 08AAA20線偏移后成圖Fig.4 Migration profile of 08AAA20 line

    與圖5(a)的水平疊加剖面對(duì)比可明顯看出,在給定參數(shù)下的Kirchhoff疊前時(shí)間偏移的效果。散射點(diǎn)、回轉(zhuǎn)波得到收斂,地層背斜明顯,偏移后橫向成像效果得到明顯提高。

    3 Kirchhoff偏移參數(shù)分析

    3.1 偏移距分組

    偏移距分組即是實(shí)現(xiàn)等效偏移距偏移(EOM)成像。將克希霍夫時(shí)間偏移的雙平方根方程(DSR),轉(zhuǎn)化為一個(gè)單平方根方程,使輸入采樣在沒有發(fā)生時(shí)移的情況下,映射到以等效偏移距為變量的共散射點(diǎn)道集(CSP)上,使散射能量按照雙曲線規(guī)律分布,之后在CSP道集上應(yīng)用Kirchhoff動(dòng)校正并疊加。

    分好組后,不同的偏移距分組通過公式(單平方根方程),可以計(jì)算出各自所對(duì)應(yīng)的共散射點(diǎn)道集(CSP),對(duì)各自同一CMP道集內(nèi)具有相同偏移距的道進(jìn)行疊加。

    通過分貝譜(見下頁(yè)圖6)和對(duì)數(shù)能量譜(見下頁(yè)圖7)可以看出,不同分組下各自剖面的差異。其中,采用五組的分組方式,在分貝譜和對(duì)數(shù)能量譜的高頻部份,會(huì)出現(xiàn)明顯的噪音干擾。

    3.2 偏移孔徑和傾角

    偏移孔徑和傾角選取是影響Kirchhoff偏移的重要參數(shù)。偏移孔徑的寬度是用所占的地震道數(shù)來度量的[4]。傾角是在確定合理孔徑后,對(duì)其相應(yīng)時(shí)間上的偏移孔徑,利用傾角進(jìn)一步控制孔徑,其目的是提高偏移質(zhì)量和節(jié)約運(yùn)行時(shí)間。

    在08AAA20研究工區(qū)上,結(jié)合偏移孔徑和傾角參數(shù)的范圍變化,分別選取了范圍在1 000 m、2 500 m、4 500 m、6 000 m下的偏移剖面(見下頁(yè)圖8)。

    偏移孔徑和傾角參數(shù)的范圍變化在1 000 m時(shí),偏移剖面的橫向分辨率十分低,背斜邊界陡傾角不能正確成像,繞射沒有很好的歸位,剖面面貌差;在2 000 m時(shí)背斜成像變得清晰,但是深部的噪音影響依舊嚴(yán)重;隨著范圍的不斷增大,在4 500 m和6 000 m時(shí),剖面整體成像效果很好。

    圖5 偏移前后構(gòu)造部位對(duì)比Fig.5 Comparison of the construction site before and aftermigration

    圖6 5組、100組、30組的分貝譜Fig.6 dB spectrum of 5 groups,100 groups and 30 groups

    圖7 100組、30組、5組的對(duì)數(shù)能量譜Fig.7 The logarithmic energy spectrum of 100 groups,30 groups and 5 groups

    圖8 不同偏移孔徑下的偏移成像Fig.8 Themigrated images under differentmigration aperturs

    3.3 偏移速度的選取

    速度是偏移計(jì)算中的又一重要參數(shù),之前進(jìn)行的偏移計(jì)算,都是基于偏移距分組之前的速度分析所提供的速度進(jìn)行計(jì)算。在利用Kirchhoff疊前時(shí)間偏移得到時(shí)間域CRP道集后,通過反動(dòng)校正交互拾取速度譜,可建立均方根速度模型。只有在滿足該點(diǎn)的均方根速度時(shí),才會(huì)出現(xiàn)拉平。否則,只會(huì)出現(xiàn)上翹(速度過小)或下彎(速度過大)[8]。

    圖9是將拾取的速度進(jìn)行0.8倍、1.2倍的速度變化,對(duì)剖面進(jìn)行偏移成像??梢钥吹剑S著速度的增加,背斜構(gòu)造頂端的同向軸由0.8倍圖像中的下彎,到1.2倍圖像中同向軸上翹。

    圖9 不同速度場(chǎng)下的剖面面貌Fig.9 The profiles with different velocity fields

    在實(shí)際運(yùn)算中,由于選取的是彎曲射線旅行時(shí),采用的是層狀速度模型,上面均方根速度的提取,是為了提取層速度參與運(yùn)算(是在整體速度選好后,對(duì)每個(gè)速度解釋點(diǎn)進(jìn)行的0.8倍、1.2倍的速度掃描)。

    3.4 旅行時(shí)的算法

    在時(shí)間偏移假設(shè)在水平層狀介質(zhì)的基礎(chǔ)上,Kirchhoff疊前時(shí)間偏移利用均方根速度模型計(jì)算旅行時(shí)。而在實(shí)際工區(qū)的研究中,由于選取的Kirchhoff彎曲射線疊前時(shí)間偏移,采用的是層狀速度模型。上面提到的均方根速度的提取,通過變換求出層速度來參與運(yùn)算。

    在生產(chǎn)中,主要采用直射線法進(jìn)行旅行時(shí)間的計(jì)算。式(2)是彎曲射線偏移的旅行時(shí)間公式[6]。

    直射線法采用了式(2)的前二項(xiàng)進(jìn)行計(jì)算,在高階的計(jì)算精度上不如彎曲射線法。

    圖10 水平層狀介質(zhì)下兩種算法的計(jì)算路徑Fig.10 The path of two algorithms for horizontal layered media

    試射法通過循環(huán)迭代的方式,統(tǒng)計(jì)每一層的旅行時(shí),較之直射線法和彎曲射線法在旅行時(shí)的計(jì)算上更為準(zhǔn)確。我們編制了試射法的C語(yǔ)言程序,并將它用到了疊前時(shí)間偏移的旅行時(shí)計(jì)算中,取得了一定效果。下面將討論三種算法在水平層狀模型中的旅行時(shí)間差。

    假設(shè)一個(gè)四層模型,速度是V1=2 000 m/s、V2=3 000 m/s、V3=4 000 m/s,到每層底部界面的雙程旅行時(shí)是1 s、2 s、3 s、4 s,代入偏移距范圍是500 m~12 000 m。

    采用試射法計(jì)算、彎曲射線法進(jìn)行旅行時(shí)間計(jì)算,通過模型正演的結(jié)果見表1

    通過表1和圖11可以看出,隨著偏移距的增大,彎曲射線旅行時(shí)計(jì)算的旅行時(shí),與試射法的旅行時(shí)差值越來越大。

    在模型的對(duì)比中也發(fā)現(xiàn),當(dāng)選取偏移距與埋深相差不大時(shí),其計(jì)算的結(jié)果也是相差無(wú)幾的。只有在偏移距遠(yuǎn)大于目的層埋深時(shí),采用試射法較之另外兩種方法更為準(zhǔn)確。

    對(duì)比圖13和圖14的剖面,兩種方法在大偏移距下,對(duì)剖面構(gòu)造主體(構(gòu)造相對(duì)復(fù)雜)的成像效果影響更為明顯。

    表1 不同旅行時(shí)間算法的計(jì)算結(jié)果Tab.1 Computation results of different travel time computingmethods

    圖11 對(duì)比試射法和彎曲射線在不同偏移距下的旅行時(shí)間差Fig.11 Comparison of testmethod and bending rays travel under different offset time difference

    圖12 對(duì)比直射線法和彎曲射線在不同偏移距下的旅行時(shí)間差Fig.12 Comparison of direct linemethod and bending rays travel under different offset time difference

    在表2中我們可以看見,彎曲射線法比直射線法采用更高階的運(yùn)算,在偏移距較小時(shí),直射線法與彎曲射線法計(jì)算的旅行時(shí)間差可以是忽略不計(jì)的;而在偏移距變大時(shí),彎曲射線法的高階項(xiàng)對(duì)曲線精度的計(jì)算也沒有太大的實(shí)際價(jià)值。因此,在實(shí)際中人們采用彎曲射線法時(shí)不愿意采用高階項(xiàng),而通常把偏移距x控制在很小的數(shù)值,以消除計(jì)算上的影響(與偏移距分組相聯(lián)系)。

    圖13 試射法在偏移距7 000 m的剖面Fig.13 The section of testmethod in the 7000m cross-section offset

    圖14 彎曲射線法在偏移距7 000 m的剖面Fig.14 The section of bending raymethod in the 7000mcross-section offset

    表2 小偏移距各種旅行算法的時(shí)間差Tab.2 Time difference of travel time algorithm in a variety small-offset

    在實(shí)際運(yùn)用中,應(yīng)結(jié)合工區(qū)實(shí)際要求情況來選取適合的旅行時(shí)間算法。

    3.5 關(guān)鍵參數(shù)的選取原則

    (1)偏移距分組原則:在實(shí)際的生產(chǎn)過程中,同一偏移距組的空道,應(yīng)控制在30%以內(nèi)。對(duì)所缺的偏移距補(bǔ)充零道,使每一個(gè)CMP道集內(nèi)偏移距分布盡量均勻、齊全。此外,若選用覆蓋次數(shù)低的小偏移距道集,偏移處理后就會(huì)填充上假的大偏移距數(shù)據(jù),這樣的道集求取的速度會(huì)影響之后速度拾取時(shí)的精確度[10]。

    (2)偏移孔徑選取原則:當(dāng)偏移孔徑取值小時(shí),由于不能含蓋成像帶范圍,將損失掉陡傾角信息,偏移剖面會(huì)出現(xiàn)畸變現(xiàn)象。偏移孔徑越大,由于含有大量來自其它反射點(diǎn)的波場(chǎng),偏移剖面噪音會(huì)加重,降低了剖面信噪比,還增加了機(jī)器的運(yùn)算時(shí)間,因此可以先根據(jù)目的層的深度選取偏移孔徑大小。在首次選取偏移孔徑大小時(shí),應(yīng)與目的層深度相同,保證目的層位信息量的獲取,再進(jìn)行實(shí)驗(yàn)驗(yàn)證,選取合適的偏移孔徑大小。

    (3)速度的選取原則:經(jīng)預(yù)處理后,疊加剖面速度的選取,基本符合地下地質(zhì)情況,偏移歸位使速度譜能量更加收斂。從能量譜入手,不斷細(xì)化速度的選取,建立更加準(zhǔn)確的速度模型。

    (4)旅行時(shí)間計(jì)算方法的選取原則:試射法作為一種高級(jí)旅行時(shí)間算法,在精度上比其它二種算法更符合地下實(shí)際地層速度分層。在地下地質(zhì)結(jié)構(gòu)相對(duì)簡(jiǎn)單,偏移距不大的情況下,算法對(duì)于剖面成像影響不大。但在地下地質(zhì)結(jié)構(gòu)復(fù)雜,采用大偏移距接收地層信息時(shí),采用試射法更符合地下地層的信息的求取。

    4 結(jié)束語(yǔ)

    通過參與四川工區(qū)08AAA20線地震資料疊前時(shí)間偏移處理工作,得出以下幾點(diǎn)體會(huì)和認(rèn)識(shí):

    (1)疊前時(shí)間偏移并不是一個(gè)獨(dú)立的處理過程,其成像效果與基礎(chǔ)處理密不可分。疊前時(shí)間偏移要求輸入經(jīng)過基礎(chǔ)處理的高質(zhì)量道集數(shù)據(jù),數(shù)據(jù)具有較高的信噪比,振幅能量均勻,且不存在靜校正影響。因此在基礎(chǔ)處理中,應(yīng)著重做好振幅恢復(fù)與補(bǔ)償、疊前去噪、子波處理、速度分析與靜校正迭代等工作,為疊前時(shí)間偏移打下良好的基礎(chǔ)[3]。

    (2)速度場(chǎng)是直接影響偏移結(jié)果的關(guān)鍵,合理的、接近地質(zhì)真實(shí)的、高精度的速度模型,是疊前時(shí)間偏移獲得高質(zhì)量效果的前提。疊前時(shí)間偏移除了其算法本身外,剩下的工作主要就是建立速度模型,而建立速度模型是一個(gè)復(fù)雜而費(fèi)時(shí)的迭代過程[9]。

    (3)對(duì)偏移孔徑的實(shí)驗(yàn),可以在偏移質(zhì)量和偏移效率中尋找到最佳的結(jié)合點(diǎn)。

    (4)在信噪比問題上,Kirchhoff積分法有其固有理論上的缺陷:即原始地震道上的任何一個(gè)樣點(diǎn)幅值,不管它是信號(hào)還是噪聲,都參與運(yùn)算。有效波信號(hào)符合反射波傳播路徑規(guī)律,會(huì)被正確地偏移歸位;而某些幅值較大的噪聲經(jīng)偏移算子改造后,放到偏移孔徑范圍內(nèi)的所有網(wǎng)格點(diǎn)上,無(wú)法相互抵消,形成畫弧現(xiàn)象,降低了疊前時(shí)間偏移的成像質(zhì)量[7]。

    [1]賀振華.反射地震資料偏移處理與反演方法[M].重慶:重慶大學(xué)出版社,1989.

    [2]路基孟.地震勘探原理下冊(cè)[M].北京:中國(guó)石油大學(xué)出版,2006.

    [3]張穎.三維地震疊前時(shí)間偏移處理技術(shù)應(yīng)用與展望[J].石油地球物理勘探,2006,36(05):174.

    [4]渥·伊爾馬茲.地震資料分析[M].劉懷山,王克斌,童思友,等譯.北京:石油工業(yè)出版社,2006.

    [5]EKOROMARIUSBLANCHARD,GU HANMING.Prestack Kirchhoff Time Migration for Imaging of Complex Structure in the Nanyang Oilfield[J].The leading edge,2007,4(2):101.

    [6]P.赫布拉爾Th.克雷.根據(jù)地震反射波旅行時(shí)測(cè)定——計(jì)算層速度[M].北京:石油工業(yè)出版社,1987年.

    [7]WANG DI,WANG HUAZHONGl,MA ZAITIANL,et al.Review of prestack time migration methods[J].Geophysics,1984,49(11):2064.

    [8]徐基祥,崔化娟.疊前時(shí)間偏移解析速度分析[J].石油地球物理勘探,2004,39(2):126.

    [9]李建明,楊品榮.地震資料疊前時(shí)間偏移處理技巧與效果——以坦南斷裂帶為例[J].石油地球物理進(jìn)展,2002,07(02):38.

    [10]麻三懷,楊長(zhǎng)春,孫福利,等.克希霍夫疊前時(shí)間偏移技術(shù)在復(fù)雜構(gòu)造帶地震資料處理中的應(yīng)用[J].石油物探,2008,54(07):447.

    P 631.4

    A

    1001—1749(2011)02—0140—07

    2010-09-01改回日期:2010-12-29

    顧雯(1986-),女,四川成都人,碩士,現(xiàn)主要從事構(gòu)造解釋儲(chǔ)層預(yù)測(cè)研究工作,同時(shí)也研究疊前時(shí)間偏移處理。

    猜你喜歡
    射線孔徑傾角
    地球軸傾角的改斜歸正
    激光傾角儀在CT引導(dǎo)下經(jīng)皮肺穿刺活檢中的應(yīng)用
    車輪外傾角和前束角匹配研究
    北京汽車(2021年3期)2021-07-17 08:25:42
    “直線、射線、線段”檢測(cè)題
    『直線、射線、線段』檢測(cè)題
    不同滲透率巖芯孔徑分布與可動(dòng)流體研究
    赤石脂X-射線衍射指紋圖譜
    中成藥(2017年3期)2017-05-17 06:09:16
    分布式孔徑相參合成雷達(dá)技術(shù)
    基于子孔徑斜率離散采樣的波前重構(gòu)
    大孔徑淺臺(tái)階控制爆破在重慶地區(qū)的應(yīng)用
    重慶建筑(2014年12期)2014-07-24 14:00:32
    国产亚洲精品av在线| 色av中文字幕| 亚洲在线自拍视频| 免费无遮挡裸体视频| 老熟妇仑乱视频hdxx| 国产精品一及| 欧美xxxx黑人xx丫x性爽| 亚洲成a人片在线一区二区| 国内精品久久久久精免费| 国产黄片美女视频| 久久精品夜夜夜夜夜久久蜜豆| 九九热线精品视视频播放| 欧美成人一区二区免费高清观看| 成年女人毛片免费观看观看9| 波多野结衣巨乳人妻| 欧美高清成人免费视频www| 国产精品一区二区三区四区免费观看 | 天美传媒精品一区二区| xxxwww97欧美| 国产免费男女视频| 日韩精品青青久久久久久| 长腿黑丝高跟| 熟女人妻精品中文字幕| 色综合婷婷激情| 22中文网久久字幕| 欧美日韩中文字幕国产精品一区二区三区| 在线免费十八禁| 最近最新中文字幕大全电影3| 婷婷亚洲欧美| 国产蜜桃级精品一区二区三区| 欧美成人性av电影在线观看| 亚洲男人的天堂狠狠| 日韩,欧美,国产一区二区三区 | 欧美国产日韩亚洲一区| 在线a可以看的网站| 草草在线视频免费看| 91久久精品国产一区二区成人| 蜜桃亚洲精品一区二区三区| 少妇熟女aⅴ在线视频| 波多野结衣高清无吗| avwww免费| 在线免费观看的www视频| 国产激情偷乱视频一区二区| 国产伦精品一区二区三区四那| 日韩欧美在线乱码| 亚洲av二区三区四区| 精品人妻1区二区| 精品久久国产蜜桃| 色精品久久人妻99蜜桃| 能在线免费观看的黄片| 久久久久性生活片| 欧美3d第一页| 级片在线观看| 亚洲av中文字字幕乱码综合| 日韩,欧美,国产一区二区三区 | 中文字幕av成人在线电影| 少妇人妻一区二区三区视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲一级一片aⅴ在线观看| www.色视频.com| 国产精品野战在线观看| 乱码一卡2卡4卡精品| 窝窝影院91人妻| 中文亚洲av片在线观看爽| 午夜影院日韩av| 小说图片视频综合网站| 日日啪夜夜撸| 亚洲乱码一区二区免费版| 国产乱人视频| 欧美+亚洲+日韩+国产| 人人妻人人澡欧美一区二区| 国产主播在线观看一区二区| 深夜精品福利| 最近视频中文字幕2019在线8| 男女视频在线观看网站免费| 国产淫片久久久久久久久| 女人被狂操c到高潮| 国产色爽女视频免费观看| 国产免费男女视频| 色综合亚洲欧美另类图片| 午夜老司机福利剧场| 午夜久久久久精精品| 亚洲人成网站在线播放欧美日韩| 亚洲内射少妇av| 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 999久久久精品免费观看国产| 五月玫瑰六月丁香| 免费在线观看成人毛片| 国产成人一区二区在线| 亚洲精品日韩av片在线观看| 成人无遮挡网站| 中文字幕人妻熟人妻熟丝袜美| 欧美区成人在线视频| 我的女老师完整版在线观看| 乱系列少妇在线播放| 国产在线男女| 亚洲精品一区av在线观看| 免费看美女性在线毛片视频| 在线播放国产精品三级| 成人毛片a级毛片在线播放| 国产精品日韩av在线免费观看| 亚洲人与动物交配视频| 欧美精品啪啪一区二区三区| 亚洲av美国av| 国产成人av教育| 亚洲国产日韩欧美精品在线观看| 日本 av在线| 国产精品久久久久久av不卡| 99热精品在线国产| 婷婷精品国产亚洲av在线| 嫩草影院精品99| 91麻豆av在线| 欧美最新免费一区二区三区| 亚洲精品色激情综合| 成人无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 真实男女啪啪啪动态图| 久久精品夜夜夜夜夜久久蜜豆| 日本精品一区二区三区蜜桃| 九九久久精品国产亚洲av麻豆| 国产av不卡久久| videossex国产| aaaaa片日本免费| 国产精品不卡视频一区二区| 亚洲色图av天堂| 午夜福利欧美成人| 亚洲狠狠婷婷综合久久图片| 俺也久久电影网| 99热这里只有精品一区| 两性午夜刺激爽爽歪歪视频在线观看| 97碰自拍视频| 九色成人免费人妻av| 日韩高清综合在线| 欧美成人免费av一区二区三区| 少妇的逼水好多| 色哟哟哟哟哟哟| 国产精品一及| 最近最新免费中文字幕在线| 如何舔出高潮| netflix在线观看网站| 国内精品美女久久久久久| 亚洲自偷自拍三级| 一区二区三区免费毛片| 91麻豆av在线| 欧美激情国产日韩精品一区| 91久久精品电影网| 高清日韩中文字幕在线| 欧美成人一区二区免费高清观看| 最近在线观看免费完整版| 国产大屁股一区二区在线视频| 成人国产综合亚洲| av天堂在线播放| 精品午夜福利在线看| 欧美性猛交╳xxx乱大交人| 久久久色成人| 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 丰满人妻一区二区三区视频av| 免费观看人在逋| 最近在线观看免费完整版| 高清在线国产一区| 国产精品乱码一区二三区的特点| 欧美+亚洲+日韩+国产| 国产精品三级大全| 国产精品一区www在线观看 | 美女大奶头视频| 亚洲av一区综合| av女优亚洲男人天堂| 亚洲不卡免费看| 欧美绝顶高潮抽搐喷水| 日本熟妇午夜| 国产精品嫩草影院av在线观看 | 成人鲁丝片一二三区免费| 精品久久久久久久久久免费视频| 日韩精品中文字幕看吧| 成人国产一区最新在线观看| www.色视频.com| 日本a在线网址| 乱系列少妇在线播放| 免费在线观看日本一区| 三级男女做爰猛烈吃奶摸视频| 午夜福利欧美成人| 搡老妇女老女人老熟妇| 韩国av在线不卡| av专区在线播放| 国产免费一级a男人的天堂| 人妻丰满熟妇av一区二区三区| 久久欧美精品欧美久久欧美| 18禁裸乳无遮挡免费网站照片| 国产aⅴ精品一区二区三区波| 97碰自拍视频| 一边摸一边抽搐一进一小说| 国产午夜精品论理片| 久久精品国产清高在天天线| 欧美3d第一页| 大型黄色视频在线免费观看| 午夜免费激情av| 午夜日韩欧美国产| 精品久久久久久,| 久久久久久大精品| 不卡一级毛片| 亚洲av二区三区四区| 最好的美女福利视频网| 亚洲一级一片aⅴ在线观看| 91在线精品国自产拍蜜月| 性欧美人与动物交配| 中文字幕av在线有码专区| 嫩草影视91久久| 久久久久久国产a免费观看| 嫁个100分男人电影在线观看| 国产精品精品国产色婷婷| 色在线成人网| 麻豆国产97在线/欧美| 国产淫片久久久久久久久| 男女那种视频在线观看| 白带黄色成豆腐渣| 男女下面进入的视频免费午夜| 国产伦人伦偷精品视频| 国产成人av教育| 99久久精品一区二区三区| 国产午夜福利久久久久久| 免费av毛片视频| 日韩高清综合在线| 又黄又爽又刺激的免费视频.| 日韩中文字幕欧美一区二区| 99热网站在线观看| 桃红色精品国产亚洲av| 男女那种视频在线观看| 村上凉子中文字幕在线| а√天堂www在线а√下载| 中出人妻视频一区二区| 免费无遮挡裸体视频| 99热只有精品国产| 久9热在线精品视频| 麻豆一二三区av精品| 午夜免费成人在线视频| 免费人成视频x8x8入口观看| 又黄又爽又刺激的免费视频.| 此物有八面人人有两片| а√天堂www在线а√下载| 国产大屁股一区二区在线视频| 97碰自拍视频| 欧美成人a在线观看| 99热精品在线国产| 此物有八面人人有两片| 91在线精品国自产拍蜜月| 乱码一卡2卡4卡精品| 国产女主播在线喷水免费视频网站 | 国产亚洲精品久久久久久毛片| 国产精品无大码| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产 | 极品教师在线免费播放| 日本撒尿小便嘘嘘汇集6| 人人妻人人看人人澡| 久久久午夜欧美精品| 美女xxoo啪啪120秒动态图| 啦啦啦啦在线视频资源| 欧美绝顶高潮抽搐喷水| 日韩精品有码人妻一区| 美女 人体艺术 gogo| 久久精品91蜜桃| 91久久精品国产一区二区成人| 中文字幕高清在线视频| 变态另类丝袜制服| 久久久久精品国产欧美久久久| 人人妻人人澡欧美一区二区| 99热6这里只有精品| 欧美三级亚洲精品| 久久天躁狠狠躁夜夜2o2o| 中文资源天堂在线| 亚洲最大成人av| 国产色婷婷99| 色吧在线观看| 亚洲天堂国产精品一区在线| 国内揄拍国产精品人妻在线| 成熟少妇高潮喷水视频| 18禁黄网站禁片免费观看直播| 国产淫片久久久久久久久| 小说图片视频综合网站| 午夜福利在线观看吧| 老司机深夜福利视频在线观看| 国产精品三级大全| 精品一区二区三区视频在线| 日本 欧美在线| 久久天躁狠狠躁夜夜2o2o| 国产精品伦人一区二区| 中文字幕久久专区| 国产精品综合久久久久久久免费| 九色成人免费人妻av| 中文字幕精品亚洲无线码一区| 欧美一级a爱片免费观看看| 老司机午夜福利在线观看视频| 淫妇啪啪啪对白视频| 国产视频一区二区在线看| 久久久久久伊人网av| 69人妻影院| 久久久国产成人精品二区| 日韩国内少妇激情av| ponron亚洲| av在线天堂中文字幕| 国产v大片淫在线免费观看| 欧美日韩综合久久久久久 | 99九九线精品视频在线观看视频| 亚洲av电影不卡..在线观看| 国产精品98久久久久久宅男小说| 丰满人妻一区二区三区视频av| 狂野欧美白嫩少妇大欣赏| 欧美+日韩+精品| 亚洲自拍偷在线| 夜夜爽天天搞| 国产色爽女视频免费观看| 亚洲 国产 在线| 日韩欧美免费精品| 午夜久久久久精精品| 91午夜精品亚洲一区二区三区 | 少妇高潮的动态图| 国产av不卡久久| 草草在线视频免费看| 亚洲av第一区精品v没综合| 亚洲avbb在线观看| 精品一区二区免费观看| 夜夜看夜夜爽夜夜摸| 国产真实伦视频高清在线观看 | 久久久久久久精品吃奶| 久久久久精品国产欧美久久久| 成人国产一区最新在线观看| 99久久无色码亚洲精品果冻| 国内精品宾馆在线| 国产精华一区二区三区| 联通29元200g的流量卡| 久久久午夜欧美精品| 国产精品自产拍在线观看55亚洲| 人妻少妇偷人精品九色| 国内精品一区二区在线观看| 欧美日韩精品成人综合77777| 禁无遮挡网站| 小蜜桃在线观看免费完整版高清| 久久国产乱子免费精品| 亚洲欧美日韩卡通动漫| av黄色大香蕉| 男女啪啪激烈高潮av片| 国产高清有码在线观看视频| 久久久精品大字幕| 久久欧美精品欧美久久欧美| 亚洲七黄色美女视频| 国产在线精品亚洲第一网站| 国产 一区精品| 精品无人区乱码1区二区| 成人毛片a级毛片在线播放| 成年人黄色毛片网站| 成年女人看的毛片在线观看| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 国产淫片久久久久久久久| 欧美三级亚洲精品| 毛片女人毛片| 日本黄大片高清| 久久久色成人| 国产高清视频在线播放一区| 简卡轻食公司| 91午夜精品亚洲一区二区三区 | 999久久久精品免费观看国产| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在 | 亚洲精品亚洲一区二区| 1024手机看黄色片| 欧美+亚洲+日韩+国产| 精品人妻熟女av久视频| 国产av在哪里看| 嫁个100分男人电影在线观看| 亚洲国产日韩欧美精品在线观看| 欧美一区二区国产精品久久精品| 午夜视频国产福利| 精品一区二区免费观看| 麻豆精品久久久久久蜜桃| 人人妻,人人澡人人爽秒播| 乱系列少妇在线播放| 男女啪啪激烈高潮av片| 麻豆精品久久久久久蜜桃| 欧美日本亚洲视频在线播放| 麻豆精品久久久久久蜜桃| 午夜福利18| 国产亚洲欧美98| 午夜影院日韩av| 久久草成人影院| 色在线成人网| 午夜福利成人在线免费观看| 色综合亚洲欧美另类图片| 直男gayav资源| 日韩欧美一区二区三区在线观看| 最近最新中文字幕大全电影3| 亚洲精品色激情综合| 啦啦啦观看免费观看视频高清| 少妇的逼好多水| 真人做人爱边吃奶动态| 亚洲av一区综合| 成人高潮视频无遮挡免费网站| 欧美激情在线99| .国产精品久久| 在线观看午夜福利视频| 成人特级黄色片久久久久久久| 久久久色成人| 男女啪啪激烈高潮av片| 又爽又黄无遮挡网站| 日韩中字成人| 亚洲精品乱码久久久v下载方式| 免费av毛片视频| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 人人妻人人看人人澡| 亚洲精品国产成人久久av| av在线老鸭窝| a级毛片免费高清观看在线播放| a级一级毛片免费在线观看| 亚洲最大成人av| 内射极品少妇av片p| 一个人免费在线观看电影| 亚洲欧美日韩东京热| 成人午夜高清在线视频| 全区人妻精品视频| 国产麻豆成人av免费视频| 波多野结衣高清无吗| 性欧美人与动物交配| 少妇丰满av| 女的被弄到高潮叫床怎么办 | 国产黄色小视频在线观看| 国国产精品蜜臀av免费| 久久午夜亚洲精品久久| 在线天堂最新版资源| 久久天躁狠狠躁夜夜2o2o| 又黄又爽又免费观看的视频| 免费大片18禁| 亚洲精华国产精华液的使用体验 | 欧美不卡视频在线免费观看| a级毛片a级免费在线| 亚洲精品乱码久久久v下载方式| 午夜影院日韩av| 久久精品国产亚洲网站| 国产三级中文精品| 久久精品国产亚洲网站| 我的老师免费观看完整版| 99在线视频只有这里精品首页| 99视频精品全部免费 在线| 免费搜索国产男女视频| 亚洲av五月六月丁香网| 亚洲七黄色美女视频| 成人永久免费在线观看视频| 99九九线精品视频在线观看视频| 美女被艹到高潮喷水动态| 一a级毛片在线观看| 蜜桃久久精品国产亚洲av| 亚洲精品一区av在线观看| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久成人| 欧美精品国产亚洲| 最近最新中文字幕大全电影3| 亚洲精品成人久久久久久| 夜夜夜夜夜久久久久| 国产精品综合久久久久久久免费| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 日本爱情动作片www.在线观看 | 日韩,欧美,国产一区二区三区 | 国内精品久久久久久久电影| 1000部很黄的大片| 春色校园在线视频观看| 国产高清不卡午夜福利| 亚洲av成人av| 中文字幕熟女人妻在线| 淫妇啪啪啪对白视频| 久久婷婷人人爽人人干人人爱| 热99re8久久精品国产| 欧美日韩黄片免| 男插女下体视频免费在线播放| 老熟妇乱子伦视频在线观看| 日本a在线网址| 日日摸夜夜添夜夜添av毛片 | 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av在线| 色播亚洲综合网| 日韩 亚洲 欧美在线| 我的老师免费观看完整版| 美女cb高潮喷水在线观看| 网址你懂的国产日韩在线| 两个人视频免费观看高清| 国产精品人妻久久久久久| 波多野结衣高清作品| 中文在线观看免费www的网站| 国产三级在线视频| 在线免费十八禁| 1000部很黄的大片| 高清毛片免费观看视频网站| 麻豆成人av在线观看| 少妇猛男粗大的猛烈进出视频 | 中文欧美无线码| 性色av一级| 亚洲精品久久久久久婷婷小说| 午夜激情福利司机影院| 亚洲av欧美aⅴ国产| 国产 精品1| 中文天堂在线官网| 国产美女午夜福利| 久久人妻熟女aⅴ| 色婷婷久久久亚洲欧美| 亚洲欧美一区二区三区黑人 | 国产精品无大码| 男女边吃奶边做爰视频| 国产大屁股一区二区在线视频| 99久久精品热视频| 亚洲美女搞黄在线观看| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 日本午夜av视频| 如何舔出高潮| 久久精品国产自在天天线| 国产精品蜜桃在线观看| 日本色播在线视频| 成人二区视频| 日本色播在线视频| 国产精品伦人一区二区| 欧美日韩精品成人综合77777| 嘟嘟电影网在线观看| 少妇丰满av| 熟妇人妻不卡中文字幕| 亚洲国产av新网站| av线在线观看网站| 九九久久精品国产亚洲av麻豆| 日日啪夜夜爽| 久久久久久久久大av| 直男gayav资源| 18禁裸乳无遮挡动漫免费视频| av在线老鸭窝| 一边亲一边摸免费视频| 久久久精品免费免费高清| 国产成人免费观看mmmm| 日韩欧美一区视频在线观看 | 狂野欧美白嫩少妇大欣赏| 久久精品国产自在天天线| 国产色婷婷99| 91久久精品国产一区二区成人| 午夜福利网站1000一区二区三区| 国产欧美日韩一区二区三区在线 | 七月丁香在线播放| 91精品一卡2卡3卡4卡| 久久午夜福利片| 国产亚洲5aaaaa淫片| 久久精品国产a三级三级三级| 国产精品爽爽va在线观看网站| 亚洲av二区三区四区| 99热6这里只有精品| 菩萨蛮人人尽说江南好唐韦庄| 成人午夜精彩视频在线观看| 在线观看免费高清a一片| 久久久久久九九精品二区国产| 欧美成人午夜免费资源| 国产精品免费大片| 久久女婷五月综合色啪小说| 午夜激情福利司机影院| 成人国产av品久久久| 亚洲成人一二三区av| 97超视频在线观看视频| 国产在线视频一区二区| 春色校园在线视频观看| 观看av在线不卡| 久热这里只有精品99| 毛片一级片免费看久久久久| 午夜日本视频在线| 久久久a久久爽久久v久久| 一级毛片久久久久久久久女| 寂寞人妻少妇视频99o| 欧美极品一区二区三区四区| 免费久久久久久久精品成人欧美视频 | 国产午夜精品久久久久久一区二区三区| 精品亚洲成国产av| 久久久久久久精品精品| 欧美日韩综合久久久久久| 在线免费十八禁| 91午夜精品亚洲一区二区三区| 五月伊人婷婷丁香| 成人二区视频| 免费播放大片免费观看视频在线观看| 国产高清三级在线| 一级毛片 在线播放| 精品亚洲成a人片在线观看 | 成人二区视频| 99久久中文字幕三级久久日本| 男女无遮挡免费网站观看| 国产精品一区二区在线不卡| 欧美老熟妇乱子伦牲交| 久久 成人 亚洲| 熟女电影av网| 国产人妻一区二区三区在| 日韩中字成人| 国产 一区 欧美 日韩| 黑人高潮一二区| 小蜜桃在线观看免费完整版高清| 国产高清有码在线观看视频| 黑丝袜美女国产一区| 熟女人妻精品中文字幕| 街头女战士在线观看网站| 精品99又大又爽又粗少妇毛片| 美女国产视频在线观看| 成人午夜精彩视频在线观看| 少妇裸体淫交视频免费看高清| 国产精品一区二区在线观看99| 亚洲精品日本国产第一区| av在线观看视频网站免费| 亚洲国产成人一精品久久久| 国产久久久一区二区三区| 日韩精品有码人妻一区| 精品亚洲成a人片在线观看 | 一级爰片在线观看| 亚洲精品国产成人久久av|