大地电磁（MT）和瞬变电磁（TEM）正演可以通过SimPEG开源框架实现，首先定义一个张量网格，比如用Mesh.TensorMesh创建10x10x10的网格，然后设置电导率模型，例如用np.ones填充低阻值1e-2，再修改深部为高阻1e-1，正演模拟MT问题时导入from simpeg import *和numpy，构建问题prob并求解二次场。反演则使用SimPEG的内置反演器，设置数据误差点、初始模型和参考模型，通过最小化目标函数迭代优化参数，支持多层含水模型、立方体和断层模型，直接运行开源代码即可模拟勘探地下隐蔽断层。

SimPEG框架正演示例

from simpeg import * import numpy as np # 定义网格 mesh = Mesh.TensorMesh([np.ones(10), np.ones(10), np.ones(10)]) # 定义电导率模型（例如含水层） sigma = np.ones(mesh.nC) * 1e-2 # 假设为低阻层 sigma[mesh.gridCC[:, 2] < 5] = 1e-1 # 深部高阻 # 定义正演问题（MT） probl…（代码可扩展到TEM，通过调整源和观测器模拟瞬变响应，反演用simpeg.inversion模块加载观测数据迭代求解模型参数，支持梯度下降或Gauss-Newton方法）。

开源代码推荐与使用建议

在进行综合大地电磁法（MT）和瞬变电磁法（TEM）勘探地下隐蔽断层的研究时，正演和反演是关键的步骤。为了完成这些任务，你需要选择合适的软件或开源代码来构建模型、进行正演模拟，并进行反演分析。以下是一些推荐的软件和开源代码资源，并附上使用建议与注意事项：SimPEG是最常用的Python开源框架，支持MT和TEM的全流程，从网格构建到正演反演，支持多层含水、立方体和断层模型。

MT正演详细步骤

使用SimPEG的MT模块：1. 安装SimPEG：pip install SimPEG；2. 创建3D网格mesh = TreeMesh(); mesh.insert_cells(...); 3. 定义物理属性model = np.log(sigma); 4. 设置源survey = magnetotellurics.sources.MTSource(...); 5. 运行模拟fields = problem.fieldsPair(model); dpred = problem.dpred(model); 对于TEM，切换到transient模块，定义线圈源和衰减曲线模拟瞬变响应。

TEM反演实现

TEM反演用SimPEG的EM1D或3D模块：加载实地数据data = np.load('tem_data.npy'); 设置反演器inv = Inversion.BaseInv(...); inv.run(m0); 针对断层模型，预设不连续界面参数，通过正则化约束如最小结构约束避免过拟合，支持多层含水层厚度、电导率联合反演，迭代直到数据残差小于阈值。

其他工具补充

除了SimPEG，推荐WSINV3D_MT用于MT 3D反演（Fortran代码，编译后运行），或Apex Inversion for TEM商业软件但有试用版；开源还有em1d用于1D TEM正演反演，简单易上手。对于毕业论文，SimPEG代码可直接复制修改，生成图表用matplotlib可视化模型和拟合曲线。

FAQ

Q: SimPEG支持哪些模型类型？

A: 支持多层含水模型、立方体异常体和断层模型，通过自定义网格和属性数组轻松构建。

Q: 正演计算时间长怎么优化？

A: 用TreeMesh自适应网格减少单元数，或并行计算fields求解器加速。

Q: 反演不收敛怎么办？

A: 调整beta值、初始模型平滑度，或增加数据权重和L2/Lp正则化。

Q: 有中文教程吗？

A: SimPEG官网文档英文为主，但GitHub有中文issue讨论和B站视频教程。