Since oil exploration is reaching several kilometers below seafloor, an adequate numerical simulation for prediction of rock behavior prior to drilling is becoming increasingly important for the oil industry due to the high operating costs to ensure well stability. One of the typical problems that oil engineers face is the instability of wells in evaporitic rocks, which are often subject to excessive deformation and have time-dependent mechanical behavior characteristics. Huge oil reserves that were found around the world in evaporitic reservoirs encouraged engineers to further investigate this mechanical time-dependent response (creep). On the other hand, the concentration of stresses around the well can make it an unviable engineering project. Destabilization due to large deformations may be involved, leading to problems such as collapse and well closure. Thus, in order to study the stability of the boreholes and design an appropriate oil production system, creep must be properly considered in the expected behavior of the material. In this research, the equations that define the elasto-viscoplastic model of Sterpi and Gioda (2007) will be implemented in the computer program FLAC 3D - Itasca, code based on the finite difference method. The implementation procedure consists of generating a Dynamic Link Library (DLL) written in C (plus plus) programming language and integrating the stress along a strain increment by using an explicit integration algorithm. The elasto-viscoplastic model has been validated with experimental results from triaxial creep tests in salt rocks available in the literature, indicating a successful implementation. Finally, a numerical simulation in FLAC 3D considering the model of Sterpi and Gioda is performed in order to study a pre-salt wellbore closure.