Abstract: Method for estimating geological properties in a subsurface region using multiple types of geophysical data (21). An initial physical properties model 22 is constructed. Some parameters in the model are frozen (23) and optionally portions of the model wave number and spatial domains (24) and the data frequency and data time domains (25), are also frozen. Then, a joint inversion (26) of the multiple data types is performed to calculate an update to the model only for the portions that are not frozen. The converged model (27) for this inversion is used as a new starting model, and the process is repeated (28), possibly several times, unfreezing more parameters and data each time until the desired spatial and parameter resolution (29) has been achieved.

Abstract: Method for adaptive weighting of geophysical data types in iterative joint inversion to speed convergence and aid escape from local minima of the penalty (objective) function. Two or more geophysical data sets (11) representing a region of interest are obtained, and are jointly inverted to infer models of the physical properties that affect the particular types of data used. The misfit for each data type is a weighted tem in the penalty function (13). The invention involves changing the weights (51) as the iteration cycles progress when the iteration convergence criteria are satisfied (15), to see if they remain satisfied (52) with the modified penalty function.

Abstract: Method for joint inversion of geophysical data to obtain 3-D models of geological parameters for subsurface regions of unknown lithology. Two or more data sets of independent geophysical data types are obtained, e.g. seismic and electromagnetic. Then they are jointly inverted, using structural coupling, to infer geophysical parameter volumes, e.g. acoustic velocity and resistivity. Regions of common lithology are next identified based on similar combinations of geophysical parameters. Then a joint inversion of the multiple data types is performed in which rock physics relations vary spatially in accordance with the now-known lithology, and 3-D models of geological properties such as shale content and fracture density are inferred. The computational grid for the last inversion may be defined by the lithology regions, resulting in average geological properties over such regions, which may then be perturbed to determine uncertainty in lithologic boundaries.

Abstract: Method for estimating a model of seismic velocity in a subsurface region from seismic data (31) recorded on 3-component instruments. The method measures the apparent incidence angle (32) of a seismic wave observed at a surface as a function of wave frequency. This apparent angle is converted to an effective velocity as a function of frequency (33), which is then inverted (34) to obtain a subsurface velocity model.

Abstract: Method for joint inversion of geophysical data to obtain 3-D models of geological parameters for subsurface regions of unknown lithology. Two or more data sets of independent geophysical data types are obtained, e.g. seismic and electromagnetic. Then they are jointly inverted, using structural coupling, to infer geophysical parameter volumes, e.g. acoustic velocity and resistivity. Regions of common lithology are next identified based on similar combinations of geophysical parameters. Then a joint inversion of the multiple data types is performed in which rock physics relations vary spatially in accordance with the now-known lithology, and 3-D models of geological properties such as shale content and fracture density are inferred. The computational grid for the last inversion may be defined by the lithology regions, resulting in average geological properties over such regions, which may then be perturbed to determine uncertainty in lithologic boundaries.

Abstract: Method for adaptive weighting of geophysical data types in iterative joint inversion to speed convergence and aid escape from local minima of the penalty (objective) function. Two or more geophysical data sets (11) representing a region of interest are obtained, and are jointly inverted to infer models of the physical properties that affect the particular types of data used. The misfit for each data type is a weighted tem in the penalty function (13). The invention involves changing the weights (51) as the iteration cycles progress when the iteration convergence criteria are satisfied (15), to see if they remain satisfied (52) with the modified penalty function.

Abstract: Method for estimating geological properties in a subsurface region using multiple types of geophysical data (21). An initial physical properties model 22 is constructed. Some parameters in the model are frozen (23) and optionally portions of the model wave number and spatial domains (24) and the data frequency and data time domains (25), are also frozen. Then, a joint inversion (26) of the multiple data types is performed to calculate an update to the model only for the portions that are not frozen. The converged model (27) for this inversion is used as a new starting model, and the process is repeated (28), possibly several times, unfreezing more parameters and data each time until the desired spatial and parameter resolution (29) has been achieved.

Abstract: Method for estimating a model of seismic velocity in a subsurface region from seismic data (31) recorded on 3-component instruments. The method measures the apparent incidence angle (32) of a seismic wave observed at a surface as a function of wave frequency. This apparent angle is converted to an effective velocity as a function of frequency (33), which is then inverted (34) to obtain a subsurface velocity model.