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

权利要求书

1.一种获取风化层的地质参数的方法，所述方法包括：

(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个点中的每个点的微测井资料的速度值，计算风化层速度的估计值。

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

3.根据权利要求2所述的方法，其中，步骤(g)包括：

利用下面的协同克里金方程组计算与m个点中的每个点的风化层速度值对应的m个加权系数和与n个点中的每个点的微测井数据的速度值对应的n个加权系数，

Σi=1nαxiCx(xi,xp)+Σj=1mβyjCc(xi,yj)+μx=Cx(0,xi),p∈[1,n]Σi=1nαxiCc(yj,xi)+Σj=1mβyjCy(yq,yj)+μy=Cy(0,yj),q∈[1,m]Σi=1nαxi=1Σj=1mβyj=0]]>

其中，C_x(x_i,x_p)为n个点中的第i点的坐标位置x_i与n个点中的第p点的坐标位置x_p对应的微测井资料的速度值的变差函数值，C_x(0,x_i)为待求位置0与n个点中的第i点的坐标位置x_i对应的微测井资料的速度值的变差函数值，C_y(y_q,y_j)为m个点中的第q点的坐标位置y_q与m个点中的第j点的坐标位置y_j对应的风化层速度值的变差函数值，C_y(0,y_j)为待求位置0与m个点中的第j点的坐标位置y_j对应的风化层速度值的变差函数值，C_c(x_i,y_j)和C_c(y_j,x_i)为关于n个点中的第i点的坐标位置x_i对应的微测井资料的速度值与m个点中的第j点的坐标位置y_j对应的风化层速度值的交叉编译函数，为协同克里金方程组中与n个点中的第i点的微测井数据的速度值对应的加权系数，为协同克里金方程组中与m个点中的第j点的风化层速度值对应的加权系数，μ_x、μ_y为计算方程组的解与原始值的误差的参数。