分形插值地震数据重建方法研究

对分形插值方法作了较详细的探讨,给出了分形插值函数的显式表达方式,同时给出了垂直比例因子的局部显式表达式,旨在提高地震道插值重建的精度及突出局部信息,并从单道地震图的角度分析其在地震道插值重建中的应用效果.利用该方法对理论模型和济阳坳陷实际地震台站资料进行了重建处理,结果表明,分形插值重建的地震道是原始地震道的良好近似,缺失道的振幅和相位都得到了很好的恢复.该法克服了随机分形插值方法必须进行多步迭代的弱点,提高了计算效率.通过对单道地震图插值重建结果的分析,说明了本文分形插值方法具有较高的精度和较高的效率,有深入研究的潜力.本文提出的显式分形插值方法既能够突出地震道数据的局部信息,又较好地保持了地震道数据的总体变化趋势.     

基于显式垂直比例因子的显式分形插值地震数据重建

在前人工作的基础上对分形插值方法作了详细的探讨,给出了分形插值函数的显式表达方式.在量纲分析的基础上给出了垂直比例因子的局部显式表达式,旨在提高地震道插值重建的精度及突出局部信息,并从单道地震图的角度分析其在地震道插值重建中的应用效果.研究了垂直比例因子的变化对分形插值精度的影响.数值实验表明,随着垂直比例因子的增大,分形垂直的误差逐渐增大,二者之间呈显出指数增长的趋势.该法克服了随机分形插值方法必须进行多步迭代的弱点,提高了计算效率.通过对理论地震道插值重建的分析,说明了本文分形插值方法的高精度和高效率.本文提出的显式分形插值方法既能够突出地震道数据的局部信息,又较好地保持了地震道数据的总体变化趋势.     

不规则地震道数据规则化重建方法研究

不规则地震数据会对地震多道处理技术的正确运行产生不良影响,降低地震资料处理质量.本文依据不规则地震数据的表现特征将其划分为四种类型并针对第三类不规则地震道数据采用抗泄露Fourier变换方法进行规则化重建.不规则采样数据会破坏Fourier基函数的正交性并产生频谱泄漏现象.抗泄露Fourier变换方法通过递归相减来压制...     

基于jitter采样和曲波变换的三维地震数据重建

传统的地震勘探数据采样必须遵循奈奎斯特采样定理,而野外数据采样可能由于地震道缺失或者勘探成本限制,不一定满足采样定理要求,因此存在数据重建问题.本文基于压缩感知理论,利用随机欠采样方法将传统规则欠采样所带来的互相干假频转化成较低幅度的不相干噪声,从而将数据重建问题转为更简单的去噪问题.在数据重建过程中引入凸集投影算法(POCS),提出采用e-√x(0≤x≤1)衰减规律的阈值参数,构建基于曲波变换三维地震数据重建技术.同时针对随机采样的不足,引入jitter采样方式,在保持随机采样优点的同时控制采样间隔.数值试验表明,基于曲波变换的重建效果优于傅里叶变换,jitter欠采样的重建效果优于随机欠采样,最后将该技术应用于实际地震勘探资料,获得较好的应用效果.     

基于泊松碟采样的地震数据压缩重建

在地震资料处理领域,数据的压缩和重建是非常重要的问题,但往往由于数据的严重缺失或采样原因而达不到理想的效果.新发展起来的压缩感知理论为重建欠采样数据提供了可能,而选用合适的采样方法是其中的关键技术之一.本文基于傅里叶变换和压缩感知理论,采用泊松碟采样,对不完整地震数据进行恢复重建.数值实验表明,与传统的单纯随机采样方法相比,泊松碟采样方法在保持采样随机性的同时,使采样点的分布更加均匀,有效地调节了采样间距,从而达到更好的恢复效果,可以有效地指导地震数据采集设计及重建.     

地震数据重建方法综述

地震勘探的目的是为了获得地下构造的精确成像.由于人为因素和环境原因,地震数据在空间方向上往往是不规则采样或缺失采样的,这可能会在数据处理时产生假象,最终导致错误的解释.因此经常需要在空间方向对缺失的地震数据进行重建.重建问题可以看作是一个反演问题,即从不完整的地震数据中重建出完整的地震波场.本文重点研究了国内外比较成熟的地震数据重建方法,如基于滤波的重建方法、基于波场延拓算子的重建方法、基于变换域的重建方法以及相干倾角插值等方法,分析了这些方法的优缺点,以及目前地震数据重建所面临的挑战.     

地震丛集的分形新方法—物理分形

