趙中平

西北大學(xué)大陸動力學(xué)國家重點(diǎn)實(shí)驗(yàn)室/地質(zhì)學(xué)系，西安 710069

0 引言

裂縫性油氣藏開發(fā)的關(guān)鍵之一是要盡可能多地鉆遇裂縫。國內(nèi)外學(xué)者研究表明，不同的井軌跡鉆遇裂縫的概率是不同的。Aguilera［1］在對某天然裂縫性儲集層的3口井的研究中發(fā)現(xiàn)：1號井（直井）是一口油氣發(fā)現(xiàn)井，其幸運(yùn)地鉆遇了高角度裂縫；隨后完鉆的2號井也是直井但未鉆遇任何裂縫，進(jìn)行壓裂增產(chǎn)措施后產(chǎn)能很低而被廢棄了；3號井是一口定向井，井軌跡垂直于裂縫走向、鉆遇了多條垂直裂縫而成為生產(chǎn)井。另外，Nelson［2］針對典型非對稱碳酸鹽巖褶皺發(fā)育的不同組系裂縫總結(jié)了相關(guān)的定向井井軌跡的優(yōu)化方案，指出所有的定向井均應(yīng)垂直或近似垂直于裂縫走向，以便鉆遇裂縫的概率最大。由此可見，直井鉆遇高角度裂縫的概率是最低的，而井軌跡垂直于裂縫走向的定向井能夠大大提高裂縫鉆遇概率。

在定向井井軌跡方位一定的情況下，井斜角也會對裂縫鉆遇率有著重要的影響，從而容易誤導(dǎo)人們對儲層裂縫的傾角、密度等特征參數(shù)的統(tǒng)計(jì)，影響對裂縫發(fā)育程度的客觀認(rèn)識。筆者以某裂縫性碳酸鹽巖油藏100余口井（其中90%以上是定向井，均大致垂直于裂縫走向）的全井眼地層微電阻率成像測井（fullbore formation micro image，F(xiàn)MI）資料為例，詳細(xì)分析鉆井井斜角對裂縫傾角、裂縫密度參數(shù)統(tǒng)計(jì)的影響，并討論其意義。

1 井斜角對裂縫傾角統(tǒng)計(jì)的影響

某油田古近系目的層巖性為坳陷期沉積的厚達(dá)500m以上的泥灰?guī)r和粒泥灰?guī)r，在后期區(qū)域擠壓作用下，受基底走滑斷裂控制形成裂縫性儲層。該套地層巖性致密，基質(zhì)幾乎不含油，已發(fā)現(xiàn)的油藏為裂縫性油藏，目前主要利用斜井和水平井進(jìn)行開發(fā)。從下往上劃分為 A、B、C、D、E、F、G、H 共8個小層（圖1）。眾多的FMI資料揭示該套地層主要發(fā)育兩組裂縫，其中：占絕對主導(dǎo)的近南北向裂縫（平均走向方位為NE15°）各層均有發(fā)育；而近東西走向的裂縫（平均走向方位NE105°）主要在F層發(fā)育，其次為G層，其他層則相對不發(fā)育（圖1）。

1.1 未考慮井軌跡影響的裂縫傾角統(tǒng)計(jì)

未考慮井軌跡影響的裂縫傾角統(tǒng)計(jì)結(jié)果見圖1，結(jié)果表明，近南北走向和近東西走向的兩組裂縫的傾角分布特征完全不同。近南北向裂縫傾角總體以高角度裂縫發(fā)育為特征，縱向上從下往上平均裂縫傾角總體逐漸減小，至G和H層裂縫傾角減至最??；底部的A、B層平均裂縫傾角也相對較小，C層平均裂縫傾角最大，達(dá)80.4°；G、H層平均裂縫傾角最小、不足55.0°，是高角度裂縫不發(fā)育造成的。而近東西走向的裂縫從E至H層平均裂縫傾角分布無明顯變化，均為58.0°±1.0°，總體以中低角度斜交裂縫發(fā)育為主，高角度裂縫相對不發(fā)育。

