This research utilizes the imbibition experiments and X-ray tomography results for modeling fluid flow in naturally fractured reservoirs. Conventional dual porosity simulation requires large number of runs to quantify transfer function parameters for history matching purposes. In this study empirical transfer functions (ETF) are derived from imbibition experiments and this allows reduction in the uncertainness in modeling of transfer of fluids from the matrix to the fracture. The application of the ETF approach is applied in two phases. In the first phase, imbibition experiments are numerically solved using the diffusivity equation with different boundary conditions. Usually only the oil recovery in imbibition experiments is matched. But with the advent of X-ray CT, the spatial variation of the saturation can also be computed. The matching of this variation can lead to accurate reservoir characterization. In the second phase, the imbibition derived empirical transfer functions are used in developing a dual porosity reservoir simulator. The results from this study are compared with published results. The study reveals the impact of uncertainty in the transfer function parameters on the flow performance and reduces the computations to obtain transfer function required for dual porosity simulation.

Naturally fractured reservoirs (NFR) represent an important percentage of the worldwide hydrocarbon reserves and production. Reservoir simulation is a fundamental technique in characterizing this type of reservoir. Fracture properties are often not available due to difficulty to characterize the fracture system. On the other hand, well test analysis is a well known and widely applied reservoir characterization technique. Well testing in NFR provides two characteristic parameters, storativity ratio and interporosity flow coefficient. The storativity ratio is related to fracture porosity. The interporosity flow coefficient can be linked to shape factor, which is a function of fracture spacing. The purpose of this work is to investigate the feasibility of estimating fracture porosity and fracture spacing from single well test analysis and to evaluate the use of these two parameters in dual porosity simulation models. The following assumptions were considered for this research: 1) fracture compressibility is equal to matrix compressibility; 2) no wellbore storage and skin effects are present; 3) pressure response is in pseudo-steady state; and 4) there is single phase flow. Various simulation models were run and build up pressure data from a producer well was extracted. Well test analysis was performed and the result was compared to the simulation input data. The results indicate that the storativity ratio provides a good estimation of the magnitude of fracture porosity. The interporosity flow coefficient also provides a reasonable estimate of the magnitude of the shape factor, assuming that matrix permeability is a known parameter. In addition, pressure tests must exhibit all three flow regimes that characterizes pressure response in NFR in order to obtain reliable estimations of fracture porosity and shape factor.

This research utilizes the imbibition experiments and X-ray tomography results for modeling fluid flow in naturally fractured reservoirs. Conventional dual porosity simulation requires large number of runs to quantify transfer function parameters for history matching purposes. In this study empirical transfer functions (ETF) are derived from imbibition experiments and this allows reduction in the uncertainness in modeling of transfer of fluids from the matrix to the fracture. The application of the ETF approach is applied in two phases. In the first phase, imbibition experiments are numerically solved using the diffusivity equation with different boundary conditions. Usually only the oil recovery in imbibition experiments is matched. But with the advent of X-ray CT, the spatial variation of the saturation can also be computed. The matching of this variation can lead to accurate reservoir characterization. In the second phase, the imbibition derived empirical transfer functions are used in developing a dual porosity reservoir simulator. The results from this study are compared with published results. The study reveals the impact of uncertainty in the transfer function parameters on the flow performance and reduces the computations to obtain transfer function required for dual porosity simulation.

