[发明专利]一种碳酸盐岩储层中火成岩侵入范围的识别方法

说明书

本发明公开了一种碳酸盐岩储层中火成岩侵入范围的识别方法，包括初步分析、建过完备字典和判断火成岩的侵入范围等步骤。本发明公开的一种碳酸盐岩储层中火成岩侵入范围的识别方法，根据稀疏分解思想，利用多道匹配追踪从过完备字典中获取的最优匹配子波，对侵入火成岩地层强反射轴进行追踪识别、分离，提高了计算的效率，更能保证原始地震信号的完整性，而且抗噪性增强，更加精确地确定火成岩的分布范围，改善地震资料的可靠性。

技术领域

本发明涉及海洋地质领域，具体涉及一种碳酸盐岩储层中火成岩侵入范围的识别方法。

背景技术

在生物礁、生屑滩等碳酸盐岩储层中，火成岩的侵入造成碳酸盐岩储层描述不精确，在一定条件下甚至会破坏碳酸盐岩储层的储集性能，因此，如何准确的刻画火成岩侵入范围对于碳酸盐岩储层储量的准确估算来说具有重要的意义。当火成岩侵入的规模和波及范围不是特别大时，由于火成岩速度较高，通常会在地震剖面上呈现出一定的反射特征，在实际应用中通常采用井震精细标定来确定火成岩的反射特征，并采用对火成岩特征敏感的典型地震属性特征开展人工地震解释，进一步结合火成岩的地震相分布特征对火成岩的侵入范围进行描述，但该方法效率低，且在火成岩厚度变化较快的区域开展人工解释时往往会产生较大的误差，导致火成岩侵入范围识别不准确。如何从地震反射信号中高精度地识别出火成岩反射特征，并确定出火成岩侵入的规模大小和准确刻画出火成岩的空间分布特征及其对储层的影响，是开展碳酸盐岩储层的精确描述和储量计算时面临的难题。

发明内容

本发明的目的在于提供一种碳酸盐岩储层中火成岩侵入范围的识别方法，用以解决开展碳酸盐岩储层的精确描述和储量计算时火成岩侵入范围识别不准确的问题。

本发明提供的一种碳酸盐岩储层中火成岩侵入范围的识别方法，所述识别方法包括以下步骤：

步骤S1：初步分析是否有火成岩侵入以及火成岩在地震剖面对应时窗处的强振幅反射特征；

步骤S2：在众多强振幅反射特征中寻找最优匹配子波，利用多个的最优匹配子波组建过完备字典；

步骤S3：利用过完备字典重构由侵入岩造成的强反射干扰地震数据，识别地震资料中强反射地震信息，据此判断火成岩的侵入范围。

优选地，所述步骤S1包括以下步骤：

步骤S11：分析地震资料的地质背景和区域构造，判定该区域是否有火成岩侵入的影响；

步骤S12：在地震资料的地质背景和区域构造的基础上，结合井旁道地震资料和测井资料开展如振幅、相位、频率等强反射地震属性分析处理工作，请从中关注火成岩在地震剖面对应时窗处的强振幅反射特征。

优选地，所述步骤S2包括以下步骤：

步骤S21：大致估算出地震信号对应的子波参数范围；

步骤S22：对原始地震信号进行稀疏分解，从过完备字典中选择与原始信号最为匹配的匹配子波。

优选地，所述步骤S22包括以下步骤：

步骤S221：对相邻的地震道信号利用匹配追踪算法进行稀疏分解，通过对过完备字典中子波进行遍历搜素，就是计算原始信号与过完备字典矩阵中的原子进行内积计算，选取内积绝对值最大的一个原子得到最优的匹配子波，其表达式如下：

其中，t为时间变量，y(t)为原始地震信号，为与相邻地震道内积值最大的子波，即在该次迭代运算中的与地震信号y(t)最匹配的子波，R^(n-1){y(t)},为y(t)和第n-1次匹配后的信号残差，每一次迭代求出原始信号与重构信号的差值；

步骤S222：在上一次循环得到的残差信号基础上，再进行循环迭代，不停循环往复，直至循环结束。