[发明专利]获取风化层的地质参数的方法

说明书

技术领域

本发明总体来说涉及石油地震资料分析处理领域，更具体地讲，涉及一种基于协同克里金技术融合微测井资料的获取风化层的地质参数的方法。

背景技术

在地震勘探中，近地表的低降速带(即，风化层)的速度和厚度不均匀变化会引起地震波传播旅行时的不均匀延迟，致使实际观测得到的反射波到达时间发生畸变，极大地影响着地震剖面的成像和构造形态。静校正就是消除地表低降速带对勘探资料的影响，对于复杂地表条件下地震资料的处理，求取正确的静校正量是非常关键和必要的。如果不能很好的解决严重的静校正问题，就会在地震剖面上产生假象或不成像，就会给地质解释带来严重困难，甚至错误的地质构造认识。

目前，在解决静校正问题过程中，对于存在稳定折射层的区域，一般采用折射静校正方法(例如，基于迭代法的折射静校正方法和广义互换静校正方法)进行计算。但这两类方法在使用时，对于风化层速度，一般参照工作经验，给定常速度值进行计算。然而，很多地区的风化层速度横向变化剧烈，实际风化层速度与前述参照工作经验给定的常速度值存在较大差距。

此外，有的方法还会对实际风化层速度进行拟合计算或者结合微测井资料进行速度反演，其效果优于给定常速度值计算的结果。但采用速度拟合的方法仍存在地震剖面的成像横向不连续现象，需要进一步处理；而进行速度反演，则计算时间长，微测井资料的好坏将直接影响速度反演的结果，不确定程度高，从而影响处理后叠加剖面成像效果，使得构造形态不准确，信噪比大大降低，不能达到资料精细处理的目的。

发明内容

本发明的目的是提出一种基于协同克里金技术融合微测井资料的获取风化层的地质参数的方法，通过变异函数这一数学工具，分别刻画出微测井资料的速度值、厚度值和由初至数据拟合出的风化层速度、厚度的变化趋势，再运用交叉编译函数得到两类数据的整体变化趋势，并加以它们的相关系数约束，最终得到符合地质分布规律的浅层速度结构。

本发明的一方面提供一种获取风化层的地质参数的方法，所述方法包括：(a)获取目标区域的标准的初至数据；(b)确定最小偏移距；(c)在所述目标区域内选择m个点，从获取的初至数据中选择m个点中的每个点在最小偏移距范围内的初至数据，并对选择的m个点中的每个点在最小偏移距范围内对应的初至数据进行拟合，确定每个点的风化层速度值，m为大于等于1的自然数；(d)根据确定的每个点的风化层速度值，计算选择的m个点中的任意两点之间的风化层速度值的变差函数值；(e)在所述目标区域内选择n个点，从所述目标区域的微测井资料中提取n个点中的每个点的微测井资料的速度值，并根据提取的每个点的微测井资料的速度值，计算选择的n个点中的任意两点之间的微测井资料的速度值的变差函数值，n为大于等于1的自然数；(f)根据确定的m个点中的每个点的风化层速度值和n个点中的每个点的微测井资料的速度值，计算m个点中的任意一点与n个点中的任意一点之间的关于速度的交叉编译函数；(g)根据步骤(d)～步骤(f)的结果获得与m个点中的每个点的风化层速度值对应的m个加权系数和与n个点中的每个点的微测井数据的速度值对应的n个加权系数；(h)根据获得的与m个点中的每个点的风化层速度值对应的m个加权系数、与n个点中的每个点的微测井数据的速度值对应的n个加权系数、确定的m个点中的每个点的风化层速度值和n个点中的每个点的微测井资料的速度值，计算风化层速度的估计值。

可选地，步骤(g)可包括：根据步骤(d)～步骤(f)的结果利用协同克里金方程组获得与m个点中的每个点的风化层速度值对应的m个加权系数和与n个点中的每个点的微测井数据的速度值对应的n个加权系数。

可选地，步骤(g)可包括：利用下面的协同克里金方程组计算与m个点中的每个点的风化层速度值对应的m个加权系数和与n个点中的每个点的微测井数据的速度值对应的n个加权系数，