基于广义S变换和JADE算法的储层识别.pdfVIP

第 5O卷第4期 石 油 物 探 Vo1．50，No．4 2011年 7月 GE0PHYSICALPROSPECTING F0RPETR0LEUM Ju1．，2011 文章编号 ：1000—1441(2011)04～0319—05 基于广义 S变换和JADE算法的储层识别 胥德平 ，郭 科 ，文晓涛。，花小琴 (1．数学地质四川省重点实验室，四川成都 610059；2．成都理工大学管理科学学院，四川成都 610059；3．成都理工大学地球物理学院，．四川成都 610059；4．江苏科技大学数理学院，江苏镇江 212003) 摘要：广义 s变换具有良好的时频聚集性 ，能够获得反映储层性质的时频谱 ，但 由此生成许多单频数据体，造成 地震数据的分析和解释工作变得繁琐。为此，将广义s变换和特征矩阵联合近似对角化(JADE)算法相结合，进 行储层识别。在时间域内对地震信号做广义S变换 ，然后在时频域内提取一些相互独立的有效频率分量，最后 根据已知井的储层发育情况，选取对储层识别有效的频段分量进行储层预测。仿真和实际地震数据应用结果表 明了该方法的有效性。 关键词：广义s变换；联合对角化；时频分析；储层识别 DOI：10．3969／j．issn．1000—1441．2011．04．001 中图分类号：P631．4 文献标识码：A 目前分析地震信号的方法很多，但都离不开傅 下，将由若干个独立源信号混合而成的混合信号分 里叶变换。傅里叶变换是全局变换，可以有效分析 解成各个独立信号。特征矩阵的联合近似对角化 平稳信号，但采用傅里叶变换分析非平稳信号缺乏 法 (joint approximative diagonalization of 物理意义。信号的局部性能需要使用时域和频域 eigenmatrix，JADE)是Cardoso提出的一种较为稳 的二维联合表示才能得到精准的描述。人们利用 健的独立分量分析算法_1¨]。该算法引入多变量 预先加窗的方法使频谱能反映时间局部特性，对傅 数据的四维累计量矩阵，利用Jacobi对角化方法 里叶变换进行推广和改造，产生一系列可对信号局 来寻找独立分量信号。 部特性进行分析的时频分析方法口]，如：短时傅里 地震信号是有一定带宽的信号[1引，在这个带 叶变换 (Shot—TimeFourierTransform，STFT)、 宽内，有效信号的频率反映各种地质信息(如油气、 Gabor变换、Wigner—Ville时频分布、小波变换 储层等)，可能某一频段的信号对油气层较为敏感， (WaveletTransform，WT)等。这些时频分析方法 而其它频段的信号对油气层不够敏感。广义 s变 各有特点，但都是基于信号在窗函数的有效持续时 换具有 良好的时频聚集性，利用地震信号的频率灵 间内是平稳的假设，因而存在一定的局限。 活选取窗函数，能够获得较为理想的时频分辨率， 1996年，Stockwell等[23提出了S变换，它综 但由此会生成许多单频数据体，使得地震数据的分 合了短时傅里叶变换和小波变换的优点，采用尺度 析和解释工作变得繁琐。我们通过对地震数据做 可以变化的局部高斯窗函数，其时频分辨率随频率 广义 S变换，利用JADE算法，在时频域 内从不同 发生变化，具有多尺度聚焦性，而且基本小波不必 的单频数据体中提取出对油气敏感的共同信息，从 满足容许性条件，这些特点在实际应用中非常有 而有效预测含油气储层。 用。然而，由于S变换的基本小波固定，因而使其 在应用中受到限制。为此，很多学者在S变换的基 1 广义S变换原理 础上，对基本小