如前所述，該油藏主要利用定向井進(jìn)行開發(fā)，并實(shí)施了部分直井。由于目的層厚達(dá)數(shù)百米，因此部分定向井在上部目的層的井軌跡也是近似垂直的。對各層具有FMI測井資料的鉆井按照井斜角0°～10°、10°～30°、30°～50°、50°～70°、70°～80°和80°～90°的區(qū)間進(jìn)行比例分布的統(tǒng)計(jì)（圖2）表明：鉆揭目的層中部的C-F層的井主要為大斜度井；鉆揭底部的A-B層和上部G-H層的井主要為直井和小斜度井。

各層統(tǒng)計(jì)的平均井斜角與近南北走向裂縫的平均裂縫傾角的變化趨勢一致（圖3a）。在C層鉆井平均井斜角最大，對應(yīng)統(tǒng)計(jì)的近南北走向裂縫傾角的平均值也最大；A-B層和G-H層的平均井斜角較小，對應(yīng)統(tǒng)計(jì)的近南北走向裂縫傾角的平均值也較小。

上述統(tǒng)計(jì)的各層裂縫平均傾角及其縱向的變化趨勢是反映了儲層真實(shí)的裂縫發(fā)育特征呢，還是由于井軌跡影響造成的統(tǒng)計(jì)假象？下面考慮井軌跡的影響、按照不同井斜角范圍統(tǒng)計(jì)裂縫傾角，分析井斜角對裂縫傾角統(tǒng)計(jì)的影響。

圖1 各層兩組裂縫的傾角分布特征Fig.1 Characteristics of fracture orientation and dip histogram in individual layers

1.2 考慮井軌跡影響的裂縫傾角統(tǒng)計(jì)

分別統(tǒng)計(jì)各層小斜度井（包括直井，井斜角范圍0°～20°）和大斜度井（井斜角范圍20°～90°）鉆遇的近南北向裂縫和近東西向裂縫的平均傾角，可以看出小斜度井鉆遇的裂縫平均傾角明顯低于大斜度井鉆遇的裂縫平均傾角（表1、2）。其中：D-F層大斜度井鉆遇的近南北向裂縫平均傾角均高出小斜度井鉆遇的裂縫平均傾角15°以上，E -F層大斜度井鉆遇的近東西向裂縫平均傾角均高出小斜度井鉆遇的裂縫平均傾角10°以上。與前述所有鉆井統(tǒng)計(jì)的各層裂縫平均傾角（圖1）比較，小斜度井鉆遇的裂縫平均傾角偏低，大斜度井鉆遇的裂縫平均傾角明顯偏高。

這一統(tǒng)計(jì)結(jié)果表明，小斜度井鉆遇高角度裂縫的幾率確實(shí)減小了，導(dǎo)致統(tǒng)計(jì)的裂縫平均傾角的減小，不能反映實(shí)際的裂縫傾角和裂縫發(fā)育程度。

2 井斜角對裂縫密度統(tǒng)計(jì)的影響

裂縫密度是描述裂縫發(fā)育程度的重要指標(biāo)，通常指沿直線測量線單位長度內(nèi)觀測到的裂縫數(shù)（線性裂縫密度，簡稱線密度），或單位面積內(nèi)的裂縫累計(jì)長度（面積裂縫密度）［3］，或單位體積內(nèi)的裂縫條數(shù)、長度、裂縫面積（均為體密度），裂縫強(qiáng)度即指單位體積內(nèi)的裂縫面積等。該油藏眾多的FMI成像測井資料為裂縫密度分析提供了良好的資料基礎(chǔ)。

圖2 各層不同井斜傾角范圍的井所占比率Fig.2 Proportion of wells with different inclination range in individual layers

與裂縫傾角統(tǒng)計(jì)一樣，首先未考慮井軌跡影響，利用所有鉆井的FMI資料統(tǒng)計(jì)各層近南北向裂縫密度，結(jié)果表明縱向上各層最大裂縫強(qiáng)度平均值、平均裂縫強(qiáng)度和線密度的變化趨勢與平均井斜角一致（圖3）。在C、D層鉆井平均井斜角最大，對應(yīng)的近南北走向裂縫的裂縫最大強(qiáng)度平均值、平均裂縫強(qiáng)度和線密度均最大；A-B層和G-H層的平均井斜角較小，對應(yīng)的近南北走向裂縫的裂縫最大強(qiáng)度平均值、平均裂縫強(qiáng)度和線密度也較小。

