LZ地区致密砂岩储层裂缝综合预测方法及应用

文章编号：1673-26（2011）03-0097-05

岩性油气藏

LITHOLOGICRESERVOIRS

Vol.23No.3Jun.2011

LZ地区致密砂岩储层裂缝综合预测方法及应用

王志萍1，秦启荣1，苏培东1，范晓丽2

（1．西南石油大学资源与环境学院；2.中国石化中原油田勘探开发科学研究院）

摘

要：LZ地区储层为低孔、低渗致密砂岩储层，裂缝控制了储层的储渗空间和井的产能。通过对LZ地

区的构造特征及演化分析，并结合野外露头资料对裂缝产状进行分期配套，认为该区主要发育横张缝、剪切缝以及断层伴生缝和派生缝等构造成因裂缝。针对以上3种构造裂缝类型，分别采用构造曲率法、古构造应力场有限元模拟法、地震不连续性检测法等对该区不同类型裂缝的分布进行了预测，并采用权重评价方法综合这3种预测成果进行裂缝的综合预测，即建立各预测方法的准确率与其影响因子之间的回归函数，再根据预测方法的准确率确定权重系数，将不同方法的预测成果进行综合权重计算，从而对研究区致密储层的裂缝发育情况进行综合预测，经钻井资料证实其预测效果较好。关键词：裂缝综合预测；构造曲率法；古构造应力场；地震不连续性检测；四川盆地中图分类号：TE122.2+21

文献标识码：A

布作出合理的评价。

储层裂缝预测是石油地质界的一个难题，至今还没有一套行之有效的技术方法来解决该问题。目前，裂缝预测方法多种多样，地质方法主要应用岩心裂缝观测描述［1］和野外露头观测描述等；测井方法主要应用裂缝的非

均质性在测井曲线上形成的异常（电阻率、密度、声波波形、地层倾角）等来实现［2］；地震方法主要应用多场信息预测技术、多方位地震属性预测裂缝技术、多分量与各向异性检测技术［3］、裂缝边缘检测技术及裂缝非线性预测技术等；数学方法主要应用构造应力场反演法［4-11］、有限单元法、分形几何法［12］、灰色评判法、统计法、神经网络法等。但由于裂缝发育受到岩性、构造应力场、岩层厚度等多种因素的影响，非均质性强，因此每种预测方法均有其局限性。针对LZ地区致密砂岩储层中裂缝特点，笔者提出了一套在地质分析基础上采用构造曲率法、古构造应力场有限元模拟法及地震不连续性检测法等3种裂缝预测方法，将不同方法的预测成果综合在一起，进行权重值计算，从而对各区块裂缝的发育分

收稿日期：2011-01-17；修回日期：2011-02-20

基金项目：国家科技重大专项（编号：2008ZX05001）资助。

1区域概况

LZ地区位于四川盆地川中隆起北部斜坡带，

是川中隆起向川北坳陷延伸的一个低幅度鼻状隆起带的一部分，与梓潼—苍溪鼻状构造带和仪陇—平昌低缓构造带毗邻（图1）。LZ地区构造主要由呈弧状弯曲的挤压性褶皱和旋扭性褶皱组成。中侏罗统千佛崖组储层埋深达2700～3000m，岩性以

0

龙山断褶带

剑阁

3060km

大巴山弧形褶皱带陷造构门

广元

旺苍

东西向构造带

南江

北川

东向北

构造

巴中

带

西坳通江

北梓

潼

苍

溪

构鼻状

造带

阆中

仪陇—平昌平缓构造带

东褶

达县

带苍溪

向川

南部

盐亭

川中隆起

蓬安

皱带开江

研究区褶皱地名构造分界线

LZ地区区域构造位置图

Fig.1RegionalstructurallocationofLZarea

图1

第一作者简介：王志萍，1981年生，女，西南石油大学在读博士研究生，研究方向主要为储层裂缝特征及预测研究。地址：（610500）四川省成

都市西南石油大学研究生部博09级。E-mail：wzp1209@163.com

岩性油气藏第23卷第3期

细粒石英砂岩、细粒长石岩屑砂岩、含钙粗粉砂岩为主，次为介壳砂质灰岩，孔隙度为0.7%～2.5%，渗透率为0.008～0.1mD，为典型的致密砂岩储层。其储集空间主要为孔隙和裂缝，其中裂缝为油气提供了储集空间，提高了储层的渗流能力。在裂缝发育带渗透率可达60mD，这对于LZ地区致密砂岩储层获得高产具有重要意义。LZ地区裂缝类型多样，既有构造缝，也有成岩缝，其中以构造缝为主，其倾角主要为15°～45°，裂缝充填程度不高，未充填缝占一半以上。

