[发明专利]一种回归-分形插值法反演测井横波时差方法

权利要求书

1.一种回归-分形插值法反演测井横波时差方法，实现的步骤如下：

第一步小层对比，建立研究地区、研究层位地层格架；

第二步利用常规测井曲线，通过多元回归建立横波时差预测初步模型；

利用阵列声波测井，结合小层对比结果，分析确定与横波时差密切相关的常规测井曲线，分层建立横波时差的多元回归预测初步模型；

第三步横波时差二维分形插值测井点循环；

利用横波时差预测初步模型，拟合得到的横波时差数据和埋深，构建数据集U，从数据集U中间隔为D，其中D≥3，抽取测井数据点构建数据集{(x_i，y_i)|i＝0，1，……，N}，对剩余的测井数据点进行二维分形插值计算；这样有N个仿射变换，第n个仿射变换的分形插值函数为：

公式(4)中的相关参数表示为：

公式(4)-(5)中，x_i为埋深，y_i为多元回归初步预测横波时差；d_n为纵向压缩比；

在确定D的基础上，利用公式(4)-(5)得到单次分形数据集Q₁，完成一次分形插值；移动插值数据点位置，移动步长为一个测井数据点间距；同样在间隔D条件下，选取未抽取过的测井数据点，构建新的数据集{(x’_i，y’_i)|i＝0，1，……，N}，依次计算得到分形数据集标记为Q₂……Q_D；在{Q₁，Q₂……Q_D}数据集中，每个测井位置对应的插值数据点取平均值得到数据Q_aver，调整d_n，使分形插值后数据集Q_aver中横波时差与真实横波时差之间的拟合系数R²达到最大，并记录对应是R²₁；

第四步横波时差二维分形插值间距循环；

依次增大间隔D，重复第三步，分别记录对应的R²₂……R²_d，当R²_d出现下降时，停止增大间隔D，记录对应的d_n；

第五步分层确定横波时差预测模型；

选取R²_d最大的模型对应的间隔D和d_n作为单层的横波时差预测模型，分别确定不同层位的横波时差预测模型，建立全区横波时差预测模型，对缺少横波时差的井进行预测。