CN101201405A - 一种提高地震数据处理分辨率的方法 - Google Patents
一种提高地震数据处理分辨率的方法 Download PDFInfo
- Publication number
- CN101201405A CN101201405A CNA2006101650182A CN200610165018A CN101201405A CN 101201405 A CN101201405 A CN 101201405A CN A2006101650182 A CNA2006101650182 A CN A2006101650182A CN 200610165018 A CN200610165018 A CN 200610165018A CN 101201405 A CN101201405 A CN 101201405A
- Authority
- CN
- China
- Prior art keywords
- deconvolution
- factor
- wavelet
- phase
- resolution
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及物探地震数据处理中采用混合相位反褶积提高地震数据处理分辨率的方法,处理步骤是采集单道地震数据,确定自相关时窗,求子波的自相关;求得各个样点值作为多项式的系数,然后求出子波因子;把所有最小相位子波的反褶积算子倒序排列,得最大,最小两个反褶积算子;再用最小、大相位反褶积算子对记录做反褶积;确定反褶积结果分辨率较高的中间结果;再对每一对因式串联处理,完成混合相位反褶积。本发明避开了子波的最小相位假设,最终实现对非最小相位子波的反褶积,能显著提高地震数据分辨率,且在一定程度上能提高信噪比。
Description
技术领域
本发明涉及地球物理勘探技术,是一种物探地震数据处理中采用混合相位反褶积提高地震数据处理分辨率的方法。
背景技术
在物探地震数据处理技术中,目前使用的反褶积方法是建立在地震子波是最小相位的基础上,利用地震记录的自相关计算反褶积算子,因为自相关不包含相位信息,只有在子波是最小相位假设条件下才能得到唯一解。所以不适应非最小相位的实际地震子波,这就限制了地震数据处理中分辨率的提高。
反褶积是高分辨率处理中的关键环节。目前使用的反褶积方法利用地震记录的自相关代替地震子波的自相关计算反褶积算子,因为自相关并不包含相位信息,只有在子波是最小相位的假设条件下才能得到唯一解。因此这种反褶积方法只适用于最小相位子波,而震源激发的地震子波在传播过程中经过大地的吸收和地震记录仪器的改造,未必是最小相位的。另外在地震数据处理中,振幅补偿可能会进一步破坏地震子波的最小相位性质。因此子波最小相位假设往往是不符合实际情况的,所以这种最小相位假设严重限制地震数据处理中分辨率的提高。
为突破这个限制,长期以来国内外对非最小相位子波的反褶积方法做过不少研究和探索,至今未得到很好解决。
发明内容
本发明目的在于提供一种混合相位反褶积方法,避开子波的最小相位假设,实现对非最小相位子波的反褶积的提高地震数据处理分辨率的方法。
本发明过程包括以下处理步骤:
(1)采集单道地震数据并用常规的地震资料处理,确定地震资料信噪比较高的时间段为自相关时窗,对选择好的时窗数据求子波的自相关;
(2)把求得的自相关各个样点值作为多项式的系数,然后对多项式作因式分解,求出自相关的各个子因式得到子波因子;
(3)用通常的方法通过子波因子求得子波因子的反褶积算子;
(4)把得到的所有最小相位子波的反褶积算子倒序排列,得到每个子波因子最大,最小两个反褶积算子;
(5)用最小相位反褶积算子对记录做反褶积;
(6)用最大相位反褶积算子对记录做反褶积;
(7)对(5),(6)两个步骤的反褶积结果进行分选,确定分辨率较高的中间结果;
(8)确定分辨率较高的中间结果的因式,得到这个因式的反褶积结果;
(9)在前一步骤因式反褶积结果基础上,用步骤(3)至(8)处理步骤(2)的每一对因式,串联进行对因式处理,完成混合相位反褶积;
(10)采用通常的物探处理得到反褶积后提高地震数据分辨率的叠加剖面。
本发明还提供以下技术方案:
所述的子波因式个数在10-20个。
所述的确定分辨率较高的中间结果采用录各样点绝对值的5次方之和。
本发明作因式分解时用单边自相关代替双边自相关,用地震数据自相关代替子波自相关。通过对自相关进行因式分解,得到子波的可能因式,然后对分解所得的子波因式逐个求反褶积算子并进行最大、最小相位筛选,从而避开了子波的最小相位假设,最终实现对非最小相位子波的反褶积。本发明的混合相位反褶积方法能显著提高地震数据分辨率,且在一定程度上能提高信噪比。
附图说明
图1(a)是理论模型子波;
图1(b)是对理论模型子波做最小相位反褶积的第一步;
图1(d)是对理论模型子波做最小相位反褶积的第二步;
图1(f)是对理论模型子波做最小相位反褶积的第三步;
图1(h)是对理论模型子波做最小相位反褶积的第四步;
图1(c)是对理论模型子波做混合相位反褶积的第一步;
图1(e)是对理论模型子波做混合相位反褶积的第二步;
图1(g)是对理论模型子波做混合相位反褶积的第三步;
图1(i)是对理论模型子波做混合相位反褶积的第四步;
图2(a),反射系数模型;
图2(b),合成地震记录,即图2(a)反射系数模型与图1(a)理论模型子波褶积结果;
图2(c),对合成地震记录做最小相位反褶积结果;
图2(d),本发明对合成地震记录做混合相位反褶积结果;
图3(a)未做反褶积的原始炮记录;
图3(b)对图3(a)做最小相位反褶积后结果;
图3(c)本发明对图3(a)做混合相位反褶积后结果;
图3(d)未做反褶积的原始炮记录的叠加剖面;
图3(e)对图3(a)做最小相位反褶积后的叠加剖面;
图3(f)本发明对图3(a)做混合相位反褶积后的叠加剖面。
具体实施方式
本发明的混合相位反褶积方法,如果设有任意N+1个点的子波,这个子波的自相关有2N个因式,其中N个因式就是子波的因式,另外N个因式的根是子波因式的根的倒数。若子波有一个最小相位因式,则自相关也有这个最小相位因式,同时还有一个最大相位因式,它的根是前一个因式的倒数。反之,若子波有一个最大相位因式,则自相关也有这个最大相位因式,同时还有一个最小相位因式,它的根是前一个因式的倒数。
因此,自相关的2N个因式有N个因式是最小相位的,另外有N个因式是最大相位的。但是不能确定那些因式是子波的因式,那些因式不是子波的因式。
如果采用子波最小相位的假设,就采用自相关的2N个因式中的N个最小相位因式当作子波的因式。在此基础上的反褶积方法就是当前使用的最小相位反褶积,它只有子波刚好是最小相位的情况下才是有效的。
本发明将自相关的2N个因式可分成N对,每一对是一个最小相位因式和一个最大相位因式,它们的根互为倒数,其中必有一个且只有一个因式是子波的因式。只要能在每一对因式中选出子波的因式,就可以得到子波的全部N个因式,即得到子波,也就得到反褶积算子。这样,反褶积就不受子波是不是最小相位的限制。
实际上每一对因式只要知道一个因式,另一个因式可以利用根的倒数关系得到。因此只要求出自相关的R(Z)在单位圆外的N个根。这些根有的是实数,有的是复数,如果是复数,必然以共轭形式成对出现,一对复根作为一组看待。一个实根可构成一个两点最小相位因式,一对复根构成一个三点最小相位因式。每个最小相位因式的样点前后倒序,就是最大相位因式。它们的根是互为倒数的。
对于一个最小相位因式(两点或三点)可以用常规方法得到该因式的反褶积算子,并用以对记录做反褶积。如果这个最小相位因式是子波的因式,那么反褶积以后的记录分辨率会有所提高。
设子波的因式是最大相位,求出最大相位因式的反褶积算子并对记录做反褶积。比较这两种反褶积结果,分辨率较高的结果就是所需要的中间结果。判断分辨率的高低可以采用例如记录各样点绝对值的5次方之和。
这样就确定了子波的一个因式,并得到了这个因式的反褶积结果。在这个因式的反褶积结果上用同样方法处理另一对因式。如此串联进行,一直到所有各对因式处理完毕,就完成了反褶积。
在此过程中,将每次选出的反褶积算子保存下来,就可以得到子波的反褶积算子。
本发明具体实现过程包括以下处理步骤:
1)采集单道地震数据并用常规的地震资料处理,确定地震资料信噪比较高的时间段为自相关时窗,对选择好的时窗数据求子波的自相关;
2)把求得的自相关各个样点值作为多项式的系数,然后对多项式作因式分解,求出自相关的各个子因式得到子波因子;
子波因式个数最好不要太多,以便影响因式分解的精度,最好把子波因式个数限制在二十个以内;
求得的自相关各个因式其实就是子波因子,其中既有最小相位子波因子,也有最大相位子波因子的对应体;
3)用通常的方法通过子波因子求得子波因子的反褶积算子;
这时反褶积算子都是最小相位的,如果是最小相位子波的因子,那么这个因子就是正确的,如果是最大相位子波对应的最小相位子波因子的对应体,那么就可以把这个最小相位子波因子反褶积算子进行倒序排列而得到最大相位子波的反褶积算子;
4)把得到的所有最小相位子波的反褶积算子倒序排列,得到每个子波因子最大,最小两个反褶积算子;
5)用最小相位反褶积算子对记录做反褶积;
如果这个最小相位因式是子波的因式,那么反褶积以后的记录分辨率会有所提高;
6)用最大相位反褶积算子对记录做反褶积;
7)对(5),(6)两个步骤的反褶积结果进行分选,确定分辨率较高的中间结果;
判断分辨率的高低可以采用例如记录各样点绝对值的5次方之和;
8)确定分辨率较高的中间结果的因式,得到这个因式的反褶积结果;
9)在前一步骤因式反褶积结果基础上,用步骤(3)至(8)试探处理步骤(2)求出的第二对子波因子,在第二对子波因子反褶积基础上进行第三对子波因子试探,如此串联进行对所有子波因子进行处理,就完成混合相位反褶积;
此时的地震数据的混合相位反褶积算子就是所试探出的所有因式反褶积算子的褶积结果;
10)采用通常的物探处理得到反褶积后提高地震数据分辨率的叠加剖面。
实施例1是本发明混合相位反褶积方法运用于理论数据的过程,图1是几种不同相位子波试验分析:图1(a)中是1 6个子波,它们具有相同的自相关,子波长度为7个采样点,图中最左边的是最小相位子波,最右边的是最大相位子波,中间14个是混合相位子波,最后一个为最大相位子波,且子波自相关一共有4个最小相位反褶积算子。图1(b),(d),(f),(h)分别是对图1(a)用四个最小相位反褶积算子串联褶积的结果;而图1(c),(e),(g),(i)是用四个算子中的最小相位反褶积算子与最大相位反褶机算子通过分辨率高低的判断串联褶积的四个步骤;不难看出最小相位反褶积只能把最小相位子波压缩成单位脉冲如图1(h),而对其它混合相位子波,最大相位子波却无能为力,本发明混合相位反褶积方法却能把最小相位,最大相位,混合相位子波都压缩成单位脉冲。如图1(i)。
实施例2是合成地震纪录试算:图2(a)是反射系数模型,所有地震道相同,中间空道。图2(b)是用图1(a)的子波与图2(a)的反射系数单道分别褶积合成的理论地震记录。图2(c),(d)是对图2(b)分别进行最小相位反褶积和本发明的混合相位反褶积结果。
由图2(c)可见,只有最左边一道最小相位合成记录得到了理想的反褶积结果,而非最小相位子波的合成记录,用最小相位反褶积不但没有取得好的反褶积效果,而且增加了高频噪音干扰。图2(d)是本发明的混合相位反褶积的结果,虽然图2(b)的合成记录各道子波不同,并且除了第一个地震道外,子波并不是最小相位的,但却通过大小相位因子反褶积后分辨率提高程度的判断,筛选出了适合混合相位子波反褶积的反褶积算子.从而使反褶积后地震数据分辨率达到了理想的结果如图2(d)。
实施例3是实际地震记录计算,从图3(b),(c)的两种反褶积,以及反褶积后叠加对比如图3(e),(f),可以看到,本发明混合相位反褶积方法的巨大优越性在于信噪比和分辨率都比最小相位反褶积高。最小相位反褶积因为不能适应非最小相位子波的实际情况,反褶积后产生了大量的噪音,掩盖了原有的有效信息,在很大程度上降低了信噪比。实际数据上应用表明,本发明的混合相位反褶积方法能显著提高地震数据分辨率,且在一定程度上能提高信噪比。
Claims (3)
1.一种提高地震数据处理分辨率的方法,其特征在于包括以下处理步骤:
1)采集单道地震数据并用常规的地震资料处理,确定地震资料信噪比较高的时间段为自相关时窗,对选择好的时窗数据求子波的自相关;
2)把求得的自相关各个样点值作为多项式的系数,然后对多项式作因式分解,求出自相关的各个子因式得到子波因子;
3)用通常的方法通过子波因子求得子波因子的反褶积算子;
4)把得到的所有最小相位子波的反褶积算子倒序排列,得到每个子波因子最大,最小两个反褶积算子;
5)用最小相位反褶积算子对记录做反褶积;
6)用最大相位反褶积算子对记录做反褶积;
7)对(5),(6)两个步骤的反褶积结果进行分选,确定分辨率较高的中间结果;
8)确定分辨率较高的中间结果的因式,得到这个因式的反褶积结果;
9)在前一步骤因式反褶积结果基础上,用步骤(3)至(8)处理步骤(2)的每一对因式,串联进行对因式处理,完成混合相位反褶积;
10)采用通常的物探处理得到反褶积后提高地震数据分辨率的叠加剖面。
2.根据权利要求1所述的提高地震数据处理分辨率的方法,其特征在于:子波因式个数在10-20个。
3.根据权利要求1所述的提高地震数据处理分辨率的方法,其特征在于:确定分辨率较高的中间结果采用录各样点绝对值的5次方之和。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2006101650182A CN101201405A (zh) | 2006-12-12 | 2006-12-12 | 一种提高地震数据处理分辨率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2006101650182A CN101201405A (zh) | 2006-12-12 | 2006-12-12 | 一种提高地震数据处理分辨率的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101201405A true CN101201405A (zh) | 2008-06-18 |
Family
ID=39516669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2006101650182A Pending CN101201405A (zh) | 2006-12-12 | 2006-12-12 | 一种提高地震数据处理分辨率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101201405A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813786A (zh) * | 2010-04-02 | 2010-08-25 | 中国石油集团西北地质研究所 | 子波处理二步法反褶积方法 |
CN102520442A (zh) * | 2011-12-02 | 2012-06-27 | 中国科学院地质与地球物理研究所 | 一种地震信号连续相干带的确定方法 |
CN103364827A (zh) * | 2012-03-30 | 2013-10-23 | 中国石油化工股份有限公司 | 基于双参数目标寻优的自适应谱模拟反褶积方法 |
CN106291679A (zh) * | 2015-05-27 | 2017-01-04 | 中国石油化工股份有限公司 | 一种针对山前带地震资料的反摺积方法 |
CN106526672A (zh) * | 2016-10-25 | 2017-03-22 | 北京天工新源科技有限公司 | 一种地震勘探特征子波整形提高分辨率方法 |
CN111736223A (zh) * | 2020-06-10 | 2020-10-02 | 中国石油天然气集团有限公司 | 地震数据的处理方法及装置 |
CN112379429A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 地震数据的振幅补偿方法及装置 |
CN112444864A (zh) * | 2019-08-28 | 2021-03-05 | 中国石油化工股份有限公司 | 薄层高分辨率反演方法及存储介质 |
CN112578436A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 一种子波提取方法及系统 |
CN112649847A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 一种提高地震垂向分辨率的调谐预测反褶积方法及系统 |
-
2006
- 2006-12-12 CN CNA2006101650182A patent/CN101201405A/zh active Pending
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813786A (zh) * | 2010-04-02 | 2010-08-25 | 中国石油集团西北地质研究所 | 子波处理二步法反褶积方法 |
CN102520442A (zh) * | 2011-12-02 | 2012-06-27 | 中国科学院地质与地球物理研究所 | 一种地震信号连续相干带的确定方法 |
CN102520442B (zh) * | 2011-12-02 | 2015-04-15 | 中国科学院地质与地球物理研究所 | 一种地震信号连续相干带的确定方法 |
CN103364827A (zh) * | 2012-03-30 | 2013-10-23 | 中国石油化工股份有限公司 | 基于双参数目标寻优的自适应谱模拟反褶积方法 |
CN106291679B (zh) * | 2015-05-27 | 2018-11-23 | 中国石油化工股份有限公司 | 一种针对山前带地震资料的反褶积方法 |
CN106291679A (zh) * | 2015-05-27 | 2017-01-04 | 中国石油化工股份有限公司 | 一种针对山前带地震资料的反摺积方法 |
CN106526672A (zh) * | 2016-10-25 | 2017-03-22 | 北京天工新源科技有限公司 | 一种地震勘探特征子波整形提高分辨率方法 |
CN112444864A (zh) * | 2019-08-28 | 2021-03-05 | 中国石油化工股份有限公司 | 薄层高分辨率反演方法及存储介质 |
CN112578436A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 一种子波提取方法及系统 |
CN112649847A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 一种提高地震垂向分辨率的调谐预测反褶积方法及系统 |
CN111736223A (zh) * | 2020-06-10 | 2020-10-02 | 中国石油天然气集团有限公司 | 地震数据的处理方法及装置 |
CN111736223B (zh) * | 2020-06-10 | 2023-12-22 | 中国石油天然气集团有限公司 | 地震数据的处理方法及装置 |
CN112379429A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 地震数据的振幅补偿方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101201405A (zh) | 一种提高地震数据处理分辨率的方法 | |
CN102109612B (zh) | 一种地震波吸收衰减补偿方法 | |
CN108267784A (zh) | 一种地震信号随机噪声压制处理方法 | |
CN102707314B (zh) | 一种多路径双谱域混合相位子波反褶积方法 | |
CN105445801B (zh) | 一种消除二维地震资料随机噪音的处理方法 | |
CN102736109B (zh) | 一种crp道集去噪、校正与叠加的方法 | |
CN102193107A (zh) | 一种地震波场分离与去噪方法 | |
CN109001813B (zh) | 一种压制多次波的方法、装置及系统 | |
CN109738951B (zh) | 一种基于地震同相轴子波谱的时变反褶积方法 | |
CN101354442B (zh) | 一种用于获取地层信息的混合相位反褶积方法及处理系统 | |
CN107144879A (zh) | 一种基于自适应滤波与小波变换结合的地震波降噪方法 | |
CN106646637A (zh) | 一种去除核磁信号中尖峰噪声的方法 | |
CN113935467A (zh) | 基于迭代多尺度注意力网络的das井中勘探数据噪声抑制方法 | |
CN100561256C (zh) | 一种避免动校拉伸的动校正叠加方法 | |
CN105044794A (zh) | 一种核磁共振回波数据的压缩方法及装置 | |
CN103630936A (zh) | 基于波束定向原理的地震单炮记录随机噪声压制方法 | |
CN113077386A (zh) | 基于字典学习和稀疏表征的地震资料高分辨率处理方法 | |
CN106950597B (zh) | 基于三边滤波的混合震源数据分离方法 | |
CN103675904A (zh) | 一种井震匹配解释性目标处理方法及装置 | |
CN105277984A (zh) | 一种时移地震互约束频率一致性处理方法 | |
CN105388519A (zh) | 一种提高地震资料分辨率的方法 | |
CN111505709A (zh) | 一种基于稀疏谱分解的衰减定性分析的方法 | |
CN113031070B (zh) | 一种深度域合成地震记录的制作方法 | |
CN104820244A (zh) | 一种提高石油勘探资料处理信噪比的方法 | |
CN105319583B (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 | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Open date: 20080618 |