频率/%Fig.2

25201510500～10LZ地区裂缝走向玫瑰花图

RosediagramoffracturestrikeinLZarea

图2

2裂缝成因分析

裂缝成因分析是在研究区构造特征及演化研

究的基础上，以裂缝构造成因理论模型为指导，利用露头资料对裂缝进行分期配套

［13］

，最终确定裂缝

的形成期次和成因类型。裂缝构造成因分析是储层裂缝综合预测的基础，为合理采用不同的裂缝预测方法提供可靠的依据［14］。

8090304060705020～～～～～～～～0000000078236541图3

LZ地区发育北东向、北西向和近东西向的褶

皱。北东向褶皱是在印支期北东东向隆起背景上发育起来的；燕山运动中后期，在南北向构造挤压作用下，LZ地区形成了一系列近东西向的褶皱；喜山运动期，构造应力的方向由近南北向转变为南偏西向，又形成了一系列北西向的褶皱。所以，LZ地区褶皱主要是在燕山后期和喜山期构造应力作用下形成的。

根据构造成因理论分析认为，伴随2期构造应力形成的裂缝系统主要与挤压褶皱作用和断层破裂作用有关。在燕山中后期南北向挤压作用下形成的裂缝，主要为北东向与北西向平面“X”型剪切缝和近东西向的剖面“X”型剪切缝；在喜山期北东向挤压力作用下形成的裂缝，主要为近南北向和近东西向的平面“X”型高角度剪切缝及北西向剖面“X”型剪切缝；在喜山期北东向应力的挤压作用下，燕山期形成的北东向剪切缝开启并形成横张缝且规模变大；燕山期北西向裂缝和喜山期北东向裂缝在构造应力的持续作用下发育为断裂。

根据露头资料对裂缝产状的分析表明，LZ地区裂缝走向以北西向为主，北东向、近南北向、近东西向裂缝次之（图2），裂缝以中—低角度缝和垂直缝为主（图3）。构造成因分析与裂缝类型进行对比认为：研究区发育燕山期北西向平面“X”型剪切缝、喜山期北西向剖面“X”型剪切缝、断裂的派生缝和伴

Fig.3

LZ地区野外裂缝倾角统计

StatisticsoffracturedipinLZarea

生缝，横张缝局部地区发育，其它类型裂缝发育程度较低。

3裂缝预测

在裂缝成因类型分析基础上，针对性地选择构

造曲率法、古构造应力场有限元数值模拟法以及地震不连续性检测法来综合预测裂缝的分布。构造曲率法主要预测横张缝分布；古构造应力场有限元数值模拟法主要预测剪切缝的分布；地震不连续性检测法主要预测断层伴生缝及派生缝的分布。

3.1构造曲率法

由于曲率值可以反映弯曲岩层面上弯曲派生

的抗张应力形成的张性缝，因此可利用构造曲率法来预测张性裂缝的分布。本次主要采用了趋势面拟合法来计算千佛崖组储层的相对曲率值，预测结果（图4a）中，Ⅰ级裂缝发育区曲率值为22～31，裂缝发育程度好；Ⅱ级裂缝发育区曲率值为13～21，裂缝发育程度中等；Ⅲ级裂缝发育区曲率值为1～13，裂缝发育程度差。千佛崖组储层曲率值总体较小且分布不均匀，大部分地区曲率值小于4，说明张性缝在大部分地区不发育，仅在局部地区发育，如石龙场构造区。这与研究区整体构造变形强度不太相符。

3.2古构造应力场有限元数值模拟法

古构造应力场有限元数值模拟法是在分析构

2011年王志萍等：LZ地区致密砂岩储层裂缝综合预测方法及应用99

造特征和演化史基础上建立地质模型，综合考虑储层的岩性、厚度、岩石力学参数等，以有限元法建立数学模型，并以摩尔-库伦准则为依据，计算岩石破坏接近程度，从而判断储层的岩体破裂程度，同时结合岩心资料和测试资料等对裂缝分布进行预测。在模拟中，用岩石破坏接近程度η值来表示岩体破裂的相对程度，也就是裂缝的发育程度。一般说来，

3.3基于局部构造熵的地震不连续性检测法基于局部构造熵的地震不连续性检测法分析

