一种利用多种测井曲线对声波时差曲线进行重构的方法

摘要

本发明涉及一种利用多种测井曲线对声波时差曲线进行重构的方法，主要包括通过对测井曲线进行储层敏感性与相关性分析，选取对储层特征响应明显的多条测井曲线；将这些测井曲线与声波曲线一起进行离散小波分解，分解层数为8层；将除声波以外的其他测井曲线各层高频分解结果分别组成矩阵，求取矩阵对应的特征值和特征向量，将不同特征值对应的特征向量作为信号重构的新分量，此时这些分量具有正交性(不相关)；各层特征向量对声波曲线小波分解对应层的高频成分进行多元回归分析，计算各向量所对应的加权系数，并返回回归显著性分析结果，判断回归质量；依据回归显著性分析结果判断利用特征向量重构声波高频成分的层数，对选定层数以上的部分选择多元回归结果作为高频成分，对选定层数以下的部分保留声波测井曲线的高频分解结果；运用声波测井曲线的低频成分与回归得到的高频成分进行曲线重构得到最终的声波重构曲线。

说明书

技术领域：

本发明涉及一种利用多种测井曲线对声波时差曲线进行重构的方法。

背景技术：

声波曲线重构是波阻抗反演当中的常用方法。当地下储层速度与泥岩速度差异性小、重合度高的问题；直接进行波阻抗反演的结果不能准确反映储层和围岩的差异，导致波阻抗反演结果与钻井结果不符。因此需要对声波测井曲线进行重构，在声波测井曲线中添加岩性信息，提高声波测井曲线对岩性的识别度，以便进行更准确的储层预测。

多元回归分析是研究多个变量之间关系的回归分析方法，按因变量和自变量的数量对应关系可划分为一个因变量对多个自变量的回归分析(简称为“一对多”回归分析)及多个因变量对多个自变量的回归分析(简称为“多对多”回归分析)，按回归模型类型可划分为线性回归分析和非线性回归分析。

F检验(F-test)，最常用的别名叫做联合假设检验(英语：joint hypotheses test)，此外也称方差比率检验、方差齐性检验。它是一种在零假设(null hypothesis，H0)之下，统计值服从F-分布的检验。其通常是用来分析用了超过一个参数的统计模型，以判断该模型中的全部或一部参数是否适合用来估计母体。

小波变换中，对于非平稳信号而言(测井信号即为非平稳信号)，需要时频窗口具有可调的性质，即要求在高频部分具有较好的时间分辨率特性，而在低频部分具有较好的频率分辨率特性。对测井曲线进行小波分析可以识别测井曲线中不同频率的曲线旋回，高频的曲线旋回对应高频的即短期的沉积旋回，低频的曲线旋回对应低频的即长期的沉积旋回，因此可以用不同频率的测井曲线旋回划分不同周期的沉积旋回，对应不同级别的层序地层单元。

发明内容：

本发明涉及一种利用多种测井曲线对声波时差曲线进行重构的方法，主要包括通过对测井曲线进行储层敏感性与相关性分析，选取对储层特征响应明显的多条测井曲线；将这些测井曲线与声波曲线一起进行离散小波分解，分解层数为8层；将除声波以外的其他测井曲线各层高频分解结果分别组成矩阵，求取矩阵对应的特征值和特征向量，将不同特征值对应的特征向量作为信号重构的新分量，此时这些分量具有正交性(不相关)；各层特征向量对声波曲线小波分解对应层的高频成分进行多元回归分析，计算各向量所对应的加权系数，并返回回归显著性分析结果，判断回归质量；依据回归显著性分析结果判断利用特征向量重构声波高频成分的层数，对选定层数以上的部分选择多元回归结果作为高频成分，对选定层数以下的部分保留声波测井曲线的高频分解结果；运用声波测井曲线的低频成分与回归得到的高频成分进行曲线重构得到最终的声波重构曲线。

附图说明

图1.利用多种测井曲线对声波时差曲线进行重构流程图

图2.特征向量多元回归显著性分析结果

图3.声波测井曲线重构结果

具体实施方式：

如图1所示，方法的实施步骤详细介绍如下：

步骤一，运用储层孔隙度、泥质含量资料与测井曲线进行相关性分析，选择与泥质含量、孔隙度相关性较高的测井曲线；

步骤二，将这些与泥质含量、孔隙度相关性较高的测井曲线进行db小波分解，分解层数为8层；

步骤三，将除声波以外的其他测井曲线各层高频分解结果分别组成矩阵，求取矩阵对应的特征值和特征向量，将不同特征值对应的特征向量作为信号重构的新分量，此时这些分量具有正交性(不相关)；

步骤四，运用这些特征向量对声波曲线小波分解高频成分进行多元回归分析，回归方程中每一项的系数就是该特征向量所对应的的权重，运用F检测判断回归质量会返回两个参数FH和FV；FH为对于多元回归结果可信度的判断；对于给定的置信度α，由F分布表查得的F_α(m，n-m-1)值与FV进行比较，当FV＞F_α(m，n-m-1)时说明多元回归结果可信，FH＝1，反之FH＝0；FV值越大，说明总体回归效果越显著(图2)；

步骤五，根据F检验返回的结果，在每一层回归结果都可信，即FH＝1的情况下，选择FV值最小的那一层，对这一层及其以上所有层的声波曲线高频分解结果用多元回归结果所代替，作为声波曲线重构的高频成分。

步骤六，运用声波测井曲线的低频成分与回归得到的高频成分进行曲线重构得到最终的声波重构曲线(图3)。