FR 2985571 A1 20130712 FR 1350081 A 20130104 Docket No.100223/0336-101 Estimation conjointe simultanée de statiques résiduelles P-P et P-S APPLICATION CONNEXE [0001] La présente invention concerne et comporte une revendication de priorité fondée sur la demande de brevet US provisoire n° 61/583 231, déposée le 5 janvier 2012, intitulée « Simultaneous Joint Estimation for P-P and P-S Residual Statics », de David Le Meur et Guillaume Poulain, dont la description est incorporée ici par voie de référence. DOMAINE TECHNIQUE [0002] Les modes de réalisation de l'objet présenté ici concernent généralement des procédés et des systèmes pour le traitement de données sismiques terrestres ou marines peu profondes et, plus particulièrement, des mécanismes et des techniques pour estimer des termes statiques pour des ensembles de données multi-composantes. CONTEXTE [0003] Des techniques d'acquisition et de traitement de données sismiques terrestres ou marines peu profondes sont utilisées pour générer un profil (image) d'une structure géophysique (sous-surface) des strates sous-jacentes. Ce profil ne fournit pas nécessairement un emplacement précis des gisements de pétrole et de gaz, mais il peut suggérer, aux hommes du métier, la présence ou l'absence de 2 9 8 5 5 7 1 Docket No. 100223/0336-101 gisements de pétrole et/ou de gaz. Ainsi, la fourniture d'une image améliorée de la sous-surface en une période de temps plus courte est un processus en cours. [0004] L'estimation de termes statiques résiduels pour un ensemble de données multi-composantes est un défi. L'estimation de statiques résiduelles en deux dimensions et trois dimensions pour des données P-P ou pour des ondes converties telles que des données P-S fait partie d'un traitement sismique dans le temps ou en profondeur. Les approches classiques d'évaluation des statiques résiduelles sur des données P-P et des données P-S sont principalement basées sur une inversion linéaire impliquant des fonctions de corrélation croisée comme décrit par J. Ronen et J. F. Claerbout dans leur article de 1985 intitulé « Surface-Consistent Residual Statics Estimation by Stack Power Maximization » publié dans Geophysics et par S. Jin, J. Li et S. Ronen dans leur article de 2004 intitulé « Robust Inversion for Converted Wave Receiver Statics » présenté au 74th Annual International Meeting of the Society of Exploration Geophysicists, ou une cohérence de trace à trace de la section de pile de récepteur commune comme décrit par P. W. Cary et W. S. Eaton dans leur article de 1993 intitulé « A Simple Method for Resolving Large Converted-Wave (P-Sv) Statics » publié dans Geophysics. [0005] En outre, il a été montré que l'estimation de statiques résiduelles de grande amplitude est mieux gérée par un système non linéaire comme décrit par Daniel H. Rothman dans son article de 1985 intitulé « Non-linear Inversion, Statistical Mechanics and Residual Statics Estimation » et son article de 1986 intitulé « Automatic Estimation of Large Residual Statics Corrections », publiés tous deux dans Geophysics, basé sur une méthode de Monte Carlo, comme décrit par D. Le 2 Docket No. 100223/0336-101 Meur et S. Merrer dans leur article intitulé « Monte Carlo Statics : The Last Frontier » présenté à la 2004 Canadian Society of Exploration Geophysicists Annual Convention, et couplé à une approche de recuit simulé comme décrit par K. Vasudevan, W. G. Wilson et W. G. Laidlaw dans leur article de 1991 intitulé « Simulated Annealing Statics Computation Using an Order-based Energy Function » publié dans Geophysics. Par conséquent, dans un traitement d'ondes converties, des statiques de récepteur d'ondes de cisaillement (transversales) sont caractérisées par une grande amplitude qui peut être deux à dix fois supérieure à des valeurs de statiques P-P ainsi que par des données d'entrée plus bruyantes que les données P-P. [0006] Il est connu dans l'art que différents procédés sont utilisés pour calculer des statiques résiduelles cohérentes avec la surface sur des données P-P et des statiques de récepteur sur des données P-S, comme décrit par D. Marsden dans son article de 1993 intitulé « Static Corrections - A Review » publié dans The Leading Edge. La plupart de ces procédés sont basés sur l'utilisation de fonctions de corrélation croisée et une solution d'un système d'équations linéaires, qui donnent très fréquemment une solution minimum locale basée sur la nature des données. Une approche non linéaire, cependant, utilisant le concept de recuit simulé, comme décrit par S. Kirkpatrick, C. D. Gelatt, Jr et M. P. Vehhi dans leur article de 1983 intitulé « Optimization by Simulated Annealing » publié dans Science, couplé à une technique Monte Carlo comme décrit par N. Metropolis, A. Rosenbluth, M. Rosenbluth, A. Teller et E. Teller dans leur article de 1953 intitulé « Equation of State Calculations by Fast Computing Machines » publié dans Journal of Chemistry 3 2 9 8 5 5 7 1 Docket No. 100223/0336-101 and Physics, permet le calcul de n'importe quel type de statiques résiduelles au minimum global. [0007] En examinant maintenant la figure 1 de l'art antérieur, un organigramme 100 pour estimer des statiques résiduelles P-P 110, 112 et des statiques résiduelles P-S 116 est représenté. Le procédé est un calcul en deux passages, le premier passage 102 agissant sur les données P-P 106 et le deuxième passage 104 agissant sur les données P-S 108 et les statiques résiduelles de source P-P 110 issues du premier passage 102. Le premier passage 102 concerne le calcul de statiques résiduelles de source P-P 110 et de statiques résiduelles de récepteur P-P 112 en appliquant un résolveur de statiques non linéaire 114 aux données P-P 106. Le deuxième passage 104 concerne le calcul de statiques résiduelles P-S 116 en appliquant les statiques résiduelles de source P-P 110 aux données P-S 108 et en appliquant ensuite le résolveur de statiques non linéaire 114 aux données P-S 108 mises à jour. La nature à deux passages de cette approche mène, de manière non souhaitable, à des longs temps de calcul basés sur des mécanismes d'accès aux données pour le grand volume de données associées aux données P-P 106 et aux données P-S 108. [0008] Par conséquent, il serait souhaitable de proposer des systèmes et des procédés qui évitent les problèmes et les inconvénients décrits précédemment, et qui améliorent la précision des images finales qui sont produites en tant que résultat de ces études sismiques. 4 Docket No. 100223/0336-101 RÉSUMÉ [0009] Selon un exemple de mode de réalisation, un procédé, mémorisé dans une mémoire et s'exécutant sur un processeur, pour estimer des statiques de réflexion résiduelles comprend l'utilisation simultanée d'une pluralité d'ensembles de données sismiques en un seul passage, dans lequel au moins une partie de ladite pluralité d'ensembles de données sismiques se superposent, l'exécution d'une ou de plusieurs simulations pondérées de la pluralité d'ensembles de données sismiques jusqu'à ce qu'une condition d'arrêt prédéterminée, basée sur une fonction de coût partagé, pour que l'estimation soit satisfaite, et la sortie (« outputting » en terminologie anglo-saxonne) des statiques de réflexion résiduelles estimées simultanément sur la base desdites une ou plusieurs simulations de la pluralité d'ensembles de données sismiques après que la condition d'arrêt a été satisfaite. [0010] Selon un autre exemple de mode de réalisation, un système pour estimer des statiques de réflexion résiduelles comprend un premier ensemble de données sismiques connecté de manière communicante à un processeur, un deuxième ensemble de données sismiques connecté de manière communicante à un processeur dans lequel au moins une partie dudit premier ensemble de données sismiques et dudit deuxième ensemble de données sismiques se superposent, un ou plusieurs processeurs configurés pour exécuter des instructions d'ordinateur, et une mémoire configurée pour mémoriser les instructions d'ordinateur, dans lequel les instructions d'ordinateur comprennent en outre un composant d'estimation pour traiter simultanément le premier ensemble de données sismiques et le deuxième ensemble de données sismiques sur la base d'une fonction de coût partagé et un5 Docket No. 100223/0336-101 composant de sortie pour sortir des statiques de réflexion résiduelles estimées simultanément. 6 Docket No. 100223/0336-101 BRÈVE DESCRIPTION DES DESSINS [0011] Les dessins joints, qui sont incorporés dans la spécification et qui font partie de celle-ci, illustrent un ou plusieurs modes de réalisation et, avec la description, expliquent ces modes de réalisation. Sur les dessins : [0012] la figure 1 est un schéma de l'art antérieur pour un système en deux passages pour estimer des statiques résiduelles P-P et P-S ; [0013] la figure 2 est un schéma de l'art antérieur montrant un système d'acquisition de données sismiques terrestres ; [0014] la figure 3 est un schéma d'un système d'estimation conjointe simultanée d'un exemple de mode de réalisation pour traiter simultanément des données pré-empilées P-P et des données pré-empilées P-S et délivrer des statiques résiduelles P-P et P-S ; [0015] la figure 4 est un organigramme illustrant un procédé pour traiter simultanément des données P-P et P-S et délivrer des statiques résiduelles P-P et P-S ; [0016] la figure 5 est un organigramme illustrant un procédé pour traiter simultanément des données P-P et P-S en utilisant un recuit simulé non linéaire ; [0017] la figure 6 est un schéma d'un système informatisé qui met en oeuvre divers procédés selon un exemple de mode de réalisation ; et 7 Docket No. 100223/0336-101 [0018] les figures 7(a)-(f) sont des schémas composites de données empilées P-P et P-S montrant des données brutes, des données estimées en cascade (deux passages) et des données empilées d'estimation conjointe simultanée. 8 Docket No. 100223/0336-101 DESCRIPTION DÉTAILLÉE [0019] La description qui suit des exemples de modes de réalisation fait référence aux dessins joints. Les mêmes numéros de référence sur différents dessins identifient les mêmes éléments ou des éléments similaires. La description détaillée qui suit ne limite pas l'invention. Au lieu de cela, l'étendue de l'invention est définie par les revendications jointes. Certains des modes de réalisation qui suivent sont examinés, par souci de simplicité, en relation avec la terminologie et la structure d'estimation de statiques résiduelles pour des données P-P et des données P-S. Cependant, les modes de réalisation qui seront examinés ensuite ne sont pas limités à ces configurations, mais peuvent être étendus à d'autres agencements, comme examiné ultérieurement. [0020] Une référence dans toute la spécification à un « mode de réalisation » signifie qu'une fonctionnalité, une structure ou une caractéristique particulière décrite en relation avec un mode de réalisation est incluse dans au moins un mode de réalisation de l'objet présenté. Ainsi, l'apparition de l'expression « dans un mode de réalisation » à différents emplacements dans toute la spécification ne fait pas nécessairement référence au même mode de réalisation. En outre, les fonctionnalités, structures ou caractéristiques particulières peuvent être combinées de n'importe quelle manière appropriée dans un ou plusieurs modes de réalisation. [0021] Afin de fournir un contexte pour les exemples de modes de réalisation suivants, une description d'aspects et de terminologie est incluse dans le présent document. Les procédés et les systèmes décrits ici génèrent et reçoivent des ondes 9 Docket No. 100223/0336-101 P et des ondes S. Une onde P est l'onde étudiée dans des données sismiques classiques et est une onde de corps élastique ou une onde sonore dans laquelle les particules oscillent dans la direction où l'onde se propage. Des ondes P frappant une interface autrement que selon une incidence normale peuvent produire des ondes S réfléchies et transmises, autrement connues en tant qu'ondes converties. Il conviendrait de noter que les données P-P sont des données sismiques associées à une onde P se propageant. [0022] Une onde S, générée par la plupart des sources sismiques terrestres et parfois en tant qu'ondes P converties, est une onde de corps élastique dans laquelle les particules oscillent perpendiculairement à la direction dans laquelle l'onde se propage. Les ondes S, connues également en tant qu'ondes de cisaillement, se propagent plus lentement que les ondes P et ne peuvent pas se propager à travers les fluides parce que les fluides ne supportent pas un cisaillement. Il conviendrait de noter que des ondes S peuvent également être converties en des ondes P. L'enregistrement d'ondes S nécessite des récepteurs couplés à la terre solide et leur interprétation peut permettre la détermination des propriétés des roches telles que la densité et l'orientation des fractures, un rapport de Poisson, et un type de roche par un tracé croisé des vitesses d'ondes P et d'ondes S et d'autres techniques. Il conviendrait de noter que les données P-S sont des données sismiques correspondant à une propagation d'ondes P descendantes, et ensuite d'ondes S montantes. [0023] Dans une autre description de la terminologie, un point de tir est l'un d'un certain nombre d'emplacements ou de stations à la surface de la terre auquel 10 Docket No. 100223/0336-101 une source sismique est activée. Une trace sismique consiste en des données sismiques enregistrées, par un canal, après que la source sismique a été déclenchée. La trace sismique représente la réponse du champ d'ondes élastiques à des contrastes de vitesse et de densité à travers des interfaces de couches de roche ou des sédiments alors que l'énergie se propage de la source sismique à travers la sous-surface vers un récepteur ou un ensemble de récepteurs. En outre, une inversion sismique est un processus de transformation des données de réflexion sismiques en une description de propriété quantitative d'une description des strates d'un emplacement souterrain et éventuellement d'un gisement contenant des ressources naturelles telles que du pétrole ou du gaz. [0024] En outre, dans le contexte des exemples de modes de réalisation suivants, il conviendrait de noter que l'estimation simultanée des statiques de réflexion résiduelles pour des sources et des récepteurs est effectuée en utilisant conjointement à la fois les données P-P et les données P-S. Selon un autre aspect d'un exemple de mode de réalisation, l'estimation conjointe simultanée (SJE « Simultaneous Joint Estimation » en terminologie anglo-saxonne) peut utiliser une approche linéaire ou non linéaire. En poursuivant avec un exemple de mode de réalisation, une estimation conjointe simultanée est un algorithme à un seul passage qui permet la déduction simultanée des statiques de réflexion résiduelles de source commune à la fois pour les données P-P et les données P-S et des statiques résiduelles de récepteurs séparés, une pour les données P-P et une pour les données P-S. Selon un autre aspect d'un exemple de mode de réalisation, le modèle susmentionné est opérationnel indépendamment de la mise en oeuvre logicielle, c'est-à-dire, d'un traitement parallèle ou non parallèle, de la taille des bases de 11 Docket No. 100223/0336-101 données, des dimensions des données, c'est-à-dire, bidimensionnelles (2D) ou tridimensionnelles (3D), et de la conception d'acquisition des données P-P et des données P-S. [0025] En poursuivant avec la description de contexte général d'un exemple de mode de réalisation, l'estimation conjointe simultanée utilise une fonction de coût unique qui garantit une cohérence pour le calcul des trois termes statiques de réflexion résiduelles. Selon un autre aspect du contexte général pour les exemples de modes de réalisation, la fonction de coût bénéficie du coefficient de pondération des données P-P plus propres pour estimer à la fois des statiques de réflexion résiduelles de sources communes et de récepteurs spécifiques pour les données P-S. En outre, dans le contexte général pour les exemples de modes de réalisation, l'approche d'estimation conjointe simultanée offre un avantage économique par une réduction du temps de traitement pour le traitement des données P-P et des données P-S pour estimer les statiques de réflexion résiduelles. [0026] Les exemples de modes de réalisation ont généralement une pluralité d'entrées de données sismiques telles que, mais sans y être limitées, des récipients, dans lesquels leurs sections de pile se superposent partiellement ou totalement, des stations de source, des stations de récepteur et n'importe quelles combinaisons des entrées de données sismiques référencées ci-dessus. En outre, les exemples de modes de réalisation ont généralement des sorties telles que, mais sans y être limitées, des statiques de réflexion résiduelles de source pour un premier ensemble de données sismiques, des statiques de réflexion résiduelles de récepteur pour un premier ensemble de données sismiques, des statiques de réflexion résiduelles de 12 2 9 8 5 5 7 1 Docket No. 100223/0336-101 source pour un deuxième ensemble de données sismiques, des statiques de réflexion résiduelles de récepteur pour un deuxième ensemble de données sismiques et n'importe quelles combinaisons des sorties de données sismiques référencées ci-dessus. En poursuivant avec une description générale des exemples de modes de réalisation, une inversion conjointe simultanée entre et traite simultanément la pluralité d'ensembles de données sismiques par une fonction de coût linéaire ou non linéaire partagé jusqu'à ce qu'une condition d'arrêt prédéterminée soit satisfaite, et des statiques de réflexion résiduelles, comme décrit précédemment, sont sorties. [0027] En examinant maintenant la figure 2, un schéma de contexte illustre les aspects décrits précédemment. Un tir est déclenché au niveau d'un point de tir 202 sur la surface 204 de la terre pour propager une série d'ondes 206 à travers différentes couches 208, 210 sous la surface 204. Les ondes sont réfléchies par les interfaces au niveau des différentes couches 208, 210 à différents instants après que le tir a été déclenché et selon différents angles. Les ondes réfléchies sont détectées par des récepteurs 212 attachés à la terre. Un dispositif d'enregistrement 214 collecte les données sismiques provenant des récepteurs et enregistre les données pour une analyse future. Il conviendrait de noter, dans l'exemple de mode de réalisation, que des calculs concernant les données enregistrées peuvent être effectués dans le dispositif d'enregistrement 214 ou ils peuvent être effectués à un autre emplacement après que les données sismiques ont été transférées. [0028] En examinant maintenant la figure 3, un système de traitement d'estimation conjointe simultanée 300 est représenté. Le système d'estimation 13 Docket No. 100223/0336-101 conjointe simultanée 300 comprend une base de données pré-empilées P-P 302 une base de données pré-empilées P-S 304, un moteur d'estimation conjointe simultanée (composant d'estimation) 306 et génère des statiques de source de réflexion pour les données P-P et P-S 308, des statiques de récepteur de réflexion pour les données P-P 310 et des statiques de récepteur de réflexion pour les données P-S 312. Il conviendrait de noter, dans l'exemple de mode de réalisation que le moteur d'estimation conjointe simultanée 306 peut utiliser une approche linéaire ou non linéaire pour déterminer les statiques de réflexion 308, 310, 312. On devrait en outre noter, dans l'exemple de mode de réalisation, que le moteur d'estimation conjointe simultanée 306 est également appelé composant d'estimation. En poursuivant avec l'exemple de mode de réalisation, le moteur d'estimation conjointe simultanée 306 lit à la fois les données pré-empilées P-P et les données pré-empilées P-S pour une opération simultanée, c'est-à-dire, une opération à un seul passage. [0029] En poursuivant avec l'exemple de mode de réalisation, le processus d'estimation conjointe simultanée tire profit de techniques informatiques à haute performance permettant l'optimisation d'opérations d'accès résultant en une minimisation du temps d'accès aux données pour des opérations, telles qu'une inversion non linéaire. Par exemple, une approche de Monte Carlo utilise une fonction de coût qui est basée sur la cohérence de la pile de données 302, 304 avec un critère robuste pour stabiliser les résultats. En poursuivant avec l'exemple, les débits d'entrée et/ou de sortie sont le goulot d'étranglement principal lorsqu'un procédé non linéaire basé sur un recuit simulé est appliqué sur de grands ensembles de données de large azimut. A chaque étape de simulation de l'exemple, 14 Docket No. 100223/0336-101 les stations sont visitées dans un ordre aléatoire et, pour chaque station, plusieurs rassemblements associés comprenant des tirs, des récepteurs, des données de point milieu commun (CMP « common mid-point » en terminologie anglo-saxonne) ou de point de conversion asymptotique (ACP « asymptotic conversion point » en terminologie anglo-saxonne) sont traités pour calculer la fonction de coût conjointe, c'est-à-dire qu'un accès à l'ensemble de données pré-empilées d'entrée doit être effectué plusieurs centaines de fois dans un ordre aléatoire. En outre, dans l'exemple, pour une étape de simulation donnée, les corrections des statiques résiduelles 308, 310, 312 ne peuvent pas être calculées de manière indépendante parce que les stations se superposent dans les domaines CMP et ACP, c'est-à-dire que l'algorithme ne peut pas être parallélisé massivement sur les stations. [0030] Ensuite, dans l'exemple de mode de réalisation, le moteur d'estimation conjointe simultanée 306 utilise une équation d'énergie globale : E(P-P, P-S) = aE(PP) + [3E(PS) avec a et 3 E [0,1] et agissant en tant que terme de pondération pour équilibrer les termes de composante d'énergie sur la base de facteurs tels que, mais sans y être limités, la qualité des données sismiques, et où E(PP) et E(PS) peuvent être écrits en tant que : E(PP) = EZ {Edh[t + + Rpp j(y,h) X {EdhY+1[t + S,(y+1,1,) + Rpp j(y+i,h) y t h h E(PS ) = EZ {Edh[t + + Rpsxy x {EdYhlt + S(y+u) + Rpsi(y+u)D y t h h où Si est le terme de source commune P-P/P-S à une position de surface I ; Rppi est le terme de récepteur P-P à une position de surface j ; et Rpsi est le terme de récepteur P-S à une position de surface j. A la fois pour l'énergie E(PP) et E(PS), y 15 2 985 5 7 1 Docket No. 100223/0336-101 est la somme sur les CDP ou les ACP ; t est la somme sur les échantillons temporels des traces ; et h est la somme sur les décalages dans un rassemblement CDP ou ACP. Les données de trace sont représentées par : dhY [t Si(y,h) R i(Y,h) où Si(y,h) et Ri(y,h) sont les statiques de tir et de récepteur pour le CDP ou l'ACP ; et y et le décalage h décalent les données de trace soit vers le haut, soit vers le bas. Il conviendrait de noter, dans l'exemple de mode de réalisation, que pour chaque station aléatoire sélectionnée qui peut être une source P-P, une source P-S, un récepteur P-P ou un récepteur P-S commun, un décalage de statiques aléatoire est choisi et une nouvelle énergie globale, c'est-à-dire, Enew(P-P, P-S), est calculée en utilisant l'équation d'énergie globale définie précédemment. En poursuivant avec l'exemple de mode de réalisation, le décalage de statiques est conservé lorsque la nouvelle valeur d'énergie globale calculée est supérieure à la valeur d'énergie globale calculée précédente. On devrait en outre noter, dans l'exemple de mode de réalisation, qu'un critère de Métropolis peut être utilisé pour s'échapper d'un minimum local potentiel et conserver le décalage de statiques sélectionné. Il conviendrait de noter, dans l'exemple de mode de réalisation, qu'une fois que toutes les stations ont été visitées, l'itération de simulation se termine. En outre, lorsque la différence entre Enew et E est dans une distance pré-configurée par rapport à zéro, le critère d'arrêt est atteint et l'estimation conjointe simultanée se termine, permettant à un composant de sortie (non montré) de sortir des statiques résiduelles de source commune à la fois pour les données P-P et les données P-S, ainsi que les statiques résiduelles de récepteur pour les données P-P et les données P-S à estimer. Selon 16 2 9 8 5 5 7 1 Docket No. 100223/0336-101 un autre aspect de l'exemple de mode de réalisation, le coefficient de pondération des données P-P et des données P-S pendant le processus d'estimation est commandé par les paramètres « a » et « r3 » alors qu'ils sont modifiés de zéro à un. [0031] En examinant maintenant la figure 4, un exemple de mode de réalisation de procédé d'une estimation conjointe simultanée (SJE) 400 est représenté. En commençant à l'étape 402, le procédé SJE 400 lit à la fois les données pré-empilées P-P et les données pré-empilées P-S. Ensuite, à l'étape 404, le procédé SJE 400 calcule, de manière itérative, l'énergie globale associée aux données jusqu'à ce qu'un critère pré-configuré, basé sur l'énergie globale, soit atteint. En continuant à l'étape 406 de l'exemple de mode de réalisation de procédé, le procédé SJE 400 sont des statiques résiduelles de sources communes à la fois pour les données P-P et les données P-S et des statiques résiduelles de récepteurs spécifiques pour les données P-P et les données P-S. [0032] En examinant maintenant la figure 5, un exemple de mode de réalisation de procédé décrit en outre les étapes itératives 500 d'une estimation conjointe. En commençant à l'étape 502 de l'exemple de mode de réalisation de procédé, une nouvelle simulation commence. Ensuite, à l'étape 504 de l'exemple de mode de réalisation de procédé, une station est sélectionnée de manière aléatoire pour un traitement supplémentaire. Il conviendrait de noter, dans l'exemple de mode de réalisation, que l'ordre de sélection de station aléatoire doit être différent pour chaque simulation. On devrait en outre noter, dans l'exemple de mode de réalisation, que, pour chaque simulation, toutes les stations, à la fois de source et de récepteur, doivent être visitées. En poursuivant à l'étape 506 de l'exemple de mode de 17 ^ Docket No. 100223/0336-101 réalisation, une énergie globale « Eold(PP, PS) » est calculé, en utilisant à la fois des données P-P et des données P-S, pour la station sélectionnée de manière aléatoire. Il conviendrait de noter, dans l'exemple de mode de réalisation, que l'énergie globale est calculée en utilisant les équations d'énergie globale précédemment définies. Ensuite, à l'étape 508, pour chaque station sélectionnée de manière aléatoire (par exemple, de source ou de récepteur), l'exemple de mode de réalisation de procédé choisit un décalage de statiques aléatoire. Ensuite, à l'étape 510 de l'exemple de mode de réalisation de procédé, une autre énergie globale « Enew(PP, PS) » est calculée, sur la base du décalage de statiques aléatoire. Il conviendrait de noter, dans l'exemple de mode de réalisation, que l'énergie globale est calculée en utilisant les équations d'énergie globale précédemment définies. En poursuivant avec l'étape 512 de l'exemple de mode de réalisation, une différence d'énergie globale (Eud et Enew) est calculée, c'est-à-dire, AE(PP, PS) = Enew - Eoid- [0033] En poursuivant à l'étape 514 de l'exemple de mode de réalisation de procédé, il est déterminé si AE est supérieur à zéro ou inférieur à zéro. Si AE est supérieur à zéro, alors l'exemple de mode de réalisation de procédé passe à l'étape 518 à laquelle le décalage de statiques sélectionné de manière aléatoire est conservé et le procédé passe à l'étape 520. Si AE est inférieur à zéro, alors l'exemple de mode de réalisation de procédé passe à l'étape 516 et il est déterminé si le calcul est un minimum local. Si le calcul est un minimum local, alors l'exemple de mode de réalisation de procédé passe de nouveau à l'étape 518 à laquelle le décalage de statiques sélectionné de manière aléatoire est conservé et le procédé passe à l'étape 520. Si le calcul n'est pas un minimum local, alors l'exemple de 18 Docket No. 100223/0336-101 mode de réalisation de procédé passe à l'étape 520. Ensuite, à l'étape 520 de l'exemple de mode de réalisation, il est déterminé si toutes les stations ont été vérifiées. Si toutes les stations n'ont pas été vérifiées, alors le procédé d'exemple de mode de réalisation retourne à l'étape 504, sélectionne de manière aléatoire une autre station, et une autre itération de la simulation commence. Si toutes les stations ont été vérifiées, alors l'exemple de mode de réalisation de procédé passe à l'étape 522 et il est déterminé si le critère d'arrêt a été atteint. Si le critère d'arrêt n'a pas été atteint, alors l'exemple de mode de réalisation de procédé retourne à l'étape 502 et commence une autre simulation. Si le critère d'arrêt a été atteint, alors l'exemple de mode de réalisation de procédé est préparé pour délivrer des statiques de source et de récepteur. [0034] Un ou plusieurs des procédés examinés ci-dessus peuvent être mis en oeuvre dans un système informatisé, comme montré sur la figure 6. Un tel système informatisé 600 peut recevoir, par l'intermédiaire de l'interface d'entrée/sortie 602, des informations pertinentes pour des données P-P et des données P-S. De plus, le système informatisé 600 peut comprendre un processeur 604 pour traiter les données indiquées ci-dessus et pour calculer, par exemple, l'énergie globale d'une station. L'interface 602 et le processeur 604 sont connectés à un bus 606. En outre, le système informatisé 600 peut comprendre une mémoire 606 pour mémoriser les données indiquées ci-dessus, un afficheur 610, une connexion 612 aux données P-P et aux données P-S et d'autres éléments communs pour un système informatisé ou un serveur, comme cela serait reconnu par les hommes du métier. 19 Docket No. 100223/0336-101 [0035] Les résultats d'un exemple de mode de réalisation de l'estimation conjointe simultanée (SJE) susmentionnée sont comparés à une approche en cascade (deux passages) de l'art antérieur et sont présentés avec les données brutes, sur les figures 7(a) à 7(f). Dans cet exemple, un ensemble de données tridimensionnelles avec mille cinq cent points de tir et quatre milles stations de récepteur est analysé. Dans l'exemple de mode de réalisation, le composite d'images empilées 700 représente les données brutes empilées P-P 702 sur la figure 7(a), les données brutes empilées P-S 704 sur la figure 7(d), les données empilées P-P estimées en cascade (deux passages) 706 de l'art antérieur sur la figure 7(b), les données empilées P-S estimées en cascade (deux passages) 708 de l'art antérieur sur la figure 7(e), les données empilées P-P estimées à un seul passage SJE 710 sur la figure 7(c) et les données empilées P-S estimées à un seul passage SJE 712 sur la figure 7(f). Il conviendrait de noter, dans l'exemple de mode de réalisation, que la distribution des tirs et des récepteurs est irrégulière sur la base des contraintes du site et a introduit une variabilité élevée de la couverture de fois à l'intérieur de chaque rassemblement CMP et ACP. On devrait en outre noter que l'estimation SJE 710, 712 est comparée à un flux en cascade classique pour une ligne bidimensionnelle extraite du volume P-P/P-S 706, 708. [0036] En poursuivant avec les résultats de l'exemple de mode de réalisation, dans le flux en cascade de contexte, les statiques résiduelles P-P pour les sources et les récepteurs utilisaient une approche Monte-Carlo non linéaire pour estimer les statiques résiduelles au cours du premier passage, le terme de source commune était ensuite appliqué aux données P-S avant d'estimer le terme de récepteur P-S au cours d'un deuxième passage. Dans le flux SJE, tous les termes résiduels sont 20 2 9 8 5 5 7 1 Docket No. 100223/0336-101 estimés au cours d'un seul passage. En examinant d'abord les données P-P, après l'application des statiques déduites du flux en cascade, on peut voir que de petites valeurs de statiques de source et de récepteur ont été déterminées, lesquelles améliorent la focalisation sur les réflecteurs principaux dans le bassin, sur la base des données brutes, comme indiqué au niveau des flèches 714, 716. En comparant la SJE de l'exemple de mode de réalisation, de meilleurs résultats sont obtenus en un temps de calcul plus court, comme illustré au niveau des flèches 718, 720 et caractérisés par une meilleure continuité des événements horizontaux dans la partie droite de la section. [0037] En examinant maintenant les données P-S, après l'application des statiques déduites du flux en cascade, la pile ACP P-S comprimée est plus focalisée 708 comparée aux données P-S brutes 704 et la continuité des réflecteurs a été partiellement rétablie, comme indiqué par les flèches 722, 724. En examinant une analyse SJE des mêmes données brutes 704, des améliorations de la continuité des événements horizontaux dans la partie supérieure de la section et sur les réflecteurs principaux, par rapport aux calculs en cascade, sont évidentes au niveau des flèches 726, 728. [0038] Pour résumer l'exemple, les problèmes de statiques P-S sont mieux résolus avec la SJE. Le calcul de statiques résiduelles de récepteur P-S bénéficie de la contribution des données P-P pendant l'estimation conjointe. En outre, une continuité latérale des événements dans les données P-P apporte de fortes contraintes, permettant un calcul plus fiable des statiques résiduelles P-S. De plus, les améliorations obtenues par rapport à la pile P-P après le flux simultané montrent 21 - Docket No. 100223/0336-101 que la contribution des données P-S n'est pas négligeable dans la zone dans laquelle une qualité des données P-P est plus faible qu'une qualité des données P-S, ce qui montre que la SJE bénéficie à la fois des ensembles de données P-P et P-S, fournissant ainsi des meilleurs résultats qu'un calcul en cascade (à deux passages). [0039] Les exemples de modes de réalisation présentés ci-dessus fournissent un système et un procédé pour calculer simultanément tous les termes statiques résiduels. On devrait comprendre que cette description n'est pas destinée à limiter l'invention. Au contraire, les exemples de modes de réalisation sont destinés à couvrir les variantes, les modifications et les équivalents, qui sont inclus dans l'esprit et l'étendue de l'invention telle que définie par les revendications jointes. En outre, dans la description détaillée des exemples de modes de réalisation, de nombreux détails spécifiques sont exposés afin de fournir une compréhension détaillée de l'invention revendiquée. Cependant, un homme du métier comprendrait que divers modes de réalisation peuvent être mis en pratique sans ces détails spécifiques. [0040] Bien que les caractéristiques et les éléments des présents exemples de modes de réalisation soient décrits dans les modes de réalisation dans des combinaisons particulières, chaque caractéristique ou élément peut être utilisé seul sans les autres caractéristiques et éléments des modes de réalisation ou dans diverses combinaisons avec ou sans d'autres caractéristiques et éléments présentés ici. En outre, on doit noter que les modes de réalisation ci-dessus peuvent être mis en oeuvre dans un logiciel, un matériel ou une combinaison de ceux-ci. On doit également noter que, bien que les exemples de modes de réalisation précédemment 22 Docket No. 100223/0336-101 décrits fassent référence à une acquisition de données sismiques terrestres, les procédés et les systèmes décrits ici sont également applicables à une acquisition de données sismiques marines, telle OBC/OBN. [0041] Cette description écrite utilise des exemples de l'objet présenté pour permettre à n'importe quel homme du métier de mettre en pratique le susdit, y compris la fabrication et l'utilisation de n'importe quels dispositifs ou systèmes et l'exécution de n'importe quels procédés incorporés. L'étendue brevetable de l'objet est définie par les revendications, et peut comprendre d'autres exemples qui apparaissent aux hommes du métier. Ces autres exemples sont destinés à être dans l'étendue des revendications.