获取风化层的地质参数的方法

摘要

本发明提供一种获取风化层的地质参数的方法，所述方法包括：(a)获取目标区域的标准的初至数据；(b)确定最小偏移距；(c)确定每个点的风化层速度值；(d)计算选择的m个点中的任意两点之间的风化层速度值的变差函数值；(e)计算选择的n个点中的任意两点之间的微测井资料的速度值的变差函数值；(f)计算m个点中的任意一点与n个点中的任意一点之间的关于速度的交叉编译函数；(g)获得与m个点中的每个点的风化层速度值对应的m个加权系数和与n个点中的每个点的微测井数据的速度值对应的n个加权系数；(h)计算风化层速度的估计值。采用本发明获取风化层的地质参数的方法可得到更符合地质分布规律的浅层速度结构。

主权项

一种获取风化层的地质参数的方法，所述方法包括：(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个点中的每个点的微测井资料的速度值，计算风化层速度的估计值，其中，步骤(c)还包括：对选择的m个点中的每个点在最小偏移距范围内对应的初至数据进行拟合，确定每个点的炮点延迟时的初始值，其中，所述方法还包括：(o)确定最大偏移距；(p)从获取的初至数据中选择m个点中的每个点在最小偏移距至最大偏移距范围内的初至数据，对选择的m个点中的每个点在最小偏移距至最大偏移距范围内对应的初至数据进行拟合，确定每个点的高速层速度的初始值；(q)根据确定的m个点中的每个点的高速层速度的初始值和m个点中的每个点的炮点延迟时的初始值，利用折射初至时间公式计算m个点中的每个点的检波点延迟时、炮点延迟时和高速层速度，以根据m个点中的每个点的风化层速度值、m个点中的每个点的高速层速度、m个点中的每个点的炮点延迟时或者检波点延迟时，计算m个点中的每个点的风化层厚度值，其中，步骤(d)还包括：根据确定的每个点的风化层厚度值，计算选择的m个点中的任意两点之间的风化层厚度值的变差函数值；步骤(e)还包括：从所述目标区域的微测井资料中提取n个点中的每个点的微测井资料的厚度值，并根据提取的每个点的微测井资料的厚度值，计算选择的n个点中的任意两点之间的微测井资料的厚度值的变差函数值；步骤(f)还包括：根据确定的m个点中的每个点的风化层厚度值和n个点中的每个点的微测井资料的厚度值，计算m个点中的任意一点与n个点中的任意一点之间的关于厚度的交叉编译函数；步骤(g)还包括：根据步骤(d)～步骤(f)的结果获得与m个点中的每个点的风化层厚度值对应的m个加权系数和与n个点中的每个点的微测井数据的厚度值对应的n个加权系数；步骤(h)还包括：根据获得的与m个点中的每个点的风化层厚度值对应的m个加权系数、与n个点中的每个点的微测井数据的厚度值对应的n个加权系数、确定的m个点中的每个点的风化层厚度值和n个点中的每个点的微测井资料的厚度值，计算风化层厚度的估计值，其中，步骤(q)包括：(q1)根据m个点中的每个点的高速层速度、炮点延迟时、初至数据中每个点对应的炮检对从激发到接收产生的时间、每个点对应的炮点与检波点之间的距离，利用折射初至时间公式计算m个点中的每个点的检波点延迟时，其中，将m个点中的每个点的高速层速度的初始值和m个点中的每个点的炮点延迟时的初始值作为初值代入；(q2)根据m个点中的每个点的高速层速度、步骤(q1)得到的m个点中的每个点的检波点延迟时、初至数据中每个点对应的炮检对从激发到接收产生的时间、每个点对应的炮点与检波点之间的距离，利用折射初至时间公式计算m个点中的每个点的炮点延迟时；(q3)根据步骤(q1)得到的m个点中的每个点的检波点延迟时、步骤(q2)得到的m个点中的每个点的炮点延迟时、初至数据中每个点对应的炮检对从激发到接收产生的时间、每个点对应的炮点与检波点之间的距离，利用折射初至时间公式计算m个点中的每个点的高速层速度；(q4)根据步骤(q1)得到的m个点中的每个点的检波点延迟时、步骤(q2)得到的m个点中的每个点的炮点延迟时、步骤(q3)得到的m个点中的每个点的高速层速度、每个点对应的炮点与检波点之间的距离，利用折射初至时间公式计算对应的初至数据中每个点对应的炮检对从激发到接收产生的时间；(q5)当所述时间小于等于设定的时间范围时，将步骤(q1)～步骤(q3)的结果作为m个点中的每个点的炮点延迟时、检波点延迟时和高速层速度，当所述时间大于设定的时间范围时，则返回执行步骤(q1)，其中，所述方法还包括：(m)根据风化层速度的估计值和风化层厚度的估计值，计算折射静校正量，其中，利用下面的公式计算折射静校正量其中，为在待求位置0处对应的炮点或检波点折射静校正量，τ为井深或检波器初至时间，k为低降速层的层数，1≤k≤f，f为大于等于1的自然数，D₀为给定的厚度，Z^*(0)为风化层速度在待求位置0处的估计值，为风化层厚度在待求位置0处的估计值，H_d为基准面高程，H_g为高速层顶界面高程，V_R为基准面校正速度。