FR2938952A1 - Estimation de proprietes lithologiques d'une zone geologique - Google Patents

Estimation de proprietes lithologiques d'une zone geologique Download PDF

Info

Publication number
FR2938952A1
FR2938952A1 FR0858021A FR0858021A FR2938952A1 FR 2938952 A1 FR2938952 A1 FR 2938952A1 FR 0858021 A FR0858021 A FR 0858021A FR 0858021 A FR0858021 A FR 0858021A FR 2938952 A1 FR2938952 A1 FR 2938952A1
Authority
FR
France
Prior art keywords
node
nodes
grid
state
target node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
FR0858021A
Other languages
English (en)
Inventor
Pierre Biver
Denis Allard
Or Dimitri D
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TotalEnergies SE
Original Assignee
Total SE
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Total SE filed Critical Total SE
Priority to FR0858021A priority Critical patent/FR2938952A1/fr
Priority to US13/131,312 priority patent/US8630832B2/en
Priority to PCT/FR2009/052301 priority patent/WO2010061135A1/fr
Publication of FR2938952A1 publication Critical patent/FR2938952A1/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling
    • G01V2210/665Subsurface modeling using geostatistical modeling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

Un procédé d'estimation de propriétés lithologiques d'une zone géologique représentée par une grille, dans lequel on cherche à assigner aux noeuds de la grille des faciès, par simulations stochastiques séquentielles. Pour un noeud cible (x ) parmi un ensemble de noeuds cibles à parcourir par itérations, on estime pour chaque faciès possible (C ) une probabilité conditionnelle ( p , , , ) d'observer ce faciès, sachant que des faciès respectifs ont déjà été assignés à d'autres noeuds du voisinage. Un tirage au sort pondéré par ces probabilités conditionnelles est ensuite effectué. L'estimation des probabilités conditionnelles fait intervenir des probabilités univariées (p ), et des probabilités bivariées (p ) d'observer deux faciès pour deux noeuds respectifs de la grille. Chaque probabilité bivariée intervenant dans l'estimation de la probabilité conditionnelle se rapporte au noeud cible (x ), et à un noeud (X ) voisin du noeud cible.

Description

ESTIMATION DE PROPRIETES LITHOLOGIQUES D'UNE ZONE GEOLOGIQUE L'invention se rapporte au domaine de l'estimation des propriétés lithologiques d'une zone géologique pour l'étude du sous-sol. II est connu de simuler une zone géologique pour estimer la structure géométrique d'un réservoir ainsi que la nature des roches qui composent le réservoir (grès, carbonates, schistes...), c'est à dire la lithologie du réservoir. On peut ainsi déterminer quelles sont les régions de la zone géologique qui ont le plus de chance de contenir des hydrocarbures, les régions dans lesquelles il peut être intéressant de forer un puits d'injection pour améliorer la récupération des hydrocarbures, le type de fluides récupérés, ou autre. II s'agit in fine d'obtenir une représentation du réservoir la plus précise possible afin de déterminer au mieux les paramètres techniques relatifs à la recherche, à l'étude ou l'exploitation de la zone. Il est crucial, dans le domaine pétrolier, que la simulation soit aussi précise que possible. En général, le géologue dispose seulement d'un petit nombre de données expérimentales (obtenues par exemple par carottage en certains points du réservoir). La simulation consiste à construire l'image d'un réservoir, en interpolant ces données expérimentales connues. Cependant, du fait du petit nombre d'informations disponibles, il est nécessaire de recourir à des procédés de construction probabilistes.
Généralement la simulation est réalisée par ordinateur, et la structure géométrique du réservoir est représentée sous la forme d'une grille. De façon conventionnelle, chacun des noeuds de la grille est identifié par ses coordonnées dans l'espace, selon les 3 axes. Par exemple, les coordonnées (xo, yo, zo) définissent un noeud, lequel est noté xo dans la présente description à des fins de clarté. On appelle variable catégorielle une fonction susceptible de prendre un nombre fini d'états. Par exemple, dans le cas d'un réservoir, si la variable catégorielle est la lithologie, les différents états sont les types de faciès susceptibles d'être observés dans le réservoir, par exemple des schistes, des argiles, des grès... Pour simuler la lithologie du réservoir, on assigne à chaque noeud de la grille un état, parmi l'ensemble des états susceptibles d'être pris par la variable catégorielle Par exemple, si 3 faciès (schistes, argiles, grès) différents peuvent être observés dans le réservoir, on aura 3 états susceptibles d'être pris par la variable catégorielle qui sont l'état c, qui correspond à des schistes, l'état c2 qui correspond à des argiles, et l'état c3 qui correspond à des grès. Il y aura en conséquence 3 assignations possibles (c,, c2, et c3) pour chaque noeud de la grille. Dans les méthodes de simulation stochastique, l'attribution d'un état ck en un noeud x; de la grille se fait par tirage au sort pondéré par les probabilités d'occurrence de l'état Ck. Lorsque les états sont des faciès on parle de façon conventionnelle de simulation de faciès.
Dans les méthodes de simulation stochastiques séquentielles, on simule les faciès en un noeud donné de la grille en fonction des faciès déjà assignés en d'autres noeuds de la grille. Au final, on parcourra séquentiellement l'ensemble de la grille. On parle de simulation séquentielle, car on utilise les noeuds précédemment simulés pour simuler les nouveaux noeuds de la grille. On appelle probabilité conditionnelle au noeud xo la probabilité d'observer un état d'une variable catégorielle, par exemple un faciès, en un noeud xo de la grille, connaissant l'état de cette variable catégorielle en d'autres noeuds de la grille, typiquement des noeuds voisins du noeud xo. Les méthodes de simulation stochastique séquentielles consistent à déterminer ces probabilités conditionnelles puis à simuler l'état d'une variable catégorielle, et ce séquentiellement pour chaque noeud de la grille. L'établissement des probabilités conditionnelles est donc une étape essentielle du procédé de simulation.
Il existe plusieurs méthodes connues pour déterminer les probabilités conditionnelles. Par exemple, dans la méthode par statistique multipoints les probabilités conditionnelles sont évaluées par dénombrement d'évènements pour des configurations données de points dans une image d'entraînement. Cette méthode est notamment décrite dans la demande WO2006023597. La méthode de simulation par champs aléatoires multinomiaux, développée par P.Bogaert (voir par exemple " Spatial Prediction of categorical variables, the Bayesian maximum entropy approach ", dans la revue : Stochastic Environmental Research and Risk Assessment, Volume 16, Numéro 6 / décembre 2002, pages 425-448, Springer Berlin / Heidelberg, DOI : 10.1007/s00477-002-0114-4), fait intervenir une étape de calcul de probabilités conjointes. Une probabilité conjointe est une probabilité d'obtenir un événement donné, l'évènement consistant en une configuration d'états, c'est-à-dire une disposition d'états d'une variable catégorielle localisés dans l'espace. Ainsi, dans le contexte de la simulation d'une zone géologique, une probabilité conjointe est la probabilité d'observer simultanément N états d'une variable catégorielle cko, ck;,..., ck(N_1) en N noeuds respectifs xo, x1,..., XN_1 de la grille. Si on ne considère qu'un couple de noeuds, on parlera de probabilité bivariée plutôt que de probabilité conjointe. A partir de ces probabilités conjointes, on peut déterminer des probabilités conditionnelles en utilisant le théorème de Bayes. Pour établir les tables de probabilités conjointes, on utilise classiquement la méthode de la maximisation de la vraisemblance. Dans l'algorithme proposé par Bogaert, la maximisation de la vraisemblance est réalisée en utilisant l'algorithme IPF (pour " Iterative Proportional Fitting voir par exemple W.E. Deming et F.F. Stephan : " On a least square adjustment of a sampled frequency table when the expected marginal totals are known ", Annals of Mathematical Statistics, vol. 11, page 427, 1940), mais en l'appliquant au contexte géostatistique. Cette méthode fait intervenir des probabilités bivariées d'observer simultanément 2 états d'une variable catégorielle ck, ck' à 2 noeuds respectifs x;, x;, de la grille. Il est ainsi nécessaire de calculer des probabilités bivariées pour chaque couple ck, ck' d'états de la variable catégorielle et pour chaque couple de noeuds x;, x;, parmi les noeuds du voisinage considéré. Cette étape est relativement coûteuse en temps de calcul et en mémoire et devient en pratique très difficile à réaliser, notamment lorsque l'on utilise plus de vingt noeuds voisins et plus de trois états de la variable catégorielle. L'invention vise à résoudre ce problème, en proposant d'estimer les probabilités conditionnelles à partir d'un nombre réduit de probabilités bivariées. Plus précisément, lors de l'estimation d'une probabilité conditionnelle au noeud xo, chaque probabilité bivariée intervenant dans cette estimation se rapporte au noeud xo et à un noeud x; voisin du noeud xo, auquel a déjà été assigné un état de la variable catégorielle. On élimine ainsi les termes croisés correspondant à des probabilités bivariées entre deux noeuds du voisinage du noeud xo. Selon un premier aspect, l'invention a pour objet un procédé d'estimation de propriétés lithologiques d'une zone géologique simulée par ordinateur ou autres moyens informatiques, comprenant a/ représenter la zone géologique à l'aide d'une grille comportant des noeuds, b/ représenter une propriété lithologique à l'aide d'une variable catégorielle, pour chaque noeud cible parmi un ensemble de noeuds cibles de la grille à parcourir itérativement : c/ pour chaque état possible de la variable catégorielle, estimer la valeur d'une probabilité conditionnelle d'observer cet état de la variable catégorielle en ce noeud cible sachant qu'un état respectif a été assigné à au moins un noeud voisin du noeud cible, à partir : de valeurs de probabilité univariée d'observer un état de la variable catégorielle pour le noeud cible, et de probabilités bivariées d'observer deux états de la variable catégorielle pour deux noeuds respectifs de la grille, dans lequel, les valeurs de probabilités bivariées intervenant dans l'estimation de la valeur de la probabilité conditionnelle se rapportent au noeud cible et à un noeud voisin du noeud cible, d/ effectuer un tirage au sort parmi états possibles de la variable catégorielle, pondéré par les valeurs de probabilités conditionnelles estimées à l'étape cl, et e/ assigner au noeud cible l'état tiré au sort à l'étape d/.
Les étapes cl, d/, e/, sont répétées itérativement, de façon à parcourir l'ensemble de noeuds à simuler, typiquement tous les noeuds de la grille auxquels aucun état n'est encore assigné. Un tel procédé requiert moins de ressources mémoire et/ou de calcul que le procédé de simulation de l'art antérieur, du fait du nombre réduit de probabilités bivariées et de l'absence de probabilités conjointes d'ordre supérieur à 2. Ce procédé, relativement simple à mettre en oeuvre, permet d'obtenir des simulations relativement fiables.
En outre, ce procédé peut permettre de calculer directement les valeurs des probabilités conditionnelles, sans calcul ni stockage préalable des tables de probabilités conjointes. La propriété lithologique représentée par la variable catégorielle peut être un faciès, un type de roche, ou un type d'environnement géologique.
Par noeuds voisins du noeud cible, on entend les noeuds de la grille suffisamment proches du noeud cible pour que les états de la variable catégorielle éventuellement assignés à ces noeuds aient une influence non négligeable sur l'état de la variable catégorielle au noeud cible. On pourra prévoir de sélectionner les noeuds à considérer pour l'étape c/ en définissant un domaine géométrique autour du noeud cible et en considérant les noeuds de ce domaine auxquels un état a déjà été assigné. Avantageusement, préalablement à l'étape d'estimation de la probabilité conditionnelle, on prévoit un modèle de probabilités bivariées donnant, pour chaque couple d'états de la variable catégorielle, la probabilité d'observer deux états de la variable catégorielle pour deux noeuds respectifs de la grille, selon les positions relatives desdits noeuds l'un par rapport à l'autre. Lors de l'étape c/ d'estimation de la probabilité conditionnelle, pour chaque état de variable catégorielle possible, et pour chaque couple de noeuds comprenant le noeud cible et un noeud du voisinage auquel a déjà été assigné un état de la variable catégorielle, on peut extraire de ce modèle la valeur de probabilité bivariée d'observer deux états de la variable catégorielle pour un couple de noeuds ayant les mêmes positions relatives l'un par rapport à l'autre que le noeud cible et le noeud du voisinage. Cette utilisation d'un modèle se base sur une hypothèse d'invariance des probabilités bivariées selon la position des noeuds cible et voisin dans la grille, le modèle étant fonction des positions relatives desdits noeuds l'un par rapport à l'autre. Le recours à cette hypothèse de stationnarité de second ordre permet d'obtenir les probabilités bivariées relativement facilement. Le modèle de probabilités bivariées est un ensemble de fonctions, qui donnent la probabilité d'observer simultanément un état d'une variable catégorielle, par exemple un faciès noté Ck en un noeud x; de la grille, et une autre état de la variable catégorielle, par exemple un faciès noté Ck' en un autre noeud x;, de la grille en fonction des positions relatives des noeuds x;, x;, l'un par rapport à l'autre, par exemple en fonction de la distance h qui les sépare, ck et Ck' pouvant être identiques ou non .
Ces fonctions sont fournies pour chaque couple (ck,ck') d'états de la variable catégorielle. Ainsi, si on considère un réservoir avec 3 faciès possibles notés c1, c2, c3, est fourni un modèle de probabilités bivariées comprenant 9 fonctions analytiques correspondant aux couples (c1, c1), (ci, C2) (Cl, C3), (C2, Cl), (C2, C2), (C2, C3), (C3, Cl), (C3, C2), (C3, C3).
Dans un premier mode de réalisation, le modèle de probabilités bivariées est stationnaire et isotrope, c'est-à-dire à dire que la probabilité d'observer deux états donnés Ck, Ck' en deux noeuds respectifs x; , x;, ne dépendra que de la distance h séparant les noeuds x; et x;,. La fonction correspondant à un couple (Ck, ck') est identique à la fonction du couple (ck', ck), de sorte que le modèle est relativement simple et occupe relativement peu de mémoire, ce qui peut être intéressant dans le cas d'un nombre d'états possibles élevé. Dans un deuxième mode de réalisation, le modèle de probabilités bivariées est stationnaire et anisotrope, c'est-à-dire à dire que la probabilité d'observer deux états donnés ck, Ck, en deux noeuds respectifs x;, x;, dépendra de la distance h séparant les noeuds x; et x;, mais aussi de l'orientation respective des noeuds x; et x;,. Dans un tel modèle, un ensemble de fonctions est fournie pour chaque couple (ck,ck') d'états de la variable catégorielle, chacune des fonctions correspondant à une direction dans l'espace. Un tel modèle permet de traduire l'asymétrie des dépôts. Le modèle de probabilités bivariées peut être obtenu : - à partir de données de puits : l'ensemble des états de la variable catégorielle par exemple les faciès, observés en différents puits sont recensés. Pour un couple d'états donnés, par exemple (c,, c2) avec ci correspondant à des argiles et c2 à des grès, observés en 2 points x; , x;, séparés d'une distance h, on recense le nombre d'occurrences de ce couple d'états pour la distance h choisie. Connaissant le nombre total de couples distants de h dans le jeu de données, on peut calculer les probabilités d'occurrence de couple de faciès (ci, c2) observés en deux points x;, x;, distants de h. En répétant cette opération pour des distances h croissantes, on obtient un nombre d'évènements expérimentaux. A partir de ces dénombrements on ajuste une fonction. On construit ainsi le modèle de probabilités bivariées pour chaque couple d'états. Dans le mode de réalisation pour lequel le modèle de probabilités bivariées est anisotrope on réalisera ces dénombrements pour plusieurs directions. On obtiendra in fine une fonction pour chaque couple d'états et pour chaque direction considérée. - À partir d'une image de référence : l'image de référence correspond à une représentation plausible du réservoir, fournie par le géologue. La méthode est la même, mais le couple de faciès (c1,c2) n'est pas recensé à partir des données observées au puits, mais à partir de l'image de référence. Un exemple d'image de référence (dunes) est donné figure 4. On peut prévoir une étape d'initialisation au cours de laquelle on assigne un état de la variable catégorielle à au moins un noeud de la grille. Ainsi, on peut utiliser cette connaissance des états d'un ou plusieurs noeuds, lors de la simulation stochastique d'un état pour le premier noeud cible de l'ensemble de noeuds cibles à parcourir itérativement. Lors de cette étape d'initialisation, on assigne à un ou plusieurs noeuds un ou des états respectifs. Ces états peuvent être par exemple des faciès observés à partir de la zone géologique réelle, et/ou des états obtenus par simulation stochastique à partir des probabilités univariées extraites du modèle de probabilités bivariées.
Par exemple, des états de la variable catégorielle correspondant à des faciès peuvent être obtenus par observation directe, par carottage par exemple, ou bien suite à l'interprétation de mesures sismiques. Ces états sont assignés aux noeuds correspondants de la grille lors de l'étape d'initialisation. Selon un autre exemple, si le géologue ne dispose d'aucun état observé, on assignera un état à un noeud de la grille, par simulation stochastique à partir des probabilités univariées extraites du modèle de probabilités bivariées. En effet, les valeurs de probabilités univariées sont facilement extraites du modèle de probabilités bivariées en considérant uniquement les couples de même état, par exemple (ci, ci) , (c2, c2), (c3, c3), et en considérant une distance nulle entre deux noeuds de la grille. Les étapes cl, dl, e/ de simulation stochastique sont répétées itérativement et permettent d'assigner des états à des noeuds cibles de la grille. L'invention n'est en rien limitée par l'ordre selon lequel les noeuds de la grille sont parcourus. Par exemple, on extrait M sous-grilles de la grille définie à l'étape a/, chaque sous-grille étant composés de noeuds deux fois moins espacés entre eux que les noeuds de la sous-grille précédente, la sous-grille de niveau M correspondant à un espacement maximal, et la sous-grille de niveau 0 correspondant à la grille définie à l'étape a/. Les étapes c/, dl, e/ sont d'abord appliquées aux noeuds de la Mième sous-grille, de façon à assigner des faciès à chaque noeud de cette Mième sous-grille. Puis, pour m courant de M-1 à 0, les étapes c/, dl, e/ sont appliquées aux noeuds de la mième sous-grille, en utilisant les états assignés aux noeuds de la sous-grille de niveau m+1. Selon un autre aspect, l'invention a pour objet un produit programme d'ordinateur destiné à être stocké dans une mémoire d'une unité centrale, et/ou stocké sur un support mémoire destiné à coopérer avec un lecteur de ladite unité centrale et/ou téléchargé via un réseau de télécommunication, caractérisé en ce qu'il comprend des instructions pour exécuter les étapes du procédé exposé ci-dessus. Selon encore un autre aspect, l'invention a pour objet un dispositif d'estimation des propriétés lithologiques d'une zone géologique, comprenant une unité de mémoire pour stocker une grille comprenant des coordonnées spatiales de noeuds, et pour stocker un nombre fini de valeurs qui correspondent à des états susceptibles d'être pris par une variable catégorielle représentant une propriété lithologique. Le dispositif comporte en outre des moyens de traitement, agencés pour itérativement, chaque itération correspondant à un noeud cible de la grille: • pour chaque état possible de la variable catégorielle, estimer une valeur de probabilité conditionnelle d'observer cet état pour le noeud cible sachant qu'un état de la variable catégorielle respectif a été affecté à au moins un noeud voisin du noeud cible, à partir : de probabilités univariées d'observer un état de la variable catégorielle pour le noeud cible, et de probabilités bivariées d'observer deux états de la variable catégorielle pour deux noeuds respectifs de la grille, dans lequel, les probabilités bivariées intervenant dans l'estimation de la probabilité conditionnelle se rapportent au noeud cible et à un noeud voisin du noeud cible, • effectuer un tirage au sort parmi les états possibles, pondéré par les valeurs de probabilités conditionnelles ainsi estimées, et assigner au noeud cible l'état ainsi tiré au sort. Le dispositif peut comprendre par exemple un ordinateur, un calculateur dédié à la simulation stochastique de zones géologiques, ou autre. L'unité de mémoire peut comprendre une ou plusieurs mémoires. Les 25 moyens de traitement peuvent par exemple comprendre un ou plusieurs processeurs. D'autres particularités et avantages de la présente invention apparaîtront dans la description détaillée ci-après, faite en référence aux dessins annexés sur lesquels : 30 - Les figures 1A et 1B montrent un exemple d'algorithme d'un procédé selon un mode de réalisation de l'invention, - Les figures 2A et 2B illustrent des procédés d'estimation de probabilités conditionnelles, respectivement suivant l'art antérieur et suivant un mode de réalisation de l'invention, - Les figures 3A à 31 montrent un exemple de modèle de probabilités bivariées, - La figure 4 montre un exemple d'image de référence utilisable pour construire un modèle de probabilités bivariées, - Les figures 5A à 5J montrent des exemples de grilles. - Les figures 6A à 6C illustrent un exemple de procédé mettant en oeuvre une approche multigrilles, selon un mode de réalisation de l'invention, - La figure 7 montre un exemple de dispositif selon un mode de réalisation de l'invention. A moins qu'il n'en soit précisé autrement, dans la description détaillée ci-dessous, la grille est bidimensionnelle. On comprendra bien que la grille est avantageusement tridimensionnelle, et que le choix d'une grille bidimensionnelle dans la description a été effectué pour faciliter la compréhension. En référence à la figure 1A, est représenté un algorithme d'un procédé d'aide à la prospection. Ce procédé permet de simuler un milieu hétérogène poreux, par exemple un réservoir d'hydrocarbures. Une représentation géométrique du réservoir est fournie lors d'une étape a/, sous la forme d'une grille. Sont conservées en mémoire les coordonnées spatiales des noeuds de la grille.
Lors d'une étape b/, on mémorise un nombre fini d'états d'une variable représentant une propriété lithologique Dans cet exemple, chaque état possible de la variable catégorielle correspond à un faciès de roche particulier. Par exemple, on mémorise trois états : l'état c, qui correspond à des schistes, l'état c2 à des argiles, et l'état c3 qui correspond à des grès. Par souci de clarté, pour la suite de la description, on appellera faciès c, l'état faciès c2 l'état c2, faciès c3 l'état c3. On cherche à assigner un faciès à chaque noeud de la grille, ce qui revient à établir la distribution des faciès sur le support géométrique.
Pour ce faire, un modèle de probabilités bivariées est défini lors d'une étape 1. Ce modèle peut par exemple être reçu d'un autre dispositif, lu dans une mémoire, ou bien construit à partir de données de puits ou d'une image de référence.
Un exemple d'image de référence est montré par la figure 4. Il s'agit d'une représentation plausible de la zone géologique, fournie par le géologue. Pour construire le modèle, l'ensemble des faciès relevés sur l'image de référence est recensé. Pour un couple de faciès donnés, par exemple (cl, c2), observés en 2 points séparés d'une longueur h, on recense le nombre d'occurrence de ce couple de faciès en fonction de la distance h séparant les 2 points. Connaissant le nombre total de couples de faciès donnés séparés d'une longueur h, on peut calculer les probabilités d'occurrence d'un couple de faciès (c,, c2) observés en 2 points xi, x2 en fonction de la distance séparant xi et x2. A partir de ces données calculées on ajuste une fonction. On construit ainsi le modèle de probabilités bivariées pour chaque couple de faciès. La figure 3 montre un exemple de modèle de probabilités bivariées, pour 3 faciès ck possibles, soit 9 couples (k,k'). Chacun des graphes représentés correspond à un couple de faciès. En abscisse figure une distance entre deux noeuds, normalisée par le pas de la grille. En ordonnée figure la probabilité que deux noeuds aient les faciès du couple correspondant au graphe en question. Le modèle de la figure 3 est stationnaire et anisotrope, c'est-à-dire que la probabilité d'observer deux états donnés ck, ck' en deux noeuds respectifs x; , x;, dépendra de la distance h séparant les noeuds x; et x;, mais aussi de l'orientation respective des noeuds x; et xe. Aussi sur chaque graphe figurent plusieurs courbes, chaque courbe correspondant à une orientation. La courbe pleine donne les probabilités pour un deuxième noeud x;, sur la même ligne horizontale de la grille et à droite du premier noeud x; (orientation à 0°), et la courbe en pointillés donne les probabilités pour un deuxième noeud x;, sur la même ligne verticale que le premier noeud et plus haut, c'est à dire plus près de la surface du sol, que ce premier noeud x; (orientation à 90°).
On peut facilement extraire de ce modèle les probabilités univariées : il s'agit des probabilités pour les couples (1,1), (2,2) et (3, 3) pour une distance entre noeuds nulle. Dans cet exemple, les probabilités univariées pour les différents faciès sont équiprobables, avec des valeurs égales à 1 / 3.
On notera qu'avec cet exemple de modèle, pour un deuxième noeud adjacent et au dessus d'un premier noeud, il est plus facile de passer du faciès c, à c2, de c2 à c3 et de c3 à ci, (les probabilités pour les couples (1,2), (2,3) et (3,1) étant autour de 0,12 pour une distance de 1), que de passer directement du faciès ci vers c3 (la probabilité pour le couple (1,3) étant autour de 0,04 pour une distance de 1). Il s'agit donc d'un modèle avec une asymétrie Nord-Sud. Pour revenir à la figure 1, au cours d'une étape d'initialisation 6, un ou plusieurs noeuds de la grille se voient affecter un faciès. Ce faciès peut être obtenu expérimentalement à partir de la zone géologique réelle ou à partir d'une simulation stochastique. Dans cet exemple, si on dispose d'observations du réservoir, certains noeuds de la grille sont renseignés en utilisant ces observations. Si l'on ne dispose d'aucune donnée observée, on assigne à un noeud un faciès par simulation stochastique .
Pour la première alternative, les observations peuvent comprendre les faciès observés expérimentalement en certains points du réservoir. Il peut s'agir d'une observation directe, par carottage par exemple, ou bien de l'interprétation d'une campagne sismique. Dans la deuxième alternative, le faciès est obtenu par tirage au sort, pondéré par des valeurs de probabilités univariées extraites du modèle défini à l'étape 1. Pour chaque noeud que l'on souhaite ainsi renseigner, on effectue un tirage au sort parmi les faciès ci, c2, c3. Avec l'exemple de modèle de la figure 3, ces probabilités univariées pl, p2, P3 sont égales, de sorte que les résultats du tirage au sort seraient ceux d'un tirage au sort non pondéré. Une fois ces étapes initiales (a), (b), 1 et 6 effectuées, on assigne à chaque noeud de la grille auquel aucune valeur n'est assignée, dit noeud cible, un faciès obtenu par simulation stochastique (étapes (c), (d), (e)). Pour chaque noeud cible, un tirage au sort est effectué à l'étape (d) parmi les faciès possibles ci, c2, c3, ce tirage au sort étant pondéré par des valeurs de probabilités estimées à l'étape (c). L'étape (e) consiste à assigner au noeud cible la valeur ainsi tirée au sort lors de l'étape (d), de sorte que le noeud cible devient ainsi un noeud renseigné dont le faciès est susceptible d'être utilisé lors de la simulation stochastique des faciès pour d'autres noeuds. Les étapes (c), (d), (e), peuvent ainsi être effectuées plusieurs fois, de façon à assigner des faciès à plusieurs noeuds respectifs, typiquement chaque noeud auquel aucun faciès n'est encore assigné. Des étapes de test et de changement de noeud cible, représentées schématiquement par les références 4 et 5, peuvent ainsi être mises en place de façon à parcourir les noeuds non renseignés, par exemple selon des procédés bien connus de l'homme du métier. La figure 1B montre de façon plus détaillée l'étape (c). Au cours de cette étape, pour un noeud cible donné, et pour chaque faciès possible ci, c2, c3, on estime une probabilité Pk0 k, kn . Il s'agit d'une probabilité conditionnelle Pko d'observer le faciès Cko au noeud cible xo sachant que certains noeuds x;,..., xn du voisinage du noeud cible xo se sont déjà vu affecter certains faciès respectifs ckl,..., Ckn.
On peut prévoir une étape non représentée au cours de laquelle on sélectionne les noeuds du voisinage du noeud cible auxquels des faciès ont déjà été assignés, pour être utilisés lors de l'étape (c). Par exemple, on définit une ou des distances, de façon à définir un domaine géométrique centré sur le noeud xo. Par exemple, trois distances vont définir un ellipsoïde.
Dans ce mode de réalisation, on recensera uniquement les noeuds pour lesquels un faciès a déjà été assigné, et qui se situent dans le domaine géométrique précédemment défini. Cette étape de sélection permet de limiter les temps de calcul lors de l'établissement des probabilités conditionnelles, puisqu'elle permet de ne considérer qu'un petit nombre de probabilités bivariées. Il s'agit in fine de réaliser une approximation lors du calcul des probabilités conditionnelles en ne considérant que les noeuds situés dans le voisinage immédiat de celui pour lequel on souhaite attribuer un faciès. Bien entendu, on pourrait alternativement prévoir d'utiliser tous les noeuds de la grille auxquels un faciès a déjà été assigné (on parle de noeuds renseignés) lors de l'estimation de la probabilité conditionnelle. On obtiendrait ainsi une estimation plus exacte de la probabilité conditionnelle, mais au prix de calculs plus longs. Lors d'une étape (cl), on extrait du modèle de probabilités bivariées des valeurs de probabilités bivariées. Pour chaque couple de noeuds ()c0, x;), xo étant le noeud cible et x; étant un noeud sélectionné, on extrait à partir des modèles de probabilités bivariées une valeur de probabilités bivariées Pk0,k, pour chaque faciès possible cko. On établit ainsi une table de probabilités bivariées indiquant la probabilité d'observer à la fois un faciès cko, au noeud xo et un faciès ck;, au noeud x; connaissant la position relative entre les noeuds xo et x;. Les faciès cko et ck; peuvent être identiques ou différents. Connaissant les positions relatives des deux noeuds xo et x;, la table de probabilités bivariées pour ce couple de noeuds est facilement extraite du modèle de probabilités bivariées fourni à l'étape 1. En effet, les fonctions analytiques donnent la probabilité d'observer simultanément un faciès ck en un noeud x; et un faciès ck, en un noeud x; (k et k' pouvant être identiques) en fonction des positions relatives des noeuds x; et xi, considérés. Par exemple, pour deux noeuds x(;=1), x(;'_2) déjà renseignés au voisinage du noeud cible xo ; le faciès c1 étant par exemple assigné au noeud x1 et le faciès c2 étant par exemple assigné au noeud x2, et pour trois faciès possibles c1, c2, c3 pour xo, la table de probabilités bivariées comporte trois valeurs pour chaque couple (xo, x1) et (xo, x2), soit 6 valeurs : P(ko=1),(ki=1) P(ko=2),(k,=1) , P(ko=3),(k1=1) P(ko=1),(k,=2) , P(ko=2),(k,=2) ' et p(ko=3),(k2=2)30 Avec un procédé selon l'art antérieur, il aurait fallu extraire en outre une valeur supplémentaire P(k,-1)(kz=2) de probabilités bivariées pour le couple (xi, x2). On conçoit que pour un nombre de noeuds du voisinage plus élevé, par exemple de l'ordre de la dizaine, le gain en calcul et en mémoire fourni par le procédé selon un mode de réalisation de l'invention est bien plus élevé encore. Les figures 2A et 2B illustrent un avantage ainsi procuré par un procédé selon un mode de réalisation de l'invention. Sur ces figures, chaque double flèche entre deux noeuds représente à une ou plusieurs valeur(s) de probabilités bivariées à estimer et se rapportant à ces deux noeuds. Le procédé selon un mode de réalisation de l'invention permet d'éviter de calculer les probabilités bivariées se rapportant à deux noeuds du voisinage de xo et distincts de xo. Seules les probabilités bivariées se rapportant à deux noeuds dont l'un est le noeud cible xo sont à estimer. Lors de l'étape (cl), on extrait également des valeurs de probabilités univariées Pk, ici trois valeurs de probabilités univariées puisque on envisage trois faciès possibles c,, c2, c3. Lors d'une étape (c2), on estime pour chaque faciès cko envisageable une probabilité conditionnelle Pk0 Iki, k,, La probabilité conditionnelle est la probabilité d'observer le faciès cko en xo connaissant les faciès Cki assignés aux n noeuds x; situés dans le voisinage de xo. Une boucle peut être mise en place pour parcourir les ne faciès possibles, ici n,=3, avec des étapes classiques d'initialisation, de test et d'incrémentation. Cette étape (c2) fait intervenir les probabilités univariées et les probabilités bivariées calculées aux étapes précédentes. Par exemple, la probabilité conditionnelle Pko Ik,,...kä d'observer un faciès cko en xo connaissant les faciès Ckl,...Ckn aux n noeuds x; voisins du noeud xo, est calculée en utilisant la formule suivante :30 n pko JJPk0,k, pk' [T pk' k. k0 =1 o 0 i=1 Selon un autre exemple, on peut utiliser la formule suivante : n (1ùn) T'T pko pko , ki i=1 pko / kl ,...,kn ne n (1ûn) pko pko ,kl k0 =1 i=1 Ces formules se calculent relativement facilement car elles font intervenir uniquement des valeurs de probabilités univariées et des valeurs de probabilités bivariées précédemment extraites des fonctions analytiques du modèle. Le nombre de facteurs pour chaque produit est du même ordre que le nombre de noeuds voisins considérés, de sorte qu'il est possible de considérer un nombre de noeuds voisins relativement élevé, sans alourdir les calculs outre mesure. Les figures 5A à 5J montrent des exemples de grilles bidimensionnelles, les couleurs des noeuds correspondant à des faciès. Les figures 5C et 5D montrent des résultats d'estimation. Pour obtenir ces résultats d'estimation, on estime pour chaque noeud cible et pour chaque faciès possible une valeur de probabilité conditionnelle, selon la méthode de P. Bogaert (figure 5C) ou selon l'étape c/ décrite ci-dessus (figure 5D). Puis, pour chaque noeud cible, on affecte à ce noeud cible le faciès correspondant à la probabilité conditionnelle la plus élevée.
On suppose certaines valeurs de noeuds connues, par exemple de données observées. Les noeuds ainsi initialement renseignés sont représentés à la figure 5B. i=1 pko / ki ,...,k,, La figure 5A montre une grille de référence, à partir de laquelle les faciès aux noeuds supposés connus ont pu être extraits. Ces estimations conduisent à des résultats relativement proches de la grille de référence (figure 5A). Ces procédés semblent donc relativement fiables pour la simulation de réservoirs inconnus. On constate que les grilles des figures 5C et 5D sont relativement semblables, c'est-à-dire que l'estimation menée selon l'étape c/ décrite ci-dessus (figure 5D) conduit à un résultat relativement similaire à celui d'une estimation menée selon la méthode de P. Bogaert (figure 5C). Pour plus de 98% des noeuds, les différences relatives de probabilités conditionnelles sont inférieures à 0,0001. En outre, le procédé conduisant aux résultats de la figure 5D peut être exécuté avec un temps de calcul relativement réduit, ici 2 minutes 30 secondes environ, à comparer aux 30 minutes nécessaires pour exécuter le procédé conduisant aux résultats de la figure 5C. Les figures 5E à 5J montrent des exemples de résultats de simulation, c'est-à-dire qu'un tirage au sort pondéré est effectué avant l'assignation. Ces résultats sont obtenus selon un procédé connu de l'art antérieur (figures 5E, 5G et 51), et selon le procédé décrit ci-dessus (figures 5F, 5H et 5J). Les figures 5E et 5F, 5G et 5H, et 51 et 5J sont à comparer deux à deux. On peut constater que les résultats sont relativement similaires. On pourra noter que le procédé décrit ci-dessus conduit à des résultats légèrement moins pixellisés que le procédé de l'art antérieur, permettant ainsi d'obtenir des régions plus cohérentes géographiquement.
Les figures 6A à 6C illustrent un exemple de procédé mettant en oeuvre une approche multi-grilles, selon un mode de réalisation de l'invention. Dans ce mode de réalisation, on extrait M sous-grilles de la grille fournie à l'étape (a), ici M=2. Chaque sous-grille est composée de noeuds deux fois moins espacés entre eux que les noeuds de la sous-grille précédente. Ainsi, la sous-grille de niveau M=2, représentée par les noeuds pleins sur la figure 6A, correspond à un pas maximal, de 2M=4 fois le pas de la grille. La sous-grille de niveau M-1=1, représentée par les noeuds pleins sur la figure 6B, correspond à un pas de 2""' = 2 fois le pas de la grille. Enfin, la sous-grille de niveau 0, représentée à la figure 6C, est la grille définie à l'étape (a). Les étapes (c), (d), (e) de simulation stochastique sont d'abord appliquées itérativement aux noeuds de la 2ème sous-grille (figure 6A), de façon à assigner un faciès à chaque noeud de cette sous-grille. On assigne alors par simulation stochastique des faciès aux noeuds de la 1 ère sous-grille (figure 6B), et ce en utilisant la connaissance des faciès assignés aux noeuds de la 2ème sous-grille. Enfin, on assigne des faciès aux noeuds de la grille initialement définie (figure 6C), en utilisant la connaissance des faciès assignés aux noeuds des sous-grilles des figures 6A et 6B. Cette approche multi-grilles permet d'obtenir des résultats de simulation de relativement bonne qualité, car on se place dans une configuration favorable au réseau en étoiles. L'ordre selon lequel les noeuds de la grille sont parcourus peut en effet jouer sur la qualité de la simulation, puisque l'on utilise les faciès précédemment assignés pour assigner un faciès à un noeud cible. On pourrait alternativement envisager d'assigner des valeurs de faciès par simulation stochastique à chaque noeud de la grille en parcourant la grille de noeud voisin en noeud voisin. Néanmoins, pour un noeud cible, les valeurs de probabilité conditionnelle risquent d'être biaisées dans la mesure où l'on disposerait d'informations de faciès très principalement pour seulement un certain voisinage (les noeuds précédents auxquels un faciès a déjà été affecté). Pour revenir à la figure 6A, lors de la simulation stochastique d'un faciès pour un noeud cible de la 2ème sous-grille, on peut prévoir d'utiliser les noeuds renseignés au voisinage du noeud cible seulement s'ils appartiennent à cette sous-grille. Alternativement et préférentiellement, on accepte d'utiliser les noeuds auxquels un faciès a été initialement assigné et qui se situent au voisinage du noeud cible, même s'ils n'appartiennent pas à cette 2ème sous-grille. Les étapes (c), (d), (e) sont ensuite appliquées aux noeuds de la l ère sous-grille auxquels aucun faciès n'a été assigné, puis aux noeuds de la grille initiale auxquels aucun faciès n'a été assigné.
La figure 7 montre un exemple de dispositif selon un mode de réalisation de l'invention. Ce dispositif comporte un ordinateur 70 relié par des moyens de communication 71 à des capteurs 72 installés dans des puits forés dans une zone géologique 73. Ces capteurs 72 permettent de fournir des données observées, à partir desquelles un processeur 75 de l'ordinateur 70 peut estimer des faciès. L'ordinateur comporte une unité de mémoire 74 pour stocker une représentation du réservoir sous forme de grille. Plus précisément, la mémoire 74 est agencée pour stocker des coordonnées spatiales des noeuds d'une grille. La mémoire 74 est agencée pour stocker un nombre fini de faciès. Le processeur 75 détermine, à partir de données de localisation des puits, à quel(s) noeud(s) de la grille affecter le ou les faciès estimés à partir des données observées. La mémoire 74 est agencée de façon à associer ce ou ces noeud(s) de la grille cette ou ces valeur(s) de faciès. Le processeur 75 est en outre agencé pour exécuter les étapes (c), (d), (e) du procédé décrit en référence aux figures 1A et 1 B. En particulier, le processeur 75 est capable d'exécuter une instruction de tirage au sort. L'ordinateur 70 comporte en outre un écran 76 pour afficher les résultats de la simulation, par exemple sous la forme de grilles similaires à celles des figures 5C et 5D. Le géologue peut étudier la grille ainsi obtenue et en tirer des conclusions quant aux localisations d'hydrocarbures. La simulation du réservoir peut ainsi constituer une aide à la prospection, et de manière plus générale une aide pour estimer l'état d'un sous-sol, par exemple pour estimer des quantités d'hydrocarbures présents. En particulier, le procédé selon un aspect de l'invention peut être mis en oeuvre pour simuler un champ déjà exploité, à des fins d'estimations des quantités d'hydrocarbures restantes et des localisations de ces hydrocarbures. f

Claims (12)

  1. REVENDICATIONS1. Dispositif (70) d'estimation de propriétés lithologiques d'une zone 5 géologique, comprenant - une unité de mémoire (74) pour stocker des coordonnées spatiales de noeuds d'une grille représentant la zone géologique, et un nombre fini de valeurs correspondant à des états d'une variable catégorielle représentant une propriété lithologique, 10 - des moyens de traitement (75), agencés pour itérativement, chaque itération correspondant à un noeud cible de la grille : • estimer, pour chaque état possible de la variable catégorielle, une valeur de probabilité conditionnelle d'observer ledit état pour le noeud cible sachant qu'un état respectif a été assigné à au moins un 15 noeud voisin du noeud cible, à partir - de probabilités univariées d'observer un état de la variable catégorielle pour le noeud cible, et - de probabilités bivariées d'observer deux états de la variable catégorielle pour deux noeuds respectifs de la grille, 20 dans lequel les valeurs de probabilités bivariées intervenant dans l'estimation de la probabilité conditionnelle se rapportent au noeud cible et à un noeud voisin du noeud cible, • effectuer un tirage au sort parmi les états possibles de la variable catégorielle, pondéré par les probabilités conditionnelles ainsi 25 estimées, et • assigner au noeud cible l'état ainsi tiré au sort.
  2. 2. Procédé d'estimation de propriétés lithologiques d'une zone géologique simulée par ordinateur, comprenant 30 a/ représenter la zone géologique à l'aide d'une grille comprenant des noeuds, b/ représenter une propriété lithologique à l'aide d'une variable catégorielle, et pour chaque un noeud cible (xo) parmi un ensemble de noeuds cibles de la grille à parcourir itérativement : c/estimer pour chaque état possible (CO de la variable catégorielle une valeur de probabilité conditionnelle (i' ko ) d'observer ledit état (cko) pour le noeud cible (xo) sachant qu'un état respectif (ckl,...ckn) a été affecté à au moins un noeud voisin (xi,...xn) du noeud cible, à partir de probabilités univariées (pk) d'observer un état (CO de la variable catégorielle pour le noeud cible (xo), et - de probabilités bivariées (Pko,k,) d'observer deux états (cko, cki) de la variable catégorielle pour deux noeuds respectifs de la grille (xo, x;), dans lequel les valeurs de probabilités bivariées (Pko,k,) intervenant dans l'estimation de la probabilité conditionnelle (Pko kl kn ) se rapportent au noeud cible (xo), et à un noeud (x;) voisin du noeud cible, d/ effectuer un tirage au sort parmi les états de la variable, pondéré par les probabilités conditionnelles estimées à l'étape cl, et e/ assigner au noeud cible (xo) l'état tiré au sort à l'étape dl.
  3. 3. Procédé selon la revendication 2, comprenant, préalablement à l'étape d'estimation de la probabilité conditionnelle : prévoir (1) un modèle donnant, pour chaque couple d'états (ck, ck') de la variable catégorielle, une valeur de probabilités bivariées d'observer lesdits deux états (ck, ck') pour deux noeuds respectifs de la grille (xi, xi+h), en fonction des positions relatives desdits noeuds l'un par rapport à l'autre, et, lors de l'étape cl d'estimation de la probabilité conditionnelle, pour chaque état possible (CkO), et pour chaque noeud du voisinage (x;) auquel est assigné un état donné (c k ), extraire du modèle la valeur de probabilité bivariée (Pko,kr) d'observer l'état possible et l'état donné pour deux noeuds ayant les mêmes positions relatives l'un par rapport à l'autre que le noeud cible (xo) et le noeud du voisinage (xi).
  4. 4. Procédé selon la revendication 3, dans lequel le modèle de probabilités bivariées est anisotrope.
  5. 5. Procédé selon l'une des revendications 2 à 4, dans lequel, à l'étape cl, pour chaque état (cko) possible de la variable catégorielle, la probabilité conditionnelle correspondante ( pko kl,...kn ) s'écrit : n pko,kl (1ùn) i=1 pko _ pko / kl ,...,kn ù n~ n (1ùn) pk0 pk0' ,ki k0 =1 i=1 Où ko est un indice de l'état possible cko pour le pixel cible xo, sont les indices des états ckl,...Ckn assignés à n noeuds respectifs voisins (x1,...xn) du pixel cible (xo), n étant supérieur ou égal à 1, ne est le nombre d'états possibles de la variable catégorielle, ne étant strictement supérieur à 1, pko est la probabilité univariée pour un noeud d'avoir l'état cko, 20 pk est la probabilité univariée pour un noeud d'avoir l'état ck.,Pko,k; est la probabilité bivariée d'observer les états C ko et C k, pour le noeud cible xo et un noeud x; situé dans le voisinage du noeud cible xo, et pko,k; est la probabilité bivariée d'observer les états c ko et c k, pour 5 le noeud cible xo et le noeud x; situé dans le voisinage du noeud cible xo.
  6. 6. Procédé selon l'une des revendications 2 à 5, dans lequel la propriété lithologique représentée par la variable catégorielle est une propriété de faciès de roche.
  7. 7. Procédé selon l'une des revendications 2 à 6, comprenant en outre une étape d'initialisation (6) au cours de laquelle on assigne un état de la variable catégorielle à au moins un noeud de la grille. 15
  8. 8. Procédé selon la revendication 7, dans lequel, à l'étape d'initialisation, pour au moins un noeud, l'état assigné audit noeud est obtenu suite à une observation de la zone géologique.
  9. 9. Procédé selon la revendication 7 ou 8, dans lequel, à l'étape 20 d'initialisation, pour au moins un noeud, l'état assigné audit noeud est obtenu par tirage au sort pondéré par des probabilités univariées.
  10. 10. Procédé selon l'une des revendications 2 à 9, comprenant extraire M sous-grilles de la grille définie à l'étape a/, chaque sous-grille 25 étant composés de noeuds deux fois moins espacés entre eux que les noeuds de la sous-grille précédente, la sous-grille de niveau M correspondant à un espacement maximal, et la sous-grille de niveau 0 correspondant à la grille définie à l'étape a/, 10 et dans lequel les étapes cl, dl, e/ sont d'abord appliquées aux noeuds de la M1ème sous-grille, puis, pour m courant de M-1 à 0, les étapes cl, dl, e/ sont appliquées aux noeuds de la mième sous-grille, en utilisant les états assignés aux noeuds de la sous-grille de niveau (m+1).
  11. 11. Procédé selon l'une des revendications 2 à 10, dans lequel on détermine le voisinage du noeud cible à considérer lors de l'étape cl en définissant un domaine géométrique centré autour dudit noeud cible, et en recensant les noeuds situés dans ledit domaine géométrique auxquels un état de variable catégorielle a déjà été assignée.
  12. 12. Programme d'ordinateur comportant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 2 à 11, lorsque lesdites instructions sont exécutées par un processeur.
FR0858021A 2008-11-26 2008-11-26 Estimation de proprietes lithologiques d'une zone geologique Withdrawn FR2938952A1 (fr)

Priority Applications (3)

Application Number Priority Date Filing Date Title
FR0858021A FR2938952A1 (fr) 2008-11-26 2008-11-26 Estimation de proprietes lithologiques d'une zone geologique
US13/131,312 US8630832B2 (en) 2008-11-26 2009-11-25 Estimation of lithological properties of a geological zone
PCT/FR2009/052301 WO2010061135A1 (fr) 2008-11-26 2009-11-25 Estimation de proprietes lithologiques d'une zone geologique.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR0858021A FR2938952A1 (fr) 2008-11-26 2008-11-26 Estimation de proprietes lithologiques d'une zone geologique

Publications (1)

Publication Number Publication Date
FR2938952A1 true FR2938952A1 (fr) 2010-05-28

Family

ID=40984695

Family Applications (1)

Application Number Title Priority Date Filing Date
FR0858021A Withdrawn FR2938952A1 (fr) 2008-11-26 2008-11-26 Estimation de proprietes lithologiques d'une zone geologique

Country Status (3)

Country Link
US (1) US8630832B2 (fr)
FR (1) FR2938952A1 (fr)
WO (1) WO2010061135A1 (fr)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102080539A (zh) * 2011-02-21 2011-06-01 黄桂芝 采用旋转交错形矿产勘查网的矿产勘查方法
WO2013156412A3 (fr) * 2012-04-16 2013-12-19 Total Sa Procédé de détermination d'une carte lithologique
WO2014181145A1 (fr) * 2013-05-07 2014-11-13 Total Sa Procédé servant à déterminer un modèle empilé décrivant des éléments architecturaux

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2795527B1 (fr) * 2012-02-10 2016-04-06 Landmark Graphics Corporation Systèmes et procédés pour la sélection de réalisations de modèle de faciès
US9164193B2 (en) 2012-06-11 2015-10-20 Chevron U.S.A. Inc. System and method for optimizing the number of conditioning data in multiple point statistics simulation
AU2013283825B2 (en) * 2012-06-26 2017-02-02 Total Sa Truncation diagram determination for a pluri-gaussian estimation
CN110927793B (zh) * 2019-11-27 2020-09-29 中国地质大学(北京) 一种基于序贯随机模糊模拟的储层预测方法及系统
CN115680645A (zh) * 2022-09-27 2023-02-03 成都理工大学 基于多源随钻信息融合的岩体特征实时预测方法和系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1930744A1 (fr) * 2006-12-11 2008-06-11 Ifp Méthode de construction d'un modèle géologique d'une formation du sous-sol contraint par des données sismiques

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5416750A (en) * 1994-03-25 1995-05-16 Western Atlas International, Inc. Bayesian sequential indicator simulation of lithology from seismic data
US7516055B2 (en) 2004-08-20 2009-04-07 Chevron U.S.A. Inc Multiple-point statistics (MPS) simulation with enhanced computational efficiency

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1930744A1 (fr) * 2006-12-11 2008-06-11 Ifp Méthode de construction d'un modèle géologique d'une formation du sous-sol contraint par des données sismiques

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
P. BOGAERT: "Spatial prediction of categorical variables: the Bayesian maximum entropy approach", STOCHASTIC ENVIRONMENTAL RESEARCH AND RISK ASSESSMENT, SPRINGER-VERLAG BERLIN / HEIDELBERG, DOI 10.1007/S00477-002-0114-4, vol. 16, no. 6, December 2002 (2002-12-01), pages 429 - 448, XP002543302 *
VARGAS-GUZMAN ET AL: "Spatial conditional simulation of facies objects for modeling complex clastic reservoirs", JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING, ELSEVIER, AMSTERDAM, NL, vol. 54, no. 1-2, 29 October 2006 (2006-10-29), pages 1 - 9, XP005839375, ISSN: 0920-4105 *
WIBRIN M A ET AL: "Combining categorical and continuous spatial information within the Bayesian maximum entropy paradigm", STOCHASTIC ENVIRONMENTAL RESEARCH AND RISK ASSESSMENT, SPRINGER, BERLIN, DE, vol. 20, no. 6, 24 February 2006 (2006-02-24), pages 423 - 433, XP019428988, ISSN: 1436-3259 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102080539A (zh) * 2011-02-21 2011-06-01 黄桂芝 采用旋转交错形矿产勘查网的矿产勘查方法
CN102080539B (zh) * 2011-02-21 2013-05-08 黄桂芝 采用旋转交错形矿产勘查网的矿产勘查方法
WO2013156412A3 (fr) * 2012-04-16 2013-12-19 Total Sa Procédé de détermination d'une carte lithologique
GB2515220A (en) * 2012-04-16 2014-12-17 Total Sa A method for determining a lithology map
GB2515220B (en) * 2012-04-16 2017-10-25 Total Sa A method for determining a lithology map
US10585199B2 (en) 2012-04-16 2020-03-10 Total Sa Method for determining a lithology map
WO2014181145A1 (fr) * 2013-05-07 2014-11-13 Total Sa Procédé servant à déterminer un modèle empilé décrivant des éléments architecturaux
GB2528011A (en) * 2013-05-07 2016-01-06 Total Sa Method for determining a stacked model describing architectural elements
US10310116B2 (en) 2013-05-07 2019-06-04 Total Sa Method for determining a stacked model describing architectural elements
GB2528011B (en) * 2013-05-07 2020-05-27 Total Sa Method for determining a stacked model describing architectural elements

Also Published As

Publication number Publication date
US20110231170A1 (en) 2011-09-22
US8630832B2 (en) 2014-01-14
WO2010061135A1 (fr) 2010-06-03

Similar Documents

Publication Publication Date Title
FR2938952A1 (fr) Estimation de proprietes lithologiques d'une zone geologique
EP2037080B1 (fr) Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité
CA2806955C (fr) Procede d'exploitation d'un gisement a partir d'une technique de selection des positions de puits a forer
EP2761340B1 (fr) Procédé d'extraction de vignette d'une image d'entraînement pour contraindre la modélisation géostatistique multipoint du sous-sol
EP2356493B1 (fr) Procede de modelisation geologique de donnees sismiques par correlation de traces
EP2963235B1 (fr) Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer
CA3043212A1 (fr) Procede pour la detection d'objets geologiques dans une image sismique
FR3045868A1 (fr) Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures
CA2821099C (fr) Procede d'exploitation d'un reservoir geologique a partir d'un modele de reservoir cale par le calcul d'une loi analytique de distribution conditionnelle de parametres incertains du modele
FR3055723A1 (fr) Modelisation basee sur un point-vecteur des proprietes de reservoir de petrole pour un modele de simulation de reservoir sans grille
CA2701524A1 (fr) Methode pour modeliser un milieu heterogene souterrain a partir de statistique multipoint
EP2873989B1 (fr) Procédé de construction d'une grille représentative de la distribution d'une propriété physique d'une formation souterraine par simulation statistique multipoint conditionnelle
FR3002270A1 (fr) Procede d'exploitation d'un reservoir geologique au moyen d'un modele de reservoir cale et coherent vis a vis des proprietes d'ecoulement
CA2711926A1 (fr) Procede, programme et systeme informatique de construction d'un modele geologique 3d
FR3029664A1 (fr) Definition de petrofacies non-lineaires pour un modele de simulation de reservoir
EP3252508B1 (fr) Procédé d'exploitation des hydrocarbures d'une formation souterraine, au moyen d'une mise à l'échelle optimisée
FR3087274A1 (fr) Étalonnage d’images sismiques répétitives pour des opérations de production
EP4252044A1 (fr) Procede pour determiner des incertitudes associees a un modele de bassin sedimentaire
FR3045097A1 (fr)

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 8

PLFP Fee payment

Year of fee payment: 9

ST Notification of lapse

Effective date: 20180731