Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. This contribution investigates the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. This contribution investigates the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. The time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes

Hydro-frac monitoring using ground time-domain electromagnetics

As motivation for considering new electromagnetic techniques for hydraulic fracture monitoring, we develop a simple financial model for the net present value offered by geophysical characterization to reduce the error in stimulated reservoir volume calculations. This model shows that even a 5% improvement in stimulated reservoir volume for a 1 billion barrel (bbl) field results in over 1 billion U.S. dollars (US$) in net present value over 24 years for US$100/bbl. oil and US$0.5 billion for US$50/bbl. oil. The application of conductivity upscaling, often used in electromagnetic modeling to reduce mesh size and thus simulation runtimes, is shown to be inaccurate for the high electrical contrasts needed to represent steel-cased wells in the earth. Fine-scale finite-difference modeling with 12.22-mm cells to capture the steel casing and fractures shows that the steel casing provides a direct current pathway to a created fracture that significantly enhances the response compared with neglecting the steel casing. We consider conductively enhanced proppant, such as coke-breeze-coated sand, and a highly saline brine solution to produce electrically conductive fractures. For a relatively small frac job at a depth of 3 km, involving 5,000 bbl. of slurry and a source midpoint to receiver separation of 50 m, the models show that the conductively enhanced proppant produces a 15% increase in the electric field strength (in-line with the transmitter) in a 10-Ωm background. In a 100-Ωm background, the response due to the proppant increases to 213%. Replacing the conductive proppant by brine with a concentration of 100,000-ppm NaCl, the field strength is increased by 23% in the 100-Ωm background and by 2.3% in the 10-Ωm background. All but the 100,000-ppm NaCl brine in a 10-Ωm background produce calculated fracture-induced electric field increases that are significantly above 2%, a value that has been demonstrated to be observable in field measurements. © 2015 European Association of Geoscientist

Distributed electric field sensing using fibre optics in borehole environments

In the past decade, rapid advances in distributed optical fibre sensing technologies have made it possible to record various geophysical data (e.g. strain, temperature and pressure) continuously in both time and space along the fibre, providing an unprecedented quantity and spatial density of data compared to traditional geophysical measurements as well as reducing data acquisition cost. To date, no distributed fibre-based electromagnetic field sensing system has been implemented although electromagnetic sensing could have a broad range of applications to geophysical imaging and monitoring in borehole environments. The goal of this paper is to provide a theoretical feasibility study regarding the design and use of an electromagnetic sensing optical fibre for geophysical applications. First, we present the sensitivity analysis of a ‘hypothetical’ optical fibre coated with polyvinylidene fluoride, a polymer that provides relatively high piezoelectric properties, yet unlike ceramics, is flexible. Using a two-dimensional electromagnetic modelling algorithm, we simulate the earth electric-field-to-fibre-strain transfer function and estimate the theoretical sensitivity of the optical fibre to electric fields. Given the state-of-the-art distributed acoustic sensing strain sensitivities in the picometres strain range, our numerical modelling analysis suggests that a perfectly coupled polyvinylidene fluoride–coated optical fibre can measure electric field values in the mV/m to V/m amplitude range. We then apply a cylindrically symmetric modelling algorithm to simulate numerical models demonstrating the applicability of such a fibre in an oilfield environment. Scenarios investigated employ an electric field source and suggest that the measurements can be used to distinguish the oil versus water ratio with a fibre mounted inside a producing steel cased oil well as well as distinguishing between brine and hydrocarbon filled reservoir zones with a fibre located outside of the casing

Magnetotelluric Investigations of the Kīlauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kīlauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kīlauea and the East and Southwest Rift zones, which has implications for understanding Kīlauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kīlauea. Melt fractions between 0.096 and 0.117 are calculated for the Kīlauea and Puʻuʻōʻō low resistivity anomalies, respectively

Seismic monitoring of well integrity

Seismic tube waves, produced by flexure of the well boundary, pressure changes in the fluid in the well, and deformation of the material immediately surrounding the well, are particularly sensitive to variations in the state of the well. We evaluate a direct approach for generating and observing tube waves as a means of detecting well damage. While we find that it can be difficult to reliably excite observable tube waves without a very strong surface source, time-frequency techniques can be employed to increase the detectability of tube wave reflections. New technologies, particularly distributed acoustic sensing, hold great promise for evaluating well integrity by monitoring tube waves, temperature changes, and seismic noise due to well deformation and fluid leakage

Enhanced Multi-Dimensional Inversion Through Target-Specific Inversion Parameter Bounds With an Application to Crosswell Electromagnetic for Sequestration Monitoring

In geophysical inversions, lower and upper model parameter bounds are a means of solution stabilization. Further, constraints that intend to let only geologically plausible inverse solutions pass are amenable to lower and upper bounds. Reliable prior information is paramount to construct such bound constraints. It is common practice to narrow and widen bound intervals for regions of, respectively, more and less certain prior information. Contrary to this practice, we experiment with widened bound intervals in zones that are poorly resolved by a given survey configuration but where prior information would suggest structural anomalies of interest. The purpose of enlarged parameter bounds that correlate spatially with predefined targets is to let the inversion explore a larger solution space, thus increasing the potential to resolve otherwise hidden anomalies. Application of the method is based on a carbon-sequestration baseline (pre-injection) crosswell electromagnetic (EM) field survey at the Containment and Monitoring Institute Field Research Station (Alberta, Canada), where impeded measurements led to generally reduced sensitivities for the interwell region. Synthetic-data proofs of concept use augmented bounds designed to boost the resolution of artificial plume targets, indicating an enhanced illumination compared to constant bounds. Comparative field data inversions with spatially variable bounds constructed from prior resistivity and velocity information highlight non-horizontal baseline structures