技术，是运用相关性原理突出相邻道之间地震信号的不连续性，进而达到检测断层和反映地质异常特征展布为目的［15］。具有相同反射特征（振幅、频率、相位）的区域呈相似性，相关值接近1，而在断层、岩性变化的突变点则呈非相似性，相关值接近0。根据相关值高低的空间变化，能快速识别出断层与裂缝的发育带。因此，地震道之间信号相关值的大小是判别断层或裂缝发育带存在的一个重要标志。对本次地震偏移资料进行不连续性处理，并针对处理后的数据体提取千佛崖组时间范围内的平均值，以检测层段内部地震信号的不连续性，从而进行裂缝预测。该方法对断层以及断层附近的裂缝检测效果较好。图5为千佛崖组不连续性平均值平面分布图，红色表示不连续性高值，为Ⅰ级裂缝发育区，裂缝发育程度好；黄色表示不连续性中值，为Ⅱ级裂缝发育区，裂缝发育程度中等；蓝—绿色表示不连续性低值，为Ⅲ级裂缝发育区，裂缝发育程度最差。从图中可以看出不连续性高值主要分布在工区中部，表明断层与裂缝发育，地震不连续性检测法的预测结果与古构造应力场有限元数值模拟法预测

η值越大，裂缝就越发育。目前，国内外通过古构造

应力场来研究储层裂缝发育程度还没有统一的判别标准，所以，η值的选取主要是根据多年来在四川盆地的研究结果，并结合LZ地区构造特征和已有的生产资料来确定。预测成果图（图4b）中，红色区域为Ⅰ级裂缝发育区，η值为1.203～1.303，岩体破裂发育，裂缝发育程度好；黄色区域为Ⅱ级裂缝发育区，η值为1.125～1.203，岩体破裂较发育，裂缝发育程度中等；浅黄色区域为Ⅲ级裂缝发育区，η值为1.022～1.125，为岩体破裂发育临界区，裂缝发育程度差。与生产井的资料进行对比可知，高产井一般位于预测结果的Ⅰ级和Ⅱ级区，而一般的生产井位于Ⅱ级和Ⅲ级区，Ⅲ级区主要是一些中低产井或干井。

0

10

20km

Cs47

Cs55石龙场Sl8

Ss1Cs42

Sl6Cs44Sl7

Sl18

曲率

28252219161310741

结果一致。

Cs60

金星场Cf50

Sl20

Cs73Sl12

Sl3

宝马Sl19

Sl16

Sl1Sl10

Sl9Sl2

Sl11Sl13老鸦Cy46

Cb54柏垭

高

Cs55

Cf50

中

Ss1

Sl19Sl18

Sl16

Sl3Sl9Sl1

（a）

-2650055-2-2650-2550Cs47

石龙场Cs55Ss1-2750金星场Cs73Cf50

Sl20

Cs60Sl3

-255050-26Sl8Cs42Sl6Cs44Sl7

Sl18

宝马Sl19

Sl16

-2650085-2-2850-2750-2650-265001020km

Cb54

柏垭

Sl12Sl2

Sl11Sl10

Sl1

低

老鸦

Sl13

Cy46

-2650-255050-26图5千佛崖组不连续性检测平均值平面分布图

-2450-2450-2450Fig.5PlanedistributionofaveragevalueofseismicdiscontinuitydetectionofQianfoyaFormation

20km

010

（b）

Ⅰ级裂缝

发育区断层

50-2450-234

4.1

裂缝综合预测及效果分析

裂缝综合预测

储层裂缝综合预测主要由综合构造曲率法、古

Ⅱ级裂缝

发育区等高线

Ⅲ级裂缝

发育区

工业油气井

Ⅳ级裂缝

发育区

非工业油气井

图4千佛崖组裂缝预测平面图

构造应力场有限元数值模拟法以及地震不连续性检测法这3种预测方法的成果得来的。基于构造曲率法原理，其预测的裂缝主要为张性缝；地震不连

Fig.4FracturedistributionpredictionofQianfoyaFormation

（a）千佛崖组顶曲率分布图；（b）古构造应力场有限元数值

模拟预测裂缝分布图

100岩性油气藏第23卷第3期

续性检测方法由于地震分辨率低，预测的裂缝多为规模相对较大的裂缝或小断层，而对规模较小的裂缝响应较差；古构造应力场有限元数值模拟法受模拟过程中岩石力学参数取值、边界条件等诸多因素的影响，主要预测的裂缝为剪切缝。因此各种方法均有一定的局限性。根据各预测成果的准确程度，设定不同的权重系数，对不同预测成果进行权重分析，从而对各区块裂缝的发育分布作出合理的评价。数学模型为：y=a1X1+a2X2+a3X3。式中，y可看作是总体测评目标；X1，X2，X3代表测评要素；a1，a2，a3为各个测评要素的权重系数。

