基于时频谱图像特征分析的油水识别方法技术

本发明专利技术提出了一种基于时频谱图像特征分析的油水识别方法，包括以下步骤：步骤S10：提取待识别油水性质的井点目的储层的地震资料，对地震资料进行频谱信号分析。根据频谱分布范围和地震主频，利用匹配追踪时频谱图像分析技术获得过井地震道的时频谱能量分布；步骤S20：对时频谱进行归一化处理；步骤S30：根据时频谱中井点目的储层所对应的时间区间中时频谱图像特征来识别储集体中的流体充填性质。采用本发明专利技术中的方法对碳酸盐缝洞储集体中的流体充填性质进行识别时，流体识别结果与实际生产井吻合率高，利用时频谱特征分析识别碳酸盐岩储层的储集体中流体充填性质的结构可靠。

【技术实现步骤摘要】
基于时频谱图像特征分析的油水识别方法
本专利技术涉及一种基于地震时频谱分析识别超深层碳酸盐岩缝洞储集体流体充填性质的方法，属于石油地球物理勘探中的油水识别领域。
技术介绍
对于含油气储层的流体识别，国内外学者在砂岩储层中进行了大量研究，但针对碳酸盐岩储层的流体识别尚处于探索研究阶段。国内学者针对BL地区长兴组地层碳酸盐岩礁滩相岩性油气藏具有埋藏深、油气分布不均匀的特征，运用黏滞弥散型波动方程模拟储层地震响应特征，利用频率衰减梯度属性和低频伴影现象识别储层含气性。还有研究人员开展了基于频谱特征差异的储层流体识别研究，以高频段振幅差异特征为依据识别洞穴与坍塌－充填洞穴的流体性质。目前基于匹配追踪算法进行地震信号的分解重构分解的技术被应用到地震属性分析研究的各个领域中，有学者利用基于雷克子波原子库的MP算法对地震信号进行稀疏分解，进而可实现对地下地质体的高精度波阻抗反演，另有将该算法用于瞬时谱识别来确定三角洲砂体尖灭线，以及利用MP谱分解算法进行峰值瞬时频率薄层反演，同时进行含气储层的烃类检测。首先，现有基于叠前地震资料的储层流体识别技术计算量大、对资料要求高。而且在碳酸盐岩缝洞型储集体中，油藏的主体是规模较大的流体充填溶洞，而大型洞穴中的流体不传播横波，无法提取像砂岩储层中辨识油、气、水的有效弹性参数，如纵、横波速度和弹性模量等，因此叠前弹性反演技术的应用受到限制。另外，基于地震属性分析的识别方法受地层埋深、岩性、厚度、温度、压力等多种因素的影响，多解性较强，运算结果的岩石物理意义不够清晰，解释难度较大。
技术实现思路
针对上述技术问题，本专利技术提出了一种基于时频谱图像特征分析的油水识别方法，包括以下步骤：步骤S10：提取待识别油水性质的井点目的储层的地震资料，对地震资料进行频谱信号分析。根据频谱分布范围和地震主频，利用匹配追踪时频谱图像分析技术获得过井地震道的时频谱能量分布；步骤S20：对时频谱进行归一化处理；步骤S30：根据时频谱中井点目的储层所对应的时间区间中时频谱图像特征来识别储集体中的流体充填性质。在一个具体的实施例中，当所述时间区间中的时频谱图像呈朝向频率减小的方向凸出的拱形，并且在此时间区间内的时频谱的频率小于主频，则所述井点目的储层的储集体中的流体含油量比含水量高。在一个具体的实施例中，当所述时间区间中的时频谱图像呈朝向频率增大的方向凸出的拱形，并且在此时间区间内的时频谱的频率大于主频，则所述井点目的储层的储集体中的流体含水量比含油量高。在一个具体的实施例中，当所述时间区间中的时频谱图像向高频方向迁移，并且在此时间区间内的时频谱的频率大于30Hz，则所述井点目的储层的储集体中不含水也不含油。在一个具体的实施例中，当所述时间区间中的时频谱图像沿平行于时间轴的方向延伸，并且在此时间区间内的时频谱的频率大于30Hz，则该储层中没有储集体发育。在一个具体的实施例中，步骤S10包括步骤S01～S03，步骤S01：对地震信号进行稀疏分解，构建字典矩阵D；步骤S02：采用匹配追踪算法，把地震信号函数y垂直投影到字典矩阵D的原子上进行重复的迭代估算，直至残值信号小于或等于预设的阈值，将进行迭代过的原子的线性和再加上最后的残差值来表示地震信号；步骤S03：采用时频谱来表示地震信号分布。在一个具体的实施例中，在步骤20中，时频谱中的每个向量||xi||＝1。在一个具体的实施例中，主频为25Hz。采用本专利技术中的方法对碳酸盐缝洞储集体中的流体充填性质进行识别时，流体识别结果与实际生产井吻合率高，利用时频谱特征分析识别碳酸盐岩储层的储集体中流体充填性质的结构可靠。附图说明在下文中将基于实施例并参考附图来对本专利技术进行更详细的描述。其中：图1为本专利技术的一个实施例中的基于时频谱图像特征分析的油水识别方法的流程图；图2为本专利技术的一个实施例中井点目的层地震资料的频谱分布图；图3为本专利技术的一个实施例中未进行归一化处理的过井地震道的时频二维谱图像；图4为本专利技术的一个实施例中完成归一化处理的过井地震道的时频二维谱图像；图5为本专利技术的一个实施例中显示了储集体中含油量较高的储层的时频谱特征的时频谱图像；图6为本专利技术的一个实施例中显示了储集体中含水量较高的储层的时频谱特征的时频谱图像；图7为本专利技术的一个实施例中显示了不含储集体的储层的时频谱特征的时频谱图像；图8为本专利技术的一个实施例中显示了储集体发育程度低的储层的时频谱特征的时频谱图像。在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。具体实施方式下面将结合附图对本专利技术作进一步说明。图1显示了本专利技术的一个实施例中的基于时频谱图像特征分析的油水识别方法的流程图，该油水识别方法用于识别碳酸盐缝洞储集体中的流体充填性质。该油水识别方法包括以下步骤：步骤S10：提取待识别油水性质的井点目的储层的地震资料，对地震资料进行频谱信号分析。根据频谱分布范围和地震主频，利用匹配追踪时频谱图像分析技术获得过井地震道的时频谱能量分布。步骤S10包括步骤S01～S03。图2中的显示了塔河油田的井点储层的频谱分布状况，塔河油田的井点储层频谱范围为5-70Hz，主频为25Hz。步骤S01：对地震信号进行稀疏分解，构建字典矩阵D。字典矩阵D为进行原始信号稀疏分解后的函数集合。构建字典矩阵D的方法为本领域的现有技术，在此不再赘述。在本实施例中，将井点目的储层地震信号函数y采用一组向量{x1,x2,x3,...xn}来表示，被表示的地震信号的长度为n，这组向量构成字典矩阵D。字典矩阵D中的每个向量x为一个原子。每个原子由其所对应的子波的振幅、中心时间、波峰频率和相位构成的四维向量表示。原子的长度与被表示的地震信号y的长度n相同。步骤S02：采用匹配追踪算法，把地震信号函数y垂直投影到字典矩阵D的原子上进行重复的迭代估算，直至残值信号小于或等于预设的阈值，将进行迭代过的原子的线性和再加上最后的残差值来表示地震信号。选择一个与原始信号匹配最佳的时频原子线性组合，构建一个稀疏逼近，并求出信号残差R，然后继续选择与信号残差最匹配的原子，反复迭代，地震信号y由这些原子来线性和，再加上最后的残差值来表示。给定匹配准则终止条件参数，即如果残差值在预设的阈值以内，则地震信号y就是该原子的线性组合。地震信号函数y的表达式为：y＝<y,xn>xn+Rmy(1)(1)式中，y为井点目的储层地震信号函数；<y,xn>为原始信号y与迭代时所选原子Xn的内积；Rmy为第m次迭代的残差值。在每次迭代过程中都是先找出与原始信号最相关的原子xn，为了使<y,xn>xn与原始信号y最逼近，残差值Rmy尽可能小，内积项<y,xn>应尽可能大。匹配追踪以其优越的自适应特征,采用具有良好时间、频率分辨率的时-频原子对地震信号不断寻找最佳匹配,真实准确地对地震信号进行时频分解,减少了频谱中交叉项的干扰。步骤S03：采用时频谱来表示地震信号分布。将一维的时间信号转换到二维的时间-频率平面上，这样能较好地反映地震信号的时频特征。步骤S20：对时频谱进行归一化处理。对每一个时间—频率谱序列进行归一化，即将时频谱中的每个向量||xi||＝1本文档来自技高网...

