[发明专利]工区地下介质频率域异常分析的方法

权利要求书

1.一种工区地下介质频率域异常分析的方法，其特征在于，包括以下步骤：A、通过测井资料的岩石物理分析，区分工区的地下介质的砂岩地层和泥岩地层；B、通过时频分析算法，求取工区的地震资料在砂岩地层和泥岩地层的频率信息资料；C、通过地震波衰减特征提取技术，对频率信息资料进行衰减特征提取，获得地震资料在砂岩地层和泥岩地层的频率衰减特征参数，并分别统计砂岩地层的频率衰减值和泥岩地层的频率衰减值；D、利用最小二乘法对砂岩地层的频率衰减值进行线性拟合或多项式拟合，得到砂岩地层的频率衰减趋势线；E、利用砂岩地层的频率衰减趋势线，分析工区地下介质频率域的异常。

2. 根据权利要求1所述的工区地下介质频率域异常分析的方法，其特征在于：所述步骤D中，利用最小二乘法对泥岩地层的频率衰减值进行线性拟合或多项式拟合，得到泥岩地层的频率衰减趋势线；所述步骤E中，以砂岩地层的频率衰减趋势线为主，以泥岩地层的频率衰减趋势线为附，分析工区地下介质频率域的异常。

3.根据权利要求1或2所述的工区地下介质频率域异常分析的方法，其特征在于：所述步骤A中，应用测井约束地震反演技术，反演得到高分辨率的波阻抗数据体，对测井资料进行岩石物理分析，分析阻抗对于砂岩、泥岩区分度的敏感性，结合岩石物理分析的结果，确定合适的纵波阻抗的阈值，通过纵波阻抗的阈值，区分工区的地下介质的砂岩地层和泥岩地层。

4.根据权利要求1或2所述的工区地下介质频率域异常分析的方法，其特征在于：所述时频分析算法，采用短时傅立叶变换算法、小波变换算法、S变换算法或者匹配追踪时频分析算法。

5.根据权利要求4所述的工区地下介质频率域异常分析的方法，其特征在于：所述匹配追踪时频分析算法，具体的算法流程如下:

用Hilbert 变换计算复地震道；

计算复地震道的瞬时包络、瞬时相位和瞬时频率；

找到包络最大值及其对应的时间位置，得到最优小波振幅的估计值；

计算相应时间处的瞬时相位和瞬时频率，得到最优小波相位角和主频的估计值；

从原地震道中减去上面确定的小波的实部；

重复至步，直到剩余的地震道幅度小于迭代误差。

6. 根据权利要求1或2所述的工区地下介质频率域异常分析的方法，其特征在于：所述衰减趋势线拟合的具体算法如下： 

分别统计砂岩和泥岩的衰减特征，建立衰减特征与深度的散点函数关系：

  ，其中，是衰减值，是深度，是样本数；

应用最小二乘法，对散点数据进行数值拟合，分别拟合出砂岩地层和泥岩地层的衰减趋势线，求取多项式（）使得

                   

满足上式的多项式称为最小二乘拟合多项式，特别的当时，称为线性拟合。