其中以各井区的裂缝发育程度为总体测评目标，以3种方法的预测成果为测评要素，以3个预测成果的准确程度为权重系数。为了满足数学模型的计算要求，需将定性分析结论进行量化以便定量计算，在预测成果中按裂缝发育程度分为好、中、差3类，分别赋予相应的权重，权重值越大，表示其裂缝发育程度越好。建立各预测方法的准确率与其影响因子之间的回归函数，从而计算出各预测方法的准确率，再将其值进行归一化处理后即得到权重系数值。曲率法预测准确率较低，权重系数值为

育程度高的区域，具有构造褶皱强、曲率高、邻近断层、破坏强度大、不连续性高等特点；权重得分在

1.73～2.3的区域为Ⅱ级裂缝发育区，这些区域其多

种方法预测结果不同，显示裂缝发育程度高或中等，具有构造褶皱强、曲率高、距断层近、破坏强度较小等特点（图6）。

-27500

25-2610

Cs47

20km

-2750-275075-28-2625050-2Cs60

-2500Cs73

金星场

Cf50Sl17

Sl20

Sl8Cs41

Ss1Sl22

Cs55Cs42Cs30

石龙场Sl6Cs43Sl21Sl7

Cb54Sl3

宝马

Cs44

Sl19

Sl16

Cs48

-2750-2625Sl18

-2375柏垭Sl1Sl2Sl10Sl12Sl15

Sl4

Sl13

老鸦

Sl11

Cy46

-2625050-2-2500-2375-2250-2250断层

工业油气井非工业油气井裂缝发育区裂缝较发育区

图6千佛崖组裂缝综合预测图

Fig.6Syntheticalfracturepredication

ofQianfoyaFormation

15%；古构造应力场有限元数值模拟法预测和地震

不连续性检测法预测的准确率较高，权重系数值分别为43%和42%。各种裂缝预测方法的测评要素值乘以其相应的权重系数，相加即得裂缝发育程度权重得分，得分越高，裂缝发育程度越好。计算结果如表1。

表1

名称

4.2预测效果分析

LZ地区主要经历了2次构造运动，构造运动

的方式和强度比较容易确定，且有野外露头验证，构造分析的结果可靠性很高，这为构造应力场地质模型的建立提供了可靠的边界条件。因此，古构造应力场有限元模拟法和地震不连续性检测法的结果相对比较可靠，LZ地区构造主要为挤压型褶皱，易形成张性缝，用曲率法来预测张性缝是可行的。裂缝综合预测成果与目前掌握的钻井、测井和岩心资料基本相符（表2）。如裂缝发育区Sl10井千二段

裂缝发育程度权重评价表

权重得分

评价

Table1Weightevaluationoffracturedevelopment

Sl7井区CS42井区CS44井区Sl19井区CS55井区Cf50井区Sl20井区Sl2井区Sl12井区Cb54井区Cy46井区

2.452.582.182.152.372.412.152.081.852.222.25

ⅠⅠⅠⅡⅠⅠⅡⅡⅡⅡⅡ

2751～2755.50m放空0.25m，说明地层中存在较

大的缝、洞。

表2

千佛崖组裂缝发育情况统计表

Table2Statisticsoffracturedevelopment

ofQianfoyaFormation

井名

井段/m

裂缝信息

岩心裂缝密度为14.29条/m岩心裂缝密度为2.23条/m

放空0.25m漏失泥浆26.7m3

岩屑中见自形—半自形石英颗粒，测井解释发育3m厚的裂

缝性油气层

Cs47井Sl11井Sl10井Sl13井Sl17井

2915.35～2918.422918.43～2918.96

2838.38～2844.012751.00～2755.50

—

权重得分在2.3～3的区域为Ⅰ级裂缝发育区，这些区域其多种方法预测的结果相同，皆为裂缝发

2857.00～2858.502908.50～2910.00

2011年

王志萍等：LZ地区致密砂岩储层裂缝综合预测方法及应用

［5］

101

5结论

（1）LZ地区主要发育燕山期北西向平面“X”型

王国刚，苏培东，秦启荣.罗家寨构造嘉陵江组嘉四2段底部裂缝预测［J］．岩性油气藏，2009，21（4）：82-86.

［6］杨凯，郭肖.裂缝性低渗透油藏三维两相黑油数值模拟研究［J］．岩性油气藏，2009，21（3）：118-121

