一种地震资料偏移成像方法技术

本发明专利技术公开一种地震资料偏移成像方法。该方法包括：获取待进行偏移成像的地震资料数据、观测系统参数、偏移速度模型及偏移参数；将偏移速度模型由深度域变换到垂向时间域；将各向同性介质一阶速度‑应力声波方程由深度域变换到垂向时间域；获取每炮对应的每一个时刻的速度‑应力优化算子；利用所得的该炮对应的每一个时刻的速度‑应力优化算子进行震源波场重建，获得该炮的每一个时刻的震源波场；获取每炮的每一个时刻的检波点波场；获得每炮对应的单炮偏移剖面；将所有炮的所述单炮偏移剖面进行叠加，获得最终的偏移剖面。本发明专利技术方法，可以有效降低常规三维逆时偏移方法的计算量和存储量，提高算效率。

【技术实现步骤摘要】
一种地震资料偏移成像方法
本专利技术属于地震勘探
，特别是涉及一种地震资料偏移成像方法。
技术介绍
地震资料逆时偏移成像技术基于双程波动理论，不受倾角和横向变速的限制，克服了传统射线类偏移和单程波偏移方法的局限性，在复杂构造油气藏成像方面具有很好的应用前景和经济效益。对于逆时偏移方法，震源波场的存储是制约该方法在实际中应用的瓶颈。为了克服该存储问题，目前业内常见的几种方法技术为：Checkpoint技术、随机边界波场重建技术、有效边界存储技术。这几种技术均有各自的局限性：Checkpoint技术计算量较大；随机边界波场重建技术引入了额外的随机噪声影响成像质量；有效边界存储技术对于三维数据，存储量依然难以接受。针对目前逆时偏移方法所存在的存储量及计算效率问题，必须建立一套新的、以“高效率、高精度、低存储”为特点的地震资料偏移成像方法。
技术实现思路
本专利技术的目的在于提供一种地震资料偏移成像方法，以提高地震资料偏移成像的效率。为实现上述目的，本专利技术提供了如下方案：一种地震资料偏移成像方法，所述方法包括：读取地震工区观测系统参数、速度模型及观测的多炮地震资料数据，确定进行偏移成像的多炮观测地震记录、观测系统参数、深度域偏移速度模型及偏移参数；利用深度域-垂向时间域坐标变换关系式，将深度域坐标变换为垂向时间域，基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型；利用深度域-垂向时间域坐标变换关系式，将深度域的各向同性介质一阶速度-应力声波方程变换为垂向时间域的各向同性介质一阶速度-应力声波方程；根据所述观测系统参数获取每炮的炮点位置坐标；在每个炮点位置坐标处分别设置震源子波；根据所述垂向时间域偏移速度模型、垂向时间域平滑偏移速度模型和偏移参数，利用交错网格优先差分数值算法求解所述的垂向时间域一阶应力-速度声波方程，实现每炮震源波场的正向延拓，获得每炮的每一个时刻的震源波场；保存每炮的每一时刻的速度-应力优化算子；根据所述观测系统参数获取每炮的检波点位置坐标；根据所述垂向时间域偏移速度模型、垂向时间域平滑偏移速度模型和偏移参数，利用数值算法求解所述的垂向时间域一阶应力-速度声波方程，对每炮的所述观测地震记录进行逆时延拓；获得每炮的每一个时刻的检波点波场；读取存储的每炮的每一个时刻的速度-应力优化算子，利用读取的优化算子以及特定的差分格式，进行震源波场重建，获得每炮的每一个时刻的重建的震源波场；在相同时刻，基于偏移成像条件，对所述的每炮的每一个时刻的重建的震源波场和检波点波场进行成像，获得每炮的单炮偏移剖面；进而获得所有炮的单炮偏移剖面，将所有炮的单炮偏移剖面进行叠加，获得垂向时间域的最终的偏移剖面；利用深度域-垂向时间域坐标反变换关系，将所述垂向时间域的最终的偏移成像剖面由垂向时间域变换至深度域，得到深度域的最终的偏移剖面。可选的，所述利用深度域-垂向时间域坐标变换关系式，将深度域坐标变换为垂向时间域，基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型，具体为：深度域-垂向时间域坐标变换关系式为：式(1)中，x，y和z分别表示笛卡尔坐标系的x，y和z方向；τ(x,y,z)是由深度域坐标变量z转换所得的垂向时间域坐标变量，它与z是一一对应的；是由所述的深度域偏移速度模型vd(x,y,z)平滑所得的深度域平滑偏移速度模型；利用式(1)即可将深度域坐标变换为垂向时间域；所述基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型，具体包括：(a)根据深度域坐标z和垂向时间域坐标τ之间的一一对应关系，所述深度域平滑偏移速度模型即为原始垂向时间域的偏移速度模型vτ1(x,y,τ)，此时垂向时间域坐标τ的坐标间隔为不均匀的；(b)以(a)中所述的垂向时间域的最大坐标间隔为准，利用三次样条插值算法对所述的垂向时间域坐标进行插值重采样，获得均匀的垂向时间域坐标，同时即可获得垂向时间域偏移速度模型vτ(x,y,τ)；同理，由深度域平滑偏移速度模型可以获得垂向时间域的平滑偏移速度模型可选的，所述利用深度域-垂向时间域坐标变换关系式，将深度域的各向同性介质一阶速度-应力声波方程变换为垂向时间域的各向同性介质一阶速度-应力声波方程，具体为：深度域各向同性介质一阶速度-应力声波方程具体为：式(2)中，ux，uy和uz分别为深度域各向同性介质声波质点速度波场在x，y和z三个方向的分量，p为深度域各向同性介质声波正应力波场；所述将深度域的各向同性介质一阶速度-应力声波方程变换为垂向时间域的各向同性介质一阶速度-应力声波方程，具体为：(a)将深度域坐标系(x,y,z)变换为垂向时间域坐标系两个坐标系的坐标变量满足如下关系：(b)引入符号fx表示：定义垂向时间域算子α＝τx和β＝τy，根据链式法则，笛卡尔坐标系的偏导数与垂向时间域坐标系的偏导数满足如下关系：那么笛卡尔坐标系的单位基矢量与曲坐标系的单位基矢量之间的关系为：(c)笛卡尔坐标系下的任意矢量其在垂向时间域对应形式为：其中，Vx，Vy和Vz为矢量V在笛卡尔坐标系在x，y和z三个方向的分量；Vξ，和Vτ为矢量V在垂向时间域坐标系ξ，和τ三个方向的分量；运用散度定理，将笛卡尔坐标系下的矢量V的梯度和散度用其垂向时间域的分量及单位基矢量进行表示，具体为：(d)将式子(8)，(9)代入式子(2)，即可得到各向同性介质垂向时域一阶速度-应力声波方程，具体为：式(10)中，Uξ，Uτ分别为垂向时间域各向同性介质声波质点速度波场在垂向时间域坐标系ξ，τ三个方向的分量，P为垂向时间域各向同性介质声波正应力波场。可选的，所述保存每炮的每一时刻的速度-应力优化算子，具体为：任意t时刻，速度-应力优化算子At定义如下：式(11)中，为t时刻模型区域外第i个网格点上的垂向时间域声波正应力波场；其中i＝1,2,…,N,表示波场延拓所用的空间差分半阶数，N表示最大半阶数；h表示空间网格间距；ai为第一组优化算子系数。可选的，所述利用读取的优化算子以及特定的差分格式，进行震源波场重建，获得每炮的每一个时刻的重建的震源波场，具体为；以正应力波场P为例，垂向时间域各向同性介质一阶应力-速度声波方程中的空间偏导数的交错网格有限差分形式为：以P在x＝0(即左边界)为例，P在x＝0右侧的第i个点的一阶偏导数可用有限差分进行表示，具体为：其中，bi,j为第二组优化算子系数，式(12)即为所述的特定差分格式，将式子(12)和式子(11)对比可知，式子(12)包含了所述的优化算子At，由式(12)式即可求出波场在左边界内N-1点的偏导数；同理，可以获得模型其余边界内N-1点的偏导数；对于模型内部的点，利用常规的交错网格有限差分计算其偏导数，从而实现震源波场重建，获得每炮的每一个时刻的重建的震源波场。可选的，所述在相同时刻，基于偏移成像条件，对所述的每炮的每一个时刻的重建的震源波场和检波点波场进行成像，获得每炮的单炮偏移剖面；进而获得所有炮的单炮偏移剖面，将所有炮的单炮偏移剖面进行叠加，获得垂向时间域的最终的偏移剖面，具体为：所述成像本文档来自技高网...

