Electrochemical reactions at the solid- electrolyte interface cause the conduction current through it and create an imperfect polarization phenomenon. This affects the magnitude of induced polarization (IP) response. The IP parameters, such as chargeability, phase difference, in-phase and imaginary conductivity are related to petrophysical, physicochemical and thermodynamic properties of porous media. So, eliminating the reactions, as it was done in classic IP models, forms the errors in the estimation of these properties, and affects the IP capability. In this paper, a complex conductivity model was introduced, and a 2D numerical model for potential and ion concentration in the electrical double layer was created by Poison-Nernst-Planck solution and considering the electrochemical reactions at the solid- electrolyte interface to simulate the dominant physicochemical and petrophysical conditions in subsurface porous media. Electrochemical IP modeling is great progress to a better understanding of the IP phenomenon and developing the IP method.