Jerry Harris | Stanford University (original) (raw)

Papers by Jerry Harris

Water Resources Research, 2020

We describe a novel framework for estimating subsurface properties, such as rock permeability and... more We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time‐lapse observed seismic data by coupling full‐waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high‐performance FWI code; (2) at the rock physics level, we used built‐in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three‐level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our nu...

SEG Technical Program Expanded Abstracts 2009, 2009

We present a 3-D finite-element time-domain (FETD) algorithm for the simulation of electromagneti... more We present a 3-D finite-element time-domain (FETD) algorithm for the simulation of electromagnetic (EM) diffusion phenomena. The algorithm simulates transient electric fields and time derivatives of the magnetic fields for a general anisotropic earth. In order to compute transient fields, the electric field wave equation is transformed into a system of ordinary differential equations (ODE) via a Galerkin method with Dirichlet boundary conditions. To ensure both numerical stability and an efficient time step size, the system of ODE is discretized in time using the implicit backward Euler scheme. The resultant FETD matrix-vector equation is solved using a sparse direct solver with a fill-in reducing algorithm. When advancing the solution in time, the algorithm adjusts the tine step by examining if or not a current step size can be doubled without affecting the accuracy of the solution. Instead of directly solving another FETD matrix-vector equation for transient magnetic fields, Faraday's law is employed to compute time-derivatives of magnetic fields only at receiver positions. The accuracy and efficiency of the FETD algorithm are demonstrated using time-domain controlled source EM (TD-CSEM) simulations.

10th International Congress of the Brazilian Geophysical Society & EXPOGEF 2007, Rio de Janeiro, Brazil, 19-23 November 2007, 2007

SEG Technical Program Expanded Abstracts 2007, 2007

Time-lapse imaging has been successfully applied by oil industry for reservoir monitoring along t... more Time-lapse imaging has been successfully applied by oil industry for reservoir monitoring along the years. Another relevant application using similar theoretical background although with different purposes, is to monitor CO 2 sequestration to detect possible leaks that may cause environmental impact. We propose to reduce the number of parameters and consequently to simplify the update surveys by two complementary key concepts: (1) keeping relevant information from previous surveys to be updated by a smaller new survey; (2) using adaptive trigonal mesh that is finer at predicted time-varying regions and coarser otherwise. To achieve our goals, we proposed an approach to integrate data from previous surveys joint with an adaptive trigonal mesh framework applied to diffraction tomography. We used these concepts to develop an inversion algorithm for CO 2 sequestration monitoring, which successfully imaged a synthetic model using VSP geometry.

Dynamic imaging provides an effective way to integrate previous time-lapse data with current data... more Dynamic imaging provides an effective way to integrate previous time-lapse data with current data in order to estimate a current model. We present two methods useful for continuous dynamic reservoir and subsurface monitoring. The first method is a stochastic method for the inversion of geophysical data. It uses the ensemble Kalman filter (EnKF) that we demonstrate using time-lapse surface seismic data. We use reflection waveform and traveltime data plus sonic log information as input to the EnKF inversion. Unlike conventional seismic methods that give relative impedance, the EnKF inversion yields absolute seismic velocity. The second method is a specific dynamic imaging algorithm we call DynaSIRT. DynaSIRT integrates data from previous surveys in an efficient way, without the reprocessing older data. It preserves state variables that store the temporally damped effective illumination of previous surveys and timestamps to track parameter updates. We show results of DynaSIRT applied t...

This proposal addresses two related aspects of the acoustic reverberation problem, namely (1) sei... more This proposal addresses two related aspects of the acoustic reverberation problem, namely (1) seismic-acoustic properties of ocean bottom/subbottom sediments, and (2) the short-range scattering of acoustic/elastic waves at the rough and heterogeneous bottom/subbottom interface zone. The primary objective of the proposed research is to quantify the importance of bottom/subbottom scattering of compressional and converted shear waves. The approach is to numerically simulate short-range scattering using realistic models for the geometry, scales, and vertical and lateral variations in material properties of the ocean bottom and especially subbottom. To treat this problem realistically, a numerical model for constrained sedimentological deposition of sands and clays over basement rock is used to describe the threedimensional geometry of sediment beds and laminae over basement rock. Models for rock properties are then used to determine the seismic velocity, attenuation, density, porosity, and permeability for the suspended and loosely consolidated granular sediment materials and basement rocks. This geological model is then used in three-dimensional seismic-acoustic simulators to predict the acoustic fields scattered by the rocks at and below the interface. Two preliminary finding of results by us to date are noteworthy: First, acoustic scattering is completely inadequate for describing the interaction, particularly near grazing angles where the elastic boundary conditions play a critical role in establishing the energy balance. Second, bathymetrically "smooth" surfaces may be acoustically very "rough", providing a mechanism for the "anomalous" scattering regions observed in reverberation data and the importance of subsurface scattering in the reverberation problem. Our simulation codes are implemented on a dedicated in-house supercomputer, thus providing a powerful and cost effective computer platform for performing the proposed research.

