基于射线密度生成层析反演模型速度界面的方法和装置制造方法及图纸

本发明专利技术提供一种基于射线密度生成层析反演模型速度界面的方法和装置，对地震初至数据进行层析反演得到三维测网的地下速度深度空间模型，该三维测网包括多个网格；根据该地下速度深度空间模型得到每个网格的深度数据及射线密度数据；确定每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度；根据每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度生成层析反演模型速度界面，即通过采用射线密度生成层析反演模型速度界面，不需要大量的表层深度数据，无需野外开展大量的微测井调查，降低采集作业成本，能够更精确的刻画界面特征，所生成的速度界面与实际界面深度误差小。所生成的速度界面与实际界面深度误差小。所生成的速度界面与实际界面深度误差小。

【技术实现步骤摘要】
基于射线密度生成层析反演模型速度界面的方法和装置


[0001]本专利技术涉及地球物理勘探
，尤其涉及一种基于射线密度生成层析反演模型速度界面的方法和装置。

技术介绍

[0002]层析反演是目前复杂地区基于初至反演表层模型的一种常用方法，在定义和解释模型时，目前主要采用提取常速度界面、折射反演界面、自定义的空变速度界面以及采用深度标定的界面。上述生成层析反演模型界面的方法在许多复杂地区的应用中受到限制，比如，基于速度生成界面的方法，其构成的表层模型精度不足，其所生成的高速顶界面与实际界面深度误差较大；基于深度标定的界面生成方法需要大量的表层深度数据，需要野外开展大量的微测井调查，增加采集作业成本。

技术实现思路

