王明飛,蘇克露,肖 偉,陳 超,劉曉晶,王靜波

(中國石化 勘探分公司，四川 成都 610000)

相對于常規(guī)油氣藏，頁巖氣的成因類型、富集機(jī)理及生產(chǎn)機(jī)制[1-3]等都具有一定的特殊性，頁巖氣產(chǎn)量高低與泥頁巖內(nèi)部天然微裂縫發(fā)育程度相關(guān)。在頁巖氣勘探突破之前，對裂縫的研究主要集中在常規(guī)油氣儲層，在裂縫發(fā)育的地質(zhì)認(rèn)識基礎(chǔ)上，發(fā)展形成了基于巖心、成像測井和疊前疊后地震資料的裂縫綜合預(yù)測技術(shù)方法，并在常規(guī)儲層勘探中取得了良好的應(yīng)用效果[4-9]。近年來，在碳酸鹽巖裂縫綜合預(yù)測中，探索了應(yīng)力場模擬約束下的裂縫綜合預(yù)測方法，該方法將應(yīng)力場預(yù)測結(jié)果作為宏觀約束條件，定性約束疊前及疊后裂縫預(yù)測成果，分別優(yōu)選單組及多組裂縫發(fā)育區(qū)[10]。在非常規(guī)泥頁巖裂縫研究方面，泥頁巖裂縫識別方法及表征、成因類型及分布規(guī)律、基本參數(shù)及其物性參數(shù)研究、儲集性能評價(jià)等方面已經(jīng)取得了重要成果，推動了頁巖氣儲層富集規(guī)律、精細(xì)評價(jià)、工程設(shè)計(jì)等研究[11-20]。在泥頁巖裂縫平面預(yù)測方面，主要應(yīng)用疊后地震屬性、聚類分析等方法開展了泥頁巖裂縫預(yù)測研究[21-22]，在泥頁巖儲層裂縫發(fā)育方向、強(qiáng)度及平面分布規(guī)律上尚未開展較全面深入的研究。

本文以涪陵頁巖氣田焦石壩區(qū)塊五峰組-龍馬溪組一段為例，在頁巖氣儲層裂縫成因類型認(rèn)識的基礎(chǔ)上，基于焦石壩疊后地震資料，開展研究區(qū)地震體曲率及應(yīng)力場研究，將歸一化應(yīng)力場及體曲率相乘，探索基于應(yīng)力場約束的泥頁巖疊后裂縫預(yù)測技術(shù)，開展泥頁巖裂縫方向及強(qiáng)度的定量-半定量綜合預(yù)測。

1 焦石壩區(qū)塊泥頁巖裂縫類型

焦石壩區(qū)塊是石柱復(fù)向斜、方斗山復(fù)背斜和萬縣復(fù)向斜等多個構(gòu)造單元的結(jié)合部。受印支、燕山和喜馬拉雅三期構(gòu)造運(yùn)動控制，發(fā)育受北東向和南北向兩組斷裂體系控制的北東向展布的菱形斷背斜[16]。背斜頂部較寬緩，呈似箱裝特征。區(qū)內(nèi)主要發(fā)育兩期斷裂，早期呈北東向，晚期呈北西向展布，背斜主體構(gòu)造穩(wěn)定，斷裂不發(fā)育(圖1)。構(gòu)造應(yīng)力控制了斷裂的發(fā)育程度及裂縫發(fā)育的主要方位。成像測井統(tǒng)計(jì)結(jié)果表明，區(qū)內(nèi)頁巖氣儲層五峰組-龍馬溪組一段內(nèi)主要發(fā)育北西向和北東向裂縫，裂縫發(fā)育的主特征方位為北東向。發(fā)育構(gòu)造成因的中-高角度的構(gòu)造縫、低角度的層面滑移縫，沉積成因的低角度層間頁理縫和有機(jī)質(zhì)演化異常壓力縫2類4種裂縫[23-24](圖2b—e)。頁巖裂縫與常規(guī)儲層裂縫相比具有明顯的差異性。以四川盆地元壩地區(qū)須家河組三段常規(guī)砂巖儲層裂縫發(fā)育特征為例，常規(guī)儲層主要以發(fā)育高角度縫為主，延伸范圍大，裂縫寬度多在毫米-厘米級(圖2a)；泥頁巖儲層裂縫以低角度滑移縫及層理縫發(fā)育為主，整體具有“小尺度，低角度”的特征，斷層伴生裂縫發(fā)育于斷褶帶局部構(gòu)造部位。在龍馬溪組一段有機(jī)碳含量較高的炭質(zhì)泥頁巖中，除大量低角度滑移縫及層理縫發(fā)育外，有機(jī)質(zhì)演化異常壓力縫普遍發(fā)育，通過氬離子拋光掃描電鏡，發(fā)現(xiàn)這種裂縫的縫面一般不規(guī)則，微裂縫寬度一般為0.02～1 μm；另外發(fā)育少量高角度縫，角度70°～90°，延伸長度一般為30～130 mm，總縫密度為0.5～3條/m，縫寬為0.1～2 mm，呈充填-半充填狀態(tài)。五峰組網(wǎng)狀縫發(fā)育(JY1井)，但多被方解石充填或半充填。龍馬溪組局部斷褶構(gòu)造頁巖氣層發(fā)育開啟的高角度斷層伴生縫(JY5井，烏江斷褶帶)(圖2e)。

