ANALYSE DE VITESSE TOMOGRAPHIQUE PAR STRUCTURE DE TENSEUR DE CONTRAINTES RÉFÉRENCE CROISÉE AUX DEMANDES CONNEXES [1] Cette demande revendique le bénéfice de la demande américaine provisoire No. 62/068,161, intitulée « STRUCTURE TENSOR CONSTRAINED TOMOGRAPHIC VELOCITY ANALYSIS », et déposée le 24 octobre 2014, qui est incorporée ici à titre de référence à toutes fins utiles. ARRIÈRE-PLAN [2] La présente invention concerne généralement l'exploration sismique et spécifiquement l'analyse de la vitesse tomographique utilisant un tenseur de structure comme contrainte. 1003] La sismologie est utilisée dans l'exploration, les études archéologiques et des projets d'ingénierie qui nécessitent des informations géologiques. La sismologie d'exploration procure des données qui, lorsqu'elles sont utilisées en association avec d'autres données géophysiques, géologiques et de trou de forage disponibles, peuvent fournir des informations sur la structure et la distribution des types de roches et de leurs contenus. De telles informations sont d'une grande aide pour la recherche de l'eau, des réservoirs géothermiques et des dépôts de minéraux tels que des hydrocarbures et des minerais. La plupart des compagnies pétrolières comptent sur la sismologie d'exploration pour sélectionner des sites dans lesquels forer des puits de prospection de pétrole. [4] La sismologie traditionnelle utilise des ondes sismiques générées artificiellement pour cartographier des structures sous la surface. Les ondes sismiques se propagent à partir d'une source vers le bas dans la terre et se réfléchissent sur les limites entre les structures sous la surface. Des récepteurs en surface détectent et enregistrent les ondes sismiques réfléchies pour des analyses ultérieures. Même si certaines structures à grande échelle peuvent souvent être perçues à partir d'un examen direct des signaux enregistrés, les signaux enregistrés sont généralement traités à l'aide d'un modèle de vitesse sous la surface afin d'éliminer les distorsions et révéler les détails les plus fins dans l'image du dessous la surface. La qualité de l'image du dessous de la surface peut dépendre de la précision du modèle de vitesse sous la surface. [5] L'analyse de la vitesse peut comprendre l'extraction des informations sur la vitesse à partir des données sismiques. Un procédé permettant l'analyse de la vitesse comprend une technique de migration en profondeur avant empilement, qui est devenue un outil attrayant pour l'analyse de la 1 vitesse, non seulement en raison de sa sensibilité au modèle de vitesse mais également en raison de sa capacité à générer des erreurs résiduelles dans le domaine post-migration. Une approche populaire à l'analyse de la vitesse par migration (MVA) est l'analyse de la courbure résiduelle sur un regroupement de point d'image commune, qui est basé sur un étalement résiduel permettant de mesurer l'erreur de vitesse. L'analyse de la courbature résiduelle dans les zones de structure complexe est couplée à un problème de migration-inversion qui peut être analysé d'un point de vue tomographique. 1006] Les procédés de traitement tomographiques MVA existants nécessitent l'étape de sélection, comprenant (1) la sélection de l'horizon dans le volume de la profondeur de l'image pour l'estimation des informations sur l'inclinaison et l'azimut local et (2) la sélection de l'étalement résiduel dans les regroupements de l'image commune migrée en profondeur pour la mesure des informations sur la profondeur résiduelles. La sélection manuelle peut être fastidieuse et chronophage, particulièrement dans le traitement itératif et les techniques d'interprétation. BRÈVE DESCRIPTION DES FIGURES [7] Une compréhension plus complète des modes de réalisation de la présente invention et des avantages de ceux-ci, peut être acquise en se référant à la description suivante prise en association avec les figures ci-jointes, dans lesquelles les numéros de référence semblables indiquent les caractéristiques semblables. [8] La Figure 1 est un diagramme illustrant une vue latérale d'un environnement d'un relevé sismique marin illustratif, selon des aspects de la présente divulgation. [9] La Figure 2 est un diagramme illustrant une vue du haut d'un environnement de relevé sismique marin illustratif, selon des aspects de la présente divulgation. [0010] La Figure 3 est un diagramme illustrant un schéma de point moyen illustratif qui provient des clichés flip-flop reçus par une chaine donnée, selon des aspects de la présente divulgation. 10011] La Figure 4 est un diagramme illustrant un système d'enregistrement de relevé sismique illustratif, selon des aspects de la présente divulgation. [0012] La Figure 5 est un diagramme illustrant un jeu de traces illustratif, selon des aspects de la présente divulgation. [0013] La Figure 6 est un diagramme illustrant un volume de données illustratif en 3D, selon des aspects de la présente divulgation. [0014] La Figure 7 est un diagramme illustrant un cliché de géométrie illustratif, selon des aspects de la présente divulgation. [0015] La Figure 8 est un organigramme illustrant un procédé de MVA tomographique, selon des aspects de la présente divulgation. [0016] La Figure 9 est un diagramme illustrant des vecteurs propres dans un tenseur de structure, selon des aspects de la présente divulgation. [0017] La Figure 10a-e sont des diagrammes illustrant l'obtention des informations sur la structure à partir d'une image migrée, selon des aspects de la présente divulgation. [0018] La Figure 11 sont des diagrammes illustrant le résultat des résiduels de profondeur automatiquement sélectionnés, selon des aspects de la présente divulgation. [0019] Les Figures 12a-b sont des diagrammes illustrant un exemple d'un angle d'inclinaison en ligne superposé sur une image migrée en temps inversé et des trajectoires de rayon superposés sur un modèle de vitesse, respectivement, selon les aspects de la présente divulgation. [0020] Les Figures 13a-b sont des diagrammes illustrant l'inclinaison et l'azimut provenant du tenseur de structure d'un ensemble de données SEAM, selon des aspects de la présente divulgation. [0021] Les Figures 14a-b sont des diagrammes illustrant une comparaison d'une couverture de densité de rayon dans une tranche de profondeur sur la surface provenant d'un ensemble de données SEAM, selon des aspects de la présente divulgation. [0022] Les Figures 15a-j sont des diagrammes illustrant des comparaisons de vitesse actualisée, d'image et de regroupement provenant d'un exemple d'ensemble de données SEAM, selon des aspects de la présente divulgation. [0023] La Figure 16 est un diagramme illustrant un système d'imagerie illustratif, selon des aspects de la présente divulgation. [0024] Alors que des modes de réalisation de cette divulgation ont été illustrés et décrits et sont définis par référence à des exemples de modes de réalisation de la divulgation, de telles références n'impliquent pas une limite sur la divulgation, et aucune limite de la sorte ne doit être déduite. L'objet de l'invention divulgué est capable de modifications, altérations et d'équivalents considérables dans la forme et dans la fonction, comme il sera apparent aux spécialistes du domaine pertinent et qui bénéficient de cette divulgation. Les modes de réalisation illustrés et décrits de cette divulgation ne sont que des exemples, et ne sont pas une description exhaustive de la divulgation.
3 DESCRIPTION DÉTAILLÉE [0025] Les modes de réalisation illustratifs de la présente invention sont décrits en détail ci-dessous. Dans un souci de clarté, toutes les caractéristiques d'un mode de réalisation réel ne sont pas décrites dans cette spécification. Il sera, bien sûr, apprécié que dans le développement d'un quelconque mode de réalisation réel, que de nombreuses décisions spécifiques à une concrétisation doivent être prises afin d'atteindre les objectifs spécifiques des développeurs, tels que la conformité avec des contraintes apparentées au système ou aux considérations monétaires, qui varieront d'une concrétisation à une autre. En outre, il sera apprécié qu'un tel effort de développement puisse être complexe et chronophage, mais serait néanmoins une entreprise de routine pour les hommes de métier qui bénéficient de la présente divulgation. [0026] Pour faciliter une meilleure compréhension de la présente divulgation, les exemples suivants de certains modes de réalisation sont donnés. En aucun cas, les exemples suivants ne doivent être interprétés comme limitant, ou définissant, la portée de l'invention. Les modes de réalisation de la présente divulgation peuvent être applicables à des puits de forage horizontaux, verticaux, déviés ou autrement non-linéaires dans un quelconque type de formation souterraine. Des modes de réalisation peuvent être applicables à des puits d'injection aussi bien que des puits de production, y compris des puits d'hydrocarbures. Les modes de réalisation peuvent être implémentés avec un outil qui est approprié pour le test, la récupération et l'échantillonnage le long des sections de la formation. [0027] Certains ou tous les aspects de la présente divulgation peuvent être implémentés dans un système de traitement d'informations ou un système informatique, les deux pouvant être utilisés de façon interchangeable ici. Des exemples de systèmes de traitement d'informations comprennent des systèmes de serveurs, des terminaux informatiques, des dispositifs informatiques de poche, des tablettes, des téléphones intelligents, etc. Pour les besoins de cette divulgation, un système de traitement d'informations peut comprendre une instrumentalité ou un agrégat d'instrumentalités fonctionnant pour calculer, classifier, traiter, transmettre, recevoir, récupérer, être à l'origine de, commuter, stocker, afficher, manifester, détecter, enregistrer, reproduire, traiter ou utiliser une quelconque forme d'informations, d'intelligence ou des données commerciales, scientifiques, de contrôle, ou pour d'autres objectifs. Par exemple, un système de traitement d'informations peut être un ordinateur personnel, un dispositif de stockage en réseau ou tout autre dispositifs appropriés et peut varier en taille, en forme, en performance, en fonctionnalité et en prix. Le système de 4 traitement d'informations peut comprendre une mémoire RAM, une ou plusieurs ressources de traitement telles qu'une unité de traitement central (CPU) ou une logique de commande de matériel ou de logiciel, une ROM et/ou d'autres types de mémoire non-volatile. Des composants additionnels du système de traitement d'informations peuvent comprendre une ou plusieurs disques, un ou plusieurs ports de réseau pour la communication avec des dispositifs externes aussi bien que divers dispositifs d'entrée et de sortie (E\S) tels qu'un clavier, une souris et un écran vidéo. Le système de traitement d'informations peut également comprendre un ou plusieurs bus fonctionnant pour transmettre des communications entre les divers composants du matériel. [0028] Pour les besoins de cette divulgation, un support lisible par ordinateur peut comprendre toute instrumentalité ou agrégat d'instrumentalités qui peuvent retenir des données et/ou des instructions pour une période de temps. Le support lisible par ordinateur peut comprendre, par exemple, sans limitation, un support de stockage tel qu'un dispositif de stockage à accès direct (par ex., un disque dur ou une disquette), un dispositif de stockage à accès séquentiel (par ex., une cassette), un disque compact, CD-ROM, DVD, RAM, ROM, une mémoire EEPROM ou une mémoire flash ; aussi bien que des supports de communication tels que des câbles, des fibres optiques, des micro-ondes, des ondes radio et d'autres supports électromagnétiques et/ou optiques ; et/ou une quelconque combinaison des éléments précédents. [0029] Telle qu'elle est décrite ici, une approche d'analyse automatique de la vitesse par migration tomographique, basée sur une grille, en 3D (MVA) peut utiliser un tenseur de structure comme contrainte. Des exemples d'approches pourraient ne pas nécessiter la sélection manuelle des horizons géologiques dans les volumes de l'image de profondeur et les étalements résiduels dans les regroupements d'image migrée commune. Le tenseur de structure peut être utile pour estimer les informations de l'inclinaison locale et de l'azimut, qui peuvent être utilisées comme contraintes pour le calcul des dérivées de Fréchet au cours de l'inversion tomographique. [0030] Les aspects de la présente divulgation peuvent être compris dans un contexte illustratif tel que le relevé sismique marin, comme celui illustré dans les figures 1-5, même si cette divulgation n'est pas limitée aux relevés marins. En mer, les navires de relevé sismique peuvent déployer une flûte derrière le navire tel que le démontre la Fig. 1. Chaque flûte 110 peut traîner derrière le navire 100 lorsque le navire se déplace vers l'avant (dans la direction de la flèche 102), et chaque flûte comprennent de multiples récepteurs également espacés 114. Chaque flûte 110 peut en outre comprendre un aiguilleur 118 programmable et des commandes de la profondeur programmables5 qui tirent la flûte à une distance de décalage de fonctionnement à partir de la trajectoire du navire (voir la Fig. 2) et jusqu'à une profondeur de fonctionnement souhaitée (Fig. 1). [0031] Les flûtes 110 peuvent faire plusieurs kilomètres de long, et sont généralement construites en sections d'une longueur de 25 à 100 mètres qui comprennent des groupes de 35 ou plus de récepteurs uniformément espacés. Chaque flûte 110 peut comprendre des câbles électriques ou en fibre optique pour interconnecter les récepteurs 114 et des équipements sismiques sur un navire 100. Les données peuvent être numérisées proches des récepteurs 114 et transmises vers le navire 110 à travers le câblage à une vitesse de 7 (ou plus) millions de bits de données par seconde. [0032] Tel que le démontre la Fig. 1, le navire de relevé sismique 100 peut également tirer une ou plusieurs sources 112. La source 112 peut être une source d'impulsion ou une source de vibration. Les récepteurs 114 utilisés dans la sismologie marine sont communément appelés des hydrophones, et sont généralement construits avec un transducteur piézoélectrique. Divers types d'hydrophones sont disponibles tels que des hydrophones à disque ou des hydrophones cylindriques. Les sources 112 et les récepteurs 114 sont généralement déployés en-dessous de la surface de l'océan 104. L'équipement de traitement, tel que des systèmes de traitement d'informations, à bord des navires de commande du fonctionnement des sources et des récepteurs et enregistre les données acquises. [0033] Les relevés sismiques peuvent procurer des données pour l'imagerie en-dessous de la surface de l'océan 104 pour révéler des structures sous la surface telles que la structure 106, qui se trouve en-dessous du plancher océanique 108. Les analystes utilisent des méthodes d'imagerie pour traiter les données et cartographier la topographie des couches sous la surface. Les données du relevé sismique révèlent diverses autres caractéristiques des couches sous la surface qui peuvent être utilisées pour déterminer les emplacements des réservoirs de pétrole et/ou de gaz. [0034] Afin de faire l'imagerie de la structure 106 sous la surface, la source 112 peut émettre des ondes sismiques 116 qui sont réfléchies lorsqu'il y a des changements dans l'impédance acoustique causées par une structure sous la surface 106 (et d'autres réflecteurs sous la surface). Les ondes réfléchies sont détectées par un schéma de récepteurs 114. En enregistrant (comme une fonction du temps) les ondes sismiques 116 entrant qui voyagent à partir d'une source 112 vers la structure 106 sous la surface 114, une image de la structure 106 sous la surface peut être obtenue après traitement approprié des données. [0035] La Fig. 2, illustre une vue du dessus (pas à l'échelle) du navire de relevé sismique 100 qui tire un jeu de flûtes 110 et deux sources 112. Lorsque le navire 100 avance, les sources 112 peuvent 6 être déclenchées alternativement dans un schéma appelé « flip-flop ». Les aiguilleurs programmables sont utilisés pour permettre un espacement plus ou moins égal entre les flûtes. Les récepteurs se trouvant à une position donnée sur les flûtes sont associés à un numéro de fichier de trace de terrain commun ou à un canal commun 202. [0036] La Fig. 3, illustre une vue du dessus des positions illustratives de la source et de récepteur pour deux clichés. Pour un premier cliché, une source est déclenchée à la position 302, et la partie illustrée du réseau de récepteurs est à la position 304 (illustrée en ligne pointillée). Pour un deuxième cliché, une source est déclenchée à la position 306 et la partie illustrée du réseau de récepteur est à la position 308 (illustrée en ligne solide). En assumant, pour le moment, que les structures réfléchissantes sous la surface sont horizontales, les ondes sismiques qui atteignent chacun des 12 récepteurs sont réfléchies à partir d'une position en-dessous du point à égale distance entre les positions de la source et du récepteur. Ainsi, le premier cliché produit des réflexions à partir du dessous des 12 points à mi-parcours 311 (illustrés en ligne pointillée avec des hachures verticales), alors que le deuxième cliché produit des réflexions du dessous des 12 points à mi-parcours 310 (illustrés en ligne solide avec une hachure horizontale). Comme un exemple, le vecteur 312 illustre une propagation de l'énergie sismique provenant d'un cliché 302 vers un point à mi-parcours 314, et un vecteur 316 de longueur égale illustre la propagation de l'énergie sismique vers une position du récepteur. Pour le deuxième cliché 306, les vecteurs 318 et 320 illustrent une trajectoire de propagation similaire. Il est à noter que le point à mi-parcours 314 est l'un des points à mi-parcours qui est « bombardé » de multiples clichés, rendant ainsi disponible plus d'énergie de signal à partir de ces zones lorsque les informations provenant des clichés sont traitées et associées. Les relevés sismiques (pour la terre et la mer) sont généralement conçus pour procurer une grille également distribués de points à mi-parcours avec un nombre moyen de coups relativement élevé pour chaque point à mi-parcours. [0037] La Fig. 4 montre un système d'enregistrement du relevé sismique illustratif ayant des récepteurs 114 couplés à un bus 402 pour communiquer des signaux numériques au circuit d'enregistrement des données 406 sur un navire de sondage 100. Des capteurs d'informations sur la position et d'autres capteurs de paramètre 404 sont couplés au circuit d'enregistrement des données 406 pour permettre au circuit d'enregistrement des données de stocker des informations additionnelles utiles pour interpréter les données enregistrées. De façon illustrative, de telles 7 informations supplémentaires peuvent inclure des informations sur l'orientation du réseau et des informations sur la vitesse. [0038] Un système de traitement de données numérique polyvalent 408, qui peut comprendre un système de traitement d'informations, est illustré comme étant couplé à un circuit d'enregistrement de données 406, et il est également illustré comme étant couplé à travers un bus 402 à des dispositifs de positionnement 410 et à des sources sismiques 112. Le système de traitement 408 configure le fonctionnement du circuit d'enregistrement 406, des dispositifs de positionnement 410 et des sources sismiques 112. Le circuit d'enregistrement 406 peut acquérir le ou les flux de données à vitesse élevée provenant des récepteurs 114 sur un support de stockage non-volatile tel qu'un réseau de stockage de disques magnétiques ou optiques. Les dispositifs de positionnement 410 (y compris des aiguilleurs programmables et des commandes de la profondeur) peuvent contrôler la position des récepteurs 114 et des sources 112. [0039] Le système d'enregistrement sismique Fig. 4 peut comprendre des composants additionnels qui ne sont pas spécifiquement illustrés ici. Par ex., chaque flûte 110 pourrait avoir un bus indépendant 402 pour le couplage avec le circuit d'enregistrement des données. Le système de traitement 408 peut comprendre une interface utilisateur possédant un écran graphique et un clavier ou un autre procédé d'acceptation de saisie par l'utilisateur, et peut également comprendre une interface de réseau permettant de communiquer les données des relevés sismiques stockées dans une centrale informatique possédant des ressources informatiques puissantes permettant de traiter les données du relevé sismique. [0040] La Fig. 5 montre des signaux sismiques illustratifs, qui peuvent être appelés des traces, détectés et échantillonnés par des récepteurs 114. Les signaux indiquent une certaine mesure de l'énergie de l'onde sismique en fonction du temps (par ex., déplacement, vélocité, accélération, pression), et ils sont numérisés à haute résolution (par ex., 24 bits) avec, une vitesse d'échantillonnage programmable. De tels signaux peuvent être regroupés de différentes façons, et lorsqu'ils sont regroupés, ils sont appelés un « regroupement ». Par ex., un « regroupement de points à mi-parcours communs » est un groupe de traces qui possèdent un point à mi-parcours dans une région définie. Un « regroupement de clichés » est un groupe de traces enregistré lors d'un déclenchement unique de la source sismique. Un « regroupement multi-clichés » est un regroupement de clichés, souvent comprenant toutes les traces enregistrées le long d'une trajectoire de navigation lors d'un relevé sismique marin. 8 [0041] Même s'il est possible de rapporter les différentes formes d'onde enregistrées dans le format illustré dans la figure 5, côte à côte avec un graphique révélant des structures à grande échelle sous la surface, de telles structures sont déformées et n'illustrent pas les structures plus fines. Dans certains modes de réalisation, les formes d'onde brutes illustrées dans la Fig. 5, peuvent être traitées pour créer un volume de l'image de la profondeur, c.-à-d., un réseau en 3D des valeurs de données telles que celles illustrées dans la Fig. 6. Le volume de l'image de la profondeur représente certains attributs sismiques à travers diverses profondeurs et orientations spatiales à l'intérieur d'une région sondée. Le réseau en 3D comprend des cellules de taille uniforme, chaque cellule ayant une valeur de données représentant l'attribut sismique pour cette cellule. Divers attributs sismiques peuvent être représentés, et dans certains modes de réalisation, chaque cellule possède de multiples valeurs de données pour représenter de multiples attributs sismiques. Des exemples d'attributs sismiques appropriés comprennent la réflectivité, l'impédance acoustique, la vitesse acoustique et la densité. Le format de données volumétriques se prête plus volontairement à une analyse informatique et à un rendement visuel, et pour cette raison, le volume de l'image de la profondeur peut être appelé une « image en 3D » de la région sondée. [0042] La Fig. 7, montre la façon dont divers paramètres sont apparentés à la géométrie d'un cliché illustratif en 2D (le cas en 3D est semblable). L'énergie sismique se propage le long d'un rayon 702 à partir d'une source sismique vers une interface cible 704 et se réfléchit vers un récepteur le long d'un rayon 706. Au point de réflexion (représenté ailleurs par des coordonnées (x,y,z) et abrégé ici comme un vecteur {flèche droite au-dessus de (x)}), la surface 704 a un vecteur normal {flèche droite au-dessus de (n)} à un angle a par rapport à la verticale. Le rayon entrant 702 et le rayon réfléchi 706 sont à des angles 0 égaux (mais opposés) « ouverts » par rapport au vecteur normal. [0043] Les données de la trace sismique initialement recueillies lors d'un relevé peuvent être acquises en fonction de l'emplacement du cliché, de l'emplacement du récepteur et du temps, c.-àd., P(s,r,t). Traditionnellement, un changement de variable est réalisé afin de placer ces données dans le domaine du point à mi-parcours/décalage/temps, c.-à-d., P(m,h,t), où le point à mi-parcours m----(s+r)12 et décalage h=ls-r112. En observant que ces données représentent les champs d'onde observés au niveau de la surface (z = 0), l'équation du champ d'onde est utilisée pour extrapoler le champ d'onde sous la surface, un procédé connu sous le nom de migration. [0044] Un exemple de technique de migration comprend les équations suivantes : P(m,h,t;z = 0)P(m,h, w;z = 0) (1) 9 P(m,h,w;z = 0)->P(m,h,w;z) (2) P(m,h,w;z)--313(m,ph,r;z) (3) P(rn,p h,r;z)-)P(m,p h,r = 0;z) (4) [0045] L'équation (1) représente une transformée de Fourier de l'ensemble des données afin de placer les données acquises au niveau de la surface (z=0) dans le domaine de point à miparcours/décalage/fréquence. L'équation (2) représente la migration de l'ensemble de données utilisant une racine carrée double (DSR) bien connue pour extrapoler un champ d'ende. L'équation (3) représente une transformée de Radon, qui peut également être appelée un fonctionnement de pile inclinée, des données dans le domaine de point à mi-parcours-p-tau. Les paramètres de rayon décalé p et tau peuvent représenter la pente et l'intersection des lignes inclinées utilisées pour empiler les données. [0046] Comme l'indique l'équation (4), le définition du tau égal à zéro procure un jeu de regroupements d'image commune de domaine d'angle, qui peuvent être visualisés comme un jeu d'images P(m,z), chaque image étant dérivée de l'énergie sismique frappant le réflecteur à un angle différent. Le paramètre de rayon décalé ph est apparenté à l'inclinaison local a et à l'angle ouvert Es par l'équation : ph= 2 * S(rn,z) * cos a * sin 0 (5) où S(m,z) représente la lenteur (l'inverse de la vitesse acoustique V(m,z)) dans le voisinage du réflecteur. [0047] La MVA tomographique peut être utilisée pour déterminer et/ou raffiner un modèle de vitesse basé sur les mésappariements de profondeur dans les regroupements d'image commune. Dans le domaine de l'angle post-migré, les données sismiques P(m,ph,z) représentent les positions de profondeur de multiples images de l'emplacement du réflecteur. L'utilisation d'un bon modèle de vitesse-profondeur V(m,z) dans la migration génère des regroupements CA1 plats dans le domaine ph-z (c'est à dire, les réflecteurs apparaissent comme des événements à une profondeur constante z, quel que soit ph). Sinon, les résiduels de profondeur sont présents sur les regroupements CAI, ce qui veut dire que la profondeur de l'événement varie avec ph). Afin d'adapter le domaine des résiduels de profondeurs ph-z à l'approche MVA tomographique, ils sont convertis en perturbations de temps de trajet At(ph), qui reflètent l'étalement résiduel d'une trajectoire de rayon spéculaire. Ayant choisi une profondeur de référence, les résiduels de profondeur Az provenant de la profondeur de référence au niveau d'un emplacement du réflecteur peut être déterminé en utilisant 10 une semblance de calcul entre des images communes calculées à différents angles dans les regroupements d'image commune de domaine d'angle. La conversion de résiduel de profondeur en perturbation de temps de trajet dans le domaine ph-z peut être exprimé sous forme de At(ph) = Az,i4S2cos2 a - ph2 (6) où S représente la lenteur locale au-dessus de la perturbation du réflecteur et a représente l'angle d'inclinaison local du réflecteur. L'équation (6) peut calculer la perturbation du temps de trajet causée par la longueur supplémentaire du trajet qu'un rayon doit emprunter à cause de la déviation de la profondeur. La dépendance de la perturbation du temps de trajet sur l'angle d'inclinaison du réflecteur est petite pour des petites inclinaisons mais devient importante pour des inclinaisons plus grandes. Il est à noter que si l'angle incident 0 est souhaité il peut être obtenu sans traçage de rayon en utilisant l'équation (5). Par conséquent, les perturbations du temps de trajet calculées à partir des regroupements CAI sont insensibles aux erreurs du trajet du rayon, permettant l'utilisation d'un algorithme plus rapide de traçage de rayons. [0048] Les procédés de traitement tomographique MVA existants, y compris le procédé de traitement tomographique MVA, peuvent nécessiter une ou plusieurs étapes dans lesquels les valeurs sont manuellement choisies par un ingénieur ou un technicien. L'étape peut comprendre, par ex., la sélection d'un horizon dans le volume de l'image de la profondeur pour l'estimation des informations de l'inclinaison locale et de l'azimut et de la sélection de l'étalement résiduel dans les regroupements de l'image commune migrée pour la profondeur pour la mesure des informations de la profondeur résiduelle. La sélection manuelle peut être fastidieuse et chronophage, particulièrement dans le traitement itératif et les techniques d'interprétation dans lesquelles les valeurs sont choisies à chaque itération. [0049] Selon des aspects de la présente divulgation, une approche automatique de migration tomographique MVA peut utiliser un tenseur de structure comme contrainte de sorte que la sélection manuelle des horizons géologiques dans les volumes d'image de profondeur et les étalements résiduels dans les recrutements de l'image commune migrée ne soient pas nécessaires. Le tenseur de structure, par ex., peut être utile pour estimer les informations de l'inclinaison locale et de l'azimut, qui peuvent être utilisées comme contraintes pour le calcul des dérivées de Fréchet au cours de l'inversion tomographique. [0050] La Figure 8 est un organigramme illustrant un procédé de MVA tomographique 800, selon des aspects de la présente divulgation. Le procédé peut commencer à l'étape 805 dans laquelle les 11 traces sismographiques sous la forme de données pré-empilées sont collectées. Les traces sismographiques peuvent être collectées, par ex., en utilisant un système de relevé sismique semblable à ou différent de celui décrit en référence à la Fig. 1. La collecte des traces sismographiques peut également comprendre la réception des informations au niveau d'un système de traitement ou d'un processeur d'un système de traitement d'informations, les traces sismographiques précédemment collectées provenant d'un support sur lequel les traces ont été précédemment stockées. Ceci peut comprendre, par ex., un dispositif de mémoire couplé à un processeur, ou à un serveur dans un centre de dépôt de données central. Les traces précédemment enregistrées peuvent être reçues, par ex., à travers un ou plusieurs canaux de communication câblés ou non. [0051] À l'étape 810, une migration en profondeur peut être réalisée basée, au moins en partie, sur le modèle de vitesse 860. La migration peut, par ex., prendre la forme de l'exemple de la technique de migration décrite ci-dessus, mais cette technique de migration n'est pas destinée à être limitée, et peut comprendre d'autres techniques de migration qui seraient comprise par un homme de métier à la lumière de cette divulgation. Dans la première itération de l'étape 810, le modèle de vitesse 860 peut comprendre un modèle de vitesse initiale. Dans les itérations ultérieures de l'étape 810, le modèle de vitesse 860 peut comprendre un modèle de vitesse actualisé provenant de l'étape 855 (présenté ci-dessous). [0052] La migration en profondeur à l'étape 810 peut être utilisée pour déterminer un volume de l'image en profondeur 815 et des regroupements d'image commune 830. Le volume de l'image en profondeur 815 peut être le résultat de la migration en profondeur de l'étape 810 et peut, mais pas nécessairement, prendre une forme semblable au volume de l'image en profondeur décrite ci-dessus en référence à la Fig. 6. Les regroupements d'image commune 830 peuvent, mais pas nécessairement, comprendre des regroupements d'image commune de domaine d'angle déterminé en utilisant le procédé décrit ci-dessus. D'autres types de regroupements d'image commune sont possibles, comme il sera compris par un homme de métier à la lumière de cette divulgation. [0053] À l'étape 820, un tenseur de structure peut être calculé à partir d'un volume d'image en profondeur 815 afin d'estimer des informations d'inclinaison structurelle et d'azimut du volume de l'image en profondeur 815. Dans certains modes de réalisation, le calcul du tenseur de structure peut, mais pas nécessairement, calculer des dérivées Gaussiennes lissées partout dans le volume de 12 l'image en profondeur 815, et ensuite la décomposition propre peut être déterminée en utilisant l'équation suivante s = .1.uuT + .1'vvr + ÀwwwT (7) [0054] où, ° 5- ÂW 212 4. La Figure 9 est un diagramme illustrant des vecteurs propres dans un tenseur de structure, selon des aspects de la présente divulgation. Comme le démontre la Fig. 9, les vecteurs propres peuvent définir un système de coordonnée binomial tangentiel au gradient d'image. Le calcul du tenseur de structure peut donner 12 des volumes de vecteurs propres (3 composants cartésiens chacun) et des valeurs propres. Les volumes d'attributs d'inclinaison peuvent être calculés à partir de la normale de la tangente U; la grandeur de l'inclinaison peut être une somme Euclidienne de 2 composants cartésiens latéraux ; et l'azimut peut être l'arc-tangente de ces composants. [0055] En revenant à la Fig. 8, à l'étape 835, la semblance peut être calculée pour les regroupements d'image commune 830, avec un tenseur de structure de l'étape 120 et les attributs de l'inclinaison calculés utilisés comme contrainte. Dans certains modes de réalisation, les composants du tenseur de structure peuvent également être utilisés comme des régions de masquage de signal faible ou d'inclinaison conflictuelle. Même si un certain nombre de semblance orientée structure et de attributs de planarité peuvent être utilisés pour les rapports des valeurs propres, en pratique, un masque de seuil 4 peut être adéquat. Pour l'illustrer, les Figs. 10a-e montrent un procédé permettant d'obtenir des informations de structure à partir d'une image migrée, selon des aspects de la présente divulgation. L'extraction de l'inclinaison peut être calculée pour un synclinal synthétique simple 1000. Comme illustré, la Fig. 10a peut comprendre une coupe à travers le milieu d'une image de profondeur du synclinal 1000. Les tenseurs de structure peuvent être construits à partir d'une image migrée de temps inversé de la Fig. 10a. La Fig. 10b illustre 4 l'attribut, la valeur propre dans la direction tangente-normale, qui peut être l'un des composants du tenseur de structure apparentée à l'amplitude et à la cohérence. La Fig. 10c illustre un profil de masque à seuil 4, qui peut être créé en utilisant la Fig. 10b. Le profil de masque de la Fig. 10c peut être créé pour éliminer un bruit de faible amplitude, non-cohérent. La Fig. 10d illustre une ligne de composant Cartésien de Uenhgne, et la Fig. 10e illustre le Uenligne masqué. Les Fig. 10d et 10e illustrent ainsi respectivement des angles d'inclinaison en ligne calculée à partir des tensions de structure avant et après l'application des profils de masques de la FIG. 10c. 13 [0056] En revenant à la Fig. 8, à l'étape 840, les résiduels en profondeur peuvent être automatiquement choisis sur les semblances résiduelles provenant de l'étape 835. Dans certains modes de réalisation, un algorithme de sélection automatique peut être utilisé qui sélectionne des fonctions permettant de maximiser la sommation à travers les valeurs de semblance dans de multiples directions simultanément basées sur une fonction de guide d'entrée et une fourchette de recherche positive et négative par rapport au guide d'entrée. Pour les étalements résiduels après la , migration, le guide d'entrée peut être un étalement zéro. L'étalement sélectionné peut être contraint pour qu'il soit lisse dans les directions verticales et toutes les directions spatiales à fin d'éviter les sélections bruitées. La sélection automatique simultanée de tous les points ensemble au moyen d'une optimisation globale peut également éviter les sélections sauvages. La Figure 11 illustre le résultat de résiduels de profondeur automatiquement sélectionnés, selon des aspects de la présente divulgation. Comme le démontre la Fig. 11, le résultat de sélection peut être un « hyper plan » de vitesse (parmi lequel, la Fig. 11 illustre la projection d'une ligne unique). [0057] De cette façon, le calcul peut être totalement automatisé et élimine le besoin d'une sélection manuelle de l'inclinaison et de l'horizon à un quelconque stade. Un contrôle peut être exercé sur ce calcul en ajustant des paramètres tels que le lissage dans le calcul du tenseur et du type de masquage utilisé. [0058] À l'étape 825, les calculs de traçage de rayon et/ou de sensibilité du noyau peuvent être effectués en utilisant le tenseur de structure de l'étape 820 comme contrainte. Dans certains modes de réalisation, les trajets de rayons peuvent être calculés en utilisant un algorithme de traçage de rayons dynamiques. Étant donné les angles de réflexion locale et de l'azimut, une paire de rayons incidente/réfléchie peut être prise en commençant au niveau du point de réflexion. Lorsque les deux rayons atteignent la surface, les emplacements de la source et du récepteur peuvent être déterminés, et le décalage source-récepteur et l'azimut de prise à la surface peuvent être obtenus. Un tel traçage de rayons permet un traitement direct des regroupements d'angle migré. Afin de traiter des regroupements de décalage migrés, pour chaque point de réflexion, l'angle de réflexion appropriée et l'angle azimut local peuvent être calculés pour correspondre au décalage en surface attendu et de l'azimut de prise jusqu'à ce que toutes les non correspondances soient minimisées à l'intérieur des tolérances données. Selon les aspects de la présente divulgation, la Fig. 12a illustre un exemple d'angle d'inclinaison en ligne superposé sur une image migrée en temps inversé, et la Fig. 12b illustre des exemples de trajets de rayons superposés sur un modèle de vitesse. 14 [00591 En revenant à la Fig. 8, à l'étape 845, une inversion tomographique peut être réalisée avec le traçage de rayons de l'étape 825 et les résiduels en profondeur automatiquement choisis de l'étape 840. Dans certains modes de réalisation, l'inversion peut être réalisée en utilisant l'équation suivante : v(r) L(h) 2 cos a[cos y / \ L(h = 0)]- As = Az(h) (8) [L(h) 1(11 = 0)] / cosi où, est un terme reflétant le traçage de rayon de l'étape 825, As est un terme reflétant la perturbation de la lenteur actualisée, et AZ ) est un terme reflétant les résiduels en profondeur automatiquement choisis de l'étape 840. Dans ce mode de réalisation, les rayons calculés ou les noyaux de sensibilité peuvent être stockés sur une matrice jacobierme clairsemée et le système d'inversion peut être résolu en utilisant le procédé conjugué-gradient. De cette façon, un modèle de vitesse actualisé peut être obtenu. [0060] Si le modèle de vitesse généré à l'étape 845 représente le modèle de vitesse fmal, il peut être la sortie du système à l'étape 865. Par ailleurs, si des itérations supplémentaires du MVA sont souhaitées pour actualiser davantage le modèle de vitesse, le modèle de vitesse 860 généré à l'étape 845 peut être fourni à l'étape 855 pour être utilisé comme le nouveau modèle de vitesse dans les itérations futures (commençant à l'étape 810 avec la migration en profondeur). [0061] Les avantages du procédé MVA tomographique décrit dans la présente divulgation peuvent être trouvés en référence à son application à un ensemble de données du « SEG Advanced Modeling Corporation (SEAM) ». Tout d'abord, une image empilée initiale peut être obtenue par migration en profondeur en utilisant un modèle de vitesse initial. Ensuite, avec le calcul des tenseurs de structure de l'image empilée initiale, l'inclinaison et l'azimut peuvent être obtenus. Dans la procédure de traçage de rayon, les informations sur l'inclinaison et l'azimut, qui associent les contraintes provenant du tenseur de structure, peuvent être saisies. En fonction du résultat du résiduel de profondeur auto-sélectionné et du traçage précis du rayon, la vitesse actualisée peut être déterminée. [0062] Les Figs. 13a-b illustrent l'inclinaison et l'azimut provenant du tenseur de structure de l'exemple de l'ensemble de données SEAM, selon des aspects de la présente divulgation. Les Figs. 13a et 13b montrent respectivement les résultats de l'inclinaison et de l'azimut provenant du tenseur 15 de structure du volume de l'image initiale. Comme le démontre les Figs. 13a-b, les résultats provenant du tenseur de structure peuvent être plus précis en comparaison avec le résultat de traçage de rayon avec une supposition d'inclinaison-azimut de zéro. [0063] Les Figs. 14a-b illustrent une comparaison d'une couverture de densité de rayon dans une tranche de profondeur sur la surface provenant d'un ensemble de données SEAM, selon des aspects de la présente divulgation. Spécifiquement, la Fig. 14a illustre la couverture de densité de rayon dans une tranche de profondeur sur la surface avec la contrainte du tenseur de structure, et la Fig. 14b illustre la couverture de densité de rayon en absence de la contrainte du tenseur de structure. La couverture de rayon améliorée avec l'utilisation de la contrainte du tenseur de structure peut être observée dans la Fig. 14a en comparaison à la Fig. 14b. [0064] Les Figures 15a-j illustrent des comparaisons de vitesse actualisée, d'image et de regroupement provenant d'un exemple d'ensemble de données SEAM, selon des aspects de la présente divulgation. Spécifiquement, la Fig. 15a illustre une vitesse initiale ; la Fig. 15b illustre une vitesse actualisée après une première itération du MVA tomographique en l'absence de la contrainte du tenseur de structure ; la Fig. 15c illustre une vitesse actualisée après une première itération du MVA tomographique en présence de la contrainte du tenseur de structure ; et la Fig. 15d illustre la vraie vitesse de l'exemple d'ensemble de données SEAM. Tel qu'illustré, la vitesse actualisée avec la contrainte du tenseur de structure (Fig. 15c) est plus proche de la vraie vitesse (Fig. 15d) après la première itération en comparaison avec la vélocité actualisée sans contrainte (Fig. 15b). [0065] La précision améliorée est également évidente dans la comparaison parmi l'image empilée (Fig. 15e-h). Spécifiquement, la Fig. 15e illustre une image de la vitesse initiale par défaut ; la Fig. 15f montre une image utilisant une vitesse actualisée après une première itération en l'absence de la contrainte du tenseur de structure ; la Fig. 15g montre une image utilisant une vitesse actualisée après une première itération en présence de la contrainte du tenseur de structure ; et la Fig. 15h montre une image de la vraie vitesse de l'exemple d'ensemble de données SEAM. Tel qu'illustré, l'image de la vitesse actualisée avec la contrainte du tenseur de structure (Fig. 15g) est plus proche de la vraie vitesse (Fig. 15h) après la première itération en comparaison avec la vélocité actualisée sans contrainte (Fig. 15f). [0066] L'amélioration de la précision est d'autant plus évidente dans deux exemples d'ensemble de données de regroupement d'image commune (Fig. 15i et Fig. 15j). La Fig. 15i illustre 4 planches 16 de regroupement d'image décalée au niveau de l'emplacement A dans la Fig. 15e. De gauche à droite, les 4 planches dans la Fig. 15i illustrent un regroupement d'images décalées avec une vitesse initiale par défaut, une vitesse actualisée en l'absence d'une contrainte de tenseur de structure, une vitesse actualisée en présence d'une contrainte de tenseur de structure, et la vraie vitesse de l'ensemble de données. La Fig. 15j illustre 4 planches similaires pour des regroupements d'image décalée au niveau d'un emplacement B dans la Fig. 15e. Comme il est illustré, dans les deux exemples de jeux de données de regroupement d'image, la vitesse actualisée avec une contrainte de tenseur de structure est plus proche de la vraie vitesse que de la vitesse actualisée en absence d'une contrainte de tenseur de structure. [0067] Certaines ou toutes les étapes du procédé illustratif décrit ci-dessus par rapport à la Fig. 8 peuvent comprendre des étapes logicielles réalisées dans un système de traitement d'informations. Un logiciel peut être caractérisé par un jeu d'instructions stocké sur un support lisible par ordinateur, qui, lorsqu'il est exécuté par un processeur, permet au processeur d'exécuter certaines fonctions. La Fig. 16 montre un système informatique illustratif 900 dans lequel le procédé illustratif peut être réalisé. Comme il est illustré, un poste de travail personnel 902 est couplé à travers un réseau local (LAN) 904 à un ou plusieurs ordinateurs multiprocesseurs 906, qui sont, à leur tour, couplés via le LAN à un ou plusieurs unités de stockage partagées 908. Le poste de travail 902 et les ordinateurs 906 peuvent comprendre des systèmes de traitement d'informations. Le poste de travail personnel 902 constitue une interface utilisateur pour le système de traitement, permettant à un utilisateur de télécharger des données du serveur vers le système, de récupérer et de visualiser des données d'images à partir du système et de configurer et de surveiller le fonctionnement d'un système de traitement. Le poste de travail personnel 902 peut prendre la forme d'un ordinateur de bureau avec un écran graphique qui montre graphiquement les données du relevé et des images en 3D de la région sondée, et avec un clavier qui permet à l'utilisateur de déplacer des fichiers et d'exécuter des logiciels de traitement. [0068] Le LAN 904 permettent une communication à haute vitesse entre les ordinateurs multiprocesseurs 906 et les postes de travail personnels 902. Le LAN 904 peut prendre la forme d'un réseau Ethernet. [0069] Le ou les ordinateurs à multiprocesseurs 906 permettent une capacité de traitement parallèle permettant une conversion relativement rapide des signaux de traces sismiques dans une image de région sondée. Chaque ordinateur 906 comprend de multiples processeurs 912, une 17 mémoire distribuée 914, un bus interne 916, est une interface LAN 920. Chaque processeur 912 fonctionne sur une partie allouée des données d'entrée pour produire une image partielle de la région sismique sondée. Un module de mémoire distribué 914 est associé à chaque processeur 912 qui stocke le logiciel de conversion et un jeu de données de travail pour l'utilisation par le processeur. Un bus interne 916 permet une communication inter-processeur et une communication avec les réseaux LAN à travers l'interface 920. La communication entre les processeurs dans les différents ordinateurs 906 peut être assurée par le LAN 904. [0070] Les unités de stockage partagées 908 peuvent être de grandes unités de stockage d'informations autonomes qui utilisent un support de disque magnétique pour un stockage non volatile des données. Afin d'améliorer la vitesse d'accès et la fiabilité, les unités de stockage partagées 908 peuvent être configurées sous forme d'un réseau de disque redondant. Les unités de stockage partagées 908 stockent un volume initial de données sur la vitesse et des regroupements de clichés provenant d'un relevé sismique. En réponse à une requête provenant d'un poste de travail 902, les données du volume d'images peuvent être récupérées par des ordinateurs 906 et alimentées vers un poste de travail pour la conversion en une image graphique qui sera affichée à un utilisateur. [0071] Un exemple de procédé pour l'analyse de la vitesse par migration tomographique peut comprendre la collecte de traces sismographiques provenant d'une formation souterraine et l'utilisation d'un modèle de vitesse initiale pour générer des regroupements d'image commune et un volume d'image de profondeur basé, au moins en partie, sur les traces sismographiques. Un tenseur de structure peut être calculé avec le volume d'image de profondeur pour l'estimation automatisée de l'inclinaison structurale et de l'azimut. Une semblance peut être générée en utilisant ladite pluralité de regroupements d'image commune et ledit tenseur de structure. Les résiduels de profondeur de l'image peuvent être sélectionnés automatiquement à partir de ladite semblance. Un calcul de traçage de rayon peut être réalisé sur lesdits modèles de vitesse en utilisant ledit tenseur de structure. Un modèle de vitesse actualisée peut être généré avec un calcul d'inversion tomographique, dans lequel ledit calcul d'inversion tomographique utilise ladite pluralité des résiduels de l'image de la profondeur et ledit calcul de traçage de rayon. [0072] Dans certains modes de réalisation décrits dans le paragraphe précédent, la collecte des traces sismographiques comprend l'émission d'au moins une onde sismique, et la réception d'une réflexion de l'au moins une onde sismique. Dans certains modes de réalisation décrits dans le présent paragraphe, l'utilisation d'un modèle de vitesse initiale pour générer la pluralité de 18 regroupement d'image commune et le volume de l'image de la profondeur basée, au moins en partie, sur les traces sismographiques comprend la réalisation d'une migration en profondeur sur les traces sismographiques. Dans certains modes de réalisation décrits dans le paragraphe précédent, la génération de la semblance utilisant ladite pluralité de regroupements d'image commune et ledit tenseur de structure comprend la génération d'une semblance utilisant le tenseur de structure comme une contrainte. Dans certains modes de réalisation décrits dans le paragaphe précédent, la sélection automatique d'une pluralité de résiduels de profondeur d'image à partir de ladite semblance comprend la sélection automatique d'une pluralité de résiduels de profondeur d'image utilisant un algorithme de sélection automatique qui maximise les valeurs de semblance dans de multiples directions basées sur une fonction de guide d'entrée et une fourchette de recherche positive et négative relativement au guide d'entrée. [0073] Dans certains modes de réalisation décrits dans le paragraphe précédent, le calcul du tenseur de structure utilisant ledit volume d'image de profondeur pour l'estimation automatisée de l'inclinaison structurale et de l'azimut comprend le calcul de dérivées Gaussiennes lissées dans le volume de l'image de profondeur. [0074] Dans certains modes de réalisation décrits dans les 3 paragraphes précédents, le procédé peut également comprendre l'utilisation du modèle de vitesse actualisée pour générer une pluralité actualisée de regroupement d'image commune et un volume d'image de profondeur. [0075] Dans certains modes de réalisation décrits dans le paragraphe précédent, le procédé peut également comprendre le calcul d'un tenseur de structure actualisé utilisant ledit volume d'image de profondeur actualisé pour l'estimation automatisée actualisée de l'inclinaison structurale et de l'azimut ; la génération d'une semblance actualisée en utilisant ladite pluralité de regroupement d'image commune et le tenseur de structure actualisée ; la sélection automatique d'une pluralité de résiduels d'image de profondeur provenant de ladite pluralité de semblance ; la réalisation d'un calcul de traçage de rayons sur ledit modèle de vitesse actualisée utilisant ledit tenseur de structure actualisé ; et la génération d'un deuxième modèle de vitesse actualisé avec le calcul de l'inversion tomographique. [0076] Dans certains modes de réalisation décrits dans les 5 paragraphes précédents, le procédé peut également comprendre la détermination d'une ou de plusieurs caractéristiques de la formation basé, au moins en partie, sur le modèle de vitesse actualisé. 19 [0077] Dans certains modes de réalisation décrits dans le paragraphe précédent, l'une ou les plusieurs caractéristiques de la formation comprennent des limites de strates de la formation. [0078] Un exemple de système peut comprendre un système de relevé sismique avec au moins une source sismique et au moins un capteur sismique, et un système de traitement d'informations comprenant un processeur et un dispositif de mémoire couplé au processeur. Le dispositif de mémoire peut comprendre un jeu d'instructions qui, lorsqu'il est exécuté par un processeur, permet au processeur de collecter des traces sismographiques provenant d'une formation souterraine, et l'utilisation d'un modèle de vitesse initiale pour générer une pluralité de regroupements d'images communes et un volume de profondeur d'image basée, au moins en partie, sur les traces sismographiques. Le jeu d'instructions peut également permettre au processeur de calculer un tenseur de structure utilisant ledit volume d'image de profondeur pour l'estimation automatisée de l'inclinaison structurelle et de l'azimut, et pour générer une semblance utilisant ladite pluralité de regroupements d'image commune et ledit tenseur de structure. Le jeu d'instructions peut également permettre au processeur de sélectionner automatiquement une pluralité de résiduels de profondeur de l'image à partir de ladite semblance ; la réalisation d'un calcul de traçage de rayon sur lesdits modèles de vitesse initiale utilisant ledit tenseur de structure ; et générer un modèle de vitesse actualisée avec un calcul d'inversion tomographique, dans lequel ledit calcul d'inversion tomographique utilise ladite pluralité de résiduels de profondeur d'image et ledit calcul du traçage de rayon. [0079] Dans certains modes de réalisation décrits dans le paragraphe précédent, les traces sismographiques comprennent au moins une onde sismique reçue au niveau de l'au moins un capteur sismique, dans lequel l'au moins une onde sismique a été générée par au moins une source sismique et réfléchie par une formation souterraine. Dans certains modes de réalisation décrits dans le présent paragraphe, le jeu d'instructions qui permet au processeur d'utiliser le modèle de vitesse initiale pour générer la pluralité de regroupements d'image commune et le volume de l'image de la profondeur basé, au moins en partie, sur les traces sismographiques permet également au processeur de réaliser une migration en profondeur sur les traces sismographiques. [0080] Dans certains modes de réalisation décrits dans le paragraphe précédent, le jeu d'instructions qui permet au processeur de calculer le tenseur de structure utilisant ledit volume d'image de profondeur l'estimation automatisée de l'inclinaison structurale et de l'azimut permet au processeur de calculer des dérivées Gaussiennes lissées dans le volume de l'image de profondeur. 20 [00$1] Dans certains modes de réalisation décrits dans les 3 paragraphes précédents, le jeu d'instructions qui permet au processeur de générer la semblance en utilisant la pluralité de regroupement d'image commune et ledit tenseur de structure permet également au processeur de générer la semblance en utilisant le tenseur de structure comme contrainte. Dans certains modes de réalisation décrits dans les 3 paragraphes précédents, le jeu d'instructions qui permet au processeur de sélectionner automatiquement une pluralité de résiduels de profondeur de l'image à partir de ladite semblance permet également au processeur de sélectionner automatiquement une pluralité de résiduels de profondeur d'image utilisant un algorithme de sélection automatique qui maximise les valeurs de semblance dans de multiples directions, basée sur une fonction de guide d'entrée et une fourchette de recherche positive et négative relativement au guide d'entrée. Dans certains modes de réalisation décrits dans les 3 paragraphes précédents, le jeu d'instructions permet également au processeur d'utiliser un modèle de vitesse actualisée pour générer une pluralité actualisée de regroupement d'image commune et un volume d'image de profondeur. [0082] Dans certains modes de réalisation décrits dans le paragraphe précédent, le jeu d'instructions permet également au processeur de calculer un tenseur de structure actualisé utilisant ledit volume d'image de profondeur actualisé pour l'estimation automatisée actualisée de l'inclinaison structurale et de l'azimut ; la génération d'une semblance actualisée en utilisant ladite pluralité de regroupements d'image commune et ledit tenseur de structure actualisée ; la sélection automatique d'une pluralité de résiduels d'image de profondeur actualisée provenant de ladite pluralité de semblance ; la réalisation d'un calcul de traçage de rayons sur ledit modèle de vitesse actualisée utilisant ledit tenseur de structure actualisé ; et la génération d'un deuxième modèle de vitesse actualisé avec le calcul de l'inversion tomographique. [0083] Dans certains modes de réalisation décrits dans les 5 paragraphes précédents, le jeu d'instructions permet également au processeur de déterminer une ou de plusieurs caractéristiques de la formation basé, au moins en partie, sur le modèle de vitesse actualisé. [0084] Dans certains modes de réalisation décrits dans le paragraphe précédent, l'une ou les plusieurs caractéristiques de la formation comprennent des limites de strates de la formation. [0085] Par conséquent, la présente invention est bien adaptée pour atteindre les finalités et obtenir les avantages mentionnés ici aussi bien que ceux qui sont inhérents à la présente description. Les modes de réalisation particuliers divulgués ci-dessus sont illustratifs seulement, étant donné que la présente invention peut être modifiée et pratiquée de façon différente mais équivalente qui sera 21 apparente aux spécialistes du domaine qui bénéficient des enseignements de la présente description. En outre, aucune limite n'est envisagée aux détails de construction ou de conception décrits ici, autres que ceux décrits dans les revendications ci-jointes. Il est donc évident que les modes de réalisation illustratifs particuliers décrits ci-dessus peuvent être altérés ou modifiés et toutes les variations de ce type sont considérées comme étant à l'intérieur de la portée ou de l'esprit de la présente invention. Mais également, les termes dans les revendications ont leur signification claire et ordinaire, sauf en cas de mention explicite et de définition claire du demandeur. Les articles indéfinis « un » ou « une », sont utilisés dans les revendications et sont tous les deux définis ici pour signifier un plusieurs éléments qu'ils introduisent.