Extended subspace method for cross-talk mitigation in multi-parameter inversion
First Claim
1. A computer-implemented method for iteratively inverting measured geophysical data to infer 3D subsurface models of N physical properties with N≥
- 2 and prospecting for hydrocarbons, comprising;
providing an initial model for each physical property, wherein a subsurface region is subdivided into discrete cells, each cell having a value of the physical property;
for each physical property and for each of a plurality of the cells, representing a search direction, indicating whether the initial model needs to be updated positively or negatively, as a linear combination of M basis vectors with M>
N, wherein;
each basis vector has its own coefficient in the linear combination, said coefficient to be determined;
the basis vectors are chosen such that their coefficients account for cross-talk between the N physical properties during inversion, wherein a mixing matrix is formed from the coefficients of the basis vectors, and degree of mixing between the model updates for the N physical properties is adjusted by scaling off-diagonal elements of the mixing matrix with scaling factors that are optimized in the inversion; and
optimal values of the coefficients are simultaneously solved for, using a computer, to minimize or maximize an objective function measuring misfit between model-simulated geophysical data and the measured geophysical data, wherein solving for the coefficients is performed at least by numerical computation constrained by the geophysical data;
generating search directions with the optimal values of the coefficients, and generating an updated model for each physical property by making changes, respectively, to the initial model for each physical property in the search directions; and
producing an image of the subsurface from the updated model, which includes subsurface reflectors, positioned with the optimal coefficients, that returned seismic energy to receivers that recorded the measured geophysical data, and prospecting for hydrocarbons according to structural features of the subsurface region.
0 Assignments
0 Petitions
Accused Products
Abstract
An extended subspace method for inverting geophysical data to infer models for two or more subsurface physical properties, using gradients of an objective function as basis vectors for forming model updates. The extended set of basis vectors provides explicit mixing between gradient components corresponding to different medium parameters, for example P-wave velocity and an anisotropy parameter. In a preferred embodiment, off-diagonal elements of the mixing matrix may be scaled to adjust the degree of mixing between gradient components. Coefficients of the basis vector expansion are determined in a way that explicitly accounts for leakage or crosstalk between different physical parameters. The same extended subspace approach may be used to make further improvement to the model updates by incorporating well constraints, where well log data are available.
-
Citations
18 Claims
-
1. A computer-implemented method for iteratively inverting measured geophysical data to infer 3D subsurface models of N physical properties with N≥
- 2 and prospecting for hydrocarbons, comprising;
providing an initial model for each physical property, wherein a subsurface region is subdivided into discrete cells, each cell having a value of the physical property; for each physical property and for each of a plurality of the cells, representing a search direction, indicating whether the initial model needs to be updated positively or negatively, as a linear combination of M basis vectors with M>
N, wherein;each basis vector has its own coefficient in the linear combination, said coefficient to be determined; the basis vectors are chosen such that their coefficients account for cross-talk between the N physical properties during inversion, wherein a mixing matrix is formed from the coefficients of the basis vectors, and degree of mixing between the model updates for the N physical properties is adjusted by scaling off-diagonal elements of the mixing matrix with scaling factors that are optimized in the inversion; and optimal values of the coefficients are simultaneously solved for, using a computer, to minimize or maximize an objective function measuring misfit between model-simulated geophysical data and the measured geophysical data, wherein solving for the coefficients is performed at least by numerical computation constrained by the geophysical data; generating search directions with the optimal values of the coefficients, and generating an updated model for each physical property by making changes, respectively, to the initial model for each physical property in the search directions; and producing an image of the subsurface from the updated model, which includes subsurface reflectors, positioned with the optimal coefficients, that returned seismic energy to receivers that recorded the measured geophysical data, and prospecting for hydrocarbons according to structural features of the subsurface region. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15)
- 2 and prospecting for hydrocarbons, comprising;
-
16. A computer-implemented method for iteratively inverting measured geophysical data to infer 3D subsurface models of N physical properties with N>
- 2 and prospecting for hydrocarbons, comprising;
providing an initial model for each physical property, wherein a subsurface region is subdivided into discrete cells, each cell having a value of the physical property; for each physical property and for each of a plurality of the cells, representing a search direction, indicating whether the initial model needs to be updated positively or negatively, as a linear combination of a plurality of basis vectors; each basis vector has its own coefficient in the linear combination, wherein;
the basis vectors are chosen such that their coefficients account for cross-talk between the N physical properties during inversion, wherein a mixing matrix is formed from the coefficients of the basis vectors, and degree of mixing between the model updates for the N physical properties is adjusted by scaling off-diagonal elements of the mixing matrix with scaling factors that are optimized in the inversion;determining coefficients of each linear combination by minimizing or maximizing, using a computer, a difference between one or more true parameters computed from well data or other known subsurface information and corresponding parameters predicted by the updated search direction; generating search directions with the determined coefficients, and generating an updated model for each physical property by making changes, respectively, to the initial model for each physical property in the search directions; and producing an image of the subsurface from the updated model, which includes subsurface reflectors, positioned with the optimal coefficients, that returned seismic energy to receivers that recorded the measured geophysical data, and prospecting for hydrocarbons according to structural features of the subsurface region. - View Dependent Claims (17, 18)
- 2 and prospecting for hydrocarbons, comprising;
Specification