楊鋒，朱春啟，王新海,3，屈海州，張福祥，牛新年

（1.油氣資源與勘探技術(shù)教育部重點實驗室（長江大學）；2.中國石油長城鉆探工程有限公司；3.中國石油大學（北京）；4.西南石油大學；5.中國石油塔里木油田公司）

0 引言

庫車坳陷位于塔里木盆地北緣，北部與天山相鄰，南面緊靠塔北隆起，整體呈北東東向展布，面積2.85×104km2，區(qū)域構(gòu)造上屬天山褶皺帶南麓前陸前淵坳陷[1-4]。庫車坳陷大北—克深地區(qū)主要含油層系為白堊系巴什基奇克組，巖性為致密砂巖，孔隙度為6%～8%。錄井和測井資料顯示該地區(qū)裂縫十分發(fā)育[5-6]，大北地區(qū)取心井裂縫傾角為70°～90°，寬0.1～5.0 mm，部分井段巖心極為破碎，大北202井日產(chǎn)高達110×104m3，大北—克深地區(qū)平均日產(chǎn)40×104m3左右。致密儲集層產(chǎn)能的高低與裂縫的發(fā)育程度密切相關(guān)，如何定量描述裂縫與產(chǎn)能的關(guān)系，目前國內(nèi)外尚無好的解決辦法。筆者根據(jù)裂縫定量描述資料，建立基于裂縫產(chǎn)能系數(shù)預測氣井產(chǎn)能的模型，研究了裂縫均勻分布和非均勻分布對產(chǎn)能的影響，為庫車地區(qū)低孔裂縫性砂巖的產(chǎn)能預測提供了一種新的思路。

1 成像測井裂縫參數(shù)定量統(tǒng)計

裂縫參數(shù)包括裂縫類型、產(chǎn)狀、長度、張開度、總縫條數(shù)、有效縫條數(shù)、縫密度、原始面縫率、有效面縫率、充填率等。庫車地區(qū)巴什基奇克組裂縫類型以構(gòu)造縫為主，少量溶蝕縫。裂縫產(chǎn)狀包括裂縫的傾向、傾角。裂縫長度是巖心或者成像測井圖像中所能觀察到的裂縫延伸長度。裂縫張開度（寬度）是兩個壁面之間的垂直距離。縫密度是裂縫條數(shù)與統(tǒng)計長度的比值。原始面縫率是巖心中裂縫分布面積占所統(tǒng)計的巖心表面積的百分數(shù)。有效面縫率是指裂縫中未被充填部分的面積占所統(tǒng)計的巖心表面積的百分數(shù)。充填率為裂縫中被充填部分的面積占所統(tǒng)計的巖心表面積的百分數(shù)。

在統(tǒng)計巖心內(nèi)的裂縫參數(shù)時，裂縫類型、產(chǎn)狀、長度、條數(shù)等肉眼均可見。裂縫張開度采用刻度放大鏡等工具讀取。面縫率則主要由網(wǎng)格法定量統(tǒng)計。網(wǎng)格法即用帶有刻度的透明方格紙包覆于巖心的表面，然后讀取裂縫在縱向、橫向上所占的線段長度（見圖1中紅線），一般在縱向、橫向上每間隔1 cm讀取數(shù)值，再用縱橫向的和除以整個網(wǎng)格總線條的長度，以定量表征裂縫的原始或有效面縫率。

圖1 網(wǎng)格法裂縫統(tǒng)計示意圖

由于巖心只能代表很小深度范圍內(nèi)裂縫的發(fā)育情況，要掌握全井段的裂縫特征，就要依靠地層微電阻率掃描成像測井資料。成像測井主要應用于巖性和沉積構(gòu)造識別、裂縫描述以及儲集層評價等領(lǐng)域[7-11]。在裂縫精細描述中，成像測井通過測量地層的微電阻率，將采集的數(shù)據(jù)進行一系列的校正、處理后，生成顏色深淺不一的直觀圖像。圖像寬度代表的實際長度為井壁的周長，圖像特征包含巖性、沉積構(gòu)造、孔洞縫等地質(zhì)信息。因此，在成像測井圖上，可以采用網(wǎng)格法對未取心井段進行裂縫參數(shù)統(tǒng)計，但受成像測井資料本身分辨率等因素的影響，需要利用巖心上裂縫參數(shù)統(tǒng)計結(jié)果對其進行校正。具體做法是在取心的深度范圍內(nèi)，將成像測井統(tǒng)計結(jié)果與巖心的定量統(tǒng)計結(jié)果進行對比，可以得出參數(shù)的一個經(jīng)驗差值?X，然后用這一數(shù)值校正該井未取心段的裂縫參數(shù)：

