一种断层破碎带宽度的确定方法及装置转让专利
申请号 : CN202010640077.0
文献号 : CN111856572B
文献日 : 2021-07-20
发明人 : 廖宗湖 , 陈硕
申请人 : 中国石油大学(北京)
摘要 :
权利要求 :
1.一种断层破碎带宽度的确定方法,其特征在于,所述方法包括:获取目标地层的平面方差属性图;
将所述平面方差属性图中,断层破碎带外方差值为预设值的位置作为基线位置;所述基线平行于所述断层破碎带的延伸方向;
以所述基线为起点作多条垂直于所述断层破碎带的测线;其中,每条测线上包括多个间隔距离相等的测点;
根据每条测线上测点的方差值随测点至基线的距离的变化,确定断层破碎带宽度;
其中,所述根据每条测线上测点的方差值随测点至基线的距离的变化,确定断层破碎带宽度,包括:
根据每条测线上测点的方差值随测点至基线的距离的变化区分围岩和断层破碎带的边界;
根据所述边界划分断层破碎带,并确定所述断层破碎带的宽度;
其中,所述根据每条测线上测点的方差值随测点至基线的距离的变化区分围岩和断层破碎带的边界,包括:
绘制每条测线上测点的方差值随测点至基线的距离的变化曲线;
将所述变化曲线中方差值的变化由缓慢增加至快速增加的点确定为围岩和断层破碎带的边界点;
将所述变化曲线中方差值等于所述边界点方差值的另一个点作为另一个边界点,其中,另一个边界点为方差值的变化趋于平缓降低的点。
2.根据权利要求1所述的方法,其特征在于,所述断层破碎带包括单核断层破碎带和多核断层破碎带。
3.根据权利要求1所述的方法,其特征在于,所述获取目标地层的平面方差属性图包括:
将目标地层的地震道数据导入petrel软件中;
根据所述petrel软件的输出结果获取目标地层的平面方差属性图。
4.根据权利要求1所述的方法,其特征在于,所述预设值为0‑0.1。
5.根据权利要求1所述的方法,其特征在于,所述每条测线与所述断层破碎带具有交点。
6.根据权利要求1所述的方法,其特征在于,各个测线之间在所述目标地层中的实际间隔距离为50‑200米。
7.根据权利要求1所述的方法,其特征在于,每条测线上各个测点在所述目标地层中的实际间隔距离为30‑60米。
8.一种断层破碎带宽度的确定装置,其特征在于,所述装置包括:获取模块,用于获取目标地层的平面方差属性图;
基线确定模块,用于将所述平面方差属性图中,断层破碎带外方差值为预设值的位置作为基线位置;所述基线平行于所述断层破碎带的延伸方向;
测线确定模块,用于以所述基线为起点作多条垂直于所述断层破碎带的测线;其中,每条测线上包括多个间隔距离相等的测点;
断层破碎带宽度确定模块,用于根据每条测线上测点的方差值随测点至基线的距离的变化,确定断层破碎带宽度;
其中,所述断层破碎带宽度确定模块包括:区分单元,用于根据每条测线上测点的方差值随测点至基线的距离的变化区分围岩和断层破碎带的边界;
划分单元,用于根据所述边界划分断层破碎带,并确定所述断层破碎带的宽度;
其中,所述区分单元具体用于:绘制每条测线上测点的方差值随测点至基线的距离的变化曲线;
将所述变化曲线中方差值的变化由缓慢增加至快速增加的点确定为围岩和断层破碎带的边界点;
将所述变化曲线中方差值等于所述边界点方差值的另一个点作为另一个边界点,其中,另一个边界点为方差值的变化趋于平缓降低的点。
说明书 :
一种断层破碎带宽度的确定方法及装置
技术领域
背景技术
宽,与断层的规模和力学性质有关。断层控制了断块油田的形成及其油气开发,特别是断层
破碎带的组成、厚度、物性和封闭性直接影响剩余油开发的方式,因此油田内部断层破碎带
的识别和描述成为亟待解决的科学和实际问题。
破碎带的外边缘(破碎带和围岩的边界)可以视为裂缝发育的程度下降到一定的值(背景
值)的位置。这种基于裂缝密度来量化断层破碎带宽度研究方法为定量地识别断层破碎带
的区域提供了有用的思路,在很多涉及断层破碎带的宽度的研究中得到了一定的应用。例
如,Brogi(Brogi A.Fault zone architecture and permeability features in
siliceous sedimentary rocks:Insights from the Rapolano geothermal area
(Northern Apennines,Italy)[J].Journal of Structural Geology,2008,30(2):237‑
256)和Kristensen(Kristensen T B,RotevatnA,Peacock D C P,et al.Structure and
flow properties of syn‑rift border faults:The interplay between fault damage
and fault‑related chemical alteration(Dombjerg Fault,Wollaston Forland,NE
Greenland)[J].Journal of Structural Geology,2016,92:99‑115)等人通过断层裂缝发
育程度的变化情况量化了破碎带的宽度,发现断层两盘的破碎带区域呈现不对称性,不仅
破碎带的宽度不相同,破碎强度也不一致。
线,然后测量其裂缝密度的统计情况,但是这种量化结果确切的说并不能是该断层的破碎
带宽度,它既没有体现断层破碎带在横向上的宽度变化,也没体现其在纵向上的宽度变化。
发明内容
为预设值的位置作为基线位置;所述基线平行于所述断层破碎带的延伸方向;以所述基线
为起点作多条垂直于所述断层破碎带的测线;其中,每条测线上包括多个间隔距离相等的
测点;根据每条测线上测点的方差值随测点至基线的距离的变化,确定断层破碎带宽度。
面方差属性图中,断层破碎带外方差值为预设值的位置作为基线位置;所述基线平行于所
述断层破碎带的延伸方向;测线确定模块,用于以所述基线为起点作多条垂直于所述断层
破碎带的测线;其中,每条测线上包括多个间隔距离相等的测点;断层破碎带宽度确定模
块,用于根据每条测线上测点的方差值随测点至基线的距离的变化,确定断层破碎带宽度。
作为基线位置;所述基线平行于所述断层破碎带的延伸方向;以所述基线为起点作多条垂
直于所述断层破碎带的测线;其中,每条测线上包括多个间隔距离相等的测点;根据每条测
线上测点的方差值随测点至基线的距离的变化,确定断层破碎带宽度。本说明书实施例提
供的断层破碎带宽度的确定方法,通过利用方差属性值量化断层破碎带宽度,由于方差属
性能够很好的体现断裂破碎带在三维空间上的变化,且方差属性值是由地震数据中计算得
到,提高了断层破碎带宽度确定的准确性。
附图说明
本说明书中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前
提下,还可以根据这些附图获得其他的附图。
具体实施方式
例。基于本说明书中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得
的所有其他实施例,都应当属于本说明书保护的范围。
可以通过测量、计算或者其他方式从地震数据中提取。基本的地震数据类型是振幅,而地震
属性可以获得隐藏在地震数据中的信息,揭示在振幅数据中不易观察到的地质现象,提高
地层的地质解释能力。目前前人已经定义了不少对断层和裂缝表征十分有用的地震属性,
比如倾角和方位角属性、曲率属性、相干属性、方差属性、相似属性等。这些属性很大程度上
取决于地下地质体的几何和构造特征,比如地震事件的形状、角度和连续性,因此可以广泛
应用于地震数据中的地质构造解释,比如断层、褶皱、裂缝网络等。在属性分析中,没有一个
属性是可以绝对正确并有效地对地下的复杂构造进行成像,因此选择正确合适的属性对地
下构造进行综合分析十分重要。
测量其裂缝密度的统计情况,但是这种量化结果确切的说并不能是该断层的破碎带宽度,
它既没有体现断层破碎带在横向上的宽度变化,也没体现其在纵向上的宽度变化。由于断
层破碎带为一个高应变区域,其结构特征十分复杂,基于裂缝密度来量化断层破碎带宽度
研究方法不能准确地体现断裂破碎带在三维空间上的变化。考虑到如果选择合适的地震属
性来量化断层破碎带宽度,则有望解决现有技术中基于裂缝密度的量化方式来确定断层破
碎带宽度的准确性不高的问题。
育或岩性突变等情况造成地层不均匀、不连续时地震波形具有差异性,通过检测这种差异
性达到提取断层等特殊构造发育信息。方差体技术的具体的运算方法是根据不同的对象选
择合适的时窗长度和采样间隔,然后计算出样点与时窗内相邻地震道的振幅均值,进而求
出方差值,最后经过加权归一化处理将方差值控制在一定的范围内。一般情况下,落差较大
(大于5m)的断层在常规地震振幅剖面上可以呈现出明显的同相轴错断、扭曲、杂乱等特征,
较容易识别。但是落差较小的断层在常规地震剖面上信息不够突出,此时运用方差技术可
将地震资料中较弱的断层信息进行放大,以检测到小断层的发育情况。所以方差属性相较
于其它属性对断层更敏感,可以解释小断层的发育情况,而且对特殊地质体的边界和几何
形态表征更清晰。
宽度,可以提高断层破碎带宽度确定的准确性。基于此,本说明书实施例一种断层破碎带宽
度的确定方法。
2
组数据的方差D 的具体计算方法如公式(1)所示。由公式(1)易知若此组数据越离散(差异
大),则方差值越大,反之越小。
形理论上就会存在比较大的差异,方差属性是通过一定的方式计算这些地震道波形之间方
差来形容道与道之间的差异程度。
按照采样间隔在其上下均取二分之一时窗长度的采样点,即当前采样点在剖面上位于一个
时窗长度的中心位置(如图2所示);计算每个时窗长度内各地震道内采样点的振幅平均值;
基于该振幅平均值算出各地震道在所选时窗长度内的方差之和;加权归一化处理后的单个
采样点的方差值,加权函数通常为三角函数,其目的是将方差值控制在0到1之间。其中,各
地震道采样点方差和加权归一化处理后的各地震道采样点方差根据以下公式计算:
数,j和t分别表示不同的采样点,w为加权函数,wj‑t为某个时窗内某采样点的三角加权函
数,其中0≤wj‑t≤1。
化趋势。
petrel软件,将目标地层的地震道数据导入petrel软件中;根据所述petrel软件的输出结
果获取目标地层的平面方差属性图。其中,所述petrel软件是由Schlumberger公司开发的,
以三维地质模型为中心的勘探开发一体化平台,属于地球物理专业软件。
核,则该断层破碎带为单核断层破碎带;若断层破碎带包括多个断层核,则该断层破碎带为
多核断层破碎带。如图3a和图3b所示,图3a和图3b为目标地层的平面方差属性图,其中,图
3a表示目标地层包括单核断层破碎带,图3b表示目标地层包括多核断层破碎带。当然,目标
地层也可以同时包括单核断层破碎带和多核断层破碎带。
方差低。因此,如果测线上各个测点之间的方差值变化不大,且接近于0,则说明各个测点所
在的区域为围岩区域;如果测线上各个测点之间的方差值变化较大,则说明各个测点处于
围岩与断层破碎带的边界区域;如果测线上各个测点之间的方差值变化不大,且接近于1,
则说明各个测点所在的区域为断层破碎带区域。因此,可以找出围岩与断层破碎带的两个
边界区域,即可根据这两个边界区域确定断层破碎带的宽度。
中,图4a为目标地层中包含单核断层破碎带的情况下,每条测线上测点的方差值随测点至
基线的距离的变化曲线;图4b为目标地层中包含多核断层破碎带的情况下,每条测线上测
点的方差值随测点至基线的距离的变化曲线。
带的边界点,该点对应的方差值为方差边界值。将所述变化曲线中,另一个方差值等于方差
边界值的点作为另一个边界点。其中,另一个边界点为方差值的变化趋于平缓降低的点。将
两个边界点所在的位置确定为围岩和断层破碎带的边界。
距离确定为所述断层破碎带的宽度。
所述基线平行于所述断层破碎带的延伸方向;以所述基线为起点作多条垂直于所述断层破
碎带的测线;其中,每条测线上包括多个间隔距离相等的测点;根据每条测线上测点的方差
值随测点至基线的距离的变化,确定断层破碎带宽度。本说明书实施例提供的断层破碎带
宽度的确定方法,通过利用方差属性值量化断层破碎带宽度,由于方差属性能够很好的体
现断裂破碎带在三维空间上的变化,且方差属性值是由地震数据中计算得到,提高了断层
破碎带宽度确定的准确性。
分次生小断层表现为单断层核形态。如图5,这里仅选择局部断层核形态较清晰的部分进行
本说明书实施例提供的断层破碎带宽度的确定方法量化断层破碎带宽度。其中L1‑L5为单
核断层破碎带的测线,L6‑L8为多核断层破碎带的测线,每条测线上测点的间隔距离为50m。
测点的方差值随测点至基线的距离的变化曲线表现出明显的单峰特征,而多核断层破碎带
每条测线上测点的方差值随测点至基线的距离的变化曲线表现出明显的多峰特征,且越位
于正中心,方差峰值相对而言越高。理论上高峰值的个数代表断层核的个数。随着离破碎最
严重的区域(峰值区)距离增加,破碎强度逐渐减弱(方差值逐渐降低)。
界值进行量化,从图6可以看出,a点为单核断层破碎带的各个测线上测点的方差值随测点
至基线的距离的变化曲线中,方差值的变化由缓慢增加至快速增加的点作为围岩和断层破
碎带的边界点,该点对应的方差值0.2为方差边界值,b点为另一个方差值等于方差边界值
的点作为另一个边界点,且b点为方差值的变化趋于平缓降低的点。将a点和b点的距离确定
为单核断层破碎带宽度,为125m。从图7可以看出,c点为多核断层破碎带的各个测线上测点
的方差值随测点至基线的距离的变化曲线中,方差值的变化由缓慢增加至快速增加的点作
为围岩和断层破碎带的边界点,该点对应的方差值0.2为方差边界值,d点为另一个方差值
等于方差边界值的点作为另一个边界点,且d点为方差值的变化趋于平缓降低的点。将c点
和d点的距离确定为单核断层破碎带宽度,为1000m。
元,用于根据所述边界划分断层破碎带,并确定所述断层破碎带的宽度。
处。尤其,对于装置实施例和设备实施例而言,由于其基本相似于方法实施例,所以描述的
比较简单,相关之处参见方法实施例的部分说明即可。
而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。
设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因
此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件
(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable Gate
Array,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员
自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作
专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻
辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,
而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言
(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL
(Advanced Boolean Expression Language)、AHDL(AlteraHardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language)等,目前最普遍使用的是VHDL(Very‑High‑Speed
Integrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也
应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路
中,就可以很容易得到实现该逻辑方法流程的硬件电路。
以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放
器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何
设备的组合。
本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件
产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机
设备(可以是个人计算机,服务器,或者网络设备等)执行本说明书各个实施例或者实施例
的某些部分所述的方法。
施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例
的部分说明即可。
置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设
备的分布式计算环境等等。
组件、数据结构等等。也可以在分布式计算环境中实践本说明书,在这些分布式计算环境
中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块
可以位于包括存储设备在内的本地和远程计算机存储介质中。
说明书的精神。