【自适应时频优化+瞬时相位梯度】基于熵准则自适应加窗时频分析和群延迟过零检测的地震反射层智能拾取方法（MATLAB）

算法提出了一种基于自适应窗口优化和瞬时相位梯度分析的地震反射层自动拾取方法。首先通过计算高斯窗宽度随频率变化的短时傅里叶变换（STFT），并依据频率域能量分布的最小熵准则自适应选择最优窗宽，构建高分辨率的时频表示；进而利用时频域能量分布估计瞬时中心频率与带宽，同时通过解析信号的瞬时相位导数计算旅行时（群延迟）参数；最后基于旅行时与观测时间的零交叉检测，自动识别反射层位置并输出二值指示序列。该方法结合了自适应时频分析、瞬时属性提取和过零检测技术，适用于非平稳信号的层位自动追踪。

算法步骤

信号生成与预处理：设定时间采样间隔、记录时长和频率参数，构造不同中心频率的Ricker子波，并与随机稀疏反射系数序列进行卷积合成地震道信号。

自适应时频变换：对输入信号进行填充并构造Hankel矩阵；遍历多个高斯窗形状参数（α），对每一列加窗后进行快速傅里叶变换（FFT），计算频率域对数能量之和作为代价函数；选择使代价函数最小的α作为该时刻的最优窗宽，构建自适应加窗的时频矩阵。

瞬时带宽估计：对时频矩阵取幅度并提升幂次以增强特征，计算每一时刻的频率加权中心（一阶矩）和频率标准差（二阶矩的开方），得到瞬时频带范围。

旅行时（群延迟）计算：对时频矩阵进行逆傅里叶变换得到时域解析信号，乘以时间对角阵后再做FFT得到时间加权频谱；利用频谱的实部与虚部及其导数的组合关系，计算瞬时相位的时间导数，即群延迟（旅行时）。

零交叉检测与层位拾取：将旅行时与观测时间作差，对差值序列进行过零检测，标记零交叉点位置，并剔除边界假象，最终输出反射层位置的二值指示序列。

结果可视化：绘制时频能量谱、旅行时剖面及零交叉点指示图，并添加相应的图题与坐标标注。