Abstract:
A method for iteratively inverting seismic data to jointly infer a model for at least P-wave velocity and attenuation parameters of the subsurface, the method including: jointly inverting the P-wave velocity and attenuation parameters with an iterative visco-acoustic full wavefield inversion process, wherein the iterative visco-acoustic full wavefield inversion process includes computing a gradient of an objective function, the objective function measuring a misfit between all or part of the seismic data and corresponding model-simulated seismic data; for each of the P-wave velocity and attenuation parameters, computing a search direction in model space from the gradient; determining line search step sizes α and β for the search directions for the P-wave velocity and attenuation parameters, respectively, wherein a ratio of the step sizes is a function of the P-wave velocity parameter; and using the step sizes α and β and the search directions for each of the P-wave velocity and attenuation parameters, computing a new search direction in model space, then performing a line search along the new search direction to arrive at a new step size, and using the new step size and the new search direction to generate an updated model for a current iteration of the iterative visco-acoustic full wavefield inversion process.
Abstract:
Embodiments described herein use stochastic inversion (460) in lower dimensions to form an initial model (458) that is to be used in higher-dimensional gradient-based inversion (466). For example, an initial model may be formed from 1.5-D stochastic inversions, which is then processed (464) to form a 3-D model. Stochastic inversions reduce or avoid local minima and may provide an initial result that is near the global minimum.
Abstract:
Method for using seismic data from earthquakes to address the low frequency lacuna problem in traditional hydrocarbon exploration methods. Seismometers with frequency response down to about 1 Hz are placed over a target subsurface region in an array with spacing suitable for hydrocarbon exploration (21). Data are collected over a long (weeks or months) time period (22). Segments of the data (44) are identified with known events from earthquake catalogs (43). Those data segments are analyzed using techniques such as traveltime delay measurements (307) or receiver function calculations (46) and then are combined with one or more other types of geophysical data acquired from the target region, using joint inversion (308-310) in some embodiments of the method, to infer physical features of the subsurface indicative of hydrocarbon potential or lack thereof (26).
Abstract:
A method for suppressing measurement system signature, or artifacts, that arise when controlled source electromagnetic survey data are inverted to obtain a resistivity image of a subsurface region. The method involves identifying regions (47) where the image has low or rapidly varying sensitivity to data acquired by a given receiver, typically regions close to and under the given receiver. Then, in the iterative inversion process where a resistivity model is updated to minimize an objective function, the model update is modified (48) to reduce the impact of such low sensitivity regions on the update.
Abstract:
A method, including: determining, with a computer, point spread functions for a plurality of parameter locations by performing at least a portion of a first iteration of an iterative full wavefield inversion process; determining at least one property for each of the point spread functions; and evaluating a candidate survey design based on the at least one property for each of the point spread functions.
Abstract:
A method for exploring for hydrocarbons, including: simulating a seismic waveform, using a computer, wherein computations are performed on a computational grid representing a subsurface region, said computational grid using perfectly matched layer (PML) boundary conditions that use an energy dissipation operator to minimize non-physical wave reflections at grid boundaries; wherein, in the simulation, the PML boundary conditions are defined to reduce computational instabilities at a boundary by steps including, representing direction of energy propagation by a Poynting vector, and dissipating energy, with the dissipation operator, in a direction of energy propagation instead of in a phase velocity direction; and using the simulated waveform in performing full waveform inversion or reverse time migration of seismic data, and using a physical property model from the inversion or a subsurface image from the migration to explore for hydrocarbons.
Abstract:
Method for efficient computation of wave equation migration angle gathers by using multiple imaging conditions. Common reflection angle or common azimuth gathers or gathers including both common reflection angles and common azimuth angles are produced as the data are migrated. In the course of either wave equation migration or reverse time migration, the pressures and particle motion velocities that need to be computed are sufficient to also compute the Poynting vector pointing in the direction of source-side (35) or receiver-side (37) wavefield propagation. From that, the reflection and azimuth angles can be computed (38). The seismic images can then be stored in the appropriate angle bins, from which common reflection angle or azimuth data volumes can be assembled (39).
Abstract:
A method for iteratively inverting seismic data to jointly infer a model for at least P-wave velocity and attenuation parameters of the subsurface, the method including: jointly inverting the P-wave velocity and attenuation parameters with an iterative visco-acoustic full wavefield inversion process, wherein the iterative visco-acoustic full wavefield inversion process includes computing a gradient of an objective function, the objective function measuring a misfit between all or part of the seismic data and corresponding model-simulated seismic data; for each of the P-wave velocity and attenuation parameters, computing a search direction in model space from the gradient; determining line search step sizes α and β for the search directions for the P-wave velocity and attenuation parameters, respectively, wherein a ratio of the step sizes is a function of the P-wave velocity parameter; and using the step sizes α and β and the search directions for each of the P-wave velocity and attenuation parameters, computing a new search direction in model space, then performing a line search along the new search direction to arrive at a new step size, and using the new step size and the new search direction to generate an updated model for a current iteration of the iterative visco-acoustic full wavefield inversion process.
Abstract:
Method for rapidly computing updates to frequency-domain seismic wave fields by utilizing a matrix perturbation approach. The method speeds up model (e.g., velocity) parameter estimation by iterative inversion of measured seismic data (21-27). The method applies to the line search where the optimal size of the model update is estimated by testing different size updates to see which one generates the minimum objective function. By treating the model update as a perturbation, perturbation theory is used to relate the model perturbation to a corresponding wavefield perturbation (35). Thus, the Helmholtz equation is solved only once per iteration cycle (22).
Abstract:
A basically time-domain method for performing full wavefield inversion of seismic data to infer a subsurface physical property model (61), where however at least one quantity required for the inversion, such as the Hessian of the cost function, is computed in the frequency domain (64). The frequency-domain quantity or quantities may be obtained at only a few discrete frequencies (62), preferably low frequencies, and may be computed on a coarse spatial grid, thus saving computing time with minimal loss in accuracy. For example, the simulations of predicted data and the broadband gradient of the objective function may be computed in the time domain (67), and the Hessian matrix, approximated by its diagonal, may be computed in the frequency domain. It may be preferable to use time-domain and the frequency-domain solvers that employ different numerical schemes, such as finite-difference method, one-way wave equation, finite-element method (63).