论文部分内容阅读
地球重力场是地球的固有物理特征,可以反映地球内部的物质分布、运动和变化的规律。重力测量是对重力场变化规律的直接反映。随着测量手段的多样化(地面测量、海洋测量、航空测量和卫星测量)和测量数据精度的提高,人们不仅可以直接测量地球重力场,而且可以对重力梯度场进行测量。重力梯度张量数据相对于传统的重力测量数据有着更高的频率信息,能更加准确、细致的研究地球浅部构造、矿产资源分布等信息。为了更加准确的处理和解释全张量重力梯度数据,需要从全张量重力梯度仪器设计原理出发,分析仪器产生噪声的原因,进而更好的从测量数据中提取反应地下异常体的真实信号进行解释。首先,本文利用理论模型对比分析了重力梯度测量相对于重力测量的优势,为文章的选题提供依据。重力梯度张量数据的获取,可以通过实际测量或计算得到。然而,计算得到的重力梯度张量数据最多能包含与测量的重力数据相同的信息成分,而且在计算过程中还会造成信息缺失,并不会增加信息成分。因此,重力梯度张量数据的获取通常需要通过重力梯度仪器测量得到。通过对比分析,重力梯度异常能反映相对波长短的高频信号,重力异常数据能反映相对波长长的低频信号。因此,联合梯度张量数据的高频信号成分和重力异常的低频信号成分可以得到加强的重力异常数据,该数据能同时保留重力异常数据及其梯度张量数据的所有信息,拓宽重力异常数据的频宽范围。这一过程主要通过余弦平方滤波完成。目前,国内外多家研究机构正在研究多种不同类型的重力梯度仪,但是,投入生产使用的只有基于旋转加速度计的重力梯度仪。因此,本文针对用于航空移动平台高精度全张量重力梯度测量系统,深入研究了全张量梯度仪12个加速度计按3个不同旋转轴圆盘形成差分组合的结构。在确认该结构具有有效抑制运动共模加速度减少外界环境干扰实现高精度探测优点的同时,着重分析了重力梯度仪测量误差的来源和影响。研究表明,主要影响包括仪器固有随机噪声和外界确定性噪声。为定量描述影响程度,笔者推导了在航空动态环境下的测量方程,并分离出加速度计的性能不匹配、平台不稳定、圆盘转速不稳定3个主要固有因素,从时间域和频率域角度定量分析固有影响因素的噪声水平。试验分析表明,通过利用Simulink仿真系统可以获得固有因素产生的噪声水平,并提出抑制方案。针对搭载环境测量误差,笔者还分析了实测飞行中姿态和质量的改变对重力梯度测量值造成的环境影响,提出了基于点质量源的自身梯度校正方法。经固有噪声和自身梯度校正后的各个梯度数据中还存在大量随机噪声,但是其噪声水平各不相同。本文利用奇偶测线网格化方法对各个梯度张量数据的噪声进行了定量估计,为后续噪声去除提供依据。全张量重力梯度张量数据相对于重力数据含有更高频的信号成分,能更好的描绘小的异常特征。然而,全张量重力梯度仪测量值的噪声成分也为高频。因此,从高频信号成分中分离出噪声将是处理重力梯度张量数据的一个挑战。通过联合全张量重力梯度仪测量的多个信息成分进行噪声处理能更好的压制随机干扰。本文利用重力梯度信号满足的拉普拉斯方程约束条件,推导了重力梯度张量的笛卡尔方程和引力位的级数解表达式,然后利用级数解去拟合测量的重力梯度张量值,拟合的部分认为是真实的梯度信号,未拟合的部分认为是噪声。在拟合过程中引入数据噪声加权矩阵和数据能量加权矩阵进行最优线性反演求解拉普拉斯方程级数解的系数,然后利用求得的系数进行正演计算得到真实的梯度信号。全张量重力梯度数据经噪声滤波后,只包含实际地下重力信息,能更加准确的进行数据的解释。重力梯度数据的解释工作通常需要从数据中获取有关场源的水平位置和深度范围。水平位置的确定通常利用边界识别方法;深度范围参数的确定需要深度计算方法来完成。边界识别方法在重力梯度数据解释中占有重要角色,它能准确且快速地确定地质体边界位置、构造水平位置而被广泛关注。一些传统的方法有的不能同时显示不同埋藏深度地质体的边界信息,有的在深部地质体边界位置确定中误差较大,且对细节信息的提取能力不足。而且,已有的方法基本上都是针对重力异常数据,专门针对重力梯度张量数据的方法还比较少。本文针对全张量重力梯度数据信息量大、信号频率高,能更好的描述小的异常特征等特点,提出了改进的水平解析信号方法、加强水平方向总水平导数方法和改进的结构张量算法进行全张量重力梯度数据的边界解释。在解释过程中,本文针对各个方法进行相应的归一化处理,使其能均衡不同埋藏深度的异常体的边界结果的振幅强度。通过模型试验和实际测量的全张量重力梯度数据验证了这些方法的可靠性、实用性和抗噪能力,并与一些传统的边界识别方法进行对比,证明了这些方法的优点。重力异常及重力梯度异常场源深度计算的方法经历了漫长的发展历史,形成了针对不同的场源类型和地质研究对象的不同深度计算方法。快速成像方法是近几年的一个发展热点,它能快速得到地下异常体分布状况,避免传统方法耗时长,内存消耗大等缺点。本文利用极大值深度计算方法DEXP进行地下异常体的埋藏深度成像,从而进行异常体的深度估计。然而,传统的DEXP方法需要事先指定异常体的构造指数,通常是根据异常形态对其进行假设,但这会对成像结果带来误差。因此,本文又利用不同阶的垂向导数的比值进行DEXP变换,有效的去除了构造指数的影响。该方法能在不知道构造指数的情况下对异常体进行深度成像,从而得到地质体的深度。利用估计的深度值在尺度函数上的对应值可以对地质体的构造指数进行估计。利用该方法对模型数据和实际测量的文顿岩丘的重力梯度数据进行分析,取得了准确的深度和构造指数估计。