基于ANSYSCSAMT法三维正演仿真模拟基础研究.docVIP

基于ANSYSCSAMT法三维正演仿真模拟基础研究 　　摘要：介绍了基于ANSYS软件的CSAMT法三维正演的建模及网格划分方法，研究模型参数设置与卡尼亚视电阻率之间的关系，提出了CSAMT三维正演时模型参数设置建议及规律。研究结果表明：合理设置过渡区可以提高网格划分精度，非均匀划分网格可以有效减少网格数量，提高计算速度；线圈高度的变化会影响磁场垂直入射效果；线圈高度一定时，模型空气层厚度会对视电阻率值产生影响；模型地层厚度会对视电阻率值产生很大影响，测量有效深度与模型地层厚度间存在一定比例关系，从而保证视电阻率值的准确性和规律性。 　　关键词：CSAMT法；三维正演；ANSYS仿真；临界比例；视电阻率 　　中图分类号：TP391 文献标识码：A 文章编号：1009-3044（2016）12-0228-04 　　Abstract： This paper makes an introduction about modeling and meshing of 3D CSAMT forward modeling by ANSYS， researching the relation between model parameters and Cagniard resistivity， coming up with a reference about setting of 3D CSAMT modeling parameters. The results indicate： meshing accuracy is improved by use of transition regionm， inhomogeneous meshing can effictively decrease quantity and enhance simulation rate. Coil height can influence the incoming of magnetic field which should be vertical. When coil height is certain， air thickness of the model influence apparent resistivity. Depth of the ground has a strong influence on apparent resistivity， and there is a proportion relation between surveying depth and ground depth. An appropriate ground depth can keep the simulated result accurate and regular. 　　Key words：CSAMT； 3D forward modeling； Ansys simulation； critical proportion； apparent resistivity 　　1 概述 　　CSAMT法是一种易于操作且成本较低的地质探测手段，它采用有限长线电流源作为激发源，在一定距离外布设测线接收电场和磁场，从而获得卡尼亚视电阻率，并据此判断地表下各地层的地质构造[1]。目前，一维、二维以及2.5维CSAMT法正、反演研究已经较为全面[2-5]，理论、模拟都很成熟。但是从实际问题考虑，特别是对于复杂局部地质体的勘探来说，三维CSAMT方法的正、反演才能更客观、准确地反映地质特征[6-9]。 　　通过理论方法求解三维正反演问题是相当困难，而有限元仿真在解决场问题上有着极大优势，采用有限元仿真软件ANSYS对三维地质模型进行正演模拟[10-11]，研究模型建立、网格划分、参数设置、仿真计算、数据提取的方法，为三维CSAMT法正演提供一种参考。 　　2 ANSYS仿真模拟原理 　　2.1 CSAMT法视电阻率计算 　　3 仿真模型建立及视电阻率计算 　　考虑到工程实际情况，假设模型需要实现探测100m以内的地层视电阻率，即模拟时的有效穿透深度[≥100m]。根据收发距要求，[r=（3～5）h有效≥300～500m]，选择收发距[r=500m]，模型长、宽为[1km×1.3km]；线圈内、外径[8m×12m]；初始空气层厚度设为300m，地层设为200m，空气层几乎不导电，电阻率设为[106Ω?m]，地层电阻率[ρ0=100Ω?m]，线圈施加环向电流，电流密度[Js=10Am2]。 　　3.1 线圈高度变化与磁场入射角 　　在线圈内、外径一定时，线圈高度的变化会影响磁场线入射地面时的角度，而卡尼亚视电阻率的计算要求是相互垂直的磁场和电场分量，所以在测量区应该让