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.

We propose a rigorous approach for well drainage volume calculations in gas reservoirs based on the flux field derived from dual porosity finite-difference simulation and demonstrate its application to optimize well placement. Our approach relies on a high frequency asymptotic solution of the diffusivity equation and emulates the propagation of a 'pressure front' in the reservoir along gas streamlines. The proposed approach is a generalization of the radius of drainage concept in well test analysis (Lee 1982), which allows us not only to compute rigorously the well drainage volumes as a function of time but also to examine the potential impact of infill wells on the drainage volumes of existing producers. Using these results, we present a systematic approach to optimize well placement to maximize the Estimated Ultimate Recovery. A history matching algorithm is proposed that sequentially calibrates reservoir parameters from the global-to-local scale considering parameter uncertainty and the resolution of the data. Parameter updates are constrained to the prior geologic heterogeneity and performed parsimoniously to the smallest spatial scales at which they can be resolved by the available data. In the first step of the workflow, Genetic Algorithm is used to assess the uncertainty in global parameters that influence field-scale flow behavior, specifically reservoir energy. To identify the reservoir volume over which each regional multiplier is applied, we have developed a novel approach to heterogeneity segmentation from spectral clustering theory. The proposed clustering can capture main feature of prior model by using second eigenvector of graph affinity matrix. In the second stage of the workflow, we parameterize the high-resolution heterogeneity in the spectral domain using the Grid Connectivity based Transform to severely compress the dimension of the calibration parameter set. The GCT implicitly imposes geological continuity and promotes minimal changes to each prior model in the ensemble during the calibration process. The field scale utility of the workflow is then demonstrated with the calibration of a model characterizing a structurally complex and highly fractured reservoir.

The oil production from well fractured carbonate reservoirs is a considerable part of the total oil production in the world. The petroleum resource base in naturally fractured reservoirs is estimated to be in the range of billions of barrels in the U.S and in addition, a multibillion- barrel international oil resource base exists in naturally fractured reservoirs. Gravity drainage is important in some of oil recovery processes, either acting as the driving force in processes using horizontal wells or altering the displacement patterns during water-flooding, chemical flooding, CO2 flooding and other EOR methods. The gravity drainage process has a major effect on oil recovery from oil reservoirs. Gravity drainage driven oil production in naturally fractured and other complex reservoirs falls into two regimes: the balk flow regime and the film flow regime. Oil recovery by gravity drainage in a fractured reservoir strongly depends on the capillary height of the porous medium. Capillarity and gravity forces are usually the major driving forces in fractured reservoirs. This PhD thesis consists of two main parts namely: 1) Experimental works on gravity drainage, and 2) Modeling and simulation of the gravity drainage processes using COMSOL® software. An appropriate design of experiment (DOE) method was selected to find the most important parameters contributing in gravity drainage and then conduct the experiments in a useful as well as economic manner. A two-dimensional experimental setup was employed to investigate free fall gravity drainage (FFGD) and controlled gravity drainage (CGD) using unconsolidated glass beads fractured porous media having various fractures configurations. Flow visualization measurements were carried out. Following the flow visualization experiments, parametric sensitivity analysis was performed considering the effects of different system parameters such as fracture aperture, matrix height, permeability, and fluid properties on the dependent variables including drainage rate, critical pumping rate, maximum drainage rate, recovery factor and so on. These experiments enabled us to capture some aspects of the recovery mechanism and the flow communication between matrix block and fracture during gravity drainage. After analyzing the experimental data for the FFGD test runs, it was found that the rate of liquid flowing from matrix to fracture is proportional to the difference of liquid levels in the matrix and in the fracture. In addition, the characteristic rate and the maximum liquid drainage rate from the fractured models were determined for such a stable gravity-dominated process. The experiments showed that the presence of fracture is more influential in lower matrix permeability systems. For a given fracture-matrix system with different initial liquid saturation conditions, it was seen that the production history can be correlated by plotting the fraction of recoverable liquid as a function of time. Furthermore, the recovery factor can be correlated using dimensionless numbers such as the Bond number and the dimensionless time. For the controlled gravity drainage (CGD) test runs conducted, the experimental results indicated that higher pumping rates cause a higher difference between the liquid levels in the fracture and in the matrix, thus the gas breakthrough happens sooner. Moreover, it was found that as long as the porous medium is drained with a constant liquid pumping rate but lower than critical rate, the height difference between the G-L interfaces in matrix and fracture remains constant. In this study, a new concept of “Critical Pumping Rate” (CPR) was defined at which each particular porous medium has recovery factor equal to the recovery factor for higher rates just before gas breakthrough. The difference between liquid levels in fracture and matrix remains unchanged at rates higher than CPR. Known this particular withdrawal rate, there are two main advantages, namely: 1) choosing a pumping rate lower than it to drain the reservoir without getting gas breakthrough; and 2) understanding the physics of pumping behaviour from fractured media and extending the concept to the real cases. In addition, the maximum liquid pumping rate from each physical model was studied and it was found that the rate depends strongly on the storage capacity of the fractures, petrophysical properties of each model as well as physical properties of test fluids. The critical rate, maximum rate, recovery factor at gas breakthrough and difference of gas liquid interface positions in matrix and fracture were correlated by dimensionless numbers such as Bond number, Capillary, and the ratio of permeabilities. Linear regression correlations presented in this study can predict production history and flow behaviour in the fractured porous media for a wide range of dimensionless numbers. The COMSOL® software was used to numerically simulate the gravity drainage processes in the two-dimensional flow experiments for fractured porous media. The parameters of the model were based on theory, as well as on the results of the two-dimensional gravity drainage experiments. The simulation results for the gravity drainage processes compared favourably with the experimental results, as a good match between the numerical solution and the experimental data was found. The simulation model developed provides a basis for further modeling of gravity drainage process in more complicated porous media.

