CN106707334A - 一种提高地震资料分辨率的方法 - Google Patents
一种提高地震资料分辨率的方法 Download PDFInfo
- Publication number
- CN106707334A CN106707334A CN201611259346.9A CN201611259346A CN106707334A CN 106707334 A CN106707334 A CN 106707334A CN 201611259346 A CN201611259346 A CN 201611259346A CN 106707334 A CN106707334 A CN 106707334A
- Authority
- CN
- China
- Prior art keywords
- frequency
- seismic data
- original
- harmonic
- seismic
- 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 36
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 abstract description 5
- 239000000284 extract Substances 0.000 abstract description 3
- 238000012731 temporal analysis Methods 0.000 abstract 1
- 238000000700 time series analysis Methods 0.000 abstract 1
- 230000001131 transforming effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 12
- 238000000465 moulding Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000002123 temporal effect Effects 0.000 description 4
- 238000012952 Resampling Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000007630 basic procedure Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种提高地震资料分辨率的方法，具体步骤如下：对原始地震记录应用连续小波变换CWT时间序列分析；利用连续小波域的多分辨率特征，通过谐频分析计算出各尺度小波的谐波信息逆变换将谐波信息变换回时间域；将得到的时间阈信息加入到原始地震记录上，有效拓展地震数据的频宽。本发明从方法原理上突破了传统褶积模型的限制，不再需要估计地震子波，不再需要假设地层反射系数是稀疏的；利用优势频段信息提取谐波分量，保证拓频结果具有较高的信噪比，具体实现采用FFT算法，计算效率高，可以同时拓展低频与高频。
Description
技术领域
本发明涉及谐频技术领域，具体是一种提高地震资料分辨率的方法。
背景技术
在地震资料处理中，信噪比和分辨率是评价地震资料品质的两个重要指标。具有高信噪比和高分辨率的地震资料，可以用来进行有效的反演和地质解释，为实现精细油藏描述奠定坚实的基础。因此，如何有效的提高地震资料的分辨率和信噪比，是地震资料处理的一个关键问题。
我们结合地震信号在连续小波域的多分辨率特征，按照谐频分析的思路，扩展地震资料的有效带宽，能够极大的提高地震资料的分辨率，这是其它方法很难达到的。
一直以来，拓展地震信号的高频端和低频端是地震资料采集和处理所追求的目标。地震信号的频带范围一般在几赫兹到几十赫兹之间，其中高频成分主要影响地震资料的分辨率，低频成分主要影响深层构造的成像精度、叠前叠后反演的稳定性和准确性，甚至还可能直接影响油气检测。目前，针对高频成分的补偿方法研究较多，而有关低频保护与补偿的研究较少。
地震波在地下介质中传播时，高频信号相对于低频信号吸收衰减更快，造成了地震信号主频低、带宽窄的特点。将地震信号分解到连续小波域，则表示高频信号的小尺度小波系数能量相对较弱。基于宽频带高分辨率地震信号振幅谱能量近似水平的假设，在连续小波域补偿因为传播损失的高频能量，以恢复地震资料高分辨率特征。下面给出应用连续小波变换扩展单道地震资料频宽步骤：
(1)对单道地震信号f(x)做连续小波变换。将该地震道分解到时间—尺度域，W(τ，s)为拓频前地震信号的小波变换系数：
连续小波变换在时频平面内的良好局部化特性和冗余特性，提供了可以在频率域预测可扩展的频率，并且利用时间频率域中的不完全信息重建信号的可能。时间频率域的联合分辨率受不确定性原理的约束，根据尺度的变化而变化，尺度增大时，时间分辨率较低，频率分辨率较高；反之尺度减小时，时间分辨率较高，频率分辨率较低。
(2)在原始地震信号的振幅谱中选择一个基准频率，以该基准频率计算需要扩展的频率信息的扩展范围，即谐波和次谐波，见图1(以高频端的扩展为例)。
在图1的频率选择中，A点表示所选择的基准频率，B＝A/2；BA之间的频率为基础频率段。扩展频谱的低端与上述类似，基准频率需要重新选择，图1虚线即为扩展低端的基准频率。基准频率定义为一个用于扩展带宽的标准，在原始信号的振幅谱上选择；带宽扩展范围的计算采用倍频程的概念，扩展的高频分量的频率信息称为谐波(一次谐波、二次谐波等)，扩展的低频分量的频率信息称为次谐波(一次次谐波、二次次谐波等)；谐波为基准频率的整数倍，次谐波为基准频率的整数倍的倒数；基础频率为基准频率到基准频率之前一个倍频程之间的频率段，用来预测谐波和次谐波的振幅谱，并且对谐波和次谐波振幅谱进行能量密度调节，从而达到带宽扩展的目的。
(3)在时频域得到带有已扩展带宽的信息对其进行连续小波变换的逆变换，根据式(1)重建已扩频的时间信号
式中，CΨ满足可允许性条件。
上述应用连续小波变换扩展单道地震资料频宽的方法存在如下缺点
1、需要提取可靠的子波；
2、地层反射系数必须满足稀疏假设；
3、原始地震记录必须具有较高的信噪比，否则拓频结果会产生较多噪声；
4、求解压缩感知问题计算量相对较大，方法效率低。
发明内容
本发明的目的在于提供一种信噪比高、计算效率高的提高地震资料分辨率的方法，以解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：
一种提高地震资料分辨率的方法，具体步骤如下：
(1)对原始地震记录应用连续小波变换CWT时间序列分析：
(2)利用连续小波域的多分辨率特征，通过谐频分析计算出各尺度小波的谐波信息
(3)逆变换将谐波信息变换回时间域：
(4)将步骤(3)得到的时间阈信息加入到原始地震记录上，有效拓展地震数据的频宽。
作为本发明进一步的方案：所述步骤(1)中的原始地震记录包括原始时间阈地震记录和原始频率阈时间记录。
作为本发明进一步的方案：所述原始时间阈地震记录是地层反射系数与地震子波的褶积。
作为本发明再进一步的方案：所述原始频率阈时间记录是反射系数谱与地震子波谱的乘积。
与现有技术相比，本发明的有益效果是：
本发明从方法原理上突破了传统褶积模型的限制，不再需要估计地震子波，不再需要假设地层反射系数是稀疏的；利用优势频段信息提取谐波分量，保证拓频结果具有较高的信噪比，具体实现采用FFT算法，计算效率高，可以同时拓展低频与高频。
附图说明
图1为现有技术中由基准频率扩展出谐波、次谐波的示意图。
图2(a)为本发明中谐波拓频的基本原理之一。
图2(b)为本发明中谐波拓频的基本原理之二。
图3为本发明中谐波拓频的基本流程示意图。
图4为本发明中高频谐波合成示意图。
图5为本发明中极值褶积法重建原始信号的示意图。
图6为本发明中谐波信号计算示例示意图。
图7为本发明中谐波信号的频谱示意图。
图8为直接拓频后的谱形态异常问题示意图。
图9为本发明中权值预估算方法示意图。
图10为本发明中加入权值预估算的拓频结果示意图。
图11为横向能量异常示意图。
图12为本发明中强轴预压制技术流程示意图。
图13为本发明中强轴预压制技术效果对比示意图。
图14为本发明中数据重采样的示意图。
具体实施方式
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
请参阅图1-14，一种提高地震资料分辨率的方法，具体步骤如下：
(1)对原始地震记录应用连续小波变换CWT时间序列分析：
(2)利用连续小波域的多分辨率特征，通过谐频分析计算出各尺度小波的谐波信息
(3)逆变换将谐波信息变换回时间域：
(4)将步骤(3)得到的时间阈信息加入到原始地震记录上，有效拓展地震数据的频宽，进而更好地刻画地下的地层反射。
所述步骤(1)中的原始地震记录包括原始时间阈地震记录和原始频率阈时间记录；所述原始时间阈地震记录是地层反射系数与地震子波的褶积，所述原始频率阈时间记录是反射系数谱与地震子波谱的乘积。
谐波拓频的基本原理如图2所示，由图2(a)中可以看出，在原始信号中加入谐波信号，可以提高原始信号的分辨率。在图2(b)中可以看出，地震子波也可以认为是由不同频率的谐波构成，缺失高频谐波地震子波的分辨率偏低，增入高频谐波地震子波的分辨率得到增强。
图4给出了一个由原始数据得到高频谐波信号的形象过程，得到低频谐波信号的过程与此类似。
本发明采用极值褶积法能够快速得到谐波：谐波计算是谐波提频技术的一个重要环节，图5给出了利用极值褶积法重建原始信号过程，首先将原始信号进行小波分解，然后提取各尺度小波信号的极值(极大值和极小值)，利用这个极值信号与对应尺度的小波进行褶积，可以很好的重新合成小波信号，变换回时间域，恢复信号与原始信号保持一致。整个过程表明，“极值信号”很好的保留了原始信号的信息，当把“极值信号”与小尺度(1/2尺度)和大尺度(2倍尺度)小波信号进行褶积时，就会得到原始小波信号的谐波信号和次谐波信息。图6给出了一个谐波信号计算的例子，图7给出了谐波信号的频谱。
本发明通过权值预估算来控制谐波谱的形态：如图8所示，利用上述的谐波法直接进行拓频，拓频后谱形态会产生一些异常，出现一个“塌陷”频段，通过权值预估算对谐波信号进行能量密度调整，这样谐波和次谐波的能量谱形状会更好。图9给出了权值预估算技术的基本流程，利用期望输出谱逐步校正原始输入谱，采用非线性优化算法，在小波尺度域预先估算各尺度小波信号的权值，进而保证输出数据的频谱逼近期望输出谱。图10给出了加入权值预估算的拓频结果。
本发明通过“强轴”预压制解决“强轴”横向能量异常：当原始地震数据中有非常强的单界面反射信息时，在拓频结果中会出现“强轴”横向能量异常，如图11所示，为了压制这种“强轴”横向能量异常，研究了“强轴”预压制技术。图12给出了“强轴”预压制技术的技术流程，从原始数据中事先估计出一个强轴压制算子，这个算子对于强轴的压制程度是可以控制的，具体参数控制后面还有讨论，对原始数据应用强轴压制算子后再进行拓频，然后在反向应用强轴压制算子，进而得到最终的拓频结果。图13是“强轴”预压制技术效果对比图，应用“强轴”预压制技术后，“强轴”横向能量异常得到了削弱。
本发明通过程序实现充分利用FFT算法来提高计算效率：整个谐波拓频技术涉及了数据重采样(如图14，数据重采样是保证算法横向连续性必要条件)、小波变换、谐波计算等算法，为了保证整个谐波拓频算法的高效性，编程实现时，充分利用了FFT算法，这也是该拓频算法走向叠前的一个必要条件。
本发明从方法原理上突破了传统褶积模型的限制，不再需要估计地震子波，不再需要假设地层反射系数是稀疏的；利用优势频段信息提取谐波分量，保证拓频结果具有较高的信噪比，具体实现采用FFT算法，计算效率高，可以同时拓展低频与高频。
上面对本专利的较佳实施方式作了详细说明，但是本专利并不限于上述实施方式，在本领域的普通技术人员所具备的知识范围内，还可以在不脱离本专利宗旨的前提下作出各种变化。
Claims (4)
1.一种提高地震资料分辨率的方法，其特征在于，具体步骤如下：
(1)对原始地震记录应用连续小波变换CWT时间序列分析：
(2)利用连续小波域的多分辨率特征，通过谐频分析计算出各尺度小波的谐波信息
(3)逆变换将谐波信息变换回时间域：
(4)将步骤(3)得到的时间阈信息加入到原始地震记录上，有效拓展地震数据的频宽。
2.根据权利要求要求1所述的提高地震资料分辨率的方法，其特征在于，所述步骤(1)中的原始地震记录包括原始时间阈地震记录和原始频率阈时间记录。
3.根据权利要求要求1或2所述的提高地震资料分辨率的方法，其特征在于，所述原始时间阈地震记录是地层反射系数与地震子波的褶积。
4.根据权利要求要求1或2所述的提高地震资料分辨率的方法，其特征在于，所述原始频率阈时间记录是反射系数谱与地震子波谱的乘积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611259346.9A CN106707334B (zh) | 2016-12-30 | 2016-12-30 | 一种提高地震资料分辨率的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611259346.9A CN106707334B (zh) | 2016-12-30 | 2016-12-30 | 一种提高地震资料分辨率的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106707334A true CN106707334A (zh) | 2017-05-24 |
| CN106707334B CN106707334B (zh) | 2018-08-10 |
Family
ID=58906251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611259346.9A Active CN106707334B (zh) | 2016-12-30 | 2016-12-30 | 一种提高地震资料分辨率的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106707334B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107436450A (zh) * | 2017-07-26 | 2017-12-05 | 西安交通大学 | 一种基于连续小波变换的地震信号带宽拓展方法 |
| CN108508489A (zh) * | 2018-03-19 | 2018-09-07 | 中国海洋石油集团有限公司 | 一种基于波形微变化匹配的地震反演方法 |
| CN109239772A (zh) * | 2018-08-31 | 2019-01-18 | 中海石油(中国)有限公司湛江分公司 | 一种特征波域的浅水多次波模型指向性预测方法 |
| CN112558158A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | 一种基于测井曲线的地震数据双向拓频方法及系统 |
| CN112684500A (zh) * | 2021-01-05 | 2021-04-20 | 南方海洋科学与工程广东省实验室(广州) | 一种获得具有丰富低频信息目标子波的方法及处理终端 |
| CN113281808A (zh) * | 2021-04-22 | 2021-08-20 | 南方海洋科学与工程广东省实验室(湛江) | 一种抗频散地震波正演方法、系统、装置及介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4490788A (en) * | 1982-09-29 | 1984-12-25 | Schlumberger Technology Corporation | Well-logging data processing system having segmented serial processor-to-peripheral data links |
| CN104914466A (zh) * | 2015-06-26 | 2015-09-16 | 中国石油大学(华东) | 一种提高地震资料分辨率的方法 |
| CN105487115A (zh) * | 2014-09-17 | 2016-04-13 | 中国石油化工股份有限公司 | 一种基于小波变换的高频延拓方法 |
-
2016
- 2016-12-30 CN CN201611259346.9A patent/CN106707334B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4490788A (en) * | 1982-09-29 | 1984-12-25 | Schlumberger Technology Corporation | Well-logging data processing system having segmented serial processor-to-peripheral data links |
| CN105487115A (zh) * | 2014-09-17 | 2016-04-13 | 中国石油化工股份有限公司 | 一种基于小波变换的高频延拓方法 |
| CN104914466A (zh) * | 2015-06-26 | 2015-09-16 | 中国石油大学(华东) | 一种提高地震资料分辨率的方法 |
Non-Patent Citations (2)
| Title |
|---|
| MICHAEL SMITH .ETL: "Extending seismic bandwidth using the continuous wavelet transform", 《LEVERAGING TECHNOLOGY》 * |
| 尹继尧 等: "基于傅里叶算法连续小波变换分析技术在提高准噶尔盆地西北缘地震资料分辨率中的应用", 《地球物理学进展》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107436450A (zh) * | 2017-07-26 | 2017-12-05 | 西安交通大学 | 一种基于连续小波变换的地震信号带宽拓展方法 |
| CN108508489A (zh) * | 2018-03-19 | 2018-09-07 | 中国海洋石油集团有限公司 | 一种基于波形微变化匹配的地震反演方法 |
| CN109239772A (zh) * | 2018-08-31 | 2019-01-18 | 中海石油(中国)有限公司湛江分公司 | 一种特征波域的浅水多次波模型指向性预测方法 |
| CN112558158A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | 一种基于测井曲线的地震数据双向拓频方法及系统 |
| CN112558158B (zh) * | 2019-09-25 | 2024-04-09 | 中国石油化工股份有限公司 | 一种基于测井曲线的地震数据双向拓频方法及系统 |
| CN112684500A (zh) * | 2021-01-05 | 2021-04-20 | 南方海洋科学与工程广东省实验室(广州) | 一种获得具有丰富低频信息目标子波的方法及处理终端 |
| CN113281808A (zh) * | 2021-04-22 | 2021-08-20 | 南方海洋科学与工程广东省实验室(湛江) | 一种抗频散地震波正演方法、系统、装置及介质 |
| CN113281808B (zh) * | 2021-04-22 | 2023-10-20 | 南方海洋科学与工程广东省实验室(湛江) | 一种抗频散地震波正演方法、系统、装置及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106707334B (zh) | 2018-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106707334B (zh) | 一种提高地震资料分辨率的方法 | |
| CN107589454A (zh) | 一种基于vmd‑tfpf压制地震勘探随机噪声方法 | |
| CN105093294B (zh) | 基于可变模态分解的地震波衰减梯度估计方法 | |
| CN102053273A (zh) | 一种对地震波信号进行反q滤波的方法 | |
| CN105388518A (zh) | 一种质心频率与频谱比联合的井中地震品质因子反演方法 | |
| CN101923176B (zh) | 一种利用地震数据瞬时频率属性进行油气检测的方法 | |
| CN103163554A (zh) | 利用零偏vsp资料估计速度和q值的自适应波形的反演方法 | |
| CN107132579A (zh) | 一种保地层结构的地震波衰减补偿方法 | |
| CN107144879A (zh) | 一种基于自适应滤波与小波变换结合的地震波降噪方法 | |
| CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
| CN106680874A (zh) | 基于波形形态特征稀疏化建模的谐波噪声压制方法 | |
| CN104199093A (zh) | 基于时频域能量自适应加权的地震信号分辨率增强方法 | |
| CN104280765A (zh) | 基于变子波反射系数反演的地震高分辨处理方法 | |
| CN107102356A (zh) | 基于ceemd的地震信号高分辨率处理方法 | |
| CN106405645A (zh) | 一种基于资料品质分析的信噪比可控的地震拓频处理方法 | |
| CN107272063A (zh) | 基于高分辨率时频分析和一致性度量的非均质性刻画方法 | |
| CN106772586A (zh) | 一种基于地震信号奇异性的隐蔽性断裂检测方法 | |
| Zoukaneri et al. | A combined Wigner-Ville and maximum entropy method for high-resolution time-frequency analysis of seismic data | |
| CN104730576A (zh) | 基于Curvelet变换的地震信号去噪方法 | |
| CN107678064A (zh) | 一种声波时差实时提取方法 | |
| CN102692647A (zh) | 一种高时间分辨率的地层含油气性预测方法 | |
| CN104932018A (zh) | 补偿变分辨率因子s变换的复时-频谱提高地震剖面分辨率的方法 | |
| CN107436450A (zh) | 一种基于连续小波变换的地震信号带宽拓展方法 | |
| CN102338886B (zh) | 一种有效衰减三分量地震记录中面波的极化滤波方法 | |
| CN106353795A (zh) | 一种高精度的地层衰减参数计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |