一种自适应射线加密式克希霍夫型束偏移地震波成像方法技术

本发明专利技术提供一种自适应射线加密式克希霍夫型束偏移地震波成像方法，涉及地震偏移成像技术领域。首先读入地震记录、偏移速度模型和参数文件；然后从炮点按照自适应射线加密方式进行射线追踪并储存射线信息；在炮点完成全部射线追踪后，计算并储存射线对应的射线束信息；再将地震单炮记录分为多个以窗中心为核心的数据体，并对每一个数据体进行局部平面波分解；针对每一个窗中心依照自适应射线加密方式进行射线追踪并储存射线信息；在窗中心完成全部射线追踪后，计算并储存射线对应的射线束信息；同时计算每条射线束相应的偏移加权系数；选取炮点和窗中心点的射线束对进行偏移成像计算，累加所有射线束对的偏移成像结果获得最终的束偏移成像结果。终的束偏移成像结果。终的束偏移成像结果。

【技术实现步骤摘要】
一种自适应射线加密式克希霍夫型束偏移地震波成像方法


[0001]本专利技术涉及地震偏移成像
，尤其涉及一种自适应射线加密式克希霍夫型束偏移地震波成像方法。

技术介绍

[0002]随着油气资源的长期持续勘探与开发，浅层油气藏的发现难度持续增加，勘探开发目的层的深度逐渐增加，深层油气资源勘探正得到越来越多的关注，但通常深部储层构造比较复杂且模型较大，给地震偏移成像技术带来了巨大挑战。常规克希霍夫型波束偏移是一种兼具成像精度与计算效率的地震成像方法，但在针对大深度区域成像时，难以保证充足的射线覆盖率，进而影响了这些区域的成像效果。
[0003]吉林大学2017年博士论文公开了《克希霍夫型动态聚焦波束偏移》，介绍了克希霍夫型波束偏移成像方法的实现细节，对其中的射线束传播算子，成像条件以及在不同介质中的应用进行了详细的讨论。并且通过多个数值模型以及实际数据对克希霍夫型波束偏移地震波成像方法进行了验证，计算结果取得了良好的效果。
[0004]201711077893.X公开了一种《一种复杂地形条件下的克希霍夫型波束偏移方法》，将克希霍夫型波束偏移方法扩展到了复杂地形条件下，并且通过采用分段式计算的波束宽度，以提高偏移成像方法在起伏地形模型浅层的成像能力。针对数值模拟数据以及实际数据的成像结果证实了这一方法的良好成像能力。
[0005]201910977970.X公开了一种《一种各向异性地震成像方法》，将克希霍夫型波束偏移方法扩展到了各向异性介质下，通过在偏移成像公式中增加了一个权系数，提高了有效信号对最终偏移成像结果的贡献比例，并且通过已经多个数值模拟数据对这一各向异性地震成像方法进行了验证，计算结果证明了该方法针对各向异性介质具有良好的成像能力。
[0006]通过以上例子可以看出，现有的克希霍夫型波束偏移方法针对常规模型、起伏地形模型以及各向异性介质都具有良好的成像能力，但是这些方法都没有能解决克希霍夫型波束偏移针对大深度区域成像时，难以保证模型深部的射线覆盖率不足问题，从而保证不了地震波成像方法在这些区域的成像效果。

技术实现思路