In fractured geologic media, flow and contaminant transport are predominantly controlled by the fractures, their distribution and connectivity. The accurate characterization of fractured geologic medium, imaging of fracture patterns and their connectivity have been a challenge for decades. Given the complexities of fractured networks in the subsurface and Dense Non Aqueous Phase Liquid (DNAPL) contamination, in this thesis, transient hydraulic tomography (THT), a recently developed tool for characterizing aquifer heterogeneity is evaluated under laboratory conditions to delineate discrete fractures. Laboratory experiments and modeling studies are also conducted to understand TCE plume behavior. A dolomite rock sample, which is 91.5 cm in length, 60.5 cm in height and 5 cm thick, was fractured in the laboratory to perform the experiments. After the fractured block was enclosed in a flow cell, flow-through and pumping tests were conducted to characterize the fractured rock block. The data from the pumping tests were then analyzed using the SSLE code developed by Zhu and Yeh [2005] and transient hydraulic tomography (THT) was conducted to image the fracture pattern and their connectivity through the delineation of K and Ss distributions (the tomograms). Synthetic pumping tests, identical in configuration to the laboratory ones were also conducted using HydroGeoSphere (HGS) [Therrien et al, 2009] in a synthetic replica of the fractured block to compare the observed and simulated drawdowns. Then synthetic THT analysis was performed utilizing the synthetic pumping test data to compare the tomograms obtained from the THT analysis of synthetic and laboratory pumping tests. Results suggest that the THT analysis of multiple laboratory pumping tests captured the fracture pattern and their connectivity quite well and they became more vivid with the additional pumping tests. The estimated high hydraulic conductivity (K) and low specific storage (Ss) zones clearly show the fractures and their connectivity. The pattern of K and Ss tomograms obtained from the analyses of synthetic and laboratory pumping tests were similar. Estimated K and Ss values for the fractures and the matrix may not exactly replicate the actual K and Ss values for the fractured rock, but the model also provides uncertainty estimates associated with the resulting K and Ss tomograms. In this study, two cases of transient hydraulic tomography (THT) analysis of the laboratory pumping tests were performed by changing the location of 2nd and 3rd pumping tests among the three to examine if there is any significant impact of these pumped location on the pattern of resulting hydraulic conductivity (K) and specific storage (Ss). The initial pumping test was the same for two cases. Results show that the patterns of estimated K and Ss tomograms obtained from these two cases are similar, although the pumped locations (2nd and 3rd tests among the three) utilized for the inversion were different for two cases suggesting that the location of these later pumping tests does not significantly impact the estimates for this fractured rock block. However, the initial test should be selected carefully as that seems to set the pattern of the tomograms. The estimated K and Ss tomograms were validated by predicting five independent pumping tests conducted in the fractured rock block. These five pumping tests were not included during the construction of the K and Ss tomograms. For most of the independent pumping tests, good correspondence between the simulated and observed drawdown was achieved. The study indicates that, it is possible to delineate discrete fractures, their pattern and connectivity by carefully applying of THT analysis of multiple pumping tests based on the inverse code SSLE [Zhu and Yeh, 2005]. In addition, hydraulic tomography seems to be a cost effective tool for characterizing fractured rock since it does not require the detailed information on fracture geometry parameters such as aperture, trace length, orientation, spatial distribution, and connectivity, which are difficult to quantify. These parameters are usually unavailable between boreholes. Therefore, THT appears to be a promising approach in delineating fractures and their connectivity in subsurface. However, it is still at the early stage as the study was conducted in the laboratory under controlled conditions. Small scale field experiments need to be conducted to validate THT as a tool for the characterization of hydraulic parameters of fractured rocks. Upon completion of the hydraulic characterization, several conservative tracer tests were conducted using bromide (Br-) as a conservative tracer to aid in the design of TCE dissolution experiment. Once the tracer experiments were completed, a known volume of pure phase TCE was injected at a known location in the flow cell to create a well-defined source zone. A constant hydraulic gradient was maintained by fixing the hydraulic heads at the two head tanks to induce steady groundwater flow through the flow cell. Water samples were obtained at a down gradient monitoring port for 3 months to obtain a long-term breakthrough curve of TCE in the aqueous phase. The purpose of this experiment was to study TCE dissolution behaviour in the fractured rock sample. Then HydroGeoSphere (HGS) was used to model the aqueous phase TCE transport using two separate approaches: 1) the Discrete Fracture Network modeling approach and 2) the stochastic continuum approach, to investigate whether they can capture the dissolution behavior. Both approaches were able to capture the pattern of the breakthrough curve in the fractured rock. The discrete fracture approach captured the observed TCE plume and the dissolution behavior quite well. On the other hand, the stochastic continuum approach, in which the fractured rock block was treated as porous medium having a heterogeneous K field obtained from THT analysis, also appeared to be promising in capturing the aqueous phase transport of TCE. Despite some early time deviation, the simulated breakthrough curve captured the overall observed concentration profile. However, the stochastic continuum approach seems to be more cost effective as it does not require detailed information about fracture aperture and their spatial distribution which are difficult if not impossible to obtain between boreholes. Note that, the studies were conducted based on a laboratory experiment conducted in a controlled environment. The experimental block was well characterized and the geometry of the experimental block as well as the flow through the system was well understood from the hydraulic and tracer experiments. Thus small scale field experiment is required to support this conclusion.

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.