2 地震裂縫預(yù)測

根據(jù)裂縫的規(guī)模及當(dāng)前的測量和預(yù)測手段，一般可以將裂縫分為兩類。一是大尺度裂縫，延伸距離長，可達(dá)幾十米到幾千米，斷距大，地震剖面可以直接解釋及識別；二是中-小尺度裂縫，延伸距離小，可達(dá)幾厘米到幾十米，斷距小，巖心和成像測井顯示為小尺度裂縫。根據(jù)裂縫性儲層地震響應(yīng)特征分析可知，當(dāng)?shù)卣鹜噍S呈現(xiàn)復(fù)波、扭曲、雜亂以及不連續(xù)等反射特征時(shí)，地層中可能發(fā)育裂縫[25]。因此，反映疊后地震同相軸橫向變化的地震屬性，如體曲率，螞蟻體和方差體等地震體屬性可以開展裂縫定性預(yù)測。針對泥頁巖裂縫自身特征，借鑒常規(guī)油氣勘探裂縫預(yù)測方法，除了對頁巖氣儲層大尺度裂縫實(shí)現(xiàn)有效預(yù)測外，如何實(shí)現(xiàn)泥頁巖儲層中小尺度、中-低角度裂縫的有效預(yù)測，目前尚未形成較系統(tǒng)的技術(shù)方法。本文基于焦石壩區(qū)塊疊后地震資料，在泥頁巖儲層地球物理特征認(rèn)識的基礎(chǔ)上[26-28]，從泥頁巖裂縫發(fā)育的控制因素入手，將構(gòu)造應(yīng)力對泥頁巖裂縫形成的影響引入泥頁巖裂縫預(yù)測中，開展了泥頁巖裂縫地震預(yù)測技術(shù)方法的探索。

2.1 疊后裂縫預(yù)測技術(shù)方法對比

目前，常規(guī)疊后地震裂縫預(yù)測技術(shù)，有相干體、地震方差體、曲率體、螞蟻?zhàn)粉櫟燃夹g(shù)[8,29]。地震相干體裂縫檢測技術(shù)，以及基于空間某點(diǎn)與相鄰點(diǎn)的相似系數(shù)或互相關(guān)函數(shù)，計(jì)算該點(diǎn)的相干屬性最后得到相干數(shù)據(jù)體，該方法檢測大尺度的斷裂效果較好。方差體技術(shù)是基于三維地震某點(diǎn)相鄰道地震信號(如振幅、相位等)之間的相似性來描述地層、巖性等的橫向變化，通過計(jì)算樣點(diǎn)的方差值，揭示三維地震數(shù)據(jù)體中的不連續(xù)信息，從而進(jìn)行斷層與裂縫的識別，該技術(shù)識別中等尺度裂縫效果較好。螞蟻?zhàn)粉櫦夹g(shù)是通過對螞蟻覓食的研究，建立函數(shù)關(guān)系，對三維地震數(shù)據(jù)體中斷裂及裂縫建立螞蟻?zhàn)粉櫂颖荆ㄟ^螞蟻?zhàn)粉櫃z測中-小尺度裂縫。這些技術(shù)方法均為常規(guī)油氣勘探中較常見的裂縫預(yù)測方法。

圖1 焦石壩區(qū)塊五峰組底界構(gòu)造圖Fig.1 Bottom structure map of Wufeng Fm in Jiaoshiba area

圖2 焦石壩區(qū)塊五峰組—龍馬溪組砂巖與泥頁巖巖心裂縫及微觀特征Fig.2 Cores and thin sections showing differences of fractures between sandstone shale in Wufeng-Longmaxi Fms in Jiaoshiba areaa.YL7井，T3x3,埋深3 605 m,深灰色鈣質(zhì)細(xì)砂巖，發(fā)育多組高角度裂縫，方解石充填；b.JY3井，龍馬溪組，埋深2 408.21～408.43 m，層理縫；c.JY井，龍馬溪組，埋深2 547.60 m，異常壓力縫；d.JY1井，五峰組， 埋深2 410 m網(wǎng)狀縫；e.JY5井，龍馬溪組，埋深3 072 m，高角度縫

