INVERSION OF LARGE, NEARLY-HOMOGENEOUS GEOBODIES VIA ULTRA LOW-DIMENSIONAL SHAPE REPRESENTATION USING ORTHOGONAL BASIS FUNCTIONS

A two-stage method for iteratively inverting geophysical data for the purpose of subsurface imaging, including: obtaining at least one geophysical dataset (400) and an initial subsurface model (401); representing subsurface such that a geometry of a geobody is defined using a set of basis functions (403), and a number of such basis functions is significantly smaller than the number of cells used in cell-based geobody representations, wherein an order of magnitude reduction is two or more for 2-D domains and 3 or more for 3-D domains; in a first stage, successively updating the initial subsurface model (413), only for the geobody, by performing iterative low-dimensional geophysical inversion based on minimizing a misfit between simulated geophysical data and the geophysical dataset, wherein the simulated geophysical data is generated from a current subsurface model at each iteration; generating a subsurface image from a final updated subsurface model (415) obtained via the low- dimensional geophysical inversion, wherein the subsurface image includes an inverted geobody; in a second stage, successively updating the subsurface model with the inverted geobody by performing iterative cell-based geophysical inversion (417) based on minimizing a misfit between simulated geophysical data and the geophysical dataset, wherein the simulated geophysical data is generated from a current subsurface model at each iteration. Procédé à deux étapes pour inverser de façon itérative des données géophysiques à des fins d'imagerie souterraine, comprenant : l'obtention d'au moins un ensemble de données géophysiques (400) et d'un modèle souterrain initial (401) ; la représentation d'une subsurface de sorte qu'une géométrie d'un corps géologique soit définie au moyen d'un ensemble de fonctions de base (403), et que le nombre de ces fonctions de base soit significativement plus faible que le nombre de cellules utilisées dans les représentations de corps géologique à base de cellules, dans lequel un ordre de grandeur de réduction est deux ou plus pour les domaines 2-D et 3 ou plus pour les domaines 3-D ; dans une première étape, la mise à jour successive du modèle souterrain initial (413), uniquement pour le corps géologique, en effectuant une inversion géophysique itérative à faible dimension basée sur la réduction au minimum d'une discordance entre les données géophysiques simulées et l'ensemble de données géophysiques, dans lequel les données géophysiques simulées sont générées à