[0007]本专利技术要解决的技术问题是针对上述现有技术的不足，提供一种自适应射线加密式克希霍夫型束偏移地震波成像方法，通过采用自适应射线加密式射线追踪方法，增强在地震波成像目标模型深部的射线照明，并在这些射线的基础上改进偏移成像公式，以此增强地震波成像方法在地震波成像目标模型深部的成像能力。
[0008]为解决上述技术问题，本专利技术所采取的技术方案是：一种自适应射线加密式克希霍夫型束偏移地震波成像方法，包括以下步骤：
[0009]步骤1：读入地震记录、偏移速度模型以及参数文件；其中参数文件中包含偏移速度模型的横向和纵向的网格点数以及横向和纵向的网格间距、目标成像区域范围、震源的
位置信息、检波器的分布和位置信息、参考频率、最大频率、地震记录每一道的采样点数以及时间采样间隔、单炮对应的接收道以及道间距、中心射线的出射角度范围以及初射角度的采样间隔；
[0010]步骤2：从炮点按照自适应射线加密方式进行射线追踪，并储存射线信息；
[0011]步骤2.1：以炮点为起始点，首先依照中心射线的出射角度范围以及初射角度的采样间隔Δθ，确定炮点初始射线的数量以及初射角度；计算目标成像区域内规则网格点对应的插入新射线的距离阈值F，其计算方式为：
[0012][0013]其中，X表示目标区域内网格点的横坐标，Z表示目标区域内网格点的纵坐标，S
x
表示炮点的横坐标，S
z
表示炮点的纵坐标，dx表示偏移速度模型横向的网格间距，dz表示偏移速度模型纵向的网格间距；
[0014]步骤2.2：从炮点沿着不同初射方向通过龙哥库塔法求解运动学射线追踪方程组进行射线追踪，得到中心射线上离散点信息；所述运动学射线追踪方程组如下公式所示：
[0015][0016]其中，x
i
、p
i
、v
i
分别表示中心射线上第i个离散点的空间位置、慢度和速度值，τ表示地震波走时；
[0017]步骤2.3：在射线追踪过程中如果同一地震波前面上相邻射线上的离散点的距离D大于离散点插入新射线的距离阈值F，在这两个离散点之间插入新的射线，新射线起始点为这两个离散点在空间上的中间位置，新插入的射线起始点的属性信息为这两个离散点各类属性信息的均值，新射线起始点之前的离散点全部记为空；
[0018]步骤3：在炮点完成全部射线追踪后，根据相应的炮点射线信息计算并储存炮点射线对应的射线束信息；
[0019]在炮点完成全部射线追踪后，射线对应的射线束的宽度w如下公式所示：
[0020][0021]其中，V0为目标成像区域内规则网格点速度的平均值，v表示中心射线上离散点的速度值，ray表示中心射线，s表示射线路径；
[0022]步骤4：将地震单炮记录按照窗间隔分为多个以窗中心为核心的数据体，并对每一个数据体进行局部平面波分解；
[0023]步骤5：针对每一个窗中心，依照自适应射线加密方式进行射线追踪并储存射线信息；
[0024]步骤6：在每个窗中心均完成全部射线追踪后，根据相应的窗中心射线信息计算并储存各窗中心射线对应的射线束信息；
[0025]步骤7：在炮点以及每一个窗中心点，依据中心射线与其相邻中心射线的出射角度，计算每条射线束覆盖范围内网格点对应的偏移加权系数；
[0026]所述每条射线束覆盖范围内网格点对应的偏移加权系数的计算方法为：
[0027]设定初射角度为θ的射线束覆盖了目标成像区域内规则网格点P，射线束到达P点走时为τ
P
，从射线初射角度向左寻找同一初射点的中心射线，按照出射角度差由小及大进行搜索，如果中心射线上走时为τ
P
的离散点不为空，结束向左搜索，记录此中心射线的初射角度为θ
left
；同理，从射线初射角度向右寻找同一初射点的中心射线，按照出射角度差由小及大进行搜索，如果中心射线上走时为τ
P
的离散点不为空，结束向右搜索，记录此中心射线的初射角度为θ
right
，此时目标成像区域内规则网格点P针对此射线束的偏移公式加权系数M如下公式所示：
[0028][0029]步骤8：选取炮点和窗中心点的射线束对按照偏移公式进行偏移成像计算，并累加所有射线束对偏移成像结果获得最终单炮地震记录数据的束偏移成像结果，如下公式所示：
[0030][0031]其中，I
s
表示单炮地震记录数据的束偏移成像值，L表示不同的窗中心，p
s
表示炮点s发出射线的慢度值，p
bc
表示窗中心发出射线的慢度值，A表示与目标成像区域内规则网格点到两条中心射线的距离有关的权函数，D
s
表示局部平面波分解结果，p表示局部平面波分解过程中使用的慢度值，M
s
表本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种自适应射线加密式克希霍夫型束偏移地震波成像方法，其特征在于：包括以下步骤：步骤1：读入地震记录、偏移速度模型以及参数文件；步骤2：从炮点按照自适应射线加密方式进行射线追踪，并储存射线信息；步骤3：在炮点完成全部射线追踪后，根据相应的炮点射线信息计算并储存炮点射线对应的射线束信息；步骤4：将地震单炮记录按照窗间隔分为多个以窗中心为核心的数据体，并对每一个数据体进行局部平面波分解；步骤5：针对每一个窗中心，依照自适应射线加密方式进行射线追踪并储存射线信息；步骤6：在每个窗中心均完成全部射线追踪后，根据相应的窗中心射线信息计算并储存各窗中心射线对应的射线束信息；步骤7：在炮点以及每一个窗中心点，依据中心射线与其相邻中心射线的出射角度，计算每条射线束覆盖范围内网格点对应的偏移加权系数；步骤8：选取炮点和窗中心点的射线束对按照偏移公式进行偏移成像计算，并累加所有射线束对偏移成像结果获得最终单炮地震记录数据的束偏移成像结果。2.根据权利要求1所述的一种自适应射线加密式克希霍夫型束偏移地震波成像方法，其特征在于：步骤1所述参数文件中包括偏移速度模型的横向和纵向的网格点数以及横向和纵向的网格间距、目标成像区域范围、震源的位置信息、检波器的分布和位置信息、参考频率、最大频率、地震记录每一道的采样点数以及时间采样间隔、单炮对应的接收道以及道间距、中心射线的出射角度范围以及初射角度的采样间隔。3.根据权利要求2所述的一种自适应射线加密式克希霍夫型束偏移地震波成像方法，其特征在于：所述步骤2的具体方法为：步骤2.1：以炮点为起始点，首先依照中心射线的出射角度范围以及初射角度的采样间隔Δθ，确定炮点初始射线的数量以及初射角度；计算目标成像区域内规则网格点对应的插入新射线的距离阈值F，其计算方式为：其中，X表示目标区域内网格点的横坐标，Z表示目标区域内网格点的纵坐标，S
x
表示炮点的横坐标，S
z
表示炮点的纵坐标，dx表示偏移速度模型横向的网格间距，dz表示偏移速度模型纵向的网格间距；步骤2.2：从炮点沿着不同初射方向通过龙哥库塔法求解运动学射线追踪方程组进行射线追踪，得到中心射线上离散点信息；所述运动学射线追踪方程组如下公式所示：其中，x
i
、p
i
、v
i
分别表示中心射线上第i个离散点的空间位置、慢度和速度值，τ表示地震波走时；
步骤2...

【专利技术属性】
技术研发人员：孙辉，李猛，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：