二维非均匀网格的差分解法计算大地电磁响应