[发明专利]一种地震数据子波相位转换处理的方法

摘要

本发明涉及石油勘探开发的一种地震数据子波相位转换处理的方法，计算地震数据零相位子波，Z变换后分解成N个相乘的因式，将特征因子按照绝对值从小到大依次排列，按绝对值的大小将特征因子分成单位圆内和单位圆外两部分，构成最小相位子波的纯相位转换因子，将零相位子波转换为最小相位子波，再波转换为零相位子波，将零相位子波地震数据转换为最小相位子波地震数据，再转换为零相位子波地震数据，实现地震子波和地震数据的纯相位转换。本发明通过反褶积处理可以有效提高地震数据分辨率，可以有效提高地震数据视分辨率，且计算简单。

主权项

1.一种地震数据子波相位转换处理的方法，其特征是包括以下步骤：1)用地震人工震源激发和采集地震数据并做预处理；2)确定地震数据子波参数和计算地震数据零相位子波；3)对子波进行Z变换；4)对子波Z变换按下式分解成N个相乘的因式：B(z)=Πi=1N(1-b^iz-1)=(1-b^1z-1)(1-b^2z-1)Λ(1-b^Nz-1)---(3)]]>且每个因式的特征因子按照绝对值从小到大依次排列，|b^1|≤|b^2|≤Λ≤|b^N|---(4)]]>5)将特征因子分成单位圆内和单位圆外两部分，按特征因子绝对值的大小，将特征因子分成两部分：第一部分是特征因子的绝对值小于或者等于1，共有K(K≥1)个，其中K是Z平面上单位圆内和单位圆上特征因子的总个数；第二部分是特征因子的绝对值大于1，共有M(M≥1)个，其中M是Z平面上单位圆外特征因子的总个数，且M+K＝N；6)根据单位圆外的特征因子，按下式构成将零相位子波转换为最小相位子波的纯相位转换因子H₁(z)：H1(z)=Πi=1M(1-b^K+i*z-1)(z-1-1b^K+i)---(6)]]>根据特征方程单位圆外的特征因子，按下式构成将最小相位子波转换为零相位子波的纯相位转换因子H₂(z)：H2(z)=Πi=1M(z-1-1b^K+i)(1-1b^K+i*z-1)---(7)]]>式中：“*”表示取共轭；7)根据最小相位子波的纯相位转换因子H₁(z)，按下式将零相位子波转换为最小相位子波：W_N(z)＝H₁(z)B_N(z)            (8)根据零相位子波的纯相位转换因子H₂(z)，按下式将最小相位子波转换为零相位子波：B_N(z)＝H₂(z)W_N(z)            (9)式中：W_N(z)表示最小相位子波的Z变换，B_N(z)表示零相位子波的Z变换；8)根据纯相位转换因子H₁(z)，按下式将零相位子波地震数据转换为最小相位子波地震数据：Y_W(z)＝H₁(z)X_B(z)            (10)式中：Y_W(z)表示最小相位子波地震数据的Z变换，X_B(z)表示零相位子波地震数据的Z变换；根据零相位子波的纯相位转换因子H₂(z)，按下式将最小相位子波地震数据转换为零相位子波地震数据：Y_B(z)＝H₂(z)X_W(z)            (11)式中：Y_B(z)表示零相位子波地震数据的Z变换，X_W(z)表示最小相位子波地震数据的Z变换；9)根据纯相位转换因子实现地震子波和地震数据的纯相位转换，绘制相位转换后的地震数据剖面和存储相位转换后的地震数据。