青藏高原东南缘多尺度重力场及构造含义

本文基于EGM2008重力场模型研究了青藏高原东南缘均衡重力异常和多尺度的布格重力异常特征, 以鲁甸和景谷地震为例, 认识其深部构造环境和动力学过程, 为该区域的构造运动和地震孕育环境研究提供依据。 结果表明, 研究区布格重力异常和均衡重力异常与地质构造格局相关性较好, 川滇地块剧烈的区域布格重力异常和非均衡状态与其强烈的地壳变形、 断裂及地震活动密切相关。 强震多分布在断裂带两侧重力异常的过渡地带和高梯度带, 断裂带两侧横向和垂向的显著介质密度差异是强震孕育的深部构造背景。 布格重力异常和均衡重力异常揭示的鲁甸、 景谷震源区深浅差异性的重力异常特征, 暗示鲁甸和景谷地震孕震环境的不同。     

利用滑动窗口的多项式拟合算法进行重力位场区域-剩余异常分离

地表观测的重力位场是地形质量、浅部地质结构产生的剩余重力异常和深部地质构造产生的区域重力异常的叠加效应。基于尺寸可变的滑动窗口的二维低阶多项式拟合算法和格网距离（到中心点）倒数的定权规则在空间域对地面观测的重力位场数据进行了不同深度层的区域-剩余异常分离。这克服了常规算法仅在水平方向上区分不同异常空间分布及垂直方向上定性分离的缺陷。并利用构建的模型重力数据和实测重力位场数据分别进行解算,数值结果验证了该方法的有效性。     

航空重力数据向下延拓的逆泊松积分算法

给出计算航空重力数据延拓的泊松积分公式,通过向上延拓评定航空重力数据的精度。基于逆泊松积分,应用快速傅里叶变换算法和迭代法求解重力异常,进行向下延拓,通过对延拓结果和地面实测数据的比较分析,评价延拓方法的精度,并讨论在向下延拓算法中引入的低通滤波器的选择。实验表明,快速傅里叶变换算法的结果较好。     

利用Extrapolation Tikhonov正则化算法进行重力梯度三维密度反演（英文）

本文利用Extrapolation Tikhonov正则化算法处理重力梯度数据三维密度反演的线性不适定问题。与Tikhonov正则化方法相比,Extrapolation Tikhonov正则化方法减小了因正则化参数的引入而带来的反演结果误差,提高了预测数据与观测数据之间的拟合精度。同时为了消除位场数据反演时位置函数快速衰减对反演结果的影响,本文提出了基于重力梯度全张量特征向量法的深度加权函数,模型试验证明了该深度加权函数能有效识别异常体密度分布特征。对澳大利亚Kauring地区实测重力梯度数据进行反演,并和已有研究成果对比分析。结果表明该反演方法能够较好的获取地下异常体的密度分布信息。     

航空重力梯度的模拟退火参数反演

针对目前地球物理问题反演时,通常利用线性或拟线性反演方法,从而导致反演过程不稳定或反演结果陷于局部最优的问题,该文采用了模拟退火算法的非线性全局反演技术,从航空重力梯度全张量数据中反演了地下密度异常体的几何参数信息。同时,利用无噪声和有噪声的两组航空重力梯度数据进行对比实验。结果表明,在重力梯度数据信噪比小于21.584 3的情况下,通过引入噪声加权的价值函数,可以提高地下密度异常体几何参数反演的准确性。     

基于插值切割位场分离技术分析2016-1-21青海门源M6.4级地震前重力变化特征（英文）

