E-mail Alert Rss
 

物探与化探, 2023, 47(3): 747-756 doi: 10.11720/wtyht.2023.2579

方法研究·信息处理·仪器研制

顺北地区碳酸盐岩断控缝洞体油气产能定量化估算技术

刘军,1,2, 黄超2, 杨林2, 张永升2, 查明1

1.中国石油大学(华东) 地球科学与技术学院,山东 青岛 266580

2.中石化西北油田分公司 勘探开发研究院,新疆 乌鲁木齐 830011

Quantitative prediction technology for the hydrocarbon production capacity of fractures and vugs in fault-controlled carbonate reservoirs in the Shunbei area

LIU Jun,1,2, HUANG Chao2, YANG Lin2, ZHANG Yong-Sheng2, ZHA Ming1

1. School of Geosciences,China University of Petroleum,Qingdao 266580,China

2. Exploration and Development Research Institute,Sinopec Northwest China Petroleum Bureau,Urumqi 830011,China

第一作者: 刘军(1982-),男,2005年获中国石油大学(华东)信息与计算科学专业学士学位,2009年获该校地球探测与信息技术专业硕士学位,目前该校博士在读,主要从事石油物探综合研究工作。Email:408842381@qq.com

责任编辑: 叶佩

收稿日期: 2022-02-22   修回日期: 2023-01-31  

基金资助: 中国石化重点科技项目“超深层碳酸盐岩地震关键技术研究与应用”(P21071)

Received: 2022-02-22   Revised: 2023-01-31  

摘要

走滑断裂带及其控制的缝洞体是顺北地区奥陶系新的勘探领域,探索断控缝洞体储层单井产量多属性量化预测技术,为顺北地区断裂带井位部署与井轨迹优选提供科学依据和技术支撑。断控缝洞体储层中,断裂是油气运移通道,其纵向延伸和宽度是断裂的通源性,主干断裂与次级断裂平面交汇结构表征断裂的汇聚性,通源性和汇聚性表征成藏断裂特征。断控缝洞体中缝洞是储集空间的主体,断控缝洞体储层规模和断裂裂缝分布密度表征缝洞特征。断控缝洞体储层的单井产量既与断裂特征有关,更与缝洞特征有关。因此,量化的断裂特征和缝洞特征是单井产量量化预测的基础。地震数据及其多种属性的断层精细解释,量化后得到断裂特征值。结构张量属性圈定断控缝洞体储层范围,异常体振幅属性表征缝洞的分布及其体量,最大似然概率属性表征断裂裂缝分布密度,3种属性融合得到缝洞特征值。最后将量化的断裂特征值和缝洞特征值融合为断控缝洞体特征值,与实测钻井的年产液量进行统计分析,得到年产液量与断控缝洞体特征值的统计关系式,实现断控缝洞体储层单井产量的量化预测。

关键词: 断控缝洞体储层; 断裂特征; 缝洞特征; 断控缝洞体特征; 单井产量量化预测; 异常体分解; 结构张量; 最大似然概率

Abstract

The strike-slip fault system and the fractures and vugs in fault-controlled reservoirs have become new Ordovician exploration targets in the Shunbei area.This study explored a multi-attribute quantitative prediction technique for the single-well production of fault-controlled reservoirs,aiming to provide a scientific basis and technical support for well deployment and well trajectory optimization in the fault zones of the Shunbei area.Faults serve as pathways for hydrocarbon migration in fault-controlled reservoirs.Their connectivity to provenance areas can be characterized by their longitudinal extension and the width of the fault zone,and their convergence can be characterized by their lateral extension and the planar intersection structure between major and secondary faults.Both characteristics of faults were used to indicate the fault characteristics in the oil source conditions for oil accumulation.The space for hydrocarbon accumulation in fault-controlled reservoirs is dominated by dissolution fractures and vugs,such as caves.The fracture-vug characteristics,which represent the space for hydrocarbon accumulation and its connectivity,were characterized by the scale of fault-controlled reservoirs,the volume of the fractures and vugs,and the density of fractures in faults.The single-well production of fault-controlled reservoirs is related to both the fault characteristics and the fracture-vug characteristics.Therefore,the quantification of fault characteristics and the fracture-vug characteristics is the basis for the quantitative prediction of single-well production.The fault characteristic values were determined through the fine-scale interpretation of faults based on seismic data and their multiple attributes.Moreover,the range of fault-controlled reservoirs was delineated based on structure tensors;the distribution and volume of caves,fractures,and vugs were characterized using the amplitude of anomalous bodies;and the density of fractures in faults can be characterized using the maximum likelihood probability.Then,the three attributes were fused to determine the fracture-vug characteristic values.Finally,the quantified fault characteristic values and fractured-vug characteristic values were fused into the characteristic value of fault-controlled reservoirs.Through the statistical analysis of the characteristic value and annual liquid production of drilled wells,the statistical relationship between the annual fluid production and the characteristic value of fault-controlled reservoirs was determined,thus achieving the quantitative prediction of the single-well production of fault-controlled reservoirs.

Keywords: fault-controlled reservoirs; fault characteristic; fracture-vug characteristic; characteristic of a fault-controlled reservoirs; quantitative prediction of single-well production; anomalous body decomposition; structure tensor; maximum likelihood probability