借鑒常規(guī)儲層裂縫預(yù)測方法，針對沉積較穩(wěn)定的泥頁巖地層中發(fā)育的裂縫，開展了基于地震方差體技術(shù)和螞蟻體技術(shù)的泥頁巖裂縫預(yù)測研究。沿五峰組-龍馬溪組底界提取方差體屬性(圖3a)及螞蟻體屬性(圖3b)?？梢钥闯觯瑑烧咚坍嫷牧芽p細(xì)節(jié)都較好，取得了一定的預(yù)測效果，在構(gòu)造翼部及大斷裂發(fā)育區(qū)，斷裂及裂縫規(guī)模較大，預(yù)測效果較好。但是在構(gòu)造穩(wěn)定區(qū)，受地震數(shù)據(jù)本身噪聲以及泥頁巖裂縫地震響應(yīng)弱敏感性的影響，裂縫預(yù)測平面顯示的裂縫的主特征方位不突出，預(yù)測平面分布趨勢與地質(zhì)統(tǒng)計(jì)規(guī)律吻合度較低。因此，為了得到更佳的泥頁巖裂縫預(yù)測效果，在常規(guī)裂縫預(yù)測方法嘗試的基礎(chǔ)上，進(jìn)一步探索了更加逼近泥頁巖裂縫地質(zhì)發(fā)育規(guī)律的裂縫預(yù)測技術(shù)方法。

圖3 焦石壩區(qū)塊常規(guī)地震裂縫預(yù)測效果對比Fig.3 Comparison of conventional seismic fracture prediction methods applied in Jiaoshiba areaa.五峰組-龍馬溪組一段方差體屬性；b.五峰組-龍馬溪組一段螞蟻體屬性

2.2 應(yīng)力場約束的泥頁巖裂縫預(yù)測

通過以上常規(guī)裂縫預(yù)測方法應(yīng)用到泥頁巖裂縫預(yù)測分析可知，由于構(gòu)造穩(wěn)定區(qū)泥頁巖儲層裂縫發(fā)育尺度較小，裂縫發(fā)育的地震響應(yīng)敏感度不及常規(guī)儲層尺度較大的中-高角度裂縫,以及研究區(qū)構(gòu)造翼部及大斷裂發(fā)育區(qū)的規(guī)模裂縫。因此，開展了增強(qiáng)裂縫地震響應(yīng)效果的地震針對性處理技術(shù)研究。在通過已鉆井成像測井資料確定裂縫主特征發(fā)育方位的基礎(chǔ)上，應(yīng)用應(yīng)力場對裂縫方位及應(yīng)變強(qiáng)度的有效預(yù)測，建立基于應(yīng)力場約束的有效裂縫模型，實(shí)現(xiàn)泥頁巖裂縫方位及強(qiáng)度的定量-半定量預(yù)測。

2.2.1 地震數(shù)據(jù)針對性處理

為了提高斷裂及裂縫成像效果，開展了焦石壩區(qū)塊疊前時(shí)間偏移地震資料的針對性處理，聯(lián)合應(yīng)用構(gòu)造平滑、振幅包絡(luò)與一階求導(dǎo)等技術(shù)方法，提高原始地震數(shù)據(jù)體對微小斷層的響應(yīng)能力。原始地震相干體屬性與針對性處理后相干體屬性對比分析表明，處理后的地震數(shù)據(jù)體有效地提高了斷層及裂縫在地震剖面上的辨識能力(圖4)，沿泥頁巖底界(即五峰組底界)提取針對性處理前后相干地震數(shù)據(jù)體屬性，可以觀測到斷裂及裂縫的更多細(xì)節(jié)(圖5)。

圖4 焦石壩區(qū)塊疊后地震數(shù)據(jù)針對性處理Fig.4 Post-stack seismic data processing in Jiaoshiba area

圖5 焦石壩區(qū)塊五峰-龍馬溪組一段相干體屬性平面圖Fig.5 Coherence attribute of the Wufeng Fm-1st member of Longmaxi Fm in Jiaoshiba areaa.針對性處理前；b.針對性處理后

2.2.2 構(gòu)造應(yīng)力場分析