式中 ?X——單井某參數(shù)的經(jīng)驗差值；Xq——取心深度內(nèi)巖心上的參數(shù)統(tǒng)計值；XFq——取心深度內(nèi)成像測井上的參數(shù)統(tǒng)計值；Xw——未取心段的參數(shù)統(tǒng)計值；XFw——未取心段成像測井上的參數(shù)統(tǒng)計值。

在進行裂縫的識別、參數(shù)統(tǒng)計時，一般采用成像測井資料的動態(tài)圖像。裂縫的傾角、傾向、視長度等參數(shù)在成像測井上可以準確識別，不需要校正，而裂縫條數(shù)、張開度、面縫率則需要利用巖心統(tǒng)計結(jié)果進行校正。此方法對單井目的層全部深度范圍內(nèi)的裂縫參數(shù)統(tǒng)計效果較好。

對大北—克深地區(qū)具有巖心描述和成像測井資料的單井進行裂縫參數(shù)定量統(tǒng)計，結(jié)果表明，巴什基奇克組的裂縫視長度多大于 50 mm，巖心裂縫密度為1.0～1.5條/m，各單井巖心上裂縫面縫率差異較大，為0.005～0.050。

以大北202井白堊系的第2筒（5 717.22～5 720.52 m）、第3筒（5 788.00～5 790.85 m）巖心為例加以說明。巖心上識別的裂縫總條數(shù)略大于成像測井上統(tǒng)計的結(jié)果，成像測井上統(tǒng)計的裂縫張開度和面縫率（下文中如不特別說明，面縫率均指有效面縫率）大于巖心的統(tǒng)計結(jié)果（見表1）。裂縫張開度、密度、面縫率的經(jīng)驗差值分別為?1.2 mm、0.2條/m、?0.9%。油田現(xiàn)場巖心裂縫描述資料相對成像測井資料要少得多，可利用大北 202井成像測井資料統(tǒng)計裂縫參數(shù)，然后利用經(jīng)驗差值對各參數(shù)進行校正，即可較準確地求得大北202井成像測井井段的各項裂縫參數(shù)值。

表1 大北202井巖心與成像測井裂縫參數(shù)統(tǒng)計對比表

2 面縫率與裂縫密度的關(guān)系

由大北—克深地區(qū)11口井測試井段裂縫精細描述結(jié)果（見表2）可知，裂縫密度為0.15～2.70條/m；面縫率為 0.4%～6.0%；米采氣指數(shù)為 0～1 081.47 m3/（d?MPa?m）。對裂縫密度（n）與面縫率（ρ）關(guān)系進行回歸分析發(fā)現(xiàn)，二者具較好的相關(guān)性（見圖2），回歸公式為：

表2 大北—克深地區(qū)巴什基奇克組測試段裂縫參數(shù)及產(chǎn)能特征統(tǒng)計表

圖2 裂縫密度與面縫率回歸分析圖

若已知裂縫密度，可由回歸公式計算面縫率。

3 裂縫產(chǎn)能系數(shù)與米采氣指數(shù)的關(guān)系

根據(jù)裂縫精細描述結(jié)果，裂縫特征參數(shù)主要有傾角、傾向、延伸長度、密度及面縫率等，通過對裂縫特征參數(shù)的篩選及與現(xiàn)場試氣結(jié)果進行相關(guān)性分析研究，發(fā)現(xiàn)裂縫密度、面縫率與產(chǎn)能密切相關(guān)[12-13]，面縫率與產(chǎn)能的相關(guān)性最好，面縫率越大，產(chǎn)能越高，因此選擇面縫率為變量構(gòu)造裂縫產(chǎn)能系數(shù)（fq），并且保證裂縫產(chǎn)能系數(shù)為0～1（0表示裂縫不發(fā)育，1為裂縫非常發(fā)育）。篩選大北—克深地區(qū)7口井進行回歸分析（見表3），將米采氣指數(shù)進行歸一化處理（各井米采氣指數(shù)除以米采氣指數(shù)最大值），使處理后的米采氣指數(shù)為0～1：產(chǎn)能最高的井為1，最低的井為0。將其與面縫率進行回歸分析，裂縫產(chǎn)能系數(shù)公式為：

表3 大北—克深地區(qū)巴什基奇克組裂縫特征參數(shù)及米采氣指數(shù)回歸分析表

米采氣指數(shù)（Jsg）與裂縫的發(fā)育程度密切相關(guān)。將米采氣指數(shù)與裂縫產(chǎn)能系數(shù)進行回歸分析（見圖3），建立裂縫產(chǎn)能系數(shù)與米采氣指數(shù)關(guān)系：

若已知裂縫產(chǎn)能系數(shù)，便可利用此系數(shù)預測米采氣指數(shù)，進而預測測試井產(chǎn)能。