PDF (6046KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘军, 黄超, 杨林, 张永升, 查明. 顺北地区碳酸盐岩断控缝洞体油气产能定量化估算技术[J]. 物探与化探, 2023, 47(3): 747-756 doi:10.11720/wtyht.2023.2579

LIU Jun, HUANG Chao, YANG Lin, ZHANG Yong-Sheng, ZHA Ming. Quantitative prediction technology for the hydrocarbon production capacity of fractures and vugs in fault-controlled carbonate reservoirs in the Shunbei area[J]. Geophysical and Geochemical Exploration, 2023, 47(3): 747-756 doi:10.11720/wtyht.2023.2579

0 引言

走滑断裂带及其控制的缝洞体形成了碳酸盐岩一种特有的断控缝洞体油气圈闭,在塔里木盆地深部碳酸盐岩油气藏的勘探开发中具有重要意义[1-2]。断控缝洞体圈闭成为顺北地区奥陶系新的勘探领域,已取得重大的油气突破。近期,重点探井顺北42X井获高产工业油气流,日产气82.2 万m3、日产油300 t,油气当量近千吨。

顺北油气田为超深油藏,深度在7 300 m以下,目前超过8 000 m的钻井有40口,是亚洲最深商业开发油气田。因此,开展地下8 000 m断裂带立体特征的精细刻画与储集体含油气性预测,强化规模储集体预测与目标评价,深化油藏特征与富集规律认识,优化井轨迹设计,努力实现少井高产,是地震预测的艰巨任务和目标。

当前,对断控缝洞体储层的地震预测技术研究取得不少进展。吴涛等[3]以塔河油田托普台北为例研究了走滑断裂对碳酸盐岩断控缝洞体油藏形成的控制作用;马乃拜等[4]依据塔里木盆地顺北地区断控缝洞体特征,进行了地震反射特征与识别的正演模拟研究;王震等[5]对塔河油田碳酸盐岩断溶体的刻画技术取得了新的认识,用结构张量取代精细相干、蚂蚁体等属性,用于刻画断溶体的轮廓。 李鹏飞等[6]综合利用 GeoEast 解释系统中的构造解释、储层预测和地震属性融合等技术对塔北地区某区块进行断溶体描述与精细刻画研究;徐红霞等[7]利用GeoEast系统在塔北地区开展断溶体多属性预测。

20世纪90年代,亿吨级塔河油田的发现,开启了碳酸盐岩缝洞储集体的地震预测方法技术的发展。唐文榜、姚姚等[8-9]指出充填流体溶洞可以形成强反射,后来被业内学者形象地称之为串珠;刘坤岩等[10]对塔河奥陶系隐蔽溶洞体进行地震精细识别;王晓梅、孙海宁等[11-12]提出用反演低波阻抗和AVO分析同时提取的强纵波和强横波振幅用于预测溶洞充填流体。李宗杰、樊佳方等[13-14]深化了缝洞体及其含油气性的预测,采用频率特征差异提取异常体属性表征洞穴的体量,基于频谱特征差异开展缝洞储集体中的油水识别与预测。与此同时,许多专家采用基于分频相干体的蚂蚁追踪技术、地震数据结构张量相干分析技术等开展裂缝预测、不同尺度裂缝预测断裂系统的精细分析和刻画[15-21]

实现断控缝洞体储层单井产量的量化预测是油田勘探开发的现实需要,前期对断控缝洞体储层的研究和历年来对碳酸盐岩缝洞型储集体的预测研究为断控缝洞体储层单井产量的量化预测奠定了基础。断控缝洞体储层单井产量既与断裂纵向延伸和横向延展有关,更与断控缝洞体规模、断裂裂缝分布密度、缝洞储集体的体量有关。充分利用现有地震断裂解释技术,进行断裂的精细解释,得到表征断裂通源性和汇聚性等成藏断裂特征;利用异常体振幅属性、结构张量属性和最大似然概率属性等3种属性分析技术,用以表征断控缝洞体储层规模和断裂裂隙分布密度,即为缝洞特征。对断裂特征值和缝洞特征进行量化处理并加以融合,得到断控缝洞体特征值,将其与实钻井的年产液量进行统计分析,得到年产液量与断控缝洞体特征值的统计关系式,即可用于断控缝洞体储层单井产量在一定精度范围内的量化预测。

1 断控缝洞体储层基本特征

断控缝洞体是碳酸盐岩特有的一种圈闭类型,为受多期次构造作用,沿深断裂带发育的局部脆性灰岩破碎带缝洞系统,在上覆泥灰岩、泥岩等盖层封堵以及侧向致密灰岩遮挡下,形成一种不规则断控缝洞体圈闭类型[5]。大型断裂带是油气疏导的有利通道,缝洞系统是油气富集成藏的有利空间,不同级次走滑断裂的变形样式和缝洞发育程度,控制了断控缝洞体油藏的空间分布及油藏形态,构成了一种与传统砂岩层状介质油藏完全不同的油藏类型。

钻井、岩心、成像测井等资料综合分析表明,顺北地区奥陶系储层类型为沿深大断裂发育的裂缝—洞穴型储层,储集空间类型可划分为洞穴、孔洞、裂缝。首先,碳酸盐岩断控缝洞体储层,既有大、中、小型断裂裂缝,也有沿大断裂分布的大小洞穴、孔洞。洞穴是储集空间的主体,大型断层和断裂裂缝是将储集空间与源岩沟通的油气运移通道。其次,断控缝洞体储层中断裂、断层的纵向延伸和横向延展直接关系到储集体与油源的沟通和汇聚;断裂类型及沿断裂分布的缝洞体等储集空间的体量决定了富集成藏规模;断裂裂缝分布密度同时决定了缝洞的连通性及油气产出能力的大小。储层的规模由大型洞穴的体量和油藏富集程度两大因素决定,断裂裂缝连通储集体产出能力强。所以,研究断裂分布特征和缝洞的体量和分布特征,将其加以量化,进行单井产量量化预测,得到断控缝洞体储层中的高产区分布特征,指导碳酸盐岩断控缝洞体油气藏规模储层预测,提高钻井成功率。

1.1 表征通源性和汇聚性的断裂特征

通源性和汇聚性是断控缝洞体成藏的断裂特征,在单井产量预测中,起到控制作用。

1)断裂的通源性。顺北地区断控缝洞体储层中的大型断裂按成因可分为挤压、拉分和平移3种类型。图1给出了挤压、拉分和平移3种类型断裂的地震反射特征:①挤压式:正花状断裂、能量强、断裂带宽度大,从浅到深的纵向延伸度大,形成良好的连通性;②拉分式:负花状断裂、能量弱、错断明显,纵向延伸度也很大,断裂带宽度也较大,加上拉张力的作用可形成较大空隙,具有良好的连通性;③平移式:以线性异常出现在平面图上,断裂带宽度小,连通性中等。

