# US 2008/0189043 A1 · Direct Time Lapse Inversion of Seismic Data > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 2008/0189043 A1 - Document ID: US20080189043A1 - Status: Published application - Kind: A1 - Application number: 11/672,000 - Filed: 2007-02-06 - Published: 2008-08-07 - Earliest claimed: 2007-02-06 - Assignees: ConocoPhillips Company - Inventors: Phil Dean Anno; Partha Sarathi Routh - Web page: https://patentgrep.com/patent/US20080189043A1 ## CPC - G01V 1/308: {Time lapse or 4D effects, e.g. production related effects to the formation (fluid flow per se E21B47/00)} - G01V 2210/612: Previously recorded data, e.g. time-lapse or 4D - G01V: GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS (means for indicating the location of accidentally buried, e.g. snow-buried, persons A63B29/02) ## Abstract A method and system of seismic data processing includes acquiring first and second prestack seismic datasets wherein the seismic datasets have a time-lapse relationship. The first prestack seismic dataset is inverted to determine a first source wavelet data set and a first reflectivity parameters data set. A perturbation data set is obtained by simultaneously minimizing an Lp-Norm between a portion of the first reflectivity parameters data set and reflectivity parameters for the second data set, an Lp-Norm between the first source wavelet data set and a source wavelet for the second data set and an Lp-Norm for data misfit between the first prestack seismic data set and the second prestack seismic data set. The perturbation dataset may be inverted for petrophysical parameters and for determining time-lapse changes in a reservoir. ## Claims ### Claim 1 A method of seismic data processing comprising: a) acquiring a first prestack seismic dataset and a second prestack seismic dataset wherein the second prestack seismic dataset has a time-lapse relationship to the first prestack seismic dataset; b) inverting the first prestack seismic dataset for determination of a first source wavelet data set and a first reflectivity parameters data set; c) obtaining a perturbation data set by using a computer to simultaneously minimize: i) an Lp-Norm between a portion of the first reflectivity parameters data set and reflectivity parameters for the second data set; ii) an Lp-Norm between the first source wavelet data set and a source wavelet for the second data set; iii) an Lp-Norm for data misfit between the first prestack seismic data set and the second prestack seismic data set. ### Claim 2 (depends on 1) The method of claim 1 further comprising inverting the perturbation dataset for at least one petrophysical parameter. ### Claim 3 (depends on 1) The method of claim 1 further comprising inverting the perturbation dataset to obtain at least one of: i) impedance, ii) Poisson's ratio, iii) density, iv) fluid saturation, v) pressure, vi) permeability, vii) porosity, viii) shear modulus, ix) Lame's constant. and x) Young's modulus. ### Claim 4 (depends on 1) The method of claim 1 further comprising determining intercept and gradient parameters using the first and second source wavelet data sets and the first and second reflectivity parameters data sets. ### Claim 5 (depends on 1) The method of claim 1 wherein the first and second reflectivity parameters Lp-Norm is an L1-Norm. ### Claim 6 (depends on 1) The method of claim 1 wherein the first and second source wavelet data sets Lp-Norm is an L2-Norm. ### Claim 7 (depends on 1) The method of claim 1 wherein the first and second prestack seismic data sets Lp-Norm for data misfit is an L2-Norm. ### Claim 8 (depends on 1) The method of claim 1 further comprising displaying the perturbation dataset. ### Claim 9 (depends on 1) The method of claim 1 wherein inverting the first prestack seismic dataset for determination of the first source wavelet data set includes inverting for a source wavelet for each of a plurality of common image point gathers and determining an average source wavelet. ### Claim 10 (depends on 1) The method of claim 1 wherein inverting the second prestack seismic dataset for determination of the second source wavelet data set includes inverting for a source wavelet for each of a plurality of common image point gathers and determining an average source wavelet. ### Claim 11 A method of seismic data processing to obtain a time-lapse data volume comprising: a) acquiring a first prestack seismic dataset and a second prestack seismic dataset wherein the second prestack seismic dataset has a time-lapse relationship to the first prestack seismic dataset; b) inverting the first prestack seismic dataset for determination of a first source wavelet data set and a first reflectivity parameters data set; c) obtaining a perturbation data set by simultaneously minimizing: i) an L1-Norm between a portion of the first reflectivity parameters data set and reflectivity parameters for the second data set; ii) an L2-Norm between the first source wavelet data set and a source wavelet for the second data set; iii) an L2-Norm for data misfit between the first prestack seismic data set and the second prestack seismic data set; and d) recording the perturbation data set as a time-lapse data volume in a form for display. ### Claim 12 (depends on 11) The method of claim 11 further comprising inverting the perturbation dataset for petrophysical parameters. ### Claim 13 (depends on 11) The method of claim 11 further comprising inverting the perturbation dataset to obtain at least one of: i) impedance, ii) Poisson's ratio, iii) density, iv) fluid saturation, v) pressure, vi) permeability, vii) porosity, viii) shear modulus, ix) Lame's constant, and x) Young's modulus. ### Claim 14 (depends on 11) The method of claim 11 further comprising determining intercept and gradient parameters using the first and second source wavelet data sets and the first and second reflectivity parameters data sets. ### Claim 15 (depends on 11) The method of claim 11 wherein inverting the first prestack seismic dataset for determination of the first source wavelet data set includes inverting for a source wavelet for each of a plurality of common image point gathers and determining an average source wavelet. ### Claim 16 (depends on 11) The method of claim 11 wherein inverting the second prestack seismic dataset for determination of the second source wavelet data set includes inverting for a source wavelet for each of a plurality of common image point gathers and determining an average source wavelet. ### Claim 17 A set of application program interfaces embodied on a computer readable medium for execution by a processor for calculating a direct time-lapse inversion data set to locate temporal changes in subsurface reservoirs comprising: a) a first interface that receives a first prestack seismic dataset and a second prestack seismic dataset wherein the second prestack seismic dataset has a time-lapse relationship to the first prestack seismic dataset; b) a second interface that receives a first source wavelet data set inverted from the first prestack seismic dataset; c) a third interface that receives a first reflectivity parameters dataset inverted from the first prestack seismic dataset; and d) a fourth interface that receives a perturbation dataset determined from a simultaneous minimization of: i) an L1-Norm between a portion of the first reflectivity parameters data set and reflectivity parameters for the second data set; ii) an L2-Norm between the first source wavelet data set and a source wavelet for the second data set; iii) an L2-Norm for data misfit between the first prestack seismic data set and the second prestack seismic data set. ### Claim 18 (depends on 17) The set of application interface programs according to claim 17 further comprising a fifth interface to receive petrophysical parameters inverted from the perturbation dataset. ### Claim 19 (depends on 17) The set of application interface programs according to claim 17 further comprising a sixth interface to receive intercept and gradient parameters determined using the first and second source wavelet data sets and the first and second reflectivity parameters data sets. ### Claim 20 (depends on 17) The set of application interface programs according to claim 17 further comprising a seventh interface to receive an average source wavelet from an inversion of a plurality of common image point gathers from one of the first prestack seismic dataset and the second prestack seismic dataset.