波场分离方法和装置

摘要

本发明提供了一种波场分离方法和装置，该方法包括：计算初始模型各向异性参数分布的变异函数值，拟合得到该初始模型的变异函数；利用所述变异函数和参考点搜索策略，从初始模型中选取N个参考模型；将空间域的矢量弹性波场变换至波数域；对选取的N个参考模型中的每个参考模型执行以下操作：截断拟微分算子，在波数域进行波场分离，将纵横波波场反变换回空间域，得到该参考模型下的波场分离结果；利用变异函数计算各个参考模型相对于初始模型的权重系数；在空间域对各个参考模型下的波场分离结果进行加权插值处理得到初始模型的波场分离结果。本发明实施达到了在保证波场分离的分离结果准确度的基础上，减少了计算量的效果。

主权项

一种波场分离方法，其特征在于，包括：计算初始模型各向异性参数分布的变异函数值，拟合得到该初始模型的变异函数；利用所述变异函数和参考点搜索策略，从初始模型中选取N个参考模型，其中，N为正整数；将空间域的矢量弹性波场变换至波数域；对选取的N个参考模型中的每个参考模型执行以下操作：根据参考模型计算得到拟微分算子，并采用自褶积窗函数截断拟微分算子；对该参考模型下的矢量弹性波场在波数域进行波场分离；将在波数域进行分离后的纵横波波场反变换回空间域，得到该参考模型下的波场分离结果；利用所述变异函数计算各个参考模型相对于初始模型的权重系数；根据计算得到的各个参考模型的权重系数，在空间域对各个参考模型下的波场分离结果进行加权插值处理；将插值处理后得到的结果作为初始模型的波场分离结果；其中，所述自褶积窗函数是按照以下方式得到的：选择主瓣和旁瓣性能高于预定阈值的窗函数作为原始窗函数；对所述原始窗函数做L次自褶积运算得到自褶积后的窗函数，其中，L为正整数；对自褶积后的窗函数与原始窗函数进行加权运算，得到所述自褶积窗函数；按照以下公式计算初始模型各向异性参数分布的变异函数值：$\mu \left(h\right)=\frac{1}{2n}\underset{i=1}{\overset{n}{\Sigma }}{[\Phi (ϵ,\delta ,\theta )-\Phi (ϵ+h_{ϵ},\delta +h_{\delta },\theta +h_{\theta })]}^2$其中，μ(h)表示变异函数值，Φ＝f[δ+2(ε‑δ)sin²(α‑θ)]sin2(α‑θ)，ε、δ表示横向各向同性TI介质各向异性的系数，α表示相角，传播方向与Z轴的夹角，代表的是各向同性部分，θ表示具有倾斜对称轴的横向各向同性TTI介质对称轴的倾角，表示插值参考点的位置增量，i＝1,2...n，其中，n表示插值参考点的个数；按照以下公式计算各个参考模型相对于初始模型的权重系数：$w^k=\frac{1/\mu ({\left({ϵ}_k-ϵ\right)}^2+{\left({\delta }_k-\delta \right)}^2+{\left({\theta }_k-\theta \right)}^2)}{\underset{k=1}{\overset{n}{\Sigma }}1/\mu ({\left({ϵ}_k-ϵ\right)}^2+{\left({\delta }_k-\delta \right)}^2+{\left({\theta }_k-\theta \right)}^2)}$其中，w^k表示权重系数，(ε_k,δ_k,θ_k)表示插值参考点。