杨怀杰,何中波,吴仙明

(核工业北京地质研究院,北京 100029)

精确的地层划分和岩性特征分析结果是岩石样品在实验室中通过化验分析得到的,不仅价格昂贵而且耗费时间,同时还需依靠经验丰富的地质学家。低成本且有效的地层划分和岩性特征分析,一般要通过测井数据综合解释来实现,因为测井数据具有高采集率和包含地下不同深度域岩石信息的特征。传统的利用测井数据分析岩性特征的方法不仅依靠经验丰富的解释人员,而且可能产生多种解释结论。为了排除主观因素的影响,国内外学者提出了许多方法,主要是多元统计分析和人工智能方法[1]。虽然这些方法能够自动地识别岩性,但是在参数的选择上,不仅需要经验丰富的解释人员,而且需要大量的测井曲线数据。

国内外的地球物理学家将信号处理技术应用于测井数据处理解释中,包括傅里叶变换[2]、沃尔什变换[3]和小波变换,特别是小波变换被广泛的应用于地球物理数据的处理解释中[4--7]。本文以黑龙江三江盆地ZKQX2 钻孔为例,对三侧向测井数据进行连续小波变换,得到不同尺度下的小波系数,绘制小波系数尺度图和高频小波系数曲线图,来识别地层的岩性界面。

1 小波变换原理

小波变换的核函数定义为[8--9]:

式中:ψ(t)称为小波函数;a为尺度因子;b为位移因子。

当尺度因子和位移因子连续变化时,称为连续小波变换,公式如下[10]:

式中:WT f(a,b)为连续小波变换系数;f(t)为待处理信号。……