It is a complicated task to evaluate the reservoirs influenced by mud-filtrate invasion. This study presents numerical simulation results of array laterolog response and the approach for mud-filtrate invasion correction. Considering differential equations and boundary conditions of the steady current field in the rock medium, finite element method is utilized to solve the Laplace equation. Then the potential value at each node is calculated along with the apparent resistivity of six measuring modes of array laterolog. Based on the numerical results, the array laterolog response is affected by various parameters including borehole radius, mud-filtrate resistivity, mud-filtrate invasion radius and invaded zone resistivity. Correction charts for borehole and mud-filtrate invasion are developed to accurately obtain the mud-filtrate invasion depth and the true resistivity of the uninvaded formation. The correction chart has been applied to the real array laterolog data of a sandstone reservoir, which the corrected resistivity is closer to the true formation resistivity.