[0003]有鉴于此，本专利技术提供一种基于射线密度生成层析反演模型速度界面的方法和装置、电子设备以及计算机可读存储介质，通过采用射线密度生成层析反演模型速度界面，不需要大量的表层深度数据，无需野外开展大量的微测井调查，降低采集作业成本，另一方面，由于利用了射线密度，能够更精确的刻画界面特征，所生成的速度界面与实际界面深度误差小，确保解释模型更接近实际表层模型，从而更好地满足复杂地区表层建模和静校正的需要。
[0004]为了实现上述目的，本专利技术采用如下技术方案：
[0005]第一方面，提供一种基于射线密度生成层析反演模型速度界面的方法，包括：
[0006]对地震初至数据进行层析反演得到三维测网的地下速度深度空间模型，该三维测网包括多个网格；
[0007]根据该地下速度深度空间模型得到每个网格的深度数据及射线密度数据；
[0008]确定每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度；
[0009]根据每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度生成层析反演模型速度界面。
[0010]进一步地，该确定每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度包括：
[0011]对该网格的预定深度范围内的深度数据及射线密度数据进行二次曲线拟合，得到深度-射线密度的二次函数；
[0012]计算该二次函数的极大值，得到该网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度。
[0013]进一步地，该二次函数为：
[0014]H＝as2+bs+c，
[0015]其中，H表示模型深度，s表示射线密度；a、b、c分别是二次函数的系数，
[0016][0017][0018][0019]其中，h
i
是网格i处的深度，i的取值范围为i＝1,2,
…
m，m是网格在垂向上的数目，s
i
是网格i处的射线密度。
[0020]进一步地，计算该二次函数的极大值，得到该网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度，采用如下公式：
[0021][0022][0023]其中，s
max
表示该网格的预定深度范围内的射线密度最大值；h表示该射线密度最大值对应的深度。
[0024]进一步地，该确定每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度包括：
[0025]遍历该网格的预定深度范围内的射线密度得到该射线密度最大值；
[0026]根据该射线密度最大值获取其对应的深度。
[0027]进一步地，基于射线密度生成层析反演模型速度界面的方法还包括：
[0028]获取地震初至数据。
[0029]进一步地，该获取地震初至数据，包括：
[0030]获取SPS文件；
[0031]对该SPS文件中的全部生产炮拾取预设偏移距范围内的地震初至时间得到该地震初至数据。
[0032]进一步地，该确定每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度之前，还包括：
[0033]根据预设的有效射线深度范围筛选每个网格的深度数据及射线密度数据，得到有
效射线深度范围内的网格深度数据和射线密度数据。
[0034]进一步地，该根据每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度生成层析反演模型速度界面，包括：
[0035]对所有网格的预定深度范围内的射线密度最大值及与该射线密度最大值进行后处理，得到后处理数据；
[0036]根据预设半径对该后处理数据进行平滑得到所有物理点位置处的深度数据；
[0037]根据该深度数据生成层析反演模型速度界面。
[0038]进一步地，该后处理包括：异常值删除、内插和/或边界效应处理。
[0039]进一步地，该对地震初至数据进行层析反演得到三维测网的地下速度深度空间模型，包括：
[0040]定义三维测网网格参数、模型深度范围、初始模型以及偏移距使用范围；
[0041]根据该三维测网网格参数、该模型深度范围、该初始模型以及该偏移距使用范围对该地震初至数据进行层析反演得到三维测网的地下速度深度空间模型。
[0042]第二方面，提供一种基于射线密度生成层析反演模型速度界面的装置，包括：
[0043]层析反演模块，对地震初至数据进行层析反演得到三维测网的地下速度深度空间模型，该三维测网包括多个网格；
[0044]网格数据获取模块，根据该地下速度深度空间模型得到每个网格的深度数据及射线密度数据；
[0045]极值获取模块，确定每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度；
[0046]速度界面生成模块，根据每个网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度生成层析反演模型速度界面。
[0047]进一步地，该极值获取模块包括：
[0048]拟合单元，对该网格的预定深度范围内的深度数据及射线密度数据进行二次曲线拟合，得到深度-射线密度的二次函数；
[0049]计算单元，计算该二次函数的极大值，得到该网格的预定深度范围内的射线密度最大值及与该射线密度最大值对应的深度。
[0050]进一步地，该极值获取模块包括：
[0051]遍历单元，遍历该网格的预定深度范围内的射线密度得到该射线密度最大值；
[0052]深度获取单元，根据该射线密度最大值获取其对应的深度。
[0053]进一步地，基于射线密度生成层析反演模型速度界面的装置还包括：
[0054]初至数据获取模块，获取地震初至数据。
[0055]进一步地，该初至数据获取模块包括：
[0056]SPS文件获取单元，获取SPS文件；
[0057]数据拾取单元，对该SPS文件中的全部生产炮拾取预设偏移距范围内的地震初至时间得到该地震初至数据。
[0058]进一步地，基于射线密度生成层析反演模型速度界面的装置还包括：
[0059]数据筛选模块，根据预设的有效射线深度范围筛本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于射线密度生成层析反演模型速度界面的方法，其特征在于，包括：对地震初至数据进行层析反演得到三维测网的地下速度深度空间模型，所述三维测网包括多个网格；根据所述地下速度深度空间模型得到每个网格的深度数据及射线密度数据；确定每个网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度；根据每个网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度生成层析反演模型速度界面。2.根据权利要求1所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述确定每个网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度包括：对所述网格的预定深度范围内的深度数据及射线密度数据进行二次曲线拟合，得到深度-射线密度的二次函数；计算所述二次函数的极大值，得到所述网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度。3.根据权利要求2所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述二次函数为：H＝as2+bs+c，其中，H表示模型深度，s表示射线密度；a、b、c分别是二次函数的系数，其中，H表示模型深度，s表示射线密度；a、b、c分别是二次函数的系数，其中，H表示模型深度，s表示射线密度；a、b、c分别是二次函数的系数，其中，h
i
是网格i处的深度，i的取值范围为i＝1,2,
…
m，m是网格在垂向上的数目，s
i
是网格i处的射线密度。4.根据权利要求3所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，计算所述二次函数的极大值，得到所述网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度，采用如下公式：
其中，s
max
表示所述网格的预定深度范围内的射线密度最大值；h表示所述射线密度最大值对应的深度。5.根据权利要求1所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述确定每个网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度包括：遍历所述网格的预定深度范围内的射线密度得到所述射线密度最大值；根据所述射线密度最大值获取其对应的深度。6.根据权利要求1所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，还包括：获取地震初至数据。7.根据权利要求6所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述获取地震初至数据，包括：获取SPS文件；对所述SPS文件中的全部生产炮拾取预设偏移距范围内的地震初至时间得到所述地震初至数据。8.根据权利要求1所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述确定每个网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度之前，还包括：根据预设的有效射线深度范围筛选每个网格的深度数据及射线密度数据，得到有效射线深度范围内的网格深度数据和射线密度数据。9.根据权利要求1所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述根据每个网格的预定深度范围内的射线密度最大值及与所述射线密度最大值对应的深度生成层析反演模型速度界面，包括：对所有网格的预定深度范围内的射线密度最大值及与所述射线密度最大值进行后处理，得到后处理数据；根据预设半径对所述后处理数据进行平滑得到所有物理点位置处的深度数据；根据所述深度数据生成层析反演模型速度界面。10.根据权利要求9所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述后处理包括：异常值删除、内插和/或边界效应处理。11.根据权利要求1所述的基于射线密度生成层析反演模型速度界面的方法，其特征在于，所述对地震初至数...

【专利技术属性】
技术研发人员：冯发全，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：