FR 2792419 A1 20001020 FR 9905288 A 19990416 La présente invention concerne de façon générale une méthode pour obtenir un modèle optimal d'une caractéristique physique dans un milieu hétérogène. La méthode selon l'invention permet notamment de modéliser l'impédance d'une zone souterraine relativement à des ondes qui s'y 0o propagent, et d'obtenir par là une image du sous-sol o les discontinuités génératrices de ces variations d'impédance, sont localisées plus précisément. Elle trouve des applications notamment pour faciliter la localisation de gisements d'hydrocarbures. Etat de la technique Il est bien connu à l'heure actuelle que la précision et la fiabilité d'un modèle de réservoir dépend beaucoup d'une intégration de données compatibles mais de nature différentes, plus ou moins reliées aux propriétés d'un gisement étudié. De nombreuses études ont par exemple démontré le rôle que peuvent jouer les données sismiques pour contraindre l'image d'un réservoir dans les espaces compris entre les puits qui y sont forés. Une stratégie courante pour intégrer des données sismiques dans un modèle de réservoir consiste à estimer les impédances acoustiques ou 2 2792419 élastiques à partir des amplitudes sismiques et ensuite à les traduire en propriétés du réservoir. Parmi les algorithmes d'inversion les techniques exploitant un modèle a priori demandent le choix de paramètres qui pondèrent les influences respectives des données et de ce modèle a priori dans l'optimisation. Par les brevets EP 354.112 (US 4,972,383) et EP 889.331 du demandeur, on connaît une méthode dite méthode Interwell pour obtenir un modèle à au moins deux dimensions représentant les variations de l'impédance acoustique d'une zone souterraine par exemple, qui s'accorde au mieux avec des données connues déduites d'observations ou de mesures: données géologiques, données sismiques après sommation, autres données obtenues dans un ou plusieurs puits traversant cette zone, etc. Elle comporte la construction d'un modèle d'impédance a priori en combinant des données de puits et des informations structurales ou stratigraphiques connues, notamment par krigeage des impédances connues le long des puits en suivant les surfaces de corrélation, à l'intérieur des unités stratigraphiques. Les horizons sismiques donnent la géométrie de ces unités stratigraphiques et les corrélations sont faites en suivant les surfaces définies en accord avec le mode de dépôt (concordant, parallèle au toit ou au mur). La distribution initiale d'impédance 2o entre les puits est ensuite mise à jour en tenant compte des amplitudes sismiques relevées. Ceci est réalisé itérativement par la minimisation d'une fonction-coût à deux termes F=Fs+Fg, o: - Fs est un terme sismique défini comme une certaine norme norme représentative de l'écart entre les traces synthétiques obtenues à partir du 3 2792419 modèle et les amplitudes des traces sismiques réelles: traces avant sommation (pre-stack) traces après, sommation (post-stack) etc., et - Fg est un terme géologique mesurant l'écart entre le modèle a priori et le modèle courant. Le modèle d'impédance optimal est un compromis entre les deux termes dépendant de leurs poids respectifs. Ces poids sont en relation avec les paramètres d'inversion c'est-à-dire, pour Fs, avec le rapport du signal au bruit et, pour Fg, en relation avec la longueur de corrélation X et la variance C2 de l'erreur sur le modèle a priori, en formant l'hypothèse que cette erreur est une io variable aléatoire avec une covariance exponentielle stationnaire de variance (2 et de longueur de corrélation A. Habituellement, on estime le poids de Fs à partir des amplitudes réelles. Inversement la longueur de corrélation X et l'écart-type î sont définis par l'utilisateur selon le degré de confiance qu'il peut avoir dans le modèle a priori. Cependant, il ne faut pas que ce dernier point soit évalué de façon trop heuristique car la confiance que l'on a du modèle a priori dépend de la qualité des données, du nombre et de l'emplacement des puits, et aussi du comportement spatial du champ d'impédance. Définition de l'invention La méthode selon l'invention permet d'obtenir, par inversion un modèle optimisé d'une grandeur physique représentative d'un milieu hétérogène stratifié, (l'impédance présentée par ce milieu relativement à des ondes qui s'y propagent, par exemple) par une caractérisation plus quantitative de paramètres d'inversion (écart-type et longueur de corrélation par exemple), qui 4 2792419 soit cohérente avec la structure des données observées ou mesurées dans le milieu. En précisant l'interdépendance de ces paramètres, elle a pour avantage de guider l'interprétateur dans la sélection des valeurs à prendre en compte pour les paramètres. Le modèle qui en résulte est amélioré par rapport à celui que l'on obtient par un choix purement heuristique. Elle comporte la construction d'un modèle a priori à partir de données connues obtenues par des mesures, des enregistrements ou des observations, en un certain nombre de points du milieu, et la construction d'un modèle optimal en minimisant une fonction-coût dépendant des écarts entre le modèle 1o recherché et les données connues, compte-tenu du modèle a priori. La méthode est caractérisée en ce que la construction du modèle a priori comporte une corrélation par krigeage avec une covariance (Cz) entre les valeurs de la grandeur physique connues en différents points du milieu suivant des discontinuités (directions de strates), les incertitudes sur les valeurs de la grandeur physique dans le modèle a priori par rapport aux valeurs correspondantes dans le milieu, en tous points le long de ces directions, respectant un modèle de covariance (CE) qui contrôle les paramètres d'inversion, lequel modèle comporte un terme de covariance stationnaire (Cz) ne dépendant que d'un vecteur d'écart (hi) entre les points et un terme de covariance non stationnaire dépendant de la position des points et de leurs écarts respectifs. Suivant un premier mode de mise en oeuvre, on détermine une covariance moyenne (CE) que l'on ajuste sur un modèle de covariance exponentielle stationnaire (d'un type connu), pour définir les paramètres de l'inversion, de façon plus précise que par la méthode heuristique classique. 2792419 Suivant un deuxième mode de mise en oeuvre, on ajuste en tous points le terme de covariance (CE) sur un modèle exponentiel stationnaire pour obtenir des valeurs locales (3x, o.).des paramètres d'inversion. Figures s D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un mode de réalisation donné à titre d'exemple non limitatif, en se référant aux dessins annexés o: - les Fig.lA, lB représentent deux grilles de covariance de l'erreur de krigeage dans un plan de corrélation autour de trois puits Wl, W2, W3, to correspondant respectivement à deux vecteurs h de distance différents, qui illustrent le caractère non stationnaire de la covariance des incertitudes; - la Fig.2 montre un exemple de modèle d'impédance obtenu avec des paramètres d'inversion usuels; et - la Fig.3 montre un exemple de modèle d'impédance obtenu par la méthode selon l'invention, avec des paramètres d'inversion mieux ajustés au comportement de la covariance de l'incertitude. Description détaillée Comme le modèle a priori est construit par krigeage, on a étudié la covariance théorique des erreurs de krigeage qui doivent en principe commander le poids du terme géologique Fg. Dans le cas d'un krigeage ordinaire, si Z désigne le champ d'impédance, et si Z d priori est le champ d'impédance krigé Z*, alors l'erreur a pour expression E(x) = Z(x) - Z*(x). 6 2792419 Si x et y désignent deux points du milieu distants de h, la covariance C(x,y) de l'erreur e(x) s'exprime alors par: C (x,y) = C. (})-'3:(x) KO(x + I) o K est la matrice de krigeage, /(x) et fi(x + h)sont les pondérations de krigeage respectivement aux abscisses x et x + I, h étant le vecteur de distance de x à y et '/3, le transposé de la pondération 3. La covariance théorique de l'erreur Ce(x,y)entre deux points x et y apparaît comme la covariance initiale C, (h) de krigeage permettant la construction du modèle a priori, covariance qui est le plus souvent isotrope et 1o ne dépend que de la distance ou écart fi, modifiée par un terme non stationnaire '/3(x)K/3(x +h)qui lui, dépend à la fois de la distance hI et de la position du point x considéré. Loin de tout puits, la covariance initiale C, (h)prédomine alors que le second terme prend une importance grandissante au fur et à mesure que l'on se rapproche d'un puits quelconque. Ce comportement est illustré par un exemple sur un champ réel avec trois puits verticaux. Pour l'unité stratigraphique du réservoir, les impédances ont été krigées le long de surfaces parallèles à sa base avec un modèle exponentiel commun (longueur de corrélation 750m, écart-type sur l'impédance égal à 1122g/cm3.m/s). La covariance des erreurs associées est calculée pour deux vecteurs de distance, comme on le voit sur les Fig. 1A, lB. La méthode comporte le choix de paramètres d'inversion compatibles quantitativement avec les données disponibles intégrées dans le processus d'inversion. L'approche la plus simple, qui ne demande aucune modification 7 2792419 de l'algorithme d'inversion, consiste à calculer une covariance moyenne d'erreur CE = {}C, qui sert à ajuster une covariance exponentielle stationnaire x y (d'un type connu des spécialistes) pour obtenir automatiquement les paramètres moyens d'inversion X =À et c6 =5. Une autre approche consiste à réaliser en tout point un ajustement local en calculant le terme correctif '/(x)K3(x + h) approprié pour la position x de ce point, qu'il faut ajouter au terme stationnaire Cz(h), et à déterminer les paramètres d'inversion locaux X =,x, cy = cx, en approchant la covariance locale par une covariance exponentielle stationnaire. La méthode permet donc facilement à l'interprétateur, d'exploiter le comportement spatial de l'erreur de krigeage pour régler le poids à accorder au terme géologique Fg. Dans l'exemple qui a été décrit ci-dessus, la méthode a été appliquée à la modélisation de l'impédance d'une formation relativement à des ondes qui i s'y propagent: ondes élastiques, ondes électromagnétiques etc. Il est bien évident qu'elle peut s'appliquer aussi bien à la modélisation d'autres grandeurs physiques de la formation: porosité, perméabilité, saturation etc., ou relatives à des essais de puits. D'une manière plus générale, la méthode peut s'appliquer à tous les modèles d'inversion de type bayesien bien connu des spécialistes. REVENDICATIONS 1) Méthode pour obtenir par inversion un modèle optimisé d'une grandeur physique représentative d'un milieu hétérogène stratifié, comportant la construction d'un modèle a priori à partir de données connues obtenues par des mesures, des enregistrements ou des observations en un certain nombre de points du milieu, et la construction d'un modèle optimal par inversion et la construction d'un modèle optimal en minimisant une fonction-coût dépendant des écarts entre le modèle recherché et les données connues, compte-tenu du 1o modèle a priori, caractérisée en ce que la construction du modèle a priori comporte une corrélation par krigeage avec une covariance (Cz) entre les valeurs de la grandeur physique connues en différents points du milieu suivant des directions de strates, les incertitudes sur les valeurs de la grandeur physique dans le modèle a priori par rapport aux valeurs correspondantes dans le milieu, en tous points le long de ces directions, respectant un modèle de covariance (CE) qui contrôle des paramètres d'inversion (., c), lequel modèle comporte le terme de covariance stationnaire (Cz) ne dépendant que d'un vecteur d'écart (h) entre les points et un terme de covariance non stationnaire dépendant de la position des points (x) et de leurs distances respectives (h). 2) Méthode selon la revendication 1, caractérisée en ce qu'elle comporte la détermination d'une covariance moyenne (C) que l'on ajuste sur un modèle de covariance exponentielle stationnaire, pour définir les paramètres d'inversion (S, cy). 3) Méthode selon la revendication 1, caractérisée en ce qu'elle comporte en tous points, un ajustement de la covariance (CE) sur un modèle 9 2792419 exponentiel stationnaire pour définir des valeurs locales des paramètres d'inversion (.x, cy) 4) Méthode selon la revendication 1, caractérisée en ce que l'on sélectionne les paramètres d'inversion pour définir les paramètres d'inversion (X,) compatibles avec un modèle de covariance d'erreur défini par: C (x,y) = C (h)-'/3(r) K/3(x+ l) o x et y désignent deux points quelconques du milieu distants de h, K est une matrice de krigeage, 3(x)et (x + Fz) sont des pondérations de krigeage respectivement aux abscisses x et x + h, et h est le vecteur d'écart entre les lo deux points.