2维和2.5维起伏地表直流电法有限差分数值模拟

起伏地表直流电场数值模拟现多采用有限元法,主要是因为其有灵活的处理曲边界的能力,然而有限元法比有限差分法要复杂,如果让有限差分法也同样具有较好的处理曲边界的能力,那数值模拟将变得更为简单.本文通过在非正则内点处采用不等距差分,在起伏地表点处直接实现边界条件,克服转移法人为改变地表形状的弊端,使得基于笛卡尔网格的有限差分法有了较强的处理曲边界的能力.通过2维和2.5维的水平地表数值解和解析解的对比,以及2维起伏地表情况下数值解和保角变换理论解的对比,均说明该方法的有效性.最后例举了一些数值计算实例.