地應(yīng)力作用形成的儲層裂縫、斷層及構(gòu)造是油氣運(yùn)移、聚集的通道和場所之一。通過地層應(yīng)力場分析，可以預(yù)測構(gòu)造成因的裂縫在研究區(qū)域的發(fā)育和分布規(guī)律，從而快速優(yōu)選有利裂縫發(fā)育方位。在焦石壩區(qū)塊，針對五峰組-龍馬溪組一段泥頁巖儲層，從構(gòu)造力學(xué)出發(fā)，利用地層的幾何信息(構(gòu)造面)、巖性信息(速度、密度)、巖石物理彈性信息(楊氏模量、泊松比)等(表1)建立泥頁巖層的地質(zhì)模型、力學(xué)模型和數(shù)學(xué)模型，基于彎曲薄板理論開展五峰組-龍馬溪組一段泥頁巖地層應(yīng)力場預(yù)測。根據(jù)薄板理論，假設(shè)z為薄板在垂直方向的厚度，則薄板內(nèi)任意一點(diǎn)(x,y)的應(yīng)力分量為[30]：

(1)

將構(gòu)造曲率和井震結(jié)合的楊氏模量、泊松比疊前反演結(jié)果代入上式可以對構(gòu)造應(yīng)力場進(jìn)行預(yù)測。

表1 焦石壩區(qū)塊JY1-JY4井彈性參數(shù)Table 1 Elastic parameters of JY1-JY4 wells in Jiaoshiba area

預(yù)測結(jié)果顯示，研究區(qū)主要發(fā)育早期北西向構(gòu)造應(yīng)力及晚期北東向構(gòu)造應(yīng)力，整體以北西向?yàn)橹?，?gòu)造應(yīng)變強(qiáng)度較大(圖6)。經(jīng)兩期構(gòu)造應(yīng)力改造，是研究區(qū)五峰組-龍馬溪組一段頁巖氣儲層發(fā)育北東向及北西向裂縫的內(nèi)因，符合該區(qū)地質(zhì)特征。

2.2.3 應(yīng)力場約束的有效裂縫預(yù)測

通過針對性處理的地震數(shù)據(jù)體的空間特征掃描，得到能夠?qū)α芽p空間分布特征進(jìn)行精細(xì)刻畫的地震曲率體。高精度曲率體呈現(xiàn)出的是地下地層的幾何屬性，缺乏對裂縫發(fā)育方向及強(qiáng)度的有效描述。構(gòu)造應(yīng)力場是控制儲層裂縫形成、分布以及發(fā)育程度的重要因素，裂縫的發(fā)育特征也在一定程度上反映出構(gòu)造應(yīng)力的狀態(tài)，二者具有較好的對應(yīng)關(guān)系。

因此，在高精度曲率體預(yù)測裂縫的基礎(chǔ)上，將應(yīng)力場特征對裂縫發(fā)育的影響引入到裂縫預(yù)測中，探索建立了應(yīng)力場約束的裂縫預(yù)測模型：

(2)

圖6 焦石壩區(qū)塊五峰組-龍馬溪組一段應(yīng)力場Fig.6 Stress field map of the Wufeng Fm-1st member of Longmaxi Fm in Jiaoshiba area

(3)

(4)

地震體曲率值與裂縫發(fā)育程度呈正相關(guān)關(guān)系，曲率值越大，裂縫發(fā)育程度越高，反之越低。應(yīng)力場的強(qiáng)度越大表明裂縫發(fā)育越強(qiáng)、越集中；反之，則裂縫發(fā)育程度較弱。本文將表征裂縫發(fā)育特征的體曲率與裂縫發(fā)育方向及強(qiáng)度的應(yīng)力場進(jìn)行歸一化，建立兩者相乘

表2 曲率體、應(yīng)力場強(qiáng)度與裂縫發(fā)育程度的關(guān)系Table 2 Correlation among seismic curvature,stress field and F(x,y) fracture intensity

