It is important to clarify activities of geological structures until earthquakes occur.In this paper, we propose a theory for numerical analysis to predict initiation and development of faults.The formulation is based upon the multiscale analysis and the faults are treated as localized deformations.Considering high order terms in the singular perturbation method, we deduce a nonlocal constitutive law, such that stresses are determined by not only strains but also strain gradients, whose effect dissolve a problem that the numerical results depend on the mesh size in finite element analysis.