重力异常多尺度反演研究及应用.pdfVIP

·338· 中国地球物理2011 (9)岩石圈结构及大陆动力学 重力异常多尺度反演研究及应用 onMulti-scaleInversionof Anomalies Study Gravity 侯遵泽1杨文采2 HouZunzel Wencai2 Yang 1．中国人民武装警察部队学院 河北廊坊065000；2．中国地质科学院地质研究所北京100037 利用重力资料进行地质解释的一项重要任务就是通过对重力异常的反演，求取地壳密度结构或密度分 布。近几年，我们研究重力场的小波变换分析方法、并结合功率谱分析和广义线性密度反演等方法，创新 了重力异常多尺度密度反演方法，用以描述地壳的密度结构。在中国大陆重力场的多尺度密度反演试验… 和塔里木盆地密度结构反演【2J的应用表明，这种集成所形成的方法技术是成功的。 采用小波多尺度分解方法对重力异常作多重分解，我们应用Mallat塔式算法，构造出专门的小波基。 无论从原理分析还是实际计算都可得至JJd,波多尺度分析的一个准则，我们称为低阶细节不变准则，它反映 了离散小波变换的优越性，此准则表述如下：设位场数据矩阵为D，经疗阶(以为整数，≥2)离散小波变 换后产生小波细节Dl’．··A和刀阶逼近A．Ⅳ，其中小波细节D1．．··A—l不随11的增大而改变。这个准则说明， 不管怎么选择甩，小波变换出来的低阶小波细节都是一样的，所不同的只是小波细节的个数和n阶逼近。 这一准则是离散小波变换特有的优点，对重力异常分解非常有利。例如，大陆地壳在许多地区可分为上、 中、下地壳三层，每一层的厚度横向变化很大。但是，大陆地壳在一些地区只分为上、下地壳两层，而且 莫霍面有明显起伏。对大区域布格重力异常的离散小波变换而言，不同人选定的不同阶数刀会造成多尺度 分解结果的不同，这种似是而非的结果是不利于科学解释的。低阶小波细节不变准则的启示就是：由于低 阶小波细节JD．…不会随总的分阶数改变，可以放心地先对低阶小波细节的分析结果和对研究目标实际要求 最后确定总分阶数n。因此，采用这种战略不需要事先人为选定阶数，l，而同样可以取得客观的重力场小 波变换多尺度分解。 在完成重力场小波变换多尺度分解之后，将重力异常各阶小波细节和逼近作二维Fourie变换，求它们 的平均功率谱，即各个方向功率谱曲线的平均，看对应的对数功率谱的直线段，其斜率是否随小波细节的 阶数增加而增大。根据位场波数域的理论，可以根据功率谱直线段斜率或它的渐近线在横轴上的切点计算 出板状场源的埋深。无论在解析上还是在统计意义上，重力场功率谱斜率都正比于场源埋藏深度。因此， 由功率谱斜率可以计算出小波细节及逼近对应的场源深度，使小波分析结果有对应的平均埋藏深度概念。 最后，用广义线性反演方法求取目的层的密度分布。包括引力场在内，位场正演问题是求下述Poisson 方程边值问题的解： 卜Ⅷ抄枷。2擘y，?，、--00X,乃“∞ (1) 【 Ul，o=V(工，y) 其中v(x。v,zO)为地表位场值，有 啦∽一击垣珠考詈转出出砂· (2) 反演问题就是给出v仅∥，由(2)求等效层体密度p6x,y,z)。 参考文献 【1】侯遵泽，杨文采，刘家琦．中国大陆地壳密度差异多尺度反演Ⅲ．地球物理学报，1998，4l(5)：651--656． 【2】侯遵泽．杨文采．塔里木盆地多尺度重力场反演与密度结构【J】中国科学：地球科学，2011．41(1)：29～39．