Brevet FR3067746A1 — publié le 20181221

ESTIMATION DES PROPRIÉTÉS MÉCANIQUES DE MILIEUX TRANSVERSALEMENT ISOTROPES Référence croisée aux applications apparentées [0001] La présente demande revendique la priorité sur la demande provisoire américaine no. de série 62/520,402 intitulée «ESTIMATION DES PROPRIÉTÉS MÉCANIQUES DE MILIEUX TRANSVERSALEMENT ISOTROPES », déposée le 15 juin 2017, dont l'intégralité est incorporée ici à titre de référence. DOMAINE TECHNIQUE [0002] La présente description concerne généralement la détermination des propriétés mécaniques, et plus particulièrement, par ex., sans limitation, les systèmes et les procédés pour l'extraction du gaz et du pétrole impliquant la détermination des propriétés mécaniques des milieux anisotropes basé sur les mesures d'une formation. Contexte DE LA DIVULGATION [0003] Les milieux anisotropes polaires, qui sont communément appelés milieux transversalement isotropes (TI), possèdent un axe de rotation à boucle fermée et un jeu infini d'axes doubles perpendiculaires à celui-ci. Un plan de symétrie existe perpendiculairement à l'axe à boucle fermée. Les types de milieux sont connus comme transversalement isotrope (TI), avec des noms alternatifs comme isotrope transverse vertical (VTI), isotrope transverse horizontal (HTI) et isotrope transverse incliné (TTI). Le VTI est un modèle approprié pour décrire des milieux stratifiés non fracturés et en particulier le schiste. Les formations schisteuse constituent environ 75 % des bassins sédimentaires, et ceci rend le VTI le modèle anisotrope le plus courant dans la sismologie exploratoire. Les schistes jouent un rôle important dans l'écoulement de fluide et l'imagerie sismique en raison de leur faible perméabilité et leurs propriétés anisotropes. [0004] Un tenseur de rigidité TI contient cinq constantes élastiquement indépendantes et fournit des informations utiles pour aider les diverses opérations de puits, telles que la détermination des emplacements pour forer un puits horizontal ou identifier des intervalles pour perforer. Malheureusement, il est difficile d'obtenir de façon fiable le tenseur de rigidité à l'aide des procédés classiques ou l'obtention nécessite des mesures qui ne sont pas disponibles à toutes les profondeurs dans un puits. 2 BRÈVE DESCRIPTION DES FIGURES [0005] La figure 1 A est un schéma d'un exemple d e système de forage qui peut utiliser les principes de la présente divulgation. [0006] La figure 1B est un schéma d'un exemple d e système câblé qui peut utiliser les principes de la présente divulgation. [0007] La figure 2 est un schéma d'un exemple d'outil de diagraphie qui peut utiliser les principes de la présente divulgation. [0008] La figure 3 est un schéma illustrant un exemple de symétrie anisotrope polaire. [0009] La figure 4 est un schéma illustrant des exemples de rapports de Poisson. [0010] Les figures 5A à 51 sont des graphiques illustrant un exemple de données d'entrée et des erreurs résultantes. [0011] Les figures 6A à 6D sont des graphiques illustrant des exemples de sensibilités de C13 à différentes vitesses. [0012] Les figures 7A à 7C sont des schémas illustrant des exemples de rapports de Poisson. [0013] La figure 8 est un graphique illustrant un exemple de tracé sous forme d'une fonction de C13. [0014] La figure 9 est un graphique illustrant des exemples de limites supérieure et inférieure. [0015] La figure 10 est un graphique illustrant des exemples de précisions de l'estimation C13. [0016] Les figures 11A à 111 sont des graphiques illustrant un exemple de données d'entrée et des erreurs résultantes. [0017] Les figures 12A à 12F sont des schémas illustrant des exemples de rapports de Poisson. [0018] Les figures 13A à 13F sont des graphiques illustrant des exemples de paramètres de Thomson. [0019] La figure 14 est un graphique illustrant un exemple d'estimation de 3 à l'aide d'une analyse de régression. [0020] La figure 15 est un organigramme illustrant un exemple de procédé pour calculer les propriétés mécaniques d'un milieu anisotrope. [0021] Dans une ou plusieurs implémentations, tous les composants illustrés dans chaque figure ne seraient pas nécessaires, et une ou plusieurs implémentations 3 peuvent comprendre des composants additionnels qui ne sont pas illustrés dans une figure. Des variations au niveau de l'agencement et du type de composants peuvent être apportées sans s'écarter de l'esprit ou de la portée de la présente divulgation. Des composants supplémentaires, des composants différents ou un nombre plus petit de composants peut être utilisé dans le cadre de la portée de la présente divulgation. DESCRIPTION DÉTAILLÉE [0022] La description détaillée présentée ci-dessous est destinée à être une description de diverses implémentations et elle n'est pas destinée à représenter seulement les implémentations dans lesquelles la technologie décrite peut être pratiquée. Comme le comprendront les spécialistes du domaine, les implémentations décrites peuvent être modifiées de diverses façons différentes, sans s'écarter de la portée de la présente divulgation. Par conséquent, les figures et la description doivent être considérées comme de nature illustrative et non pas restrictive. [0023] La présente description concerne généralement la détermination des propriétés mécaniques, et plus particulièrement, par ex., sans limitation, les systèmes et les procédés pour l'extraction de gaz et de pétrole impliquant la détermination des propriétés mécaniques des milieux anisotropes, tel que les milieux transversalement isotrope, orthotrope et/ou orthotrombique, basée sur les mesures d'une formation. Systèmes et procédés permettant de déterminer les propriétés mécaniques des milieux anisotropes. [0024] Les principes de cette divulgation peuvent être utilisés, par ex., dans divers types de systèmes de puits, tels que des systèmes de forage, des systèmes de complétion et des systèmes câblés pour des puits qui peuvent être utiles pour l'extraction du gaz et/ou du pétrole. La figure 1 A est un schéma d'un exemple de système de forage 100 qui peut utiliser les principes de la présente divulgation, selon un ou plusieurs modes de réalisation. Tel qu'illustré, le système de forage 100 peut comprendre une plateforme de forage 102 positionnée au niveau de la surface de la terre et un puits de forage 104 qui se prolonge de la plateforme de forage 102 vers une ou plusieurs formations souterraines 106. Dans d'autres modes de réalisation, tels que dans une opération de forage offshore ou sous-marine, un volume d'eau peut séparer la plateforme de forage 102 et le puits de forage 104. [0025] Le système de forage 100 peut comprendre un derrick 108 supporté par la plateforme de forage 102 et comportant une moufle mobile 110 permettant de remonter et de descendre la colonne de forage 112. Un kelly 114 peut soutenir la colonne de forage 112 lorsqu'elle est descendue à travers une table tournante 116. Un trépan de forage 118 4 peut être couplé à la colonne de forage 112 et alimenté par un moteur au fond du puits et/ou par rotation de colonne de forage 112 par la table tournante 116. Lorsque le trépan de forage 118 tourne, il crée le puits de forage 104, qui pénètre les formations souterraines 106. Une pompe 120 peut faire circuler du fluide de forage à travers un tuyau d'alimentation 122 et le kelly 114, vers le fond du trou à l'intérieur de la colonne de forage 112, à travers des orifices dans le trépan de forage 118, et de retour vers la surface à travers l'anneau défini autour de la colonne de forage 112, et dans un bassin de rétention 124. Le fluide de forage refroidit le trépan de forage 118 au cours du fonctionnement et transporte les déblais du puits de forage 104 jusque dans le bassin de rétention 124. [0026] Le système de forage 100 peut également comprendre un ensemble de fond de trou (BHA) couplé à la colonne de forage 112 à proximité du trépan de forage 118. Le BHA peut comprendre divers outils de mesure de fond de trou tels que, sans limitation, des outils de mesure pendant le forage (MWD) et de diagraphie pendant le forage (LWD), qui peuvent être conçus pour prendre des mesures de fond de trou des paramètres apparentés au puits, tels que les conditions de forage, les propriétés de la formation, etc. Les outils MWD et LWD peuvent comprendre au moins un outil de diagraphie 126, qui peut comprendre un ou plusieurs capteurs ayant une pluralité d'éléments de détection capable de collecter des mesures ou des données de journal utiles pour la détermination des paramètres du puits de forage, y compris les propriétés mécaniques de la formation. Dans certains modes de réalisation, l'outil de diagraphie 126 est implémenté sous forme de ou autrement comprend un outil de diagraphie sonique ayant un ou plusieurs émetteurs acoustique et/ou un ou plusieurs récepteurs acoustique conçus pour mesurer et collecter des données qui peuvent être analyser pour fournir des mesures de la vitesse de l'onde élastique d'un milieu anisotrope dans la formation 106. Dans le présent contexte, un «outil de diagraphie sonique» décrit un quelconque outil de diagraphie conçu pour collecter des données de forme d'onde élastique basé sur les acoustiques, qui peuvent être dans des fréquences audibles et/ou inaudibles. Les émetteurs et/ou récepteurs peuvent être conçus pour émettre et/ou recevoir des ondes se propageant dans de multiples directions différentes (par ex., verticale, horizontale et/ou hors axe) et/ou se propageant avec de multiples modes différents (par ex., monopôle, dipôle, quadrupôle). L'analyse des différents modes à l'aide des procédés standards peut fournir les mesures de l'onde élastique souhaitées, telles que des vitesses de compression et de cisaillement, à travers la formation d'intérêt. [0027] Lorsque le trépan de forage 118 prolonge le puits de forage 104 à travers les formations 106, l'outil de diagraphie 126 peut recueillir des mesures qui peuvent être utilisées pour estimer les propriétés mécaniques des formations 106. L'outil de diagraphie 126 et d'autres capteurs des outils MWD et LWD peuvent être couplés en communication à un module de télémétrie 128 utilisé pour transférer des mesures et des signaux à partir du BHA vers un récepteur en surface (non illustré) et/ou pour recevoir des commandes provenant du récepteur en surface. Le module de télémétrie 128 peut englober un quelconque moyen connu de communication de fond de trou comprenant, sans limitation, un système de télémétrie par impulsion dans la boue, un système de télémétrie acoustique, un système de communication sur fil, un système de communication sans fil, ou une quelconque combinaison de ceux-ci. Dans certains modes de réalisation, le module de télémétrie 128 peut être omis et la colonne de forage 112 peut comprendre à la place un tube de forage câblé ou un tube enroulé câblé utilisé pour transférer des données à travers des conducteurs câblés vers un récepteur en surface. Dans certains modes de réalisation, certaines ou toutes les mesures prises par l'outil de diagraphie 126 peuvent être stockées à l'intérieur de l'outil de diagraphie 126 ou du module de télémétrie 128 pour une récupération ultérieure au niveau de la surface lors de la récupération de la colonne de forage 112. 100281 À divers moments au cours du procédé de forage, la colonne de forage 112 peut être enlevée du puits de forage 104 tel qu'il est illustré dans la figure 1B, pour réaliser des opérations de mesure/de diagraphie. Plus particulièrement, la figure 1B illustre un schéma d'un exemple de système de forage câblé 200 qui peut utiliser les principes de la présente divulgation, selon un ou plusieurs modes de réalisation. Les chiffres semblables utilisés dans les figures 1 A et 1B se rapportent aux mêmes composants ou éléments et, par conséquent, pourraient ne pas être de nouveau décrits en détail. Tel qu'illustré, le système câblé 200 peut comprendre une sonde d'instrument câblé 202 qui peut être suspendue dans le puits de forage 104 par un moyen de transport 204. Même si le moyen de transport 204 est illustré sous forme d'un câble dans la figure 1B, dans divers modes de réalisation, un moyen de transport peut comprendre, par ex., une ligne câblée, un câble lisse, un tube de forage, un tubage enroulé, un tracteur de fond de trou, ou une combinaison de ceux-ci. La sonde d'instrument câblé 202 peut comprendre l'outil de diagraphie 126, qui peut être couplé en communication au moyen de transport 204. Dans divers modes de réalisation, le moyen de transport 204 peut transporter la télémétrie et/ou le courant ou pas. Par exemple, le moyen de transport 204 peut comprendre des conducteurs pour transporter du courant à la sonde d'instrument câblé 202 et également pour faciliter la communication entre la surface et la sonde d'instrument câblé 202. Une unité de diagraphie 206, illustrée dans la figure 1B sous forme d'un camion, peut recueillir les mesures provenant de l'outil de diagraphie 126, et peut comprendre des composants informatiques 208 pour commander, traiter, stocker et/ou visualiser les mesures recueillies par l'outil de diagraphie 126. Les composants informatiques 208 peuvent être couplés en communication à l'outil de diagraphie 126 avec le moyen de transport 204. 6 [0029] Même si les composants informatiques 208 sont illustrés en haut du trou dans la figure 1B, dans divers modes de réalisation de la présente divulgation, les procédés décrits ici peuvent être implémentés en haut du trou, au fond du trou, ou les deux. Par exemple, un système de traitement comprenant un ou plusieurs processeurs et/ou mémoires conçues pour implémenter l'un quelconque des procédés pour calculer les propriétés mécaniques décrites ici peut être placé au fond du trou dans le puits de forage 104 (par ex., dans l'outil de diagraphie 126 ou la sonde d'instrument 202), en haut du trou dans l'unité de diagraphie 206, ou une combinaison de ceux-ci, peut être conçu pour implémenter ces procédés à l'aide de techniques de traitement réparties. [0030] La figure 2 est un schéma d'un exemple d'outil de diagraphie 126 plus en détail, selon certains modes de réalisation. L'outil de diagraphie 126 illustré dans la figure 2 comprend généralement un ou plus émetteurs TX et un ou plusieurs récepteurs RX. Les émetteurs TX sont généralement conçus pour émettre des ondes dans un milieu anisotrope d'intérêt, par ex., la formation 106 (voir les figures 1 A à 1B), alors que les récepteurs RX sont généralement conçus pour mesurer une réponse correspondante aux ondes émises après interaction des ondes avec le milieu d'intérêt. Comme exemple, le ou les émetteurs TX peuvent comprendre un ou plusieurs émetteurs acoustiques (par ex., émetteurs piézoélectriques) qui agissent comme une source de sons audibles et/ou inaudibles pour générer des ondes élastiques dans le milieu anisotrope, et le ou les récepteurs RX peuvent comprendre un ou plusieurs récepteurs acoustiques (par ex., récepteurs piézoélectriques) qui captent les ondes retournées émises par le ou les émetteurs acoustiques. Dans certains modes de réalisation, les émetteurs TX comprennent une pluralité d'émetteurs de deux ou plusieurs types différents (par ex., choisis parmi les émetteurs monopôles, les émetteurs dipôles et les émetteurs quadripôles) pour permettre la détermination des diverses vitesses différentes divulguées ici. Même si seulement un émetteur TX est illustré dans la figure 2, dans divers modes de réalisation l'outil de diagraphie 126 peut comprendre un quelconque nombre approprié d'émetteurs. Par exemple, l'outil de diagraphie 126 peut comprendre 2, 3, 4, 5 ou plus d'émetteurs. Même si seulement deux récepteurs RX sont illustrés dans la figure 2, dans divers modes de réalisation l'outil de diagraphie 126 peut comprendre un quelconque nombre approprié de récepteurs. Par exemple, l'outil de diagraphie 126 peut comprendre 2, 3, 4, 5 ou plus de récepteurs. [0031] Comme le démontre la figure 2, l'outil de diagraphie 126 peut généralement avoir la forme d'un composant allongé définissant un axe longitudinal 127, qui peut être placé au fond du trou dans un puits de forage 104 (voir les figures 1 A à 1B). Le ou les émetteurs TX sont généralement axialement espacés le long de l'axe longitudinal par rapport aux récepteurs RX pour permettre la capture des données souhaitées. Dans l'exemple illustré, de 7 multiples récepteurs sont compris, qui sont également axialement espacées les uns par rapport aux autres, pour permettre la capture d'informations de temps de retard provenant des ondes émises. [0032] Dans certains modes de réalisation, les données' provenant de l'outil de diagraphie 126 peuvent être utilisées pour prédire les propriétés mécaniques de la formation. La prédiction des valeurs justes des rapports de Poisson (y) et du module de Young (E), qui peut être utilisée pour calculer le gradient de la fracture et la contrainte horizontale minimale, est utile pour la géométrie mécanique, la fracturation hydraulique et la complétion. La contrainte horizontale minimale est utile pour sélectionner l'endroit où il faut s'installer et perforer. Avec des prédictions améliorées des modules d'élasticité, la contrainte horizontale minimale et maximale peut être mieux évaluée. Un profil de contrainte versus profondeur associé à la friabilité de la roche prédite à partir des profils des modules d'élasticité est utile pour aider à choisir les «points idéaux» pour forer un puits horizontal, et également pour déterminer les intervalles de perforation. [0033] Dans certains modes de réalisation, un procédé est proposé pour le calcul précis des modules de Young, du rapport de Poisson et d'autres propriétés mécaniques et sismiques du milieu avec superposition (ou stratification ou fractures). Dans certains modes de réalisation, les procédés décrits ici permettent d'avoir des estimations hautement précises tout en réduisant le nombre d'entrées, qui peut être difficile à avoir même à partir des mesures par carottage. Dans certaines applications, les carottes ne sont pas acquises à toutes les profondeurs ; par conséquent, il existe un besoin pour un procédé qui est moins affecté par le manque de données de carottage. [0034] Les approches classiques du calcul de ces propriétés utilisent des mesures à des angles de 0, 90 et 45 degrés (le troisième peut également être un angle arbitraire hors axe, entre 0 et 90) par rapport à l'axe de symétrie du milieu. Dans certains modes de réalisation, les procédés décrits ici ne demandent pas les mesures à 45° (ou autres mesures hors axe). [0035] La caractérisation des propriétés élastiques du milieu VTI a été d'un grand intérêt durant les trois 3 dernières décennies pour plusieurs applications telles que l'imagerie sismique, le forage et les géomécaniques de complétion. Même s'il y a eu un bon progrès dans la prise en compte des modèles TI pour les applications sismiques, ces modèles sont toujours faiblement compris pour d'autres applications telles que la géomécanique. Pour des applications géomécaniques et de forage, les suppositions idéalisées pour les modèles TI tels que l'homogénéité, l'élasticité et la dépendance à l'échelle, sont violées plus que dans les applications sismiques. Pour les applications sismiques, plusieurs de ces problèmes sont moins 8 importants, en partie en raison des longueurs d'ondes sismiques plus grandes qui peuvent ignorer les hétérogénéités à petite échelle et en partie parce que les roches restent dans la plage élastique, contrairement aux applications géomécaniques. [0036] Lors de la prise en compte des modèles TI, selon au moins certains modes de réalisation divulgués ici, les soi-disant propriétés mécaniques dynamiques (par ex. ,modules de Young et rapports de Poisson) dérivées des mesures de la vitesse d'onde dans les roches peuvent être apparentées aux propriétés mécaniques statiques, qui sont plus pertinentes aux applications d'ingénierie. Ceci peut être appelé un problème de correction dynamique-statique. Le fait d'ignorer l'anisotropie causée par la stratification et une organisation en couches fines (VTI) peut entraîner des erreurs et des problèmes importants dans la compréhension de la vraie relation entre les propriétés dynamique et statique. Ce problème semble être plus grave pour les rapports de Poisson étant donné qu'il est défini comme une fraction et les valeurs dynamiques obtenues peuvent être très différentes des valeurs réelles. Comme il est décrit ici, de petites incertitudes au niveau des vitesses mesurées peuvent entraîner de très grandes erreurs dans le calcul des rapports de Poisson. [0037] Dans certains modes de réalisation, le fait d'ignorer le modèle VTI pour les roches qui ont réellement une symétrie VTI affectera l'investigation de la propriété dynamique-statique. Il existe une différence importante entre le VTI et les propriétés mécaniques isotropes, ce qui aurait pu empêcher l'observation antérieure d'une quelconque relation significative pour les rapports de Poisson (et peut-être les modules de Young). La seule propriété, qui semble la moins affectée par le fait d'ignorer le VTI, est l'estimation du module de Young Ev, qui semble être proche de ce qui serait observé par la supposition isotrope. [0038] Le concept du rapport de Poisson (définition isotrope et généralement calculé à partir des vitesses d'onde) a été un sujet controversé dans la géophysique, certains suggérant qu'il ne doit pas être utilisé pour des applications géophysiques en raison du fait qu'il est calculé à l'aide de vitesses alors que la vraie valeur du rapport de Poisson doit être mesurée par la charge mécanique réelle et l'utilisation des jauges de contrainte. Cependant, le rapport de Poisson peut être un outil utile pour caractériser les propriétés mécaniques de la roche à travers, par exemple, la friabilité. Dans certains modes de réalisation, une mesure soigneuse des rapports de Poisson dynamiques peut préserver les relations relatives entre les rapports de Poisson VTI et même procurer des résultats dynamiques, qui sont proches des mesures statiques. Une supposition fondamentale ici est que différents effets tels que la supposition de dispersion et d'élasticité seront annulés lors du calcul des rapports de Poisson. [0039] Selon certains modes de réalisation, en utilisant les relations qui doivent être maintenues pour les rapports de Poisson dans le milieu VTI, une fourchette de 9 variation pour C13 peut être définie. Les limites pour C13 sont relativement étroites pour plusieurs échantillons de roches. Pour les mesures dans lesquelles les relations entre les rapports de Poisson VTI sont violées, ceci peut être attribué à des erreurs dans C13 étant donné que les autres quatre composants de la matrice de rigidité VTI sont mesurés directement alors que C13 est calculé à partir d'autres vitesses et implique des mesures à 45°. Comme il est décrit ici, C13 est très sensible aux erreurs dans d'autres vitesses particulièrement les mesures de vitesse à 45°. [0040] La notation de Thomsen et les paramètres c, y et 6 ont été utilisés pour plusieurs applications, en raison de leur simplicité. La notation de Thomsen et les paramètres sont décrits, par ex., dans Thomsen, L. ,1986, "Weak Elastic Anisotropy," Geophysics 51(10), 1954-1966. Même s'il est largement accepté que s et y sont étroitement liés (y dans plusieurs cas légèrement plus élevés que c), il n'y a pas eu une quelconque relation significative entre 6 et les deux autres paramètres. Tel qu'il est décrit ici, cette ambiguïté réside dans les complications au niveau de la mesure des vitesses hors axe (et ultérieurement C13) utilisées pour calculer 8 mais également en raison de la définition compliquée de 8 elle-même. Si Cnest calculé avec précision, 8 semble avoir une relation significative à E et y. Théorie [0041] L'élasticité linéaire peut décrire le comportement d'une diversité de matériaux, tels que les milieux constituant les schistes et d'autres formations. Les relations stress-contrainte pour un matériau élastique linéaire en 3D sont complexes étant donné que le matériau peut être soumis à une diversité de différentes contraintes compressives et de cisaillement dans différentes directions autour d'un quelconque point donné. [0042] La loi de Hooke généralisée capte cette complexité en modélisant le stress (u) et la contrainte (E) comme des tenseurs de second ordre (ai» Eii) dans lesquels chaque composant de stress dépend linéairement de chaque composant de contrainte. Cette relation entre le stress et la contrainte est caractérisée par une constante de matériau appelée la rigidité (C) qui peut être représentée par un tenseur de quatrième ordre Cipa qui définit le mappage linéaire entre les deux tenseurs de second ordre correspondant au stress et à la contrainte (o-ii = Cijkleii). Il est à noter que l'inverse de la rigidité est appelé l'élasticité (S), et il doit être compris qu'un quelconque calcul de l'élasticité peut généralement être considéré comme équivalent au calcul de la rigidité dans le cadre de cette description. [0043] Mathématiquement, les tenseurs de stress et de contrainte peuvent être représentés dans un système de coordonnées cartésiennes comme ce qui suit, dans lequel les indices 1, 2, 3 correspondent aux 3 axes dans le système de coordonnées (par ex., x = 1, y = 2, z = 3) : au = e11 cr12 an ; Ei = En En £13 (1) 021 F#531 az2 "23 i £21 £22 E23 0-32 C1-33 £31 E32 £33 [0044] Comme mappage entre des matrices 3 par 3, le tenseur de rigidité Ciiki peut être représenté par une matrice 3 par 3 par 3 par 3 ayant 81 composants. Les symétries inhérentes de ces propriétés mécaniques permettent de simplifier considérablement ces équations et relations. La notation de Voigt procure un mappage standard pour les indices de tenseur et permet la réduction des tenseurs symétriques pour le stress, la contrainte et la rigidité en deux tenseurs de premier ordre et un tenseur de second ordre, respectivement. Ce mappage est illustré dans l'équation suivante, qui omet les composants de stress et de contrainte redondants en raison de la symétrie (o-23 = (1329E13 = E31, - C11 .. C13 C14 C15 C16 C26 C36 C46 C56 C66 -E1 = E11E2 = £22 (2) -el = ell- (3-2 = a22 = C21 C31 C12 C23 C24 C25 £3 = £33 0-3 = 0-33 C41 L41 C22 C33 C34 C35 £4 = £23 a 4 = 1723 as = (713 (76 = C512- C51 C32 C43 C53 C44 C45 £ E 61 C42 C63 C54 C55 C65 5 = l3 C52 C64 -E6 = E12- C62 36 composants sont illustrés pour la matrice de rigidité précédente, qui sont quelquefois appelés coefficients de rigidité ou constantes élastiques. Alors que 36 composants de rigidité sont illustrés, la symétrie de la matrice de rigidité veut dire que seulement 21 de ses composants sont indépendants pour le cas le plus général de l'élasticité anisotrope. [0045] La matrice de rigidité peut davantage être simplifiée pour une diversité de types d'isotropie et d'anisotropie. Par exemple, pour les milieux isotropes la matrice de rigidité peut être réduite à seulement deux composants indépendants correspondant aux changements en volume et les déformations de cisaillement. Pour divers types d'anisotropie (isotropie transverse, orthotropique, orthothrombique, et autres), la matrice de rigidité peut être simplifiée à moins de 21 composants indépendants. [0046] Les milieux anisotropes polaires, qui sont communément appelés milieux transversalement isotropes (TI), possèdent un axe de rotation à boucle fermée et un jeu infini d'axes doubles perpendiculaires à celui-ci. Un plan de symétrie existe perpendiculairement à l'axe à boucle fermée. Les types de milieux sont connus comme transversalement isotrope (TI), avec des noms alternatifs comme isotrope transverse vertical (VTI), isotrope transverse horizontal (HTI) et isotrope transverse incliné (TTI). [0047] La figure 3 illustre un exemple de milieu TI 351 dans lequel l'axe de rotation à boucle fermée correspond à l'axe de symétrie 353 et se prolonge dans une direction verticale, perpendiculaire à la direction de stratification 355 du milieu. Un tenseur de rigidité TI contient cinq constantes élastiquement indépendantes. Le VTI est un modèle approprié pour 11 décrire des milieux stratifiés non fracturés et en particulier le schiste. Les formations schisteuse constituent environ 75 % des bassins sédimentaires, et ceci rend le VTI le modèle anisotrope le plus courant dans la sismologie exploratoire. Les schistes jouent un rôle important dans l'écoulement de fluide et l'imagerie sismique en raison de leur faible perméabilité et leurs propriétés anisotropes. Il sera compris que même si des exemples sont décrits ici en référence au milieu VTI, les principes de cette divulgation peuvent être étendus à divers autres types de milieux, s'il y a lieu. [0048] La matrice de rigidité pour le milieu VTI dans un système de coordonnées dans lequel les directions xl et x2 sont dans un plan horizontal et x3 est dans la direction verticale (voir la figure 3) a la forme de : CVTI = C11 Cl , - 2 C66 C11 C13 C13 C33 0 0 0 0 0 0 0 0 C55 0 0 (3) C11 - 2G66 Ci 3 0 0 Ci 3 0 0 0 0 0 C' 0 0 0 0 0 0 0 C66 [0049] Les composants de la matrice de rigidité VTI, C11, C33, C55 et C66 peuvent être définis en termes de vitesses de compression et de cisaillement verticales et horizontales. Il est à noter que C44 = C55 pour le milieu VTI, et, par conséquent, la matrice précédente peut de façon équivalente être exprimée à l'aide de C44. C13, d'autre part nécessite généralement des mesures de vitesse hors axe. Pour des raisons de simplicité, la vitesse de l'onde de compression à 45° est utilisée dans l'équation pour calculer (4) C13. C13 = -\./(2pV45 - C11 - C55) ( 2 p1145 - C33 - C55) - C55 (4) [0050] Les paramètres de Thomsen sont également définis en termes de E = 2C33 C66 - C55 Y - 2 C55 (C13 +C55)2 -(C33 -055)2 8 = (7) 2c33(c33-c55) [0051] Les rapports de Poisson anisotropes (vi1) pour le milieu VTI sont définis par l'équation générale suivante : (5) (6) (8) [0052] En remplaçant le marquage des axes de 1,2 par H,h et 3 par V, les axes définissent : 12 vv = VHH = VHV = [0053] Les trois rapports de Poisson dans un milieu VTI sont illustrés dans la figure 4, dans chaque cas, la première flèche indique la direction du stress appliqué et la flèche suivante illustre la direction de l'expansion, orthogonale au stress appliqué. Il est à noter que dans le cas VHH, Eh est utilisé pour distinguer les directions 1 et 2 parce que les valeurs de contrainte horizontale sont différentes lors de l'évaluation VHH, sinon on obtiendrait VHH = 1 qui n'est pas juste. [0054] En termes de Cii, les rapports de Poisson VTI aussi bien que le module de Young vertical (Ev) et le module de Young horizontal (EH) sont, 2C13 C66 VHV - Cil C33 -CL C33 (Cil - 2 C66) - CL VHH Ci1C33 -C12.3 C13 VV 2(C11-C66) E C33 (Ci -C66)- CL = V Cil -C66 4C66 (C33 (C11 - C66) - CL) Cil C33 -CL [0055] Comme le démontre les équations précédentes, C13 entre dans le calcul de toutes ces propriétés mécaniques aussi bien que le paramètre de Thomsen ô. Par conséquent, les incertitudes dans C13 affecteront toutes ses propriétés. Le calcul C13 à l'aide de l'équation (4) est très sensible aux incertitudes dans les vitesses, particulièrement à la vitesse à 45°. [0056] En se reportant aux figures 5A à 51, une simulation de Monte- Carlo simple révèle les effets des petites erreurs dans les mesures de vitesse pour le calcul de C13, Ô, et des rapports de Poisson. Les figures 5A à 51 illustrent des exemples de graphiques 500a-500i pour démontrer les sensibilités basées sur les mesures provenant d'un jeu de données pour un échantillon de roche. Les graphiques 500a-500i illustrent un exemple de la façon dont les paramètres et les propriétés mécaniques sont sensibles aux variations ou aux imprécisions dans les vitesses. Les graphiques 500a à 500i illustrés dans les figures 5A à SI illustrent les résultats d'une expérience réalisée à l'aide de données provenant d'un échantillon de roche, mais il doit être noté que cette expérience a été répétée pour beaucoup d'autres échantillons de roche et que les mêmes observations ont été faites pour toutes les expériences. EH (9) eV (10) EH Eh eV EH EH= (12) (13) (14) (15) (16) 13 [0057] Les graphiques 500a-d illustrés dans les figures 5A à 5D contiennent des histogrammes basés sur quatre valeurs de vitesse mesurées pour l'échantillon, respectivement (c. -à-d. , vitesse de compression verticale Vp0=3350 m/s, vitesse de compression horizontale Vp90-5533 m/s, vitesse de compression à 450 Vp45-4360 m/s et vitesse de cisaillement horizontale Vs90=3246 m/s). Les barres horizontales 590a-d illustrent ces valeurs de vitesse mesurées. 100 valeurs de bruit ont été générées pour chaque vitesse, le bruit maximal ajouté étant égal à 1 % de chaque grandeur de vitesse. Ces valeurs sont illustrées sous forme de barres verticales dans les figures 5A à 5D et ont été prises comme entrées pour l'analyse. [0058] Dans les figures 5E à 5I, les équations (4, 7, 12, 13 et 14) ont été utilisées pour calculer C13, 8 et les rapports de Poisson VTI pour chacune des 100 permutations. Les résultats de ces calculs sont rapportés dans les graphiques 500e-i. Dans la figure 5E, le graphique 500e illustre les valeurs de C13 calculées pour chaque permutation sous forme de barres verticales et illustre la valeur de C13 attendue sous forme de barre horizontale 590e. Les résultats sont représentés de la même façon dans les figures 5F à SI dans les graphiques 500f-i pour 8 et les rapports de Poisson VTI. Dans la figure 5F, le graphique 500f illustre les valeurs de 8 calculées pour chaque permutation sous forme de barres verticales et illustre la vraie valeur de 8 sous forme de barre horizontale 590f. Dans la figure 5G à 5I, les graphiques 500g-i illustrent les valeurs de viiv, vHH calculées, et les valeurs vv. pour chaque permutation sous forme de barres verticales et illustrent les vraies valeurs de vHv, vHH, et les valeurs de vv sous forme de barres horizontales 590g-i, respectivement. [0059] Comme le démontre la figure 5E, même une petite erreur maximale de 1 % peut entraîner des valeurs de C13 qui sont différentes de 10 gigapascals (GPa). La différence pour 8 peut être aussi élevée que 0,5 (voir la figure 5F), et d'environ 0,3 pour les rapports de Poisson (voir les figures 5G à 51). Le même taux d'erreur que C13, peut etre observé pour les modules de Young dynamiques (non illustré). Ces différences peuvent donner des résultats et des conclusions complètement différents concernant les propriétés élastiques de la roche. Dans certaines applications, les mesures à 45° sont les principales causes d'erreur dans l'estimation duC13, et par conséquent, les propriétés mécaniques de 8 et VTI. [0060] Les figures 6A à 6D illustrent un exemple de la façon dont les paramètres et les propriétés mécaniques sont sensibles à différentes vitesses. En particulier, les figures 6A à 6D illustrent des graphiques 600a-d contenant des exemples de résultats d'une analyse de sensibilité démontrant la façon dont les erreurs dans la vitesse hors axe Vp45 peuvent avoir un effet plus grand sur les C13 valeurs en comparaison à d'autres vitesses. 14 [0061] Dans cet exemple, les données ultrasoniques mesurées sur environ 450 échantillons de roche ont été utilisées pour l'analyse. Une fourchette de variation pour les vitesses a été définie pour chaque échantillon. La fourchette de la variation pour le nième échantillon était de [V' - 600 V' + 6001 mètre par seconde m/s. En particulier, pour chaque mesure de vitesse, 1201 valeurs ont été générées parmi lesquelles la valeur moyenne dans la série générée a été la valeur mesurée provenant des données. L'axe vertical illustre le numéro de l'échantillon de roche pour l'ensemble de données qui a été utilisé. Il y a environ 450 échantillons de roche (par conséquent, l'axe vertical va de 1 à environ 450). En allant le long d'une ligne horizontale de gauche à droite, la première valeur est V' - 600 m/s et la dernière valeur indiquée est V' + 600 m/s. Les valeurs mesurées pour les vitesses sont données par la courbe blanche qui passe à travers le milieu de chaque nuage. [0062] En se référant au graphique 600a illustré dans la figure 6A, pour voir l'effet de Vp0 sur C13 dans cet exemple, C13 a été calculé à l'aide de l'équation (4) en faisant varier Vp0 pour chaque échantillon dans la fourchette [V130 - 600 à V ip0 + 600], tout en n'utilisant que les valeurs mesurées pour le restant des vitesses (V1)90, Vp45 et Vs90) pour l'échantillon ième (dans lequel i est le numéro de l'échantillon de roche). Ceci a été répété pour les graphiques 600b (figure 6B), 600c (figure6C) et 600d (figure 6D) pour calculer les sensibilités de C13 to Vip90, Vp45 et Vs90, respectivement, dans chaque cas définissant une fourchette de variation pour la vitesse donnée, tout en maintenant inchangé les autres vitesses. Afin d'être en mesure de voir des tendances significatives, pour chaque tracé, les données ont été triées en se basant sur la vitesse donnée qui était variée. Par exemple, dans le C13 : le tracé Vp0 illustré dans la figure 6A, C13 est trié en augmentant Vp0. [0063] Les hachures montrent la valeur de C13 dans cette analyse de sensibilité. Par conséquent, une variation plus latérale (si elle est cohérente pour la majorité des 450 échantillons) dans chaque tracé signifie que C13 est plus sensible que le paramètre spécifique. Les résultats dans cet exemple montrent que C13 est plus sensible à Vp45 que d'autres vitesses parce que C13 change le plus lors de la variation de la Vp45. Ceci peut être observé dans la figure 6B, en allant de gauche à droite, un changement beaucoup plus abrupt dans les valeurs de C13 pour la Vp45 est observé. Ceci est cohérent pour presque tous les 450 échantillons. C13 a également une sensibilité relativement élevée à Vs0 ou C55. Même si elles ne sont pas aussi apparentes, ces conclusions peuvent également être tirées de l'équation (4). [0064] Dans certaines applications, Vp45 est la vitesse la plus problématique pour la mesure pour plusieurs raisons telles que les difficultés dans la préparation de l'échantillon avec un angle exact de 45° et dans certains cas des problèmes de vitesse de rayon versus groupe. Ajoutant à la complexité, les figures 6A à 6D demontrent que C13 peut être grandement affecté par de petites erreurs dans Vp45. La figure 5A à SI montre que ces erreurs in C13 peuvent donner des propriétés élastiques de la roche qui sont complètement différentes de la réalité. Par conséquent, dans certaines applications, à moins que toutes les vitesses (particulièrement les 45°) soient mesurées très précisément, des erreurs majeures dans le calcul des propriétés de la roche et de la caractérisation de l'anisotropie (ô) peuvent se produire. Les problèmes associés aux mesures de la Vp45, et conséquemment C13, sont quelquefois attribués aux problèmes de vitesse de groupe/de phase, particulièrement lorsque des transducteurs à point sont utilisés. Cependant, dans certaines applications, la correction de la vitesse de phase à la vitesse de groupe ne règle pas toujours le problème des mesures imprécises C13. Les limites des rapports de Poisson et C13 [0065] Les relations parmi les rapports de Poisson VTI peuvent être utilisées pour définir des limites supérieure et inférieure pratiques pour C13. Un jeu d'exemples de limites supérieure et inférieure pour C13 est illustré dans l'inégalité (17). La limite inférieure vient du fait que :1) C13 doivent être positifs, et 2) 0 >> EH. En outre, lorsque vv est calculé, la contrainte EH verticale Ey va au dénominateur (vv = - -), rendant vv très petit. Ceci montre que vv doit être Ey le rapport de Poisson le plus petit dans les roches VTI. Par conséquent, ce qui reste dans la recherche de la relation entre les trois rapports de Poisson VTI est de déterminer si VHH VHy. [0069] La figure 7B illustre une représentation du rapport de Poisson statique VHH dans le milieu VTI 351. Lors du calcul du VHH, le milieu VTI 351 est comprimé dans la direction horizontale h. h est utilisée comme la première direction horizontale pour le distinguer de H qui est la seconde direction horizontale. Il est important de faire cette distinction parce que lorsque la roche est comprimée dans la direction h,Eh et EH seront différents, sinon VHH = 1, ce qui n'a aucun sens physique. Lors du calcul de VHH, la roche est moins élastique dans la direction horizontale parce que cette fois-ci les roches ne sont pas poussées contre les plans de faiblesse. Par conséquent, une contrainte moyenne dans la direction h (Eh) est observée, et une contrainte encore plus petite (en comparaison à h) dans la direction H (EH) est observée. Par conséquent, Eh >> EH. Lors du calcul de VHH, Eh va au dénominateur (VHH = - -E2-), mais il Eh n'est pas aussi grand que ev lors du calcul de V. Par conséquent, VHH sera plus grand que vv. [0070] La figure 7C illustre une représentation du rapport de Poisson statique vHvdans le milieu VTI 351. Lors du calcul de vHv, si le milieu VTI 351 est comprimé dans la direction H, EH sera plus grand que ev (pour avoir un rapport de Poisson inférieur à 1). Cependant, ev, dans ce cas, est relativement grand (et va au numérateur) en raison des propriétés du milieu VTI et il est facile de les dilater dans la direction verticale lorsqu'on les comprime dans la direction horizontale. Par conséquent, vHv doit être le rapport de Poisson le plus grand dans un milieu VTI, et il devient évident que VHH 0 (20) dans laquelle les coefficients b, c, d peuvent comprendre les composants de rigidité, comme suit b = -2(C11 - C66) (21) 17 c = (22) d = 2C33 (C11 C66)(C11 2C66) (23) [0073] Si l'inégalité (18) est traitée comme une équation : + bq3 + cC1.3 + d = 0, (24) [0074] En se référant maintenant à la figure 8, pour voir le comportement de l'inégalité (20), les données provenant d'un échantillon de roche ont été prises, et tout en maintenant les autres Cii inchangés, on a fait varier, C13. La figure 8 est un graphique 800 illustrant les résultats de cette analyse, qui montre l'inégalité (20) rapportée en fonction de C13. La ligne verticale gauche 871 est la valeur mesurée de C13, et la ligne verticale 873 à la droite est la seconde racine de l'équation (24). [0075] En référence à la figure 8, comme un exemple, la fourchette acceptable pour C13 semble toujours être entre la première et la seconde racine de l'équation (24) et les valeurs mesurées de C13 semblent être toujours très proches de la seconde racine (à la gauche de la seconde racine). Cette expérience a été répétée sur plusieurs échantillons de roche et le même comportement a été observé. [0076] Selon certains modes de réalisation, ceci fournit un procédé pour déterminer une limite supérieure pratique ou estimer une valeur pour C13 basé sur l'équation (24). Par exemple, les trois racines de la fonction cubique (24) peuvent être résolues en utilisant, parmi d'autres procédés, le procédé trigonométrique. La solution générale a la forme suivante : tk = 7Icos (-cos-1 (3q \I 3) 2n-k 3 pour k = 0,1,2 3 2p p - 3 (25) où 3ac-b2 P = 3a2 2b3-9abc+27a2d = 27a3 Les solutions sont : C13 t k a [0077] La limite supérieure acceptable de C13 tel que mentionné précédemment, semble toujours se produire sur la gauche de la seconde racine (k=1) ; par conséquent, une limite supérieure pratique de C13 peut être définie comme : C113 = - 3a (29) (26) (27) (28) 18 [0078] Tel qu'il a été précédemment discuté, la limite inférieure pour C13 peut être basée sur l'équation (17) ; par conséquent une limite inférieure pratique de C13 peut être définie comme : Cf3 = N/C33C12 ± CL - C66 (30) [0079] La figure 9 illustre un exemple de graphique 900 dans lequel ces limites pour C13 sont déterminées pour l'exemple de jeu de données provenant des figures 5A à 51 correspondant à un échantillon de roche. Dans la figure 9, les limites déterminées sont rapportées avec les valeurs mesurées de C13. Dans la figure 9, la courbe 947 montre les valeurs mesurées de C13 provenant de l'ensemble de données. La courbe 949 représente la limite inférieure calculée à l'aide de l'équation (30) et la courbe 943 représente la limite supérieure nouvellement définie calculée à l'aide de l'équation (29). La courbe 941 représente la limite supérieure définie dans l'équation (17). Comme on peut le voir dans la figure 9, toutes les mesures C13 dans cet exemple, à l'exception d'une mesure (le 15' point de données), se trouvent à l'intérieur des limites 943 et 949 définies dans les équations (29 et 30). De façon avantageuse, les limites supérieure et inférieure (courbes 943 et 949) sont très proches l'une de l'autre. Par conséquent, ceci fournit un procédé alternatif pour l'estimation de C13 (et d'autres propriétés telles que Ô) simplement en calculant la moyenne des limites supérieure et inférieure, qui n'ont pas besoin des mesures de la Vp45. Cette moyenne est représentée par représentée par l'équation suivante : cL+ c U c = 13 13 2 [0080] Dans certains modes de réalisation, ce procédé de calcul C13 n'est utilisé que si le C13 mesuré (par ex., calculé à partir de la Vp45 à l'aide de l'équation 4) ne se trouve pas à l'intérieur des limites définies (équations 29 et 30). Selon au moins certains modes de réalisation divulgués ici, un procédé de contrôle de la qualité des calculs C13 peut également ou alternativement être utilisé comme un procédé pour l'estimation de C13, sans avoir à mesurer la Vp45 et utiliser l'équation (4). Par exemple, la nouvelle limite supérieure correspondant à l'équation (29) peut fournir un paramètre pour le contrôle de la qualité du C13 mesuré (calculé à partir de la Vp45 à l'aide de l'équation 4), et la même limite supérieure peut en outre ou par ailleurs être utilisée pour estimer C13 sans utilisation de la Vp45 mesurée pour calculer l'estimation à l'aide de l'équation (31). Résultats et applications [0081] La figure 10 est un graphique 1000 contenant un exemple d'histogramme pour démontrer les précisions de l'estimation C13 à l'aide de l'équation (31). Le (31) 19 graphique 1000 est basé sur environ 450 échantillons de roche provenant des données utilisées dans les figures 6A à 6D. Afin de générer le graphique 1000, C13 a été estimé à l'aide de l'équation (31) pour chacun des environs 450 échantillons de roche. Chaque barre correspond à une précision donnée du C13 estimé en comparaison à la valeur mesurée provenant de l'ensemble de données, alors que la hauteur de chaque barre montre le nombre de points de données (c. -à-d. , le nombre de valeurs C13 estimées) ayant la précision donnée. Comme on peut le voir dans la figure 10, une précision de ± 0,5 GPa est observée pour 135 des environs 450 échantillons de roche, une précision de ± 1 GPa est observée pour 100 échantillons, une précision de ± 1,5 GPa est observée pour 125 échantillons et une précision de ± 2 GPa est observée pour 40 échantillons. Par conséquent, dans cet exemple, C13 peut être estimé en calculant la moyenne des limites, avec une précision de ± 2 GPa pour 90 % des données et ± 1,5 GPa pour 80 % des données. [0082] Les figures 11A à 111 illustrent un autre exemple de jeu de graphiques 1100a-i basé sur la simulation de Monte-Carlo. Les graphiques 1100a-i démontrent la façon dont l'estimation de C13 en calculant la moyenne des limites, comme il est décrit ci-dessus, peut fournir un procédé qui est robuste par rapport aux erreurs dans les vitesses. En particulier, les graphiques 1100a-i montrent les valeurs déterminées de C13, de ô, et les valeurs de propriété mécanique qui sont moins éparpillées que ce qui est montré dans les figures 5A à SI. [0083] Les résultats montrés dans les figures 11E à 11I sont basés sur les mêmes entrées que celles utilisées dans les figures 5E à SI. Par conséquent, les graphiques 1100a-100d illustrés dans les figures 11A à 11D contiennent des histogrammes qui sont les mêmes graphiques 500a-500d illustrés dans les figures 5A à 5D. Les barres horizontales 590a-d illustrent les mêmes valeurs de vitesse mesurées que les barres horizontales dans les figures 5A à 5D. Les mêmes valeurs de bruit ont 100 été générées pour chaque vitesse, le bruit maximal ajouté étant égal à 1 % de chaque grandeur de vitesse, comme le démontre les barres verticales dans les figuresl 1 A à 11D et sont prises comme entrées pour l'analyse. [0084] Dans les figures 11E, l'équation (31) a été utilisée pour estimer C13 pour chacune des 100 permutations. Dans les figures 11F à 11I, les équations (7, 12, 13 et 14) ont été utilisées pour calculer 6 et les rapports de Poisson VTI pour chacune des 100 permutations, basé sur les valeurs estimées de C13 provenant de la figure 11E. Les valeurs attendues basées sur les données pour C13, 6 et les rapports de Poisson VTI sont illustrées comme des barres verticales 1190e à 1190i. Comme on peut le voir dans la figure 11E à 11I, des résultats de calcul plus robustes pour toutes les propriétés sont observés. Même s'il n'est pas illustré dans les figures 11A à 111, le même comportement (résultats de calcul plus robustes à l'aide de l'équation 31) a été observé pour les autres données provenant de différents échantillons de roche. [0085] Afin d'étudier davantage la performance de l'exemple de procédé de calcul de la moyenne précédent («procédé de calcul de la moyenne» dans le Tableau 1) il a été comparé à un certain nombre d'autres procédés. Les données provenant du même échantillon de roche des figures 5A à 5! ont été utilisées comme référence. C13 on suppose que 5, et les rapports de Poisson dynamiques mesurés à partir de ces données sont des valeurs justes. [0086] L'un des procédés (« l' procédé alt. » dans le Tableau 1) qui a été utilisé pour la comparaison, provient de la théorie de glissement linéaire C13 peut être calculé en solutionnant l'équation (32). C11C33-C132 = 2C66(C33+C13) (32) [0087] Le procédé suivant (« 26' procédé alt. » dans le Tableau 1) est un procédé empirique récent pour estimer C13. En outre, le soi-disant « Modified ANNIE » (« 3è1e procédé alt. » dans le Tableau 1) a été utilisé pour la comparaison. Après le calcul de C13 à l'aide de ces procédés, Ô, et les rapports de Poisson ont été calculés. Les statistiques de performance pour C13 et 8 sont donnés dans le Tableau 1. Comme on peut le voir dans les statistiques de performance, le procédé de calcul de la moyenne démontre une meilleure performance que les autres procédés dans la prédiction de tous les paramètres. La même performance est retrouvée pour les rapports de Poisson. A C13 Coefficient de Erreur RMS Coefficient de Erreur RMS corrélation corrélation Procédé. moy. 0,9648 0,0289 0,9958 0,9668 l' procédé alt. 0,8749 0,0898 0,9916 2,9604 211'e procédé alt. 0,6780 0,0867 0,9418 3,2092 3ème procédé alt. 0,8464 0,01594 0,9306 3,6264 Tableau 1 :Coefficient de corrélation et erreurs dans la prédiction de la moyenne quadratique (RMS) C13 et de Ô. [0088] Les figures 12A à 12F illustrent un autre exemple de jeu de graphiques 1200a-f. Dans cet exemple, les données provenant des mêmes 450 échantillons de roche ont été utilisées dans les figures6A à 6D ont été corrigées à l'aide de l'équation (31) à chaque fois que le C13 mesuré s'est retrouvé à l'extérieur des limites données dans les équations (29 et 30). Par conséquent, pour ces points de données, en sus du C13, et les rapports de Poisson ont également été recalculés. Cette correction, cependant, n'a pas affecté les paramètres 21 de Thomsen c et y parce que C13 n'est pas utilisé dans leur définition. Les données pour les échantillons de roche qui ont respecté les limites dans les équations (29 et 30) n'ont pas été touchées dans cette expérience. [0089] Les résultats de cette analyse, avant la correction, sont donnés dans les graphiques 1200a-c des figures 12A à 12C, et les résultats après correction sont donnés dans les graphiques 1200d-f des figures 12D à 12F. Pour tous les tracés illustrés, les hachures illustrent la valeur des données de Thomsen s, qui n'ont pas été modifiées par la correction. Comme on peut le voir dans la figure 12B, le fait que vHv a la valeur la plus élevée parmi les rapports de Poisson VTI et vv a les valeurs les plus faibles (impliquant vv REVENDICATIONSMODIFIÉES Ce qui est revendiqué est : 1. Procédé permettant de déterminer les propriétés mécaniques d'un milieu anisotrope, le procédé comprenant : l'obtention des données de journal du milieu anisotrope, les données de journal correspondant aux mesures du milieu anisotrope recueillies avec un outil de diagraphie ; la détermination des valeurs pour une pluralité de premiers composants de rigidité d'une matrice de rigidité basée sur les vitesses horizontale et verticale dérivées des données de journal ; la détermination d'une limite supérieure pour un second composant de rigidité de la matrice de rigidité basée sur les valeurs pour la pluralité des premiers composants de rigidité ; l'estimation d'une valeur pour le second composant de rigidité basée sur la limite supérieure déterminée ; la détermination d'une propriété mécanique du milieu anisotrope basée sur la valeur estimée du second composant de rigidité ; et la fourniture de la propriété mécanique déterminée. 2. Procédé permettant de déterminer les propriétés mécaniques d'un milieu anisotrope, le procédé comprenant : l'obtention des données de mesure du milieu anisotrope ; la détermination des valeurs pour une pluralité de premiers composants de rigidité d'une matrice de rigidité basée sur les vitesses horizontale et verticale dérivées des données de mesure ; la détermination d'une valeur mesurée pour un second composant de rigidité de la matrice de rigidité basée sur une vitesse hors axe dérivée des données de mesure ; la détermination d'une limite supérieure et d'une limite inférieure pour le second composant de rigidité basée sur les valeurs pour la pluralité des premiers composants de rigidité ; la comparaison de la valeur mesurée pour le second composant de rigidité à la limite supérieure et à la limite inférieure déterminées pour le second composant de rigidité ; 37 la détermination d'une propriété mécanique du milieu anisotrope basée sur la valeur mesurée du second composant de rigidité lorsque la valeur mesurée est déterminée comme étant à l'intérieur de la limite supérieure et de la limite inférieure ; et la fourniture de la propriété mécanique déterminée. 3. Procédé selon la revendication 1, comprenant également la collecte des données de journal en émettant une pluralité d'ondes dans le milieu anisotrope et la mesure d'une réponse correspondante, dans lequel la pluralité d'ondes implique une pluralité de directions de propagation différentes et une pluralité de modes de propagation différents. 4. Procédé selon la revendication 3, comprenant également le placement de l'outil de diagraphie au fond du trou dans un puits de forage, dans lequel l'outil de diagraphie comprend un ou plusieurs émetteurs conçus pour émettre la pluralité d'ondes et un ou plusieurs récepteurs conçus pour mesurer la réponse correspondante, dans lequel l'un ou les plusieurs émetteurs sont axialement espacés de l'un ou des plusieurs récepteurs. 5. Procédé selon la revendication 4, dans lequel l'un ou les plusieurs émetteurs sont conçus pour émettre des ondes acoustiques, dans lequel l'un ou les plusieurs récepteurs sont conçus pour mesurer la réponse correspondante après interaction des ondes acoustiques émises avec le milieu anisotrope, et dans lequel les données de journal comprennent des données de forme d'onde acoustique. 6. Procédé selon la revendication 2, dans lequel les données de mesure comprennent des données de carottage correspondant à un échantillon du milieu anisotrope obtenues à l'aide d'un outil de carottage, et dans lequel la vitesse hors axe est dérivée des données de carottage. 7. Procédé selon la revendication 6, comprenant également la collecte des données de carottage en extrayant l'échantillon avec l'outil de carottage et en mesurant l'échantillon extrait. 8. Procédé selon la revendication 1 ou 2, dans lequel la propriété mécanique déterminée comprend un rapport de Poisson, dans lequel le procédé comprend également la détermination d'une pression de fracturation basée sur le rapport de Poisson. 38 9. Procédé selon la revendication 1 ou 2, dans lequel la limite supérieure est déterminée en se basant sur une racine d'une fonction cubique, dans lequel la pluralité des premiers composants de rigidité est incluse dans les coefficients de la fonction cubique. 10. Procédé selon la revendication 1 ou 2, dans lequel la pluralité des premiers composants de rigidité correspond aux composants C11, C33, C55 et C66 de la matrice de rigidité, et dans lequel le second composant de rigidité correspond à un composant C13 de la matrice de rigidité. 11. Procédé selon la revendication 10, dans lequel la limite supérieure est déterminée en se basant sur Clu3 = t1 - -, dans lequel : C113 représente la limite supérieure, 3a b = -2(C11 - C66), c = -C 11C33, -33, d = 2C33 (C11 - C66)(C11 - 2C66), et t1 est la seconde racine de q3 + bq3 + cC13 + d = 0. 12. Procédé selon la revendication 1 ou 2, comprenant également : l'estimation de la valeur du second composant de rigidité basée sur une moyenne de la limite supérieure et d'une limite inférieure pour le second composant de rigidité. 13. Système permettant de déterminer les propriétés mécaniques d'un milieu anisotrope, le système comprenant : un outil de diagraphie conçu pour collecter les mesures du milieu anisotrope ; et un système de traitement ayant un processeur et une mémoire, le système de traitement étant conçu pour implémenter le procédé selon l'une quelconque des revendications 1, 3 à 5 ou 8 à 12, dans lequel les données de journal proviennent de l'outil de diagraphie. 14. Système selon la revendication 13, dans lequel l'outil de diagraphie comprend : un ou plusieurs émetteurs conçus pour émettre une pluralité d'ondes dans le milieu anisotrope, la pluralité d'ondes étant conçue pour impliquer une pluralité de directions de propagation différentes et une pluralité de modes de propagation différents ; et 39 un ou plusieurs récepteurs axialement séparés de l'un ou des plusieurs émetteurs et conçus pour mesurer une réponse correspondante par rapport à la pluralité d'ondes émises. 15. Système permettant de calculer les propriétés mécaniques d'un milieu anisotrope, le système comprenant : un outil de carottage conçu pour collecter un échantillon du milieu anisotrope ; et un système de traitement ayant un processeur et une mémoire, le système de traitement étant conçu pour implémenter le procédé selon l'une quelconque des revendications 2 ou 6 à 12, dans lequel les données de mesure proviennent de l'outil de carottage.