基于上述裂縫預(yù)測模型[公式(2)]對焦石壩區(qū)塊五峰組-龍馬溪組一段頁巖氣儲層進(jìn)行了應(yīng)力場約束的裂縫預(yù)測(圖7)。與方差體、螞蟻體(圖3)裂縫預(yù)測效果相比，應(yīng)力場約束的裂縫預(yù)測優(yōu)勢明顯，該方法既能實(shí)現(xiàn)對有效裂縫及方位的較為準(zhǔn)確的預(yù)測，又能實(shí)現(xiàn)有效裂縫發(fā)育強(qiáng)度的平面預(yù)測。同時(shí)，通過對有效裂縫發(fā)育方位預(yù)測結(jié)果分析可知(顏色代表不同的方位)，焦石壩區(qū)塊主要發(fā)育北東向及北西向有效裂縫，以北東向?yàn)橹?，符合地質(zhì)統(tǒng)計(jì)規(guī)律。從有效裂縫發(fā)育強(qiáng)度預(yù)測分析可知，在焦石壩似箱狀背斜構(gòu)造主體區(qū)，構(gòu)造相對穩(wěn)定，應(yīng)力強(qiáng)度中等，裂縫程度相對較低，強(qiáng)度中-弱為主(裂縫強(qiáng)度：0～0.125)，適中的微裂縫發(fā)育，有利于頁巖氣的富集及保存，是該區(qū)頁巖氣井富集高產(chǎn)的因素之一(JY1井測試產(chǎn)量20.3×104m3/d)。而在構(gòu)造翼部或大型斷裂(褶)發(fā)育區(qū)，構(gòu)造應(yīng)變強(qiáng)度較大，斷裂發(fā)育，中高強(qiáng)度、中高角度斷裂(褶)伴生裂縫集中發(fā)育(裂縫強(qiáng)度：0.125～0.75)，大型斷裂帶與局部規(guī)模裂縫疊加，造成頁巖氣藏局部泄壓，頁巖氣井產(chǎn)量相對下降。JY5井距離烏江斷層1 200 m，位于烏江斷褶帶，斷裂及裂縫疊加，測試產(chǎn)量與其它勘探開發(fā)井相比較低，為4.5×104m3/d。

勘探實(shí)踐表明，靠近裂縫發(fā)育帶的鉆井產(chǎn)量相對較低，鉆井過程中易發(fā)生裂縫性漏失。以實(shí)鉆井為例，JY3HF井及JY17-3HF井分別位于有效裂縫弱發(fā)育區(qū)及有效裂縫中強(qiáng)發(fā)育區(qū)，在鉆井過程中JY3HF水平井段并未發(fā)生鉆井液漏失現(xiàn)象，但是JY17-3HF井水平井段在五峰組-龍馬溪組一段泥頁巖層段內(nèi)發(fā)生多段裂縫性漏失(表3)。通過與本次有效裂縫預(yù)測結(jié)果剖面標(biāo)定來看，JY17-3HF井漏失點(diǎn)位置與有效裂縫預(yù)測地震剖面中裂縫發(fā)育點(diǎn)具有一一對應(yīng)關(guān)系(圖8)，進(jìn)一步說明應(yīng)用該方法開展泥頁巖裂縫預(yù)測的可靠性，對頁巖氣儲層評價(jià)及工程施工具有一定的指導(dǎo)意義。

圖7 焦石壩區(qū)塊五峰組-龍馬溪組一段基于應(yīng)力場約束的裂縫綜合預(yù)測Fig.7 Comprehensive prediction of fractures in the Wufeng Fm-1st member of Longmaxi Fm in Jiaoshiba area under the constraint of stress fielda.五峰組-龍馬溪組一段有效裂縫及方位預(yù)測;b.五峰組-龍馬溪組一段有效裂縫發(fā)育強(qiáng)度預(yù)測

井深/m層位漏失井段/m厚度/m漏失層位漏失量/m3平均漏速/(m3·h-1)漏失段主要巖性2600．00龍馬溪組2600．00～2655．3755龍馬溪組25．642．0～3．0灰黑色炭質(zhì)頁巖2865．14龍馬溪組2865．14～2880．0015龍馬溪組17．603．0～4．0灰黑色炭質(zhì)頁巖2948．10龍馬溪組2948．10～2958．4210龍馬溪組61．001．6～48．0灰黑色炭質(zhì)頁巖3268．00龍馬溪組3268．00～3355．0087龍馬溪組18．601．6灰黑色炭質(zhì)頁巖3538．00龍馬溪組3538．00～3618．0080龍馬溪組15．461．3灰黑色炭質(zhì)頁巖3708．00龍馬溪組3708．00～3727．0019龍馬溪組13．801．2灰黑色炭質(zhì)頁巖4288．00龍馬溪組4280．00龍馬溪組16．342．5灰黑色炭質(zhì)頁巖

圖8 焦石壩區(qū)塊JY17-3HF井與JY3HF井裂縫預(yù)測成果剖面對比Fig.8 Comparison of fracture prediction profiles between JY3HF and JY17-3HF wells in Jiaoshiba areaa.過JY17-3HF井裂縫預(yù)測+地震疊合剖面；b.過JY3HF井裂縫預(yù)測+地震疊合剖面

3 結(jié)論

1) 應(yīng)用常規(guī)裂縫預(yù)測方法難以對地震響應(yīng)較弱的“中小尺度，低角度”的泥頁巖裂縫實(shí)現(xiàn)有效預(yù)測，通過對疊后地震資料開展針對性處理，可以增強(qiáng)裂縫地震響應(yīng)特征，改善預(yù)測效果。

