Chebyshev低通滤波器在位场数据波数域导数换算中的应用

导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。     

位场高阶导数的波数域迭代法

位场导数换算是位场数据处理和解释的基本方法,但其不稳性限制了高阶导数在实际资料处理中的应用。针对于此,提出了位场高阶导数的波数域迭代法。在导数算子中引入一个包含求导阶次和高阶导数算子本身相关的低通滤波器,然后通过迭代法进行逐次逼近,进而得到高阶导数波数域迭代法的递推公式;同时证明了该方法的收敛性,对比分析了迭代法与常规方法导数算子的滤波特性。在高阶导数的模型检验和实例应用中,迭代法与常规方法的对比分析结果表明,迭代法具有稳定性强和计算精度高的优点。     

计算重力异常垂向二阶导数的DCT法——模型实验

基于离散余弦变换（DCT）的重力异常垂询二阶导数的计算方法是笔者提出的新方法。以无限长水平圆柱体为例研究了DCT法的计算规律：模型实验证实,该法计算的无限长水平圆柱体重力异常垂向二阶导数的精度与采样闻隔、剩余密度及圆柱体半径大小无关。与圆柱体轴心埋藏深度成正比。通过分析给出了基于DCT的合适的线性滤波方式,滤波后的重力异常垂向二阶导数与理论二阶导数拟合效果好,具有较高的计算精度。     

局部波数在磁异常解释应用中的方法技术

研究了提高频率域中导数计算结果的精度以及局部波数反演结果的稳定性和精度,对比了余弦扩边方法和最小曲率扩边方法在频率域中计算导数的精度。通过对比表明,最小曲率扩边结果的精度高于余弦扩边结果的精度,提出了向上延拓某一高度来进行局部波数的反演技术。此技术不但能够解决高频干扰,提高导数计算精度,而且利用不同延拓高度的反演结果能够得到一个平稳的反演结果,从而提高了局部波数反演的精度。理论模型试算验证了方法技术的精度,实际资料处理验证了方法技术的稳定性,达到了满意的效果。     

在波数域计算一维重磁异常导数的Matlab语言算法

利用Matlab内建的快速傅氏变换函数可以方便地在波数域计算重磁异常导数。介绍了基于Matlab语言的波数域求导的算法,给出了程序源代码,讨论了一些有助于提高计算精度的编程技巧。通过模型试验和数据分析,发现在计算垂向导数时波数域求导算法的精度比傅氏级数的精度有明显改善;而水平导数的计算,2种方法的精度相当。在某区钾盐勘探中,用该方法处理高精度重力剖面数据,取得了较好效果。     

应用计算机“推导”位场圆滑滤波及二次导数公式的方法

根据台劳级数展开式,将窗口内网格结点位场数据在计算点展开,且对应点两两相加,由最小二乘法求得窗口内网格结点的圆滑滤波算子和二次导数算子,将此算子与观测的位场数据褶积,就可得到位场圆滑后的数据及水平、垂直二次导数。     

基于常Q模型的分数阶粘弹介质数值模拟方法

基于常Q模型的解耦分数阶拉普拉斯算子粘滞波动方程,可以分开模拟振幅衰减和相位错动。但该方程拉普拉斯算子的阶数是随空间变化的,因此数值求解存在一定困难。这里基于截断的泰勒展开,经过一系列近似,推导出拉普拉斯算子的阶数与空间无关的解耦分数阶粘滞弹性波动方程。采用中心差分计算时间导数,使用交错网格伪谱法计算空间导数。数值算例表明,新的方程在处理非均匀介质时具有精度高,计算简便的优点。     

弹性波场和稳恒电流场的Frchet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的偏导数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度△Ｇ求得。把任意场源正演问题的解和Ｆｅｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解。文中以稳恒电流场问题为例，给出了直流电法视电阻率偏导数的数值结果。     

弹性波场和稳恒电流场的Frechet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的编数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度ΔＧ求得。把任意场源正演问题的解和Ｆｒｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解     

应用计算机“推导”位场水平导数公式的方法及程序

根据台劳级数展开式，将窗口内网格结点位场数据在计算点展开，且对应点两两相减，由最小二乘法求得窗口内网格结点的水平导数算子，将此算子与观测的位场数据褶积，就可求得位场的水平导数。