Brevet FR2904982A1 — publié le 20080222

La présente invention concerne une méthode pour optimiser l'exploitation d'un milieu poreux hétérogène dans le cadre de récupération assistée d'un fluide en place, par détermination de la position du front séparant un fluide de balayage et le fluide en place. En particulier la méthode selon l'invention permet de déterminer la position d'un front séparant deux fluides immiscibles en mouvement, sans remise à jour systématique du champ 15 de pression, de façon à obtenir une simulation des écoulements polyphasiques dans le milieu poreux suffisamment rapide et précise pour obtenir des informations quantitatives permettant une exploitation optimale du milieu. La méthode trouve des applications notamment pour l'exploitation de gisements de pétrole ou de gaz, ou l'exploitation de stockages souterrains de gaz par exemple. Présentation de l'art antérieur Toutes les notations utilisées pour décrire l'art antérieur et l'invention sont définies en fin de description. Savoir décrire et simuler les écoulements polyphasiques dans les réservoirs souterrains 25 est à la base du métier des ingénieurs de gisement dans les compagnies pétrolières ou gazières (et de façon analogue dans des compagnies d'adduction d'eau). La présence des hétérogénéités du sous-sol, et la complexité croissante des systèmes de drainage rendent impossible toute solution analytique simple, ce qui impose le développement de solutions numériques utilisant un modèle maillé. La simulation des écoulements polyphasiques en milieu poreux hétérogène 20 2 2904982 peut alors requérir d'importantes ressources informatiques, en particulier lorsque le modèle numérique du milieu considéré est fortement détaillé. C'est le cas notamment en ingénierie de réservoirs, dans le domaine pétrolier. Ce coût provient principalement de la résolution des systèmes linéaires de grande taille issus de l'équation gouvernant la pression, qu'il convient de réactualiser en suivant le déplacement des fluides, si l'on désire disposer d'une solution possédant une bonne précision. Dans le cadre d'un écoulement diphasique dans un milieu poreux hétérogène, on considère le déplacement d'un fluide en place (par exemple l'huile) sous l'effet de l'injection d'un autre fluide (par exemple de l'eau). La loi de Darcy généralisée gouvernant le mouvement des fluides s'écrit, moyennant des hypothèses et notations standard, de la façon suivante : unw = ù/tn YPnw, 2nw = rnw TTrT~ /unw il/S uw -2vpw~ nw ù rw w @1w u=unw+uw où les indices nw et w note respectivement les fluides non mouillant et mouillant. La différence des pressions entre les deux fluides est notée pc (S) = pnw ù pw . On a la relation : Snw + Sw = 1 (4) De ce fait, on peut écrire toutes les fonctions dépendant de la saturation S en terme de la saturation Sw. En négligeant la pression capillaire, on a P. = Pw = P (5) La vitesse totale prend alors la forme u = ùa,(Sw)V/p (6) avec D.0 = 0 (7) Ce système d'équations (6) et (7) définit l'équation en pression. Dans les équations (1) et (2), les fonctions kr sont les perméabilités relatives. On a k,. = kr (Sw) . Les fonctions /Lnw, 2w sont les mobilités respectives des fluides, et 2(Sw) = 2nw (Sw) + 2 (Sw) est la mobilité totale, qui dépend donc explicitement de Sw. (1) (2) (3) 3 2904982 Les lois de conservation de la masse relatives à chacun des fluides s'écrivent aatw + Q • unw =: 0 a~+V uw=0 Ce système d'équations (8) et (9) définit l'équation en saturation. Ici est la porosité 5 (supposée uniforme dans le réservoir), t le temps, et u n, nw est la vitesse du fluide considéré. Une des difficultés pour la résolution de ces équations vient du couplage entre les équations (6) à (9) qui couplent la saturation au champ de pression. Ces effets sont bien connus et pilotent le développement d'éventuelles instabilités visqueuses. En particulier, dans le cas 1 D, en considérant le balayage d'un milieu initialement saturé en huile par de l'eau, ce 10 système d'équations peut être résolu par la méthode des caractéristiques, dont les solutions sont caractérisées par l'existence d'oscillations, correspondant à des saut de saturation se propageant dans le milieu. La discontinuité est caractérisée par une saturation en eau passant de Swi à Sf, la saturation au front dont la valeur est donnée par une condition de type Rankine-Hugoniot. 15 Dans le contexte pétrolier, la construction de la tangente de Welge permet de déterminer Sf graphiquement. En arrière de ce front, une "onde de raréfaction" assure la transition entre Sf et la saturation en eau à l'injection. On connaît plusieurs techniques pour déterminer la simulation des écoulements polyphasiques en milieu poreux hétérogène. Ces méthodes reposent toutes sur la résolution du 20 système d'équations suivant : Équation en pression : u = -2(Sw)Op V.u=0 Équation en saturation : 25 0 nw + 0 •un . = 0 a~+Q uw=0 Les techniques se différencient par la méthode utilisée pour résoudre ce système. Dans le cas général, en milieu hétérogène complexe, ce type de système d'équations doit être résolu 30 numériquement. (8) (9) (10) (11) (12) (13) 4 2904982 Les compagnies pétrolières utilisent le plus souvent les techniques de discrétisation classiques, de type volumes finis (Aziz, Kb., Settary, A.: Petroleum reservoir simulation. Applied science publishers, London, 1979). Concernant la discrétisation temporelle, les choix sont multiples : si on privilégie la robustesse numérique, on peut être totalement implicite en 5 pression comme en saturation. Si on préfère éviter des oscillations de la solution, et une meilleure précision, on choisira un schéma implicite en pression, explicite en saturation. On constate en général l'existence de fronts caractérisés par la même discontinuité de la saturation dont la valeur passe brutalement de la saturation irréductible en eau, &d, à la saturation du front, notée Sf. Ces discontinuités sont dues au caractère hyperbolique de l'équation de 10 transport en saturation, et sont bien décrites analytiquement dans le cas du déplacement 1D : c'est la théorie de Buckley-Leverett dont on trouvera une description détaillée dans l'ouvrage de Charles-Michel Marle. Les écoulements polyphasiques en milieu poreux. Seconde édition revue et augmentée 1972. Editions Technip, Paris. Quelque soit la méthode retenue, des remises à jour régulières de l'équation en pression 15 sont nécessaires, de façon à estimer les vitesses avec précision. Dans une autre classe de techniques, qui deviennent de plus en plus employées dans les compagnies opératrices, on trouve les techniques dites de lignes de courant, dont une référence est : R. P. Batycky, M.J. Blunt, and M.R. Thiele, A 3d field-scale streamline-based reservoir simulator, SPERE, 11, 246-254 (1997). 20 Dans ces méthodes les saturations sont mises à jour en suivant le déplacement des fluides dans leur mouvement le long des lignes de courant, courbes paramétrées définies par dx/dt = u(x,t). Un changement de variable permet de se ramener à la théorie de Buckley-Leverett précédemment évoquée. Toutefois, afin de bien calculer la vitesse u(x, t) la pression est résolue autant de fois que la précision l'exige, sur une grille cartésienne. Ces méthodes 25 sont plus rapides que les méthodes classiques. Elles n'en possèdent pas non plus la versatilité et la robustesse pour traiter des problèmes compositionnels complexes, où la méthode traditionnelle reste la seule issue possible. Les mises à jour onéreuses du champ de pression ne sont pas évitées. Enfin, il existe les techniques de suivi de front (R Juanes a,d KA Lie "A front tracking 30 method for efficient simulation of miscible gas injection." paper SPE 92298, Houston 31 jan. to 2 feb. 2005), dont l'idée est de suivre le front par une approche lagrangienne, qui nécessite les mêmes mises à jour de la pression. Dans ce type de méthode, la description des 5 2904982 écoulements polyphasiques dans les réservoirs souterrains consiste à déterminer la position du front (appelé aussi interface) séparant deux fluides immiscibles en mouvement : un fluide en place (de l'huile par exemple) et un fluide de balayage (de l'eau par exemple), appelé également fluide injecté 5 L'évolution du front dans le réservoir au cours de l'écoulement est considérablement influencée par le couplage (appelé couplage visqueux par les ingénieurs réservoirs) entre le champ de pression et le champ de saturation (Saffman P.G. and G. Taylor. The penetration of a fluid into a porous medium or hele-shaw cell containing a more viscous liquid. Proc. Royal Society of London, A245:312-329, 1958, Noetinger B., V. Artus, and L.Ricard. Dynamics of 10 the water-oil front for two-phase, immiscible flows in heterogeneous porous media. 2isotropie media. Transport in porous media, 56:305-328, 2004). En particulier, lorsque le fluide injecté est moins visqueux et par conséquent plus mobile au niveau du front que le fluide en place, les instabilités visqueuses favoriseront toujours l'écoulement des fluides dans les couches les plus perméables du milieu. Le temps de 15 percée au travers de celles-ci est beaucoup plus rapide que dans le reste du réservoir. Au contraire, si le fluide injecté est moins mobile, le couplage visqueux peut le ralentir dans les couches initialement plus rapides, compensant ainsi les différences de perméabilité dues à la stratification. Un front stationnaire apparaît alors (Artus V., B. Noetinger, and L. Ricard. Dynamics of the water-oil front for two-phase, immiscible flows in heterogeneous porous 20 media. 1-stratified media. Transport in Porous Media, 56:283-303, 2004). Il est connu que la stabilité du front suite à de petites perturbations est déterminée par le rapport des viscosités des fluides. Si le fluide le moins visqueux est déplacé par le fluide le plus visqueux, alors le front est stable, et vice versa. Si la situation est instable, alors des phénomènes, appelés doigts visqueux , se développent avec le temps. 25 Les effets de la perméabilité relative sont pris en compte en considérant les déplacements de type Bucley-Leverett dans un milieu homogène. De tels déplacements sont décrits par une équation en saturation hyperbolique. Une solution à cette équation par la méthode des caractéristiques fournis des solutions qui ont des oscillations caractérisées par une condition de Rankine Hugoniot. Il a été montré que la stabilité du front où la saturation a 30 un choc est déterminée par un rapport de mobilité frontal, noté Mf : si Mf 0). 11 2904982 On peut utiliser une discrétisation de type volumes finis cinq points standard pour résoudre cette équation à l'aide du simulateur : la pression est estimée au centre des mailles, les flux entre mailles faisant intervenir des transmissivités calculées comme moyenne harmonique des perméabilités des deux mailles. La loi de Darcy permet alors de déterminer 5 les flux et les vitesses. La carte de perméabilité K(x, y) peut être générée à l'aide d'un générateur de champs aléatoire, comme décrit dans le document suivant par exemple : Le Ravalec, M., Noetinger, B., and Hu, L.Y. [2000] The FFT moving average (FFTMA) generator: An efficient numerical method for generating and conditioning Gaussian 10 simulations. Mathematical Geology 32(6), 701-723. Définition d'une relation décrivant laposition du front en milieu hétérogène L'idée de base consiste donc à éviter de recalculer à chaque itération en temps l'écoulement dans tout le domaine de calcul, en se restreignant à la description du mouvement 15 du front, décrit par la fonction bh(y, t). Ceci est rendu possible par une approximation permettant de relier la modification du champ de pression bp(x, y, t), et donc des vitesses, à la distorsion du front bh(y, t). On peut ainsi éliminer la pression du problème et on est alors amenés à la résolution d'un système différentiel impliquant bh(y, t). On a vu que dans un milieu homogène, l'évolution de la transformée de Fourier des 20 fluctuations de la position du front, fluctuations notées bh(q, t), est donnée par : M -1 a 8h(q, t) = co I q I Mj .8h(q, t) (14) j +1 _ uû Jw(sj)ùfw(swr) ua f w Co \sj) (Y sf swr q5 Avec : 25 Mf : le rapport de mobilité frontal uo : la vitesse de filtration moyenne le long de la direction de l'axe x q5 : la porosité où 12 2904982 la fonction de Buckley-Leverett représentant le flux fractionnaire en eau. Saturation en eau au front Saturation maximale en eau vecteur d'onde dans l'espace de Fourier 5 En milieu hétérogène, le problème est plus complexe. Les hétérogénéités du champ de perméabilité initient des perturbations du front, c'est-à-dire que le front se déforme en fonction de la peiinéabilité des strates. Ainsi, ces perturbations peuvent grossir ou diminuer en fonction du rapport des viscosités. C'est le problème du couplage visqueux. Le problème 10 est non linéaire et une analyse mathématique s'avère très compliquée. Selon l'invention, on propose de tenir compte des hétérogénéités, d'une part en s'affranchissant du couplage visqueux à l'aide de la théorie des perturbations, et d'autre part en prenant en compte les fluctuations transversales de la vitesse, c'est-à-dire les fluctuations dans un axe Y perpendiculaire à l'axe X d'avancée du front. Découplage par la technique des perturbations La méthode selon l'invention estime le couplage entre les perturbations de vitesses Sic et la 'fluctuation de la position du front bh de façon serai analytique. En ce sens, la méthode utilise la théorie des perturbations. 20 D'un point de vue heuristique, la théorie des perturbations est une méthode mathématique générale qui permet de trouver une solution approchée d'une équation mathématique (El) dépendante d'un paramètre 1. lorsque la solution de l'équation (Eo), correspondant à la valeur X = 0, est connue exactement. L'équation mathématique (El) peut être une équation algébrique, une équation différentielle, une équation aux valeurs propres, . 25 La méthode consiste à chercher la solution approchée de l'équation (El) sous la forme d'un développement en série des puissances du paramètre 1, cette solution approchée étant supposée être une approximation d'autant meilleure de la solution exacte, mais inconnue, que la valeur absolue du paramètre 1 est plus petite . Ainsi, selon l'invention, la théorie des perturbations conduit à considérer localement la 30 perturbation de vitesse comme une somme de deux perturbations : Sf Swr q 15 13 2904982 Su(x, y, t) = Sua (x, y, t) + Sub (x, y, t) (15) Le premier terme, Sua décrit la perturbation due à la variation de saturation, et le second terme Su b décrit la perturbation due à l'hétérogénéité du milieu. Selon cette décomposition, on peut décrire la position du front au cours de son avancée 5 dans un milieu poreux hétérogène par l'équation suivante : a,Sh(q, t) = coA I q Sh(q, t) +a Subx (xf, q, t) (16) uo Mfù1 Avec: A= Mf +1 les fluctuations de la vitesse de filtration dans la direction d'avancée du front (axe X) l'abscisse de la position du front sur l'axe X Ainsi, à l'aide de cette formulation de la position du front, il est possible de découpler le problème. De cette façon, pour des champs de perméabilité à petites variances, la stabilité du front subissant de petites perturbations est déterminée par le rapport de mobilité frontal. Prise en compte des fluctuations transversales de la vitesse Selon l'invention, on simule la propagation du front dans un milieu hétérogène, c'est-à-dire dans un champ de vitesse hétérogène, en modélisant l'influence des effets visqueux de façon analytique. Le champ de vitesse hétérogène a été obtenu en utilisant une simulation 20 d'écoulement monophasique (étape 2 de la méthode). Pour simplifier les explications, on se place en deux dimensions. Le front séparant les deux fluides immiscibles est alors décrit par une fonction 2D, x(y, t) = h(y, t). On peut écrire cette équation sous la forme : F(x, y, t) = 0 (17) 25 En la résolvant par rapport à x, on a : FI (x, y,t) = h(y,t)ùx= 0 (18) SUbx 10 xi 15 14 2904982 Puis, en utilisant le fait que la surface en question est matérielle, i.e. composée d'éléments de fluide qui suivent le mouvement local du fluide, et en utilisant l'expression pour la vitesse de propagation du front pour les déplacements de type Buckley-Leverett, on obtient : 5 ath(y, t) = co (u . D çP) 19=o (19) En considérant maintenant l'hétérogénéité du champ de vitesse, on représente la vitesse comme une somme de deux termes. Le premier représente la moyenne, et le second les fluctuation dues aux hétérogénéités : u(r) = uo + 8u(r) 10 Avec : r = {x, y}, uo = {uo,0} et 8u={8ux,8uy} Les perturbations du champ de vitesse provoquent les perturbations du front h(y, t) = ho (t) + 8h(y, t) . En substituant cette décomposition dans l'équation (19) et en extrayant le terme moyen en utilisant la relation a,ho (t) = co, qui est valable pour les déplacement de type Buckley-Leverett, on obtient l'équation suivante, décrivant les 15 perturbations de la position du front provoquées par les hétérogénéités du milieu : a,8h(y,t)=û (8ux(xf,y)û8u~,(xf,y)ay8h(y,t))(20) 0 Ainsi, par application de la théorie des perturbations, de la théorie des surfaces matérielles, et en utilisant l'expression de la vitesse de propagation du front pour les déplacement de type Buckley-Leverett, on peut décrire la position du front dans un milieu 20 poreux hétérogène de la façon suivante (8h a été substituée par h dans cette expression par soucis de clarté) : arh(y,t) = coAJ 2TC j q I h(q,t)e`gy + û (8ux(xf,y)û8uy(xf,y)ayh(y,t)) (21) 0 Le premier terme (à gauche) est l'expression du couplage visqueux et peut être calculé avantageusement dans l'espace de Fourier. Il n'y a donc plus lieu d'estimer ce couplage par 25 une coûteuse résolution de système linéaire. Dans tout ce qui suit, les perturbations de vitesses 3u,, et 8uy peuvent maintenant être supposées découplées de l'équation en saturation. Le second terme (à droite) est obtenu à l'aide d'un simulateur d'écoulement monophasique (étape 2 de l'invention). 15 2904982 L'étape suivante de la méthode consiste à proposer une technique de résolution de l'équation du front complet. Estimation de la position du front 5 Un champ de vitesse hétérogène dans le cadre d'un écoulement monophasique, a été établit à l'étape 2 de la méthode selon l'invention par résolution numérique (à l'aide d'un simulateur d'écoulement monophasique) de l'équation en pression. Ensuite, l'équation (21) est discrétisée sur la grille discrétisant le milieu. Puis cette équation discrétisée est utilisée pour obtenir les perturbations du front pour chaque pas de 10 temps. On utilise alors un schéma numérique explicite, décrit ci-après. On introduit l'indexe i pour numéroter les mailles selon la direction de l'axe X, l'indexe j pour numéroter les mailles selon l'axe Y, et n pour numéroter les pas de temps. On note Ax la taille d'une maille selon la direction X, Ay la taille d'une maille selon la direction Y, et At le pas de temps. 15 La méthode de résolution de l'équation (21) discrétisée se présente alors comme suit. 1- Tout d'abord on effectue une transformée de Fourier rapide pour les perturbations du front. En effet, l'approximation par un développement en série de perturbation possède une formulation explicite simple dans l'espace de Fourier. Ainsi, selon la méthode, on effectue le calcul dans le domaine où celui-ci est le plus simple, puisque le coût d'une Transformée de 20 Fourier Rapide (FFT) est considéré comme négligeable. Un algorithme de Transformée de Fourier Rapide (FFT) en y permet d'estimer rapidement le premier membre de l'équation (21). 5h Ny 8h . ex 27i (k -1)(n -1) (22) p ù ) J=1 Ny / 2- Puis on multiplie les modes de Fourier des fluctuations de la position du front par un module du vecteur d'onde dans l'espace de Fourier, et l'on utilise une transformée de Fourier 25 inverse pour revenir dans l'espace réel. Pour obtenir une approximation de la dérivée de la vitesse selon y, on utilise la différence centrale, bien connue des spécialistes. On obtient donc le schéma explicite suivant : 16 2904982 Ny 8h;+' = an +c0AAt 2n-(k -1) Shk exp Ny k=1 NyAY Sh;+~ -.8h;_1 où l'index i est défini par la position du front au pas de temps précèdent n. Ici le terme Sh;"k représente la transformée de Fourier du front à la fréquence k, Su représente la 5 fluctuation de la vitesse longitudinale selon l'axe X (Sur), et Sv représente la fluctuation de la vitesse longitudinale selon l'axe Y (8uy). La transformation de Fourier discrète et son inverse peuvent être implémentée en utilisant un algorithme de transformée de Fourier rapide (FFT) pour accélérer les calculs. Ayant calculé les fluctuations du front pour un pas de temps donné Nt, on reconstruit la 10 position du front dans la grille de discrétisation du milieu à l'aide de la formule suivante : hN' = AtN co + 8hN' (24) On rappelle que l'approximation de la vitesse u a été réalisée une fois pour toute, puisque le couplage visqueux est modélisé via le terme de convolution (premier terme de 15 l'équation (21)). On peut donc déterminer de façon explicite la position du front pour chaque pas de temps, après avoir déterminer une seule fois, à l'instant t = 0, les vitesses des fluides par résolution de l'équation en pression à l'aide d'un simulateur d'écoulement. On notera que la technique proposée n'a besoin que d'une seule résolution numérique 20 complète de l'équation en pression, contrairement aux techniques du marché qui en demande autant qu'il y a de pas de temps. Une comparaison entre la méthode selon l'invention (figure 1A) et une simulation utilisant un simulateur à ligne de courant du marché (figure 1B) illustre que malgré l'approximation réalisée pour gagner en temps de calcul, la précision des résultats restent comparables aux autres méthodes. La figure lA montre le front de l'interface simulé 25 avec la méthode, tandis que la figure 1B montre le front de l'interface simulé à l'aide d'un + co At 8u"1ù 8v, i 2Ay Ny ( û 2'ri (k -1)(n -1) 8hk = E Sh' exp =1 Ny "2Tri(k -1)(n -1)1 Ny (23) 17 2904982 simulateur à ligne de courant du marché. Ces figures représentent une coupe du sous-sol, l'axe des abscisses correspond à une coordonnée géographique horizontale x, tandis que l'axe des ordonnées représente la profondeur y. 5 Optimisation de la récupération d'huile par injection A partir de la position du front, il est connu de déterminer les saturations en chaque maille du milieu et pour chaque pas de temps. On peut par exemple utiliser la méthode des tubes de courant (bien connue des spécialistes) loin du front, et interpoler les saturations près du front dont la position est maintenant bien déterminée. 10 En suivant l'évolution géographique et temporelle des saturations on peut déterminer à quel moment et où le fluide injecté pour faciliter la récupération va atteindre le puits. Lorsque l'on utilise de l'eau pour pousser l'huile en place d'un réservoir, la connaissance des saturations permet de déterminer le temps de percée de l'eau, qui une donnée clef de l'exploitation d'un champ pétrolier. 15 Du fait d'une détermination rapide des saturations, l'homme du métier est en mesure : 1 D'échantillonner rapidement l'espace des paramètres d'étude sur lesquels porte la modélisation, de façon à optimiser un ou plusieurs critères économiques d'exploitation (décision d'exploitation, choix de la position des puits, d'un procédé de récupération etc.). 2 De modifier de façon cohérente le modèle géologique afin de caler au mieux les 20 données de production observées, y compris des données de sismiques répétées. Ceci peut permettre de forer en s'adaptant aux hétérogénéités géologiques, de localiser les fluides, et donc de mieux contrôler le scénario de récupération. 3 D'estimer des incertitudes en effectuant des simulations de Monte Carlo, de façon à tester le rôle d'hétérogénéités, ou de paramètres mal connus caractérisant le sous-sol. Ceci 25 permet de quantifier le niveau de risque lié à l'exploitation d'un réservoir et donc de redéfinir des paramètres d'exploitation voire les paramètres économiques. Une étude exhaustive incorporant ces divers aspects (étude de sensibilité, calage aux données et estimation d'incertitudes) nécessite potentiellement de nombreux appels à un 30 simulateur d'écoulement (logiciel dédié). De ce fait, les solutions existantes ne permettent pas d'obtenir des réponses quantitatives à des questions pratiques que l'ingénieur pétrolier se pose 18 2904982 dès lors que l'on procède à des études de simulation de réservoir. D'où l'intérêt de la méthode selon l'invention qui propose un d'outil de calcul rapide, de précision compatible avec la précision des données d'entrée. 19 2904982 Les notations suivantes sont utilisées au cours de la description de l'invention : Vecteur vitesse du fluide non-mouillant Vecteur vitesse du fluide mouillant Vitesse totale des fluides. Mobilité du fluide non-mouillant Mobilité du fluide mouillant Pression du fluide non-mouillant Pression du fluide mouillant Tenseur perméabilité absolue du milieu Perméabilité relative du fluide non-mouillant Perméabilité relative du fluide mouillant Viscosité du fluide non-mouillant Viscosité du fluide mouillant Pression capillaire Saturation du fluide non-mouillant Saturation du fluide mouillant Pression commune des deux fluides, si pc(S) = O. Porosité (supposée uniforme dans le réservoir) Temps Saturation irréductible en eau Saturation en eau au front Rapport de mobilité totale au front. vecteur d'onde dans l'espace de Fourier Vitesse du front non perturbé Vitesse constante de référence Saturation maximale en eau Dérivée du flux fractionnaire en eau Position de référence du front Position du front Fluctuation de la position du front. Vecteur normal au front Vitesse totale moyenne. Perturbation du vecteur vitesse du fluide non-mouillant Perturbation du vecteur vitesse du fluide mouillant Perturbation du vecteur vitesse totale des fluides. Perturbation de la composante x de la vitesse totale. Perturbation de la composante y de la vitesse totale. F Equation du front. 40 FI Equation du front sous forme canonique LX, Ly Taille en x et en y du domaine. Vn Vitesse normale au front f(S) Flux fractionnaire en eau A~ Aire d'injection 45 Q Débit total injecté. V Volume du milieu poreux x Coordonnée longitudinale. y Coordonnée transverse. k Vecteur d'onde discrétisé 50 tD Temps adimensionné par le temps de percée moyen. Pnw Pw 10 K krnw krw rw A, 15 Pc(S) Snw Sn, p 0 20 t Swr Sf Mf q 25 co uo Swr Jw' ho 30 h(y, t) bh(y, t) n uo bunw buw bu 8ux buy..FT: METHODE POUR OPTIMISER LA RECUPERATION ASSISTEE D'UN FLUIDE EN PLACE DANS UN MILIEU POREUX PAR SUIVI DE FRONT. REVENDICATIONS 1) Méthode pour optimiser l'exploitation d'un milieu poreux hétérogène par injection dans ledit milieu d'un fluide de balayage pour provoquer un écoulement d'un fluide en place contenu dans ledit milieu, lesdits fluides étant immiscibles, dans laquelle on discrétise ledit milieu en une grille constituée d'un ensemble de mailles, la méthode étant caractérisée en ce qu'elle comporte les étapes suivantes : - on détermine une vitesse et une direction d'écoulement d'au moins un desdits fluides, en utilisant un simulateur d'écoulement pour résoudre une équation en pression ; - on définit une relation décrivant une position d'un front séparant lesdits fluides, en s'affranchissant d'un couplage visqueux à l'aide de la théorie des perturbations, et en prenant en compte des fluctuations de vitesse dans la direction d'avancée du front et des fluctuations de vitesse dans la direction perpendiculaire à celle d'avancée du front ; et pour différents pas de temps, - on reconstruit la position du front dans ladite grille à l'aide d'une discrétisation de ladite relation et d'une transformée de Fourier rapide ; et - on optimise l'injection dudit fluide de balayage en fonction de la position dudit front. 2) Méthode selon la revendication 1, dans laquelle ladite relation comporte un premier terme représentant un couplage visqueux décrivant la perturbation due à la variation de saturation, et un second terme décrivant les perturbations de la position du front provoquées par les hétérogénéités du milieu. 3) Méthode selon la revendication 2, dans laquelle ledit premier terme est obtenu en considérant ledit milieu homogène et les déplacements dudit front de type Bucley-Leverett. 4) Méthode selon la revendication 2, dans laquelle ledit second terme est obtenu en utilisant le fait que le front est une surface matérielle, et en représentant la vitesse comme une somme d'une vitesse moyenne avec des fluctuations de vitesses dues aux hétérogénéités. 5) Méthode selon la revendication 3, dans laquelle ledit premier terme est fonction d'au moins les paramètres suivants : un rapport de mobilité frontal, une vitesse de filtration moyenne le long d'une direction d'avancée du front, une porosité du milieu, une fonction de 2904982 21 Buckley-Leverett représentant un flux fractionnaire en eau, une saturation en eau au front, une saturation maximale en eau, le vecteur d'onde dans l'espace de Fourier. 6) Méthode selon la revendication 4, dans laquelle ledit second terme est fonction d'au moins les paramètres suivants : une vitesse du front non perturbé, une vitesse totale moyenne, 5 des perturbations des composantes de la vitesse totale dans la direction d'avancée du front et perpendiculaire à cette direction. 7) Méthode selon l'une des revendications précédentes, dans laquelle on optimise l'injection en déterminant en chaque maille et pour différents pas de temps la saturation d'au moins un desdits fluides à partir de la position dudit front.