Fault zones are geological structures usually present in subsurface. In the oil and gas industry, such structures can have different impacts on geomechanical behavior and on reservoir production. In general, fault zones are composed by a core and a damage zone. The core is responsible for the compartmentalization of reservoirs, acting usually as a barrier for fluid flow. The damage zone is the region of rock deformed adjacent to the core that may directly impact on production, creating preferential flow paths or barriers, depending on the geological features triggered in the deformation process. While the core width can be in the order of centimeters to a few meters, the damage zone width is uncertain. This work presents a methodology for numerical modeling of the generation of damage zones. The methodology is based on finite element method, elastoplastic constitutive models representative of the geomechanical behavior of rocks, and on the application of prescribed displacements to simulate the fault formation process. Several scenarios are analyzed considering the impact of the geomechanical parameters of rocks and the relative distances between two faults to characterize the damage zones. Then, the methodology is used to analyze regions favorable to the trajectory of a production well located between three geological faults. The obtained results demonstrate that the proposed methodology can serve as a basis for characterizing damage zones in geological faults.