[发明专利]一种基于谱分解的煤层火成岩侵入位置直接指示方法

权利要求书

1.一种基于谱分解的煤层火成岩侵入位置直接指示方法，其特征是：火成岩位置直接指示方法：1)通过人机交互，在常规三维地震数据体上拾取煤层反射波同相轴；2)根据煤层反射波的特征，沿煤层反射波同相轴提取煤层反射波数据，生成子数据体；3)利用谱分解算法，对子数据体中的每一个地震道进行谱分解，获得相应的时频谱；4)对于每一个时频谱，搜索其各个频率曲线的最大响应值，将最大响应值对应的时间和频率绘制成交会图；5)利用线性回归分析方法，拟合交会图中的散点，获得截距和梯度值；6)将所有地震道的截距和梯度值绘制成交会图，根据每一个地震道在交会图中的位置，直接指示火成岩的侵入位置。

2.根据权利要求1所述的一种基于谱分解的煤层火成岩侵入位置直接指示方法，其特征是：所述的火成岩位置直接指示方法，具体步骤如下：

1)通过人机交互，在常规三维地震数据体上拾取煤层反射波同相轴；

具体方法如下：①收集采区内已有的测井资料，制作合成地震记录；②将合成地震记录与井旁道进行对比，识别出煤层反射波的位置和相位；③将识别出的煤层反射波位置和相位作为控制点，利用人机交互在全区拾取煤层反射波同相轴；

2)根据煤层反射波的特征，沿煤层反射波同相轴提取煤层反射波数据，生成子数据体；

具体方法如下：①根据煤层反射波的特征，确定数据提取时的时间窗上界、时间窗下界和时间窗长度；②将提取的时间窗函数确定为钟形高斯窗；③沿交互拾取的煤层反射波同相轴，依次提取每个地震道的煤层反射波数据，生成子数据体；

3)利用谱分解算法，对子数据体中的每一个地震道进行谱分解，获得相应的时频谱；

具体方法如下：①提取子数据体中每一道的煤层反射波数据；②根据煤层反射波的特征，合理选择时频谱的起始频率、终止频率和步长等参数；③选取合适的谱分解算法，为小波变换或S变换等；④计算数据体中每一个地震道的时频谱；

4)对于每一个时频谱，搜索其各个频率曲线的最大响应值，将最大响应值对应的时间和频率绘制成交会图；

具体方法如下：①搜索时频谱中各个频率曲线的最大响应值；②读取各个频率曲线最大响应值对应的时间；③以频率为横坐标、最大响应值对应的时间为纵坐标，绘制交会图；

5)利用线性回归分析方法，拟合交会图中的散点，计算截距和梯度值；

具体方法如下：①根据交会图中散点分布的特点，选取线性回归的数据范围，起始频率选为20Hz，终止频率选为100Hz；②根据最小二乘法原理，对数据范围内的散点进行线性回归；③根据线性回归所获得的直线方程，计算拟合直线的截距和梯度值；

6)将所有地震道的截距和梯度值绘制成交会图，根据每一个地震道在交会图中的位置，直接指示火成岩的侵入位置；

具体方法如下：①将所有地震道所对应的截距和梯度依次排序；②绘制截距和梯度的交会图；③根据交会图中散点的分布特点，定义阈值线，初步指示火成岩侵入位置；④将初步指示结果与已知钻孔揭露情况进行对比，验证定义的阈值线是否合理；⑤如果合理，则初步指示火成岩位置为煤层火成岩侵入位置的直接指示；如果不合理，则重新定义阈值线，重复步骤④直到获得满意结果为止。