Search

We present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

We present a synthetic multi-scale, multi-physics dataset constructed from the Kimberlina 1.2 CO2 reservoir model based on a potential CO2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir-simulation scenario produces hydrologic-state models at the onset and after 20 years of CO2 injection. Subsequently, these models were transformed into geophysical properties, including P- and S-wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO2, and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time-lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real-world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time-lapse geophysical responses and well logs forms a multi-scale, multi-physics testbed for designing and testing geophysical CO2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

Multiphysics inversion exploits different types of geophysical data that often complement each other and aims to improve overall imaging resolution and reduce uncertainties in geophysical interpretation. Despite the advantages, traditional multiphysics inversion is challenging because it requires a large amount of computational time and intensive human interactions for preprocessing data and finding trade-off parameters. These issues make it nearly impossible for traditional multiphysics inversion to be applied as a real-time monitoring tool for geological carbon storage. In this paper, we present a deep learning (DL) multiphysics network for imaging CO2 saturation in real time. The multiphysics network consists of three encoders for analysing seismic, electromagnetic and gravity data and shares one decoder for combining imaging capabilities of the different geophysical data for better predicting CO2 saturation. The network is trained on pairs of CO2 label models and multiphysics data so that it can directly image CO2 saturation. We use the bootstrap aggregating method to enhance the imaging accuracy and estimate uncertainties associated with CO2 saturation images. Using realistic CO2 label models and multiphysics data derived from the Kimberlina CO2 storage model, we evaluate the performance of the deep learning multiphysics network and compare its imaging results to those from the deep learning single-physics networks. Our modelling experiments show that the deep learning multiphysics network for seismic, electromagnetic, and gravity data not only improves the imaging accuracy but also reduces uncertainties associated with CO2 saturation images. Our results also suggest that the deep learning multiphysics network for the non-seismic data (i.e., electromagnetic and gravity) can be used as an effective low-cost monitoring tool in between regular seismic monitoring.

Electromagnetics has been shown to be a viable tool to monitor CO2 plumes embedded in saline reservoirs. However, the majority of studies focus on measuring the electric field, which requires cumbersome equipment in the field and it is difficult to install permanent measurement stations. Magnetic field receivers offer an opportunity to reduce the form factor of the survey and increase the mobility by utilizing upcoming technologies, such as drones. We explore the use of frequency-domain electric dipole sources, and time-domain loops with a focus on measuring the secondary magnetic field at the surface for a conceptual injection scenario based on the Kemper CarbonSAFE site. We find that electric dipole sources give a response above the sensitivity of current sensor technology and, therefore, be a viable tool for CO2 monitoring. The time-domain loop source does provide fields that are useful for determining the location of the CO2 plume, however the field magnitude is below the sensitivity of the current generation of instruments. To explore the use of a potential borehole receiver we generate a map of the magnetic field at depth to explore potential borehole placement for monitoring efforts. Finally, we limit the spatial extent of the electric dipole survey to a single parcel of land to help understand how the fields change with survey geometry. We find that the shape of the secondary fields change slightly with the small transmitter, but are still measurable provided that the cultural noise at the site is low. Thus, we conclude that at the Kemper site a frequency-domain cross-dipole source with magnetometer receivers is suitable to monitor the expansion of the CO2 plume in the saline reservoir, even with a limited transmitter footprint on the surface.

The basin environment is a widely studied subject in both geology and geophysics for its economic significance in energy and mineral explorations. However, the estimation of the basement depth is often a challenging task given the complexity of the basement relief and lateral physical property change. Previous works simplify the problem by only inverting for the depth to the basement, and more recent studies have suggested the need to incorporate the variation of physical properties to improve basement structure imaging. In this study, we develop an inversion method with the associated workflow to simultaneously recover both the depth to a magnetic basement and a laterally varying magnetic susceptibility in the basement rock. To achieve this, we employ a set of constraints on the inverse problem. Particularly, both the recovered susceptibility and basement depth models are bounded below a possible maximum value, and the depth model is guided by a few depth points obtained from the resistivity models that are obtained from the one‐dimensional blocky inversions of magnetotelluric (MT) data. In addition, we apply the fuzzy C‐means (FCM) clustering to the susceptibility model during the inversion and use the inverted cluster centers to differentiate for different geological units in the basement. To show the effectiveness of our work, we compare the existing approaches and our method using two test inversions on one synthetic model resembling the basin–basement environment before demonstrating our method on a field data example with magnetic data collected by the U.S. Geological Survey (USGS) over the Illinois Basin. Our results show improved recovery in both basement relief and susceptibility in the basement rock, and inversion with field data is able to identify three different susceptibility zones in basement rock below the Illinois Basin. [ABSTRACT FROM AUTHOR]

We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO2 reservoir model based on a potential CO2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO2, and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website. [ABSTRACT FROM AUTHOR]

