同步压缩小波变换分频谱蓝化在薄层识别中的应用_刘庆文.pdf

刘庆文，孙永壮，李琴 同步压缩小波变换分频谱蓝化在薄层识别中的应用 石油物探，（）：，（）：收稿日期：。第一作者简介：刘庆文（），男，硕士，工程师，现主要从事储层地球物理研究工作。：基金项目：中国海油“七年行动计划”东海专项课题（）和中国海油“十四五”重大科技项目“海上深层超深层油气勘探技术”（）共同资助。“”（）“”（）同步压缩小波变换分频谱蓝化在薄层识别中的应用刘庆文，孙永壮，李琴（中海石油（中国）有限公司上海分公司研究院，上海 ）摘要：传统谱蓝化拓频技术由于谱蓝化算子的单一性，蓝化拓频后频谱分布异常，从而导致地震资料高频段存在不合理的现象；同时，蓝化拓频后频谱低频段容易出现震荡，导致拓频效果不佳。为此，提出基于同步压缩小波变换的分频谱蓝化拓频技术。该技术首先利用同步压缩小波变换对频率的分辨率优势，将地震数据精确地分成多个分频数据体；然后，在分频段计算谱蓝化算子过程中，考虑不同频段振幅谱特征，提出一种基于不同频带能量差异加权的改进谱蓝化算子；再依据不同分频段振幅谱占比确定其对应的改进谱蓝化算子权重；最后，将改进谱蓝化算子与地震反射系数褶积，得到优化后的拓频数据。模型试算及实际地震资料测试结果表明：相较于传统谱蓝化技术，联合同步压缩小波变换分频及谱蓝化算子优化进行地震拓频，较大程度提高了地震资料的高频信息，拓频后地震数据分辨率更高，在保证薄层精细刻画下，其低频段信息也得到有效保留，视觉假分辨率得到一定程度的压制。关键词：谱蓝化；同步压缩小波变换；时频分析；算子改进；薄层识别中图分类号：文献标识码：文章编号：（）：，（.，.，）：，第 卷第期 年月石油物探 ，：，随着海上成熟油气区块勘探开发的不断深入，薄互层岩性油气藏逐渐成为新增探明储量的重要研究领域，但中深层地震资料频带宽度有限、分辨率低，薄层刻画难度较大，为了解决上述问题，资料处理人员多采用地震数据高分辨率处理技术。常用的高频恢复技术包括谱白化、反褶积、反滤波等，这些技术都很好地改善了实际地震资料分辨率，但存在各自的适用条件或缺陷。谱白化方法存在纵向时间、横向空间补偿不均衡问题；反褶积方法对地层反射系数及子波均有一定假设条件，且时变子波获取难度大；反滤波在高频补偿问题上需要估算值。谱蓝化拓频首先对测井数据反射系数与地震信号进行频谱分析，得到两者匹配的谱蓝化算子，然后将其与原始地震反射系数褶积得到谱蓝化拓频数据，其中，谱蓝化算子的计算是核心步骤。认为谱蓝化拓频技术可以提高叠后地震资料分辨率，并结合有色反演进行了叠后地震高分辨率处理及薄层识别。等对比分析了叠后及叠前谱蓝化拓频技术在 区的二氧化碳预测效果，论证了叠前谱蓝化压制边缘震荡效应的优势，证明了其拓频效果更佳。近年，国内学者杨瑞召等、刘建伟等、陈文雄等、李贤兵等 利用谱蓝化技术拓宽地震数据中高频段信息，实现地震数据的高分辨率处理，较好地解决了薄煤层识别以及薄储层刻画中的一些问题，但这些方法技术中的谱蓝化算子均基于地震整体频谱计算，未考虑地震不同频带信息的特点，导致了拓频后的地震数据出现分辨率不合理、低频震荡等问题。对此，陈平等 提出基于谱反演的谱蓝化拓频技术，通过压缩子波影响，拓宽地震低频与高频，提高地震资料分辨能力；张璐等 提出地震分频谱蓝化思想，借助变换对地震数据分频，然后分别求取谱蓝化算子，较大程度地改善了薄层刻画的精度，但变换存在时频分辨率精度不够等问题。为了优化谱蓝化地震拓频效果，解决薄层刻画难等问题，本文提出了一种基于同步压缩小波变换的分频谱蓝化拓频技术，利用同步压缩小波变换时频分辨率优势完成地震数据的分频，在分频数据体谱蓝化算子计算过程中，利用不同频带频谱能量差异改进谱蓝化算子，通过不同分频数据体振幅谱占比确定其对应的改进谱蓝化算子权重，并将改进后的谱蓝化算子与地震数据进行褶积，得到精度更高、更合理的拓频地震数据，进而将其用于高分辨率薄层识别。同步压缩小波变换方法原理同步压缩小 波 变 换（，）由 等 在 年提出，其基本思想是在频率域压缩小波变换的时频结果，使能量尽量集中到信号真正的频率附近，使得计算结果的分辨率和能量聚集性更高。该技术的实现主要包括以下步。）求取信号的连续小波变换结果，获取复小波系数谱（，）。对于信号（），其连续小波变换公式计算如下：（，）（）（）（）式中：（，）为待求信号的复小波系数谱；为时间；，分别为尺度因子、时间平移因子；（）为小波基函数的复共轭。结合 定理，得到频率域复小波系数谱为：（，）?（）?（）（）式中：为角频率；?（），?（）分别是（），（）的傅里叶变换结果。）求取信号（）的瞬时频率（，）。等 的测试分析表明，信号（）的小波系数在时间 尺度域容易发生扩散，进而降低小波系数的可读性，而信号的相位具有稳定性，不随值变化而改变。因此，计算其相位偏导数就能求得瞬时频率（，），如下：（，）（，）（，）（，）（，）（）石油物探第 卷）在频率方向对小波系数压缩重组，得到同步压缩小波变换量值（，）。将复小波系数从时间 尺度域（，）转换到时间 频率域（，），在频率域将其压缩，此时频率变量和尺度因子被离散化。对复小波系数（，），的任一中心频率附近，的值进行压缩重组，这样就将小波系数从频率方向发散的区域压缩到中心频率附近，从而大幅度提高频率分辨率。由于实际信号处理过程中得到的均为时间域离散值，所以要将尺度离散化，同步压缩变量（，）可表示为：（，）（）：（，）（，）（）式中；是第个离散频率值；为第个离散尺度因子；。采用上述技术，可以使各时频脊线细化，从而减弱频率混叠现象。分析可知，同步压缩小波变换以小波变换为基础，对相位求偏导数得到瞬时频率，然后将小波系数进行压缩重组。在时间 尺度域到时间频率域转变的过程中，能量聚焦性得到提升，从而提高了时频分析结果的分辨率。模型测试为了验证同步压缩小波变换时频分析高分辨率特性，建 立 如 图 所 示 的地震 信号，模 型 长 度 为 ，采样间隔，纵向包含个有效信号，分别在 ，信号振幅分别为，其中 与 信号时间间隔较小。对选取的地震信号分别进行短时傅里叶变换（）、广义变换（）及同步压缩小波变换（），结果如图、图 和图 所示，可以看出短时傅里叶变换能量聚焦性差，对 ，的两个信号区分度不够。广义变换后，时频域分辨率明显提升，纵向时间分辨率与信号匹配性高，但当信号间隔较小时（），频率域能量聚焦存在一定误差。同步压缩小波变换效果最佳，不仅时间分辨率高，可以较好地区分 ，的两个相邻信号，且频率能量聚焦在 附近，时频域分辨率高。为了进一步验证该本文所述方法的时频分析的稳定性，设计如图所示的添加 高斯噪声的模型，图 为含噪声地震信号，分别应用上述种时频分析方法，结果如图，图 和图 所示。分析可知，相较于其它方法，同步压缩小波变换抗噪性更强，能量聚集性依然最好，时频域分辨率最高。图不同时频分析方法的效果地震信号；短时傅里叶变换；广义变换；同步压缩小波变换第期刘庆文等 同步压缩小波变换分频谱蓝化在薄层识别中的应用图不同时频分析方法的效果对比（含 噪声）含噪地震信号；短时傅里叶变换；广义变换；同步压缩小波变换 分频谱蓝化拓频技术技术原理及流程传统谱蓝化拓频技术利用地震频率与振幅的正相关“蓝谱”特征，先将地震和测井的反射系数振幅谱进行匹配以求取谱蓝化算子，再利用井控的拓频思路来提高地震数据的分辨率。传统谱蓝化算子由地震数据整体平均频谱与测井数据对比计算得到，因此缺乏对地震数据不同频带信息的充分匹配，导致拓频后的地震数据出现低频信息损失、精度不够等问题。本文在传统谱蓝化拓频基础上，提出基于同步压缩小波变换的分频谱蓝化拓频技术，其具体流程如下。）对地震信号进行 分频，并计算不同频段的地震平均频谱（）。（）（）（）式 中：（）为 第个 分 频 段 的 地 震 平 均 频 谱；（）为第个频段中第道地震振幅谱；为该分频段的地震道数。）测井反射系数频谱计算。目标井反射系数表达式为：（）（）（）式中：为密度；为速度；为井上纵向第个采样点；为第层反射系数。利用傅里叶变换，可以进一步得到：（）（）（）（）式中：（）为测井反射系数频率域信号；（）为测井反射系数振幅谱；（）为相位谱。）不同分频段的谱蓝化算子计算。计算不同分频段与测井反射系数频谱匹配的谱蓝化算子（），使得两者差异最小化：（）（）（）（）（）（）式中：（）为锥度函数；（）为测井反射系数振幅谱；（）为地震平均振幅谱；（）为正则化函数。依据（）式计算每一个分频段的谱蓝化算子。石油物探第 卷）分频段谱蓝化权重因子计算。首先对分频段振幅谱能量进行归一化处理，具体如下：（）（）（）式中：（）为第个分频段振幅的谱能量；，分别为原始地震谱能量最大值和最小值；（）为第个分频段归一化后振幅谱能量。根据不同分频段能量占比，确定其谱蓝化权重因子。（）（）（）式中：为分频段个数。）对不同分频段的谱蓝化算子（）进行加权求和，得到优化后的谱蓝化算子 。（）（）计算优化拓频数据。先利用循环反褶积技术对地震数据进行重采样，生成稀疏的反射系数体，然后将优化后的谱蓝化分频算子与进行褶积，得到拓频数据：（）因为求取了不同频段的谱蓝化算子，故优化后的谱蓝化算子更能表征不同频段信号。利用分频谱蓝化算子拓频后的地震数据体精度高，该技术不仅有效地拓宽了高频信息，提高薄层刻画精度，同时，有效解决了传统谱蓝化技术低频段损失问题，使得优化后的拓频数据对储层的刻画更加合理。模型测试为了验证该技术高分辨率处理的优势，设计如图所示的薄互层模型并添加 高斯噪声，黄色填充为砂岩，灰色为泥岩，套砂岩均为 薄层，沿着 测线方向泥岩隔层厚度逐渐增大，模型参数为实钻数据，砂岩速度为 ，砂岩密度为 ，泥岩速度为 ，泥岩密度为 ，子波选用 的零相位雷克子波。图为根据该模型正演得到的地震记录、传统谱蓝化拓频和 分频谱蓝化拓频结果。分析可知，泥岩隔层为 处（图 红色箭头），第一套砂岩底出现波峰，而泥岩隔层厚度约为 处（图 黑色箭头），第二套砂岩底出现波峰，且薄层的波形厚图薄互层模型（含 噪声）图模型正演地震记录（）、传统谱蓝化拓频结果（）与 分频谱蓝化（）技术拓频结果度明显超过实际厚度，薄层分辨率整体不高。相较于原始地震记录，图 中第一套砂底波峰出现在隔层厚度处（图 红色箭头），传统谱蓝化拓频方法虽然无法有效区分第一套与第二套砂体，但第二套砂岩底波峰较原始地震记录出现更早，泥岩隔层厚度约（图 黑色箭头），整体薄层分辨能力得到了较大程度的提升。对比图 和图 可知，分频谱蓝化的第一套砂岩底波峰相较于传统谱蓝化拓频结果有了进一步提升（图 红色箭头），泥岩隔层厚度约为（图 黑色箭头）处第二套砂岩底出现波峰，薄层识别精度更高。相较于传统谱蓝化拓频，分频谱蓝化技术的抗噪性更好，在泥岩隔层偏薄的情况下，薄层识别精度受噪声影响较小，且厚度与实际模型匹配度更高（图 蓝色箭头）。图为原始地震记录频谱与两种谱蓝化拓频的频谱，蓝色为原始地震记录频谱，绿色为传统谱蓝化拓频得到的频谱，红色为 分频谱蓝化拓频得到的频谱。对比可知，利用传统谱蓝化拓频方法得到的地震资料高频段信息得到了一定程度的拓宽，但 以第期刘庆文等 同步压缩小波变换分频谱蓝化在薄层识别中的应用图原始地震记录频谱与两种谱蓝化拓频技术得到的频谱下的低频信息存在损失，有效频带（信噪比）范围为 ；而利用 分频谱蓝化拓频技术得到的结果中，中、高频段信息更丰富，且有效地保留了原始地震低频段信息，有效频带范围为 ，故该技术对地震资料拓频效果更好。应用实例为了验证本文技术的有效性，以 凹陷中深层某气田 井为例进行分析。图为过该井的原始地震剖面，该井在 附近钻遇两套优质薄气层 、，厚度分别为 和，两层之间泥岩隔层为。受限于凹陷中深层地震资料品质，无法有效区分两套薄互层，地震剖面上 砂顶表现为弱波峰，与下覆 表现为典型复波叠置特征。对此，首先应用 分频谱蓝化技术进行分频及优化谱蓝化算子求取。为了提高靶区拓频效果，此次选用 井进行测井谱与不同频段地震谱匹配算子求取，其中，地震谱计算时采用 井邻道目的层段地震数据。图为中心频率为 的分频段地震频谱与井谱蓝化算子求取过程，图 为井目的层段反射系数谱，图 为分频段地震数图过井的原始地震剖面图中心频率为 的分频段地震频谱与井谱蓝化算子求取过程 井目的层段反射系数谱；分频段地震数据在 井邻道目的层段地震频谱；分频段谱蓝化算子匹配；时间域分频段谱蓝化算子石油物探第 卷据在井邻道目的层段的地震频谱，蓝色为平滑后的振幅谱，图 为分频段谱蓝化算子匹配结果，其中蓝色为地震平均谱，绿色为测井反射系数谱，红色曲线代表该分频段的频率域谱蓝化算子，图 为该算子的时间域形态。类似地，对其它分频段分别求取分频段谱蓝化算子，最后利用（）式和（）式分别计算优化后的谱蓝化算子及优化后的拓频地震数据。分别应用传统谱蓝化及 分频谱蓝化两种技术对 段、段地震资料进行拓频处理，得到的剖面如图所示。不难发现，两种技术均能提高 两套薄互层地震分辨率，且高频信息有所提高，但在原始地震低频信息的保真上，传统谱蓝化拓频在一定程度上破坏了低频段信息，且横向分辨率存在损失（图 蓝色箭头），而该技术的局限性也使得 段振幅能量失真，相较于原始地震剖面，其振幅能量过强（图 黑色箭头）。图对 ，段地震资料应用两种分频谱蓝化拓频技术得到的剖面传统谱蓝化拓频；分频谱蓝化拓频图为分别应用两种谱蓝化拓频技术在 段得到的地震频谱，黑色为原始地震频谱，蓝色为传统谱蓝化拓频后的地震频谱，红色为 分频谱蓝化拓频后的地震频谱。对比可知，对于中、高频信息，两种技术相较于原始地震信号均有一定程图应用两种谱蓝化拓频技术得到的地震频谱（）度的拓宽，其中，采用 分频谱蓝化拓频技术得到的结果中高频信息更丰富；对于低频信息，采用传统谱蓝化拓频技术造成 以下的信息有较大损失，而采用 分频谱蓝化拓频技术可以有效地保留原始地震资料的低频段信息，保幅性更好。为了更直观地对比两种谱蓝化拓频处理得到的地震资料对薄储层的预测精度，对拓频后的地震资料进行 相移，将地震界面信息判别转为更直观的岩性厚度识别。图 为 相移后的原始地震资料岩性预测结果，图 和图 分别为传统谱蓝化拓频后及 分频谱蓝化拓频后 相移地震资料岩性预测结果。分析可知，受限于地震资料的分辨率，原始地震资料对薄气层段 预测效果较差，仅表现为虚弱反射（图 黑色箭头）；传统谱蓝化拓频技术可以明显改善对 段的岩性预测结果，但由于拓频造成地震横向分辨率损失，预测的砂体展布连续第期刘庆文等 同步压缩小波变换分频谱蓝化在薄层识别中的应用广泛，故对 段下覆岩性假象有一定放大作用，实钻泥岩段预测为砂岩（图 黑色箭头）；利用 分频谱蓝化拓频技术得到的地震资料对 段岩性预测更合理，在保留原始地震横向分辨率情况下，该技术较好地提高了 段纵向分辨率，且对 段下覆岩性假象有一定压制效果（图 蓝色箭头）。图 、图 和图 分别显示原始地震资料以及分别采用传统谱蓝化及 分频谱蓝化拓频得到的过 井 段砂体横向展布情况。相较于原始地震数据，两种谱蓝化技术均可以较好地表征图 相移后的原始地震资料（）、传统谱蓝化拓频后 相移地震资料（）和 分频谱蓝化拓频后 相移地震资料（）岩性预测结果图 采用两种谱蓝化拓频技术得到的过 井 段砂体平面展布对比原始地震资料；传统谱蓝化拓频；分频谱蓝化拓频石油物探第 卷井点 段砂体，但其平面展布形态存在差异，地质分析认为 井 段砂体为受潮汐改造、断坡遮挡形成的砂坝，图 中 分频谱蓝化拓频后砂体展布与上述认识一致，在北侧断层尖灭处砂体不发育。图 中 段砂体表现为北东向连片展布，在没有断坡遮挡的区域发育砂体，这与实际地质认识不符。对比分析表明，分频谱蓝化拓频技术对薄互层砂体纵横向刻画精度更高。结论谱蓝化地震高分辨处理技术的核心是谱蓝化算子。针对传统谱蓝化算子计算方法的单一性，本文联合同步压缩小波变换时频分析对分频谱蓝化算子进行优化改进，获取与实际更匹配的谱蓝化算子，进而完成地震数据的高分辨处理，提高地震薄储层识别能力。分频谱蓝化拓频技术较好地克服了传统谱蓝化中高频带异常、低频损失等不足，在高分辨处理的原则下，较好地保留了低频段信息，且信号保幅性更好。中、深层地震资料普遍存在主频低、信噪比不高等问题，进而影响了储层纵、横向的精细刻画，分频谱蓝化拓频技术在凹陷某气田应用表明，其较好地提高了薄储层纵、横向分辨率，砂体刻画结果与地质认识相符，对薄互储层高分辨率处理具有较好的应用推广价值。参考文献李庆忠走向精确勘探的道路高分辨率地震勘探系统工程剖析北京：石油工业出版社，：，：王江，涂国田，王杰基于载波调制的高分辨率地震双向拓频技术及其应用石油物探，（）：，（）：范小冬，曾濒，刘益强 地震资料时空变谱白化处理 石油地球物理勘探，（）：，（）：，：郭廷超，曹文俊，陶长江，等 时变谱模拟反褶积方法研究 石油物探，（）：，（）：，（）：姬战怀 小波变换在叠后地震数据处理中的应用研究西安：西北工业大学，：，：，：，：，（）：杨瑞召，赵争光，马彦龙，等利用谱蓝化和有色反演分辨薄煤层天然气地球科学，（）：，（）：刘建伟，高秋菊，师涛 谱蓝化技术在大王庄油田储层预测中的应用复杂油气藏，（）：，（）：陈文雄渤海西南部新近系超薄储层定量预测技术研究与应用地球物理学进展，（）：，（）：李贤兵，赵俊杰，晋剑利，等 叠前谱蓝化提频技术在乍得 油田储层预测中的应用石油地球物理勘探，（）：，：，（）：陈平，谭辉煌，秦童，等基于谱反演的谱蓝化拓频方法研究及应用东北石油大学学报，（）：，（）：张璐，汪毓铎 地震数据分频谱蓝化算子计算方法及应用 北京信息科技大学学报，（）：，（）：，：，（）：（编辑：朱珠）第期刘庆文等 同步压缩小波变换分频谱蓝化在薄层识别中的应用