<p>Owing to the lack of knowledge pertaining to the fate and transport of microorganisms in fractured aquifers, the research presented in this thesis was designed to improve the mechanistic understanding of particulate transport in fractures by conducting tracer experiments in natural and epoxy replica fractures. This research demonstrated that particulate retention within fractures is heavily dependent on the equivalent mass balance aperture, followed by the coefficient of variation of the aperture field, and then by the flow conditions. It was also shown that the fracture aperture field alone, not the flow rate or the matrix properties, determines the number of fracture volume flushes required to achieve a 2-3 log decrease in effluent concentration. Moreover, a statistical model was developed that identifies the most important factors affecting particulate retention as the ratio of the ionic strength of solution to the charge of the collector, the ratio of the particle to collector charges, and the Peclet number. The model is able to reasonably predict particulate retention. Finally, tracer experiments conducted in a natural fracture and an epoxy replica of that fracture isolated the effects of matrix properties on attachment, and hence, retention. The transparent nature of the replica fracture was exploited to capture images of <em>E. coli</em> RS2-GFP transport. These images reveal preferential transport within the fracture, and also show that the preferential pathway broadens slightly under increasing flow conditions. This broadening is likely due to higher fluid pressures associated with larger specific discharges. In the groundwater field, there is so little fracture-specific information available that coupling the understanding of a critical environmental setting (fractures) with high-quality particulate tracer experiments and associated modeling represents a significant contribution to the body of science.</p> / Doctor of Philosophy (PhD)

<p>A series of experiments were carried out to gain a better understanding of the mechanisms governing the transport of biological and non-biological particles through single, saturated dolomite fractures at the laboratory scale. Fracture apertures and general roughness were characterized using hydraulic and conservative solute tracer experiments.</p> <p>The effects of particle size, surface characteristics and ionic strength of carrying solution were all evaluated. Particulate material studied included MS2, <em>E. coli</em> and two sizes of carboxylated microspheres. To elucidate the effect of ionic strength on particulate transport, the ionic strength of the carrying solution was altered during each experiment. All particulate experiments were completed at a specific discharge of 15 m/day to facilitate comparisons.</p> <p>Recovery of biological particulate material was found to be much less relative to the carboxylated microspheres, even though the energy profiles predicted similar interactions with the fracture surface. This suggests that the biological surface has a significant impact on retention within the fracture. Further, altering the ionic strength of the carrying solution did not spur significant elution of additional particulate material, regardless of surface characteristics. Therefore, it was determined that retention within the secondary energy minimum was negligible under these operating conditions.</p> <p>With respect to carboxylated microspheres, increased retention was observed within the less variable fracture. This suggests that increased variability within a fracture results in increased eddying within the aperture field. This eddying effectively reduces the aperture region available for particle transport, lessening the particle/fracture interaction. Overall, while mean residence times were similar, recovery of biological particles was poorly replicated by microspheres.</p> / Master of Applied Science (MASc)

Pump and treat remediation is the most commonly used method to remediate contaminated aquifers, but the effectiveness decreases when heterogeneities are introduced. Fractures within the matrix cause large differences in hydraulic conductivity. The low hydraulic conductivity of the matrix acts as an area of storage for contaminant, allowing for attenuation of the plume. The attenuation of the plume causes the effectiveness of the system to decrease and cost of remediation to increase. In order to understand what parameters enhance contaminant storage in the matrix, rapid transport in fractures, and both of their influences on the efficiency of the pumping system, a hypothetical model was developed to simulate the release and remediation of a plume using pumping. The code used was HydroGeoSphere, which allowed for the interpretation of parameters influencing contaminant storage during the withdrawal phase of the pump and treat remediation by allowing transport of contaminant within both the matrix and the fractures. Matrix parameters of porosity and hydraulic conductivity influenced the effectiveness of the withdrawal system most. For instance, the difference in percent mass extracted between porosity values of 0.01 and 0.4 was 23.75%, while the difference between fracture lengths of 200 and 400 m was 5.59%. Fracture pattern influenced where the stored contaminant was located within the matrix. Downgradient of the source, six different fracture patterns resulted in a difference in relative concentration of 0.4 at the start of the withdrawal phase. Evaluation of remediation included both percent extraction of contaminant and finer scale remediation of the contaminant specifically within the matrix. Multiple length-scale observations helped determine how fracture and matrix parameters influence remediation in dual permeability media. / Geology

