[发明专利]一种确定地震子波的方法

说明书

技术领域

本发明涉及石油勘探开发技术，具体是地震资料处理、解释中的一种确定地震子波的方法。

背景技术

地震子波的确定是地震勘探中一个非常关键的问题。在正演问题中，需要通过波动方程或褶积模型结合地震子波来形成正演模拟地震数据。在反演和反褶积问题中，也需要通过地震道提取一个地震子波，不同的地震子波得到不同的反演结果。

地震子波提取的基础是褶积模型，它包括地震子波、反射系数、含有噪声的地震道数据。地震子波提取方法是利用测井资料计算出反射系数，然后结合地震道由褶积模型求出地震子波。目前常用的井震联合子波确定方法是基于合成道与地震道残差能量（以下简称地震道残差能量）最小原理：

Σt=t1t=t2[x(t)-Στ=-pτ=qw(τ)r(t-τ)]2→min---(1)]]>

公式(1)中r(t)是由测井资料计算得到的时间域反射系数，其时窗范围为：t₁≤t≤t₂,x(t)为地震道振幅，其时窗范围为在t₁+p≤t≤t₂-q，w(t)为要确定的地震子波，其起至时间范围为-p到q，将公式（1）转化为矩阵运算形式可得到如下表达式：

E＝||RW-X||₂→min            (2)

由∂E∂W=0]]>得：

W＝R^-1X                     （3）

公式（3）即为目前常用的井震联合确定地震子波的计算公式，R为反射系数褶积矩阵。设定地震子波样点个数为n（ｎ为奇数），反射系数样点个数为m，反射系数褶积矩阵R为(m-n+1)×n阶矩阵，具体可表达为：