直流电法勘探中,三维点源正演计算复杂、计算量大,虽进行傅里叶变换可转化为二维点源,但相比二维线源计算还是复杂。基于Matlab编程,采用第一类边界条件,给出了线源直流电法二维有限元数值解法,通过与均匀半空间理论值比较,验证了计算结果的正确性;再建立高阻、低阻地电模型,并采用单边供电和双边供电分别对其响应进行比较和分析,得到了相应的视电阻率断面图。结果表明单边供电可以准确反映异常体在纵向上的位置,而双边供电更能准确反应异常体的横向位置和形态。论文研究为生产实际提供了理论依据。
2021-12-08 22:05:08
616KB
1
基于MATLAB编写的大地电磁二维有限单元法正演程序，矩形单元剖分，线性插值，
使用说明：
1. 运行MT2DMODEL.m构建正演模型
2. 运行TMmodel.m和TEmodel.m正演计算，计算结果自动保存
3. 程序进行了优化，并将BICGSTAB（稳定双共轭梯度法）用于求解线性方程组，求解效率提高，一般情况下正演计算仅需1~2秒
2021-05-17 15:20:15
57KB
1