【技术保护点】
1.一种地震资料偏移成像方法，其特征在于，所述方法包括：读取地震工区观测系统参数、速度模型及观测的多炮地震资料数据，确定进行偏移成像的多炮观测地震记录、观测系统参数、深度域偏移速度模型及偏移参数；利用深度域‑垂向时间域坐标变换关系式，将深度域坐标变换为垂向时间域，基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型；利用深度域‑垂向时间域坐标变换关系式，将深度域的各向同性介质一阶速度‑应力声波方程变换为垂向时间域的各向同性介质一阶速度‑应力声波方程；根据所述观测系统参数获取每炮的炮点位置坐标；在每个炮点位置坐标处分别设置震源子波；根据所述垂向时间域偏移速度模型、垂向时间域平滑偏移速度模型和偏移参数，利用交错网格优先差分数值算法求解所述的垂向时间域一阶应力‑速度声波方程，实现每炮震源波场的正向延拓，获得每炮的每一个时刻的震源波场；保存每炮的每一时刻的速度‑应力优化算子；读取存储的每炮的每一个时刻的速度‑应力优化算子，利用读取的优化算子以及特定的差分格式，进行震源波场重建，获得每炮的每一个时刻的重建的震源波场；根据所述观测系统参数获取每炮的检波点位置坐标；根据所述垂向时间域偏移速度模型、垂向时间域平滑偏移速度模型和偏移参数，利用数值算法求解所述的垂向时间域一阶应力‑速度声波方程，对每炮的所述观测地震记录进行逆时延拓；获得每炮的每一个时刻的检波点波场；在相同时刻，基于偏移成像条件，对所述的每炮的每一个时刻的重建的震源波场和检波点波场进行成像，获得每炮的单炮偏移剖面；进而获得所有炮的单炮偏移剖面，将所有炮的单炮偏移剖面进行叠加，获得垂向时间域的最终的偏移剖面；利用深度域‑垂向时间域坐标反变换关系，将所述垂向时间域的最终的偏移成像剖面由垂向时间域变换至深度域，得到深度域的最终的偏移剖面。...

