FR 2951555 A1 20110422 FR 0905050 A 20091021 10 La présente invention concerne le domaine de l'exploitation d'un milieu poreux souterrain, tel que la production d'hydrocarbures issus d'un gisement pétrolier, ou l'injection de gaz acides dans le sous-sol. 15 Que ce soit dans le cadre de production d'hydrocarbures, ou dans le cadre d'injection de gaz, il est indispensable pour réussir cette exploitation, de pouvoir suivre l'évolution des fluides au sein du milieu poreux. Par exemple, il est nécessaire de déterminer où se trouve les hydrocarbures dans le gisement exploité, puis une fois en production, de suivre 20 l'évolution des hydrocarbures pour déterminer comment améliorer la production. Par exemple, si une partie des hydrocarbures reste bloquée, il peut être décidé de forer d'autres puits de production ou des puits d'injection pour améliorer la récupération de ces hydrocarbures. Dans le cadre de l'injection de gaz acides, tels que du CO2, il est indispensable de 25 suivre l'évolution du gaz injecté pour déterminer si le stockage s'effectue correctement, à l'endroit souhaité, évaluer la quantité effectivement stockée, et vérifier l'intégrité du réservoir en s'assurant qu'il n'existe aucune fuite. A partir de ces informations, on peut également déterminer si d'autres puits d'injection sont nécessaires. Pour suivre l'évolution des fluides produits ou injectés au sein d'un milieu poreux, de 30 nombreuses techniques ont été développées par les industriels. Parmi ces techniques, la sismique répétitive, dite sismique 4D, est utilisée dans l'industrie (pétrolière ou environnementale). Une telle technique consiste à réaliser différentes campagnes sismiques, à des moments différents (en général les campagnes sont espacées d'au moins une année). Ainsi, le spécialiste peut suivre l'évolution des fluides du réservoir en production ou du site géologique de stockage. On détermine alors si la variation de vitesses sismiques (ou impédances) observées sur les différentes sismiques acquises à différents temps calendaires, est principalement due à une variation de pression (fluide et confinement) ou à une substitution de fluide. En d'autres termes, on détermine comment les répartitions de fluides dans le réservoir varient avec le temps. Cette évolution des fluides produits ou injectés au sein du milieu poreux se caractérise, à partir de deux jeux de données sismiques, par deux paramètres : la variation du module d'incompressibilité du fluide saturant, et la variation de la pression de pore (ou pression de fluide). Le premier paramètre indique le type de fluide en un point donné du milieu, le second indique l'augmentation de pression induite. Ces deux paramètres sont généralement déterminés à partir des données sismiques inversées (vitesses ou impédances) : la loi de Biot-Gassmann relie le module d'incompressibilité du fluide saturant aux impédances ou vitesses sismiques, et la loi de Hertz-Mindlin relie la pression de pore aux impédances ou vitesses sismiques. La loi de Hertz-Mindlin est une loi en puissance de type : V = k.Ph, où V représente la vitesse, P la pression, k un réel, et h un exposant appelé coefficient de Hertz. Le coefficient de Hertz peut avoir des valeurs très diverses selon la nature des roches et la porosité. Les valeurs sont généralement comprises entre 0,005 et 0,05 pour les calcaires, tant pour les ondes P que pour les ondes S, tandis que pour les grès les valeurs croissent avec la porosité de 0,04 à 0,06 pour les ondes P et de 0,06 à 0,1 pour les ondes S. Le rapport des coefficients de Hertz pour les ondes S et P tourne autour de 1,5 pour les grès. Ce coefficient est déterminé à partir de mesures ultrasonores (typiquement de quelques centaines de kilohertz au mégahertz) sur échantillons de roche au laboratoire, ou à partir de données de diagraphies acoustiques (typiquement de quelques kilohertz à quelques dizaines de kilohertz). Les documents suivants décrivent cette procédure : Rasolofosaon, P.N.J, Zinszner, B.E., 2006, Complete Pressure and Fluid Dependence of Seismic Properties of Reservoir Rocks, 68th EAGE conference and exhibition, paper PO46 Rasolofosaon P., Zinszner B., 2004, Laboratory petroacoustics for seismic monitoring feasibility study. The Leading Edge, v. 23; no. 3, p. 252-258. - Rasolofosaon, P.N.J. and Zinszner, B.E., 2007. The unreasonable success of Gassmann's theory ... Revisited, Journal of seismic Exploration, Volume 16, Number 2-4, 281-301. Cependant, la bande de fréquence des données sismiques est typiquement de quelques hertz à quelques dizaines de hertz. Cette différence de fréquences entre les ondes utilisées pour déterminer la loi vitesse/pression et les ondes utilisées lors des campagnes sismiques, peuvent être à l'origine d'interprétations erronées, comme l'illustre la figure 3 (décrite ultérieurement). 10 L'objet de l'invention est une méthode alternative pour exploiter un milieu poreux souterrain, dans laquelle on interprète des données de sismique 4D au moyen d'une loi de dépendance des vitesses sismiques (ou impédances) en fonction de la pression de pore, prenant en compte la bande de fréquences des données sismiques. La méthode y parvient notamment en corrigeant le coefficient de Hertz au moyen d'une estimation de l'atténuation 15 des ondes utilisées pour définir ce coefficient de Hertz. '- La méthode selon l'invention L'invention concerne une méthode pour exploiter un milieu poreux souterrain par injection ou production d'un fluide, dans laquelle on acquiert au moins deux jeux de données 20 sismiques en réalisant au moins deux campagnes sismiques réalisées à des moments différents par émission d'ondes sismiques dans le milieu. On détermine deux jeux de vitesses des ondes sismiques en inversant les jeux de données, et on considère une loi de puissance de type Hertz-Mindlin reliant au moyen d'une puissance appelée coefficient de Hertz, une vitesse de propagation d'ondes à une pression de pore. La méthode comporte les 25 étapes suivantes: on détermine ledit coefficient de Hertz au moyen de mesures sur des échantillons dudit milieu en utilisant des ondes de fréquences f1; on corrige ledit coefficient de Hertz pour tenir compte d'une différence de fréquence entre lesdites ondes de fréquence f1 et lesdites ondes sismiques, au moyen d'une 30 mesure d'atténuation des ondes de fréquence f1 dans lesdits échantillons; on détermine en tout point du milieu poreux une variation de pression de pore au moyen de ladite loi utilisant ledit coefficient de Hertz corrigé et desdits jeux de données inversés, et on suit l'exploitation dudit milieu poreux souterrain par injection ou production d'un fluide, en déduisant desdites variations de pression de pore une localisation dudit fluide et éventuellement une substitution dudit fluide au sein du milieu poreux. On peut corriger le coefficient de Hertz au moyen de la formule suivante : n (f2,Pdff)-n (fi,Pâ=ff) Y(Pdiff)ln(f fJ + [no(.f2)s(f2)-no(fi)S(.fi)lxPdff avec : f2 : fréquence centrale des enregistrements sismiques; Pdl ff : pression différentielle à l'état X (=A ou B); n(f,Pdff) : coefficient de Hertz à la fréquence f et à la pression différentielle Pdiff n (f2,Pd .) étant le coefficient de Hertz corrigé et n (f,Pd .) le coefficient de Hertz non corrigé; Y atténuation des ondes P et S aux fréquences fj; un réel; no : coefficient de Hertz aux pressions différentielles tendant vers zéro, no (f2) étant le coefficient de Hertz corrigé et no (f,) le coefficient de Hertz non corrigé aux pressions différentielles tendant vers zéro. Selon un mode de réalisation préféré, on corrige le coefficient de Hertz au moyen de la formule suivante : no(f2)ùno(f) çr(Pdff)ln(f ' ) avec : f2 fréquence centrale des enregistrements sismiques; no : coefficient de Hertz aux pressions différentielles tendant vers zéro, no (f2) étant le coefficient de Hertz corrigé et no(f) le coefficient de Hertz non corrigé; P. : pression différentielle dans ledit milieu; y : atténuation des ondes P et S aux fréquences fi; : un réel. Enfin, selon l'invention, les fréquences f1 peuvent être des fréquences ultrasonores ou 5 soniques. D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Présentation succincte des figures La figure 1 illustre, sur une échelle bi-logarithmique, des valeurs expérimentales de l'atténuation y de l'onde P (losanges noirs) et de l'onde S (triangles noirs), en fonction de la pression différentielle (Pdff) mesurées au laboratoire. Les deux ajustements linéaires 15 sont également représentés en trait continu pour l'onde P et en trait discontinu pour l'onde S. Cas du grès de Fontainebleau. - La figure 2 illustre les vitesses expérimentales ultrasonores (courbes repérées US pour 'ultrasons') et soniques (courbes repérées SO pour 'sonic') de l'onde P (carré et triangle) et de l'onde S (losange et rond) en fonction de la pression différentielle (Pd ff), en échelle 20 bi-logarithmique. Les prédictions théoriques des vitesses P et S soniques (en traits discontinus) sont comparées aux vitesses soniques expérimentales. Cas du grès de Fontainebleau. - La figure 3 illustre les résultats d'interprétation selon la méthode. 25 Description détaillée de la méthode L'invention concerne une méthode pour exploiter un milieu poreux souterrain, dans laquelle on acquiert et interprète des données de sismique 4D. L'interprétation permet de suivre l'évolution des fluides produits ou injectés au sein d'un milieu poreux. L'idée consiste à calculer une variation de vitesses sismiques (ou impédances) observées sur les différentes 30 sismiques acquises à différents temps calendaires, puis de déterminer si cette variation est principalement due à une variation de pression (fluide et confinement) ou à une substitution 10 de fluide. En d'autres termes, on détermine comment les répartitions de fluides dans le réservoir varient avec le temps. Pour ce faire, on utilise la loi de Biot-Gassmann, et une loi de dépendance des vitesses sismiques (ou impédances) en fonction de la pression de pore, prenant en compte la bande de fréquences des données sismiques. La méthode y parvient notamment en corrigeant le coefficient de Hertz au moyen d'une estimation de l'atténuation des ondes utilisées pour définir ce coefficient de Hertz. La méthode peut comporter les étapes suivantes: 1. Acquisition et inversion de sismique répétitive pour imager le réservoir souterrain. 2. Établissement de lois de dépendance vitesses sismiques / pression de pore 2a. Détermination de coefficient de Hertz à partir d'ondes ultrasonores ou soniques 2b. Détermination d'une relation entre coefficient de Hertz et fréquence des ondes 2c. Utilisation de la relation entre coefficient de Hertz et fréquence des ondes 3. Interprétations des images sismiques 1. Acquisition et inversion de sismique répétitive pour imaqer le réservoir souterrain. La sismique-réflexion a pour but d'obtenir une image du sous-sol. Elle consiste à créer un ébranlement à la surface du sol et à enregistrer les très petits déplacements à la surface du sol entraînés par les réflexions de l'onde sur les surfaces de discontinuité des caractéristiques mécaniques (dureté ou compressibilité) des roches. La sismique répétitive, dite sismique 4D, est utilisée dans l'industrie (pétrolière ou environnementale). Une telle technique consiste à réaliser différentes campagnes sismiques, à des moments différents (en général les campagnes sont espacées d'au moins une année). Ainsi, le spécialiste peut suivre l'évolution des fluides du réservoir en production ou du site géologique de stockage, par comparaison des images issues de chaque campagne. Selon l'invention, on acquiert au moins deux jeux de données sismiques en réalisant au moins deux campagnes sismiques réalisées à des moments différents par émission d'ondes sismiques dans le milieu. Puis on transforme ces jeux d'amplitudes sismiques en impédances sismiques ou en vitesses sismiques. Pour ce faire on utilise une technique bien connue des spécialistes : l'inversion stratigraphique. On peut de façon avantageuse utiliser une inversion avant sommation. 2. Établissement de lois de dépendance vitesses sismiques / pression de pore 2a. Détermination de coefficient de Hertz à partir d'ondes ultrasonores ou soniques Les lois utilisées sont les lois de Hertz-Mindlin. Il s'agit d'une loi en puissance de type : V = k.Ph, où V représente la vitesse, P la pression, k un réel, et h un exposant appelé coefficient de Hertz. Le coefficient de Hertz peut avoir des valeurs très diverses selon la nature des roches et la porosité. L'établissement de lois de dépendance vitesses sismiques / pression de pore consiste à calibrer le coefficient de Hertz pour différentes valeurs de saturation en fluide. De façon discrète, on peut avoir une courbe pour un échantillon totalement saturée en eau, une autre pour un échantillon totalement saturée en gaz ou huile. Les documents suivants décrivent cette procédure : Rasolofosaon, P.N.J, Zinszner, B.E., 2006, Complete Pressure and Fluid Dependence of Seismic Properties of Reservoir Rocks, 68th EAGE conference and exhibition, paper PO46 Rasolofosaon P., Zinszner B., 2004, Laboratory petroacoustics for seismic monitoring feasibility study. The Leading Edge, v. 23; no. 3, p. 252-258. Rasolofosaon, P.N.J. and Zinszner, B.E., 2007. The unreasonable success of Gassmann's theory ... Revisited, Journal of seismic Exploration, Volume 16, Number 2-4, 281-301. A partir de mesures ultrasonores de laboratoire, on mesure sur les échantillons réservoir considérés, les vitesses/impédances avec différents fluides en fonctions de la pression de pore Pp (cf. Rasolofosaon et Zinszner, 2004, 2006 et 2007). On trace alors des courbes de calibration telles que celles de la figure 3. Les courbes de calibration peuvent être également obtenues par des mesures de diagraphies acoustiques, aux fréquences soniques et non plus aux fréquences ultrasonores comme au laboratoire. Les courbes de calibrations sonic sont évidemment moins détaillées qu'au laboratoire, car on ne dispose que de quelques mesures de vitesses P et S, sous différentes pressions de fluide dans le réservoir qui sont classiquement mesurées simultanément par d'autres diagraphies (par exemple PLT etc...). Ces courbes de calibration sont donc obtenues soit aux fréquences soniques, soit aux fréquences ultrasonores, et non aux fréquences sismiques. Cependant, une courbe de calibration dépend de manière importante de la fréquence (courbes en traits continus comparées aux courbes en traits discontinus sur la figure 3). 2951555 2b. Détermination d'une relation entre coefficient de Hertz et fréquence des ondes Pour corriger cet effet de fréquence, selon l'invention, on modifie ces courbes par une mesure de l'atténuation ultrasonore des ondes P et S, dans le cas de calibration aux fréquences ultrasonores au laboratoire, et de l'atténuation sonique des ondes P et S, dans le 5 cas de calibration aux fréquences soniques avec les diagraphies. Ces mesures conventionnelles d'atténuations sont assez connues dans la littérature et détaillées dans différentes références (voir par exemple, pour les mesures d'atténuations ultrasonores: Bourbié, T., Coussy, O., Zinszner, B., 1987, Acoustics of porous media, Editions Technip, Paris, et pour les mesures d'atténuation soniques : Serra, O., 2008, The Weil Logging 10 Handbook, Editions Technip, Paris). Concrètement, on définit une relation, permettant de modifier le coefficient de Hertz déterminé précédemment (étape 2a) en fonction de l'atténuation des ondes utilisées pour déterminer ce coefficient (de fréquence f1) et de la fréquence sismique (f2). Selon un mode préféré de réalisation, la formule utilisée est : 15 no (f2) - no (A) Ç Y(Pd1Î) ln\.f2 fl / Les formules de passage des vitesses aux fréquences ultrasonores aux vitesses aux fréquences sismiques, utilisant les atténuations P et S, sont détaillées ci-après. Les notations suivantes sont utilisées dans la description de l'invention. 20 Pd~ ff : pression différentielle dans un réservoir, différence entre la pression de confinement et la pression de pore. Pion f : pression de confinement, pression des sédiments pression de pore, pression régnant dans le fluide contenu dans le réservoir: : fréquence 25 P/ h& : pression différentielle à l'état X (=A ou B), n (f, Pdii f ) exposant de la loi de Hertz à la fréquence f et à la pression différentielle Pd .. De plus nous supposons que l'exposant de Hertz n(f,Pdff) ne Ppore f dépend que faiblement de la pression différentielle, soit: n(f,Pdiff ) = no(f)(1+8(f)xPdff) 8 une fonction de la fréquence, indépendante de la pression différentielle Pdiff et homogène à l'inverse d'une pression. y atténuation (yp et ys , respectivement de l'onde P et de l'onde S) y=rùtan _1(ù) Q facteur de qualité indépendant de la fréquence (modèle à Q-constant ), exposant de la loi puissance supposée, reliant l'atténuation y et la pression différentielle Pdiff . APde accroissement de pression différentielle entre deux états A et B B A APdiff = Pdiff Pdiff fl, f2 deux fréquences quelconques entre lesquels vont être exploitées les équations d'atténuation des vitesses. Par exemple si l'on s'intéresse au passage des vitesses soniques aux vitesses sismiques, la fréquence fi désigne la fréquence centrale de l'outil sonique (typiquement autour du kilohertz) et la fréquence f2 la fréquence centrale des enregistrements sismiques (typiquement autour de quelques hertz à quelques dizaines de hertz). no : coefficient de Hertz aux pressions différentielles tendant vers zero (Pd1 ff =0) La pression différentielle Pd ff dans un réservoir est la différence entre la pression de confinement Pconf, approximativement égale à la pression géostatique due au poids des sédiments sus-jacents, et la pression de pore PpOYe , ou pression régnant dans le fluide contenu dans le réservoir: Pd= Pconf ûP iff pore II est clairement établi que les propriétés sismiques, vitesse et atténuation, sont essentiellement fonctions de la pression différentielle Pd ff , et non des valeurs particulières (1) individuelles de la pression de confinement Pconf ou de la pression de pore Ppore (par exemple, Rasolofosaon et Zinszner, 2006). Considérons la loi puissance la plus simple reliant la vitesse P ou S et la pression différentielle Pdiff à la fréquence f , ou loi de Hertz-Mindlin (Rasolofosaon et Zinszner, 5 2004): A \n(f,Pdiif ) fiff Pdiff (2) où P d ff la pression différentielle à l'état X (=A ou B), et n( f, Pd ff) désigne l'exposant de la loi de Hertz à la fréquence f et à la pression différentielle Pdiff. De plus, nous supposons que l'exposant de Hertz n(f,Pdiff) ne dépend que faiblement de la pression 10 différentielle, soit: (3) n(f,Pdiff )=no(.f)(1+8(.f)xPdiff) où les dépendances en fréquence et en pression de n ont été découplées, et où l'hypothèse de faible dépendance de n vis-à-vis de Pdiff impose: 15 (4) I8(f)xPdiff « 8 une fonction de la fréquence, indépendante de la pression différentielle Pdiff et homogène à l'inverse d'une pression. Par ailleurs nous introduisons une loi reliant l'atténuation y = - tan-1(') [Q étant le facteur de qualité, modèle à Q-constant , indépendant de la fréquence (Bourbié et al., 1987) 20 et la pression différentielle. On suppose qu'il s'agit aussi d'une loi puissance caractérisée par un exposant y(PA "Pdiff \S A diff B B 7(Pdiff) Pdiff ~ La figure 1 illustre, sur. une échelle bi-logarithmique, des valeurs expérimentales de l'atténuation y de l'onde P (losanges noirs) et de l'onde S (triangles noirs), en fonction de la 25 pression différentielle (Pdie) mesurées au laboratoire. Les deux ajustements linéaires sont (5) également représentés en trait continu pour l'onde P et en trait discontinu pour l'onde S. Cas du grès de Fontainebleau. Par ailleurs si l'on admet que le modèle d'atténuation est à Q-constant (Bourbié et al., 1987) nous avons la relation d'atténuation suivante à la pression différentielle Pdff V (Pdi ) 11 Y (Pdiff ) ff ~(Pdiff . f~iJ f1) L'équation (3) reportée dans (2) conduit alors à: (7) A no(f)(1+5(f)xPdiff) Pdff B (Pdff)-no(f)8(fxd B iff Pdiff où APdiff =Pdff ùPâiff désigne l'accroissement de pression différentielle supposé 10 faible , c'est-à-dire: (8) APdi ff «Pâiff De plus le rapport membre à membre des deux équations d'atténuation (6) écrites aux pressions différentielles Pdff et Pdi ff se réduit à: 15 (9) V(Pâff,.f2) âff,.f2) V(Pdff,.fl) Les équations (7) et (9) entrainent: / A \no(f2)(1+5(f2)xPdff)-no(fl)1+S(fl)xPdiff) Pdiff (Pd Bff )-lno(f2)S(f2) no(fl)S(fl))x~dff B Pd1.f , 'A B = f2 )r(Pdiff1 p d e') (6) (PdzAff )no(f)(1+5(f)xPdff) (Pdff,.f) (B)no(f)(1+6(f)xPff J 20 (10) Soit en prenant le logarithme de chaque membre: /PA \ [no(f2)(1+S(.Î2)xPl)-no(.f1(1+8(.Î1)x dff)] ln B Pdiff / ù[no(f2)S(fz)ùno(fi)8(fi)]xAPdtff = [Y(Pd=f)ùr(Pf)] 1n / J En supposant l'accroissement de pression différentielle APdff =PP~ff ùPdff faible, PPZff et Pâ ff = Pd ff +APdI . ne sont pas trop différents l'un de l'autre, et par conséquent nous avons: ( A = ln Pdzff Paff B PdlA ff + OPd1i. de même, sous la même hypothèse, en posant y(Pd ff) = y(Pdlff)+Dy la différentielle logarithmique de l'équation (5) fournit: (13) Dy S APdiff ee' Y(Pdff) Pdtff Les équations (12) et (13) reportées dans l'équation (11) conduisent à: OP rr . ù[no(f2)(l+a(f2)xPff)ùno(fi)(1+«fi)xPdlfrf)] A Pdt ff ) OP. ù[no(f2)8(f2)ùno(fi)S(fi)]xAPdiff ùç 7( \ Aff Pdiff (12) ln Aide ù PdAtff Soit, en divisant chaque membre par A Pd`ff (15) n (.f2,Pdff)-n (fi,Pff ''-'' 4- ~ y(PdA1f) lnrffi/ [no(f2)S(f2)ùno(f1)S(f1)]xPff Cette équation est la relation recherchée, liant l'atténuation contenue dans le premier terme du membre de gauche en inf fil et l'effet de pression, couplé à une dépendance en fréquence, contenu dans le second terme du membre de gauche. Les deux membres de l'égalité sont fonctions de la fréquence et de la pression différentielle au premier ordre. Plus précisément, l'hypothèse de perturbation au premier ordre inclue d'une part l'hypothèse de faible dépendance de n vis-à-vis de Pdiff [condition (4)], et d'autre part, l'hypothèse de faible accroissement de pression différentielle APdff = Pdiff ùPdiff [condition (8)]. A l'ordre zéro, du fait de la condition (4), l'équation (15) se simplifie notablement en se réduisant à: (16) no(f2)ùno(fi)='r(Pd .) 1nl f/ En d'autre termes à l'ordre zéro, la dépendance de y en fonction de la pression différentielle peut être négligée. Il est toujours possible d'utiliser l'équation complète (15). 2c. Utilisation de la relation entre coefficient de Hertz et fréquence des ondes A partir de cette relation, on détermine de nouveaux coefficients de Hertz, et on définit donc une nouvelle loi de dépendance vitesses sismiques / pression de pore. Cette loi est très bien vérifiée expérimentalement comme l'illustre la figure 2. Cette figure illustre les vitesses expérimentales ultrasonores (courbes repérées US pour 'ultrasons') et soniques (courbes repérées SO pour 'sonic') de l'onde P (carré et triangle) et de l'onde S (losange et rond) en fonction de la pression différentielle (Pdiff), en échelle bi-logarithmique. Les prédictions théoriques des vitesses P et S soniques (en traits discontinus) sont comparées aux vitesses soniques expérimentales. Les ajustements linéaires de ces données expérimentales sont excellents, validant une nouvelle fois la loi de Hertz. De plus les prédictions théoriques des vitesses P et S soniques (en traits discontinus) résumées par la formule (16) sont tout-à-fait comparables aux vitesses soniques expérimentales. Ainsi, d'une part le coefficient de Hertz de la loi puissance entre la vitesse et la pression différentielle (équation 2) n'est pas la même dans les différentes bandes de fréquences et, d'autre part cela est vérifié expérimentalement et prédit par la formule simple (16). Nous disposons alors de nouvelles courbes de dépendances vitesses sismiques pression de pore aux fréquences sismiques. . Interprétations des imaqes sismiques L'interprétation des images sismiques, issues de la sismique 4D et de l'inversion, a pour objectif de déterminer si la variation de vitesses sismiques (ou impédances) observées sur les différentes sismiques acquises à différents temps calendaires, est principalement due à une variation de pression (fluide et confinement) ou à une substitution de fluide. Cette évolution des fluides produits ou injectés au sein du milieu poreux se caractérise, à partir de deux jeux de données sismiques, par deux paramètres : la variation du module d'incompressibilité du fluide saturant, et la variation de la pression de pore. Le premier paramètre indique le type de fluide en un point donné du milieu, le second indique l'augmentation de pression induite. Ces deux paramètres sont généralement déterminés à partir des données sismiques inversées (vitesses ou impédances) : la loi de Biot-Gassmann relie, entre autres, le module d'incompressibilité du fluide saturant aux impédances ou vitesses sismiques, et la loi de Hertz-Mindlin relie la pression de pore aux impédances ou vitesses sismiques. Selon l'invention, on utilise la loi de Hertz-Mindlin avec un coefficient de Hertz modifié selon l'étape 2. La figure 3 illustre les résultats d'interprétation selon la méthode. Sur cette figure sont représentées les variations des vitesses V en fonction de la pression de pore Pp, dans une roche réservoir clastique typique (grès), saturée d'eau (W), d'huile (0) ou de gaz (G), et pour trois types d'exploitation : le balayage à l'eau (repères A), la production de gaz (repères B) et l'injection de gaz (repères C) telle l'injection de CO2 dans un aquifère salin. Les courbes des vitesses sismiques selon l'invention sont indiquées par des courbes en traits continus et les vitesses ultrasonores "classiques" par des courbes en traits discontinus. La vitesse sismique augmente avec la pression différentielle Pdz ff comme l'illustre la figure 2 et l'équation (2). Comme la pression de confinement PConf, égale à la pression géostatique due au poids des sédiments sus-jacents, reste approximativement constante lors de l'exploitation d'un réservoir, seule la pression de pore Ppo,.e varie substantiellement. Aussi, à une augmentation de la pression de pore correspond une diminution de la pression différentielle et donc, une diminution de la vitesse ou de l'impédance sismique, ce qui est observé sur la figure 3. D'après l'équation (16) cette effet est beaucoup plus marqué aux basses fréquences, par exemple aux fréquences sismiques (typiquement de quelques hertz 14 à quelques dizaines de hertz), qu'aux hautes fréquences, par exemple aux fréquences ultrasonores du laboratoire (typiquement de quelques centaines de kilohertz au mégahertz). Considérons ainsi une injection de gaz (repères C), telle l'injection de CO2 dans un aquifère salin. Si l'on se base sur des données ultrasonores au laboratoire, lors d'une injection de gaz (repère C) dans un réservoir la pression de pore augmente (état C2 à état C3) et induit une diminution d'impédance notablement moins importante que celle induite par la substitution de l'eau par du gaz (état C, à état C2) lors du même mode d'exploitation. En revanche, toujours lors du même mode d'exploitation, dans la bande de fréquence sismique (courbes discontinues) l'augmentation de pression de pore (état C', à état C'2) peut induire une diminution d'impédance tout-à-fait comparable à celle induite par la substitution de l'eau par du gaz (état C'2 à état C'3). Ainsi une interprétation des données de sismique 4D basées uniquement sur la loi vitesses-pression ultrasonore obtenue au laboratoire conclurait, de manière erronée, à une influence majeure des substitutions de fluides par rapport aux variations de pression de pore. Le phénomène est encore plus marqué lors de la production de gaz (repère B). Si l'on se base sur des données ultrasonores au laboratoire, lors d'une production de gaz dans un réservoir la pression de pore diminue (état B, à état B2) et induit par conséquent une augmentation d'impédance notablement moins importante que celle induite par la substitution du gaz de l'huile (état B2 à état B3) lors du même mode d'exploitation. En revanche, toujours lors du même mode d'exploitation, dans la bande de fréquence sismique (courbes discontinues) la diminution de pression de pore (état B', à état B'2) peut induire une variation d'impédance encore plus importante que lors de l'injection de gaz (état C'2 à état C'3) et tout-à-fait comparable à celle induite par la substitution de gaz par de l'huile (état B'2 à état B'3). Ainsi une interprétation des données de sismique 4D basées uniquement sur la loi vitesses-pression ultrasonore obtenue au laboratoire conclurait une nouvelle fois, de manière erronée, à une influence majeure des substitutions de fluides par rapport aux variations de pression de pore. Dans le cas de récupération assistée par balayage à l'eau (repère A) le phénomène est beaucoup moins important. En effet l'influence de la pression de pore mesurée par ultrasons (état A2 à état A3) semble tout-à-fait comparable à celle mesurée aux fréquences sismiques (état A'2 à état A'3). Toutefois l'influence de la substitution fluide semble un peu plus marquée aux fréquences sismiques (état A', à état A'2) qu'aux fréquences ultrasonores (état A, à état A2). L'invention permet donc de contribuer à l'interprétation correcte des sismiques 4D dans des opérations d'exploitation du sous-sol, telles l'injection de gaz, la production de gaz ou le balayage à l'eau (figure 1), que ce soit dans le domaine environnemental (stockage géologique de CO2) ou dans le domaine de l'industrie pétrolière ou gazière. L'invention trouve des applications dans le domaine de l'exploitation des réservoirs souterrains ('monitoring' sismique), dans le domaine du forage (problème sécurité-évaluation des pressions de fluides dans le sous-sol), et dans l'évaluation des systèmes pétroliers (modélisation des migrations de fluides). En effet, ce brevet permet de faire plus correctement la part entre les variations de pressions (fluide et confinement) et des substitutions fluides dans les modifications observées sur les différentes sismiques observées à différents temps calendaires, et par conséquent de fournir un modèle de réservoir plus réaliste loin des puits et un meilleur calage de production ("history matching"), et donc une meilleure prévision de l'exploitation du réservoir. Les principaux "bénéfices" attendus sont pour le management de réservoir. En effet, la méthode, permettant de mieux décrire comment les répartitions de fluides dans le réservoir varient avec le temps. Ces "bénéfices" peuvent se regrouper autour des problématiques suivantes: Compartimenter le réservoir: la méthode permet de repérer les différents compartiments du réservoir et, par conséquent, d'estimer les volumes d'huile et de gaz non récupérés (dans le domaine de l'industrie pétrolière), ou de repérer et estimer les volumes de CO2 effectivement stockés dans le sous-sol et contrôler l'intégrité des couvertures (dans le domaine environnemental). Connectivité des différents compartiments: la méthode, permettant de repérer les différents compartiments du réservoir et leur balayage par les fluides, permet d'aider à prendre les meilleures décisions pour l'exploitation du réservoir (emplacement des nouveaux puits par exemple dans le domaine pétrolier, intégrité des couvertures pour le stockage de CO2: fuites de CO2 ou pas...). - Répartition des contacts fluides: la méthode, permettant de mieux suivre l'évolution spatiale des contacts fluides lors de la production, permet de confirmer ou, au contraire d'invalider le modèle d'écoulement, et par conséquent permet au producteur pétrolier de modifier son programme de récupération. Dans le domaine de l'environnement, cela permet de suivre plus distinctement le front de CO2 stocké dans le sous-sol.35 REVENDICATIONS1 Méthode pour exploiter un milieu poreux souterrain par injection ou production d'un fluide, dans laquelle on acquiert au moins deux jeux de données sismiques en réalisant au moins deux campagnes sismiques réalisées à des moments différents par émission d'ondes sismiques dans ledit milieu, on détermine deux jeux de vitesses desdites ondes sismiques en inversant lesdits jeux de données, et on considère une loi de puissance de type Hertz-Mindlin reliant au moyen d'une puissance appelée coefficient de Hertz, une vitesse de propagation d'ondes à une pression de pore, caractérisé en ce que l'on réalise les étapes suivantes: - on détermine ledit coefficient de Hertz au moyen de mesures sur des échantillons dudit milieu en utilisant des ondes de fréquences f1; on corrige ledit coefficient de Hertz pour tenir compte d'une différence de fréquence entre lesdites ondes de fréquence f1 et lesdites ondes sismiques, au moyen d'une mesure d'atténuation des ondes de fréquence f1 dans lesdits échantillons; - on détermine en tout point du milieu poreux une variation de pression de pore au moyen de ladite loi utilisant ledit coefficient de Hertz corrigé et desdits jeux de données inversés, et - on suit l'exploitation dudit milieu poreux souterrain par injection ou production d'un fluide, en déduisant desdites variations de pression de pore une localisation dudit fluide et éventuellement une substitution dudit fluide au sein du milieu poreux. 2. Méthode selon la revendication 1, dans laquelle on corrige ledit coefficient de Hertz au moyen de la formule suivante : n (.f2,Pdff)ùn (fi,P ff)'(Pdlf) In(fYi) + [no(f2)S(f2)ùno(fi)S(fi)]xPdf avec : f2 fréquence centrale des enregistrements sismiques; Pdff pression différentielle à l'état X (=A ou B); n(f,Pd ff) : coefficient de Hertz à la fréquence f et à la pression différentielle Pdtff n (f2,Pâf.) étant le coefficient de Hertz corrigé et n (f,,PdA,ff ) le coefficient de Hertz non corrigé; y : atténuation des ondes P et S aux fréquences f1; : un réel; no : coefficient de Hertz aux pressions différentielles tendant vers zéro, no(f2) étant le coefficient de Hertz corrigé et no(f) le coefficient de Hertz non corrigé aux pressions différentielles tendant vers zéro. 3. Méthode selon la revendication 1, dans laquelle on corrige ledit coefficient de Hertz au moyen de la formule suivante : no (f2 ) ù no (f) C7(1 )ln(% ) avec : f2 : fréquence centrale des enregistrements sismiques; no : coefficient de Hertz aux pressions différentielles tendant vers zéro, no(f2) étant le coefficient de Hertz corrigé et no(fi) le coefficient de Hertz non corrigé; Pdiff : pression différentielle dans ledit milieu; y : atténuation des ondes P et S aux fréquences fi; : un réel. 4. Méthode selon l'une des revendications 1 à 3, dans laquelle les fréquences f1 sont des fréquences ultrasonores. 5. Méthode selon l'une des revendications 1 à 3, dans laquelle les fréquences f1 sont des fréquences soniques.