地震时间丛集的传统分形分析方法把地震事件看作是时间轴上无质量的数学点。这些方法得到的结果主要反映了大量小地震的活动特征。本文发展了一个新的物理分形分析的筛子方法，采用一系列震级筛子来研究不同震级地震的分布特征。对华北和南加洲的研究结果表明，同一地区地震目录中不同震级范围的子集合具有近似相等的分维数，华北Ｄ＝０．３，南加洲Ｄ＝０．４。     

饱和砂土地震液化判别的分形插值模型

借鉴分形基本理论,提出了基于分形插值模型的饱和砂土地震液化判别方法.该方法首先选取影响饱和砂土地震液化判别的7个主要因素,根据分类标准,采用在每级标准中随机内插的方法,得到40个标准样本,用于构建饱和砂土地震液化判别的分形插值模型;其次根据最大似然分类原则确定每个饱和砂土地震液化判别指标的评价分维数;然后利用加权求和法计算样本的综合评价值,并根据样本综合评价值与经验等级之间的关系建立分形插值评价模型;最后,进行了实例分析结果表明:该模型的评价结果合理、客观,计算得到的每个样本具体得分值,即使对属于同一级的样本也可以给出其地震液化程度的顺序,为饱和砂土地震液化评价工作提供了一种新的研究方法与思路.     

最小加权范数插值与调制升频结合的地震数据重建研究

提出一种新的反假频地震数据重建的两步算法,将最小加权范数插值(MWNI)方法与调制升频方法有效地结合起来.首先利用MWNI方法构建数据谱的低频部分,为了提高计算效率,在低频重建算法中,引入了预条件共轭梯度法求解反问题方程,并使用了与频率有关的变波数带宽技术;然后,基于重建的低频数据,采用调制升频方法重构数据的高频部分.调制升频方法灵活,简便,能有效地从低频资料中恢复出高频成分,克服了以往AR模型预测高频走不远的限制,当数据存在严重的空间假频时,亦能获得较好的重建效果.该两步算法不仅可用于规则地震数据的内插重建,也可用于含空道地震数据的重建.理论模型和实际地震数据重建试验表明,该方法效率高,精度高,反假频能力强,重建剖面波形连续、自然,与正确完整的地震剖面相似程度高,具有良好的实用价值和应用前景.     

地震数据Bregman整形迭代插值方法

人工地震方法由于受到野外观测系统和经济因素等的限制,采集的数据在空间方向总是不规则分布.但是,许多地震数据处理技术的应用（如：多次波衰减,偏移和时移地震）都基于空间规则分布条件下的地震数据体.因此,数据插值技术是地震数据处理流程中关键环节之一.失败的插值方法往往会引入虚假信息,给后续处理环节带来严重的影响.迭代插值方法是目前广泛应用的地震数据重建思路,但是常规的迭代插值方法往往很难保证插值精度,并且迭代收敛速度较慢,尤其存在随机噪声的情况下,插值地震道与原始地震道之间存在较大的信噪比差异.因此开发快速的、有效的迭代数据插值方法具有重要的工业价值.本文将地震数据插值归纳为数学基追踪问题,在压缩感知理论框架下,提出新的非线性Bregman整形迭代算法来求解约束最小化问题,同时在迭代过程中提出两种匹配的迭代控制准则,通过有效的稀疏变换对缺失数据进行重建.通过理论模型和实际数据测试本文方法,并且与常规迭代插值算法进行比较,结果表明Bregman整形迭代插值方法能够更加有效地恢复含有随机噪声的缺失地震信息.

     

显式分形插值在有限元叠前逆时偏移成像中的应用

以分形理论为基础,借助于压缩映射原理和不动点理论及实变函数理论,给出了分形插值函数的显式表达形式,同时给出了垂直比例因子的局部显式表达式.利用该显式分形插值方法对稀疏采样的地震记录进行了插值加密处理,并对原始记录和分形插值重建记录进行了叠前逆时偏移处理.数值重建实验表明,显式分形插值重建的剖面与原始地震剖面非常类似,单道地震记录则表明显式分形插值重建的地震道是原始地震道的良好近似,缺失道的振幅和相位都得到了很好的估计.叠前偏移结果则表明,分形插值叠前偏移剖面是原始记录偏移剖面的良好近似,二者的分辨率几乎一样.     

地震前兆观测资料的分形处理方法

本文从自仿射分形的角度探讨地震兆汾观测资料的分形处理方法。在对自仿射分形的分维量测方法进行对比分析的基础上，给出了处理地震前兆观测资料较为适用的分形方法和满足计算需要的样本点数，从对部分水氡和地电阻率资料进行了处理分析中，得到了有益的启示：地震前兆时序观测资料可用分形方法处理分析，分维值可反映观测资料的结构变化特征，其作为一个无量纲量为综合对比分析各类前兆现象，探讨地震前兆场的复杂性特征提供了有效     

