[发明专利]子波处理二步法反褶积方法

摘要

本发明涉及到石油地球物理勘探中地震资料处理，是一种子波处理二步法反褶积方法，按以下步骤进行：1.将采集于不同施工年代、不同施工参数的三维地震资料按常规方法进行分块预处理，包括静校正处理、振幅衰减补偿处理、噪音压制处理等；2.地震子波频宽处理，3.将地震子波频宽处理后的资料进行常规的速度分析与剩余静校正、叠加等处理；4.求取相位校正反褶积算子，利用式求相位校正反褶积算子m(t)，5.用求取的相位校正反褶积算子m(t)与需要调整相位区块的叠前地震道进行褶积，即可完成相位校正子波处理。它将不同区块的子波统一起来，以便连片处理，方法简单而有效。获取的相位校正反褶积算子稳定性高，对地震资料的保幅性较好。

主权项

1.一种子波处理二步法反褶积方法，其特征在于按以下步骤进行：第一步、将采集于不同施工年代、不同施工参数的三维地震资料按常规方法进行分块预处理，包括静校正处理、振幅衰减补偿处理、噪音压制处理等；第二步、地震子波频宽处理假设同一位置的不同时间采集的不同信号，包含震源不同信号，是反射系数与地震子波的褶积，第一块三维采集的地震信号为：s₁(t)＝w₁(t)*r(t)第二块三维采集的地震信号为：s₂(t)＝w₂(t)*r(t)……第N块三维采集的地震信号为：s_n(t)＝w_n(t)*r(t)            (1)其中s(t)是用检波器接收到的地震子波，w(t)为激发的震源子波，r(t)是地层反射系数，假定同一位置的地震反射系数r(t)是不变的，只是子波受到不同采集参数的影响；首先，不考虑子波相位的差异，认为子波只是频宽不同；根据子波的特征，采用反褶积算子a(t)，对地震记录s₁(t)和s₂(t)进行处理，相当于对地震子波w₁(t)、w₂(t)……w_n(t)进行处理，得到处理后的w_a1(t)、w_a2(t)……w_an(t)：s_a1(t)＝a(t)*s₁(t)＝a(t)*(r(t)*w₁(t))＝r(t)*(a(t)*w₁(t))＝r(t)*w_a1(t)s_a2(t)＝a(t)*s₂(t)＝a(t)*(r(t)*w₂(t))＝r(t)*(a(t)*w₂(t))＝r(t)*w_a2(t)……s_an(t)＝a(t)*s_n(t)＝a(t)*(r(t)*w_n(t))＝r(t)*(a(t)*w_n(t))            (2)＝r(t)*w_an(t)将反褶积算子a(t)设计成零相位反褶积算子，进行地震子波频宽处理；第三步、将地震子波频宽处理后的资料进行常规的速度分析与剩余静校正、叠加等处理；第四步、求取相位校正反褶积算子设s_a1(t)，s_a2(t)为经过第二步地震子波频宽处理后的任意地理位置上相邻的两个区块的叠前或叠后地震记录信号，s_a1(t)＝m(t)*s_a2(t)                   (3)式中，m(t)为相位校正反褶积算子；再利用式(4)求相位校正反褶积算子m(t)：Q=Σt=0K+M[sa1(t)-Στ=0mm(t)sa2(t-τ)]2---(4)]]>Q为趋于零的极小值，根据函数求极值的原理，令s＝0，1，2…，M  (5)  可列出M+1个方程组，组成托布利兹矩阵方程组：Στ=0Mraa(τ-s)sa1(τ)=rab(s)---(6)]]>式中：r_aa(τ-s)是时间延迟为τ-s的地震信号s_a1(t)自相关数列；r_ab(s)是时间延迟为s的地震信号s_a1(t)与s_a2(t)的互相关数列；用递推法解托布利兹矩阵方程(6)，得到相位校正反褶积算子m(t)；第五步、用求取的相位校正反褶积算子m(t)与需要调整相位区块的叠前地震道进行褶积，即可完成相位校正子波处理。