一种基于模型分割的波场重建方法

摘要

本发明涉及石油物探确定地下油气藏分布的波场重建方法，利用野外采集的叠加和偏移剖面和离散地层属性模型重建每一时刻各个离散网格点上的波场值，并将每一时刻检波点所在位置波场作为该时刻的波场重建记录，处理得到波场重建的叠加和偏移剖面，将波场重建得到的叠加和偏移剖面与野外采集数据处理得到的叠加和偏移剖面对比确定地下油气藏的分布。本发明适用于复杂构造波场重建，精度高、效率高和频散低，是确定地下油气藏的重要依据。

主权项

1.一种基于模型分割的波场重建方法，其特征在于采用以下步骤：1)利用野外观测系统采集地震数据；2)对野外采集的数据处理得到叠加和偏移剖面；3)利用步骤2得到的处理结果建立反映油气藏分布的离散地层属性模型；4)利用步骤1的野外观测系统和步骤3的离散地层属性模型，利用以下公式重建每一时刻各个离散网格点上的波场值，并将每一时刻检波点所在位置波场作为该时刻的波场重建记录；$u(t+\mathrm{\Delta t})=2u\left(t\right)-u(t-\mathrm{\Delta t})+{\left(\mathrm{v\Delta t}\right)}^2(\frac{{\partial }^{2}u}{{\partial x}^2}+\frac{{\partial }^{2}u}{{\partial y}^2}+\frac{{\partial }^{2}u}{{\partial z}^2})---\left(2\right)$$+\frac{1}{12}{\left(\mathrm{v\Delta t}\right)}^4(\frac{{\partial }^{4}u}{{\partial x}^4}+\frac{{\partial }^{4}u}{{\partial z}^4}+\frac{{\partial }^{4}u}{{\partial }^{4}y}+2\frac{{\partial }^{4}u}{{\partial y}^2{\partial z}^2}+2\frac{{\partial }^{4}u}{{\partial y}^2{\partial x}^2}+2\frac{{\partial }^{4}u}{{\partial x}^2{\partial z}^2})+O\left({\mathrm{\Delta t}}^6\right)$公式(2)中O(Δt⁶)是所述的公式在时间方向的差分截断误差。公式(2)中的所有空间导数都转化为二阶导数，并通过下式的10阶差分来近似：$\frac{{\partial }^{2}f}{{\partial s}^2}=\frac{1}{{\mathrm{\Delta s}}^2}\underset{m=1}{\overset{5}{\Sigma }}{C}_m^{\left(5\right)}[f(s+\mathrm{m\Delta x})-2f\left(s\right)+f(s-\mathrm{m\Delta s})]---\left(3\right)$公式(3)中差分系数通过下式求得：$\left[\begin{array}{ccccc}{1}^2& 2^2& 3^2& 4^2& 5^2\\ 1^4& 2^4& 3^4& 4^4& 5^4\\ 1^6& 2^6& 3^6& 4^6& 5^6\\ 1^8& 2^8& 3^8& 4^8& 5^8\\ 1^{10}& 2^{10}& 3^{10}& 4^{10}& 5^{10}\end{array}\right]\left[\begin{array}{c}{C}_1^{\left(5\right)}\\ C_2^{\left(5\right)}\\ C_3^{\left(5\right)}\\ C_4^{\left(5\right)}\\ C_5^{\left(5\right)}\end{array}\right]=\left[\begin{array}{c}1\\ 0\\ 0\\ 0\\ 0\end{array}\right]$5)将得到的记录按照步骤2)进行相同的处理，得到波场重建的叠加和偏移剖面；6)将波场重建得到的叠加和偏移剖面与步骤2)野外采集数据处理得到的叠加和偏移剖面对比，确定地下油气藏的分布。