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

    玉北地區(qū)碳酸鹽巖縫洞型儲(chǔ)層精細(xì)成像方法應(yīng)用研究

    2015-06-27 05:54:49劉志成劉小民楊瑞娟
    石油物探 2015年4期
    關(guān)鍵詞:縫洞波場(chǎng)碳酸鹽巖

    李 鵬,劉志成,楊 楠,劉小民,楊瑞娟

    (中國(guó)石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

    玉北地區(qū)碳酸鹽巖縫洞型儲(chǔ)層精細(xì)成像方法應(yīng)用研究

    李 鵬,劉志成,楊 楠,劉小民,楊瑞娟

    (中國(guó)石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

    針對(duì)新疆塔里木盆地玉北地區(qū)目標(biāo)儲(chǔ)層埋藏深、非均質(zhì)性強(qiáng)以及地震資料信噪比低、波場(chǎng)復(fù)雜的特點(diǎn),以逆時(shí)疊前深度偏移成像為核心、深度域速度建模為關(guān)鍵、保真的弱信號(hào)補(bǔ)償為基礎(chǔ),建立了一套面向玉北碳酸鹽巖縫洞型儲(chǔ)層的高精度地震成像技術(shù)流程。將該技術(shù)流程應(yīng)用于玉北6井區(qū)474km2實(shí)際三維地震資料處理,所得結(jié)果與傳統(tǒng)疊前時(shí)間偏移結(jié)果對(duì)比分析表明,該技術(shù)流程對(duì)于復(fù)雜構(gòu)造帶成像、“串珠”成像、縫洞儲(chǔ)層刻畫等具有明顯的優(yōu)勢(shì),有效提高了研究區(qū)奧陶系碳酸鹽巖縫洞型目標(biāo)儲(chǔ)層的成像精度。

    逆時(shí)偏移(RTM);速度建模;弱信號(hào)補(bǔ)償;縫洞儲(chǔ)層

    玉北地區(qū)古地貌特征復(fù)雜,由于長(zhǎng)期的風(fēng)化淋濾及巖溶作用,風(fēng)化面附近形成了良好的儲(chǔ)層,為油氣聚集提供了優(yōu)越的儲(chǔ)集空間,成為塔里木盆地具有勘探前景的油氣聚集區(qū)帶。但玉北地區(qū)奧陶系巖溶縫洞型目標(biāo)儲(chǔ)層埋藏深、類型復(fù)雜、儲(chǔ)層非均性強(qiáng)的地質(zhì)特點(diǎn)造成資料信噪比低、地震波場(chǎng)復(fù)雜等一系列地球物理難題,對(duì)后續(xù)的速度建模和地震成像帶來更高要求?;谌暡ǚ匠痰哪鏁r(shí)偏移(RTM)技術(shù)克服了偏移傾角和偏移孔徑的限制,具有相位準(zhǔn)確、成像精度高、保幅性能好、對(duì)縱橫向劇烈速度變化和高陡傾角構(gòu)造適應(yīng)性強(qiáng)、甚至可以利用回轉(zhuǎn)波、多次波等特殊類型波場(chǎng)正確成像等優(yōu)點(diǎn),是目前高度復(fù)雜介質(zhì)高精度地震成像的理想技術(shù)手段。

    Whitemore[1]早在20世紀(jì)80年代就提出了RTM的概念,但受限于海量計(jì)算和存儲(chǔ)以及特有低頻噪聲等技術(shù)瓶頸,RTM技術(shù)一直沒有能夠?qū)崿F(xiàn)實(shí)用化。近年來隨著計(jì)算機(jī)軟硬件技術(shù)的飛速發(fā)展,RTM技術(shù)在實(shí)際生產(chǎn)中得到了越來越廣泛的應(yīng)用。借鑒求解波動(dòng)方程的三類主流算法:有限差分法、偽譜法和有限元法,Wu等[2]、劉紅偉等[3]、劉定進(jìn)等[4]發(fā)展了高階有限差分RTM技術(shù),孫文博等[5]、薛東川等[6]和李信富等[7]分別發(fā)展了偽譜法RTM技術(shù)和有限元RTM技術(shù),而董淵等[8]和王童奎等[9]根據(jù)這三類方法的優(yōu)缺點(diǎn)又分別實(shí)現(xiàn)了RTM的有限元—有限差分混合解法和有限元—偽譜混合解法;為了解決RTM海量計(jì)算、海量存儲(chǔ)、低頻噪聲等技術(shù)瓶頸問題,李博等[10]、趙磊等[11]在GPU平臺(tái)上實(shí)現(xiàn)了RTM,劉紅偉等[12]深入分析了RTM產(chǎn)生的特有低頻噪聲和存儲(chǔ)需求并提出了解決方案;隨著RTM的實(shí)用化程度的提高,一些學(xué)者又將該技術(shù)拓展到新的領(lǐng)域,張美根等[13]實(shí)現(xiàn)了各向異性RTM,底青云[14]、陳可洋[15]、張會(huì)星[16]等人實(shí)現(xiàn)了彈性波RTM,Zhang等[17-18]實(shí)現(xiàn)了真振幅RTM,劉紅偉等[19]實(shí)現(xiàn)了起伏地表RTM。

    本文根據(jù)玉北地區(qū)地震資料的特點(diǎn)和處理要求,從影響地震成像精度的3個(gè)要素:①正確反映地下地質(zhì)情況的品質(zhì)較好的地震數(shù)據(jù);②與地下地質(zhì)情況匹配較好的速度模型;③高精度偏移成像方法出發(fā),建立了一套以RTM疊前深度偏移成像為核心、深度域速度建模為關(guān)鍵、保真的弱信號(hào)補(bǔ)償為基礎(chǔ)的面向玉北地區(qū)碳酸鹽巖縫洞型儲(chǔ)層的高精度地震成像技術(shù)流程。

    1 地震資料的弱信號(hào)補(bǔ)償

    玉北地區(qū)地震資料一方面受地表沙丘的影響;

    另一方面由于目標(biāo)儲(chǔ)層埋藏深、波阻抗差異小等原因,奧陶系內(nèi)部信噪比低,有效波信號(hào)湮沒在噪聲中,常規(guī)的疊前去噪技術(shù)很難在去噪過程中對(duì)弱信號(hào)保護(hù)進(jìn)行有效監(jiān)控,容易損失有效信號(hào)而影響到最終的成像效果。弱信號(hào)補(bǔ)償技術(shù)在原始單炮上對(duì)有效弱信號(hào)進(jìn)行補(bǔ)償,在此基礎(chǔ)上再對(duì)不同類型的噪聲進(jìn)行壓制,在壓制過程中能夠較好地實(shí)施質(zhì)量控制,最大限度地保留有效信號(hào),進(jìn)而提高成像品質(zhì)。其基本思路是利用原始數(shù)據(jù)CMP道集構(gòu)建宏面元,通過Hilbert變換(公式(1))得到正弦、余弦相位(公式(2))道集和瞬時(shí)能量(公式(3))道集,經(jīng)子波自識(shí)別和能量、相位以及頻率恢復(fù)后實(shí)現(xiàn)保真重構(gòu)(公式(4)),輸出信號(hào)補(bǔ)償?shù)兰糜诤罄m(xù)處理。其優(yōu)點(diǎn)是將振幅和相位進(jìn)行分離,對(duì)相位處理的同時(shí),能保護(hù)振幅的原始特性。

    (1)

    (2)

    (3)

    (4)

    圖1為弱信號(hào)補(bǔ)償前、后的單炮記錄對(duì)比,可以看到,經(jīng)過弱信號(hào)補(bǔ)償后單炮層間尤其是深層有效信息得到較明顯的補(bǔ)償,而在補(bǔ)償前后單炮記錄差值中也幾乎看不到有效信息,進(jìn)一步印證了本方法的振幅保真性。圖2 為弱信號(hào)補(bǔ)償前、后的速度譜和超道集對(duì)比,可以看到,補(bǔ)償后速度譜的品質(zhì)更高,有效信號(hào)能量團(tuán)聚集得非常好,為后續(xù)速度分析與提高建模精度提供了保障。圖3是弱信號(hào)補(bǔ)償前、后疊加剖面對(duì)比,補(bǔ)償后剖面信噪比明顯提高,尤其碳酸鹽巖目標(biāo)儲(chǔ)層(圖中綠線框)的波場(chǎng)豐富,地層之間的接觸關(guān)系更加清晰。而從信噪比模型對(duì)比圖上(圖4)可以看到,補(bǔ)償后信噪比剖面較為均衡,信噪比明顯提高,這有利于后續(xù)的地震偏移成像處理。

    圖1 弱信號(hào)補(bǔ)償前(a)、后(b)單炮記錄及單炮記錄補(bǔ)償前后的差值(c)

    圖2 弱信號(hào)補(bǔ)償前(a)、后(b)的速度譜及超道集

    圖3 弱信號(hào)補(bǔ)償前(a)、后(b)的疊加剖面

    圖4 弱信號(hào)補(bǔ)償前(a)、后(b)信噪比模型對(duì)比

    2 遞進(jìn)式深度域速度分析與建模

    對(duì)于任何偏移方法,速度模型的正確性都是決定構(gòu)造成像質(zhì)量的關(guān)鍵,尤其是RTM算法對(duì)速度模型的依賴性更強(qiáng)、反應(yīng)更敏感、要求更高。而玉北地區(qū)碳酸鹽巖裂縫、溶蝕孔洞型目標(biāo)儲(chǔ)層的儲(chǔ)集空間類型多樣,空間配置復(fù)雜,而且規(guī)模、尺度不同,因而地震響應(yīng)特征復(fù)雜,增加了速度分析與建模的難度。面向RTM的速度建模過程是一個(gè)處理解釋一體化的過程,主要包括初始模型建立和模型優(yōu)化迭代兩個(gè)步驟。圖5為面向RTM的遞進(jìn)式深度域速度分析與建模流程。

    2.1 初始模型的建立

    初始模型的建立非常關(guān)鍵,建立一個(gè)比較準(zhǔn)確的初始模型直接決定了后期速度模型修正的精度和迭代效率。由于疊前時(shí)間偏移對(duì)速度的敏感度較小,我們首先進(jìn)行基于疊前時(shí)間偏移的速度分析,并采用剩余速度自動(dòng)拾取的方法修正偏移速度場(chǎng),經(jīng)多次迭代后得到最優(yōu)的時(shí)間域瞬時(shí)層速度模型,為疊前深度偏移的初始模型的建立奠定良好的基礎(chǔ);然后以地質(zhì)模型為約束條件,沿解釋的層位從瞬時(shí)層速度體上提取各層的層速度,再經(jīng)過平滑處理形成初始層速度;利用初始層速度,將時(shí)間域?qū)嶓w模型轉(zhuǎn)化為深度域?qū)嶓w模型;將層速度充填到深度域?qū)嶓w模型中就得到了一個(gè)比較理想的初始深度域速度模型。圖6分別為獲得的深度域初始速度模型和對(duì)應(yīng)的疊前深度偏移剖面與共成像點(diǎn)道集,可以看到剖面大套地層結(jié)構(gòu)清晰,地震資料信噪比高,古生界以上地震反射資料品質(zhì)好,說明建立的初始速度模型整體精度較高;但同時(shí)可以看到在中深層(尤其是奧陶系碳酸鹽巖目標(biāo)儲(chǔ)層)成像品質(zhì)不高、道集同相軸不平,說明該層位速度模型精度還有進(jìn)一步提高的必要。

    圖5 面向RTM的遞進(jìn)式深度域速度分析與建模流程

    圖6 深度域初始速度模型(a)、疊前深度偏移剖面(b)和共成像點(diǎn)道集(c)

    2.2 速度模型的迭代與修正

    與疊前時(shí)間偏移相比,疊前深度偏移對(duì)速度的敏感性更高。為了獲得精度更高的可靠速度模型,我們借助疊前深度偏移方法,采用遞進(jìn)式思路對(duì)速度模型進(jìn)行進(jìn)一步細(xì)化。首先利用垂向速度分析技術(shù)進(jìn)行第一輪迭代,經(jīng)過多次迭代可獲得垂向變化的精度較高的連續(xù)介質(zhì)模型;然后利用基于層位的層析成像技術(shù)進(jìn)行第二輪迭代,這種技術(shù)引入了更多的地質(zhì)解釋方案,并考慮了聯(lián)絡(luò)線方向的閉合,由此獲得的橫向和垂向上均有變化的速度模型更加符合地質(zhì)發(fā)育史的連續(xù)沉積和地層壓實(shí)作用,能夠保證成像結(jié)果在主要地質(zhì)界面(即大套地層模型中的解釋層)成像清楚;由于基于層位的層析成像主要考慮大套層位的平均層速度,對(duì)層速度具有平均效應(yīng),對(duì)層內(nèi)大部分同相軸是合適的,而對(duì)那些層間速度異常或?qū)娱g速度變化較大的區(qū)域卻可能還有一些局部速度誤差,因此我們?cè)倮没诰W(wǎng)格層析成像技術(shù)進(jìn)行第三輪迭代,對(duì)速度模型進(jìn)行精細(xì)修正,最終得到與地下地質(zhì)構(gòu)造匹配度較高的可靠速度模型,為RTM應(yīng)用提供支撐。圖7和圖8分別為由初始速度模型和三輪迭代修正的最終速度模型獲得的疊前深度偏移剖面與共成像點(diǎn)道集,可以看到經(jīng)三輪迭代修正后,剖面(圖8)的質(zhì)量得到進(jìn)一步提高,尤其是奧陶系碳酸鹽巖目標(biāo)儲(chǔ)層的成像品質(zhì)得到較大幅度改進(jìn),一些弱小縫洞異常體刻畫得更清楚,成像道集也校得更平,說明速度模型的精度更高。圖9 所示的最終速度模型與疊前深度偏移剖面疊合圖也表明速度模型與地下的地質(zhì)情況吻合非常好。

    圖7 初始速度模型疊前深度偏移剖面(a)與共成像點(diǎn)道集(b)

    圖8 經(jīng)3輪迭代后的最終速度模型疊前深度偏移剖面(a)與共成像點(diǎn)道集(b)

    圖9 最終速度模型與疊前深度偏移剖面疊合顯示

    3 保幅的逆時(shí)偏移成像

    玉北地區(qū)地震資料的整體品質(zhì)較高,傳統(tǒng)的疊前時(shí)間偏移剖面上大套地層結(jié)構(gòu)清晰,信噪比高,古生界以上地震反射資料品質(zhì)好。但奧陶系內(nèi)幕目標(biāo)儲(chǔ)層在信噪比、縫洞成像精度、地層之間接觸關(guān)系等方面無法滿足儲(chǔ)層精細(xì)描述的需求。究其原因,一方面在于碳酸鹽巖儲(chǔ)層速度橫向變化劇烈、地震波場(chǎng)復(fù)雜,基于地下橫向均勻理論的疊前時(shí)間偏移不再適用,需要理論更先進(jìn)的成像方法提供支撐;另一方面小規(guī)模、小尺度縫洞儲(chǔ)層對(duì)振幅保真性提出了更高的要求?;陔p程波動(dòng)理論的逆時(shí)偏移技術(shù)遵循了波的繞射、反射和折射定律,能夠較好地處理復(fù)雜介質(zhì)中的復(fù)雜波場(chǎng)問題,能自然地處理多路徑問題以及由速度變化引起的聚焦或焦散效應(yīng),并具有很好的振幅保持特性,是玉北地區(qū)碳酸鹽巖縫洞儲(chǔ)層高精度地震成像的理想手段。

    3.1 緊致差分波場(chǎng)延拓算子

    選取一種計(jì)算精度好、效率高的正演算法是RTM技術(shù)的核心?,F(xiàn)有的地震波正演手段大致可分為3類:有限差分法、有限元法和偽譜法。相對(duì)于另兩種方法,有限差分法憑借其算法簡(jiǎn)單快速、對(duì)速度場(chǎng)任意變化適應(yīng)性強(qiáng)的理論優(yōu)勢(shì)仍然占據(jù)工業(yè)應(yīng)用的主要位置。

    從全聲波方程出發(fā):

    (5)

    式中:c(x,y,z)表示地下介質(zhì)速度;p(x,y,z,t)為地表接收波場(chǎng),為了使波場(chǎng)延拓精度更高,平衡模擬精度與計(jì)算量,本文引入隱式的緊致差分格式對(duì)方程(5)進(jìn)行求解,形成基于優(yōu)化的五對(duì)角緊致差分RTM波場(chǎng)延拓算子,(6)式為該偏移波場(chǎng)延拓算子方程:

    βp″i-2+αp″i-1+p″i+αp″i+1+βp″i+2=

    (6)

    式中:p為計(jì)算波場(chǎng);i為離散化為空間坐標(biāo);h是空間采樣間隔;α,β,a,b和c都是緊致差分系數(shù)。

    圖10為傳統(tǒng)高階有限差分和本文緊致差分算子在均勻介質(zhì)中的脈沖響應(yīng),空間網(wǎng)格參數(shù)Δx和Δy都為20m,深度步長(zhǎng)Δz為10m,均勻介質(zhì)的速度為c=2000m/s,所用雷克子波主頻為fpeak=30Hz。圖10a為高階有限差分RTM算子(時(shí)間2階、空間8階)的脈沖響應(yīng),可以看到頻散比較嚴(yán)重,如果在此基礎(chǔ)上提高精度就得減小空間網(wǎng)格參數(shù),這樣計(jì)算量和存儲(chǔ)需求就會(huì)急劇上升。而本文緊致差分算子的脈沖響應(yīng)頻散壓制比較干凈,說明了緊致差分格式能夠在相同的基架上構(gòu)造出更高精度的差分格式,可以保證波場(chǎng)在傳播過程中頻散更小和遞推過程更加穩(wěn)定(圖10b)。

    3.2 動(dòng)力學(xué)成像條件

    成像條件是波動(dòng)類地震偏移成像方法的另外一個(gè)核心要素。本文采用帶阻尼因子的動(dòng)力學(xué)成像條件,其理論公式可以表示為:

    (7)

    式中:tmax是最大記錄時(shí)間;Ss(x,y,z,t)為正向外推的震源波場(chǎng);Rs(x,y,z,t)為反向外推的記錄波場(chǎng);I(x,y,z)為點(diǎn)(x,y,z)的成像結(jié)果。

    RTM波場(chǎng)延拓算子具有非常好的振幅保真性能,動(dòng)力學(xué)成像條件也具有比較高的振幅保真性,兩者結(jié)合形成了構(gòu)造成像精度高、振幅保真性好的RTM技術(shù)。這種振幅保真性高的RTM技術(shù)系列不但有助于碳酸鹽巖縫洞型特殊儲(chǔ)層的精確成像,而且可以為后續(xù)的AVO等屬性分析提供更真實(shí)的地震信息,促進(jìn)更深層次的勘探開發(fā)。

    圖10 RTM算子在均勻介質(zhì)中的脈沖響應(yīng)

    4 實(shí)際資料應(yīng)用效果分析

    本研究對(duì)塔里木盆地玉北地區(qū)玉北6井區(qū)474km2實(shí)際三維地震資料進(jìn)行了常規(guī)疊前時(shí)間偏移處理(PSTM)和逆時(shí)疊前深度偏移處理(RTM),野外采集面元25m×25m,總炮數(shù)27929炮,采用的計(jì)算孔徑縱向13600m和橫向7000m,頻帶范圍4~70Hz,時(shí)間延拓步長(zhǎng)為0.5ms,最大延拓時(shí)間為8.0s。相干屬性分析技術(shù)和振幅變化率屬性分析技術(shù)被廣泛應(yīng)用于塔里木盆地碳酸鹽巖縫洞型儲(chǔ)層預(yù)測(cè)[20-24],本研究利用地震剖面對(duì)比、地震相干屬性分析和振幅變化率屬性分析等技術(shù)進(jìn)行了PSTM和RTM偏移成像效果的對(duì)比。

    PSTM和RTM偏移成像結(jié)果對(duì)比分析發(fā)現(xiàn),逆時(shí)疊前深度偏移成像在高陡構(gòu)造帶的斷面成像更加清晰、地震剖面中“串珠”成像更精細(xì)、“串珠”的形態(tài)更規(guī)整、“串珠”能量與圍巖的差異更加明顯(圖11),平面相干屬性分析圖中(圖12),逆時(shí)偏移對(duì)大斷裂的成像更加清楚,局部小斷裂成像更為清晰;通過振幅變化率屬性分析(圖13),在工區(qū)南部,逆時(shí)偏移結(jié)果中“串珠”的個(gè)數(shù)明顯增多(圖13b)。逆時(shí)偏移技術(shù)明顯提高了玉北地區(qū)奧陶系碳酸鹽巖縫洞儲(chǔ)層成像的精度。

    圖11 PSTM成像結(jié)果(a)與RTM成像結(jié)果(b)對(duì)比

    圖12 鷹山組頂面沿層相干屬性分析剖面

    圖13 鷹山組沿層振幅變化率屬性分析

    5 結(jié)束語

    針對(duì)玉北地區(qū)碳酸鹽巖目標(biāo)儲(chǔ)層埋藏深、類型復(fù)雜、儲(chǔ)層非均性強(qiáng)的地質(zhì)特點(diǎn)造成的資料信噪比低、地震波場(chǎng)復(fù)雜等一系列地球物理難題,本文采用弱信號(hào)補(bǔ)償技術(shù)對(duì)碳酸鹽巖目標(biāo)儲(chǔ)層中弱有效信號(hào)的能量、相位以及頻率實(shí)現(xiàn)保真重構(gòu),有效地解決了埋藏深、波阻抗差異小導(dǎo)致的低信噪比問題,為后續(xù)的速度建模與地震成像奠定了良好的數(shù)據(jù)基礎(chǔ);采用遞進(jìn)式深度域速度建模技術(shù)為后續(xù)的RTM成像建立了與地下地質(zhì)結(jié)構(gòu)匹配較好的可靠速度模型;采用基于緊致差分的保幅RTM地震成像技術(shù)有效解決了碳酸鹽巖縫洞儲(chǔ)層非均性強(qiáng)、波場(chǎng)復(fù)雜問題,從而建立了一套以RTM逆時(shí)偏移技術(shù)為核心的高精度地震成像技術(shù)流程,并將該技術(shù)流程應(yīng)用于玉北6井區(qū)474km2實(shí)際三維地震資料處理,取得了較好的應(yīng)用效果。不僅提高了研究區(qū)小“串珠”、小斷裂和奧陶系內(nèi)幕地層成像的精度,同時(shí),“串珠”的空間位置更加準(zhǔn)確,為該區(qū)后續(xù)儲(chǔ)層預(yù)測(cè)、井位部署等提供了高品質(zhì)的地震資料。

    [1] Whitemore N D.Iterative depth imaging by backward time propagation[J].Expanded Abstracts of 53rdAnnual Internat SEG Mtg,1983,382-384

    [2] Wu W,Lines L R,Lu H.Analysis of higher-order finite-difference schemes in 3-D reverse-time migration[J].Geophysics,1996,61(3):845-846

    [3] 劉紅偉,李博,劉洪,等.地震疊前逆時(shí)偏移高階有限差分算法及GPU實(shí)現(xiàn)[J].地球物理學(xué)報(bào),2010,53(7):1725-1733 Liu H W,Li B,Liu H,et al.The algorithm of high order finite difference pre-stack reverse time migration and GPU implementation[J].Chinese Journal of Geophysics,2010,53(7):1725-1733

    [4] 劉定進(jìn),楊勤勇,方伍寶,等.疊前逆時(shí)深度偏移成像的實(shí)現(xiàn)與應(yīng)用[J].石油物探,2011,50(6):545-549 Liu D J,Yang Q Y,Fang W B,et al.Realization and practices of prestack reverse-time depth migration[J].Geophysical Prospecting for Petroleum,2011,50(6):545-549

    [5] 孫文博,孫贊東.基于偽譜法的VSP逆時(shí)偏移及其應(yīng)用研究[J].地球物理學(xué)報(bào),2010,53(9):2196-2203 Sun W B,Sun Z D.VSP reverse time migration based on the pseudo-spectral method and its applications[J].Chinese Journal of Geophysics,2010,53(9):2196-2203

    [6] 薛東川,王尚旭.波動(dòng)方程有限元疊前逆時(shí)偏移[J].石油地球物理勘探,2008,43(1):18-21 Xue D C,Wang S X.Wave-equation finite-element prestack reverse-time migration[J].Oil Geophysical Prospecting,2008,43(1):18-21

    [7] 李信富,張美根.顯式分形插值在有限元疊前逆時(shí)偏移成像中的應(yīng)用[J].地球物理學(xué)進(jìn)展,2008,23(5): 1406-1411 Li X F,Zhang M G.Application of explicit fractal interpolation in pre-stack finite element reverse-time migration[J].Progress in Geophysics,2008,23(5):1406-1411

    [8] 董淵,楊慧珠,杜啟振.有限元-有限差分法二維波動(dòng)逆時(shí)偏移初探[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,27(6):25-29 Dong Y,Yang H Z,Du Q Z.Study on reverse-time migration of two-dimensional wave equation with finite element finite difference method[J].Journal of China University of Petroleum(Edition of Natural Science),2003,27(6):25-29

    [9] 王童奎,付興深,朱德獻(xiàn),等.譜元法疊前逆時(shí)偏移研究[J].地球物理學(xué)進(jìn)展,2008,23(3):681-685 Wang T K,Fu X S,Zhu D X,et al.Spectral-element method for prestack reverse-time migration[J].Progress in Geophysics,2008,23(3):681-685

    [10] 李博,劉紅偉,劉國(guó)峰,等.地震疊前逆時(shí)偏移算法的CPU/GPU實(shí)施對(duì)策[J].地球物理學(xué)報(bào),2010,53(12):2938-2943 Li B,Liu H W,Liu G F,et al.Computational strategy of seismic pre-stack reverse time migration on CPU/GPU[J].Chinese Journal of Geophysics,2010,53(12):2938-2943

    [11] 趙磊,王華忠,劉守偉.逆時(shí)深度偏移成像方法及其在CPU/GPU異構(gòu)平臺(tái)上的實(shí)現(xiàn)[J].巖性油氣藏,2010,22(7):36-41 Zhao L,Wang H Z,Liu S W.Reverse-time depth migration imaging and its application at CPU/GPU platform[J].Lithologic Reservoirs,2010,22(7):36-41

    [12] 劉紅偉,劉洪,鄒振,等.地震疊前逆時(shí)偏移中的去噪與存儲(chǔ)[J].地球物理學(xué)報(bào),2010,53(9):2171-2180 Liu H W,Liu H,Zhou Z,et al.The problems of denoise and storage in seismic reverse time migration[J].Chinese Journal of Geophysics,2010,53(9):2171-2180

    [13] 張美根,王妙月.各向異性彈性波有限元疊前逆時(shí)偏移[J].地球物理學(xué)報(bào),2001,44(5):711-719 Zhang M G,Wang M Y.Prestack finite element reverse-time migration for anisotropic elastic waves[J].Chinese Journal of Geophysics,2001,44(5):711-719

    [14] 底青云,王妙月.彈性波有限元逆時(shí)偏移技術(shù)研究[J].地球物理學(xué)報(bào),1997,40(4):570-579 Di Q Y,Wang M Y.The study of finite element reverse-time migration for elastic wave[J].Chinese Journal of Geophysics,1997,40(4):570-579

    [15] 陳可洋.各向異性彈性波動(dòng)方程多分量聯(lián)合疊后逆時(shí)偏移[J].內(nèi)陸地震,2009,23(4):455-460 Chen K Y.Multi-component combined post-stack reverse-time migration of anisotropic elastic wave equation[J].Inland Earthquake,2009,23(4):455-460

    [16] 張會(huì)星,寧書年.彈性波動(dòng)方程疊前逆時(shí)偏移[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào),2002,31(5):371-375 Zhang H X,Ning S N.Pre-stack reverse time migration of elastic wave equation[J].Journal of China University of Mining & Technology,2002,31(5):371-375

    [17] Zhang Y,Sun J,Gray S.Reverse-time migration:amplitude and implementation issues[J].Expanded Abstracts of 77thAnnual Internat SEG Mtg,2007,2145-2149[18] Zhang Y,James S.Practical issues of reverse time migration:true amplitude gathers,noise removal and harnlonic source encoding[J].First Break,2009,26(1):29-35

    [19] 劉紅偉,劉洪,李博,等.起伏地表疊前逆時(shí)偏移理論及GPU加速技術(shù)[J].地球物理學(xué)報(bào),2011,54(7):1883-1892 Liu H W,Liu H,Li B,et al.Pre-stack reverse time migration for rugged topography and GPU acceleration technology[J].Chinese Journal of Geophysics,2011,54(7):1883-1892

    [20] 楊江峰,劉定進(jìn),朱文斌,等,逆時(shí)偏移技術(shù)在碳酸鹽巖縫-洞儲(chǔ)層成像中的應(yīng)用[J].高校地質(zhì)學(xué)報(bào),2014,20(2):303-308 Yang J F,Liu D J,Zhu W B,et al.Technique and application of reverse time migration of fractured-vuggy carbonate reservoir[J].Geological Journal of China Universities,2014,20(2):303-308

    [21] 楊江峰,關(guān)達(dá),朱文斌,等.塔里木盆地塔中卡1三維工區(qū)碳酸鹽巖儲(chǔ)層預(yù)測(cè)技術(shù)與應(yīng)用[J].石油與天然氣地質(zhì),2012,33(1):61-69 Yang J F,Guan D,Zhu W B,et al.Carbonate reservoir prediction technology and its application in Ka-1 3D acreage in Tazhong area Tarim Basin[J].Oil & Gas Geology,2012,33(1):61-69

    [22] 王振卿,王宏斌,龔洪林,等.塔中地區(qū)碳酸鹽巖儲(chǔ)層預(yù)測(cè)技術(shù)[J].天然氣技術(shù),2009,3(1):23-26 Wang Z Q,Wang H B,Gong H L,et al.Reservoir prediction technologies of carbonate rock in central Tarim Basin[J].Natural Gas Technology,2009,3(1):23-26

    [23] 閆相賓,管路平,王世星.塔里木盆地碳酸鹽巖縫洞系統(tǒng)的地震響應(yīng)特征及預(yù)測(cè)[J].石油與天然氣地質(zhì),2007,28(6):828-835 Yan X B,Guan L P,Wang S X.Seismic response characteristics and prediction of carbonate rock fractured-vuggy systems in the Tarim Basin[J].Oil & Gas Geology,2007,28(6):828-835

    [24] 楊子川.塔河油田碳酸鹽巖儲(chǔ)層預(yù)測(cè)技術(shù)與應(yīng)用[J].勘探地球物理進(jìn)展,2004,27(6):432-438 Yang Z C.Carbonate reservoir prediction technology and its application in Tahe oilfield[J].Progress in Exploration Geophysics,2004,27(6):432-438

    (編輯:朱文杰)

    Application of fine target reservoir imaging in fractured-vuggy carbonate reservoir Yubei area

    Li Peng,Liu Zhicheng,Yang Nan,Liu Xiaomin,Yang Ruijuan

    (SinopecGeophysicalResearchInstitute,Nanjing211103,China)

    The target reservoir of Yubei area in the Tarim Basin has complex characteristics,including deep buried,strong heterogeneity,complex reflected wavefield and low signal-to-noise ratio for seismic data.In view of this,we created a high resolution reverse time migration imaging technical workflow for Yubei fractured-vuggy carbonate reservoir based on depth domain velocity model building and high fidelity weak signal compensation.We applied the technology on an area of 474km23D seismic data in Yubei 6 well region.Compared with the conventional prestack time migration,our technique technology has obvious advantages in the complicated structure imaging,“moniliform” imaging and fractured-vuggy reservoir characterization,and effectively improve the imaging precision of Ordovician fractured-vuggy carbonate target reservoir.

    reverse time migration,velocity model building,weak signal compensation,fractured-vuggy reservoir

    2014-12-28;改回日期:2015-04-15。

    李鵬(1979—),男,工程師,主要從事地震資料處理方法研究工作。

    P631

    A

    1000-1441(2015)04-0443-09

    10.3969/j.issn.1000-1441.2015.04.011

    猜你喜歡
    縫洞波場(chǎng)碳酸鹽巖
    適用于碳酸鹽巖儲(chǔ)層的固體酸解堵體系
    碳酸鹽巖裂縫描述七大難點(diǎn)
    碳酸鹽巖縫洞儲(chǔ)集體分尺度量化表征
    哈拉哈塘奧陶系縫洞型成巖圈閉及其成因
    大數(shù)據(jù)統(tǒng)計(jì)在碳酸鹽巖油氣田開發(fā)中的應(yīng)用
    彈性波波場(chǎng)分離方法對(duì)比及其在逆時(shí)偏移成像中的應(yīng)用
    交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析
    基于Hilbert變換的全波場(chǎng)分離逆時(shí)偏移成像
    考慮酸蝕蚓孔的碳酸鹽巖儲(chǔ)層酸壓產(chǎn)能預(yù)測(cè)
    縫洞型介質(zhì)結(jié)構(gòu)對(duì)非混相氣驅(qū)油采收率的影響
    久久久久视频综合| 黄色毛片三级朝国网站| 久久国产精品影院| tube8黄色片| 欧美日本中文国产一区发布| 香蕉国产在线看| 人人妻人人澡人人看| 日韩欧美在线二视频 | 午夜福利欧美成人| 亚洲熟妇熟女久久| 少妇裸体淫交视频免费看高清 | av超薄肉色丝袜交足视频| 国产极品粉嫩免费观看在线| 看片在线看免费视频| 在线观看免费视频日本深夜| 欧美老熟妇乱子伦牲交| 大香蕉久久成人网| xxx96com| 国产乱人伦免费视频| av中文乱码字幕在线| 黄色视频不卡| 99久久国产精品久久久| 亚洲午夜理论影院| 国产高清视频在线播放一区| 1024香蕉在线观看| 成人黄色视频免费在线看| 亚洲欧洲精品一区二区精品久久久| 757午夜福利合集在线观看| 亚洲成人免费av在线播放| 亚洲av日韩精品久久久久久密| 欧美日韩视频精品一区| xxxhd国产人妻xxx| 亚洲五月色婷婷综合| 国产淫语在线视频| 一级,二级,三级黄色视频| 久久久国产一区二区| 亚洲自偷自拍图片 自拍| 午夜免费成人在线视频| 亚洲av电影在线进入| 国产精品美女特级片免费视频播放器 | 日本欧美视频一区| cao死你这个sao货| 在线播放国产精品三级| 亚洲欧美激情在线| 亚洲av成人av| 亚洲人成电影免费在线| 一区二区三区国产精品乱码| 日韩欧美国产一区二区入口| 日韩视频一区二区在线观看| 少妇的丰满在线观看| 女人爽到高潮嗷嗷叫在线视频| 嫩草影视91久久| 成年人免费黄色播放视频| 欧美日韩视频精品一区| 99精品在免费线老司机午夜| 亚洲专区中文字幕在线| 俄罗斯特黄特色一大片| 亚洲免费av在线视频| 精品国产国语对白av| 99久久人妻综合| 成人亚洲精品一区在线观看| 久久久国产成人精品二区 | 无遮挡黄片免费观看| 免费高清在线观看日韩| 成年女人毛片免费观看观看9 | 久久精品人人爽人人爽视色| 国产97色在线日韩免费| 香蕉久久夜色| 搡老岳熟女国产| 狠狠婷婷综合久久久久久88av| 国产区一区二久久| 欧美不卡视频在线免费观看 | 夜夜爽天天搞| 黑人欧美特级aaaaaa片| 国产精品综合久久久久久久免费 | 999久久久国产精品视频| 欧美激情极品国产一区二区三区| 色婷婷av一区二区三区视频| 岛国在线观看网站| 亚洲成人国产一区在线观看| 久久性视频一级片| 久久国产亚洲av麻豆专区| 老熟妇仑乱视频hdxx| 亚洲熟妇熟女久久| 精品久久蜜臀av无| 亚洲午夜精品一区,二区,三区| 国产xxxxx性猛交| 18禁美女被吸乳视频| 久久久久国产一级毛片高清牌| 高清欧美精品videossex| 天天躁狠狠躁夜夜躁狠狠躁| 少妇裸体淫交视频免费看高清 | 国产片内射在线| 欧美精品亚洲一区二区| 欧美黄色片欧美黄色片| 亚洲熟女精品中文字幕| 丝袜美足系列| 欧美精品高潮呻吟av久久| cao死你这个sao货| 国产精品一区二区免费欧美| 久久精品国产a三级三级三级| 老司机影院毛片| 两人在一起打扑克的视频| 成人永久免费在线观看视频| 日韩欧美一区二区三区在线观看 | 99国产精品一区二区三区| 巨乳人妻的诱惑在线观看| 狠狠婷婷综合久久久久久88av| av网站在线播放免费| 亚洲色图av天堂| 91九色精品人成在线观看| 亚洲国产精品sss在线观看 | 色婷婷av一区二区三区视频| 欧美日韩av久久| 在线永久观看黄色视频| tocl精华| 亚洲国产欧美网| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一区二区三区欧美精品| 啦啦啦视频在线资源免费观看| 大型黄色视频在线免费观看| 免费女性裸体啪啪无遮挡网站| 精品久久蜜臀av无| 丰满饥渴人妻一区二区三| 国产欧美日韩精品亚洲av| 久久香蕉激情| 两性午夜刺激爽爽歪歪视频在线观看 | 一边摸一边做爽爽视频免费| 亚洲精品美女久久av网站| 黄色片一级片一级黄色片| 王馨瑶露胸无遮挡在线观看| 日韩免费av在线播放| 国产精品自产拍在线观看55亚洲 | 亚洲国产看品久久| 久久午夜亚洲精品久久| 精品久久久精品久久久| 国产欧美日韩一区二区三| 视频在线观看一区二区三区| 中文字幕av电影在线播放| 嫩草影视91久久| 狠狠狠狠99中文字幕| 日韩欧美在线二视频 | 中文字幕人妻熟女乱码| 成人永久免费在线观看视频| 成年动漫av网址| 桃红色精品国产亚洲av| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品sss在线观看 | 欧美丝袜亚洲另类 | 亚洲午夜理论影院| 亚洲精品一二三| 一进一出抽搐动态| 久久久久国产精品人妻aⅴ院 | 国产精品电影一区二区三区 | 午夜成年电影在线免费观看| 亚洲专区中文字幕在线| 亚洲精品美女久久av网站| 亚洲av日韩精品久久久久久密| 黑人猛操日本美女一级片| 美女扒开内裤让男人捅视频| 国产精品欧美亚洲77777| 看片在线看免费视频| 亚洲熟女毛片儿| av欧美777| 狠狠婷婷综合久久久久久88av| 一进一出抽搐动态| 日日夜夜操网爽| 欧美亚洲日本最大视频资源| 18禁国产床啪视频网站| 啪啪无遮挡十八禁网站| 亚洲精品乱久久久久久| 国产高清videossex| 欧美不卡视频在线免费观看 | 国产精品98久久久久久宅男小说| 天天影视国产精品| 狠狠婷婷综合久久久久久88av| 欧美黄色淫秽网站| 天天躁日日躁夜夜躁夜夜| 亚洲av成人不卡在线观看播放网| 19禁男女啪啪无遮挡网站| 在线观看日韩欧美| 亚洲精品久久午夜乱码| 日本黄色视频三级网站网址 | 免费久久久久久久精品成人欧美视频| 在线国产一区二区在线| 热re99久久精品国产66热6| 久久亚洲精品不卡| 国产av精品麻豆| 人人妻人人澡人人看| av国产精品久久久久影院| 水蜜桃什么品种好| a级毛片黄视频| 纯流量卡能插随身wifi吗| 亚洲成人免费av在线播放| 国产精品免费一区二区三区在线 | 亚洲精品国产一区二区精华液| 看片在线看免费视频| 精品一区二区三区视频在线观看免费 | 性少妇av在线| 日日爽夜夜爽网站| 国产精品1区2区在线观看. | 国产精品 国内视频| 久久这里只有精品19| 国产在线精品亚洲第一网站| 日本黄色视频三级网站网址 | 国产99久久九九免费精品| 久久久精品区二区三区| av欧美777| 最近最新免费中文字幕在线| 国产精品永久免费网站| 在线观看66精品国产| 久热这里只有精品99| 久久99一区二区三区| 夜夜夜夜夜久久久久| 黄色丝袜av网址大全| 欧美国产精品va在线观看不卡| 精品人妻1区二区| 欧美日韩中文字幕国产精品一区二区三区 | 黄频高清免费视频| 一级a爱片免费观看的视频| 久久久精品免费免费高清| 夜夜爽天天搞| 黄色丝袜av网址大全| 国产亚洲精品久久久久久毛片 | 丰满人妻熟妇乱又伦精品不卡| 老司机在亚洲福利影院| 午夜福利一区二区在线看| 午夜福利在线观看吧| 久久ye,这里只有精品| 亚洲熟女毛片儿| 亚洲片人在线观看| av超薄肉色丝袜交足视频| 日韩欧美在线二视频 | 成人三级做爰电影| 亚洲专区中文字幕在线| tocl精华| 侵犯人妻中文字幕一二三四区| 久久精品人人爽人人爽视色| 欧美大码av| 精品亚洲成a人片在线观看| 高清毛片免费观看视频网站 | 下体分泌物呈黄色| 俄罗斯特黄特色一大片| 黄色a级毛片大全视频| 一级a爱片免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 下体分泌物呈黄色| 一夜夜www| 日日爽夜夜爽网站| 99国产精品免费福利视频| 欧美一级毛片孕妇| 99国产极品粉嫩在线观看| 国产成人影院久久av| 国产在视频线精品| 免费高清在线观看日韩| 黄色毛片三级朝国网站| 亚洲第一青青草原| 亚洲综合色网址| 久久久久久久国产电影| 青草久久国产| 国产精品1区2区在线观看. | 真人做人爱边吃奶动态| 欧美黑人欧美精品刺激| 一级片'在线观看视频| 中文字幕制服av| 18在线观看网站| 亚洲一区二区三区不卡视频| 美女高潮到喷水免费观看| 国产人伦9x9x在线观看| 久久精品国产亚洲av香蕉五月 | 日韩欧美一区二区三区在线观看 | 99久久人妻综合| 久久性视频一级片| 午夜两性在线视频| 国产成人精品久久二区二区91| 啦啦啦视频在线资源免费观看| 国产不卡av网站在线观看| 精品久久蜜臀av无| 天天躁狠狠躁夜夜躁狠狠躁| 午夜老司机福利片| 国产成人欧美| 嫁个100分男人电影在线观看| 久久精品亚洲熟妇少妇任你| 免费观看人在逋| 久久久久久久久免费视频了| 在线观看免费视频网站a站| 中文字幕人妻熟女乱码| 天堂俺去俺来也www色官网| 18禁美女被吸乳视频| 亚洲国产欧美网| 天堂√8在线中文| 亚洲欧洲精品一区二区精品久久久| 午夜免费成人在线视频| 99国产精品一区二区蜜桃av | 亚洲免费av在线视频| 69精品国产乱码久久久| 免费久久久久久久精品成人欧美视频| 又黄又爽又免费观看的视频| 国产主播在线观看一区二区| 欧美乱妇无乱码| 在线观看舔阴道视频| 欧美老熟妇乱子伦牲交| 18禁裸乳无遮挡免费网站照片 | 高清毛片免费观看视频网站 | 成人18禁在线播放| 丝袜在线中文字幕| 91av网站免费观看| 欧洲精品卡2卡3卡4卡5卡区| 天堂俺去俺来也www色官网| 露出奶头的视频| 亚洲avbb在线观看| 动漫黄色视频在线观看| 老司机影院毛片| 国产精品 欧美亚洲| 精品国产超薄肉色丝袜足j| 视频区欧美日本亚洲| 色综合欧美亚洲国产小说| 国产一区二区三区视频了| 中文字幕另类日韩欧美亚洲嫩草| 多毛熟女@视频| 免费黄频网站在线观看国产| 女人被狂操c到高潮| www.自偷自拍.com| 视频区图区小说| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩视频精品一区| 老司机靠b影院| 十八禁人妻一区二区| 9热在线视频观看99| 狠狠婷婷综合久久久久久88av| 9191精品国产免费久久| 精品久久久久久久毛片微露脸| 精品少妇久久久久久888优播| 黄片小视频在线播放| 午夜成年电影在线免费观看| 在线永久观看黄色视频| 91大片在线观看| 99re在线观看精品视频| 午夜久久久在线观看| tube8黄色片| 亚洲免费av在线视频| 精品人妻1区二区| 久久这里只有精品19| 亚洲专区中文字幕在线| 精品久久久久久电影网| 国产在线一区二区三区精| 在线观看免费午夜福利视频| 国产精品综合久久久久久久免费 | 成人免费观看视频高清| 亚洲av片天天在线观看| 91成年电影在线观看| 亚洲 欧美一区二区三区| 久久精品国产清高在天天线| xxx96com| 亚洲男人天堂网一区| 精品第一国产精品| 男女床上黄色一级片免费看| 怎么达到女性高潮| 国产一区二区三区综合在线观看| 一级毛片女人18水好多| 久久久水蜜桃国产精品网| 我的亚洲天堂| 人妻一区二区av| 黄片播放在线免费| 国产午夜精品久久久久久| 夜夜躁狠狠躁天天躁| 亚洲精品美女久久久久99蜜臀| 精品国产乱子伦一区二区三区| 日韩欧美三级三区| 国产精品电影一区二区三区 | 色婷婷久久久亚洲欧美| 一级片免费观看大全| 色老头精品视频在线观看| 亚洲欧美一区二区三区黑人| 黑人猛操日本美女一级片| 80岁老熟妇乱子伦牲交| 久久久久精品国产欧美久久久| 他把我摸到了高潮在线观看| 91大片在线观看| 人妻一区二区av| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久av网站| 国精品久久久久久国模美| 亚洲精品美女久久av网站| 一级毛片精品| 国产精品一区二区在线观看99| 一级毛片高清免费大全| 成在线人永久免费视频| 亚洲一区二区三区不卡视频| 男女下面插进去视频免费观看| 亚洲色图av天堂| 久久精品91无色码中文字幕| 又紧又爽又黄一区二区| 成人手机av| 亚洲精品成人av观看孕妇| 一级片免费观看大全| 女同久久另类99精品国产91| 国产精品成人在线| 12—13女人毛片做爰片一| 日本一区二区免费在线视频| 午夜免费成人在线视频| 可以免费在线观看a视频的电影网站| 巨乳人妻的诱惑在线观看| 亚洲精品国产一区二区精华液| 久久久精品免费免费高清| 中文字幕最新亚洲高清| 91精品三级在线观看| 日日爽夜夜爽网站| 黄色视频,在线免费观看| 咕卡用的链子| 亚洲欧美一区二区三区黑人| 侵犯人妻中文字幕一二三四区| 桃红色精品国产亚洲av| 午夜精品国产一区二区电影| 精品久久久久久电影网| 国产一区在线观看成人免费| 亚洲精品国产一区二区精华液| 少妇被粗大的猛进出69影院| 中文字幕最新亚洲高清| 亚洲精品粉嫩美女一区| 亚洲欧美一区二区三区久久| 精品国产乱码久久久久久男人| 又黄又粗又硬又大视频| 色精品久久人妻99蜜桃| 国产99白浆流出| 69精品国产乱码久久久| 国产野战对白在线观看| 中文字幕色久视频| 久久天堂一区二区三区四区| 国精品久久久久久国模美| 国产男女超爽视频在线观看| 国产精品乱码一区二三区的特点 | 日本a在线网址| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| 国产男靠女视频免费网站| 91麻豆av在线| 黄片小视频在线播放| 国精品久久久久久国模美| 免费高清在线观看日韩| 下体分泌物呈黄色| 亚洲人成77777在线视频| 精品久久久精品久久久| 亚洲国产欧美一区二区综合| 激情视频va一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲国产毛片av蜜桃av| 久久久国产一区二区| 91麻豆精品激情在线观看国产 | 亚洲全国av大片| 99热网站在线观看| 久久久精品免费免费高清| 黄片播放在线免费| 黄色视频,在线免费观看| a级毛片黄视频| 欧美日韩成人在线一区二区| 久久久久久久国产电影| 777米奇影视久久| 午夜免费观看网址| 啦啦啦在线免费观看视频4| 老司机影院毛片| 亚洲一区中文字幕在线| 亚洲五月婷婷丁香| 大型黄色视频在线免费观看| 色综合欧美亚洲国产小说| a级片在线免费高清观看视频| xxxhd国产人妻xxx| 国产亚洲精品久久久久5区| 黄色成人免费大全| 在线观看免费午夜福利视频| 日本欧美视频一区| 淫妇啪啪啪对白视频| 国产男女超爽视频在线观看| 欧美精品亚洲一区二区| 中文字幕人妻熟女乱码| 啦啦啦免费观看视频1| 色尼玛亚洲综合影院| 99在线人妻在线中文字幕 | 99riav亚洲国产免费| 12—13女人毛片做爰片一| 欧美黄色淫秽网站| 国产精品免费一区二区三区在线 | 深夜精品福利| 好看av亚洲va欧美ⅴa在| 国产精品久久视频播放| 欧美亚洲日本最大视频资源| 亚洲精品乱久久久久久| 午夜福利,免费看| 一级毛片高清免费大全| 啪啪无遮挡十八禁网站| 中文字幕制服av| 久久久久国产一级毛片高清牌| a级片在线免费高清观看视频| 成人精品一区二区免费| 欧美精品人与动牲交sv欧美| 美女视频免费永久观看网站| 欧美丝袜亚洲另类 | av在线播放免费不卡| 久久人妻av系列| 日韩视频一区二区在线观看| 侵犯人妻中文字幕一二三四区| 日韩有码中文字幕| 久久精品熟女亚洲av麻豆精品| 精品第一国产精品| 久久中文字幕一级| 中文亚洲av片在线观看爽 | avwww免费| 免费观看精品视频网站| 免费女性裸体啪啪无遮挡网站| 亚洲三区欧美一区| 国产97色在线日韩免费| 黄片小视频在线播放| 欧美日韩国产mv在线观看视频| 久久精品91无色码中文字幕| 丝瓜视频免费看黄片| 校园春色视频在线观看| 久久香蕉国产精品| 欧美日韩精品网址| 黄色片一级片一级黄色片| 在线观看免费视频网站a站| 脱女人内裤的视频| 国产成人精品无人区| 亚洲国产欧美日韩在线播放| 婷婷丁香在线五月| 亚洲熟妇中文字幕五十中出 | 免费少妇av软件| 9色porny在线观看| 欧美日韩亚洲综合一区二区三区_| 一级黄色大片毛片| 久久久精品国产亚洲av高清涩受| 国产精品久久久人人做人人爽| 一级毛片精品| 久久人妻福利社区极品人妻图片| 99国产精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 两个人看的免费小视频| 91精品三级在线观看| 怎么达到女性高潮| 五月开心婷婷网| 亚洲国产欧美网| 老司机影院毛片| 又黄又粗又硬又大视频| 在线看a的网站| 水蜜桃什么品种好| 中国美女看黄片| 亚洲熟女毛片儿| tocl精华| 欧美+亚洲+日韩+国产| 老司机午夜十八禁免费视频| 男女午夜视频在线观看| 91麻豆av在线| 热99久久久久精品小说推荐| 成人18禁在线播放| 欧美黄色片欧美黄色片| 啦啦啦在线免费观看视频4| 日韩视频一区二区在线观看| 国产精品久久久久成人av| 性色av乱码一区二区三区2| 亚洲午夜精品一区,二区,三区| 麻豆av在线久日| 日本黄色视频三级网站网址 | 三上悠亚av全集在线观看| 国产国语露脸激情在线看| 欧美精品高潮呻吟av久久| 亚洲精品久久午夜乱码| 亚洲黑人精品在线| 亚洲伊人色综图| 久久ye,这里只有精品| 91av网站免费观看| 免费在线观看影片大全网站| www.精华液| 久久精品人人爽人人爽视色| 80岁老熟妇乱子伦牲交| 色播在线永久视频| 久久性视频一级片| 精品国产一区二区三区四区第35| 99精品在免费线老司机午夜| a级毛片在线看网站| 免费看a级黄色片| 久久ye,这里只有精品| 午夜亚洲福利在线播放| 日本黄色视频三级网站网址 | 男女床上黄色一级片免费看| 久久久久国产精品人妻aⅴ院 | 日本五十路高清| 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲综合一区二区三区_| 婷婷成人精品国产| 高清av免费在线| 久久国产精品大桥未久av| 精品一区二区三区视频在线观看免费 | 欧美精品一区二区免费开放| 搡老岳熟女国产| 一区二区三区国产精品乱码| 欧美日韩成人在线一区二区| 久久国产精品男人的天堂亚洲| 99国产精品99久久久久| 老鸭窝网址在线观看| 中文字幕制服av| 搡老乐熟女国产| 叶爱在线成人免费视频播放| 黄网站色视频无遮挡免费观看| 国产成人欧美在线观看 | 欧美性长视频在线观看| 69精品国产乱码久久久| 他把我摸到了高潮在线观看| 欧美 亚洲 国产 日韩一| 亚洲熟女精品中文字幕| 人妻 亚洲 视频| 亚洲精品久久午夜乱码| 丝袜美腿诱惑在线| 欧美不卡视频在线免费观看 | 极品教师在线免费播放| 悠悠久久av| 成人影院久久|