大北A井位于新疆阿克蘇地區(qū)拜城縣城西約26.5 km、大宛齊油田北約9.7 km處，是庫車坳陷克拉蘇構(gòu)造帶西段大北1氣田大北201斷背斜東翼上的1口評價井。該井于2010年5月6日開鉆，2011年8月2日完鉆，完鉆井深6 170 m。試油測試井段為5 937～6 058 m，測試厚度 121 m，有效厚度 84.5 m，平均孔隙度8.9%，測試層位為白堊系，巖性為砂巖。成像測井顯示裂縫密度1.2條/m，面縫率2.1%，在10 MPa壓差下利用上述方法預測產(chǎn)量為69.16×104m3/d，實際試油結(jié)果為72.02×104m3/d，說明產(chǎn)能預測結(jié)果與實際試油結(jié)果相符。

圖3 裂縫產(chǎn)能系數(shù)與米采氣指數(shù)回歸分析圖

4 裂縫均勻分布與非均勻分布模型

在裂縫精細定量描述方法中，對某一井段的裂縫密度、面縫率給出的都是井段的平均值，而忽略了各部分裂縫發(fā)育程度的非均質(zhì)性。裂縫均勻分布與非均勻分布對產(chǎn)能的影響目前沒有理論研究結(jié)果可以借鑒?；谏鲜鲇懻摚⒘芽p均勻分布與非均勻分布的理論模型，并分析其對產(chǎn)能的影響。假設測試地層有效厚度為300 m，每100 m為1個單元（見圖4）。模型a表示3個單元的裂縫密度全部為1條/m；模型b表示從上到下3個單元的裂縫密度分別為1.5，0.5，1條/m；模型c表示3個單元的裂縫密度分別為2，0，1條/m；模型d表示3個單元的裂縫密度分別為3，0，0條/m；4個模型中裂縫的總條數(shù)相同。

圖4 裂縫均勻分布與非均勻分布模型

已知裂縫密度，可根據(jù)（3）—（5）式分別求出各單元面縫率、裂縫產(chǎn)能系數(shù)以及產(chǎn)能，然后對各單元產(chǎn)能求和即可得到總產(chǎn)能。假設壓差為10 MPa，模型d裂縫分布最不均勻，產(chǎn)能最大，達447×104m3/d；模型a裂縫均勻分布，產(chǎn)能最小，為156×104m3/d（見圖 5）。

圖5 裂縫均勻分布與非均勻分布產(chǎn)能對比

5 結(jié)論及建議

對于未取心井段，可利用成像測井資料統(tǒng)計裂縫參數(shù)，經(jīng)校正后，可較準確定量描述裂縫延伸長度、傾角、張開度、面縫率及密度等參數(shù)。大北—克深地區(qū)氣井測試段裂縫密度、面縫率與產(chǎn)能有較好的相關(guān)性，據(jù)此得出裂縫產(chǎn)能系數(shù)計算公式，根據(jù)裂縫產(chǎn)能系數(shù)與米產(chǎn)氣指數(shù)的回歸公式即可預測氣井產(chǎn)能，實例分析驗證了方法的有效性。

建立裂縫均勻分布與非均勻分布的理論模型，基于裂縫產(chǎn)能系數(shù)進行產(chǎn)能分析，結(jié)果表明：裂縫性砂巖地層中，裂縫總條數(shù)一定時，裂縫集中分布于某一段時產(chǎn)能更大。

裂縫精細定量描述方法中，裂縫密度的平均化處理直接影響產(chǎn)能預測精度，在現(xiàn)場實際應用時，建議對測試井段進行分段，考慮每一小段的裂縫密度，再利用裂縫產(chǎn)能系數(shù)模型預測產(chǎn)能，可減小誤差。

[1]張榮虎, 姚根順, 壽建峰, 等.沉積、成巖、構(gòu)造一體化孔隙度預測模型[J].石油勘探與開發(fā), 2011, 38(2): 145-151.Zhang Ronghu, Yao Genshun, Shou Jianfeng, et al.An integration porosity forecast model of deposition, diagenesis and structure[J].Petroleum Exploration and Development, 2011, 38(2): 145-151.

[2]曾聯(lián)波.庫車前陸盆地喜馬拉雅運動特征及其油氣地質(zhì)意義[J].石油與天然氣地質(zhì), 2004, 25(2): 175-179.Zeng Lianbo.Characteristics and petroleum geological significance of Himalayan orogeny in Kuqa foreland basin[J].Oil & Gas Geology,2004, 25(2): 175-179.

