对应 的波场；
[0106] 检波点差分模拟模块806,用于通过金字塔网格的親合二阶偏微分方程对所述一 炮数据进行有限差分波场模拟得到检波点对应的波场；
[0107] 成像模块807,用于应用互相关成像条件对所述炮点对应的波场和所述检波点对 应的波场进行成像；
[0108] 插值模块808,用于将成像结果插值至规则网格，将插值至规则网格的成像结果作 为所述炮点的单炮成像结果。
[0109] 在一个实施方式中，第一变换模块802和第二变换模块803具体可以用于按照以下 公式将坐标由规则网格变换至金字塔网格：
[0110]
[0111] 其中，（x，y，z)表示规则网格的坐标系，（Xp，yp，zp)表示金字塔网格的坐标系，α表 示坐标变换系数。
[0112] 在一个实施方式中，金字塔网格的耦合二阶偏微分方程为：
[0113]
[0114]
[0115]
[0116] 其中
Ρ表不Ρ波 波场，q表示q波波场，νρζ表示Ρ波的速度，vsz表示S波的速度。
[0117] 在一个实施方式中，炮点差分模拟模块805在通过金字塔网格的耦合二阶偏微分 方程对激发产生的波场进行有限差分波场模拟得到所述炮点对应的波场的过程中，还按照 预定时间间隔对所述炮点对应的波场进行波场压缩，并将压缩后的结果存储在本地磁盘 中；成像模块807具体可以用于从所述本地磁盘中解压缩并读取所述炮点对应的波场;应用 互相关成像条件对所述检波点对应的波场和解压缩出来的所述炮点对应的波场，进行成 像。
[0118] 在一个实施方式中，插值模块808具体可以用于将所述成像结果在X方向、Y方向和 Z方向插值至所述规则网格中。
[0119] 从以上的描述中，可以看出，本发明实施例实现了如下技术效果:在本例中，通过 将地震数据由规则网格变换至金字塔网格，通过金字塔网格的TTI地震各向异性介质进行 逆时偏移成像，且基于金字塔网格的耦合二阶偏微分方程实现波动方程的差分求解，因此， 能够在不影响成像效果的前提下，显著提高计算效率、降低内存消耗，并最终解决速度急剧 变化的三维复杂构造成像问题。通过上述方式解决了现有技术中无法在水平方向和垂直方 向同时采用变化的网格进行成像的技术问题，达到了在水平方向和垂直方向同时采用变化 的网格进行成像的目的。
[0120] 显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以 用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置 所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它 们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执 行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个 模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬 件和软件结合。
[0121] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技 术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的 任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种逆时偏移成像方法，其特征在于，包括： SI:获取一炮数据； S2:读取所述一炮数据的炮点的成像空间对应的深度域速度场，并将所述深度域速度 场由规则网格变换至金字塔网格； S3:读取所述成像空间对应的各向异性参数S和e，并将所述各向异性参数S和e的坐标 由规则网格变换至金字塔网格； S4:读取所述成像空间对应的对称轴倾角参数0、对称轴方位角参数d)，并将所述对称 轴倾角参数e和对称轴方位角参数4的坐标由规则网格变换至金字塔网格； S5:通过放置一个子波对所述炮点进行激发，通过金字塔网格的禪合二阶偏微分方程 对激发产生的波场进行有限差分波场模拟得到所述炮点对应的波场； S6:通过金字塔网格的禪合二阶偏微分方程对所述一炮数据进行有限差分波场模拟得 到检波点对应的波场； S7:应用互相关成像条件对所述炮点对应的波场和所述检波点对应的波场进行成像； S8:将成像结果插值至规则网格，将插值至规则网格的成像结果作为所述炮点的单炮 成像结果； 对多个炮重复上述Sl至S8得到所述多个炮的单炮成像结果，将所述多个炮的单炮成像 结果进行叠加得到逆时偏移成像剖面。2. 如权利要求1所述的方法，其特征在于，按照W下公式将坐标由规则网格变换至金字 塔网格：其中，（x，y，z)表示规则网格的坐标系，（xp，yp，zp)表示金字塔网格的坐标系，a表示坐 标变换系数。3. 如权利要求2所述的方法，其特征在于，金字塔网格的禪合二阶偏微分方程为： 其中，其中，xdip = tan(目）*cos( 4 )，ydip = tan(目）*sin( 4 ) ，p表示P波波场， q表不q波波场，Vpz表不P波的速度，Vsz表不S波的速度。4. 如权利要求1至3中任一项所述的方法，其特征在于： 在通过金字塔网格的禪合二阶偏微分方程对激发产生的波场进行有限差分波场模拟 得到所述炮点对应的波场的过程中包括:按照预定时间间隔对所述炮点对应的波场进行波 场压缩，并将压缩后的结果存储在本地磁盘中； 所述应用互相关成像条件对所述炮点对应的波场和所述检波点对应的波场进行成像， 包括:从所述本地磁盘中解压缩并读取所述炮点对应的波场;应用互相关成像条件对所述 检波点对应的波场和解压缩出来的所述炮点对应的波场，进行成像。5. 如权利要求1至3中任一项所述的方法，其特征在于，将成像结果插值至规则网格包 括： 将所述成像结果在X方向、Y方向和Z方向插值至所述规则网格中。6. -种逆时偏移成像装置，其特征在于，包括： 获取模块，用于获取一炮数据； 第一变换模块，用于读取所述一炮数据的炮点的成像空间对应的深度域速度场，并将 所述深度域速度场由规则网格变换至金字塔网格； 第二变换模块，用于读取所述成像空间对应的各向异性参数S和e，并将所述各向异性 参数S和e的坐标由规则网格变换至金字塔网格； 第=变换模块，用于读取所述成像空间对应的对称轴倾角参数9、对称轴方位角参数 4，并将所述对称轴倾角参数e和对称轴方位角参数(I)的坐标由规则网格变换至金字塔网 格； 炮点差分模拟模块，用于通过放置一个子波对所述炮点进行激发，通过金字塔网格的 禪合二阶偏微分方程对激发产生的波场进行有限差分波场模拟得到所述炮点对应的波场； 检波点差分模拟模块，用于通过金字塔网格的禪合二阶偏微分方程对所述一炮数据进 行有限差分波场模拟得到检波点对应的波场； 成像模块，用于应用互相关成像条件对所述炮点对应的波场和所述检波点对应的波场 进行成像； 插值模块，用于将成像结果插值至规则网格，将插值至规则网格的成像结果作为所述 炮点的单炮成像结果。7. 如权利要求6所述的装置，其特征在于，所述第一变换模块和所述第二变换模块具体 用于按照W下公式将坐标由规则网格变换至金字塔网格：其中，（X，y，Z)表示规则网格的坐标系，（Xp，yp，Zp)表示金字塔网格的坐标系，a表示坐 标变换系数。8. 如权利要求7所述的装置，其特征在于，金字塔网格的禪合二阶偏微分方程为： 其中，其中，xdip = tan(目）*cos( 4 )，ydip = 1:an(目）*sin( 4 )，，:P表示P波波场， q表示q波波场，Vpz表示P波的速度，Vsz表示S波的速度。9. 如权利要求6至8中任一项所述的装置，其特征在于： 所述炮点差分模拟模块在通过金字塔网格的禪合二阶偏微分方程对激发产生的波场 进行有限差分波场模拟得到所述炮点对应的波场的过程中，还按照预定时间间隔对所述炮 点对应的波场进行波场压缩，并将压缩后的结果存储在本地磁盘中； 所述成像模块具体用于从所述本地磁盘中解压缩并读取所述炮点对应的波场;应用互 相关成像条件对所述检波点对应的波场和解压缩出来的所述炮点对应的波场，进行成像。10. 如权利要求6至8中任一项所述的装置，其特征在于，所述插值模块具体用于将所述 成像结果在X方向、Y方向和Z方向插值至所述规则网格中。
【专利摘要】本发明提供了一种逆时偏移成像方法和装置，其中，该方法包括：获取一炮数据；并将深度域速度场由规则网格变换至金字塔网格；将各向异性参数δ和ε、对称轴倾角参数θ和对称轴方位角参数φ的坐标由规则网格变换至金字塔网格；通过金字塔网格的耦合二阶偏微分方程对激发产生的波场进行有限差分波场模拟得到炮点对应的波场和检波点对应的波场；应用互相关成像条件对炮点对应的波场和检波点对应的波场进行成像；将成像结果插值至规则网格，将插值至规则网格的成像结果作为炮点的单炮成像结果。本发明解决了现有技术中无法在水平方向和垂直方向同时采用变化的网格进行成像的技术问题，达到了在水平方向和垂直方向同时采用变化的网格进行成像的目的。
【IPC分类】G01V1/34
【公开号】CN105652320
【申请号】
【发明人】戴南浔, 吴晓丰, 武威, 章晓
【申请人】中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司
【公开日】2016年6月8日
【申请日】2015年12月30日