SEG Technical Program Expanded Abstracts 2009, 2009

A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DAR... more A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DARS) measurements. In these experiments the bulk modulus of a sample is determined by the change in the resonance frequency of a tube due to the introduction of the sample. DARS experiments on porous samples confirm the perturbation theory by Morse and Ingard. For sealed porous samples, good agreement between the Biot-Gassmann modulus and the DARS bulk modulus was obtained. For open-pore samples the DARS bulk modulus is governed by the relative fluid motion at the outer wall. Measurements based on more realistic in-situ conditions (pressurized samples with multi-constituent properties) should lead to revisions of the theory, and enhance the use of DARS for reservoir characterization.

SEG Technical Program Expanded Abstracts 2001, 2001

Single-well seismic imaging has recently emerged as a valuable tool for delineating well-parallel... more Single-well seismic imaging has recently emerged as a valuable tool for delineating well-parallel features within structurally complex regions. Forward modeling can provide crucial guidance in both the design of single-well experiments and in the later phases of processing and ...

A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DAR... more A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DARS) measurements. In these experiments the bulk modulus of a sample is determined by the change in the resonance frequency of a tube due to the introduction of the sample. DARS experiments on porous samples confirm the perturbation theory by Morse and Ingard. For sealed porous samples, good agreement between the Biot-Gassmann modulus and the DARS bulk modulus was obtained. For open-pore samples the DARS bulk modulus is governed by the relative fluid motion at the outer wall. Measurements based on more realistic in-situ conditions (pressurized samples with multi-constituent properties) should lead to revisions of the theory, and enhance the use of DARS for reservoir characterization.

5th International Congress of the Brazilian Geophysical Society

Brazilian Journal of Geophysics

ABSTRACT. Seismic tomography is used in reservoir geophysics as an important method for high-reso... more ABSTRACT. Seismic tomography is used in reservoir geophysics as an important method for high-resolution imaging. The classical Born approach, which is used in single-frequency diffraction tomography under the condition of weak scattering, is limited by the requirement to know the background velocity in advance. We propose tomographic inversion approaches within matrix formalism and the Born approximation conditions. These approaches are iterative (in the sense that the background velocity field is updated at each iteration) and do not require knowledge of the true background velocity. In the first approach, a single-frequency that is kept constant is used. In the second approach, several frequencies are also kept constant and are used simultaneously. In the third approach, in addition to the background velocity, the working frequency is also updated. Finally, in the last approach, the multiple frequencies used simultaneously are updated throughout the iteration. The proposed approac...

Attenuation-difference, borehole-radar tomography was used to monitor a series of sodium chloride... more Attenuation-difference, borehole-radar tomography was used to monitor a series of sodium chloride tracer injection tests conducted within the FSE wellfield at the U.S. Geological Survey Fractured-Rock Hydrology Research Site in Grafton County, New Hampshire, USA. Borehole-radar tomography surveys were conducted by using the sequential-scanning and injection method in three boreholes that form a triangular prism of adjoining tomographic image planes.

Ph D Thesis California Inst of Tech Pasadena, 1980

True 4-D seismic subsurface monitoring for CO2 storage in Coal is examined. Our GCEP project is d... more True 4-D seismic subsurface monitoring for CO2 storage in Coal is examined. Our GCEP project is divided into two tasks: (1) Innovative scenarios for continuous seismic monitoring, and (2) development of a laboratory measurement for the signature of CO2 in coal. Under task #1, we have carried out a comprehensive simulation study that involved reservoir simulation, seismic simulation, and simulated

Journal of Geophysical Research: Solid Earth, 2015

DOI to the publisher's website. • The final author version and the galley proof are versions of t... more DOI to the publisher's website. • The final author version and the galley proof are versions of the publication after peer review. • The final published version features the final layout of the paper including the volume, issue and page numbers. Link to publication General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal. If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the "Taverne" license above, please follow below link for the End User Agreement:

Dynamic imaging provides an effective way to integrate previous time-lapse data with current data... more Dynamic imaging provides an effective way to integrate previous time-lapse data with current data in order to estimate a current model. We present two methods useful for continuous dynamic reservoir and subsurface monitoring. The first method is a stochastic method for the inversion of geophysical data. It uses the ensemble Kalman filter (EnKF) that we demonstrate using time-lapse surface seismic

