用高分辨率双曲线Radon变换实现波场分离

【摘要】：多波多分量地震勘探比传统的单一纵波勘探具有许多优势,同时采集的纵、横波资料可以分别形成叠加剖面,有利于构造成像。更可贵的是它可以用来提取更多的岩性、储层、裂隙分布的地球物理参数以及应力场信息。波场分离是多波多分量地震资料处理的重要环节之一。 Radon变换已广泛应用于地震资料的处理和解释中,在地震同相轴识别、波场分离、压制多次波、速度分析等方面具有良好的效果。无论是线性Radon变换、抛物线Radon变换,还是双曲线Radon变换,以往都采用最小二乘反演方法实现。由于采用最小二乘法在分辨率和计算效率方面较低,因此目前在这一方面仍值得进一步研究。 本文重点讨论了影响Radon变换分辨率的因素,以及解决这一因素的办法,给出了高分辨率Radon变换的实现方法。选取了一种分辨率较高的双曲线Radon变换方法,为了克服其计算速度慢的缺点,采用共轭梯度法来求解稀疏矩阵,这样既不影响分辨率又提高了计算速度。根据纵波和转换波的传播速度和截距时间的不同,采用高分辨率双曲线Radon变换法实现波场分离,并分析了变换中参数值的选取以及端点效应现象,比如,频率越高,道间距越大,越容易出现空间假频,则分辨率越低。文中给出了模型及应用实例,研究表明,采用共轭梯度法求解的高分辨率双曲线Radon变换算法的分辨率和计算效率都较高,可以更好地实现波场分离。