FR2757957A1 - Methode pour simplifier la modelisation d'un milieu geologique poreux traverse par un reseau irregulier de fractures - Google Patents
Methode pour simplifier la modelisation d'un milieu geologique poreux traverse par un reseau irregulier de fractures Download PDFInfo
- Publication number
- FR2757957A1 FR2757957A1 FR9616331A FR9616331A FR2757957A1 FR 2757957 A1 FR2757957 A1 FR 2757957A1 FR 9616331 A FR9616331 A FR 9616331A FR 9616331 A FR9616331 A FR 9616331A FR 2757957 A1 FR2757957 A1 FR 2757957A1
- Authority
- FR
- France
- Prior art keywords
- medium
- transposed
- transfer function
- equivalent
- geological
- 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.)
- Granted
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Abstract
Méthode permet une modélisation simplifiée d'un milieu géologique hétérogène poreux (tel qu'un réservoir traversé par un réseau irrégulier de fissures par exemple) sous la forme d'un milieu transposé ou équivalent de manière à ce que le milieu transposé soit équivalent au milieu d'origine, relativement à un type déterminé de fonction de transfert physique (connu pour le milieu transposé). Elle comporte: a) la formation d'une image en au moins deux dimensions du milieu géologique d'origine sous la forme d'une série de pixels et l'association à chaque pixel de la série d'une valeur initiale particulière pour cette fonction, b) la détermination par étapes de valeurs à attribuer pour la fonction de transfert physique à chaque pixel de la série (telle que la distance minimale séparant le pixel de la fissure la plus proche) par référence à des valeurs de la fonction attribuées à des pixels voisins de l'image, et c) la détermination d'une propriété physique du milieu transposé ou équivalent en identifiant les valeurs de la fonction de transfert connue pour le milieu transposé avec les valeurs de la fonction de transfert du milieu d'origine déterminées par étapes. La fonction de transfert physique peut représenter des variations entre différentes parties du milieu géologique, par exemple des distances ou des transmissivités ou des transferts thermiques (entre un réservoir et un puits traversant le réservoir par exemple), etc. Application à la détermination d'un milieu transposé offrant la même récupération d'un fluide au cours d'un processus d'imbibition capillaire que le milieu réel.
Description
L'invention concerne une méthode destinée à simplifier la modélihation d'un milieu géologique poreux traversé par un réseau irrégulier de fissures qui simplifie la mise en relation de modèles de caractérisation de réservoirs fissurés et de modèles double porosité. La méthode peut être mise en oeuvre par exemple dans le domaine de la production pétrolière par des ingénieurs de gisement en vue d'obtenir des prédictions d'écoulement fiables.
Les réservoirs fissurés constituent un type extrême de réservoirs hétérogènes comportant deux milieux contrastés, un milieu matriciel contenant la plus grande part de l'huile en place et présentant une faible perméabilité, et un milieu fissuré représentant moins de I % de l'huile en place et hautement conducteur. Le milieu fissuré lui-même peut être complexe, avec différents ensembles de fissures caractérisés par leur densité, longueur, orientation, inclinaison et ouverture respectives. Les images en 3D de réservoirs fissurés ne peuvent pas être utilisées directement sous forme de données d'entrée de simulation de réservoir. La représentation d'un réseau de fissures dans des simulateurs d'écoulement a été longtemps considérée comme irréaliste car la configuration du réseau est en partie inconnue et à cause des limitations numériques liées à la juxtaposition de nombreuses cellules présentant des dimensions et des propriétés extrêmement contrastées. C'est pourquoi la modélisation simplifiée mais réaliste de tels milieux présente un grand intérêt pour les ingénieurs de gisement.
L"'approche double porosité" telle qu'elle est enseignée, par exemple. par Warren J.E. et al dans "The Behavior of Naturally Fractured Reservoirs", SPE Journal (septembre 1963). 245-255, est connue dans l'art pour interpréter le comportement d'un écoulement monophasique observé en testant un réservoir fissuré. Selon ce modèle de base, tout volume élémentaire du réservoir fissuré est modélisé sous la forme d'un ensemble de blocs parallélépipédiques identiques limité par un système orthogonal de fissures uniformes continues oricntées dans la direction de l'un des trois principaux sens d'écoulement. L'écoulement des fluides à l'échelle du réservoir s'effectue à travers le milieu de fissure seulement et des échantes de fluides surviennent localement entre les fissures et les blocs matriciels.
De nombreux simulateurs de réservoirs fissurés ont été développés en utilisant un tel modèle. avec des améliorations spécifiques relatives à la modélisatioti des échanges de fluides entre matrice et fissure régis par des forces capillaires gravitationnelles, visqueuses et des mécanismes compositionnels. ainsi que des échanges matrice-matrice (simulateurs double perméabilité double porosité). De nombreux exemples de techniques de l'art antérieur sont cités dans les références indiquées ci-après.
- Thomas, L.K. et al : "Fractured Reservoir Simulation", SPE Journal (février 1983), 42-54;
- Quandalle, P. et al "Typical Features of a New Multipurpose Reservoir
Simulator", SPE 16007 présenté au 9ème Symposium sur la simulation de réservoirs de
San Antonio, Texas, 1-4 février 1987;
- Coats. K.H.: "Implicit Compositional Simulation of Single-Porosity and Dual
Porosity Reservoirs", SPE 18427 présenté au Symposium SPE sur la simulation de réservoirs de Houston, Texas, 6-8 février 1989.
- Quandalle, P. et al "Typical Features of a New Multipurpose Reservoir
Simulator", SPE 16007 présenté au 9ème Symposium sur la simulation de réservoirs de
San Antonio, Texas, 1-4 février 1987;
- Coats. K.H.: "Implicit Compositional Simulation of Single-Porosity and Dual
Porosity Reservoirs", SPE 18427 présenté au Symposium SPE sur la simulation de réservoirs de Houston, Texas, 6-8 février 1989.
L'un des problèmes que rencontrent les ingénieurs de gisement consiste à doter ce modèle de base de paramètres en vue d'obtenir des prédictions d'écoulement fiables. En particulier, les perméabilités de fissures équivalentes et la taille des blocs matriciels doivent être connues pour chaque cellule du simulateur d'écoulement. Alors que la perméabilité matricielle peut être évaluée à partir de carottes. les perméabilités du réseau de fissures contenu dans la cellule, c'est-à-dire les perméabilités de fissures équivalentes. ne peuvent pas être évaluées simplement et nécessitent la prise en compte de la géométrie et des propriétés du réseau de fissures réel. Une méthode permettant de déterminer les perméabilités de fissures équivalentes d'un réseau de fissures est présentée dans la demande de brevet parallèle EN. 96/16330.
Il existe une procédure de référence connue permettant de détcrmincr les dimensions a.b de chaque bloc d'une section traversée par un maillage régulier de fissures Fet qui est équivalente à la section d'un milieu multicouche naturel fissuré traversé par un rcsc;lu tic @issures FN le long d'un plan dc rèlèrence parallèle aux couches (ce plan étant habituellement horizontal ou sensiblement horizontal). Pour chaquc couche du volume de roche fissurée étudié (ligure I ). les dimensions "horizontales" a. b des blocs de la section équivalente sont déterminées de manière itérative au moyen de calculs et par comparaison des fonctions de récupération d'huile en fonction du temps Rit) et Req t) lespectivement dans la section réelle RE du volume de roche fissurée étudié et dans la section EQ découpée en "morceaux de sucre" (sugar lumps) de taille identique équivalent à la distribution des blocs réels. Cette méthode conventionnelle nécessite un simulateur d'écoulement polyphasique simple porosité discrétisant les blocs matriciels et les fissures de manière à pouvoir comparer les courbes de récupération. Une telle procédure est extrêmement coûteuse dans la mesure où la discrétisation de la section réelle peut comporter un nombre très élevé de cellules. En fait, la forme réelle des blocs doit être représentée au moyen de cellules de fissures fines le long des limites de chaque bloc. La matrice doit également être discrétisée avec un nombre suffisant de cellules en vue d'obtenir une fonction de transfert d'imbibition bloc-fissure précise.
On peut citer différentes techniques de l'art antérieur mises en oeuvre dans ce domaine, par exemple:
- Bourbiaux, B. et al : "Experimental Study of Cocurrent and Countercurrent Flows in
Natural Porous Media", SPE Reservoir Engineering (août 1990), 361-368,
- Cuiec, L. et al : "Oil Recovery by Imbibition in Low-Permeability Chalk", SPE
Formation Evaluation (septembre 1994), 200-208.
- Bourbiaux, B. et al : "Experimental Study of Cocurrent and Countercurrent Flows in
Natural Porous Media", SPE Reservoir Engineering (août 1990), 361-368,
- Cuiec, L. et al : "Oil Recovery by Imbibition in Low-Permeability Chalk", SPE
Formation Evaluation (septembre 1994), 200-208.
Cependant, les caractéristiques spécifiques d'imbibition n'ont encore jamais été utilisées pour trouver les dimensions du bloc équivalent dans les modèles double porosité. Ainsi, les ingénieurs de gisement manquent d'un outil systématique permettant de calculer les dimensions des blocs parallélépipédiques équivalents à des écoulements polyphasiques pour déterminer la distribution réelle des blocs dans chaque zone de réservoir fissurée.
On connaît également dans l'art antérieur des techniques permettant d'intégoer des données de lissuration naturelle dans des modèles de réservoirs fissurés. Les données dc tissu rat ion sont notamment des données de nature géométrique et elles comportent des mesures de densité, dc longueur, d'azimut et d'inclinaison des plans dc tissure soit observés sur des affleurements. dans des galeries de mines. sur des carottes. soit déduits à partir de diagraphies. Différentes ensembles de fissures peuvent étire distingués et caractérisés par différentes distributions statistiques des caractères des fissures. Les configurations des tissures Line fois caractérisées, on peut réaliser des réseaux numériques de ces ensembles de tissures en utilisant un procédé stochastique respect lut les distributions statistiques des paramètres de fissures. De tels procédés sont présentés par exemple dans les brevets FR-A-2.725.8 14: 2.725.794 ou 2.733.073 du demandeur.
La méthode selon l'invention permet une modélisation simplifiée d'un milieu géologique d'origine hétérogène poreux (tel qu'un réservoir traversé par un réseau irrégulier de fissures par exemple) sous la forme d'un milieu transposé ou équivalent de telle manière que le milieu transposé soit équivalent au milieu d'origine relativement à un type déterminé de fonction de transfert physique connue pour le milieu transposé. La méthode comporte
- la formation d'une image à au moins deux dimensions du milieu géologique sous la forme d'une série de pixels et l'association à chaque pixel de la série d'une valeur initiale particulière pour ladite fonction,
- la détermination par étapes de la valeur à attribuer pour la fonction de transfert physique à chaque pixel de ladite série, en référence à des valeurs de la fonction attribuées à des pixels voisins de l'image, et
- la détermination d'une propriété physique du milieu transposé ou équivalent par identification des valeurs de la fonction de transfert connue pour le milieu transposé (simplifié) avec la valeur de la fonction de transfert déterminée par étapes pour le milieu d'origine.
- la formation d'une image à au moins deux dimensions du milieu géologique sous la forme d'une série de pixels et l'association à chaque pixel de la série d'une valeur initiale particulière pour ladite fonction,
- la détermination par étapes de la valeur à attribuer pour la fonction de transfert physique à chaque pixel de ladite série, en référence à des valeurs de la fonction attribuées à des pixels voisins de l'image, et
- la détermination d'une propriété physique du milieu transposé ou équivalent par identification des valeurs de la fonction de transfert connue pour le milieu transposé (simplifié) avec la valeur de la fonction de transfert déterminée par étapes pour le milieu d'origine.
La fonction de transfert physique peut représenter des variations entre différentes parties du milieu géologique. par exemple des variations de distances, de transmissivités ou de chaleur (tels que des transferts thermiques entre un réservoir et un puits traversant ce réservoir). etc.
La méthode peut être appliquée par exemple pour déterminer. à partir d'une image d'un milieu géologique poreux traversé par un réseau irrégulier de fissures. Lin milieu transposé comportant un ensemble de blocs disposés de manière régulière et séparés pal- un maillage régulier de fissures lequel milieu transposé donne sensiblement une meme récupération de fluide au cours d'un processus d'imbibition capillaire que Ic milieu réel.
Dans ce cas, la méthode comporte
- la formation d'une image à au moins deux dimensions du milieu réel sous la forrne d'une série de pixels.
- la formation d'une image à au moins deux dimensions du milieu réel sous la forrne d'une série de pixels.
- la détermination. pour chaque pixel. de la distance minimale séparant le pixel de la tissure la plus proche.
- la formation d'une distribution de nombres de pixels en fonction des distances minimales aux fissures et la détermination, à partir de cette distribution. de la fonction de récupération (R) dudit ensemble de blocs, et
- la détermination de dimensions (a, b) des blocs réguliers équivalents de l'ensemble à partir de la fonction de récupération (R) et de la fonction de récupération (Req) de l'équivalent (en utilisant par exemple une procédure d'identification desdites fonctions de récupération).
- la détermination de dimensions (a, b) des blocs réguliers équivalents de l'ensemble à partir de la fonction de récupération (R) et de la fonction de récupération (Req) de l'équivalent (en utilisant par exemple une procédure d'identification desdites fonctions de récupération).
Avec la méthode définie ci-dessus utilisant une représentation du milieu par pixels, de nombreuses fonctions de transfert différentes appliquées à tout type de milieu hétérogène peuvent être facilement et rapidement calculées.
La méthode géométrique par exemple permet de trouver des dimensions de blocs équivalents permettant une très bonne adéquation du comportement d'imbibition du bloc réel ou de la distribution des blocs réels, quelle que soit la forme (ou les formes) des blocs considérés. La courbe de récupération d'huile calculée sur la section de bloc équivalent, bien que simplifiée par rapport aux méthodes antérieures, se révèle très proche de celle calculée sur la section de bloc réel.
D'autres caractéristiques et avantages de la méthode selon l'invention apparaîtront plus clairement à la lecture de la description ci-après de modes de réalisation donnés à titre d'exemples non limitatifs, en se référant aux schémas figurant en annexe parmi lesquels
- la figure I illustre une procédure connue permettant de déterminer un milieu présentant des fissures régulières équivalent au milieu fissuré réel,
- la figure 2 illustre la procédure selon l'invention permettant de déterminer un milieu présentant des fissures régulières équivalent à un milieu fissuré réel,
- la figure 3 présente un exemple de pixels voisins servant au calcul de la valeur attribuée à un pixel.
- la figure I illustre une procédure connue permettant de déterminer un milieu présentant des fissures régulières équivalent au milieu fissuré réel,
- la figure 2 illustre la procédure selon l'invention permettant de déterminer un milieu présentant des fissures régulières équivalent à un milieu fissuré réel,
- la figure 3 présente un exemple de pixels voisins servant au calcul de la valeur attribuée à un pixel.
- la figure 4 présente un histogramme d'une distribution possible de pixels par rapport à la distance aux fissures.
- la figure 5 présente une variation possible d'une zone envahie normée en fonction de la distance aux fissures,
- la figure 6 montre un autre exemple de pixels voisins dans trois plans différents Sk1, Sk et Sk+ I utilisés pour un calcul 3D de valeurs attribuées à un pixel,
- la figure 7 présente une extension possible de pixels voisins permettant d'améliorer le calcul des valeurs attribuées à un pixel, et
- la figure 8 montre, à des fins de validation, la bonne adéquation entre deux courbes de récupération d'huile OR(t) déterminées en utilisant d'une part un bloc "en forme de peigne" réel et d'autre part un bloc rectangulaire équivalent.
- la figure 6 montre un autre exemple de pixels voisins dans trois plans différents Sk1, Sk et Sk+ I utilisés pour un calcul 3D de valeurs attribuées à un pixel,
- la figure 7 présente une extension possible de pixels voisins permettant d'améliorer le calcul des valeurs attribuées à un pixel, et
- la figure 8 montre, à des fins de validation, la bonne adéquation entre deux courbes de récupération d'huile OR(t) déterminées en utilisant d'une part un bloc "en forme de peigne" réel et d'autre part un bloc rectangulaire équivalent.
Une nouvelle procédure simplifiée permettant de calculer les dimensions de la section de bloc équivalente à la section "horizontale" d'un milieu fissuré naturel est présentée ciaprès.
En premier lieu, il convient d'indiquer que, selon l'hypothèse de fissures "verticales", c'est-à-dire perpendiculaires aux plans de stratification, le milieu matriciel est continu d'une couche géologique à une autre, et le problème lié à la découverte de dimensions de blocs équivalents devient un problème bidimensionnel. Ainsi, le problème auquel on est confronté ici consiste à déterminer la section carrée ou rectangulaire équivalente de blocs matriciels numériques pour chaque couche ou groupe de couches présentant des propriétés de fissuration similaires.
Deuxièmement. l'équivalence d'un modèle double porosité par rapport à un réservoir fissuré doit être établie pour le comportement de l'écoulement, Les écoulements au sein de réservoirs d'huile fissurés sont essentiellement polyphasiques au cours de l'exploitation du gisement. avec deux principaux mécanismes de drainage pour la récupélatioll de l'huile dans la matrice. l'inibibition capillaire et le drainage par gravité. Les ettcts de ces deux mécanismes se conjuguent dans le cas de procédés de récupélation eau-huile qui demeurent une stratégie prédominante dans le cadre du développement de nombreux résel\olrs fissurés. Des mécanismes compositionnels tels que la diffusion Intel-\iennent également dans les procédés de récupération de gaz. C'est pourquoi la méthode géométrique décrite ci-après pour déterminer les blocs équivalents est basée sur des conccpts d'écoulement polyphasique.
Un mode de réalisation de l'invention est décrit ci-dessous en relation avec la figure 2. qui consiste sensiblement en la mise en adéquation de la fonction de récupération d'huile
R(t) (du milieu fissuré naturel) obtenue par la méthode de référence citée et de la fonction de récupération connue Req(t) pour le milieu transposé, pour un processus d'imbibition eau-huile diphasique (au cours d'un mécanisme de drainage par imbibition capillaire eau-huile). Cette mise en adéquation est réalisée pour chaque couche du milieu fissuré, puis pour des ensembles de n couches. Dans ce cas, la fonction de récupération
R(t) obtenue est la somme des différentes fonctions Rn(t) des n couches pondérées par les épaisseurs correspondantes Hn. Les fissures étant verticales, seules les dimensions horizontales du bloc équivalent sont déterminées. La concordance des fonctions R(t) et
Req(t) est donc un problème bidimensionnel.
R(t) (du milieu fissuré naturel) obtenue par la méthode de référence citée et de la fonction de récupération connue Req(t) pour le milieu transposé, pour un processus d'imbibition eau-huile diphasique (au cours d'un mécanisme de drainage par imbibition capillaire eau-huile). Cette mise en adéquation est réalisée pour chaque couche du milieu fissuré, puis pour des ensembles de n couches. Dans ce cas, la fonction de récupération
R(t) obtenue est la somme des différentes fonctions Rn(t) des n couches pondérées par les épaisseurs correspondantes Hn. Les fissures étant verticales, seules les dimensions horizontales du bloc équivalent sont déterminées. La concordance des fonctions R(t) et
Req(t) est donc un problème bidimensionnel.
1 Formulation géométrique
Les fissures étant définies par les coordonnées de leurs points limites sur une coupe bidimensionnelle XY d'une couche, le processus d'imbibition par lequel de l'eau est présente dans les fissures et de l'huile est présente dans les blocs matriciels doit être déterminé. On suppose que l'invasion de la matrice par de l'eau est du type à piston. On admet que la fonction x=f(t) qui relie l'avancée du front d'eau au temps est la même pour tous les blocs matriciels, quelle que soit leur forme, et pour tous les blocs élémentaires.
Les fissures étant définies par les coordonnées de leurs points limites sur une coupe bidimensionnelle XY d'une couche, le processus d'imbibition par lequel de l'eau est présente dans les fissures et de l'huile est présente dans les blocs matriciels doit être déterminé. On suppose que l'invasion de la matrice par de l'eau est du type à piston. On admet que la fonction x=f(t) qui relie l'avancée du front d'eau au temps est la même pour tous les blocs matriciels, quelle que soit leur forme, et pour tous les blocs élémentaires.
Par conséquent, la mise en adéquation des fonctions R(t) et Req(t) équivaut à la mise en adéquation des fonctions R(x) et Req(x). Ces fonctions définissent physiquement des zones normées envahies par l'eau en fonction de l'avancée du front dtimbibiìioll dans le milieu fissuré.
où a et b sont les dimensions du bloc rectangulaire ou carré équivalent (a et b > 0)
La fonction R(x) n'a pas d'expression analytique. Elle est calculée à partir d'une discrétisation de la coupe XY de la couche étudiée suivant l'algorithme défini ci-après.
La fonction R(x) n'a pas d'expression analytique. Elle est calculée à partir d'une discrétisation de la coupe XY de la couche étudiée suivant l'algorithme défini ci-après.
2 Algorithme de calcul de la fonction R(x)
La coupe XY de la couche étudiée est considérée comme une image dont chaque pixel représente un élément de surface. Ces pixels sont espacés régulièrement d'un pas dx dans la direction X et dy dans la direction Y. L'algorithme mis en oeuvre vise à déterminer, pour chaque pixel de cette image, la distance minimale qui le sépare de la fissure la plus proche.
La coupe XY de la couche étudiée est considérée comme une image dont chaque pixel représente un élément de surface. Ces pixels sont espacés régulièrement d'un pas dx dans la direction X et dy dans la direction Y. L'algorithme mis en oeuvre vise à déterminer, pour chaque pixel de cette image, la distance minimale qui le sépare de la fissure la plus proche.
L'image est traduite par un tableau de nombres réels à deux dimensions : Pict[0 :nx+1,0 :ny+l] où nx et ny sont les nombres de pixels de l'image dans les directions X et
Y. En pratique, le nombre total de pixels (nx.ny) est par exemple de l'ordre du million.
Y. En pratique, le nombre total de pixels (nx.ny) est par exemple de l'ordre du million.
Les valeurs des éléments du tableau Pict sont les distances recherchées.
Initialisation : Tous les pixels par lesquels passe une fissure sont à une distance nulle de la fissure la plus proche. Pour ces pixels, le tableau Pict est donc initialisé à la valeur 0. Ceci est fait par un algorithme connu en soi (I'algorithme de Bresline par exemple) auquel on donne les coordonnées des pixels correspondants aux deux extrémités d'une fissure considérée comme un segment de droite et qui initialise (à 0 dans le cas présent) les pixels les plus proches. Les autres éléments de Pict sont initialisés à une valeur supérieure à la plus grande distance existant entre deux pixels de l'image. Cette valeur est par exemple nx.dx+ny.dy.
Calcul : Pour un pixel donné. le calcul de la distance recherchée à la fissure la plus proche se fait à partir des valeurs de distance déjà calculées pour les pixels V()iSII. On lui affecte une valeur qui. si elle s'avère inférieure à la valeur qui lui a été assignée initialement. est le minimum des valeurs des pixels voisins auxquelles on ajoute la distance de ces pixels à celui considéré.
Ce calcul est réalisé en deux phases successives. Lors de la passe descendante. on parcourt l'image ligne par ligne. de haut en bas et de gauche à droite (de Pict @@. I . à @ à Pict [nx.nyj). Les pixels dont on tient compte sont différents selon que l'on est dans une passe descendante ou une passe montante. Comme le montre la figure 3. les pixels en noir et en grisé sont ceux que l'on prend en compte respectivement durant les passes descendantes et les passes ascendantes. pour le pixel Px.
L'écart oblique dxy étant défini comme
l'algorithme s'écrit
pour j= I àny i pouri=l ânx
Pict[ij]= min Pict[i-1,j]+dx, :passe descendante i I Pict[i-l,j-l]+dxy,
Pict[i,j-1]+dy,
Pict[i+1,j-1]+dxy,
I l Pict[i,j] I fin de boucle sur i fin de boucle surj pourj=ny à 1,
pour i=1x à 1,
Pict[i,j]=min Pict[i+1,j]+dx, :passe descendante
Pict[i+1,j+1]+dxy,
Pict[i,j+1]+dy, i I Pict[i-l,j+I] + dxy,
Pict[i,j] I fin de boucle sur i fin de boucle surj
Histogramme : A partir du tableau Pict ainsi calculé, on peut tracer un histogramme en classant les valeurs non nulles (celles affectées aux pixels hors des fissures) par ordre croissant.
l'algorithme s'écrit
pour j= I àny i pouri=l ânx
Pict[ij]= min Pict[i-1,j]+dx, :passe descendante i I Pict[i-l,j-l]+dxy,
Pict[i,j-1]+dy,
Pict[i+1,j-1]+dxy,
I l Pict[i,j] I fin de boucle sur i fin de boucle surj pourj=ny à 1,
pour i=1x à 1,
Pict[i,j]=min Pict[i+1,j]+dx, :passe descendante
Pict[i+1,j+1]+dxy,
Pict[i,j+1]+dy, i I Pict[i-l,j+I] + dxy,
Pict[i,j] I fin de boucle sur i fin de boucle surj
Histogramme : A partir du tableau Pict ainsi calculé, on peut tracer un histogramme en classant les valeurs non nulles (celles affectées aux pixels hors des fissures) par ordre croissant.
Le cumulé de cet histogramme donne, pour toute distance délimitant deux intervalles de l'histogramme. le nombre de pixels non nuls dont la valeur est inférieure à cette distance. Dans l'application décrite à un milieu poreux fissuré où cette distance correspond à l'avancée du front d'eau, le cumulé de l'histogramme indique donc l'aire envahie par l'eau. La courbe R(x) est obtenue en divisant ce cumulé par le nombre total de pixels non nuls (pour la normer). Le nombre d'intervalles utilises en abscisse pour l'histogramme correspond au nombre de points de discrétisation de la courbe R(x). On le choisit égal à 500 par exemple.
3 Recherche des dimensions du bloc équivalent
A ce stade. on connaît la fonction R(x) et l'on cherche les paramètres (#,#) (dimensions du bloc équivalent qui minimisent la fonctionnelle):
A ce stade. on connaît la fonction R(x) et l'on cherche les paramètres (#,#) (dimensions du bloc équivalent qui minimisent la fonctionnelle):
où N est le nombre de points de discrétisation de R(x) et (xi) sont les abscisses de ces points de discrétisation.
Discrétisation suivant l'ordonnée de R(x)
Pour donner autant de poids à tous les volumes d'huile récupérés lors de l'imbibition, la courbe R(x) est re-discrétisée suivant un pas constant sur l'axe des ordonnées (figure 5). La suite (xi) utilisée par la fonctionnelle est déduite de cette discrétisation.
Pour donner autant de poids à tous les volumes d'huile récupérés lors de l'imbibition, la courbe R(x) est re-discrétisée suivant un pas constant sur l'axe des ordonnées (figure 5). La suite (xi) utilisée par la fonctionnelle est déduite de cette discrétisation.
Minimisation de la fonctionnelle:
Comme a et b jouent des rôles symétriques dans l'expression Req(a,b,x), on utilise en fait la fonctionnelle
Comme a et b jouent des rôles symétriques dans l'expression Req(a,b,x), on utilise en fait la fonctionnelle
<tb> <SEP> Re <SEP> q(u, <SEP> =uxx+vxx <SEP> u=2x(i+11
<tb> ÌL'ec <SEP> Req(u,v,x) <SEP> < <SEP> I <SEP> l <SEP> -4
<tb> <SEP> axb
<tb>
Niinimiser cette fonctionnelle revient à trouver le couple (u.v) pour lequel J(#.#) = O. Ceci est fait par un algorithme de Newton.
<tb> ÌL'ec <SEP> Req(u,v,x) <SEP> < <SEP> I <SEP> l <SEP> -4
<tb> <SEP> axb
<tb>
Niinimiser cette fonctionnelle revient à trouver le couple (u.v) pour lequel J(#.#) = O. Ceci est fait par un algorithme de Newton.
Ensuite, le couple (a.b) recherché est déduit de )-u.v). Trois cas peuvent se présenter
@) # v > 0 signifie qu'une des valeurs du couple (a.b) est négative. ce qui n'a aucun sens physique. On pose alors v=0 dans l'expression de Req( u. v x ) . ce qui implique que les fissures sont parallèles. L'opération est recommencée et le couple : a. b) est calculé comme suit
@) # v > 0 signifie qu'une des valeurs du couple (a.b) est négative. ce qui n'a aucun sens physique. On pose alors v=0 dans l'expression de Req( u. v x ) . ce qui implique que les fissures sont parallèles. L'opération est recommencée et le couple : a. b) est calculé comme suit
2) Le cas u2 + 4v < 0 est dépourvu de tout sens physique également puisqu'il signifie que (a,b) ne sont pas réels. On pose alors u2 + 4v = 0, ce qui impose que le bloc élémentaire recherché a la forme d'un carré (a=b). Après minimisation, le couple (a.b) est calculé comme suit:
3) Pour les autres valeurs du couple (#,#), on a
3) Pour les autres valeurs du couple (#,#), on a
4 Validation de la détermination d'un bloc équivalent
La méthode géométrique indiquée ci-avant pour une fonction de transfert d'imbibition a été validée par rapport à une méthode de référence conventionnelle d'un coût très élevé qui nécessite un simulateur d'écoulement polyphasique simple porosité discrétisant les blocs matriciels et les fissures de telle manière que les courbes de récupération puissent être comparées. Des simulations d'écoulement diphasique conventionnelles ont été réalisées pour valider les solutions fournies par la méthode géométrique décrite plus haut. La validation peut comporter les étapes suivantes
a) Calcul de la fonction de récupération d'huile R(t) pour la section géologique (réelle) avec la méthode conventionnelle (solution dc référence);
b) application de la méthode géométrique à la section réelle. ce qui donne une solution ( t:
c) nouvelle utilisation de la méthode conventionnelle. calcul de la fonction de récupération d'huile Req(t) sur la section de bloc équivalent de dimensions ta.bt préalablement déterminées. et comparaison avec la fonction de récupération d'huile de référence R(t).
La méthode géométrique indiquée ci-avant pour une fonction de transfert d'imbibition a été validée par rapport à une méthode de référence conventionnelle d'un coût très élevé qui nécessite un simulateur d'écoulement polyphasique simple porosité discrétisant les blocs matriciels et les fissures de telle manière que les courbes de récupération puissent être comparées. Des simulations d'écoulement diphasique conventionnelles ont été réalisées pour valider les solutions fournies par la méthode géométrique décrite plus haut. La validation peut comporter les étapes suivantes
a) Calcul de la fonction de récupération d'huile R(t) pour la section géologique (réelle) avec la méthode conventionnelle (solution dc référence);
b) application de la méthode géométrique à la section réelle. ce qui donne une solution ( t:
c) nouvelle utilisation de la méthode conventionnelle. calcul de la fonction de récupération d'huile Req(t) sur la section de bloc équivalent de dimensions ta.bt préalablement déterminées. et comparaison avec la fonction de récupération d'huile de référence R(t).
La méthode géométrique donne des dimensions de bloc équivalent permettant une bonne adéquation avec le comportement à l'imbibition du bloc réel. quelle que soit la forme de bloc considérée. La courbe de récupération d'huile calculée sur la section de bloc équivalent est toujours très proche de celle calculée sur la section de bloc réel comme le montre la figure 8.
Une telle validation a également été réalisée avec succès pour une distribution de blocs présentant différentes autres dimensions et formes.
5 Autres applications de la méthode
Précision du calcul des distances aux fissures
Dans l'algorithme de calcul des distances des pixels aux fissures à partir de l'image 2D. on peut améliorer la précision du calcul en prenant en compte un plus grand nombre de voisins du pixel considéré. La figure 7 montre les pixels voisins dont on tient compte quand on étend l'influence d'un pixel Px à deux lignes et deux colonnes précédant ou suivant le pixel Px. De la même façon, les pixels en noir et en grisé sont ceux que l'on prend en compte respectivement durant les passes descendantes et les passes ascendantes, ceux indiqués par une croix étant éliminés pour cause de redondance.
Précision du calcul des distances aux fissures
Dans l'algorithme de calcul des distances des pixels aux fissures à partir de l'image 2D. on peut améliorer la précision du calcul en prenant en compte un plus grand nombre de voisins du pixel considéré. La figure 7 montre les pixels voisins dont on tient compte quand on étend l'influence d'un pixel Px à deux lignes et deux colonnes précédant ou suivant le pixel Px. De la même façon, les pixels en noir et en grisé sont ceux que l'on prend en compte respectivement durant les passes descendantes et les passes ascendantes, ceux indiqués par une croix étant éliminés pour cause de redondance.
Pour augmenter encore la précision des calculs, la zone d'influence des pixels peut être agrandie davantage (à 3 lignes et 3 colonnes ou plus). Cependant. pour l'utilisatioi présentée plus haut, une telle extension n'apporte pas dans la pratiquc d'amélioration notable sur les résultats finaux.
Extension de la méthode à un objet tridimensionnel
L'algorithme présenté ci-dessus peut être appliqué à un volume. Dans ce cas. chaque pixel représente un élément de volume. Le tableau Pict est remplacé par un tableau à 3 dimensions pict3D[() :nu+l.() :ny+1.0 :nz+l] où nx, ny et nz sont les nombres de pixels en X. Y et Z. Pour le calcul en un pixel Px du plan horizontal nulllélo k. les pixels voisins dont on tient compte lors des passes descendante et ascendante sont représentés sur la figure 6.
L'algorithme présenté ci-dessus peut être appliqué à un volume. Dans ce cas. chaque pixel représente un élément de volume. Le tableau Pict est remplacé par un tableau à 3 dimensions pict3D[() :nu+l.() :ny+1.0 :nz+l] où nx, ny et nz sont les nombres de pixels en X. Y et Z. Pour le calcul en un pixel Px du plan horizontal nulllélo k. les pixels voisins dont on tient compte lors des passes descendante et ascendante sont représentés sur la figure 6.
Extension à une fonction quelconque
Dans l'exemple qui a été développé d'une étude d'un phénomène d'imbibition diphasique (eau-huile par exemple), on a cherché à déterminer la taille des blocs en relation avec les distances de points à la fissure la plus proche. La méthode géométrique selon l'invention peut être aussi utilisée pour d'autres types de transferts entre deux milieux contrastés comme par exemple les transferts thermiques entre un puits et un réservoir. Mais surtout, la fonction "distance entre pixels" utilisée dans l'algorithme précédent peut être remplacée par une fonction quelconque reliant les points de l'image.
Dans l'exemple qui a été développé d'une étude d'un phénomène d'imbibition diphasique (eau-huile par exemple), on a cherché à déterminer la taille des blocs en relation avec les distances de points à la fissure la plus proche. La méthode géométrique selon l'invention peut être aussi utilisée pour d'autres types de transferts entre deux milieux contrastés comme par exemple les transferts thermiques entre un puits et un réservoir. Mais surtout, la fonction "distance entre pixels" utilisée dans l'algorithme précédent peut être remplacée par une fonction quelconque reliant les points de l'image.
Il faut alors connaître, pour tout pixel de l'image, la valeur de cette fonction entre ce pixel et ses voisins pris en compte pour le calcul. Cette fonction peut. par exemple. exprimer les valeurs de transmissivités entre les mailles d'un réservoir dont les centres sont les pixels de l'image et servir au calcul du volume de drainage d'un puits foré dans ce réservoir.
Dans un tel cas, les deux passes ascendante et descendante effectuées dans l'algorithme peuvent s'avérer insuffisantes pour trouver une valeur minimale en tout pixel de l'image. On recommence alors jusqu'à ce que les valeurs calculées ne changent plus.
En reprenant les notations présentées plus haut et en supposant que la fonction
F(i,j.k,l) renvoie la valeur de la fonction entre les pixels (i.j) et (k,l), l'algorithme 2D
temp=Pict[i,j]
Pict[i,j]=min(Pict[i-1,j]+F(i,j,i-1,j), : passe descendante
Pict[i-1,j-1]+F(i,j,i-1,j-1)
Pict[i,j-1]+F(i,j,i,j-1)
Pict[i+1,j-1]+F(i,j,i+1,j-1)
Pict[i,j]
si (Pict[i.j]#temp) change = vrai
fin boucle sur i fin boucle surj
pour j=ny à 1.-1 pour i=nx à 1.-1 temp=Pict[i,j]
Pict[i,j] = min (Pict[i+lj] + F(i,j,i+1,j), : passe ascendante
Pict[i+1,j+1]+F(i,j,i+1,j+1)
Pict[i,j+1]+F(i,j,i,j+1)
Pict[i-1,j+1] + F(i,j,i- I ,j+ I)
I I Pict[i,jj si (Pict[i,j]#temp) change=vrai t fin boucle sur i fin boucle surj fin tant~que
F(i,j.k,l) renvoie la valeur de la fonction entre les pixels (i.j) et (k,l), l'algorithme 2D
temp=Pict[i,j]
Pict[i,j]=min(Pict[i-1,j]+F(i,j,i-1,j), : passe descendante
Pict[i-1,j-1]+F(i,j,i-1,j-1)
Pict[i,j-1]+F(i,j,i,j-1)
Pict[i+1,j-1]+F(i,j,i+1,j-1)
Pict[i,j]
si (Pict[i.j]#temp) change = vrai
fin boucle sur i fin boucle surj
pour j=ny à 1.-1 pour i=nx à 1.-1 temp=Pict[i,j]
Pict[i,j] = min (Pict[i+lj] + F(i,j,i+1,j), : passe ascendante
Pict[i+1,j+1]+F(i,j,i+1,j+1)
Pict[i,j+1]+F(i,j,i,j+1)
Pict[i-1,j+1] + F(i,j,i- I ,j+ I)
I I Pict[i,jj si (Pict[i,j]#temp) change=vrai t fin boucle sur i fin boucle surj fin tant~que
Claims (7)
- - la détermination d'une propriété physique du milieu transposé ou équivalent par identification de valeurs de la fonction de transfert connue pour le milieu transposé par rapport à la valeur de la fonction de transfert déterminée par étapes pour le milieu d'origine.- la détermination par étapes de la valeur à attribuer pour la fonction de transfert physique à chaque pixel de ladite série, en référence à des valeurs de la fonction attribuées à des pixels voisins de l'image, et- la formation d'une image à au moins deux dimensions (2D) du milieu géologique sous la forme d'une série de pixels et l'association à chaque pixel de la série d'une valeur initiale particulière pour ladite fonction,Méthode permettant de simplifier la modélisation d'un milieu géologique poreux hétérogène d'origine sous la forme d'un milieu transposé ou équivalent afin que le milieu transposé soit équivalent au milieu d'origine relativement à un type déterminé de fonction de transfen physique connu pour Ic milieu transposé, la méthode comportantREVENDICATIONS
- 2) Méthode selon la revendication 1, caractérisée en ce que le milieu géologique hétérogène est traversé par un réseau irrégulier de fissures toutes géométriquement définies dans des blocs de formes et de tailles irrégulières.
- 3) Méthode selon l'une des revendications I ou 2, caractérisée en ce que ladite fonction de transfert physique représente une distance entre différentes parties du milieu géologique.
- 4) Méthode selon l'une des revendications I ou 2, caractérisée en ce que ladite fonction de transfert physique représente des transmissivités entre différentes parties du milieu géologique.
- 5) Méthode selon l'une des revendications I ou 2. caractérisée en ce que ladite fonction de transfert physique représente des transferts thermiques entre différentes panies du niilicu géologique tels que des transferts thermiques entre un réservoir et uii pu@ts traversant ce réservoir.
- 6) Méthode selon l'une des revendications I ou 2. caractérisée en ce que ladite fonction de transfert physique représente tout transfert de masse ou de tlux entre différentes parties du milieu géologique.
- 7) Méthode selon la revendication 2, permettant de déterminer, à partir d'une image d'un milieu géologique poreux réel traversé par un réseau irrégulier de fissures. un milieu transposé comportant un ensemble de blocs régulièrement disposés et séparés par un maillage régulier de fissures ledit milieu transposé donnant sensiblement la meme récupération (Req) de fluide lors d'un processus d'imbibition capillaire que le milieu réel. la méthode comportant- la formation d'une image à au moins deux dimensions (2D) du milieu réel sous la forme d'une série de pixels,- la détermination, pour chaque pixel, de la distance minimale séparant le pixel de la fissure la plus proche,- la formation d'une distribution du nombre de pixels par rapport à la distance minimale au milieu fissuré et la détermination, à partir de cette distribution, de la fonction de récupération (R) dudit ensemble de blocs. et- la détermination de dimensions (a.b) des blocs réguliers équivalents du milieu transposé à partir de la récupération (R) et de la récupération (Req) du bloc équivalent
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9616331A FR2757957B1 (fr) | 1996-12-30 | 1996-12-30 | Methode pour simplifier la modelisation d'un milieu geologique poreux traverse par un reseau irregulier de fractures |
GB9727231A GB2322949B (en) | 1996-12-30 | 1997-12-23 | Method of simplifying the task of modelling a porous geological medium crossed by an irregular fracture network |
NO19976129A NO326402B1 (no) | 1996-12-30 | 1997-12-29 | Fremgangsmate for a forenkle modelldannelse av et geologisk, porost medium som krysses av et uregelmessig nett av sprekker |
NL1007916A NL1007916C2 (nl) | 1996-12-30 | 1997-12-29 | Werkwijze ter vereenvoudiging van de modelvorming van een poreus geologisch medium dat doorsneden wordt door een onregelmatig netwerk van breuken. |
US09/000,767 US6064944A (en) | 1996-12-30 | 1997-12-30 | Method for simplifying the modeling of a geological porous medium crossed by an irregular network of fractures |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9616331A FR2757957B1 (fr) | 1996-12-30 | 1996-12-30 | Methode pour simplifier la modelisation d'un milieu geologique poreux traverse par un reseau irregulier de fractures |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2757957A1 true FR2757957A1 (fr) | 1998-07-03 |
FR2757957B1 FR2757957B1 (fr) | 1999-01-29 |
Family
ID=9499399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR9616331A Expired - Lifetime FR2757957B1 (fr) | 1996-12-30 | 1996-12-30 | Methode pour simplifier la modelisation d'un milieu geologique poreux traverse par un reseau irregulier de fractures |
Country Status (5)
Country | Link |
---|---|
US (1) | US6064944A (fr) |
FR (1) | FR2757957B1 (fr) |
GB (1) | GB2322949B (fr) |
NL (1) | NL1007916C2 (fr) |
NO (1) | NO326402B1 (fr) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1008720A1 (fr) * | 1998-12-11 | 2000-06-14 | Institut Francais Du Petrole | Procédé pour créer un modèle des écoulements de fluides dans un milieu poreux multi-couches fissuré |
EP1158312A1 (fr) * | 2000-05-26 | 2001-11-28 | Institut Francais Du Petrole | Méthode pour modéliser des écoulements dans un milieu fracturé traversé par de grandes fractures |
FR2850773A1 (fr) * | 2003-01-31 | 2004-08-06 | Inst Francais Du Petrole | Methode pour modeliser des flux de fluides compressibles dans un milieu poreux multicouches fracture |
FR2858444A1 (fr) * | 2003-07-29 | 2005-02-04 | Inst Francais Du Petrole | Methode pour modeliser les transferts compositionnels et/ou polyphasiques entre la matrice poreuse et les fractures d'un milieu poreux multicouches |
FR2923930A1 (fr) * | 2007-11-19 | 2009-05-22 | Inst Francais Du Petrole | Methode de modelisation d'un milieu geologique poreux traverse par un reseau de fractures. |
EP2072752A1 (fr) | 2007-12-20 | 2009-06-24 | Ifp | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures |
EP2453106A1 (fr) | 2010-11-10 | 2012-05-16 | IFP Energies Nouvelles | Methode pour caracteriser le réseau de fractures d'un gisement fracture et méthode pour l'exploiter |
EP2530493A1 (fr) | 2011-06-01 | 2012-12-05 | IFP Energies Nouvelles | Méthode pour construire un maillage d'un réseau de fractures à partir de diagramme de voronoi |
EP2581767A1 (fr) | 2011-10-12 | 2013-04-17 | IFP Energies Nouvelles | Méthode pour construire un maillage d'un réservoir fracturé avec un nombre limité de noeuds dans le milieu matrice |
EP3144468A1 (fr) | 2015-09-15 | 2017-03-22 | IFP Energies nouvelles | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter |
EP3181804A1 (fr) | 2015-12-17 | 2017-06-21 | IFP Énergies nouvelles | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures |
EP3199749A1 (fr) | 2016-01-26 | 2017-08-02 | IFP Énergies nouvelles | Procede pour l'exploitation d'un gisement de fluide traverse par des fractures au moyen d'une simulation d'ecoulement basee sur un flux d'echange et un facteur correctif |
EP4112875A1 (fr) | 2021-07-01 | 2023-01-04 | IFP Energies nouvelles | Procédé pour simuler les écoulements d'un fluide dans une formation souterraine comportant un réseau de fractures |
FR3127250A1 (fr) | 2021-09-21 | 2023-03-24 | IFP Energies Nouvelles | Procédé pour déterminer la viscosité apparente d’une mousse à injecter en vue de la récupération d’une phase organique présente dans une formation souterraine |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2801710B1 (fr) * | 1999-11-29 | 2002-05-03 | Inst Francais Du Petrole | Methode pour generer un maillage hybride permettant de modeliser une formation heterogene traversee par un ou plusieurs puits |
US6516080B1 (en) | 2000-04-05 | 2003-02-04 | The Board Of Trustees Of The Leland Stanford Junior University | Numerical method of estimating physical properties of three-dimensional porous media |
FR2842321B1 (fr) * | 2002-07-11 | 2008-12-05 | Inst Francais Du Petrole | Methode pour contraindre un champ de permeabilite heterogene representant un reservoir souterrain par des donnees dynamiques |
GB2387000B (en) * | 2002-03-20 | 2005-06-01 | Inst Francais Du Petrole | Method for modelling fluid flows in a multilayer porous medium crossed by an unevenly distributed fracture network |
US7565278B2 (en) | 2006-12-04 | 2009-07-21 | Chevron U.S.A. Inc. | Method, system and apparatus for simulating fluid flow in a fractured reservoir utilizing a combination of discrete fracture networks and homogenization of small fractures |
AU2009308037B2 (en) | 2008-10-24 | 2015-04-30 | Exxonmobil Upstream Research Company | Tracking geologic object and detecting geologic anomalies in exploration seismic data volume |
AU2010229934A1 (en) * | 2009-03-24 | 2011-09-29 | Chevron U.S.A. Inc. | A system and method for characterizing fractures in a subsurface reservoir |
US8781806B2 (en) * | 2009-04-30 | 2014-07-15 | Schlumberger Technology Corporation | Determining elastic and fluid flow properties of a fractured reservoir |
US8301427B2 (en) * | 2009-06-05 | 2012-10-30 | Schlumberger Technology Corporation | Fracture network characterization method |
US8583411B2 (en) | 2011-01-10 | 2013-11-12 | Saudi Arabian Oil Company | Scalable simulation of multiphase flow in a fractured subterranean reservoir as multiple interacting continua |
US10417354B2 (en) * | 2013-12-17 | 2019-09-17 | Schlumberger Technology Corporation | Model order reduction technique for discrete fractured network simulation |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2725814A1 (fr) * | 1994-10-18 | 1996-04-19 | Inst Francais Du Petrole | Methode pour cartographier par interpolation, un reseau de lignes, notamment la configuration de failles geologiques |
FR2725794A1 (fr) * | 1994-10-18 | 1996-04-19 | Inst Francais Du Petrole | Methode pour modeliser la distribution spatiale d'objets geometriques dans un milieu, tels que des failles dans une formation geologique |
FR2733073A1 (fr) * | 1995-04-12 | 1996-10-18 | Inst Francais Du Petrole | Methode pour modeliser un milieu geologique stratifie et fracture |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5038378A (en) * | 1985-04-26 | 1991-08-06 | Schlumberger Technology Corporation | Method and apparatus for smoothing measurements and detecting boundaries of features |
US4745550A (en) * | 1985-08-16 | 1988-05-17 | Schlumberger Technology Corporation | Processing of oriented patterns |
US4926394A (en) * | 1989-10-27 | 1990-05-15 | Western Atlas International, Inc. | Monte Carlo method for estimating lithology from seismic data |
US5416750A (en) * | 1994-03-25 | 1995-05-16 | Western Atlas International, Inc. | Bayesian sequential indicator simulation of lithology from seismic data |
US5539704A (en) * | 1995-06-23 | 1996-07-23 | Western Atlas International, Inc. | Bayesian sequential Gaussian simulation of lithology with non-linear data |
US5671136A (en) * | 1995-12-11 | 1997-09-23 | Willhoit, Jr.; Louis E. | Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects |
-
1996
- 1996-12-30 FR FR9616331A patent/FR2757957B1/fr not_active Expired - Lifetime
-
1997
- 1997-12-23 GB GB9727231A patent/GB2322949B/en not_active Expired - Lifetime
- 1997-12-29 NO NO19976129A patent/NO326402B1/no not_active IP Right Cessation
- 1997-12-29 NL NL1007916A patent/NL1007916C2/nl not_active IP Right Cessation
- 1997-12-30 US US09/000,767 patent/US6064944A/en not_active Expired - Lifetime
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2725814A1 (fr) * | 1994-10-18 | 1996-04-19 | Inst Francais Du Petrole | Methode pour cartographier par interpolation, un reseau de lignes, notamment la configuration de failles geologiques |
FR2725794A1 (fr) * | 1994-10-18 | 1996-04-19 | Inst Francais Du Petrole | Methode pour modeliser la distribution spatiale d'objets geometriques dans un milieu, tels que des failles dans une formation geologique |
FR2733073A1 (fr) * | 1995-04-12 | 1996-10-18 | Inst Francais Du Petrole | Methode pour modeliser un milieu geologique stratifie et fracture |
Non-Patent Citations (10)
Title |
---|
B.J.BOURBIAUX: "Experimental Study of Cocurrent and Countercurrent Flows in Natural Porous Media", SPE JOURNAL, no. 18283, August 1990 (1990-08-01), DALLAS, TEXAS, USA, pages 361 - 368, XP002041400 * |
CUIEC: "Oil Recovery by Imbibition in Low-Permeability Chalk", SPE JOURNAL, no. 20259, September 1994 (1994-09-01), DALLAS, TEXAS, USA, pages 200 - 208, XP002041401 * |
J.C.L.LONG: "A Model for Steady Fluid Flow in Random Three-Dimentional Networks of Disc-Shaped Fractures", WATER RESOURCES RESEARCH, 1105-1115, vol. 21, no. 8, August 1985 (1985-08-01), XP002041398 * |
J.E.WARREN: "The Behavior of Naturally Fractured Reservoirs", SPE JOURNAL, no. 426, September 1963 (1963-09-01), DALLAS, TEXAS, USA, pages 245 - 255, XP002039628 * |
K.H.COATS: "Implicit Compositional Simulation of Single-Porosity and Dual-Porosity Reservoirs", SPE JOURNAL, no. 18427, February 1989 (1989-02-01), DALLAS, TEXAS, USA, XP002039631 * |
L.K.THOMAS: "Fractured Reservoir Simulation", SPE JOURNAL, no. 9305, February 1983 (1983-02-01), DALLAS, TEXAS, USA, pages 42 - 54, XP002039629 * |
M.C.CACAS: "Modeling Fracture Flow With a Stochastic Discrete Fracture Network: Calibration and Validation. 1. The Flow Model", WATER RESOURCES RESEARCH, vol. 26, no. 3, March 1990 (1990-03-01), XP002041399 * |
N.E.ODLING: "Permeability of natural and simulated fracture patterns", STRUCTURAL AND TECTONIC MODELLING AND ITS APPLICATION TO PETROLEUM GEOLOGY, NORWEGIAN PETROLEUM SOCIETY (NPF) SPECIAL PUBLICATION 1, 365-380, ELSEVIER, 1992, AMSTERDAM, XP002041397 * |
P.QUANDALLE: "Typical Features of a Multipurpose Reservoir Simulator", SPE JOURNAL, no. 16007, February 1987 (1987-02-01), DALLAS, TEXAS, USA, pages 475 - 480, XP002039630 * |
S.M.REYNOLDS: "Tools assist in mapping fractured reservoirs", OIL AND GAS JOURNAL, vol. 88, no. 23, 4 June 1990 (1990-06-04), TULSA, OKLAHOMA, USA, pages 106 - 111, XP000137794 * |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1008720A1 (fr) * | 1998-12-11 | 2000-06-14 | Institut Francais Du Petrole | Procédé pour créer un modèle des écoulements de fluides dans un milieu poreux multi-couches fissuré |
FR2787219A1 (fr) | 1998-12-11 | 2000-06-16 | Inst Francais Du Petrole | Methode pour modeliser les flux de fluides dans un milieu poreux multi-couches fissure et les interactions correlatives dans un puits de production |
EP1158312A1 (fr) * | 2000-05-26 | 2001-11-28 | Institut Francais Du Petrole | Méthode pour modéliser des écoulements dans un milieu fracturé traversé par de grandes fractures |
FR2809494A1 (fr) * | 2000-05-26 | 2001-11-30 | Inst Francais Du Petrole | Methode pour modeliser des ecoulements dans un milieu fracture traverse par de grandes fractures |
FR2850773A1 (fr) * | 2003-01-31 | 2004-08-06 | Inst Francais Du Petrole | Methode pour modeliser des flux de fluides compressibles dans un milieu poreux multicouches fracture |
FR2858444A1 (fr) * | 2003-07-29 | 2005-02-04 | Inst Francais Du Petrole | Methode pour modeliser les transferts compositionnels et/ou polyphasiques entre la matrice poreuse et les fractures d'un milieu poreux multicouches |
NL1026723C2 (nl) * | 2003-07-29 | 2005-11-01 | Inst Francais Du Petrole | Werkwijze voor het modelleren van samenstellings- en/of veelfaseoverdracht tussen het poreuze medium en de breuken van een veellagig poreus medium. |
US7373288B2 (en) | 2003-07-29 | 2008-05-13 | Institut Francais Du Petrole | Method for modelling compositional and/or multiphase transfers between the porous matrix and the fractures of a multilayer porous medium |
FR2923930A1 (fr) * | 2007-11-19 | 2009-05-22 | Inst Francais Du Petrole | Methode de modelisation d'un milieu geologique poreux traverse par un reseau de fractures. |
WO2009098366A2 (fr) * | 2007-11-19 | 2009-08-13 | Ifp | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures |
WO2009098366A3 (fr) * | 2007-11-19 | 2010-03-25 | Ifp | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures |
US8688424B2 (en) | 2007-11-19 | 2014-04-01 | Ifp | Method of modelling a porous geologic medium traversed by a network of fractures |
EP2072752A1 (fr) | 2007-12-20 | 2009-06-24 | Ifp | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures |
FR2925726A1 (fr) * | 2007-12-20 | 2009-06-26 | Inst Francais Du Petrole | Methode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'echange geologique et transitoire entre blocs matriciels et fractures |
US8386225B2 (en) | 2007-12-20 | 2013-02-26 | Ifp | Method of optimizing the development of a fluid reservoir by taking into account a geologic and transient exchange term between matrix blocks and fractures |
EP2453106A1 (fr) | 2010-11-10 | 2012-05-16 | IFP Energies Nouvelles | Methode pour caracteriser le réseau de fractures d'un gisement fracture et méthode pour l'exploiter |
US8983818B2 (en) | 2010-11-10 | 2015-03-17 | IFP Energies Nouvelles | Method for characterizing the fracture network of a fractured reservoir and method for developing it |
EP2530493A1 (fr) | 2011-06-01 | 2012-12-05 | IFP Energies Nouvelles | Méthode pour construire un maillage d'un réseau de fractures à partir de diagramme de voronoi |
EP2581767A1 (fr) | 2011-10-12 | 2013-04-17 | IFP Energies Nouvelles | Méthode pour construire un maillage d'un réservoir fracturé avec un nombre limité de noeuds dans le milieu matrice |
US9665537B2 (en) | 2011-10-12 | 2017-05-30 | IFP Energies Nouvelles | Method for generating a fractured reservoir mesh with a limited number of nodes in the matrix medium |
EP3144468A1 (fr) | 2015-09-15 | 2017-03-22 | IFP Energies nouvelles | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter |
US10288544B2 (en) | 2015-09-15 | 2019-05-14 | IFP Energies Nouvelles | Method for characterizing the fracture network of a fractured reservoir and method for exploiting it |
EP3181804A1 (fr) | 2015-12-17 | 2017-06-21 | IFP Énergies nouvelles | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures |
EP3199749A1 (fr) | 2016-01-26 | 2017-08-02 | IFP Énergies nouvelles | Procede pour l'exploitation d'un gisement de fluide traverse par des fractures au moyen d'une simulation d'ecoulement basee sur un flux d'echange et un facteur correctif |
US10872182B2 (en) | 2016-01-26 | 2020-12-22 | IFP Energies Nouvelles | Method for the development of a fluid deposit traversed by fractures by means of a flow simulation based on an exchange flow and a corrective factor |
EP4112875A1 (fr) | 2021-07-01 | 2023-01-04 | IFP Energies nouvelles | Procédé pour simuler les écoulements d'un fluide dans une formation souterraine comportant un réseau de fractures |
FR3124869A1 (fr) | 2021-07-01 | 2023-01-06 | IFP Energies Nouvelles | Procédé pour simuler les écoulements d'un fluide dans une formation souterraine comportant un réseau de fractures |
FR3127250A1 (fr) | 2021-09-21 | 2023-03-24 | IFP Energies Nouvelles | Procédé pour déterminer la viscosité apparente d’une mousse à injecter en vue de la récupération d’une phase organique présente dans une formation souterraine |
Also Published As
Publication number | Publication date |
---|---|
NL1007916C2 (nl) | 1998-07-01 |
NO976129L (no) | 1998-07-01 |
GB9727231D0 (en) | 1998-02-25 |
US6064944A (en) | 2000-05-16 |
GB2322949B (en) | 2002-02-20 |
NO976129D0 (no) | 1997-12-29 |
NO326402B1 (no) | 2008-11-24 |
GB2322949A (en) | 1998-09-09 |
FR2757957B1 (fr) | 1999-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2757957A1 (fr) | Methode pour simplifier la modelisation d'un milieu geologique poreux traverse par un reseau irregulier de fractures | |
EP1158312B1 (fr) | Méthode pour modéliser des écoulements dans un milieu fracturé traversé par de grandes fractures | |
EP2212722B1 (fr) | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures | |
EP2072752B1 (fr) | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures | |
EP1008720B1 (fr) | Procédé pour créer un modèle des écoulements de fluides dans un milieu poreux multi-couches fissuré | |
EP3144468B1 (fr) | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter | |
EP2037080B1 (fr) | Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité | |
EP2253797B1 (fr) | Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide | |
EP2530493B1 (fr) | Méthode pour construire un maillage d'un réseau de fractures à partir de diagramme de voronoi | |
FR2757947A1 (fr) | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches | |
EP2581767A1 (fr) | Méthode pour construire un maillage d'un réservoir fracturé avec un nombre limité de noeuds dans le milieu matrice | |
EP1837683B1 (fr) | Utilisation de l'état stationnaire d'automates cellulaires échangeant de l'énergie pour modéliser des architectures sédimentaires | |
EP3181804A1 (fr) | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures | |
FR2967200A1 (fr) | Methode pour caracteriser le reseau de fractures d'un gisement fracture et methode pour l'exploiter | |
EP3199749A1 (fr) | Procede pour l'exploitation d'un gisement de fluide traverse par des fractures au moyen d'une simulation d'ecoulement basee sur un flux d'echange et un facteur correctif | |
EP2813668B1 (fr) | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures | |
FR3000579A1 (fr) | Fractures a segments multiples | |
EP2365359B1 (fr) | Méthode de calage d'historique d'un modèle géologique | |
EP4112875B1 (fr) | Procédé pour simuler les écoulements d'un fluide dans une formation souterraine comportant un réseau de fractures | |
EP3502753B1 (fr) | Procede de modelisation d'un bassin sedimentaire | |
CN111027780B (zh) | 油势场模拟方法及装置 | |
EP3502751B1 (fr) | Procédé de modelisation d'un bassin sédimentaire | |
Lopez et al. | Investigation into the reliability of predictive pore-scale modeling for siliciclastic reservoir rocks | |
Tilke et al. | Characterizing flow in oil reservoir rock using smooth particle hydrodynamics | |
FR2837592A1 (fr) | Methode pour modeliser des flux de fluides dans un milieu poreux multicouches traverse par un reseau de fractures inegalement reparties |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CD | Change of name or company name | ||
PLFP | Fee payment |
Year of fee payment: 20 |