CN104597485A - 一种微小断层检测方法及断层检测装置 - Google Patents
一种微小断层检测方法及断层检测装置 Download PDFInfo
- Publication number
- CN104597485A CN104597485A CN201310532574.9A CN201310532574A CN104597485A CN 104597485 A CN104597485 A CN 104597485A CN 201310532574 A CN201310532574 A CN 201310532574A CN 104597485 A CN104597485 A CN 104597485A
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- amplitude
- inter
- module
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001514 detection method Methods 0.000 claims abstract description 38
- 238000012937 correction Methods 0.000 claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims description 42
- 238000005070 sampling Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000007781 pre-processing Methods 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 230000001502 supplementing effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000011161 development Methods 0.000 abstract description 4
- 230000001427 coherent effect Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 9
- 230000005012 migration Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种在高分辨率的地震数据中检测微小断层的方法及断层检测装置。所述微小断层检测方法包括:计算装置将地震数据处理成高分辨率的共中心点道集叠加数据;计算装置为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理;计算装置计算所述时间窗内地震数据的道间固有时差;算装置将振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据;所述相邻的两个地震道振幅相减前做道间固有时差校正;计算装置对相邻地震道振幅相减获得的道间差分数据做切片分析。利用本发明,可以实现相干体断层识别技术无法实现的微小断层检测,为小构造油气藏的发现和开发奠定了理论上的基础。
Description
技术领域
本方法涉及地震勘探中地震数据解释领域,特别涉及一种微小断层检测的方法及断层检测装置。
背景技术
为了勘探开发出地层中更多的油气资源,往往需要寻找出更多更隐蔽、更细小的构造油气藏。通常油气资源分布在断层中。断层既是油气藏的边界,也是油气资源运移、聚集的通道,因此断层的识别是油气勘探开发中非常重要的工作。
识别断层的首要步骤是获得高分辨率地震数据。高分辨率地震数据主要通过在野外加大地震数据的采集密度、对采集的地震数据做高分辨率处理等过程获得。
其次是分析上述获得的高分辨率地震数据,通过对地震数据振幅、相位、波形等信息的分析,尽可能在地震数据解释中准确识别断层。
目前在高分辨率地震数据中检测断层的主要方法为相干体断层识别技术。这种方法是利用在出现断层、地层岩性突变、特殊地质体的地带,地震道之间的波形特征发生变化,进而导致局部的道与道之间相关性的发生突变,相干值较小,求得所有样点与其周边样点的相关值形成一个相干数据体,通过对相干数据体作水平切片来揭示断层的地质现象。相干体断层识别技术的实施过程主要包括三个步骤,第一步是将获得的高分辨率三维地震数据做偏移处理。偏移处理能够提高地震剖面分辨率,恢复子波波形和振幅特征,有利于地震数据解释。第二步是在经过偏移获得的三维数据体中对每一个样点求得与周围样点数据的相干值,形成一个相干数据体。第三步是对第二步形成的相干数据体做切片分析,进行构造和岩性解释。这种方法能够识别落差较大的断层,在地震数据信噪比和分辨率较高的情况下也能识别出断层厚度较小的断层。通常断层厚度小于1/4地震波的波长时(通常小于10米),该方法就很难识别断层,尤其是当对断距不超过5米的微小断层的则无法识别。当断层厚度小于1/8地震波的波长时则相干体断层识别技术就无法识别断层了。
相干体断层识别技术是利用地震数据的波形特性的相干性去识别断层,而微小断层两侧波形、振幅、相位等特征并没有发生显著的特征改变。即使有微小断层存在,微小断层处地震波特性的相干值变化不大,尤其是对断层厚度小于1/8地震波波长的微小断层则找不到微小断层的特征,也就无法识别微小断层,影响了微小断层所在的油气藏的勘探和开发。
发明内容
本发明目的在于提供一种地震资料中微小断层的检测方法及断层检测装置,把共中心点道集叠加的高分辨率地震数据进行振幅调整处理,然后将相邻地震道振幅相减获得道间差分数据,从道间差分数据体的切片上找到微小断层的特征。本发明解决了相干体断层识别技术通过波形相干性无法识别微小断层特性的技术问题。
本发明的实现方法包括以下处理步骤:
S1:计算装置将地震数据处理成高分辨率的共中心点道集叠加数据;
S2:计算装置为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理;
S3:计算装置计算所述时间窗内地震数据的道间固有时差;
S4:计算装置将振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据;所述相邻的两个地震道振幅相减前做道间固有时差校正;
S5:计算装置对相邻地震道振幅相减获得的道间差分数据做切片分析。
上述所述的微小断层检测方法,其优选方案为,还包括为S4中所述振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据设定一个门槛值,将道间差分数据中振幅绝对值小于门槛值的数据用零代替。
上述所述的微小断层检测方法,其优选方案为,在S2中所述将地震数据的振幅调整到同一振幅级别,S2中对所述时间窗内的地震数据做振幅均衡处理,包括以下处理步骤:
S301:在所述的时间窗内,将地震记录上的总共N个地震道平均分成K段,每段内有L个地震道,其中N,K,L为正整数;
如果S301中N个地震道不能平均分成K段,还包括进行如下处理:
将N个地震道加上使其可以平均分成K段的P个地震道数,加上的P个地震道数在分成K段的任意一段内进行删、补,其中N、K为正整数,P为整数。
S302:计算S301中平均分成的K段每段内L个地震道内所有样点振幅绝对值的平均值,并以该振幅绝对值的平均值的倒数作为该段的权系数;
S303:将每段内样点的振幅值与其相应段的权系数相乘,并将乘积结果作为该样点的新振幅值。
上述所述的微小断层检测方法,其优选方案为,S1中所述的高分辨率的共中心点道集叠加数据至少经过以下处理步骤:静校正,反褶积运算,动校正,叠加。
上述所述的微小断层检测方法,其优选方案为,S2中所述选定的时间窗T范围为:S≤T≤N,其中S表示该地震数据记录中包含需要识别的微小断层的最小时间跨度,N为该地震数据记录的总时间长度,T表示所选定的时间窗宽度。
本发明提供使用本发明方法进行断层检测的断层检测装置,所述断层检测装置,包括预处理模块,振幅均衡模块,道间固有时差计算模块,相邻地震道相减模块,切片分析模块;其中,
预处理模块,用于将地震数据处理成高分辨率的共中心点道集叠加数据;
振幅均衡模块,用于为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理;
道间固有时差计算模块,用于计算地震数据的道间固有时差;
相邻地震道相减模块,用于振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据;并在相减前对所述相邻的两个地震道做道间固有时差校正;
切片分析模块,用于对所述相邻地震道相减获得的道间差分数据做切片分析。
上述所述的断层检测装置,其优选方案为,所述振幅均衡模块还包括时间窗获得模块,均等分模块,权系数计算模块,振幅加权模块;其中,
时间窗获得模块,用于获取时间窗信息;
均等分模块,用于将地震数据记录上的N个地震道平均分成K等份,每段L个地震道,N、K、L为正整数;
权系数计算模块,用于计算每段内L个地震道内所有样点振幅绝对值的平均值,并以该振幅绝对值的平均值的倒数作为该段的权系数;
振幅加权模块,用于将每段内样点的振幅值与其相应段的权系数相乘,并将乘积结果作为该样点的新振幅值。
上述所述的断层检测装置,其优选方案为,所述相邻地震道相减模块还包括固有时差校正模块,振幅相减模块;其中,
固有时差校正模块,用于根据道间固有时差对相减的两个相邻地震道数据做道间固有时差校正;
振幅相减模块,用于将相邻地震道振幅相减。
本发明将高分辨率的共中心点道集叠加数据做振幅均衡处理,然后将相邻地震道振幅相减获得道间差分数据。利用道间差分数据的振幅对断层两侧地震波的相位差非常敏感的特性,将微小断层的特征在地震数据中突显出来。再对道间差分数据做切片分析,通过微小断层的特征结合地震解释资料以及当地实际地质情况识别出微小断层。
附图说明
图1是本发明实施例1中室内模拟室外的含微小断层的高分辨率共中心点道集叠加数据示意图;
图2是本发明实施例1中图1相邻地震道振幅相减后的道间差分数据示意图;
图3是本发明实施例1中图2用门槛值处理后的道间差分数据示意图;
图4是本发明提供的一种微小断层检测方法的流程图;
图5是本发明实施例2提供的一种断层检测装置的模块结构示意图;
图6是本发明实施例2提供的一种断层检测装置中振幅均衡模块结构示意图;
图7是本发明实施例2提供的一种断层检测装置中相邻地震道相减模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都应当属于本发明保护的范围。
实施例1是本发明提供的一种在高分辨率地震数据背景下检测微小断层的方法,图4是本发明方法检测微小断层的流程图。如图4所示,一种微小断层检测方法,其实施方式包括以下处理步骤:
S1:计算装置将地震数据处理成高分辨率的共中心点道集叠加数据。
计算装置对高密度数据采集、高分辨率处理后获得的高分辨率地震数据进行共中心点道集叠加,形成高分辨率的共中心点道集叠加数据。计算装置可以先对采集后的地震数据做高分辨率处理,然后进行共中心点道集叠加。计算装置也可以先对采集的地震数据进行共中心点道集叠加后再进行高分辨率处理。高分辨率处理是指地震数据处理中一系列消除干扰、拓宽有效信号频谱、压缩地震子波等提高采集的地震数据分辨率的处理过程,其主要过程包括:静校正、反褶积运算、动校正、叠加等。例如所述的静校正是按静校正量所确定的移动量,把地震道的振幅离散值进行整体移动,消除地面山谷、地表风化层和低速带等地表异常的影响。检波器接收到的地震信号可以看成是反射系数序列与地震子波的褶积,子波的存在会降低分辨率,反褶积是把接收到的地震信号进行子波压缩,最后在地震道上仅保留反射系数,提高地震资料分辨率。动校正是通过动校正公式计算出动校正量来消除地震波到达各不同检波点的正常时差。叠加是将经动、静校正处理后的各地震道上序号相同的采样值取算术平均值,组成叠加输出道,目的是压制多次波和随机干扰。
在地震数据采集中,每个检测点上的地震记录称为地震道或道。在地震数据处理时,通常根据需要把地震数据按照一定规律重新排列组成新的地震数据集合。本发明中是把地震数据记录中具有公共中心点的道抽取出来,重新排列形成共中心点道集,也称为CMP道集。上述CMP道集中的公共中心点指的是炮点与其相应的检测点的中点。
计算装置对采集的地震数据完成高分辨率处理、CMP道集叠加,形成高分辨率的CMP道集叠加数据。
S2:计算装置为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理。
上述所述选定的时间窗是限定所需处理数据的范围,该时间窗可以为该地震数据记录的总时间长度,即该地震数据记录中的所有数据参与计算。所述时间窗也可以是至少把需要识别的微小断层包括在内的时间窗。本实施例中以时间窗是该地震数据记录的总时间长度为例进行说明。
地震波的振幅随着传播距离和传播时间的增大而不断衰减,同时又随着频率的不同而变化。在地震数据采集中,检波器到震源激发点的距离为炮检距。同一张地震数据记录中炮检距较小的地震道对应的振幅较大,炮检距较大的距地震道对应的振幅较小,因此同一张记录上的各地震道之间的能量是不均衡的。虽然在地震资料处理时进行了振幅补偿校正,但是各地震道之间还会存在残余能量不均衡的情况。振幅均衡处理就是把各道之间残余能量不均衡的地震道的误差控制在允许范围内,同时可以将地震记录上不同地震道的振幅调整到同一个级别,将强振幅相对减弱,将弱振幅相对放大,缩小振幅横向变化的差距,提高识别微小断层的准确度。
所述振幅均衡处理的方法,可以采用平均振幅法,即取地震道中数据采样点与周边数据采样点振幅的平均值作为采样点的新振幅值。所述振幅均衡处理的方法也可以采用加权系数法。本实施例中振幅均衡处理以加权系数法为例进行说明。所述加权系数法的处理步骤可以包括以下步骤:第一步,在上述所选定的时间窗内,把地震记录上的总共N个地震道平均分成K段,每段包含L个地震道,N,K,L均为正整数。若N个地震道不能平均分成K段,则将N个地震道加上使其可以平均分成K段的P个地震道数后再平均分成K段,与N个地震道相差的P个地震道数可以在分成K段的任意一段内删、补,P为整数。一般加上P个地震道数在最后一段地震道进行减去或增补。第二步,计算每段内L个地震道内所有样点振幅绝对值的平均值,并以该振幅绝对值的平均值的倒数作为该段的权系数。第三步,把每段地震道内的所有样点都进行振幅加权,即把每段内的样点振幅值与其对应的权系数的乘积作为相应样点的新振幅值。完成上述所选时间窗内所有样点的振幅加权即完成了上述所选时间窗内的数据振幅均衡处理。
这样,计算装置对叠加后的共中心点道集地震数据,在选定的时间窗内完成不同地震道的振幅调整处理。
S3:计算装置计算所述时间窗内地震数据的道间固有时差。
在地震数据记录中,由于炮检距不同,同一地震波到达不同检波点的时间也不同,相邻两个检波点接收到同一地震波时间的差值称为道间时差。通常情况下同一条测线上检波点之间是等间距的,所以道间时差不大,当道间时差一致的时候这个差值被称为道间固有时差。道间固有时差不为零时,在下一步相邻地震道相减前要做道间固有时差校正,目的是使相减的相邻地震道振幅样点是同一个同相轴的样点,这样可以获得准确的道间差分数据。
道间固有时差的获得,通常是在上述S2所述时间窗内的同相轴的一个方向上从左到右或从右到左依次扫描出每两个相邻地震道的道间时差,依次记为t1,t2,t3,…,tk。将所有道间时差按照降序或升序的顺序排列,然后取所有道间时差的中值(若k为偶数,则取中间两个道间时差的平均值),这个中值就被认为是沿这个同相轴方向上的道间固有时差。地震记录上各道振动相位相同的振幅极值(俗称波峰或波谷)的连线称为同相轴。
S4:计算装置将振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据;所述相邻的两个地震道振幅相减前做道间固有时差校正。
相邻的两个地震道振幅相减前若道间固有时差不为零,则将相减的两个相邻地震道做固有时差校正后再相减。若道间固有时差为零,此时道间固有时差的校正量为零,实质上未做道间固有时差校正。其校正方法为:设两个相邻地震道的道间固有时差为T,相减前,若T为正值,相减地震道的数据整体向上(时间值减小的方向)移动T毫秒;与此相反,相减前,若T为负值,则被减地震道的数据整体向上(时间值减小的方向)移动T毫秒。
相邻地震道振幅可以依次相减,可以是第一道振幅减去第二道振幅形成第一个差分地震道,第二道振幅减去第三道振幅形成第二个差分地震道,…,依次相邻地震道振幅相减直到整个地震数据的倒数第二道振幅减去最后一道振幅形成最后一个差分地震道为止。也可以是第二道振幅减去第一道振幅形成第一个差分地震道,第三道振幅减去第二道振幅形成第二个差分地震道,…,依次相邻地震道振幅相减直到整个地震数据的最后一道振幅减去倒数第二道振幅形成最后一个差分地震道为止。相邻地震道振幅依次相减完成后获得道间差分数据。
道间差分数据的获得是相邻地震道的振幅相减的过程,也就是逐个样点计算相邻地震道的振幅差值。地震波在跨越断层时,断层点振幅、相位会发送突变,所以相邻地震道波形特征就不一样。在道间固有时差校正后,如果相邻地震道波形相同,则道间差分为零;如果相邻地震道波形不相同,则道间差分不为零。虽然相邻地震道振幅、相位变化很小,但相减后道间差分数据会有明显的表现。设两个相邻地震道内的地震谐波为A(t)和B(t),A(t)和B(t)具有相位差φ,最大振幅值为1,则两个地震谐波可以表示为:
A(t)=sinωt,B(t)=sin(ωt-φ),t是传播时间,ω为周期,φ为相位差。
二者相减得到谐波差分,也就是道间差分:
D(t)=A(t)-B(t)=sinωt-sin(ωt-φ)=2sin(φ/2)cos(ωt-φ/2)。
由上式此可知:道间差分D(t)与A(t)和B(t)是同频率的谐波,其幅度为2sin(φ/2),与A(t)和B(t)的相位差有关,谐波的频率不变,但是正弦函数变成了余弦函数。当没有遇到断层时φ=0,D(t)的幅度是0。如果落差为谐波的1/12周期(即地震波的1/12波长),即φ=π/6,则D(t)的幅度为0.518,此时的振幅差值已经可以被所属技术领域人员识别出来。如果落差为1/6周期,即φ=π/3,则D(t)的幅度达到1。两个谐波的相位差越大,D(t)的幅度就越大,道间差分数据就越明显。
将幅调整到同一振幅级别的地震数据中每两个相邻地震道的振幅依次相减,获得道间差分数据。
此时可以为获得的道间差分数据设定一个门槛值,将道间差分数据中振幅绝对值小于门槛值的数据用零代替。
由于计算误差或者其它干扰导致相邻地震道相减后振幅不为零,道间差分数据存在干扰。为了排除这些不是由微小断层引起的非零道间差分数据的影响,增加微小断层特性的显著性,可以为道间差分数据设置一个振幅门槛值,将门槛值以下的数据用零代替。该门槛值的选定可以是一个分选处理的过程,地震解释人员可以先选定若干预选门槛值进行门槛值处理,将处理后的地震数据进行比较,选择最能显出微小断层特性的门槛值。地震解释人员也可以根据对数据处理要求自行设定门槛值大小。
S5:计算装置对相邻地震道振幅相减获得的道间差分数据做切片分析。
计算装置对S4相邻地震道振幅相减获得的道间差分数据做切片分析,本领域技术人员可以通过切片上微小断层振幅特性体现出的同相轴找到微小断层。切片是沿着三维地震数据体的一个方向,以平面或曲面抽取的具有地质物理意义的二维空间的数据体。常用切片分析方法为水平切片。水平切片是一种等时平面图,是地下不同层位的地质信息在同一时间平面内的反映。断层在水平切片上的反映主要表现在振幅呈规律性分布。由断层形成的差分数据在水平切片上主要表现为规律性出现的一段同相轴。而使用本发明方法后,1/12地震波波长的微小断层所体现出来的振幅特性在切片上表现出来的同相轴已经可以辨识,再结合其它地质资料和分析方法可以识别相干体断层识别技术无法找到的微小断层。
对S4获得的道间差分数据做切片分析,根据切片上的微小断层体现出来的特征识别出微小断层。
图1至图3是本发明使用理论模拟数据突出微小断层特征的过程。
首先,室内模拟一段野外采集的地震数据,对该地震数据进行高分辨率处理并进行共中心点道集叠加,形成一段高分辨率的共中心点道集叠加数据,如图1所示。图1中横轴为地震道号,对每个地震道进行的编号为地震道号,如图1中32到172。纵轴为记录地震道的时间,单位为毫秒。其次,本实施1中将所有地震数据参与计算,因此时间窗选择整个地震数据的最大时间长度300毫秒。对300毫秒长的时间窗内的地震数据进行振幅均衡处理。第三步,根据振幅的能量在这个300毫秒时间窗内扫描出一个同相轴的道间时差,扫描方式为同相轴上右边的地震道减去左边的地震道获得该相邻地震道的道间时差。获得这个时间窗内的地震数据道间时差的中值为6毫秒,这个6毫秒就被认为是沿这个同相轴方向的道间固有时差。第四步,进行相邻地震道振幅依次相减。用第一道振幅减去第二道振幅前,第二个地震道内的数据向上(也就是时间值减小的方向)移动6毫秒,完成道间时间差校正,然后第一道振幅再减去第二道移动后的振幅,得到第一个道间差分数据。接着第三道地震内的数据向上移动6毫秒,用移动前的第二道振幅减去移动后第三道振幅得到第二个道间差分数据,依次相减直到道间固有时差校正后最后一道减去固有时差校正前的倒数第二道振幅得到最后一个道间差分数据,这样下去就完成了道间差分处理,获得道间差分数据,如附图2所示。本实施例1门槛值设置为0.4,将差分数据振幅绝对值小于0.4的数据用零代替,形成门槛值处理后的道间差分数据,如图3所示。由图3可以看出,在92号地震道附近微小断层的振幅特性在道间差分数据上已经体现的很明显,最后再对其做切片分析就可以辨识出微小断层。
本发明方法通过模拟数据的分析、理论的推理证明了本发明方法可以识别出相干体断层识别技术无法识别的微小断层的特征,进行识别出微小断层,为微小断层的油气藏开发奠定了理论基础。
实施例2为本发明提供的一种使用本发明方法进行断层检测的装置。图5示出了本发明装置的模块结构示意图。所述断层检测装置包括预处理模块1,振幅均衡模块2,道间固有时差计算模块3,相邻地震道相减模块4,切片分析模块5。其中,
预处理模块1,可以用于将地震数据处理成高分辨率的共中心点道集叠加数据。
振幅均衡模块2,可以用于为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理。
图6是本实施例提供的一种振幅均衡模块结构示意图。如图6所示,该振幅均衡模块2包括时间窗获得模块201,均等分模块202,权系数计算模块203,振幅加权模204。其中,
时间窗获得模块201,可以用于获得操作人员输入的时间窗信息;
均等分模块202,可以用于将地震数据记录上的N个地震道平均分成K等份,每一份内有L个地震道。若N个地震道不能平均分成K等份,可以将N个地震道加上P个地震道数后使其可以平均分成K等份。相差的P个地震道数在最后一段地震道内删、补。所述N、K、L为正整数,P为整数;
权系数计算模块203,可以用于计算每段内L个地震道内所有样点振幅绝对值的平均值,并以该振幅绝对值的平均值的倒数作为该段的权系数;
振幅加权模块204,可以用于将每段内样点的振幅值与其相应段的权系数相乘,并将乘积结果作为该样点的新振幅值。
道间固有时差计算模块3,可以用于计算振幅均衡模块3中选定时间窗内地震数据的道间固有时差。该模块计算地震数据同相轴的任意一个方向上所有相邻地震道的道间时差,将道间时差按照升序或降序排列,然后取所有道间时差的中值为该时间窗内数据的道间固有时差。
相邻地震道相减模块4,可以用于将相邻地震道的振幅相减,并在相减前对相邻的两个地震道做道间固有时差校正。
相邻地震道相减模块4,可以将相邻地震道的振幅相减,获得道间差分数据。图7是本实施例提供的一种相邻地震道相减模块结构示意图。如图7所示,相邻地震道模块分为固有时差校正模块401,振幅相减模块402。其中,
固有时差校正模块401,可以用于根据道间固有时差对相减的两个相邻地震道数据做道间固有时差校正;
振幅相减模块402,可以用于将相邻地震道振幅相减,获得相应的道间差分数据。
切片分析模块5,可以用于对相邻地震道相减获得的道间差分数据做切片分析。
所述的使用本发明方法的断层检测装置,可以显示出1/12地震波波长的微小断层体现出的振幅特性,再通过对差分数据的切片分析进一步识别出微小断层。该断层检测装置实现了一般利用相干体断层识别技术的装置无法实现的微小断层检测。
Claims (9)
1.一种微小断层检测方法,其特征在于,包括以下处理步骤:
S1:计算装置将地震数据处理成高分辨率的共中心点道集叠加数据;
S2:计算装置为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理;
S3:计算装置计算所述时间窗内地震数据的道间固有时差;
S4:计算装置将振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据;所述相邻的两个地震道振幅相减前做道间固有时差校正;
S5:计算装置对相邻地震道振幅相减获得的道间差分数据做切片分析。
2.如权利要求1所述的微小断层检测方法,其特征在于,还包括为S4中所述振幅均衡处理后地震数据的每两个相邻地震道的振幅相减获得的道间差分数据设定一个门槛值,将所述道间差分数据中振幅绝对值小于门槛值的数据用零代替。
3.如权利要求1或2所述的微小断层检测方法,其特征在于,S2中对所述时间窗内的地震数据做振幅均衡处理,包括以下处理步骤:
S301:在所述的时间窗内,将地震记录上的总共N个地震道平均分成K段,每段内有L个地震道,其中N,K,L为正整数;
S302:计算S301中平均分成的K段每段内L个地震道内所有样点振幅绝对值的平均值,并以该振幅绝对值的平均值的倒数作为该段的权系数;
S303:将每段内样点的振幅值与其相应段的权系数相乘,并将乘积结果作为该样点的新振幅值。
4.如权利要求3所述的微小断层检测方法,其特征在于,如果S301中N个地震道不能平均分成K段,还包括进行如下处理:
将N个地震道加上使其可以平均分成K段的P个地震道数,加上的P个地震道数在分成K段的任意一段内进行删、补,其中N、K为正整数,P为整数。
5.如权利要求1或2所述的微小断层检测方法,其特征在于,S1中所述的高分辨率的共中心点道集叠加数据至少经过以下处理步骤:静校正,反褶积运算,动校正,叠加。
6.如权利要求1或2所述的微小断层检测方法,其特征在于,S2中所述选定的时间窗T的范围为:S≤T≤N,其中S表示该地震数据记录中包含需要识别的微小断层的最小时间跨度,N为该地震数据记录的总时间长度,T表示所选定的时间窗宽度。
7.一种断层检测装置,其特征在于,所述断层检测装置,包括预处理模块,振幅均衡模块,道间固有时差计算模块,相邻地震道相减模块,切片分析模块;其中,
预处理模块,用于将地震数据处理成高分辨率的共中心点道集叠加数据;
振幅均衡模块,用于为高分辨率的共中心点道集叠加数据选定一个时间窗,对所述时间窗内的地震数据做振幅均衡处理;
道间固有时差计算模块,用于计算地震数据的道间固有时差;
相邻地震道相减模块,用于振幅均衡处理后地震数据的每两个相邻地震道的振幅相减,获得道间差分数据;并在相减前对所述相邻的两个地震道做道间固有时差校正;
切片分析模块,用于对所述相邻地震道相减获得的道间差分数据做切片分析。
8.如权利要求7所述的断层检测装置,其特征在于,所述的振幅均衡模块包括时间窗获得模块,均等分模块,权系数计算模块,振幅加权模块;其中,
时间窗获得模块,用于获取时间窗信息;
均等分模块,用于将地震数据记录上的N个地震道平均分成K等份,每段L个地震道,N、K、L为正整数;
权系数计算模块,用于计算每段内L个地震道内所有样点振幅绝对值的平均值,并以该振幅绝对值的平均值的倒数作为该段的权系数;
振幅加权模块,用于将每段内样点的振幅值与其相应段的权系数相乘,并将乘积结果作为该样点的新振幅值。
9.如权利要求7所述的断层检测装置,其特征在于,所述相邻地震道相减模块包括固有时差校正模块,振幅相减模块;其中,
固有时差校正模块,用于根据道间固有时差对相减的两个相邻地震道数据做道间固有时差校正;
振幅相减模块,用于将相邻地震道振幅相减。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310532574.9A CN104597485B (zh) | 2013-10-31 | 2013-10-31 | 一种微小断层检测方法及断层检测装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310532574.9A CN104597485B (zh) | 2013-10-31 | 2013-10-31 | 一种微小断层检测方法及断层检测装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104597485A true CN104597485A (zh) | 2015-05-06 |
CN104597485B CN104597485B (zh) | 2017-03-08 |
Family
ID=53123390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310532574.9A Active CN104597485B (zh) | 2013-10-31 | 2013-10-31 | 一种微小断层检测方法及断层检测装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104597485B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873030A (zh) * | 2017-02-07 | 2017-06-20 | 中国石油天然气股份有限公司 | 断层封堵性的确定方法和装置 |
CN107843919A (zh) * | 2017-10-27 | 2018-03-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种微地震监测数据振幅均衡方法 |
CN110646848A (zh) * | 2019-10-15 | 2020-01-03 | 中国石油集团东方地球物理勘探有限责任公司 | 一种断层的地震解释方法 |
CN113050161A (zh) * | 2021-03-29 | 2021-06-29 | 中国矿业大学 | 一种基于断层交面线的断裂走向断距图解计算方法 |
CN113253345A (zh) * | 2020-02-10 | 2021-08-13 | 中国石油天然气股份有限公司 | 相位差异属性获取方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1460339A1 (ru) * | 1987-04-13 | 1989-02-23 | В. С. Маевский, М. И. Большинский, О. Г. Р ременев и Л. А. Вайнщтейн | Способ определени местоположени мелкоамплитудного геологического нарушени при разработке газоносного угольного пласта |
-
2013
- 2013-10-31 CN CN201310532574.9A patent/CN104597485B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
Non-Patent Citations (3)
Title |
---|
侯水平等: "《利用相对振幅曲线解释小断层的应用》", 《河北煤炭》 * |
杨文强等: "《时间剖面上小断层的计算机自动检测》", 《煤田地质与勘探》 * |
王彦君等: "《小断层识别技术研究及应用》", 《勘探地球物理进展》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873030A (zh) * | 2017-02-07 | 2017-06-20 | 中国石油天然气股份有限公司 | 断层封堵性的确定方法和装置 |
CN107843919A (zh) * | 2017-10-27 | 2018-03-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种微地震监测数据振幅均衡方法 |
CN107843919B (zh) * | 2017-10-27 | 2020-02-07 | 中国石油天然气集团有限公司 | 一种微地震监测数据振幅均衡方法 |
CN110646848A (zh) * | 2019-10-15 | 2020-01-03 | 中国石油集团东方地球物理勘探有限责任公司 | 一种断层的地震解释方法 |
CN113253345A (zh) * | 2020-02-10 | 2021-08-13 | 中国石油天然气股份有限公司 | 相位差异属性获取方法及系统 |
CN113050161A (zh) * | 2021-03-29 | 2021-06-29 | 中国矿业大学 | 一种基于断层交面线的断裂走向断距图解计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104597485B (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100501449C (zh) | 一种地震数据拆分处理方法 | |
Aghaei et al. | Crustal thickness and Moho character of the fast‐spreading East Pacific Rise from 9° 42′ N to 9° 57′ N from poststack‐migrated 3‐D MCS data | |
CN103954992B (zh) | 一种反褶积方法及装置 | |
US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
CN104280777B (zh) | 一种压制陆上地震资料多次波干扰的方法 | |
CN104597485B (zh) | 一种微小断层检测方法及断层检测装置 | |
Al-Hagan et al. | Iterative supervirtual refraction interferometry | |
CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
WO2015104637A2 (en) | Device and method for deghosting seismic data using sparse tau-p inversion | |
CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
Colombo et al. | Near-surface full-waveform inversion in a transmission surface-consistent scheme | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
Duret et al. | Near-surface velocity modeling using a combined inversion of surface and refracted P-waves | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
Qiu et al. | Denoising surface waves extracted from ambient noise recorded by 1‐D linear array using three‐station interferometry of direct waves | |
Ikeda et al. | Advanced surface-wave analysis for 3D ocean bottom cable data to detect localized heterogeneity in shallow geological formation of a CO2 storage site | |
EA030770B1 (ru) | Система и способ адаптивной сейсмической оптики | |
Bouchaala et al. | Azimuthal Investigation of a Fractured Carbonate Reservoir | |
Oliveira et al. | Estimation of quality factor based on peak frequency-shift method and redatuming operator: Application in real data set | |
Zhao et al. | Target-oriented diversity stacking for virtual-source imaging and monitoring | |
Shekar et al. | Structural information derived from ambient noise tomography over a hydrocarbon‐producing region in the Cachar fold belt, lower Assam, northeast India | |
Almutlaq et al. | Surface-consistent matching filters for time-lapse seismic processing | |
Santos et al. | Diffraction velocity analysis in a single-channel seismic survey in the Joetsu Basin | |
CN103760599B (zh) | 一种微小断层检测方法及断层检测装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |