三维地电模型电磁响应有限元模拟.pdfVIP

三维地电模型电磁响应的有限元模拟 徐志锋 吴小平 中国科学技术大学地球与空间科学学院，合肥，230026 tO3-D model Modelingelectromagneticresponsegeoelectrical ’· finiteelementmethods by Xu Wu zhifeng xiaoping School and and science,UniversityofsciencetechnologyofChina,Heifei ofearthspace ThreedimensionalGalerkinfilliteelementmethodhasbeenusedtosolvecontrolledsourceaudio in media．Theisbasedon frequencymagnetotellurics(CSAMT)problemsinhomogeneousconducting solution theweaksolutionofthe for linear Wassolvedthe equationmagneticpotential．Thesparse equation by Bi-CGSTAB at can the fieldfinite values method；aftergottenpotential verynodes，weget electromagneticby differencemethod．Itshowsthattherelativeelrorislessthan calculatedresults． quotient 4％by 可控音频率大地电磁法(CSAMT)在资源劫查、深部构造研究中应用的广度和深度都在增加。早期， 可控音频大地电磁法的理论电磁响应模型都是基于层状地电模型的响应，随着计算机的发展，人们利用 计算机开展了三维地电模型电磁响应的数值模拟，特别是二维地电模型的电磁响应：但是目前我国金属 矿勘探状况是，面临在地形起伏、构造复杂的山区进行工作，则必须直接面对三维地电模型电磁响应的 问题，这就给可控音频率大地电磁法的三维地电模型正演计算提出了更高的要求，以适应当前的情况。 设时间因子的形式如下： expl甜 (1) 则频率域的Maxwell方程组形式如下： v．E：丝 ￡ V×E—f翻B=0 V×H=erE+，，一f枷 V．B=0 H=V×A 设： E=f掣A—V矽 则可以得到： VV·A-V2A=k2A-(a-ioJe)Vq， V·(faczA—V咖=0 其中： k2=彩2e／z—ia乒zo 为了保证数值稳定，在式(4)第一式的左边加上一项一VV