The Unconventional Resource Plays consisting of the lowest tier of resources (large volumes and most difficult to develop) have been the main focus of US domestic activity during recent times. Horizontal well drilling and hydraulic fracturing completion technology have been primarily responsible for this paradigm shift. The concept of drainage volume is being examined using pressure diffusion along streamlines. We use diffusive time of flight to optimize the number of hydraulic fracture stages in horizontal well application for Tight Gas reservoirs. Numerous field case histories are available in literature for optimizing number of hydraulic fracture stages, although the conclusions are case specific. In contrast, a general method is being presented that can be used to augment field experiments necessary to optimize the number of hydraulic fracture stages. The optimization results for the tight gas example are in line with the results from economic analysis. The fluid flow simulation for Naturally Fractured Reservoirs (NFR) is performed by Dual-Permeability or Dual-Porosity formulations. Microseismic data from Barnett Shale well is used to characterize the hydraulic fracture geometry. Sensitivity analysis, uncertainty assessment, manual & computer assisted history matching are integrated to develop a comprehensive workflow for building reliable reservoir simulation models. We demonstrate that incorporating proper physics of flow is the first step in building reliable reservoir simulation models. Lack of proper physics often leads to unreasonable reservoir parameter estimates. The workflow demonstrates reduced non-uniqueness for the inverse history matching problem. The behavior of near-critical fluids in Liquid Rich Shale plays defies the production behavior observed in conventional reservoir systems. In conventional reservoirs an increased gas-oil ratio is observed as flowing bottom-hole pressure is less than the saturation pressure. The production behavior is examined by building a compositional simulation model on an Eagle Ford well. Extremely high pressure drop along the multiple transverse hydraulic fractures and high critical gas saturation are responsible for this production behavior. Integrating pore-scale flow modeling (such as Lattice Boltzmann) to the field-scale reservoir simulation may enable quantifying the effects of high capillary pressure and phase behavior alteration due to confinement in the nano-pore system.

Understanding economic success of unconventional production from shales requires an explanation of the relationship between induced seismicity and hydraulic fracturing. This thesis deals with observing and analyzing synthetic and real microseismic monitoring data acquired during hydraulic fracturing. The thesis is based on observation and analyses of source mechanisms of induced microseismic events that have recently become regularly inverted and interpreted in the oil and gas industry. The results of analyses are interpreted with the geomechanical model of the relationship between hydraulic fracturing and induced seismicity. The study of source mechanisms starts with detailed analyses of spatial distribution of full moment tensor inversion stability. It was mapped based on synthetically computed condition numbers in the vicinity of different monitoring arrays including dense arrays at the surface and sparse arrays with sensors in the boreholes. Stability of inversion was tested under several conditions, mainly dependency on size and geometry of monitoring array and level of noise in the data. In this part of the thesis it is shown that dense surface arrays may provide very stable inversion of source mechanisms which may be interpreted. The study shows that an increasing percentage of non-shear...

By projecting observed microseismic data backward in time to when fracturing occurred, it is possible to locate the fracture events in space, assuming a correct velocity model. In order to achieve this task in near real-time, a robust computational system to handle backward propagation, or Reverse Time Migration (RTM), is required. We can then test many different velocity models for each run of the RTM. We investigate the use of a Graphics Processing Unit (GPU) based system using Compute Unified Device Architecture for C (CUDA-C) as the programming language. Our preliminary results show a large improvement in run-time over conventional programming methods based on conventional Central Processing Unit (CPU) computing with Fortran. Considerable room for improvement still remains. / text

Les données sismologiques enregistrées sur les îles océaniques offrent l'opportunité d'analyser la houle via sa signature dans le "bruit" microsismique. Nous avons d'abord analysé les sources de bruit « secondaires », qui sont générées par l'interaction entre des vagues de même période dans une tempête, un cyclone ou par le phénomène de réflexion des vagues sur la côte. L'analyse des « microséismes secondaires », à l'échelle du bassin océanique, permet d'assurer un suivi spatio-temporel de la source qui les génère, même si elle est distante de plusieurs milliers de km des stations sismiques d'enregistrement. À plus long terme, leur étude permet d'assurer un suivi climatique global de l'activité des vagues dans une région donnée. Nous avons également étudié les sources de bruit « primaires », qui sont générées par l'interaction directe de la houle avec la côte. L'analyse des « microséismes primaires » permet de caractériser la houle localement par son amplitude, sa période et sa direction de propagation. Il est alors possible d'utiliser un capteur sismique comme substitut d'un houlographe. Dans le contexte des îles Éparses et de La Réunion, nous avons étudié plusieurs épisodes de houles extrêmes qui les touchent fréquemment et qui peuvent avoir d'importantes répercutions sociales ou environnementales. Enfin, en utilisant un réseau temporaire de stations sismologiques à La Réunion, nous avons analysé les variations du niveau de bruit microsismique pour caractériser l'impact des houles extrêmes sur les différentes façades de l'île. Cette étude permet d'identifier les sites les plus exposés aux vagues qui présentent un intérêt pour la récupération de cette énergie renouvelable encore trop sous-exploitée. / Seismic data recorded on oceanic islands can be used to analyze the swell through its signature in the microseismic noise. We first analysed the "secondary" noise sources, which are generated by the interaction of ocean waves with similar period within a storm, a cyclone or by the reflection phenomena off the coast. By analyzing secondary microseisms at the scale of the oceanic basin, we first performed a spatio-temporal tracking of the source, even localized thousands km off the recording seismic stations. Secondly, in the long-term, they can be used to follow the global climate change related to the ocean waves activity in a specific region. We also studied the "primary" seismic noise sources, which are created by the direct interaction of ocean waves with the coast. By analyzing these primary microseisms, we succeeded to characterize ocean waves locally in terms of amplitude, period, and, sometimes, direction of propagation. This showed that it is possible to use a seismic station as an ocean waves gauge to get precious swell data, particularly in remote and poorly instrumented areas. In the context of the Iles Éparses and of La Réunion Island, we studied a few extreme swells, which occur there frequently and can have strong social and environmental consequences. Finally, by using a temporary network of seismic stations installed in La Reunion, we studied the spatial variations of the seismic noise level across the island to characterize the swell impact on the different coasts. This study enables us to identify the most exposed spots to the swell, which may be attractive for generating renewable energy with this powerful resource yet underexploited.

The master's thesis consists of two parts. The first part contains the azimuth -- slowness analysis for period 3s<=T<=6s, which was conducted by the FK analysis for 10 stations of the WEBNET network. In the second part there is the surface wave group velocity measurement for every station pair of 10 WEBNET stations. The interstation group velocity was obtained by the cross--correlation of microseismic noise records. Local group velocities are connected with the instantaneous period data. The Love wave group velocities are visualised by the 2D tomography calculation. The thesis also contains the additional Love waves dispersion analysis for just one station pair. Powered by TCPDF (www.tcpdf.org)