CN107589448A - 一种多道地震记录反射系数序列同时反演方法 - Google Patents
一种多道地震记录反射系数序列同时反演方法 Download PDFInfo
- Publication number
- CN107589448A CN107589448A CN201710572136.3A CN201710572136A CN107589448A CN 107589448 A CN107589448 A CN 107589448A CN 201710572136 A CN201710572136 A CN 201710572136A CN 107589448 A CN107589448 A CN 107589448A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msup
- reflection coefficient
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种多道地震记录反射系数序列同时反演方法。该方法充分考虑地层反射系数序列的稀疏特征以及地震道的空间横向连续性，利用参数化指数函数来度量反射系数序列的稀疏性，建立相应的约束优化目标函数。首先，基于褶积模型的特点和反射系数的稀疏性假设，在不丢失任何信息的前提下，设计降维预处理矩阵，降低约束目标函数的维度；然后，利用信赖域序列二次规划算法对优化问题进行高效求解。该算法具有全局收敛性和超线性收敛率，计算效率高，且具有较高反演精度。相比于传统的单道地震记录反褶积算法，本发明算法由于考虑了地震道的空间相关性，反演结果更稳定、可靠。
Description
技术领域
本发明属于地球物理勘探技术领域，涉及一种地震记录反射系数序列反演方法，特别涉及一种用于多道地震记录反射系数序列同时反演方法。
背景技术
在地层－岩性油气藏中，薄互层型储层是主要储层类型之一。由于激发地震所用人工震源频带范围有限，薄层产生的反射波走时差别小，从而使得反射波相互叠加干涉，接收到的地震记录无法直接分辨出薄互层结构，为薄互层型储层识别带来诸多困难。因此，有必要提高反射地震记录的分辨率，实现对薄互层地层结构的高精度刻画。
反褶积类方法通过压缩地震子波，进而获取地层反射系数序列，是地震资料处理中提高纵向分辨率、精细刻画地层结构常用方法之一。反褶积的理论基础一般是建立在满足若干基本假设的褶积模型之上。根据对地震资料的不同假设，地震记录反褶积方法种类较多，比如维纳反褶积、最小熵反褶积、最小平方反褶积、预测反褶积、互信息率反褶积、统计反褶积和稀疏脉冲反褶积等。
稀疏脉冲反褶积是建立在地层反射系数满足稀疏性假设的条件之上的，即大部分反射系数幅值趋近于零，剩余的大幅值反射系数对应着主要的岩性变化界面。基于对大量的实际测井数据分析，人们认识到在一定的地质沉积背景下，地层反射系数可用稀疏脉冲序列进行表示。统计上，可用 Bernoulli-Gaussian过程来描述服从稀疏分布的地层反射系数序列。在此基础之上，已经发展了一系列稀疏脉冲反褶积方法。Taylor等于1979年率先探讨了利用lp范数来实现地震稀疏脉冲反褶积方法。Oldenbug等利用l1范数反演层状地质结构的反射系数序列，取得较好的效果。Kormylo和Mendel在稀疏脉冲假设条件下，提出最大似然反褶积技术。Lavielle提出了基于 Bernoulli-Gaussian过程的Bayesian反褶积方法。Sachhi提出利用重加权策略实现地震反褶积。Kaaresen基于后验概率极大准则，提出迭代窗口校正函数极大化的稀疏脉冲反褶积方法。Velis提出统计稀疏脉冲反褶积方法，利用模拟退火法来确定反射层所在的位置，然后利用线性最小二乘方法估计反射系数的振幅。Gholami和Sacchi将基追踪去噪模型引入地震反褶积，提出一种从含噪地震记录中估计稀疏脉冲反射系数序列的方法。
以上现有技术具有如下缺点：
(1)常规的反射系数序列反演算法通常是基于梯度或共轭梯度构造下降方向，其收敛率是线性的，收敛速度较慢；
(2)常规的反射系数序列反演通常是逐道反演的，忽略了地震记录的横向连续性，进而使得后续叠加阻抗低频噪声严重。
发明内容
本发明的目的在于克服上述现有技术的缺点，提供一种新的用于多道地震记录反射系数序列同时反演方法。该方法基于地震记录褶积模型，充分考虑临近地震道之间的相关性，并利用近似零范数度量反射系数序列的稀疏特征，建立一种新的多道地震记录反射系数同时反演模型。在数值实现上，基于褶积模型特点和反射系数稀疏性特征，利用预处理技术降低求解维度，然后采用信赖域二次规划算法对本发明建立的反演模型进行求解，具有较高的计算效率和全局收敛性。该技术方案易于实现，可操作性强。
本发明的目的是通过以下技术方案来解决的：
该种用于多道地震记录反射系数序列同时反演方法，包括以下步骤：
1)采集原始地震资料，然后对采集到的地震资料进行预处理，包括静校正处理、噪音压制处理、真振幅恢复处理，得到叠后地震数据，称这个地震数据为观测地震数据，记为S(x,t)，其中x表示CDP号，t表示时间轴；
2)基于统计性方法提取地震子波w(t)；
3)利用快速衰减的指数函数fσ(·)来逐渐逼近零范数，并给定初始值σ；
4)确定正则化参数μ，估计地震记录噪声能量ε2，记反射系数剖面为R，记优化目标函数为h(R,σ,μ)，记约束条件为c(R)，建立多道地震记录反射系数序列同时反演优化目标函数；
5)利用阻尼最小二乘法估计反射系数初值R0；
6)构造预处理降维矩阵A；
7)确定当前迭代步的零范数逼近参数σk，构造优化问题对应的Lagrange 函数L(R(k),λ(k))，计算目标函数的梯度gk，以及Lagrange函数的海森阵的正定近似矩阵Bk；
8)利用信赖域技术求解二次规划子问题Q(d,λ)，寻求当前迭代点R(k)的修正量d(k)并改进相应的Lagrange乘子λk；
9)进行迭代k＝k+1；判断是否满足终止条件，若不满足，则返回步骤7)；否则，执行步骤10)；
10)输出最终的反射系数剖面R*。
上述步骤3)中，利用快速衰减的指数函数来逐渐逼近零范数，该函数的表达式如下
式中，σ为高斯函数的标准差，当σ趋近于0时，fσ(x)趋近于零范数。
上述步骤4)中，优化目标函数h(R,σ,μ)的表达式为
式中，M为地震道总道数，N为时间采样点数，i为时间离散指标，j为道离散指标，μ为正则化参数，由GCV准则确定，ε2为噪声能量水平，R为反射系数剖面，Ri,j为第j道第i个时间采样点反射系数。
约束条件c(R)的表达式为
式中，W为地震子波w(t)张成的Toepolitz矩阵，S为叠后地震记录。
进而，多道地震记录反射系数序列同时反演优化目标函数可简写为如下形式
上述步骤6)中，基于高斯随机分布构造降维矩阵A，该矩阵的行数远小于列数，然后对其进行谱范数归一化。记降维处理后的地震记录和Toeplitz 矩阵分别为
上述步骤7)中，经过降维处理后，约束条件可改写为
此时，优化问题对应的Lagrange函数L(R(k),λ(k))的表达式为
L(R(k),λ(k))＝h(R,σk,μ)-λc(R)
式中，k表示迭代步数。目标函数的梯度Bk由BFGS校正公式生成。
上述步骤8)中，二次规划子问题Q(d,λ)如下所示
式中，d为当前迭代点R(k)的修正量，ρk为信赖域半径。
本发明具有以下有益效果：
本发明从地震记录褶积模型出发，充分考虑了临近地震道之间的空间相关性和反射系数的稀疏特征，建立了一种新的多道地震记录反射系数同时反演的约束优化目标函数。利用参数化的光滑指数函数逐渐逼近零范数，以度量反射系数的稀疏特征。基于褶积模型的特点和反射系数的稀疏性假设，在不丢失任何信息的前提下，设计降维预处理矩阵，降低约束目标函数的维度。在算法实现上，利用信赖域序列二次规划算法对优化问题进行高效求解。该算法具有全局收敛性和超线性收敛率，计算效率高。相比于传统的单道地震记录反褶积算法，本发明算法考虑了地震道的空间相关性，反演结果更稳定、可靠。
附图说明
图1是本发明流程示意图；
图2是地震子波估计示意图；
图3是某陆上实际叠后地震资料。
图4是地震子波估计结果。
图5是本发明反射系数反演剖面。
具体实施方式
下面结合附图对本发明做进一步详细描述：
反射系数序列反演方法通过压缩地震子波，进而获取地层反射系数序列，是地震资料处理中提高纵向分辨率、精细刻画地层结构常用方法之一。本发明基于反射系数稀疏性假设，充分考虑地震道的空间相关性，建立多道地震记录反射系数同时反演模型。在算法实现上，利用信赖域二次规划算法对该模型进行高效求解，反演结果稳定、可靠。
本发明的物质基础是通过野外地震采集设备采集到的大量地震数据。本发明的理论基础是地震记录褶积模型。
在一定的假设条件之下，地震记录可以看做是地震子波与地层反射系数序列褶积结果。若假设地震子波在一定时空区域内是平稳的，即地震子波与空间和时间无关，此时，离散形式的多道褶积模型可写为
S＝WR+E
式中，R∈RN×M为反射系数序列，共M道，N个时间采样点；若地震子波的长度为l，则W∈RL×N为地震子波构成的Toeplitz矩阵，L＝N+l-1；E∈RL×M为加性噪声；S∈RL×M为观测到的叠后地震记录。
本发明的多道地震记录反射系数序列同时反演框架，具体步骤分别为：
1)采集原始地震资料，然后对采集到的地震资料做常规预处理，包括静校正处理、噪音压制处理、真振幅恢复处理。处理后得到叠前共中心点道集或叠后地震数据，称这个地震数据为观测地震数据，记为S(x,t)，其中x表示CDP号，t表示时间轴。
2)基于统计性方法提取地震子波w(t)。先对叠后地震记录S(x,t)逐道做傅里叶变换，然后计算多道平均振幅谱，利用多项式拟合地震子波的振幅谱，然后利用峭度最大化准则估计地震子波的相位，结合估计得到的地震子波振幅谱和相位，通过逆傅里叶变换得到地震子波w(t)。提取地震子波的流程如图2所示。
3)构造逼近零范数的光滑函数fσ(·)，并给定初始值σ。一个序列的零范数定义为该序列中非零元素的个数。尽管零范数是度量反射系数序列稀疏特征最有力的工具，但是零范数度量非凸非光滑，求解难度大，计算量高。为了后续迭代算法设计的需要，需要选取合适的函数对零范数进行近似，这里选用如下所示的高斯函数
式中，σ为高斯函数的标准差，当σ趋近于0时，fσ(x)趋近于零范数。σ在后续迭代中可动态选取，这里设定其初值为0.5。
4)首先，通过广义交叉确认法，估计正则化参数μ。地震记录噪声能量ε2可通过对地震记录做小波变换分析得到。记反射系数剖面为R，为了后续叙述方便，记优化目标函数为h(R,σ,μ)，其表达式，
式中，i为时间离散指标，j为道离散指标，μ为正则化参数，由GCV准则确定，ε2为噪声能量水平，Ri,j为第j道第i个时间采样点反射系数。
类似地，为了后续叙述方便，记约束条件为c(R)，其表达式为
最终，本发明建立多道地震记录反射系数序列同时反演约束优化目标函数如下所示
5)利用阻尼最小二乘法估计反射系数R0，作为后续迭代的初值。阻尼最小二乘法估计公式为
R0＝(WTW+δI)-1WTS
式中，T为矩阵转置符号，I为单位矩阵，δ为阻尼参数，为了保证初值估计的稳定性，可取为0.05。
6)构造预处理降维矩阵A。由于褶积模型是一个超定矩阵，即方程个数大于未知量的维度，有必要对褶积模型进行降维处理。此外，由于反射系数序列满足稀疏性假设，在不丢失任何信息的前提下，可对褶积模型 S＝WR+E两边同时乘以预处理矩阵A，其中，A∈RK×N是行满秩的，其维度满足K＜N＜L，且||A||2＝1。记降维处理后的地震记录和Toeplitz矩阵分别为则褶积模型的维度得到有效降低。此外，由于||A||2＝1，则||AE||2≤||A||2||E||2≤||E||2。
7)确定当前迭代步的零范数逼近参数σk，构造优化问题对应的Lagrange 函数L(R(k),λ(k))，计算目标函数的梯度gk，以及Lagrange函数的海森阵的正定近似矩阵Bk。经过降维处理后，约束条件可改写为
此时，优化问题对应的Lagrange函数L(R(k),λ(k))的表达式为
L(R(k),λ(k))＝h(R,σk,μ)-λc(R)
式中，k表示迭代步数。目标函数的梯度Bk由BFGS校正公式生成
式中，yk＝gk-gk-1，ηk＝R(k)-R(k-1)，trc(·)表示矩阵的迹。当k＝1时，B1＝I 为单位矩阵。
8)利用信赖域技术求解二次规划子问题Q(d,λ)，寻求当前迭代点R(k)的修正量d(k)并改进相应的Lagrange乘子λk；
二次规划子问题Q(d,λ)如下所示
式中，d为当前迭代点R(k)的修正量，ρk为信赖域半径。利用预处理共轭梯度法求取信赖域子问题Q(d,λ)，获取最优的修正量d(k)。
9)进行迭代R(k+1)＝R(k)+d(k)，k＝k+1；判断是否满足终止条件
式中，ξ为误差容许阈值，可取为0.0001。
若不满足，则返回步骤7)；否则，R*＝R(k+1)，执行步骤10)；
10)输出最终的反射系数剖面R*。
实际地震资料测试
本节利用实际地震资料来验证本发明提出的多道地震记录反射系数同时反演算法的有效性。图3为某陆上实际叠后地震资料中的一条测线。该测线共1751 道，道间距20m，截取时间范围为1400～2200ms，时间采样间隔为2ms。测井位于CDP 1020处。可见，地震资料中采集脚印干扰严重，信噪比较差。利用图2给出流程估计地震子波，估计结果如图4所示。图5为本发明提出的多道地震记录反射系数同时反演算法反演得到的反射系数剖面。可见，反射系数剖面结构清晰、空间连续性好，且具有较高的纵向分辨率，实现对地质结构得到精细刻画。
Claims (6)
1.一种多道地震记录反射系数序列同时反演方法，其特征在于，包括以下步骤：
1)采集原始地震资料，对采集到的地震资料进行预处理，包括静校正处理、噪音压制处理、真振幅恢复处理，得到叠后地震数据，称这个地震数据为观测地震数据，记为S(x,t)，其中x表示CDP号，t表示时间轴；
2)基于统计性方法提取地震子波w(t)；
3)利用快速衰减的指数函数fσ(·)来逐渐逼近零范数，并给定初始值σ；
4)确定正则化参数μ，估计地震记录噪声能量ε2，记反射系数剖面为R，记优化目标函数为h(R,σ,μ)，记约束条件为c(R)，建立多道地震记录反射系数序列同时反演的优化目标函数；
5)利用阻尼最小二乘法估计反射系数初值R0；
6)构造预处理降维矩阵A；
7)确定当前迭代步的零范数逼近参数σk，构造优化问题对应的Lagrange函数L(R(k),λ(k))，计算目标函数的梯度gk，以及Lagrange函数的海森阵的正定近似矩阵Bk；
8)利用信赖域技术求解二次规划子问题Q(d,λ)，寻求当前迭代点R(k)的修正量d(k)并改进相应的Lagrange乘子λk；
9)进行迭代k＝k+1；判断是否满足终止条件，若不满足，则返回步骤7)；否则，执行步骤10)；
10)输出最终的反射系数剖面R*。
2.根据权利要求1所述的一种多道地震记录反射系数序列同时反演方法，其特征在于，所述步骤3)中，利用快速衰减的指数函数来逐渐逼近零范数，该函数的表达式如下
<mrow>
<msub>
<mi>f</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</msup>
</mrow>
式中，σ为高斯函数的标准差，当σ趋近于0时，fσ(x)趋近于零范数。σ的初值可去较大数值，可取值为0.5，以保证目标函数的光滑性。
3.根据权利要求1所述的一种多道地震记录反射系数序列同时反演方法，其特征在于，上述步骤4)中，优化目标函数h(R,σ,μ)的表达式为
<mrow>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>,</mo>
<mi>&sigma;</mi>
<mo>,</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>&mu;</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
式中，M为地震道总道数，N为时间采样点数，i为时间离散指标，j为道离散指标，μ为正则化参数，由GCV准则确定，ε2为噪声能量水平，R为反射系数剖面，Ri,j为第j道第i个时间采样点反射系数；
约束条件c(R)的表达式为
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>|</mo>
<mo>|</mo>
<mi>S</mi>
<mo>-</mo>
<mi>W</mi>
<mi>R</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
式中，W为地震子波w(t)张成的Toepolitz矩阵，S为叠后地震记录；
进而，多道地震记录反射系数序列同时反演优化目标函数可简写为如下形式
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>R</mi>
</munder>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>,</mo>
<mi>&sigma;</mi>
<mo>,</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
4.根据权利要求1所述的一种多道地震记录反射系数序列同时反演方法，其特征在于，上述步骤6)中，基于高斯随机分布构造降维矩阵A，该矩阵的行数远小于列数，然后对其进行谱范数归一化，记降维处理后的地震记录和Toeplitz矩阵分别为
5.根据权利要求1所述的一种多道地震记录反射系数序列同时反演方法，其特征在于，上述步骤7)中，经过降维处理后，约束条件可改写为
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>|</mo>
<mo>|</mo>
<mover>
<mi>S</mi>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<mover>
<mi>W</mi>
<mo>&OverBar;</mo>
</mover>
<mi>R</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
此时，优化问题对应的Lagrange函数L(R(k)，λ(k))的表达式为
L(R(k),λ(k))＝h(R,σk,μ)-λc(R)
式中，k表示迭代步数，目标函数的梯度Bk由BFGS校正公式生成。
6.根据权利要求1所述的一种多道地震记录反射系数序列同时反演方法，其特征在于，上述步骤8)中，二次规划子问题Q(d,λ)如下所示
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>d</mi>
</munder>
</mtd>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>g</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<mi>d</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>d</mi>
<mi>T</mi>
</msup>
<msub>
<mi>B</mi>
<mi>k</mi>
</msub>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>d</mi>
<mi>T</mi>
</msup>
<mo>&dtri;</mo>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>d</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>&le;</mo>
<msub>
<mi>&rho;</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中，d为当前迭代点R(k)的修正量，ρk为信赖域半径。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710572136.3A CN107589448B (zh) | 2017-07-13 | 2017-07-13 | 一种多道地震记录反射系数序列同时反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710572136.3A CN107589448B (zh) | 2017-07-13 | 2017-07-13 | 一种多道地震记录反射系数序列同时反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107589448A true CN107589448A (zh) | 2018-01-16 |
| CN107589448B CN107589448B (zh) | 2019-02-05 |
Family
ID=61042971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710572136.3A Active CN107589448B (zh) | 2017-07-13 | 2017-07-13 | 一种多道地震记录反射系数序列同时反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107589448B (zh) |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108387935A (zh) * | 2018-02-08 | 2018-08-10 | 中国地质调查局油气资源调查中心 | 一种地震叠前反演的方法及装置 |
| CN108957404A (zh) * | 2018-05-18 | 2018-12-07 | 中国科学技术大学 | 一种地震定位方法、装置及系统 |
| CN109212602A (zh) * | 2018-09-05 | 2019-01-15 | 湖南科技大学 | 一种改进地震数据分辨率的反射系数反演方法 |
| CN109901216A (zh) * | 2019-02-22 | 2019-06-18 | 山东省科学院海洋仪器仪表研究所 | 一种检测海水放射性核素的寻峰方法 |
| CN110308483A (zh) * | 2019-05-23 | 2019-10-08 | 中国石油天然气股份有限公司 | 基于多任务贝叶斯压缩感知的反射系数求取方法及装置 |
| CN110471113A (zh) * | 2019-08-01 | 2019-11-19 | 中国石油大学(北京) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 |
| CN110618453A (zh) * | 2019-08-07 | 2019-12-27 | 成都理工大学 | 一种基于改进阻尼最小二乘法的波阻抗反演方法 |
| CN111856568A (zh) * | 2019-04-30 | 2020-10-30 | 中国石油天然气股份有限公司 | 基于mwv模型的频率域多道反射系数联合反演方法及系统 |
| CN112415587A (zh) * | 2019-08-21 | 2021-02-26 | 中国石油化工股份有限公司 | 储层地震波衰减特性分析方法及储层反射系数的反演方法 |
| CN112444846A (zh) * | 2019-08-27 | 2021-03-05 | 中国石油天然气集团有限公司 | 地震资料中浅水层影响优化方法及装置 |
| CN112526599A (zh) * | 2019-09-17 | 2021-03-19 | 中国石油化工股份有限公司 | 基于加权l1范数稀疏准则的子波相位估计方法及系统 |
| CN112578439A (zh) * | 2019-09-29 | 2021-03-30 | 中国石油化工股份有限公司 | 一种基于空间约束的地震反演方法 |
| CN113341463A (zh) * | 2021-06-10 | 2021-09-03 | 中国石油大学(北京) | 一种叠前地震资料非平稳盲反褶积方法及相关组件 |
| CN114167492A (zh) * | 2020-09-10 | 2022-03-11 | 中国石油天然气股份有限公司 | 基于自适应地质结构约束的多道稀疏反褶积方法及装置 |
| CN114397700A (zh) * | 2022-01-26 | 2022-04-26 | 西安交通大学 | 一种基于图信号约束的节点地震仪叠前地震数据插值方法、装置、设备及存储介质 |
| CN115184986A (zh) * | 2022-06-28 | 2022-10-14 | 吉林大学 | 不依赖震源的全局包络互相关全波形反演方法 |
| CN115373024A (zh) * | 2022-08-09 | 2022-11-22 | 中国科学院南海海洋研究所 | 基于地层记录沉降反演被动陆缘地壳结构的方法及装置 |
| CN116027406A (zh) * | 2023-02-10 | 2023-04-28 | 成都理工大学 | 提高反演分辨率的多道同时反演识别方法、装置及介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004095074A2 (en) * | 2003-04-01 | 2004-11-04 | Conocophillips Company | Application for letters patent for simultaneous inversion for source wavelet and avo parameters from prestack seismic data |
| CN103018773A (zh) * | 2012-11-15 | 2013-04-03 | 中国石油天然气股份有限公司 | 测井信号驱动的地震记录有色补偿处理方法 |
| CN103792571A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 点约束贝叶斯稀疏脉冲反演方法 |
| CN103954992A (zh) * | 2014-04-03 | 2014-07-30 | 中国石油天然气股份有限公司 | 一种反褶积方法及装置 |
| CN104062680A (zh) * | 2013-03-22 | 2014-09-24 | 中国石油化工股份有限公司 | 一种计算波阻抗反演目标函数梯度的方法 |
-
2017
- 2017-07-13 CN CN201710572136.3A patent/CN107589448B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004095074A2 (en) * | 2003-04-01 | 2004-11-04 | Conocophillips Company | Application for letters patent for simultaneous inversion for source wavelet and avo parameters from prestack seismic data |
| CN103792571A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 点约束贝叶斯稀疏脉冲反演方法 |
| CN103018773A (zh) * | 2012-11-15 | 2013-04-03 | 中国石油天然气股份有限公司 | 测井信号驱动的地震记录有色补偿处理方法 |
| CN104062680A (zh) * | 2013-03-22 | 2014-09-24 | 中国石油化工股份有限公司 | 一种计算波阻抗反演目标函数梯度的方法 |
| CN103954992A (zh) * | 2014-04-03 | 2014-07-30 | 中国石油天然气股份有限公司 | 一种反褶积方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 汪玲玲 等: "基于矩阵Toeplitz 稀疏分解的相对波阻抗反演方法", 《地球物理学报》 * |
Cited By (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108387935A (zh) * | 2018-02-08 | 2018-08-10 | 中国地质调查局油气资源调查中心 | 一种地震叠前反演的方法及装置 |
| CN108957404A (zh) * | 2018-05-18 | 2018-12-07 | 中国科学技术大学 | 一种地震定位方法、装置及系统 |
| CN109212602A (zh) * | 2018-09-05 | 2019-01-15 | 湖南科技大学 | 一种改进地震数据分辨率的反射系数反演方法 |
| CN109212602B (zh) * | 2018-09-05 | 2019-11-08 | 湖南科技大学 | 一种改进地震数据分辨率的反射系数反演方法 |
| CN109901216A (zh) * | 2019-02-22 | 2019-06-18 | 山东省科学院海洋仪器仪表研究所 | 一种检测海水放射性核素的寻峰方法 |
| CN111856568A (zh) * | 2019-04-30 | 2020-10-30 | 中国石油天然气股份有限公司 | 基于mwv模型的频率域多道反射系数联合反演方法及系统 |
| CN110308483A (zh) * | 2019-05-23 | 2019-10-08 | 中国石油天然气股份有限公司 | 基于多任务贝叶斯压缩感知的反射系数求取方法及装置 |
| CN110471113B (zh) * | 2019-08-01 | 2020-08-04 | 中国石油大学(北京) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 |
| CN110471113A (zh) * | 2019-08-01 | 2019-11-19 | 中国石油大学(北京) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 |
| CN110618453A (zh) * | 2019-08-07 | 2019-12-27 | 成都理工大学 | 一种基于改进阻尼最小二乘法的波阻抗反演方法 |
| CN112415587A (zh) * | 2019-08-21 | 2021-02-26 | 中国石油化工股份有限公司 | 储层地震波衰减特性分析方法及储层反射系数的反演方法 |
| CN112444846A (zh) * | 2019-08-27 | 2021-03-05 | 中国石油天然气集团有限公司 | 地震资料中浅水层影响优化方法及装置 |
| CN112444846B (zh) * | 2019-08-27 | 2023-11-28 | 中国石油天然气集团有限公司 | 地震资料中浅水层影响优化方法及装置 |
| CN112526599A (zh) * | 2019-09-17 | 2021-03-19 | 中国石油化工股份有限公司 | 基于加权l1范数稀疏准则的子波相位估计方法及系统 |
| CN112526599B (zh) * | 2019-09-17 | 2024-04-09 | 中国石油化工股份有限公司 | 基于加权l1范数稀疏准则的子波相位估计方法及系统 |
| CN112578439B (zh) * | 2019-09-29 | 2023-06-02 | 中国石油化工股份有限公司 | 一种基于空间约束的地震反演方法 |
| CN112578439A (zh) * | 2019-09-29 | 2021-03-30 | 中国石油化工股份有限公司 | 一种基于空间约束的地震反演方法 |
| CN114167492A (zh) * | 2020-09-10 | 2022-03-11 | 中国石油天然气股份有限公司 | 基于自适应地质结构约束的多道稀疏反褶积方法及装置 |
| CN114167492B (zh) * | 2020-09-10 | 2023-08-22 | 中国石油天然气股份有限公司 | 基于自适应地质结构约束的多道稀疏反褶积方法及装置 |
| CN113341463A (zh) * | 2021-06-10 | 2021-09-03 | 中国石油大学(北京) | 一种叠前地震资料非平稳盲反褶积方法及相关组件 |
| CN113341463B (zh) * | 2021-06-10 | 2023-05-26 | 中国石油大学(北京) | 一种叠前地震资料非平稳盲反褶积方法及相关组件 |
| CN114397700B (zh) * | 2022-01-26 | 2023-08-22 | 西安交通大学 | 一种基于图信号约束的节点地震仪叠前地震数据插值方法、装置、设备及存储介质 |
| CN114397700A (zh) * | 2022-01-26 | 2022-04-26 | 西安交通大学 | 一种基于图信号约束的节点地震仪叠前地震数据插值方法、装置、设备及存储介质 |
| CN115184986B (zh) * | 2022-06-28 | 2023-09-29 | 吉林大学 | 不依赖震源的全局包络互相关全波形反演方法 |
| CN115184986A (zh) * | 2022-06-28 | 2022-10-14 | 吉林大学 | 不依赖震源的全局包络互相关全波形反演方法 |
| CN115373024B (zh) * | 2022-08-09 | 2023-04-18 | 中国科学院南海海洋研究所 | 基于地层记录沉降反演被动陆缘地壳结构的方法及装置 |
| CN115373024A (zh) * | 2022-08-09 | 2022-11-22 | 中国科学院南海海洋研究所 | 基于地层记录沉降反演被动陆缘地壳结构的方法及装置 |
| CN116027406A (zh) * | 2023-02-10 | 2023-04-28 | 成都理工大学 | 提高反演分辨率的多道同时反演识别方法、装置及介质 |
| CN116027406B (zh) * | 2023-02-10 | 2023-08-18 | 成都理工大学 | 提高反演分辨率的多道同时反演识别方法、装置及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107589448B (zh) | 2019-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107589448B (zh) | 一种多道地震记录反射系数序列同时反演方法 | |
| CN113610945B (zh) | 一种基于混合神经网络的地应力曲线预测方法 | |
| CN103954992B (zh) | 一种反褶积方法及装置 | |
| CN104122588A (zh) | 基于谱分解提高叠后地震资料分辨率的方法 | |
| CN101206264A (zh) | 高分辨率非线性地震波阻抗反演方法 | |
| CN103454677B (zh) | 基于粒子群与线性加法器结合的地震数据反演方法 | |
| CN110895348B (zh) | 一种地震弹性阻抗低频信息提取方法、系统及存储介质 | |
| CN104422969B (zh) | 一种减小电磁测深反演结果非唯一性的方法 | |
| CN110187384A (zh) | 贝叶斯时移地震差异反演方法及装置 | |
| CN111722283B (zh) | 一种地层速度模型建立方法 | |
| CN109738950A (zh) | 基于三维稀疏聚焦域反演的噪声型数据一次波反演方法 | |
| CN103217715B (zh) | 多尺度规则网格层析反演静校正方法 | |
| CN116859478B (zh) | 一种基于瞬变电磁法成像的地下水模拟方法及系统 | |
| CN111914609B (zh) | 井震联合叠前地质统计学弹性参数反演方法及装置 | |
| CN112100906B (zh) | 数据驱动的大尺度密度建模方法、计算设备及存储介质 | |
| CN111273346B (zh) | 去除沉积背景的方法、装置、计算机设备及可读存储介质 | |
| CN116068644B (zh) | 一种利用生成对抗网络提升地震数据分辨率和降噪的方法 | |
| CN114942473B (zh) | 基于注意力门神经网络的叠前地震速度反演方法 | |
| Malvić et al. | Neural networks in petroleum geology as interpretation tools | |
| CN116011338A (zh) | 一种基于自编码器和深度神经网络的全波形反演方法 | |
| CN107678065B (zh) | 提高地震分辨率的保构造井控空间反褶积方法和装置 | |
| CN111175824B (zh) | 岩相驱动下的时频联合域地震反演方法 | |
| Wu et al. | A spatially constrained divisive hierarchical k-means clustering to capture prior features from migration velocity model to build training model set for deep-learning LSRTM | |
| CN117388921B (zh) | 弹性参数的叠前反演方法、装置和电子设备 | |
| CN118295016B (zh) | 一种多元参数耦合数据－知识双驱动地震反演方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |