FR 2898382 A1 20070914 FR 0602139 A 20060310 La présente invention concerne une méthode pour modéliser l'acidification au sein d'un milieu poreux suite à l'injection d'un produit chimique, tel qu'un acide. 10 En particulier l'invention permet d'optimiser des paramètres d'injection d'acide, tels que le débit et les zones à traiter, dans le cadre de stimulation acide de puits en contexte carbonaté. Dans l'industrie pétrolière, la production d'un puits peut être fortement réduite 15 suite à des endommagements au voisinage du puits. Les endommagements se traduisent par une altération de la perméabilité et de la nature de la roche autour du puits. Les opérations endommageant le puits sont multiples : forage, tubage, cimentation, exploitation, complétion et traitement. La conséquence de ces endommagements est le colmatage de la formation, et donc la réduction voir l'arrêt de la production 20 d'hydrocarbure. Il est donc très important pour l'industrie pétrolière, d'identifier d'une part le type d'endommagement, et d'autre part les zones endommagées, pour décider et élaborer un traitement adapté. L'un des traitements couramment utilisés dans l'industrie pétrolière est l'injection d'acide autour d'un puits. Cette injection permet de réduire les endommagements, et 25 donc, d'améliorer la production du puits. L'objectif premier de la stimulation acide est de réduire la résistance à l'écoulement des fluides du réservoir due aux endommagements. L'acide injecté dissout le matériel dans la matrice du réservoir et crée des chenaux qui augmentent la perméabilité de la matrice du réservoir. Ces chenaux sont d'autant plus fréquents dans les roches carbonatées, c'est-à-dire des roches contenant plus de 50% de minéraux carbonés (calcite, dolomie), telles que les calcaires. L'efficacité de ce procédé dépend du type d'acides utilisé, des conditions d'injections, de la structure du milieu, du transfert entre le fluide et le solide, du taux de réactions, etc. Alors que la dissolution augmente la perméabilité, l'augmentation relative de la perméabilité pour l'injection d'un volume d'acide donnée est observé comme étant fortement fonction des conditions d'injections. Dans les réservoirs gréseux, les fronts de réaction tendent à être uniforme et les chenaux d'écoulement ne sont pas observés. Dans les réservoirs carbonatés, en fonction des conditions d'injection, de multiples percées vermiculaires dans la roche, communément appelées wormholes , peuvent être créées. Présentation de l'art antérieur Il est donc très important pour l'industrie pétrolière, d'identifier d'une part le type d'endommagement, et (l'autre part les zones endommagées, pour optimiser les paramètres de la stimulation acide, dans le but de produire des percées vermiculaires ( whormholes ) avec une densité et une profondeur de pénétration dans la formation optimales. L'étude de la formation et du comportement des percées vermiculaires pour déterminer les paramètres de l'injection acide peut se faire selon quatre échelles différentes : - L'échelle du pore qui est l'échelle à laquelle nous décrivons les mécanismes de la réaction chimiques : - L'échelle de la carotte à laquelle apparaît l'instabilité des wormholes - L'échelle du puits qui est l'échelle à laquelle nous pouvons apprécier la compétition entre les wormholes et l'impacte des hétérogénéités à cette échelle. L'échelle du réservoir à laquelle l'effet de la stimulation est mesuré par le skin. Les figures lA à 1D, où le milieu a représente la roche et le milieu 13 l'eau et l'acide, illustre ces différentes échelles impliquées dans la stimulation acide : Figure lA : échelle du pore ( m - mm) Figure 1B : échelle de la carotte (mm - cm) Figure 1C : échelle du puits (cm-m) Figure 1D : échelle du réservoir (m -km) De nombreux modèles, tels que celui présenté par Wang, Y., Hill, A. D., et Schechter, R. S., The Optimum Injection Rate for MAtrix Acidizing of Carbonate Formations , papier SPE 26578, SPE ATCE, Houston, 1993, ont déjà été proposés pour étudier l'effet de fuite de fluides, de réaction cinétiques, etc. sur le taux de propagation des wormholes et l'effet des wormholes voisins sur le taux de croissance des wormholes dominant. La structure simple de ces modèles offre l'avantage d'étudier en détail la réaction, les mécanismes de diffusion et de convection au sein des percées vermiculaires. Ces modèles, cependant, ne peuvent pas être utilisé pour étudier l'initialisation de wormholes et les effets sur les hétérogénéités de la formation. Des modèles décrivant la dissolution lors d'injection d'acide ont été pour la première fois mis en plhce pour décrire ce phénomène à l'échelle du pore. Une telle méthode est décrite par exemple dans Hoefner, M. L., Fogler, H., S., "Pore Ecolution and Channel Formation During Flow and Reaction in Porous Media", AIChE J, 34, 45-54 (1998). Cependant la simulation à l'échelle de la carotte à partir de ces modèles est difficile et nécessite une grande capacité de calcul. Or c'est à cette échelle que les instabilités dues aux wormholes apparaissent. Le premier modèle à l'échelle de la carotte capable de reproduire totalement les mécanismes de la dissolution à été proposé par Golfier,F. et al., "A discussion on a Darcy-Scale modelling of porous media dissolution in homogeneous systems", Computational Methods in Water Resources, 2, 1195-1202 (2002). Ce modèle simple milieu est construit à partir d'une prise de moyenne volumique des équations à l'échelle du pore. Cette modélisation à été reprise dans la demande internationale WO 03/102362, qui a élargi le modèle au cas d'une dissolution limitée par la cinétique de réaction. Ces modèles se basent sur une description de la physique à l'échelle de la carotte, ce qui nécessite des tailles de mailles de l'ordre du millimètre. Cependant, un procédé d'injection acide est un procédé à l'échelle du puits. Il apparaît donc nécessaire de modéliser la formation et le comportement des wormholes à cette échelle. D'autant plus que, dans l'industrie pétrolière, la généralisation de puits horizontaux a engendré une augmentation des quantités d'acides injectées dans un seul puits. Les besoins en moyen de simulation afin d'augmenter les chances de réussite du traitement se sont accrus. Or, les modélisations décrites précédemment ne permettent pas de simuler l'acidification sur un domaine représentant la section d'un puits et ses environs (1 à 3m). Des modèles destines à simuler le traitement acide à une échelle supérieure à celle de la carotte ont déjà été proposés. On peut citer par exemple : - Buisje,M.A. Urderstanding Wormholing Mechanisms Can Improve Acid Treatments in Carbonate Formations. (SPE 38166). 1997. SPE European Formation Damage Conference. Buisje,M.A. & Glasbergen,G. (SPE 96892). 2005. SPE Annual Technical conference and Exhibition. Gdanski,R. A Fundamentally New Model of Acid Wormholing in Carbonates. (SPE 54719). 1999. SPE European Formation Damage Conference. Ces méthodes reposent sur des considérations empiriques basées sur des observations faites en laboratoire, très éloignées des conditions et des dimensions réelles. La méthode selon l'invention est une méthode pour modéliser à l'échelle métrique l'acidification au sein d'un milieu poreux suite à l'injection d'acide, permettant de répondre aux besoins des ingénieurs réservoir, pour définir un scénario adapté de stimulation acide de puits dans le cadre de réservoirs carbonatés. La méthode selon l'invention L'invention concerne une méthode pour modéliser l'acidification au sein d'un milieu poreux suite à l'in- ection d'un acide, dans laquelle on représente ledit milieu par un modèle double milieu, caractérisée en ce que la méthode comporte les étapes suivantes : a) on construit ledit modèle double milieu - en considérant un premier sous-milieu propice à des percées de dissolution, et un second sous-milieu r.on propice à des percées de dissolution ; - en réalisant, pour chacun desdits sous-milieux, une description à l'échelle métrique du transport de l'acide, de la conservation de la masse dudit sous-milieu et de la conservation de la masse dudit acide ; - en décrivant un transfert d'acide d'un sous-milieu à l'autre sous-milieu ; b) on initialise ledit modèle double milieu à partir de calages expérimentaux ; c) on modélise, à l'aide dudit modèle double milieu, ladite acidification, en déterminant des pararntres physiques représentatifs dudit milieu poreux et des paramètres physiques relatifs à l'acide injecté. Les paramètres physiques représentatifs du milieu poreux peuvent être choisis, pour chacun des sous-m:.lieux, parmi les paramètres suivants : la porosité moyenne, la perméabilité à l'échelle métrique et la pression totale moyenne. Les paramètres physiques relatifs à l'acide peuvent être choisis, pour chacun des sous-milieux, parmi les paramètres suivants : la concentration moyenne de l'acide, la vitesse de Darcy moyenne. Selon l'invention, la description peut être réalisée à l'aide d'équations obtenues en réalisant une moyenne volumique à l'échelle métrique d'équations décrivant la propagation d'un acide dans un modèle simple milieu à une échelle centimétrique. Ces équations comportent alors de préférence un terme de dissolution. Ce dernier peut être défini comme le produit d'une concentration moyenne de l'acide à l'échelle métrique par un coefficient dépendant d'une vitesse locale de l'acide. Il peut également être défini comme le produit d'un paramètre par la divergence d'un produit entre une concentration de l'acide. une fonction de flux fractionnaire et un vecteur vitesse. Le paramètre de ce dernier terme de dissolution peut dépendre d'une norme d'une vitesse locale de l'acide et de la ,porosité moyenne à l'échelle métrique. Selon l'invention, la procédure de calage peut être soit basée sur des simulations à des échelles inférieures à l'échelle métrique, soit basée sur des études d'injection de l'acide à débit constant dans un échantillon du milieu. Selon une réalisation, le milieu poreux peut être un réservoir carbonaté traversé par un puits, l'injection d'acide étant réalisée pour stimuler une production d'hydrocarbure par ledit puits, et dans laquelle on détermine des paramètres optimaux de l'injection acide, en réalisant les étapes suivantes : a) on construit un maillage dudit puits et de son voisinage ; b) on définit des paramètres initiaux de l'injection d'acide ; c) on détermine, en modélisant l'acidification due à l'injection de l'acide, au moins les paramètres physiques représentatifs dudit réservoir suivants : une porosité et une perméabilité dudit réservoir après injection de l'acide ; d) on simule la production du puits en fonction de ladite porosité et de ladite perméabilité à l'aide d'un simulateur réservoir ; e) on modifie lesdits paramètres initiaux et l'on réitère à l'étape c) jusqu'à obtenir un maximum de production. Selon cette réalisation, les paramètres initiaux peuvent être choisis parmi l'un au moins des paramètres suivants : le débit d'injection de l'acide, la vitesse initiale d'injection, le volume d'acide injecté, la concentration de l'acide utilisé lors de la stimulation, les zones à traiter. Présentation sommaire des figures - Les figures 1A à 1 D montrent les différentes échelles impliquées dans la stimulation acide. - La figure 2 illustre les différentes étapes de la méthode selon l'invention. - La figure 3 illustre les différentes étapes de calage du paramètre x. - La figure 4 montre la répartition des volumes VH et VM dans l'approche double milieu selon l'invention. - Les figures 5A à 5C montrent la représentation simplifiée de la répartition des volumes VH et VM utilisée dans la modélisation du terme d'échange. - La figure 6 illustre des résultats en utilisant les termes de dissolution Elle montre la chute de pression dans l'échantillon au cours du temps. - La figure 7 illustre des résultats en utilisant les termes de dissolution Elle montre la porosité de l'échantillon en fonction de l'abscisse. - La figure 8 illustre des résultats en utilisant les termes de dissolution g2M . Elle montre la chute de pression dans l'échantillon au cours du temps. - La figure 9 illustre des résultats en utilisant les termes de dissolution g2H et g2M . Elle montre la porosité de l'échantillon en fonction de l'abscisse. Description détaillée de la méthode La méthode selon l'invention permet de modéliser l'acidification d'un milieu poreux due à l'injection d'un produit chimique, tel qu'un acide. L'acidification comprend plusieurs phénomènes, dont les principaux sont : la dissolution du milieu par l'acide et le transport (propagation) de l'acide au sein du milieu. Pour ce faire on construit un modèle dotble milieu permettant de modéliser ces phénomènes à une échelle métrique. La description de l'invention est réalisée dans le cadre de la stimulation acide de puits de production. Cette stimulation consiste à injecter de l'acide autours d'un puits pour en augmenter la production d'hydrocarbure. La méthode, une fois appliquée à un domaine maillé représentant les abords d'un puits à stimuler, permet de simuler l'évolution de la porosité et de la perméabilité de la roche, et ainsi optimiser les paramètres de la stimulat: on acide tels que le débit d'injection et la zone de traitement, pour définir le scénario d'acidification optimum pour ce puits. La figure 2 illustre les différentes étapes de la méthode appliquée à l'injection acide autours d'un puits : gl H gui et g,M . et g,M. et g2H 1- Maillage du puits et de son voisinage (MAI) 2- Modélisation de l'acidification due à la stimulation acide (MOD -+ e, K) 3- Simulation de la Droduction du puits (SIM ---~ Prod) 4- Optimisation des paramètres d'injection d'acide (OPT) 5- Stimulation acide optimisée du puits pour augmenter sa productivité (STIM) 1- Maillage du puil:s et de son voisinage Pour permettre de modéliser les effets de la stimulation acide sur un puits et son environnement direct, on discrétise cet espace (puits + voisinage) à l'aide d'un maillage structuré de type radial. Ce type de maillage, bien connu des spécialistes, permet de prendre en compte les directions radiales des écoulements autours des puits, et donc d'améliorer la précision des calculs. 2- Modélisation de l'acidification due à la stimulation Cette étape nécessite dans un premier temps la définition d'un modèle de dissolution et propagation (acidification) permettant de modéliser la formation et le comportement de toutes les figures de dissolution à l'échelle du puits : front compact, wormhole conique, wormhole dominant, wormhole ramifié et dissolution uniforme. 3.1- Modèle d'acidification à l'échelle du puits Pour modéliser les phénomènes de dissolution et de propagation de l'acide à l'échelle du puits, on construit un modèle basé sur des équations de la mécanique des fluides et des lois de conservation de la masse de roche, de fluide et d'acide. Selon l'invention, ce modèle est un modèle double milieu, construit à partir d'une prise de moyenne volumique à l'échelle du puits des équations décrivant la propagation d'acide dans un modèle simple milieu à l'échelle de la carotte (cm-mm). Ces équations à l'échelle de la carotte ont été développées par Golfier,F.. et al., On The ability of a Darcy-scale mode/ to tapure wormhole formation during the dissolution of a porous media. , Journal of Fluid Mechanic. 547, 213-254. 2002. Le modèle selon l'invention permet ainsi d'utiliser un maillage radial dont l'extension radiale est de l'ordre du centimètre et du mètre. On peut donc simuler l'acidification à une échelle suffisamment grande pour reproduire l'ensemble d'un traitement acide et évaluer l'augmentation de la perméabilité aux abords du puits. La simulation de l'évolution. de la perméabilité permet ensuite de simuler la production puis d'optimiser les paramètres d'injection acide. Le modèle double milieu est défini en considérant que la roche du réservoir est constituée de deux milieux, notés H et M, de volumes respectifs VH et Vu, caractérisés par deux régimes de dissolutions différents. La figure 4 illustre ces deux milieux où Vsection représente le volume d'une section d'un puits, et la courbe noire représente les percées vermiculaires ( wormholes ). Selon l'invention, les deux régimes associés aux deux milieux sont définis de la manière suivante : - le volume V1 contient la forte densité de wormholes dont la croissance s'achève rapidement et qui sont de petites tailles (mm-cm) (courbe noire sur la figure 4). Ce milieu est représenta:if d'un régime compact dans lequel les wormholes ont une croissance courte ; - le volume VH contient les wormholes dominants, c'est-à-dire les wormholes pour lesquels la compétition s'étend sur des distances (dm-m) et des temps longs (figure 4). Ce milieu est propice au développement de wormholes et par conséquent est caractérisé par un front de dissolution rapide. Ainsi, le milieu H est un milieu propice à la formation de percées de dissolution : l'injection d'acide dans un tel milieu provoque la formation de percées (wormholes) de grandes tailles, en générale supérieures au décimètre. Le milieu M est un milieu non propice à la formation de percées de dissolution : l'injection d'acide dans un tel milieu ne provoque pas la formation de percées (wormholes) de grandes tailles, et permet, au mieux, la formation de percées de petites tailles, en générale inférieures au décimètre. Pour décrire les phénomènes d'acidification à l'échelle du puits, on réalise une prise de moyenne volumique à l'échelle du puits des équations décrivant la propagation d'acide dans un modèle simple milieu à l'échelle de la carotte. Ainsi, pour une variable quelconque 9, permettant de décrire l'acidification à l'échelle de la carotte, l'application du théorème de prise de moyenne permet de calculer les variables suivantes, permettant de décrire l'acidification à l'échelle du puits : {rP}H = V 1 f çpdV et {çp}M = V 1 J~dV Section V,1 Section VM {H 1 1 rp}H = JçpdV et {q'}MM = `o ;_ V VHVä VMVM Y'H VH et OM = ù~M = , soit OH = 1ù çbM VSection V Section {rp}H = 0H {p}H et {rp}A, _ 0M {çp}M avec Owh = {v} =OH {e}H +0M {r}, Le modèle d'acidification double milieu à l'échelle du puits selon l'invention 5 comporte pour chacun des deux milieux M et H : • une équation de transport de l'espèce acide composée : - de termes convectifs contenant une fonction de type flux fractionnaire permettant de reproduire en partie la dispersion liée au wormholing 10 - de te: mes réactifs - d'un terme d'accumulation • une équation de Darcy • une équation de conservation de la masse de roche • une équation de conservation de la masse fluide. 15 • Une équation de fermeture du système reliant la peunéabilité et la porosité. ( 1) ( 2) (3) Vwh VSection ( 4) ( 5) Le modèle double milieu selon l'invention s'écrit ainsi : Milieu H Équation du transport de l'espèce acide dans le milieu H H ac'H rH f rH - VrCH_M rM ù rH He at + v.( 0H V ) (P P )=gH' Equation de Darcy appliquée au milieu H VIN = ùK'x pp'H Équation de conservation de la masse de roche dans le milieu H at LV acgh' Équation de conservation de la masse fluide dans le milieu H v. (çby' !) ù w prM ùPH), o Équation reliant la perméabilité et la porosité dans le milieu H (eH ~ KH =K +(.ùK ) ù 1ùe i Milieu M Équation du transport de l'espèce acide dans le milieu M MeM aC'M }V.(~MVrMffM)+VCH (prM--P'H) ù MgM at Équation de Darcy appliquée au milieu M V'M = ùK'M.7P'M Équation de conservation de la masse de roche dans le milieu M aeM at NacgA1 Equation de conservation de la masse fluide dans le milieu M v.(IoMV1M ) + v/r( prM ù prH) = o Équation reliant la perméabilité et la porosité dans le milieu M rem ^x KM = K +(K. O. ù K ) 1ùe aeH (6) (7) (8) (9) (1o) ( 12) ( 13) ( 14) (15) 12 Avec les variables sJivantes utilisées pour adimensionnaliser le système : Km = 2898382 1H f Afl H H H CAR ù+(1- CAR H)/HH Co Co V t' = ù ùt L H C' H = CAR H Co f,H = f Co P'H pVoL )H K'H KH Ko g H L gH =- VoCo N _ ac Vo c'm = CAR Co 1 ~r VIM = IV~M_ Vo f M \ICAR ~M / {CAR 1M Co f fM Co P M _ L ter o ter K 0 H Pa ( 16) ( 17) (18) 1M ICAR M M ( 19) +Gù CAR m)/Hm Co (19') ( 20) (21) (22) ( 23) 13 Les sorties du modèle d'acidification sont : {CAJ f M = concentration moyenne d'acide dans le milieu M(Kg/m3) M {AQ ~H C = concentration moyenne d'acide dans le milieu H (Kg/m3) M { VA } = vitesse de Darcy moyenne dans le milieu M (m/s) {Vfl }H = vitesse de Darcy moyenne dans le milieu H (m/s) {P}H = pression totale moyenne dans le milieu M (Pa) {P}HI = pression totale moyenne dans le milieu H (Pa) {~ }M = M porosité moyenne dans le milieu M {e }H = eH porosité moyenne dans le milieu H KM == perméabilité dans le milieu M à l'échelle de la section (m2) KH == perméabilité dans le milieu M à l'échelle de la section (m2) Les entrées du modèle d'acidification sont : Vo = vitesse initiale d'injection de l'acide Co = concentration de l'acide utilisé lors de la stimulation eo = porosité initia: e = perméabilité initiale L longueur caractéristique du problème (rayon de la zone acidifié) ~c viscosité cinématique (Pa.$) v = coefficient stoechiométrique de la réaction de dissolution p0- = masse volumique de la roche(Kg/m3) Kf perméabilité dans le milieu dissout (m) ( 24) ( 25) (26) ( 27) ( 28) Ces données sont obtenues à partir de mesures diagraphiques, mesures sur carotte ou mesures en laboratoire. Ces paramètres peuvent également provenir de la connaissance géologique du spécialiste ou de simulations. La porosité et la perméabilité initiales sont ensuite opÉmisées au cours d'un processus d'optimisation basé sur la modélisation de l'acidification suite à une injection d'acide dans le puits. Kf correspond à la perméabilité dans le wormhole et sa valeur est donc très grande. Elle est calculée par analogie avec un écoulement de Poiseuille dans un wormhole. En prenant comme rayon b caractéristique du wormholes égale à 1 millimètre, nous obtenons : z Kf = b2 = 8,33.10-8 m2 en 2D (écoulement entre deux plaques parallèles) 2 (29) Kf = b = 1,25.107 m2 en 3D (écoulement dans un tube) 8 La précision de la valeur affectée à Kf n'a que peu d'importance dès l'instant où elle est très supérieure à Ko. Certains paramètres du modèle d'acidification doivent être déterminés avant la modélisation de l'acidification. Concernant les coefficient de dissolution gm et gH , deux formulations différentes peuvent être construites à travers deux approches différentes : g1 , g1H d'une part et g2M, g2H d'autre part. L'une est fonction de la concentration et de la porosité, et l'autre est dépendante de la vitesse, de la porosité et du bilan local du flux d'acide. Terme de dissolution g1H Lgim, La prise de moyenne volumique donne des termes de dissolution gH et gu non-linéaire qu'il est donc nécessaire de modéliser. Une première approche consiste à linéariser ces termes. On obtient ainsi les termes g1 et gui qui ne dépendent que du paramètre A.25 H glH == al y S( lCAp H M glM = alM {CAP )M Avec au, =A(1-H)2/3 z/3 (31) am A(1ùem) gui est le terme de dissolution pour le milieu H et gjj celui du milieu M. Cette expression des coefficients am et am a pour rôle de prendre en compte l'évolution de la surface de réaction par le biais de la variation de porosité. Après adimensionalisation on obtient : glH == V CyA(1ù y) Calage des paramètres des coefficients de dissolution Les paramètres utilisés dans notre modèle sont déterminés par une procédure de calage par rapport à des résultats de références couvrant une large gamme de débit. Ces débits doivent être choisis dans la gamme des débits pour lesquels des wormholes se forment. Ces résultats de références sont d'une part la porosité exacte à l'échelle de la carotte, moyennée à l'échelle du puits, et d'autre part le champ de pression, noté AP(t). Ceux-ci peuvent provenir soit d'expériences en laboratoire, telles que des études d'injection à débit constant dans un échantillon de roche, soit des simulations simple milieu à l'échelle de la carotte sur des domaines de petites tailles (échelle de la carotte). La détermination des paramètres à l'échelle du puits, permettant de reproduire les résultats obtenus à l'échelle de la carotte, est réalisée à travers une méthode d'inversion utilisant un algorithme de Levenberg-Maquart (K. Madsen, H.B. Nielsen, O. Tingleff, Methods for Non-Linear Least Squares Problems. 2004. Informatics and Mathematical Modelling. Technical University of Denmark. Cette procédure permet de déterminer les paramètres A, 0H HM, Hi Ay qui apparaissent dans le modèle double milieu utilisant les termes de dissolution glu et g,,, . De la même manière cette procédure permet de déterminer les paramètres ni, n2, y, 0H, HM le Ay qui apparaissent dans le modèle double milieu utilisant les termes de dissolution g2M et g2t/ . Ces déterminations de paramètres sont réalisées par rapport à des résultats de références couvrant une .arge gamme de débit. Pour un débit autre, on détermine la valeur affectée à chaque paramètre pour le modèle à l'échelle du puits, par une interpolation linéaire réalisée en comparant la vitesse à l'échelle de la section moyennée sur le volume Vsectjon, avec la vitesse d'injection des simulations simple milieu à l'échelle de la carotte. L"interpolation de ces valeurs en fonction du débit nous permet de réaliser des simulations à grande échelle sur un domaine de grande taille (échelle du puits). Calage du paramètre ,y Ce paramètre est également déterminé par une procédure de calage par rapport à des résultats de références couvrant une large gamme de débit choisis dans la gamme des débits pour lesquels des wormholes se forment. Ceux-ci peuvent provenir soit d'expériences en laboratoire, telles que des études d'injection à débit constant dans un échantillon de roche, soit des simulations simple milieu à l'échelle de la carotte sur des domaines de petites tailles(échelle de la carotte). Cette méthode de calage est illustrée sur la figure 3. Pour déterminer le coefficient x, on peut utiliser un calage par rapport à des résultats de simulations à l'échelle de la carotte à un débit constant (SimuC). Le débit imposé à un échantillon de la taille d'une carotte doit être choisi dans la gamme des débits pour lesquels des wormholes se forment. De ces résultats obtenus à l'échelle de la carotte, la porosité et le champ de pression (AP(t)) sont extraits, servant ainsi de résultats de références. Pour différents instants de la dissolution, la moyenne à l'échelle du puits, de la porosité obtenue à l'échelle de la carotte, est calculée. Cette nouvelle porosité eexact est appliquée à la relation reliant la perméabilité et la porosité décrite précédemment, K(s, x) : èxact ù 0 1ù ea ~ La perméabilité ainsi calculée représente la perméabilité moyenne (K) à l'échelle 15 du puits. Le champ de pression (P) induit par cette perméabilité est résolu en utilisant la relation suivante : O.(KVP) = 0 (38) Pour simplifier la détermination du paramètre x, on peut également utiliser des essais sur un écoulemen: linéaire dans un milieu homogène de manière à pouvoir résoudre l'équation (41) en ID de façon analytique en utilisant la relation suivante : P(L) -- P(0) = U,u 1ù dx K(x) 20 U correspond à la vi:esse d'injection. Ainsi on obtient ainsi la différence de pression entre les limites du domaine à différents instants, c'est-à-dire le gradient de pression aux bornes de l'échantillon (AP(t)exp). On calcule ensuite la différence entre ce gradient et le résultat de référence K =Ko+(KfùKo) (37) (39) AP(t). En fonction de l'erreur alors mesurée, le paramètre x est modifié (Ax) par conséquent. I)e manière itérative on obtient une valeur optimum de x qui minimise cette erreur. Cette valeur optimale de x est obtenue pour un débit donné. On répète l'opération pour différer.ts débits, choisi dans la gamme des débits pour lesquels des wormholes se forment. On utilise alors la relation K(s) ainsi paramétrée pour tous les débits lors des simulations à l'échelle du puits. Si le débit donné n'a pas servi à l'évaluation, on réalise une interpolation de valeur pour les débits utilisés encadrant le débit donné. Calage du paramètre d'échange yr entre les deux milieux H et M Les deux milieux If et M interagissent par le biais d'un terme d'échange dépendant de la différence de pression entre ces deux milieux. Ce terme permet de modéliser le détournement du flux d'acide vers les wormholes dominants au détriment de ceux qui se trouve dans le milieu M. Pour déterminer le terme d'échange yr on applique le modèle à un cas particulier. On représente le volume par un milieu dans lequel on injecte linéairement de l'acide. Des wormholes de forme; cylindriques s'y développent, disposés de manière périodique en fonction de leur taire. L'équivalence entre cette représentation et la réalité est assurée par un paramètre 4y qui doit être déterminé par calage. Ay défini ici la distance entre deux wormholes dominants. Les figures 5A, 5B et 5C montrent la représentation simplifiée de la répartition des volumes VH et Vy utilisé dans la modélisation du terme d'échange : la figure 5A montre la figure de dissolution réelle, la figure 5B illustre la représentation simplifiée, et la figure 5C illustre le motif de base. Cette représentation périodique permet de représenter les termes d'échanges pour tout le domaine à partir de sa description sur un motif de base. Dans cette description le volume Vsection contient n fois le motif de base. V 1 fVf.nMds= 1 1 j ûK .OP .nMds section Ay_y section A"_M 1 Keg y aP Ax.n.Ay , f2 ay On évalue le gradient de pression à l'interface en divisant l'écart des pressions moyennes des deux milieux par la hauteur dy/2 du motif de base (Figure 5C). La perméabilité équivalente Keu, est une variable calculée en faisant une moyenne harmonique des perméal: ilités transversales (Ky H et Ky m) des deux milieux. 1 4Keg _ y \{P}M ù IP}N ) V ù f V/3.nMds û z section AH_M y Avec Ky H •Ky M K = eg_y Y'MKyH +OHKY_H Il faut maintenant définir les perméabilités transversales. Pour cela on utilise une représentation idéale de chaque milieu, en les modélisant comme des blocs traversés par une certaine quantité de wormholes de section constante. En appliquant la loi de Darcy à cette représentation pour déterminer Ky m et Ky H,, on obtient : K0 .K f EH û s0 Ky û AH )K f. + AH Ko ' AH =0- (1 1û so (40) 2.n.Ax (41) (42) (43) K K 0' f (1û AM )K f + AM Ko , AM Sm ù E0 1 û e0 Ky_M (44) On peut enfin écrire le terme r/sous la forme suivante, en fonction du paramètre dy. K eq_y PQyz Calage de la concentration à l'interface entre les deux milieux CH-M En ce qui concerne la concentration CHM utilisée avec le terme d'échange dans le modèle double milieu, on utilise soit la concentration C'M soit la concentration C'H 5 selon les valeurs de P'M et P'H : si p 'm > p 'N alors CH_M = C'M si P M Exemple d'application Selon un exemple d'application de la méthode selon l'invention, on simule à l'échelle du puits une irjection d'acide à débit constant sur un échantillon de roche d'une longueur de 2m, 4Ccm de largeur et 40cm de hauteur. Maillage et initialisation : Après avoir maillé l'échantillon à l'aide d'un maillage cartésien (dans ce cas d'application le maillage est cartésien et non radial comme dans le cas d'un puits par exemple), on détermine ou définit les données d'entrée du modèle : - Vitesse initiale d'injection Vo = 1.10.4 m/s - Concentration initiale Co = 210 Kg/m3 Porosité initiale s0 = 0,36 Perméabilité initiale Ko = 2.318.10-12 m2 -Viscosité cinématique ,u = 1.10 Pals - Masse volumique de la roche pa = 2160 Kg/m3 - Perméabilité dans le milieu dissout Kf= 8,331.10-8 m2 Coefficient stoechiométrique massique v = 1 - Longueur caractéristique du problème L = 0,1 m Détermination expérimentale des paramètres : On détermine des résultats de références à l'aide de simulations à l'échelle de la carotte grâce au modèle développé par Golfier,F. et al., On The ability of a Darcyscale model to tapure wormhole formation during the dissolution of a porous media. , Journal of Fluid Mecharic. 547, 213-254. 2002. Ces simulations encadrent le débit utilisé ultérieurement pour la simulation à l'échelle du puits. Une série de simulations est ainsi réalisée en utilisant un modèle simple milieu à l'échelle de la carotte sur un domaine qui représente une petite portion du domaine à simuler, sur une gamme de débit suffisamment large pour reproduire les différents types de figures de dissolution possibles (front compact, wormhole conique, wolinhole dominant, wormholes ramifiés 23 15 et dissolution uniforme). Les dimensions du domaine sont 25 cm de longueur, 40cm de hauteur et 1mm de hauteur. Pour déterminer le coefficient x lié à la relation perméabilité/porosité, on utilise les résultats de pression 2t de porosité issus des simulations à l'échelle de la carotte auxquels on applique la procédure détaillée précédemment. On obtient ainsi la valeur optimale x = 3,08. On utilise ensuite le modèle à l'échelle du puits sur un domaine équivalent à celui utilisé à petite échelle, el appliquant les mêmes conditions d'injection. Pour chaque simulation à l'échelle de la carotte à un débit donné, on utilise un algorithme d'optimisation afin de déterminer les paramètres équivalents utilisés par notre modèle à l'échelle du puits pour chacun de ces débits. Les résultats sont présentés dans les tableaux 1 et 2. Le tableau 1 montre les valeurs des paramètres du modèle utilisant les termes de dissolution g1 et gui pour différentes vitesses d'injections. Vo (mis) A OH Ay (m) Hm HH 9,27E-08 0,0023 0,5 0,2 1 1 4,64E-06 0,001 0,05 0,0041 1,1 1,48 9,27E-06 0,0023 0,08459 0,005 1,1939 1,461 2,32E-05 0,04 0,08459 0,005884 1,16373 1,2 4,64E-05 0,04 0,0769 0,006 1,13 1,207 9,27E-05 0,04 0,07 0, 006 1,13 1,3 1,85E-04 0,04773 0,0697 0,01 1 1,226 9,27E-04 0,3936 0,136 0,1 1,0012 2,99 9,27E-03 7 0,2 0,1366 1,5105 3 Le tableau 2 montre les valeurs des paramètres du modèle utilisant les termes de dissolution g2M et g211 pour différentes vitesses d'injections. Vo (mis) ni OH Ay (m) HM HH y na 9,27E-08 0 0,5 10 1 1 0,667 0 4,64E-06 0,265 0,499 0,1783 3,8 4,3 0,755 0,83 9,27E-06 0,345 0,499 0,143 2,9722 3,962 0,71 0,8 2,32E-05 0,345 0,5 0,0938 2,5 3,5 0,68 0,55 4,64E-05 0,339 0,5 0,04859 2,47 3,35 0,667 0,5 9,27E-05 0,3 0,5 0,0344 2,479 _ 3,3 0,705 0,45 1,85E-04 0,2918 0,48 0, 0678 2,522 _ 3,15 0,8 0,35 9,27E-04 0,149 0,4537 0,1157 2,5742 2,8547 0,8858 0,2902 9,27E-03 0,0454 0,4148 0,1852 2,797 2,623 0,902 0,2728 Pour déterminer les paramètres à l'échelle du puits à utiliser pour l'injection d'acide, une interpolation des valeurs obtenues dans l'étape précédente est réalisée. Pour une vitesse d'injection de 1.10-4 mis, on obtient les paramètres suivants : - Avec les termes de dissolution glm et gui : A = C.0406 ÇH = 0.07 HM = 1.192 HH = 1.447 Ay = 0.006314 Avec les termes de dissolution g2M et g2H ni = 0.299 n2 = 0.4421 = 0.7124 OH = 0.498 Hm = 2.482 HH = 3.288 Ay = 0.037 Modélisation de l'acidification à l'échelle du puits : Pour modéliser l'acidification, on utilise le modèle double milieu selon l'invention (équations 6 à 15). L'échantillon est homogène en porosité et perméabilité à l'échelle de la section. Pour cette raison le modèle à l'échelle de la section est appliqué sur une seule dimension, dans le sens et la direction de l'injection. Les figures 6 et 8 montrent la différence de pression entre l'entrée et la sortie de l'échantillon, AP exprimé en pascal, en fonction du temps t exprimé en seconde. Les figures 7 et 9 montrent la porosité E de l'échantillon en fonction de l'abscisse x de l'échantillon en mètre. Les figures 6 et 7 illustrent les résultats pour le modèle utilisant les termes de dissolution g1 et g,H . Le temps de percée est de 4 heures et 2 minutes. Le volume d'acide injecté est de 2.677.10-1 m3. Les figures 8 et 9 illustrent les résultats pour le modèle utilisant les termes g2M et g2H . Le temps de percée est de 3 heures et 53 minutes. Le volume d'acide injecté est de 2.244.10-1 m3. Les deux modèles montrent une forte chute de la pression pour une faible augmentation de la porosité, ce qui est caractéristique du wonnholing. Ils montrent également qu'il faut environ 4 heures d'injection à la vitesse d'injection 1.10-4 mis pour obtenir un wormhole d'une longueur de deux mètres, longueur caractéristique d'une stimulation de puits par injection d'acide. REVENDICATIONS 1) Méthode pour modéliser l'acidification au sein d'un milieu poreux suite à l'injection d'un acide, da:as laquelle on représente ledit milieu par un modèle double milieu, caractérisée en ce que la méthode comporte les étapes suivantes : a) on construit ledit modèle double milieu - en considérant un premier sous-milieu propice à des percées de dissolution, et un second sous-milieu non propice à des percées de dissolution ; - en réalisant, pour chacun desdits sous-milieux, une description à l'échelle métrique du transport de l'acide, de la conservation de la masse dudit sous-mileu et de la conservation de la masse dudit acide ; - en décrivant un transfert d'acide d'un sous-milieu à l'autre sous-milieu ; b) on initialise ledit modèle double milieu à partir de calages expérimentaux ; c) on modélise, à l'aide dudit modèle double milieu, ladite acidification, en déterminant des paramètres physiques représentatifs dudit milieu poreux et des paramètres physiques relatifs à l'acide injecté. 2) Méthode selon la revendication 1, dans laquelle lesdits paramètres physiques représentatifs dudit milieu poreux sont choisis, pour chacun desdits sous-milieux, parmi les paramètres suivants : la porosité moyenne, la perméabilité à l'échelle métrique et la pression totale moyenne. 3) Méthode selon l'une des revendications précédentes, dans laquelle lesdits paramètres physiques relatifs à l'acide sont choisis, pour chacun desdits sous-milieux, parmi les paramètres suivants : la concentration moyenne de l'acide, la vitesse de Darcy moyenne. 4) Méthode selon :'une des revendications précédentes, dans laquelle ladite description est réalisée à l'aide d'équations obtenues en réalisant une moyenne volumique à l'échelle métrique d'équations décrivant la propagation d'un acide dans un modèle simple milieu à uni échelle centimétrique. 5) Méthode selon la revendication 4, dans laquelle lesdites équations comportent un terme de dissolution. 6) Méthode selon la revendication 5, dans laquelle ledit terme de dissolution est défini comme le produit d'une concentration moyenne de l'acide à l'échelle métrique par un coefficient dépendant d'une vitesse locale de l'acide. 7) Méthode selon la revendication 5, dans laquelle ledit terme de dissolution est défini comme le produit d'un paramètre par la divergence d'un produit entre une concentration de l'acide, une fonction de flux fractionnaire et un vecteur vitesse. 8) Méthode selon la revendication 5, dans laquelle ledit paramètre dudit terme de dissolution dépend d'une norme d'une vitesse locale de l'acide et de la porosité moyenne à l'échelle métrique. 9) Méthode selon l'une des revendications précédentes, dans laquelle ladite procédure de calage est basée sur des simulations à des échelles inférieures à l'échelle métrique. 10) Méthode selon l'une des revendications précédentes, dans laquelle ladite procédure de calage est basée sur des études d'injection de l'acide à débit constant dans un échantillon dudit milieu. 11) Méthode selon l'une des revendications précédentes, dans laquelle ledit milieu poreux est un réservoir ca:-bonaté traversé par un puits, l'injection d'acide étant réalisée pour stimuler une production d'hydrocarbure par ledit puits, et dans laquelle on détermine des paramètres optimaux de l'injection acide, en réalisant les étapes suivantes: a) on construit un maillage dudit puits et de son voisinage ; b) on définit des paramètres initiaux de l'injection d'acide ; c) on détermine, er. modélisant l'acidification due à l'injection de l'acide, au moins les paramètres physiques représentatifs dudit réservoir suivants : une porosité et une perréabilité dudit réservoir après injection de l'acide ; d) on simule la production du puits en fonction de ladite porosité et de ladite perméabilité à l'aide d'un simulateur réservoir ;e) on modifie lesdits paramètres initiaux et l'on réitère à l'étape c) jusqu'à obtenir un maximum de production. 12) Méthode selon la revendication 11, dans laquelle lesdits paramètres initiaux sont choisis parmi l'un au moins des paramètres suivants : le débit d'injection de l'acide, la vitesse initiale d'injection, le volume d'acide injecté, la concentration de l'acide utilisé lors de la stimulation, les zones à traiter.