一种多路径双谱域混合相位子波反褶积方法

摘要

公开了一种多路径双谱域混合相位子波反褶积方法。所述方法包括如下步骤：从采集到的反射地震记录中选取n道地震记录；根据选取的n道地震记录确定子波的振幅谱；根据选取的n道地震记录确定子波的相位谱；根据确定的子波的振幅谱和相位谱，进行反傅里叶变换来获得时域的混合相位子波；根据用户基于确定的子波的振幅谱输入的期望的输出子波的振幅谱，确定希望输出的零相位子波；以混合相位子波为输入以及零相位子波为输出，进行匹配滤波来确定反褶积滤波器；基于反褶积滤波器和地震记录，获得分辨率提高的地震记录。

主权项

1.一种多路径双谱域混合相位子波反褶积方法，其特征在于，包括如下步骤：a)从采集到的反射地震记录中选取n道地震记录x(ix,t)，其中，n为大于1的整数且ix=1,2,…,n；b)根据选取的n道地震记录确定子波的振幅谱w(f)；c)根据选取的n道地震记录确定子波的相位谱φ(ω)；d)根据确定的子波的振幅谱w(f)和相位谱φ(ω)，按如下的数学式进行反傅里叶变换来获得时域的混合相位子波w(t)，$w\left(t\right)=\underset{-\infty }{\overset{+\infty }{\int }}w\left(f\right)e^{\mathrm{i\phi t}}\mathrm{df};$e)根据用户基于确定的子波的振幅谱w(f)输入的期望的输出子波的振幅谱w′(f)，利用如下的数学式来确定希望输出的零相位子波w′(t)$w^{\prime }\left(t\right)=\underset{-\infty }{\overset{+\infty }{\int }}{w}^{\prime }\left(f\right)e^{i2\mathrm{\pi f}2}\mathrm{df};$f）以混合相位子波w(t)为输入以及零相位子波w′(t)为输出，进行匹配滤波来确定反褶积滤波器d(t)；g)基于反褶积滤波器d(t)和地震记录x(ix,t)，按如下的数学式来获得分辨率提高的地震记录y(ix,t)，y(ix,t)=d(t)*x(ix,t)，其中，*表示褶积，其中，在所述步骤b）中，对选取的n道地震记录进行傅里叶变换来获得每个地震记录的振幅谱x(ix,f)；计算n道地震记录的的振幅谱的平均值确定在最小二乘意义下使下面的数学式具有最小值时的a_m，从而确定子波的振幅谱w(f)，$Q=\underset{f_1}{\overset{f_2}{\Sigma }}{\left(\right|\bar{x}\left(f\right)|-|w\left(f\right)\left|\right)}^2,$其中，m的数值范围为，L为大于1的整数，其中，在所述步骤c）中，利用选取的n道地震记录和三阶累积量方式，确定双谱的相位谱ψ(ω₁,ω₂)，根据如下的数学式确定子波的相位谱φ(ω)，ψ(ω₁,ω₂)=φ(ω₁)+φ(ω₂)-φ(ω₁+ω₂)，其中，设置l=1，且c-1)按如下数学式计算其中，k=1,2,...2N-l,这里N是大于1的整数，这里，PV[x]是相位主算子，表示将[]中的值限定在模为π的(-π,π]之内；c-2)按如下数学式计算${\phi }_p^{l+1}(l+1)=\frac{1}{2N-2l+1}\underset{j=l}{\overset{2N-2l+1}{\Sigma }}{\phi }_p^l\left(j\right);$c-3）对l进行加1处理并判断l是否等于N；c-4）若l不等于N，则重复上述步骤c-1）、c-2）和c-3）；c-5）若l等于N，则按下述数学式计算φ(ω)，其中，ω为不小于0的整数，若ω大于1且不大于N，则若ω大于N且不大于2N-1，则