【技术特征摘要】
1.一种地震资料偏移成像方法，其特征在于，所述方法包括：读取地震工区观测系统参数、速度模型及观测的多炮地震资料数据，确定进行偏移成像的多炮观测地震记录、观测系统参数、深度域偏移速度模型及偏移参数；利用深度域-垂向时间域坐标变换关系式，将深度域坐标变换为垂向时间域，基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型；利用深度域-垂向时间域坐标变换关系式，将深度域的各向同性介质一阶速度-应力声波方程变换为垂向时间域的各向同性介质一阶速度-应力声波方程；根据所述观测系统参数获取每炮的炮点位置坐标；在每个炮点位置坐标处分别设置震源子波；根据所述垂向时间域偏移速度模型、垂向时间域平滑偏移速度模型和偏移参数，利用交错网格优先差分数值算法求解所述的垂向时间域一阶应力-速度声波方程，实现每炮震源波场的正向延拓，获得每炮的每一个时刻的震源波场；保存每炮的每一时刻的速度-应力优化算子；读取存储的每炮的每一个时刻的速度-应力优化算子，利用读取的优化算子以及特定的差分格式，进行震源波场重建，获得每炮的每一个时刻的重建的震源波场；根据所述观测系统参数获取每炮的检波点位置坐标；根据所述垂向时间域偏移速度模型、垂向时间域平滑偏移速度模型和偏移参数，利用数值算法求解所述的垂向时间域一阶应力-速度声波方程，对每炮的所述观测地震记录进行逆时延拓；获得每炮的每一个时刻的检波点波场；在相同时刻，基于偏移成像条件，对所述的每炮的每一个时刻的重建的震源波场和检波点波场进行成像，获得每炮的单炮偏移剖面；进而获得所有炮的单炮偏移剖面，将所有炮的单炮偏移剖面进行叠加，获得垂向时间域的最终的偏移剖面；利用深度域-垂向时间域坐标反变换关系，将所述垂向时间域的最终的偏移成像剖面由垂向时间域变换至深度域，得到深度域的最终的偏移剖面。2.根据权利要求1所述一种地震资料偏移成像方法，其特征在于，所述利用深度域-垂向时间域坐标变换关系式，将深度域坐标变换为垂向时间域，基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型，具体为：深度域-垂向时间域坐标变换关系式为：式(1)中，x，y和z分别表示笛卡尔坐标系的x，y和z方向；τ(x,y,z)是由深度域坐标变量z转换所得的垂向时间域坐标变量，它与z是一一对应的；是由所述的深度域偏移速度模型vd(x,y,z)平滑所得的深度域平滑偏移速度模型；利用式(1)即可将深度域坐标变换为垂向时间域；所述基于所述深度域偏移速度模型，利用三次样条插值算法，获得垂向时间域的偏移速度模型，进而获得垂向时间域的平滑偏移速度模型，具体包括：(a)根据深度域坐标z和垂向时间域坐标τ之间的一一对应关系，所述深度域平滑偏移速度模型即为原始垂向时间域的偏移速度模型vτ1(x,y,τ)，此时垂向时间域坐标τ的坐标间隔为不均匀的；(b)以(a)中所述的垂向时间域的最大坐标间隔为准，利用三次样条插值算法对所述的垂向时间域坐标进行插值重采样，获得均匀的垂向时间域坐标，同时即可获得垂向时间域偏移速度模型vτ(x,y,τ)；同理，由深度域平滑偏移速度模型可以获得垂向时间域的平滑偏移速度模型3.按照权利要求1所述一种地震资料偏移成像方法，其特征在于：所述利用深度域-垂向时间域坐标变换关系式，将深度域的各向同性介质一阶速度-应力声波方程变换为垂向时间域的各向同性介质一阶速度-应力声波方程，具体为：深度域各向同性介质一阶速度-应力声波方程具体为：式(2)中，ux，uy和uz分别为深度域各向同性介质声波质点速度波场在x，y和z三个方向的分量，p为深度域各向同性介质声波正应力波场；所述将深度域的各向同性介质一阶速度-应力声波方程变换为...

【专利技术属性】
技术研发人员：谷丙洛，段沛然，李振春，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东,37