Water Resources Research, 2020

We describe a novel framework for estimating subsurface properties, such as rock permeability and... more We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time‐lapse observed seismic data by coupling full‐waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high‐performance FWI code; (2) at the rock physics level, we used built‐in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three‐level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our nu...

SEG Technical Program Expanded Abstracts 2009, 2009

We present a 3-D finite-element time-domain (FETD) algorithm for the simulation of electromagneti... more We present a 3-D finite-element time-domain (FETD) algorithm for the simulation of electromagnetic (EM) diffusion phenomena. The algorithm simulates transient electric fields and time derivatives of the magnetic fields for a general anisotropic earth. In order to compute transient fields, the electric field wave equation is transformed into a system of ordinary differential equations (ODE) via a Galerkin method with Dirichlet boundary conditions. To ensure both numerical stability and an efficient time step size, the system of ODE is discretized in time using the implicit backward Euler scheme. The resultant FETD matrix-vector equation is solved using a sparse direct solver with a fill-in reducing algorithm. When advancing the solution in time, the algorithm adjusts the tine step by examining if or not a current step size can be doubled without affecting the accuracy of the solution. Instead of directly solving another FETD matrix-vector equation for transient magnetic fields, Faraday's law is employed to compute time-derivatives of magnetic fields only at receiver positions. The accuracy and efficiency of the FETD algorithm are demonstrated using time-domain controlled source EM (TD-CSEM) simulations.

10th International Congress of the Brazilian Geophysical Society & EXPOGEF 2007, Rio de Janeiro, Brazil, 19-23 November 2007, 2007

SEG Technical Program Expanded Abstracts 2007, 2007

Time-lapse imaging has been successfully applied by oil industry for reservoir monitoring along t... more Time-lapse imaging has been successfully applied by oil industry for reservoir monitoring along the years. Another relevant application using similar theoretical background although with different purposes, is to monitor CO 2 sequestration to detect possible leaks that may cause environmental impact. We propose to reduce the number of parameters and consequently to simplify the update surveys by two complementary key concepts: (1) keeping relevant information from previous surveys to be updated by a smaller new survey; (2) using adaptive trigonal mesh that is finer at predicted time-varying regions and coarser otherwise. To achieve our goals, we proposed an approach to integrate data from previous surveys joint with an adaptive trigonal mesh framework applied to diffraction tomography. We used these concepts to develop an inversion algorithm for CO 2 sequestration monitoring, which successfully imaged a synthetic model using VSP geometry.

Dynamic imaging provides an effective way to integrate previous time-lapse data with current data... more Dynamic imaging provides an effective way to integrate previous time-lapse data with current data in order to estimate a current model. We present two methods useful for continuous dynamic reservoir and subsurface monitoring. The first method is a stochastic method for the inversion of geophysical data. It uses the ensemble Kalman filter (EnKF) that we demonstrate using time-lapse surface seismic data. We use reflection waveform and traveltime data plus sonic log information as input to the EnKF inversion. Unlike conventional seismic methods that give relative impedance, the EnKF inversion yields absolute seismic velocity. The second method is a specific dynamic imaging algorithm we call DynaSIRT. DynaSIRT integrates data from previous surveys in an efficient way, without the reprocessing older data. It preserves state variables that store the temporally damped effective illumination of previous surveys and timestamps to track parameter updates. We show results of DynaSIRT applied t...

This proposal addresses two related aspects of the acoustic reverberation problem, namely (1) sei... more This proposal addresses two related aspects of the acoustic reverberation problem, namely (1) seismic-acoustic properties of ocean bottom/subbottom sediments, and (2) the short-range scattering of acoustic/elastic waves at the rough and heterogeneous bottom/subbottom interface zone. The primary objective of the proposed research is to quantify the importance of bottom/subbottom scattering of compressional and converted shear waves. The approach is to numerically simulate short-range scattering using realistic models for the geometry, scales, and vertical and lateral variations in material properties of the ocean bottom and especially subbottom. To treat this problem realistically, a numerical model for constrained sedimentological deposition of sands and clays over basement rock is used to describe the threedimensional geometry of sediment beds and laminae over basement rock. Models for rock properties are then used to determine the seismic velocity, attenuation, density, porosity, and permeability for the suspended and loosely consolidated granular sediment materials and basement rocks. This geological model is then used in three-dimensional seismic-acoustic simulators to predict the acoustic fields scattered by the rocks at and below the interface. Two preliminary finding of results by us to date are noteworthy: First, acoustic scattering is completely inadequate for describing the interaction, particularly near grazing angles where the elastic boundary conditions play a critical role in establishing the energy balance. Second, bathymetrically "smooth" surfaces may be acoustically very "rough", providing a mechanism for the "anomalous" scattering regions observed in reverberation data and the importance of subsurface scattering in the reverberation problem. Our simulation codes are implemented on a dedicated in-house supercomputer, thus providing a powerful and cost effective computer platform for performing the proposed research.

