FR2863052A1 - Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse - Google Patents
Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse Download PDFInfo
- Publication number
- FR2863052A1 FR2863052A1 FR0314211A FR0314211A FR2863052A1 FR 2863052 A1 FR2863052 A1 FR 2863052A1 FR 0314211 A FR0314211 A FR 0314211A FR 0314211 A FR0314211 A FR 0314211A FR 2863052 A1 FR2863052 A1 FR 2863052A1
- Authority
- FR
- France
- Prior art keywords
- sample
- liquid
- quantities
- tensor
- flows
- 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
- 239000007788 liquid Substances 0.000 title claims abstract description 16
- 230000035699 permeability Effects 0.000 title claims description 31
- 238000012360 testing method Methods 0.000 title description 2
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000005259 measurement Methods 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 4
- 239000011435 rock Substances 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000011144 upstream manufacturing Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- 240000002234 Allium sativum Species 0.000 description 1
- 244000000626 Daucus carota Species 0.000 description 1
- 235000002767 Daucus carota Nutrition 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 235000004611 garlic Nutrition 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000005303 weighing Methods 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
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Complex Calculations (AREA)
Abstract
- Méthode pour déterminer l'ensemble des composantes d'un tenseur de perméabilité absolue d'un échantillon de milieu poreux à partir des mesures de perméabilité, obtenues par exemple en plaçant l'échantillon dans un perméamètre.- On impose une différence de pression ΔP entre les faces d'entrée et de sortie d'une "carotte" de laboratoire, avec des conditions de flux nul sur les bords parallèles à l'écoulement moyen obtenu en la confinant dans une gaine sous pression. En partant de mesures de perméabilité conventionnelles dans trois directions, et des mesures des deux composantes des forces visqueuses transverses à l'échantillon (soit en tout 9 mesures à 3D), on peut "inverser" le tenseur de perméabilité k par la résolution numérique du problème aux limites correspondant. Comme ces quantités peuvent alternativement être obtenues à partir de résolutions numériques du même écoulement en milieu hétérogène, la méthode peut être utilisée aussi en outil de changement d'échelle (upscaling) dans un simulateur de gisement.- Applications à la détermination rapide des anisotropies de perméabilité des roches, à la détection d'hétérogénéité internes.
Description
La présente invention concerne une méthode pour déterminer l'ensemble des
composantes 2D ou 3D d'un tenseur de perméabilité absolue d'un échantillon de milieu poreux etlou permettre une mise à l'échelle de champs de perméabilité hétérogènes, à partir de données mesurées en laboratoire ou issues de cartes de perméabilité d'une zone
souterraine fournies par les géologues.
ETAT DE LA TECHNIQUE
La mesure en laboratoire de la perméabilité des roches est un étape incontournable des études de réservoirs pétroliers, d'aquifères, en génie civil, voire de catalyseurs employés par l'industrie chimique. Si le milieu présente des directions d'anisotropie évidentes, l'expérimentateur tâchera de les respecter lors de sa campagne de mesures. En d'autre termes, il se placera d'emblée suivant les axes propres du tenseur de perméabilité. Il y a bien des cas où cette direction n'est pas donnée à l'avance, et il faut donc a priori déterminer le tenseur sans idée préconçue. Il peut aussi s'avérer que les échantillons de roche n'ont pas été carottés dans les bonnes directions; de plus les axes propres du tenseur ne sont pas forcement alignés avec le litage que peuvent présenter certaines carottes.
Il est donc particulièrement utilé de disposer de mesures directes de l'ensemble des composantes du tenseur, ou d'au moins en estimer l'éventuelle anisotropie.
Une méthode de détermination d'un tenseur de perméabilité d'un échantillon de roche à partir de mesures permettant d'estimer avec suffisamment de précision les pressions moyennes autour d'un échantillon placé dans un perméamètre adapté à ces mesures, est décrit par exemple dans la publication suivante: - Renard P. "Laboratory détermination of the full permeability tensor " JGR 106, B11 2001 pp26 443-26 451.
Une autre méthode connue de détermination d'un tenseur de perméabilité dans laquelle on crée dans l'échantillon des écoulements rendus volontairement tortueux par des modifications des conditions aux limites des faces d'injection, et on résout un problème inverse, est décrit par'exemple dans la publication suivante: - Bernabé Y.: "On the measurement of perméability in anisotropic rocks " in "Fault mechanics and transport properties of rocks " edited by B.Evans and TF Wong pp 147-167, Academic San Diego 1992 Pour être mises en oeuvre, les méthodes connues exigent une modification substantielle des conditions d'écoulement qui les rendent onéreuses et difficilement 10 applicables en pratique.
Dans le contexte d'un changement d'échelle (upscaling) bien connu des ingénieurs de réservoir, on est amené à remplacer l'expérimentation de laboratoire par les résultats d'une simulation numérique "fine" sur le milieu hétérogène, afin de le remplacer par un milieu homogène équivalent dont le tenseur de perméabilité que l'on note par exemple Keif. Dans le cas d'un milieu anisotrope, la simulation fine de référence se fait en général en considérant des conditions aux limites périodiques. La plupart des auteurs considèrent à tort que ce sont les seules conditions permettant d'obtenir l'ensemble des éléments du tenseur Keff.
Comme va le montrer la description qui va suivre, une interprétation adaptée d'une simulation numérique du même perméamètre est à même de fournir cette information. Les conditions aux limites imposées dans un perméamètre, peuvent être plus réalistes dans le cas où l'écoulement à grande échelle est contraint par des barrières argileuses. De plus, la comparaison entre les deux tenseurs résultat peuvent fournir des indications sur l'échelle du Volume Elémentaire Représentatif (VER).
LA METHODE SELON L'INVENTION La méthode selon l'invention permet de déterminer, à partir des mesures connues de peiméabilité, toutes les composantes d'un tenseur de perméabilité absolue (Keq,) d'un échantillon de milieu poreux anisotrope. Elle comporte essentiellement les étapes suivantes: - on mesure ou on calcule en imposant deux écoulements de liquide 'successivement suivant deux (ou trois) directions orthogonales (Ox, Oy) (et Oz) au travers de l'échantillon et sous gradient de pression imposé, les débits de liquide et les forces transversales générées dans l'échantillon du fait du passage du liquide, ces quantités étant ramenées à des gradients de pression unitaires; et - on utilise un simulateur pour déterminer de proche en proche, à partir d'une solution numérique a priori d'un écoulement confiné dans un milieu homogène, les composantes du tenseur équivalent (k) par une méthode inverse quelconque tel que les dites quantités correspondant s'ajustent sensiblement aux quantités mesurées ou calculées.
Suivant un premier mode de mise en oeuvre, on mesure les quantités (débits et forces transversales) au moyen d'un perméamètre et de moyens de mesure de forces transversales appliquées à l'échantillon en réponse à des écoulements imposés au travers de l'échantillon.
Ces forces transversales appliquées à l'échantillon, en réponse à des écoulements imposés à travers lui, peuvent être mesurées par exemple en mesurant les variations de poids de l'échantillon dans le perméamètre quand on alterne le sens de passage des écoulements de liquide imposés au travers de l'échantillon. En répétant cette procédure après avoir effectué une rotation de l'échantillon de 90 autour d'un axe parallèle à l'écoulement, on peut mesurer l'autre composante.
Suivant un deuxième mode de mise en oeuvre, on détermine les mêmes quantités par simulation d'écoulements dans un milieu, à partir de données de perméabilité connues ou estimées telles que des cartes de perméabilité fournies par exemple par des géologues dans les applications en géophysique.
La mise en oeuvre de la méthode selon l'invention est avantageuse en ce qu'elle permet de simplifier et d'accélérer la détermination du tenseur de perméabilité d'échantillons.
4 PRESENTATION SUCCINCTE DES FIGURES Les caractéristiques et avantages de la méthode selon l'invention, apparaîtront plus clairement à la lecture de la description ci-après d'un exemple non limitatif de mise en oeuvre, en se référant aux dessins annexés où : - la figure 1 montre, un organigramme, permettant une détermination itérative du tenseur équivalent Keq(n) d'un échantillon de roche, à partir de différentes mesures de débit et de forces effectuées en le confinant dans un perméamètre; - les figures 2 à 4 illustrent différentes exemples de simulation montrant la pertinence de la méthode, qui seront commentées dans le cours de la description; et -la figure 5 montre un exemple schématique de perméamètre convenant pour effectuer des mesures sur un échantillon de roche.
DESCRIPTION DETAILLEE DE LA METHODE
Les mesures de perméabilité utiles à la mise en oeuvre de la méthode qui va être décrite peuvent être obtenues au moyen d'un perméamètre de type connu. Un tel perméamètre comporte par exemple (fig.5), une cellule de confinement 1 pour l'échantillon. La paroi latérale de la cellule est une gaine déformable. A une première extrémité, la cellule 1 est connectée à un circuit 2 d'injection de fluide à une pression d'injection P1. A son extrémité opposée, la cellule est connectée à un circuit 3 d'évacuation du fluide issu de l'échantillon maintenu à une pression P2. La cellule 1 est placée dans une enceinte extérieure 4 connectée à un dispositif 5 pour appliquer une pression de confinement Pc à l'échantillon, de façon que le flux de liquide sur le pourtour latéral de l'échantillon soit nul (bords imperméables). Les embouts opposés de la cellule 1 sont tels que les pressions P1 et P2 respectivement exercées sur les faces amont et aval de l'échantillon sont uniformes. On note dPx,x = P2 Pi De façon à simplifier la présentation, on considère tout d'abord un milieu poreux 2D, carré de côté L. La méthode, transposée aux cas 3 D avec des proportions quelconques sera exposée plus loin. Dans toute cette partie, on considère que l'échantillon de roche poreuse est homogène mais avec un tenseur de perméabilité k dont les axes propres ne coïncident pas nécessairement avec les axes du perméamètre où il est placé (fig. 1).
Dans un premier temps, on applique le flux de liquide à l'échantillon orienté suivant une première direction A par rapport à l'axe de la cellule et on fait un ensemble de mesures que l'on va préciser plus loin. Dans un deuxième temps, on recommence l'expérience après avoir orienté l'échantillon selon une direction B orthogonale à la première et l'on fait les mêmes mesures selon le même protocole.
Sur les faces de l'échantillon parallèles à l'écoulement imposé, on a les conditions standard de flux nul de Neumann, soit n.k.Vp(r) = 0 (1) où r désigne le vecteur position dans l'espace, de composantes x, y (ou xi, x2) n est le 10 vecteur unitaire orthogonal à la frontière au point considéré, Vp(r) est le gradient de pression. Composantes par composantes, cette formule s'écrit: E n..k... ap(r) = 0, ou bien de façon plus concise encore en adoptant à partir de t,J=1,D axf maintenant la convention d'Einstein de sommation sur les indices répétés: k ap(r) n =0 axi - Equation en pression: L'équation en pression est donnée classiquement, où le point (.) correspond à des contractions de tenseurs: k 0. .Op(r) =0 ^ i soit, in extenso, toujours avec la convention d'Einstein: (2) (3) k. ap(r) axi a ax; L'équation en pression et les conditions aux limites définissent un problème noté problème (P) bien posé possédant une solution unique. Dans le cas particulier où la direction x coïncide avec un axe propre de k, la solution linéaire p(x, y, z) = P + APX,X É L est l'unique solution du problème. On obtient alors dans ce cas la relation classique Qx,x = kx,x x'x dont on tire kx x si Qx x et OPx x ont été simultanément mesurés.
Dans le cas anisotrope général, la solution p(x, y, z) ne possède pas de structure simple et il n'y a pas de solution analytique à notre connaissance. On se contente donc de supposer que l'on est dans la situation précédemment présentée et l'on fournit donc en pratique des estimations des valeurs propres du tenseur k.
Nous allons montrer maintenant que à condition de mesurer quelques quantités supplémentaires, il est possible de déterminer tous les éléments k du tenseur k. Pour cela, on remarque que l'équation (4) peut se réécrire sous la forme: fera axi ap(r) axi f (r) = O Vf (r) (5) Avec r désignant le vecteur position à D composantes, et dDr l'élément d'intégration. En choisissant alors la fonction f(r) de façon convenable (f(r) =x, y puis z) à 3D, en intégrant par parties l'égalité précédente on obtient le jeu d'égalités suivantes à 2D: xx = kxx Lz f L dy(p(x = L, y) p(x = 0, y)) + kxy Lz f L dx(p(x, y = L) - p(x, y = 0)) soit, avec des notations plus condensées n -k api+k ap r1 xz xx ax x xY ay (7 (6) ap 8x = qxx kYx + ax x (7) où qxx correspond au débit normalisé par le volume de l'échantillon q,x,x = Q,,x /V. Les {... }X désignent une moyenne spatiale et l'indice x indique ici la direction du gradient 20 de pression imposé. Les quantités Byx et Sxy sont définies par: f L y kxxj ôp (L, y) - ap (0, y) ]dy - fo y kxs[ âx (L, y) --- (0, y) ]dy Syx Vgxx fo kxx ax (0, y) dy Cette quantité correspond à la différence des coordonnées y des barycentres des débits en x locaux sur les faces aval et amont de l'échantillon lors d'un écoulement imposé selon x. De même on a: I x k , [ ap (x, L) ô y (x,0) l dx Jo x kYY [ â (x, L) ap (x,0) 1dx &y = J L o kYY â (x,0) Vg3Y y Des formules analogues peuvent être démontrées par la même technique à 3 D, introduisant les trois débits qxx, qyy, qzz, les six quantités (8 x Szx Sxy Sz,, Sxz et (5yz), et les six autres grandeurs additionnelles ( ap ap) ay. ' {al)} azx' axy' {al)} azy' ax lapl ay Z Les conditions de Dirichlet sur les faces amont et aval fournissent bien sûr la relation {a j '. . En considérant une expérience effectuée dans une autre direction ax x L orthogonale, y on obtient: ap 1 aP 4u q = kYx { ax} + k yy a y, y, (8) Y /iSxy qyy =kxx1ap1
Y (9)
De façon à alléger les notations, nous convenons de travailler avec des gradients de pression imposés unitaires ap = PxxouYY -1. En inversant ces égalités, on obtient dans axJ 0 L ce cas les relations suivantes: Ç a gxx gYY SxY la y x k,,_ 1 1ax1 Ylayfx agYY rg x( y ( 11) 1 lay xapl ly (10) kYY a kxy = {--P ax kxx + 8xyq
Y
ap l kyx = 1aY x }k +8yxgxx On peut rassembler ces relations sous une forme symbolique: (12) (14) k = F(qxx, qyy> ap Sap l ayJx'laxJy' x, axy) ' (15) Pour résumer, dans le cas homogène 2D, en mesurant les quantités qxx, gyy,lapl, ax y f ap 1, SyX, Sxy, on peut retrouver directement le tenseur k. Comme le tenseur est défini layJx par 3 nombres à 2D, et 6 à 3D (mais dans ce cas on détermine 15 quantités additionnelles), ceci signifie que des relations implicites existent entre les grandeurs pl ap x 8x. En théorie, si on retient au moins trois valeurs, on doit donc qxx, qyy, x {ax}y' y pouvoir encore retrouver k.
Nous allons retenir le sous ensemble qxx, qyy, {ap}x, {â J y (16) Le principal intérêt de ces données est de posséder un sens physique immédiat et partant d'être directement mesurables: que, qyy correspondent aux débits usuels (ramenés 15 au volume de l'expérience). Les quantités apx' azysont mesurables car: (ap = 1 fLdx(p(x, y = L) p(x, y = 0)) n'est autre qu'une résultante de la force jl ay x LZ Jo transverse par unité de volume imposée à l'échantillon. de plus, en utilisant la relation générale: f d rVp(1-) V fad -'r p(r)n, on voit que cette interprétation du gradient de
V
pression moyen en terme de force surfacique exercée sur l'échantillon subsiste.
Toutes ces quantités ou forces peuvent être mesurées directement en laboratoire. Pour mesurer les forces transversales, on peut utiliser par exemple des dynamomètres ou des jauges de contraintes appliquées localement contre les parois latérales de" l'échantillon dans sa gaine (cf. fig.5).
Pour mesurer plus commodément la résultante des forces transversales, il est possible également de mesurer les variations de poids de l'échantillon selon que les forces transversales se rajoutent ou se soustraient à son poids statique. Toutes choses égales par ailleurs, on provoque un écoulement au travers de l'échantillon successivement suivant un premier sens et suivant le sens opposé. Dans un premier cas, la force transversale s'ajoute par exemple au poids de la cellule et dans l'autre cas elle se retranche. En soustrayant les mesures, on peut obtenir directement la mesure de la force transversale. Pour la mise en oeuvre, on pose l'enceinte 4 (cf. fig.5) à l'horizontale sur un dispositif de pesée (non représenté). On mesure la variation de poids pour deux positions à 90 l'une de l'autre correspondant aux directions A et B fléchées sur la figure 1.
Exemple 2D
On suppose que l'on a obtenu comme indiqué le quadruplet que, qyy, {aF' 1ax yde mesures de débits et de forces transversales appliquées à l'échantillon dans les positions successives qu'on lui fait prendre dans le perméamètre (4 mesures globales de forces et de débits à 2D, 9 à 3D).
A partir de maintenant, les milieux poreux considérés peuvent être hétérogènes, c'est à dire caractérisés par une carte fine de tenseurs de perméabilité dépendant de la position, k(r) qui peut être fournie par le géologue ou tout autre technique.
De façon à remonter au tenseur k que l'on notera Keq à partir de maintenant pour rappeler qu'il s'agit d'un tenseur effectif à partir de la connaissance de 25}ap}X' {aply, il est nécessaire de disposer d'une solution numérique de l'équation de Laplace anisotrope. Toute méthode connue d'inversion peut être utilisée. L'utilisateur peut choisir par exemple une méthode d' "optimisation" (pour minimiser l'écart entre les valeurs de ap mesurées et les prédictions du modèle numérique), une méthode par q,g" 'laP ax}Y essai et erreur, par dichotomie, etc. Un exemple est donné ci-après.
De façon à trouver le tenseur Keq, associé au milieu homogène équivalent, on va utiliser un simulateur numérique d'un type connu disponible sur le marché, opérant par éléments ou volumes finis, pour résoudre le problème aux limites 2D (ou 3D). La procédure consiste par exemple à construire une suite de tenseurs Keq(1),..., Keq(n), etc. définie par son premier terme représentant une première estimation de Keq obtenue en utilisant la relation (15) : Keq(1) F(q.,'YY,ap x,apY,O,o), Comme les moments SyXet 8xy sont inconnus, car non mesurés, ceci ne peut être qu'une approximation puisque comme on l'a vu, il existe une relation implicite entre les moments et q q ap. Cette première estimation ne la prend pas en compte et YY Li ay x {ail aY doit donc être corrigée. On modifie l'estimation de Keq(n) par exemple par la relation de récurrence suivante: Keq(n+1)=F(q q) , }x, 1'1y Sx;, ,x).
Ici, les quantités Sxn, fjy fl sont calculées à l'aide du simulateur numérique: on y, résout le problème aux limites (P) en prenant comme tenseur de perméabilité l'estimation précédente Keq(n). On arrête les itérations une fois que les quantités f p a Y, } x a calculées par le simulateur avec l'estimation correspondante Keq(n) sont à e de la valeur 20 désirée, étant un critère d'arrêt fixé par l'utilisateur par rapport à une norme choisie à l'avance. L'organigramme est indiqué à la figure 1.
Le quadruplet de mesure utilisé pour déterminer cômplètement le tenseur k peut également avoir été calculé à l'aide d'une "simulation fine" de l'expérience sur un milieu hétérogène dont on souhaite définir le tenseur de perméabilité équivalent par cette méthode. La perméabilité de ce milieu est connue en tout ou en partie, en se basant par exemple sur une carte de perméabilité fournie par le géologue. Durant cette simulation, on applique les lois de Darcy aux écoulements dans le milieu et on en déduit le quadruplet de valeurs applicable.
Validation Par simulation, nous avons vérifié l'existence de forces transversales, mesurables et non négligeables dans certains cas, apparaissant lorsqu'un écoulement traverse un milieu confiné, due à la présence d'hétérogénéités internes ou à une anisotropie du milieu, et validé la possibilité d'identifier toutes les composantes du tenseur de perméabilité équivalent grâce à des mesures simples de forces sur échantillons obtenues au moyen d'un perméamètre. En plus, cette étude nous a permis d'établir l'erreur commise sur le comportement d'un milieu équivalent (défini par une méthode ou une autre) par rapport au comportement du milieu réel à homogénéiser. C'est sur cette base d'un meilleur ajustement entre les comportements respectifs que la méthode selon l'invention a été développée.
Sur la fig.2B, on voit l'allure du champ de pression pour un milieu homogène anisotrope avec un rapport d'anisotropie de 10 et dont les axes principaux du tenseur diagonal (fig.2A) sont tournés de 45 par rapport à l'horizontale (0 = 45 ). Sur la fig.2C, on montre l'évolution du rapport de la différence cumulée de pression sur les bords imperméables (bords perpendiculaire aux champs de pression imposé) sur la différence de pression appliquée sur l'entrée-sortie du milieu (bords parallèles aux pressions imposées) pour 0 (0 <_0_<180 ).
Sur la fig.3A, on montre une carte de perméabilité log-normale (moyenne géométrique Kg=121, a2=2) générée avec le code de génération stochastique bien connu dit FFTMA. Localement, la perméabilité est isotrope (scalaire), en revanche sa fonction de corrélation est spatialement anisotrope, orientée à 45 . Sur la fig.3B,on a représenté l'évolution de la pression aux bords haut et bas (bords imperméables pour un écoulement confiné horizontal) pour le milieu réel et son équivalent homogène anisotrope déterminé par la méthode du perméamètre. Sur les fig.3C, 3D respectivement, on montre les cartographies du champ de pression pour le milieu réel (à gauche) et le milieu équivalent (à droite). Après inversion, on estime les composantes du tenseur de perméabilité du milieu équivalent valent: KXa=122.05; kyy=119.86; kxy=kyx=58.08. Les valeurs propres de ce tenseur valent: K',,,=179.05; K'yy=62.86. Le système d'axes principaux est bien tourné de 44.46 par rapport à la direction de l'écoulement imposé. Ceci correspond bien à la direction de l'anisotropie de départ.
Sur la fig.4A, on montre le réseau de fractures dont la fraction volumique des fractures vaut 10%, la perméabilité de la matrice est de l'ordre de 1 Darcy alors que celle des fractures est de l'ordre de 100 Darcy's (à gauche). La fig.4B montre l'évolution de la pression aux bords haut et bas (bords imperméables pour un écoulement confiné horizontal) pour le milieu réel et son équivalent homogène anisotrope déterminé par la méthode du perméamètre. Sur les fig.4C, 4D respectivement, on montre les cartographies du champ de pression pour le milieu réel et son équivalent. Les composantes du tenseur de perméabilité du milieu équivalent valent: KXX=5.79 D; kyy=4.37 D; kXy=ky,=0.26 D. Les valeurs propres de ce tenseur valent: K'X,F5.83 D; K'yy=4.32 D et dont le système d'axes principaux sont tournés à 9.83 par rapport au système d'axe formé par la direction horizontale et celle verticale.
Claims (4)
1) Méthode pour déterminer, à partir des mesures connues de perméabilité, toutes les composantes d'un tenseur de perméabilité absolue (Keq,) d'un échantillon de milieu 5 poreux anisotrope, dans laquelle: on mesure ou bien on calcule à partir de deux (ou trois) écoulements selon Ox, Oy (et Oz) sous gradient de pression imposé les quantités les débits de liquide et les composantes des forces transversales générées dans l'échantillon du fait du passage du liquide, ces quantités étant ramenées à des gradients de pression unitaires; et - - on utilise un simulateur pour déterminer de proche en proche, à partir d'une solution numérique a priori d'un écoulement confiné dans un milieu homogène, les composantes du tenseur équivalent (k) par une méthode inverse quelconque tel que les dites quantités correspondant s'ajustent sensiblement aux quantités mesurées ou calculées.
2) Méthode selon la revendication 1, dans laquelle on mesure les dites quantités au moyen d'un perméamètre et de moyens de mesure de forces transversales appliquées à l'échantillon en réponse à des écoulements imposés au travers de l'échantillon.
3) Méthode selon la revendication 1., dans laquelle on détermine les dites quantités par simulation d'écoulements dans un milieu, à partir de données de perméabilité connues 20 ou estimées.
4) Méthode selon la revendication 2, dans laquelle on détermine les forces transversales appliquées à l'échantillon en réponse à des écoulements imposés au travers de l'échantillon, eh mesurant les variations de poids apparent de l'échantillon dans le perméamètre quand on alterne le sens de passage des écoulements de liquide au travers de l'échantillon.
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0314211A FR2863052B1 (fr) | 2003-12-02 | 2003-12-02 | Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse |
NO20044943A NO20044943L (no) | 2003-12-02 | 2004-11-12 | Fremgangsmate for bestemmelse av komponentene i en absolutt permeabilitets-tensor for en poros bergart |
GB0426320A GB2409730B (en) | 2003-12-02 | 2004-12-01 | Method of determining the components of an effective permeability tensor of a porous rock |
US11/001,035 US7191071B2 (en) | 2003-12-02 | 2004-12-02 | Method of determining the components of an effective permeability tensor of a porous rock |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0314211A FR2863052B1 (fr) | 2003-12-02 | 2003-12-02 | Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2863052A1 true FR2863052A1 (fr) | 2005-06-03 |
FR2863052B1 FR2863052B1 (fr) | 2006-02-24 |
Family
ID=34043798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0314211A Expired - Fee Related FR2863052B1 (fr) | 2003-12-02 | 2003-12-02 | Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse |
Country Status (4)
Country | Link |
---|---|
US (1) | US7191071B2 (fr) |
FR (1) | FR2863052B1 (fr) |
GB (1) | GB2409730B (fr) |
NO (1) | NO20044943L (fr) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1927846A1 (fr) * | 2006-12-01 | 2008-06-04 | Ifp | Méthode pour caractériser la distribution de la perméabilité absolue d'un échantillon hétérogène |
RU2549216C2 (ru) * | 2010-01-22 | 2015-04-20 | Тоталь Са | Измерение параметров, связанных с прохождением текучих сред в пористом материале |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7363161B2 (en) * | 2005-06-03 | 2008-04-22 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
US8140309B2 (en) * | 2007-07-18 | 2012-03-20 | Council Of Scientific & Industrial Research | Method of predicting the dynamic behavior of water table in an anisotropic unconfined aquifer having a general time-varying recharge rate from multiple rectangular recharge basins |
US8744817B2 (en) * | 2007-12-21 | 2014-06-03 | Schlumberger Technology Corporation | Method for upscaling a reservoir model using deep reading measurements |
US8738341B2 (en) * | 2007-12-21 | 2014-05-27 | Schlumberger Technology Corporation | Method for reservoir characterization and monitoring including deep reading quad combo measurements |
EP2564309A4 (fr) | 2010-04-30 | 2017-12-20 | Exxonmobil Upstream Research Company | Procédé et système de simulation de flux par volumes finis |
EP2599032A4 (fr) | 2010-07-29 | 2018-01-17 | Exxonmobil Upstream Research Company | Procédé et système de modélisation d'un réservoir |
WO2012015515A1 (fr) | 2010-07-29 | 2012-02-02 | Exxonmobil Upstream Research Company | Procédés et systèmes pour une simulation de flux par apprentissage automatique |
CA2803315A1 (fr) | 2010-07-29 | 2012-02-02 | Exxonmobil Upstream Research Company | Procedes et systemes pour une simulation de flux par apprentissage automatique |
US9187984B2 (en) | 2010-07-29 | 2015-11-17 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
AU2011289788B2 (en) | 2010-08-09 | 2015-03-19 | Conocophillips Company | Reservoir upscaling method with preserved transmissibility |
CA2807300C (fr) | 2010-09-20 | 2017-01-03 | Exxonmobil Upstream Research Company | Formulations souples et adaptatives pour des simulations de gisements complexes |
AU2011320352B2 (en) | 2010-10-29 | 2015-04-23 | Schlumberger Technology B.V. | Model based inversion of seismic response for determining formation properties |
WO2013039606A1 (fr) | 2011-09-15 | 2013-03-21 | Exxonmobil Upstream Research Company | Opérations matricielles et vectorielles optimisées dans des algorithmes à instructions limitées qui effectuent des calculs eos |
EP2901363A4 (fr) | 2012-09-28 | 2016-06-01 | Exxonmobil Upstream Res Co | Suppression des failles dans des modèles géologiques |
CN102937561B (zh) * | 2012-10-17 | 2014-07-02 | 西北工业大学 | 一种正交非织造三维矩形织物渗透率的确定方法 |
US9146333B2 (en) * | 2012-10-23 | 2015-09-29 | Schlumberger Technology Corporation | Systems and methods for collecting measurements and/or samples from within a borehole formed in a subsurface reservoir using a wireless interface |
US9416642B2 (en) | 2013-02-01 | 2016-08-16 | Halliburton Energy Services, Inc. | Modeling subterranean rock blocks in an injection treatment simulation |
US9798042B2 (en) | 2013-02-01 | 2017-10-24 | Halliburton Energy Services, Inc. | Simulating an injection treatment of a subterranean zone |
US9297250B2 (en) | 2013-03-14 | 2016-03-29 | Halliburton Energy Services, Inc. | Controlling net treating pressure in a subterranean region |
US9217318B2 (en) | 2013-03-14 | 2015-12-22 | Halliburton Energy Services, Inc. | Determining a target net treating pressure for a subterranean region |
CN103235888B (zh) * | 2013-04-27 | 2016-04-06 | 杭州电子科技大学 | 一种精确计算双各向异性介质球电磁散射的方法 |
CN103335928B (zh) * | 2013-05-30 | 2015-03-04 | 中国石油天然气集团公司 | 一种测量孔隙岩石渗透率的方法和装置 |
US9239407B2 (en) | 2013-08-27 | 2016-01-19 | Halliburton Energy Services, Inc. | Injection treatment simulation using condensation |
EP3175265A1 (fr) | 2014-07-30 | 2017-06-07 | ExxonMobil Upstream Research Company | Procédé de génération de maillage volumétrique dans un domaine ayant des propriétés de matériau hétérogènes |
EP3213126A1 (fr) | 2014-10-31 | 2017-09-06 | Exxonmobil Upstream Research Company | Gestion de discontinuité de domaine dans un modèle de grille de sous-surface à l'aide de techniques d'optimisation de grille |
CA2963092C (fr) | 2014-10-31 | 2021-07-06 | Exxonmobil Upstream Research Company | Procedes de gestion de discontinuite dans la construction d'espace de conception de modele de subsurface faillee a l'aide de moindres carres mobiles |
CN104809275A (zh) * | 2015-04-13 | 2015-07-29 | 中国石油大学(北京) | 碳酸盐岩微观流动模拟方法及装置 |
CN104866706A (zh) * | 2015-04-13 | 2015-08-26 | 中国石油大学(北京) | 碳酸盐岩渗透率确定方法及装置 |
CN105572012B (zh) * | 2015-12-22 | 2018-03-27 | 山东大学 | 一种充填贯通裂隙岩石渗流性能检测方法 |
CN106053319A (zh) * | 2016-07-29 | 2016-10-26 | 中国电建集团华东勘测设计研究院有限公司 | 一种各向异性岩体渗透梯度测试装置及方法 |
US10839114B2 (en) | 2016-12-23 | 2020-11-17 | Exxonmobil Upstream Research Company | Method and system for stable and efficient reservoir simulation using stability proxies |
US10801943B2 (en) * | 2017-06-26 | 2020-10-13 | China University Of Petroleum-Beijing | Apparatus and method for measuring apparent permeability of tight rock core |
CN110671100B (zh) * | 2019-10-10 | 2022-08-30 | 东北石油大学 | 一种制造模拟岩石非均质性装置中棋盘状仿真体的方法 |
CN112485096B (zh) * | 2020-11-16 | 2022-01-07 | 中国石油大学(北京) | 用于地热试验的具有非均质岩样的模拟地层的制作方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2862765B1 (fr) * | 2003-11-20 | 2006-10-27 | Inst Francais Du Petrole | Methode pour former un modele stochastique optimal d'une zone souterraine heterogene, cale sur des donnees dynamiques par parametrisation de distributions continues |
-
2003
- 2003-12-02 FR FR0314211A patent/FR2863052B1/fr not_active Expired - Fee Related
-
2004
- 2004-11-12 NO NO20044943A patent/NO20044943L/no not_active Application Discontinuation
- 2004-12-01 GB GB0426320A patent/GB2409730B/en not_active Expired - Fee Related
- 2004-12-02 US US11/001,035 patent/US7191071B2/en not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
BERNABE Y: "ON THE MEASUREMENT OF PERMEABILITY IN ANISOTROPIC ROCKS", FAULT MECHANICS AND TRANSPORT PROPERTIES OF ROCKS, ISBN 0-12-243780-2, B. EDITED BY EVANS AND T.F. WONG ,ACADEMIC SAN DIEGO, XP008032906 * |
DURLOFSKY L J: "NUMERICAL CALCULATION OF EQUIVALENT GRID BLOCK PERMEABILITY TENSORS FOR HETEROGENEOUS POROUS MEDIA", WATER RESOURCES RESEARCH, vol. 27, no. 5, May 1991 (1991-05-01), pages 699 - 708, XP008033632 * |
RENARD P ET AL: "LABORATORY DETERMINATION OF THE FULL PERMEABILITY TENSOR", JOURNAL OF GEOGRAPHICAL RESEARCH, RICHMOND, VA, US, vol. 106, no. B11, 10 November 2001 (2001-11-10), pages 26,443 - 26,452, XP008032904, ISSN: 0148-0227 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1927846A1 (fr) * | 2006-12-01 | 2008-06-04 | Ifp | Méthode pour caractériser la distribution de la perméabilité absolue d'un échantillon hétérogène |
FR2909448A1 (fr) * | 2006-12-01 | 2008-06-06 | Inst Francais Du Petrole | Methode pour caracteriser la distribution de la permeabilite absolue d'un echantillon heterogene |
US7693677B2 (en) | 2006-12-01 | 2010-04-06 | Institut Francais Du Petrole | Method of characterizing the distribution of the absolute permeability of a heterogeneous sample |
RU2549216C2 (ru) * | 2010-01-22 | 2015-04-20 | Тоталь Са | Измерение параметров, связанных с прохождением текучих сред в пористом материале |
Also Published As
Publication number | Publication date |
---|---|
US7191071B2 (en) | 2007-03-13 |
NO20044943L (no) | 2005-06-03 |
GB2409730B (en) | 2007-01-31 |
GB2409730A (en) | 2005-07-06 |
FR2863052B1 (fr) | 2006-02-24 |
GB0426320D0 (en) | 2005-01-05 |
US20050229680A1 (en) | 2005-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2863052A1 (fr) | Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse | |
Saxena et al. | Estimating permeability from thin sections without reconstruction: Digital rock study of 3D properties from 2D images | |
Sochi | Flow of non‐Newtonian fluids in porous media | |
Zou et al. | Experimental and theoretical evidence for increased ganglion dynamics during fractional flow in mixed‐wet porous media | |
Joekar‐Niasar et al. | Trapping and hysteresis in two‐phase flow in porous media: A pore‐network study | |
Walderhaug et al. | Prediction of permeability in quartz-rich sandstones: examples from the Norwegian continental shelf and the Fontainebleau sandstone | |
Berg et al. | Non-uniqueness and uncertainty quantification of relative permeability measurements by inverse modelling | |
EP1927846A1 (fr) | Méthode pour caractériser la distribution de la perméabilité absolue d'un échantillon hétérogène | |
FR2823877A1 (fr) | Methode pour contraindre par des donnees dynamiques de production un modele fin representatif de la repartition dans le gisement d'une grandeur physique caracteristique de la structure du sous-sol | |
Behbahani et al. | Analysis of imbibition in mixed-wet rocks using pore-scale modeling | |
CA2688842A1 (fr) | Methode pour determiner l'evolution de proprietes petrophysiques d'une roche au cours de la diagenese | |
WO2002006794A1 (fr) | Methode pour modeliser des deplacements de fluides dans un milieu poreux tenant compte d'effets d'hysteresis | |
EP2037080A1 (fr) | Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité | |
EP2400320B1 (fr) | Procédé d'exploitation d'un gisement pétrolier à partir d'un modèle de réservoir déformé graduellement au moyen de cosimulations | |
Durlofsky | Use of higher moments for the description of upscaled, process independent relative permeabilities | |
Al-Shalabi et al. | Single-well chemical-tracer modeling of low-salinity-water injection in carbonates | |
Herbas et al. | Comprehensive micromodel study to evaluate polymer EOR in unconsolidated sand reservoirs | |
Abadpour et al. | Integrated geo-modeling and ensemble history matching of complex fractured carbonate and deep offshore turbidite fields, generation of several geologically coherent solutions using ensemble methods | |
EP1889999B1 (fr) | Méthode pour optimiser la récupération assistée d'un fluide en place dans un milieu poreux par suivi de front | |
Leray et al. | Multi-objective assisted inversion of chemical EOR corefloods for improving the predictive capacity of numerical models | |
Di Stefano et al. | Estimating flow resistance in steep slope rills | |
EP2770162A1 (fr) | Procédé d'exploitation d'un réservoir géologique au moyen d'un modèle de réservoir calé et cohérent vis à vis des propriétés d'écoulement | |
Han et al. | Effect of the pore structure on resistivity index curves | |
Newgord et al. | Wettability quantification in mixed-wet rocks using a new NMR-based method: Experimental model verification | |
Goodfield et al. | New coreflood interpretation method for relative permeabilities based on direct processing of in-situ saturation data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CD | Change of name or company name | ||
ST | Notification of lapse |
Effective date: 20150831 |