FR 2792080 A1 20001013 FR 0004147 A 20000331 L'invention concerne de façon générale le traitement de données sismiques, et plus particulièrement un procédé de compensation des effets d'un échantillonnage de données spatiales irrégulières de champs d'ondes sismiques et d'une illumination irrégulière de réflecteurs sismiques souterrains par des champs d'ondes sismiques. Dans l'exploration sismique, un champ d'onde acoustique (un tir) est généré par une source acoustique. Le champ d'onde se propage à travers la terre depuis la position d'une source. Le champ d'onde est réfléchi par des couches de terrain situées au-dessous de la surface, d'o il revient vers la surface. Plusieurs détecteurs sismiques sont répartis sur, ou à proximité de, la surface de la terre, à distance de la position de la source, suivant des lignes de prospection ou dans de grandes zones de surface. Les détecteurs, qui constituent des postes discrets d'échantillonnage du champ d'onde, sont de préférence répartis uniformément dans l'espace afin que le champ d'onde puisse être échantillonné uniformément à la fois en surface et dans le temps. Les données échantillonnées sont quantifiées et archivées en vue d'un traitement ultérieur. Le mouvement, ou la "migration", des données à partir de la position o elles sont enregistrées jusqu'à la position d'o partent les réflexions en profondeur dans le sol, est fondamental pour la plupart des procédés de formation d'images de réflexions en profondeur. La migration de Kirchhoff est actuellement l'algorithme de migration le plus communément utilisé pour une migration en profondeur avant sommation en trois dimensions. L'intégrale de Kirchhoff est basée sur une approximation de fréquence élevée du champ d'onde et est basée sur des calculs de temps de parcours associés. Pour cette raison, elle ne permet pas de traiter convenablement les situations les plus complexes. Dans tous les cas, dans les applications de l'art antérieur, l'intégrale est évaluée comme si le champ d'onde était échantillonné de façon régulière. En réalité, tel est loin d'être le cas pour un certain nombre de raisons. Premièrement, l'échantillonnage des données à la surface n'est pas uniforme du fait d'obstacles tels que des bâtiments, des routes ou autres cultures. Dans le cas d'une exploration marine, la répartition des positions des récepteurs peut être irrégulière du fait d'erreurs dans les positions prises par les détecteurs par suite d'une dérive de câble sous l'effet de courants et du vent ou par suite de la présence de plates-formes de forage et de production. Souvent, non seulement les données sismiques sont localement sous-échantillonnées, mais elles peuvent également être échantillonnées localement de façon excessivement dense. Un second problème avec l'évaluation de l'intégrale de Kirchhoff est que lorsque la vitesse de propagation sous la surface n'est pas uniforme, comme c'est presque toujours le cas, l'illumination des réflecteurs par le champ d'onde sismique n'est pas non plus uniforme, alors que, lorsque l'intégrale de Kirchhoff est évaluée, l'illumination des réflecteurs est habituellement supposée comme étant uniforme. Des procédés de traitement sismique bien connus tels que l'addition, le filtrage multicanal, la correction du delta-t dû au pendage (DMO), la migration avant sommation, les analyses de vitesse, les études d'anisotropie, la migration et l'extrapolation du champ d'onde supposent tous que les données sont échantillonnées uniformément. Cependant, comme indiqué ci-dessus, les données regroupées peuvent être échantillonnées de façon irrégulière. Cette irrégularité peut être due à des obstacles comme expliqué précédemment ou à des tirs manqués ou à des détecteurs ou récepteurs inopérants. Lorsque ces données résultant d'un échantillonnage spatial irrégulier ou impropre ne sont pas corrigées, il peut en résulter des artefacts de calcul non voulus qui se superposent aux données de sortie traitées. Des sous-programmes de traitement d'équations d'ondes tels que des échantillons DMO et de mouvements de données, par exemple par repositionnement du fait d'une migration de pendage, peuvent être affectés par un échantillonnage spatial irrégulier. En fait, l'algorithme devant améliorer les données de sortie dégrade en réalité ces données du5 fait qu'un sous-échantillonnage ou un sur- échantillonnage spatial laisse des restes des opérateurs de traitement de données dans les données de sortie. Un sous-échantillonnage grossier peut évidemment aboutir à un repliement du spectre de pendage. Le bruit, qui dépend du déport horizontal, n'est que partiellement éliminé dans l'opération d'addition en présence d'un échantillonnage spatial localement clairsemé. Le terme "opérateur" sera souvent utilisé dans cette description. Il est défini comme signifiant un élément spécifique impliqué dans une opération de traitement de données. Ainsi, un opérateur DMO est une expression spécifique utilisée dans l'application d'une correction à un delta-t normal dû au pendage. Un opérateur peut être exprimé sous la forme d'un symbole indiquant une opération devant être effectuée et lui-même peut être l'objet d'une manipulation mathématique. Divers auteurs ont abordé le problème de l'échantillonnage d'un champ d'onde spatiale localement clairsemé. Dans un article intitulé "Wave-equation Trace Interpolation", (Geophysics, vol. 52, no. 7, juillet 1987, pages 973 à 984) J. Ronen indique qu'une séquence de traitement dans laquelle on traite des données manquantes comme étant des données zéro et on effectue une migration partielle avant l'addition est équivalente à l'application de la transposée de l'opérateur d'addition qui doit en réalité être inversé. Ronen indique que l'inverse de l'opérateur ne peut pas être déterminé de façon unique, mais qu'il peut être estimé en utilisant un équilibrage spectral spatial dans un schéma itératif de gradients conjugués. La première itération est simplement une migration partielle avant sommation. Dans le cas o un repliement de spectre spatial est présent, plusieurs itérations additionnelles sont nécessaires. R. G. Williams et collaborateurs, dans un article présenté au 5lème congrès annuel de EAEG, mai 1989, intitulé "Model-constrained Anti- alias Filtering for Improved DMO", applique un filtre anti- repliement aux données afin que le sous-échantillonnage azimutal de l'opérateur DMO puisse être réduit par l'utilisation d'un modèle de pendage pour la prospection de façon à limiter l'ouverture de l'opérateur d'une manière dépendant de l'azimut. Dans un article publié sous la forme de l'abrégé développé n 1144, au 59ème congrès international de la SEG, 1989, intitulé "Effect of Irregular Sampling on Prestack DMO", J. Black et collaborateurs expliquent que des distributions irrégulières du point milieu et de l'azimut des données sismiques engendrent des artefacts dans le signal de sortie DMO si ces distributions irrégulières sont ignorées. Il est affirmé que ceci pose un problème en particulier pour la correction DMO avant sommation avec des géométries de collecte de données terrestres. En utilisant la géométrie réelle de collecte de données, la réponse DMO peut être calculée pour des événements plats et peut être utilisée pour concevoir des corrections de grandes longueurs d'ondes et des éditions sélectives qui minimisent l'impact des artefacts d'origine DMO. Ce processus peut être utilisé avec toute mise en oeuvre de la correction DMO, mais il est particulièrement approprié pour une correction DMO en trois dimensions utilisant le procédé de Kirchhoff. Le brevet des États-Unis d'Amérique n 5 206 837 délivré à Beasley et collaborateurs décrit un procédé basé sur la décomposition d'un regroupement d'un champ d'onde sismique en ses composantes constitutives de pendage et de déport horizontal source-détecteur, qui tient compte de l'effet d'un échantillonnage irrégulier des données spatiales en appliquant un opérateur d'échantillonnage inverse au regroupement. Le procédé peut être appliqué à tout algorithme DMO à deux dimensions ou trois dimensions dans une opération de regroupement de champ d'onde multicanal. Dans le brevet de Beasley et collaborateurs, un champ d'onde se propage depuis la position d'une source. Le champ d'onde est échantillonné dans l'espace en plusieurs postes discrets d'échantillonnage de données qui sont répartis sur une zone éloignée de la position de la source. Les échantillons du champ d'onde sont combinés en un regroupement brut auquel un opérateur de traitement choisi est appliqué dans une étape de traitement pour produire un regroupement traité. Un opérateur d'échantillonnage inverse est ensuite appliqué au regroupement traité pour produire un regroupement de champ d'onde égalisé. L'opérateur d'échantillonnage inverse comprend les étapes consistant à décomposer le regroupement traité en composantes constitutives de pendage. Le mouvement et la mise à l'échelle de caractères choisis des composantes, résultant de l'étape de traitement, sont analysés. Le mouvement et la mise à l'échelle de ces composantes font l'objet d'une correction portant sur des artefacts de traitement indésirés dus à un échantillonnage spatial localement, relativement clairsemé, ou relativement dense. Selon un aspect de cette invention, l'opérateur de traitement peut être une forme de DMO. Bien que le procédé d'égalisation dans le brevet de Beasley puisse être appliqué à tout régime de traitement de données multicanal tel que des études d'anisotropie, des opérateurs de filtrage, une migration et une extrapolation de champ d'onde, il n'y a pas d'exemples explicites de l'utilisation du procédé d'égalisation pour compenser l'échantillonnage irrégulier du champ d'onde aux interfaces réfléchissantes par suite d'irrégularités de vitesse. La présente invention montre comment le procédé d'égalisation décrit dans le brevet de Beasley précité peut être utilisé pour corriger un échantillonnage irrégulier du champ d'onde sismique par suite d'irrégularités de vitesse, donnant ainsi une migration de Kirchhoff améliorée dans laquelle les amplitudes sont mieux préservées. Les particularités nouvelles considérées comme étant caractéristiques de l'invention, à la fois en ce qui concerne l'organisation et des procédés de mise en ouvre, ainsi que leurs objets et avantages, ressortiront mieux de la description détaillée suivante et des dessins annexés à titre nullement limitatif et sur lesquels: la figure lA est un schéma illustrant un opérateur de migration dans des milieux à vitesse variable; la figure lB est un schéma illustrant la cartographie d'une zone sur le plan d'enregistrement en une zone d'angle solide dans la sphère unité au point image; la figure 2A est un organigramme illustrant une exécution par ordinateur du procédé de l'invention lorsqu'il est utilisé dans l'application d'un opérateur de migration avant sommation aux données; la figure 2B est un organigramme illustrant une exécution par ordinateur du procédé de l'invention lorsqu'il est utilisé dans l'examen de caractères d'illumination de position du signal de sortie du sous-sol et de problèmes potentiels basés sur une géométrie d'acquisition source- récepteur avec une structure de vitesse connue; la figure 2C est un organigramme illustrant une exécution par ordinateur du procédé de l'invention lorsqu'il est utilisé dans l'application d'un filtre basé sur l'invention avec un opérateur de traitement tel qu'un opérateur de migration; les figures 3A et 3B montrent une comparaison de collectes d'images migrées en profondeur avant sommation pour une géométrie de couverture souterraine irrégulière sans égalisation (figure 3A) et avec égalisation (figure 3B) des amplitudes en utilisant l'invention; la figure 4A illustre une collecte de données montrant des amplitudes sans égalisation au voisinage d'une masse salée; et la figure 4B illustre le résultat d'une égalisation des amplitudes portant sur la migration de données à partir de la figure 4A au voisinage d'une masse salée. Dans un souci de clarté et d'explication, le procédé de l'invention sera décrit à titre d'exemple, nullement limitatif, portant sur une migration avant sommation dans trois dimensions de champs d'ondes sismiques. Il est bien entendu que le procédé préféré peut être appliqué à tout régime multicanal de traitement de données tel que des études d'anisotropie, des opérateurs de filtrage, une extrapolation de champ d'onde et un traitement DMO comme mentionné précédemment. L'efficacité de la formation d'images sismiques se manifeste fondamentalement dans l'interaction de l'échantillonnage d'acquisition et de la distorsion du champ d'onde pendant que le son se propage à travers la terre. Pour qu'une expérience sismique forme une image satisfaisante d'un réflecteur situé en profondeur dans le sol, le son revenant par réflexion du réflecteur doit arriver avec une amplitude appréciable au récepteur. Dans l'approximation de haute fréquence, on peut considérer que la propagation et la réflexion du son sont approximées grossièrement par des rayons qui se propagent dans un domaine de vitesse complexe, et se réfléchissent en un point de formation d'image sur le réflecteur conformément à la loi de Snell. Intuitivement, on ne peut former l'image du réflecteur que s'il y a des paires source-récepteur dans la géométrie d'acquisition pour lesquelles la loi de Snell au point de réflexion peut être satisfaite. On imagine à présent comment ceci se manifeste de lui-même dans un algorithme simple de formation d'image de Kirchhoff qui travaille simplement en étalant et superposant des échantillons le long d'isochrones de temps de parcours (c'est-à-dire des opérateurs de migration). Une interférence constructive apparait lorsque la s tangente à l'opérateur de migration s'adapte au pendage du réflecteur au point de formation d'image, et une interférence destructive apparaît autrement, ce qui aboutit à une image du réflecteur. Si la tangente correspond au pendage, des rayons allant de la source au point image au récepteur pour cet opérateur satisfont essentiellement à la loi de Snell au point image. Dans une approximation acoustique simple qui néglige les pertes par transmission, trois effets s'associent à présent pour produire l'amplitude d'image finale dans le processus de sommation. Premièrement, le son a traversé un domaine de vitesse complexe, et son amplitude a donc chuté conformément à l'étalement géométrique. Un simple empâtement d'échantillon le long d'une isochrone suivi d'une somme peut éliminer une grande partie de cet effet, mais non nécessairement d'une manière qui est compatible avec une image réelle de réflectivité. Deuxièmement, si de nombreuses paires source- récepteur dans la géométrie d'acquisition satisfont à la loi de Snell au point de réflexion, on peut s'attendre à une grande amplitude par suite du processus d'interférence constructive. Inversement, si seulement peu de paires satisfont à la loi, on peut s'attendre à une chute d'amplitude. Il en résulte une variation d'amplitude dans l'image du réflecteur, qui est indépendante de celle due à des variations de la réflectivité: ceci résulte de la façon dont l'échantillonnage dans la géométrie d'acquisition se traduit en un échantillonnage de pendage au réflecteur en passant par des chemins de rayons complexes. Enfin, troisièmement, il s'agit de la réflectivité elle- même qui a un effet évident sur l'amplitude sismique. La façon d'éliminer l'interaction entre l'étalement géométrique et l'échantillonnage effectif au point image pour obtenir seulement la réflectivité est loin d'être évidente, mais est un élément essentiel du maintien d'une amplitude fidèle dans la migration de Kirchhoff. Récemment, des progrès importants ont été réalisés dans la dérivation de formules explicites d'inversion de Kirchhoff pour résoudre le problème de la réflectivité (Jaramillo, 1998, Bleistein, 1987, et les références citées). Les formules d'inversion qui sont actuellement les plus intéressantes dans la formation d'images de profondeur avant sommation à grande échelle sont celles qui tiennent convenablement compte des facteurs d'amplitude et de phase, mais négligent les pertes par transmission et la conversion de mode. Une telle inversion laisse intacte la réflectivité dépendante de l'angle, et est habituellement appelée inversion acoustique d'amplitude vraie. On voit le plus aisément les effets de la discrétisation en transformant les intégrales d'inversion, qui sont exprimées à l'origine sous forme d'intégrales sur une configuration source- récepteur à la surface, en intégrales sur le pendage de l'opérateur au point de formation d'image. De cette manière, la géométrie d'acquisition cartographie naturellement le pendage de l'opérateur dans la position de formation d'image sensiblement de la même manière que celle décrite dans le simple exemple de sommation de Kirchhoff ci-dessus. Des progrès notables sur l'effet de l'échantillonnage du pendage ont été réalisés dans le problème de la production DMO, et beaucoup s'appliquent au problème de l'inversion. Une technique qui a amélioré les résultats dans le cas du DMO a consisté à renormaliser les amplitudes par la distribution de pendages participants dans la position du signal de sortie, une technique appelée égalisation. Ceci a été décrit dans le brevet n 5 206 837 de Beasley, précité. On peut obtenir la réflectivité dans une inversion d'amplitude vraie en évaluant une intégrale de la forme (1) r(x)= If(, x)h(ó,x)d2 Ici, l'intégration s'effectue sur une configuration source- récepteur fixe à la surface, (3) h(; x) =6(.X Ip(x)' et hb(,x) est le déterminant de Beylkin, as,r sont les amplitudes des fonctions de Green de la source et du récepteur, et - a, (ó,x)t(), +r est la dérivée des données soumises à une convolution avec la fonction delta unilatérale, l'ensemble évalué à la somme des temps de parcours depuis la source jusqu'au point de sortie et du récepteur jusqu'au point de sortie, et p(x) est la somme des vecteurs lenteur le long des rayons de la source et du récepteur dans la position de sortie. Bien que cette équation semble impressionnante, un grand nombre des termes ont des interprétations géométriques. Les amplitudes des fonctions de Green peuvent être interprétées comme étant la racine carrée du rapport d'éléments de petite aire sur le front d'onde, l'un à proximité de la source et le second après que le front d'onde s'est déplacé jusqu'à la position de sortie. Similairement, le terme dérivé par rapport au temps impliquant les données est essentiellement le filtre rho connu appliqué aux données, qui est évalué ensuite à la somme des temps de parcours dans le processus de sommation de Kirchhoff. L'un des facteurs les plus intéressants est le facteur h impliquant le déterminant de Beylkin. Ce facteur est conceptuellement similaire à un rapport de petites aires de la manière suivante. On considère un point dans la configuration source- récepteur. Ceci peut être considéré comme définissant une trace sismique avec une position fixe de la source et du récepteur à la surface. Étant donné une position de sortie x, il peut exister un opérateur de migration, c'est-à-dire une courbe de temps de parcours source-récepteur constant, qui passe par le point de sortie. Si cet opérateur de migration existe, la normale à l'opérateur dans la position de sortie définit alors une direction qui est essentiellement le pendage de l'opérateur dans cette position. Étant donné que le pendage consiste en deux angles dans un espace à trois dimensions, il peut être considéré comme définissant un point v sur une sphère unité. Lorsque le point % varie, le pendage de l'opérateur passant par le point x varie aussi, en sorte qu'une petite aire dans le plan de surface autour du point est représentée par une petite aire sur la sphère unité autour du point v. Le facteur h de Beylkin est essentiellement le rapport de ces deux aires; c'est le jacobien de l'application de % à v. Ceci est illustré schématiquement sur la figure lA pour un milieu à vitesse variable. Une trajectoire sismique 12 est représentée en coupe depuis un point de tir 10 sur la surface de la terre 5, laquelle, après réflexion à partir d'un point 32 dans la profondeur du sol, atteint la position 20 d'un récepteur comme illustré par la trajectoire sismique 12'. L'opérateur de migration 30 est une courbe (surface dans trois dimensions) qui est le lieu géométrique de la totalité de ces trajectoires 12 - 12'. Il s'agit également d'une courbe (surface) à temps de parcours source-récepteur constant. La figure 1A représente également une sphère unité 34 centrée au point de réflexion 32 et une normale 36 à l'opérateur de migration 30 au point de réflexion. Le signal de sortie après migration de Kirchhoff est obtenu dans la position 32 de sortie en intégrant les contributions provenant de la totalité des combinaisons possibles de tir et de récepteurs. Dans un domaine échantillonné spatialement, cette intégrale est approximée par sommation sur un nombre fini de paires de tir et de récepteur. Étant donné plusieurs combinaisons de tir 10 et de positions de récepteurs 20, le signal de sortie migré au point 32 est obtenu par sommation des données provenant de la totalité de ces combinaisons de tir et de récepteurs. Pour exécuter convenablement une telle migration, y compris une préservation des amplitudes du point de réflexion, l'intégration doit être effectuée sur la sphère unité 34 centrée dans la position de sortie. En pratique, du fait d'une ouverture finie des positions de tir et de récepteurs sur le plan d'enregistrement 40 (voir figure lB), il est évidemment impossible de "couvrir" la sphère unité. De plus, en présence d'espaces dans la couverture, avec des sources ou des récepteurs, la couverture de la sphère unité 34 n'est pas uniforme. Cet effet a été convenablement traité dans le brevet n 5 206 837 précité. La figure lB montre également une autre source d'erreur due à la courbure des trajectoires sismiques 12 et 12' (de la figure 1A) lorsque la distribution de vitesse dans le sous- sol n'est pas uniforme. L'effet en est que des points se trouvant dans un petit élément d'aire 42 sur le plan d'acquisition s'appliquent de façon non uniforme à une région 44 sur la sphère unité. Une migration préservant l'amplitude nécessite que l'intégrale de Kirchhoff soit obtenue sur des éléments uniformes sur la sphère unité 34. Ceci peut être réalisé en définissant de petits éléments sur la sphère unité qui sous-tendent des angles solides égaux au point 34 de réflexion. Par conséquent, une migration préservant l'amplitude nécessite que des données collectées sur des aires sur le plan d'acquisition 40 qui correspondent à des angles solides égaux sur la sphère unité 34 soient affectées de poids égaux. Ceci est contraire aux procédés de l'art antérieur dans lesquels l'intégration est effectuée sur de petites aires sur le plan d'acquisition ayant des dimensions Ax et Ay dans les directions x et y et les contributions de ces petites aires font l'objet d'une sommation. Etant donné cette représentation graphique, l'équation de la réflectivité (1) peut être écrite sous la forme r(x)= Jf(x,V)d2V (4) o le terme h est absent en comparaison avec l'équation (1). Une discrétisation de l'équation (4) a pour résultat r(x)= _ f(ó,x)h(ó,x)A24 = f (V,x)A2V (5) i i En général, même si la géométrie d'acquisition est régulière, la distribution de points Vi sur la sphère unité n'est pas régulière, en particulier si le domaine des vitesses au-dessus du point de sortie est complexe. Un calcul de la sommation dans l'équation (5) nécessite de calculer les amplitudes des fonctions de Green qui font partie de f. Dans la présente invention, un calcul du facteur h de Beylkin n'est pas réalisé pour les raisons indiquées ci-dessous. Il existe des expressions analytiques pour le facteur h de Beylkin en termes d'amplitudes de fonctions de Green pour la configuration de tir commune, mais une expression analytique simple exprimée par des amplitudes des fonctions de Green pour le cas d'un départ horizontal commun n'existe pas. Cependant, au lieu d'utiliser une expression analytique pour h, il est beaucoup plus simple en pratique d'utiliser le processus suivant pour évaluer le facteur. On suppose qu'on fixe une géométrie d'acquisition donnée. Si l'on divise le plan d'acquisition en cases suffisamment petites, on peut alors réaliser un comptage du nombre d'échantillons qui tombent dans une case donnée. En divisant le nombre compté par la dimension de l'échantillon, on obtient le nombre d'échantillons dans la géométrie d'acquisition par unité d'aire, ce qu'on désigne n%. Ensuite, étant donné une petite aire élémentaire A2, le nombre d'échantillons dans cette aire élémentaire est donné par N = niA2t. Dans l'application à la sphère unité, ces échantillons s'appliquent à l'aire élémentaire A2v, en sorte que le nombre d'échantillons par unité d'aire sur la sphère unité est nv = N/A2v. La réflectivité finale devient alors r(x)= _ _(,x)/n.(v) (6) Une particularité notable de l'équation (6) est qu'elle implique seulement f et le nombre d'échantillons qui s'appliquent à la sphère unité de pendage par aire pour la géométrie d'acquisition donnée; elle n'implique aucune aire différentielle sur la surface d'acquisition, et le facteur de Beylkin est calculé de façon implicite. Par conséquent, aucune expression analytique pour le facteur de Beylkin n'est nécessaire; une division d'amplitude par le nombre compté d'échantillons par aire sur la sphère unité suffit. La figure 2A montre l'application du procédé de l'invention à un opérateur de traitement de données tel qu'une migration de Kirchhoff avant sommation. Les échantillons d'entrée de champ d'onde sismique bruts sont regroupés en 100 et les positions source-récepteurs sont identifiées en 102 par rapport aux échantillons de champ d'onde respectifs pour déterminer des déports horizontaux source-récepteurs. Pour un modèle de vitesse connu du sous- sol, un opérateur de migration est déterminé en 104 pour chaque point de sortie pour la géométrie appropriée de source et récepteurs. Comme indiqué précédemment, le procédé de l'invention est indépendant de la géométrie tir- récepteur, en sorte que les données peuvent être dans un tir commun, dans un récepteur commun, dans un déport horizontal commun ou dans tout autre regroupement convenable. Pour chaque point de sortie, l'échantillonnage des données est analysé en termes d'échantillons par angle solide unité 106. A des fins de calcul, une approximation est réalisée dans la détermination de la densité d'échantillons par angle solide unité: premièrement, la sphère unité est divisée en segments sur la base d'un échantillonnage en pendage et en azimut et on calcule une densité pour ces segments. Ces valeurs sont ensuite corrigées par un facteur trigonométrique approprié pour transformer les valeurs de densité en une densité par angle solide unité. Un opérateur d'échantillonnage inverse correspondant est déterminé en 108. Une opération de conditionnement est appliquée pour éliminer toutessingularités, c'est-à-dire des valeurs nulles dues à un échantillonnage spatial localement incorrect. L'opérateur de migration est appliqué avec l'opérateur d'échantillonnage inverse aux données sismiques 112. Ceci produit des données migrées avec une correction d'amplitude pour compenser l'effet d'une illumination non uniforme. Une opération de sommation facultative est appliquée en 114 et la section migrée est produite. Cette sommation facultative est effectuée si un nombre d'images partielles est produit par le traitement des données en 112 et améliore le rapport signal/bruit de la section de sortie. L'image de réflecteurs dans le sous-sol peut être formée du fait que de l'énergie sismique en est revenue par réflexion. Dans certains cas, l'image formée d'un emplacement dans le sous-sol peut être mauvaise du fait que la géométrie d'acquisition de position source- récepteurs et la structure de vitesse dans le sous-sol sont telles qu'il ne revient par réflexion que peu ou aucune énergie sismique d'une position de sortie. Ce problème est souvent appelé problème d'illumination, c'est-à-dire que l'emplacement dans le sous-sol (position de sortie) est faiblement illuminé. Pour comprendre de tels problèmes, on peut évaluer les caractères d'illumination en considérant une géométrie d'acquisition de positions source-récepteurs et un modèle de vitesse. Par exemple, le caractère d'illumination d'une densité d'échantillons par angle solide unité peut être affiché pour déterminer quelles sont les positions de sortie qui risquent d'être sujettes à des problèmes de bruit dus à des variations avec le pendage du réflecteur dans l'illumination. Normalement, une analyse d'illumination était réalisée dans l'art antérieur par un lancer de rayons depuis la géométrie spécifique dans le sous-sol et par une analyse de la quantité de rayons revenant par réflexion de réflecteurs souterrains spécifiques. Cependant, bien que le lancer de rayons montre la densité de couverture en un réflecteur, il ne convient pas à une analyse de tous les réflecteurs possibles, de tous les pendages illuminés possibles et de toutes les positions de sortie, contrairement à la présente invention. En examinant l'illumination de tous les pendages dans des positions de sortie, on peut comprendre si l'on peut mettre en évidence un bruit cohérent au prix d'une structure souterraine réelle. Par conséquent, comprenant quelque peu la structure de vitesse souterraine, on peut évaluer des géométries d'acquisition et quels types de caractères d'illumination liés à l'échantillonnage et de problèmes de bruit potentiel vont ensuite apparaître pour toutes les positions de sortie. On peut utiliser le lancer de rayons à des fins d'analyse uniquement lorsque des réflecteurs coïncident avec des positions de sortie. La présente invention peut être utilisée en pratique pour une analyse d'illumination pour tous les caractères de position de sortie dépendant du pendage dans le sous-sol. Toute sensibilité des positions de sortie du sous-sol aux problèmes de bruit liés à l'échantillonnage apparaît souvent à l'examen des caractéristiques affichées d'illumination des caractères des positions de sortie dépendant du pendage, telles qu'une densité d'échantillons par angle solide unité en des points de sortie. Cette information de densité d'échantillons peut être combinée pour être affichée avec d'autres caractères d'illumination liés à la géométrie d'acquisition et à la structure de vitesse. La figure 2B montre le procédé de l'invention appliqué à l'examen de la façon dont des réflecteurs, des structures et d'autres particularités du sous-sol sont illuminés. Pour la géométrie source-récepteurs 200 et un modèle de vitesse du sous-sol 201, on calcule en 202 un opérateur de traitement, par exemple un opérateur de migration de Kirchhoff. Ces caractéristiques d'illumination peuvent ensuite être calculées en chaque point de sortie 204 et affichées en 206. Les caractéristiques affichées montrent aux spécialistes de la technique si des particularités souterraines peuvent être bien ou mal illuminées et quels types de problèmes de bruit peuvent interférer avec l'illumination du sous-sol, et quels pendages peuvent poser des problèmes de bruit. L'application d'une migration de profondeur égalisée, comme indiqué sur la figure 2A, peut être considérée comme étant un mécanisme généralisé de filtrage dépendant du pendage commandé par l'illumination. En variante, cet aspect de filtrage de l'invention peut être dirigé ou polarisé pour renforcer ou supprimer des aspects d'illumination du sous-sol que le spécialiste souhaite. Ceci est obtenu en ajustant la pondération en une position5 de sortie du sous-sol pour un pendage spécifique afin de renforcer ou d'affaiblir l'amplitude à ce pendage dans l'image finale. Les pendages qui sont affectés peuvent varier d'une position de sortie à une autre, et peuvent dépendre de la façon dont le technicien souhaite polariser la formation d'images. Par exemple, en plus d'une égalisation des contributions des signaux à partir de tous les angles de pendage présents dans toute position de sortie, les données en une position de sortie de point image quelconque peuvent être polarisées sur la base de la structure de vitesse du sous-sol, de la structure géologique, de toutes caractéristiques d'énergie dominante ou exactement des angles de pendage que le spécialiste choisit de renforcer ou de supprimer. La figure 2C montre l'application du procédé de l'invention à un opérateur de traitement de données combiné à un filtre pour renforcer ou supprimer des aspects de données choisis. Les échantillons d'entrée de champ d'onde sismique bruts 300 sont regroupés et les positions source- récepteurs 302 sont identifiées par rapport aux échantillons de champ d'onde respectifs pour déterminer des déports horizontaux source-récepteurs. Pour un modèle de vitesse connu du sous- sol, un opérateur de traitement (par exemple un opérateur de migration de Kirchhoff) est déterminé en 304 pour chaque point de sortie pour les géométries source-récepteurs appropriées. Pour chaque point de sortie, l'échantillonnage des données est analysé en termes d'échantillons par angle solide unité 306 comme expliqué en référence à la figure 2A ci-dessus. Un opérateur de filtrage est ensuite choisi en 308, lequel renforce ou supprime des données à des pendages souhaités par le technicien. Une opération de conditionnement 310 est appliquée pour éliminer toutes singularités si le filtre ou la donnée n'est pas déjà bien conditionné. L'opérateur de traitement et cet opérateur de filtrage sont appliqués en 312 aux données sismiques. Ceci produit des données sismiques traitées avec des caractères d'illumination renforcés et/ou supprimés. Une sommation facultative des autres regroupements de données traitées 314 peut être réalisée, laquelle peut améliorer le signal rapport/bruit de la section traitée de sortie. Les données pour la figure 3A ont été acquises le long de plusieurs lignes de prospection sismique radiales et circonférentielles. Ces données ont une plage de déport horizontal de 900 à 960 mètres et la densité de la couverture superficielle est irrégulière. La figure 3A montre le résultat d'une migration de Kirchhoff des données acquises dans cette plage de déport horizontal en utilisant des procédés de l'art antérieur. On représente quatre regroupements migrés séparés à point-milieu commun (CMP). La figure 3B montre les mêmes quatre regroupements migrés CMP utilisant une migration de Kirchhoff avec la correction d'amplitude de la présente invention en suivant le procédé décrit en regard de la figure 2A. L'amélioration due à la migration apparaît sous la forme d'un comportement d'amplitude amélioré d'un déport horizontal à un autre et d'une image plus continue. La figure 4A montre des données utilisant une migration de Kirchhoff classique obtenues en un emplacement o il n'y a pas d'irrégularités importantes dans l'échantillonnage de surface des données. Cependant, la section géologique comprend ici une couche de sel qui produit de graves distorsions dans les champs d'ondes du fait d'un contraste de vitesse important avec des sédiments adjacents. La figure 4B montre la migration obtenue en utilisant la présente invention et en compensant l'illumination non uniforme du sous- sol par suite de variations de vitesse dans celui-ci. L'invention a été décrite dans une application particulière à une migration de Kirchhoff à titre d'exemple nullement limitatif de l'invention. Comme indiqué précédemment, l'invention est une généralisation du procédé5 d'égalisation décrit dans le brevet n 5 206 837 précité pour des opérations DMO (correction du delta-t dû au pendage). Les concepts développés ici peuvent être appliqués à tout procédé de traitement ou de regroupement de champ d'onde multicanal, par exemple à des opérateurs10 DMO plus exotiques conçus pour traiter des analyses de vitesse, une anisotropie, un filtrage et d'autres processus liés à une migration et à une extrapolation de champ d'onde. L'invention n'est en aucune manière limitée à une migration de Kirchhoff d'un type spécifique quelconque et15 elle n'est limitée que par le cadre des revendications annexées. REVENDICATIONS 1. Procédé de traitement de champs d'ondes sismiques propagés depuis une position de source (10), caractérisé en ce qu'il comprend:5 (a) l'échantillonnage spatial du champ d'onde en une pluralité préalablement choisie de postes d'échantillonnage discrets répartis sur une aire éloignée de la position de source; (b) la combinaison des échantillons de champ d'onde en un regroupement d'entrée de champs d'ondes bruts; (c) pour un opérateur de traitement choisi, la détermination d'un nombre d'échantillons par angle solide unité du regroupement traité en plusieurs positions de sortie (20) afin de donner plusieurs valeurs de densité d'échantillons, (d) la détermination d'un opérateur de mise à l'échelle inverse pour la pluralité de positions de sortie; et (e) l'application de l'opérateur de traitement choisi et de l'opérateur de mise à l'échelle inverse au regroupement d'entrée de champs d'ondes bruts et la production d'une section sismique à la pluralité de positions de sortie. 2. Procédé selon la revendication 1, caractérisé en ce que l'opérateur de traitement est un opérateur de migration de Kirchhoff. 3. Procédé selon la revendication 1, caractérisé en ce que la détermination de l'opérateur de mise à l'échelle inverse comprend en outre le renforcement ou la suppression de données sismiques en des positions de sortie pour des pendages prédéterminés associés aux valeurs de densité d'échantillons. 4. Procédé selon la revendication 3, caractérisé en ce que la suppression ou le renforcement des données sismiques en des positions de sortie pour des pendages prédéterminés est basé sur au moins l'un des points suivants: (i) une égalisation de données sismiques à tous les angles de pendage des positions de sortie; ou (ii) une égalisation de données sismiques à tous les angles de pendage des positions de sortie o des données sont présentes; ou (iii) une structure de vitesse du sous-sol; ou (iv) une structure géologique du sous-sol; ou (v) des caractéristiques de pendage à énergie dominante; ou (vi) des caractéristiques d'énergie de signal prédéterminées; ou (vii) des caractéristiques d'énergie de bruit prédéterminées. 5. Procédé de compensation de l'effet d'une illumination irrégulière de réflecteurs du sous-sol dans une prospection sismique multicanal, caractérisé en ce qu'il comprend: (a) la propagation d'un champ d'onde sismique depuis une position de source (10) en utilisant une source de ce champ d'onde; (b) l'échantillonnage spatial du champ d'onde en une pluralité préalablement choisie de postes d'échantillonnage discrets (20) répartis sur une aire éloignée de ladite position de la source; (c) la combinaison des échantillons de champ d'onde en un regroupement de champ d'onde brut; (d) le calcul, pour un opérateur de traitement sélectionné, d'un opérateur d'échantillonnage qui met en relation l'échantillonnage spatial à ladite pluralité discrète de postes d'échantillonnage avec un échantillonnage local associé à ladite illumination irrégulière en plusieurs points image; et (e) l'application d'un opérateur d'échantillonnage inverse et de l'opérateur de traitement sélectionné au regroupement de champ d'onde brut pour produire un regroupement traité dont l'illumination irrégulière fait l'objet d'une compensation. 6. Procédé selon la revendication 5, caractérisé en ce que l'opérateur de traitement est un opérateur de migration de Kirchhoff. 7. Procédé d'analyse des caractères d'illumination de points de sortie du sous-sol par rapport à une position de source dans une prospection sismique multicanal, caractérisé en ce qu'il comprend: (a) pour plusieurs postes d'échantillonnage discrets (20) sélectionnés, répartis à distance de ladite source (10) et un domaine de vitesse connu dans le sous-sol, la détermination d'au moins un caractère de position de sortie, dépendant du pendage, d'un regroupement de données devant être traitées en plusieurs positions de sortie, afin de donner plusieurs valeurs de position de sortie dépendant du pendage pour un opérateur de traitement choisi; (b) la dérivation de caractères d'illumination à partir des valeurs de position de sortie dépendant du pendage, de géométries d'acquisition des postes d'échantillonnage discrets et du domaine de vitesse du sous-sol; (c) la combinaison et l'affichage des caractères d'illumination en plusieurs positions de sortie. 8. Procédé selon la revendication 7, caractérisé en ce que le caractère de position de sortie dépendant du pendage est la densité d'échantillon par angle solide unité. 9. Procédé selon la revendication 7, caractérisé en ce que l'opérateur de traitement est un opérateur de migration de Kirchhoff. 10. Procédé de filtrage de données sismiques pour l'illumination d'une structure du sous-sol dans une prospection sismique multicanal, caractérisé en ce qu'il comprend: (a) la propagation d'un champ d'onde sismique à partir d'une position de source (10) en utilisant une source de celui-ci; (b) l'échantillonnage spatial du champ d'onde en une pluralité préalablement choisie de postes d'échantillonnage discrets (20) répartis sur une aire éloignée de la position de la source; (c) la combinaison des échantillons de champ d'onde dans un regroupement de champ d'onde brut; (d) le calcul, pour un opérateur de traitement choisi, d'un opérateur de filtrage et d'échantillonnage qui met en relation l'échantillonnage spatial par angle solide unité à la pluralité discrète de postes d'échantillonnage avec une fonction de polarisation prédéterminée, renforçant ou supprimant ainsi l'illumination à divers pendages; (e) l'application de l'opérateur de filtrage calculé et de l'opérateur de traitement choisi au regroupement de champ d'onde brut pour produire un regroupement traité qui renforce ou supprime les caractéristiques d'illumination. 11. Procédé selon la revendication 10, caractérisé en ce que le calcul de l'opérateur de filtrage comprend en outre le renforcement ou la suppression de données sismiques en des positions de sortie pour des pendages sur la base d'une fonction de polarisation prédéterminée basée sur au moins l'un des points suivants: (i) une égalisation de données sismiques à tous les angles de pendage des positions de sortie; ou (ii) une égalisation des données sismiques à tous les angles de pendage des positions de sortie o des données sont présentes; ou (iii) une structure de vitesse du sous-sol; ou (iv) une structure géologique du sous-sol; ou (v) des caractéristiques de pendage d'énergie dominante; ou (vi) des caractéristiques d'énergie de signal prédéterminée; ou (vii) des caractéristiques d'énergie de bruit prédéterminée. 12. Procédé selon la revendication 10, caractérisé en ce que l'opérateur de traitement est un opérateur de migration de Kirchhoff.