基于重力异常导数反演垂直断层模型参数.pdfVIP

第 23 卷第 6 期 中 国 科 学 院 研 究 生 院 学 报 Vol. 23 No. 6 2006 年 11 月 Journal of the Graduate School of the Chinese Academy of Sciences November 　2006 ( ) 文章编号 2006 简 报 基于重力异常导数反演垂直断层模型参数 陈　石　张　健 ( 中国科学院研究生院计算地球动力学实验室, 北京 100049) ( ) 2006 年 1 月 11 日收稿; 2006 年 4 月 10 日收修改稿 Chen S , Zhang J. Inversion of the vertical fault model parameters based on the gravity anomaly derivative. J ournal of the Graduate School of the Chinese Academy of Sciences ,2006 ,23( 6) :808～814 摘 　要 　断层的定量解释模型一般是一个无限延伸的二度体模型. 但当断层的深度加大或断 层延伸范围变小时 ,实际的模型与二度体模型的假设差别较大 ,致使反演结果误差较大 ,难以 适应实际问题的需要. 本文从六面体模型出发 ,将二度体模型修改为似二度体模型 ,建立了一 个带约束的重力 - 垂直断层反演模型 ,经过理论模型试验 ,证明效果要优于二度体模型的反演 结果 ,在此基础上对美国内华达州的圣安德烈斯地区的布格重力异常资料做了定量解释. 关键词 　重力异常 ,垂直断层 ,模型反演 中图分类号 　P3 1 　引言 断层或断裂带是地球表面最受关注的地质现象之一. 那些规模大、断距宽的区域性断层构造 ,能引 起明显的重力异常 ,因此重力异常被作为识别地下隐伏断层的重要标志. 重力异常的梯级带常对应断层 带的分布 , 以重力梯级带的中轴线或水平方向总梯度来划定识别断层带是传统的常用方法. 但这种传统 方法很难具体定量化确定断层产状参数. 近年来 , 国内外大多利用二维无限延伸的二度体模型来构造重力异常定量计算断层参数的反演算 法. 例如 ,S. Chacko 和 B. Bhattacharya[1 ] 提出了对斜台阶重力异常的水平二阶导数进行 Fourier 变换 ,用频 谱来确定斜台阶的水平位置、台阶倾角、上顶深度以及界面两侧的密度差 ,从理论上给出了由地质界面 [2 ] [3 ] 引起重力异常的计算公式和频率域的反演方法. 蔡宗熹 在此基础上利用侯重初 的上延叠加场思想 [4 ] 做二维倾斜断层重力异常的反演 ,提高了反演精度. Mohan 等 提出了利用 Hilbert 积分变换反演倾斜台 [5 ] 阶的理论和方法. 另外 ,顾心如 提出了利用重力异常的一阶和二阶导数来求解断层参数 ,基于时间域 的导数计算公式构造了反演公式. 但上述方法 ,在实际的应用上也存在一定的不足 ,主要表现在模型的 抗干扰能力不足 ,导致在实际资料处理中 ,难以提取断层构造信息. 为克服这种不足 ,我们从六面体模型 出发 ,将二度体模型修改为似二度体模型 ,建立了一个带约束的重力 —断层反演模型 ,经过理论模型试 验和实