[发明专利]一种结合层位信息的微测井方位加权插值建模方法

权利要求书

1.一种结合层位信息的微测井方位加权插值建模方法，其特征在于，包括以下步骤：

步骤一，确定微测井的主要信息，微测井的主要信息包括微测井速度曲线、解释层位、工区范围、反演所需的微测井数、反演的网格、约束权重，然后输入微测井的主要信息；

步骤二，统计各微测井的距离分布，确定参考距离；

步骤三，基于微测井的参考距离和层位，插值出层位面；

步骤四，统计各层位面的平均值，并将微测井的速度曲线变换为对应层位平均值的速度曲线；

步骤五，基于微测井的参考距离和水平位置，计算径向权重和方位权重；

步骤六，以微测井的方位系数作为未知量，构建反演方程组，以变换后的微测井速度作为右端项，并添加反演约束项；

步骤七，求解反演方程组，获得方位系数。

步骤八，逐点进行方位插值，并将各点的插值结果进行层位变换，获得表层速度模型。

2.根据权利要求1所述的结合层位信息的微测井方位加权插值建模方法，其特征在于：所示步骤二中，对于各微测井，统计出其它微测井与该微测井的距离，基于输入的反演所需微测井数，确定与之相对应的距离，并作为该微测井的参考距离。

3.根据权利要求2所述的结合层位信息的微测井方位加权插值建模方法，其特征在于，在步骤三中，插值层位面是通过径向基函数方法计算得到，具体计算方法如下：某层位面在n口微测井处测得的层位深度为s_i,(i＝1,…,n)，以作为目标函数，式中λ_j为第j口微测井的权重系数，s_j为第j口微测井的层位深度，||x_i-x_j||₂表示第i口微测井与第j口微测井的水平距离，||.||₂为二范数运算，x_i、x_j分别为第i、j口微测井的位置；R为径向权重函数，其形式为：式中，为第i口微测井和第j口微测井的距离与参考距离的比值，r_ref为参考距离；结合最小二乘思想，求解得到各口微测井的权重，并进行全区的层面插值，插值公式如下：其中，s(x)为插值得到的层位面，x为插值点的水平位置。

4.根据权利要求1所述的结合层位信息的微测井方位加权插值建模方法，其特征在于，在步骤四中，以各层位面的平均值为基准，对微测井速度所对应的深度进行线性变换，使对应同一层位的深度统一为该层位平均深度，层位间的深度由上下层位的平均深度线性计算得到，具体计算公式为：其中，H为变换后的深度，h为原始深度，h_i、h_i+1为h对应的上下层位面深度，为相应的层位面的平均深度。

5.根据权利要求1所述的结合层位信息的微测井方位加权插值建模方法，其特征在于，所述步骤五中，每口微测井的径向权重的计算公式为：式中，R(r)为每口微测井的径向权重函数，为第i口微测井和第j口微测井的距离与参考距离的比值，x_i、x_j分别为第i、j口微测井的位置，r_ref为参考距离。

6.根据权利要求1所述的结合层位信息的微测井方位加权插值建模方法，其特征在于，在步骤五中，基于微测井的空间位置，对于各微测井，计算该微测井与其它微测井之间的方位角，并以正北、东北、正东、东南、正南、西南、正西、西北这八个主方位系数加权的形式表示相应方位角的权重，方位系数加权所用函数为二阶三角B样条函数，具体形式如下：其中，t∈[0,1]，为由方位角计算得到的参数，计算公式为：其中θ为方位角，{...}为提取t的小数部分的数学运算符。