Kirchhoff积分时间偏移的两种走时计算及并行算法

走时计算的可行性在很大程度上决定了Kirchhoff积分法偏移的效果,而基于PC-Cluster的并行计算是实施叠前偏移数据计算量大的基本工具,本文针对这两个问题的研究,阐明两种旨在提高复杂情况叠前时间偏移聚焦能力的弯曲射线对称走时及非对称走时的计算方法,以及旨在提高叠前偏移计算效率的"数据并行"和"成像点并行"策略,并分别对其结果进行了对比分析.     

基于GPU的任意三维复杂形体重磁异常快速计算

提出了基于图形处理单元的任意三维复杂形体的重磁异常快速正演计算方法.将地下半空间剖分为大小相等规则排列的一组长方体单元,任意三维复杂形体可以表示成很多不同体积和密度(磁性)的长方体的近似组合.用解析方法计算出所有这些长方体在计算点的重力(磁力)异常,并累加求和,就可以得到整个模型体在计算点引起的重(磁)异常值.为了提高近似程度,需将地下半空间剖分得很细,用传统的CPU串行程序计算相当耗时.GPU在处理能力和存储器带宽上相对CPU有明显优势,采用GPU并行算法,可大大提高计算速度.相关试验结果表明,用GPU实现的正演快速算法计算结果正确,效率明显提高,为重磁异常三维物性反演提供了基础.     

一个基于多GPU的TTI介质逆时偏移方法（英文）

地震数据TTI介质逆时偏移计算量巨大,阻碍了其深度应用。NVIDIA? GPU以及其通用计算具有高度并行特性,为加速该类高密度计算提速提供了可能,可解决由于地震处理成像空间增加所带来的计算量大难题。本文提出了一种高效率的多GPU并行计算策略,以解决TTI介质逆时偏移庞大的数据处理问题。该策略流程,以GPU及其统一设备架构CUDA为前提,拥有多核心同时计算、作业流、点对点GPU直接传输等一些列运算特点。其核心是将GPU间边界数据的数据传输时间和不同区域的差分计算时间重叠。由于逆时偏移计算强度主要与差分计算空间相关。重叠计算后,数据传输时间可以忽略,因而计算效率随GPU数量的增加呈线性提高。并用于TTI逆时偏移成像的处理,以验证本文提出的高计算效率的多GPU并行计算策略的正确性。对比试验表明,利用本文提出的多GPU策略可大大提供高密度数据成像计算的效率及实现多GPU计算时效率呈线性增加,提升了该计算的延展特性。     

被动源反射地震数据起伏地表静校正问题分析

主动源地震中采用静校正技术将地震数据校正到平面上进行常规处理,被动源反射地震数据处理中同样面临该问题,其中涉及到了两种方法,一种是生成拟炮集记录后,采用主动源地震方法进行静校正处理;第二种方法是在背景噪声上,根据检波点高程将数据校正到平地表上再进行拟炮集记录计算.本文主要分析了第二种静校正计算方法的原理、以及实现步骤,并通过数值模拟起伏被动源地震数据,对两种静校正方法校方法进行了综合对比.实验结果表明：在被动源原始数据上静校正和在互相关生成的拟炮集上静校正都能实现起伏地表校正为平地表,但在被动源原始数据上静校正再相关会导致静校正量发生改变,需对反射记录非因果部分进行两倍炮点校正,因果部分进行两倍检波点校正才能叠加生成拟炮集.     

地球物理不规则分布数据的空间网格化法

不规则分布数据的网格化处理是地球物理数据处理和解释的基础问题,是保证许多地球物理数据处理方法得以成功实施的前提.文中介绍了线性插值法、多元二次函数法、普通克里格法、反插值法等4种空间网格化法的原理,并利用这些网格化方法进行了理论模型和实际航磁数据的试验分析.     

逆时偏移去噪方法研究进展

逆时偏移(RTM)作为精度较高的深度域成像方法之一越来越受到人们的重视..基于双程波动方程的逆时偏移具有原理简单、无倾角限制、可以用于复杂区域和陡倾角构造成像等优点.本文概述了国内外逆时偏移的发展历程和现状,介绍了应用随机边界条件后逆时偏移基本原理以及存在的问题并分析了成像噪声的产生机制,重点给出成像后压制逆时偏移低频噪声的几种方法及特征并给出总体认识.     

地震叠前逆时偏移算法的CPU/GPU实施对策

相较于单程波偏移算法而言,逆时偏移成像方法以其物理基础为依托优势,几十年来一直备受国内外地球物理学家的青睐.目前的逆时偏移(RTM)若直接采用双程波动方程进行延拓,尽管可以回避上下行波的分离处理,然就已有算法而言,其计算量和I/O(输入/输出)量却是最大的.针对此问题,本文在分析现行逆时偏移的多种算法基础上,提出利用CPU/GPU(中央处理器/图形处理器)作为数值计算核心,建立随机边界模型,从而克服存储I/O难题和提高计算效率.在实际的数据测试中,本文的方法可以大幅度的提高计算效率和减少存储单元,从而促使其高效地应用于生产实际.     

浅层地震在薄层采空区探测技术研究实践

济南东部城市采空区是该地区城市化进程的隐患.由于采空区厚度在1 m左右,常规方式难于直接分辨.地震勘探因其具有高分辨率的优势被采用,但是同样因为分辨率问题,需对地震成像数据进行提频处理后进行采空区特征的识别.本文采用基于同步挤压小波变换高频重建方法和正演计算相结合的方式,就该区薄采空区在地震剖面上的特征进行了识别.在高频重建后的地震剖面上,采空区因为低速填充的存在,出现煤层反射下拉的现象.通过钻井验证,该方法在浅层薄煤层采空区识别中效果良好.本文方法直接、稳定,综合应用拓展了浅层地震在城市地下空间探测的应用领域,对济南东部采空区识别具有重要实践意义.     

基于剩余异常相关成像的重磁物性反演方法

将场源区剖分成长方体单元,通过采集的重磁数据反演出这些单元的密度或者磁化率变化,勾画出场源的分布图像,这种方式是重磁三维反演的重要方向.重磁相关成像通过计算测量的重磁异常与地下各点在测区上的重磁异常的归一化相关,显示出异常地质体的空间赋存状态和等效剩余重磁物性.该方法计算速度快,方法简单、稳定,但是反演的结果只是在-1到+1之间的等效物性,不能够直接反演剩余密度或者磁化率,并且无法引入已知的地质约束.本文通过对物性模型的正演和实测结果的残差进行相关成像,迭代更新物性模型实现对物性参数的反演过程.模型实验证明该方法相对相关成像不仅能提高分辨率,还能够得到真正的物性参数.