基于球冠的卫星重磁数据处理和分析方法.pdfVIP

中国地球物理2011 ·337· (9)岩石圈结构及大陆动力学 基于球冠的卫星重磁数据处理与分析方法 SatelHte and Data GravityMagneticProcessing and MethodsBasedon Analysis SphericalCap 杜劲松+陈超+陈波王秋革 Du Chao。ChenBo Jinsong。Chen WangQiuge 中国地质大学(武汉)地球物理与空间信息学院武汉430074 随着卫星重力观测技术的发展，重力场的精度与分辨率都大幅度地提高，卫星重磁被越来越多地应用 于地球、月球与火星等内部物质分布、结构与构造的研究中。然而，区域性的重磁研究急需在球坐标系中 进行处理。 本文基本思想是在球冠内进行卫星重磁的数据处理与分析，研究不同点时，需要进行球坐标系旋转、 球坐标变换与模型重构，将新球坐标系中需重构点的球坐标反变换到原球坐标系中，利用已有卫星重磁场 模型直接解算。将球冠内的数据体建立后，可以进行如下几个方面的研究： (1)重力位、重力异常与重力梯度张量正演计算，用于重力地形校正、均衡校正、“剥层法”分离异常 及全球重力异常场合成计算。该方法的优点是采用球冠层重力正演的严格解析解公式，计算精度高，只取 决于地形或者界面模型精度，其缺点是由于需要进行坐标旋转，因此不适用于三维反演。 (2)GRACE时变重力场去噪，建立好球冠内的数据体后，可以设计各种空间域滤波器，本文首先设 计各向同性高斯滤波器，与球谐域内的高斯滤波比较验证了本文方法的正确性，进而设计各向异性滤波器 与自适应高斯滤波器。 (3)异常分离，在一定范围内区域异常的变化可以看成是线性的，而局部异常则可以视为随机扰动， 因此可以采用统计分析方法，即滑动平均方法(包括圆周平均与网格平均)与小子域方法，该方法分离异 常的效果取决于滑动平均的窗口大小，在实际应用中，最佳平滑窗口大小可以通过先验知识与试验方法估 计。在异常分离中出现的“虚假异常”，可以通过多次迭代消除。模型试验及实际应用均表明，该方法优于 直接利用截断重力场阶次分离异常的方法，能够反映实际地质情况。 (4)地质体边界信息提取，通过球谐与球冠谐计算球冠内各个网格点的各次导数与梯度等，进而采用 余钦范等(1994)的搜索水平梯度最大模的方法提取线性构造，结合位场延拓，可以构建地质体边界的拟三 维结构；基于球冠内的数据体，可以计算各个方向的各次导数。 (5)重磁对应分析，该方法是重磁异常联合定量解释的一种方法，它根据泊松定理，通过研究重磁异 常的相关性、斜率和截距，寻找异常与地质体赋存状态之间的关系，其难点在于如何确定窗口及球冠张角 大小：现在卫磁数据分析方法相对欠缺，进行卫星重磁对应分析研究区域磁异常是可行的。 (6)有效弹性厚度您计算，死代表岩石圈综合强度，通过对死观测值的比较，可以揭示一个地区 是否经历过热重建或其他显著构造事件，还可区分不同的地壳组成。往往其计算方法在谱域进行，本文采 用空间域基于球冠的相干法计算Te-球冠内数据体之间的间隔与球冠张角对应的球面弧长代表数据的最短 波长与最长波长，采用不同大小的球冠计算重力异常与地形的相关性。得到累计相关系数后，根据特征波 长计算岩石圈弹性厚度。 基于球冠的卫星重磁数据处理与分析方法为球面上区域性地质构造研究提供了一种方法，其最大优点 在于：一是基于球坐标，更加适合实际地质情况：二是不仅适于全球而且适于局部。 学基金联合资助。