2) 通過基于針對性處理疊后地震數(shù)據(jù)的高精度地震體曲率屬性研究以及基于彎曲波板理論的應(yīng)力場模擬，獲得了體現(xiàn)裂縫發(fā)育程度的體曲率屬性及裂縫發(fā)育方向與強(qiáng)度的應(yīng)力場屬性。將體曲率屬性與應(yīng)力場屬性進(jìn)行歸一化處理，依據(jù)裂縫發(fā)育程度與地震體曲率屬性及應(yīng)力場屬性之間的相互關(guān)系，建立了歸一化體曲率與應(yīng)力場相乘的裂縫預(yù)測模型。通過該模型可以實(shí)現(xiàn)頁巖氣儲層裂縫發(fā)育方向及強(qiáng)度的定量-半定量有效預(yù)測，預(yù)測結(jié)果符合區(qū)塊內(nèi)地質(zhì)裂縫發(fā)育規(guī)律。

3) 應(yīng)力場約束的裂縫預(yù)測結(jié)果表明，在焦石壩地區(qū)箱狀背斜構(gòu)造主體，五峰組-龍馬溪組一段地層較穩(wěn)定，構(gòu)造形變較弱，發(fā)育北東為主的中-弱程度的裂縫，是該局部區(qū)域富集高產(chǎn)的重要因素之一；而在大型斷裂(褶)發(fā)育區(qū)，構(gòu)造應(yīng)變強(qiáng)度較大，中強(qiáng)程度裂縫集中發(fā)育，規(guī)模斷裂與規(guī)模裂縫疊加是頁巖氣局部富集保存及鉆井工程的相對不利因素。

4) 不同裂縫強(qiáng)度發(fā)育區(qū)井震標(biāo)定成果表明，實(shí)鉆井鉆井工程裂縫性漏失點(diǎn)與應(yīng)力場約束的裂縫預(yù)測成果具有一一對應(yīng)關(guān)系，進(jìn)一步驗(yàn)證了該方法的有效性和實(shí)用性，對頁巖氣儲層評價(jià)及工程施工具有一定的指導(dǎo)意義。

[1] 王祥，劉玉華，張敏，等.頁巖氣形成條件及成藏影響因素研究[J].天然氣地球科學(xué)，2010，21(2)：350-356.

Wang Xiang,Liu Yuhua,Zhang Min,et al.Conditions of formation and accumulation for shale gas[J].Natural Gas Geoscience,2010,21(2):350-356.

[2] 張金川，薛會，張德明，等.頁巖氣及其成藏機(jī)理[J].現(xiàn)代地質(zhì)，2003，17(4)：466.

Zhang Jinchuan,Xue Hui,Zhang Deming,et al.Shale gas and its accumulation mechanism[J].Geoscience,2003,17(4):466.

[3] 張金川，金之鈞，袁明生.頁巖氣成藏機(jī)理和分布[J].天然氣工業(yè)，2004，24(7)：15-18.

Zhang Jinchuan,Jin Zhijun,Yuan Mingsheng.Reservoiring mechanism of shale gas and its distribution[J].Natural Gas Industry,2004,24(7):15-18.

[4] 蘇培東，秦啟榮，黃潤秋.儲層裂縫預(yù)測研究現(xiàn)狀與展望[J].西南石油學(xué)院學(xué)報(bào)，2005，(5)：14-17.

Su Peidong,Qin Qirong,Huang Runqiu.Prospectsand statusfor the study on reservoir fractures[J].Journal of Southwest Petroleum Institute,2005,(5):14-17.

[5] 童亨茂.成像測井資料在構(gòu)造裂縫預(yù)測和評價(jià)中的應(yīng)用[J].天然氣工業(yè)，2006，(9)：58-61.

Tong Hengmao.Application of imaging well logging data in prediction of structural fracture[J].Natural Gas Industry,2006,(9):58-61.

[6] 孫松領(lǐng)，李琦，李娟，等.低滲透砂巖儲層構(gòu)造裂縫預(yù)測及開啟性分析[J].特種油氣藏，2007，(1)：30-33.

Sun Songling,Li Qi,Li Juan,et al.Structural fracture prediction and aperture analysis for low permeability sandstone reservoir[J].Special Oil and Gas Reservoirs,2007,(1):30-33.

[7] 時(shí)華星，王秀玲，季玉新，等.構(gòu)造正反演裂縫預(yù)測方法及應(yīng)用實(shí)例[J].石油物探，2004，(4)：337-340.

Shi Hoaxing,Wang Xiuling,Ji Yuxin,et al.Structurally forward and inverse fracture prediction and its application[J].Geophysical Prospection for Petroleum,2004,43(4):337-340.

[8] 劉振峰，曲壽利，孫建國，等.地震裂縫預(yù)測技術(shù)研究進(jìn)展[J].石油物探，2012，(2)：191-198.