This thesis is aimed to calculate the effective permeability tensor and to simulate the fluid flow in naturally fractured reservoirs. This requires an understanding of the mechanisms of fluid flow in naturally fractured reservoirs and the detailed properties of individual fractures and matrix porous media. This study has been carried out to address the issues and difficulties faced by previous methods; to establish possible answers to minimise the difficulties; and hence, to improve the efficiency of reservoir simulation through the use of properties of individual fractures. The methodology used in this study combines several mathematical and numerical techniques like the boundary element method, periodic boundary conditions, and the control volume mixed finite element method. This study has contributed to knowledge in the calculation of the effective permeability and simulation of fluid flow in naturally fractured reservoirs through the development of two algorithms. The first algorithm calculates the effective permeability tensor by use of properties of arbitrary oriented fractures (location, size and orientation). It includes all multi-scaled fractures and considers the appropriate method of analysis for each type of fracture (short, medium and long). In this study a characterisation module which provides the detail information for individual fractures is incorporated. The effective permeability algorithm accounts for fluid flows in the matrix, between the matrix and the fracture and disconnected fractures on effective permeability. It also accounts for the properties of individual fractures in calculation of the effective permeability tensor. The second algorithm simulates flow of single-phase fluid in naturally fractured reservoirs by use of the effective permeability tensor. This algorithm takes full advantage of the control volume discretisation technique and the mixed finite element method in calculation of pressure and fluid flow velocity in each grid block. It accounts for the continuity of flux between the neighbouring blocks and has the advantage of calculation of fluid velocity and pressure, directly from a system of first order equations (Darcy???s law and conservation of mass???s law). The application of the effective permeability tensor in the second algorithm allows us the simulation of fluid flow in naturally fractured reservoirs with large number of multi-scale fractures. The fluid pressure and velocity distributions obtained from this study are important and can considered for further studies in hydraulic fracturing and production optimization of NFRs.

Naturally fractured reservoirs (NFRs) hold a significant amount of the world's hydrocarbon reserves. Compared to conventional reservoirs, NFRs exhibit a higher degree of heterogeneity and complexity created by fractures. The importance of fractures in production of oil and gas is not limited to naturally fractured reservoirs. The economic exploitation of unconventional reservoirs, which is increasingly a major source of short- and long-term energy in the United States, hinges in part on effective stimulation of low-permeability rock through multi-stage hydraulic fracturing of horizontal wells. Accurate modeling and simulation of fractured media is still challenging owing to permeability anisotropies and contrasts. Non-physical abstractions inherent in conventional dual porosity and dual permeability models make these methods inadequate for solving different fluid-flow problems in fractured reservoirs. Also, recent approaches for discrete fracture modeling may require large computational times and hence the oil industry has not widely used such approaches, even though they give more accurate representations of fractured reservoirs than dual continuum models. We developed an embedded discrete fracture model (EDFM) for an in-house fully-implicit compositional reservoir simulator. EDFM borrows the dual-medium concept from conventional dual continuum models and also incorporates the effect of each fracture explicitly. In contrast to dual continuum models, fractures have arbitrary orientations and can be oblique or vertical, honoring the complexity and heterogeneity of a typical fractured reservoir. EDFM employs a structured grid to remediate challenges associated with unstructured gridding required for other discrete fracture models. Also, the EDFM approach can be easily incorporated in existing finite difference reservoir simulators. The accuracy of the EDFM approach was confirmed by comparing the results with analytical solutions and fine-grid, explicit-fracture simulations. Comparison of our results using the EDFM approach with fine-grid simulations showed that accurate results can be achieved using moderate grid refinements. This was further verified in a mesh sensitivity study that the EDFM approach with moderate grid refinement can obtain a converged solution. Hence, EDFM offers a computationally-efficient approach for simulating fluid flow in NFRs. Furthermore, several case studies presented in this study demonstrate the applicability, robustness, and efficiency of the EDFM approach for modeling fluid flow in fractured porous media. Another advantage of EDFM is its extensibility for various applications by incorporating different physics in the model. In order to examine the effect of pressure-dependent fracture properties on production, we incorporated the dynamic behavior of fractures into EDFM by employing empirical fracture deformation models. Our simulations showed that fracture deformation, caused by effective stress changes, substantially affects pressure depletion and hydrocarbon recovery. Based on the examples presented in this study, implementation of fracture geomechanical effects in EDFM did not degrade the computational performance of EDFM. Many unconventional reservoirs comprise well-developed natural fracture networks with multiple orientations and complex hydraulic fracture patterns suggested by microseismic data. We developed a coupled dual continuum and discrete fracture model to efficiently simulate production from these reservoirs. Large-scale hydraulic fractures were modeled explicitly using the EDFM approach and numerous small-scale natural fractures were modeled using a dual continuum approach. The transport parameters for dual continuum modeling of numerous natural fractures were derived by upscaling the EDFM equations. Comparison of the results using the coupled model with that of using the EDFM approach to represent all natural and hydraulic fractures explicitly showed that reasonably accurate results can be obtained at much lower computational cost by using the coupled approach with moderate grid refinements. / text

