三维混合延拓一步法波动方程正演模拟

提出三维相位移加有限差分混合延拓法,以此为基础,实现速度纵横向变化,构造任意复杂的三维介质一步法正演模拟,为三维复杂波场的地震,地质解释提供了精确,实用的方法。对于二维介质,相位移延拓方法对陡倾斜地层成像精度高,稳定性好,计算快,但实现的速度的横向变化困难,而频率、空间域的有限差分法算法简单,稳定性好,能适应速度的纵横向任意变化,但偏移陡倾地层存在很多问题,用快速的４５°有限差分对相位移延拓作补充     

基于频率域高阶有限差分法的正演模拟及并行算法

在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。     

波动方程有限差分正演模拟误差来源分析

有限差分波动方程正演数值模拟主要是解决差分格式、震源选择、边界条件处理等问题.这些问题的合理解决对得到正确和准确的波动方程数值解有着十分重要的作用.笔者通过对波动方程计算机有限差分正演模拟过程的叙述,分析了在做正演模拟时的误差来源以及如何有效地减小误差.     

波动方程显式短算子三维正演模拟

通过对三维波动方程的显式差分法的研究,并对模小于1.0的稳定性约束作线性化处理时引入加权因子,提出了一个基于二次规划方法的一维显式短算子设计方法.与以往的显式算子设计方法的不同之处,在于对传播区内算子的模也引入了最大值限制,这样可更严格地保证深层波场延拓的稳定性.同时,给出了基于一维显式算子及McClellan变换的三维波场深度延拓方法,大大减少了计算量,提高了计算速度,并能更好地处理介质横向非均匀性.     

可控震源粘弹性波动方程有限差分模拟

在一般的地震波计算中,地球介质可以当作各向同性的完全弹性体来对待。而对地表土壤和精细观测,粘弹性介质模型比完全弹性模型更符合实际。可控震源是近年来发展起来的精确可控、安全有效的新型震源,但国内对可控震源传播动力学的研究较少。对可控震源激发的地震波在粘弹性介质中的传播,用有限差分方法进行了数值模拟,给出了Kelvin一阶粘弹性波动方程组、吸收边界条件和震源处理。提出了检波器接收信号加速度和传播距离的三次多项式指数关系式,比线性指数关系更准确。使用该公式可以预测采集系统有效接收距离,指导测线布置。     

基于球坐标系下有限差分的地磁测深三维正演

为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。     

地震波的有限差分法正演模拟

基于波动方程的地震波正演模拟能够准确揭示波在介质中传播的振幅、频率与相位变化,真实地反映出波的动力学特征。有限差分法因其精度高、计算效率高而成为波场数值模拟的一个重要方法。根据地震波的近似方程———声波方程,采用二阶有限差分网格数值算法,利用透明边界条件,设计典型地质模型,进行地面波场的地震记录模拟。通过稳定性条件分析、震源子波优选以及频散抑制,成功合成了倾斜地层模型与正断层模型的地面共炮点波场记录,直观地显示出直达波、反射波、绕射波等丰富的波场信息,获取各模型的波场快照。结果表明,有限差分法结合透明边界条件对于获取介质丰富的波场信息具有很大的优势。     

SeisModTM波动方程正演模拟软件在圈闭识别中的应用

准噶尔盆地石西凸起侏罗系八道湾组具很大勘探潜力,现有地震资料对低幅度构造背斜圈闭不能进行有效识别.为满足圈闭识别需要,在不同解释方案基础上,利用SeisModTM波动方程正演模拟软件构建3种地质模型,正演模拟出不同地层构造模式下地震响应特征.通过对正演模拟结果对比分析,为背斜圈闭解释方案提供理论依据,对该区侏罗系八道湾...     

基于波动方程正演模拟分析薄砂层含不同流体的AVO特征

塔河油田碎屑岩储层厚度较小,薄层AVO特征受薄层内多次波反射以及薄层调谐效应和频率滤波效应影响,导致薄层AVO特征与厚层AVO特征差异较大。为此,笔者设计了不同含油气性薄层砂体模型来进行波动方程正演模拟,得到了CRP道集记录,提取了目的层反射振幅,基于入射波振幅计算了反射系数,从而获得了薄层地震正演AVO曲线,分析了塔河油田碎屑岩储层含不同流体时的AVO特征,找出了实际地质条件下的含油气性与AVO属性特征的对应关系,以指导地震剖面上的含流体预测工作。     

