[发明专利]一种提高测井横波速度预测精度的方法及系统

权利要求书

1.一种提高测井横波速度预测精度的方法，其特征在于：所述方法包括：统计工区内的特征井的孔隙结构参数与孔隙度、泥质含量的回归关系，利用该回归关系预测该工区内的待预测井的横波速度；

所述特征井为该工区内具有横波测井资料的测井；

所述待预测井为该工区内不具有横波测井资料的测井，即需要预测横波速度的测井。

2.根据权利要求1所述的提高测井横波速度预测精度的方法，其特征在于：所述方法包括：

(1)输入特征井的测井数据，包括孔隙度、泥质含量、孔隙结构参数，统计特征井的孔隙结构参数与孔隙度、泥质含量的回归关系；

(2)输入待预测井的孔隙度、泥质含量、含水饱和度和纵波速度，将待预测井的孔隙度、泥质含量代入到步骤(1)得到的回归关系中求得该待预测井的孔隙结构参数计算值L；

(3)根据所述待预测井的孔隙结构参数计算值L设定待预测井的孔隙结构参数的线性扫描范围[(1-a)·L,(1+a)·L]，同时设定步长；

(4)在所述线性扫描范围内对待预测井的孔隙结构参数进行线性扫描，找到待预测井的最优孔隙结构参数：

(5)输出结果。

3.根据权利要求2所述的提高测井横波速度预测精度的方法，其特征在于：所述步骤(3)中的a取值为5％。

4.根据权利要求2所述的提高测井横波速度预测精度的方法，其特征在于：所述步骤(4)的操作具体包括：

(41)设孔隙结构参数的初始值为(1-a)·L；

(42)利用所述孔隙结构参数计算得到纵波速度；

(43)计算步骤(42)获得的纵波速度与步骤(2)输入的待预测井的纵波速度之间的误差fo：

(44)判断fo≤ε是否成立，如果是，则将使得fo≤ε成立的孔隙结构参数作为最优孔隙结构参数，然后转入步骤(5)；如果否，则转入步骤(45)；

(45)判断孔隙结构参数是否大于(1+a)·L，如果是，则将计算得到的所有误差进行排序，找到排序中的最小值，该最小值对应的孔隙结构参数即为最优孔隙结构参数，然后转入步骤(5)，如果否，则将孔隙结构参数加上一个步长作为新的孔隙结构参数，然后返回步骤(42)。

5.根据权利要求4所述的提高测井横波速度预测精度的方法，其特征在于：所述步骤(42)的操作包括：

利用所述孔隙结构参数、输入的待预测井的孔隙度、泥质含量、含水饱和度和岩石物理模型正演计算得到纵波速度RPM(α_s,φ,Vsh,Sw)和横波速度，其中，α_s,φ,Vsh,Sw依次为孔隙结构参数、孔隙度、泥质含量、含水饱和度。

6.根据权利要求5所述的提高测井横波速度预测精度的方法，其特征在于：所述步骤(43)的操作包括：

采用下式计算得到误差fo：

fo＝||RPM(α_s,φ,Vsh,Sw)-VPm|| (1)

所述公式(1)中的VPm为输入的待预测井的纵波速度。

7.根据权利要求6所述的提高测井横波速度预测精度的方法，其特征在于：所述步骤(44)中的ε是设定的阈值，ε的取值为0.1。

8.根据权利要求4所述的提高测井横波速度预测精度的方法，其特征在于：所述步骤(5)的操作包括：

将所述最优孔隙结构参数输出，同时输出由该最优孔隙结构参数正演计算得到的横波速度，该横波速度即为待预测井的横波速度。