This thesis is aimed to calculate the effective permeability tensor and to simulate the fluid flow in naturally fractured reservoirs. This requires an understanding of the mechanisms of fluid flow in naturally fractured reservoirs and the detailed properties of individual fractures and matrix porous media. This study has been carried out to address the issues and difficulties faced by previous methods; to establish possible answers to minimise the difficulties; and hence, to improve the efficiency of reservoir simulation through the use of properties of individual fractures. The methodology used in this study combines several mathematical and numerical techniques like the boundary element method, periodic boundary conditions, and the control volume mixed finite element method. This study has contributed to knowledge in the calculation of the effective permeability and simulation of fluid flow in naturally fractured reservoirs through the development of two algorithms. The first algorithm calculates the effective permeability tensor by use of properties of arbitrary oriented fractures (location, size and orientation). It includes all multi-scaled fractures and considers the appropriate method of analysis for each type of fracture (short, medium and long). In this study a characterisation module which provides the detail information for individual fractures is incorporated. The effective permeability algorithm accounts for fluid flows in the matrix, between the matrix and the fracture and disconnected fractures on effective permeability. It also accounts for the properties of individual fractures in calculation of the effective permeability tensor. The second algorithm simulates flow of single-phase fluid in naturally fractured reservoirs by use of the effective permeability tensor. This algorithm takes full advantage of the control volume discretisation technique and the mixed finite element method in calculation of pressure and fluid flow velocity in each grid block. It accounts for the continuity of flux between the neighbouring blocks and has the advantage of calculation of fluid velocity and pressure, directly from a system of first order equations (Darcy???s law and conservation of mass???s law). The application of the effective permeability tensor in the second algorithm allows us the simulation of fluid flow in naturally fractured reservoirs with large number of multi-scale fractures. The fluid pressure and velocity distributions obtained from this study are important and can considered for further studies in hydraulic fracturing and production optimization of NFRs.

As easy target reservoirs are depleted around the world, the need for intelligent enhanced oil recovery (EOR) methods increases. The first part of this work is focused on modeling aspects of novel chemical EOR methods for naturally fractured reservoirs (NFR) involving wettability modification towards more water wet conditions. The wettability of preferentially oil wet carbonates can be modified to more water wet conditions using alkali and/or surfactant solutions. This helps the oil production by increasing the rate of spontaneous imbibition of water from fractures into the matrix. This novel method cannot be successfully implemented in the field unless all of the mechanisms involved in this process are fully understood. A wettability alteration model is developed and implemented in the chemical flooding simulator, UTCHEM. A combination of laboratory experimental results and modeling is then used to understand the mechanisms involved in this process and their relative importance. The second part of this work is focused on modeling surfactant/polymer floods using a fully implicit scheme. A fully implicit chemical flooding module with comprehensive oil/brine/surfactant phase behavior is developed and implemented in general purpose adaptive simulator, GPAS. GPAS is a fully implicit, parallel EOS compositional reservoir simulator developed at The University of Texas at Austin. The developed chemical flooding module is then validated against UTCHEM. / text

The primary objective of this study is to improve understanding of the mechanics of hydraulic fracturing in naturally-fractured reservoirs. The study focuses on enhancing the interpretation of hydraulic fracture-induced microseismic data using an S-wave Gaussian-beam method and numerical modeling techniques for interpretation. The S-wave Gaussian-beam method was comprehensively calibrated by synthetic and real data sets with different recording networks, and this showed the potential to retrieve additional microseismic data from hydraulic fracturing with linear receiver arrays. This approach could enhance current practice because a large number of induced events in these environments have very strong S-waves with P-wave amplitudes similar, or less than, background noise levels. The numerical study using the distinct element methods PFC2D and PFC3D was used to validate the understanding of the hydraulic fracturing mechanisms induced in laboratory and field fluid treatments in naturally-fractured reservoirs. This was achieved through direct comparison with the results of the geometry of hydraulic fractures and seismic source information (locations, magnitudes, and mechanisms) from both laboratory experiments and field observations. A suite of numerical models with fully-dynamic and hydro-mechanical coupling has been used to examine in detail the interaction between natural and induced fractures with the variations of the differential stresses and the orientations of the pre-fractures, and the relationship between the fluid front, the fracture tip, and the induced seismicity. The numerical results qualitatively agreed with the laboratory and field observations of the geometry of hydraulic fractures, confirmed the possible mechanics of new fracture development and their interactions with natural fractures, and illustrated the possible relationship between the fluid front and the fracture tip. The validated model could therefore help track the potential extent of induced fracturing in naturally-fractured reservoirs and the extent to which it can be detected by a microseismic monitoring array in order to assess the effectiveness of a hydraulic fracturing project.