图1

图1   断溶体储层中不同类型断裂特征的地震响应特征

a—挤压断裂;b—拉分断裂;c—平移断裂

Fig.1   Seismic response characteristics of different fault types

a—compressional fault;b—pull-apart fault;c—translation fault


根据断裂规模将断裂细化分为五类:主干1级、主干2级、分支、次级(包括褶曲)、隐伏断裂。在地震数据上采用断裂精细解释与最大似然概率属性等多属性分析相结合的方法,对断裂进行精细解释,确定出主干断裂,分支断裂,并进行平面分段的活动期次的厘定。同时加强小断裂的多属性识别与描述,在用相干属性和趋势面识别较大尺度断裂基础上,联合振幅变化率、曲率、AFE、倾角、蚂蚁及照明属性,以及AFE属性与能量类属性联合对低序级断裂进行识别,确定断裂的纵向延伸、横向延展、断裂带的规模。

钻井资料分析表明,挤压和拉分式断裂是组成主干断裂的主要类型,这两类断裂的断裂带宽度大、纵向延伸度大,是良好的通道。顺北地区寒武系之下(即T80反射之下)存在深部油源,若断裂向深部延伸穿至震旦系(即T90反射之下),表示该断裂具有良好通源性。

2)断裂的汇聚性。顺北地区存在大小断裂相交的平面结构特征,使更大范围的油气汇聚在相交部位。图2展示的顺北 1 号断裂为NE向主干断裂,有 7 条NNE向分支断裂和1条NNW向分支断裂与其相交,形成断裂的平面结构特征。断裂彼此相交,表示不仅有深部油源通向于此,还有横向更大范围的油源在此汇聚,从而形成富集区。除了主干断裂和分支断裂相交之外,还可见到主干断裂有与多条平行小断裂相交的情况,如 3 号断裂。

图2

图2   断裂平面相交的结构特征

Fig.2   Structural characteristics plane of intersecting fracture


综上所述,断裂的规模和纵向延伸表示断控缝洞体储层内部具有良好的通源性,断裂的平面交叉结构特征表示断控缝洞体储层具有更大范围的油气汇聚性。通源性和汇聚性是断控缝洞体储层富集的重要条件。

1.2 表征断控缝洞体规模和缝洞体体量的缝洞特征

断控缝洞体规模和缝洞体体量等缝洞特征是成藏的空间条件。选择结构张量、异常体振幅、最大似然概率3种属性来表征断控缝洞体储层规模、缝洞体体量和断裂分布密度。

结构张量属性与钻时变化相结合已经用于确定断控缝洞体边界,圈定断控缝洞体储层有利范围和规模[5]。异常体分解技术是一种从层状反射波场中提取纵、横向尺度有限、非均匀分布的非规则洞穴反射波场的技术,异常体振幅属性可以表征洞穴的体量[14]图3为偏移剖面及相应的异常体振幅属性、结构张量属性和最大似然概率属性剖面,图4是不同层位的异常体、结构梯度张量和最大似然概率3种属性平面。

图3

图3   偏移剖面及其3种属性剖面

a—偏移剖面;b—异常体振幅属性剖面;c—结构梯度张量属性剖面;d—最大似然概率属性剖面

Fig.3   The migration and its three attributes section

a—migration;b—abnormal body amplitude attribute;c—structural gradient tensor attribute;d—maximum likelihood probability attribute


图4

图4   不同层位异常体振幅、结构梯度张量和最大似然概率3种属性平面

a—T76-T78 异常体振幅属性;b—T76-T78 结构梯度张量属性;c—T76-T78最大似然概率属性

Fig.4   Plans of attributes of abnormal body amplitude,structural gradient tensor and maximum likelihood probability at different horizons

a—T76-T78 abnormal body amplitude attribute;b—T76-T78 structural gradient tensor attribute;c—T76-T78 maximum likelihood probability attribute


图3图4可以看到:结构张量属性具有成片分布特征,描述了断控缝洞体有利储层的分布范围和规模;最大似然概率属性以非规则线状特征表示大小断裂的分布,反映了缝洞体之外的断裂裂缝的细节,使偏移剖面上出现的断裂更加清晰,可提高地震断裂解释的精度,也可表示断裂裂缝分布密度。异常体振幅属性则表征单个尺度较大洞穴,这些洞穴随机分散分布或沿大断裂纵向分布,从图3可见异常体表征的洞穴主要分布在奥陶侵蚀面(T74之下)至寒武系中部(T82之上)。3种属性从不同侧面描绘出断控缝洞体储层的范围和内部结构,它们分别反映了断控缝洞体整体分布特征(包括断裂的尺度特征),大小断裂和洞穴的分布和断裂裂缝分布密度的特征;既反映了通源性、汇聚性等特征,也反映了富集成藏的条件。

2 断控缝洞体储层基本特征的量化表征

将通过精细地震断裂解释用以表征断裂通源性和汇聚性(“二性”)的特征称之为断裂特征,将异常体振幅、结构梯度张量和最大似然概率3种地震属性(“三属性”)描述的断控缝洞体储层规模、储集空间体量、断裂裂缝分布密度的特征称之为缝洞特征。断裂特征和缝洞特征可作为断控缝洞体储层油藏富集程度和富集区分布范围的量化分析基础。将断控缝洞体储层中由“二性”表示的断裂特征和由“三属性”表征的缝洞特征量化之后融合为断控缝洞体特征值,将其与单井产量进行统计分析,得到单井产量与断控缝洞体特征值的统计关系式,实现量化预测。

2.1 断裂特征的量化分析

图5图6为典型的高产井的剖面属性和平面属性。图5是过高产井的偏移剖面及其异常体属性、结构梯度张量属性和最大似然属性剖面,偏移剖面和最大似然属性展示出井旁有一条断裂向深部延伸到T90之下,表明此断层沟通了深部油源,成为深部油源运移通道。结构梯度张量属性圈定的断控缝洞体向下延伸到T90之下,表明断控缝洞体储层是下延伸到T90之下的大型断控缝洞体储层。异常体属性表示在T74之下和T81T82之间存在多个大型洞穴,具有良好的储集空间。

