本项目旨在利用MATLAB开发环境，实现一种高精度的地震子波提取工具，主要解决地震资料数字处理中由于子波未知或时变导致的盲反褶积难题。传统的基于二阶统计量（如自相关函数）的方法通常只能提取最小相位子波，无法保留相位信息。本项目利用高阶累积量（主要是四阶累积量）对高斯噪声不敏感且保留系统相位信息的特性，构建基于各阶累积量匹配的目标函数。为了解决该目标函数高度非线性且存在多个局部极值的问题，项目采用模拟退火算法（Simulated Annealing）作为核心优化引擎。该系统具备以下详细功能：1. 数据预处理模块，负责读取SEGY或MAT格式地震数据并进行去噪；2. 统计量计算模块，能够准确估计地震记录的二阶及四阶累积量；3. 模拟退火优化模块，通过模拟热力学退火过程（包含温度初始化、扰动产生新解、Metropolis接受准则、降温策略等），在参数空间内全局搜索最优子波参数，有效避免陷入局部最优解；4. 子波重构与反褶积模块，根据优化结果重构混合相位子波，并进行维纳滤波或反滤波操作以获取反射系数序列；5. 结果可视化模块，提供时域波形对比、频域谱分析以及收敛曲线展示，验证提取结果的准确性与鲁棒性。
