台阶（断层）模型重力异常正演计算的一种新算法.pdfVIP

第33卷 第2期 世 界 地 质 Vo1．33 No．2 2014年 6月 GL0BAL GE0L0GY Jun．2Ol4 文章编号：1004— 5589 2014 02—0477—07 台阶 断层 模型重力异常正演计算的一种新算法 冯v—兰—…天 吉林大学 地球探测科学与技术学院，长春 130026 摘要：传统的台阶模型重力异常正演公式会随着台阶倾斜面与水平面夹角的变化而产生畸变现象。为 克服这一局限性，在前人推导台阶模型重力异常正演表达式的基础上，重新对台阶模型重力异常正演 公式进行 了理论推导，得 出全新的台阶模型重力异常正演公式，并证明传统公式仅为新公式的特殊情 况。模型实验证明新的正演公式有效克服 了传统正演公式产生畸变的问题。 关键词：台阶模型；重力场；正演；积分运算 ；断层 中图分类号 ：P631．1 文献标识码 ：A doi：10．3969／j．issn．1004—5589．2014．02．026 A new algorithm forgravityanomalyforward calculation ofstepped fault model FENG Lan．tian CollegeofGeo—explorationScienceandTechnology，JilinUniversity，Changchun130026，China Abstract：Traditionalfl~rwardformulaofsteppedmodelforgravity anomalymaycausedistortion asthestep includedangleofitsinclinedplaneandhorizontalplanechanged．Toovercomethislimitation，basedontheprevi— OUSderivationaboutthegravityanomalyforwardexpression，theauthorderivesthisformularenewedlyandgetsa newformandfurthermoreprovesthatthetraditionalformulaisjustonespecialcaseofthenew．Themodelexperi— mentsprovethatthenew equationcanconquerthedistortionproblem efficiently． Keywords： steppedmodel；gravityfield；forward；integraloperation；fault 用 j。 0 引言 目前国内编写出版的大部分重力勘探类教材 重力场的正演在重力异常解释 中有着重要的意 中，关于台阶模型重力异常的理论表达式都是相 同 义，是研究重力场源与重力场分布特征之问关系的 的。在进行正演模拟过程中存在一定的局限性，主 基础。台阶模型是重力异常正演 中常用的基本模 要反映在台阶倾斜面与水平面夹角 简称倾斜角， 型，同时也是有限延伸板状体、无限长三棱柱体等 下同 取值变化到一定角度或台阶厚度与台阶上 常见二度体模型的基础。自然界中的断层、接触带 顶面深度比例发生变化时，所得到的重力异常会发 和超覆岩层等实际地质体可以抽象成台阶模型，断 生畸变。针对这一局限性 ，笔者推导了一种新的台 裂构造可以抽象成两组台阶模型的组合。因此，在 阶模型重力异常正演公式，克服 了现有公式的局 断裂构造位置识别、深度反演等重力异常地质解释 限性。 中，台阶模型及其重力 Ⅱ演公式有着广泛 的应 收稿 日期：2013-o3-o5；改回日期：2013-124 2 基金项 目：国家级大学生创新性实验项 目 2012A62156 478 世 界 地 质 第33卷 1 台阶模型重力异常 △g 2J Ⅳ 戋 根据万有引力定律和定积分基本原理 ，如 推导过程实质是对上式进行积分过程。为了推 图1所示密度均匀的台阶所产生的重力场表达式 导方便，首先给出推导过程中所用公式 ¨。】： 为： Jf ：arefan+c Zig 22rL 舞 1 14-2一 “ “ 。 3 r dx 2 2ax -4-b ． 式中：G为万有引力常数；P为场源与围岩的 J “ 密度差； ， ， 为观测点坐标； ，叩， 4 为台阶场源点坐标；h为台阶上顶面埋深；H为台 阶下顶面深度 ； h，日 为 的积分限；Ol为台阶 』 n 【+ +c一b2+音一篆] 倾斜角 0 叮T ，而 一~xota，∞ 为 的积 分限。 一 』 ㈩ 对 1 式进行积分求解，即可得到台阶模型 重力异常的理论表达式，现有大部分教