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

    開發(fā)地震技術(shù)在蘇里格氣田的應(yīng)用

    2016-09-23 08:14:00龐崇友李大軍張亞東武鐵嶺
    物探化探計(jì)算技術(shù) 2016年4期
    關(guān)鍵詞:里格橫波砂層

    龐崇友 , 李大軍 , 章 雄 , 張亞東 , 夏 銘, 武鐵嶺

    (1.中國石油集團(tuán) 川慶鉆探工程有限公司地球物理勘探公司,成都 610213;2.中國石油 長慶油田分公司勘探部,西安 710018)

    ?

    開發(fā)地震技術(shù)在蘇里格氣田的應(yīng)用

    龐崇友1, 李大軍1, 章雄1, 張亞東2, 夏銘1, 武鐵嶺2

    (1.中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司,成都610213;2.中國石油長慶油田分公司勘探部,西安710018)

    蘇里格氣田盒8氣藏屬低孔低滲氣藏,已有的地震預(yù)測方法與技術(shù)已難以滿足當(dāng)前水平井開發(fā)的需求。針對(duì)這一現(xiàn)狀,形成了以地震屬性分析與巖石物理分析技術(shù)為基礎(chǔ)、以高分辨率疊前反演技術(shù)與高精度時(shí)深轉(zhuǎn)換技術(shù)為核心的開發(fā)地震預(yù)測技術(shù)系列。首先,通過地震屬性分析宏觀預(yù)測主砂帶的展布,進(jìn)行橫波速度曲線預(yù)測與巖石物理分析;在此基礎(chǔ)上,開展高分辨率疊前反演、高精度時(shí)深轉(zhuǎn)換,準(zhǔn)確落實(shí)有效砂層的空間分布及其深度。應(yīng)用結(jié)果表明,均方根振幅屬性能宏觀反映盒8主砂帶的空間展布;Xu-White模型法是預(yù)測橫波速度曲線的有效方法,泊松比是準(zhǔn)確識(shí)別砂巖、含氣砂巖與泥巖的敏感參數(shù);高分辨率疊前反演技術(shù)極大地提高了預(yù)測結(jié)果的縱向分辨率,有效解決了盒8段砂泥巖薄互層的單砂體地震識(shí)別難題;高精度時(shí)深轉(zhuǎn)換技術(shù)與三維可視化技術(shù)相結(jié)合,為水平井位的部署、井軌跡的設(shè)計(jì)提供了直觀、可靠的依據(jù)。這一技術(shù)系列是陸相致密碎屑巖性油氣藏進(jìn)行儲(chǔ)層精細(xì)描述與高效開發(fā)的有效技術(shù)系列。

    開發(fā)地震技術(shù); 蘇里格氣田; 主砂帶; 橫波曲線預(yù)測; 巖石物理分析

    0 引言

    開發(fā)地震是在勘探地震的基礎(chǔ)上,充分利用針對(duì)油氣藏的觀測方法和信息處理技術(shù),緊密結(jié)合鉆井、測井、巖石物理、油田地質(zhì)和油藏工程等多學(xué)科資料,對(duì)油氣藏特征進(jìn)行橫向預(yù)測,在油氣田開發(fā)和開采過程中對(duì)油氣藏做出完整描述和進(jìn)行動(dòng)態(tài)監(jiān)測的一門新興學(xué)科,其主要任務(wù)包括:①確定油氣藏的圈閉形態(tài)與分布范圍、斷裂展布;②精細(xì)描述油藏參數(shù),包括不同巖性體的空間展布特征及其連通性、儲(chǔ)層的厚度、孔隙度與滲透率以及含油氣飽和度、孔隙流體壓力等;③監(jiān)測水驅(qū)、氣驅(qū)、注水、壓裂等增產(chǎn)措施的實(shí)施效果,完善油氣藏地質(zhì)模型,優(yōu)化井位部署和開發(fā)方案[1]。我國的開發(fā)地震自20世紀(jì)60年代末提出以來,受到了油氣工作者的廣泛關(guān)注并為處于開發(fā)后期老油氣田的剩余油氣開采與提高采收率發(fā)揮了重要的作用[2-10]。近年來,隨著蘇里格氣田等大批大中型低滲透油氣田的陸續(xù)發(fā)現(xiàn)并投入開發(fā),如何實(shí)現(xiàn)此類氣藏的經(jīng)濟(jì)有效開發(fā)已引起油氣生產(chǎn)與科研人員的高度重視并取得了可喜的成效[11-15]。

    1 開發(fā)現(xiàn)狀及對(duì)地震技術(shù)的需求

    蘇里格氣田是我國目前已發(fā)現(xiàn)的氣田中最大的天然氣氣田,主力產(chǎn)氣層為二疊系下石盒子組盒8段,儲(chǔ)層巖性主要為三角洲平原分流河道砂泥巖薄互層沉積的中-粗粒石英砂巖、砂巖,儲(chǔ)層的分布受巖性和物性的雙重控制:①產(chǎn)層段砂體縱向多期疊置、橫向復(fù)合連片,累計(jì)砂體厚度大,但單砂體規(guī)模小且多呈透鏡狀,厚度小于5 m的單砂層占78%以上;橫向上物性、含氣性變化快,儲(chǔ)層平均孔隙度為2%~12%,平均滲透率為0.06 mD~2 mD,有效砂體的橫向分布局限、連通性差,有效砂層厚為3 m~10 m,厚度小于5 m的有效單砂層更是占88%以上;②氣藏大面積含氣,但儲(chǔ)量豐度低,鉆井結(jié)果往往表現(xiàn)為井井有氣而單井產(chǎn)量低的特點(diǎn)[11,13-14]。因此,如何提高單井產(chǎn)量與氣田的采收率、實(shí)現(xiàn)經(jīng)濟(jì)有效開發(fā)是蘇里格氣田面臨的重大難題。

    蘇里格氣田自2001年發(fā)現(xiàn)并建成投產(chǎn)以來,截止2012年底已建成了年產(chǎn)300×108m3以上的天然氣生產(chǎn)能力,當(dāng)年產(chǎn)氣量達(dá)290×108m3,已經(jīng)成為了我國重要的天然氣生產(chǎn)基地[15]。地震技術(shù)在蘇里格氣田的前期開發(fā)中,實(shí)現(xiàn)了儲(chǔ)層的定性識(shí)別,并通過疊前疊后儲(chǔ)層綜合預(yù)測方法有效預(yù)測了砂巖儲(chǔ)層含氣富集區(qū),為蘇里格氣田的開發(fā)做出了重要的貢獻(xiàn)。盡管目前的開發(fā)井網(wǎng)已達(dá)到一定密度,但由于主力產(chǎn)層具有極強(qiáng)的非均質(zhì)性,僅僅依靠完鉆井資料的對(duì)比就對(duì)井間的地質(zhì)及其含氣特征進(jìn)行推斷是遠(yuǎn)遠(yuǎn)不夠的,其結(jié)果也是非常不確定的;而地震資料具有橫向采樣率高、可追蹤的優(yōu)勢,地震技術(shù)是目前唯一可對(duì)井間信息進(jìn)行直接量化描述的技術(shù),因此,必須充分利用地震信息與技術(shù)才能真實(shí)再現(xiàn)地下目標(biāo)地質(zhì)體的三維空間展布。近年來,蘇里格氣田的開發(fā)方式從初期以直井開發(fā)為主逐漸轉(zhuǎn)變?yōu)閰彩骄?、水平井開發(fā)并重,而隨著大斜度水平井的規(guī)模應(yīng)用,已有的地震預(yù)測方法與技術(shù)已越來越難以滿足水平井開發(fā)的需求,存在的主要問題有:①受地震帶寬的影響,地震預(yù)測結(jié)果縱向分辨率低,難以精細(xì)描述不同巖性體、儲(chǔ)層的厚度、孔隙度、滲透率、含氣飽和度等氣藏參數(shù)的空間展布特征,也無法滿足水平井開發(fā)對(duì)于單砂層、特別是有效單砂層的識(shí)別要求,準(zhǔn)確確定水平實(shí)施段的目標(biāo)層位置難度極大;②地震無法提供水平井開發(fā)所需的高精度構(gòu)造成果與深度域儲(chǔ)層成果,設(shè)計(jì)目標(biāo)層深度與實(shí)鉆深度往往偏差很大,也難以準(zhǔn)確把握入靶的深度與水平段軌跡,常常造成提前或推后于實(shí)際目標(biāo)層位置實(shí)施水平段而未達(dá)到預(yù)期效果。因此,為滿足大斜度水平開發(fā)井規(guī)模應(yīng)用的需求,必須建立精確的速度場,準(zhǔn)確刻畫目的層砂體的深度,進(jìn)一步提高儲(chǔ)層預(yù)測精度,精細(xì)識(shí)別單砂體、特別是有效單砂層的空間展布特征,為水平井的設(shè)計(jì)與實(shí)施提供準(zhǔn)確的目標(biāo)層、精確的入靶深度與有效的水平段軌跡。

    通過持續(xù)的攻關(guān)研究,逐步形成了一套以地震屬性分析與巖石物理分析為基礎(chǔ),以高分辨率疊前反演技術(shù)與高精度時(shí)深轉(zhuǎn)換技術(shù)為核心的薄互儲(chǔ)層開發(fā)地震預(yù)測技術(shù)系列,應(yīng)用效果明顯,為蘇里格氣田的經(jīng)濟(jì)有效開發(fā)提供了重要的技術(shù)保障。

    2 開發(fā)地震技術(shù)在蘇里格氣田的應(yīng)用

    2.1基于地震屬性分析的主砂帶預(yù)測技術(shù)

    地震屬性是蘊(yùn)含在地震反射波中的有關(guān)地震波的幾何形態(tài)、運(yùn)動(dòng)學(xué)特征、動(dòng)力學(xué)特征和統(tǒng)計(jì)學(xué)特征的信息,它能從多方面反映地下特殊地質(zhì)現(xiàn)象、沉積現(xiàn)象、流體特征等。因此,地震屬性分析已成為進(jìn)行儲(chǔ)層描述的常用手段之一,也是勘探地震與開發(fā)地震之間的橋梁。對(duì)于蘇里格地區(qū)石盒子組盒8段這樣的非均質(zhì)性強(qiáng)的薄儲(chǔ)層,雖然常規(guī)的地震分辯率很難達(dá)到對(duì)單一河道砂體的準(zhǔn)確識(shí)別,但地震屬性分析可以客觀地反映儲(chǔ)層巖性、物性的差異性特征,宏觀反映出河道發(fā)育帶的展布情況。

    前期的勘探開發(fā)揭示,蘇里格氣田盒8段沉積時(shí)期為辮狀河三角洲平原沉積環(huán)境,發(fā)育辮狀河道、心灘、天然堤、決口扇、泛濫平原等微相類型。一般認(rèn)為盒8段的辮狀河道控制了河道砂體的規(guī)模,心灘控制了粗粒砂巖的分布,單一河道或心灘砂巖厚度往往較薄,而多期河道或心灘砂巖的疊置可以形成厚度較大的連續(xù)粗巖相沉積,從而發(fā)育良好的儲(chǔ)層,因此通過尋找主砂帶可以預(yù)測有效儲(chǔ)層發(fā)育的有利區(qū)。

    大量井震標(biāo)定結(jié)果表明,Tp8h的波峰反射是盒8段砂泥巖的綜合地震響應(yīng),而Tp8h反射的振幅能量橫向變化與巖性、巖相密切相關(guān),中強(qiáng)振幅區(qū)主要反映盒8段辮狀河河道或心灘厚層砂巖沉積區(qū),而弱振幅或空白反射區(qū)一般為河道間泛濫平原的細(xì)粒沉積(圖1)。盒8沉積期,辮狀河道雖然頻繁改道,但主河道仍分布在一定范圍內(nèi),因此通過分析Tp8h反射振幅的強(qiáng)弱,可以預(yù)測盒8段主河道的分布范圍。通過多種振幅屬性的比較分析,認(rèn)為均方根振幅屬性能較好地反映盒8段河道砂體的空間展布特征。

    圖1 蘇里格氣田A井區(qū)石盒子組盒8段地震相-沉積相解釋剖面Fig.1 The seismic facies and sedimentary facies interpretation section of the 8th member of Shihezi formation in A well area of Sulige gas field

    圖2是蘇里格氣田A井區(qū)盒8段的均方根振幅平面圖,圖2中的中強(qiáng)振幅區(qū)(黃色和紅色)呈近西北、南東向的條帶狀展布,辮狀河主河道分布形態(tài)清晰可見;中強(qiáng)振幅區(qū)內(nèi)的完鉆井在盒8段鉆遇厚砂巖的比例達(dá)82%。

    2.2巖石物理分析技術(shù)

    2.2.1橫波速度曲線預(yù)測

    橫波速度是開展巖石物理分析、疊前彈性參數(shù)反演的關(guān)鍵巖石物理參數(shù),對(duì)于儲(chǔ)層巖性的識(shí)別、儲(chǔ)層物性及含油氣性的預(yù)測具有至關(guān)重要的作用,然而目前蘇里格氣田的橫波測井曲線資料卻十分匱乏。因此,準(zhǔn)確預(yù)測盒8段的橫波速度曲線,是蘇里格氣田開發(fā)地震所迫切需要解決的問題。目前,橫波速度曲線的預(yù)測方法主要有經(jīng)驗(yàn)公式法和巖石物理模型法兩大類。經(jīng)驗(yàn)公式方法以縱波測井?dāng)?shù)據(jù)為基礎(chǔ),通過擬合縱波速度、巖石礦物成分、密度、泥質(zhì)含量等參數(shù)與橫波速度的關(guān)系來求取橫波速度,其中Greenberg-Castagna公式法是常用的經(jīng)驗(yàn)公式法[16-18];巖石物理理論模型法是通過建立橫波速度與巖石物理參數(shù)之間的理論關(guān)系,由一條或多條測井曲線預(yù)測橫波速度曲線的方法(Gassmann方程法、Xu-White模型法)。利用經(jīng)驗(yàn)公式法預(yù)測橫波速度較為方便、快捷,但由于各地區(qū)的經(jīng)驗(yàn)參數(shù)差異較大而難以準(zhǔn)確求取,得到的橫波速度往往與實(shí)測橫波速度存在較大的誤差。Xu-White模型法是建立在Gassmann方程基礎(chǔ)上的一種理論巖石物理模型,理論計(jì)算公式較為繁雜,在此不做闡述,可參閱文獻(xiàn)[19-20]。利用Xu-White模型法計(jì)算橫波速度的理論公式雖然相對(duì)繁雜,輸入?yún)?shù)也較多,但該方法是一種完全的理論模型,預(yù)測結(jié)果能夠同時(shí)獲得橫波速度和縱波速度,因此可以通過不斷調(diào)整模型參數(shù)來減小預(yù)測縱波速度與實(shí)測縱波速度的誤差,從而達(dá)到約束、提高橫波速度的預(yù)測精度的目的,故其計(jì)算精度較高,是目前較為常用的一種橫波速度預(yù)測方法。

    圖2 A井區(qū)盒8段均方根振幅平面圖Fig.2 The RMS amplitude plan of the 8th member of Shihezi formation in A well area

    X井是蘇里格氣田A井區(qū)內(nèi)唯一具有實(shí)測橫波測井曲線的井,圖3是X井利用Greenberg - Castagna公式法與Xu-White模型法預(yù)測的橫波速度結(jié)果。由圖3可以看出,利用Greenberg-Castagna公式法預(yù)測的橫波速度曲線(Vs_GC)與實(shí)測的橫波速度曲線(DTS)存在較大的誤差,而基于Xu-White模型法預(yù)測的橫波速度曲線(Vs_Xu)與實(shí)測的橫波速度曲線具有很好的一致性,吻合程度較高,誤差較小。實(shí)驗(yàn)表明,Xu-White模型法是預(yù)測該區(qū)砂泥巖地層橫波速度曲線的有效技術(shù)方法。在此基礎(chǔ)上,根據(jù)X井測試獲得的各項(xiàng)巖石物理模型輸入?yún)?shù),利用Xu-White模型法預(yù)測了其余井的橫波速度曲線,預(yù)測精度均較高,為后續(xù)的巖石物理分析、疊前反演提供了可靠的橫波速度數(shù)據(jù)。

    第一道GR是自然伽馬曲線;第二道DEN是密度曲線;第三道AC是聲波時(shí)差曲線轉(zhuǎn)換而來的縱波速度曲線;第四道的DTS為實(shí)測橫波速度曲線,Vs_GC為利用Greenberg-Castagna公式法預(yù)測的橫波速度曲線,Vs_Xu是利用Xu-White模型法預(yù)測的橫波速度曲線;第五道是砂泥巖剖面,CARB為泥質(zhì)含量,SAND為砂巖含量,POR為孔隙度(巖性剖面中的空白部分)圖3 X井橫波速度曲線預(yù)測圖Fig.3 The shear-wave curve prediction map of X well

    2.2.2巖石物理分析

    巖石物理分析是建立油氣藏特征參數(shù)與地震(反演)數(shù)據(jù)之間關(guān)系的橋梁,是利用地震反演結(jié)果開展巖性識(shí)別、儲(chǔ)層精細(xì)描述與流體預(yù)測的前提和基礎(chǔ),是地震勘探識(shí)別精度從勘探階段的定性向開發(fā)階段的半定量到定量轉(zhuǎn)變的重要手段。通過盒8段的測井巖石彈性參數(shù)交會(huì)分析可以看出(圖4),無論是縱波阻抗還是橫波阻抗,盒8段的砂巖與泥巖間存在很大的重疊空間,難以直接將兩者很好地分開;而泊松比與縱波阻抗、橫波阻抗的交匯,可以很好地識(shí)別砂巖、含氣砂巖與泥巖,通常砂巖的泊松比小于0.18;泊松比與縱波阻抗的交會(huì)識(shí)別含氣砂巖效果更好。因此,采用常規(guī)的疊后縱波阻抗反演技術(shù)難以準(zhǔn)確識(shí)別砂巖與泥巖,只有通過疊前反演技術(shù),求取縱波阻抗、橫波阻抗,進(jìn)而獲取縱橫波速度比、泊松比等巖石彈性參數(shù),才能將砂巖、含氣砂巖與泥巖進(jìn)行有效地區(qū)分。

    孔隙度、滲透率是蘇里格氣田開發(fā)中儲(chǔ)層預(yù)測與描述的重要物性參數(shù),利用地震資料往往無法直接獲取這些物性參數(shù),但由于地震資料具有覆蓋面廣、橫向分辨率高的優(yōu)勢,利用基于地震資料反演結(jié)果獲得的地震巖石物理參數(shù)與測井解釋孔隙度、滲透率之間的數(shù)學(xué)統(tǒng)計(jì)關(guān)系,進(jìn)而估算這些參數(shù)值的方法是目前應(yīng)用較為普遍、效果也較好的方法。盒8段的多種巖石物理參數(shù)交會(huì)分析表明,有效砂巖的縱波阻抗與孔隙度、孔隙度與滲透率之間存在良好的線性關(guān)系(圖5),其數(shù)學(xué)表達(dá)式如式(1)、式(2)所示。

    φ=86.1701-0.0144862×ZP+

    (1)

    K=0.00120045+0.0453542×φ-

    0.0104677×φ2+0.00125738×φ3

    (2)

    式中:φ為孔隙度,%;ZP為縱波阻抗,g·cm-3·m·s-1;K為滲透率,mD。

    2.3高分辨率疊前反演技術(shù)

    疊后波阻抗反演、疊前同時(shí)反演或彈性阻抗反演等確定性反演技術(shù),是蘇里格氣田目前常用的儲(chǔ)層地震反演技術(shù),在蘇里格氣田的勘探階段以及開發(fā)前期發(fā)揮了重要的作用。但由于蘇里格氣田盒8段砂巖與泥巖互層沉積,橫向相變快,單砂體、有效砂體厚度薄,儲(chǔ)層孔隙度、滲透率低,屬典型的低孔低滲型儲(chǔ)層,且砂巖的波阻抗特征與泥巖疊置嚴(yán)重,利用常規(guī)的確定性反演開展優(yōu)質(zhì)儲(chǔ)層識(shí)別難度較大、預(yù)測精度較低,尤其是預(yù)測結(jié)果受到地震帶寬的限制,縱向分辨率低,無法滿足氣田開發(fā)對(duì)識(shí)別單砂層特別是有效單砂層的需求,而準(zhǔn)確識(shí)別盒8段單砂層的規(guī)模及其縱橫向展布并優(yōu)選出高孔砂層,是開發(fā)中進(jìn)行定向井和大斜度水平井設(shè)計(jì)進(jìn)而提高單井產(chǎn)能的前提和關(guān)鍵。近年來,基于地質(zhì)統(tǒng)計(jì)學(xué)思想的地震隨機(jī)反演受到了國內(nèi)、外越來越多油氣工作者的重視[21-23],將疊前確定性反演技術(shù)與隨機(jī)反演技術(shù)相結(jié)合的高分辨率疊前反演技術(shù)更是倍受關(guān)注[24-26]。

    圖4 蘇里格氣田A井區(qū)盒8段巖石物理參數(shù)交會(huì)圖Fig.4 The petrophysical parameter cross-plot of the 8th member of Shihezi formation in A well area of Sulige gas field(a)泊松比與橫波阻抗的交會(huì)圖;(b)泊松比與縱波阻抗的交會(huì)圖;(c)橫波阻抗與縱波阻抗的交會(huì)圖

    圖5 盒8段有效砂巖的物性參數(shù)交會(huì)圖Fig.5 The physical parameters cross-plot of the effective sandstone in the 8thmember of Shihezi formation(a)縱波阻抗與孔隙度交會(huì)圖; (b)孔隙度與滲透率交會(huì)圖

    高分辨率疊前反演技術(shù)是以測井、地震、地質(zhì)信息為基礎(chǔ),采用嚴(yán)格的馬爾科夫鏈-蒙特卡羅算法(Markov Chain- Monte Carlo,MCMC),將疊前確定性反演與地質(zhì)統(tǒng)計(jì)學(xué)隨機(jī)模擬相結(jié)合的同時(shí)AVA隨機(jī)反演方法[26]。在高分辨率疊前反演中,以地質(zhì)統(tǒng)計(jì)學(xué)分析結(jié)果為基礎(chǔ),以疊前地震數(shù)據(jù)為約束,包括隨機(jī)模擬過程和反演過程。模擬和反演均以井點(diǎn)出發(fā),通過MCMC采樣方法對(duì)概率密度函數(shù)進(jìn)行充分采樣,從而對(duì)確定性反演中喪失的高頻數(shù)據(jù)信息進(jìn)行了補(bǔ)償,而井間嚴(yán)格遵從原始地震數(shù)據(jù),每次模擬所對(duì)應(yīng)的合成地震記錄必須與實(shí)際的地震數(shù)據(jù)具有很高的相似性,因此反演的結(jié)果既獲得了超過地震帶寬的、與測井?dāng)?shù)據(jù)接近的縱向分辨率,同時(shí)也保證了與地震數(shù)據(jù)完全相同的空間分布趨勢。

    圖6為蘇里格氣田A井區(qū)過X1、X2、X3井的高分辨率疊前反演結(jié)果。由圖6可以看出:①泊松比反演結(jié)果具有很高的縱向分辨率,縱向上預(yù)測結(jié)果的各砂層(泊松比小于0.18)位置、厚度與伽馬曲線解釋結(jié)果吻合較好,X1井砂薄(13.6 m)、X2井砂厚(29.9 m)、X3井較厚(27.8 m)的橫向展布特征刻畫清晰,河道、心灘砂體表現(xiàn)出明顯的縱向疊置、橫向連片的展布特征,盲井(X3井)的反演結(jié)果與實(shí)鉆結(jié)果的一致性高;②預(yù)測的盒8段砂體孔隙度、滲透率與測井解釋結(jié)果吻合良好,孔隙度、滲透率的橫向分布變化較大、非均質(zhì)性極強(qiáng),儲(chǔ)層孔隙度主要分布在4%~12%,滲透率主要分布在0.2 mD~1.2 mD。利用預(yù)測的孔隙度、滲透率結(jié)果,可以清晰刻畫出高孔隙度、高滲透率的有效單砂層的空間展布,為后續(xù)的儲(chǔ)層精細(xì)描述、儲(chǔ)量計(jì)算,以及為定向井和大斜度水平井的鉆井軌跡和靶點(diǎn)的優(yōu)選提供了可靠的依據(jù)。

    圖6 蘇里格氣田A井區(qū)盒8段高分辨率疊前反演結(jié)果Fig.6 The high-resolution pre-stack inversion result of the 8th member of Shihezi formation in A well area, Sulige gas field(a)泊松比剖面;(b)孔隙度剖面;(c)滲透率剖面

    2.4高精度時(shí)深轉(zhuǎn)換技術(shù)

    在勘探階段,基于大套層段的層控變層速度時(shí)深轉(zhuǎn)換技術(shù),能夠很好滿足目的層系頂、底界面構(gòu)造成果的精度需求[27-28],但已遠(yuǎn)遠(yuǎn)無法滿足處于開發(fā)階段的蘇里格低滲透氣田需要精確定位大套地質(zhì)層段內(nèi)部有效單砂層位置及其深度的精度要求?,F(xiàn)場技術(shù)人員常常通過相鄰已完鉆井的目標(biāo)層海拔、地層厚度、巖性等進(jìn)行對(duì)比來確定擬鉆井的目標(biāo)層入靶點(diǎn)位置、深度[29],由于盒8砂層尤其是有效砂層在橫向上變化較快,導(dǎo)致預(yù)測精度往往不高。在蘇里格氣田開發(fā)地震解釋中,提出一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的高精度時(shí)深轉(zhuǎn)換方法,其關(guān)鍵是建立精確的三維空間時(shí)深轉(zhuǎn)換速度場,主要技術(shù)思路與流程為:①通過對(duì)已完鉆井VSP、聲波測井的縱向速度結(jié)構(gòu)分析與精細(xì)的井震標(biāo)定,確定合適數(shù)量的區(qū)域性速度控制層的時(shí)間域地震反射界面并對(duì)其進(jìn)行精細(xì)的對(duì)比解釋,在此基礎(chǔ)上建立由區(qū)域性速度控制層反射界面控制的全區(qū)三維地震層位框架模型;②利用已完鉆井的VSP、聲波測井等速度數(shù)據(jù)建立各速度控制層段的速度概率密度函數(shù)、速度縱向變差函數(shù),根據(jù)各速度控制層段的地震屬性、沉積相、巖相等研究成果建立各速度控制層段的橫向變差函數(shù);③在全區(qū)三維地震層位框架模型的控制下,利用各速度控制層段的速度概率密度函數(shù)、速度縱向變差函數(shù)和速度橫向變差函數(shù)采用序貫高斯協(xié)模擬算法進(jìn)行速度地質(zhì)統(tǒng)計(jì)學(xué)模擬,得到多個(gè)速度模型實(shí)現(xiàn),并對(duì)多個(gè)速度模型實(shí)現(xiàn)進(jìn)行統(tǒng)計(jì)計(jì)算,建立全區(qū)三維時(shí)深轉(zhuǎn)換速度場,該方法仍然采用了大套層段進(jìn)行層控的變層速度建模這一基本思路,但井點(diǎn)位置的大套層段內(nèi)部采用VSP、聲波測井等實(shí)測速度數(shù)據(jù),井點(diǎn)與井點(diǎn)之間的速度通過地質(zhì)統(tǒng)計(jì)學(xué)模擬得到,這樣既充分考慮了巖相、沉積相等地質(zhì)因素對(duì)速度的控制作用,其速度求取算法也明顯優(yōu)于常規(guī)的線性插值方法,結(jié)果更符合地質(zhì)規(guī)律、精度更高;④在此基礎(chǔ)上,利用已建立的高精度三維空間時(shí)深轉(zhuǎn)換速度場,將時(shí)間域的地震數(shù)據(jù)(疊加數(shù)據(jù)、偏移數(shù)據(jù))、高分辨率疊前反演成果數(shù)據(jù)(縱波阻抗數(shù)據(jù)、橫波阻抗數(shù)據(jù)、縱橫波速度比數(shù)據(jù)、泊松比數(shù)據(jù)、孔隙度數(shù)據(jù)、滲透率數(shù)據(jù)等)轉(zhuǎn)換為深度域,該方法在蘇里格氣田A井區(qū)的應(yīng)用表明,利用轉(zhuǎn)換得到的深度域數(shù)據(jù),不僅可以準(zhǔn)確落實(shí)各大套層段頂、底界面的構(gòu)造起伏形態(tài),而且能精細(xì)刻畫目的層段盒8內(nèi)部各砂層的深度位置及其空間展布,從而為定向井和大斜度水平井準(zhǔn)確定位目標(biāo)砂層的位置、入靶深度,并確定水平段穿越的砂層個(gè)數(shù)、水平段長度及鉆井軌跡(圖7),極大地提高了有效儲(chǔ)層鉆遇率和井口產(chǎn)能,鉆探效果達(dá)到了預(yù)期目的。

    圖7 A井區(qū)盒8段的深度域泊松比剖面與孔隙度剖面Fig.7 The poisson's ratio section and porosity section in depth field of the 8thmember of Shihezi formation in A well area(a)泊松比剖面;(b) 孔隙度剖面

    2.5三維可視化技術(shù)

    三維可視化技術(shù)是利用三維地震數(shù)據(jù)體顯示、描述和理解地下諸多地質(zhì)現(xiàn)象和特征的一種圖像顯示工具,已被廣泛應(yīng)用于油氣勘探與開發(fā)的各個(gè)領(lǐng)域,成為地震資料解釋及儲(chǔ)層預(yù)測與精細(xì)描述的一項(xiàng)關(guān)鍵技術(shù)。它利用圖形工作站的先進(jìn)顯示技術(shù),通過改變可視化參數(shù)對(duì)不同的地震數(shù)據(jù)體進(jìn)行全方位的透視,使三維數(shù)據(jù)體內(nèi)的地質(zhì)特征信息立體地顯示出來。

    蘇里格氣田盒8段儲(chǔ)層雖然非均質(zhì)性極強(qiáng),但通過利用三維可視化技術(shù)精細(xì)刻畫了砂體、有效砂體的空間展布特征,為開展定向井、大斜度水平井的井軌跡設(shè)計(jì)提供了快速的三維顯示,大大節(jié)省了地震解釋與開發(fā)人員的決策時(shí)間,提高了鉆井成功率。①通過對(duì)常規(guī)的振幅數(shù)據(jù)體以及經(jīng)過提取處理的各種屬性(如均方根振幅屬性、反射強(qiáng)度屬性、分頻屬性)等數(shù)據(jù)體進(jìn)行三維可視化,快速識(shí)別有利沉積相帶、厚砂層分布區(qū);②利用反演得到的縱波阻抗、橫波阻抗、縱橫波速度比、泊松比、孔隙度以及滲透率等數(shù)據(jù)體進(jìn)行三維可視化透視顯示,精細(xì)雕刻了蘇里格氣田A井區(qū)盒8段砂層、有效砂層的空間展布特征(圖8),為開發(fā)井位的部署、水平井軌跡的設(shè)計(jì)提供了直觀、可靠的依據(jù)。

    3 結(jié)論

    1)均方根振幅屬性能宏觀反映盒8段主砂帶的空間展布特征,中強(qiáng)振幅主要反映河道或心灘厚層砂巖沉積,弱振幅或空白反射一般為泥巖或薄砂巖沉積。

    2)Xu-White模型法是該區(qū)盒8砂泥巖段預(yù)測橫波速度的有效技術(shù)方法,泊松比是準(zhǔn)確識(shí)別盒8段砂巖、含氣砂巖與泥巖的敏感參數(shù),砂巖的縱波阻抗與孔隙度、孔隙度與滲透率之間存在良好的線性關(guān)系。

    3)高分辨率疊前反演技術(shù)極大地提高了預(yù)測結(jié)果的縱向分辨率,有效解決了砂泥巖薄互層的單砂體地震識(shí)別難題,是陸相致密碎屑巖性油氣藏進(jìn)行儲(chǔ)層精細(xì)描述與高效開發(fā)的有效技術(shù)方法。

    4)高精度時(shí)深轉(zhuǎn)換技術(shù)與三維可視化技術(shù)相結(jié)合,準(zhǔn)確落實(shí)并精細(xì)雕刻了各砂層、有效砂層的深度位置及其空間展布特征,為水平開發(fā)井位的部署、井軌跡的設(shè)計(jì)提供了直觀、可靠的依據(jù)。

    圖8 盒8段高孔砂體的三維空間展布Fig.8 The 3D distribution of high- porosity sandstone in the 8th member of Shihezi formation(a)整個(gè)三維空間; (b)兩條不同方向剖面

    [1]劉雯林.油氣田開發(fā)地震技術(shù)[M].北京:石油工業(yè)出版社,1996.

    LIU W l.Development seismic technology in oil & gas field[M].Beijing:Oil Industry Publishing Company,1996.(In Chinese)

    [2]劉振武,撒利明,張昕,等.中國石油開發(fā)地震技術(shù)應(yīng)用現(xiàn)狀和未來發(fā)展建議[J].石油學(xué)報(bào),2009,30(5):711-716.

    LIU Z W,SA L M,ZHANG X,et al.Current application and future development of production seismology in PetroChina[J].Acta Petrolei Sinica,2009,30(5):711-716.(In Chinese)

    [3]撒利明,甘利燈,黃旭日,等.中國石油集團(tuán)油藏地球物理技術(shù)現(xiàn)狀與發(fā)展方向[J].石油地球物理勘探,2014,49(3):611-626.

    SA L M,GAN L D,HUANG X R,et al.The Status and developmnet direction of reservoir geophysics in CNPC[J].Oil Geophysical Prospecting,2014,49(3):611-626.(In Chinese)

    [4]馬在田.關(guān)于油氣開發(fā)地震學(xué)的思考[J].天然氣工業(yè),2004,24(6):43-46.

    MA Z T.Some thoughts on petroleum development seismology[J].Natural Gas Intustry,2004,24(6):43-46.(In Chinese)

    [5]康德權(quán).開發(fā)地震技術(shù)在大慶外圍東部地區(qū)勘探開發(fā)中的應(yīng)用[J].資源與產(chǎn)業(yè),2009,11(2):140-145.

    KANG D Q.Application of development seismic technology in oil exploration in eastern Daqing oilfiled[J].Resources & Industries,2009,11(2):140-145.(In Chinese)

    [6]李棟明,劉群星.地震技術(shù)在油氣田開發(fā)中的應(yīng)用[J].內(nèi)蒙古石油化工,2008(4):97-99.

    LI D M,LIU Q X.The application of seismic technology in oil &gas field development[J].Inner Mongolia petrochemical industry,2008(4):97-99.(In Chinese)

    [7]劉薇薇,周家雄,馬光克,等.東方1-1氣田開發(fā)地震技術(shù)的應(yīng)用[J].天然氣工業(yè),2012,32(8):22-26.

    LIU W W,ZHOU J X,MA G K,et al.Application of seismic techniques to the development of the Dongfang 1-1 gas field,Yingahai Basin[J].Natural Gas Intustry,2012,32(8):22-26.(In Chinese)

    [8]魏小東,涂小仙,曹麗麗,等.高精度三維地震資料精細(xì)解釋技術(shù)在吐哈丘陵油田開發(fā)中的應(yīng)用[J].石油地球物理勘探,2008,43(增刊1):92-97.

    WEI X D,TU X X,CAO L L,et al.Application of fine interpretation technique of 3-D high precision seismic data to development of Tu-ha hill oilfield[J].Oil Geophysical Prospecting,2008,43(Supplement 1):92-97.(In Chinese)

    [9]羅士利,高興友,彭承文.開發(fā)地震技術(shù)在扶余油層分支水平井地質(zhì)設(shè)計(jì)中的應(yīng)用[J].吉林大學(xué)學(xué)報(bào):自然科學(xué)版,2006,37(專輯):109-116.

    LUO S L,GAO X Y,PENG C W.The application of production seismology techniques in branched horizontal well design of Fuyu oil layer[J].Journal of Jilin University:Earth Science Edition,2006,37(Supplement):109-116.(In Chinese)

    [10]任芳祥,孫巖,樸永紅,等.地震技術(shù)在油田開發(fā)中的應(yīng)用問題探討.地球物理學(xué)進(jìn)展,2010,25(1):4-8.

    RRN F X,SUN Y,PIAO Y H,et al.Discussion on application of seismic technology to oilfield development[J].Progress in Geophysics,2010,25(1):4-8.(In Chinese)

    [11]冉新權(quán),李安琪.蘇里格氣田開發(fā)論[M].北京:石油工業(yè)出版社,2008.

    RAN X Q,LI A Q.Sulige gas field development approach[M].Beijing:Oil Industry Publishing Company,2008.(In Chinese)

    [12]劉振武,撒利明,張研,等.中國天然氣勘探開發(fā)現(xiàn)狀及物探技術(shù)需求[J].天然氣工業(yè),2009,29(1):1-7.

    LIU Z W,SA L M,ZHANG Y,et al.Status quo of natural gas exploration and development and technical demand for geophysical exploration in China[J].Natural Gas Intustry,2009,29(1):1-7.(In Chinese)

    [13]張明祿,吳正,樊友宏,等.鄂爾多斯盆地低滲透氣藏開發(fā)技術(shù)及開發(fā)前景[J].天然氣工業(yè),2011,31(7):1-4.

    ZHANG M L,WU Z,FAN Y H,et al.Technologies and prospects of the development of low permeability gas reservoirs in the Ordos Basin[J].Natural Gas Intustry,2011,31(7):1-4.(In Chinese)

    [14]楊華,劉新社,楊勇.鄂爾多斯盆地致密氣勘探開發(fā)形勢與未來發(fā)展展望[J].中國工程科學(xué),2012,14(6):40-48.

    YANG H,LIU X S,YANG Y.Status and prospects of tight gas exploration and develoment in the Ordos basin[J].Engineering Sciences,2012,14(6):40-48.(In Chinese)

    [15]張明祿,樊友宏,何光懷,等.長慶氣區(qū)低滲透氣藏開發(fā)技術(shù)新進(jìn)展及攻關(guān)方向[J].天然氣工業(yè),2013,33(8):1-7.

    ZHANG M L,FAN Y H,HE G H,et al.Latest progress in development technologies for low-permeability gas reservoirs in the Changqing gas zone[J].Natural Gas Intustry,2013,33(8):1-7.(In Chinese)

    [16]李慶忠.巖石的縱、橫波速度規(guī)律[J].石油地球物理勘探,1992,27(1):1-12.

    LI Q Z.Velocity regularities of P and S-waves in formations[J].Oil Geophysical Prospecting,1992,27(1):1-12.(In Chinese)

    [17]CASTAGNA J P,BATZLE M L,EASTWOOD R I.relationship between compressional wave and shear wave velocities in clastic rocks[J].Geophysics,1985,50(5):571-581.

    [18]GREENBERG M L,CASTAGNA J P.Shearwave velocity estimation in porous rocks:theoretical formulation,preliminary verification and application[J].Geophysical Prospecting,1992,40(2):195-209.

    [19]XU S,WHITE R E.A new velocity model for clay-sand mixtures[J].Geophysical Prospecting,1995,43(1):91-118.

    [20]XU S,WHITE R E.A physical model for shear-wave velocity prediction[J].Geophysical Prospecting,1996,44(4):687-717.

    [21]孫思敏,彭仕宓.地質(zhì)統(tǒng)計(jì)學(xué)反演及其在吉林扶余油田儲(chǔ)層預(yù)測中的應(yīng)用[J].物探與化探,2007,31(1):51-54.

    SUN S M,PENG S M.Geostatistical inversion method and its application to reservoir prediction of Fuyu oilfield in Jilin province[J].Geophysical and Geochemical Exploration,2007,31(1):51-54.(In Chinese)

    [22]韓翀,臧殿光,李建華.地質(zhì)統(tǒng)計(jì)學(xué)反演在四川盆地L地區(qū)儲(chǔ)層預(yù)測中的應(yīng)用[J].物探化探計(jì)算技術(shù),2010,32(3):310-315.

    HAN C,ZANG D G,LI J H.The application of geostatistical inversion to reservoir prediction of L area in Sichuan basin[J].Computing Techniques for Geophysical and Geochemical Exploration,2010,32(3):310-315.(In Chinese)

    [23]孫煒,王彥春,李玉鳳,等.基于地質(zhì)統(tǒng)計(jì)學(xué)反演的碳酸鹽巖孔縫洞預(yù)測研究[J].現(xiàn)代地質(zhì),2012,26(6):1258-1263.

    SUN W,WANG Y C,LI Y F,et al.Fracture and solution pore prediction in carbonate reservoirs based on geostatistics inversion[J].Geoscience,2012,26(6):1258-1263.(In Chinese)

    [24]姜文龍,楊鍇.巖石物理參數(shù)高分辨率地質(zhì)統(tǒng)計(jì)學(xué)反演[J].石油物探,2012,51(6):639-647.

    JIANG W L ,YANG K.High-resolution geostatistical petrophysical-parameter inversion[J].Geopthysical Prospecting for Petroleum,2012,51(6):639-647.(In Chinese)

    [25]姜巖,徐立恒,張秀麗,等.疊前地質(zhì)統(tǒng)計(jì)學(xué)反演方法在長垣油田儲(chǔ)層預(yù)測中的應(yīng)用[J].地球物理學(xué)進(jìn)展,2013,28(5):2579-2585.

    JIANG Y,XU L H,ZHANG X L,et al.Prestack geostatistical inversion method and its application on the reservoir prediction of Changyuan oil field[J].Progress In Geophysics,2013,28(5):2579-2585.(In Chinese)

    [26]錢玉貴,葉泰然,張世華,等.疊前地質(zhì)統(tǒng)計(jì)學(xué)反演技術(shù)在復(fù)雜儲(chǔ)層量化預(yù)測中的應(yīng)用[J].石油與天然氣地質(zhì),2013,34(6):834-839.

    QIAN Y G,YE T R,ZHANG S H,et al.Application of pre-stack geostatistics inversion technology in quantitative prediction of complex reservoirs[J].Oil & Gas Geology,2013,34(6):834-839.(In Chinese)

    [27]張華軍,王海蘭,肖富森,等.基于反射層的變層速度模型時(shí)深轉(zhuǎn)換方法[J].天然氣工業(yè),2003,23(1):36-38.

    ZHANG H J,WANG H L,XIAO F S,et al.Time-depth conversion method in view of variable interval velocity model of reflection horizon[J].Natural Gas Intustry,2003,23(1):36-38.(In Chinese)

    [28]張健,汪少華,陳家琪,等.變速成圖技術(shù)在準(zhǔn)噶爾盆地南緣沖斷帶的應(yīng)用[J].天然氣工業(yè),2010,30(4):35-37.

    ZHANG J,WANG S H,CHEN J Q,et al.Application of varying velocity mapping in the thrust belt in the southern margin of the Junggar Basin[J].Natural Gas Intustry,2010,30(4):35-37.(In Chinese)

    [29]李建奇,楊志倫,陳啟文,等.蘇里格氣田水平井開發(fā)技術(shù)[J].天然氣工業(yè),2011,31(8):1-8.

    LI J Q,YANG Z L,CHEN Q W,et al.Horizontal well technology for the development of the Sulige gas Field[J].Natural Gas Intustry,2011,31(8):1-8.(In Chinese)

    The application of development seismic technology in Sulige gas field

    PANG Chong-you1, LI Da-jun1, ZHANG Xiong1, ZHANG Ya-dong2, XIA Ming1, WU Tie-ling2

    (1. Geophysical Prospecting Company, Chuanqing Drilling Engoneering Ltd Co., CNPC, Chengdu610213, China;2. Department of Exploration, Changqing Oilfield Company, PetroChina Company Limited, Xi’an710018, China)

    The present seismic prediction methods and technologies cannot produce the demand for the development of horizontal wells in the Member 8 gas reservoir of Shihezi formation in Sulige gas field with low porosity and low permeability. Contrapose to this situation, the development seismic prediction technology series based on seismic attributes analysis technique and petrophysical analysis technique, and that is centering on high-resolution pre-stack incersion technology and high-precision time-depth conversion method is provided. Firstly, the distribution of host sand band is macroscopically forecasted according to seismic attributes analysis and the shear-wave velocity curves are calculated, and then the petrophysical analysis is developed. Based on the above situations, the spatial distribution and depth of the effective sand layer are accurately confirmed after high-resolution pre-stack incursion and high-precision time-depth conversion. The application results show that the RMS amplitude attribute can macroscopically response the spatial distribution of host sand band in the member 8 of Shihezi formation. The Xu-White model method is a valuable method for predicting shear-wave velocity curve and poisson's ratio is the sensitive parameter for identifying sandstone, gas sandstone and mudstone accurately. Besides, the high-resolution pre-stack incersion technology greatly improves the vertical resolution of seismic prediction results, thus the problem for identifying single sand body in thin alternating beds of sandstone and mudstone by seismic method has been solved effectively in the 8th member of Shihezi formation. In addition, the high-precision time-depth conversion method that combined with 3D visualization techniques provides intuitive and reliable basis for deploying horizontal well and designing well trajectory. Therfore, this technology series is a valuable technology series for fine reservoir description and high efficiency development in the continental tight clastic lithological oil-gas reservoir.

    development seismic technology; Sulige gas field; host sand band; shear-wave curve prediction; petrophysical analysis

    2015-05-19改回日期:2015-06-23

    國家重大科技專項(xiàng)(2011ZX05013-001-02)

    龐崇友(1979-),男,碩士,工程師,研究方向?yàn)橛蜌獾卣鸬刭|(zhì)綜合,E-mail:pangchy2015@163.com。

    1001-1749(2016)04-0530-10

    P 631.4

    A

    10.3969/j.issn.1001-1749.2016.04.15

    猜你喜歡
    里格橫波砂層
    橫波技術(shù)在工程物探中的應(yīng)用分析
    第四系膠結(jié)砂層水文地質(zhì)鉆進(jìn)技術(shù)研究初探
    西安主城區(qū)砂層工程地質(zhì)特征分析
    蘇里格氣田總產(chǎn)量突破1000×108m3
    蘇里格氣田致密砂巖氣層識(shí)別難點(diǎn)及方法評(píng)述
    汽車制動(dòng)檢驗(yàn)臺(tái)滾筒粘砂層脫落的原因及維護(hù)方法
    揚(yáng)眉一顧,妖嬈橫波處
    橫波一顧,傲殺人間萬戶侯
    火花(2015年1期)2015-02-27 07:40:24
    蘇里格氣田東區(qū)組合式排水采氣應(yīng)用淺析
    橫波淺層地震在城市勘探中的應(yīng)用
    久久香蕉国产精品| 中文字幕av在线有码专区| 欧美最黄视频在线播放免费| 久久国产精品影院| 草草在线视频免费看| 国产精品日韩av在线免费观看| 亚洲精品粉嫩美女一区| 黄片小视频在线播放| 99热只有精品国产| 黄色女人牲交| 最新中文字幕久久久久 | 又黄又爽又免费观看的视频| 亚洲va日本ⅴa欧美va伊人久久| 日本免费一区二区三区高清不卡| 国产三级黄色录像| 一级作爱视频免费观看| 99国产精品99久久久久| 波多野结衣高清无吗| 亚洲av成人av| 国产真实乱freesex| 久久热在线av| 欧美日韩黄片免| 精品国产美女av久久久久小说| 亚洲精品中文字幕一二三四区| 日本黄色片子视频| 一个人看视频在线观看www免费 | 岛国在线免费视频观看| 人妻久久中文字幕网| 好看av亚洲va欧美ⅴa在| 神马国产精品三级电影在线观看| 天天添夜夜摸| 成人鲁丝片一二三区免费| 国产成人精品久久二区二区免费| 久久欧美精品欧美久久欧美| 听说在线观看完整版免费高清| 18禁黄网站禁片午夜丰满| 两个人看的免费小视频| av在线蜜桃| 丰满人妻一区二区三区视频av | 欧美日韩一级在线毛片| 99久久精品一区二区三区| 亚洲一区二区三区色噜噜| 欧美午夜高清在线| 久久性视频一级片| 国产精品电影一区二区三区| 日韩欧美国产一区二区入口| 亚洲九九香蕉| 亚洲av电影不卡..在线观看| 一个人看的www免费观看视频| 精华霜和精华液先用哪个| 久久中文字幕人妻熟女| 亚洲第一欧美日韩一区二区三区| 曰老女人黄片| 日韩欧美三级三区| 亚洲成人久久爱视频| 日日夜夜操网爽| 老汉色∧v一级毛片| 88av欧美| 国产成人aa在线观看| 日韩大尺度精品在线看网址| 久久精品人妻少妇| 香蕉丝袜av| 美女高潮的动态| 久久久久免费精品人妻一区二区| 成人三级黄色视频| 久久亚洲真实| 国产高清三级在线| 日韩人妻高清精品专区| 亚洲九九香蕉| 日本三级黄在线观看| 国产精品1区2区在线观看.| 最新中文字幕久久久久 | 久久精品亚洲精品国产色婷小说| 国产一区在线观看成人免费| 免费看十八禁软件| 午夜久久久久精精品| 精品不卡国产一区二区三区| 一二三四社区在线视频社区8| 男人和女人高潮做爰伦理| 91老司机精品| www.自偷自拍.com| 91在线精品国自产拍蜜月 | 又黄又爽又免费观看的视频| 免费电影在线观看免费观看| 国产成人aa在线观看| 在线国产一区二区在线| 精华霜和精华液先用哪个| 99国产精品一区二区蜜桃av| 制服人妻中文乱码| 91字幕亚洲| 亚洲无线在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产极品精品免费视频能看的| www.999成人在线观看| 精品熟女少妇八av免费久了| 国产伦一二天堂av在线观看| 丰满的人妻完整版| 99久久国产精品久久久| 99国产精品一区二区三区| 午夜福利欧美成人| 国产亚洲av嫩草精品影院| 一级毛片女人18水好多| 日韩欧美国产一区二区入口| 老鸭窝网址在线观看| 亚洲天堂国产精品一区在线| 人妻夜夜爽99麻豆av| 亚洲精品粉嫩美女一区| 色综合欧美亚洲国产小说| 国产黄片美女视频| 在线免费观看的www视频| 亚洲欧美日韩高清专用| 国产亚洲精品一区二区www| 日韩成人在线观看一区二区三区| 这个男人来自地球电影免费观看| 国产av麻豆久久久久久久| 嫩草影视91久久| 丁香欧美五月| 亚洲欧美一区二区三区黑人| 国产精品 欧美亚洲| 亚洲国产高清在线一区二区三| 亚洲国产欧洲综合997久久,| 国产97色在线日韩免费| 老熟妇乱子伦视频在线观看| 国产三级黄色录像| 国产成人啪精品午夜网站| 亚洲电影在线观看av| 亚洲熟妇熟女久久| 99热6这里只有精品| 久久精品亚洲精品国产色婷小说| 18禁国产床啪视频网站| e午夜精品久久久久久久| 国产精品乱码一区二三区的特点| 国产主播在线观看一区二区| 99精品久久久久人妻精品| 91九色精品人成在线观看| 少妇丰满av| www.www免费av| 高潮久久久久久久久久久不卡| 好男人电影高清在线观看| 免费大片18禁| 国产高清三级在线| 大型黄色视频在线免费观看| 亚洲黑人精品在线| 国产av麻豆久久久久久久| 亚洲欧美日韩无卡精品| 一本久久中文字幕| 又黄又粗又硬又大视频| 美女高潮喷水抽搐中文字幕| 国产精品亚洲av一区麻豆| 手机成人av网站| 日本精品一区二区三区蜜桃| 99久久精品国产亚洲精品| 免费无遮挡裸体视频| 欧美黄色片欧美黄色片| 中文字幕精品亚洲无线码一区| 亚洲自偷自拍图片 自拍| 国产伦精品一区二区三区视频9 | 中出人妻视频一区二区| 亚洲一区二区三区色噜噜| 嫩草影院精品99| 一卡2卡三卡四卡精品乱码亚洲| 夜夜看夜夜爽夜夜摸| 亚洲精品粉嫩美女一区| 精品日产1卡2卡| 日本黄色视频三级网站网址| 观看美女的网站| 一卡2卡三卡四卡精品乱码亚洲| 亚洲专区字幕在线| 网址你懂的国产日韩在线| 美女高潮的动态| 日本熟妇午夜| 成人三级黄色视频| 一个人观看的视频www高清免费观看 | 天堂动漫精品| 欧美成人性av电影在线观看| 一本久久中文字幕| 精品久久久久久久久久久久久| 欧美乱色亚洲激情| 国产三级中文精品| 淫秽高清视频在线观看| 国产成人aa在线观看| 看片在线看免费视频| 久久欧美精品欧美久久欧美| 黄色视频,在线免费观看| 一进一出抽搐动态| 欧美不卡视频在线免费观看| 久久久久久久久中文| 国产精品国产高清国产av| 久久久久久九九精品二区国产| 欧美日韩乱码在线| 制服人妻中文乱码| 一a级毛片在线观看| 久9热在线精品视频| 国产精品亚洲美女久久久| 美女黄网站色视频| 99久国产av精品| 久久精品91蜜桃| 99国产综合亚洲精品| 亚洲第一电影网av| 午夜福利18| 高清毛片免费观看视频网站| 免费看光身美女| 美女黄网站色视频| 女人被狂操c到高潮| 全区人妻精品视频| 俄罗斯特黄特色一大片| 99精品在免费线老司机午夜| 亚洲美女视频黄频| 成年免费大片在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产伦人伦偷精品视频| 国产精品一区二区免费欧美| 国产精品香港三级国产av潘金莲| 国产麻豆成人av免费视频| 亚洲av五月六月丁香网| 最近视频中文字幕2019在线8| www.熟女人妻精品国产| 婷婷精品国产亚洲av在线| 欧美xxxx黑人xx丫x性爽| 亚洲中文日韩欧美视频| 久久中文看片网| 日本黄色片子视频| 成人国产一区最新在线观看| 禁无遮挡网站| 精品福利观看| 欧美日韩乱码在线| 亚洲国产欧美网| 他把我摸到了高潮在线观看| 欧美极品一区二区三区四区| 男女午夜视频在线观看| 99精品久久久久人妻精品| h日本视频在线播放| 精品久久久久久久毛片微露脸| 国产久久久一区二区三区| 大型黄色视频在线免费观看| 国模一区二区三区四区视频 | 一个人免费在线观看电影 | 日韩欧美免费精品| 丁香六月欧美| 久久久久久久久久黄片| 法律面前人人平等表现在哪些方面| 国内少妇人妻偷人精品xxx网站 | 国产精品影院久久| 99在线视频只有这里精品首页| 国内精品一区二区在线观看| 欧美性猛交╳xxx乱大交人| 欧美大码av| 日韩成人在线观看一区二区三区| 国内精品久久久久久久电影| 老司机在亚洲福利影院| a级毛片a级免费在线| 午夜久久久久精精品| 精品久久久久久,| 真人一进一出gif抽搐免费| 欧美激情久久久久久爽电影| 听说在线观看完整版免费高清| 淫妇啪啪啪对白视频| 亚洲精品一区av在线观看| 国内毛片毛片毛片毛片毛片| 成人一区二区视频在线观看| 国产精品,欧美在线| 这个男人来自地球电影免费观看| 国产精品99久久久久久久久| 九九热线精品视视频播放| 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆 | 久久精品91蜜桃| 亚洲aⅴ乱码一区二区在线播放| 亚洲中文日韩欧美视频| 国产精品九九99| 国产又黄又爽又无遮挡在线| 久久久成人免费电影| 久久久国产成人精品二区| 国内少妇人妻偷人精品xxx网站 | 国产成人系列免费观看| 成年免费大片在线观看| 精品人妻1区二区| 欧美成狂野欧美在线观看| 天堂av国产一区二区熟女人妻| 久久精品国产亚洲av香蕉五月| 国产69精品久久久久777片 | 欧美极品一区二区三区四区| 亚洲成人中文字幕在线播放| 欧美黄色淫秽网站| 亚洲色图av天堂| 国产私拍福利视频在线观看| 日韩三级视频一区二区三区| 国产精品电影一区二区三区| 色老头精品视频在线观看| 国产av一区在线观看免费| 麻豆久久精品国产亚洲av| 久久精品aⅴ一区二区三区四区| 99久久无色码亚洲精品果冻| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 久久久久久久午夜电影| 日本三级黄在线观看| 国产精品1区2区在线观看.| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全免费视频| 最近最新中文字幕大全电影3| av天堂在线播放| 一本久久中文字幕| 免费看光身美女| 美女免费视频网站| av视频在线观看入口| 特大巨黑吊av在线直播| 国产精品98久久久久久宅男小说| 丁香欧美五月| 免费观看的影片在线观看| 亚洲av中文字字幕乱码综合| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人一区二区免费高清观看 | 999久久久国产精品视频| 亚洲无线在线观看| 精品久久久久久成人av| 亚洲电影在线观看av| 天天添夜夜摸| 丰满人妻熟妇乱又伦精品不卡| 国产高清视频在线播放一区| 国产高清有码在线观看视频| 亚洲第一电影网av| 国产黄片美女视频| 欧美激情在线99| 欧美日韩国产亚洲二区| 国产精品自产拍在线观看55亚洲| 亚洲成av人片在线播放无| 国产亚洲欧美在线一区二区| 成人精品一区二区免费| 女生性感内裤真人,穿戴方法视频| 国产又黄又爽又无遮挡在线| 亚洲,欧美精品.| av黄色大香蕉| 91麻豆精品激情在线观看国产| 亚洲avbb在线观看| 久久精品aⅴ一区二区三区四区| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 老司机午夜福利在线观看视频| 日本与韩国留学比较| 国产成人aa在线观看| 每晚都被弄得嗷嗷叫到高潮| 一进一出抽搐动态| 欧美三级亚洲精品| 亚洲片人在线观看| 黄片大片在线免费观看| 国产97色在线日韩免费| 国产伦人伦偷精品视频| 久久久久九九精品影院| 中文字幕精品亚洲无线码一区| 国产熟女xx| 一二三四社区在线视频社区8| 老司机午夜福利在线观看视频| 成人国产综合亚洲| 999久久久精品免费观看国产| 国产免费av片在线观看野外av| 后天国语完整版免费观看| 午夜福利在线观看免费完整高清在 | 免费无遮挡裸体视频| 欧美激情在线99| 又黄又爽又免费观看的视频| 国产精品 国内视频| 天天一区二区日本电影三级| 99久久无色码亚洲精品果冻| 女同久久另类99精品国产91| 久久久久九九精品影院| 色哟哟哟哟哟哟| 青草久久国产| 国产99白浆流出| 又大又爽又粗| 亚洲电影在线观看av| 男插女下体视频免费在线播放| 又紧又爽又黄一区二区| 两性夫妻黄色片| 国产极品精品免费视频能看的| 亚洲成av人片免费观看| 国产爱豆传媒在线观看| 三级毛片av免费| 好男人电影高清在线观看| 麻豆国产av国片精品| 中文亚洲av片在线观看爽| 久久天躁狠狠躁夜夜2o2o| 又紧又爽又黄一区二区| 亚洲精华国产精华精| 婷婷精品国产亚洲av| 国产精品 国内视频| 三级男女做爰猛烈吃奶摸视频| 最近最新免费中文字幕在线| 日本黄色视频三级网站网址| 国产真实乱freesex| 18禁观看日本| 亚洲一区高清亚洲精品| 亚洲成人中文字幕在线播放| 国内久久婷婷六月综合欲色啪| 久久婷婷人人爽人人干人人爱| 日本三级黄在线观看| 国产欧美日韩一区二区三| 黄色丝袜av网址大全| 免费av毛片视频| 黑人欧美特级aaaaaa片| 长腿黑丝高跟| 看黄色毛片网站| 叶爱在线成人免费视频播放| 色老头精品视频在线观看| 夜夜爽天天搞| 欧美三级亚洲精品| tocl精华| 好男人电影高清在线观看| 中文字幕精品亚洲无线码一区| 亚洲av成人不卡在线观看播放网| 很黄的视频免费| 国产精品久久久久久亚洲av鲁大| 久久草成人影院| av国产免费在线观看| 99热这里只有精品一区 | 很黄的视频免费| 欧美+亚洲+日韩+国产| 久久久久性生活片| 美女午夜性视频免费| 精品久久久久久成人av| 国产成人欧美在线观看| 国产真人三级小视频在线观看| 亚洲精品乱码久久久v下载方式 | a级毛片在线看网站| 变态另类成人亚洲欧美熟女| 国产欧美日韩一区二区精品| 亚洲真实伦在线观看| 女人被狂操c到高潮| 少妇丰满av| 亚洲天堂国产精品一区在线| 免费高清视频大片| 人人妻人人看人人澡| 观看美女的网站| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久人妻精品电影| 老熟妇乱子伦视频在线观看| 国产亚洲av嫩草精品影院| 老熟妇仑乱视频hdxx| 非洲黑人性xxxx精品又粗又长| 日本与韩国留学比较| 亚洲avbb在线观看| 国产蜜桃级精品一区二区三区| 久久午夜综合久久蜜桃| 亚洲精品粉嫩美女一区| 搡老熟女国产l中国老女人| 18禁美女被吸乳视频| 亚洲午夜精品一区,二区,三区| 18禁裸乳无遮挡免费网站照片| 小说图片视频综合网站| 操出白浆在线播放| av国产免费在线观看| 午夜成年电影在线免费观看| 久久精品国产综合久久久| 亚洲性夜色夜夜综合| 一本综合久久免费| 18美女黄网站色大片免费观看| 国内精品美女久久久久久| 女警被强在线播放| 少妇人妻一区二区三区视频| 日韩欧美国产在线观看| 久久久久久国产a免费观看| 久久精品亚洲精品国产色婷小说| 中国美女看黄片| 一区二区三区高清视频在线| 久久精品国产99精品国产亚洲性色| 啦啦啦免费观看视频1| 午夜久久久久精精品| 久久香蕉国产精品| 欧美精品啪啪一区二区三区| 亚洲五月婷婷丁香| 日韩人妻高清精品专区| 色av中文字幕| 亚洲精品色激情综合| 亚洲av成人一区二区三| 啪啪无遮挡十八禁网站| 99久久久亚洲精品蜜臀av| xxxwww97欧美| 他把我摸到了高潮在线观看| 99国产精品99久久久久| 亚洲国产精品成人综合色| 亚洲中文字幕一区二区三区有码在线看 | 色老头精品视频在线观看| 亚洲五月婷婷丁香| 免费看a级黄色片| 精品久久久久久成人av| 我要搜黄色片| 国产欧美日韩精品一区二区| 亚洲,欧美精品.| 成人18禁在线播放| 男人舔奶头视频| 久久精品91无色码中文字幕| 免费看光身美女| 亚洲av美国av| 亚洲成a人片在线一区二区| 最好的美女福利视频网| 中国美女看黄片| 99热精品在线国产| 99精品欧美一区二区三区四区| 国产综合懂色| 久久精品国产亚洲av香蕉五月| 色av中文字幕| 欧美av亚洲av综合av国产av| 亚洲欧美激情综合另类| 久久精品人妻少妇| 女警被强在线播放| 此物有八面人人有两片| 午夜激情欧美在线| 亚洲av成人不卡在线观看播放网| 热99re8久久精品国产| 午夜福利高清视频| 熟女人妻精品中文字幕| 国产三级中文精品| 免费在线观看亚洲国产| 神马国产精品三级电影在线观看| 欧美日韩一级在线毛片| 999精品在线视频| 真人一进一出gif抽搐免费| 夜夜看夜夜爽夜夜摸| 一个人观看的视频www高清免费观看 | 看片在线看免费视频| 综合色av麻豆| 国产激情偷乱视频一区二区| 99热只有精品国产| 国产综合懂色| 亚洲,欧美精品.| 亚洲国产日韩欧美精品在线观看 | 久久久久精品国产欧美久久久| 麻豆国产97在线/欧美| 亚洲国产精品合色在线| 后天国语完整版免费观看| 亚洲成人中文字幕在线播放| 看免费av毛片| 狠狠狠狠99中文字幕| 两个人的视频大全免费| 亚洲成a人片在线一区二区| 1024香蕉在线观看| 免费一级毛片在线播放高清视频| 69av精品久久久久久| 亚洲熟女毛片儿| 成人鲁丝片一二三区免费| 在线免费观看的www视频| 欧美乱色亚洲激情| 久久久色成人| 少妇丰满av| 看片在线看免费视频| 国产精品 国内视频| 此物有八面人人有两片| 中文在线观看免费www的网站| 久久中文字幕一级| 国产午夜福利久久久久久| 亚洲精品456在线播放app | 不卡av一区二区三区| 国产亚洲av嫩草精品影院| 亚洲国产精品sss在线观看| 欧美一级a爱片免费观看看| 欧美成人性av电影在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲熟妇熟女久久| 色老头精品视频在线观看| 国产精品美女特级片免费视频播放器 | 一进一出抽搐动态| 国产精品久久久久久精品电影| 亚洲成av人片在线播放无| 嫁个100分男人电影在线观看| 亚洲精品一区av在线观看| 国产在线精品亚洲第一网站| 51午夜福利影视在线观看| 久久久久久人人人人人| 色尼玛亚洲综合影院| 天天添夜夜摸| 国产99白浆流出| 成人亚洲精品av一区二区| 成人特级av手机在线观看| 黄色成人免费大全| x7x7x7水蜜桃| 国产成年人精品一区二区| 亚洲国产精品合色在线| 国产成人一区二区三区免费视频网站| 男女做爰动态图高潮gif福利片| 91老司机精品| 露出奶头的视频| 91九色精品人成在线观看| 偷拍熟女少妇极品色| 亚洲av五月六月丁香网| 午夜福利在线观看吧| 亚洲电影在线观看av| bbb黄色大片| 亚洲真实伦在线观看| 国产69精品久久久久777片 | 老司机福利观看| 精品国产三级普通话版| 国产aⅴ精品一区二区三区波| 久久久水蜜桃国产精品网| 欧美高清成人免费视频www| 日本a在线网址| 国产亚洲av高清不卡| 一本一本综合久久| 男人舔奶头视频| 老司机午夜福利在线观看视频| 啦啦啦免费观看视频1| www.熟女人妻精品国产| 香蕉av资源在线| 国产伦精品一区二区三区四那| 国产精品女同一区二区软件 | 色哟哟哟哟哟哟| 香蕉国产在线看| 真实男女啪啪啪动态图| 亚洲美女黄片视频| 欧美大码av| 听说在线观看完整版免费高清| 日韩精品中文字幕看吧| 青草久久国产| 国产精品九九99| 欧美乱色亚洲激情| 欧美xxxx黑人xx丫x性爽|