转换横波共炮记录单程声波方程正演模拟

本文分别把分离的纵波和横波方程分裂为各自的两个单向旁轴方程（Ｊ．ＦＣｌｅａｒ－ｂｏｕｔ）,然后用单程纵波和横波方程,结合相位移加有限差分混合波场延拓法（Ｄ．Ｒｉｓｔｏｗ）,在共炮记录叠前深度偏移原理的基础上,借鉴几何地震学的射线法两点问题原理,提出合成转换Ｐ—ＳＶ波共炮记录的时间一致性原理和逆过程方法,从而推导出精确的转换Ｐ—ＳＶ波共炮记录正演模拟算法,由于单程延拓波场,可完全解决多次反射波问题。     

三维定位原理与三维反射波场模拟

在三维正演模拟定位原理基础上 ,以单程声波方程作为波场外推手段 ,实现三维叠前反射波场模拟 ,完全克服了全波方程模拟的 3个问题。在波场外推中 ,使用波动方程的混合解法 ,即以波动方程解析解 (相位移法 )为主 ,在介质速度横向变化的地方 ,用数值解 (45°有限差分 )作补充 ,使得整个算法精确、快速 ,既能处理复杂构造 ,又能适应速度的任意变化。     

直流电阻率三维正演的预条件共轭梯度法研究

三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。     

三维任意形体高密度电阻率法数值模拟方法

为研究高密度电阻率法对各种三维地电体的反映规律，本文介绍了正演程序编制的方法并给出了计算实例。     

WINSURF二次开发在瑞雷波频散曲线正演模拟中的应用

WINSURF二次开发是物化探数据处理的一种强有力工具,以其在瑞雷波频散曲线正演模拟中的应用为例,阐述了WINSURF二次开发的基本原理、方法和技巧,并给出了源代码。正演模拟试算结果表明,该程序稳定性好、通用性强,可直接作为工具来使用,从而大大缩短了程序开发周期,避免了重复缟程和资源浪费,提高了科研和工作效率。     

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

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

成矿过程三维并行模拟技术探讨

成矿作用的多阶段多过程，成矿因素的多变量，成矿空间的多维多分辨率，导致成矿过程模拟的复杂性和在单机上实施模拟成矿过程的困难性，利用由微机组成的以太网构成的机群并行环境和内在共享,任务分配的并行处理策略，可以实施模拟三维空间的成矿过程，采用Client-Server结构，插件和除算法技术来构建的模拟系统原型，在由8台微机构成的并行环境中具有线性加速比。     

三维井地电阻率有限元数值模拟及反演

油气资源和矿产资源勘探中,常需要较大的勘探深度和高的深部分辨能力,充分利用现有钻孔的三维井地电阻率成像技术能够适应于这种要求。通过设计一个高阻水平板状体模型,利用Ansys有限元软件对模型进行三维井地电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律,其结论和认识具有一定实用价值。     

台吉营子金矿三维建模及储量估算

三维矿山地质模型能够完整准确地表述各种地质现象,直观展现地质空间分布及相互关系并为矿山动态管理和合理利用资源提供依据。利用3DMine矿业工程软件,整理并数字化台吉营子金矿大量地质勘探资料,探索并建立了矿区地表模型、1号矿体模型、探矿工程模型,实现了矿区三维数据的可视化。在构建矿体模型基础上,对1号主矿体进行探索性储量估算。     

二阶弹性波动方程高精度交错网格波场分离数值模拟

给出了一种等价的二阶弹性波动方程,以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波波场问题.应用高阶交错网格有限差分法求解该波动方程,并使用通量校正技术(FCT)进一步压制频散,采用均匀介质模型和层状介质模型进行波场分离数值试验,精确得到了混合波场、完全分离的纯纵波及纯横波波场.数值结果分析表明,本文方法在均匀介质情况下准确可靠,在分离后的纯纵、横波波场中可观察到较为丰富的能量转换信息,这对认识复杂弹性波的传播规律及弹性波理论具有重要意义.