The primary objective of this study is to improve understanding of the mechanics of hydraulic fracturing in naturally-fractured reservoirs. The study focuses on enhancing the interpretation of hydraulic fracture-induced microseismic data using an S-wave Gaussian-beam method and numerical modeling techniques for interpretation. The S-wave Gaussian-beam method was comprehensively calibrated by synthetic and real data sets with different recording networks, and this showed the potential to retrieve additional microseismic data from hydraulic fracturing with linear receiver arrays. This approach could enhance current practice because a large number of induced events in these environments have very strong S-waves with P-wave amplitudes similar, or less than, background noise levels. The numerical study using the distinct element methods PFC2D and PFC3D was used to validate the understanding of the hydraulic fracturing mechanisms induced in laboratory and field fluid treatments in naturally-fractured reservoirs. This was achieved through direct comparison with the results of the geometry of hydraulic fractures and seismic source information (locations, magnitudes, and mechanisms) from both laboratory experiments and field observations. A suite of numerical models with fully-dynamic and hydro-mechanical coupling has been used to examine in detail the interaction between natural and induced fractures with the variations of the differential stresses and the orientations of the pre-fractures, and the relationship between the fluid front, the fracture tip, and the induced seismicity. The numerical results qualitatively agreed with the laboratory and field observations of the geometry of hydraulic fractures, confirmed the possible mechanics of new fracture development and their interactions with natural fractures, and illustrated the possible relationship between the fluid front and the fracture tip. The validated model could therefore help track the potential extent of induced fracturing in naturally-fractured reservoirs and the extent to which it can be detected by a microseismic monitoring array in order to assess the effectiveness of a hydraulic fracturing project.

