本发明属于石油天然气勘探技术领域，涉及一种基于深度域地震波频散分析的储层含气性评价方法。
背景技术：
现有基于地震波频散分析的储层含气性评价方法技术在时间域进行(时频分析)，因为地质介质地震波传播速度的可变性，使得识别的含气层的空间位置确定存在不确定性，同时也无法确定体效应。
技术实现要素：
本发明的目的在于提供一种基于深度域地震波频散分析的储层含气性评价方法。
为实现上述技术目的，达到上述技术效果，其技术方案具体为：
包括以下步骤；
步骤1.将时域地震记录转换为深度域地震记录；利用测井曲线进行层位划分并建立深-时关系，根据多层介质反射系数的波动方程推导公式，求取频率域的反射系数；对时间域的井旁道记录利用同态反褶积方法提取地震子波，得到频率域的子波谱；利用频率域褶积模型计算合成记录，作反傅里叶变换得到时间域合成地震记录，再根据时-深关系求出每个时域采样点对应的深度点，按照所要求的深度采样间隔进行重采样，得到深度域合成地震记录；
步骤2.在深度域进行深频分析，获得频率成分随深度的变化；
深度域地震记录进行预处理。测井曲线的预处理包括去奇异值、环境校正和重采样等，对预处理后的数据(信号)在深度域进行小波变换和频谱特征分析，获得频率成分随深度的变化；
步骤3.计算各深度的频散，识别含气层。
进一步，深时和时深转换在求取地震子波的过程中，需要把井旁道记录从深度域转换到时间域，可以用离散曲线的连续化公式，求取对应测井曲线深度点的记录值，然后重采样得到等时间隔的记录；或者根据每个测井曲线的分层层速度，求逐点的时间值，然后重采样就得到等时间隔的记录。
进一步，步骤1中，所述测井曲线为声波、密度测井曲线。
进一步，步骤1中，所述测井曲线为等间隔采样。
本发明具有以下有益效果：本发明在深度域进行地震波的频散分析，则能直接确定含气层的空间深度位置，并能定量估算储层的体效应。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。
一种基于深度域地震波频散分析的储层含气性评价方法，包括以下步骤：
步骤1.将时域地震记录转换为深度域地震记录；利用测井曲线进行层位划分并建立深-时关系，根据多层介质反射系数的波动方程推导公式，求取频率域的反射系数；对时间域的井旁道记录利用同态反褶积方法提取地震子波，得到频率域的子波谱；利用频率域褶积模型计算合成记录，作反傅里叶变换得到时间域合成地震记录，再根据时-深关系求出每个时域采样点对应的深度点，按照所要求的深度采样间隔进行重采样，得到深度域合成地震记录；
步骤2.在深度域进行深频分析，获得频率成分随深度的变化；
深度域地震记录进行预处理。测井曲线的预处理包括去奇异值、环境校正和重采样等，对预处理后的数据(信号)在深度域进行小波变换和频谱特征分析，获得频率成分随深度的变化；
步骤3.计算各深度的频散，识别含气层。
深时和时深转换在求取地震子波的过程中，需要把井旁道记录从深度域转换到时间域，可以用离散曲线的连续化公式，求取对应测井曲线深度点的记录值，然后重采样得到等时间隔的记录；或者根据每个测井曲线的分层层速度，求逐点的时间值，然后重采样就得到等时间隔的记录。
步骤1中，所述测井曲线为声波、密度测井曲线。
步骤1中，所述测井曲线为等间隔采样。
以上所述，仅为本发明较佳的具体实施方式，本发明的保护范围不限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。