Liu Zhenfeng,Qu Shouli,Sun Jianguo,et al.Progress of seismic fracture characterization technology[J].Geophysical Prospection for Petroleum,2012,(2):191-198.

[9] 蔡希源，唐建明，陳本池.三維多分量地震技術(shù)在川西新場地區(qū)深層致密砂巖裂縫檢測中的應(yīng)用[J].石油學(xué)報(bào)，2010，(5)：737-743.

Cai Xiyuan,Tang Jianming,Chen Benchi.A 3D multi-component seismic exploration technique applied to the detection of deep-seated tight sandstone fractures in the Xinchang area of the western Sichuan Basin[J].Acta Petrolei Sinica,2010,(5):737-743.

[10] 狄貴東，孫贊東，龐雄奇，等．應(yīng)力場模擬約束下的碳酸鹽巖裂縫綜合預(yù)測[J].石油物探，2016,55(1):150-156

Di Guidong,Sun Zandong,Pang Xiongqi,et al.Comprehensive fracture prediction technology constrained by stress field simulation:a case study from ZG8 area of central Tarim basin[J].Geophysical Prospection for Petroleum,2016,55(1):150-156.

[11] Curtis J B.Fractured shale-gas systems[J].AAPG Bulletin,2002,86(11):1921-1938.

[12] Montgomery S L,Jarvie D M,Bowker K A,et al.Mississippi-an Barnett shale,Fort Worth basin,north-central Texas:Gas-shale play with multi-trillion cubic foot potential[J].AAPG Bulletin,2005,89(2):155-175.

[13] Bowker K A.Barnett shale gas production,Fort Worth basin:Issues and discussion[J].AAPG Bulletin,2007,91(4):523-533.

[14] Zeng W,Zhang J,Ding W,et al.Fracture development in Paleozoic shale of Chongqing area (South China).Part one:Fracture characteri-stics and comparative analysis of main controlling factors[J].Journal of Asian Earth Sciences,2013,75(8):251-266.

[15] Zeng W,Ding W,Zhang J,et al.Fracture development in Paleozoic shale of Chongqing area (South China).Part two:Numerical simulation of tectonic stress field and prediction of fractures distribution[J].Journal of Asian Earth Sciences,2013,75(8):267-279.

[16] 郭彤樓，張漢榮.四川盆地焦石壩頁巖氣田形成與富集高產(chǎn)模式[J].石油勘探與開發(fā)，2014，41(1)：28-36.

Guo Tonglou,Zhang Hanrong.Formation and enrichment mode of Jiaoshiba shale gas field,Sichuan Basin[J].Petroleum Exploration and Development,2014,41(1):28-36.

[17] 郭旭升.南方海相頁巖氣“二元富集”規(guī)律——四川盆地及周緣龍馬溪組頁巖氣勘探實(shí)踐認(rèn)識[J].石油勘探與開發(fā)，2014，88(7)：209-1218.

Guo Xusheng.Rules of two-factor enrichiment for marine shale gas in southern China——understanding from the Longmaxi Formation shale gas in Sichuan Basin and its surrounding area[J].Petroleum Exploration and Development,2014,88(7):209-1218.

[18] 曾維特，張金川，丁文龍，等.延長組頁巖儲層納米級孔隙特征及影響因素——以鄂爾多斯盆地柳坪171井為例[J].煤炭學(xué)報(bào),2014,39(6):1118-1126.

Zeng Weite,Zhang Jinchuan,Ding Wenlong,et al.Characteristics and influence factors of nanopores in Yanchang shale reservoir:A case study of Liuping-171 Well in Erdos Basin[J].Journal of China Coal Society,2014,39(6):1118-1126.

[19] 丁文龍，曾維特，王濡岳，等.頁巖儲層構(gòu)造應(yīng)力場模擬與裂縫分布預(yù)測方法及應(yīng)用[J].地學(xué)前緣,2016,23(2):63-74.

Ding Wenlong,Zeng Weite,Wang Ruyue,et al.Methods and application of tectonic stress field numerical simulation and fractures distribution prediction for shale reservoir[J].Earth Science Frontier,2016,23(2):63-74.

[20] 曾維特，丁文龍，張金川，等.渝東南-黔北地區(qū)下寒武統(tǒng)牛蹄塘組頁巖裂縫有效性研究[J].地學(xué)前緣,2016,23(1):96-106.

Zeng Weite,Ding Wenlong,Zhang Jinchuan,et al.Research on the fracture effectiveness of Lower Cambrian Niutitang Shale in Southeast of Chongqing and North of Guizhou area (South China)[J].Earth Science Frontier,2016,23(1):96-106.