图5

图5   过高产井的偏移剖面(a)、异常体振幅属性(b)、结构梯度张量属性(c)、最大似然概率属性(d)剖面

Fig.5   Sections of migration(a),abnormal body amplitude attribute(b),structural gradient tensor attribute(c),maximum likelihood probability attribute(d) across a high productivity well


图6

图6   主干断裂和分支断裂交汇部位的高产井

Fig.6   High-productivity wells at the intersection of main and branch faults


图6是顺北1区T74~T74+500 ms时窗范围内的异常体属性切片。图中有一NE向主干断裂,该主干断裂的东北段存在大小不等的分支断裂。沿主干断裂钻有多口井,其中高产井S3、S7、 S11、S2、S1在分支断裂与主干断裂的交汇部位。说明主干断裂与分支断裂的交汇部位具有良好的油气汇聚性,有利于油气富集成藏。

将断裂的通源性量化为A1,汇聚性量化为A2,将二者之和定义为断裂特征值A:

A=A1+A2

式中通源性数值A1和汇聚性数值A2都可以分为两级。断裂的纵向延伸深度在T90 以下的A1为I级,延伸深度在T90 以上的A1为II级;主干断裂与分支断裂的交汇部位的A2为I级,无分支断裂交汇部位的A2则为II级。将I级赋值为0.5,II级赋值为0,于是断裂特征值取值为0、0.5或1.0。A值大表示有利于形成高产储层。表1给出了通源性值A1和汇聚特征值A2的量化分级数值。

表1   断裂的通源性和汇聚性的量化分级

Table 1  Quantitative classification of nature to communicate oil sources and nature to convergence more oil sources

断裂特征断裂通源性A1断裂汇聚性A2
延伸至
T90以下(I)
延伸在
T90以上(II)
主干+分支
交汇部位(I)
仅有主
干断裂(II)
量化值0.500.50

新窗口打开| 下载CSV


2.2 缝洞体特征的量化

断控缝洞体储层中,断裂裂缝是运移通道,也是储集空间,洞穴则是储集空间的主体,因此洞穴的体量是确定断控缝洞体储集空间体量的关键。

为说明异常体振幅属性值大小用以表征洞穴体量的大小,进行了正演模拟分析(图7)。图7a是宽度为 60 m、高度为0~40 m的充填流体洞穴的正演模拟偏移剖面,洞穴形成了串珠状的短反射。图7c给出了串珠振幅与高度的关系曲线。可以看到:洞穴高度从2 m增高到40 m,串珠的振幅随洞穴高度变化大致呈线性增大。图7b7d分别是高度6 m、宽度10~ 300 m的洞穴正演模拟偏移剖面及串珠的振幅随宽度变化曲线。可以看出串珠振幅总体上随宽度的变化分为两个区间:在宽度小于100 m 范围内,振幅随宽度的加宽大体按线性增加,但在这一宽度范围内,串珠的宽度随洞穴宽度的加大变化很小;当宽度大于100 m,随洞穴宽度的加大,振幅则基本保持不变,但偏移剖面上串珠的宽度随洞穴宽度的加大而变宽。

图7

图7   串珠振幅与高度和宽度相关关系的正演模拟

a—宽度60 m、高度0~40 m洞穴正演模拟偏移剖面;b—高度6 m、宽度 10~300 m 洞穴正演模拟偏移剖面;c—宽度60 m、高度0~40 m洞穴串珠振幅变化曲线;d—高度6 m、宽度 10~300 m 洞穴串珠振幅变化曲线

Fig.7   Forward modeling of the correlation between the amplitude of beads and the height and width of caves

a—the cave forward modeling migration profile with a width of 60 m and a height of 0~40 m;b—the cave forward modeling migration profile with height of 6 m and width of 10~300 m;c—the amplitude variation curve of the beads with a width of 60 m and a height of 0~40 m;d—the amplitude variation curve of the beads with a height of 6 m and a width of 10~300 m


模拟试验表明,洞穴高度增大时,串珠振幅与高度变化大体成正比。洞穴宽度增大时,宽度在100 m范围内串珠振幅与洞穴宽度变化成正比;洞穴宽度大于100 m时,虽然振幅不随洞穴宽度发生变化,但串珠宽度随洞穴宽度的加大而加大,对于一个完整串珠,串珠振幅之和会随洞穴加宽而加大,即串珠振幅之和与洞穴的体量(长、宽、高)的变化呈正相关关系。

由于断控缝洞体储层具有一定的空间体积,因此用一定体积范围内的属性值,即体属性值,可以表示断控缝洞体储层规模。体属性值按下述步骤求取。

1)对地震偏移数据体进行结构梯度张量、异常体分解和最大似然概率分析,得到结构张量、异常体振幅和最大似然概率3个属性数据体;

2)在给定层位的每一个CRP上求取一定时窗的均方根属性值,得到该点一定厚度内的均方根属性值;

3)以每一个CRP为中心的圆周内,求取圆内多个相邻CRP的均方根属性值之和,作为该中心点的属性值,即为该点的体属性值;

4)将异常体振幅的体属性值记为B1,将结构梯度张量属性值记为B2,将最大似然概率属性值记为B3,求3种属性值B1B2B3之和,得到复合属性B,用以表征断控缝洞体储层的规模和断裂密度,称之为缝洞特征值:

B=B1+B2+B3

2.3 断控缝洞体特征值

根据地震断裂精细解释结果得到的通源性和汇聚性,按照表1所给出的特征量化值得到断裂特征值A,根据结构张量、异常体振幅和最大似然概率的体属性值按照式(2)计算得到表征断控缝洞体规模和断裂密度的缝洞特征值B,即可按照式(3)得到断控缝洞体特征值 EFK:

EFK=(1+A)×B

3 断控缝洞体储层单井产量量化预测

