Increasing of the oil recovery is an essential task of the reservoir engineers. Many highly fractured carbonate reservoirs with low matrix porosity and low matrix permeability in the world have been depleted with low ultimate recovery due to the improper management. To obtain higher oil recovery and an optimum oil production, a three dimensional geological model for the reservoir is needed. The reservoir model should be calibrated via the history matching process which makes the model reliable for forecasting and development planning. An investigation has been carried out to find the most sensitive parameter which affects the matching between the real production data and the simulated production data. After the calibration of the model, some prediction scenarios have been run to realize the future performance of the reservoir. Three wells have been suggested and included in the simulation as producers in some prediction cases and as gas injectors in another other case. The results show that the three producers with low production rate will assist in increasing the recovery and by converting those to gas injectors will assist more especially when the oil production rate is higher. An oil field in Kurdistan region has been selected for this research.

The anisotropic behavior of a rock mass with persistent and planar joint sets is mainly governed by the geometrical and mechanical characteristics of the joints. The aim of the study is to develop a continuum-based approach for simulation of multi jointed geomaterials. There are two available numerical techniques for the strain-stress analysis of rock masses: continuum-based methods and discontinuum based methods. Joints are simulated explicitly in discontinuous methodology. This technique provides a more accurate description for the behavior of a rock mass. However, in some projects, the explicit definition becomes impractical, especially with increasing number of joints. Besides, the calculation efficiency will be significant reduced as the number of joints increases within the model. Considering the above mentioned shortcomings of the discontinuous method, the continuum-based approach is widely used in rock mechanics. Within the continuum methods, the discontinuities are regarded as smeared cracks in an implicit manner and all the joint parameters are incorporated into the equivalent constitutive equations. A new equivalent continuum model, called multi-joint model, is developed for jointed rock masses which may contain up to three arbitrary persistent joint sets. The Mohr-Coulomb yield criterion is used to check failure of the intact rock and the joints. The proposed model has solved the issue of multiple plasticity surfaces involved in this approach combined with multiple failure mechanisms. The multi-joint model is implemented into FLAC and is verified against the distinct element method (UDEC), analytical solutions, and experimental data. Uniaxial compression tests with artificial rock-like material (gypsum) are carried out in the laboratory in order to verify the developed constitutive model and to investigate the behavior of jointed specimen. Samples with two crossing joints covering more than 20 angle configurations and two different property sets were prepared and tested. Simulation results are in good agreement with experimental observations. The developed model is applied to two potential practical applications: the stability analysis of a slope and a tunnel under different stress conditions. Finally, the main achievements of the whole PhD study are summarized and future research work is proposed.