FR 2932574 A1 20091218 FR 0803385 A 20080617 10 La présente invention concerne le domaine de la caractérisation de gisement souterrain. En particulier l'invention concerne une méthode pour évaluer de façon quantitative les pressions de fluide dans le sous-sol. On peut utiliser la méthode dans le domaine pétrolier pour la détection des zones de surpression pour le forage. 15 La présence de zones de surpression en exploration peut avoir de lourdes conséquences financières, et parfois humaines, pour le forage en cas de méconnaissance des pressions de fluide. Prédire la présence de zones en surpression et plus généralement arriver à évaluer quantitativement les surpressions est devenu 20 une priorité pour les compagnies pétrolières. En effet, dans le domaine de l'exploration, la pression de fluide peut approcher la contrainte minimale principale et induire la réouverture de fracture ou éventuellement initier une fracturation hydraulique. Dans le domaine du forage, il est important de connaître l'écart de pression entre la pression de fluide et la contrainte minimale en place pour le design 25 des � casings � de puits et pour prévoir le poids de boue de manière à prévenir les éruptions (blow-out) en forage en � underbalanced � ou les pertes de boue en forage en � overbalanced �. Enfin la déplétion dans les zones de surpression peut induire des réagencements notables des contraintes avec éventuellement des conséquences importantes sur la productivité des réservoirs. Aussi une bonne évaluation quantitative des pressions de fluide et de leurs liens avec les variations de contrainte est également importante dans le domaine de la production. Présentation de l'art antérieur Il existe de nombreuses méthodes permettant à partir de mesures physiques (et non à partir de modélisations) une évaluation quantitative des pressions de fluide : tests des formations, vitesse d'avancement de l'outil de forage, mesures de densité des argiles, Indices de gaz, débit-métrie, diagraphies différées, etc. Parmi ces méthodes, les méthodes géophysiques, et plus particulièrement les méthodes sismiques, à plus grande résolution spatiale que ses concurrentes (gravimétrie par exemple), sont les seules à pouvoir y arriver loin des puits. Il est donc primordial d'exploiter au mieux les données sismiques. Toutefois les traitements sismiques conventionnellement utilisées, type analyse de vitesse, ont une efficacité limitée, principalement du fait d'une part de leur résolution spatiale trop faible pour être efficacement utilisées pour le forage, et d'autre part de la faible prise en compte des variations lithologiques, souvent critiques dans les phénomènes de surpression : Reynolds, E.B, 1970, Predicting overpressured zones with seismic data: World oil, 171, 78-82. Pour comprendre ces principaux problèmes techniques, on décrit sommairement la procédure classique d'évaluation quantitative des surpressions à partir des données sismiques. Les différentes étapes sont les suivantes: - obtention d'un modèle de vitesses sismiques aussi précis que possible par analyse fine de vitesse ; - déduction d'une courbe de compaction de référence (vitesse sismique en fonction de la profondeur) dite de "compaction normale" (correspondant à la répartition hydrostatique de la pression fluide) - interprétation des écarts entre la courbe de compaction mesurée par la sismique et la courbe de compaction normale en terme d'anomalies de pression fluide. Les anomalies (ou écarts par rapport à l'hydrostatique) peuvent être positives (surpressions) ou négatives. Les principaux problèmes de ces méthodes sont d'une part la faiblesse relative de la résolution spatiale des méthodes conventionnelles, les rendant difficilement utilisable pour les opérations de forage. Le second problème est l'hypothèse implicite d'attribuer tout changement anomalique de vitesse à une surpression, en excluant par exemple des causes telles le changement de lithologie (Reynolds, 1970, par exemple). La vérification lithologique vient a posteriori dans la méthode classique. En d'autres termes on vérifie après tous les traitements que les anomalies de pressions ne sont pas dues à une variation lithologique. On connaît également, par le brevet FR 2.893.421, une méthode basée sur l'inversion avant sommation des données sismiques. Cette méthode est difficile à mettre en oeuvre, en particulier l'analyse lithosismique des cubes d'impédances sismiques. Cette analyse lithosismique, particulièrement délicate et longue à mettre en oeuvre, est une interprétation des cubes 3D d'impédances sismiques en terme de faciès sismique au regard des différentes lithologies rencontrées aux puits. L'objet de l'invention est une méthode alternative pour évaluer des pressions de fluide dans une zone du sous-sol à partir de données de puits et de données sismiques. La méthode permet de surmonter les difficultés de l'art antérieur en fournissant un cube des pressions de fluide à une échelle suffisamment précise pour le foreur tout en prenant en compte la lithologie explicitement dans le traitement. La lithologie est alors réduite à une différenciation entre lithologie argileuse et lithologie non argileuse. La méthode selon l'invention L'invention concerne une méthode pour évaluer des pressions de fluide dans une zone du sous-sol à partir de données de puits et de données sismiques. Elle comporte les étapes suivantes : on construit un cube de vitesse d'onde sismique P et un cube de vitesse d'onde sismique S au moyen d'une inversion stratigraphique desdites données sismiques, et on déduit desdites vitesses un cube de lithologies identifiant des lithologies argileuses et des lithologies non argileuses ; - on détermine, à partir des données de puits et pour chacune des deux lithologies, une relation permettant d'estimer la pression de fluide à partir des vitesses d'onde sismique P ; on évalue les pressions de fluide dans ladite zone du sous-sol en construisant un cube de pressions de fluide en appliquant lesdites relations audit cube de vitesse d'onde sismique P en fonction dudit cube de lithologies. Selon un mode de réalisation, les données sismiques comportent au moins un cube sismique discrétisant la zone en volumes élémentaires repérés par leurs coordonnées horizontales (x, y) et verticale en temps (t) ; on construit alors le cube de lithologies en déterminant un cube de vitesse d'onde sismique S au moyen de l'inversion stratigraphique, et en appliquant une valeur seuil de rapport des vitesses d'onde sismique P et S, de façon à ce que l'on affecte une lithologie argileuse aux volumes élémentaires ayant un rapport supérieur à ce seuil, et une lithologie non argileuse aux autres volumes. Le seuil peut par exemple être égal à 2. Les relations peuvent être de la forme suivante : Pe (z) = P" (z . VP (z) pore pore ) VP (z) 4 avec : VP (z) : vitesse des ondes sismiques P mesurée aux puits à une profondeur z VP (z) : vitesse des ondes sismiques P estimée dans une hypothèse où il n'existe aucune surpression à la profondeur z Pore (z) : pression de fluide mesurée aux puits à une profondeur z PPo,e (z) : pression de fluide estimée dans une hypothèse où il n'existe aucune surpression à la profondeur z Dans les lithologies argileuses, on peut déterminer VP (z) au moyen des étapes suivantes : - on identifie un intervalle en profondeur où la pression de fluide est voisine de la pression hydrostatique, et on définit ladite relation Vp (z) par une relation linéaire sur cet intervalle. Dans les lithologies non argileuses, on peut déterminer au moyen des étapes suivantes : on mesure la pression de fluide au moyen de diagraphie en un nombre réduit de profondeurs z, et pour les mêmes profondeurs z, on calcule une pression de fluide égale à p.g.z , avec p 1030kg / m3 et g 9,81m / s2 ; on mesure la vitesse sismique des ondes P à ces mêmes profondeurs ; on déduit V; (z) en interpolant une droite entre des profondeurs où la pression de pore a pu être mesurée. Selon le méthode, on peut déterminer des zones de surpressions au sein de la zone, en construisant un cube de pression de confinement Pconf(x, y, t) et en détectant les zones de suppression lorsque la pression de fluide est supérieur à a*Pconf(x, y, t), où a est un seuil préalablement choisi. Ce seuil peut être égal à 0,9 par exemple. 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 sommaire des figures - La figure 1 fournit une illustration générale des différentes étapes de la méthode. - La figure 2 représente en détail la partie de la figure 1 correspondant à la boucle d'analyse sismique (SAL). - La figure 3 illustre le rapport des vitesses sismiques Vp et Vs en fonction du 10 coefficient de Poisson (y) pour différentes lithologies (courbe) ; les lithologies argileuses (Arg) se distinguent des lithologies non argileuses (Narg) par un rapport Vp/Vs très élevé, typiquement supérieur à 2 (droite horizontale). - Les figures 4 et 5 représentent des variations de la vitesse (log Vp) des ondes sismiques P en fonction de la profondeur (z) dans un certain nombre de puits de 15 calibration. La figure 4 illustre le cas des lithologies argileuses, et la figure 5 celui des lithologies non argileuses. Description détaillée de la méthode La méthode permet d'évaluer les pressions de fluide d'une zone du sous-sol à 20 partir de données de puits (telles que des diagraphies) et de données sismiques. Elle comporte principalement les trois étapes suivantes : 1) Construction de cubes de vitesses sismiques et d'un cube de lithologies argileuses, par traitement des données sismiques ; 2) Détermination d'une relation entre la vitesse des ondes sismiques P et la 25 pression de fluide pour chacune des deux lithologies ; 3) Construction d'un cube de pression de fluide (PpOYe(x, y, t)). Les différentes étapes de la méthode sont représentées schématiquement sur les figures 1 et 2. La figure 1 fournit une illustration générale des différentes étapes de la méthode. Ces étapes se composent de deux sous ensembles d'étapes, à savoir, sur le côté gauche de la figure 1, les étapes correspondant au traitement des données de puits, et sur le côté droit de la figure 1, les étapes correspondant au traitement des données sismiques (détaillées sur la figure 2). Suivant la convention habituelle, les rectangles contiennent les données d'entrées, ou les résultats obtenus à une certaine étape du traitement, l'étape étant repérée par un chiffre suivi éventuellement d'une lettre romaine minuscule (la par exemple). Ces rectangles sont reliés par des flèches descendantes pour la plupart, orientées dans le sens séquentiel des traitements, depuis l'entrée des données (en haut de la figure), jusqu'à la sortie finale des résultats (en bas de la figure). Pour clarifier la description, des rectangles sont parfois accompagnés du descriptif très sommaire de la technique (SI, ...) permettant de passer d'un résultat au résultat suivant. Selon la méthode, les pressions de fluide sont évaluées dans une zone du sous-sol sous la forme d'un cube discrétisant la zone à étudier. Cette discrétisation consiste à découper la zone en volumes élémentaires repérés par leurs coordonnées horizontales (x, y) et verticale soit en temps (t), soit en profondeur (z). 1) Construction d'un cube de vitesse sismique et d'un cube de lithologies argileuses Il est bien connu que les milieux argileux ou argilo-sableux, typiquement de perméabilités faibles, sont plus favorables au développement des surpressions. Par contre les milieux sableux, souvent de lithologie grossière et surtout suffisamment perméable, facilitent l'écoulement des fluides contenus dans leurs pores, ce qui empêche le développement des surpressions. Ainsi toute méthode ne prenant pas en compte d'emblée ces réalités géologiques sont biaisées à la base et ne peuvent conduire qu'à des résultats erronés au moins partiellement avant correction a posteriori en tenant compte d'une manière ou d'une autre les variations de lithologie. La prise en compte de cette lithologie argileuse peut se faire (voir figure 3) en utilisant le lien existant entre la lithologie et le rapport des impédances sismiques P (Ip) et S (Is) , ou de manière équivalente le rapport des vitesses des ondes sismiques P (Vp) et S (Vs), car : IP VP Is Vs En effet, les lithologies argileuses se distinguent par des rapports Vp/Vs très élevés, typiquement entre 1,9 et 3, contrairement aux autres lithologies des bassins sédimentaires, à savoir principalement les sables/grès avec 1,6 Classiquement, on découpe la zone en intervalles d'analyse en temps. A partir des données sismiques, on identifie des horizons, encore appelés � marqueurs sismiques �. Ces horizons indiquent des discontinuités sismiques, lithologiques ou non, caractérisées par une variation de l'impédance sismique. On considère donc généralement que la partie du sous-sol comprise entre deux horizons est homogène d'un point de vue de ses propriétés pétro-élastiques. Ainsi, on divise la zone souterraine en plusieurs intervalles d'analyse en temps, délimités par des horizons sismiques, dans le but d'obtenir une précision accrue dans les résultats. Chaque intervalle d'analyse en temps est ainsi traité séparément pour identifier des propriétés bien spécifiques (ondelette, relation entre la lithologie et la sismique, etc.) et successivement, pour fournir un résultat global, décrit ultérieurement. Généralement on choisit ces intervalles d'analyse en temps inférieurs à 500ms, typiquement de l'ordre de 300ms à 400ms. On débute ensuite l'analyse avec un premier intervalle d'analyse en temps, TAI (étape 2a), et l'on tronque les cubes 3D correspondant à chacune des classes d'angles choisies pour se réduire à ce premier intervalle d'analyse en temps (étape 3a). Puis, à partir de ces cubes tronqués (TSDA1(x, y, t), TSDA2(x, y, t),...), on effectue une inversion stratigraphique (SI) avant sommation avec utilisation d'information a priori géologique. Cette technique est bien connue des spécialistes et l'on peut par exemple utiliser les techniques proposées par : - Brac J.P. et al., 1988, Inversion with A Priori Information: An Approach to Integrated Stratigraphic Interpretation, Reservoir Geophysics R.E. Sheriff ed. Investigation in Geophysics, 7, SEG, Tulsa. - T. Tonellot, D. Macé, V. Richard, 1999, Prestack elastic waveform inversion using a priori information, 69th Ann. Internat. Mtg: Soc. of Expl. Geophys., paper 0231, p.800-804. - Lucet, N., Déquirez, P. -Y. and Cailly, F., 2000, Weil to seismic calibration: A multiwell analysis to extract one single wavelet, 70th Ann. Internat. Mtg: Soc. of Expl. Geophys., 1615-1618. Ce type d'inversion comprend deux phases. La première phase (WE), selon la méthode décrite par Lucet et al. (2000), consiste à extraire pour chaque cube tronqué, c'est-à-dire pour chaque classe d'angle, la meilleure ondelette (wl(t), w2(t), ...) cohérente avec les données observées au puits (étape 4a). La seconde phase (MB), décrite par Tonellot et al. (1999), consiste à construire un modèle 3D a priori (étape 4b) nécessaire pour initier et contraindre l'inversion dans la prochaine étape. Il s'agit principalement de deux cubes 3D d'impédances sismiques, à savoir le cube a priori d'impédance des ondes P, noté Ip,,n(x, y, t), et le cube a priori d'impédance des ondes S, noté Is,,n(x, y, t). Les coordonnées x et y sont les deux coordonnées horizontales liées à l'acquisition, typiquement � on-line � et � cross-line �. La troisième dimension n'est pas la profondeur z mais le temps d'enregistrement t, directement lié à la mesure sismique. Enfin l'inversion (SI) à proprement parler est effectuée. Plus précisément, la connaissance des ondelettes et du modèle a priori sur l'intervalle d'analyse en temps choisi, permet d'inverser simultanément tous les cubes 3D (TSDA1(x, y, t), TSDA2(x, y, t),...), via une inversion stratigraphique avant sommation selon la méthode décrite par : - Tonellot, T., Macé, D. and Richard, V., 2001, Joint stratigraphie inversion of angle-limited stacks, 71 st Ann. Internat. Mtg: Soc. of Expl. Geophys., 227-230. Cette inversion produit deux cubes 3D d'impédances sismiques, à savoir le TAI cube d'impédances sismiques des ondes P, noté I P (x, y't) , et le cube d'impédances 'TAI sismiques des ondes S, noté s (x, y, t) , ainsi qu'un cube de densité pTA' (x, y, t) (étape 5a). A l'issu de cette inversion, l'ensemble des cubes IP I (x, y, t) , I T 2 (x, y, t) ., ITA` (x, y, t) , ITA2 (x, y, t) , etc., permettent de former deux cubes d'impédances et un cube de densité représentatifs de l'ensemble de la zone étudiée (étape 7a) : - IP(x,y,t) - Is(x,y,t) p(x, y, t) En divisant les cubes d'impédance par le cube de densité, on obtient deux 10 cubes de vitesse sismique, à savoir le cube de vitesses sismiques des ondes P, noté Vp (x, y, t) , et le cube de vitesses sismiques des ondes S, noté Vs (x, y, t) , Par ailleurs, à partir des deux cubes d'impédances, on construit un troisième cube Vp, s (x, y, t) , correspondant, en chaque point de discrétisation, au rapport des deux impédances, ou de manière équivalente au rapport des deux vitesses sismiques : V,,s(x,y,t) V(x,Y,t) Ip(x,Y,t) Vs(x,y,t) /s(x,Y,t) 15 Enfin, à partir de ce cube Vp, s (x, y, t) , on construit (étape 7a) un cube de lithologie argileuse Arg(x, y, t) : ce cube binaire indique la localisation des faciès argileux et celle des faciès non argileux. Pour ce faire on choisit un seuil de rapport Vp/Vs. Selon un exemple, sans que 20 ce soit limitatif dans la méthode parce parfaitement modulable, on utilise la valeur seuil 2 comme valeur limite inférieure pour le rapport Vp/Vs dans les lithologies argileuses. Par conséquent, par convention, tous les milieux caractérisés par un rapport Vp/Vs inférieur à 2, correspondent à des lithologies non argileuses. A l'issu de cette boucle sismique (SAL), on dispose des deux cubes suivants : Vp(x,y,t) et Arg(x,y,t). 2) Détermination d'une relation entre la vitesse des ondes sismiques P mesurée et la pression de fluide, pour chacune des deux lithologies (étape 7b) Ces relations sont établies par traitement des données de puits en tenant compte des deux types de lithologie. Ceci tient au fait que dans les lithologies non argileuses les mesures de pression (par exemple par diagraphie MDT) sont possibles, alors que les pressions de fluides ne peuvent être qu'estimées dans les lithologies argileuses. Selon un exemple, les données de puits, WD(z), acquises au cours d'une étape lb, comprennent principalement : les vitesses sismiques issues des diagraphies acoustiques : vitesse des ondes P notée Vp(z), et éventuellement vitesse des ondes S notée Vs(z) ; la pression de fluide, notée Ppore(z) ; - le type de lithologie argile ou non argile, noté Arg(z), où z désigne la profondeur. Le détail de cette nouvelle méthode de calibration est illustré par les figures 4 et 5 pour les deux types de lithologie. Plus précisément, sur les figures 4 et 5 sont portées les variations de la vitesse (log Vp) des ondes sismiques P en fonction de la profondeur (z) dans un certain nombre de puits de calibration. Selon un exemple de réalisation, sans que ce soit limitatif dans la méthode parce que parfaitement interchangeable avec une autre relation reliant la vitesse sismique et la pression de pore, on utilise la relation suivante : pore Y P(Z) Pe (z) = Pn (z) V . (z) pore (1) avec : Vp �z) : vitesse des ondes sismiques P mesurée aux puits à une profondeur z Vp "(z) vitesse des ondes sismiques P estimée dans une hypothèse où il n'existe aucune surpression à la profondeur z P4 ore (z) : pression de fluide (pore) mesurée aux puits à une profondeur z P" (z) : pression de fluide (pore) estimée dans une hypothèse où il n'existe pore aucune surpression à la profondeur z La pression de fluide "normale" Pp°re (z) à la profondeur z considérée, est donnée par pwgZ , où pw 1030 kg / m3 g 9,81m / s 2 et Z désignent respectivement la masse volumique de l'eau de mer, l'accélération de la pesanteur et la profondeur. La vitesse des ondes sismiques P mesurée aux puits à une profondeur z VP (z) est connue à l'issue de l'inversion stratigraphique (cube VP (x, y, t) ) Ainsi, on peut utiliser la relation suivante, dans laquelle il reste à déterminer Vp" (z) pour chaque lithologie : Vm(z) P;°re (z) = p.g.z. Vp �z) Estimation de Vp" (z) On estime la vitesse des ondes sismiques P dans l'hypothèse où il n'existe aucune surpression à la profondeur z. La technique consiste à définir une relation entre Vp et z par traitement des données de puits. A partir de ces données de puits, on trace pour chacune des deux lithologies (argileuse ou non argileuse), une courbe représentant le logarithme de la vitesse des ondes P en fonction de la (2) profondeur. Pour les lithologies argileuses, un exemple est illustré sur la figure 5, et pour les lithologies non argileuses, un exemple est illustré sur la figure 5. Ensuite, pour les lithologies argileuses (figure 4), on identifie un intervalle en profondeur où la pression de fluide est voisine de la pression hydrostatique. Cet intervalle (NT) caractérise un comportement dit �normal �, c'est-à-dire, sans surpression dans le sous-sol. Cette information est par exemple fournie par le foreur, par le géologue, ou constitue une information régionale connue par ailleurs. Sur l'exemple de la figure 4 cet intervalle va de la surface jusqu'à une profondeur inférieure à 2500m. Puis, on définit une relation linéaire log(Vp)=f(z) sur cet intervalle. Cette relation définit une tendance normale dans les argiles, notée VP,arg)(z), c'est-à-dire qu'elle traduit l'évolution de la vitesse des ondes P dans des argiles du sous-sol en absence de pression de pore anormalement élevées. Ensuite, pour les lithologies non-argileuses (figure 5), la technique de détermination d'une tendance normale, notée V' g) (z) , comprend trois étapes : on mesure la pression de fluide, au moyen de diagraphie MDT par exemple, en un nombre réduit de profondeurs z, et pour les mêmes profondeurs z, on calcule la pression de fluide dite � normale �, c'est-à-dire p.g.z . on mesure la vitesse sismique des ondes P à ces mêmes profondeurs ; - on déduit VP�narg)(z) en utilisant la relation (1) pour calculer la vitesse dite � normale � aux mêmes profondeurs, puis en interpolant une droite entre les points où la pression de pore a pu être mesurée. Ainsi, dans les argiles on utilise : V(z) Pore (z) = p.g.z. V n arg) (z) P et dans les lithologies non argileuse, on utilise :25 Vm(z) P;ore(z) = p.g.z. Vn(narg) (z) P 3) Construction d'un cube de pressions de fluide PpOYe(x, y, t) (étape 8) Pour construire le cube de pressions de fluide PpOYe(x, y, t) on scanne le cube de lithologie argileuse Arg(x, y, t) , ainsi que le cube de vitesse VP (x, y, t) , pour obtenir en chaque point x, y, t une valeur de lithologie et de vitesse P. En fonction de la lithologie, on applique la relation définie à l'étape 2 pour attribuer une valeur de pression de pore au point x, y, t. Si les relations on été établies en profondeur z, il est nécessaire de convertir les cubes Arg(x, y, t) et Vp (x, y, t) en profondeur. Une telle conversion temps ù profondeur est une technique classique pour l'homme du métier. On peut également convertir les données de puits, pour obtenir des relations directement en temps : Ainsi, dans les argiles on utilise : V m t Ppore(t) = Ppore(t). Vn(ar) (t) P et dans les lithologies non argileuse, on utilise : V1(t) Ppore (t) = Pp re (t)' 7p n(n arg) (t) P On peut ainsi construire un cube de pressions de fluide en temps Ppore(x, Y, t), ou en profondeur PpOYe(x, y, z) : Pour un point de coordonnées x, y et z : - si Arg(x,y,z) indique une lithologie argileuse, alors : Vp (x, y, z) Ppore (x, y, z) = p.g.z. Vn(arg) (z) - si Arg(x,y,z) indique une lithologie non argileuse, alors : V. (x,Y,z) Ppore(x,Y,Z) = p.g.z. Vn(narg)(z) P 4) Détermination de zones de surpressions En temps A partir de la méthode selon l'invention, il est possible de prédire, avec une bonne résolution spatiale, les zones éventuelles de surpressions pouvant présenter un danger lors d'un forage pétrolier par exemple. En effet, un cube de pression de confinement Pconj(x, y, t) peut être obtenue par l'équation suivante : Pconf (x,Y,t) = 2 p(x,y,t) VP(x,Y,t) g dt Pconf (x, y, t) = 1 IP (x, y, t)g dt Ainsi les cubes 3D de pression de fluide PpOYe(x, y, t) et de pression de confinement Pconf(x, y, t), étant connus, il suffit d'appliquer un critère de seuil sur la pression de pore au choix de l'utilisateur. On peut par exemple, sans que ce soit limitatif dans la méthode parce que parfaitement modulable, utiliser le seuil égale à 0,9*Pconf(x, y, t). En profondeur Selon la méthode, il est également possible d'obtenir le même type de résultat en profondeur z et non en temps t, ce qui peut être capital pour définir les conditions de forage. En effet, le problème du passage des cubes en coordonnées temporelles à des cubes en coordonnées exprimées en profondeur est un problème général bien connu en traitement sismique, et toute méthode ayant fait les preuves de son efficacité (stretch vertical, �map migration �, etc.) peuvent être appliquées. En effet une telle méthode de conversion permet simplement de passer de � l'espace temps � à � l'espace profondeur � en utilisant la fonction de conversion de la variable temps t en la variable profondeur z. On obtient alors sans difficulté un cube en profondeur des pressions de fluide PpOYe(x, y, z). Cette évaluation quantitative des pressions de fluide dans le sous-sol permet évidemment de repérer, en profondeur cette fois, les zones de surpressions anormalement élevées et pouvant présenter un danger pour un forage pétrolier. A partir de cette information, on peut modifier la trajectoire du puits de forage pour éviter ces zones de surpression, ou encore, on peut modifier la pression d'injection des fluides de forage pour compenser les surpressions des fluides du sous-sol. Selon un autre mode de réalisation, les surpressions sont détectées directement sans construire un cube de pression de pore. En effet, les tendances normales, fournissent pour chaque lithologie une valeur de vitesse d'onde sismique P en fonction de la profondeur. Tout écart entre la valeur de vitesse Vp, issue de l'inversion stratigraphique et la vitesse Vp donnée par cette tendance normale est interprétée comme une surpression. Avantages La méthode selon l'invention permet donc d'estimer les pressions de fluide dans une zone du sous-sol, ainsi que les zones de surpression, en temps ou en profondeur, et même pour des profondeurs non encore atteintes. La méthode est caractérisée par une grande résolution spatiale par rapport aux méthodes conventionnelles basées sur des analyses de vitesse. Elle permet de définir les conditions de forage (trajectoire, pression de fluide de forage, etc.), car elle donne des résultats très précis en tenant compte de la lithologie dès le début du traitement de manière quantitative, et non de manière qualitative et a posteriori comme c'est le cas dans les approches conventionnelles. Enfin, la méthode exploite au maximum des données acquises au niveau de l'échelle sismique, proche de la mesure directe, pour éviter les problèmes de changement d'échelle (géologique, réservoir et, sismique). Il faut également noter que pour simplifier la description, l'exemple particulier est illustré à partir de données particulières qui ne limitent pas l'invention. D'autres données de puits ou sismiques peuvent être utilisées, comme par exemple des données sismiques multi composantes. REVENDICATIONS1. Méthode pour évaluer des pressions de fluide dans une zone du sous-sol à partir de données de puits et de données sismiques, caractérisée en ce qu'elle comporte les étapes suivantes : on construit un cube de vitesse d'onde sismique P et un cube de vitesse d'onde sismique S au moyen d'une inversion stratigraphique desdites données sismiques, et on déduit desdites vitesses un cube de lithologies identifiant des lithologies argileuses et des lithologies non argileuses ; on détermine, à partir des données de puits et pour chacune des deux lithologies, une relation permettant d'estimer la pression de fluide à partir des vitesses d'onde sismique P ; on évalue les pressions de fluide dans ladite zone du sous-sol en construisant un cube de pressions de fluide en appliquant lesdites relations audit cube de vitesse d'onde sismique P en fonction dudit cube de lithologies. 2. Méthode selon la revendication 1, dans laquelle lesdites données sismiques comportent au moins un cube sismique discrétisant ladite zone en volumes élémentaires repérés par leurs coordonnées horizontales (x, y) et verticale en temps (t), et dans laquelle on construit le cube de lithologies en déterminant un cube de vitesse d'onde sismique S au moyen de ladite inversion stratigraphique, et en appliquant une valeur seuil de rapport des vitesses d'onde sismique P et S, de façon à ce que l'on affecte une lithologie argileuse aux volumes élémentaires ayant un rapport supérieur audit seuil, et une lithologie non argileuse aux autres volumes. 3. Méthode selon la revendication 2, dans laquelle le seuil est égale à 2. 4. Méthode selon l'une des revendications précédentes, dans laquelle on détermine lesdites relations sont de la forme suivante : Pe (z) = P" (z . n ( V z) pore pore ) VP (z) avec : Vp (z) : vitesse des ondes sismiques P mesurée aux puits à une profondeur z VP (z) : vitesse des ondes sismiques P estimée dans une hypothèse où il n'existe aucune surpression à la profondeur z Pore (z) : pression de fluide mesurée aux puits à une profondeur z Ppore(z) : pression de fluide estimée dans une hypothèse où il n'existe aucune surpression à la profondeur z 5. Méthode selon la revendication 4, dans laquelle on détermine VP (z) dans les lithologies argileuses au moyen des étapes suivantes : - on identifie un intervalle en profondeur où la pression de fluide est voisine de la pression hydrostatique, et on définit ladite relation VP (z) par une relation linéaire sur cet intervalle. 6. Méthode selon l'une des revendications 4 et 5, dans laquelle on détermine dans les lithologies non argileuses au moyen des étapes suivantes : on mesure la pression de fluide au moyen de diagraphie en un nombre réduit de profondeurs z, et pour les mêmes profondeurs z, on calcule une pression de fluide égale à p.g.z , avec p 1030kg / m3 et g 9,81m/ s2 ; on mesure la vitesse sismique des ondes P à ces mêmes profondeurs ;on déduit VP (z) en interpolant une droite entre des profondeurs où la pression de pore a pu être mesurée. 7. Méthode selon l'une des revendications précédentes, dans laquelle on détermine des zones de surpressions au sein de ladite zone, en construisant un cube de pression de confinement Pconf(x, y, t) et en détectant les zones de suppression lorsque la pression de fluide est supérieur à a*Pconf(x, y, t), où a est un seuil préalablement choisi. 8. Méthode selon la revendication 7, dans laquelle le seuil a est égal à 0,9.