选择了单井年产液量作为产量指标与断控缝洞体特征值进行统计分析,得到产量指标与断溶体特征值之间相关关系。

3.1 单井年产液量的量化表达式

将断控缝洞体特征值与单井年累计产液量进行对比分析。表2列出了16口样本井单井年产液量、断裂特征值A、缝洞特征值B,以及由此计算出的断控缝洞体特征值EFK

表2   16口样本井年产液量SWOC与断控缝洞体特征值EFK统计

Table 2  Annual fluid production SWOC and fault-controlled reservoirs features value EFK of 16 sample wells were calculated

井名单井年产液量/ t断裂特征值缝洞特征值断控缝洞体特征值
SWOCA1A2ABEFK=(1+A)×B
S1770000.50.513484.1846968.368
S2540000.50.513546.9037093.806
S34200007605.8667605.866
S4270000.500.54052.3426078.515
S52500001985.6221985.622
S62400002696.0852696.085
S72400001354.2541354.254
S82200002408.2702408.27
S92100002810.1412810.141
S102100002207.3042207.304
S11190000521.6521.6
S121800003565.8603565.860
S131200003936.6283936.628
S14700001609.0801609.080
S15500001503.8011503.801
S16100001406.61406.6

新窗口打开| 下载CSV


图8是单井年产液量SWOC与断控缝洞体特征值EFK的变化趋势曲线,可以看到二者的数值变化有着相似的趋势。因此,将两种数值进行拟合得到图9所示的线性拟合曲线,其拟合关系式为

SWOC=0.8366EFK×895.48

图8

图8   单井年产液量与断控缝洞体特征值变化趋势曲线

Fig.8   Variation trend curve of single well annual fluid yield and fault-controlled reservoirs features value


图9

图9   单井年产液量与断控缝洞体特征值的拟合关系曲线

Fig.9   The fitting relation curve between the annual liquid production of a single well and the fault-controlled reservoirs features value


SWOCEFK的拟合度接近65%,即

R2=0.6448

按照式(4),根据每一个点的断控缝洞体特征值EFK,即可估计出每一个点以年累计产液量表示的单井年产量SWOC,实现断控缝洞体储层单井产量的量化预测。

3.2 断控缝洞体储层单井产量多属性量化预测使用

根据单井年产液量与断控缝洞体特征值的拟合曲线式(4),预测得到的顺北1区的单井年产液量,按年产液量划分出高产、中产和低产区域,高于3 万t为高产区,低于1 万t为低产区,介于二者之间为中产区。作出预测产量平面分布(图10),图中展示出沿断控缝洞体(主要为主干断裂)的产量分布,色标的不同颜色表示不同的产量,红色表示高产区,橙色表示中产区,蓝色表示低产区。在图上同时标注出高产井、中产井和低产井。可以看到在顺北1区,1号断裂带北段和中段产量较南段高。高产井主要集中在主干断裂北段与分支断裂交汇部位。

图10

图10   顺北1区预测年累计产液量平面

Fig.10   Plan of predicted annual liquid production in Shunbei 1 area


为清晰起见,将顺北1区实测与预测单井年累计产液量列于表3中。按照高、中、低三级划分的产量进行对比,20口钻井中预测符合13口,符合率为65%,与式(4)的拟合度64.48% 基本一致。表明以年累计产液量为产量指标的多属性量化预测,其符合率可达到65% 。

表3   顺北1区年产液量预测符合率

Table 3  Coincident rate of predicted annual liquid production in Shunbei 1 area

序号井名年累计产
液量/万t
按产液
量分级
预测产
量分级
预测
符合情况
1S17.7高产井高产
2S25.4高产
3S34.2高产
4S43.6低产×
5S52.7中产井低产×
6S62.5中产
7S72.4高产×
8S82.4中产
9S92.2低产×
10S102.1高产×
11S112.1中产
12S121.9中产
13S131.3中产
14S141.2低产×
15S151.1中产
16S161中产
17S170.5低产低产
18S180.1低产
19S190.1低产
20S200中产×
总计
预测符合率
137
65%35%

新窗口打开| 下载CSV


4 结论与讨论

针对顺北地区奥陶系断控缝洞体储层特征,开展了单井产量的量化预测技术研究,提出了由“二性”“三属性”组成的断控缝洞体油气产能定量化估算技术。该技术的实现由四部分组成:①在走滑断裂带精细解释基础上得到表示断裂纵向延伸的通源性和断裂平面展布结构的汇聚特征,并将两种特征量化为断裂特征值;②以异常体、结构张量、最大似然属性等3种属性融合,得到表征断控缝洞体规模和断裂裂缝分布密度的断控缝洞体征值;③将断裂特征值与缝洞特征值融合为断控缝洞体特征值;④将断控缝洞体特征值与实测单井年累计产液量进行统计分析,得到年累计产液量与断控缝洞体特征值的拟合表达式,实现断控缝洞体储层单井产量的量化预测,预测符合率达到65%。

断控缝洞体储层单井产量的量化预测仅仅走出第一步,其有效性、稳定性还有待于深入研究。用于表征产量指标的选择,使用年累计产液量,还是单位压降产量,或者其他指标,还需进行对比试验。断裂特征的优选和断裂特征值的确定,空间每一个点上异常体、结构梯度张量和最大似然概率等3种体属性值计算的时窗大小、横向范围的选取也还要深入试验。尽管如此,目前研究中发现的断裂特征值与缝洞特征值融合得到的断控缝洞体特征值与年累计产液量之间变化趋势的相似性,为产量量化预测的统计分析工作打下了基础,为技术方法的发展提供了线索。通过进一步深入研究之后得到一种或多种有效的方法来预测断控缝洞体的产量,为勘探开发提供技术支撑。

参考文献

鲁新便, 胡文革, 汪彦, .

塔河地区碳酸盐岩断溶体油藏特征与开发实践

[J]. 石油与天然气地质, 2015, 36(3):347-355.

[本文引用: 1]

Lu X B, Hu W G, Wang Y, et al.

Characteristics and development practice of fault-karst carbonate reservoirs in Tahe area,Tarim Basin

