基于CRUST2.0地壳模型高精度地形重力改正研究.docVIP

基于CRUST2.0地壳模型高精度地形重力改正研究 　　【摘 要】 对实测的重力数据进行重力归算时需要了解地表的密度结构，实算中常假设地壳密度为一常数。这一假设引起的误差将直接影响大地水准面的确定。为了提高地形改正的精度，本文提出利用CRUST2.0模型中的地壳密度分布计算地形改正，以得到大地水准面上更为准确的重力异常值。计算结果表明，由于忽略地壳密度分布引起的误差可达几个毫伽，地壳密度不均匀的影响在厘米级大地水准面的确定中应当于以考虑。 　　【关键词】 CRUST2.0；重力归算；地形改正；密度横向分布 　　【中图分类号】 【文献标识码】A 【文章编号】 　　1引言 　　 地球重力场是地球内部质量分布的直接映像,它是研究地球内部结构、地球动力学的重要基础信息[3]。根据Stokes理论计算大地水准面时需已知大地水准面上的重力值，且要求大地水准面以外没有质量分布，这就要求对测???的地表重力值进行重力归算。重力地形改正的计算需要了解地表到大地水准面间的密度结构，由于地质资料的不完备，实践中通常采用2.67g/cm3近似表示地壳密度。然而真实的密度分布与之差异可达10%甚至更多。对于高精度的大地水准面确定，这一近似引起的误差是不可忽略的[7]。本文分析了影响地形改正的各种误差源，对青藏高原地区的地形改正进行了计算，比较了不同计算模型的误差。应用CRUST2.0的地壳密度分布计算得到更符合实际的地形改正，结果表明密度横向不均匀的影响在该地区可达几个毫伽。 　　2 地形改正基本原理 　　 按Stokes理论确定大地水准面时不仅要求已知大地水准面上的重力值，还需假定大地水准面外部没有质量，所以需通过一定的重力归算方法消去大地水准面外部质量的影响。根据牛顿引力公式，大地水准面外质量对空间一点P所产生的引力位和径向方向上的引力分别为式（1），（2）： 　　(1) 　　 (2) 　　 其中，G为是地球引力常数，为流动点密度，为流动点高程，为计算点与流动点间的距离，为积分球冠。 　　 根据（1）、（2）式，当已知地表高程以及地壳密度时，即可精确的计算出大地水准面外部地形质量对一点的重力影响，从而求出其对大地水准面的影响。但由于地形与地壳密度信息的不完备以及计算模型的近似处理，在进行重力归算时存在着各种误差。 　　3 主要误差源 　　3.1 计算模型 　　 目前常用的方法包括严格积分法、近似法和组合法。这三种方法大都是以平面模型为基础。当假设局部范围内高程起算面为平面时，(2)式可表示为： 　　(3) 　　 由于不可能获得完全的地形数据，只能以一定精度，一定分辨率的DEM模型表示，假定每个DEM格网为一个高程不变、密度相同的棱柱体，如图1所示。则利用严格积分法计算地形改正的公式为： 　　 　　 　　 　　 图1 棱柱模型 　　 　　 严格积分法可以精确求定平面直角坐标系下棱柱体对一点的重力影响。但这种方法计算速度慢，对于距离较远的柱体计算结果不稳定[8]。在一定精度前提下，为了提高计算速度国内外学者研究了各种近似方法，通常采用的有点模型法、线模型法，以及McMillan法。本文算例中采用的是McMillan法，McMillan法是根据棱柱积分的球谐展开式所推导而得到的[8]。其具体公式如下。 　　 　　 （5） 　　 　　 （6） 　　 式中为棱柱体在三个方向的边长。 　　 近似模型法可以转换成卷积公式利用二维FFT方法计算，计算效率将大幅提高。但是为了保证级数展开式的收敛性，一般要求地形坡度小于45°，且近似误差所引起的影响在高山区可达到50mGal。李建成教授提出的数值积分法与快速Fourier变换技术的严密卷积联合既完全消除了数值计算误差又提高了计算效率。 　　3.2 密度横向不均匀的影响 　　 一直以来，由于难以获得精细岩石圈密度结构，在做重力归算时，通常设定地壳密度为常数2.67g/cm3。然而这一近似的误差，有时可达10%以上，对于大地水准面的影响可达分米级[3]。假设在高程为1km的区域，存在着两个相距4km、宽1km与周岩密度差分别为-0.1g/cm3和0.2 g/cm3的两块密度异常体。其对地形改正的结果如下图所示。 　　 　　 　　 　　图2 密度对地形改正的影响 　　 通常情况下，将地形改正分两步计算，局部地形改正和布格改正。而布格改正正是基于对布格片密度为一常数的假设推导而得的。 　　 　　 假设密度误差为，地形改正的误差，对于高程分别为1000米的地区,由于地壳密度误差引起的地形改正计算误差将分别达到11.19mgal，其对大地水准面的影响约11cm。可见地壳横向密度的影响在重力归算中是不可忽视的。 　　3.3 积分半径 　　 积分半径的选