适于Kirchhoff叠前深度偏移的地震走时李代数积分算法

本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率.     

沿圆曲线的Radon变换数值解

研究具有紧支集且在支集内连续的二元函数沿上半圆曲线的Radon变换反演问题。基于对投影函数的Fourier变换,反演问题可以归结为具有弱奇性及震荡核的Abel积分方程的求解。我们证明了当圆曲线中心及半径在一定范围内变化时,在已知沿上半圆曲线的Radon变换情况下,这个积分方程的解具有唯一性,并给出了消除Abel积分方程弱奇性的数值方法。在考虑投影数据噪声的情况下,给出了多次加权改善系数矩阵条件数稳定的数值方法,并通过数值模拟验证所提出方法的有效性。     

用积分方程法模拟各向异性地层中三维电性异常体的电磁响应

本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律.     

等效偏移距方法及应用

当勘探目标地质构造复杂、地层倾角较陡或岩性横向突变时,在地表接收到的反射信号较弱、信噪比较低,甚至接收不到反射波,常规偏移成像方法对这类数据很难精确成像,因此可以考虑基于散射理论的等效偏移距方法来成像.该方法分两步进行：第一步是抽取共散射点道集;第二步是对共散射点道集进行克希霍夫积分偏移.文中数值模拟及实际资料的实例表明,该方法在复杂地质目标及低信噪比地震资料的成像处理中具有一定的实用价值.     

实际井眼条件下过套管电阻率测井响应的传输线方程正演算法

过套管电阻率测井通过测量套管壁电势实现测量地层视电阻率,不同的地层模型对过套管电阻率测井会有不同的测井响应. 在传输线方程系数中本文考虑了径向含多个界面地层的影响,对轴向（井轴方向）呈层状、径向阶跃变化地层模型,给出了电势分布递推公式,计算了过套管电阻率测井响应,绘出了测井响应曲线. 算例表明,该方法有较强的轴向地层界面分辨能力;对径向均匀地层计算结果非常接近地层模型电导率,能较好地反映实际地层特征;对高电导地层,低电导水泥环有明显的测井响应;当径向存在多个地层界面时计算结果则为径向各层电导率的综合反映,不是某一径向层的电导率;算例显示该方法有较快（算例运行时间在1s内）的计算速度.     

O(T2)精度下椭球界面Dirichlet边值问题的积分解

研究了边界是参考椭球面的Laplace方程Dirichlet边值问题的求解,在O(ε4·T)精 度下给出了参考椭球界面上扰动重力位Dirichlet外问题的积分解式. 该结果理论上优于目 前常用的球近似下的积分解式,从而为研究物理大地测量中边值问题的求解提供了新的依据     

基于积分图像的面状注记配置方法研究

提出了一种基于积分图像的面状图形注记自动配置算法,设计并实现了基于积分图像的面状注记配置的核心算法。实验表明,该算法能快速准确地定位待标记注记的位置,合理而又美观地实现面状图形注记的自动配置。该算法已成功应用于深圳市地籍测绘大队的房产测绘系统中。     

理论地震图计算方法

近二十多年来,随着理论地震学和计算技术的发展,理论地震图的计算技术有了飞速发展.目前已能对不同的介质模型和震源计算各种体波、面波、地球自由震荡和静态位移场.它们在研究天然地震震源过程、地球内部结构、近场强地面运动、核爆作监测以及地震勘探等领域中发挥了越来越大的作用. 目前常用的理论地震图计算方法主要包括积分变换法、离散数值方法和射线方法几大类.本文对这些方法分类进行了简单介绍,并评论了各种方法的特点和各自的适用范围.     

海洋重力大地水准面的确定

美国海洋卫星测高仪的出现,使应用Hotine积分确定海洋大地水准面成为现实。本文通过对Hotine积分及垂线偏差的计算公式进行改进,较好地改善了求和项的收敛性,减小了截断误差影响,并提出了利用Hotine函数和重力异常确定海洋大地水准面的方法。实际计算表明:海洋重力大地水准面的精度在1米以内;卫星测高大地水准面间存在0.5米系统差;它和海底地形有一定的相关性,能较好地反映出海底地形的宏观特性。     

极低频电磁法超高速振荡Bessel积分数值计算的解决方案——积分外推法

极低频电磁法需要考虑电离层和位移电流的影响,Bessel积分核函数超高速振荡,超过了现有滤波系数的截止频率,收敛慢,甚至不收敛,导致数值计算极为困难,即使采用现有的、最高采样密度的Hankel滤波系数进行数值计算,高频时的计算误差大,不能满足需求.针对以上问题,本文将积分外推算法应用到极低频电磁法的数值计算中.首先通过计算CSAMT在均匀半空间的数值解与解析解验证算法的准确性;然后研究了极低频电磁法三层模型中给定收发距的波场频率响应曲线及给定频率的波场结构特征,积分外推法的低频结果与高采样密度Hankel滤波计算结果一致,高频结果优于Hankel滤波计算结果.本文研究表明,积分外推算法可有效解决超高速振荡Bessel积分数值计算难题.