提高蚂蚁追踪裂缝预测精度的应用研究
Application study of improving the precision of the ant-tracking-based fracture prediction technique
责任编辑: 叶佩
收稿日期: 2021-04-26 修回日期: 2021-06-17
基金资助: |
|
Received: 2021-04-26 Revised: 2021-06-17
作者简介 About authors
谢清惠(1985-),女,博士,地质工程师,从事地震解释研究工作。Email:
四川盆地外YS1井区经历多期构造演化,发育的断裂具有多期次、多走向等特征,在已钻井过程中发生了严重的井漏。利用常规的地震解释方法,较难达到精细刻画断裂的要求,也难以捋清断裂期次和展布。针对研究区实际情况,引进蚂蚁追踪技术来进行裂缝预测,并通过降噪、增强连续性、突出断层边界、剔除假构造等手段提高裂缝预测的精度。首先对输入的地震数据体进行构造导向滤波和不连续探测等预处理,之后通过对比优选出蚂蚁追踪参数设置,最后进行产状控制下的蚂蚁追踪,得到一套高精度的蚂蚁追踪体。使用该数据体不但能辅助精细刻画断层展布,还为接下来的断裂期次发育分析和水平井位部署提供数据基础。
关键词:
The YS1 area in Sichuan basin has undergone multiple stages of tectonic evolution.The faults are characterized by multiple stages and multiple strikes,and serious mud losses occurred in the drilled well.It is difficult to finely characterize fracture and analyze structural development when the conventional seismic interpretation method is used.In this study,according to the actual situation,ant tracking technology was used to predict fractures,whose accuracy was improved by noise reduction,continuity enhancement,fault boundaries height,and fake structures elimination.The process is as follows:firstly,the input seismic data volume is pre-processed by structure-oriented filtering and discontinuity detection,then the ant tracking parameter settings are optimized,and finally the ant tracking under the occurrence control is performed.The data can not only be used to assist fault interpretation but also provide a data basis for subsequent fault development analysis and horizontal well location deployment.
Keywords:
本文引用格式
谢清惠, 蒋立伟, 赵春段, 王仲达, 唐协华, 罗瑀峰.
XIE Qing-Hui, JIANG Li-Wei, ZHAO Chun-Duan, WANG Zhong-Da, TANG Xie-Hua, LUO Yu-Feng.
0 引言
目前,随着油气勘探的发展,断层和裂缝预测技术也在不断地革新,针对于大、中型尺度的断层识别方法多且较为成熟,但小尺度的断层和裂缝预测技术仍然是地球物理勘探领域的难点[1,2]。裂缝的地震预测技术包括叠前预测和叠后预测两大类,前者主要是利用快、慢横波的差异或纵波的方位各向异性来预测裂缝的方位、发育程度甚至有效性,该方法要求有昂贵的三维地震采集和处理技术相配套,例如多波多分量地震采集与处理、高密度宽方位地震采集与处理等;后者主要利用叠后三维地震属性对裂缝进行识别和预测,其中很多振幅类、频率类和相位类的属性已被广泛应用,而较为精细的裂缝地震属性分析则主要是围绕地震反射波形的突变(不连续性)来开展,例如相干分析和蚂蚁追踪等边界强化探测技术[3,4,5]。
研究区四川盆地外YS1井区经历了多期构造演化,发育多旋回沉积。由于历次构造运动叠加,区内断裂复杂多样,具有多期次、多尺度、多类型、多走向等特征[13]。目的层存在着可钻性差、井壁稳定性差、井漏、固井难度大等钻完井技术难点,在施工过程中频繁发生井塌、井漏等井壁失稳现象。针对研究区实际情况,利用常规的地震解释方法,较难达到精细刻画断层的要求,也难以捋清断裂期次和展布;故引进蚂蚁追踪技术,但若仅使用常规的蚂蚁追踪流程,存在对中、小裂缝预测精度不够的问题。因此,急需研究一套适用于本区的高精度蚂蚁追踪流程,提高裂缝预测精度,捋清裂缝发育特征并查明井漏原因,为下一步井位部署提供数据基础。本文正是针对研究区的地质特点,结合现有资料条件,围绕如何提高裂缝预测精度进行探讨,最终通过提高分辨率、降噪、增强连续性、突出断层边界、剔除假构造等手段提高精度,并在实际应用中取得良好效果。
1 方法流程简介
常规蚂蚁追踪主要包括以下几个步骤:
1)构造平滑:对原始地震数据体进行构造平滑、降低噪声影响,增强地震有效反射的连续性;
2)不连续性检测:对地震数据不连续点进行探测,并对这一不连续进行强化;
3)蚂蚁追踪参数设置:蚂蚁追踪的效果受控于蚂蚁追踪参数设置,主要有6个参数,包括初始蚂蚁边界、追踪偏差、搜索步长、允许的非法步长、必须合法步长以及搜索中的门槛值等。
4)蚂蚁追踪:得到裂缝预测结果。
针对研究区面临的问题,为了提高裂缝预测精度,依托Petrel软件,笔者探索出一套适用于研究区的蚂蚁追踪技术流程(图1)。首先,对资料进行构造导向滤波,取代了常规蚂蚁追踪第一步中的构造平滑,目的是提高信噪比,增强不连续性;第二步,分别进行了最大曲率、混沌体以及方差体属性计算,对比并分析了3种地震属性的不连续性识别能力,选择最优的不连续性探测方法;第三步,通过对比并优选一套适用于研究区的蚂蚁追踪参数得到初始的蚂蚁追踪结果;与常规蚂蚁追踪对比,为了剔除虚假构造,本文增加了第四步中的产状控制下的蚂蚁追踪,并将不同走向的断层分离开,便于下一步的构造期次分析。
图1
1.1 构造导向滤波
图2
图2
YS1井区构造导向滤波效果
a—原始地震剖面;b—构造导向滤波之后的地震剖面
Fig.2
Structure-oriented filtering maps in YS1 area
a—original seismic profile;b—seismic profile after structure-oriented filtering
1.2 不连续性检测
在构造导向滤波之后,断层的分辨率和清晰度得到了一定的改善,但一些断层依然存在边界不清晰、连续性不强等问题,因此接下来将使用不连续性检测技术来增强断裂的连续性,以突出断层边界。目前常用的方法种类繁多,针对不同研究区的适用性不同,效果也是良莠不齐。其中,较为广泛适用且效果较好的断层边界识别方法有方差体属性、曲率属性、混沌体属性等。方差体技术是一种基于概率方差分析的地层不连续性检测,其通过计算相邻地震道之间的方差来表示各个地震道反射特征的差异,从而达到识别断层的作用。早期的面曲率分析技术主要用于宏观的断裂描述和应力场分析[17,18],近期出来的最大正曲率和最小负曲率大幅提高了曲率研究精度,也常被用于裂缝预测中。混沌属性则是通过对地震数据体中层间地震反射波倾角相似性计算,描述地层内地震反射波的反射特征及反射结构的连续性变化,可有效突出特殊地质体的边界,较好地刻画特殊岩性体和小尺度断层[19,20]。
图3a、3b和3c分别为方差体、最大曲率以及混沌体的沿目的层切片。相比之下,方差体切片仅为较大尺度的NE向断层;最大曲率属性切片比方差体切片多了近SN向断层;而混沌体切片对近SN向断裂刻画得更为清晰。经过对比,认为混沌体属性能更有效地识别YS1井区的断层边界。
图3
图3
YS1井区地震属性沿层切片
a—方差体切片;b—最大曲率切片;c—混沌体切片
Fig.3
Seismic attributes slices in YS1 area
a—variance slice;b—max curvature slice;c—chaos slice
1.3 受参数控制的蚂蚁追踪
蚂蚁追踪的效果受控于蚂蚁追踪参数设置,主要有6个参数,各参数含义和特点见表1。在实际应用中,若对各参数进行分别试验、对比和优选,可能会组合出成千上万种参数组合结果,需要耗费大量的人力和时间。根据试验对比和经验设置,Petrel软件中已为解释人员提供了两种默认的追踪参数组合模式:Aggressive主动模式和Passive被动模式,具体参数设置见表1,二者的区别主要在于初始蚂蚁边界和搜索中的门槛值。主动模式的初始蚂蚁边界较小,而搜索中的门槛值较大,因此,主动蚂蚁追踪算法类似“勤劳的蚂蚁”,更善于挖掘断层,但由于其“主动性”较强,致使噪声较清晰。被动模式的初始蚂蚁边界较大,而搜索中的门槛值较小,因此被动蚂蚁追踪算法即“懒惰的蚂蚁”,其倾向于追踪极强信号,放弃较弱信号,因此有助于压制噪声,体现大断层趋势,但是易导致局部不清晰的低级别小断层呈断续特征,进而破坏其连续性[12]。
表1 蚂蚁追踪方法的参数设置及含义
Table 1
参数 | 含义 | 特点 | 主动模式 | 被动模式 |
---|---|---|---|---|
初始蚂蚁边界 | 定义蚂蚁的初始分布范围 | 越小越利于识别小尺度断层 | 5 | 7 |
追踪偏差 | 设置蚂蚁可搜索方向的偏离范围 | 越大越利于弯曲断层的识别 | 2 | 2 |
搜索步长 | 蚂蚁每次搜索的步长 | 越大搜索能力越大,但可能忽略细节 | 3 | 3 |
允许的非法步长 | 允许搜索超越定义步长的范围 | 越大越连续 | 2 | 1 |
合法步长 | 搜索路径中必须包含的合法步数 | 越小搜索越受限,越不连续 | 2 | 3 |
搜索中的门槛值 | 追踪过程中允许非法步长百分比 | 越大搜索能力越强 | 10 | 5 |
图4a和4b分别为Aggressive主动蚂蚁追踪和Passive被动蚂蚁追踪结果的沿目的层切片。相比之下,被动蚂蚁追踪可以识别大、中型断裂;而主动模式更善于追踪中、小型裂缝,其沿层切片断裂数量大增,甚至对于裂缝细节刻画得更为清晰。因此选取Aggressive主动蚂蚁追踪算法为研究区裂缝预测的优选方法。
图4
图4
YS1井区蚂蚁追踪结果沿层切片
a—被动蚂蚁追踪结果;b—主动蚂蚁追踪结果;c—产状控制下的主动蚂蚁追踪结果
Fig.4
Ant tracking slices in YS1 area
a—passive ant tracking;b—aggressive ant tracking;c—occurrence-controlled aggressive ant tracking
1.4 产状控制的蚂蚁追踪
在优选了蚂蚁追踪的最终参数组合之后,进行初步的蚂蚁追踪。此时的蚂蚁追踪结果搜索到了全方位的裂缝,而实际工区的断裂发育并非在全方位都发育,因此往往存在虚假断裂。产状控制的蚂蚁追踪正是去掉虚假断裂信息的手段。具体做法是:根据生成的地震属性以及初始蚂蚁结果,结合前人研究成果,统计出大断层走向、倾向和倾角等产状参数信息,得到裂缝发育空间展布形态,确定过滤原则。同时考虑层位痕迹一般为低倾角响应,因此需滤去较小的倾角。实际开展蚂蚁追踪时,通过Petrel软件内的赤平投影图来控制蚂蚁追踪的搜索方向(图5):将要过滤掉的断裂倾角和方位选为灰色,则该部分断裂将不会被蚂蚁追踪结果记录到;反之,白色部分则为蚂蚁追踪结果所追踪和记录的。
图5
图5
Petrel内用于产状控制的赤平投影
Fig.5
Stereonet map using by occurrence controlled in Petrel
根据前人研究表明,研究区早期受到近SN向挤压,主要形成NE向断层;晚期受到NE向挤压,主要形成近SN向断层。从图3的几个地震属性中可以看出,研究区发育有NE向和近SN向断层,但初步的蚂蚁追踪结果存在其他走向的断层(图4b),因此需要对此进行虚假构造的剔除,提高裂缝预测精度。结合YS1井区实际地质情况,设定参与计算的断层倾角为45°~90°(区域内发育走滑断层,近乎垂直,选取最大倾角为90°)。对于方位角设置,仅将NE向和近SN向的断裂留下参与计算。具体参数设置见如图4c中的赤平投影图,计算得到的沿层切片见如4c,与图4b蚂蚁追踪切片相比,NW向断裂被剔除,基本只保留了NE向和近SN向断裂,更符合研究区断裂发育的实际情况。但是,此时得到的结果是各时期断裂全部交杂在一起。为了捋清各期次的断裂发育特征和展布规律,可将两个不同走向断裂分离开,便于接下来的断裂期次发育分析。图6a为NE向断裂为主的沿层切片,在YS1井区中南部发育几条平行排列的大断层,北部以短、小断裂发育为主。图6b为近SN向断裂为主的沿层切片,发育中、小尺度裂缝为主,集中在YS1井区西北部。
图6
图6
YS1井区分走向蚂蚁追踪沿层切片
Fig.6
Ant tracking slices with different trends in YS1 area
2 应用效果
图7
图7
YS1井区断层解释多边形与地震属性的叠合
a—与方差体的叠合图;b—与蚂蚁追踪结果的叠合体
Fig.7
Overlapping maps of fault polygons with seismic attributes
a—overlapping map with variance;b—overlapping map with ant tracking
图8为过YS1-3水平井段拉取的剖面,图上可见水平井段上的蚂蚁追踪结果与泥浆漏失点对应关系良好,说明了裂缝发育是泥浆漏失的主要原因,接下来的水平井部署应避开裂缝发育带。
图8
图8
过YS1-3井拉取的剖面
a—原始地震剖面;b—蚂蚁追踪结果
Fig.8
Seismic profile crossing YS1-3
a—original seismic profile;b—ant tracking
3 结论
1)不同于常规的蚂蚁追踪流程,本文方法用构造导向滤波取代了普通的构造平滑,意在增强断层不连续性。除此之外,还增加了产状控制下的蚂蚁追踪,不仅从整体上对追踪结果进行产状控制、剔除虚假构造,还将不同走向的断层分离开进行计算。该方法尤其适用于发育多期次、多走向的复杂构造地区。
2)该方法通过降低噪声、增强连续性、突出断层边界、剔除虚假构造等方面提高裂缝预测精度,使得裂缝刻画更加精细,断层解释更为准确,并且定位了泥浆漏失的主要原因。这不仅为地震解释提供参考,也能为构造期次分析提供依据,还能为接下来的水平井部署提供数据基础,减小钻井风险。
3)蚂蚁追踪结果受控于蚂蚁追踪参数设置,但参数数量甚多,不同的解释员进行不同的参数设置会产生成千上万个结果,因此需谨慎选择参数。软件中常设置有主动蚂蚁追踪和被动追踪两种不同的参数组合模式,主动追踪善于挖掘裂缝,大、小裂缝信息都丰富;被动模式善于追踪较强信号,主要预测大断层。主动追踪适用于精细裂缝预测,被动追踪可用于预测大断层以辅助断层解释。
参考文献
不同尺度裂缝的叠后地震预测技术研究
[J].
Poststack seismic prediction techniques for fractures of different scales
[J].
塔中东部走滑断裂带特征及油气地质意义
[J].
Characteristics of strike-slip fault and its hydrocarbon geological significance in the eastern of central Tarism Basin
[J].
地震裂缝预测技术研究进展
[J].
Progress of seismic fracture characterization technology
[J].
页岩储层地震尺度断裂系统分析及其石油地质意义——以四川盆地长宁地区宁201井区为例
[J].
Analysis of seismic-scale fracture system of shale reservoir and its petroleum significance:A case study of well Ning 201 area of Changning Block, Sichuan Basin
[J].
改进的蚂蚁追踪裂缝检测算法及其应用研究
[J].
Improved ant tracking crack detection algorithm and its application
[J].
Ant system:Optimization by a colony of cooperating agents
[J].
蚂蚁追踪技术在裂缝预测中的应用——以青西油田下沟组为例
[J].
Application of ant tracking technology to fracture prediction:A case study from Xiagou Formation in Qingxi Oilfield
[J].
基于蚂蚁体各向异性的裂缝表征方法
[J].
Fracture characterization based on azimuthal anisotropy of ant-tracking attribute volumes
[J].
Fault and fracture distribution within a tight-gas sandstone reservoir:Mesaverde Group,Mamm Creek Field,Piceance Basin,Colorado,USA
[J].DOI:10.1144/petgeo2011-093 URL [本文引用: 1]
Integration of surface/subsurface techniques reveals faults in Gulf of Suez oilfields
[J].DOI:10.1144/1354-079310-014 URL [本文引用: 1]
基于蚂蚁算法的断裂追踪技术在乌夏地区的应用
[J].
Application of fault tracking technology based on ant colony algorithm in Wuxia Area
[J].
利用高清蚂蚁体精细解释复杂断裂带
[J].
3D seismic fine structural interpretation in complex fault zones based on the high-definition ant-tracking attribute volume
[J].
中扬子地区晚奥陶世—早志留世古地理演化及页岩沉积模式
[J].
Paleogeographic evolution of the late Ordovician-early Silurian in upper and middle Yangtze regions and depositional model of shale
[J].
基于构造导向滤波的高斯束层析速度建模方法及其应用
[J].
Gaussian beam tomography with structure-filtering and its applications
[J].
利用构造导向滤波技术识别复杂断块圈闭
[J].
Complex fault-block traps identification with structure-oriented filter
[J].
Fast Structural interpretation with structure-oriented filtering
[J].
构造曲率与振幅曲率在地震资料解释中的应用
[J].
The application structural curvature and amplitude curvature attribute to seismic interpretation
[J].
地震属性在微断层解释中的应用
[J].
Application of seismic attributes in micro-fault interpretation
[J].
混沌属性预测泌阳凹陷陡坡带小型砂砾岩体
[J].
Application of chaos attributes to predict the small-scale glutenite bodies in the steep-slope zone of Biyang depression
[J].
/
〈 |
|
〉 |
