Brevet FR2990519A1 — publié le 20131115

Procédé, dispositif et algorithme de traitement pour l'élimination de multiples et de bruit de données sismiques marines CONTEXTE DOMAINE TECHNIQUE [1] 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 marines et, plus particulièrement, des mécanismes et des techniques pour retirer des multiples et un bruit de données sismiques marines. EXAMEN DU CONTEXTE [2] L'acquisition et le traitement de données sismiques marines génèrent un profil (image) de la structure géophysique sous les fonds océaniques. Bien que ce profil ne fournisse pas un emplacement précis des gisements de pétrole et de gaz, il suggère, aux hommes du métier, la présence ou l'absence de ceux-ci. Ainsi, la fourniture d'une image de haute résolution de la sous-surface est un processus en cours. [3] Généralement, une source sismique est utilisée pour générer un signal sismique qui se propage dans la terre, et qui est au moins partiellement réfléchi par divers réflecteurs sismiques dans la sous-surface. Les ondes réfléchies sont enregistrées par des récepteurs sismiques. Les récepteurs sismiques peuvent être placés sur le fond de l'océan, à proximité du fond de l'océan, au-dessous d'une surface de l'eau, à proximité de la surface de l'eau, à proximité de la surface de la terre, ou dans des trous de forage dans la terre. Les ensembles de données sismiques enregistrés, par exemple, un temps de propagation, peuvent être traités pour produire des informations concernant l'emplacement des réflecteurs de sous-surface et les propriétés physiques des formations de sous-surface, par exemple, pour générer une image de la sous-surface. [004] De nombreux ensembles de données marines souffrent de niveaux élevés de bruit, qui rendent la tâche de traitement et d'interprétation difficile. Cela est plus prononcé pour les ensembles de données à faibles fois. Les ensembles de données à fois élevées d'un capteur unique moderne peuvent également présenter des niveaux de bruit élevés du fait d'un couplage médiocre et d'une onde de surface ou de boue. Pour ces ensembles de données, il peut être plus pratique de réduire le niveau de bruit plutôt que d'interpoler encore plus densément. [005] En outre, les ensembles de données marines sont également affectés par des multiples. Les multiples apparaissent lorsqu'une énergie sismique est réfléchie entre la surface de la mer et le fond de l'océan ou entre la surface de la mer et diverses structures de réflecteur de la sous-surface. Le concept des multiples est illustré sur les figures 1A à 1C. La figure 1A montre une source 10 générant une onde sismique 12 qui se propage vers le bas, vers le fond de l'océan 14. L'onde sismique 12 se propage plus loin dans la sous-surface et est réfléchie au niveau d'un réflecteur 16. Ensuite, l'onde réfléchie 18 (réflexion principale) se propage vers le haut, vers la surface de la mer 22. Un récepteur 20, situé dans l'eau, est configuré pour enregistrer l'onde réfléchie 18. On doit noter que l'onde de réflexion principale ne subit qu'une seule réflexion. Une analyse du signal sismique généré par la réflexion principale fournit des 4.% 2990519 3 informations concernant la caractéristique géologique (réflecteur) responsable de la réflexion de l'onde sismique 12. [6] Cependant, en plus des ondes sismiques montrées sur la figure 1A, il existe d'autres ondes qui ne suivent pas les trajets décrits ci-dessus 5 et qui forment les multiples. Par exemple, la figure 1B montre une onde sismique 24 qui se propage vers le haut et qui est réfléchie par la surface de la mer 22. Ensuite, l'onde réfléchie 26 se propage vers le fond de l'océan 14 et est de nouveau réfléchie. Ce processus peut se poursuivre avec l'onde réfléchie par le réflecteur 16 avant d'être enregistrée par le récepteur 20. D'autres trajets qui sont 10 différents de ceux montrés sur les figures 1A et 1B et qui forment les multiples peuvent être imaginés. Un exemple supplémentaire est montré sur la figure 1C sur laquelle une onde sismique 28 est réfléchie par d'autres réflecteurs 30 et 32 et également par la surface de la mer 22 avant d'être enregistrée par le récepteur 20. Ainsi, dans ce cas, l'onde sismique enregistrée 34 est une onde descendante 15 alors que, précédemment, l'onde enregistrée était une onde montante. [7] Les ondes sismiques enregistrées par le récepteur 20 sur les figures 1B et 1C ont un temps de propagation de la source au récepteur différent de celui de l'énergie qui suit le trajet principal de la figure 1A. Par ailleurs, un multiple subit plusieurs réflexions. Ainsi, ces ondes sismiques sont une source 20 indésirable de contamination des données sismiques parce qu'elles ont tendance à obscurcir l'interprétation des données produites par la réflexion principale. [8] On sait que le bruit et les multiples sont le plus souvent présents dans les données sismiques enregistrées. Il existe des procédés classiques pour éliminer le bruit ou les multiples des données sismiques enregistrées. Par exemple, un tel procédé peut reposer sur le concept de la séparation des champs d'ondes, c'est-à-dire, la séparation des composantes montantes (l'onde 18 sur la figure 1A) des composantes descendantes (l'onde 34 sur la figure 10). Cependant, les résultats de la séparation des champs d'ondes sont parfois affectés par des niveaux élevés de bruit sur une composante verticale Z, tandis qu'une composante de pression P est généralement de bonne qualité. On doit noter que la composante verticale Z est enregistrée par un géophone ou un accéléromètre, tandis que la composante de pression P est enregistrée par un hydrophone, et ces dispositifs forment, ensemble, le récepteur 20. [009] La composante Z est utilisée pour un traitement sismique au fond de l'océan (OBS) et, en général, lors d'un traitement de données de multiples composants, comprenant des données de flûte à multiples composants. Par exemple, Zabihi et d'autres, (« Enhanced wavefield separation of OBC data », 73rd EAGE conference and exhibition, expanded abstract, 2011) indiquent que «Z est nécessaire pour réaliser une séparation des champs d'ondes de pré- empilage complète et également pour commander des processus tels qu'une formation d'image miroir et une déconvolution haut-bas ». [0010] En utilisant la composante Z, trois stratégies d'atténuation de bruit sont possibles. Un bruit peut être retiré avant une séparation des champs d'ondes, pendant une séparation des champs d'ondes, ou après une séparation des champs d'ondes. Plusieurs algorithmes d'atténuation de bruit connus dans l'art et mis en pratique dans des centres de traitement de données sismiques tentent une atténuation du bruit avant une séparation des champs d'ondes (par exemple, Craft, « Geophone noise attenuation and wavefield separation using multi-dimensional decomposition technique », 70th EAGE conference and exhibition, expanded abstract G037, 2008). [0011] Le document de Zabihi et d'autres propose un procédé pour atténuer le bruit pendant une séparation des champs d'ondes, en utilisant le fait qu'une séparation en des ondes montantes et descendantes simplifie le champ d'ondes enregistré et permet un fonctionnement amélioré des algorithmes de traitement de signal. [0012] En ce qui concerne les multiples, un procédé (voir Lokhstanov, « Suppression of sea-floor effects from multicomponent seafloor data », 62nd EAGE conference and exhibition, expanded abstract L52, 2000), propose un procédé adaptatif pour le retrait des effets de la surface libre du champ d'ondes montant. Ce procédé implique la déduction d'un modèle de multiples en utilisant des résultats de déconvolution haut-bas et en exécutant ensuite une soustraction adaptative de ce modèle du champ d'ondes montant. Cependant, parce que le champ d'ondes montant et le modèle sont tous deux généralement affectés par un bruit, cette approche n'effectue pas une atténuation du bruit. [0013] Ainsi, on peut voir que les procédés existants retirent soit le bruit, soit les multiples pendant ou avant une séparation des champs d'ondes. Par conséquent, il est nécessaire de proposer un nouveau mécanisme ou procédé qui effectue conjointement une atténuation de bruit et de multiples. RESUME [0014] Selon un exemple de mode de réalisation, un procédé pour retirer un bruit et atténuer des multiples simultanément dans des données sismiques enregistrées par des récepteurs sismiques est proposé. Le procédé comprend une étape de réception des données sismiques, où les données sismiques comprennent une composante de pression P et une composante verticale Z; une étape de séparation des données sismiques en des champs d'onde U montants et des champs d'ondes D descendants ; une étape de calcul d'une déconvolution haut-bas R sur la base des champs d'ondes U montants et D descendants ; une étape de génération d'un modèle de multiples M sur la base de la déconvolution haut-bas R; et une étape de soustraction adaptative du modèle de multiples M de la composante de pression P pour obtenir une composante de pression corrigée P - free de laquelle le bruit et les multiples ont été retirés simultanément. Un multiple est un champ d'ondes qui se propage d'une source vers un récepteur et qui subit plusieurs réflexions. [0015] Selon un autre exemple de mode de réalisation, il existe un dispositif informatique pour retirer le bruit et atténuer les multiples simultanément dans des données sismiques enregistrées par des récepteurs sismiques. Le dispositif informatique comprend une interface configurée pour recevoir les données sismiques, où les données sismiques comprennent une composante de pression P et une composante verticale Z; et un processeur connecté à l'interface. Le processeur est configuré pour séparer les données sismiques en des champs d'onde U montants et des champs d'ondes D descendants, calculer une déconvolution haut-bas R sur la base des champs d'ondes U montants et D descendants, générer un modèle de multiples M sur la base de la déconvolution haut-bas R, et soustraire de manière adaptative le modèle de multiples M de la composante de pression P pour obtenir une composante de pression corrigée Pfree de laquelle le bruit et les multiples ont été retirés simultanément. Un multiple est un champ d'ondes qui se propage d'une source vers un récepteur et qui subit plusieurs réflexions. [0016] Selon encore un autre exemple de mode de réalisation, il existe un support pouvant être lu par un ordinateur configuré pour mémoriser en permanence des instructions d'ordinateur de manière à, lorsqu'un ordinateur exécute les instructions, simultanément, retirer le bruit et atténuer les multiples dans des données sismiques enregistrées par des récepteurs sismiques. Les instructions exécutent les étapes de procédé indiquées ci-dessus. BREVE DESCRIPTION DES DESSINS [0017] Les dessins joints, qui sont incorporés dans la description, et qui constituent une 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 : [0018] les figures 1A à 1C illustrent une composante principale et des multiples enregistrés par un récepteur sismique ; [0019] la figure 2 est un organigramme illustrant un procédé pour retirer le bruit et retirer les multiples de données sismiques selon un exemple de mode de réalisation ; [0020] la figure 3 illustre des données synthétiques enregistrées par un hydrophone ; [0021] la figure 4 illustre des données synthétiques enregistrées par un géophone ; [0022] les figures 5 et 6 illustrent des composantes montantes et descendantes séparées des données des figures 3 et 4 selon un exemple de mode de réalisation ; [0023] la figure 7 illustre une déconvolution haut-bas calculée selon un exemple de mode de réalisation ; [0024] la figure 8 illustre une composante de pression de laquelle les multiples et le bruit ont été retirés selon un exemple de mode de réalisation ; [0025] la figure 9 est un schéma d'un algorithme original pour retirer un bruit et des multiples de données sismiques selon un exemple de mode de réalisation ; et [0026] la figure 10 est un schéma d'un dispositif informatique pour exécuter l'algorithme original selon un exemple de mode de réalisation. DESCRIPTION DETAILLEE [0027] 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 les 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. Les modes de réalisation qui suivent sont examinés, par souci de simplicité, en relation avec des données sismiques d'OBC dont le bruit a été retiré et les multiples ont été atténués après la séparation des champs d'ondes. Cependant, les modes de réalisation qui seront examinés ensuite ne sont pas limités à l'OBC, étant donné qu'ils peuvent être étendus ou utilisés avec n'importe quelles données marines, par exemple, des données de flûte. [0028] Une référence dans toute la description à « 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 » à divers emplacements dans toute la description 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. [0029] Selon un exemple de mode de réalisation, un procédé original atténue un bruit et des multiples conjointement après une séparation des champs d'ondes. Ainsi, dans un mode de réalisation, des données sismiques brutes sont reçues par une application d'ordinateur et/ou un dispositif informatique approprié et elles sont séparées en ondes montantes (U) et descendantes (D). La séparation peut avoir lieu au-dessus ou au-dessous du fond de l'océan. Le procédé effectue également une déconvolution haut-bas pour obtenir une réponse R de la terre, où R est une estimation du champ d'ondes principal. Plus spécifiquement, R peut être compris, si on utilise un exemple à une dimension (unidimensionnel) de la terre, comme étant la déconvolution de U et de D, où D correspond à l'énergie qui entre dans la terre pendant la durée d'étude et U correspond à l'énergie qui sort de la terre en tant que réponse de la terre à l'énergie correspondant à D. En utilisant le langage de traitement de signal, D correspond à l'énergie qui va dans un filtre, R est le filtre, et U correspond à la sortie du filtre. [0030] Sur la base de la déconvolution haut-bas R et de la composante de pression P des données sismiques enregistrées, un modèle M pour les multiples est défini. Le modèle M peut également comprendre un terme de repositionnement par rapport au fond de l'océan. Le modèle est ensuite soustrait de manière adaptative de la composante de pression P pour obtenir les données sismiques à bruit retiré et multiples atténués. [0031] L'algorithme de traitement original indiqué ci-dessus est maintenant examiné plus en détail. Selon un organigramme illustré sur la figure 2, l'algorithme original comprend une étape 200 de réception de données sismiques. Les données sismiques reçues comprennent une composante Z et une composante P. Par exemple, la figure 3 montre une composante P synthétique et la figure 4 montre une composante Z synthétique. Les données sismiques reçues sont enregistrées, par exemple, par les récepteurs 20. Chaque récepteur peut comprendre un hydrophone pour enregistrer la composante P et un géophone pour enregistrer la composante Z. Ces données comprennent un bruit et des multiples. A l'étape 202, diverses étapes de traitement peuvent être appliquées avant de séparer les ondes montantes et descendantes. Par exemple, les bulles peuvent être retirées des composantes P et Z (voir Wood et d'autres « The debubbling of marine source signatures », Geophysics, 1978), c'est-à-dire qu'un opérateur de retrait de bulles est appliqué aux données pour le retrait des oscillations des bulles des données. Les oscillations des bulles sont une conséquence de la nature de la source sismique, c'est-à-dire, des canons à air qui produisent une grande bulle qui oscille pendant une quantité de temps limitée. [0032] Ensuite, à l'étape 204, la séparation des champs d'ondes est appliquée aux données sismiques pour obtenir les ondes U montantes et D descendantes. La figure 5 montre la composante U et la figure 6 'montre la composante D. La séparation des champs d'ondes peut être effectuée, soit d'une manière classique, soit simultanément par l'application d'une atténuation de bruit initiale pour retirer partiellement un bruit tout en préservant un contenu de signal. [0033] Les champs d'ondes montant et descendant obtenus sont affectés par le bruit présent sur la composante Z. La séparation des champs d'ondes peut être effectuée juste au-dessus ou au-dessous du fond de l'océan. On doit noter que les champs d'ondes U et D sont tous deux affectés à la fois par le bruit et les multiples. En option, à l'étape 205, une extraction de signature directionnelle (voir G. A. Roberts et N. R. Goulty, 1990, « Directional deconvolution of marine seismic reflection data, North Sea example », Geophysical Prospecting) peut être appliquée à la composante P. Cette étape peut être exécutée en déconvoluant la composante P avec une version silencieuse du champ d'ondes D qui comprend des données de l'arrivée directe à juste au-dessus du multiple de fond marin du premier ordre uniquement, ou par des procédés classiques. [0034] A l'étape 206, une déconvolution haut-bas est effectuée en utilisant, par exemple, une estimation imparfaite des énergies montante et descendante (par exemple, les énergies correspondant aux composantes U et D) corrompues par le bruit ou par un bruit résiduel. Une manière d'effectuer une déconvolution haut-bas est décrite par L. Amundsen dans « Wavenumber-based filtering of marine point source data », Geophysics, vol. 58, n° 9, 1993, dont le contenu entier est inclus ici par voie de référence. [0035] En rappelant que l'énergie associée à la composante D peut être considérée comme l'énergie entrant dans la terre pendant l'étude sismique, et que R est la réponse de la terre à l'énergie correspondant à la composante D, l'énergie correspondant à la composante U peut être écrite en tant que U = D*R, où l'opération «* » est comprise comme une convolution moyenne. Le résultat de la déconvolution haut-bas R est ensuite défini par R = U/D, où «I » désigne ici l'opération de déconvolution. Le résultat de cette étape est illustré sur la figure 7. Cependant, il existe d'autres manières pour calculer R, comme cela sera reconnu par les hommes du métier. [0036] Ensuite, un modèle de multiples M est généré à l'étape 208. Il existe diverses manières de déterminer le modèle de multiples M. Par exemple, il est possible d'obtenir le modèle M sur la base d'une modélisation ou sur la base du modèle d'élimination de multiples associés à une surface (SRME). En variante, une manière originale de générer le modèle de multiples M est maintenant examinée. Selon un exemple de mode de réalisation, il est possible d'utiliser la déconvolution haut-bas R, la composante P et un terme de repositionnement pour générer le modèle de multiples M. Ainsi, dans ce cas, le modèle M est donné par F(R, P, AZ), où F est une fonction, et AZ concerne le repositionnement, et ce terme exprime une distance verticale entre l'emplacement où le modèle M est calculé et le fond de l'océan. Une définition concrète du modèle M est donnée, par exemple, dans le domaine f-k (fréquence- onde) comme étant M = R * P * On doit noter qu'une autre définition est vraie si un autre domaine est utilisé. Le modèle M peut être exprimé non seulement dans le domaine f-k, mais également dans d'autres domaines qui sont coutumier dans l'art, par exemple, tau-p. [0037] Le modèle de multiples M étant défini, l'algorithme original soustrait de manière adaptative, à l'étape 210, le modèle M de la composante de pression P. La composante P est choisie étant donné que cette composante peut être considérée comme étant exempte de bruit. Ainsi, l'étape 210 effectue en même temps une atténuation de bruit et une atténuation de multiples. Le résultat de cette étape est une estimation du champ de pression P - free exempt de multiples et de bruit côté source et récepteur, et ce résultat est illustré sur la figure 8. On doit noter que la composante P free illustrée sur la figure 8 est meilleure que la composante P d'origine illustrée sur la figure 3. [0038] Comme indiqué ci-dessus, l'étape 208 peut être exécutée dans un grand nombre de domaines. C'est également vrai pour l'étape 210 et pour d'autres étapes. Ainsi, les données sismiques à n'importe quelle étape peuvent être transformées dans ces domaines, soit sur une base rassemblement par rassemblement, soit localement dans des fenêtres glissantes sélectionnées par l'utilisateur. En outre, divers outils mathématiques peuvent être utilisés pour soustraire de manière adaptative le modèle M de la composante P. Il est connu dans le domaine du traitement de signal que la soustraction adaptative d'un modèle M d'un ensemble de données implique le calcul d'un ensemble d'opérateurs déduits des données qui font correspondre le modèle aux données en utilisant un ou plusieurs d'un grand nombre de critères. Par exemple, une norme L1 ou L2 peut être utilisée. Les données sismiques peuvent être transformées dans n'importe quel domaine de manière appropriée pour le problème traité, avec ou sans contraintes de rareté. La soustraction adaptative peut être effectuée en utilisant un grand nombre de critères de fenêtrage et un grand nombre de dimensions de filtre. Par exemple, des filtres unidimensionnels, bidimensionnels ou tridimensionnels peuvent être utilisés. Les paramètres de soustraction adaptative devraient être choisis de sorte que les multiples et le bruit soient atténués alors que l'énergie principale est préservée. [0039] D'autres mises en oeuvre du procédé proposé sont possibles. Par exemple, l'estimation R du champ d'ondes principal peut être obtenue en utilisant un procédé intégral tel que proposé par Amundsen (Amundsen « Elimination of free-surface related multiples without need of a source wavelet », Geophysics, 2001) au lieu d'une déconvolution haut-bas. En plus, un modèle de multiples M peut être calculé en utilisant un processus de repositionnement plus complexe que celui spécifié ci-dessus, pour prendre en compte, par exemple, des variations de la vitesse de l'eau et une topographie du fond de la mer. Divers procédés de repositionnement sont examinés dans la littérature (voir, par exemple, Berryhill, « Submarine Canyons-Velocity replacement by waveequation datuming before stack », Geophysics, 1986). [0040] Un concept original du procédé proposé est d'utiliser à la fois les enregistrements de composante de pression P et de composante verticale Z pour déduire le modèle de multiples M, et ensuite de soustraire de manière adaptative le modèle M obtenu de la composante de pression P pour des suppressions de multiples et de bruit simultanées en supposant que la composante P est un signal de bonne qualité, tandis que la composante verticale est corrompue par le bruit. [0041] Sur la base de la pression Pfree calculée, une image finale de la sous-surface étudiée peut être générée à l'étape 212. Bien entendu, la pression Pfree calculée peut être utilisée à d'autres fins pendant le traitement des données sismiques. [0042] Un aspect original de l'algorithme examiné ci-dessus est qu'il évite les efforts prenant du temps que nécessitent les procédés classiques pour une atténuation de bruit avant séparation. Autrement dit, si un bruit est atténué avant séparation, l'estimation de la qualité de l'atténuation du bruit nécessite l'exécution d'une séparation des champs d'ondes, d'une déconvolution haut-bas et d'un empilage. Plusieurs tentatives de séparation, de déconvolution haut-bas et d'empilage peuvent être nécessaires pour évaluer un ensemble d'autres paramètres dans des scénarios difficiles. [0043] En utilisant le procédé original, une séparation et une déconvolution haut-bas peuvent être exécutées en appliquant une atténuation de bruit initiale, alors qu'un retrait de bruit précis se produit après séparation et après déconvolution haut-bas, lorsque l'estimation des paramètres de soustraction adaptative est plus facile parce que les données sont prêtes pour un empilage. Comme examiné ci-dessus, le procédé original peut également être appliqué à des données de flûte. En outre, bien que le procédé suppose qu'un bruit n'est pas présent sur la composante d'hydrophone P alors que le bruit corrompt la composante de géophone verticale Z, le procédé peut également être modifié pour inclure le cas dans lequel un bruit est présent dans la composante P et pas dans la composante Z. [0044] Un algorithme plus complexe pour retirer le bruit et atténuer les multiples conjointement est maintenant examiné en relation avec la figure 9. La figure 9 montre schématiquement les étapes à effectuer pour obtenir la composante de pression P - free- On doit noter que la majeure partie de ces étapes ont déjà été examinées ci-dessus et, ainsi, elles sont simplement mentionnées dans la description qui suit. [0045] L'algorithme original 900 comprend une partie 902 qui détermine les composantes montante et descendante U et D à partir des données sismiques d'entrée. Ces composantes sont fournies à un dispositif informatique (qui sera examiné ultérieurement) pour exécuter la déconvolution haut-bas R dans une partie 904 de l'algorithme original. La déconvolution haut-bas R résultante est fournie, avec la composante de pression P (ou la composante de pression dont les bulles sont retirées) à une partie 906 pour déterminer le modèle de multiples M. Une partie 908 de l'algorithme original reçoit la composante descendante D (par exemple, de l'instant zéro jusqu'à un multiple du premier ordre) et la composante de pression P (ou la composante de pression de laquelle les bulles ont été retirées) et calcule la composante P - design de laquelle la signature a été extraite de manière directionnelle. [0046] La composante P - design de laquelle la signature a été extraite de manière directionnelle et le modèle de multiples M sont tous deux fournis à une partie 910 de l'algorithme original qui effectue la soustraction adaptative du modèle M de la composante P design de laquelle la signature a été extraite de manière directionnelle. Le résultat de ce calcul est la composante de pression Pfree de laquelle les multiples, la signature et le bruit ont été retirés. On doit noter que les parties 902 à 910 indiquées ci-dessus peuvent être mises en oeuvre dans un logiciel, un matériel ou une combinaison de ceux-ci. Par exemple, chaque partie peut être mise en oeuvre dans un module correspondant d'un dispositif informatique ou dans un même processeur du dispositif informatique. [0047] Un exemple d'un dispositif informatique représentatif capable d'exécuter les opérations selon les exemples de modes de réalisation examinés ci-dessus est illustré sur la figure 10. Un matériel, un micrologiciel, un logiciel ou une combinaison de ceux-ci peut être utilisé pour exécuter les diverses étapes et opérations décrites ici. [0048] L'exemple de dispositif informatique 1000 approprié pour exécuter les activités décrites dans les exemples de modes de réalisation peut comprendre un serveur 1001. Un tel serveur 1001 peut comprendre une unité centrale (CPU) 1002 couplée à une mémoire vive (RAM) 1004 et à une mémoire à lecture seule (ROM) 1006. La mémoire ROM 1006 peut également être un support de mémorisation d'autres types pour mémoriser des programmes, telle qu'une mémoire ROM programmable (PROM), une PROM effaçable (EPROM), etc. Le processeur 1002 peut communiquer avec d'autres composants internes et externes par l'intermédiaire d'éléments de circuit d'entrée-sortie (E/S) 1008 et d'un système de bus 1010, pour fournir des signaux de commande et similaires. Le processeur 1002 effectue un grand nombre de fonctions comme cela est connu dans l'art, telles que dictées par des instructions de logiciel et/ou de micrologiciel. [0049] Le serveur 1001 peut également comprendre un ou plusieurs dispositifs de mémorisation de données, comprenant des lecteurs de disque dur 1012, des lecteurs de CD-ROM 1014, et un autre matériel capable de lire et/ou de mémoriser des informations, par exemple un DVD, etc. Dans un mode de réalisation, un logiciel pour exécuter les étapes examinées ci-dessus peut être mémorisé et distribué sur un CD-ROM ou un DVD 1016, un support amovible 1018 ou une autre forme de support capable de mémoriser de manière portable des informations. Ces supports de mémorisation peuvent être insérés dans des dispositifs tels que le lecteur de CD-ROM 1014, le lecteur de disque 1012, etc., et lus par ceux-ci. Le serveur 1001 peut être couplé à un afficheur 1020, qui peut être n'importe quel type d'afficheur ou d'écran de présentation connu, tel que des afficheurs LCD ou à DEL, des afficheurs à plasma, des tubes à rayons cathodiques (CRT), etc. Une interface d'entrée d'utilisateur 1022 est prévue, comprenant un ou plusieurs mécanismes d'interface utilisateur, tels qu'une souris, un clavier, un microphone, un pavé tactile, un écran tactile, un système de reconnaissance vocale, etc. [0050] Le serveur 1001 peut être couplé à d'autres dispositifs informatiques par l'intermédiaire d'un réseau. Le serveur peut être une partie d'une configuration de réseau plus grande comme dans un réseau global (GAN) tel qu'Internet 1028. [0051] Comme cela sera apprécié par un homme du métier, les exemples de modes de réalisation peuvent être mis en oeuvre dans un dispositif de communication sans fil, un réseau de télécommunication, en tant que procédé ou dans un produit-programme d'ordinateur. Par conséquent, les exemples de modes de réalisation peuvent prendre la forme d'un mode de réalisation entièrement matériel ou d'un mode de réalisation combinant des aspects matériels et logiciels. En outre, les exemples de modes de réalisation peuvent prendre la forme d'un produit-programme d'ordinateur mémorisé sur un support de mémorisation pouvant être lu par un ordinateur comportant des instructions pouvant être lues par un ordinateur mises en oeuvre sur le support. N'importe quel support pouvant être lu par un ordinateur approprié peut être utilisé, comprenant des disques durs, des CD-ROM, des disques polyvalents numériques (DVD), des dispositifs de mémorisation optiques, ou des dispositifs de mémorisation magnétiques, tels qu'un disque amovible ou une bande magnétique. D'autres exemples non limitatifs de supports pouvant être lus par un ordinateur comprennent des mémoires de type mémoire flash ou d'autres types connus de mémoires. [0052] Les exemples de modes de réalisation présentés fournissent un appareil et un procédé pour le retrait de bruit et la réduction de multiples, simultanément, dans des données sismiques. 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. [0053] Bien que les caractéristiques et é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 autres caractéristiques et éléments présentés ici. [0054] 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, comprenant la réalisation 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. REVENDICATIONS1. Procédé pour le retrait de bruit et l'atténuation de multiples simultanément, dans des données sismiques enregistrées par des récepteurs sismiques (20), le procédé comprenant : la réception des données sismiques, dans lequel les données sismiques comprennent une composante de pression P et une composante verticale Z; la séparation des données sismiques en des champs d'onde U montants et des champs d'ondes D descendants ; le calcul d'une déconvolution haut-bas R sur la base des champs d'ondes U montants et D descendants ; la génération d'un modèle de multiples M sur la base de la déconvolution haut-bas R; et la soustraction de manière adaptative du modèle de multiples M de la composante de pression P pour obtenir une composante de pression corrigée Pfree de laquelle le bruit et les multiples ont été retirés simultanément, dans lequel un multiple est un champ d'ondes qui se propage d'une source vers un récepteur et qui subit plusieurs réflexions. 2. Procédé selon la revendication 1, dans lequel la composante de pression P est mesurée par un hydrophone et la composante verticale Z est mesurée par un géophone ou un accéléromètre et les deux composantes comprennent un bruit et des multiples. 3. Procédé selon la revendication 1, dans lequel l'étape de calcul de la déconvolution haut-bas R comprend : la déconvolution des champs d'ondes U montants et D descendants de sorte que R = U/D. 4. Procédé selon la revendication 1, dans lequel l'étape de génération du modèle de multiples M comprend : la définition du modèle de multiples M, dans un domaine fréquence-nombre d'onde, en tant que produit de la déconvolution haut-bas R et de la composante de pression P. 5. Procédé selon la revendication 4, comprenant en outre : le repositionnement du produit de la déconvolution haut-bas R et de la composante de pression P par une exponentielle qui comprend un nombre d'onde k, et une distance verticale Az, dans lequel la distance verticale définit une distance entre un fond de l'océan et un point auquel le modèle de multiples M est calculé. 6. Procédé selon la revendication 1, comprenant en outre : l'application d'un opérateur de retrait de bulles à la composante de pression P et à la composante verticale Z avant de séparer les données sismiques ; et l'application d'une extraction de signature directionnelle à la composante de pression P avant de calculer la déconvolution haut-bas R.25 7. Procédé selon la revendication 6, dans lequel l'extraction de signature directionnelle comprend : la déconvolution de la composante de pression P avec une version silencieuse de la composante descendante D; et la soustraction de manière adaptative du modèle de multiples M de la composante de pression de laquelle la signature a été extraite pour obtenir la pression corrigée P free. 8. Dispositif informatique (1000) pour le retrait de bruit et l'atténuation de multiples, simultanément, dans des données sismiques enregistrées par des récepteurs sismiques (20), le dispositif informatique (1000) comprenant : une interface (1010) configurée pour recevoir les données sismiques, dans lequel les données sismiques comprennent une composante de pression P et une composante verticale Z ; et un processeur (1002) connecté à l'interface (1010) et configuré pour : séparer les données sismiques en des champs d'onde U montants et des champs d'ondes D descendants, calculer une déconvolution haut-bas R sur la base des champs d'ondes U montants et D descendants, générer un modèle de multiples M sur la base de la déconvolution haut- bas R, et soustraire de manière adaptative le modèle de multiples M de la composante de pression P pour obtenir une composante de pression corrigée Pfree de laquelle le bruit et les multiples ont été retirés simultanément,dans lequel un multiple est un champ d'ondes qui se propage d'une source vers un récepteur et qui subit plusieurs réflexions. 9. Dispositif informatique selon la revendication 8, dans lequel le processeur est en outre configuré pour : déconvoluer les champs d'ondes U montants et D descendants de sorte que R = U/D ; définir le modèle de multiples M, dans un domaine fréquence-nombre d'onde, en tant que produit de la déconvolution haut-bas R et de la composante de pression P; repositionner le produit de la déconvolution haut-bas R et de la composante de pression P par une exponentielle qui comprend un nombre d'onde kz et une distance verticale Az, dans lequel la distance verticale définit une distance entre un fond de l'océan et un point auquel le modèle de multiples M est calculé ; appliquer un opérateur de retrait de bulles à la composante de pression P et à la composante verticale Z avant de séparer les données sismiques ; appliquer une extraction de signature directionnelle à la composante de pression P avant de calculer la déconvolution haut-bas R; déconvoluer la composante de pression P avec une version silencieuse de la composante descendante D ; et soustraire de manière adaptative le modèle de multiples M de la composante de pression de laquelle la signature a été extraite pour obtenir la pression corrigée P - free. , 2990519 24 10. Support pouvant être lu par un ordinateur configuré pour mémoriser en permanence des instructions d'ordinateur, de manière à, lorsqu'un ordinateur exécute les instructions, retirer le bruit et atténuer les multiples, simultanément, des données sismiques enregistrées par des récepteurs sismiques (20), les 5 instructions comprenant : la réception des données sismiques, dans lequel les données sismiques comprennent une composante de pression P et une composante verticale Z; la séparation des données sismiques en des champs d'onde U montants et des champs d'ondes D descendants ; 10 le calcul d'une déconvolution haut-bas R sur la base des champs d'ondes U montants et D descendants ; la génération d'un modèle de multiples M sur la base de la déconvolution haut-bas R; et la soustraction de manière adaptative du modèle de multiples M de la 15 composante de pression P pour obtenir une composante de pression corrigée Pfree de laquelle le bruit et les multiples ont été retirés simultanément, dans lequel un multiple est un champ d'ondes qui se propage d'une source vers un récepteur et qui subit plusieurs réflexions.