[3]杜金虎, 王招明, 胡素云, 等.庫車前陸沖斷帶深層大氣區(qū)形成條件與地質(zhì)特征[J].石油勘探與開發(fā), 2012, 39(4): 385-393.Du Jinhu, Wang Zhaoming, Hu Suyun, et al.Formation and geological characteristics of deep giant gas provinces in the Kuqa foreland thrust belt,Tarim Basin, NW China[J].Petroleum Exploration and Development,2012, 39(4): 385-393.

[4]韓登林, 李忠, 壽建峰.背斜構(gòu)造不同部位儲集層物性差異: 以庫車坳陷克拉 2氣田為例[J].石油勘探與開發(fā), 2011, 38(3):282-286.Han Denglin, Li Zhong, Shou Jianfeng.Reservoir property difference between structural positions in the anticline: A case study from Kela-2 gas field in the Kuqa Depression, Tarim Basin, NW China[J].Petroleum Exploration and Development, 2011, 38(3): 282-286.

[5]李軍, 張超謨, 李進福, 等.庫車前陸盆地構(gòu)造壓實作用及其對儲集層的影響[J].石油勘探與開發(fā), 2011, 38(1): 47-51.Li Jun, Zhang Chaomo, Li Jinfu, et al.Tectonic compaction and its influence on reservoirs in the Kuqa foreland basin, Tarim[J].Petroleum Exploration and Development, 2011, 38(1): 47-51.

[6]曾聯(lián)波, 周天偉.塔里木盆地庫車坳陷儲層裂縫分布規(guī)律[J].天然氣工業(yè), 2004, 24(9): 23-25.Zeng Lianbo, Zhou Tianwei.Reservoir fracture distribution law of Kuche Depression in Talimu Basin[J].Natural Gas Industry, 2004,24(9): 23-25.

[7]?；⒘? 胡欣, 徐志強, 等.基巖油氣藏裂縫性儲層的成像測井評價及裂縫預測[J].石油學報, 2010, 31(2): 264-269.Niu Hulin, Hu Xin, Xu Zhiqiang, et al.Evaluation of imaging logging and fracture prediction in fractured basement reservoirs[J].Acta Petrolei Sinica, 2010, 31(2): 264-269.

[8]閆建平, 蔡進功, 趙銘海, 等.電成像測井在砂礫巖體沉積特征研究中的應用[J].石油勘探與開發(fā), 2011, 38(4): 444-451.Yan Jianping, Cai Jingong, Zhao Minghai, et al.Application of electrical image logging in the study of sedimentary characteristics of sandy conglomerates[J].Petroleum Exploration and Development,2011, 38(4): 444-451.

[9]曾焱, 楊麗兵.川西坳陷 DY構(gòu)造須家河組單井裂縫綜合評價[J].成都理工大學學報: 自然科學版, 2010, 37(3): 237-243.Zeng Yan, Yang Libing.Evaluation of single well fractures of Upper Triassic Xujiahe Formation in DY structure of West Sichuan depression, China[J].Journal of Chengdu University of Technology:Science and Technology Edition, 2010, 37(3): 237-243.

[10]景建恩, 梅忠武, 李舟波.塔河油田碳酸鹽巖縫洞型儲層的測井識別與評價方法研究[J].地球物理學進展, 2003, 18(2): 336-341.Jing Jianen, Mei Zhongwu, Li Zhoubo.Investigation of well logging evaluation and identification methods for carbonate reservoirs with fracture and vug in Tahe oil field[J].Progress in Geophysics, 2003,18(2): 336-341.

[11]肖麗, 范曉敏.利用成像測井資料標定常規(guī)測井資料裂隙發(fā)育參數(shù)的方法研究[J].吉林大學學報: 地球科學版, 2003, 33(4):559-563.Xiao Li, Fan Xiaomin.Study on the method of evaluating fracture intensity from conventional logs calibrated imaging logs[J].Journal of Jilin University: Earth Science Edition, 2003, 33(4): 559-563.

[12]甘秀娥.低孔低滲砂泥巖儲層裂縫發(fā)育程度與產(chǎn)能關(guān)系[J].天然氣工業(yè), 2003, 23(5): 41-43.Gan Xiue.Relation between the productivity and the degree of development of the fractures in low porosity and low permeability sandshale reservoirs[J].Natural Gas Industry, 2003, 23(5): 41-43.

[13]禹定孺, 劉春雷.庫車坳陷大北地區(qū)白堊系裂縫及產(chǎn)能特征[J].重慶科技學院學報: 自然科學版, 2011, 13(5): 16-18, 34.Yu Dingru, Liu Chunlei.On the characters of productivity and crack at Dabei area in Kuqa depression of the Cretaceous stage[J].Journal of Chongqing University of Science and Technology: Natural Science Edition, 2011, 13(5): 16-18, 34.