[21] 唐誠.儲層裂縫表征及預(yù)測研究進(jìn)展[J].科技導(dǎo)報(bào)，2013，31(21)：74-79.

Tang Cheng.Progress in fracture characterization and prediction[J].Science and Technology Review,2013,31(21):74-79.

[22] 陸明華，駱璞，姜傳芳，等.地震屬性技術(shù)在頁巖裂縫預(yù)測中的應(yīng)用[J].石油天然氣學(xué)報(bào)，2013，35(8)：62-64.

Lu Minghua,Luo Pu,Jiang Chuanfang,et al.Application of seismic attribute technology in shale fracture prediction[J].Journal of Oil and Gas Technology,2013,35(8):62-64.

[23] 郭旭升，胡東風(fēng)，文治東，等.四川盆地及周緣下古生界海相頁巖氣富集高產(chǎn)主控因素——以焦石壩地區(qū)五峰組—龍馬溪組為例[J].中國地質(zhì)，2014，41(3)：893-901.

Guo Xusheng,Hu Dongfeng,Wen Zhidong,et al.Major factors controlling the accumulation and high productivity in marine shale gas in the Lower Paleozoic of Sichuan Basin and its periphery:A case study of the Wufeng-Longmaxi Formation of Jiaoshiba area[J].Geology in China,2014,41(3):893-901.

[24] 郭旭升，李宇平，劉若冰，等.四川盆地焦石壩地區(qū)龍馬溪組頁巖微觀孔隙結(jié)構(gòu)特征及其控制因素[J].天然氣工業(yè)，2014，34(6)：9-16.

Guo Xusheng,Li Yuping,Liu Ruobing,et al.Characteristics and controlling factors of micro-pore structures of Longmaxi Shale Play in the Jiaoshiba area,Sichuan Basin[J].Natural Gas Industry,2014,34(6):9-16.

[25] 仇鵬，牟中海，蔣裕強(qiáng)，等.裂縫性儲層的地震響應(yīng)特征—以川中地區(qū)致密砂巖儲層為例[J].油氣田開發(fā)，2011，29(3)：49-53.

Qiu Peng,Mu Zhonghai,Jiang Yuqiang,et al.Seismic response of fractured reservoirs—taking the tight sandstone reservoirs in the central region of Sichuan for example [J].Oil and Gas Field Development,2011,29(3):49-53.

[26] 陳祖慶.海相頁巖TOC地震定量預(yù)測技術(shù)及其應(yīng)用——以四川盆地焦石壩地區(qū)為例[J].天然氣工業(yè),2014,34(6):24-29.

Chen Zuqin.Quantitative seismic prediction technique of marine shale TOC and its application:A case from the Longmaxi Shale Play in the Jiaoshiba area,Sichuan Basin[J].Natural Gas Industry,2014,34(6):24-29.

[27] 李金磊,尹正武.四川盆地焦石壩地區(qū)頁巖氣儲層地震定量預(yù)測方法[J].石油物探,2015,54(3):324-330.

Li Jinlei,Yin Zhengwu.Seismic quantitative prediction method of shale gas reservoirs in the Jiaoshiba Area,Sichuan Basin[J].Geophysical Prospection for Petroleum,2015,54(3):324-330.

[28] 王明飛，陳超，屈大鵬，等.涪陵頁巖氣田焦石壩區(qū)塊五峰組—龍馬溪組一段頁巖氣儲層地球物理特征分析[J].石油物探，2015，54(5)：613-620.

Wang Mingfei,Chen Chao,Qu Dapeng et al.The geophysical characteristics of shale gas reservoir fron Wufeng member to Longmaxi member in Jiaoshiba block of Fulin shale gasfield[J].Geophysical Prospection for Petroleum,2015,54(5):613-620.

[29] 劉雷頌，高軍，代雙河，等.“螞蟻?zhàn)粉櫋绷芽p預(yù)測技術(shù)在中東AD油田開發(fā)中的應(yīng)用[J].石油天然氣學(xué)報(bào)，2013，35(8)：57-61.

Liu Leisong,Gao Jun,Dai Shuanghe,et al.Application of “ant tracki-ng” fracture prediction technique in development of AD oilfild[J].Journal of Oil and Gas Technology,2013,35 (8):57-61.

[30] 余夫，金衍，陳勉，等.基于薄板理論的碳酸鹽巖地層壓力檢測方法探討[J].石油鉆探技術(shù)，2014，42(5)：57-61.

Yu Fu,Jin Yan,Chen Mian,et al.Discussion on a formation pore pressure detection method for carbonate rocks based on the thin plate theory[J].Petroleum Drilling Techniques,2014,42 (5):57-61.