The oil production from well fractured carbonate reservoirs is a considerable part of the total oil production in the world. The petroleum resource base in naturally fractured reservoirs is estimated to be in the range of billions of barrels in the U.S and in addition, a multibillion- barrel international oil resource base exists in naturally fractured reservoirs. Gravity drainage is important in some of oil recovery processes, either acting as the driving force in processes using horizontal wells or altering the displacement patterns during water-flooding, chemical flooding, CO2 flooding and other EOR methods. The gravity drainage process has a major effect on oil recovery from oil reservoirs. Gravity drainage driven oil production in naturally fractured and other complex reservoirs falls into two regimes: the balk flow regime and the film flow regime. Oil recovery by gravity drainage in a fractured reservoir strongly depends on the capillary height of the porous medium. Capillarity and gravity forces are usually the major driving forces in fractured reservoirs. This PhD thesis consists of two main parts namely: 1) Experimental works on gravity drainage, and 2) Modeling and simulation of the gravity drainage processes using COMSOL® software. An appropriate design of experiment (DOE) method was selected to find the most important parameters contributing in gravity drainage and then conduct the experiments in a useful as well as economic manner. A two-dimensional experimental setup was employed to investigate free fall gravity drainage (FFGD) and controlled gravity drainage (CGD) using unconsolidated glass beads fractured porous media having various fractures configurations. Flow visualization measurements were carried out. Following the flow visualization experiments, parametric sensitivity analysis was performed considering the effects of different system parameters such as fracture aperture, matrix height, permeability, and fluid properties on the dependent variables including drainage rate, critical pumping rate, maximum drainage rate, recovery factor and so on. These experiments enabled us to capture some aspects of the recovery mechanism and the flow communication between matrix block and fracture during gravity drainage. After analyzing the experimental data for the FFGD test runs, it was found that the rate of liquid flowing from matrix to fracture is proportional to the difference of liquid levels in the matrix and in the fracture. In addition, the characteristic rate and the maximum liquid drainage rate from the fractured models were determined for such a stable gravity-dominated process. The experiments showed that the presence of fracture is more influential in lower matrix permeability systems. For a given fracture-matrix system with different initial liquid saturation conditions, it was seen that the production history can be correlated by plotting the fraction of recoverable liquid as a function of time. Furthermore, the recovery factor can be correlated using dimensionless numbers such as the Bond number and the dimensionless time. For the controlled gravity drainage (CGD) test runs conducted, the experimental results indicated that higher pumping rates cause a higher difference between the liquid levels in the fracture and in the matrix, thus the gas breakthrough happens sooner. Moreover, it was found that as long as the porous medium is drained with a constant liquid pumping rate but lower than critical rate, the height difference between the G-L interfaces in matrix and fracture remains constant. In this study, a new concept of “Critical Pumping Rate” (CPR) was defined at which each particular porous medium has recovery factor equal to the recovery factor for higher rates just before gas breakthrough. The difference between liquid levels in fracture and matrix remains unchanged at rates higher than CPR. Known this particular withdrawal rate, there are two main advantages, namely: 1) choosing a pumping rate lower than it to drain the reservoir without getting gas breakthrough; and 2) understanding the physics of pumping behaviour from fractured media and extending the concept to the real cases. In addition, the maximum liquid pumping rate from each physical model was studied and it was found that the rate depends strongly on the storage capacity of the fractures, petrophysical properties of each model as well as physical properties of test fluids. The critical rate, maximum rate, recovery factor at gas breakthrough and difference of gas liquid interface positions in matrix and fracture were correlated by dimensionless numbers such as Bond number, Capillary, and the ratio of permeabilities. Linear regression correlations presented in this study can predict production history and flow behaviour in the fractured porous media for a wide range of dimensionless numbers. The COMSOL® software was used to numerically simulate the gravity drainage processes in the two-dimensional flow experiments for fractured porous media. The parameters of the model were based on theory, as well as on the results of the two-dimensional gravity drainage experiments. The simulation results for the gravity drainage processes compared favourably with the experimental results, as a good match between the numerical solution and the experimental data was found. The simulation model developed provides a basis for further modeling of gravity drainage process in more complicated porous media.

