JOINT INVERSION WITH UNKNOWN LITHOLOGY

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. La présente invention concerne un procédé d'inversion conjointe de données géophysiques qui permet d'obtenir des modèles tridimensionnels de paramètres géophysiques pour des régions du sous-sol ayant une lithologie inconnue. Au moins deux ensembles de données ayant des types de données géophysiques indépendants sont obtenus, par exemple des données sismiques et électromagnétiques. Celles-ci sont ensuite conjointement inversées au moyen d'un couplage structurel, afin d'en déduire des volumes de paramètres géophysiques, par exemple la vitesse acoustique et la résistivité. Des régions présentant une lithologie commune sont ensuite identifiées sur la base de combinaisons semblables de paramètres géophysiques. Une inversion conjointe des multiples types de données est ensuite effectuée, lors de laquelle des relations concernant la physique des roches varient spatialement en fonction de la lithologie maintenant connue, et des modèles tridimensionnels de propriétés géologiques telles que la teneur en schistes et la densité de fracture en sont déduits. La grille informatique utilisée pour la dernière inversion peut être définie par les régions lithologiques, ce qui conduit à des propriétés géologiques moyennes sur ces régions, qui peuvent ensuite être perturbées pour déterminer l'incertitude affectant les frontières lithologiques.