A significant amount of Viscous Oil (e.g., heavy oil, extra heavy oil, and bitumen) is trapped in Naturally Fractured Carbonate Reservoirs also known as NFCRs. The word VO endowment in NFCRs is estimated at ~ 2 Trillion barrels mostly reported in Canada, the USA, Russia, and the Middle East. To date, contributions to the world daily oil production from this immense energy resource remains negligible mainly due to the lack of appropriate production technologies. Implementation of a VO production technology such as steam injection is expensive (high capital investment), time-consuming, and people-intensive. Hence, before selecting a production technology for detailed economic analysis, use of cursory or broad screening tools or guides is a convenient means of gaining a quick overview of the technical feasibility of the various possible production technologies applied to a particular reservoir. Technical screening tools are only available for the purpose of evaluation of the reservoir performance parameters in oil sands for various thermal VO exploitation technologies such as Steam Assisted Gravity Drainage (SAGD), Cyclic Steam Stimulation (CSS), Horizontal well Cyclic steam Stimulation (HCS), and so on. Nevertheless, such tools are not applicable for VO NFCRs assessment without considerable modifications due to the different nature of these two reservoir types (e.g., presence and effects of fracture network on reservoir behavior, wettability, lithology, fabric, pore structure, and so on) and also different mechanisms of energy and mass transport. Considering the lack of robust and rapid technical reservoir screening tools for the purpose of quick assessment and performance prediction for VO NFCRs under thermal stimulation (e.g., steamflooding), developing such fast and precise tools seems inevitable and desirable. In this dissertation, an attempt was made to develop new screening tools for the purpose of reservoir performance prediction in VO NFCRs using all the field and laboratory available data on a particular thermal technology (vertical well steamflooding). Considering the complex and heterogeneous nature of the NFCRs, there is great uncertainty associated with the geological nature of the NFCRs such as fracture and porosity distribution in the reservoir which will affect any modeling tasks aiming at modeling of processes involved in thermal VO production from these types of technically difficult and economically unattractive reservoirs. Therefore, several modeling and analyses technqiues were used in order to understand the main parameters controlling the steamflooding process in NFCRs and also cope with the uncertainties associated with the nature of geologic, reservoir and fluid properties data. Thermal geomechanics effects are well-known in VO production from oil sands using thermal technologies such as SAGD and cyclic steam processes. Hence, possible impacts of thermal processes on VO NFCRs performance was studied despite the lack of adequate field data. This dissertation makes the following contributions to the literature and the oil industry: Two new statistical correlations were developed, introduced, and examined which can be utilized for the purpose of estimation of Cumulative Steam to Oil Ratio (CSOR) and Recovery Factor (RF) as measures of process performance and technical viability during vertical well steamflooding in VO Naturally Fractured Carbonate Reservoirs (NFCRs). The proposed correlations include vital parameters such as in situ fluid and reservoir properties. The data used are taken from experimental studies and also field trials of vertical well steamflooding pilots in viscous oil NFCRs reported in the literature. The error percentage for the proposed correlations is < 10% for the worst case and contains fewer empirical constants compared with existing correlations for oil sands. The interactions between the parameters were also considered. The initial oil saturation and oil viscosity are the most important predictive factors. The proposed correlations successfully predicted steam/oil ratios and recovery factors in two heavy oil NFCRs. These correlations are reported for the first time in the literature for this type of VO reservoirs. A 3-D mathematical model was developed, presented, and examined in this research work, investigating various parameters and mechanisms affecting VO recovery from NFCRs using vertical well steamflooding. The governing equations are written for the matrix and fractured medium, separately. Uncertainties associated with the shape factor for the communication between the matrix and fracture is eliminated through setting a continuity boundary condition at the interface. Using this boundary condition, the solution method employed differs from the most of the modeling simulations reported in the literature. A Newton-Raphson approach was also used for solving mass and energy balance equations. RF and CSOR were obtained as a function of steam injection rate and temperature and characteristics of the fractured media such as matrix size and permeability. The numerical solution clearly shows that fractures play an important role in better conduction of heat into the matrix part. It was also concluded that the matrix block size and total permeability are the most important parameters affecting the dependent variables involved in steamflooding. A hybrid Artificial Neural Network model optimized by co-implementation of a Particle Swarm Optimization method (ANN-PSO) was developed, presented, and tested in this research work for the purpose of estimation of the CSOR and RF during vertical well steamflooding in VO NFCRs. The developed PSO-ANN model, conventional ANN models, and statistical correlations were examined using field data. Comparison of the predictions and field data implies superiority of the proposed PSO-ANN model with an absolute average error percentage < 6.5% , a determination coefficient (R2) > 0.98, and Mean Squared Error (MSE) < 0.06, a substantial improvement in comparison with conventional ANN model and empirical correlations for prediction of RF and CSOR. This indicates excellent potential for application of hybrid PSO-ANN models to screen VO NFCRs for steamflooding. This is the first time that the ANN technique has been applied for the purpose of performance prediction of steamflooding in VO NFCRs and also reported in the literature. The predictive PSO-ANN model and statistical correlations have strong potentials to be merged with heavy oil recovery modeling softwares available for thermal methods. This combination is expected to speed up their performance, reduce their uncertainty, and enhance their prediction and modeling capabilities. An integrated geological-geophysical-geomechanical approach was designed, presented, and applied in the case of a NFCR for the purpose of fracture and in situ stresses characterization in NFCRs. The proposed methodology can be applied for fracture and in situ stresses characterization which is beneficial to various aspects of asset development such as well placement, drilling, production, thermal reservoir modeling incorporating geomechanics effects, technology assessment and so on. A conceptual study was also conducted on geomechanics effects in VO NFCRs during steamflooding which is not yet well understood and still requires further field, laboratory, and theoretical studies. This can be considered as a small step forward in this area identifying positive potential of such knowledge to the design of large scale thermal operations in VO NFCRs.