During HVDC earth return operation systems, a high magnitude current will be injected into soil through earth electrode, the potential on the surface would change widely and produce unfavorable effects on the AC systems around. This paper presents an effective finite element method (FEM) coupling electric field with thermal field to evaluate the electrical field induced by the injected DC current. Firstly, owe to the characteristic of FEM, this method can consider arbitrary soil and earth electrode structure. Secondly, by setting the electrical and thermal parameters of soil as a function of temperature at the same time, the dynamic coupling process of electric field and thermal field is simulated accurately. Thirdly, to deal with the singular point in FEM subdivision and the huge computation in traditional three-dimensional FEM, the FEM coupling 2-D earth electrode with 3-D soil based on "shell" theory is introduced. Finally, based on the suggested method, the effect of abnormal resistance region (ARR) near DC earth electrode on electric field distribution is analyzed.