Hydraulic fracturing is a controversial practice because of concerns about environmental impacts due to its widespread use in recovering unconventional petroleum and natural gas deposits. However, water-only hydraulic fracturing has been used safely and successfully for many years to increase the permeability of aquifers used for drinking and irrigation water supply. This process extends and widens existing bedrock fractures, allowing groundwater storage to increase. Researchers have studied the behavior of fractured-rock aquifers for decades, but little has been published on the hydraulic and mechanical properties of hydraulically enhanced fractures. In this study, a multi-faceted approach consisting of aquifer and tracer testing is used to estimate the transmissivity and storage coefficient of a hydraulically induced fracture and observe its behavior as a contaminant flow pathway. The results of the aquifer tests indicated a decrease in both the transmissivity and storage coefficient of the fracture of three orders of magnitude after cessation of pumping. The aquifer temporarily experienced incomplete recovery following pumping tests, likely due to slow recharge. After complete recovery occurred, subsequent tests showed that these hydraulic properties returned to their original values, indicating elastic compression of the fracture during periods of applied stress. The results of the tracer test indicated rapid, uniform, one-dimensional flow through the fracture, with average fluid velocity approaching 1 km/day in an induced steady flow field of 6 x 10-5 m3/s (1 gal/min) and a fracture volume of 0.238 m3 (63 gal). The complex heterogeneity of fractured-rock aquifers necessitates the use of multiple lines of testing in order to arrive at a detailed description of the behavior of these systems. This study demonstrates one effective method of investigating a single fracture that can uncover information about the behavior of a hydraulically enhanced aquifer that is otherwise difficult to obtain. / Master of Science

No description available.

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

A large percentage of oil and gas reservoirs in the most productive regions such as the Middle East, South America, and Southeast Asia are naturally fractured reservoirs (NFR). The major difference between conventional reservoirs and naturally fractured reservoirs is the discontinuity in media in fractured reservoir due to tectonic activities. These discontinuities cause remarkable difficulties in describing the petrophysical structures and the flow of fluids in the fractured reservoirs. Predicting fluid flow behavior in naturally fractured reservoirs is a challenging area in petroleum engineering. Two classes of models used to describe flow and transport phenomena in fracture reservoirs are discrete and continuum (i.e. dual porosity) models. The discrete model is appealing from a modeling point of view, but the huge computational demand and burden of porting the fractures into the computational grid are its shortcomings. The affect of natural fractures on the permeability anisotropy can be determined by considering distribution and orientation of fractures. Representative fracture permeability, which is a crucial step in the reservoir simulation study, must be calculated based on fracture characteristics. The diagonal representation of permeability, which is customarily used in a dual porosity model, is valid only for the cases where fractures are parallel to one of the principal axes. This assumption cannot adequately describe flow characteristics where there is variation in fracture spacing, length, and orientation. To overcome this shortcoming, the principle of the full permeability tensor in the discrete fracture network can be incorporated into the dual porosity model. Hence, the dual porosity model can retain the real fracture system characteristics. This study was designed to develop a novel approach to integrate dual porosity model and full permeability tensor representation in fractures. A fully implicit, parallel, compositional chemical dual porosity simulator for modeling naturally fractured reservoirs has been developed. The model is capable of simulating large-scale chemical flooding processes. Accurate representation of the fluid exchange between the matrix and fracture and precise representation of the fracture system as an equivalent porous media are the key parameters in utilizing of dual porosity models. The matrix blocks are discretized into both rectangular rings and vertical layers to offer a better resolution of transient flow. The developed model was successfully verified against a chemical flooding simulator called UTCHEM. Results show excellent agreements for a variety of flooding processes. The developed dual porosity model has further been improved by implementing a full permeability tensor representation of fractures. The full permeability feature in the fracture system of a dual porosity model adequately captures the system directionality and heterogeneity. At the same time, the powerful dual porosity concept is inherited. The implementation has been verified by studying water and chemical flooding in cylindrical and spherical reservoirs. It has also been verified against ECLIPSE and FracMan commercial simulators. This study leads to a conclusion that the full permeability tensor representation is essential to accurately simulate fluid flow in heterogeneous and anisotropic fracture systems. / text

