Taking eliminating terrain effect as the main target, this paper adopts the 2D earth and 3D source terrain model and uses the weighted cosine integral method to carry out finite-element calculation in wave-number domain. To simulate complex terrain the author selected cross-symmetrical triangular grid to realize forward calculation of 2D CSAMT data for the electrical equatorial dipole array. On the basis of the 2D forward calculation, a set of rolling terrain 2D CSAMT inversion technique was developed based on the OCCAM method. Through theoretical calculation and real data processing, the terrain effects could be effectively corrected.