FR2757947A1 - Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches - Google Patents
Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches Download PDFInfo
- Publication number
- FR2757947A1 FR2757947A1 FR9616330A FR9616330A FR2757947A1 FR 2757947 A1 FR2757947 A1 FR 2757947A1 FR 9616330 A FR9616330 A FR 9616330A FR 9616330 A FR9616330 A FR 9616330A FR 2757947 A1 FR2757947 A1 FR 2757947A1
- Authority
- FR
- France
- Prior art keywords
- crack
- network
- cracks
- permeability
- medium
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000012530 fluid Substances 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims description 5
- 238000012512 characterization method Methods 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract description 3
- 239000011159 matrix material Substances 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 9
- 210000004027 cell Anatomy 0.000 description 9
- 238000004088 simulation Methods 0.000 description 7
- 238000005336 cracking Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- RVRCFVVLDHTFFA-UHFFFAOYSA-N heptasodium;tungsten;nonatriacontahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W] RVRCFVVLDHTFFA-UHFFFAOYSA-N 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000010200 validation analysis Methods 0.000 description 3
- 235000002767 Daucus carota Nutrition 0.000 description 2
- 244000000626 Daucus carota Species 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000010420 art technique Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000001739 density measurement Methods 0.000 description 1
- 210000005045 desmin Anatomy 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
Méthode de détermination de la perméabilité de fissure équivalente d'un réseau de fissures dans un milieu souterrain multicouche à partir d'une représentation connue de ce réseau. La perméabilité de fissure équivalente d'un réseau de fissures dans un milieu multicouche souterrain est déterminée en discrétisant, à l'aide d'une procédure spécifique, chaque fissure (F) du réseau de fissures en éléments de fissure (R) (tels que des rectangles par exemple) et en définissant des noeuds N représentant des éléments de fissure interconnectés dans chaque couche du milieu et déterminant des écoulements de fluides (par exemple des écoulements en régime stable) à travers le réseau discrétisé en imposant des conditions de pression limites et des transmissivités de fluide à chaque couple de noeuds voisins. La méthode permet de relier de manière systématique des modèles de caractérisation de réservoir fissuré à des simulateurs double porosité en vue de réaliser une modélisation plus réaliste d'une structure géologique souterraine fissurée. Application par exemple en production pétrolière par des ingénieurs de gisement en vue d'obtenir des prédictions d'écoulement fiables.
Description
L'invention concerne une niéthode destinée à déterminer la perméabilité de fissure équivalente d'un réseau fisurè dans un milieu multicouche souterrain t'issOk. permettant de réaliser une modélisation plus réaliste d'une structure géologique souterraine fissurée.
La méthode peut être mise en oeuvle par exemple par des ingénieurs de gisement en ue d'obtenir des prédictions d'écoulement d'huile 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 de fissure représentant moins de I ic de l'huile en place et hautement conducteur. Le milieu de fissure lui-même peut être complexe, avec différents ensembles de fissures caractérisés par leur densité, longueur, orientåtion, 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 du réseau de fissures au sein de simulateurs d'écoulement de réservoir 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és par un système orthogonal de fissures uniformes continues orientées dans la direction de l'un des trois principaux sens de l'écoulement. L'écoulement des fluides à l'échelle du réservoir effectue à travers le milieu de fissure seulement et des échanges de fluide surviennent localement entre les fissures et les blocs matriciels.
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és par un système orthogonal de fissures uniformes continues orientées dans la direction de l'un des trois principaux sens de l'écoulement. L'écoulement des fluides à l'échelle du réservoir effectue à travers le milieu de fissure seulement et des échanges de fluide 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élioratlons spécifiques relatives à la modélisation des échanges de fluide matrice-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é'rier 1953) 42-54:
Quandalle, P. et aI : "Typical Features of a New Multipurpose Reservoir Simulatol-".
Quandalle, P. et aI : "Typical Features of a New Multipurpose Reservoir Simulatol-".
SPE 16007. présenté au 9ème Symposium sur la simulation de réservoir 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éservoir de Houston, Texas, 6-8 février 1989.
Porosity Reservoirs", SPE 18427 présenté au Symposium SPE sur la simulation de réservoir 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 propriétés pétrophysiques de base des fissures et des matrices, ainsi que 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.
la perméabilité du réseau de fissures contenu dans la cellule, c'est-à-dire la perméabilité de fissure équivalente, ne peut pas être évaluée simplement et nécessite la prise en compte de la géométrie et des propriétés du réseau de fissures.
Une méthode directe connue permet de déterminer l'écoulement en régime stable dans un réseau de fissures. Elle comprend l'utilisation de mailles fines et régulières conventionnelles discrétisant à la fois les fissures et les blocs matriciels du volume de roche fissurée parallélépipédique considéré. Pour différentes raisons. cette méthode connue ne fournit pas de résultats fiables, sauf si le volume de roche fissurée est discrétisé en utilisant une maille comportant un nombre extrêmement élevé de cellules. ce qui nécessitc des moyens de calcul très importants.
D'autres modèles spécifiques servant à calculer les perméabilités équivalentes de réseaux fissurés en 2D ou 3D sont également connus. par exemple par
- Odling. N.E. ''Perme;lbility of Natural and Simulated Fracture Patterns". Structural and Tectonic Modelling and its Application to Petroleum Geology. NPF Special
Publication 1,365-380. Elsevier. Norwegian Petroleum Society (N'PF) 199':
- Long. J.C.S. et al : "A Modes for Steady Fluid Flou in Random Three-Dimensional Networks of Disc-Shaped Fractures", Water Resources Research (août 1985 > . vol.9 I. n 8. 1105-1115:
- Cacas. M.C. et al : ",NIodeling Fracture Flow with a Stochastic Discrete Fracture
Network : Calibration and Validation. 1. The Flow Model". Water Resources Research (mars 1990) vol.26, n 3;
- Billaux, D. : Hydrogéologie des milieux fracturés. Géométrie, connectivité et comportement hydraulique". Thèse de doctorat préentée à l'école Nationale Supérieure des Mines de Paris, document du BRGM n 186, Editions du BRGM, 1990;
- Robinson, P.C. "Connectivity, Flow and Transport in networks Models of
Fractured Media", Thèse de doctorat, St Catherine's College, Oxford University, Réf.
- Odling. N.E. ''Perme;lbility of Natural and Simulated Fracture Patterns". Structural and Tectonic Modelling and its Application to Petroleum Geology. NPF Special
Publication 1,365-380. Elsevier. Norwegian Petroleum Society (N'PF) 199':
- Long. J.C.S. et al : "A Modes for Steady Fluid Flou in Random Three-Dimensional Networks of Disc-Shaped Fractures", Water Resources Research (août 1985 > . vol.9 I. n 8. 1105-1115:
- Cacas. M.C. et al : ",NIodeling Fracture Flow with a Stochastic Discrete Fracture
Network : Calibration and Validation. 1. The Flow Model". Water Resources Research (mars 1990) vol.26, n 3;
- Billaux, D. : Hydrogéologie des milieux fracturés. Géométrie, connectivité et comportement hydraulique". Thèse de doctorat préentée à l'école Nationale Supérieure des Mines de Paris, document du BRGM n 186, Editions du BRGM, 1990;
- Robinson, P.C. "Connectivity, Flow and Transport in networks Models of
Fractured Media", Thèse de doctorat, St Catherine's College, Oxford University, Réf.
TP1072, mai 1984.
La présente invention concerne une méthode de détermination de la perméabilité de fissure équivalente d'un réseau de fissures dans un milieu multicouche souterrain.
La méthode est caractérisée en ce qu'elle comporte les étapes suivantes
- discrétisation du réseau de fissures en éléments de fissure (tels que des rectangles par exemple) et définition de noeuds représentant les éléments de fissure interconnectés dans chaque couche du milieu, et
- détermination des écoulements de fluides à travers le réseau discrétisé en imposant des conditions de pression limites et des transmissivités de fluide à chaque couple de noeuds voisins.
- discrétisation du réseau de fissures en éléments de fissure (tels que des rectangles par exemple) et définition de noeuds représentant les éléments de fissure interconnectés dans chaque couche du milieu, et
- détermination des écoulements de fluides à travers le réseau discrétisé en imposant des conditions de pression limites et des transmissivités de fluide à chaque couple de noeuds voisins.
Plus précisément, la méthode est caractérisée en ce que:
- I'on divise le milieu cn une série de couches parallèles s'étendant chacune dais un plan de référence perpendiculaire à l'axe de référence et définies chacune par une coordonnée le long dudit axe,
- on divise chaque fissure en une série de rectangles limités le long dudit axe de rétëlence par deux couches adjacentes et on catalogue les rectangles en leur associant des caractères géométriques et physiques tels que coordonnées et tailles des rectangles et conductivités hydrauliques des fissures,
- on positionne des noeuds dans chaque couche pour toutes les fissures interconnectées, et
- pour tous les couples de noeuds voisins, on calcule des facteurs de transmissivité et résolution des équations d'écoulement en vue de déterminer les perméabilités équivalentes du milieu dans trois directions orthogonales.
- I'on divise le milieu cn une série de couches parallèles s'étendant chacune dais un plan de référence perpendiculaire à l'axe de référence et définies chacune par une coordonnée le long dudit axe,
- on divise chaque fissure en une série de rectangles limités le long dudit axe de rétëlence par deux couches adjacentes et on catalogue les rectangles en leur associant des caractères géométriques et physiques tels que coordonnées et tailles des rectangles et conductivités hydrauliques des fissures,
- on positionne des noeuds dans chaque couche pour toutes les fissures interconnectées, et
- pour tous les couples de noeuds voisins, on calcule des facteurs de transmissivité et résolution des équations d'écoulement en vue de déterminer les perméabilités équivalentes du milieu dans trois directions orthogonales.
Dans un mode de réalisation préféré, la perméabilité équivalente du milieu comprend ia détermination directe du tenseur d'anisotropie de perméabilité équivalente et l'étalonnage de valeurs absolues de perméabilité à partir de résultats d'essais de puits.
La méthode ainsi résumée permet de relier de manière systématique des modèles de caractérisation de réservoirs fissurés et des simulateurs double porosité en vue de réaliser une modélisation plus réaliste d'une structure géologique souterraine fissurée. La méthode peut être mise en oeuvre par exemple en production pétrolière par des ingénieurs de gisement en vue d'obtenir des prédictions d'écoulement fiables.
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 montre par exemple une image en 3D d'un réseau de fissures stochastiquement obtenue à partir d'observations et de mesures prises sur un affleurement de gres.
- la figure I montre par exemple une image en 3D d'un réseau de fissures stochastiquement obtenue à partir d'observations et de mesures prises sur un affleurement de gres.
- les figures 7. 3 montrent une fissure divisée en une série de rectangles R.
- la figure 4 montre la structure des données d'entrée spécifiant les caractères de la fissure.
- la figure 5 montre le mode préféré de discrétisation d'un plan de fissure.
- les figures 6. 7. 8 et 9 illustrent schématiquement le calcul des facteurs de transmissivité pour différentes positions de noeuds les uns par rapport aux autres ou avec une limite.
Les perméabilités équivalentes d'un réseau de fissures 3D sont déterminées ci-après en utilisant une technique numérique basée sur la méthode connue de "réseau de résistance" décrit par exemple dans le document d'Odling N.E. cité ci-avant. Dans la présente méthode, on admet que la matrice est imperméable afin qu'elle soit compatible avec l'approche double porosité. Dans les simulateurs de réservoir, les écoulements matricefissure et matrice-matrice sont en fait calculés séparément par rapport aux écoulements au sein des fissures.
On admet que le réseau de fissures en 3D considéré représente. dans un volume égal à une cellule de réservoir, la distribution réelle des fissures donnée par intégration des caractères de fissure du champ dans un modèle de caractérisation. Le principal objectif de calculs d'écoulements monophasiques dans le réseau de fissures en 3D est d'évaluer l'anisotropie de perméabilité équivalente (Kv/Kh et Ky/Kx) de la cellule de fissure considérée, qui constitue un paramètre important régissant le comportement des écoulements polyphasiques dans les réservoirs. Les valeurs de perméabilité équivalente obtenues par ces calculs pourraient en pratique être comparées aux résultats d'essais de puits en vue d'étalonner les caractères de fissures tels que les conductivités hydrauliques de fissure (ou les ouvertures hydrauliques équivalentes). qui ne peuvent être définies que de manière médiocre a priori.
De plus. les résultats dc perméabilité équivalente peuvent être utilisés en vue dc déterminer un tenseur de perméabilité dont les principales directions permettent tine orientation optimale de la maille du modèle de réservoir. Cependant. des conditions limites spécifiqucs sont necessaires pour obtenir ces informations Des contlittons limites sans écoulement latéral imposées sur les quatre faces latérales du volume parallèlépipédique étudié ne donnent pas accès aux termes non diagonaux du tenseur de perméabilité équivalente. alors que des potentiels (ou des pressions; à 'aviation linéale SUI les faces latérales permettent d'imposer la direction du gradient de potentiel au selon du milieu anisotrope et de déduire diirctement les termes de perméabilité non diagonaux des débits latéraux.
Les techniques d'intégration de données de fissuration naturelle dans des modèles de réservoirs fissurés sont bien connues. Les données de fissuration sont essentiellement dc nature géométrique et comportant des mesures de densité. de longueur. d'azimut et d'inclinaison des plans de fissure soit observés sur des affleurements. dans des galeries de mines, sur des carottes. soit déduits à partir de diagraphies. Différents ensembles de fissures peuvent être distingués et caractérisés par différentes distributions statistiques de leurs caractères de fissure. Une fois les schémas de fissuration caractérisés. on peut créer des réseaux numériques de ces ensembles de fissures en utilisant un procédé stochastique respectant les distributions statistiques des paramètres de fissures. De tels procédés sont décrits par exemple dans les brevets FR-A-2.725.814; 2.725.794 ou 2.733.du demandeur.
La méthode selon l'invention s'applique à des images de structures géologiques fissurées présentant différentes tailles ou volumes et/ou situées en différents endroits.
créées par un générateur de modèle de fissure. La figure 1 présente une telle image.
DONNEES D'ENTREE
Avant de développer les procédures recommandées en vue de déterminer les paramètres hydrauliques équivalents d'images de fissures en 3D, une étape importante consiste à définir une structure commune de données d'entrée pqur ces images afin qu'elles puissent être traitées indépendamment de l'outil de traitement utilisé pour leur création.
Avant de développer les procédures recommandées en vue de déterminer les paramètres hydrauliques équivalents d'images de fissures en 3D, une étape importante consiste à définir une structure commune de données d'entrée pqur ces images afin qu'elles puissent être traitées indépendamment de l'outil de traitement utilisé pour leur création.
Comme le montrent les figures 2. 3, on admet que les fissures F sont sensiblement verticales (c'est-à-dire perpendiculaires aux limites de la couche). Cependant. une nlcllle structure de données peut être appliquée aux fissures s'écartant légèrement de la verticale. L'image 3D est discrCtisée verticalement en respectant la stratification géologique réelle lorsque cette information est disponible. Dans le cas cintraire. on applique une discrétisation arbitraire à l'image. Chaque couche horizontale L est caractérisée par sa coordonnée verticale zL dans le système de référence de coordonnées > ON. OY. OZ).
Une série de rectangles R doit être définie pour chaque couche L. Chaque rectangle consiste en un élément de plan de fissure compris entre les limites d'une couche donnée.
Ainsi. chaque fissure naturelle consiste en une série de rectangles superposés R et on lui attribue une origine (origine de la fissure). Chaque rectangle est défini par:
- les trois coordonnées (xO, yO, zO) de l'origine du rectangle O. Pour une fissure naturelle donnée, tous les points d'origine des rectangles constitutifs sont situés sLir la même ligne verticale (ou plus fort pendage) tracée depuis l'origine de la fissure:
- les coordonnées du vecteur d'unité horizontal i (xH, yH) et du vecteur d'unité vertical j (xV, yV) définissant l'orientation du rectangle dans le système de référence de coordonnées, avec xVertical et yVertical égaux à zéro en cas de fissures verticales mais considérés comme des données d'entrée afin de pouvoir traiter des fissures non verticales;
- les deux longueurs horizontales algébriques 1- et 1+ séparant l'origine du rectangle et les deux limites latérales (verticales) de ce rectangle;
- la hauteur h du rectangle, c'est-à-dire la longueur du rectangle dans le sens j. qui est l'épaisseur de la couche si la discrétisation en sens j correspond à la géologie;
- la conductivité hydraulique c dérivée de l'application de la loi de Darcy relative à
AP l'écoulement dans les fissures (pour un gradient de pression , le débit au sein de la
ch A P fissure de hauteur h est -. , tt étant la viscosité du fluide). La conductivité c est donnée par la relation c=k.a. où k = a2/l 2 (en utilisant la représentation idéalisée clefs fissures de Poiseuille) est la perméabilité intrinsèque de la fissure et a son oLIvcltulc hydraulique équivalente. La conductivité hydraulique c est une valeur dc rêlèrence donnée pour une direction de la contrainte horizontale maximale parallèle au sens de la fisstire:
- les deux rectangles sLipérieLir et inférieur voisins l UR. LR;
- la série de fissures FS à laquelle appartient le rectangle considéré:
- l'angle d'orientation ao de la direction de la contrainte horizontale maximale prise à partir de l'axe (OX) dans le système de référence de coordonnées:
- pour chaque série de fissures, une table de corrélation corrélant I ) l'angle compris entre la direction de la contrainte horizontale maximale et la direction de la fissure (azimut) avec 2) la conductivité hydraulique c ou l'ouverture hydraulique équivalente a précédemment définies. Les termes "horizontal" et "vertical" utilisés dans ce contexte se rapportent aux directions respectivement parallèles et perpendiculaires aux limites de couches considérées ici comme étant horizontales. Les limites de couche discrétisent les plans de fissure dans le sens "vertical" Il convient de noter que les données d'entrée mentionnées ci-dessus I ) conviennent pour tous les outils logiciels existants utilisés pour caractériser et générer des fissures et 9) qu'elles pourraient ête utilisées pour discrétiser un réseau de fissures légèrement non verticales, c'est-à-dire qui ne sont pas perpendiculaires aux limites de couche.
- les trois coordonnées (xO, yO, zO) de l'origine du rectangle O. Pour une fissure naturelle donnée, tous les points d'origine des rectangles constitutifs sont situés sLir la même ligne verticale (ou plus fort pendage) tracée depuis l'origine de la fissure:
- les coordonnées du vecteur d'unité horizontal i (xH, yH) et du vecteur d'unité vertical j (xV, yV) définissant l'orientation du rectangle dans le système de référence de coordonnées, avec xVertical et yVertical égaux à zéro en cas de fissures verticales mais considérés comme des données d'entrée afin de pouvoir traiter des fissures non verticales;
- les deux longueurs horizontales algébriques 1- et 1+ séparant l'origine du rectangle et les deux limites latérales (verticales) de ce rectangle;
- la hauteur h du rectangle, c'est-à-dire la longueur du rectangle dans le sens j. qui est l'épaisseur de la couche si la discrétisation en sens j correspond à la géologie;
- la conductivité hydraulique c dérivée de l'application de la loi de Darcy relative à
AP l'écoulement dans les fissures (pour un gradient de pression , le débit au sein de la
ch A P fissure de hauteur h est -. , tt étant la viscosité du fluide). La conductivité c est donnée par la relation c=k.a. où k = a2/l 2 (en utilisant la représentation idéalisée clefs fissures de Poiseuille) est la perméabilité intrinsèque de la fissure et a son oLIvcltulc hydraulique équivalente. La conductivité hydraulique c est une valeur dc rêlèrence donnée pour une direction de la contrainte horizontale maximale parallèle au sens de la fisstire:
- les deux rectangles sLipérieLir et inférieur voisins l UR. LR;
- la série de fissures FS à laquelle appartient le rectangle considéré:
- l'angle d'orientation ao de la direction de la contrainte horizontale maximale prise à partir de l'axe (OX) dans le système de référence de coordonnées:
- pour chaque série de fissures, une table de corrélation corrélant I ) l'angle compris entre la direction de la contrainte horizontale maximale et la direction de la fissure (azimut) avec 2) la conductivité hydraulique c ou l'ouverture hydraulique équivalente a précédemment définies. Les termes "horizontal" et "vertical" utilisés dans ce contexte se rapportent aux directions respectivement parallèles et perpendiculaires aux limites de couches considérées ici comme étant horizontales. Les limites de couche discrétisent les plans de fissure dans le sens "vertical" Il convient de noter que les données d'entrée mentionnées ci-dessus I ) conviennent pour tous les outils logiciels existants utilisés pour caractériser et générer des fissures et 9) qu'elles pourraient ête utilisées pour discrétiser un réseau de fissures légèrement non verticales, c'est-à-dire qui ne sont pas perpendiculaires aux limites de couche.
MODES OPERATOIRES
Les modes opératoires et les tests de validation de la méthode permettant de calculer l'anisotropie de perméabilité du réseau de fissures dans son ensemble sont décrits ci-après à partir de l'image en 3D ainsi codifiée. La procédure numérique destinée à calculer les perméabilités équivalentes d'un réseau de fissures 3D est également décrite.
Les modes opératoires et les tests de validation de la méthode permettant de calculer l'anisotropie de perméabilité du réseau de fissures dans son ensemble sont décrits ci-après à partir de l'image en 3D ainsi codifiée. La procédure numérique destinée à calculer les perméabilités équivalentes d'un réseau de fissures 3D est également décrite.
Le problème consiste à trouver la distribution des débits dans le réseau pour les conditions frontières suivantes aux limites du volume parallélépipédique étudié, c'est-à- dire les pressions fixes imposées sur deux faces opposées et les pressions variant de manière linéaire sur les quatre faces latérales (entre les valeurs inlposees sur les deux autres faces)
Les principales étapes sont résumées ci-dessous
1) Discrétisation du réseau
En utilisant les définitions données pour la structure des données d'entrée. le résealLI cie tissure est discrétlsé sous la forme d'une série de "noeuds" N. chaque noeud étant placé au centre des segments d'intersection IS de deux rectangles R (c'est-à-dire de deux plans de fissure au sein d'une couche donnée). Comme le montre la figure 5. des noeuds supplémentaires AN sont placés au-dessus et au-dessous des noeuds précédents N afin de représenter d'autres rectangles discrétisant les fissures et de minimiser les longueurs d'écoulement au sein d'une fissure donnée. Sur la figure 5, BL est une limite latérale de deux cellules de fissure voisines.
Les principales étapes sont résumées ci-dessous
1) Discrétisation du réseau
En utilisant les définitions données pour la structure des données d'entrée. le résealLI cie tissure est discrétlsé sous la forme d'une série de "noeuds" N. chaque noeud étant placé au centre des segments d'intersection IS de deux rectangles R (c'est-à-dire de deux plans de fissure au sein d'une couche donnée). Comme le montre la figure 5. des noeuds supplémentaires AN sont placés au-dessus et au-dessous des noeuds précédents N afin de représenter d'autres rectangles discrétisant les fissures et de minimiser les longueurs d'écoulement au sein d'une fissure donnée. Sur la figure 5, BL est une limite latérale de deux cellules de fissure voisines.
Une fois le réseau discrétisé, une procédure de tri est appliquée à ce réseau de fissures en vue d'éliminer les noeuds isolés ou des groupes de noeuds sans lien avec l'une des limites latérales FL du volume 3D étudié, car les fissures ainsi "triées" ne contribuent pas au transport du fluide et peuvent entraver les procédures de résolution mises en oeuvre en vue de trouver les pressions aux noeuds de fissures lors d'un écoulement en régime stable à travers le réseau.
2) Calcul des transmissivités
Un facteur de transmissivité T est calculé pour chaque paire de noeuds associés en utilisant la relation
ch kah
T = - = -- i i
où c est la conductivité hydraulique de la fissure, k la perméabilité intrinsèque de la fissure, a l'ouverture de la fissure, h la hauteur de la fissure et I la distance entre deux noeuds d'une fissure.
Un facteur de transmissivité T est calculé pour chaque paire de noeuds associés en utilisant la relation
ch kah
T = - = -- i i
où c est la conductivité hydraulique de la fissure, k la perméabilité intrinsèque de la fissure, a l'ouverture de la fissure, h la hauteur de la fissure et I la distance entre deux noeuds d'une fissure.
Différentes situations doivent être envisagées selon la position respective des deux noeuds. Pour des noeuds contenus au sein de la même couche (fig.6 > . le facteurs' dc transmissivité horizontale T est obtenu directement sous la forme de la distance (11+12 > séparant les deux noeuds dans le sens de l'écoulement (fi.7). Pour des noeuds SitLiés dans deux couches difi;e'rentes (Fie ())* le facteur de transrllissivité horizonLale est la somme arithmétique des facteurs de transmissivité (T'+T") relatifs aux deux éléments du plan de fissure des cellules de fissure superposées. II comporte une longueur d'écoulement égale à la demi-somme des deux épaisseurs de couche hl et hl. Pour les noeuds supplémentaires tels que dètinîs précédemment. reliés pal un seul plan de fissure.
un seul facteur de transmissivité est calculé pour cet élément du plan de fissure.
Le facteur de transmissivité T entre un noeud et une limite du volume 7D étudié est exprimé comme dans le cas de deux noeuds. selon les deux cas suivants.
Pour une limite verticale latérale, le facteur de transmissivité T peut être exprimé directement pour un seul élément du plan de fissure (fig.8), et sous la forme de la somme de deux transmissivités si deux plans de fissure relient le noeud et la limite.
Pour une limite inférieure ou supérieure horizontale. le facteur de transmissivité verticale peut être exprimé en considérant une longueur d'écoulement égale à la moitié de l'épaisseur de la couche (fig.9).
3) Equations d'écoulement
En régime stable, un écoulement monophasique incompressible à travers le réseau de fissure est déterminé en résolvant un série de n équations, une pour chaque noeud, comme il est connu dans l'art. Chaque équation exprime le fait que le débit total est nul à chaque noeud de fissure. Pour calculer un tenseur de perméabilité, on considère qu'une pression constante est imposée à chacune des limites aval et amont. Une pression variant de manière linéaire en fonction de la position entre les limites aval et amont est imposée.
En régime stable, un écoulement monophasique incompressible à travers le réseau de fissure est déterminé en résolvant un série de n équations, une pour chaque noeud, comme il est connu dans l'art. Chaque équation exprime le fait que le débit total est nul à chaque noeud de fissure. Pour calculer un tenseur de perméabilité, on considère qu'une pression constante est imposée à chacune des limites aval et amont. Une pression variant de manière linéaire en fonction de la position entre les limites aval et amont est imposée.
La matrice de perméabilité équivalente (Kij) déterminée précédemment est diagonalisée pour calculer les directions principales de l'écoulement avec les perméabilités équivalentes respectives dans ces directions.
En pratique, le problème se limite souvent à la découverte des principales directions horizontales de l'écoulement U et V dans la mesure où la direction perpendiculaire aux limites de couche t ènéralement verticales) est toujours prise comme l'axe des z. Dans un tel cas, seuls les termes extra-diagonaux Kxy et Kyx doivent être calculés: ils peuvent être obtenus avec les conditions limites mixtes suivantes
- les écoulements horizontaux sont calculés avec des faces Ints-rieui-es et supérieures mperméahles et des pressions à variation linéaire sur les faces verticale parallèles au .sens de l'écoulement.
- les écoulements horizontaux sont calculés avec des faces Ints-rieui-es et supérieures mperméahles et des pressions à variation linéaire sur les faces verticale parallèles au .sens de l'écoulement.
- I'écoulement vertical est calculé avec toutes les faces latérales imperméables.
Ainsi. un tenseur de perméabilité simplifié est obtenu. à partir duquel on peut facilement déduire les principales directions horizontales de l'écoulement U et V
Validation
La méthode a été validée avec succès par rapport aux calculs d'écoulement monophasique de référence mentionnés ci-avant, réalisés avec un simulateur de gisement conventionnel. Les calculs de référence ont été obtenus sur de fines mailles régulières discrétisant les fissures, ainsi que les blocs matriciels du volume de roche fissurée parallélépipédique considéré. Pour un sens d'écoulement donné, une pression d'injection et une pression de production fixes ont été imposées sur les faces d'entrée et de sortie et le débit résultant a été calculé pour des conditions exemptes d'écoulement latéral.
La méthode a été validée avec succès par rapport aux calculs d'écoulement monophasique de référence mentionnés ci-avant, réalisés avec un simulateur de gisement conventionnel. Les calculs de référence ont été obtenus sur de fines mailles régulières discrétisant les fissures, ainsi que les blocs matriciels du volume de roche fissurée parallélépipédique considéré. Pour un sens d'écoulement donné, une pression d'injection et une pression de production fixes ont été imposées sur les faces d'entrée et de sortie et le débit résultant a été calculé pour des conditions exemptes d'écoulement latéral.
On a réalisé trois étapes pour valider le calcul de
- la perméabilité verticale équivalente d'un volume de roche traversé par une seule fissure, cette dernière étant représentée par plusieurs noeuds correspondant aux intersections avec de petites fissures isolées;
- les perméabilités horizontales équivalentes (dans une géométrie d'écoulement en 2D > et les principaux sens d'écoulement:
- les perméabilités équivalentes et l'anisotropie de perméabilité dans Lin seul éseaLi présentant une géométrie d'écoulement en 3D.
- la perméabilité verticale équivalente d'un volume de roche traversé par une seule fissure, cette dernière étant représentée par plusieurs noeuds correspondant aux intersections avec de petites fissures isolées;
- les perméabilités horizontales équivalentes (dans une géométrie d'écoulement en 2D > et les principaux sens d'écoulement:
- les perméabilités équivalentes et l'anisotropie de perméabilité dans Lin seul éseaLi présentant une géométrie d'écoulement en 3D.
Les résultats obtenus pour la troisième étape (pour une géométrie d'écoulement en 3D > sont donnés dans le tableau ci-après. Pour les directions horizontales d'écoulement.
une solution analytique de référence peut également être calculée puisque la géométrie de l'écoulement est un écoulement 2D dans ces directions (la géométrie d'écoulement en 3D concerne la direction z).
<tb>
<SEP> Perméabilités <SEP> Simulation <SEP> METHODE <SEP> Solution
<tb> équivalentes <SEP> (md) <SEP> MAILLE <SEP> FINE <SEP> PRESENTEE <SEP> ANALYTIQUE
<tb> <SEP> Kx <SEP> 0,119 <SEP> 0,120 <SEP> 0,120
<tb> <SEP> Ky <SEP> 0,2240,227 <SEP> 0,226 <SEP>
<tb> <SEP> Kz <SEP> 0,255 <SEP> 0,267
<tb> <SEP> Anisotropie
<tb> <SEP> Kz/(KxKy)0-5 <SEP> 1,56 <SEP> 1,62
<tb>
Il apparaît clairement que les résultats obtenus avec la méthode présentée sont très proches des valeurs correspondantes obtenues avec la solution analytique et la simulation à maille fine pour les directions X et Y.
<tb> équivalentes <SEP> (md) <SEP> MAILLE <SEP> FINE <SEP> PRESENTEE <SEP> ANALYTIQUE
<tb> <SEP> Kx <SEP> 0,119 <SEP> 0,120 <SEP> 0,120
<tb> <SEP> Ky <SEP> 0,2240,227 <SEP> 0,226 <SEP>
<tb> <SEP> Kz <SEP> 0,255 <SEP> 0,267
<tb> <SEP> Anisotropie
<tb> <SEP> Kz/(KxKy)0-5 <SEP> 1,56 <SEP> 1,62
<tb>
Il apparaît clairement que les résultats obtenus avec la méthode présentée sont très proches des valeurs correspondantes obtenues avec la solution analytique et la simulation à maille fine pour les directions X et Y.
De plus, la différence entre les valeurs de perméabilité équivalente verticale concernant l'écoulement 3D reste acceptable. Ainsi, le rapport d'anisotropie, égal à 1,6, est prédit de manière satisfaisante par la méthode avec un nombre très limité de cellules.
La méthode selon l'invention, qui fournit une représentation facilement transposable d'un réseau de fissures naturel, se prête bien au calcul d'écoulements au sein de fissures.
Elle peut également être utile pour améliorer l'image d'origine du réseau de fissures. Une telle image est en fait obtenue à partir d'un générateur de fissure stochastique utilisant comme données d'entrée les résultats de l'intégration de données de fissuration archivées dans un modèle de caractértsation de fissure comme il est décrit dans les brevets FR-A 2.725.814: 2.725.794 ou ' 733.073 du demandeur déja mentionnés. De telles Images.
une fois discrétisées par la méthode selon l'invention. peuvent être facilement modifiées
afin d'être adaptées aux règles géologiques. Par exemple. l'interruption systématiqLie d'une fissure donnée contre une autre série de fissures peut être prise en compte dans l'image d'origine en annulant les éléments du plan de fissure d'une série donnée 'étendant au-delà des fissures entrecoupées de l'autre série.
afin d'être adaptées aux règles géologiques. Par exemple. l'interruption systématiqLie d'une fissure donnée contre une autre série de fissures peut être prise en compte dans l'image d'origine en annulant les éléments du plan de fissure d'une série donnée 'étendant au-delà des fissures entrecoupées de l'autre série.
Claims (3)
- REVENDICATiONS- détermination d'écoulements de fluides à travers le réseau de fissures en imposant des conditions limites de pression et de transmissivités de fluides pour chaque couple de noeuds voisins.- discrétisation de chaque fissure (F) du réseau de fissures en éléments de fissure et de définition de noeuds (N) représentant des éléments de fissure interconnectés dans chaque couche du milieu, et1) Méthode permettant de déterminer la perméabilité de fissure êqui'alente d'lll réseau de fissures dans un milieu multicouche souterrain fissuré. à partir diane représentation prédéterminée dudit réseau. comportant les étapes de
- 2) Méthode selon la revendication I, caractérisée en ce que:- I'on divise le milieu en une série de couches parallèles s'étendant chacune dans un plan de référence perpendiculaire à l'axe de référence et définies chacune par une coordonnée le long dudit axe,- on divise chaque fissure en une série de rectangles limités le long dudit axe de référence par deux couches adjacentes et on catalogue les rectangles en leur associant des caractères géométriques et physiques tels que coordonnées et tailles des rectangles et conductivités hydrauliques des fissures,- on positionne des noeuds dans chaque couche pour toutes les fissures interconnectées, et- pour tous les couples de noeuds voisins, on calcule des facteurs de transmissivité et résolution des équations d'écoulement en vue de déterminer les perméabilités équivalentes du milieu dans trois directions orthogonales.
- 3) Méthode selon l'une des revendications précédentes, caractérisée en ce que la perméabilité équivalente du milieu comprend la détermination directe d'un tenseur de perméabilité équivalent et l'étalonnage de valeurs absolues de perméabilité à partir de résultats d'essais de puits.A) Méthode selon l'une des revendications précédentes. caractérisée cn ce que lesdits caractères géométriques et phvsiques sont les coordonnées et les dimensions des rectangles et les conductivités hydrauliques des fissures.
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9616330A FR2757947B1 (fr) | 1996-12-30 | 1996-12-30 | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches |
GB9727227A GB2322948B (en) | 1996-12-30 | 1997-12-23 | Method of determining the equivalent permeability of a fracture network in a sub-surface, multi-layered medium |
NL1007917A NL1007917C2 (nl) | 1996-12-30 | 1997-12-29 | Werkwijze voor het bepalen van de equivalente breukpermeabiliteit van een breuknetwerk in een gebroken meerlaags ondergronds milieu. |
NO19976127A NO326401B1 (no) | 1996-12-30 | 1997-12-29 | Fremgangsmate for a bestemme ekvivalent fraktur-permeabilitet for et fraktur-nettverk i et flerlags undergrunnsmedium |
US09/000,972 US6023656A (en) | 1996-12-30 | 1997-12-30 | Method for determining the equivalent fracture permeability of a fracture network in a subsurface multi-layered medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9616330A FR2757947B1 (fr) | 1996-12-30 | 1996-12-30 | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2757947A1 true FR2757947A1 (fr) | 1998-07-03 |
FR2757947B1 FR2757947B1 (fr) | 1999-01-29 |
Family
ID=9499398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR9616330A Expired - Fee Related FR2757947B1 (fr) | 1996-12-30 | 1996-12-30 | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches |
Country Status (5)
Country | Link |
---|---|
US (1) | US6023656A (fr) |
FR (1) | FR2757947B1 (fr) |
GB (1) | GB2322948B (fr) |
NL (1) | NL1007917C2 (fr) |
NO (1) | NO326401B1 (fr) |
Cited By (13)
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 |
FR2918179A1 (fr) * | 2007-06-29 | 2009-01-02 | Inst Francais Du Petrole | Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite |
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 |
CN108756859A (zh) * | 2018-05-03 | 2018-11-06 | 中国石油天然气股份有限公司 | 气藏多层合采边水水侵机理的模拟实验装置及方法 |
CN109932129A (zh) * | 2017-12-15 | 2019-06-25 | 核工业北京地质研究院 | 一种表面接触式水压力传感器的快速标定方法 |
CN113533157A (zh) * | 2021-07-02 | 2021-10-22 | 浙江大学 | 一种用于可视化实验的变开度可拆卸裂隙装置 |
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 (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6853921B2 (en) | 1999-07-20 | 2005-02-08 | Halliburton Energy Services, Inc. | System and method for real time reservoir management |
US7283941B2 (en) * | 2001-11-13 | 2007-10-16 | Swanson Consulting Services, Inc. | Computer system and method for modeling fluid depletion |
US7584165B2 (en) * | 2003-01-30 | 2009-09-01 | Landmark Graphics Corporation | Support apparatus, method and system for real time operations and maintenance |
CN101371164A (zh) * | 2006-01-20 | 2009-02-18 | 界标制图有限公司 | 动态生产系统管理 |
US7925482B2 (en) * | 2006-10-13 | 2011-04-12 | Object Reservoir, Inc. | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
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 |
FR2930350B1 (fr) * | 2008-04-17 | 2011-07-15 | Inst Francais Du Petrole | Procede pour rechercher des hydrocarbures dans un bassin geologiquement complexe,au moyen d'une modelisation de bassin |
EP2411912A2 (fr) * | 2009-03-24 | 2012-02-01 | Chevron U.S.A. Inc. | Système et procédé pour identifier des fractures dans un réservoir souterrain |
US9085957B2 (en) | 2009-10-07 | 2015-07-21 | Exxonmobil Upstream Research Company | Discretized physics-based models and simulations of subterranean regions, and methods for creating and using the same |
CN101871876B (zh) * | 2010-06-09 | 2011-08-31 | 中国矿业大学 | 可视化多功能裂隙渗流模拟实验台 |
WO2012021292A1 (fr) * | 2010-08-09 | 2012-02-16 | Conocophillips Company | Procédé de mise à l'échelle supérieure de réservoir avec une transmissibilité préservée |
US9262561B2 (en) * | 2013-02-15 | 2016-02-16 | Chevron U.S.A. Inc. | Modeling of non-newtonian fluids in subterranean reservoirs |
US20140290937A1 (en) * | 2013-03-27 | 2014-10-02 | Baker Hughes Incorporated | Shale fracture flow simulation apparatus |
FR3005988B1 (fr) | 2013-05-21 | 2015-05-15 | IFP Energies Nouvelles | Procede d'exploitation d'un milieu fracture a partir d'un modele de reservoir cale pour des puits choisis au moyen d'un modele de transmissivite equivalente |
CN103410502B (zh) * | 2013-08-05 | 2016-03-16 | 西南石油大学 | 一种网络状缝洞型油藏的三维渗透率场的获取方法 |
US10408028B2 (en) * | 2013-10-15 | 2019-09-10 | Halliburton Energy Services, Inc. | Optimization of engine emissions from equipment used in well site operations |
US11125912B2 (en) * | 2013-11-25 | 2021-09-21 | Schlumberger Technology Corporation | Geologic feature splitting |
CN105114065A (zh) * | 2015-08-06 | 2015-12-02 | 中国石油天然气股份有限公司 | 水力压裂裂缝延伸过程的模拟方法 |
US10337294B2 (en) | 2015-10-01 | 2019-07-02 | Baker Hughes, A Ge Company, Llc | Reservoir permeability upscaling |
US10329904B2 (en) | 2016-03-18 | 2019-06-25 | Baker Hughes, A Ge Company, Llc | Determining the robustness of discrete fracture network permeability estimates |
CN106050224B (zh) * | 2016-05-26 | 2019-05-17 | 中国石油集团渤海钻探工程有限公司 | 一种基于吸水剖面的油藏非均质性试井分析方法 |
US10612356B2 (en) | 2017-03-01 | 2020-04-07 | Proptester, Inc. | Fracture fluid and proppant transport testing systems and methods of using same |
US11365626B2 (en) | 2017-03-01 | 2022-06-21 | Proptester, Inc. | Fluid flow testing apparatus and methods |
EP3394379B1 (fr) | 2017-03-10 | 2020-02-05 | Allmetal, Inc. | Construction d'entretoise pour vitrage isolant |
CN108804750B (zh) * | 2018-04-19 | 2021-11-12 | 中国石油化工股份有限公司 | 适用于油藏数值模拟的三维矢量渗透率的获取方法 |
CN111798328A (zh) * | 2019-03-22 | 2020-10-20 | 中国石油化工股份有限公司 | 一种计算低渗透油藏五点井网瞬时产量的方法 |
CN110261282B (zh) * | 2019-07-23 | 2021-09-28 | 西南石油大学 | 一种剪切缝页岩岩心脉冲衰减渗透率测试方法 |
US11782182B2 (en) | 2019-08-10 | 2023-10-10 | Research Institute Of Petroleum Exploration And Development | Systems and methods for combined physical and numerical simulation of subterranean characteristics |
CN113643603B (zh) * | 2021-08-05 | 2022-11-29 | 山东科技大学 | 分离式巷道物理实验模型结构、制作辅具以及制作方法 |
CN113866066B (zh) * | 2021-10-26 | 2023-11-03 | 贵州乌江水电开发有限责任公司东风发电厂 | 一种三维可视化水电站地下厂房洞室渗水监控方法及系统 |
CN117386349B (zh) * | 2023-10-31 | 2024-04-30 | 西南石油大学 | 基于产液剖面的致密油压裂水平井人工裂缝参数反演方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
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 |
FR2733073A1 (fr) * | 1995-04-12 | 1996-10-18 | Inst Francais Du Petrole | Methode pour modeliser un milieu geologique stratifie et fracture |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4442710A (en) * | 1982-03-05 | 1984-04-17 | Schlumberger Technology Corporation | Method of determining optimum cost-effective free flowing or gas lift well production |
US4797821A (en) * | 1987-04-02 | 1989-01-10 | Halliburton Company | Method of analyzing naturally fractured reservoirs |
US5070457A (en) * | 1990-06-08 | 1991-12-03 | Halliburton Company | Methods for design and analysis of subterranean fractures using net pressures |
-
1996
- 1996-12-30 FR FR9616330A patent/FR2757947B1/fr not_active Expired - Fee Related
-
1997
- 1997-12-23 GB GB9727227A patent/GB2322948B/en not_active Expired - Fee Related
- 1997-12-29 NL NL1007917A patent/NL1007917C2/nl not_active IP Right Cessation
- 1997-12-29 NO NO19976127A patent/NO326401B1/no not_active IP Right Cessation
- 1997-12-30 US US09/000,972 patent/US6023656A/en not_active Expired - Lifetime
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
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 |
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 (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
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 |
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 |
US6922662B2 (en) | 2000-05-26 | 2005-07-26 | Institut Francais Du Petrole | Method for modelling flows in a fractured medium crossed by large 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 |
FR2918179A1 (fr) * | 2007-06-29 | 2009-01-02 | Inst Francais Du Petrole | Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite |
EP2037080A1 (fr) | 2007-06-29 | 2009-03-18 | Ifp | Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité |
US8078405B2 (en) | 2007-06-29 | 2011-12-13 | Ifp | Method of estimating the permeability of a fracture network from a connectivity analysis |
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 |
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 |
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 |
CN109932129A (zh) * | 2017-12-15 | 2019-06-25 | 核工业北京地质研究院 | 一种表面接触式水压力传感器的快速标定方法 |
CN108756859A (zh) * | 2018-05-03 | 2018-11-06 | 中国石油天然气股份有限公司 | 气藏多层合采边水水侵机理的模拟实验装置及方法 |
CN113533157A (zh) * | 2021-07-02 | 2021-10-22 | 浙江大学 | 一种用于可视化实验的变开度可拆卸裂隙装置 |
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 |
---|---|
FR2757947B1 (fr) | 1999-01-29 |
NL1007917C2 (nl) | 1998-07-01 |
NO976127D0 (no) | 1997-12-29 |
GB2322948B (en) | 2002-01-09 |
GB9727227D0 (en) | 1998-02-25 |
US6023656A (en) | 2000-02-08 |
NO326401B1 (no) | 2008-11-24 |
GB2322948A (en) | 1998-09-09 |
NO976127L (no) | 1998-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2757947A1 (fr) | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches | |
Manzocchi et al. | Fault transmissibility multipliers for flow simulation models | |
Kuchuk et al. | Fractured-reservoir modeling and interpretation | |
Seifert et al. | Using sequential indicator simulation as a tool in reservoir description: Issues and uncertainties | |
Bourbiaux et al. | A rapid and efficient methodology to convert fractured reservoir images into a dual-porosity model | |
EP1158312A1 (fr) | Méthode pour modéliser des écoulements dans un milieu fracturé traversé par de grandes fractures | |
WO2009098366A2 (fr) | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures | |
Al-Shalabi et al. | Effect of Pore‐Scale Heterogeneity and Capillary‐Viscous Fingering on Commingled Waterflood Oil Recovery in Stratified Porous Media | |
Mullins et al. | The dynamics of reservoir fluids and their substantial systematic variations | |
Kløv et al. | Pore-to-field scale modeling of WAG | |
Milliotte et al. | Well-data-based discrete fracture and matrix modelling and flow-based upscaling of multilayer carbonate reservoir horizons | |
Fisher et al. | Quantifying Reservoir Heterogeneity Through Outcrop Characterization: 1. Architecture, Lithology, and Permeability Distribution of a Seaward-Stepping Fluvial-Deltaic Sequence, Ferron Sandstone (Cretaceous), Central Utah | |
Slatt | Geologic controls on reservoir quality | |
Boya Ferrero et al. | Predicting water in the crest of a giant gas field: Ormen Lange Hydrodynamic Aquifer Model | |
Barros-Galvis et al. | Analytical modeling and contradictions in limestone reservoirs: Breccias, vugs, and fractures | |
Jin et al. | Capillary pressure prediction from rock models reconstructed using well log data | |
Bihani et al. | On the role of spatially correlated heterogeneity in determining mudrock sealing capacity for CO2 sequestration | |
EP2912582B1 (fr) | Système, procédé et produit de programme informatique permettant d'évaluer et de classer des corps géologiques à l'aide d'une caractéristique d'euler | |
Rodriguez et al. | Recognition of Oil-Water Contact Architecture in Water Drive Reservoirs by Rock Typing Implementation, Sacha Field of Ecuador | |
Costa Gomes et al. | An open access carbonate reservoir benchmarking study for reservoir characterisation, uncertainty quantification & history matching | |
Thai et al. | Applying the evaluation results of porosity-permeability distribution characteristics based on hydraulic flow units (HFU) to improve the reliability in building a 3D geological model, GD field, Cuu Long Basin | |
Galvis | Geomechanics, fluid dynamics and well testing, applied to naturally fractured carbonate reservoirs: Extreme naturally fractured reservoirs | |
de Joussineau et al. | Organization, flow impact and modeling of natural fracture networks in a karstified carbonate bitumen reservoir: An example in the Grosmont Formation of the Athabasca Saleski leases, Alberta, Canada | |
Barros Galvis et al. | Phenomenology and contradictions in carbonate reservoirs | |
Khashman et al. | 3D Geological Modeling with Petrel Software: Estimating Original Oil in Place of Lower Cretaceous Y-Formation: A Case Study in a Selected Oilfield |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CD | Change of name or company name | ||
ST | Notification of lapse |
Effective date: 20150831 |