一种类RNN的改进ISTA稀疏脉冲反褶积

稀疏脉冲反褶积方法对提高地震资料分辨率有着重要作用,迭代阈值收缩算法(ISTA)是其核心算法,首先利用地震数据提取子波,再利用ISTA求解反射系数。当地震子波提取不准确时,反褶积效果不理想。为此,在ISTA基础上,结合循环神经网络(RNN)中反向传播(BPTT)的思想,研究形成了一种类RNN的改进ISTA稀疏脉冲反褶积方法。该算法首先使用常规手段从实际地震数据中提取地震子波,构建反褶积的子波字典;然后将构建的地震子波字典作为已知的初始条件,结合ISTA求取的反射系数;再根据BPTT算法思想,将求取的反射系数与子波褶积并与实际数据进行比较,反向修改地震子波;最终,经过多次迭代修改获得合理的地震子波字典,并利用该地震子波字典求解实际地震数据的反射系数序列。为验证算法的有效性,采用不同信噪比的理论地震记录,给定存在较大误差的初始子波,进行了反褶积计算。采用传统的ISTA和类RNN的改进ISTA进行对比处理,结果表明,改进ISTA具有较好的抗噪能力和子波自适应能力,可使实测地震资料的有效频带拓展约1.5倍,能够较好地适应实际地震资料的反褶积处理。     

基于地质统计学的沙漠地区逐点时深曲线静校正方法

在我国西部沙漠地区,起伏剧烈的沙丘造成了严重的静校正问题。在沙丘成因单一、物性变化不大的沙漠区,通常应用时深曲线静校正方法来解决静校正问题,但在沙丘成因多样化或者沙漠边缘区域,应用常规的时深曲线静校正技术误差较大。基于微测井测量数据的地质统计学分析结果,提出了一种逐点时深曲线静校正方法。基于数据正态分布检验、变异函数分析等地质统计学手段,明确了沙漠地区沙层厚度与地震波垂直传播到地表的时间之间具有空间分布相关性,在此基础上利用克里金插值方法对微测井测量数据进行插值,得到各炮点和检波点处的时间-深度散点数据,再利用最小二乘法对所有散点进行拟合,获得炮点和检波点的时深曲线量板,计算静校正量。以准噶尔盆地S1X工区野外64口微测井的测量数据为控制点,采用上述方法计算得到各个炮点和检波点的静校正量,并对计算结果进行了交叉验证,该方法获得的静校正量与微测井获得的静校正量基本一致。应用该静校正量进行了叠加成像,并与常规静校正方法处理后的共炮点道集和叠加剖面进行了对比,结果显示,逐点时深曲线静校正方法的应用效果明显优于常规时深曲线静校正方法。