• <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ū)油采收率的影響
    国产亚洲一区二区精品| 高清欧美精品videossex| 久久久久精品性色| 一级爰片在线观看| 亚洲色图av天堂| 91午夜精品亚洲一区二区三区| 免费播放大片免费观看视频在线观看| 在线观看一区二区三区| 综合色丁香网| 欧美性感艳星| 成人漫画全彩无遮挡| 99九九线精品视频在线观看视频| av卡一久久| 97超视频在线观看视频| 99久久中文字幕三级久久日本| 男人舔奶头视频| 久久久亚洲精品成人影院| 午夜日本视频在线| 一区二区三区精品91| 亚洲精品自拍成人| 国产精品福利在线免费观看| 久久国产乱子免费精品| 哪个播放器可以免费观看大片| 自拍偷自拍亚洲精品老妇| 精华霜和精华液先用哪个| 黄片无遮挡物在线观看| 色网站视频免费| 看免费成人av毛片| 亚洲无线观看免费| 欧美 日韩 精品 国产| 国产淫片久久久久久久久| 夫妻午夜视频| 精品久久久久久久久av| 国内揄拍国产精品人妻在线| 国产伦精品一区二区三区视频9| 亚洲综合色惰| 国产成人免费观看mmmm| 一本色道久久久久久精品综合| 欧美国产精品一级二级三级 | 一本色道久久久久久精品综合| 大码成人一级视频| 99热网站在线观看| 国产午夜精品久久久久久一区二区三区| 成人免费观看视频高清| 亚洲国产日韩一区二区| 日韩电影二区| 在线观看免费视频网站a站| 久久毛片免费看一区二区三区| 熟女av电影| 成人高潮视频无遮挡免费网站| 黑丝袜美女国产一区| 天堂中文最新版在线下载| 国产女主播在线喷水免费视频网站| 交换朋友夫妻互换小说| 自拍偷自拍亚洲精品老妇| 色哟哟·www| 男人爽女人下面视频在线观看| 在线观看一区二区三区激情| 在线天堂最新版资源| 精品一区二区免费观看| 97热精品久久久久久| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 男女边吃奶边做爰视频| 亚洲aⅴ乱码一区二区在线播放| 久久久a久久爽久久v久久| 91精品国产国语对白视频| 国精品久久久久久国模美| 最后的刺客免费高清国语| 下体分泌物呈黄色| 热re99久久精品国产66热6| 亚洲av不卡在线观看| 精品亚洲乱码少妇综合久久| av线在线观看网站| 久久毛片免费看一区二区三区| av福利片在线观看| 色吧在线观看| 国产黄色视频一区二区在线观看| 国产在线男女| 日本色播在线视频| 亚洲婷婷狠狠爱综合网| 国产免费又黄又爽又色| 只有这里有精品99| 美女福利国产在线 | 欧美人与善性xxx| 欧美激情极品国产一区二区三区 | 2018国产大陆天天弄谢| 色吧在线观看| 一级a做视频免费观看| 丰满迷人的少妇在线观看| 超碰av人人做人人爽久久| 国模一区二区三区四区视频| 最近的中文字幕免费完整| 男人狂女人下面高潮的视频| 国产精品久久久久久精品古装| av女优亚洲男人天堂| 精品一区二区三区视频在线| 夜夜骑夜夜射夜夜干| 乱系列少妇在线播放| 天堂俺去俺来也www色官网| 日日啪夜夜爽| 少妇的逼水好多| 少妇高潮的动态图| 免费观看av网站的网址| 亚洲精品久久久久久婷婷小说| 99久久精品热视频| 国产成人精品婷婷| 大片免费播放器 马上看| 2021少妇久久久久久久久久久| 国产精品一区二区在线不卡| 日韩三级伦理在线观看| 菩萨蛮人人尽说江南好唐韦庄| 成人国产av品久久久| 身体一侧抽搐| 哪个播放器可以免费观看大片| 国产精品人妻久久久影院| 在线观看国产h片| 国产av国产精品国产| 亚洲人成网站在线播| 成人亚洲欧美一区二区av| 中国三级夫妇交换| a级毛色黄片| 不卡视频在线观看欧美| 97超视频在线观看视频| 欧美丝袜亚洲另类| 中文资源天堂在线| 国产色爽女视频免费观看| 一本色道久久久久久精品综合| 这个男人来自地球电影免费观看 | 国产精品秋霞免费鲁丝片| 777米奇影视久久| 校园人妻丝袜中文字幕| 毛片女人毛片| 久久精品国产亚洲网站| 中文字幕人妻熟人妻熟丝袜美| 国产精品不卡视频一区二区| 欧美日韩一区二区视频在线观看视频在线| 日日啪夜夜撸| 99久久中文字幕三级久久日本| 亚洲婷婷狠狠爱综合网| 观看免费一级毛片| 男人舔奶头视频| 最近中文字幕高清免费大全6| 国产深夜福利视频在线观看| 我的老师免费观看完整版| 免费av中文字幕在线| 久热久热在线精品观看| 身体一侧抽搐| 91精品国产国语对白视频| 亚洲av成人精品一区久久| 22中文网久久字幕| 精品少妇久久久久久888优播| 国产91av在线免费观看| 亚洲精品乱码久久久久久按摩| 国产综合精华液| 内地一区二区视频在线| 这个男人来自地球电影免费观看 | 一级毛片久久久久久久久女| 嘟嘟电影网在线观看| 国产男女超爽视频在线观看| 在线观看三级黄色| 欧美成人精品欧美一级黄| 亚洲国产高清在线一区二区三| 伦理电影免费视频| 国产爽快片一区二区三区| 免费观看av网站的网址| 一区二区三区四区激情视频| 国产精品久久久久久久久免| 国产精品一区www在线观看| 91aial.com中文字幕在线观看| 国产 一区精品| 一级毛片 在线播放| 国产av国产精品国产| 麻豆乱淫一区二区| 国产成人a区在线观看| 亚洲中文av在线| 亚洲av国产av综合av卡| 久久久久视频综合| 久久精品久久精品一区二区三区| 干丝袜人妻中文字幕| 小蜜桃在线观看免费完整版高清| av国产免费在线观看| 中文字幕精品免费在线观看视频 | 亚洲内射少妇av| 婷婷色av中文字幕| 制服丝袜香蕉在线| 毛片一级片免费看久久久久| 美女中出高潮动态图| av在线观看视频网站免费| 午夜福利视频精品| 国产毛片在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美激情国产日韩精品一区| 十分钟在线观看高清视频www | 丝瓜视频免费看黄片| 国产一区二区三区综合在线观看 | 国产精品蜜桃在线观看| 久久久久久久久久久丰满| 人人妻人人爽人人添夜夜欢视频 | 国产av一区二区精品久久 | 精品午夜福利在线看| 建设人人有责人人尽责人人享有的 | 人人妻人人澡人人爽人人夜夜| 亚洲国产欧美在线一区| 成年av动漫网址| 国产高潮美女av| 精品少妇久久久久久888优播| 在线观看一区二区三区激情| 哪个播放器可以免费观看大片| 日本欧美国产在线视频| 看免费成人av毛片| av一本久久久久| 亚洲美女搞黄在线观看| 国产精品福利在线免费观看| 日本欧美国产在线视频| 国产在线免费精品| 国产精品精品国产色婷婷| 精品人妻一区二区三区麻豆| 日本黄大片高清| 欧美一区二区亚洲| 免费人成在线观看视频色| 精品久久久久久久久亚洲| 亚洲欧美清纯卡通| 99re6热这里在线精品视频| 国产精品人妻久久久影院| 蜜桃亚洲精品一区二区三区| 国产精品一区二区性色av| 嫩草影院新地址| 男女无遮挡免费网站观看| 精品国产乱码久久久久久小说| 成年人午夜在线观看视频| 亚洲国产日韩一区二区| 亚洲欧洲日产国产| 久久久a久久爽久久v久久| 欧美日韩视频精品一区| 最近的中文字幕免费完整| 大陆偷拍与自拍| 久久人人爽av亚洲精品天堂 | 我要看日韩黄色一级片| 久久精品熟女亚洲av麻豆精品| 国产真实伦视频高清在线观看| 99久国产av精品国产电影| 伊人久久精品亚洲午夜| 国产高清国产精品国产三级 | 激情 狠狠 欧美| videos熟女内射| 亚洲欧美精品专区久久| 水蜜桃什么品种好| 99精国产麻豆久久婷婷| 亚洲色图综合在线观看| 中国美白少妇内射xxxbb| 大码成人一级视频| 国产亚洲av片在线观看秒播厂| 国产真实伦视频高清在线观看| 在线免费十八禁| 老司机影院毛片| 天堂8中文在线网| 在线观看免费视频网站a站| 少妇人妻久久综合中文| 女的被弄到高潮叫床怎么办| 国产乱人视频| 美女cb高潮喷水在线观看| 欧美激情国产日韩精品一区| 欧美少妇被猛烈插入视频| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩视频精品一区| 性色av一级| 国产女主播在线喷水免费视频网站| 天天躁日日操中文字幕| 国产在线视频一区二区| 丰满迷人的少妇在线观看| 久久国产精品男人的天堂亚洲 | 赤兔流量卡办理| 亚洲精品乱久久久久久| videos熟女内射| 国产一区二区在线观看日韩| 色视频www国产| 边亲边吃奶的免费视频| 99久久精品国产国产毛片| 九草在线视频观看| 人人妻人人添人人爽欧美一区卜 | 爱豆传媒免费全集在线观看| 九九爱精品视频在线观看| 久久99精品国语久久久| 一级黄片播放器| 亚洲欧美精品自产自拍| 久久国内精品自在自线图片| 成人免费观看视频高清| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 亚洲国产欧美人成| tube8黄色片| 最新中文字幕久久久久| 如何舔出高潮| 欧美日韩精品成人综合77777| 免费在线观看成人毛片| 日韩av在线免费看完整版不卡| 日本欧美国产在线视频| 成年美女黄网站色视频大全免费 | 一二三四中文在线观看免费高清| 18+在线观看网站| 亚洲国产欧美在线一区| 亚洲第一av免费看| 777米奇影视久久| 涩涩av久久男人的天堂| 国产精品一及| a级毛色黄片| 亚洲经典国产精华液单| 国产精品福利在线免费观看| 久久久久久久久久人人人人人人| 2022亚洲国产成人精品| 少妇裸体淫交视频免费看高清| 欧美精品亚洲一区二区| 久久久精品免费免费高清| 日日啪夜夜撸| 免费观看性生交大片5| 国产在线一区二区三区精| 两个人的视频大全免费| 91在线精品国自产拍蜜月| av在线老鸭窝| 亚洲av国产av综合av卡| 成人国产av品久久久| 妹子高潮喷水视频| 欧美精品一区二区免费开放| 99热这里只有是精品在线观看| 两个人的视频大全免费| 老司机影院成人| 免费观看av网站的网址| 亚洲av免费高清在线观看| 熟女电影av网| 亚洲四区av| 在线天堂最新版资源| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 亚洲,一卡二卡三卡| 人妻一区二区av| 久久久久人妻精品一区果冻| 另类亚洲欧美激情| 成人18禁高潮啪啪吃奶动态图 | 99热全是精品| 免费人妻精品一区二区三区视频| 一级黄片播放器| 日韩av在线免费看完整版不卡| 久久久久精品性色| 国产乱人视频| 三级国产精品片| freevideosex欧美| 日本黄色片子视频| 成人黄色视频免费在线看| 久久99热6这里只有精品| 2022亚洲国产成人精品| 国产中年淑女户外野战色| 日韩一本色道免费dvd| 亚洲国产精品成人久久小说| 纵有疾风起免费观看全集完整版| 丝袜脚勾引网站| 两个人的视频大全免费| 最近2019中文字幕mv第一页| 22中文网久久字幕| 韩国av在线不卡| 亚洲精品日韩在线中文字幕| 麻豆精品久久久久久蜜桃| 亚洲精品日韩在线中文字幕| 国产精品国产三级国产专区5o| 男女啪啪激烈高潮av片| 久久精品国产亚洲网站| 欧美日韩视频精品一区| 中文字幕免费在线视频6| 五月伊人婷婷丁香| 午夜精品国产一区二区电影| 男女下面进入的视频免费午夜| 晚上一个人看的免费电影| 成人黄色视频免费在线看| 欧美丝袜亚洲另类| 日本午夜av视频| 美女视频免费永久观看网站| 国产一区二区在线观看日韩| 亚洲不卡免费看| 精品久久久久久电影网| tube8黄色片| 国产一区亚洲一区在线观看| 免费少妇av软件| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 免费av中文字幕在线| 国产午夜精品一二区理论片| 国产高清国产精品国产三级 | 欧美日韩国产mv在线观看视频 | 韩国av在线不卡| 一区二区三区精品91| 妹子高潮喷水视频| 免费播放大片免费观看视频在线观看| 亚洲国产精品国产精品| 大码成人一级视频| 亚洲av中文字字幕乱码综合| 亚洲av电影在线观看一区二区三区| 蜜桃在线观看..| 搡老乐熟女国产| 伦精品一区二区三区| 免费观看性生交大片5| 国产成人精品福利久久| 精品久久久精品久久久| 亚洲av日韩在线播放| 国产人妻一区二区三区在| 亚洲一区二区三区欧美精品| 欧美3d第一页| 少妇人妻 视频| 高清在线视频一区二区三区| 一级毛片 在线播放| 亚洲av成人精品一区久久| 久久人人爽人人爽人人片va| 日韩一区二区视频免费看| 欧美另类一区| 内射极品少妇av片p| 在线观看美女被高潮喷水网站| 国产欧美亚洲国产| 王馨瑶露胸无遮挡在线观看| 欧美xxⅹ黑人| 韩国av在线不卡| 亚洲欧美日韩无卡精品| 内地一区二区视频在线| av又黄又爽大尺度在线免费看| 成年女人在线观看亚洲视频| 噜噜噜噜噜久久久久久91| 国产精品99久久久久久久久| 国产精品久久久久久av不卡| 免费观看在线日韩| 亚洲精品第二区| 深夜a级毛片| 内地一区二区视频在线| av又黄又爽大尺度在线免费看| 日韩一区二区视频免费看| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品乱码久久久v下载方式| 亚洲av男天堂| 国产淫语在线视频| 舔av片在线| 国产精品99久久久久久久久| 精品久久久久久久末码| 国产黄频视频在线观看| 最近的中文字幕免费完整| 国产精品福利在线免费观看| 国产久久久一区二区三区| 97超视频在线观看视频| 午夜免费男女啪啪视频观看| 国产黄片视频在线免费观看| 99精国产麻豆久久婷婷| 国产成人aa在线观看| 国产一区二区三区av在线| 内地一区二区视频在线| 又粗又硬又长又爽又黄的视频| a级一级毛片免费在线观看| 舔av片在线| 大香蕉久久网| 老司机影院成人| 久久久午夜欧美精品| 特大巨黑吊av在线直播| 亚州av有码| 国产午夜精品一二区理论片| 日日摸夜夜添夜夜爱| 日本免费在线观看一区| 最近的中文字幕免费完整| 日本色播在线视频| tube8黄色片| av一本久久久久| 亚洲最大成人中文| 一本久久精品| 男女边摸边吃奶| 女的被弄到高潮叫床怎么办| 国产欧美另类精品又又久久亚洲欧美| 亚洲精华国产精华液的使用体验| 女人十人毛片免费观看3o分钟| 大陆偷拍与自拍| 精品人妻熟女av久视频| 美女主播在线视频| 亚洲人成网站高清观看| 欧美性感艳星| 欧美区成人在线视频| 久久久久精品性色| 亚洲成人中文字幕在线播放| 亚洲人成网站高清观看| 大码成人一级视频| 国产精品麻豆人妻色哟哟久久| 日韩大片免费观看网站| 丝瓜视频免费看黄片| 欧美精品一区二区大全| 大香蕉97超碰在线| 久久国产精品大桥未久av | 亚洲av综合色区一区| 99热网站在线观看| 国产精品国产av在线观看| 高清不卡的av网站| 亚洲欧美精品自产自拍| 干丝袜人妻中文字幕| 久久久久久久久大av| av视频免费观看在线观看| 黄色怎么调成土黄色| 波野结衣二区三区在线| 国产精品一区www在线观看| 亚洲欧美日韩无卡精品| 又黄又爽又刺激的免费视频.| a级毛片免费高清观看在线播放| 成人高潮视频无遮挡免费网站| 久久久久久伊人网av| 赤兔流量卡办理| 在线观看免费高清a一片| 1000部很黄的大片| 晚上一个人看的免费电影| 国产欧美另类精品又又久久亚洲欧美| 国产伦精品一区二区三区四那| 午夜免费鲁丝| 日本-黄色视频高清免费观看| 国产亚洲一区二区精品| 十分钟在线观看高清视频www | av福利片在线观看| 少妇 在线观看| 搡女人真爽免费视频火全软件| 精品一品国产午夜福利视频| 日日摸夜夜添夜夜爱| 国产女主播在线喷水免费视频网站| 一级av片app| 亚洲精品色激情综合| 高清视频免费观看一区二区| 99久久精品国产国产毛片| 日本与韩国留学比较| 老女人水多毛片| 美女xxoo啪啪120秒动态图| 国产高清不卡午夜福利| 三级国产精品片| 高清日韩中文字幕在线| 日韩在线高清观看一区二区三区| 日韩欧美精品免费久久| 国产精品偷伦视频观看了| 最近最新中文字幕大全电影3| 丰满迷人的少妇在线观看| 美女脱内裤让男人舔精品视频| 韩国av在线不卡| 日日撸夜夜添| 久久亚洲国产成人精品v| 亚洲va在线va天堂va国产| 自拍偷自拍亚洲精品老妇| videos熟女内射| 18禁在线播放成人免费| 久久人妻熟女aⅴ| 国产精品久久久久久精品古装| 亚洲精品色激情综合| 国产探花极品一区二区| 精品一区二区三区视频在线| 久久精品久久久久久噜噜老黄| 亚洲丝袜综合中文字幕| 天天躁日日操中文字幕| 国产黄色免费在线视频| 乱码一卡2卡4卡精品| 深爱激情五月婷婷| 亚洲欧美精品自产自拍| 中文字幕亚洲精品专区| 日韩av免费高清视频| 性高湖久久久久久久久免费观看| 国产成人aa在线观看| 久久久色成人| 2022亚洲国产成人精品| 特大巨黑吊av在线直播| 国产精品不卡视频一区二区| 亚洲av福利一区| 久久这里有精品视频免费| 久久精品久久精品一区二区三区| av.在线天堂| 又粗又硬又长又爽又黄的视频| 中文资源天堂在线| 久久6这里有精品| 久久久久国产网址| 天天躁日日操中文字幕| 国产免费一级a男人的天堂| 亚洲精品第二区| 哪个播放器可以免费观看大片| 熟妇人妻不卡中文字幕| 国产一区二区在线观看日韩| 国产亚洲精品久久久com| 九草在线视频观看| 2021少妇久久久久久久久久久| 两个人的视频大全免费| 婷婷色av中文字幕| 亚洲国产av新网站| 亚洲三级黄色毛片| 黄色怎么调成土黄色| 一级毛片黄色毛片免费观看视频| 51国产日韩欧美| 看非洲黑人一级黄片| 国产 一区 欧美 日韩| 91久久精品国产一区二区三区| 熟女av电影| 黑人猛操日本美女一级片| 99re6热这里在线精品视频| 六月丁香七月| 身体一侧抽搐| 三级国产精品片| 亚洲av男天堂| 国产一区二区三区综合在线观看 | 亚洲欧美精品自产自拍| 99热这里只有精品一区| 国产视频首页在线观看| av.在线天堂| 99热这里只有精品一区| 偷拍熟女少妇极品色| 91久久精品电影网| 尤物成人国产欧美一区二区三区| 久久综合国产亚洲精品| 久久久久久九九精品二区国产| 久久精品国产鲁丝片午夜精品| 黄色怎么调成土黄色| 男女啪啪激烈高潮av片| 99热这里只有是精品50| 天天躁日日操中文字幕| 欧美日韩在线观看h| 亚洲丝袜综合中文字幕|