【技术保护点】
1.一种基于时频谱图像特征分析的油水识别方法，包括以下步骤：步骤S10：提取待识别油水性质的井点目的储层的地震资料，对地震资料进行频谱信号分析，根据频谱分布范围和地震主频，利用匹配追踪时频谱图像分析技术获得过井地震道的时频谱能量分布；步骤S20：对时频谱进行归一化处理；步骤S30：根据时频谱中井点目的储层所对应的时间区间中时频谱图像特征来识别储集体中的流体充填性质。

【技术特征摘要】
1.一种基于时频谱图像特征分析的油水识别方法，包括以下步骤：步骤S10：提取待识别油水性质的井点目的储层的地震资料，对地震资料进行频谱信号分析，根据频谱分布范围和地震主频，利用匹配追踪时频谱图像分析技术获得过井地震道的时频谱能量分布；步骤S20：对时频谱进行归一化处理；步骤S30：根据时频谱中井点目的储层所对应的时间区间中时频谱图像特征来识别储集体中的流体充填性质。2.根据权利要求1所述的油水识别方法，其特征在于，当所述时间区间中的时频谱图像呈朝向频率减小的方向凸出的拱形，并且在此时间区间内的时频谱的频率小于主频，则所述井点目的储层的储集体中的流体含油量比含水量高。3.根据权利要求1所述的油水识别方法，其特征在于，当所述时间区间中的时频谱图像呈朝向频率增大的方向凸出的拱形，并且在此时间区间内的时频谱的频率大于主频，则所述井点目的储层的储集体中的流体含水量比含油量高。4.根据权利要求1所述的油水识别方法，其特征在于，当所述时间区间中的时频谱图像向高频方向迁移，并且在此时间区间内的时频谱的频率大于30...

【专利技术属性】
技术研发人员：解丽慧，邬兴威，司朝年，王萍，刘坤岩，韩东，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京,11