[J]. Oil & Gas Geology, 2015, 36(3):347-355.

[本文引用: 1]

焦方正.

塔里木盆地顺北特深碳酸盐岩断溶体油气藏发现意义与前景

[J]. 石油与天然气地质, 2018, 39(2):208-216.

[本文引用: 1]

Jiao F Z.

Significance and prospect of ultra-deep carbonate fault-karst reservoirs in Shunbei area,Tarim Basin

[J]. Oil & Gas Geology, 2018, 39(2):208-216.

[本文引用: 1]

吴涛, 戴少康, 曹飞, .

走滑断裂系对碳酸盐岩“断溶体”油藏形成的控制作用——以塔河油田托普台北为例

[J]. 地球科学前沿, 2017, 7(5):681-694.

[本文引用: 1]

Wu T, Dai S K, Cao F, et al.

Strike slip fault system and it's controlling on the formation of "reservoir of fault controlling dissoulution":A case study of Toufutai Area,Tahe Oil Fields

[J]. Advances in Geosciences, 2017, 7(5):681-694.

DOI:10.12677/AG.2017.75069      URL     [本文引用: 1]

马乃拜, 金圣林, 杨瑞召, .

塔里木盆地顺北地区断溶体地震反射特征与识别

[J]. 石油地球物理勘探, 2019, 54(2):398-403.

[本文引用: 1]

Ma N B, Jin S L, Yang R Z, et al.

Seismic reservoir characteristics and identification of fault karst reservoir in Shunbei area,Tarim Basin

[J]. Oil Geophysical Prospecting, 2019, 54(2):398-403.

[本文引用: 1]

王震, 文欢, 邓光校, .

塔河油田碳酸盐岩断溶体刻画技术研究与应用

[J]. 石油物探, 2019, 58(1):149-154.

DOI:10.3969/j.issn.1000-1441.2019.01.017      [本文引用: 3]

碳酸盐岩断溶体是近年来发现的一种重要的碳酸盐岩圈闭,常规断裂检测属性如精细相干、蚂蚁体等方法以表征断溶体内部断裂分布特征为主,不能很好地刻画断溶体圈闭的轮廓特征。引入梯度结构张量方法用于刻画断溶体的轮廓,并对其单一特征值及组合特征值进行断溶体轮廓刻画效果的对比分析。结果表明,梯度结构张量第二特征值的高值区域与断溶体的反射异常区域较为一致,其门槛值可利用实钻井的钻时曲线进行标定。在断溶体轮廓刻画的基础上,融合反演孔隙度数据体及蚂蚁数据体属性,表征其内部溶洞及断裂发育特征,由此形成一套刻画断溶体外部轮廓及内部结构的技术。塔河油田实际地震资料应用表明,该技术在刻画断溶体方面具有很好的效果,可直观体现断溶体在空间上发育的差异性,亦能刻画其内部缝洞结构。

Wang Z, Wen H, Deng G X, et al.

Fault-karst characterization technology in the Tahe Oilfield,China

[J]. Geophysical Prospecting for Petroleum, 2019, 58(1):149-154.

DOI:10.3969/j.issn.1000-1441.2019.01.017      [本文引用: 3]

Fault-karst is an important type of trap that was recently discovered.Conventional fracture detection attributes,such as fine coherence and ant body,are mainly used to characterize fracture distribution inside the fault-karst instead of structural characteristics.This study uses gradient structure tensor to characterize the structural characteristics of fault-karst.The application effect of single eigenvalues and combined eigenvalues were analyzed in this paper.The results of contrastive analysis indicated that the high value region of the second eigenvalue is consistent with the reflection anomaly region of fault-karst,and its threshold value can be calibrated by using the drilling time curve.Based on the structural characterization using the proposed technique,porosity inversion and ant attribute can be used to characterize the caves and fractures inside the fault-karst.Therefore a set of techniques for describing the structure and internal structure of fault-karst were developed.Testing on field data showed that the proposed technique could reflect the structural spatial difference and the internal fracture hole.

李鹏飞, 崔德育, 田浩男.

塔里木盆地塔北地区X区块断溶体刻画方法与效果

[J]. 石油地球物理勘探, 2017, 52(s1):189-194.

[本文引用: 1]

Li P F, Cui D Y, Tian H N.

Fault karst carbonate reservoir description with GeoEast interpretation sub-system in the Tabei area,Tarim Basin

[J]. Oil Geophysical Prospecting, 2017, 52(s1):189-194.

[本文引用: 1]

徐红霞, 沈春光, 李斌, .

多属性分析技术在碳酸盐岩断溶体预测中的应用

[J]. 石油地球物理勘探, 2017, 52(s2):158-163.

[本文引用: 1]

Xu H X, Shen C G, Li B, et al.

Fault-karst carbonate reservoir prediction with comprehensive multi-attribute analysis

[J]. Oil Geophysical Prospecting, 2017, 52(s2):158-163.

[本文引用: 1]

唐文榜, 韩革华, 刘来祥, .

溶洞充填物判识的频率差异分析技术

[J]. 石油与天然气地质, 2002, 23(1):41-44.

[本文引用: 1]

Tang W B, Han G H, Liu L X, et al.

Analysis technique of difference for discrimination of cavity fillers

[J]. Oil & Gas Geology, 2002: 23(1):41-44.

[本文引用: 1]

姚姚, 唐文榜.

深层碳酸盐岩岩溶风化壳洞缝型油气藏可检测性的理论研究

[J]. 石油地球物理勘探, 2003, 38(6):623-629.

[本文引用: 1]

Yao Y, Tang W B.

The theoretical researching of detecting for deeply carbonate grotto mantlerrock cavity reservoir

[J]. Oil Geophysical Prospecting, 2003, 38(6):623-629.

[本文引用: 1]

刘坤岩, 许杰.

塔河奥陶系隐蔽溶洞体地震精细识别

[J]. 石油地球物理勘探, 2019, 54(5):1106-1114.

[本文引用: 1]

Liu K Y, Xu J.

