基于二次场的二维大地电磁有限元法数值模拟.docVIP

第35卷第8 期 2007年8月同济大学学报(自然科学版)JOURNALOFTONGJIUNIVERSITY(NATURALSCIENCE)Vol.35No.8　Aug.2007 基于二次场的二维大地电磁有限元法数值模拟 刘小军,王家林,于　鹏 (同济大学海洋地质国家重点实验室,上海　200092) 摘要:引入基于计算二次场响应的方法,将其用于二维大地电磁数值模拟.导出了二维大地电磁二次场的表达式,在使用有限元法解偏微分方程时,利用二次场偏微分方程与总场偏微分方程的相似性,直接采用总场法中的系数矩阵作为二次场的合成矩阵.将根据此原理所编制的二次场有限元法模拟软件,对几个典型模型进行了试算,并与传统总场模拟法作了比较分析.结果表明,直接计算二次场的结果精度更高,更接近真实解. 关键词:大地电磁;有限元法;总场;二次场;数值模拟 中图分类号:P631　　　　　　　文献标识码:A　　　　　　　　文章编号:0253-374X(2007)08-1113-05 SecondaryField2BasedTwo2NumericalModelingbyunnUPeng (Geology,TongjiUniversity,Shanghai200092,China) Abstract:Thispaperpresentsanefficientmethodnamedsecondaryfiniteelementmethod(FEM)fortwo2dimensionalmagnetotelluricnumericalmodeling,whichprovidessolutionforsecondaryvariationsinthefield.Todeducetherationaleofthenewprocedure,thetotalfieldisdividedintoprimaryandsecondaryfield.Onlythesecondaryfieldequationswhichneedtobesolvedfortheprimaryfieldcanbeobtainedeasily.Theequationsofthesecondaryfieldcanbesolvedeffortlesslybyfiniteelementmethodwhentheboundaryconditionsaresimple.TheeffectofresolutionisanalysedtostudytheresponseofsyntheticmodelbythesecondaryfieldandacomparisonisconductedbetweentheresultsofsecondaryFEMandthetotalfieldFEM.Itisconcludedthatthenewarithmeticofnumericalmodelingcanbeusedtoderivemorepreciseresultsincalculatingthesecondaryfieldandisalmostnotaffectedbythedensityofmesh.Itisexpectedtobeusedwidelyfortheinversionofmagnetotelluricdata. Keywords:magnetotelluric;finiteelementmethod;totalfield;secondaryfield;numericalmodeling　　有限元法作为一种解偏微分方程的数值方法, 已广泛应用于各种地球物理方法的正演模拟. Coggon[1]最早提出将有限元法应用电磁场与激发的模拟效果,随后Wannamaker等人[2]将其应用到大地电磁场的模拟中,并提出了带地形的模拟方法.但有限元法是用非无限小的单元和有限网格来模拟 无限空间,本身就是一种近似.因此,在对二维大地极化场的模拟,并讨论了二维带地形的激发极化场 收稿日期:2006-06-06 基金项目:中科院知识创新工程资助项目(KZCX1-SW-18-01) 作者简介:刘小军(1979-),男,湖北潜江人,博士生.E2mail:joneslxj@163.com 　　　1114 同济大学学报(自然科学版)第35卷　 电磁场模拟时有一些固有缺点[3].国内很多学者提出一些改善计算精度和提高计算速度的方法,黄俊革等人[4]提出了计算异常电位的有限元激发极化数值模拟,采用异常电位法计算节点电位,保证了计算精度;欧阳联华等人[5]通过改变插值函数和剖分方式来提高解的精确度,并且详细讨论了选择不同插值函数和不同网格密度计算时对电磁场值的影响;陈小斌等人[6]提出直接迭代算法,直接在计算网格的基本结构上,由泛函求极值过程中直接