中心频率法估算地层吸收参数.docx

中心频率法估算地层吸收参数魏文;李红梅;穆玉庆;王树刚;王红【摘要】本文针对储层含气性检测，提出了利用中心频率估算地层吸收参数(地层Q值)的方法。该方法首先利用S变换对地震信号进行时频域转换，得到能量分布函数;然后基于该函数，计算中心频率及其与地层Q值之间的关系，获得地层等效Q值，进而得到单套地层Q值。模型正演表明，含气层中心频率变小,在地层Q值参数剖面上表现为强吸收异常，从而验证了应用中心频率法估算地层Q值检测气层的可行性。将该方法应用于东营凹陷永安地区实际资料中，估算的地层Q值比较真实地反映了地层吸收特征，提高了气层检测的精度。 【期刊名称】《石油地球物理勘探》【年(卷)，期】2012(047)005【总页数】5页(P735-739)【关键词】品质因子;吸收衰减;S变换;中心频率;气层检测【作者】魏文;李红梅;穆玉庆;王树刚;王红【作者单位】中国石化胜利油田分公司物探研究院,山东东营257022;中国石化胜利油田分公司物探研究院仙东东营257022;中国石化胜利油田西部新区研究中心,山东东营257022;中国石化胜利油田分公司物探研究院仙东东营257022;中国石化胜利油田分公司物探研究院,山东东营257022【正文语种】中文【中图分类】P631实际地层介质并非是理想的完全弹性介质。当地震波在地层中传播时弓I起的能量衰减方式有两类：一类是与地震波传播特性有关的衰减，包括球面扩散、散射及透射损失等，这部分能量在地震资料处理时已对其进行补偿；另一类是地层本身的固有衰减，即地层的吸收，这种介质本身所固有的吸收特性通常用品质因子Q来描述［1,2］,Q值与介质的结构、孔隙度以及孔隙流体的性质密切相关。 在由固态、液态、气态物质构成的多相介质中，对吸收性质影响最显著的是气态物质，在岩石孔隙饱和液中掺入少量气态物质，可以明显提高对纵波能量的吸收，导致地震波的高频成分很快衰减，因此地层吸收参数是检测储层含气性的一个重要参数［3~6］。 当储层含气后，地震资料的振幅和频率都有一定的反映，振幅有异常，在频率特征上通常都有反映，而振幅没有异常，在频率特征上也会有所反映，因此利用地震资料检测储层含气性时，频率特征比振幅特征的应用范围更广泛。本文通过S变换，基于地震频谱信息，利用中心频率估算地层吸收参数（地层Q值），实现对储层含气性的检测。 假设震源为理想脉冲源，地层品质因子Q为常数，同时忽略散射引起的衰减，那么在一维黏弹性介质中，地震波U（3,z）z轴增大方向以角频率3传播时，其传播方程［7］为式中：3是角频率；Z是传播距离；C（3）是相速度。 由于z=tc（3）,3=2nf,则式（1）可表示为将式（2）进行S变换［8］，得到地震能谱中心频率定义为功率谱的频率加权累加和与总能量的比值，它反映了信号能量按频率分配的特点，高频能量强，则中心频率高；高频能量弱，低频部分能量强，则中心频率低，即此时的高频成分吸收衰减严重。其计算公式为式（7）表明中心频率fm与地层Q值成正比关系。文中通过地震资料中心频率估算的Q值是地层的等效品质因子，单套地层的品质因子可通过式（8）求取[12]，即式中：t1、t2分别为地震波通过目的层顶、底界面时所用的传播时间；分别为t1、t2时刻的等效品质因子；Q2为t1与t2之间的品质因子。 在对地震数据处理时，可以通过以下步骤从浅到深估算每个时间采样点的等效品质因子，实现单套地层品质因子的转换。 （1）对地震数据逐道进行S变换，得到时频域能量分布。 （2）假设地震数据从浅到深时间采样点分别为t1,t2，…，tn-1,tn（n为正整数），根据时频域能量分布和中心频率的定义，计算拾取t1,t2，…，tn-1,tn对应的中心频率，并利用式（7）估算出t1,t2，…，tn-1,tn对应的等效品质因子。 （3）将t1,t2时刻等效品质因子代入式（8），计算得到t1,t2之间地层品质因子，依次类推，将tn-1,tn时刻等效品质因子代入式（8）,计算得到tn-1,tn之间地层品质因子，并按此过程依次得到两两时刻之间的单套地层品质因子。 （4）在地层品质因子参数剖面上，记录t1,t2，…，tn-1,tn时间采样点品质因子为对应的单套地层品质因子值，如：tn时刻品质因子记为tn-1,tn之间地层品质因子。 通常情况下，地层含气与不含气所引起的地震波吸收衰减特征存在明显差异。为了研究地层含气与不含气引起的地震波吸收衰减变化特征，并验证中心频率法估算的地层Q值检测气层的可行性，设计一个三层水平层状模型，其地层参数如图1所示，其中上覆及下伏地层的Q值均为300，中间砂岩地层含气部位的Q值从5开始以增量10逐渐增大至155，含水部位的Q值为300。 利用单程波法地震正演模拟制作反射合成记录。图2为含气层Q值为5时对应的正演模拟记录及吸收衰减分析。现选取气层