简介：野外空间采样密度的提高将增加室内数据分析的工作量，常规的基于地震数据的点与线分析方法有一定的局限性。本文简要说明空间子集的抽取方法，列举了正交子集和斜交子集的特点，并通过三维可视化展示了子集数据的空间特性。提出在数据处理中利用子集的时间切片：（1）分析地震道空间分布的均匀性和规则性；（2）研究面波及规则干扰的空间分布特点；（3）检测叠前数据中的异常信息；（4）监控叠前去噪的效果。实际数据的应用结果表明基于空间子集的分析方法是一种独特有效的地震数据分析方法，从另一种视角观察地震数据，可以发现其中某些新的特征，以提高处理人员对数据的洞察能力。

简介：针对有限差分数值模拟的频散问题，本文将交错网格技术和紧致差分格式相结合，推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式；对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项，并利用Fourier误差分析方法分析了上述三种差分格式的近似精度；在此基础上，分别采用上述三种差分格式进行了波场数值模拟。结果表明，当差分方程阶数相同时，紧致交错网格差分格式截断误差最小，数值频散最弱，差分精度最高，证实了该方法的有效性。

简介：角度域共成像点道集（ADCIGS）是偏移速度分析和振幅随角度变化分析（AVA）的基础数据。传统Kirchhoff叠前深度偏移（KPSDM）按偏移距组织数据,能方便的输出偏移距域共成像点道集（ODCIGS）,其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为：（1）利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;（2）根据旅行时场的梯度方向计算反射张角;（3）在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是：（1）按照炮数据来组织输入数据;（2）旅行时场的输入与单炮覆盖范围相联系以节省内存;（3）多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：传统的f-x域经验模态分解法（Empiricalmodedecomposition,EMD）能够有效地对主要由水平同相轴构成的地震记录进行随机噪声衰减。然而,当同相轴倾斜时,f-x域经验模态分解法在衰减随机噪声的同时去除大部分有效信号。本文提出了一种基于f-x域经验模态分解法的改进算法。我们通过局部相似度对所去除的噪声信号中的有效信号进行提取。局部相似度可以用来检测噪声信号中的有效信号点并用来构造一权重算子进行信号提取。新方法与f-x域经验模态分解法、f-x域预测滤波法以及f-x域经验模态分解预测滤波法相比能够在衰减随机噪声的同时保留更多的有用信号。数值模拟实验以及实际地震资料处理结果均表明该方法能更为有效地去噪。

简介：角度域共成像点道集(ADCIGS)是偏移速度分析和振幅随角度变化分析(AVA)的基础数据。传统Kirchhoff叠前深度偏移(KPSDM)按偏移距组织数据,能方便的输出偏移距域共成像点道集(ODCIGS),其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为:(1)利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;(2)根据旅行时场的梯度方向计算反射张角;(3)在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是:(1)按照炮数据来组织输入数据;(2)旅行时场的输入与单炮覆盖范围相联系以节省内存;(3)多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：位场边界检测在位场数据处理中占有重要的地位,大多数的检测方法都是基于位场的梯度计算,因此算法易受干扰、稳定性较低。本文基于数理统计理论,不需要进行位场梯度计算,提出了利用各方向均方差相关系数进行位场边界检测,并且对算法及其合理性进行了详细的阐述与分析。在模型试验中,分别做了单一模型、组合模型以及加入随机噪声的组合模型试验,验证了方法的可靠性,并进一步与其它边界识别方法作了比较。各方向均方差相关系数法的特点为:算法简单稳定,结果辨识度较高,能同时对不同埋深的地质体边界都有较好的检测效果,能较好地保留边界形态,对噪声敏感度较低。最后将方法应用于老挝万象附近某地实测布格重力异常的处理中,利用研究区遥感解译的构造格架作为佐证,说明了各方向均方差相关系数法在实际应用中的可行性,为进一步判读区内构造展布提供了依据。

简介：地震波场正演模拟是地震资料处理、解释中最为重要的技术之一。地震波场正演模拟在大时间步长、长时程的波场延拓中,存在计算不稳定的问题。本文基于声波方程的Hamilton表述,在波动方程求解中用辛差分格式进行时间网格离散,用傅里叶有限差分进行空间网格离散,提出一种新的保结构地震波场正演模拟方法一辛格式傅里叶有限差分法,在保证计算精度的同时提高计算的稳定性。利用声学近似处理空间-波数混合域的积分算子,将该方法推广至各向异性介质。给出各向同性和各向异性条件下的地震正演模拟的计算流程,并将本文方法用于BP盐丘、BPTTI等模型的波场正演模拟。数值算例表明本文开发的方法适用于速度变化剧烈的复杂介质地震波场正演模拟,计算精度高,数值频散小,在各向异性介质正演中能够有效避免qSV波残余,在大时间步长的迭代计算中稳定性好。本文为在辛算法的框架下实现高精度地震正演模拟提供了一种新的选择。

简介：对VSP资料进行偏移成像可提高井附近地下结构的成像分辨率。本文给出了一种基于局域化相空间波场分解的VSP偏移成像方法。此方法采用了基于Gabor-Daubechies紧标架的延拓算子(G-D延拓算子)及其高频渐近形式对相空间波场进行延拓;基于局部平面假设,提出了一种局部角度域相关成像条件。合成和实际VSP资料的偏移成像结果表明,在满足渐近展开的条件下利用G-D延拓算子的高频近似式能够有效的减少计算时间;采用局部角度域相关成像条件能够在不增加计算量的同时,有效减弱VSP成像剖面上的偏移假象。更多还原

简介：常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI（ATI）条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的＂精确走时＂结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。

简介：边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原