细胞自动机地震波模拟的并行化算法

采用区域分解的办法对格子气自动机地震波传播模拟的并行化算法进行了研究,提出了两种区域分解算法,并在SGIPowerChallenge并行计算机上进行了实际地质模型中平面波反射地震响应的计算,获得了较高的并行计算效率．     

几类叠前深度偏移方法的研究现状

叠前深度偏移是解决复杂地质构造成像的有力手段,关于叠前深度偏移方法的研究一直是勘探地球物理研究的热点和难点。本文在收集、整理国内外有关研究资料的基础上,有选择地对若干类型方法研究背景及现状、基本原理及方法特点等进行分析,目的在于作为研究工作的起点,把握叠前深度偏移方法的发展趋势,以期生成适合我国复杂油气储集特点的地震成像方法及技术手段。     

基于螺旋线上谱因式分解的地震波场隐式辛算法

均匀介质、复杂各向同性介质和各向异性介质中的地震波传播过程,可用统一形式的标量声波方程描述.考虑到在无损耗条件下,地震波方程描述了地震波场这一个无穷维的哈密顿体系随时间的演化过程,该过程为一个单参数连续辛变换,因而可以在其哈密顿形式表述下导出其辛格式.与显式辛算法相比,隐式辛格式对应的隐式辛几何算法具有无条件稳定的特点,可以允许较大的计算步长.但是由于隐式算法不可避免地面临高阶矩阵的求逆,其每一步的计算速度较慢.为实现矩阵快速求逆,文中采用了螺旋边界条件下谱因式分解的方法.在螺旋边界条件下,需要求逆的矩阵化为带状矩阵,而且其各列非零元素的位置和大小具有非常好的相似性,因而可以采用谱因式分解的方法实现快速LU分解.文中采用二阶精度的隐式蛙跳辛格式和谱因式分解方法,计算了常速度、层状介质和Marmousi模型中的波场.计算表明,隐式辛算法不失为波场计算的一种好方法.     

"九五油储"基金项目分析

本文概述由国家自然科学基金委员会与大庆油田联合资助的"九五油储"重大项目的立项背景、研究目标和学科内涵、实验基地及地质目标、组织措施、具体成果概述、实际效果、人才培养和"九五油储"后续的重大成果产出等相关内容.     

一种适于近场理论地震图宽频域计算的矩阵分解算法

理论地震图方法在地震震源过程的研究中得到了广泛的应用。为了研究震源过程的细节,必须利用近场地震资料的高频信息,但是在应用Haskell矩阵法计算近场理论地震图时,格林函数的高频成分的数值不稳定性是计算宽频域理论地震图的一个基本困难。因此,目前仅有不超过10Hz的算例。本研究采用Haskell矩阵的一种新的分解组合形式,在数值计算时有效地避免了在计算过程中数值结果的溢出,实现了近场理论地震图的宽频域计算。同时,由于利用了矩阵运算中的解析关系,减少了运算的次数,从而提高了计算的速度。根据本算法建立的Fortran程序,在Univoc-1100计算机上进行了数值检验。结果表明,计算得到的格林函数至少在0—40Hz频率域范围内仍保持良好的数值稳定性。本文给出的这种Haskell矩阵分解组合形式的矩阵元素简单,且具有一定的对称性。对于这种分解所包含的物理意义,将有待进一步深入研究。     

三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

地震静校正全局最优化问题的求解

针对地震静校正存在的非线性和多参数性等问题，本文提出了一种全局优化的剩余静校正方法.在模型更新规则上，把模拟退火方法和均匀优化设计方法结合在一起，使模型的更新更为合理，从而加快了搜寻全局最优解的速度.同时提出了模型的不同分量的温度参数和退火过程的选取方式，把温度参数的选取与地震剖面的能量联系起来，使温度参数的选取具有自适应的特点.该方法克服了常规模拟退火方法所具有的寻优空间不均匀以及退火参数需通过多次试验选取的缺陷.通过实际地震资料的计算证明，本文提出的方法对地震静校正问题合理而有效.     

共反射元弧叠加方法在火山岩成像中的应用

常规叠加方法假设叠加的地震道集的反射波来自地下同一个反射点,而地表观测到 的反射波场在空间上小于菲涅尔体范围是不可分辨的. 本文根据火山岩体特征,用元弧描述 火山岩的反射界面. 将来自菲涅尔体范围内的反射波沿着元弧界面求和,进而扩大叠加次数 ,达到增加反射信号能量提高地震反射波的叠加成像质量的目的. 文中给出了反射波走时参 数与地下元弧反射界面之间的联系和叠加实现的步骤. 计算结果表明,元弧叠加方法能够提 高剖面的信噪比,增加反射波同相轴的连续性.     

三维相干切片断层多边形检测

作者在本文中将三维切片视为二维图像，利用图像边缘检测技术，检测出相干切片上断层多边形，并将其直接投影到测网底图上，用于构造成图。该项技术不仅提高了用相干切片解释断层的精度，而且提高了断层解释效率，使相干数据体技术在三维解释中的作用得到更大的发挥。实际应用表明，三维相干切片断层多边形检测技术，是全三维高精度解释的有力工具，提高了解释效率，缩短了解释周期。     

透射波波场重建和介质层厚反演研究

本文以P波和SV波波场比值,提出层状介质情况下透射波波场的反向重建算法原理,并据此给出在已知介质速度值时,求介质层厚度的反演方法。文中还以深源远震记录波形资料给出反演算例,证实本文的算法是可行的。