其次考慮井軌跡影響，利用具有FMI資料的大斜度井統(tǒng)計(jì)裂縫密度信息，并與小斜度井資料統(tǒng)計(jì)結(jié)果對比，可分析井斜角對裂縫密度統(tǒng)計(jì)的影響。按照0.5m采樣間隔和1.0m的采樣窗口生成近南北走向裂縫的裂縫強(qiáng)度曲線，并統(tǒng)計(jì)各層段所有單井最大裂縫強(qiáng)度的平均值、平均裂縫強(qiáng)度以及線密度（表3）。結(jié)果表明，總體上根據(jù)大斜度井鉆遇的裂縫統(tǒng)計(jì)的裂縫密度明顯高于小斜度井的統(tǒng)計(jì)結(jié)果。其中A、B和H層大斜度井較少，結(jié)果可能不具統(tǒng)計(jì)意義。

表1 根據(jù)不同斜度鉆井統(tǒng)計(jì)的各層近南北走向裂縫傾角Table 1 Average dip of near NS trending fractures derived from wells with different inclination range in individual layers

表2 根據(jù)不同斜度鉆井統(tǒng)計(jì)的各層近東西走向裂縫傾角Table 2 Average dip of near EW trending fractures derived from wells with different inclination range in individual layers

3 討論

實(shí)例統(tǒng)計(jì)結(jié)果表明，當(dāng)應(yīng)用鉆井資料評價裂縫性儲層時應(yīng)相當(dāng)謹(jǐn)慎，否則由于未鉆遇裂縫很容易漏掉許多有潛力的好的儲集層，正如Aguilera［1］研究中的某裂縫性儲層，如果先鉆2號井，則可能漏掉了這個有潛力的油氣儲集層。在本文某裂縫性碳酸鹽巖儲層中，C-F層是產(chǎn)層（A、B層為水層）且C層產(chǎn)能最高，H和G層幾無產(chǎn)能或產(chǎn)能很低，各層產(chǎn)能差異與未考慮井軌跡影響以及按大斜度井的統(tǒng)計(jì)結(jié)果反映的縱向上裂縫發(fā)育程度差異均一致。然而，考慮到井斜角對裂縫鉆遇率的影響，有必要重新認(rèn)真審視H和G層的潛力。該油藏主要利用定向井開發(fā)，由于目的層厚達(dá)數(shù)百米，因此部分定向井在上部目的層的井軌跡也是近似垂直的。鉆揭上部G-H層的井主要為直井和小斜度井，井斜角大于20°的鉆井比率?。▓D2）；其總體鉆遇高角度裂縫的幾率必然很小。目前H層井斜角大于20°的鉆井僅5口，統(tǒng)計(jì)意義較差；對G層井斜角大于20°的12口鉆井按井斜角20°～40°和40°以上進(jìn)一步分組統(tǒng)計(jì)裂縫密度（表4），結(jié)果表明井斜角40°以上井鉆遇的裂縫密度更大。因此，G-H層實(shí)際裂縫發(fā)育程度可能遠(yuǎn)高于統(tǒng)計(jì)結(jié)果，今后應(yīng)部署水平井等大斜度井進(jìn)一步評價其潛力。

圖3 井斜角與近南北向裂縫統(tǒng)計(jì)的裂縫傾角、強(qiáng)度和線密度的關(guān)系Fig.3 Relationship between the wellbore inclination angle and the average fracture dip and fractured intensity

表3 各層近南北走向裂縫按不同斜度鉆井分組統(tǒng)計(jì)的裂縫強(qiáng)度和線密度Table 3 Density of near NS trending fractures derived from wells with different inclination range in individual layers

表4 G層近南北走向裂縫按大斜度鉆井分組統(tǒng)計(jì)的裂縫強(qiáng)度和線密度Table 4 Density of near NS trending fractures derived from wells with different inclination range in Layer G