Identification of Ordovician subtle cave reservoirs in Tahe on seismic data

[J]. Oil Geophysical Prospecting, 2019, 54(5):1106-1114.

[本文引用: 1]

王晓梅, 唐文榜, .李学著,

碳酸盐岩溶洞型储层反演方法

[J]. 物探与化探, 2005, 29(1):44-46.

[本文引用: 1]

Wang X M, Tang W B, Li X Z, et al.

An inversion techniques for water-eroded case type carbonate reservoirs

[J]. Geophysical and Geochemical Exploration, 2005, 29(1):44-46.

[本文引用: 1]

孙海宁, 王晓梅, 刘来祥.

AVO技术在识别充填流体溶洞中的应用

[J]. 物探与化探, 2008, 32(4):397-399.

[本文引用: 1]

Sun H N, Wang X M, Liu L X.

The application of AVO to the predication of water-eroded caves filled with liquids for carbonate reservoirs

[J]. Geophysical and Geochemical Exploration, 2008, 32(4):397-399.

[本文引用: 1]

李宗杰, 刘群, 李海英, .

塔河油田缝洞储集体油水识别的谐频特征分析技术应用研究

[J]. 石油物探, 2014, 53(4):484-490.

DOI:10.3969/j.issn.1000-1441.2014.04.015      [本文引用: 1]

针对塔河油田下奥陶统碳酸盐岩缝洞型储集体的油水识别难题,开展了以频谱特征差异为基础的瞬时谐频特征(Harmonic Frequence Characters,HFC)分析技术应用研究。充填流体溶洞、坍塌-充填洞穴模型的正演模拟数据和井旁地震道产层段的频谱分析表明,主频段振幅是区分充填流体洞穴与坍塌-充填洞穴的主要依据,30 Hz以上的高频段振幅的差异是区分油和水的主要标志,由此形成缝洞型储集体的油水识别模式。对塔河油田十区和托普台地区各200 km<sup>2</sup>三维地震资料开展HFC油水识别应用研究,将识别结果用43口盲井进行验证,其识别符合率达到81.4%。

Li Z J, Liu Q, Li H Y, et al.

Application of HFC technique for hydrocarbon identification in fracture-cave reservoirs of lower-Ordovician Carbonate in Tahe oilfield

[J]. Geophysical Prospecting for Petroleum, 2014, 53(4):484-490.

[本文引用: 1]

樊佳方, 李宗杰, 韩革华, .

基于属性差异薄储集层(体)流体预测技术

[C]// 北京: SPG/SEG 南京2020国际地球物理会议论文集, 2020:723-726.

[本文引用: 2]

Fan J F, Li Z J, Han G H, et al.

Fluid prediction techniques for thin reservoirs based on seismic attribute differences

[C]// Beijing: SPG/SEG Nanjing 2020 International Geophysical Conference, 2020:723-726.

[本文引用: 2]

马艺璇, 李慧莉, 刘坤岩, .

基于分频相干体的蚂蚁追踪技术在塔河油田断裂刻画中的应用

[J]. 石油物探, 2020, 59(2):258-266.

DOI:10.3969/j.issn.1000-1441.2020.02.012      [本文引用: 1]

利用全频带常规蚂蚁追踪技术刻画断裂易出现横向连续性不强或平面断裂格局不明显的现象,造成断裂刻画与井间动态响应匹配度低。提出了基于匹配追踪频谱分解的分频相干体技术,并改进常规蚂蚁追踪技术流程,以提高断裂刻画精度。首先,以原始全频带地震数据为基础,采用匹配追踪频谱分解技术生成一系列分频体,并利用倾角导向滤波提高分频体的信噪比;然后,采用本征相干计算方法获得不同中心频率的分频相干体,对照不同分频剖面及目的层相干切片优选出某一分频相干体;最后对二次去噪后的分频相干体应用蚂蚁追踪技术并结合井区动静态资料,反复调整追踪参数,获得分频蚂蚁体。在此过程中,合理的频谱分解和去噪处理可提高中小尺度断裂的分辨率,本征相干计算和精调参数的蚂蚁追踪可提高中小尺度断裂的识别效果。将此技术应用于塔河油田T井区碳酸盐岩储层的中小尺度断裂识别,结果表明,相较于常规技术,该技术对断裂的刻画精度高,刻画结果与井间动态连通关系匹配性好,与地质认识吻合度高。

Ma Y X, Li H L, Liu K Y, et al.

Application of an ant-tracking technique based on spectral decomposition to fault characterization

[J]. Geophysical Prospecting for Petroleum, 2020, 59(2):258-266.

DOI:10.3969/j.issn.1000-1441.2020.02.012      [本文引用: 1]

The conventional ant-tracking technique applied to fault characterization often results in weak transverse continuity or unobvious plane fracture patterns.It also leads to a poor matching relationship between the characterized faults and the inter-well dynamic response.To solve these issues,we proposed an improved ant-tracking technique based on matching-tracking spectral decomposition and improved the conventional workflow of ant-tracking.Firstly,the full spectral-band seismic data were decomposed into a series of frequency-divided data,and the dip steering filtering method was used to improve the signal-to-noise ratio of the initial seismic data.Then,the intrinsic coherence of seismic data in different frequency bands were calculated and extracted.We compared the profiles and the coherent slices of the target layer from different single-frequency data,so as to optimize the required coherent data.In addition,the data after further noise reduction were tracked by artificial ants,and the appropriate ant-tracking parameters were repeatedly adjusted based on dynamic or static data.Finally,dividing-frequency ant data were obtained by the aforementioned ant-tracking process.The resolution for medium-and small-scale fracture identification was improved through denoising and reasonable spectral decomposition.Furthermore,the effect of fault identification was enhanced by the ant-tracking parameters determination and intrinsic coherence calculation.This method was applied to the identification of medium-and small-scale faults in carbonate reservoirs in well area T in Tahe Oilfield.The results showed that the so-identified faults were consistent with the inter-well dynamic connectivity and geological interpretation.

陈双全, 季敏.