Data acquisition using equal spacing has been a standard practice in geophysics. The dense uniform sampling derived from Nyquist–Shannon sampling includes redundant samples, and it is sufficient but not necessary to adequately record target signals. We propose an ergodic sampling, which avoids the redundant samples in the dense uniform sampling and possesses the ability to capture the sufficiently similar information content as does Nyquist sampling. Ergodicity means that a key part of the system can represent the average performance of the entire system. Our ergodic samples are a critical subset of dense uniform samples and can represent the full uniform Nyquist samples. To find such a critical subset, we first examine the properties of different sampling patterns, including the sampling interval distribution, sampling angle distribution, areal sample density, and resolution in the spectral domain. We define the information sampling ability of sampling patterns based on these properties. The concept of information sampling ability that we have proposed serves as the criterion to compare different sampling patterns and assess their sampling performances. The sampling patterns with the same information sampling ability have the same capability to gather information, even though the appearance of sampling patterns may be different. We formulate an optimization problem to find this critical subset of sample locations, which has the fewest number of samples but has a similar information sampling ability as that of the desired dense uniform samples. This critical subset is irregularly located, has the optimized properties and forms the ergodic sampling pattern. We define this process of sampling design and associated understanding as the ergodic sampling. Ergodic sampling can be applied to gain two major benefits in practice. First, this approach can save a significant number of samples. We demonstrate ergodic sampling using one‐dimensional synthetic data and a two‐dimensional field geophysical dataset. The simulations confirm that, compared with other sampling strategies, ergodic sampling can use fewer samples to acquire the same amount of information, so that we can save cost. Alternatively, with the same budget, we can use the same number of samples through ergodic sampling to acquire more information. The new ergodic sampling can lead to a new generation of economic and efficient geophysical data acquisition, which could assist in increasing the discovery rate in resource exploration, tackling more earth science problems with a limited budget and can also benefit the environment in the process by reducing the invasiveness in potentially sensitive regions. [ABSTRACT FROM AUTHOR]

The field of geophysics faces the daunting task of monitoring complex reservoir dynamics and imaging carbon dioxide storage up to several decades into the future. This presents numerous challenges, including sensitivity to parameter changes, resolution of obtained results and the cost of long‐term deployment. To effectively store CO2 subsurface, it is necessary to monitor and account for the injected CO2. The gravity method provides several advantages for CO2 monitoring, as changes in fluid saturation correspond directly and uniquely to observed density changes. Three‐axis borehole gravity has demonstrated significant promise as a next‐generation tool for reliably monitoring reservoir dynamics across a range of depths and sizes. However, the gravity inverse problem is highly ill‐posed, necessitating regularization that incorporates prior knowledge. To address this issue, we propose using a feed‐forward neural network, a machine learning method, to invert time‐lapse three‐axis borehole gravity data and monitor CO2 movement within a reservoir. By training the neural network on models that analyse changes in density and corresponding gravity responses resulting from perturbations made to the reservoir model, we can create scenarios that train the algorithm to identify unexpected CO2 migration in addition to the normal movement of CO2. Our method is demonstrated using reservoir models for the Johansen formation in offshore Norway. We convert reservoir saturation models into density changes and generate their corresponding three‐axis gravity data in a set of boreholes. Our results show that the developed machine learning inversion algorithm has high reliability and resolution for imaging density change associated with CO2 plumes, as demonstrated in the Johansen reservoir models utilized by the simulator. We also investigate machine learning inversion using regularization parameters and show that it is robust, with a strong tolerance for higher levels of noise. Our study demonstrates that the developed machine learning algorithm is a powerful tool for inverting three‐axis borehole gravity data and monitoring the migration and long‐term storage of injected CO2. [ABSTRACT FROM AUTHOR]

Multiphysics inversion exploits different types of geophysical data that often complement each other and aims to improve overall imaging resolution and reduce uncertainties in geophysical interpretation. Despite the advantages, traditional multiphysics inversion is challenging because it requires a large amount of computational time and intensive human interactions for preprocessing data and finding trade‐off parameters. These issues make it nearly impossible for traditional multiphysics inversion to be applied as a real‐time monitoring tool for geological carbon storage. In this paper, we present a deep learning (DL) multiphysics network for imaging CO2 saturation in real time. The multiphysics network consists of three encoders for analysing seismic, electromagnetic and gravity data and shares one decoder for combining imaging capabilities of the different geophysical data for better predicting CO2 saturation. The network is trained on pairs of CO2 label models and multiphysics data so that it can directly image CO2 saturation. We use the bootstrap aggregating method to enhance the imaging accuracy and estimate uncertainties associated with CO2 saturation images. Using realistic CO2 label models and multiphysics data derived from the Kimberlina CO2 storage model, we evaluate the performance of the deep learning multiphysics network and compare its imaging results to those from the deep learning single‐physics networks. Our modelling experiments show that the deep learning multiphysics network for seismic, electromagnetic, and gravity data not only improves the imaging accuracy but also reduces uncertainties associated with CO2 saturation images. Our results also suggest that the deep learning multiphysics network for the non‐seismic data (i.e., electromagnetic and gravity) can be used as an effective low‐cost monitoring tool in between regular seismic monitoring. [ABSTRACT FROM AUTHOR]