同時，上述統(tǒng)計(jì)表明對于類似裂縫性儲層，只有利用井軌跡垂直于裂縫走向和主要裂縫面的鉆井鉆遇的裂縫來分析統(tǒng)計(jì)才能較準(zhǔn)確地統(tǒng)計(jì)裂縫傾角并更客觀地評價儲層裂縫發(fā)育程度。另外，雖然由于地下實(shí)際裂縫資料的不完全性［4－8］而應(yīng)綜合利用各種裂縫資料開展相關(guān)研究，但該實(shí)例統(tǒng)計(jì)結(jié)果也啟示人們需重新審視不加選擇地利用所有測井資料獲得的裂縫信息在地震資料預(yù)測裂縫密度［9－14］或多學(xué)科綜合預(yù)測裂縫密度［15－19］中應(yīng)用的合理性。

4 結(jié)論

1）井斜角對裂縫傾角和裂縫密度的統(tǒng)計(jì)有重要影響。對高角度裂縫而言，直井和小斜度井鉆遇高角度裂縫的幾率明顯低于大斜度井或水平井。因此利用直井和小斜度井統(tǒng)計(jì)的裂縫傾角和密度較實(shí)際值偏小，利用大斜度井和水平井統(tǒng)計(jì)的裂縫傾角和裂縫密度更能較真實(shí)地反映地下儲層裂縫的傾角和裂縫發(fā)育程度。

2）為較準(zhǔn)確地統(tǒng)計(jì)高角度裂縫性儲層的裂縫傾角并更客觀地評價裂縫密度，應(yīng)盡量利用井軌跡垂直于裂縫走向和主要裂縫面的水平井鉆遇的裂縫來分析統(tǒng)計(jì)。對于同類高角度裂縫性油藏的勘探評價和開發(fā)，應(yīng)盡量采用大斜度井和水平井并確保井軌跡垂直于主要裂縫面。

3）對于本文某裂縫性碳酸鹽巖油藏實(shí)例中的G-H層，考慮到井斜角對裂縫密度統(tǒng)計(jì)的影響，有必要重新認(rèn)真審視其潛力，應(yīng)部署水平井等大斜度井進(jìn)一步評價其潛力。

4）應(yīng)用鉆、測井資料評價裂縫發(fā)育程度和勘探開發(fā)潛力時，尤其當(dāng)初期鉆井未鉆遇良好的裂縫性儲層時，應(yīng)充分考慮井軌跡對裂縫鉆遇率的影響，以免漏掉有潛力的良好儲集層。

（References）：

［1］Aguilera R.Naturally Fractured Reservoirs［M］.2nd ed.Tulsa：Pennwell Publishing Company，1995：50－60.

［2］Nelson R A.Geologic Analysis of Naturally Fractured Reservoirs［M］.2nd ed.Houston：Gulf Professional Publishing，2001：332－335.

［3］劉建中，孫慶友，徐國明，等.油氣田儲層裂縫研究［M］.北京：石油工業(yè)出版社，2008.Liu Jianzhong，Sun Qingyou，Xu Guoming，et al.Study on Fracture in Oil－Gas Reservoirs［M］.Beijing：Petroleum Industry Press，2008.

［4］Ortega O J，Marrett R A，Laubach S E.A Scale－Independent Approach to Fracture Intensity and Average Spacing Measurement［J］.AAPG Bulletin，2006，90（2）：193－208.

［5］Narr W.Fracture Density in the Deep Subsurface：Techniques with Application to Point Arguello Oil Field［J］.AAPG Bulletin，1991，75（8）：1300－1312.

［6］Laubach S E.Practical Approaches to Identifying Sealed and Open Fractures［J］.AAPG Bulletin，2003，87（4）：561－579.

［7］Laubach S E，Reed R M，Olson J E，et al.Coevolution of Crack－Seal Texture and Fracture Porosity in Sedimentary Rocks： Cathodoluminescence Observations of Regional Fractures［J］.Journal of Structural Geology，2004，26（5）：967－982.

［8］Barthélémy J F，Guiton M L E，Daniel J M.Estimates of Fracture Density and Uncertainties from Well Data［J］.International Journal of Rock Mechanics and Mining Sciences，2009，46（3）：590－603.