基于非均匀Fourier变换的地震数据重建方法研究

不规则采样地震数据会对地震数据的多道处理造成严重影响,将非均匀Fourier变换和贝叶斯参数反演方法相结合,对不规则空间带限地震数据进行反演重建.对每一个频率依据最小视速度确定出重建数据的带宽,然后从不规则地震数据中估计出重建数据的空间Fourier系数.将不规则地震数据重建视为信息重建的地球物理反演问题,运用贝叶斯参数反演理论来估计Fourier系数.在反演求解时,使用共轭梯度算法,以保证求解的稳定性,加快解的收敛速度.理论模型和实际资料处理验证了本方法的有效性和实用性.     

澜沧地震的分形研究

从现代物理学观点看地震成因,在漫长的孕震过程中,这种大的相互作用系统总是朝着某种临界状态不断地自组织(Bak P.et al,1989)。在临界状态下,小事件引起的链锁反应很可能导致大系统失稳,产生大地震。某些大震前中、小地震活动密度增大且速率增强,很可能说明了这一点。     

分割地震剖面的分形特征研究

本文针对地球物理中地震剖面图象的分割,进行特征提取和计算的探讨．分析了纹理结构的自相似性特性,应用分形原理进行分形维数的估测,探索了用分形特征进行图象分割的可行性．     

中国地震分形研究的进展

80年代中期以来地震分形研究在中国兴起,并发展为一个十分活跃的科学分支。本文综述了近年中国地震分形研究的概况,及取得的重要进展,并提出了急需解决的主要问题。     

连续缺失地震数据的高阶流式预测滤波插值方法

由于观测系统实施以及经济因素的限制,采集到的勘探地震数据在空间方向上总是不规则分布的,并且往往会出现数据的大范围连续缺失情况.许多后续地震数据处理方法（例如：多次波压制和波动方程偏移等）都需要空间上规则分布的数据.插值技术是一种解决地震数据缺失问题的有效手段,但是传统的数据插值方法在进行连续缺失数据重建时往往会出现失效的情况,尤其在处理非平稳地震同相轴时精度不高,并且大多数的方法需要迭代计算,在处理高维大规模数据时效率较低.针对连续缺失地震数据的快速插值问题,本文提出了一种非迭代的时空域高阶流式预测滤波插值方法,通过使用高阶限制条件来提高连续缺失数据的滤波器估计精度,提高局部约束条件的稳定性,改善低阶流式计算由于滤波器系数无法连续更新所造成的插值失效情况.同时,空间非因果滤波器和蛇形插值路径的设计方案可以有效减小大范围连续缺失数据和数据边界对于预测滤波器的计算误差,本方法能够有效处理包括近炮检距缺失情况在内的连续缺失数据插值重建.通过与工业标准傅里叶凸集投影（POCS）方法进行比较,理论模型和实际数据处理结果表明,本文提出的高阶流式预测滤波插值方法对高维连续缺失地震数据有较好的重建效果,在插值精度和计算效率两个方面有更好的平衡性.

     

基于压缩感知的Curvelet域联合迭代地震数据重建

由于野外采集环境的限制,常常无法采集得到完整规则的野外地震数据,为了后续地震处理、解释工作的顺利进行,地震数据重建工作被广泛的研究.自压缩感知理论的提出,相继出现了基于该理论的多种迭代阈值方法,如CRSI方法（Curvelet Recovery by Sparsity-promoting Inversion method）、Bregman迭代阈值算法（the linearized Bregman method）等.CSRI方法利用地震波形在Curvelet的稀疏特性,通过一种基于最速下降的迭代算法在Curvelet变换域恢复出高信噪比地震数据,该迭代算法稳定,收敛,但其收敛速度慢.Bregman迭代阈值法与CRSI最大区别在于每次迭代时把上一次恢复结果中的阈值前所有能量都保留到本次恢复结果中,从而加快了收敛速度,但随着迭代的进行重构数据中噪声干扰越来越严重,导致最终恢复出的数据信噪比低.综合两种经典方法的优缺点,本文构造了一种新的联合迭代算法框架,在每次迭代中将CRSI和Bregman的恢复量加权并同时加回本次迭代结果中,从而加快了迭代初期的收敛速度,又避免了迭代后期噪声干扰的影响.合成数据和实际数据试算结果表明,我们提出的新方法不仅迭代快速收敛稳定,且能得到高信噪比的重建结果.