Shale oil (tight oil) is oil trapped in low permeability shale or sandstone. Shale oil is a resource with great potential as it is heavily supplementing oil production in the United States (U.S. Energy Information Administration, 2013). The shale rock must be stimulated using hydraulic fracturing before the production of shale oil. When the hydrocarbons are produced from fractured systems, the resulting flow is influenced by the fracture, the stimulated rock, and the matrix rock. The production decline rates from shale oil reservoirs experience flow regimes starting with fracture linear flow (fracture dominated), then bilinear flow (fracture and stimulated rock dominated), then formation linear flow (stimulated rock dominated), and finally pseudo-radial flow (unstimulated matrix rock dominated) (Cinco-Ley 1982). In this thesis, daily production rates from a shale oil reservoir are modeled using a simple spreadsheet-based, finite difference serial flow simulator that models the single-phase flow of a slightly-compressible oil. This simulator is equivalent to flow through multiple tanks (subsequent part of the thesis will call these cells) through which flow passes serially through one tank into the other. The simulator consists of 11 tanks. The user must specify the compressibility-pore volume product of each tank and the transmissibility that governs flow from one tank to another. The calculated rate was fitted to the given data using the Solver function in Excel. The fitted matches were excellent. Although we can adjust all 22 parameters (2 per cell) to affect the simulation results, we found that adjusting only the first three cells nearest to the well was sufficient. In many cases, only two cells were enough. Adjusting 4 or more cells resulted in non-unique matches. Furthermore, the properties of the very first cells proved insensitive to the matches when using the 3 cells to match the data. The cells in the 2 cell model represent the stimulated zone and the unstimulated rock. Likewise, the cells in the 3 cell model represent the hydraulic fracture, the stimulated zone, and the unstimulated rock. The accessed pore volume and transmissibility were responsive to the injected sand mass and fluid volume up to approximately 10⁶ kg and 7000 m³ respectively; injecting more sand and fluids than this caused negligible increases in the accessed pore volume and transmissibility. This observation suggests that the sand does not migrate far into the fractures. Similarly, it was observed that the number of stages was positively correlated with cell transmissibility and pore volume up to 20 stages. These results suggest that fracture treatments were significantly over designed and injecting less sand and water in fewer stages would optimize the economics of similar projects. To our knowledge this is the first work to analyze the results of fracture treatments by matching with pore volumes and transmissibility in a simple serial cell flow. / text

Pour étudier le comportement mécanique des massifs rocheux, en prenant en compte le réseau des discontinuités au sein de la roche intacte, cette recherche a pour objectif la représentation du massif par des modèles géométriques basés sur des relevés de terrain et l'analyse de ces modèles par l'utilisation d'outils informatiques adaptés pour les milieux granulaires. Le premier chapitre fait l'état de l'art des roches fracturées, des méthodes numériques de la mécanique des roches et des approches du calcul de structure d'un tunnel. Ces études conduisent à la proposition d'une méthodologie depuis les recherches in situ jusqu'à la modélisation et l'analyse mécanique, présentée dans le deuxième chapitre. Le massif rocheux est d'abord représenté géométriquement par la distribution de ses discontinuités, et l'utilisation du logiciel RESOBLOK basé sur la méthode du Réseau de Fractures Discrètes. Les modèles mécaniques de massifs rocheux sont ensuite présentés à partir des données sur les études de l'histoire du massif, et des mesures faites sur site et en laboratoire. Les modèles numériques en 3D sont analysés par l'utilisation du logiciel LMGC90 basé sur la méthode de la Dynamique des Contacts Non Réguliers. Les premières applications de la méthodologie sont exposées : la création d'une roche numérique pour simuler un essai de compression triaxiale, et la simulation d'une excavation multi phases d'un tunnel au rocher. La méthodologie proposée a été appliquée sur le marbre blanc de Saint Béat (Haute Garonne, France) et les résultats préliminaires sont donnés dans le chapitre trois. Les réponses mécaniques de la roche numérique sont analysées et son comportement est caractérisé. / Aimed at studying the mechanical behaviour of rock mass and considering the presence of the discontinuity network in the intact rock, this research concentrates on how the rock can be represented in suitable geometrical models, on the basis of site measurements, and then appropriately analysed using computer tools developed for the study of granular media. The first chapter deals with a bibliographical study on fractured rock and tunnel engineering. Different computational methods of rock mechanics are introduced. Simultaneously, three principal approaches for tunnel structural design are recalled. These studies lead to the proposition of a methodology from the in situ investigation to in-door modelling and mechanical analysis, presented in the second chapters. The rock mass is first geometrically represented through the distribution of discontinuities in the rock mass and the use of the RESOBLOK code based on the Discrete Fracture Network method. Mechanical models of rock mass are then presented from the data of historical studies on the rock mass and from laboratory and in situ measurements. The 3D computational models are analysed using the LMGC90 based on the Non Smooth Contact Dynamics method. The first two applications of the methodology are introduced: the generation of the numerical rock for the simulation of the triaxial compression test, and the simulation of multi-phase excavation of rock tunnel. The proposed methodology has been applied on the white marble of Saint Béat (Haute Garonne, France) and the initial results are given in the third chapter. The mechanical responses of the numerical rock mass are analysed and the bulk behaviour of the rock is evaluated.