地震数据结构张量相干计算方法

[J]. 石油物探, 2012, 51(3):234-238.

[本文引用: 1]

Chen S Q, Ji M.

Structure tensor coherence computation method of seismic data

[J]. Geophysical Prospecting for Petroleum, 2012, 51(3):234-238.

[本文引用: 1]

梁志强.

不同尺度裂缝的叠后地震预测技术研究

[J]. 石油物探, 2019, 58(5):766-772.

DOI:10.3969/j.issn.1000-1441.2019.05.016      [本文引用: 1]

分析了大尺度、中尺度和小尺度3种不同尺度断裂裂缝的类型、成因和地质特点,总结并提出了这3种不同尺度裂缝的地球物理特征及其综合预测方法。多尺度裂缝地震预测技术主要利用多尺度相干、多尺度曲率以及Likelihood属性。研究认为:利用断面增强、多尺度相干、结构张量等多种属性方法可实现大尺度断裂裂缝的边界刻画和内幕识别;利用多尺度曲率加蚂蚁体的计算方法可获得不同级别的中尺度裂缝发育体的预测结果;利用叠后Likelihood属性可预测小尺度裂缝空间展布。最终形成了一套针对不同级别裂缝的叠后地震预测技术系列,并展示了应用该技术得到的3种不同尺度裂缝地震预测结果。

Liang Z Q.

Post stack seismic prediction techniques for fractures of different scales

[J]. Geophysical Prospecting for Petroleum, 2019, 58(5):766-772.

[本文引用: 1]

刘振峰, 曲寿利, 孙建国, .

地震裂缝预测技术研究进展

[J]. 石油物探, 2012, 51(2):191-196.

DOI:10.3969/j.issn.1000-1441.2012.02.013      [本文引用: 1]

地下裂缝识别和预测对于油气勘探具有重要的实践价值。相关的理论研究和现场试验已经证明裂缝诱发的各向异性会对岩石的地震波传播 特征产生影响,通过地震资料识别和表征地下裂缝的走向、密度及分布范围是可行的。地震识别裂缝方法具有费用低、覆盖范围广和探测深度 大的优点,主要包括多分量转换波裂缝检测、纵波方位各向异性检测和叠后地震属性分析三大类。从应用情况来看,由于数据基础、数据要求 、技术完备程度3方面的原因,不同方法的应用程度及效果存在差异。随着采集和处理等技术的不断进步,多分量转换波裂缝检测、纵波方位各 向异性检测方法会得到更为广泛的应用;而在传统三维叠后数据基础上发展起来的属性分析技术也会因为数据对象类型的多样化而获得更为广 阔的发展空间。

Liu Z F, Qu S L, Sun J G, et al.

Progress of seismic fracture characterization technology

[J]. Geophysical Prospecting for Petroleum, 2012, 51(2):191-196.

DOI:10.3969/j.issn.1000-1441.2012.02.013      [本文引用: 1]

Fracture identification and prediction is significant in hydrocarbon exploration. Corresponding theory and experiments have <BR>shown that the anisotropy caused by fracture will affect seismic reflection features during wave propagation. For this <BR>reason,it is feasible to characterize fracture trend,intensity and distribution by using seismic data. For fracture <BR>prediction,seismic technology is characterized by low cost,extensive coverage area and large probing depth compared with <BR>other methods,which can be divided into three types,including multi-component converted-wave fracture <BR>characterization,azimuth anisotropy detection of P-wave response and post-stack seismic attribute analysis. Various <BR>techniques have different popularity and efficiency because of the difference from data source,data requirement and technique <BR>maturity. With the development of surveying and processing techniques,multi-component converted-wave fracture <BR>characterization and azimuth anisotropy detection of P-wave response will be more widely used in fracture characterization. <BR>For the same reason,seismic attribute anomaly detection will also be more important in using various dataset to recognize <BR>fractured reservoirs.

朱成宏, 黄国骞, 秦瞳.

断裂系统精细分析技术

[J]. 石油物探, 2002, 41(1):42-48.

[本文引用: 1]

断裂系统精确解释是精确构造和岩性分析的基础。概要介绍了利用地震资料断裂系统解释的多种方法的演变过程和相互联系,分析了它们的优缺点,阐述了三代相干体、倾角、方位角、方差及混沌等增强不连续的属性计算方法,给出了属性计算结果细线化的具体算法与应用实例。

Zhu C H, Huang G Q, Qin T.

Methods for detailed fracture system description

[J]. Geophysical Prospecting for Petroleum, 2002, 41(1):42-48.

[本文引用: 1]

Detailed fracture system analysis lays the foundation for precise structure and stratigraphy description. In this paper several fracture system description methods using seismic data are briefly discussed. Their mutual relation, evolution, advantages and weaknesses are also addressed. The approaches to computing three-generation coherence cubes, to extracting enhancing attributes of discontinuity including dip, azimuth, variance, chaos are summarized. A new graph thinning method adapted from computer vision processing is introduced and illustrated with a real data example. Finally, it is pointed out that the quality of input data is vital to automatic seismic data interpretation comparing to various fracture delineation methods.

黄苇, 周捷, 高利君, .

基于同步挤压改进短时傅里叶变换的分频蚂蚁追踪在断裂识别中的应用

[J]. 物探与化探, 2021, 45(2):432-439.

[本文引用: 1]

Huang W, Zhou J, Gao L J, et al.

The application of frequency division ant tracking based on synchronous extrusion improvement of short time Fourier transform in crack detection

[J]. Geophysical and Geochemical Exploration, 2021, 45(2):432-439.

[本文引用: 1]

马玉歌, 苏朝光, 张健, .

低序级断层结构导向坎尼属性边缘检测识别方法

[J]. 物探与化探, 2020, 44(3):698-703.

[本文引用: 1]

Ma Y G, Su C G, Zhang J, et al.

Low-order fault structure-oriented Canny property edge detection and recognition method

[J]. Geophysical and Geochemical Exploration, 2020, 44(3):698-703.

[本文引用: 1]

/

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com , whtbjb@163.com