剪切缝、喜山期北西向剖面“X”型剪切缝以及断裂的派生缝和伴生缝，其它类型裂缝发育程度较低。

（2）在裂缝构造成因分析的基础上，针对性地选用了构造曲率法、古构造应力场有限元模拟法与地震不连续性检测法，分别对裂缝进行了预测，再根据预测方法的准确率确定权重系数，将不同方法的预测成果进行综合权重计算，从而对LZ地区侏罗系千佛崖组致密储层的裂缝发育情况进行了综合预测，已钻井资料表明预测效果较好。

参考文献：

［1］

［10］［7］

王孟华，崔永谦，张锐峰，等.泥灰岩裂缝储层预测方法研究———以束鹿凹陷为例［J］．岩性油气藏，2007，19（3）：110-119.

［8］葛善良，周文，谢润成，等.复杂构造变形地区致密岩石力学性质实验分析［J］．岩性油气藏，2010，22（1）：59-.

［9］MaertenL，MaertenF.Chronologicmodelingoffaultedandfracturedreservoirsusinggeomechanicallybasedrestoration:Techniqueandindustryapplications［J］.AAPGBulletin，2006，90（8）：1201-1226.

陈波，田崇鲁.储层构造裂缝数值模拟技术的应用实例［J］.石油学报，1998，19（4）：50-54.

［11］李虎，秦启荣，王志萍，等.大邑构造须二段储层裂缝预测［J］.岩性油气藏，2010，22（2）：107-115.

AmeenMS，BuhidmaIM，RahimZ.Thefunctionoffracturesandin-situstressesintheKhuffreservoirperformance，onshorefields，SaudiArabia［J］.AAPGBulletin，2010，94（1）：27-60.

［12］

黄捍东，魏修成，叶连池，等．分形边缘检测在裂缝预测中的应用［J］.石油地球物理勘探，2002，37（1）：65-68.

［2］［3］［4］

许同海.致密储层裂缝识别的测井方法及研究进展［J］.油气地质与采收率，2005，12（3）：75-78.

杨克明，张虹．地震三维三分量技术在致密砂岩裂缝预测中的应用［J］．石油与天然气地质，2008，29（5）：683-6.

周新桂，邓宏文，操成杰．储层构造裂缝定量预测研究及评价技术［J］．地球学报，2003，24（2）：175-180.

［13］徐开礼，朱志澄．构造地质学［M］.北京：地质出版社，1998：

111-118.［14］

童亨茂．储层裂缝综合预测方法在GBEIBE油田的应用［J］．中国石油勘探，2007，（3）：77-80.［15］

周艳辉，高静怀．局部结构熵算法在地震数据不连续性检测中的应用［J］．煤田地质与勘探，2007，35（1）：71-73.

Predictionmethodoffractureintightsandstonereservoir

anditsapplicationinLZarea

WANGZhi-ping1，QinQi-rong1，SUPei-dong1，FANXiao-li2

（1.SchoolofResourcesandEnvironment，SouthwestPetroleumUniversity，Chengdu610500，China；

2.ResearchInstituteofExplorationandDevelopment，ZhongyuanOilfieldCompany，Sinopec，Puyang457001，China）

Abstract：ThereservoirinLZareaistightsandstonereservoir,andfracturescontrolthereservoirspacesandwellcapacity.Basedontheanalysisofstructuralfeaturesandevolution,combinedwiththefracturedividedperiodcoordinationonthebasisofoutcropdata,thetectonicoriginoffracturesisanalyzed.Itisconsideredthatthefracturesdevelopedinthestudyareamainlyareshearfracture,tensionalfracture,associatedandderivedfractureswithfaults.Theconstructionalcurvaturemethod,paleotectonicstressfieldfiniteelementmodelingmethodandseismicdiscontinuitydetectionmethodareappliedtopredictthedifferentfracturesrespectively.Theweightisadoptedtoevaluatethepredictionresults.Itistoestablishtheregressionfunctionbetweentheaccuracyofthepredictionmethodanditsinfluencingfactors,andthendefinetheweightcoefficientthroughthepredictionaccuracyandcalculatetheweightsofthedifferentpredictionresults.Finallyitcanbeusedtopredictthefracturedevelopmentintightsandstonereservoirinthestudyarea.Thewelldrillingdataprovethegoodpredictioneffects.

Keywords：syntheticalfracturepredication；constructionalcurvaturemethod；paleotectonicstressfield；seismicdiscontinuitydetection；SichuanBasin

（本文编辑：杨琦）