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

    油氣藏蓋層封閉性研究現(xiàn)狀及未來(lái)發(fā)展趨勢(shì)

    2018-06-22 05:50:58付曉飛呂延防柳少波盧明旭東北石油大學(xué)非常規(guī)油氣成藏與開(kāi)發(fā)省部共建國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地黑龍江大慶1618黑龍江省高??萍紕?chuàng)新團(tuán)隊(duì)斷層變形封閉性及與流體運(yùn)移黑龍江大慶1618中國(guó)石油勘探開(kāi)發(fā)研究院石油地質(zhì)實(shí)驗(yàn)研究中心北京10008中國(guó)石油大慶油田有限責(zé)任公司第四采油廠黑龍江大慶16511
    石油與天然氣地質(zhì) 2018年3期
    關(guān)鍵詞:蓋層脆性泥巖

    付曉飛,吳 桐,呂延防,柳少波,田 華,盧明旭(1.東北石油大學(xué) 非常規(guī)油氣成藏與開(kāi)發(fā)省部共建國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地,黑龍江 大慶 1618;2.黑龍江省高??萍紕?chuàng)新團(tuán)隊(duì)“斷層變形、封閉性及與流體運(yùn)移”,黑龍江 大慶 1618;.中國(guó)石油 勘探開(kāi)發(fā)研究院 石油地質(zhì)實(shí)驗(yàn)研究中心,北京 10008;.中國(guó)石油 大慶油田有限責(zé)任公司 第四采油廠,黑龍江 大慶 16511)

    油氣藏蓋層封閉性研究是伴隨石油地質(zhì)理論的發(fā)展而逐漸發(fā)展起來(lái)的,大致可分為3個(gè)階段(圖1)。油氣封蓋概念建立階段(1966年以前),油氣勘探由露頭區(qū)轉(zhuǎn)入覆蓋區(qū),背斜聚油理論指導(dǎo)油氣勘探,認(rèn)識(shí)到浮力是油氣運(yùn)移的主要?jiǎng)恿?,泥巖蓋層在油氣成藏中起到重要的作用,以1966年Smith提出蓋層和斷層封閉性概念為標(biāo)志,開(kāi)始探索蓋層封閉性和斷層封閉性問(wèn)題[1]。蓋層封閉機(jī)理探索階段(1966—1996),構(gòu)造圈閉的復(fù)雜性和巖性圈閉的隱蔽性導(dǎo)致鉆探失利,促使勘探家考慮油氣運(yùn)聚后的保存問(wèn)題,以1987年Watts提出蓋層毛細(xì)管封閉機(jī)理和水力封閉機(jī)理為標(biāo)志[2],開(kāi)始系統(tǒng)研究蓋層和斷層封閉機(jī)理并進(jìn)行定性評(píng)價(jià)。1983美國(guó)地質(zhì)學(xué)家協(xié)會(huì)年(AAPG)召開(kāi)蓋層和油氣保存方面的專題研討會(huì),Downey認(rèn)為蓋層研究要從微觀和宏觀兩方面進(jìn)行[3];Grunau認(rèn)為頁(yè)巖和蒸發(fā)巖是有效蓋層[4],沉積環(huán)境、擴(kuò)散速率和裂縫形成演化是蓋層評(píng)價(jià)的重要參數(shù)。1996年呂延防等人又提出了超壓封閉機(jī)理和烴濃度封閉機(jī)理[5],針對(duì)不同封閉機(jī)理建立了利用排替壓力、有效應(yīng)力和抗張強(qiáng)度、異??紫读黧w壓力和擴(kuò)散系數(shù)評(píng)價(jià)封閉能力方法,依據(jù)蓋層宏觀發(fā)育特征和微觀封閉能力,建立了蓋層封閉性評(píng)價(jià)行業(yè)標(biāo)準(zhǔn)。蓋層與斷層封閉性定量評(píng)價(jià)階段(1997—2013),斷層附近大量有效的經(jīng)濟(jì)目標(biāo)發(fā)現(xiàn)和鉆探失利形成的矛盾,使地質(zhì)學(xué)家開(kāi)始考慮如何定量評(píng)價(jià)斷層和蓋層封閉性。以呂延防出版的“油氣藏封蓋研究”為標(biāo)志,建立了蓋層封閉能力定量評(píng)價(jià)方法[5],初步建立了蓋層完整性評(píng)價(jià)思路[6-11],指出蓋層完整性破壞因素有構(gòu)造破裂(斷裂錯(cuò)斷、亞地震斷層與砂體匹配、構(gòu)造裂縫)和水力破裂作用。認(rèn)識(shí)到蓋巖脆-韌性特征是控制蓋層完整性的關(guān)鍵因素,提出了利用密度-應(yīng)變[7,10]、BRI和OCR[8-9,12]、Byerlee摩擦定律和Goetze’s準(zhǔn)則定量判斷泥頁(yè)巖和膏鹽巖脆-韌性轉(zhuǎn)換定量評(píng)價(jià)標(biāo)準(zhǔn)[13]。蓋層封閉能力動(dòng)態(tài)演化研究階段(2014—),伴隨常規(guī)油氣由中淺層轉(zhuǎn)向深層,多旋回盆地油氣多期成藏和改造迫切需要開(kāi)展蓋層封閉能力動(dòng)態(tài)演化研究,加之非常規(guī)油氣勘探開(kāi)發(fā),在蓋巖脆-韌性定量表征基礎(chǔ)上,金之鈞提出了泥頁(yè)巖蓋層封閉能力動(dòng)態(tài)演化過(guò)程評(píng)價(jià)方法[14]。此外,地球物理新技術(shù)在蓋層完整性評(píng)價(jià)中也逐漸得到應(yīng)用[15]。迄今為止,針對(duì)蓋層封閉能力的研究在封閉機(jī)理和評(píng)價(jià)方法上都取得了長(zhǎng)足的進(jìn)展,勘探目標(biāo)也由常規(guī)油氣藏向非常規(guī)和深層油氣藏?cái)U(kuò)展,未來(lái)對(duì)蓋層封閉能力的研究必定從定性轉(zhuǎn)為定量評(píng)價(jià),從靜態(tài)轉(zhuǎn)為動(dòng)態(tài)評(píng)價(jià),從過(guò)去僅依據(jù)宏觀發(fā)育特征和微觀封閉機(jī)理評(píng)價(jià)轉(zhuǎn)為從潤(rùn)濕性改變、構(gòu)造破裂和水力破裂等對(duì)不同成巖程度和不同脆韌性蓋層完整性的評(píng)價(jià)。

    圖1 油氣藏蓋層封閉性理論研究進(jìn)展Fig.1 Theoretical research progress in reservoir caprock sealing properties

    1 蓋層發(fā)育特征及封閉能力

    1.1 巖性

    統(tǒng)計(jì)表明,大中型油氣田蓋層主要巖性為泥巖、膏鹽巖和碳酸鹽巖,泥頁(yè)巖最為常見(jiàn)[16-17](圖2)。不同巖性蓋層封閉的最大烴柱高度存在明顯差異[6](圖3a),膏巖封閉烴柱高度最大,其次為泥巖和碳酸鹽巖,砂巖最小。我國(guó)陸相盆地不同巖性蓋層排替壓力也存在差異,封閉能力依然遵循膏巖、泥巖、碳酸鹽巖和砂巖依次變差的規(guī)律(圖3b)。但膏巖、泥巖和碳酸鹽巖均能封閉住幾千米的烴柱高度,因此,巖性不是蓋層封閉性定量評(píng)價(jià)的關(guān)鍵參數(shù)。

    1.2 蓋層分布范圍

    蓋層按分布范圍可分為區(qū)域性蓋層和局部性蓋層[5]。區(qū)域性蓋層遍布在含油氣盆地或坳陷的大部分地區(qū),厚度大、面積廣且分布穩(wěn)定;局部性蓋層指分布在局部構(gòu)造上的蓋層。深湖-半深湖相、潟湖相通常發(fā)育質(zhì)純、厚度較大、分布面積廣的區(qū)域性蓋層[5],河流相、三角州相、濱淺湖相等通常形成砂-泥薄互層、厚度較小、分布面積小的局部性蓋層。因此說(shuō)沉積環(huán)境是決定蓋層宏觀發(fā)育程度的關(guān)鍵因素。僅從蓋層和烴源巖分布范圍看,蓋層可分為三級(jí):好的蓋層分布范圍大于烴源巖的范圍,中等的蓋層分布范圍與烴源巖的面積相當(dāng),差的蓋層分布范圍遠(yuǎn)小于烴源巖面積。

    1.3 蓋層厚度

    從毛細(xì)管封閉機(jī)理講,蓋層厚度與其封閉能力無(wú)直接的函數(shù)關(guān)系[2,18]。但一些研究者認(rèn)為蓋層厚度與烴柱高度之間存在某些聯(lián)系[19],即蓋層厚度影響其空間分布[5],厚度越大,沉積環(huán)境越穩(wěn)定,均質(zhì)性越好,橫向分布的面積越廣,欠壓實(shí)并導(dǎo)致超壓形成的可能性越大,斷裂越不容易導(dǎo)致蓋層滲漏。從厚度影響分布的角度看,厚度越大,越有利于油氣的保存。因此,無(wú)論是氣田還是油田,直接蓋層的厚度與可封蓋的最大烴柱高度總體呈現(xiàn)正相關(guān)關(guān)系(圖4)。

    1.4 單層厚度和蓋地比

    單層厚度和蓋地比是反映蓋層“純度”的參數(shù),同一層位的蓋層單層厚度越大,蓋地比越高,縱向連續(xù)分布蓋層的厚度越大,封閉能力越強(qiáng)[20]。King利用逾滲理論研究了疊置砂體間連通性問(wèn)題(不考慮斷裂),發(fā)現(xiàn)存在一個(gè)砂地比特征門(mén)限值,低于該門(mén)限值砂體之間基本不連通[21]。Jackson等人提出一個(gè)薄層砂泥巖互層的儲(chǔ)層模型[22],估計(jì)水平的逾滲閾值為0.28(C0),垂直的逾滲閾值為0.5(C)。這意味著泥質(zhì)巖蓋層的蓋地比為35%~50%時(shí),垂向不連通;統(tǒng)計(jì)中西部前陸盆地氣藏之上蓋層的泥地比發(fā)現(xiàn),封住天然氣所需最小的泥地比為40%~60%[20],當(dāng)有小斷層和裂縫存在的時(shí)候,最小的泥地比可能升高。根據(jù)蓋地比蓋層品質(zhì)可分為三級(jí)(圖5):蓋地比小于1-C的蓋層為封閉能力差的蓋層,蓋地比介于1-C和1-C0時(shí)為封閉能力中等的蓋層,蓋地比大于1-C0時(shí)為封閉能力強(qiáng)的蓋層。

    1.5 蓋層的巖石力學(xué)特征

    泥巖脆性指數(shù)(BRI)、超固結(jié)比(OCR)、巖石內(nèi)摩擦角和巖石密度是蓋層脆-韌性評(píng)價(jià)的主要參數(shù)。脆性指數(shù)為BRI=UCS/UCSNC,UCSNC是正常固結(jié)巖石的單軸抗壓強(qiáng)度,UCS是超固結(jié)巖石的單軸抗壓強(qiáng)度,當(dāng)BRI大于2時(shí),巖石更容易發(fā)生脆性破裂[9](圖6)。巖石在地質(zhì)歷史時(shí)期受到過(guò)的最大有效應(yīng)力稱為前固結(jié)應(yīng)力,而前固結(jié)應(yīng)力與現(xiàn)在有效應(yīng)力之比定義為超固結(jié)比。Nygard等人利用OCR對(duì)北海盆地的泥巖進(jìn)行評(píng)價(jià),當(dāng)OCR大于2.5時(shí),巖石發(fā)生滲漏的風(fēng)險(xiǎn)增加[12]。巖石比表面積是指單位重量?jī)?nèi)巖石內(nèi)表面積和外表面積之和,它是巖石粒徑、孔隙度、壓實(shí)程度和膠結(jié)程度的綜合反映,巖石比表面積增加,內(nèi)摩擦角減小,更易表現(xiàn)為韌性,表面積大于300 m2/g指示了膠結(jié)作用較弱,為韌性巖石。大多數(shù)沉積盆地內(nèi)的頁(yè)巖密度小于2.1 g/cm3,只發(fā)生韌性變形;若密度大于2.1 g/cm3,在足夠的應(yīng)變下將發(fā)生脆性破裂[6]。脆性巖石在變形過(guò)程中易于形成裂縫,與蓋層內(nèi)薄砂體相配合,形成連續(xù)的滲漏通道。研究表明,當(dāng)砂泥薄互層蓋層中裂縫的數(shù)量超過(guò)泥巖層數(shù)量的5倍時(shí),由砂巖層和裂縫構(gòu)成的“之字型”運(yùn)移通道導(dǎo)致油氣穿越蓋層垂向運(yùn)移[9]。

    圖3 不同巖性蓋層封閉能力評(píng)價(jià)參數(shù)對(duì)比[6]Fig.3 Parameter comparison for various lithological caprock sealing capacitiesa.世界不同巖性蓋層封閉烴柱高度對(duì)比;b.中國(guó)不同巖性蓋層排替壓力對(duì)比[6]

    圖4 直接蓋層厚度與封閉烴柱高度的關(guān)系Fig.4 Relationship between the thickness of direct caprock and the height of hydrocarbon columna.油田;b.氣田

    圖5 渤海灣盆地王家崗油田砂地比與砂體幾何連通概率關(guān)系模型[23]Fig.5 Probability model of relationship between sand-stratum ratio and sand body geometric connectivity of Wang jiagang oilfield,Bohai Bay Basin[23]

    圖6 頁(yè)巖脆性指數(shù)及發(fā)生破裂的臨界條件[9]Fig.6 Brittle indexes and fracture critical conditions of shale[9]

    2 蓋層封閉機(jī)理及封閉能力動(dòng)態(tài)演化過(guò)程

    2.1 毛細(xì)管封閉能力及定量評(píng)價(jià)

    蓋層多為水潤(rùn)濕或水飽和的,且具有較高的毛細(xì)管壓力,因此能封閉住一定的烴柱高度[2]。毛細(xì)管壓力定義為在儲(chǔ)-蓋層接觸面上油壓力和水壓力的差異。只有當(dāng)油氣柱產(chǎn)生的浮壓超過(guò)毛細(xì)管進(jìn)入壓力時(shí),毛細(xì)管封閉失效,油氣突破蓋層垂向運(yùn)移[2](圖7)。

    表征毛細(xì)管封閉能力的參數(shù)包括孔隙度、滲透率、孔喉半徑和排替壓力[5]。泥頁(yè)巖原始孔隙度可達(dá)60%~80%,隨著埋深的增加而逐漸減小,埋深為0.5~1 km后下降到35%,埋深為3.5~4.0 km下降到5%左右(圖8a)。Khanin測(cè)試泥質(zhì)巖滲透率為(10-6~10-2)×10-3μm,滲透率隨著埋深、有效應(yīng)力增加或泥質(zhì)含量的增加而逐漸降低[24](圖8b)。滲透率和孔隙度之間沒(méi)有明顯的正相關(guān)關(guān)系[25-26]。北海盆地白堊蓋層孔喉半徑為5~160 nm,且隨著埋藏深度增加而逐漸減小,超過(guò)3 000 m普遍小于30 nm[27]。松遼盆地深層和中淺層泥巖蓋層孔喉半徑存在差異[28](圖9),中淺層泥巖孔喉半徑為0.8~200 nm,峰值普遍為1~2 nm,深層泥巖孔喉半徑為0.8~40 nm,峰值范圍為1~6 nm。排替壓力定義為潤(rùn)濕性流體被非潤(rùn)濕性流體驅(qū)替所需的最小壓力[5],可以通過(guò)實(shí)驗(yàn)法、壓汞法、測(cè)井和地震資料等手段獲得。排替壓力與孔隙度[14]、滲透率[8]和泥質(zhì)含量[5]之間存在定量關(guān)系,是定量預(yù)測(cè)蓋層封閉能力演化過(guò)程的基礎(chǔ)[14]。

    影響蓋層毛細(xì)管封閉能力的因素主要有4方面。一是成巖程度影響,泥巖壓實(shí)歷經(jīng)快速壓實(shí)、穩(wěn)定壓實(shí)、突變壓實(shí)和緊密壓實(shí)四個(gè)階段[28],封閉能力在突變壓實(shí)和緊密壓實(shí)階段最強(qiáng)。深層泥質(zhì)巖蓋層處于緊密壓實(shí)階段,對(duì)應(yīng)中成巖晚期和晚成巖階段,化學(xué)膠結(jié)作用強(qiáng)烈,孔喉半徑很小,封閉能力增強(qiáng),但泥質(zhì)巖脆性增強(qiáng),裂縫導(dǎo)致滲漏的可能性增大[2,8,24]。二是流體性質(zhì)影響,氣的界面張力隨深度增加的速率與油不同,相同蓋層封閉氣可能多于油。深層處于高溫高壓環(huán)境,多以凝析油和氣為主,對(duì)蓋層封閉能力影響需要進(jìn)一步探討。三是界面張力影響,油-水界面張力隨著溫度增加而減小,隨著壓力增加變化不大;氣-水界面張力隨溫度和壓力升高而降低,深層蓋層封閉油氣能力受界面張力影響而降低。四是潤(rùn)濕性影響,蓋層能夠封閉住油氣是因其為水潤(rùn)濕的,但原油中的極性化合物可改變巖石的潤(rùn)濕性[29-31],如顆粒表面吸附的有機(jī)金屬?gòu)?fù)合物導(dǎo)致油潤(rùn)濕,去掉后變?yōu)樗疂?rùn)濕[31]。實(shí)驗(yàn)證實(shí)束縛水可以通過(guò)油層進(jìn)入蓋層[32-33],酸性化合物從油中分離出進(jìn)入殘留水中,水?dāng)y帶酸性化合物進(jìn)入蓋層中,導(dǎo)致水膜破裂,吸附到礦物表面,導(dǎo)致局部變?yōu)橛蜐?rùn)濕,減小毛細(xì)管壓力,從而形成不連續(xù)的微滲漏空間。五是儲(chǔ)層超壓的影響,當(dāng)儲(chǔ)層存在超壓時(shí),蓋層封閉能力減小,相當(dāng)于蓋層毛細(xì)管壓力要抵消儲(chǔ)層超壓。1998年Bj?rkum等人提出不同觀點(diǎn):認(rèn)為連續(xù)水相可以從油層向上覆蓋層中運(yùn)移[34]。油層中水相和蓋層中水相壓力差是極小的,因此封閉層之下烴類并不受儲(chǔ)層超壓的影響,在水潤(rùn)濕的儲(chǔ)層中的水相是流動(dòng)的,烴類界面之下水層的超壓可以傳遞到儲(chǔ)層-封閉層的界面,因此超壓水潤(rùn)濕性儲(chǔ)層并不比正常壓力儲(chǔ)層更容易發(fā)生毛細(xì)管滲漏。

    圖7 蓋層封閉機(jī)理模式[2]Fig.7 Schematic diagram of caprock sealing mechanism pattern[2]

    圖8 不同類型泥巖孔隙度(a)與滲透率(b)隨深度的變化規(guī)律Fig.8 Variation of porosity and permeability with depth for different mudstones

    泥質(zhì)巖蓋層在埋藏-抬升過(guò)程中,物性和封閉能力呈現(xiàn)動(dòng)態(tài)變化規(guī)律,埋藏過(guò)程中,泥質(zhì)巖蓋層孔隙度逐漸降低,封閉能力逐漸增強(qiáng)[14];抬升過(guò)程中由于裂縫產(chǎn)生[35],孔隙度變化不大,但滲透率明顯增加?;谂盘鎵毫涂紫抖?、滲透率之間的定量關(guān)系,Jin等人建立了泥質(zhì)巖蓋層封閉能力動(dòng)態(tài)變化過(guò)程定量評(píng)價(jià)方法[14](圖10)。

    2.2 水力封閉機(jī)理及封閉能力演化過(guò)程

    水力破裂指由于孔隙流體壓力增加導(dǎo)致的巖石破裂作用,包括了孔隙流體對(duì)完整巖石和原有裂縫的水力作用[2,36],且水力破裂是由孔隙流體導(dǎo)致的一種張性破裂[36-37]。根據(jù)有效差應(yīng)力與巖石抗張強(qiáng)度之比的不同,將巖石脆性破裂分為:張性破裂、張性剪切破裂和剪切破裂,不同類型的脆性破裂需要不同的孔隙壓力[37]。

    圖9 松遼盆地泥質(zhì)巖蓋層孔喉半徑分布[28]Fig.9 Pore throat radius distribution of argillaceous caprock in the Songliao Basin[28]

    圖10 泥巖封閉能力動(dòng)態(tài)變化過(guò)程定量評(píng)價(jià)方法[14]Fig.10 Quantitative evaluation method for the change of caprock sealing properties [14]a.泥巖毛細(xì)管力和孔隙度的關(guān)系;b.泥巖毛細(xì)管力和滲透率的關(guān)系;c.泥巖毛細(xì)管力隨地質(zhì)歷史時(shí)間的變化關(guān)系

    在孔隙流體壓力增加的過(guò)程中,有效應(yīng)力降低,有效差應(yīng)力不變,表現(xiàn)為應(yīng)力圓向靠近破裂包絡(luò)線方向移動(dòng),位移量等于孔隙流體壓力的增加量(Δp=p2-p1),但應(yīng)力圓的大小不變。當(dāng)應(yīng)力圓與破裂包絡(luò)線相切于(-T,0)時(shí),發(fā)生水力破裂,形成與最小主應(yīng)力σ3方向垂直的水力裂縫。因此,水力破裂準(zhǔn)則表示為:

    (1)

    也就是說(shuō),當(dāng)孔隙流體壓力大于最小主應(yīng)力與巖石抗張強(qiáng)度之和,即當(dāng)最小有效主應(yīng)力表現(xiàn)為張力且大于巖石抗張強(qiáng)度時(shí),發(fā)生水力破裂[9,12,37-38]。對(duì)墨西哥灣超壓泥巖和挪威Snorre油田研究發(fā)現(xiàn),當(dāng)孔隙流體壓力大于最小主應(yīng)力,即接近靜巖壓力的85%和82%時(shí),泥巖發(fā)生水力破裂,形成垂直的水力裂縫[39-40]。因此,在對(duì)水力破裂機(jī)理研究的基礎(chǔ)上,提出了評(píng)價(jià)水力破裂風(fēng)險(xiǎn)的定量方法[41-42]。Gaarenstroom定義最小水平主應(yīng)力σh與孔隙流體壓力p之差為保持力[41],孔隙流體壓力越大則保持力就越小,相應(yīng)水力破裂的風(fēng)險(xiǎn)就越大。研究發(fā)現(xiàn)當(dāng)保持力低于7 MPa時(shí),水力破裂的風(fēng)險(xiǎn)就會(huì)明顯增加。

    Leach研究表明,油氣主要分布在超壓帶頂界及稍下部位,頂部之下勘探成功率逐漸減小[43]。隨著埋藏深度增加,盡管地層壓力明顯增大,但超壓在深層并不容易引起蓋層破裂,淺層很小超壓即可引起蓋層破裂,超壓引起的蓋層破裂可能不是超壓盆地深層油氣勘探潛力的主要制約因素。按著超壓演化、地層破裂壓力和油氣充注時(shí)間關(guān)系,可將超壓圈閉按含油氣性分為四種情況(圖11):充滿型圈閉、欠充滿圈閉、散失型圈閉和未充注圈閉[11]。

    3 蓋層脆-韌性轉(zhuǎn)化定量表征

    無(wú)論是泥質(zhì)巖還是膏鹽巖蓋層,隨著埋藏深度增加,成巖程度、物性及溫壓環(huán)境發(fā)生改變,蓋層均發(fā)生力學(xué)性質(zhì)變化,變形歷經(jīng)3個(gè)過(guò)程:脆性、脆-韌性和韌性(圖12a)。不同變形階段蓋層封閉能力及破裂方式不同,蓋層頂部封閉能力存在差異[13]。

    3.1 蓋層脆-韌性變形影響因素

    膏鹽巖脆韌性變形的主要影響因素是圍壓[44-50]和溫度[49],其次是加載速率[44,50-51]、流體壓力[46,51]、組成和內(nèi)部結(jié)構(gòu)[45]。保持單一變量原則,隨圍壓增加膏鹽巖向脆-韌性域甚至向韌性域轉(zhuǎn)變;同理,溫度越高膏鹽巖韌性變形特征越明顯。實(shí)驗(yàn)證實(shí)加載速率足夠大(大于5×10-9/s)時(shí)鹽巖就會(huì)發(fā)生脆性變形[52],但在大多數(shù)鹽構(gòu)造背景下應(yīng)變速率很少大于10-12/s,只有在孔隙流體超壓條件下,鹽巖才會(huì)發(fā)生脆性破裂[51],否則,孔隙流體對(duì)膏鹽巖的脆韌性影響可以忽略[46]。由于鹽巖的可塑性大于膏巖,因此發(fā)生脆韌性變形的難易程度為:鹽巖>含鹽夾層的膏巖>純膏巖。此外,粗粒硬石膏樣品比細(xì)粒樣品更易發(fā)生脆韌性變形,平行于面理方向壓縮比垂直于面理方向壓縮更易發(fā)生脆韌性變形[45]。

    圖11 超壓圈閉類型模式[11]Fig.11 Schematic diagrams showing the types of overpressured reservoirs[11].a.充滿型圈閉;b.欠充滿圈閉;c.散失型圈閉;d.未充注圈閉

    圖12 巖石脆韌性轉(zhuǎn)化應(yīng)力-應(yīng)變曲線及不同圍壓下實(shí)驗(yàn)變形曲線特征Fig.12 Stress-strain curves for the rock brittle-ductile transition and the variation of the curves under varying confining pressures in experimentsa.不同脆韌性巖石應(yīng)力-應(yīng)變曲線模式[55];b,c.啟莫里階泥巖和含油頁(yè)巖應(yīng)力-應(yīng)變曲線(數(shù)值為有效垂直固結(jié)應(yīng)力)[12];d,e,f.庫(kù)車坳陷純膏 巖、含鹽泥巖和含泥鹽巖應(yīng)力-應(yīng)變曲線(數(shù)值為圍壓)

    對(duì)于泥頁(yè)巖而言,除以上因素外,礦物成分[53]、成巖作用[9-10]和流體性質(zhì)[12]也影響脆-韌性轉(zhuǎn)化。泥頁(yè)巖中石英為典型的脆性礦物,其含量越高,泥頁(yè)巖脆性程度越高,Jarvie等人提出了利用巖石礦物分析法判斷巖石脆性程度[54]。頁(yè)巖在未固結(jié)-半固結(jié)成巖階段密度很低,在大多數(shù)沉積盆地的圍壓范圍內(nèi)只發(fā)生韌性變形而不會(huì)發(fā)生脆性破裂。如果頁(yè)巖的密度大于2.1 g/cm3,在足夠的應(yīng)變下將發(fā)生脆性破裂,為此,提出了利用密度定量判斷泥頁(yè)巖脆韌性變形定量評(píng)價(jià)方法[6,9-10]。對(duì)比研究表明,含油泥巖比純泥巖更容易發(fā)生韌性變形(圖12b,c)[12-13],在15 MPa圍壓條件下即發(fā)生韌性變形。

    3.2 蓋層脆-韌性轉(zhuǎn)化階段表征

    3.2.1 膏鹽巖脆-韌性定量表征

    影響膏鹽巖脆-韌性的主要因素是圍壓和溫度,依據(jù)三軸壓縮實(shí)驗(yàn)可以獲得不同圍壓和溫度條件下的峰值強(qiáng)度和殘余強(qiáng)度[55](圖12d,f)。巖石力學(xué)中,脆性破裂的表征方法有庫(kù)侖破裂準(zhǔn)則、格里菲斯準(zhǔn)則、修正格里菲斯準(zhǔn)則以及摩爾-庫(kù)侖破裂準(zhǔn)則。摩爾-庫(kù)侖破裂包絡(luò)線一般為二次曲線[56],為巖石不同溫度壓應(yīng)力-應(yīng)變曲線上峰值強(qiáng)度(剪切破裂強(qiáng)度)的二次擬合曲線。Byerlee摩擦定律是一個(gè)與巖石類型和滑動(dòng)面粗糙度等完全無(wú)關(guān),適用于自然界各種摩擦滑動(dòng)現(xiàn)象[57]。在差應(yīng)力(σ1-σ3)和圍壓σ3坐標(biāo)中,低圍壓下巖石體現(xiàn)出純脆性破裂,峰值破裂曲線是直線,斜率與破裂后的滑動(dòng)摩擦曲線相近。隨著圍壓增大,巖石中部分礦物體現(xiàn)出塑性變形,巖石的內(nèi)摩擦角及滑動(dòng)摩擦角都減??;當(dāng)巖石破裂強(qiáng)度與圍壓恰好滿足Byerlee定律時(shí),即摩爾-庫(kù)侖破裂包絡(luò)線與Byerlee摩擦滑動(dòng)曲線相交時(shí),巖石開(kāi)始從脆性向脆-韌性轉(zhuǎn)變。因此,Byerlee摩擦滑動(dòng)定律標(biāo)志著脆性破裂的結(jié)束(圖13)。應(yīng)力降是指巖石破裂峰值強(qiáng)度與殘余強(qiáng)度的差值[58],微觀上脆性和脆韌性變形巖石普遍存在微破裂和聲發(fā)射,從應(yīng)力-應(yīng)變曲線上看具有明顯的應(yīng)力降。隨著圍壓的增加,巖石塑性成分逐漸增多,應(yīng)力降逐漸減小,當(dāng)圍壓增加到某一臨界值時(shí),應(yīng)力降為零,即巖石殘余強(qiáng)度破裂包絡(luò)線與摩爾庫(kù)侖破裂包絡(luò)線相交時(shí),巖石開(kāi)始轉(zhuǎn)變?yōu)轫g性變形,不發(fā)生脆性破裂。Goetze研究表明,當(dāng)應(yīng)力降為零時(shí),大部分?jǐn)?shù)據(jù)表明所施加圍壓約與破裂強(qiáng)度相近時(shí),標(biāo)志著半脆性向塑流過(guò)渡的轉(zhuǎn)變[58]。

    圖13 厘定巖石脆韌性變形階段模式(a)、脆韌性轉(zhuǎn)化臨界圍壓(b)及庫(kù)車坳陷庫(kù)姆格列木組膏巖應(yīng)力-應(yīng)變曲線(c)Fig.13 Schematic diagram showing the brittle-ductile deformation stage of rock,the stress-strain curves and brittle-ductile transition critical confining pressure of Kumugeliemu Formation gypsolyte in the Kuqa Depression

    利用此方法,付曉飛厘定了庫(kù)車膏巖脆-韌性轉(zhuǎn)化圍壓[13]。結(jié)果表明,膏巖脆韌性轉(zhuǎn)換臨界圍壓明顯低于含鹽泥巖,即相同圍壓條件下,泥巖表現(xiàn)為脆性時(shí),膏巖可能已經(jīng)轉(zhuǎn)變?yōu)榇囗g性甚至韌性。因此,以膏巖脆韌性轉(zhuǎn)換臨界圍壓作為庫(kù)姆格列木組區(qū)域蓋層的脆韌性轉(zhuǎn)換的臨界條件,分別為46 MPa和90 MPa。

    3.2.2 泥頁(yè)巖脆-韌性定量表征

    隨著埋藏深度增加,泥巖經(jīng)歷更復(fù)雜的脆-韌性轉(zhuǎn)化過(guò)程。在大多數(shù)沉積盆地的圍壓范圍內(nèi),頁(yè)巖的密度在大致小于2.1 g/cm3的情況下,只發(fā)生韌性變形而不會(huì)發(fā)生脆性破裂[9,13]。如果頁(yè)巖的密度大于2.1 g/cm3,在足夠的應(yīng)變下將發(fā)生脆性破裂。Corcoran和Dore利用密度和破裂時(shí)應(yīng)變定量判斷泥頁(yè)巖脆-韌性轉(zhuǎn)化過(guò)程[10](圖14)。脆性階段,密度大于2.5 g/cm3,破裂前應(yīng)變小于3%。過(guò)渡階段,密度介于2.5~2.25 g/cm3,應(yīng)變介于5%~8%。韌性階段,密度小于2.25 g/cm3,應(yīng)變大于8%。但Nyg?rd等人研究表明,固結(jié)泥巖在15 MPa時(shí)開(kāi)始從脆性變形向韌性變形轉(zhuǎn)化[12]。綜合考慮這些因素,泥頁(yè)巖脆-韌性轉(zhuǎn)化可能分為多個(gè)階段。淺埋未固結(jié)-半固結(jié)成巖階段,密度普遍小于2.1 g/cm3,主要發(fā)生韌性變形。固結(jié)成巖階段,受圍壓和溫度影響經(jīng)歷脆、過(guò)渡和韌性轉(zhuǎn)化(圖12b)[9-10];化學(xué)成巖作用導(dǎo)致泥頁(yè)巖在高圍壓條件下,不經(jīng)歷韌性變形而直接發(fā)生脆性破裂;抬升過(guò)程,泥頁(yè)巖為超固結(jié)成巖階段,伴隨壓力釋放和圍壓降低,泥巖變脆,易于發(fā)生破裂。Ingram 1997年提出利用超固結(jié)比定量判斷破裂發(fā)生的條件,認(rèn)為OCR大于2.5時(shí)蓋層破裂導(dǎo)致滲漏的風(fēng)險(xiǎn)增加[8,12,14,59]。

    圖14 利用泥巖密度和破裂時(shí)應(yīng)變判斷脆韌性[10]Fig.14 Ductility and brittleness distinguished by means of mudstone density(g/cm3)and fracturing strain[10]

    4 蓋層完整性評(píng)價(jià)

    蓋層完整性破壞通常由潤(rùn)濕性改變和破裂作用形成微滲漏空間引起的[6],分為五種類型(圖15)。

    4.1 潤(rùn)濕性改變

    蓋層內(nèi)孔隙規(guī)模流體流動(dòng)實(shí)驗(yàn)表明[60],垂向運(yùn)移的殘留水進(jìn)入油層和蓋層的接觸面,油中酸性化合物從油中分離出進(jìn)入殘留水中,水?dāng)y帶酸性化合物進(jìn)入蓋層中,這種流動(dòng)發(fā)生在不連續(xù)、狹窄的通道中,酸性極性化合物形成有機(jī)-金屬化合物,貼附在礦物表面,導(dǎo)致水膜破裂,有機(jī)-金屬化合物導(dǎo)致局部潤(rùn)濕性改變,變?yōu)橛蜐?rùn)濕,減小毛細(xì)管壓力,促使局部油和氣進(jìn)入蓋層。

    4.2 斷裂破壞作用

    斷裂在不同脆-韌性蓋層中變形機(jī)制不同,封閉能力評(píng)價(jià)方法也不同。脆性階段,脆性域的蓋層由于破裂作用易被斷層破壞,隨著斷距的增大,裂縫的密度逐漸增大,直至互相連通后油氣就會(huì)通過(guò)斷層發(fā)生垂向運(yùn)移[61]。野外觀察證實(shí)脆性的膏鹽巖一般會(huì)形成貫通性的大斷層,其斷裂帶通常被斷層角礫巖和軟的斷層泥充填,如意大利亞平寧山脈北部三疊系蒸發(fā)巖內(nèi)的脆性斷層[62]。基于裂縫垂向連通程度受控于蓋層厚度和斷距,提出斷接厚度來(lái)定量表征裂縫垂向連通性,即平行于斷面的蓋層厚度與斷層位移的差值[63]。該值越大,裂縫垂向?qū)芰υ讲?。處于脆性域的蓋層存在臨界的斷接厚度值,該值越大封閉能力越強(qiáng)。

    脆-韌性階段,蓋巖發(fā)育典型的涂抹結(jié)構(gòu),只要涂抹保持連續(xù)性,斷層垂向就是封閉的。通過(guò)泥巖涂抹系數(shù)SSF(Shale smear factor),即斷距與泥巖厚度的比率[64],可以預(yù)測(cè)涂抹的發(fā)育程度,實(shí)現(xiàn)斷層封閉性定量評(píng)價(jià)。多數(shù)學(xué)者認(rèn)為[64-66],泥巖涂抹的連續(xù)性受控于斷距與泥巖厚度的比率,小規(guī)模斷層(斷距小于15 m)[66]泥巖涂抹連續(xù)SSF值變化范圍為1~50,通常泥巖層厚度為0~10 cm,斷距為3 cm至幾分米[64]。對(duì)于規(guī)模較大的斷層(斷距大于15m),泥巖涂抹保持連續(xù)性的臨界值較小,一般為4~8[65-67],但也有SSF值為6~20的斷層保持連續(xù)的泥巖涂抹[67]。Takahashi通過(guò)高溫高壓物理模擬表明[68],有效正應(yīng)力為30 MPa,當(dāng)SSF大于4.9時(shí),粉砂巖形成的涂抹失去連續(xù)性;有效應(yīng)力提高到40 MPa,涂抹保持連續(xù)性的臨界SSF值為6.6。因此,泥巖隨著埋深增加,泥巖涂抹越發(fā)育,且容易保持連續(xù)性。

    圖15 蓋層完整性影響因素模式Fig.15 Schematic diagram of factors affecting the caprock integrity

    韌性階段,鹽巖具有流動(dòng)特征,伴隨斷裂逆沖滑動(dòng),鹽巖沿著斷面流動(dòng),并在逆沖帶前鋒擠出,此階段的斷層垂向封閉。但是在差異壓實(shí)作用下,韌性的鹽巖向構(gòu)造高部位(低應(yīng)力區(qū))流動(dòng)和局部集中,使得上覆地層發(fā)生隆起,形成鹽構(gòu)造;同時(shí)在相對(duì)薄的或厚度變化劇烈的鹽巖區(qū),鹽巖的蠕動(dòng)作用會(huì)形成一些塑性斷層斷穿鹽層,形成鹽巖缺失區(qū),即“鹽窗”[69]。鹽窗的形成為在鹽巖區(qū)域性蓋層之下的油氣向鹽上及鹽間運(yùn)移調(diào)整提供了通道。從巴西坎波斯盆地和濱里海盆地Sagizski區(qū)塊的剖面圖中可以看出,鹽上油氣藏的分布與鹽窗的發(fā)育位置具有極好的相關(guān)性[69-70]。

    4.3 亞地震斷層與砂體連通

    通過(guò)二維和三維地震解釋,可以清楚地識(shí)別出大尺度斷裂作用,小尺度斷裂作用也能夠通過(guò)空間上孤立的一維井資料來(lái)進(jìn)行識(shí)別。但是對(duì)于中尺度的斷裂作用(位移大約在幾分米到30 m之間),常稱為亞地震斷層,通常既不能從地震數(shù)據(jù)上識(shí)別,也不能從井資料上識(shí)別[71]。然而,這種中尺度的斷裂作用是控制油氣沿蓋層滲漏的主要因素[72]。蓋層內(nèi)亞地震斷裂與砂體連通可形成有效的滲漏通道[8],數(shù)值模擬證實(shí)砂泥互層蓋層垂向封閉性取決于裂縫數(shù)量和泥巖層層數(shù)?;跀嗔焉L(zhǎng)冪率規(guī)律和斷層長(zhǎng)度與斷距關(guān)系,結(jié)合最大單層泥巖厚度分布,可以定量厘定蓋層垂向滲漏最小斷距和斷層頻數(shù)。Ingram提出當(dāng)斷層條數(shù)(斷距大于最大單層泥巖厚度)是泥巖層層數(shù)5倍時(shí),蓋層垂向滲漏風(fēng)險(xiǎn)極大[8]。

    4.4 構(gòu)造裂縫連通

    泥頁(yè)巖蓋層常含有大量的微裂縫,其在正常應(yīng)力狀態(tài)下往往是孤立與閉合的,彼此不連通,蓋層可以很好地封閉油氣。當(dāng)應(yīng)力狀態(tài)發(fā)生改變時(shí),就會(huì)形成新的微裂縫或閉合裂縫重新張開(kāi),從而使微裂縫相互連通形成微滲漏空間。引起這種應(yīng)力發(fā)生變化的機(jī)制主要有構(gòu)造擠壓、構(gòu)造抬升以及異常高壓流體等[73]。例如,在層狀泥頁(yè)巖中發(fā)育有大量的與層面高角度相交或近垂直的構(gòu)造裂縫,裂縫的形成與分布受巖層的控制,裂縫通常在巖層內(nèi)發(fā)育,并終止層面上,很少穿越巖層界面,裂縫高度一般等于裂隙化巖層的厚度[74-75]。在一定的巖層厚度范圍內(nèi),裂縫的平均間距與裂隙化巖層的平均厚度呈較好的線性關(guān)系,其比值稱作裂縫間距比率(FSR),裂縫間距比率與巖石的應(yīng)變呈正比。根據(jù)大量野外露頭實(shí)測(cè),當(dāng)FSR<1.17時(shí),在巖層中只發(fā)育這種層控裂縫,它們是孤立的、不連通的。當(dāng)FSR>1.17時(shí),開(kāi)始出現(xiàn)貫穿裂縫,且貫穿裂縫的頻率與FSR呈正比,說(shuō)明應(yīng)變的增加導(dǎo)致裂縫發(fā)生連通,從而形成滲漏通道。構(gòu)造抬升作用沒(méi)有改變最大主應(yīng)力,但是使最小主應(yīng)力降低,造成應(yīng)力莫爾圓變大與包絡(luò)線相交,致使巖石中容易產(chǎn)生破裂,使裂縫發(fā)生連通,從而形成滲漏空間。

    4.5 水力破裂

    當(dāng)儲(chǔ)層孔隙流體壓力超過(guò)圍巖的最小圍壓時(shí),蓋層破裂形成裂縫或原有裂縫張開(kāi),構(gòu)成了由高壓沉積物到上覆地層的高滲透性通道[38]。大量的流體通過(guò)裂縫進(jìn)行運(yùn)移,發(fā)生流體排放活動(dòng),隨著流體經(jīng)裂縫流出,孔隙流體壓力降低,引起礦物膠結(jié)沉淀,裂縫的寬度減小,滲透率降低,最終裂縫閉合,流體排放活動(dòng)終止[8-9]。一旦裂縫閉合,蓋層會(huì)一直保持非滲透性,直到深部運(yùn)移來(lái)的流體使孔隙流體壓力增加至破裂為止,新的排放活動(dòng)開(kāi)始從而周期性地循環(huán)。對(duì)于相對(duì)堅(jiān)硬的蓋層,裂縫會(huì)在20~50年內(nèi)保持張開(kāi)狀態(tài)[76]。Caillet以挪威Snorre油田為例,研究發(fā)現(xiàn)蓋層中含有與儲(chǔ)層內(nèi)油氣同源的油氣。根據(jù)漏失實(shí)驗(yàn)數(shù)據(jù),當(dāng)儲(chǔ)層頂部的孔隙流體壓力達(dá)到靜巖壓力的82%,近似等于最小水平應(yīng)力時(shí),發(fā)生水力破裂導(dǎo)致泥巖的滲透性增加,油氣從儲(chǔ)層中逃逸出來(lái)。

    5 蓋層研究存在的主要問(wèn)題

    盡管蓋層封閉性研究取得了長(zhǎng)足進(jìn)展,但仍有許多問(wèn)題需要深入探索,主要包括以下幾方面。

    不同級(jí)別蓋層與油氣保存。對(duì)鄂爾多斯盆地深部油氣成藏研究表明,石炭系-二疊系煤系地層為區(qū)域性蓋層,上組合直接蓋層為本溪組鋁土質(zhì)泥巖,下組合直接蓋層為石炭系煤系地層,儲(chǔ)層內(nèi)部還發(fā)育膏鹽巖和泥云巖隔層,三級(jí)蓋層對(duì)油氣保存均起作用;研究表明,對(duì)裂縫改造的致密儲(chǔ)層斷層油氣藏而言,油氣富集受控于區(qū)域性蓋層,油氣緊鄰區(qū)域性蓋層分布[20];對(duì)常規(guī)儲(chǔ)層斷層油氣藏而言,油氣富集受控于直接蓋層,區(qū)域性蓋層起到油氣充注的頂部封閉作用;厘定三級(jí)蓋層分布,結(jié)合典型油氣藏解剖,剖析不同級(jí)別蓋層對(duì)油氣保存作用,對(duì)油氣勘探才能起到指導(dǎo)作用。

    影響蓋層封閉能力的宏觀地質(zhì)因素,蓋層厚度、巖性組合和蓋地比對(duì)封閉能力的影響缺少數(shù)據(jù)支撐;按著毛細(xì)管封閉機(jī)理,蓋層厚度并非是封閉的主控因素[2],但一些盆地依然顯示油氣藏分布在厚度較大蓋層范圍[19-20],可能與蓋層巖性組合有關(guān)。依據(jù)巖性組合可將泥質(zhì)巖蓋層分為3類:一是質(zhì)純泥巖,如松遼盆地青一段和嫩一二段泥巖蓋層,單層厚度可達(dá)百米,含砂量不超過(guò)10%[28];二是砂泥互層泥巖,如柴達(dá)木盆地區(qū)域性蓋層[20];三是致密砂礫巖,如松遼盆地徐家圍子斷陷火山巖氣藏頂部營(yíng)城組四段砂礫巖蓋層。蓋地比是蓋層質(zhì)量評(píng)價(jià)的重要指標(biāo)[6],應(yīng)用砂體連通性原理和砂體-亞地震斷層連通性去定量評(píng)價(jià)蓋層品質(zhì)是未來(lái)蓋層研究的重要領(lǐng)域之一[9,77]。

    泥頁(yè)巖巖石物理特征及成巖演化階段,泥頁(yè)巖蓋層礦物組成、孔隙結(jié)構(gòu)、物性特征和力學(xué)性質(zhì)等因素決定的巖石物理特征及成巖演化階段缺少有效數(shù)據(jù)。受巖心資料和實(shí)驗(yàn)手段限制,對(duì)泥頁(yè)巖蓋層測(cè)試難度大,缺少統(tǒng)一的行業(yè)標(biāo)準(zhǔn),極大限制了對(duì)蓋層封閉能力動(dòng)態(tài)變化規(guī)律的認(rèn)識(shí)[14]。

    蓋層封閉能力動(dòng)態(tài)評(píng)價(jià),盡管Jin等人建立了蓋層毛細(xì)管封閉能力動(dòng)態(tài)評(píng)價(jià)方法[14],但埋藏過(guò)程中泥質(zhì)巖蓋層的孔隙度在一定范圍內(nèi)變化不大,能否與毛細(xì)管壓力是正相關(guān)關(guān)系,還需要利用不同埋藏深度蓋層樣品去驗(yàn)證。三軸壓縮實(shí)驗(yàn)證實(shí)抬升過(guò)程中滲透率增大[78],但不是數(shù)量級(jí)變化,導(dǎo)致抬升過(guò)程中滲透率變化的核心地質(zhì)因素還是裂縫。如何定量評(píng)價(jià)需要通過(guò)實(shí)驗(yàn)驗(yàn)證,僅利用OCR大于2.5去判斷仍缺乏理論依據(jù)[12]。Hao等人建立了水力封閉動(dòng)態(tài)演化模型[11],但并沒(méi)得到油田驗(yàn)證,因此需要開(kāi)展巖石破裂條件和古壓力恢復(fù)研究,才能建立能應(yīng)用到油田的地質(zhì)模型。

    蓋層脆-韌性轉(zhuǎn)化過(guò)程及破裂條件定量表征,埋藏過(guò)程中,膏鹽巖脆-韌性變形轉(zhuǎn)化分3個(gè)階段,即脆性階段、脆-韌性階段和韌性階段。抬升階段與埋藏過(guò)程呈現(xiàn)對(duì)稱式反方向變化[13]。但泥頁(yè)巖脆-韌性轉(zhuǎn)化影響因素很多,當(dāng)密度低于2.1 g/cm3時(shí)發(fā)生韌性變形[6,9-10],固結(jié)成巖后受圍壓和溫度影響較大,也要經(jīng)歷脆、脆-韌性和韌性變形[12];特定礦物組成和化學(xué)成巖作用也可能導(dǎo)致韌性變形階段缺失。抬升過(guò)程,伴隨壓力釋放,超固結(jié)的泥巖會(huì)變脆。影響泥頁(yè)巖脆-韌性變形的因素多而復(fù)雜,加之泥頁(yè)巖樣品難以制備,目前研究相對(duì)薄弱。

    深層斷層對(duì)油氣的輸導(dǎo)及保存,經(jīng)歷半個(gè)多世紀(jì)的研究,對(duì)斷層輸導(dǎo)和保存作用有了一定認(rèn)識(shí)[79,80],但針對(duì)蓋層內(nèi)斷裂帶結(jié)構(gòu)研究還缺乏系統(tǒng)認(rèn)識(shí)[80],特別是斷裂在不同脆-韌性蓋層內(nèi)變形機(jī)制、斷裂帶內(nèi)部結(jié)構(gòu)和封閉性評(píng)價(jià)方法還不全面[81]。付曉飛等人深入研究了不同脆-韌性膏鹽巖中斷裂帶結(jié)構(gòu)及封閉能力評(píng)價(jià)方法[13],Holland解剖了超固結(jié)泥巖斷裂帶內(nèi)部結(jié)構(gòu)及封閉能力演化規(guī)律,對(duì)脆-韌性泥巖涂抹發(fā)育規(guī)律開(kāi)展了系統(tǒng)研究[64-65,67-68,81],但均缺少合理的定量評(píng)價(jià)方法,對(duì)斷裂輸導(dǎo)和保存作用動(dòng)態(tài)演化過(guò)程研究仍較薄弱。

    蓋層完整性定量評(píng)價(jià),斷裂、構(gòu)造裂縫和水力裂縫是導(dǎo)致蓋層完整性破壞的主因,目前僅從機(jī)理探討了對(duì)蓋層完整性影響[8-9],缺少能夠應(yīng)用到地下的定量評(píng)價(jià)方法。

    6 蓋層研究趨勢(shì)

    6.1 斷裂和蓋層耦合控藏機(jī)理和控藏模式

    斷裂和蓋層是油氣保存條件的關(guān)鍵因素[5]。由于斷裂對(duì)儲(chǔ)層具有反向改造作用[82],致密儲(chǔ)層中斷裂帶為高滲透性油氣運(yùn)移通道[83],油氣聚集側(cè)向封閉靠巖性對(duì)接;對(duì)盤(pán)若為區(qū)域性蓋層才能形成規(guī)模性聚集,因此斷裂和區(qū)域性蓋層是致密儲(chǔ)層油氣保存關(guān)鍵因素[82]。高孔性儲(chǔ)層中斷裂帶為低滲透性的[79,84-87],油氣聚集側(cè)向封閉靠斷層巖[85,88],因此斷裂和直接蓋層、局部蓋層、區(qū)域性蓋層均能有效封閉油氣。伴隨油氣精細(xì)勘探,在大尺度沉積微相研究的基礎(chǔ)上,深入開(kāi)展油氣藏解剖,進(jìn)而明確直接蓋層、局部性蓋層、區(qū)域性蓋層和斷裂耦合配置及控藏機(jī)理和控藏模式。

    6.2 不同類型蓋巖成巖演化及封閉能力動(dòng)態(tài)演化過(guò)程

    成巖演化階段劃分是蓋層巖石封閉能力動(dòng)態(tài)演化的基礎(chǔ)。常規(guī)油氣勘探階段針對(duì)中淺層蓋層取心相對(duì)較多。伴隨非常規(guī)油氣勘探深入,深層烴源巖成為重要的勘探目的層,取心越來(lái)越多。具備測(cè)試不同埋藏深度和不同抬升幅度的蓋巖物性、孔隙結(jié)構(gòu)、比表面積、排替壓力、地球化學(xué)和巖石力學(xué)等參數(shù)條件。依據(jù)各種參數(shù)對(duì)埋藏和抬升的響應(yīng)[14],考慮將成巖階段在傳統(tǒng)的埋藏物理壓實(shí)、化學(xué)成巖、變質(zhì)階段和抬升階段基礎(chǔ)上進(jìn)一步細(xì)分,結(jié)合斷裂和裂縫形成演化及對(duì)蓋層破壞,深入研究蓋巖封閉能力動(dòng)態(tài)演化過(guò)程。

    6.3 蓋層脆-韌性轉(zhuǎn)化過(guò)程定量表征及斷裂和裂縫形成演化

    在埋藏和抬升過(guò)程中蓋層巖石力學(xué)發(fā)生明顯變化[13],但在什么條件下蓋層發(fā)生脆-韌性轉(zhuǎn)化依然是個(gè)難題,需要通過(guò)三軸壓縮試驗(yàn)、聲發(fā)射特征、微觀變形機(jī)理和蓋層內(nèi)構(gòu)造變形特征研究來(lái)建立蓋層脆-韌性轉(zhuǎn)化過(guò)程定量表征方法?;诼额^區(qū)斷裂帶內(nèi)部結(jié)構(gòu)解剖,分析斷裂和裂縫在不同巖石力學(xué)特征蓋層段變形機(jī)制及形成演化規(guī)律,才能搞清斷裂和裂縫對(duì)蓋層破壞作用。

    6.4 蓋層完整性控制因素及定量評(píng)價(jià)

    由界面張力、潤(rùn)濕性和儲(chǔ)層超壓引起蓋層毛細(xì)管封閉能力降低導(dǎo)致油氣散失并沒(méi)得到有效證實(shí),而由于構(gòu)造變動(dòng)產(chǎn)生的斷裂和裂縫引起油氣藏破壞實(shí)例很多[89],蓋層完整性評(píng)價(jià)成為保存條件研究的核心[75]。影響蓋層完整性的因素包括斷裂、構(gòu)造裂縫、水力裂縫、亞地震斷層和砂體的連通?;趲r層能干性差異形成的力學(xué)界面對(duì)裂縫和斷裂擴(kuò)展影響和斷裂源于裂縫遞進(jìn)發(fā)展原理,合理解釋斷裂和有效預(yù)測(cè)亞地震及裂縫分布;基于原地應(yīng)力場(chǎng)特征和破裂機(jī)理,有效預(yù)測(cè)斷裂和裂縫啟閉性[90],成為解決蓋層完整性的有效手段。

    7 結(jié)論

    1) 蓋層封閉機(jī)理主要包括毛細(xì)管封閉和水力封閉機(jī)理。蓋層毛細(xì)管封閉能力主要受成巖程度、流體性質(zhì)、界面張力、潤(rùn)濕性和儲(chǔ)層超壓的影響。多期構(gòu)造演化的盆地蓋層封閉能力是變化的,埋藏過(guò)程中伴隨壓實(shí)成巖和化學(xué)成巖程度增強(qiáng),蓋層孔隙度逐漸降低,毛細(xì)管壓力逐漸增大;抬升過(guò)程由于裂縫的產(chǎn)生,滲透率明顯增加,毛細(xì)管壓力減小,利用孔隙度、滲透率與毛細(xì)管壓力之間的定量關(guān)系可以分別判斷蓋層在埋藏和抬升階段封閉能力變化。蓋層水力破裂主要受最小水平主應(yīng)力與孔隙流體壓力的控制,應(yīng)用保持力定量評(píng)價(jià)水力破裂風(fēng)險(xiǎn),保持力越小,風(fēng)險(xiǎn)越高。

    2) 蓋層脆韌性變形階段厘定是研究蓋層頂部封閉能力的基礎(chǔ)。對(duì)于膏鹽巖蓋層,應(yīng)用摩爾-庫(kù)侖準(zhǔn)則、拜爾利摩擦定律和應(yīng)力降規(guī)律定量判斷脆韌性轉(zhuǎn)換的臨界圍壓。而泥巖脆-韌性轉(zhuǎn)化過(guò)程更為復(fù)雜,在淺埋未固結(jié)-半固結(jié)成巖階段,密度普遍小于2.1 g/cm3,主要發(fā)生韌性變形;固結(jié)成巖階段,受圍壓和溫度影響經(jīng)歷脆性、脆韌性和韌性變形;化學(xué)成巖階段泥頁(yè)巖在高圍壓條件下不經(jīng)歷韌性變形而直接發(fā)生脆性破裂;抬升過(guò)程,泥頁(yè)巖為超固結(jié)成巖階段,泥巖變脆,易于發(fā)生破裂,目前普遍認(rèn)為OCR等于2.5是蓋層破裂滲漏的臨界。

    3) 影響蓋層完整性的因素主要包括5個(gè)方面。一是潤(rùn)濕性改變形成的微滲漏空間。二是斷裂的破壞,野外露頭證實(shí)斷裂在脆性蓋層中變形產(chǎn)生大量裂縫,利用斷接厚度判斷斷裂垂向封閉性;斷裂在脆-韌性蓋層中變形產(chǎn)生涂抹,利用SSF判斷斷裂垂向封閉性。三是亞地震斷層與砂體形成的滲漏空間,當(dāng)亞地震斷層條數(shù)為砂巖層數(shù)5倍時(shí)形成連續(xù)性滲漏通道。四是構(gòu)造裂縫連通形成的滲漏空間,當(dāng)FSR>1.17時(shí),開(kāi)始出現(xiàn)貫穿裂縫,F(xiàn)SR越大,裂縫連通程度越大。五是水力破裂形成的間歇性滲漏空間。

    4) 隨著勘探向非常規(guī)油氣藏和深層油氣藏?cái)U(kuò)展,蓋層對(duì)油氣的保存顯得更為重要,需要建立一套完善的適用于不同類型盆地、不同巖石類型蓋層的成巖演化階段和脆韌性階段劃分方法,明確不同成巖程度和脆韌性蓋層封閉能力評(píng)價(jià)參數(shù),從而開(kāi)展蓋層封閉能力動(dòng)態(tài)演化定量評(píng)價(jià)。此外,還需要開(kāi)展油氣藏精細(xì)解剖,明確不同級(jí)別蓋層和斷裂耦合配置及控藏機(jī)理和控藏模式;合理解釋斷裂、有效預(yù)測(cè)亞地震斷層和裂縫分布,基于原地應(yīng)力場(chǎng)和破裂機(jī)理準(zhǔn)確恢復(fù)古應(yīng)力,有效預(yù)測(cè)斷裂和裂縫啟閉性,從而定量評(píng)價(jià)蓋層完整性。

    參 考 文 獻(xiàn)

    [1] Smith D A.Theoretical consideration of sealing and nonsealing faults[J].AAPG Bulletin,1966,50:363-374.

    [2] Watts N L.Theoretical aspects of cap-rock and fault seals for single-and two-phase hydrocarbon columns[J].Marine and Petroleum Geology,1987,4(4):274-307.

    [3] Downey M W.Evaluating seals for hydrocarbon accumulations[J].AAPG Bulletin,1984,68:1752-1763.

    [4] Grunau H R.A worldwide look at the cap-rock problem[J].Journal of Petroleum Geology,1987,10(3):245-266.

    [5] 呂延防,付廣,高大嶺,等.油氣藏封蓋研究[M].北京:石油工業(yè)出版社,1996:118-120.

    Lv Yanfang,Fu Guang,Gao Daling,et al.Study on the cap rock of reservoir[M].Beijing:Petroleum Industry Press,1996:118-120.

    [6] Skerlec G M.頂部封閉與斷層封閉性評(píng)價(jià)[C]∥劉德來(lái),薛良清.油氣圈閉勘探.北京:石油工業(yè)出版社,2002:395-444.

    Skerlec G M.Evaluation top and fault seal[C]∥Liu Delai,Xue Liangqing.Hydrocarbon traps exploration.Beijing:PetroleumIndustry Press,2002:395-444.

    [7] Jin Zhijun,Yuan Yusong,Liu Quanyou,et al.Controls of Late Jurassic-Early Cretaceous tectonic event on source rocks and seals in marine sequences,South China[J].Science China Earth Sciences,2013,56(2):228-239.

    [8] Ingram G M,Urai J L,Naylor M A.Sealing processes and top seal assessment[J].Norwegian Petroleum Society Special Publications,1997,7:165-174.

    [9] Ingram G M,Urai J L.Top-seal leakage through faults and fractures:the role of mudrock properties[J].Geological Society,1999,158(1):125-135.

    [10] Corcoran D V,Dore A G.Top seal assessment in exhumed basin settings-some insights from Atlantic Margin and borderland basins[C]∥Koestler A G,Hunsdale R.Hydrocarbon Seal Quantification.NPF Special Publication,2002:89-107.

    [11] Hao F,Zhu W,Zou H,et al.Factors controlling petroleum accumulation and leakage in overpressured reservoirs[J].AAPG Bulletin,2015,99(5):831-925.

    [12] Nyg?rd R,Gutierrez M,Bratli R K,et al.Brittle-ductile transition,shear failure and leakage in shales and mudrocks[J].Marine and Petroleum Geology,2006,23(2):201-212.

    [13] 付曉飛,賈茹,王海學(xué),等.斷層-蓋層封閉性定量評(píng)價(jià)-以塔里木盆地庫(kù)車坳陷大北-克拉蘇構(gòu)造帶為例[J].石油勘探與開(kāi)發(fā),2015,42(3):300-309.

    Fu Xiaofei,Jia Ru,Wang Haixue,et al.Quantitative evaluation on the fault-caprock sealing capacity of Dabei-kelasu structural belt in Kuqa depression,Tarim basin,NW China[J].Petroleum Exploration and Deveopment,2015,42(3):300-309.

    [14] Jin Zhijun,Yuan Yusong,Sun Dongsheng,et al.Models for dynamic evaluation of mudstone/shale cap rocks and their applications in the lower Paleozoic sequences,Sichuan basin,SW China[J].Marine and Petroleum Geology,2014,49:121-128.

    [15] Baranova V,Mustaqeem A,Brouwer F.Application of seismic stratigraphy,multi-attribute analysis and neural networks to mitigate risk in new exploration frontiers-west newfoundland example[J].GeoCanada,2010.

    [16] 付廣,胡欣蕾.蓋層封閉天然氣有效性定量評(píng)價(jià)方法及應(yīng)用[J].山東科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,34(4):28-36.DOI:10.16452/j.cnki.sdkjzk.2015.04.004.

    Fu Guang,Hu Xinlei.Quantitative evaluation method of sealing gas efficirncy of caprock and its application[J]Journal of Shandong University of Science and Technology(Natural Science),2015,34(4):28-36.DOI:10.16452/j.cnki.sdkjzk.2015.04.004.

    [17] Ovcharenko A V,Ermakov B V,Myatchin K M,Shlezinger A E.Caprocks in hydrocarbon fields[J].Lithology and Mineral Resources,2007,42(2):180-191.

    [18] Zieglar D L.Hydrocarbon columns,buoyancy pressures,and seal efficiency:comparisons of oil and gas accumulations in California and the Rocky mountain area[J].AAPG Bulletin,1992,76(4):501-508.

    [19] 蔣有錄.油氣藏蓋層厚度與所封蓋烴柱高度關(guān)系問(wèn)題探討[J].天然氣工業(yè),1998,18(2):20-23.

    Jiang Youlu.A discussion on the relation between the cap rock thickness of oil & gas reservoir and its hydrocarbon column height[J].Natural Gas Industry,1998,18(2):20-23.

    [20] 付曉飛,劉小波,宋巖,等.中國(guó)中西部前陸沖斷帶蓋層品質(zhì)與油氣成藏[J].地質(zhì)論評(píng),2008,54(1):82-93.

    Fu Xiaofei,Liu Xiaobo,Song Yan,et al.Caprock quality and hydrocarbon accumulation in the basin of foreland thrust belts,central and western China[J].Geological Review,2008,54(1):82-93.

    [21] King P R.The connectivity and conductivity of overlapping sand bodies[C]∥North Sea Oil and Gas Reservoirs-Ⅱ.Springer Netherlands,1990:353-362.

    [22] Jackson M D,Yoshida S,Muggeridge A H,et al.Three-dimensional reservoir characterization and flow simulation of heterolithic tidal sandstones[J].AAPG Bulletin,2005,89(4):507-528.

    [23] Lei Yuhong,Luo Xiaorong,Zhang Likuan,et al.A quantitative method for characterizing transport capability of compound hydrocarbon carrier system[J].Journal of Earth Science,2013,24:328-342.

    [24] Dewhurst D N,Aplin A C,Sarda J P.Compaction driven evolution of porosity and permeability in natural mudstones:an experimental study[J].Journal of Geophysical Research,1998a,103(B1),651-661.

    [25] Dewhurst D N,Yang Y L,Aplin A C.Permeability and fluid flow in natural mudstones[C]∥In:Aplin A C,Fleet A J,Macquaker J H S.Muds and mudstones physical and fluid flow properties.Geological Society,London,Special Publications,1999,158,23-43.

    [26] Schl?mer S,Krooss B M.Experimental characterisation of the hydrocarbon sealing efficiency of cap rocks[J].Marine and Petroleum Geology,1997,14(5):565-580.

    [27] Mallon A J,Swarbrick R E.Diagenetic characteristics of low permeability,non-reservoir chalks from the Central North Sea[J].Marine and Petroleum Geology,2008,25(10):1097-1108.

    [28] 高瑞祺,蔡希源.松遼盆地油氣田形成條件與分布規(guī)律[J].北京:石油工業(yè)出版社,1997.

    Gao Ruiqi,Cai Xiyuan.Formation conditions and distribution rules for oil-gas fields in Songliao basin[J].Beijing:Petroleum Industry Press,1997.

    [29] Sayyouh M H,Hemeida A M,Al-Blehed M S,et al.Role of polar compounds in crude oils on rock wettability[J].Journal of Petroleum Science and Engineering,1991,6(3):225-233.

    [30] Bennett B,Buckman J O,Bowler B F J,LARTER S R.Wettability alteration in petroleum systems:the role of polar non-hydrocarbons.Petroleum Geoscience,2004,10:271-277.

    [31] Ruesl?tten H G,Hjelmeland O,Selle O M.Wettability of reservoir rocks and the influence of organo-metallic compounds.North Sea Oil and Gas Reservoirs-Ⅲ.Norwegian Institute of Technology(NTH),1994:317-324.

    [32] Melrose J C.Valid capillary pressure data at low wetting-phase saturations[J].SPE Reservoir Engineering,1990:95-99.

    [33] Rodgers S.Discussion:‘Physical constraints on hydrocarbon leakage and trapping revisited’by P.A.Bj?rkum-further aspects[J].Petroleum Geoscience,1999,5(4):421-423.

    [34] Bj?rkum P A,Walderhaug O,Nadeau P.Physical constraints on hydrocarbon leakage and trapping revisited[J].Petroleum Geoscience,1998,4:237-239.

    [35] 袁玉松,范明,劉偉新,等.蓋層封閉性研究中的幾個(gè)問(wèn)題[J].石油實(shí)驗(yàn)地質(zhì),2011,33(4):336-340.

    Yuan Yusogn,Fan Ming,Liu Weixin,et al.Several discussions of sealing capacity studies of caprock[J].Petroleum Geology & Experiment,2011,33(4):336-340.

    [36] Engelder T,Lacazette A.Natural hydraulic fracturing[J].Rock joints:Rotterdam,AA Balkema,1990:35-44.

    [37] Sibson R H.Structural permeability of fluid-driven fault-fracture meshes[J].Journal of Structural Geology,1996,18(8):1031-1042.

    [38] Lacazette A.Mechanics of natural hydraulic fracturing-review and introduction[C]∥AAPG Hedberg Research conference structural diagenesis:fundamental adrances and new application from a Holistie view of mechanical and chemical Processes,2004.

    [39] Anderson R N,He W,Hobart M A,et al.Active fluid flow in the Eugene island area,offshore Louisiana[J].The Leading Edge,1991,10(4):12-17.

    [40] Caillet G.The caprock of the Snorre field,Norway:a possible leakage by hydraulic fracturing[J].Marine and Petroleum Geology,1993,10(1):42-50.

    [41] Gaarenstroom L,Tromp R A J,Brandenburg A M.Overpressures in the central North Sea:implications for trap integrity and drilling safety[C]∥Geological Society,London,Petroleum Geology Conference series,Geological Society of London,1993,4:1305-1313.

    [42] Mildren S D,Hillis R R.FAST:a new approach to risking fault reactivation and related seal breach[J].AAPG Bulletin,2002:101-103.

    [43] Leach W G.Maximum hydrocarbon window determination in south Louisiana[J].Oil and Gas Journal,1993,91(13):81-84.

    [44] Fuenkajorn K,Sriapai T,Samsri P.Effects of loading rate on strength and deformability of Maha Sarakham salt[J].Engineering Geology,2012,135:10-23.

    [45] De Paola N,Faulkner D R,Collettini C.Brittle versus ductile deforma-tion as the main control on the transport properties of low-porosity anhydrite rocks[J].Journal of Geophysical Research:Solid Earth(1978-2012),2009:114(B6)258-266.

    [46] Hangx S J T,Spiers C J,Peach C J.Mechanical behavior of anhydrite caprock and implications for CO2sealing capacity[J].Journal of Geophysical Research:Solid Earth(1978-2012),2010:115,B07402.

    [47] Brantut N,Schubnel A,Guéguen Y.Damage and rupture dynamics at the brittle-ductile transition:The case of gypsum[J].Journal of Geophysical Research:Solid Earth(1978-2012),2011:116,B11404.

    [48] Popp T,Kern H,Schulze O.Evolution of dilatancy and permeability in rock salt during hydrostatic compaction and triaxial deformation[J].Journal of Geophysical Research:Solid Earth(1978-2012),2001,106(B3):4061-4078.

    [49] 高小平,楊春和,吳文,等.溫度效應(yīng)對(duì)鹽巖力學(xué)特性影響的試驗(yàn)研究[J].巖土力學(xué),2005,26(11):1775-1778.

    Gao Xiaoping,Yang Chunhe,Wu Wen,et al.Experimental studies on temperature effect of mechanical properties of rock salt[J].Rock and Soil Mechanics,2005,26(11):1775-1778.

    [50] Hamami M.Simultaneous effect of loading rate and confining pressure on the deviator evolution in rock salt[J].International Journal of Rock Mechanics and Mining Sciences,1999,36(6):827-831.

    [51] Davison I.Faulting and fluid flow through salt[J].Journal of the Geological Society,2009,166(2):205-216.

    [52] Thorel L,Ghoreychi M.Rock salt damage-experimental results and interpretation[J].Series on Rock and Soil Mechanics,1996,20:175-190.

    [53] Alqahtani A A.Effect of mineralogy and petrophysical characteristics on acoustic and mechanical properties of organic rich shale[C]∥Unconventional Resources Technology Conference(URTEC),2013:199-411.

    [54] Jarvie D M,Hill R J,Pollastro R M.Assessment of the gas potential and yields from shales:The Barnett shale model[J].Oklahoma Geological Survey Circular,2005,110:9-10.

    [55] Faulkner D R,Mitchell T M,Rutter E H,et al.On the structure and mechanical properties of large strike-slip faults[J].Geological Society,London,Special Publications,2008,299(1):139-150.

    [56] Kohlstedt D L,Evans B,Mackwell S J.Strength of the lithosphere:constraints imposed by laboratory experiments[J].Journal of Geophysical Research:Solid Earth(1978-2012),1995,100(B9):17587-17602.

    [57] Byerlee J D.Friction of rocks[J].Pure and Applied Geophysics,1978,116(4-5):615-626.

    [58] Goetze C.High temperature rheology of Westerly granite[J].Journal of Geophysics Research,1971,76:1223-1230.

    [59] Nyg?rd R,Gutierrez M,Gautam R,et al.Compaction behaviour of argillaceous sediments as function of diagenesis[J].Marine and Petroleum Geology,2004,21(3):349-362.

    [60] Teige G M G,Hermanrud C,Thomas W L H,et al.Capillary resistance and trapping of hydrocarbons:a laboratory experiment[J].Petroleum Geoscience,2005,11:125-129.

    [61] Bolton A,Maltman A.Fluid-flow pathways in actively deforming sediments:the role of pore fluid pressures and volume change[J].Marine and Petroleum Geology,1998,15(4):281-297.

    [62] De Paola N,Collettini C,Faulkner D R,et al.Fault zone architecture and deformation processes within evaporitic rocks in the upper crust[J].Tectonics,2008,27(4):1156-1178.

    [63] 呂延防,萬(wàn)軍,沙子萱,等.被斷裂破壞的蓋層封閉能力評(píng)價(jià)方法及其應(yīng)用[J].地質(zhì)科學(xué),2008,43(1):162-174.

    Lv Yanfang,Wan Jun,Sha Zixuan,et al.Evaluation method for seal ability of cap rock destructed by faulting and its application[J].Chinese Journal of Geology,2008,43(1):162-174.

    [64] Lindsay N G,Murphy F C,Walsh J J,Watterson J.Outcrop studies of shale smear on fault surface[J].International Association of Sedimentologists Special Publication,1993,15:113-123.

    [65] Yielding G,Freeman B,Needham D T.Quantitative fault seal prediction[J].AAPG Bulletin,1997,81(6):897-917.

    [66] Faerseth R B.Shale smear along large faults:continuity of smear and the fault seal capacity[J].Journal of the Geological Society,2006,163:741-751.

    [67] Doughty P T.Clay smear seals and fault sealing potential of an exhumed growth fault,Rio Grande rift,New Mexico[J].AAPG Bulletin,2003,3:427-444.

    [68] Takahashi M.Permeability change during experimental fault smearing[J].Journal of Geophysical Research,2003,108(B5):1-15.

    [69] 馬中振,謝寅符,張志偉,等.南美東緣鹽巖發(fā)育特征及其與油氣聚集的關(guān)系[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2013,43(2):360-370.

    Ma Zhongzhen,Xie Yinfu,Zhang Zhiwei,et al.Salt development characteristics and its controlling on hydrocarbon accumulation in eastern margin of south America[J].Journal of Jinlin University(Earth Science Edition),2013,43(2):360-370.

    [70] 古俊林,朱桂生,李永林.濱里海盆地Sagizski區(qū)塊鹽上層系成藏條件及分布規(guī)律研究[J].中國(guó)石油勘探,2012,17(2):57-61.

    Gu Junlin,Zhu Guisheng,Li Yonglin.Study on reservoir formation conditions and distribution in post-salt strata of Block Sagizski,Pre-Caspian Basin[J].China Petroleum Exploration,2012,17(2):57-61.

    [71] Lohr T,Krawczyk C M,Tanner D C,et al.Prediction of subseismic faults and fractures:integration of three-dimensional seismic data,three-dimensional retrodeformation,and well data on an example of deformation around an inverted fault[J].AAPG Bulletin,2008,92(4):473-485.

    [72] 付曉飛,郭雪,朱麗旭,等.泥巖涂抹形成演化與油氣運(yùn)移及封閉[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào),2012,41(1):52-63.

    Fu Xiaofei,Guo Xue,Zhu Lixu,et al.Formation and evolution of clay smear and hydrocarbon migration and sealing[J].Journal of China University of Mining & Technology,2012,41(1):52-63.

    [73] Zeng Lianbo,Qi Jiafu,Wang Chenggang,et al.The influence of tectonic stress on fracture formation and fluid flow[J].Earth Science Frontiers,2008,15(3):292-298.

    [74] Narr W,Suppe J.Joint spacing in sedimentary rocks[J].Journal of Structural Geology,1991,13(9):1037-1048.

    [75] Lamarche J,Lavenu A P C,Gauthier B D M,et al.Relationships between fracture patterns,geodynamics and mechanical stratigraphy in carbonates(south-east basin,France)[J].Tectonophysics,2012,581:231-245.

    [76] Roberts G,Stewart I.Uplift,deformation and fluid involvement within an active normal fault zone in the Gulf of Corinth,Greece[J].Journal of the Geological Society,1994,151(3):531-541.

    [77] Luo X R,Lei Y H,Zhang L K,et al.Characterization of carrier formation for hydrocarbon migration:concepts and approaches[J].Acta Petrolei Sinica,2012,33(9):428-436.

    [78] 李雙建,周雁,孫冬勝.評(píng)價(jià)蓋層有效性的巖石力學(xué)實(shí)驗(yàn)研究[J].石油實(shí)驗(yàn)地質(zhì),2013,35(5):574-578.

    Li Shuangjian,Zhou Yan,Sun Yongsheng.Rock mechanic experiment study of evaluation on cap rock effectiveness[J].Petroleum Geology & Experiment,2013,35(5):574-578.

    [79] Fossen H.Structural Geology[M]∥Structural geology.Cambridge University Press,2010:1-4.

    [80] Wibberley C A J,Yielding G,Toro G D.Recent advances in the understanding of fault zone internal structure:a review[C]∥Wibberley C A J,Kurz W,Imber J,et al.The internal structure of fault zones,implications for mechanical and fluid-flow properties,2008,299:5-33.

    [81] Holland M,Urai J L,Van der Zee W,et al.Fault gouge evolution in highly overconsolidated claystones[J].Journal of Structural Geology,2006,28:323-332.

    [82] 賈茹,付曉飛,孟令東,等.斷裂及其伴生微構(gòu)造對(duì)不同類型儲(chǔ)層的改造機(jī)理[J].石油學(xué)報(bào),2017,38(3):286-296.

    Jia Ru,Fu Xiaofei,Meng Lingdong,et al.Transform mechanism of fault and its associated microstructures to different types of reservoirs[J].Acta Petrolei Sinica,2017,38(3):286-296.

    [83] 付曉飛,陳哲,閆百泉,等.海拉爾-塔木察格盆地中部斷陷帶油氣富集主控因素分析——斷層和蓋層雙控模式[J].中國(guó)科學(xué):地球科學(xué),2013,43(8):1338-1351.

    Fu Xiaofei,Chen Zhe,Yan Baiquan,et al.Analysis of main controlling factors for hydrocarbon accumulation in central rift zones of the Hailar-Tamtsag basin using a fault-caprock dual control mode[J].Science China:Earth Sciences,2013,43(8):1338-1351.

    [84] Aydin A,Johnson A M.Analysis of faulting in porous sandstones[J].Journal of Structural Geology,1983,5:19-31.

    [85] Fossen H,Bale A.Deformation bands and their influence on fluid flow[J].AAPG Bulletin,2007,91(12):1685-1700.

    [86] Bernard X D,Labaume P,Darcel C,et al.Cataclastic slip band distribution in normal fault damage zones,Nubian sandstones,Suez rift[J].Journal of Geophysical Research Atmospheres,2002,107(7):1-12.

    [87] Bernard X D,Eichhubl P,Aydin A.Dilation bands:A new form of localized failure in granular media[J].Geophysical Research Letters,2002,29(24):2176.

    [88] Knipe R J.Juxtaposition and seal diagrams to help analyze fault seals in hydrocarbon reservoirs[J].AAPG Bulletin,1997,81(2):187-195.

    [89] 吳元燕,平俊彪,付建林,等.中國(guó)油氣藏破壞類型及分布[J].地質(zhì)論評(píng),2002,48(4):377-383.

    Wu Yuanyan,Ping Junbiao,Fu Jianlin,et al.Destructive style and distribution of oil and gas reservoirs in China[J].Geological Review,2002,48(4):377-383.

    [90] 張立寬,羅曉容,宋國(guó)奇,等.油氣運(yùn)移過(guò)程中斷層啟閉性的量化表征參數(shù)評(píng)價(jià)[J].石油學(xué)報(bào),2013,34(1):92-100.

    Zhang Likuan,Luo Xiaorong,Song Guoqi,et al.Quantitative evaluation of parameters to characterize fault opening and sealing during hydrocarbon migration[J].Acta Petrolei Sinica,2013,34(1):92-100.

    猜你喜歡
    蓋層脆性泥巖
    一種零件制造過(guò)程工序脆性源評(píng)價(jià)方法
    風(fēng)化泥巖地質(zhì)斷層水疏排工藝探討
    考慮初始損傷的脆性疲勞損傷模型及驗(yàn)證
    基于能量耗散的頁(yè)巖脆性特征
    高孔低滲泥巖滲流-損傷耦合模型與數(shù)值模擬
    區(qū)域性泥巖蓋層阻止油氣沿輸導(dǎo)斷裂運(yùn)移機(jī)制及其判別方法
    高強(qiáng)度厚壁鋼的回火脆性研究
    大型鑄鍛件(2015年1期)2016-01-12 06:33:06
    CO2咸水層封存中組合蓋層密封效果的影響因素
    塔里木盆地巴-麥地區(qū)古生界油氣蓋層動(dòng)態(tài)演化評(píng)價(jià)
    時(shí)間域瞬變電磁在玄武巖蓋層區(qū)煤田勘探中的應(yīng)用
    观看av在线不卡| 熟女电影av网| a级一级毛片免费在线观看| 观看免费一级毛片| 久久久久久久久久久丰满| 欧美丝袜亚洲另类| 精品人妻熟女av久视频| 97热精品久久久久久| 在线观看美女被高潮喷水网站| 国国产精品蜜臀av免费| 我的女老师完整版在线观看| 日韩精品有码人妻一区| 精品人妻视频免费看| 亚洲精品自拍成人| 免费观看在线日韩| av.在线天堂| 夜夜看夜夜爽夜夜摸| h日本视频在线播放| 国产成人精品福利久久| 十八禁网站网址无遮挡 | tube8黄色片| 在线观看人妻少妇| av在线老鸭窝| 久久精品国产a三级三级三级| 免费黄网站久久成人精品| 亚洲精品456在线播放app| 人妻少妇偷人精品九色| 国产人妻一区二区三区在| 欧美xxxx性猛交bbbb| 亚洲欧美中文字幕日韩二区| 亚洲欧美中文字幕日韩二区| 在线免费观看不下载黄p国产| 最近手机中文字幕大全| 日韩免费高清中文字幕av| 2021少妇久久久久久久久久久| 2021少妇久久久久久久久久久| 国产免费一区二区三区四区乱码| 高清不卡的av网站| 最近中文字幕高清免费大全6| 久久精品久久精品一区二区三区| 最近中文字幕高清免费大全6| 校园人妻丝袜中文字幕| 欧美日韩视频高清一区二区三区二| 啦啦啦中文免费视频观看日本| 国产精品女同一区二区软件| 青春草亚洲视频在线观看| 91aial.com中文字幕在线观看| 欧美亚洲 丝袜 人妻 在线| 91午夜精品亚洲一区二区三区| 午夜福利高清视频| 午夜日本视频在线| 深夜a级毛片| 精品久久久久久久久亚洲| 日本色播在线视频| 国产真实伦视频高清在线观看| 国产精品国产三级专区第一集| 国产乱人视频| 婷婷色麻豆天堂久久| 日韩av不卡免费在线播放| 国产黄色免费在线视频| 亚洲成人中文字幕在线播放| 看非洲黑人一级黄片| 大码成人一级视频| 啦啦啦在线观看免费高清www| 免费av中文字幕在线| 91久久精品国产一区二区三区| 欧美日韩视频高清一区二区三区二| 欧美日韩国产mv在线观看视频 | 美女脱内裤让男人舔精品视频| 亚洲电影在线观看av| 寂寞人妻少妇视频99o| 日韩av不卡免费在线播放| freevideosex欧美| 一区在线观看完整版| 十分钟在线观看高清视频www | 日韩人妻高清精品专区| av播播在线观看一区| 国产精品一及| 日本猛色少妇xxxxx猛交久久| 青青草视频在线视频观看| videossex国产| 人妻系列 视频| 国产爽快片一区二区三区| 亚洲三级黄色毛片| 一个人看视频在线观看www免费| 99久久精品一区二区三区| 日韩在线高清观看一区二区三区| 99re6热这里在线精品视频| 精品国产露脸久久av麻豆| 午夜福利高清视频| 少妇丰满av| 在线观看三级黄色| 久久ye,这里只有精品| 欧美zozozo另类| 亚洲精品国产av成人精品| 亚洲人与动物交配视频| 久热久热在线精品观看| videos熟女内射| a级一级毛片免费在线观看| 欧美日本视频| 黄色日韩在线| 国内少妇人妻偷人精品xxx网站| 国产亚洲5aaaaa淫片| 一区二区av电影网| 麻豆国产97在线/欧美| 精品少妇久久久久久888优播| 久久国产亚洲av麻豆专区| 老司机影院成人| 亚洲精品中文字幕在线视频 | 好男人视频免费观看在线| 国内精品宾馆在线| 在线免费十八禁| 久久久久久久国产电影| h视频一区二区三区| 黄色一级大片看看| 国产女主播在线喷水免费视频网站| 亚洲真实伦在线观看| av免费在线看不卡| 国产亚洲91精品色在线| 偷拍熟女少妇极品色| 男男h啪啪无遮挡| 在线 av 中文字幕| 久久精品夜色国产| 国产欧美日韩一区二区三区在线 | 午夜激情久久久久久久| 国产精品精品国产色婷婷| 成人无遮挡网站| 成年免费大片在线观看| av黄色大香蕉| 亚洲人成网站在线观看播放| 国产在线免费精品| 大话2 男鬼变身卡| 免费人妻精品一区二区三区视频| 国产成人91sexporn| 精品一区在线观看国产| 成人漫画全彩无遮挡| 美女cb高潮喷水在线观看| 久久精品久久久久久久性| 精品亚洲乱码少妇综合久久| 免费人成在线观看视频色| a级一级毛片免费在线观看| 国产在线视频一区二区| 中文字幕av成人在线电影| 男女下面进入的视频免费午夜| 亚洲av.av天堂| 欧美性感艳星| 五月天丁香电影| 欧美3d第一页| 中文在线观看免费www的网站| 亚洲欧洲日产国产| 久久久成人免费电影| 国产精品三级大全| 一区二区三区四区激情视频| 五月天丁香电影| 国产91av在线免费观看| 午夜免费观看性视频| .国产精品久久| 亚洲美女黄色视频免费看| 18禁裸乳无遮挡动漫免费视频| 精品久久久久久久久av| 国产高清有码在线观看视频| 日韩欧美精品免费久久| 一级黄片播放器| 九九爱精品视频在线观看| 在线观看免费视频网站a站| 欧美日韩综合久久久久久| freevideosex欧美| 爱豆传媒免费全集在线观看| 卡戴珊不雅视频在线播放| 久久久久国产精品人妻一区二区| 直男gayav资源| 久久鲁丝午夜福利片| 这个男人来自地球电影免费观看 | 国产亚洲最大av| 18禁裸乳无遮挡免费网站照片| 国产69精品久久久久777片| 亚洲人成网站高清观看| 五月开心婷婷网| 欧美亚洲 丝袜 人妻 在线| 国产亚洲91精品色在线| 啦啦啦视频在线资源免费观看| 少妇 在线观看| 国产精品麻豆人妻色哟哟久久| av视频免费观看在线观看| 久久人妻熟女aⅴ| 黄片wwwwww| 激情 狠狠 欧美| 人妻少妇偷人精品九色| 久久精品国产鲁丝片午夜精品| 99re6热这里在线精品视频| 久久久亚洲精品成人影院| 亚洲精品aⅴ在线观看| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| 国产欧美日韩一区二区三区在线 | 国产男女内射视频| 国精品久久久久久国模美| 熟女人妻精品中文字幕| 在线观看人妻少妇| 日韩,欧美,国产一区二区三区| 91午夜精品亚洲一区二区三区| 亚洲精品乱久久久久久| 天天躁日日操中文字幕| 性色av一级| 91精品一卡2卡3卡4卡| 夜夜看夜夜爽夜夜摸| 亚洲国产av新网站| 91精品国产国语对白视频| 天堂俺去俺来也www色官网| av专区在线播放| 免费观看av网站的网址| 99热国产这里只有精品6| 欧美日韩视频高清一区二区三区二| 熟女人妻精品中文字幕| 又爽又黄a免费视频| 岛国毛片在线播放| av在线观看视频网站免费| 国产精品免费大片| 亚洲va在线va天堂va国产| 高清毛片免费看| 亚洲性久久影院| 亚洲国产高清在线一区二区三| 精品久久久久久电影网| 毛片一级片免费看久久久久| 亚洲精品中文字幕在线视频 | 99热这里只有是精品在线观看| tube8黄色片| 99精国产麻豆久久婷婷| 免费播放大片免费观看视频在线观看| 亚洲av电影在线观看一区二区三区| 色婷婷久久久亚洲欧美| 中文字幕人妻熟人妻熟丝袜美| 99国产精品免费福利视频| 国产精品偷伦视频观看了| 久久久精品94久久精品| 色婷婷久久久亚洲欧美| 少妇精品久久久久久久| 波野结衣二区三区在线| 亚洲美女黄色视频免费看| 成人漫画全彩无遮挡| 国产乱来视频区| 80岁老熟妇乱子伦牲交| 欧美精品一区二区大全| 亚洲国产欧美在线一区| 国语对白做爰xxxⅹ性视频网站| 99热这里只有精品一区| 97超碰精品成人国产| 国产女主播在线喷水免费视频网站| 爱豆传媒免费全集在线观看| 久久国产亚洲av麻豆专区| 纯流量卡能插随身wifi吗| 80岁老熟妇乱子伦牲交| 日韩欧美 国产精品| 美女国产视频在线观看| 精品少妇久久久久久888优播| 免费在线观看成人毛片| 99热这里只有是精品50| 人人妻人人澡人人爽人人夜夜| 日韩不卡一区二区三区视频在线| 高清不卡的av网站| 蜜桃久久精品国产亚洲av| 五月开心婷婷网| 少妇丰满av| 免费观看性生交大片5| 九色成人免费人妻av| 大片免费播放器 马上看| 老熟女久久久| 一个人看的www免费观看视频| 欧美bdsm另类| 91在线精品国自产拍蜜月| 狠狠精品人妻久久久久久综合| 久久久久国产网址| av国产免费在线观看| 一级毛片久久久久久久久女| 天堂俺去俺来也www色官网| 99国产精品免费福利视频| 亚洲天堂av无毛| 天堂8中文在线网| 国产精品99久久99久久久不卡 | 国产成人精品一,二区| 欧美xxxx性猛交bbbb| 成人国产av品久久久| 丝袜脚勾引网站| 国产成人精品一,二区| 亚洲av日韩在线播放| 中文精品一卡2卡3卡4更新| 高清欧美精品videossex| 免费人成在线观看视频色| 熟女电影av网| 亚洲av国产av综合av卡| 国产亚洲午夜精品一区二区久久| 欧美一区二区亚洲| 简卡轻食公司| 日韩在线高清观看一区二区三区| 免费高清在线观看视频在线观看| 久久久精品94久久精品| 免费黄色在线免费观看| 亚洲在久久综合| 国内揄拍国产精品人妻在线| 一区二区三区四区激情视频| 亚洲图色成人| 最近中文字幕2019免费版| 少妇的逼好多水| 青春草国产在线视频| 国产在线男女| 国产 一区精品| 中国国产av一级| 成年av动漫网址| 亚洲国产高清在线一区二区三| 国产亚洲欧美精品永久| 免费黄色在线免费观看| 熟女av电影| 国产一区二区三区av在线| 久久97久久精品| 99热这里只有精品一区| 99精国产麻豆久久婷婷| 成人国产av品久久久| 亚洲av中文av极速乱| 丰满人妻一区二区三区视频av| 最近手机中文字幕大全| 成人美女网站在线观看视频| 国产女主播在线喷水免费视频网站| 欧美精品国产亚洲| 韩国av在线不卡| 欧美激情极品国产一区二区三区 | 亚洲自偷自拍三级| 精品人妻熟女av久视频| 国产伦在线观看视频一区| 亚洲美女视频黄频| 一级毛片我不卡| 在线观看免费视频网站a站| 日韩av在线免费看完整版不卡| 午夜免费鲁丝| 老熟女久久久| 国产免费福利视频在线观看| 国产精品欧美亚洲77777| 美女高潮的动态| 免费av不卡在线播放| 九色成人免费人妻av| 水蜜桃什么品种好| 国产伦在线观看视频一区| 精品久久久久久电影网| 国产精品嫩草影院av在线观看| 国产成人午夜福利电影在线观看| 欧美亚洲 丝袜 人妻 在线| 免费观看无遮挡的男女| 亚洲国产毛片av蜜桃av| 91午夜精品亚洲一区二区三区| 校园人妻丝袜中文字幕| 最近最新中文字幕大全电影3| 国产高潮美女av| 91久久精品国产一区二区三区| 欧美成人精品欧美一级黄| 亚洲国产高清在线一区二区三| 2018国产大陆天天弄谢| 最近的中文字幕免费完整| 看免费成人av毛片| 国产精品一及| 日韩 亚洲 欧美在线| 久久久久久久久久成人| 久久99热这里只有精品18| 亚洲人与动物交配视频| 久久久久久久久久久免费av| 日韩三级伦理在线观看| 下体分泌物呈黄色| 18+在线观看网站| 一区二区av电影网| 亚洲国产色片| 人妻一区二区av| 亚洲av.av天堂| 久久久久精品性色| 人妻制服诱惑在线中文字幕| 日本vs欧美在线观看视频 | 国产人妻一区二区三区在| 久久久久久久久久成人| 搡女人真爽免费视频火全软件| 国产淫片久久久久久久久| 亚洲无线观看免费| 日本vs欧美在线观看视频 | 亚洲精品成人av观看孕妇| 国产高清三级在线| 久久女婷五月综合色啪小说| 免费不卡的大黄色大毛片视频在线观看| 国产黄色视频一区二区在线观看| 日本黄色日本黄色录像| 黑人高潮一二区| 国产精品一二三区在线看| 亚洲电影在线观看av| 久久综合国产亚洲精品| 又爽又黄a免费视频| 深夜a级毛片| 国产美女午夜福利| 黄色日韩在线| 久久99热这里只频精品6学生| 赤兔流量卡办理| 日韩制服骚丝袜av| 最近2019中文字幕mv第一页| 亚洲无线观看免费| 久久久久视频综合| 久久99热6这里只有精品| 夫妻性生交免费视频一级片| 有码 亚洲区| 国产欧美日韩一区二区三区在线 | 一区二区三区乱码不卡18| 18+在线观看网站| 日本vs欧美在线观看视频 | 最近的中文字幕免费完整| 欧美最新免费一区二区三区| 一级a做视频免费观看| 国产v大片淫在线免费观看| 午夜福利网站1000一区二区三区| av在线播放精品| 精品久久久噜噜| 哪个播放器可以免费观看大片| 亚洲欧美清纯卡通| 久久av网站| 日韩三级伦理在线观看| 欧美成人精品欧美一级黄| 国产女主播在线喷水免费视频网站| 亚洲,一卡二卡三卡| 午夜免费观看性视频| 免费播放大片免费观看视频在线观看| 亚洲av免费高清在线观看| 黄色配什么色好看| 久久久精品94久久精品| 日日啪夜夜爽| 国产黄频视频在线观看| 国内精品宾馆在线| 大香蕉久久网| 国产精品三级大全| 久久精品国产a三级三级三级| 亚洲真实伦在线观看| 日韩免费高清中文字幕av| 99国产精品免费福利视频| 一边亲一边摸免费视频| 在线亚洲精品国产二区图片欧美 | 国产精品一二三区在线看| 波野结衣二区三区在线| 久久久久人妻精品一区果冻| 99久久人妻综合| 日本与韩国留学比较| 国产老妇伦熟女老妇高清| 亚洲成色77777| 青春草亚洲视频在线观看| 精品久久久噜噜| 亚洲美女视频黄频| 久久精品国产亚洲网站| 成年美女黄网站色视频大全免费 | 亚洲成色77777| 久久ye,这里只有精品| 亚洲国产毛片av蜜桃av| 日本午夜av视频| 少妇 在线观看| 制服丝袜香蕉在线| 国产免费视频播放在线视频| 另类亚洲欧美激情| 亚洲精品乱久久久久久| 免费黄网站久久成人精品| 亚洲av国产av综合av卡| 色婷婷久久久亚洲欧美| .国产精品久久| 亚洲色图综合在线观看| 免费看不卡的av| 啦啦啦在线观看免费高清www| 国产精品久久久久久久久免| .国产精品久久| 人人妻人人添人人爽欧美一区卜 | 男女边摸边吃奶| 中文字幕精品免费在线观看视频 | 在线观看免费日韩欧美大片 | 精品人妻视频免费看| 人妻少妇偷人精品九色| h日本视频在线播放| 天堂8中文在线网| 久热这里只有精品99| 18禁裸乳无遮挡动漫免费视频| 美女高潮的动态| 午夜日本视频在线| 国产片特级美女逼逼视频| 久久女婷五月综合色啪小说| 视频中文字幕在线观看| 久久久欧美国产精品| 十分钟在线观看高清视频www | 久久韩国三级中文字幕| 欧美+日韩+精品| 青春草亚洲视频在线观看| 日韩一区二区视频免费看| 精品国产三级普通话版| 香蕉精品网在线| 中文欧美无线码| 97精品久久久久久久久久精品| 国产精品久久久久久久电影| 美女福利国产在线 | 亚洲精品,欧美精品| 日韩 亚洲 欧美在线| 寂寞人妻少妇视频99o| 能在线免费看毛片的网站| 伦理电影大哥的女人| 狂野欧美激情性xxxx在线观看| 99热这里只有是精品50| 国产精品秋霞免费鲁丝片| 国产精品久久久久久精品古装| 亚洲av中文av极速乱| 熟女人妻精品中文字幕| 亚洲国产高清在线一区二区三| 最后的刺客免费高清国语| 亚洲国产日韩一区二区| 18禁在线播放成人免费| 亚洲人成网站在线观看播放| 一级a做视频免费观看| 亚洲三级黄色毛片| 日日撸夜夜添| 中文字幕人妻熟人妻熟丝袜美| 国产精品av视频在线免费观看| 亚洲av不卡在线观看| 成人黄色视频免费在线看| 亚洲综合精品二区| 欧美激情极品国产一区二区三区 | 能在线免费看毛片的网站| 女性生殖器流出的白浆| 日日摸夜夜添夜夜爱| 人妻制服诱惑在线中文字幕| 视频中文字幕在线观看| 91精品一卡2卡3卡4卡| 啦啦啦在线观看免费高清www| 国产色爽女视频免费观看| 爱豆传媒免费全集在线观看| 精品亚洲成国产av| 日韩三级伦理在线观看| 少妇精品久久久久久久| 高清日韩中文字幕在线| 亚洲精品乱码久久久v下载方式| 国产一区二区在线观看日韩| 香蕉精品网在线| 日本午夜av视频| 极品教师在线视频| 亚洲天堂av无毛| 嘟嘟电影网在线观看| 亚洲熟女精品中文字幕| 超碰97精品在线观看| 自拍偷自拍亚洲精品老妇| 麻豆乱淫一区二区| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 国内少妇人妻偷人精品xxx网站| 99久国产av精品国产电影| 国产精品久久久久久久电影| 日韩免费高清中文字幕av| 亚洲欧美精品自产自拍| 国产精品久久久久成人av| 日本黄大片高清| 亚洲国产色片| 精品人妻偷拍中文字幕| 水蜜桃什么品种好| 久久久亚洲精品成人影院| 熟女人妻精品中文字幕| 亚洲精品亚洲一区二区| 国产视频内射| 日韩中字成人| 中文字幕免费在线视频6| 少妇人妻久久综合中文| 乱码一卡2卡4卡精品| 各种免费的搞黄视频| 精品一区二区三卡| 亚洲欧洲日产国产| 在线播放无遮挡| 国产精品一区二区在线不卡| 91精品一卡2卡3卡4卡| 偷拍熟女少妇极品色| 99re6热这里在线精品视频| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩东京热| 男男h啪啪无遮挡| 99精国产麻豆久久婷婷| 中文字幕精品免费在线观看视频 | 少妇人妻久久综合中文| 老女人水多毛片| 国产男女内射视频| 人人妻人人看人人澡| 欧美日韩视频高清一区二区三区二| 欧美人与善性xxx| 啦啦啦视频在线资源免费观看| 黑丝袜美女国产一区| 成年人午夜在线观看视频| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品古装| 99热国产这里只有精品6| 中文字幕人妻熟人妻熟丝袜美| 少妇人妻 视频| 成人二区视频| 亚洲欧美精品专区久久| 国产欧美日韩一区二区三区在线 | 黄色欧美视频在线观看| 多毛熟女@视频| 亚洲精华国产精华液的使用体验| 在线观看美女被高潮喷水网站| 色综合色国产| 一区二区三区乱码不卡18| 国产免费一级a男人的天堂| 亚洲欧美日韩另类电影网站 | 啦啦啦中文免费视频观看日本| 黄色日韩在线| 男人和女人高潮做爰伦理| 日韩一本色道免费dvd| 亚洲成人一二三区av| 成人国产麻豆网| 国产成人精品一,二区| 免费av中文字幕在线| 在线精品无人区一区二区三 | 亚洲色图av天堂| 最近手机中文字幕大全| 亚洲av欧美aⅴ国产| 啦啦啦啦在线视频资源| 国产精品秋霞免费鲁丝片| 夜夜爽夜夜爽视频|