[发明专利]地震建模中的浅中深层速度融合方法

权利要求书

1.地震建模中的浅中深层速度融合方法，其特征在于：包括下列步骤：

步骤1：取近地表速度模型；获取中、深层速度模型；

步骤2：选取适当的窗口，将所述近地表速度模型和所述中、深层速度模型融合为一个统一的浅中深层速度模型；

融合时，对于近地表和中、深层速度函数在深度域上出现重叠区域时，通过加权函数合并融合浅中深层速度模型,对于两速度函数在深度域上出现空缺区域时，通过插值补缺，建立统一的全速度函数；

步骤3：输出该统一的浅中深层速度模型；

所述步骤2的具体步骤包括：

第一步：融合数据准备，主要包括近地表层析速度模型和叠前时间偏移获取的中、深层速度模型d=0对应的物理界面、各道大地坐标信息和采样深度间隔的统一，以及两个速度模型重叠或者缺失的区域确定；

第二步：融合起始和终止界面选取

如果地表相对光滑，直接从地表开始进行融合计算；如果地表局部变化剧烈，对剧烈起伏地表进行适度平滑或选取靠近地表的某一平缓界面开始融合计算；

如果地下出现速度随深度变化较大的强反射面、低速带底界或者折射面，选取这一界面为融合基准面，分别向上、向下开窗口进行融合计算；

第三步：针对融合计算的深度区域，选取适当的窗口大小将两速度模型相应深度的数据分段，窗口大小根据每道进行融合的点数确定；

第四步：首先依据近地表速度模型修改中、深层速度模型值，计算每个窗口下两组速度模型数据的均值及其均值比；窗口与窗口之间数据可以交叉，也可以不交叉；

第五步：根据各窗口下计算的均值比与1的相对大小关系，将均值比分为大于1，小于1和等于1三组数据组；

选出均值比大于1数据组中的最小值和均值比小于1的数据组中的中的最大值；

将上述最小值和最大值作为修改中、深层速度模型的修改均值比；

第六步：按深度采样点比较每一深度处近地表速度模型值和中、深层速度模型值的相对大小，根据大小关系利用第五步中选择的修改均值比修改中、深层速度模型值；

第七步：将近地表速度模型和修改后的中、深层速度模型通过加权得到靠近浅层速度模型的融合速度模型；

第八步：依据中、深层速度模型值修改近地表速度模型值，得到靠近中、深层速度模型的融合速度模型，具体步骤为：

计算每个窗口下两组速度模型数据的均值及其均值比；窗口与窗口之间数据可以交叉，也可以不交叉；

根据各窗口下计算的均值比与1的相对大小关系，将均值比分为大于1，小于1和等于1三组数据组；

选出均值比大于1数据组中的最小值和均值比小于1的数据组中的中的最大值；

将上述最小值和最大值作为修改近地表速度模型的修改均值比；

按深度采样点比较每一深度处近地表速度模型值和中、深层速度模型值的相对大小，根据大小关系利用选择的修改均值比修改近地表速度模型；

第九步：将第七步和第八步得到的两个融合速度模型加权输出最终的浅中深层融合速度模型；

第十步：将更新后的融合速度模型选择适当的平滑因子，进行速度模型平滑，输出可用于叠前深度偏移的初始全速度模型。

2.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：在步骤1中，所述近地表速度模型通过基于初至波的近地表层析反演技术获取。

3.如权利要求2所述的地震建模中的浅中深层速度融合方法，其特征在于：在步骤1中，所述取中、深层速度模型通过基于叠前时间偏移获得。

4.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：第七步所述加权为反距离加权或sin函数加权。

5.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：当所述第二步存在融合基准面时，第七步所述加权为斜坡加权。

6.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：第八步所述修改近地表速度模型值是指：计算中、深层速度模型值与近地表速度模型值的均值比，并利用选择的修改均值比修改近地表速度模型。

7.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：第九步所述加权为反距离加权或sin函数加权。

8.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：当所述第二步存在融合基准面时，第九步所述加权为斜坡加权。

9.如权利要求1所述的地震建模中的浅中深层速度融合方法，其特征在于：在步骤3中还包括将建立的浅中深层速度模型进行平滑的步骤。