O estudo da hidrogeologia em aqüífero cristalino fraturado apresenta como maior desafio a imprevisibilidade do meio, uma vez que o fluxo ocorre exclusivamente pelas fraturas. Neste sentido informações de atitude, abertura e conectividade de fraturas tornam-se essenciais no entendimento do fluxo da água subterrânea. O presente trabalho busca aplicar o conhecimento da tectônica rúptil para caracterizar o meio fraturado, com intuito de auxiliar o estudo hidrogeológico do aqüífero. A área de estudo está localizada na porção leste do Município de Cotia, SP, no contexto do Complexo Embu, de idade proterozóica, na porção central da Faixa Ribeira. Ocorrem predominantemente rochas granitóides, com permeabilidade primária desprezível, rochas cataclásticas relacionadas às grandes zonas de cisalhamento de Taxaquara e Caucaia do Alto, metassedimentos e coberturas aluviais cenozóicas, estas ao longo das principais drenagens. Os principais eventos tectônicos de caráter rúptil da região estiveram relacionados à reativação mesozóico-cenozóica, inicialmente com a ruptura continental e abertura do Oceano Atlântico Sul e, posteriormente, a formação do Rift Continental do Sudeste do Brasil (RCSB), e seus principais eventos de deformação. Os estudos foram desenvolvidos em duas escalas de trabalho. A primeira de semi-detalhe, em 1:50.000, e a segunda de detalhe, em 1:5.000, esta última em duas porções distintas dentro da área de semi-detalhe. Os resultados obtidos em escala de semi-detalhe, a partir da análise da tectônica rúptil (morfometria, fotogeologia e levantamentos de campo), serviram para direcionar as investigações em escala de detalhe, onde foram usados dados provenientes de sondagens, métodos BHTV (Bore Hole Television), levantamento de perfis de eletrorresistividade e dados hidrogeológicos de poços de monitoramento instalados no aqüífero em meio fraturado. Os estudos de escala de semi-detalhe permitiram delinear lineamentos de direções N-S, E-W, NE-SW e NW-SE, originados pelas diversas fases de deformação rúptil. Esses dados de lineamentos foram correlacionados com dados de atitude de falhas e juntas medidas em afloramento, e, posteriormente, associados às diferentes fases de deformação descritas para o RCSB. Foram definidas seis famílias de juntas: a) família NNW-SSE, com mergulho alto para NE); b) família NE-SW com mergulho sub-vertical; c) família E-W com mergulho alto para NW); d) família WNW-ESE com mergulho alto para NE); e) família NESW com mergulho médio para NW; e f) família NE-SW com mergulho médio para SE. Os resultados dos estudos desenvolvidos em escala de semi-detalhe, somados aos condicionantes geológicos e geomorfológicos das áreas selecionadas para estudos de detalhe, orientaram a instalação de poços de monitoramento verticais e inclinados nessas últimas. Foram instalados 11 poços na área de detalhe A e quatro poços na área de detalhe B. Os poços foram perfilados por televisionamento acústico e as fraturas identificadas nesses poços foram separadas por famílias e subfamílias: Família 1 (sub-horizontal), Família 2A (N-S com mergulho médio para W) e 2B (N-S com mergulho alto para E), Família 3A (NE-SW com mergulho médio para NW) e 3B (NE-SW com mergulho médio para SE) e Família 4 (WNW-ESE sub-vertical). Na maioria dos poços foram realizados ensaios hidráulicos para determinação da condutividade hidráulica (K), que variou de 1,55x10-8 a 2x10-5 cm/s, com média de 2,98x10-6 e desvio padrão de 5,60 x10-6 cm/s. Ensaios de bombeamento indicaram fluxo de água subterrânea através das estruturas N-S, com boa conexão com as estruturas NE-SW. A abertura média das fraturas foi de 1,13x10-3 cm, calculada com base nos valores de condutividade hidráulica correlacionados aos dados estruturais. A partir dos dados hidrogeológicos e estruturais foi proposto um método para se estimar a transmissividade e condutividade hidráulica de um poço instalado no aqüífero fraturado. / The main challenge in a hydrogeological study of a fractured-bedrock aquifer is the unpredictability of media, because the ground-water flows only through the fractures. Informations about strike, dip, opening and connectivity among fractures become essential in order to understand the ground-water flow. This work is an application of the concepts of the brittle tectonics to characterize a fractured media in order to help the hydrogeological study of a fractured-bedrock aquifer. The study area is located at the eastern part of the City of Cotia, State of São Paulo, Southeastern Brazil, in the context of the Proterozoic Embu Complex, in the central part of Ribeira Belt. This area encompass predominantly granitoid rocks (with negligible primary permeability), cataclastic rocks related to the Taxaquara and Caucaia do Alto shear zones, and metasedimentary rocks. Quaternary alluvial deposits occur along the main drainages. The main brittle tectonic events in the region are related to the Mesozoic-Cenozoic reactivation, associated with the continental break-up, the opening of the South Atlantic Ocean, and the formation and deformation of the basins of the Continental Rift of the Southeastern Brazil (CRSB). The study was carried out in two working scales. First in a semi-detail scale, at 1:50,000, of a large area, were brittle tectonic analysis (morphometry, photogeology and field survey in outcrops) was performed. The achieved results were used to select two small areas, within the semi-detail area, for a detailed investigation, on scale of 1:5,000, in which data from rock-boring, BHTV (Bore Hole Television), eletroresistivity and hydrogeology from bedrock monitoring wells were obtained. At a semi-detail scale it was possible to delineate lineaments of N-S, E-W, NE-SW and NW-SE-trending directions, probably originated during several brittle deformational events. Data from lineaments were first correlated with data of fault and joint measured in outcrops and latterly associated with different deformation phases described in CRSB. Six joint families were defined: a) steeply northeastward dipping NNW-SSE-trending family; b) sub-vertically dipping NE-SW trending family; c) steeply northwestward dipping E-W-trending family; d) steeply northeastward dipping WNW-ESE-trending family); e) NE-SW-trending family with moderate dipping towards the NW; and f) NE-SW-trending family with moderate dipping towards the SE. Results of studies in the semi-detail investigation area allowed the selection of sites for drilling of 11 wells in detail area A and 4 wells in the detail area B. All the wells were surveyed by acoustic teleview and the identified fractures were classified in four families and two subfamilies: family 1 (sub-horizontal to gently dipping); family 2A (N-S-trending with moderate dipping towards the W) and 2B (NS-trending with steep dipping towards the E); family 3A (NE-SW-trending with moderate dipping towards the NW) and 3B (NE-SW-trending with moderate dipping towards the SE); and family 4 (WNW-ESE-trending with sub-vertical dipping). Hydraulic conductivity (K) measured in the wells ranged from 1.55x10-8 to 2x10-5 cm/s, with a mean of 2.98x10-6 and standard deviation of 5.60 x10-6 cm/s. Pumping tests showed the groundwater flow through N-S-trending structures and good connection with NE-SW-trending structures. The average opening of fractures attained 1.13x10-3 cm and was calculated on the basis of the amount of hydraulic conductivity correlated with structural data. Based on structural and hydrogeological data it was proposed a method for estimating the transmissivity and hydraulic conductivity of a well installed in the fractured-bedrock aquifer.