基于离散余弦变换的磁位谱分析及磁异常导数计算方法

针对提高磁异常导数的计算精度,提出磁位离散余弦变换谱的分析方法. 根据重磁位场的泊松公式,利用余弦变换给出磁位与磁场分量间的余弦变换谱关系,推导出磁异常n阶导数的余弦变换谱公式. 利用余弦变换法计算的无限长水平圆柱体磁异常水平和垂向一阶导数的最大误差分别为-0.28 nT/m、0.47 nT/m;水平一阶导数的误差一般在-3.57%～3.27%之间,垂向一阶导数的误差一般在-1.94%～1.88%之间;计算的磁异常一阶导数值与理论值大致重合,而且不受有效磁化倾角的影响. 而Fourier变换法计算的水平和垂向一阶导数最大误差分别为-10.62 nT/m、14.42 nT/m,计算曲线与理论曲线偏离大,受磁化倾角的影响也较大. 这说明与Fourier变换法相比,余弦变换法计算的异常导数精度高,而且具有良好的稳定性.     

重力位谱分析及重力异常导数换算新方法——余弦变换

为了提高重力异常导数换算的精度,真实有效地反映地质体的异常特征,提出用余弦变换计算异常导数的新方法. 给出并证明了两个定理,利用它们推导出重力位余弦谱一般表达式以及重力异常各阶导数计算公式,建立了位场余弦谱分析理论. 模型实验中发现,用Fourier变换计算的一阶导数与理论导数偏差很大,而余弦变换计算的导数与理论异常导数拟合效果非常好,除边界几个数据因重力异常的有限截断产生的吉布斯效应残留使误差较大外,数据的计算精度均很高,误差为-0.09%～5%.     

基于离散余弦变换（DCT）的化磁极方法

针对提高磁异常化磁极的质量问题,提出基于离散余弦变换(DCT)的化磁极方法.以位场余弦变换谱分析为基础,从理论上推导了基于DCT的二度和三度体总场磁异常化磁极转换公式.在倾斜板状体模型实验中,化磁极误差小于0.001 nT,具有较高的精度;在单球体及多球体模型实验中,采用基于DCT的化磁极方法在5.倾斜磁化时就可以取得较好的化磁极结果,15°时化磁极的效果更加明显,其等值线的形态、幅值以及所反映的磁性体的水平位置都得到较好的恢复,这说明,采用本文方法进行化磁极时,可以取得较好的效果.     

利用余弦变换计算重力异常的向上延拓

利用余弦变换计算重力异常的向上延拓是一种新方法.根据余弦变换的基本性质,推导了二度、三度体异常向上延拓余弦变换谱理论公式,采用离散余弦变换实现了该法的数值计算;研究了无限长水平圆柱体的补偿因子中主频段的特性,给出了二度体的线性补偿方式;补偿后的理论模型异常向上延拓具有较高的计算精度,除边部几个数据因数据的离散和有限截断使误差较大外（最大误差为6.23%）,其余数据的误差均在1%以内,理论值和计算值曲线基本重合.这说明,与Fourier变换相比,离散余弦变换在数值计算中,受非周期性深度因子的影响小,补偿方式易于选择,其计算方法优于Fourier变换.     

由重力异常计算重力全梯度张量的余弦变换法(英文)

在傅里叶变换的基础上,提出在只有重力异常的情况下,利用余弦变换计算重力全梯度张量,并探讨了重力异常信号中含噪声时,余弦变换法计算的梯度张量的精度问题。通过对二维无限长均匀水平圆柱体及三维"Y"型岩脉模型梯度张量的计算表明,在原始重力异常中加入噪声与否的情况下,余弦变换法的计算结果均要优于传统的傅里叶变换法,特别是从均方差上看,余弦变换法的精度提高了近2'3倍。在实际应用中,实现了对黑龙江虎林盆地的重力梯度张量计算,计算结果表明,由余弦变换获得的重力梯度张量与傅里叶变换法的结果具有较高的一致性,但余弦变换法获得的梯度张量异常更加平滑,表明其受噪声的影响更小,能够更好地反映断裂的分布特征。     

