_D地质体重磁异常无解析奇点正演计算研究.docxVIP

22009 年 8 月第 44 卷 第 4 期·非地震 ·21 5D 地质体重磁异常无解析 奇点正演计算研究骆 遥 3 ①② 姚长利 ① 薛典军 ② 金玉洁 ③( ①地质过程与矿产资源国家重点实验室和地下信息探测技术与仪器教育部重点实验室 (中国地质大学 ·北京) ,北京 100083 ; ②中国国土资源航空物探遥感中心 ,北京 100083 ; ③中国矿业大学 ,北京 100083)骆遥 ,姚长利 ,薛典军 , 金玉洁. 21 5D 地质体重磁异常无解析奇点正演计算研究. 石油地球物理勘探 , 2009 , 44( 4) :487～493摘要 本文明确指出现有文献中有关 215D 重磁异常正演模型在上半无源空间存在某些点的场值无法计算的 解析“奇点”问题 ,长方体Δ T 场及其梯度场正演计算理论表达式在上半无源空间存在场值无法计算的解析“奇 点”问题仅是 215D 体的一种特例情况 。通过详尽的理论推导 ,建立了 215D 体重磁异常在上半无源空间无解析 “奇点”正演计算模型 ,并在此基础上导出了新的长方体Δ T 场在上半无源空间无解析奇点理论表达式 。上述建 立的模型或表达式是普遍适用的 ,同现有模型或表达式相比 ,具有形式统一 、正演计算快速的优点 ,一定程度上 满足了正反演中快速计算的需求 。对比计算表明 ,所建立的 215D 体重磁异常无解析“奇点”正演计算模型和新 导出的Δ T 场无解析“奇点”理论表达式是正确的 。关键词 215D 体 重磁异常 正演 解析“奇点” 长方体1 引言21 5D 体 (即二度半体) 是指沿走向截面位置 、形 状和物性参数不变的有限延伸地质体 。限于三维地 质体建模以及重磁场正 、反演问题的复杂性[ 1～5 ] ,21 5D 体不仅可以模拟真实地质体用于地球物理重、生的原因或条件 ,并进一步找出相应的解决办法是 非常必要的 ,具有理论和实际意义。2 21 5D 重磁异常表达式及其“奇点” 问题二度半重磁异常正演模型由二度体重磁异常正[ 12～14 ]磁资料处理及解释 ,而且可以将 21 5D 体模型方便演模型修改而得,Rasmussen 等系统地建立了地加载到地震剖面上进行重力 、磁法和地震联合解 释 ,以提高地球物理资料解释的可信度 。基于 21 5D 体及 21 5D 组合体的重磁异常人机交互解释已比较 普及[ 6～11 ] ,并出现了 GM2SYS 等可视化功能很强的 21 5D 重磁场建模及正反演解释的人机交互系统。 上述正反演及实用化系统的核心问题在于 21 5D 体 的重 、磁异常场的正演计算 , 在利用现有文献中的 21 5D 体的重 、磁异常模型进行正演计算时 ,发现这 些正演计算理论表达式均存在未曾被人发现或重视 的在上半无源空间某些点的场值无法直接计算的解 析“奇点”问题 ,因此深入分析和研究解析“奇点”产基于二度半多边形截面棱柱体重磁异常模型及正演理论表达式[ 13 ] ,后来许多学者[ 6～11 ] 建立的正反演方 法及人机交互系统都以此为基础 ,其使用的正演公 式本质上具有相同的解析形式。下面对其加以说 明 ,并对计算的“奇点”问题进行分析 。设 21 5D 体棱柱模型 (图 1) 的密度为ρ, 磁化强 度为 M , 参考现有教科书[ 15 ] , 空间任一点 p ( r0 ) 引 起的重力异常和磁力异常正演公式可以表示为 重力异常N∑Δgz ( r0 ) = Gρ co sφi [ Ii ( Y2 - y0 ) - Ii ( Y1 - y0 ) ]∑i = 1(1)3 北京市海淀区学院路 29 号中国地质大学 ,100083本文于 2008 年 10 月 29 日收到 ,修改稿于 2009 年 6 月 7 日收到。本项研究得到国家 863 计划项目(2007AA06 Z134) 、(2006AA06A201) ,高等学校学科创新引智计划 (B07011) ,长江学者和创新团队发展计 划资助 。 488 石 油 地 球 物 理 勘 探 2009 年i上述式中 : G 为引力常数 ; M 为磁化强度 ; i 为棱柱 体角点标号 ; N 为棱柱体的边数 ; I0 、D0 分别为地磁 场的倾角 、偏角 ; I 、D为磁化强度方向的倾角 、偏角 ; ui = x i co sφi + z i sinφi , ui+1 = xi+1 co sφi + z i+1 sinφi ;iiwi = - xi sinφi + z i co sφ