SEG Technical Program Expanded Abstracts 2009, 2009

A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DAR... more A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DARS) measurements. In these experiments the bulk modulus of a sample is determined by the change in the resonance frequency of a tube due to the introduction of the sample. DARS experiments on porous samples confirm the perturbation theory by Morse and Ingard. For sealed porous samples, good agreement between the Biot-Gassmann modulus and the DARS bulk modulus was obtained. For open-pore samples the DARS bulk modulus is governed by the relative fluid motion at the outer wall. Measurements based on more realistic in-situ conditions (pressurized samples with multi-constituent properties) should lead to revisions of the theory, and enhance the use of DARS for reservoir characterization.

SEG Technical Program Expanded Abstracts 2001, 2001

Single-well seismic imaging has recently emerged as a valuable tool for delineating well-parallel... more Single-well seismic imaging has recently emerged as a valuable tool for delineating well-parallel features within structurally complex regions. Forward modeling can provide crucial guidance in both the design of single-well experiments and in the later phases of processing and ...

A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DAR... more A set of 23 oil-saturated samples were used for Differential Acoustic Resonance Spectroscopy (DARS) measurements. In these experiments the bulk modulus of a sample is determined by the change in the resonance frequency of a tube due to the introduction of the sample. DARS experiments on porous samples confirm the perturbation theory by Morse and Ingard. For sealed porous samples, good agreement between the Biot-Gassmann modulus and the DARS bulk modulus was obtained. For open-pore samples the DARS bulk modulus is governed by the relative fluid motion at the outer wall. Measurements based on more realistic in-situ conditions (pressurized samples with multi-constituent properties) should lead to revisions of the theory, and enhance the use of DARS for reservoir characterization.

5th International Congress of the Brazilian Geophysical Society

Brazilian Journal of Geophysics

ABSTRACT. Seismic tomography is used in reservoir geophysics as an important method for high-reso... more ABSTRACT. Seismic tomography is used in reservoir geophysics as an important method for high-resolution imaging. The classical Born approach, which is used in single-frequency diffraction tomography under the condition of weak scattering, is limited by the requirement to know the background velocity in advance. We propose tomographic inversion approaches within matrix formalism and the Born approximation conditions. These approaches are iterative (in the sense that the background velocity field is updated at each iteration) and do not require knowledge of the true background velocity. In the first approach, a single-frequency that is kept constant is used. In the second approach, several frequencies are also kept constant and are used simultaneously. In the third approach, in addition to the background velocity, the working frequency is also updated. Finally, in the last approach, the multiple frequencies used simultaneously are updated throughout the iteration. The proposed approac...

Attenuation-difference, borehole-radar tomography was used to monitor a series of sodium chloride... more Attenuation-difference, borehole-radar tomography was used to monitor a series of sodium chloride tracer injection tests conducted within the FSE wellfield at the U.S. Geological Survey Fractured-Rock Hydrology Research Site in Grafton County, New Hampshire, USA. Borehole-radar tomography surveys were conducted by using the sequential-scanning and injection method in three boreholes that form a triangular prism of adjoining tomographic image planes.

Ph D Thesis California Inst of Tech Pasadena, 1980

True 4-D seismic subsurface monitoring for CO2 storage in Coal is examined. Our GCEP project is d... more True 4-D seismic subsurface monitoring for CO2 storage in Coal is examined. Our GCEP project is divided into two tasks: (1) Innovative scenarios for continuous seismic monitoring, and (2) development of a laboratory measurement for the signature of CO2 in coal. Under task #1, we have carried out a comprehensive simulation study that involved reservoir simulation, seismic simulation, and simulated

Journal of Geophysical Research: Solid Earth, 2015

DOI to the publisher's website. • The final author version and the galley proof are versions of t... more DOI to the publisher's website. • The final author version and the galley proof are versions of the publication after peer review. • The final published version features the final layout of the paper including the volume, issue and page numbers. Link to publication General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal. If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the "Taverne" license above, please follow below link for the End User Agreement:

Dynamic imaging provides an effective way to integrate previous time-lapse data with current data... more Dynamic imaging provides an effective way to integrate previous time-lapse data with current data in order to estimate a current model. We present two methods useful for continuous dynamic reservoir and subsurface monitoring. The first method is a stochastic method for the inversion of geophysical data. It uses the ensemble Kalman filter (EnKF) that we demonstrate using time-lapse surface seismic