The presented article demonstrates the probabilistic method based modeling of the 2D chloride ingress into reinforced concrete structures with respect to concrete heterogeneity and epoxy-coated steel reinforcement. Spatial change of concrete diffusion is assessed through the investigation of random variation of the ability of concrete to resist chloride ingress. Time-dependent chloride concentration at the reinforcement level in both homogeneous and heterogeneous models is comparatively considered taking into account of the influence of reinforcement protection as well as the defects and holidays of the coating. Expansion optimal linear estimation method is exploited to generate a random field for the structure at the mesoscale and correlation length is employed to simplify the modeling process. Preliminary analyses of the built model are conducted in both deterministic and probabilistic solutions under the scheme of the finite element method. Thus, possibility of such analyses is exploited.