Study of time lapse data processing for dynamic hydrologic conditions
There are no author-identified significant results in this report.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
There are no author-identified significant results in this report.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Accurate estimation of the spatio-temporal variation in soil saturation and salt concentration is valuable for calibrating and evaluating Earth system models as well as detecting important eco-hydrological interactions. Such representative models play a vital role in investigating the stability of ecosystems facing changing hydrologic disturbance regimes. Ecosystem-scale field experiments are a useful way to capture hydrological disturbance and gain a deeper understanding of the complex mechanisms driving environmental changes. However, previous studies have not explored the use of quantitative imaging at a large spatial scale to reproduce breakthrough curves and inform system response and recovery following hydrologic disturbance events. Therefore, in this study, non-invasive geophysical methods were used to capture the spatio-temporal variation in subsurface saturation and salt concentrations during an ecosystem-scale coastal forest flooding experiment.
Seismic characterization of flow properties is a difficult prospect, due to the indirect relationship between permeability and seismic velocity and attenuation. At best, seismic time-lapse changes can detect the effects of saturation and pressure changes in a reservoir due to fluid flow. Even in this case, the interpretation of the seismic observation in terms of the state of the reservoir depends intimately on the properties of the rock physics model, which are usually poorly known. The onset time, the calendar time of geophysical changes, provides an alternative datum for characterizing properties such as reservoir permeability. The main advantage of an onset time is that it is sensitive to the flow properties of the reservoir yet insensitive to the details of the rock physics model. Two examples of the utility of onset times are discussed: The use of travel time shifts induced by the injection of carbon dioxide between two wells, and the time-lapse time shifts for elastic waves propagating through a reservoir undergoing enhanced oil recovery. In both examples, the onset times are mapped into permeability estimates using a trajectory-based approach akin to seismic tomography.
Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.
With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.
We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.
This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.