本文基于覆盖2016年1月21日中国青海省门源县M6.4级地震的河西地区2011年至2015年每年两期的流动重力网观测资料,根据不同的拟合插值算法进行了区域重力位场值网格化方法优选,并通过引人应用于资源勘探领域中的插值切割位场分离算法(改进了边缘点和角点的处理)和图像处理领域的自适应噪声滤除技术,解算了门源M6.4级地震发生前,在震源深度层附近的区域重力场时空演化特征及研究区域重力观测点的半年尺度变化量,解算结果发现,在震源深度范围(-11.12 km)分离的重力位场自2015年下半年,重力变化达到-400~+150x l0-8m/s~2量级左右,大于2011年至2014年同期的变化结果(±30xl0-8m/s~2);半年尺度(2014年9月至2015年5月以及2015年5月至2015年9月)重力观测点值变化趋势也验证了门源地震前重力加速变化的特征,这可能在一定程度上显示了门源地震发震前地壳内部物质运移等地质活动在地表重力场上的反映。     

基于改进的2D多项式拟合及加权迭代算法的重力异常垂向分离

基于稳健的二维低阶多项式拟合及加权迭代算法，分别对模型重力位场数据和川滇地区实测地表重力异常数据进行地下不同深度层地质源的位场分离，并将所得结果与基于插值切割算法的位场分离结果以及川滇地区地质构造图进行比对。结果表明，该算法简单可靠，其结果可以用于川滇地区地质结构的精细反演和解释，也为其他更加深入的研究奠定了基础。     

基于实测重力异常和地形数据确定重力梯度的研究

本文利用澳大利亚北领地West Arnhem Land 地区实测重力异常数据并联合DEM(9")和SRTM3(3")地形高程数据,使用移去-恢复技术和Stokes积分方法计算了该地区两条剖面的重力梯度及其功率谱密度,使用FFT方法解算了整个地区的重力梯度值,结果证明了联合重力异常数据和高分辨率地形高程数据能有效地提高重力梯度的解算精度;功率谱密度的计算结果与国外成熟的重力梯度功率谱密度模型相吻合,表明高于0.3 Hz频率范围的功率谱密度可看做噪声,为重力梯度数据处理中噪声的辨别和剔除提供了借鉴,另外对重力梯度辅助导航基准图的构建以及重力梯度测量系统的标定提供了有益的探索.     

利用重力异常数据构建重力梯度场的方法比较

为实现大范围、高精度基准重力梯度数据库的构建,考虑到重力梯度场对地形质量的敏感效应,一般利用恒密度数字高程模型来求取重力梯度值,从而忽略了地形密度变化以及水准面以下密度异常对重力梯度的影响。根据重力位理论中求解边值问题的数值应用方法,直接利用重力异常数据求取重力梯度场,弥补了密度变化和密度异常在重力梯度上的反映。根据模型算例和实测重力异常数据求取了剖面重力梯度值,结果表明,限于重力数据空间分辨率的影响,利用重力异常数据可恢复中长波段重力梯度场。该方法与地形数据求取重力梯度和卫星重力梯度测量等方法技术相结合,对重力梯度数据库的建设具有实际应用价值。     

融合多源数据构建区域航空重力梯度扰动全张量

通过联合全球重力位模型(EGM2008)、航空重力扰动数据和剩余地形模型(RTM)数据,基于频谱域(二维FFT变换)和空间域(Stokes数值积分)算法对毛乌素测区GT-2A航空重力测量系统采集的空中测线后处理重力扰动数据进行解算,构建了该地区的航空重力梯度扰动全张量.(1)残余航空重力扰动延拓结果表明:残余航空重力扰动经向下延拓至大地水准面,再向上延拓至航空高度后与原数据差值的标准差为1.0078 mGal,考虑边缘效应后,内缩计算范围得到的差值标准差减小至0.1269 mGal.(2)基于残余重力扰动数据(原航空高度数据及向下延拓数据),通过不同方案解算得到的梯度扰动结果表明:两种方案得到的研究区域重力梯度扰动各分量之差的最大标准差为6.4798E(Γ_(yz)分量),最小标准差为2.6968E(Γ_(xy)分量),内缩计算范围后得到的差值标准差最大值为1.8307E(Γ_(zz)分量),最小值为0.7223E(Γ_(yz)分量).本文的思路和方法可为未来我国自主构建航空重力梯度标定场提供参考.