Both onsite septic disposal systems and private drinking water wells are commonly utilized in rural areas of Canada. The coexistence of septic systems and drinking water wells has the potential to greatly impact the quality of water obtained in these settings. Human enteric viruses have been recognized as a potential source of groundwater borne disease, although the level of risk they pose and the processes responsible for their transport are poorly understood. As a result of thin overburden, low storage capacity, and high groundwater velocities, fractured rock aquifers are potentially at highest risk to viral contamination. However, only limited research has been conducted to explore this concern. The current study was conducted to investigate both the rate of occurrence of human viruses in fractured rock aquifers and the transport mechanisms acting in these settings. A survey was conducted to identify the prevalence of human enteric viruses in three fractured rock aquifers located across Canada. A total of 61 samples were collected from 28 wells drilled in aquifers in Ontario, Newfoundland, and British Columbia. Molecular PCR techniques were utilized to determine virus presence. Results showed that 37.7% of samples and 58.1% of wells were at some time positive for viruses. Virus presence was found to increase with housing density and viruses were found to travel distances of at least 40 meters. Poor correlation was found between the presence of viruses and traditional bacterial indicators. A field-scale viral infiltration experiment was conducted to investigate viral transport behavior. The bacteriophage ф-X174 and the fluorescent dye Lissamine FF were utilized as viral and solute tracers, respectively. Tracers were applied to an exposed rock outcrop exhibiting fractures with known connection to two nearby wells. Breakthrough was extremely rapid and the colloidal processes of decreased dispersion and slow-release kinetic sorption were identified. This study has provided concrete evidence that viral contamination poses a significant threat to fractured groundwater aquifers in rural areas where onsite septic disposal practices are utilized. The results observed in this study suggest that current set back distances and monitoring techniques may be inadequate to prevent exposure to human viruses. / Thesis (Master, Civil Engineering) -- Queen's University, 2010-11-09 23:07:31.595

The aim of this study was to assess the performance of thermal heating for the removal of chlorinated solvents from fractured rock. The study included a laboratory experimental program, a field pilot study demonstration and a mathematical modeling component. In the laboratory experimental program, thermal heating parameters, such as operational temperature, heating duration, and the corresponding degree of contaminant removal, were evaluated through a series of heating tests. To evaluate the effect of heating temperature and heating duration on the degree of contaminant mass removal, two different heating profiles were utilized during the experiments. Additionally, seven types of rock and two common contaminants were selected to evaluate the effect of thermal heating on different geological media impacted with different chlorinated compounds. In general, results showed that heating duration had the most significant effect on the degree of contaminant mass removal in post-remedy samples. Results showed that a higher porosity in combination with a lower organic content facilitates the removal of chlorinated solvents from the rock matrix. A Thermal Conductive Heating (TCH) pilot test was implemented by TerraTherm, Inc. at the former Naval Air Warfare Center (NAWC) in West Trenton, NJ to assess the performance of TCH for the removal of trichloroethylene (TCE) and daughter products (i.e cis-1,2-dichloroethylene (DCE) and vinyl chloride (VC)) from fractured bedrock. Results showed that treatment removed 318.5 kg of TCE, DCE and VC, from the treatment zone, of which 62.6 kg were recovered from the rock matrix. A total of 63 % TCE, 65.8 % of DCE and 90.4% of VC were removed during heating. Finally, Semi-analytical solutions were derived to evaluate back diffusion in a fractured bedrock environment where the initial condition comprises a spatially uniform, non-zero matrix concentration throughout the domain. It was concluded that the time required to reach a desired fracture pore water concentration is a function of the distance between the point of compliance and the upgradient face of the domain where clean groundwater is inflowing. Hence, shorter distances correspond to reduced times required to reach compliance. / Thesis (Ph.D, Civil Engineering) -- Queen's University, 2012-09-24 11:30:16.52

No description available.