茅山地区地壳深部构造与溧阳地震关系新见——重力延拓和莫霍面计算

本文据茅山地区1∶20万大比例尺布格重力异常及其延拓结果,并结合爆破地震和天然地震等资料,分析了茅山断裂带的延伸规模及切割深度.认为茅山断裂带属地壳断裂;茅山断裂带与近东西向的石桥—无锡隐伏地壳断裂相交汇构成了溧阳地震的深部构造条件.

重力场的局部虚拟向下延拓以及利用EGM96模型的模拟实验检验

阐述了基于虚拟压缩恢复法导出的虚拟向下延拓法,该法从理论上解决了重力场的向下延拓问题。为了检验理论的正确性,选取EGM96模型进行了模拟实验检验,实验结果支持局部虚拟向下延拓法。同时,对观测值加进随机误差的情形也进行了模拟实验检验,结果表明,局部虚拟向下延拓几乎不损失精度。

重力异常向上延拓严密改化模型及向下延拓应用

重力异常向上延拓全球积分模型在航空重力测量数据质量评估和向下延拓迭代计算等领域具有广泛的应用。为了消除积分核函数奇异性影响,需要对该模型进行基于积分恒等式的移去-恢复转换及全球积分域的分区改化处理。在此过程中,传统改化处理方法往往忽略了全球积分过渡到局域积分引起的积分恒等式偏差影响,从而导致不必要的计算模型误差,最终影响向上延拓计算结果的可靠性,甚至影响向下延拓迭代解算结果的稳定性。针对此问题,本文开展了重力异常向上延拓积分模型改化及向下延拓应用分析研究,依据实测数据保障条件和积分恒等式适用条件要求,导出了重力异常向上延拓积分模型的分步改化公式,提出了补偿传统改化模型缺陷的修正公式,并将最终的严密改化模型应用于重力异常向下延拓迭代解算。使用超高阶地球位模型EGM2008作为标准位场开展数值计算检验,分别对重力异常向上延拓分步改化模型的计算精度及在向下延拓迭代解算中的应用效果进行了检核评估,验证了采用严密改化模型的必要性和有效性。

对比评价两种重力位场延拓方法

本文针对两种快捷实用的重力位场延拓方法——数值积分法位场转换、直接解拉式方程法位场转换法,阐述了其方法原理。通过建立重力理论模型、数值模拟分析,进行了研究对比。

GT-2A航空重力仪精度评定

选择地形平坦的实验区A和地形起伏较大的实验区B,采用GT-2A航空重力测量数据和高精度地面重力测量数据,依据测线网交叉点平差方法评定GT-2A航空重力仪的内符合精度,通过地面重力与航空重力延拓结果的比对实验评定其外符合精度。结果表明,GT-2A重力扰动测量结果的内符合精度优于1mGal,外符合精度在实验区A优于2mGal,在实验区B优于3mGal。

海空重力测量及应用技术研究进展与展望(四):数值模型构建与综合应用技术

分析总结了海空重力测量数据向上和向下延拓技术的研究现状及发展方向,简要论述了海空重力测量多源数据融合处理技术的研究动态及发展前景,分析讨论了海空重力测量数据应用于地球外部重力场赋值和大地水准面精化技术的研究进展及发展思路,为海空重力测量数值建模与数据综合应用技术的未来发展提供参考。

虚拟压缩恢复法在向下延拓问题中的应用

阐述了基于虚拟压缩恢复法的重力场向下延拓,构造了点质量处于内部球内部、内部球外部和在内部球内外各有一个点质量三种点质量异常模型,试验结果支持基于虚拟压缩恢复法的向下延拓。

利用高阶径向导数带限模型进行重力向下延拓计算

泰勒级数展开是实施位场向下延拓解算的主要方法之一。该方法的有效性主要取决于位场延拓参量各阶垂向(或径向)偏导数的求取精度及其可靠性。为了避免使用封闭解析核函数在球边界面出现奇异性带来的不确定性问题,依据各类重力观测经滤波处理后均表现为一类有限频谱带宽信号的特点,提出将地球外部重力异常泊松积分式的核函数表示为球谐级数展开式,并将其截断为与重力观测值频谱范围相一致的带限求和式,进而通过直接求导方法,推导得到一组与带限核函数相对应的重力异常高阶径向导数带限计算公式,同时对该组公式进行了实用性改化,并将其应用于重力异常向下延拓泰勒级数展开式计算。使用超高阶地球位模型EGM2008设计两个阶段的数值计算检验方案,分别对重力异常高阶径向偏导数带限模型及其泰勒级数展开向下延拓模型的计算精度进行了检核评估,表明新模型具有良好的可靠性和有效性,在解算稳定性和计算精度两个方面都优于其他同类模型。

差分技术在海洋测量中的应用

系统分析了差分概念在无线电定位、卫星定位、平均海面和深度基准面传递、水位改正、重力场和磁场向下延拓等方面的应用模式。差分技术在卫星导航定位中的应用非常普遍,在海洋测量中举一反三,还可以推广到更多的应用领域,有效消除系统误差,解决很多的实际问题,提高信息获取与处理精度或者降低问题的复杂性。

对《基于离散余弦变换的磁位谱分析及磁异常导数计算方法》一文的评论

被评论的论文(简称原文)中应用离散余弦变换(DCT法)计算磁异常导数,并与离散Fourier变换法(DFT法)计算结果进行对比,得到"DCT法明显优于DFT法"的结论.我们对原文进行了验证.原文磁异常导数波数响应计算的错误导致DFT法计算的错误,最终导致结论的错误.原文相关两篇论文中应用DCT法分别进行重力异常的导数换算和重力异常的向上延拓,也存在同样的错误.

应用任意采样点数FFT算法时离散频率计算

快速傅里叶变换(FFT)算法在位场、电磁场和波场等地球物理场的高效数值模拟及其数据处理中发挥着重要作用.FFT算法本质上是一种实现离散傅里叶变换的计算方法,目前多种类型FFT算法组合,能够实现任意采样点数的离散傅里叶变换.离散频率计算是应用任意采样点数FFT算法求解地球物理场数值模拟和数据处理等问题的关键环节.本文从离散傅里叶变换作为傅里叶变换的一种数值逼近的观点出发,通过推导和分析任意采样点数离散傅里叶变换数学表达式,给出了离散频率的计算公式.以重力场向上延拓问题为例,通过理论模型数据实验,检验了本文给出的离散频率计算公式的正确性.