磁异常模量垂向一阶导数的特征及应用

磁异常转换模量(magnitude magnetic transforms,简称MMTs)相比传统的数据处理方法如化极、总梯度模等,其产生的异常更接近实际磁性异常体的水平位置,以便进行磁法解释.本文基于转换模量,提出了磁异常模量的垂向一阶导数的处理方法.进一步分析模量垂向一阶导数的特征可知,对于2D异常体其完全不依赖于...     

航磁资料数字处理——垂直换算、垂直二阶导数及延拓计算

本文通过实例对航磁资料进行垂直换算、垂直二阶导数和延拓计算的数字处理方法作了介绍。本方法以偶函数(cos级数)为基础,比一般富氏级数方法提高了计算精度,并改进了边界效应。文章最后给出的计算程序是采用了快速富氏变换算法和递推算法相结合,不但计算速度快,而且对计算点数也没有严格限制。这一计算方法对于找矿、地质构造与深部地壳分析都是十分有用的     

基于离散余弦变换的集合Kalman滤波方法对河流相油藏渗透率场的反演

针对河流相油藏参数场的非高斯性及油藏生产数据与地质模型之间的非线性关系,提出了耦合离散余弦变换的半迭代集合Kalman滤波方法对河流相油藏进行生产历史拟合及反演.将河流相油藏非高斯渗透率场利用离散余弦变换变为类似高斯分布的DCT(Discrete cosine transform)系数,利用能量压缩特性对DCT系数进行截断,使得反演渗透率场的分布变得光滑,较好地解决了河流相油藏的非高斯分布问题.半迭代集合Kalman滤波法很好地解决了油藏生产数据与地质模型之间的非线性问题.理论证明了半迭代集合Kalman滤波与EnRML方法的等效性.比较了标准EnKF、耦合DCT的EnKF及耦合DCT的半迭代EnKF法对河流相油藏渗透率场的反演.讨论了DCT系数保留模式及系数保留数量对耦合DCT的半迭代集合Kalman滤波方法反演结果的影响,结果表明三角形系数保留模式可以用更少的系数保留参数场的地质特征,当保留2.5%的DCT系数时,反演结果已经充分反映了河流相油藏的地质特征.     

重力异常四次导数的计算及应用

利用泰勒级数展开式推导了二度空间水平四次导数及三度空间垂直四次导数计算公式。文中对四次导数若干理论问题进行了探讨,并对某些地区实际重力资料进行了计算,取得了较好的地质效果。     

重力异常四次导数的计算及应用

利用泰勒级数展开式推导了二度空间水平四次导数及三度空间垂直四次导数计算公式。文中对四次导数若干理论问题进行了探讨,并对某些地区实际重力资料进行了计算,取得了较好的地质效果。     

一种高效的基于逆Stokes公式的海域重力异常计算方法

利用卫星测高数据解算重力场是获取大范围海域重力异常的主要途径.然而,由于涉及到全球范围的多重积分,传统串行计算方法效率低下.针对上述问题,本文研究提出一种高效的海域重力异常计算方法.该快速方法利用逆Stokes公式进行重力场解算,同时针对其中不同点位多重循环计算的效率低下问题提出了一种基于多核CPU技术的并行化改进方法,该方法通过将程序合理地划分为CPU不同线程,在运算过程中合理调配计算机资源以完成海域重力异常的高效并行计算.文章利用EGM2008模型提供的大地水准面数据进行了快速算法测试,并将解算出的重力异常与船测重力异常进行了比对,表明该快速算法能够准确高效地求取出研究区重力异常,计算效率相比较于传统串行算法提升30%以上.本文提出的并行化方法对解算大范围海域重力异常和精化地球重力场模型具有积极意义.