［9］胡光岷，賀振華，黃德濟(jì)，等.利用縱波資料反演裂縫發(fā)育密度和方向［J］.成都理工學(xué)院學(xué)報，2000，27（2）：145－150.Hu Guangmin，He Zhenhua，Huang Deji，et al.Inversion of Fracture Density and Azmuth Using PWave Seismic Data［J］.Journal of Chengdu University of Technology，2000，27（2）：145－150.

［10］賴生華，劉文碧，李德發(fā)，等.利用變參數(shù)回歸方法反演泥質(zhì)巖裂縫密度［J］.巖相古地理，1998，18（3）：50－54.Lai Shenghua，Liu Wenbi，Li Defa，et al.The Fracture Density Model for Pelite Reservoirs Based on the Regression Analysis of Variable Parameters［J］.Sedimentary Facies and Palaeogeography，1998，18（3）：50－54.

［11］李瓊，李勇，李正文，等.基于GA－BP理論的儲層視裂縫密度地震非線性反演方法［J］.地球物理學(xué)進(jìn)展，2006（2）：465－471.Li Qiong，Li Yong，Li Zhengwen，et al.A Seismic Nonlinear Inversion for Apparent Fracture Density of Hydrocarbon Reservoir Based on GA－BP Theory［J］.Progress in Geophysics，2006（2）：465－471.

［12］田鋒.裂縫體的彈性模量和裂縫密度［J］.地質(zhì)學(xué)報，2007，81（10）：1338－1344.Tian Feng.Elasticity Modulus and Fracture Density in the Rocks with Crevices［J］.Acta Geologica Sinica，2007，81（10）：1338－1344.

［13］魏永佩，張之一.油氣田裂縫密度定量研究方法［J］.西安地質(zhì)學(xué)院學(xué)報，1991，13（1）：37－40.Wei Yongpei， Zhang Zhiyi.The Quantitative Analysis Method of the Fracture Density in the Oil ＆Gas Field［J］.Journal of Xi’an College of Geology，1991，13（1）：37－40.

［14］吳朝容，段文燊，李正文.致密碎屑巖裂縫視密度高分辨率非線性反演［J］.西南石油大學(xué)學(xué)報，2007，29（2）：89－91.Wu Chaorong，Duan Wenshen，Li Zhengwen.A High Resolution Nonliear Inversion of Apparent Density of Fracture in Tight Clastic Rock［J］.Journal of Southwest Petroleum University，2007，29（2）：89－91.

［15］魯紅英，肖思和，鄧林.信息融合技術(shù)在裂縫密度綜合預(yù)測中的應(yīng)用［J］.天然氣工業(yè)，2007，27（7）：49－51.Lu Hongying，Xiao Sihe，Deng Lin.Information Fusion Technology Applied to Integrative Prediction of Fracture Intensity［J］.Natural Gas and Industry，2007，27（7）：49－51.

［16］Wong P M.A Novel Technique for Modeling Fracture Intensity：A Case Study from the Pinedale Anticline in Wyoming［J］.AAPG Bulletin，2003，87（11）：1717－1727.

［17］McLennan J A，Allwardt P F，Hennings P H，et al.Multivariate Fracture Intensity Prediction：Application to Oil Mountain Anticline，Wyoming［J］.AAPG Bulletin，2009，93（11）：1585－1595.

［18］Biadu Fred K.Inversion of Fracture Density from Field Seismic Velocities Using Artificial Neural Networks［J］.Geophysics，1998，63（2）：534－545.

［19］王軍，崔紅莊，戴俊生，等.接觸變質(zhì)帶中冷凝收縮縫裂縫參數(shù)定量研究［J］.吉林大學(xué)學(xué)報：地球科學(xué)版，2012，42（1）：58－65.Wang Jun，Cui Hongzhuang，Dai Junsheng，et al.Quantitative Study on Parameters of Condensate Contraction Joint Fracture of in Contact Metamorphic Belt［J］.Journal of Jilin University：Earth Science Edition，2012，42（1）：58－65.