CN105929447B - 考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法 - Google Patents
考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法 Download PDFInfo
- Publication number
- CN105929447B CN105929447B CN201610377985.9A CN201610377985A CN105929447B CN 105929447 B CN105929447 B CN 105929447B CN 201610377985 A CN201610377985 A CN 201610377985A CN 105929447 B CN105929447 B CN 105929447B
- Authority
- CN
- China
- Prior art keywords
- radon
- seismic
- mrow
- hyperbolic
- msubsup
- 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.)
- Active
Links
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 79
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000000694 effects Effects 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 17
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 10
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 12
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 238000011426 transformation method Methods 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 4
- 239000000523 sample Substances 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000004454 trace mineral analysis Methods 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 abstract 1
- 238000000844 transformation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 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/282—Application of seismic models, synthetic seismograms
-
- 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
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
本发明公开了一种考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法。该方法利用地震同相轴的倾角及其导数来计算双曲线同相轴的顶点位置和速度，并利用聚类算法来确定最终的双曲线同相轴个数，降低所估计的同相轴参数个数；利用所计算的同相轴顶点位置和速度，采用开窗分割地震道的方法来构建变顶点双曲线Radon变换算子，消除Radon反变换产生的地震子波拉伸效应；然后在优化问题中引入Radon域模型和数据拟合项的L1范数最小化约束，来提高Radon域模型的稀疏性和对非高斯拟合误差的适应性；另外，本发明方法利用乘数交替方向方法求解L1‑L1范数最小化优化问题，在整个迭代过程中只需一次矩阵求逆，计算复杂度较低。
Description
技术领域
本发明属于地震勘探技术中的地震信号处理领域，具体涉及一种考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法。
背景技术
Radon变换广泛应用于地震信号处理中，如多次波压制、地震数据规则化和波场分离等，其主要思想是将地震数据从时间-偏移距域映射到Radon域。Radon变换按照积分轨迹区分为三大类：线性Radon变换、抛物线Radon变换和双曲线Radon变换。
其中，双曲线Radon变换可直接应用于未做动校正的共中心点道集或共炮点道集，并利用双曲线轨迹来构造变换算子，能够更好地表示地震数据。
由于双曲线Radon变换的时变特性，计算双曲线Radon变换需要对大型稀疏双曲线Radon变换算子进行矩阵求逆，导致时间域双曲线Radon变换过程中的巨大计算量。另外，传统的双曲线Radon变换将顶点位置限制在零偏移距处，利用不同的速度参数来表征地震道集中的双曲线轨迹。当双曲线同相轴的顶点位于非零偏移距位置，如道集中存在微曲多次波、绕射多次波，传统的双曲线Radon变换进行多次波压制或数据规则化时，其效果会变差。
双曲线Radon变换的应用效果依赖于变换算子是否能有效表征道集中存在的双曲线同相轴。对于双曲线同相轴的顶点沿偏移距移动的情况，变顶点双曲线Radon变换方法通过将Radon域模型的维度空间从2D(顶点时间和速度)扩展为3D(顶点时间、速度和顶点偏移距)，能够更有效地表示顶点位于非零偏移距的双曲线同相轴。由于Radon域模型的空间维度增加一维，导致在Radon域模型中需要求解的未知参数个数变多，增加了Radon正反变换的计算量。另外，变顶点双曲线Radon变换方法对应的Radon变换方程组为欠定方程组，需要在相应的优化问题中引入Radon域模型的稀疏性约束来计算Radon域模型。传统的变顶点稀疏双曲线Radon变换方法采用迭代重加权最小二乘算法求解L1范数最小化优化问题，在每一步迭代中需要计算一次矩阵求逆，计算复杂度较高。因此，传统变顶点稀疏双曲线Radon变换的计算量较大，这限制了该变换方法的实用性。另外，传统的双曲线Radon变换在构建双曲线Radon变换算子进行正反变换时，并没有考虑地震子波的存在，在利用Radon反变换进行地震数据重构过程中，会出现地震子波拉伸效应，造成地震信号失真。
发明内容
针对现有技术中存在的上述技术问题，本发明提出了一种考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法，其采用如下技术方案：
考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法，包括如下步骤：
a、设置变量初始值：需要设置初始值的变量包括阻尼因子β、乘数交替方向方法的阈值χ和φ、乘数交替方向方法的最大迭代次数N；
b、从所有地震道集中选取某一个地震道集，利用复地震道分析方法估计地震同相轴的倾角pr，并计算倾角导数prr；
c、利用下式计算双曲线同相轴的顶点位置(τ,x0)和速度v：
x0＝xr-tprv2，其中，τ表示顶点位置中的时间位置，x0表示顶点位置中的空间位置，t表示时间，xr表示检波点的空间位置；
d、在利用步骤c得到双曲线同相轴的顶点位置和速度后，利用聚类算法来检测双曲线同相轴的个数，确定最终的同相轴顶点位置和速度；
e、利用步骤d得到的同相轴顶点位置和速度，确定双曲线同相轴的轨迹，然后，利用开窗分割地震道的方法沿着垂直该轨迹的方向选取数据窗口，该数据窗口的中心与双曲线同相轴的轨迹重合，并且在垂直双曲线同相轴轨迹的方向，数据窗口的长度选为地震子波长度，最后，采用该数据窗口所对应的多条双曲线轨迹来构建变顶点双曲线Radon变换算子L；
f、利用变顶点双曲线Radon变换算子L得到如下的数学模型：d＝Lm；
其中，d表示时间域的原始地震数据，m表示时间域的Radon域模型；
对Radon域模型施加L1范数最小化约束，同时，对拟合误差施加L1范数最小化约束，对应的L1-L1范数最小化优化问题为：其中λ表示正则化因子，用来均衡对Radon域模型的L1范数最小化约束和对数据拟合项的L1范数最小化约束；
g、计算逆矩阵其中LH表示矩阵L的转置矩阵；并计算矩阵利用乘数交替方向方法来求解上述的L1-L1范数最小化优化问题，得到Radon域模型的估计结果；
h、判断所有地震道集是否处理完毕；如果否，返回执行上述步骤b至步骤g；如果全部处理完毕，则输出处理结果。
所述步骤g中利用乘数交替方向方法求解L1-L1范数最小化优化问题的计算过程如下：
g1、在每一步迭代中，Radon域模型的计算结果为：
其中，和为中间变量，并且
g2、判断迭代次数n是否达到最大迭代次数N；如果没达到，继续执行步骤g3；如果达到，输出Radon域模型的估计结果m(n+1)；
g3、利用距离算子计算中间变量
其中，定义距离算子为：定义距离算子
v＝{vi,j}，C＝max(|vi,j|)，
g4、对中间变量和计算如下：
并返回到步骤g1。
本发明具有如下优点：
本发明方法利用地震同相轴的倾角及其导数来计算双曲线同相轴的顶点位置和速度，并利用聚类算法来确定最终的双曲线同相轴个数，降低所估计的同相轴参数个数；利用所计算的同相轴顶点位置和速度，采用开窗分割地震道的方法来构建变顶点双曲线Radon变换算子，消除Radon反变换产生的地震子波拉伸效应；然后在优化问题中引入模型空间(Radon域模型)和数据拟合项(Radon反变换残差)的L1范数最小化约束，来提高Radon域模型的稀疏性和对非高斯拟合误差的适应性；另外，本发明方法利用乘数交替方向方法(Alternating Direction Method of Multipliers)求解优化问题，乘数交替方向方法利用距离算子求解L1范数最小化问题，该方法在整个迭代过程中只需一次矩阵求逆，计算复杂度较低。
附图说明
图1为本发明中考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法的流程图；
图2为本发明中缺失近偏移距道的3个满足双曲线时距曲线关系的同相轴示意图；
图3为利用本发明方法估计的倾角示意图；
图4为利用本发明方法所估计倾角的导数示意图；
图5为本发明中聚类方法检测到的双曲线同相轴的个数示意图；
图6为本发明中Radon正变换后的结果示意图；
图7为本发明中Radon反变换后的结果示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法，包括如下步骤：
a、设置变量初始值：需要设置初始值的变量包括阻尼因子β、乘数交替方向方法的阈值χ和φ、乘数交替方向方法的最大迭代次数N；
b、从所有地震道集中选取某一个地震道集，利用复地震道分析方法估计地震同相轴的倾角pr，并计算倾角导数prr；
c、利用下式计算双曲线同相轴的顶点位置(τ,x0)和速度v：
x0＝xr-tprv2，其中，τ表示顶点位置中的时间位置，x0表示顶点位置中的空间位置，t表示时间，xr表示检波点的空间位置；
d、在利用步骤c得到双曲线同相轴的顶点位置(τ,x0)和速度v后，利用聚类算法来检测双曲线同相轴的个数，确定最终的同相轴顶点位置和速度；
e、利用步骤d得到的同相轴顶点位置和速度，确定双曲线同相轴的轨迹，然后，利用开窗分割地震道的方法沿着垂直该轨迹的方向选取数据窗口，该数据窗口的中心与双曲线同相轴的轨迹重合，并且在垂直双曲线同相轴轨迹的方向，数据窗口的长度选为地震子波长度，最后，采用该数据窗口所对应的多条双曲线轨迹来构建变顶点双曲线Radon变换算子L，来消除Radon反变换重构地震数据时出现的地震子波拉伸效应；
f、利用变顶点双曲线Radon变换算子L得到如下的数学模型：d＝Lm；
其中，d表示时间域的原始地震数据，m表示时间域的Radon域模型；
为提高Radon域模型的稀疏性，对Radon域模型施加L1范数最小化约束。同时，为有效提高本发明方法对非高斯拟合误差的适应性，对拟合误差施加L1范数最小化约束。对应的L1-L1范数最小化优化问题为：其中，λ表示正则化因子，用来均衡对Radon域模型的L1范数最小化约束和对数据拟合项的L1范数最小化约束；
g、为了求解上式中给定的优化问题，首先计算逆矩阵其中，阻尼因子β用来保证矩阵求逆的稳定性，LH表示矩阵L的转置矩阵，并计算矩阵然后利用乘数交替方向方法来求解上述的L1-L1范数最小化优化问题，得到Radon域模型的估计结果；
该步骤g中利用乘数交替方向方法求解L1-L1范数最小化优化问题的计算过程如下：
g1、在每一步迭代中，Radon域模型的计算结果为：
其中，和为中间变量，并且
g2、判断迭代次数n是否达到最大迭代次数N；如果没达到，继续执行步骤g3；如果达到，输出Radon域模型的估计结果m(n+1)；
g3、利用距离算子计算中间变量和
其中，定义距离算子为：定义距离算子
v＝{vi,j}，C＝max(|vi,j|)，
g4、对中间变量和计算如下：
并返回到步骤g1；
h、判断所有地震道集是否处理完毕；如果否，返回执行上述步骤b至步骤g；如果全部处理完毕，则输出处理结果。
在仿真实验中，利用模型数据验证本发明方法的有效性：
图2为3个满足双曲线时距曲线关系的同相轴，并且缺失近偏移距，横坐标TraceNumber表示道号，单位为米(m)，纵坐标Time Sample Number表示时间采样点个数，采样间距为2毫秒(ms)，道间距为4.5米(m)。在该例子中，利用本发明方法来对缺失的近偏移距进行外推，来验证本发明中Radon变换的有效性。图3为利用复地震道分析方法估计的局部倾角，图4为对应的倾角导数。利用倾角及其倾角导数进行同相轴顶点位置和速度估计，并利用聚类算法检测同相轴的个数，图5为所检测到的3个同相轴顶点位置。图6为利用所构建的Radon变换算子进行Radon正变换所得到的结果。图7为利用图6中的Radon变换结果进行反变换后得到的同相轴，可以看到近偏移距得到有效的重建。注意到本发明实例中的双曲线同相轴的顶点并非完全位于零偏移距处，本发明方法能有效重建顶点不在零偏移距处的同相轴。
当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。
Claims (2)
1.考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法，其特征在于，包括如下步骤：
a、设置变量初始值：需要设置初始值的变量包括阻尼因子β、乘数交替方向方法的阈值χ和φ、乘数交替方向方法的最大迭代次数N；
b、从所有地震道集中选取某一个地震道集，利用复地震道分析方法估计地震同相轴的倾角pr，并计算倾角导数prr；
c、利用下式计算双曲线同相轴的顶点位置(τ,x0)和速度v：
x0＝xr-tprv2，其中，τ表示顶点位置中的时间位置，x0表示顶点位置中的空间位置，t表示时间，xr表示检波点的空间位置；
d、在利用步骤c得到双曲线同相轴的顶点位置和速度后，利用聚类算法来检测双曲线同相轴的个数，确定最终的同相轴顶点位置和速度；
e、利用步骤d得到的同相轴顶点位置和速度，确定双曲线同相轴的轨迹，然后，利用开窗分割地震道的方法沿着垂直该轨迹的方向选取数据窗口，该数据窗口的中心与双曲线同相轴的轨迹重合，并且在垂直双曲线同相轴轨迹的方向，数据窗口的长度选为地震子波长度，最后，采用该数据窗口所对应的多条双曲线轨迹来构建变顶点双曲线Radon变换算子L；
f、利用变顶点双曲线Radon变换算子L得到如下的数学模型：d＝Lm；
其中，d表示时间域的原始地震数据，m表示时间域的Radon域模型；
对Radon域模型施加L1范数最小化约束，同时，对拟合误差施加L1范数最小化约束，对应的L1-L1范数最小化优化问题为：其中λ表示正则化因子，用来均衡对Radon域模型的L1范数最小化约束和对数据拟合项的L1范数最小化约束；
g、计算逆矩阵其中LH表示矩阵L的转置矩阵；并计算矩阵利用乘数交替方向方法来求解上述的L1-L1范数最小化优化问题，得到Radon域模型的估计结果；
h、判断所有地震道集是否处理完毕；如果否，返回执行上述步骤b至步骤g；如果全部处理完毕，则输出处理结果。
2.根据权利要求1所述的考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法，其特征在于，
所述步骤g中利用乘数交替方向方法求解L1-L1范数最小化优化问题的计算过程如下：
g1、在每一步迭代中，Radon域模型的计算结果为：其中，和为中间变量，并且
g2、判断迭代次数n是否达到最大迭代次数N；如果没达到，继续执行步骤g3；如果达到，输出Radon域模型的估计结果m(n+1)；
g3、利用距离算子计算中间变量和
<mrow>
<msubsup>
<mi>z</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>prox</mi>
<mi>&chi;</mi>
<mn>1</mn>
</msubsup>
<mo>&lsqb;</mo>
<msup>
<mi>Lm</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<mi>d</mi>
<mo>+</mo>
<msubsup>
<mi>b</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>&rsqb;</mo>
<mo>,</mo>
<msubsup>
<mi>z</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>prox</mi>
<mi>&phi;</mi>
<mn>1</mn>
</msubsup>
<mo>&lsqb;</mo>
<msup>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msubsup>
<mi>b</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
1
其中，定义距离算子为：定义距离算子
v＝{vi,j}，C＝max(|vi,j|)，
g4、对中间变量和计算如下：
并返回到步骤g1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610377985.9A CN105929447B (zh) | 2016-05-31 | 2016-05-31 | 考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610377985.9A CN105929447B (zh) | 2016-05-31 | 2016-05-31 | 考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105929447A CN105929447A (zh) | 2016-09-07 |
| CN105929447B true CN105929447B (zh) | 2017-11-21 |
Family
ID=56832828
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610377985.9A Active CN105929447B (zh) | 2016-05-31 | 2016-05-31 | 考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105929447B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107783191B (zh) * | 2018-01-26 | 2019-02-15 | 吉林大学 | 多维空间时空时频峰值滤波消减地震勘探随机噪声的方法 |
| CN113687416B (zh) * | 2021-08-19 | 2022-08-30 | 吉林大学 | 一种针对倾斜金属矿体成像方法 |
| CN116736378B (zh) * | 2023-05-05 | 2024-02-02 | 广东海洋大学 | 一种基于Lq范数约束的用于地震勘探的速度分析方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102288994A (zh) * | 2011-07-01 | 2011-12-21 | 中国海洋石油总公司 | Radon谱约束下高维地震数据规则化方法 |
| CN102879824A (zh) * | 2012-09-07 | 2013-01-16 | 清华大学 | 一种基于迭代收缩的快速稀疏Radon变换方法 |
| CN103926622A (zh) * | 2014-05-06 | 2014-07-16 | 王维红 | 一种基于l1范数多道匹配滤波压制多次波的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6889142B2 (en) * | 2003-05-30 | 2005-05-03 | Pgs Exploration (Uk) Limited | Method of correcting for time shifts in seismic data resulting from azimuthal variation |
| US7733741B2 (en) * | 2007-04-05 | 2010-06-08 | Pgs Geophysical As | Method for reducing 3-D migration operator to 2-D migration operator for inhomogeneous media |
-
2016
- 2016-05-31 CN CN201610377985.9A patent/CN105929447B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102288994A (zh) * | 2011-07-01 | 2011-12-21 | 中国海洋石油总公司 | Radon谱约束下高维地震数据规则化方法 |
| CN102879824A (zh) * | 2012-09-07 | 2013-01-16 | 清华大学 | 一种基于迭代收缩的快速稀疏Radon变换方法 |
| CN103926622A (zh) * | 2014-05-06 | 2014-07-16 | 王维红 | 一种基于l1范数多道匹配滤波压制多次波的方法 |
Non-Patent Citations (4)
| Title |
|---|
| "Latest views of the sparse Radon transform";Daniel Trad 等;《GEOPHYSICS》;20031231;第68卷(第1期);第386–399页 * |
| "基于稀疏反演算法的高分辨率Radon变换及其在多次波压制中的应用";范景文 等;《物探与化探》;20151231;第39卷(第6期);第1245-1250页 * |
| "基于高阶稀疏Radon变换的预测多次波自适应相减方法";薛亚茹 等;《中国石油大学学报（自然科学版）》;20150228;第39卷(第1期);第43-49页 * |
| "高分辨率Radon变换方法及其在地震信号处理中的应用";刘喜武 等;《地球物理学进展》;20040131;第19卷(第1期);第008-015页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105929447A (zh) | 2016-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102854533B (zh) | 一种基于波场分离原理提高地震资料信噪比的去噪方法 | |
| CN102831588B (zh) | 一种三维地震图像的降噪处理方法 | |
| CN105929447B (zh) | 考虑地震子波拉伸效应的变顶点稀疏双曲线Radon变换方法 | |
| CN106483559B (zh) | 一种地下速度模型的构建方法 | |
| CN105277985A (zh) | 一种基于图像处理的ovt域地震数据规则化方法 | |
| CN110471113A (zh) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 | |
| CN111666721B (zh) | 全波形反演方法、装置及电子设备 | |
| CN109375265A (zh) | 一种基于变相位雷克子波匹配追踪的理想地震谱分解方法 | |
| CN110046323B (zh) | 同步压缩变换与重构的快速计算方法 | |
| CN103412328A (zh) | 基于交错网格有限差分算法的波数域保幅波场分离方法 | |
| CN105259575B (zh) | 快速3d自由表面多次波预测方法 | |
| CN106772585A (zh) | 一种基于弹性波解耦方程的优化拟解析方法及装置 | |
| CN104749625B (zh) | 一种基于正则化技术的地震数据倾角估计方法及装置 | |
| CN105005073A (zh) | 基于局部相似度和评价反馈的时变子波提取方法 | |
| CN102162858A (zh) | 利用非对称走时进行动校正速度分析的方法 | |
| CN110244353A (zh) | 一种基于稀疏范数优化算法的地震数据规则化方法 | |
| CN111694056B (zh) | 一种压制地震资料异常噪声的方法、存储介质及设备 | |
| CN106842297A (zh) | 井约束非稳态相位校正方法 | |
| CN113156509A (zh) | 基于饱和介质精确Zeoppritz方程的地震振幅反演方法及系统 | |
| CN116068644B (zh) | 一种利用生成对抗网络提升地震数据分辨率和降噪的方法 | |
| CN105259574B (zh) | 基于一次波稀疏约束的多道预测反褶积方法 | |
| Trad | Five dimensional seismic data interpolation | |
| CN109856673B (zh) | 一种基于优势频率迭代加权的高分辨Radon变换数据分离技术 | |
| CN112630840A (zh) | 基于统计特征参数的随机反演方法及处理器 | |
| CN110687605A (zh) | 基于改进的k-svd算法在地震信号处理的算法分析应用 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |