Carbonate reservoirs contain a significant portion of the world's oil reserves. The uptake of oil from these reservoirs is difficult due to the complex hydro-thermo-geomechanical-geochemical processes. In recent years many efforts have been made to model these processes in fractured rocks. A few studies have been done on the stochastic nature of these processes. Geomechanical models need to take into account the reservoir rock deformations that are the result of the pores and fracture pressure, while geochemical models deal with the chemical reactions that take place into fracture and matrix with their principal processes such dissolution and precipitation. In this paper we will make the modeling from the stochastic point of view of these processes by making a model to deal with such parameters and coefficients that are stochastic, seeing them as random variables. In our knowledge this is the first work of this kind. From the geomechanical point of view we will use the stochastic nature of the Biot linear model with uncertain coefficients, while for the stochastic geochemical model we will use the Langrazhian formulation.