CN118131330A - 深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备 - Google Patents
深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备 Download PDFInfo
- Publication number
- CN118131330A CN118131330A CN202410476754.8A CN202410476754A CN118131330A CN 118131330 A CN118131330 A CN 118131330A CN 202410476754 A CN202410476754 A CN 202410476754A CN 118131330 A CN118131330 A CN 118131330A
- Authority
- CN
- China
- Prior art keywords
- attribute
- crack
- seismic data
- fracture
- sensitive
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000002708 enhancing effect Effects 0.000 claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims description 40
- 238000001914 filtration Methods 0.000 claims description 29
- 238000011161 development Methods 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 17
- 238000007670 refining Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005755 formation reaction Methods 0.000 description 30
- 238000004364 calculation method Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 9
- 239000011435 rock Substances 0.000 description 6
- 238000007476 Maximum Likelihood Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- -1 compound hydrocarbon Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本申请涉及一种深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备,涉及裂缝预测技术领域,其中方法包括:对初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,第一敏感属性体为对裂缝最为敏感的裂缝属性体,裂缝属性体为反映断裂及裂缝的属性体;基于第一敏感属性体,对初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据;从印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从第二敏感属性体中提取深层碳酸盐岩地层中目的层段的裂缝平面属性,得到目的层段的裂缝预测结果。本申请具有提高深层碳酸盐岩地层中裂缝预测的精度的效果。
Description
技术领域
本申请涉及裂缝预测技术领域,具体涉及一种深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备。
背景技术
深层碳酸盐岩地层通常是指埋藏深度在4500米到6000米之间的碳酸盐岩储层,并且深层碳酸盐岩地层中可能存在多种类型的油气藏,包括高孔高渗常规油气藏、缝洞复合型油气藏等,故而在实际生产生活中,深层碳酸盐岩地层的勘探和开发对于提高能源供应的稳定性具有重要意义。在对深层碳酸盐岩地层的勘探和开发过程中,由于裂缝是碳酸盐岩储层中油气储存和运移的重要通道,因此对裂缝预测有助于确定最佳的钻井位置,有助于提高深层碳酸盐岩地层的勘探成功率。
目前,常用的裂缝预测方法主要依赖叠后地震数据同相轴的波形相似性变化进行裂缝的预测,裂缝预测的精度与波形连续性密切相关,如果叠后地震数据的品质较低,波形杂乱,同相轴连续性差时,预测裂缝往往受到噪音干扰,尤其是在深层碳酸盐岩地层中,其对应的叠后地震数据品质较差时,直接基于常用的裂缝预测方法,以叠后地震数据同相轴错断或者弯曲程度来预测裂缝存在较多假象,导致对深层碳酸盐岩地层中裂缝预测的精度较低。
发明内容
为了提高深层碳酸盐岩地层中裂缝预测的精度,本申请提供一种深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备。
在本申请的第一方面提供了一种深层碳酸盐岩地层的裂缝预测方法,具体包括:
获取深层碳酸盐岩地层的初始叠后地震数据;
对所述初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从所述增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,所述第一敏感属性体为对裂缝最为敏感的裂缝属性体,所述裂缝属性体为反映断裂及裂缝的属性体;
基于所述第一敏感属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据;
从所述印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从所述第二敏感属性体中提取所述深层碳酸盐岩地层中目的层段的裂缝平面属性,得到所述目的层段的裂缝预测结果。
通过采用上述技术方案,对初始叠后地震数据进行断裂增强处理,从而提高识别小断裂和裂缝带的能力,接着从增强后叠后地震数据对应的裂缝属性体中确定对裂缝最为敏感的第一裂缝属性体,从而有效表征和识别裂缝。进一步地,对初始叠后地震数据中第一敏感属性体对应的裂缝印迹进行增强处理,得到印迹增强地震数据,从而实现对初始叠后地震数据中裂缝发育处的同相轴进行打断,增强裂缝发育处同相轴的不连续性,从而增强深层碳酸盐岩地层中低品质的层段的叠后地震数据的裂缝空间上的印迹,降低裂缝预测时低品质的叠后地震数据对预测精度的影响。最后从第二敏感属性体中提取目的层段的裂缝平面属性,得到裂缝预测结果,进而实现对深层碳酸盐岩地层中裂缝预测的精度的提高。
可选的,所述从所述增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,具体包括:
基于所述增强后叠后地震数据,计算第一预设个数不同的第一裂缝属性体;
分别提取各所述第一裂缝属性体的第一平面属性,并基于所述第一平面属性,从各所述第一裂缝属性体中确定第一敏感属性体。
通过采用上述技术方案,增强后叠后地震数据确定后,基于增强后叠后地震数据计算能够表征和识别深层碳酸盐岩地层中裂缝的第一裂缝属性体,接着提取每个第一裂缝属性体的第一平面属性,从而确定每个第一裂缝属性体在二维平面上裂缝分布情况,最后基于各个第一平面属性,将裂缝分布情况较好的第一面属性对应的第二裂缝属性体确定为第一敏感属性体,从而对此第一敏感属性体对应的裂缝印迹进行增强,提高裂缝预测精度。
可选的,所述基于所述第一敏感属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据,具体包括:
对所述第一敏感属性体的地震属性的目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体;
基于所述合并过滤后属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据。
通过采用上述技术方案,对目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体,从而得到噪音较低、可信度以及连续性较高的合并过滤后属性体。进一步地,对初始叠后地震数据中此合并过滤后属性体对应的裂缝印迹进行增强,从而实现对初始叠后地震数据中裂缝发育处的同相轴进行打断,增强裂缝发育处同相轴的不连续性,进而提高基于初始叠后地震数据进行裂缝预测的精度。
可选的,所述对所述第一敏感属性体的地震属性的目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体,具体包括:
在所述第一敏感属性体的地震属性的目标剖面中,沿代表裂缝发育的属性高值对裂缝进行细化,得到细化后属性体,所述细化后属性体表征与预设的地质认识相一致的裂缝;
将所述细化后属性体的水平切片中属性值高于预设门槛值的裂缝确定为目标裂缝,并将方向一致且裂缝间距小于间距阈值的各所述目标裂缝进行合并,得到合并后裂缝;
对裂缝长度小于长度阈值的合并后裂缝进行过滤处理,得到合并过滤后属性体。
通过采用上述技术方案,沿代表裂缝发育的属性高值对目标剖面中的裂缝进行细化处理,得到细化后属性体,从而较好地降低第一敏感属性体中的噪音,同时为后续三维空间裂缝追踪和合并提供数据基础。进一步地,将水平切片中属性值高于预设门槛值的裂缝确定为目标裂缝,并将所有目标裂缝中方向一致并且裂缝间距小于间距阈值的目标裂缝进行合并,得到合并后裂缝。接着对裂缝长度小于长度阈值的合并后裂缝进行过滤处理或删除处理,最终得到合并过滤后属性体,从而有助于提高此第一敏感属性体的可信度以及三维空间的连续性。
可选的,所述基于所述合并过滤后属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据,具体包括:
确定所述合并过滤后属性体对应的三维地质空间位置;
从所述初始叠后地震数据中选取所述三维地质空间位置对应的裂缝印迹,并按照预设的比例减小公式,对所述裂缝印记进行减小,得到印迹增强地震数据。
通过采用上述技术方案,对初始叠后地震数据中此三维地质空间位置对应的叠后地震数据,即,裂缝印记进行缩小处理,得到印迹增强地震数据,从而实现对初始叠后地震数据中的同相轴进行打断,进而提高基于初始叠后地震数据进行裂缝预测的精度。
可选的,所述比例减小公式为:
;
式中,B表示印迹增强地震数据,A表示初始叠后地震数据,a表示合并过滤后属性体,(1,0)表示合并过滤后属性体正常的属性值范围。
通过采用上述技术方案,当表示合并过滤后属性体中的属性值不处于正常的属性值范围内,说明合并过滤后属性体无效,不能被直接用来识别和预测裂缝,那么将初始叠后地震数据确定为印迹增强地震数据。当/>表示合并过滤后属性体中存在处于正常的属性值范围内的有效值,说明合并过滤后属性体有效,那么通过B=(1-a)*A的方式,对裂缝印记进行减小,得到印迹增强地震数据,从而增强深层碳酸盐岩地层中低品质的层段的叠后地震数据的裂缝空间上的印迹,降低裂缝预测时低品质的叠后地震数据对预测精度的影响。
可选的,所述从所述印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,具体包括:
基于所述印迹增强地震数据,计算第二预设个数不同的第二裂缝属性体;
分别提取各所述第二裂缝属性体的第二平面属性,并基于所述第二平面属性,从各所述第二裂缝属性体中确定第二敏感属性体。
通过采用上述技术方案,印迹增强地震数据确定后,基于印迹增强地震数据计算能够表征和识别深层碳酸盐岩地层中裂缝的第二裂缝属性体,接着提取每个第二裂缝属性体的第二平面属性,从而确定每个第二裂缝属性体在二维平面上裂缝分布情况,最后基于各个第二平面属性,将裂缝分布情况较好的第二平面属性对应的第二裂缝属性体确定为第二敏感属性体,从而方便后续提高裂缝预测的精度。
在本申请的第二方面提供了一种深层碳酸盐岩地层的裂缝预测装置,具体包括:
数据获取模块,用于获取深层碳酸盐岩地层的初始叠后地震数据;
属性优选模块,用于对所述初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从所述增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,所述第一敏感属性体为对裂缝最为敏感的裂缝属性体,所述裂缝属性体为反映断裂及裂缝的属性体;
印迹增强模块,用于基于所述第一敏感属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据;
裂缝预测模块,用于从所述印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从所述第二敏感属性体中提取所述深层碳酸盐岩地层中目的层段的裂缝平面属性,得到所述目的层段的裂缝预测结果。
通过采用上述技术方案,数据获取模块获取到初始叠后地震数据后,属性优选模块从增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,接着由印迹增强模块对裂缝印迹进行增强处理,得到印迹增强地震数据。最后裂缝预测模块确定第二敏感属性体,并基于第二敏感属性体得到目的层段的裂缝预测结果。
在本申请的第三方面提供了一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机程序,当所述计算机程序被处理器加载并执行时,执行如第一方面中任意一项所述的方法步骤。
在本申请的第四方面提供了一种电子设备,具体包括:
处理器、存储器及存储在存储器中并能够在处理器上运行的计算机程序,所述处理器用于加载并执行所述存储器中存储的计算机程序,以使所述电子设备执行如第一方面中任意一项所述的方法。
综上所述,本申请包括以下至少一种有益技术效果:对初始叠后地震数据进行断裂增强处理,从而提高识别小断裂和裂缝带的能力,接着从增强后叠后地震数据对应的裂缝属性体中确定对裂缝最为敏感的第一裂缝属性体,从而有效表征和识别裂缝。进一步地,对初始叠后地震数据中第一敏感属性体对应的裂缝印迹进行增强处理,得到印迹增强地震数据,从而实现对初始叠后地震数据中裂缝发育处的同相轴进行打断,增强裂缝发育处同相轴的不连续性,从而增强深层碳酸盐岩地层中低品质的层段的叠后地震数据的裂缝空间上的印迹,降低裂缝预测时低品质的叠后地震数据对预测精度的影响。最后从第二敏感属性体中提取目的层段的裂缝平面属性,得到裂缝预测结果,进而实现对深层碳酸盐岩地层中裂缝预测的精度的提高。
附图说明
图1是本申请实施例提供的一种深层碳酸盐岩地层的裂缝预测方法的流程示意图;
图2是本申请实施例提供的另一种深层碳酸盐岩地层的裂缝预测方法的流程示意图;
图3是本申请实施例提供的一种深层碳酸盐岩地层的裂缝预测装置的结构示意图。
附图标记说明:11、数据获取模块;12、属性优选模块;13、印迹增强模块;14、裂缝预测模块。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。
在本申请实施例的描述中,“示性的”、“例如”或者“举例来说”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示性的”、“例如”或者“举例来说”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言,使用“示性的”、“例如”或者“举例来说”等词旨在以具体方式呈现相关概念。
在本申请实施例的描述中,术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,单独存在B,同时存在A和B这三种情况。另外,除非另有说明,术语“多个”的含义是指两个或两个以上。例如,多个系统是指两个或两个以上的系统,多个屏幕终端是指两个或两个以上的屏幕终端。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”,除非是以其他方式另外特别强调。
参见图1,本申请实施例公开了一种深层碳酸盐岩地层的裂缝预测方法的流程示意图,可依赖于计算机程序实现,也可运行于基于冯诺依曼体系的深层碳酸盐岩地层的裂缝预测装置上。该计算机程序可集成在应用中,也可作为独立的工具类应用运行,具体包括:
S101:获取深层碳酸盐岩地层的初始叠后地震数据。
具体的,叠后地震数据是指地震数据经过动校正和静校正叠加偏移处理后形成的地震数据体,其可以为三维或者二维,在地震勘探中,叠后偏移是一个重要的处理手段,它通过将地震记录进行时间和空间的叠加,来提高地震资料的分辨率和解释质量。叠后地震数据有助于勘探人员更清晰地观察地下构造,为油气勘探和地质灾害预测提供有力支持。在本申请实施例中,初始叠后地震数据为还未经过任何处理的叠后地震数据。一种可行的获取初始叠后地震数据的方式为:通过人工地震的方式产生井位部署备选区域对应的地震数据,最后通过预设的地震仪设备对产生的地震数据进行收集,其中,井为部署备选区域为需要判定是否适合设置井位的区域,通过人工地震的方式产生井位部署备选区域对应的地震数据的具体过程为:通过机械撞击或气爆震源对井位部署备选区域制造振动,进而产生地震数据。进一步地,通过预设的地震数据处理软件GeoSoong软件对地震数据进行叠后偏移处理,最终获取到此井位部署备选区域的深层碳酸盐岩地层的初始叠后地震数据。
在其他实施例中,也可以通过终端将提前准备好的初始叠后地震数据记载到裂缝预测客户端中,最终服务器获取到层碳酸盐岩地层的初始叠后地震数据。其中,终端可以为平板电脑或者个人计算机,终端中安装有裂缝预测客户端,服务器可以为此裂缝预测客户端的后台服务器,具体为独立的物理服务器,也可以是多个物理服务器构成的服务器集群或者分布式系统。
S102:对初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体。
具体的,断裂增强处理是一种旨在提高叠后地震数据中断裂构造的可见性和清晰性。断裂构造是地下岩层中重要的地质特征,对于油气勘探、地质灾害预测等方面具有重要意义。因此,通过断裂增强处理,可以更有效地识别和分析这些断裂构造,为地质解释和决策提供更为准确的信息。裂缝属性体为反映断裂及裂缝的属性体,包含表征地下裂缝特征的属性信息集合,能描述裂缝的方位、开度、密度等属性,并在空间中展现裂缝的分布和走向。换言之,其为在三维空间中对裂缝特征进行表征的三维空间裂缝体。第一敏感属性体为对裂缝最为敏感的裂缝属性体。
在本申请实施例中,初始叠后地震数据确定后,通过构造导向滤波以及分频聚焦等处理,对初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,从而增强此深层碳酸盐岩地层同相轴连续性,并保留大尺度断裂处的地震特征。进一步地,从增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,一种可实现的实施方式中:
基于增强后叠后地震数据,进行断裂属性的计算,计算第一预设个数不同的第一裂缝属性体,在本申请实施例中,各个第一裂缝属性体可以为相干体、曲率体、方差体以及相似性体。在其他实施例中,计算的第一裂缝属性体还可以包括倾角方位角体和最大似然体。其中,相干体是通过对地震道之间相似性的数学描述来揭示深层碳酸盐岩地层的不连续性,具体通过第二代Semblance相干计算公式Semblance=来计算;t0表示样点时间;T表示时窗长度;Xj表示每个样点相似系数。相似性体用于反映增强后叠后地震数据之间的不连续性,具体通过公式Similirity=/>来计算;Xi表示X方向欧式距离,Yi表示Y方向欧式距离。方差体用于反映增强后叠后地震数据的分散或波动的程度,具体通过σ=/> 来计算;j表示计算时窗,i表示地震例数,ω表示方差体的三角加权函数,Xij表示采样点振幅值,/>表示一组采样点的振幅均值。曲率体用于评价裂缝的发育程度,具体通过2D曲率计算公式为:k2d=/>来计算;z表示切线转角,x表示两点之间弧长。倾角方位角体用于确定断层或地层的走向和倾向,有助于揭示地下结构的特征。最大似然体为一种反映断层发育位置的属性体,具体通过最大似然估计(Maximum Likelihood Estimation, MLE)来确定,其是一种在统计学中常用的参数估计方法。需要说明的是,各个第一裂缝属性体的计算过程均可通过预设的Petrel软件确定,此为现有技术,在此不再赘述。
进一步地,通过预设的Petrel软件及其内置的蚂蚁体技术和地震属性分析工具,从各个第一裂缝属性体中提取对应的第一平面属性,其中,第一平面属性为对应的第一裂缝属性体中能裂缝在平面上分布特征的属性,用于识别和评估裂缝发育区域。进一步地,基于提取出的各个第一平面属性,从各个第一裂缝属性体中确定第一敏感属性体的可行方式为:基于人员经验判定最能准确清晰地反映裂缝或断裂的第一平面属性,并通过终端发送选择指令,进而获取终端发送的第一平面属性的选择指令,将选择指令对应的第一裂缝属性体确定为第一敏感属性体。在本申请实施例中,第一敏感属性体可以为相似性体。
S103:基于第一敏感属性体,对初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据。
在一个可实现的实施方式中,对第一敏感属性体的地震属性的目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体;
基于合并过滤后属性体,对初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据。
具体的,地震属性指的是有关地震波的几何学、运动学、动力学或统计学特征的子集,用于地震相分析等。地震属性对应的目标剖面为展示有地下介质的不同物理性质或地质特征的剖面,能够显示地下岩石的物理和几何特性。第一敏感属性体确定后,从第一敏感属性体中提取地震属性对应的目标剖面,然后对目标剖面中的裂缝依次进行细化处理、合并处理和滤波处理,一种可行的实现方式为:
在第一敏感属性体的地震属性的目标剖面中,沿代表裂缝发育的属性高值对裂缝进行细化,得到细化后属性体,细化后属性体表征与预设的地质认识相一致的裂缝;
将细化后属性体的水平切片中属性值高于预设门槛值的裂缝确定为目标裂缝,并将方向一致且裂缝间距小于间距阈值的各目标裂缝进行合并,得到合并后裂缝;
对裂缝长度小于长度阈值的合并后裂缝进行过滤处理,得到合并过滤后属性体。
具体的,属性高值指的是目标剖面中高于属性值阈值的地震属性的属性值。在本申请实施例中,属性高值的地震属性一定程度上能够代表裂缝发育,即,岩石或地质体中裂缝或裂隙系统的发育程度较高。进一步地,通过预设的Petrel软件,沿代表裂缝发育的属性高值对目标剖面中的裂缝进行细化处理,得到细化后属性体,从而较好地降低第一敏感属性体中的噪音,细化处理的具体过程为:去除地震属性中的水平展布属性高值,即消除掉异常的属性高值,异常的属性高值可能由于水平展布效应或其它非地质因素导致,同时通过分析水平应力的方法,确定目标剖面中裂缝的走向方位角度,基于走向方位角度,保留与预设的地质认识相一致的裂缝,从而为后续三维空间裂缝追踪和合并提供数据基础,最终得到细化后属性体。其中,地质认识相为特定环境和条件下形成的沉积岩(物)特征的综合认识和理解,不仅反映了沉积岩的特征,还揭示了沉积环境的信息。
进一步地,通过预设的Petrel软件提取细化后属性体的水平切片,其中,水平切片指的是沿水平方向对裂缝属性体进行的切片操作,即垂直于地面的切片,水平切片能够展示裂缝在某一深度或水平层面上的分布和形态。进一步地,将水平切片中属性值高于预设门槛值的裂缝确定为目标裂缝,并将所有目标裂缝中方向一致并且裂缝间距小于间距阈值的目标裂缝进行合并,得到合并后裂缝。接着将各个合并后裂缝的裂缝长度与长度阈值进行对比,对裂缝长度小于长度阈值的合并后裂缝进行过滤处理或删除处理,最终得到合并过滤后属性体,从而有助于提高此第一敏感属性体的可信度以及三维空间的连续性。
进一步地,基于此合并过滤后属性体,对深层碳酸盐岩地层的初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据,其中,裂缝印迹指的是岩石或其他材料表面由于受力而形成的裂缝所留下的痕迹或特征,在本申请实施例中,裂缝印迹为初始叠后地震数据中合并过滤后属性体对应的叠后地震数据。另外,裂缝印迹增强处理指的是利用特定的技术手段提高初始叠后地震数据中裂缝特征的识别能力和预测精度。
一种裂缝印迹进行增强处理的可行方式为:确定合并过滤后属性体对应的三维地质空间位置;
从初始叠后地震数据中选取三维地质空间位置对应的裂缝印迹,并按照预设的比例减小公式,对裂缝印记进行减小,得到印迹增强地震数据。
具体的,通过预设的Petrel软件确定合并过滤后属性体在深层碳酸盐岩地层中的三维地质空间位置,在其他实施例中,也可以通过预设的Fraca软件确定合并过滤后属性体在深层碳酸盐岩地层中的三维地质空间位置。进一步地,通过预设的比例缩小公式,对初始叠后地震数据中此三维地质空间位置对应的叠后地震数据,即,裂缝印记进行缩小处理,得到印迹增强地震数据,从而实现对初始叠后地震数据中的同相轴进行打断,进而提高基于初始叠后地震数据进行裂缝预测的精度。其中,同相轴为地震记录上各道振动相位相同的极值的连续,其反映地层连续性的地震反射事件,而在裂缝发育区域,这种连续性往往会受到破坏。通过对同相轴的分析,可以识别出这些不连续的特征,从而帮助预测裂缝的存在。在其他实施例中,也可以通过OpendTect解释软件平台得到印迹增强地震数据。
其中,比例缩小公式为:;
式中,B表示印迹增强地震数据,A表示初始叠后地震数据,a表示合并过滤后属性体,(1,0)表示合并过滤后属性体正常的属性值范围。需要说明的是,表示合并过滤后属性体中的属性值不处于正常的属性值范围内,说明合并过滤后属性体无效,不能被直接用来识别和预测裂缝,那么将初始叠后地震数据确定为印迹增强地震数据。/>表示合并过滤后属性体中存在处于正常的属性值范围内的有效值,说明合并过滤后属性体有效,那么通过B=(1-a)*A的方式,对裂缝印记进行减小,得到印迹增强地震数据,从而增强深层碳酸盐岩地层中低品质的目的层段的叠后地震数据的裂缝空间上的印迹,降低裂缝预测时低品质的叠后地震数据对预测精度的影响。其中,目的层段指的是地质勘探和油气开发工作中,深层碳酸盐岩地层中需要重点研究的层段,即,裂缝发育程度较高、具有特殊的油气储集和运移特性。
S104:从印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从第二敏感属性体中提取深层碳酸盐岩地层中目标层段的裂缝平面属性,得到目的层段的裂缝预测结果。
具体的,印迹增强地震数据确定后,基于印迹增强地震数据,重新进行裂缝属性体的计算,即,进行裂缝属性的计算,主要涉及相干体、相似性体、方差体和曲率体等裂缝属性体的计算,具体可参见步骤S102,最终从印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,即,印迹增强地震数据中对裂缝最为敏感的裂缝属性体。进一步地,通过预设的Petrel软件从第二敏感属性体中提取深层碳酸盐岩地层的目的层段的裂缝平面属性,最终得到表征裂缝分布的平面图,即,目的层段的裂缝预测结果。其中,裂缝平面属性指的是目的层段的二维平面上裂缝的几何形态、分布特征以及物理性质等信息。
参见图2,本申请实施例公开了另一种深层碳酸盐岩地层的裂缝预测方法的流程示意图,可依赖于计算机程序实现,也可运行于基于冯诺依曼体系的深层碳酸盐岩地层的裂缝预测装置上。该计算机程序可集成在应用中,也可作为独立的工具类应用运行,具体包括:
S201:获取深层碳酸盐岩地层的初始叠后地震数据。
S202:对初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体。
S203:基于第一敏感属性体,对初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据。
具体的,可参考步骤S101-S103,在此不再赘述。
S204:基于印迹增强地震数据,计算第二预设个数不同的第二裂缝属性体。
S205:分别提取各第二裂缝属性体的第二平面属性,并基于第二平面属性,从各第二裂缝属性体中确定第二敏感属性体。
具体的,印迹增强地震数据确定后,计算第二预设个数不同的第二裂缝属性体,在本申请实施例中,第二预设个数可以为6个,各个第二裂缝属性体可以为相干体、曲率体、方差体、倾角方位角体、最大似然体和相似性体。其中,第二裂缝属性体的计算方式可以参见步骤S102。进一步地,通过预设的Petrel软件提取各个第二裂缝属性体的第二平面属性,基于第二平面属性,从各个第二裂缝属性体中确定第二敏感属性体。在其他实施例中,在第二敏感属性体确定后,可以对第二敏感属性体进行常规的增强或过滤处理,去除其中的噪音,从而有助于凸显裂缝在平面上的规律性。
S206:从第二敏感属性体中提取深层碳酸盐岩地层中目标层段的裂缝平面属性,得到目的层段的裂缝预测结果。
具体的,可参考步骤S104,在此不再赘述。在其他实施例中,从增强或过滤处理后的第二敏感属性体中提取裂缝平面属性,最终得到目的层段的裂缝预测结果。
本申请实施例深层碳酸盐岩地层的裂缝预测方法的实施原理为:对初始叠后地震数据进行断裂增强处理,从而提高识别小断裂和裂缝带的能力,接着从增强后叠后地震数据对应的裂缝属性体中确定对裂缝最为敏感的第一裂缝属性体,从而有效表征和识别裂缝。进一步地,对初始叠后地震数据中第一敏感属性体对应的裂缝印迹进行增强处理,得到印迹增强地震数据,从而实现对初始叠后地震数据中裂缝发育处的同相轴进行打断,增强裂缝发育处同相轴的不连续性,从而增强深层碳酸盐岩地层中低品质的层段的叠后地震数据的裂缝空间上的印迹,降低裂缝预测时低品质的叠后地震数据对预测精度的影响。最后从第二敏感属性体中提取目的层段的裂缝平面属性,得到裂缝预测结果,进而实现对深层碳酸盐岩地层中裂缝预测的精度的提高。
下述为本申请装置实施例,可以用于执行本申请方法实施例。对于本申请装置实施例中未披露的细节,请参照本申请方法实施例。
请参见图3,为本申请实施例提供的深层碳酸盐岩地层的裂缝预测装置的结构示意图。该应用于深层碳酸盐岩地层的裂缝预测装置可以通过软件、硬件或者两者的结合实现成为装置的全部或一部分。该装置包括数据获取模块11、属性优选模块12、印迹增强模块13和裂缝预测模块14。
数据获取模块11,用于获取深层碳酸盐岩地层的初始叠后地震数据;
属性优选模块12,用于对初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,第一敏感属性体为对裂缝最为敏感的裂缝属性体,裂缝属性体为反映断裂及裂缝的属性体;
印迹增强模块13,用于基于第一敏感属性体,对初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据;
裂缝预测模块14,用于从印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从第二敏感属性体中提取深层碳酸盐岩地层中目的层段的裂缝平面属性,得到目的层段的裂缝预测结果。
可选的,属性优选模块12,具体用于:
基于增强后叠后地震数据,计算第一预设个数不同的第一裂缝属性体;
分别提取各第一裂缝属性体的第一平面属性,并基于第一平面属性,从各第一裂缝属性体中确定第一敏感属性体。
可选的,印迹增强模块13,具体用于:
对第一敏感属性体的地震属性的目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体;
基于合并过滤后属性体,对初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据。
可选的,印迹增强模块13,具体还用于:
在第一敏感属性体的地震属性的目标剖面中,沿代表裂缝发育的属性高值对裂缝进行细化,得到细化后属性体,细化后属性体表征与预设的地质认识相一致的裂缝;
将细化后属性体的水平切片中属性值高于预设门槛值的裂缝确定为目标裂缝,并将方向一致且裂缝间距小于间距阈值的各目标裂缝进行合并,得到合并后裂缝;
对裂缝长度小于长度阈值的合并后裂缝进行过滤处理,得到合并过滤后属性体。
可选的,印迹增强模块13,具体还用于:
确定合并过滤后属性体对应的三维地质空间位置;
从初始叠后地震数据中选取三维地质空间位置对应的裂缝印迹,并按照预设的比例减小公式,对裂缝印记进行减小,得到印迹增强地震数据。
可选的,裂缝预测模块14,具体用于:
基于印迹增强地震数据,计算第二预设个数不同的第二裂缝属性体;
分别提取各第二裂缝属性体的第二平面属性,并基于第二平面属性,从各第二裂缝属性体中确定第二敏感属性体。
需要说明的是,上述实施例提供的一种深层碳酸盐岩地层的裂缝预测装置在执行深层碳酸盐岩地层的裂缝预测方法时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的一种深层碳酸盐岩地层的裂缝预测装置与一种深层碳酸盐岩地层的裂缝预测方法实施例属于同一构思,其体现实现过程详见方法实施例,这里不再赘述。
本申请实施例还公开一种计算机可读存储介质,并且,计算机可读存储介质存储有计算机程序,其中,计算机程序被处理器执行时,采用了上述实施例的一种深层碳酸盐岩地层的裂缝预测方法。
其中,计算机程序可以存储于计算机可读介质中,计算机程序包括计算机程序代码,计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间件形式等,计算机可读介质包括能够携带计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM)、随机存取存储器(RAM)、电载波信号、电信信号以及软件分发介质等,需要说明的是,计算机可读介质包括但不限于上述元器件。
其中,通过本计算机可读存储介质,将上述实施例的一种深层碳酸盐岩地层的裂缝预测方法存储于计算机可读存储介质中,并且,被加载并执行于处理器上,以方便上述方法的存储及应用。
本申请实施例还公开一种电子设备,计算机可读存储介质中存储有计算机程序,计算机程序被处理器加载并执行时,采用了上述一种深层碳酸盐岩地层的裂缝预测方法。
其中,电子设备可以采用台式电脑、笔记本电脑或者云端服务器等电子设备,并且,电子设备包括但不限于处理器以及存储器,例如,电子设备还可以包括输入输出设备、网络接入设备以及总线等。
其中,处理器可以采用中央处理单元(CPU),当然,根据实际的使用情况,也可以采用其他通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等,通用处理器可以采用微处理器或者任何常规的处理器等,本申请对此不做限制。
其中,存储器可以为电子设备的内部存储单元,例如,电子设备的硬盘或者内存,也可以为电子设备的外部存储设备,例如,电子设备上配备的插接式硬盘、智能存储卡(SMC)、安全数字卡(SD)或者闪存卡(FC)等,并且,存储器还可以为电子设备的内部存储单元与外部存储设备的组合,存储器用于存储计算机程序以及电子设备所需的其他程序和数据,存储器还可以用于暂时地存储已经输出或者将要输出的数据,本申请对此不做限制。
其中,通过本电子设备,将上述实施例的一种深层碳酸盐岩地层的裂缝预测方法存储于电子设备的存储器中,并且,被加载并执行于电子设备的处理器上,方便使用。
以上所述者,仅为本公开的示例性实施例,不能以此限定本公开的范围。即但凡依本公开教导所作的等效变化与修饰,皆仍属本公开涵盖的范围内。本领域技术人员在考虑说明书及实践这里的公开后,将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未记载的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的范围和精神由权利要求限定。
Claims (10)
1.一种深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述方法包括:
获取深层碳酸盐岩地层的初始叠后地震数据;
对所述初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从所述增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,所述第一敏感属性体为对裂缝最为敏感的裂缝属性体,所述裂缝属性体为反映断裂及裂缝的属性体;
基于所述第一敏感属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据;
从所述印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从所述第二敏感属性体中提取所述深层碳酸盐岩地层中目的层段的裂缝平面属性,得到所述目的层段的裂缝预测结果。
2.根据权利要求1所述的深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述从所述增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,具体包括:
基于所述增强后叠后地震数据,计算第一预设个数不同的第一裂缝属性体;
分别提取各所述第一裂缝属性体的第一平面属性,并基于所述第一平面属性,从各所述第一裂缝属性体中确定第一敏感属性体。
3.根据权利要求1所述的深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述基于所述第一敏感属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据,具体包括:
对所述第一敏感属性体的地震属性的目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体;
基于所述合并过滤后属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据。
4.根据权利要求3所述的深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述对所述第一敏感属性体的地震属性的目标剖面中的裂缝依次进行细化处理、合并处理和过滤处理,得到合并过滤后属性体,具体包括:
在所述第一敏感属性体的地震属性的目标剖面中,沿代表裂缝发育的属性高值对裂缝进行细化,得到细化后属性体,所述细化后属性体表征与预设的地质认识相一致的裂缝;
将所述细化后属性体的水平切片中属性值高于预设门槛值的裂缝确定为目标裂缝,并将方向一致且裂缝间距小于间距阈值的各所述目标裂缝进行合并,得到合并后裂缝;
对裂缝长度小于长度阈值的合并后裂缝进行过滤处理,得到合并过滤后属性体。
5.根据权利要求3所述的深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述基于所述合并过滤后属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据,具体包括:
确定所述合并过滤后属性体对应的三维地质空间位置;
从所述初始叠后地震数据中选取所述三维地质空间位置对应的裂缝印迹,并按照预设的比例减小公式,对所述裂缝印记进行减小,得到印迹增强地震数据。
6.根据权利要求5所述的深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述比例减小公式为:
;
式中,B表示印迹增强地震数据,A表示初始叠后地震数据,a表示合并过滤后属性体,(0,1)表示合并过滤后属性体正常的属性值范围。
7.根据权利要求1所述的深层碳酸盐岩地层的裂缝预测方法,其特征在于,所述从所述印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,具体包括:
基于所述印迹增强地震数据,计算第二预设个数不同的第二裂缝属性体;
分别提取各所述第二裂缝属性体的第二平面属性,并基于所述第二平面属性,从各所述第二裂缝属性体中确定第二敏感属性体。
8.一种深层碳酸盐岩地层的裂缝预测装置,其特征在于,包括:
数据获取模块(11),用于获取深层碳酸盐岩地层的初始叠后地震数据;
属性优选模块(12),用于对所述初始叠后地震数据进行断裂增强处理,得到增强后叠后地震数据,并从所述增强后叠后地震数据对应的裂缝属性体中确定第一敏感属性体,所述第一敏感属性体为对裂缝最为敏感的裂缝属性体,所述裂缝属性体为反映断裂及裂缝的属性体;
印迹增强模块(13),用于基于所述第一敏感属性体,对所述初始叠后地震数据中对应的裂缝印迹进行增强处理,得到印迹增强地震数据;
裂缝预测模块(14),用于从所述印迹增强地震数据对应的裂缝属性体中确定第二敏感属性体,并从所述第二敏感属性体中提取所述深层碳酸盐岩地层中目的层段的裂缝平面属性,得到所述目的层段的裂缝预测结果。
9.一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机程序,其特征在于,所述计算机程序被处理器加载并执行时,采用了权利要求1-7中任一项所述的方法。
10.一种电子设备,包括存储器、处理器及存储在存储器中并能够在处理器上运行的计算机程序,其特征在于,所述处理器加载并执行计算机程序时,采用了权利要求1-7中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410476754.8A CN118131330A (zh) | 2024-04-19 | 2024-04-19 | 深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410476754.8A CN118131330A (zh) | 2024-04-19 | 2024-04-19 | 深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN118131330A true CN118131330A (zh) | 2024-06-04 |
Family
ID=91235842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410476754.8A Pending CN118131330A (zh) | 2024-04-19 | 2024-04-19 | 深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118131330A (zh) |
-
2024
- 2024-04-19 CN CN202410476754.8A patent/CN118131330A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8358561B2 (en) | Bayesian DHI for seismic data | |
US20190080122A1 (en) | Subsurface reservoir model with 3d natural fractures prediction | |
RU2594373C1 (ru) | Идентификация преобладающих ориентаций трещин | |
EP2697734B1 (en) | Systems and methods for hydraulic fracture characterization using microseismic event data | |
EP2962135B1 (en) | System and method for detecting a fracture in a rock formation using an electromagnetic source | |
WO2012085848A1 (en) | System and method for reconstructing microseismic event statistics from detection limited data | |
CN108254784B (zh) | 一种基于二维地震数据的断层识别方法、装置及系统 | |
CN108375785B (zh) | 裂缝带位置校正方法及装置 | |
EP3264335A1 (en) | Automatic dip picking in borehole images | |
BRPI1002186A2 (pt) | mÉtodo para determinaÇço de abertura dinÂmica para eliminaÇço de méltiplas relacionadas À superfÍcie tridimensional | |
US10705234B2 (en) | Line and edge detection and enhancement | |
CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
US9829591B1 (en) | Determining seismic stratigraphic features using a symmetry attribute | |
CN113970785B (zh) | 暗河裂缝发育预测方法、系统、存储介质以及电子设备 | |
US8340912B2 (en) | Seismic attributes for structural analysis | |
CN118131330A (zh) | 深层碳酸盐岩地层的裂缝预测方法、装置、介质及设备 | |
US11740374B2 (en) | System and method for randomness measurement in sesimic image data using vectorized disorder algorithm | |
CN114876454A (zh) | 一种大斜度井复杂岩性水淹层识别方法及系统 | |
CN111502647A (zh) | 一种钻井地质环境因素确定方法、装置及存储介质 | |
WO2016118223A1 (en) | Adaptive structure-oriented operator | |
CN113138407A (zh) | 一种深层页岩气的多尺度断裂地震预测方法及系统 | |
US20220291405A1 (en) | System and method for storage and retrieval of subsurface rock physical property prediction models using seismic interpretation | |
CN114415230A (zh) | 一种线性断裂提取方法及装置 | |
CN114114390A (zh) | 油气勘探方法、装置、设备以及存储介质 | |
CN118070654A (zh) | 海底浅地层建模方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination |