摘要:
Method for reducing the time needed to perform geophysical inversion by using simultaneous encoded sources in the simulation steps of the inversion process. The geophysical survey data are prepared by encoding (3) a group of source gathers (1), using for each gather a different encoding signature selected from a set (2) of non-equivalent encoding signatures. Then, the encoded gathers are summed (4) by summing all traces corresponding to the same receiver from each gather, resulting in a simultaneous encoded gather. (Alternatively, the geophysical data are acquired from simultaneously encoded sources.) The simulation steps needed for inversion are then calculated using a particular assumed velocity (or other physical property) model (5) and simultaneously activated encoded sources using the same encoding scheme used on the measured data. The result is an updated physical properties model (6) that may be further updated (7) by additional iterations.
摘要:
A method for efficient inversion of measured geophysical data from a subsurface region to prospect for hydrocarbons. Gathers of measured data (40) are encoded (60) using a set of non-equivalent encoding functions (30). Then all data records in each encoded gather that correspond to a single receiver are summed (60), repeating for each receiver to generate a simultaneous encoded gather (80). The method employs iterative, local optimization of a cost function to invert the encoded gathers of simultaneous source data. An adjoint method is used to calculate the gradients of the cost function needed for the local optimization process (100). The inverted data yields a physical properties model (110) of the subsurface region that, after iterative updating, can indicate presence of accumulations of hydrocarbons.
摘要:
Method for converting seismic data to obtain a subsurface model of, for example, bulk modulus or density. The gradient of an objective function is computed (103) using the seismic data (101) and a background subsurface medium model (102). The source and receiver illuminations are computed in the background model (104). The seismic resolution volume is computed using the velocities of the background model (105). The gradient is converted into the difference subsurface model parameters (106) using the source and receiver illumination, seismic resolution volume, and the background subsurface model. These same factors may be used to compensate seismic data migrated by reverse time migration, which can then be related to a subsurface bulk modulus model. For iterative inversion, the difference subsurface model parameters (106) are used as preconditioned gradients (107).
摘要:
A method for efficient inversion of measured geophysical data from a subsurface region to prospect for hydrocarbons. Gathers of measured data (40) are encoded (60) using a set of non-equivalent encoding functions (30). Then all data records in each encoded gather that correspond to a single receiver are summed (60), repeating for each receiver to generate a simultaneous encoded gather (80). The method employs iterative, local optimization of a cost function to invert the encoded gathers of simultaneous source data. An adjoint method is used to calculate the gradients of the cost function needed for the local optimization process (100). The inverted data yields a physical properties model (110) of the subsurface region that, after iterative updating, can indicate presence of accumulations of hydrocarbons.
摘要:
Method for converting seismic data to obtain a subsurface model of, for example, bulk modulus or density. The gradient of an objective function is computed (103) using the seismic data (101) and a background subsurface medium model (102). The source and receiver illuminations are computed in the background model (104). The seismic resolution volume is computed using the velocities of the background model (105). The gradient is converted into the difference subsurface model parameters (106) using the source and receiver illumination, seismic resolution volume, and the background subsurface model. These same factors may be used to compensate seismic data migrated by reverse time migration, which can then be related to a subsurface bulk modulus model. For iterative inversion, the difference subsurface model parameters (106) are used as preconditioned gradients (107).
摘要:
Method for simultaneous full-wavefield inversion of gathers of source (or receiver) encoded geophysical data to determine a physical properties model for a subsurface region, especially suitable for surveys where fixed receiver geometry conditions were not satisfied in the data acquisition. First, a shallow time window of the data (202) where the fixed receiver condition is satisfied is inverted by simultaneous encoded (203) source inversion (205). Then, the deeper time window of the data (208) is inverted by sparse sequential source inversion (209), using the physical properties model from the shallow time window (206) as a starting model (207). Alternatively, the shallow time window model is used to simulate missing far offset data (211) producing a data set satisfying the stationary receiver assumption, after which this data set is source encoded (212) and inverted by simultaneous source inversion (214).
摘要:
Method for reducing the time needed to perform geophysical inversion by using simultaneous encoded sources in the simulation steps of the inversion process. The geophysical survey data are prepared by encoding (3) a group of source gathers (1), using for each gather a different encoding signature selected from a set (2) of non-equivalent encoding signatures. Then, the encoded gathers are summed (4) by summing all traces corresponding to the same receiver from each gather, resulting in a simultaneous encoded gather. (Alternatively, the geophysical data are acquired from simultaneously encoded sources.) The simulation steps needed for inversion are then calculated using a particular assumed velocity (or other physical property) model (5) and simultaneously activated encoded sources using the same encoding scheme used on the measured data. The result is an updated physical properties model (6) that may be further updated (7) by additional iterations.
摘要:
A method, including: storing, in a computer storage device, geophysical seismic data that has been separated into a multiple-free component and a multiple contaminated component; performing, with a processor, a first full wavefield inversion process on the multiple-free component of the seismic data, wherein a first subsurface physical property model is generated; determining, with a processor, an extended target reflectivity, wherein the extended target reflectivity includes a reflectivity for each of a plurality of shots; separately performing, with a processor, a second full wavefield inversion process with the multiple contaminated component of the seismic data for each of the plurality of shots using the reflectivity corresponding to each of the plurality of shots, wherein a second subsurface physical property model is generated; and generating, with a processor, multiple-free final subsurface physical property model by combining the first subsurface physical property model and the second subsurface physical property model.
摘要:
An improved method for reducing the accuracy requirements on the starting model when performing multi-scale inversion of seismic data (65) by local objective function optimization (64). The different scales of inversion are brought about by incorporating a low-pass filter into the objective function (61), and then decreasing the amount of high-frequency data that is filtered out from one scale to the next. Moreover, the filter is designed to be time varying, wherein the filter's low-pass cutoff frequency decreases with increasing traveltime of the seismic data being filtered (62). The filter may be designed using Pratt's criterion for eliminating local minima, and performing averages (or other statistical measure) of the period and the traveltime error only with respect to source and receiver location but not traveltime (63).
摘要:
Method for reducing instability and increasing computational efficiency in tomographic inversion for velocity model building. A system of tomographic equations is developed for a uniform grid. A non-uniform parameterization is found for which a linear mapping exists between the space of the uniform grid and the space of the non-uniform grid. The matrix that relates velocity to the tomographic data in the non-uniform representation is then given by the matrix product of the corresponding matrix in the uniform grid representation and the mapping matrix. Inversion can then be performed for the non-uniform parameterization on a smaller, more stable matrix.