摘要:
A multi-stage FWI workflow uses multiple-contaminated FWI models to predict surface-related multiples. A method embodying the present technological advancement, can include: using data with free surface multiples as input into FWI; generating a subsurface model by performing FWI with the free-surface boundary condition imposed on top of the subsurface model; using inverted model from FWI to predict multiples; removing predicted multiples from the measured data; using the multiple-free data as input into FWI with absorbing boundary conditions imposed on top of the subsurface model; and preparing a multiple free data set for use in conventional seismic data processing.
摘要:
A multi-stage FWI workflow uses multiple-contaminated FWI models to predict surface-related multiples. A method embodying the present technological advancement, can include: using data with free surface multiples as input into FWI; generating a subsurface model by performing FWI with the free-surface boundary condition imposed on top of the subsurface model; using inverted model from FWI to predict multiples; removing predicted multiples from the measured data; using the multiple-free data as input into FWI with absorbing boundary conditions imposed on top of the subsurface model; and preparing a multiple free data set for use in conventional seismic data processing.
摘要:
A method, including performing, with a computer, up/down separation of geophysical data, which produces an approximate up-going wavefield and an approximate down-going wavefield; creating an areal source based at least in part on the down-going wavefield; and performing, with a computer, a full wavefield inversion process with the areal source, and an objective function measuring a misfit between modeled up-going wavefields and recorded up-going wavefields, wherein the full wavefield inversion process generates a final subsurface physical property model.
摘要:
A method, including performing, with a computer, up/down separation of geophysical data, which produces an approximate up-going wavefield and an approximate down-going wavefield; creating an areal source based at least in part on the down-going wavefield; and performing, with a computer, a full wavefield inversion process with the areal source, and an objective function measuring a misfit between modeled up-going wavefields and recorded up-going wavefields, wherein the full wavefield inversion process generates a final subsurface physical property model.
摘要:
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.
摘要:
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.
摘要:
Method for improving convergence in gradient-based iterative inversion of seismic data (101), especially advantageous for full wavefield inversion. The method comprises decomposing the gradient into two (or more) components (103), typically the migration component and the tomographic component, then weighting the components to compensate for unequal frequency content in the data (104), then recombining the weighted components (105), and using the recombined gradient to update (106) the physical properties model (102).
摘要:
A computer-implemented method for generating missing frequencies within geophysical data, said method including: decomposing, with a computer, geophysical data into frequencies and slowness or wavenumber along time and one or more spatial dimensions; estimating, with a computer, a filter across known frequencies within the geophysical data for each time, spatial dimension, and slowness or wavenumber sample; extrapolating or interpolating, with a computer, the missing frequencies from the known frequencies by applying the filter; recomposing, with a computer, the known frequencies and the missing frequencies back into time and the one or more spatial dimensions; and prospecting for hydrocarbons with geophysical data that includes the missing frequencies.
摘要:
Method for reducing computational time in inversion of geophysical data to infer a physical property model (91), especially advantageous in full wavefield inversion of seismic data. An approximate Hessian is pre-calculated by computing the product of the exact Hessian and a sampling vector composed of isolated point diffractors (82), and the approximate Hessian is stored in computer hard disk or memory (83). The approximate Hessian is then retrieved when needed (99) for computing its product with the gradient (93) of an objective function or other vector. Since the approximate Hessian is very sparse (diagonally dominant), its product with a vector may therefore be approximated very efficiently with good accuracy. Once the approximate Hessian is computed and stored, computing its product with a vector requires no simulator calls (wavefield propagations) at all. The pre-calculated approximate Hessian can also be reused in the subsequent steps whenever necessary.
摘要:
A computer-implemented method for updating subsurface models including: using an offset continuation approach to update the model, and at each stage defining a new objective function where a maximum offset for each stage is set, wherein the approach includes, performing a first stage iterative full wavefield inversion with near offset data, as the maximum offset, to obtain velocity and density or impedance models, performing subsequent stages of iterative full wavefield inversion, each generating updated models, relative to a previous stage, wherein the subsequent stages include incrementally expanding the maximum offset until ending at a full offset, wherein a last of the stages yields finally updated models, the subsequent stages use the updated models as starting models, and the full wavefield inversions include constraining scales of the velocity model updates at each stage of inversion as a function of velocity resolution; and using the finally updated models to prospect for hydrocarbons.