获取风化层的地质参数的方法
【专利摘要】本发明提供一种获取风化层的地质参数的方法，所述方法包括：(a)获取目标区域的标准的初至数据；(b)确定最小偏移距；(c)确定每个点的风化层速度值；(d)计算选择的m个点中的任意两点之间的风化层速度值的变差函数值；(e)计算选择的n个点中的任意两点之间的微测井资料的速度值的变差函数值；(f)计算m个点中的任意一点与n个点中的任意一点之间的关于速度的交叉编译函数；(g)获得与m个点中的每个点的风化层速度值对应的m个加权系数和与n个点中的每个点的微测井数据的速度值对应的n个加权系数；(h)计算风化层速度的估计值。采用本发明获取风化层的地质参数的方法可得到更符合地质分布规律的浅层速度结构。
【专利说明】获取风化层的地质参数的方法
【技术领域】
[0001] 本发明总体来说涉及石油地震资料分析处理领域，更具体地讲，涉及一种基于协 同克里金技术融合微测井资料的获取风化层的地质参数的方法。
【背景技术】
[0002] 在地震勘探中，近地表的低降速带（S卩，风化层）的速度和厚度不均匀变化会引起 地震波传播旅行时的不均匀延迟，致使实际观测得到的反射波到达时间发生畸变，极大地 影响着地震剖面的成像和构造形态。静校正就是消除地表低降速带对勘探资料的影响，对 于复杂地表条件下地震资料的处理，求取正确的静校正量是非常关键和必要的。如果不能 很好的解决严重的静校正问题，就会在地震剖面上产生假象或不成像，就会给地质解释带 来严重困难，甚至错误的地质构造认识。
[0003] 目前，在解决静校正问题过程中，对于存在稳定折射层的区域，一般采用折射静校 正方法（例如，基于迭代法的折射静校正方法和广义互换静校正方法）进行计算。但这两 类方法在使用时，对于风化层速度，一般参照工作经验，给定常速度值进行计算。然而，很多 地区的风化层速度横向变化剧烈，实际风化层速度与前述参照工作经验给定的常速度值存 在较大差距。
[0004] 此外，有的方法还会对实际风化层速度进行拟合计算或者结合微测井资料进行速 度反演，其效果优于给定常速度值计算的结果。但采用速度拟合的方法仍存在地震剖面的 成像横向不连续现象，需要进一步处理；而进行速度反演，则计算时间长，微测井资料的好 坏将直接影响速度反演的结果，不确定程度高，从而影响处理后叠加剖面成像效果，使得构 造形态不准确，信噪比大大降低，不能达到资料精细处理的目的。
【发明内容】
[0005] 本发明的目的是提出一种基于协同克里金技术融合微测井资料的获取风化层的 地质参数的方法，通过变异函数这一数学工具，分别刻画出微测井资料的速度值、厚度值和 由初至数据拟合出的风化层速度、厚度的变化趋势，再运用交叉编译函数得到两类数据的 整体变化趋势，并加以它们的相关系数约束，最终得到符合地质分布规律的浅层速度结构。
[0006] 本发明的一方面提供一种获取风化层的地质参数的方法，所述方法包括：(a)获 取目标区域的标准的初至数据；(b)确定最小偏移距；（c)在所述目标区域内选择m个点， 从获取的初至数据中选择m个点中的每个点在最小偏移距范围内的初至数据，并对选择的 m个点中的每个点在最小偏移距范围内对应的初至数据进行拟合，确定每个点的风化层速 度值，m为大于等于1的自然数；（d)根据确定的每个点的风化层速度值，计算选择的m个 点中的任意两点之间的风化层速度值的变差函数值；(e)在所述目标区域内选择η个点，从 所述目标区域的微测井资料中提取η个点中的每个点的微测井资料的速度值，并根据提取 的每个点的微测井资料的速度值，计算选择的η个点中的任意两点之间的微测井资料的速 度值的变差函数值，η为大于等于1的自然数；（f)根据确定的m个点中的每个点的风化层 速度值和η个点中的每个点的微测井资料的速度值，计算m个点中的任意一点与η个点中 的任意一点之间的关于速度的交叉编译函数；（g)根据步骤（d)?步骤（f)的结果获得与 m个点中的每个点的风化层速度值对应的m个加权系数和与η个点中的每个点的微测井数 据的速度值对应的η个加权系数；(h)根据获得的与m个点中的每个点的风化层速度值对 应的m个加权系数、与η个点中的每个点的微测井数据的速度值对应的η个加权系数、确定 的m个点中的每个点的风化层速度值和η个点中的每个点的微测井资料的速度值，计算风 化层速度的估计值。
[0007] 可选地，步骤（g)可包括：根据步骤（d)?步骤（f)的结果利用协同克里金方程组 获得与m个点中的每个点的风化层速度值对应的m个加权系数和与η个点中的每个点的微 测井数据的速度值对应的η个加权系数。
[0008] 可选地，步骤（g)可包括：利用下面的协同克里金方程组计算与m个点中的每个点 的风化层速度值对应的m个加权系数和与η个点中的每个点的微测井数据的速度值对应的 η个加权系数，
【权利要求】
1. 一种获取风化层的地质参数的方法，所述方法包括： (a) 获取目标区域的标准的初至数据； (b) 确定最小偏移距； (c) 在所述目标区域内选择m个点，从获取的初至数据中选择m个点中的每个点在最小 偏移距范围内的初至数据，并对选择的m个点中的每个点在最小偏移距范围内对应的初至 数据进行拟合，确定每个点的风化层速度值，m为大于等于1的自然数； (d) 根据确定的每个点的风化层速度值，计算选择的m个点中的任意两点之间的风化 层速度值的变差函数值； (e) 在所述目标区域内选择η个点，从所述目标区域的微测井资料中提取η个点中的 每个点的微测井资料的速度值，并根据提取的每个点的微测井资料的速度值，计算选择的η 个点中的任意两点之间的微测井资料的速度值的变差函数值，η为大于等于1的自然数； (f) 根据确定的m个点中的每个点的风化层速度值和η个点中的每个点的微测井资料 的速度值，计算m个点中的任意一点与η个点中的任意一点之间的关于速度的交叉编译函 数； (g) 根据步骤（d)?步骤（f)的结果获得与m个点中的每个点的风化层速度值对应的 m个加权系数和与η个点中的每个点的微测井数据的速度值对应的η个加权系数； (h) 根据获得的与m个点中的每个点的风化层速度值对应的m个加权系数、与η个点中 的每个点的微测井数据的速度值对应的η个加权系数、确定的m个点中的每个点的风化层 速度值和η个点中的每个点的微测井资料的速度值，计算风化层速度的估计值。
2. 根据权利要求1所述的方法，其中，步骤（g)包括：根据步骤（d)?步骤（f)的结果 利用协同克里金方程组获得与m个点中的每个点的风化层速度值对应的m个加权系数和与 η个点中的每个点的微测井数据的速度值对应的η个加权系数。
3. 根据权利要求2所述的方法，其中，步骤（g)包括： 利用下面的协同克里金方程组计算与m个点中的每个点的风化层速度值对应的m个加 权系数和与η个点中的每个点的微测井数据的速度值对应的η个加权系数，
其中,Cx (Xi,χρ)为η个点中的第i点的坐标位置Xi与η个点中的第ρ点的坐标位置χρ 对应的微测井资料的速度值的变差函数值，Cx (0, Xi)为待求位置O与η个点中的第i点的坐 标位置Xi对应的微测井资料的速度值的变差函数值，Cy (y,，0为m个点中的第q点的坐标 位置!与m个点中的第j点的坐标位置yj对应的风化层速度值的变差函数值,Cy (0,yj为 待求位置0与m个点中的第j点的坐标位置t对应的风化层速度值的变差函数值，CJbyj) 和Cc (yj，Xi)为关于η个点中的第i点的坐标位置Xi对应的微测井资料的速度值与m个点 中的第j点的坐标位置&对应的风化层速度值的交叉编译函数，为协同克里金方程组中 与η个点中的第i点的微测井数据的速度值对应的加权系数，/>\,为协同克里金方程组中与 m个点中的第j点的风化层速度值对应的加权系数，μχ、μy为计算方程组的解与原始值的 误差的参数。
4. 根据权利要求1所述的方法，其中，步骤（h)包括： 利用下面的公式计算风化层速度的估计值Z(O)，
其中，Zi(O)为风化层速度在待求位置0处的估计值，Zx(Xi)为η个点中的第i点的坐 标位置Xi对应的微测井资料的速度值，A为与η个点中的第i点的微测井数据的速度值对 应的加权系数，Zy (yp为m个点中的第j点的坐标位置&对应的风化层速度值，Py,.为与m 个点中的第j点的风化层速度值对应的加权系数。
5. 根据权利要求1所述的方法，其中，步骤（c)还包括：对选择的m个点中的每个点在 最小偏移距范围内对应的初至数据进行拟合，确定每个点的炮点延迟时的初始值， 其中，所述方法还包括： (〇)确定最大偏移距； (P)从获取的初至数据中选择m个点中的每个点在最小偏移距至最大偏移距范围内的 初至数据，对选择的m个点中的每个点在最小偏移距至最大偏移距范围内对应的初至数据 进行拟合，确定每个点的高速层速度的初始值； (q)根据确定的m个点中的每个点的高速层速度的初始值和m个点中的每个点的炮点 延迟时的初始值，利用折射初至时间公式计算m个点中的每个点的检波点延迟时、炮点延 迟时和高速层速度。
6. 根据权利要求5所述的方法，其中，步骤（q)包括： (ql)根据m个点中的每个点的高速层速度、炮点延迟时、初至数据中每个点对应的炮 检对从激发到接收产生的时间、每个点对应的炮点与检波点之间的距离，利用折射初至时 间公式计算m个点中的每个点的检波点延迟时，其中，将m个点中的每个点的高速层速度的 初始值和m个点中的每个点的炮点延迟时的初始值作为初值代入； (q2)根据m个点中的每个点的高速层速度、步骤（ql)得到的m个点中的每个点的检波 点延迟时、初至数据中每个点对应的炮检对从激发到接收产生的时间、每个点对应的炮点 与检波点之间的距离，利用折射初至时间公式计算m个点中的每个点的炮点延迟时； (q3)根据步骤（ql)得到的m个点中的每个点的检波点延迟时、步骤（q2)得到的m个 点中的每个点的炮点延迟时、初至数据中每个点对应的炮检对从激发到接收产生的时间、 每个点对应的炮点与检波点之间的距离，利用折射初至时间公式计算m个点中的每个点的 高速层速度； (q4)根据步骤（ql)得到的m个点中的每个点的检波点延迟时、步骤（q2)得到的m个 点中的每个点的炮点延迟时、步骤（q3)得到的m个点中的每个点的高速层速度、每个点对 应的炮点与检波点之间的距离，利用折射初至时间公式计算对应的初至数据中每个点对应 的炮检对从激发到接收产生的时间； (q5)当所述时间小于等于设定的时间范围时，将步骤（ql)?步骤（q3)的结果作为m个点中的每个点的炮点延迟时、检波点延迟时和高速层速度，当所述时间大于设定的时间 范围时，则返回执行步骤（ql)。
7. 根据权利要求6所述的方法，其中，折射初至时间公式为，
其中，为初至数据中m个点中的第j点对应的炮检对从激发到接收产生的时间，TA(yj)为m个点中的第j点的坐标位置y』对应的炮点延迟时，TB(yj)为m个点中的第j点 的坐标位置yj对应的检波点延迟时，/iZ?/为m个点中的第j点对应的炮点与检波点之间的 距离，S(yj)为m个点中的第j点的坐标位置yj对应的高速层速度。
8. 根据权利要求5所述的方法，所述方法还包括： (k)根据m个点中的每个点的风化层速度值、m个点中的每个点的高速层速度、m个点 中的每个点的炮点延迟时或者检波点延迟时，计算m个点中的每个点的风化层厚度值； 其中，步骤（d)还包括：根据确定的每个点的风化层厚度值，计算选择的m个点中的任 意两点之间的风化层厚度值的变差函数值； 步骤（e)还包括：从所述目标区域的微测井资料中提取η个点中的每个点的微测井资 料的厚度值，并根据提取的每个点的微测井资料的厚度值，计算选择的η个点中的任意两 点之间的微测井资料的厚度值的变差函数值； 步骤（f)还包括：根据确定的m个点中的每个点的风化层厚度值和η个点中的每个点 的微测井资料的厚度值，计算m个点中的任意一点与η个点中的任意一点之间的关于厚度 的交叉编译函数； 步骤（g)还包括：根据步骤（d)?步骤（f)的结果获得与m个点中的每个点的风化层 厚度值对应的m个加权系数和与η个点中的每个点的微测井数据的厚度值对应的η个加权 系数； 步骤（h)还包括：根据获得的与m个点中的每个点的风化层厚度值对应的m个加权系 数、与η个点中的每个点的微测井数据的厚度值对应的η个加权系数、确定的m个点中的每 个点的风化层厚度值和η个点中的每个点的微测井资料的厚度值，计算风化层厚度的估计 值。
9. 根据权利要求8所述的方法，其中，步骤（k)包括： 利用下面的公式计算m个点中的每个点的风化层厚度值，
其中，Dy (yp为m个点中的第j点的坐标位置t对应的风化层厚度值，T(yp为m个点 中的第j点的坐标位置yj对应的炮点延迟时或检波点延迟时，Zy (y」）为m个点中的第j点 的坐标位置y」对应的风化层速度值，S(y」）为m个点中的第j点的坐标位置y」对应的高速 层速度。
10.根据权利要求8所述的方法，所述方法还包括：(m)根据风化层速度的估计值和风 化层厚度的估计值，计算折射静校正量。
【文档编号】G01V1/36GK104316961SQ201410612194
【公开日】2015年1月28日 申请日期:2014年11月4日 优先权日:2014年11月4日
【发明者】张亨, 陈爱萍, 何光明, 刘奇琳, 曹中林, 周强, 漆韬, 李若禹 申请人:中国石油集团川庆钻探工程有限公司地球物理勘探公司