FR 2826734 A1 20030103 FR 0108703 A 20010629 PERFECTIONNEMENTS AUX TECHNIQUES DE POINTE DE SURFACE La présente invention concerne des perfectionnements aux techniques de pointé de surfaces utilisées en traitement sismique. Dans le domaine des sciences de la terre, pour imager la structure du sous-sol, le géophysicien utilise les techniques de sismique réflexion . Celles-ci consistent à émettre à la surface du sol des signaux acoustiques et à les enregistrer après leurs réflexions successives aux limites des couches géologiques. Ces couches constituent le sous-sol dont on cherche à créer un modèle géométrique. Pour ce faire, l'un des procédés le plus couramment utilisé, notamment dans la recherche de gisements pétroliers est le procédé d'acquisition 3D . Celui-ci consiste à répartir des émetteurs et des récepteurs dans un plan horizontal en surface le long d'une grille de coordonnées X, Y. Le résultat d'un tel procédé est une image tridimensionnelle du sous-sol organisée sous la forme d'enregistrements échantillonnés en fonction du temps ou de la profondeur, noté Z. Chacun de ces enregistrements, appelé trace , est un représentatif de la verticale à l'aplomb d'un point PI) de coordonnées X,, Yj. Une telle acquisition tridimensionnelle permet d'obtenir un volume dont les dimensions sont deux axes géométriques X, Y et un axe des temps ou de profondeur d'acquisition Z. Outres ces trois dimensions spatiales (X, Y et axe vertical Z), bien d'autres grandeurs ont un intérêt pour le physicien, et notamment : - La distance source/récepteur ou l'angle d'incidence qui a servi à générer le volume. Jusqu'à il y a peu cette distance (déport) était le plus souvent ramené à 0, les volumes étaient obtenus par sommation de contributions à distance constante de façon à augmenter le rapport signal/bruit. La sommation partielle par classes de déport ou d'angle permet désormais au géophysicien de faire des analyses par classe de distance ou d'angle. - Le temps. Plusieurs acquisitions sismiques peuvent être effectuées successivement dans le temps avec un même réseau de récepteurs/émetteurs ; la comparaison des données acquises permet, par exemple, de comprendre l'évolution d'un régime de fluide dans le cas d'une exploitation de réservoir. - Le type d'ondes sismiques. Les capteurs peuvent enregistrer simultanément les ondes de compression (ondes P) et les ondes de cisaillement (ondes S) ; le rapport entre les temps d'arrivée de ces deux types d'ondes est un lien direct avec la nature du milieu où les ondes se propagent. - La vitesse. La migration temps ou profondeur consiste à déplacer les points du volume d'acquisition en tenant compte d'un champ de vitesse. La vitesse utilisée dans ce processus peut être perturbée et constituer une dimension d'analyse complémentaire. Par exemple, pour mettre à jour un modèle de vitesse, la demanderesse a proposé une approche consistant à distribuer de manière régulière les facteurs de perturbation sur un modèle de vitesse, et de s'intéresser pour une position X, Y, à la valeur de l'échantillon associé à un même réflecteur dans chacun des cubes migrés. - etc. La prise en compte de ces autres grandeurs est fondamentale dans nombre de domaines des géo-sciences, de l'analyse de vitesse à la caractérisation de réservoir. Elle est aujourd'hui accessible grâce à l'augmentation de la puissance de calcul informatique et des capacités de stockage ; il reste à disposer des outils d'analyse et d'interprétation adaptés à cette nouvelle orientation. Aujourd'hui, lorsque l'on prend en compte avec les outils existants une grandeur en plus des trois dimensions spatiales dans lesquelles les volumes sismiques sont généralement représentés, l'interprétation des données se ramène à la somme des interprétations réalisées de façon indépendantes sur chacun des volumes sismiques 3D obtenus pour différentes valeurs de ladite grandeur. Les informations que l'on peut retirer d'un tel mode d'interprétation sont nécessairement limitées. L'invention propose quant à elle une contribution efficace et originale à la résolution de ce problème. Une technique classiquement utilisée pour l'analyse de volumes sismiques 3D est celle du pointé de surface, dont un exemple est donné dans le brevet FR 2 657 695. Lorsque l'on prend en compte une grandeur supplémentaire en plus des dimensions spatiales qui constituent les volumes 3D (X, Y et axe vertical) et que l'on fait varier ladite grandeur, une difficulté tient en ce que la nature et la forme du réflecteur à interpréter peut changer d'un volume à l'autre le long de la dimension que constitue ladite grandeur. La cohérence entre les différents volumes 3D n'est pas forcément inscrite dans l'échantillon cartographié, ce qui conduit à des problèmes de choix quant aux événements interprétés par le géophysicien. Cela peut entraîner des discontinuités et des aberrations dans le pointé de la surface. L'invention permet de pallier ces insuffisances. Elle permet des interprétations prenant en compte l'ensemble des paramètres à analyser et tenant compte des lois qui les caractérisent. L'outil qu'elle propose permet en outre de valider la cohérence des analyses réalisées par rapport aux lois géophysiques dont dépendent les grandeurs que l'on prend en compte. A cet effet, il est proposé un procédé de traitement de données sismiques, selon lequel, en vue de leur analyse en fonction d'au moins quatre dimensions, on met en oeuvre un pointé et/ou une propagation de pointé sur au moins un sous-ensemble de ces données qui correspond à une collection de ces données échantillonnées selon trois ou deux dimensions, la ou les grandeurs qui correspondent à la ou aux autres dimensions étant alors de valeurs fixes, puis on applique à au moins un des points mis en évidence par ce pointé et/ou cette propagation de pointés une loi de projection donnée afin de déterminer un point qui lui correspond dans un autre sous-ensemble correspondant à une collection de données échantillonnées selon les mêmes dimensions, pour une valeur différente d'au moins une grandeur qui correspond à au moins une autre dimension. Par loi de projection, on entend une loi destinée à relier les points d'un des sous-ensembles à ceux d'un autre, cette loi étant fonction des modifications que l'on apporte à la ou aux grandeurs correspondant aux dimensions autres que celles dans lesquelles lesdits sous-ensembles sont représentés. Des exemples de loi de projection sont donnés plus loin de façon détaillée. Avantageusement, on met en oeuvre une propagation de pointé à partir d'au moins un point projeté dans un tel autre sous-ensemble. Notamment, de façon préférée, on met en oeuvre une succession d'alternance de propagations de pointés et de projections. De façon avantageuse également, on applique à au moins un des points déterminé par propagation de pointé dans un tel autre sousensemble une loi de projection inverse afin de déterminer un point qui lui correspond dans au moins un sous-ensemble à partir duquel on a mis en oeuvre une projection. Egalement, on met préférentiellement en oeuvre sur un point projeté un test de cohérence pour le valider. En particulier, on affiche une cartographie des points pointés et/ou projetés et/ou déterminés par propagation de pointés dans les différents sous-ensembles considérés. Une grandeur que l'on fait varier pour définir une autre dimension est avantageusement un déport source/récepteur et/ou un angle d'incidence et/ou un temps d'acquisition et/ou un type d'ondes sismiques et/ou une vitesse de migration et/ou un paramètre prenant en compte au moins l'un de ces paramètres. De façon particulièrement préférée, le nombre de dimensions prises en compte est supérieur à quatre. L'invention concerne également un procédé de détermination cartographique sur des données sismiques, caractérisé en ce qu'on met en oeuvre le procédé précité pour déterminer une cartographie selon au, moins quatre dimensions. Avantageusement notamment, on détermine à partir d'un point projeté et/ou pointé et/ou déterminé par propagation de pointés un voisinage de points comportant en plus des points dans un voisinage prédéterminé dudit point dans le sous-ensemble de dimensions deux ou trois auquel ledit point appartient, un ou des points projetés à partir dudit point dans d'autres sous-ensembles et en ce qu'on détermine sur l'ensemble des points de ce voisinage celui de ces points qui optimise un critère donné. Autrement formulé, il est proposé un procédé d'interprétation des données sismiques dans une espace de dimension supérieure à trois, les trois premières dimensions étant avantageusement, mais non limitativement, les dimensions spatiales (X, Y et axe vertical). Une caractéristique de ce procédé est d'alterner des phases de propagation, c'est à dire de cartographie des surfaces dans l'espace 3D (X, Y et axe vertical par exemple) et des phases de projection durant lesquelles on cherche à déplacer l'échantillon suivant une autre dimension, toutes choses égales Par ailleurs. une fois trouvé dans l'espace à N dimensions une position acceptable pour l'échantillon--c'est à dire remplissant un ensemble de critères de validation eu égard à l'échantillon d'origine --, la propagation dans l'espace des dimensions spatiale peut reprendre, les (N-3) dimensions restantes étant de nouveau figées. L'ensemble des échantillons validés durant ces alternances de phases de propagation et de projection est cartographié au fur et à mesure que progresse l'algorithme de cartographie de l'hyper-surface dans l'espace à N dimensions. D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des figures annexées parmi lesquelles : - la figure 1 est une représentation d'un volume tridimensionnel dans lequel on a pointé différentes traces sismiques obtenues, - la figure 2 représente des volumes tridimensionnels V1, V2, et V3 correspondant à différentes valeurs d'un paramètre Xd, - la figure 3 illustre un exemple de guide de pointé de surface, - les figures 4 et 5 illustrent des exemples de propagation de pointé dans trois volumes tridimensionnels V1, V2, V3, - les figures 6 à 8 illustrent des exemples de guides de pointé de surface. On se place dans la suite du texte dans l'hypothèse où un ensemble de traces sismiques a été acquis à l'occasion d'une ou plusieurs campagnes d'enregistrement, par exemple d'enregistrement 3D. Ces enregistrements sont traités pour être exploités sous forme de données (par exemple, mais non nécessairement les valeurs d'amplitudes des ondes enregistrées) dépendant d'un nombre n de paramètres qui peut être bien supérieur à 3. Ces paramètres sont appelés Xi, X2, Xn dans la suite du texte. Ils peuvent correspondre aux dimensions géométriques X, Y, Z (position et profondeur ou temps) et/ou à toutes les autres grandeurs que l'on souhaite prendre en compte pour analyser les enregistrements dont on dispose : distance source/récepteur, angle d'incidence, temps entre acquisition, type d'onde, incrément de modèle de vitesse, etc. Dans une première étape on choisit trois de ces paramètres (xi, x, Xn par exemple) pour définir trois dimensions d'un sous espace tridimensionnel sur lequel il est possible de réaliser une cartographie 3D classique, si les autres paramètres sont fixés. Les données que l'on considère sont en général des collections de données sismiques 3D générées chacune pour une valeur donnée d'un ou plusieurs autres paramètres, que l'on fait varier pour définir des dimensions supplémentaires. Pour chaque volume 3D ainsi obtenu, il est possible de réaliser des cartographies, par exemple par pointé de surfaces. La figure 1 illustre un exemple de pointé de surface (surface Si) réalisé sur un tel volume Vi. Par ailleurs, on définit également une ou plusieurs relations dites de projection ou de guides de pointé de surface destinées à relier les points d'un de ces volumes à ceux d'un autre, en fonction des modifications que l'on apporte aux paramètres autres que ceux (xi, X2... xn) qui correspondent aux dimensions dans lesquelles ces volumes sont représentés, pour passer de l'un à l'autre de ces deux volumes. Ces guides de pointé peuvent être déterminés de différentes façons et notamment outre en fonction des valeurs des échantillons considérés, en fonction également de descriptions dynamiques ou cinématiques de l'évolution des volumes dans les directions correspondant aux paramètres modifiés. Les guides de pointé peuvent également être des lois que l'on se donne de façon empirique. Les figures 2 et 3 illustrent un exemple de ce que peut être un guide de pointé. V1, V2 et V3 correspondent à trois volumes 3D donnant des collections d'échantillons en fonction des trois mêmes paramètres Xi, X, 2 Xn pour trois valeurs différentes d'un autre paramètre Xd (le terme valeur étant ici à comprendre de façon large, Xd pouvant tout aussi bien correspondre à un paramètre quantifiable, qu'à un paramètre exprimant une condition, autre qu'une condition quantifiée (type d'ondes ondes P, ondes S, etc.)). S1 étant une surface pointée autour du point P1 dans le volume V1, un guide de pointé de surfaces est par exemple une courbe (Il sur la figure 3) qui dans le plan (P, Xd, X1) permet d'associer à un point P1 d'un pointé de surface dans le volume V1 des points P2 et P3 qui dans les volumes V2 et V3 décrivent théoriquement le même événement, pour les mêmes valeurs de X1, X2. Considérons une surface S associe au volume de dimension n. Elle est définie par : v (..., ) [l,.., N,] x... x [l,.., J S ... -.) i =ou Ud silepointn'existe pasou n'est pasdéfini Soit P (i1,..., in) un point de la surface S et d une direction autre que celles formant le prédéfini volume 3D. Un guide de pointé de surfaces associé à P dans une direction d est une application ni qui permet de déterminer une partie des points p appartenant à S se trouvant au voisinage de P dans le plan défini par (P, Xd, Xn). Si on sait définir l'application n pour un point P dans une direction d, on sait étendre la surface S à partir de P dans cette direction. On augmente ainsi la dimension du voisinage, dans le procédé de propagation, par le nombre de direction, autres que celles formant le prédéfini volume 3D, où une application n peut être définie. L'application FI décrit une courbe représentant la projection de l'image décrite par S dans le plan (P, Xd, Xn). La courbe est définie par des segments reliant les extrema corrélant de S d'une colonne à l'autre dans le plan (P, Xd, Xn). Notons R l'application reliant deux points voisins de S par un segment. On applique alors la relation inverse R-1 de façon à transformer le segment en éléments de surface élémentaire. L'application n permet ainsi de passer d'un volume 3D à un autre, pour pointer dans un volume puis dans un autre des événements en correspondance eu égard aux critères de propagation sélectionnés. Des tests de cohérence sont avantageusement prévus pour valider les points que l'on retient dans les volumes dans lesquels on se projette à l'aide du guide ri : calcul de semblance et/ou de cross-corrélation par exemple. On comprend que la technique qui vient d'être décrite permet successivement de mettre en oeuvre un pointé de surface dans un premier volume, de passer dans un deuxième volume par un guide de pointé, puis de reprendre un pointé de surface dans ce deuxième. En multipliant ce type de combinaisons, on augmente considérablement le nombre de cheminements possibles à partir d'un échantillon pointé, pour parvenir à la cartographie complète de l'événement sur l'ensemble de l'espace. La figure 4 illustre en ce sens, dans le cas d'une analyse prenant en compte 4 dimensions, des combinaisons de pointés de surfaces 2D ou 3D classiques (petites flèches sur les grilles) avec des guides de propagation permettant de passer d'un volume 3D à un autre. Dans cet exemple, on a illustré un cas où l'on prend en compte quatre paramètres, l'un de ces paramètres étant celui que l'on fait varier pour passer d'un volume à un autre. Il est bien entendu possible de prendre en compte un nombre de paramètres beaucoup plus important et de faire varier plusieurs paramètres. Il est également possible de changer les paramètres qui définissent les dimensions des volumes 3D dans lesquels les techniques classiques de pointé de surface sont mises en oeuvre. Comme l'illustre la figure 5, les techniques proposées sont par exemple utilisées pour contourner une discontinuité (référencée par D) empêchant de poursuivre un pointé de surface dans un volume. Ainsi, dans l'exemple illustré sur ladite figure 5, le pointe de surface est réalisé dans un volume de dimension Xi, X2, X4. Une discontinuité D empêche de réaliser ce pointé de surface sur l'ensemble de ce volume. On projette alors un ou plusieurs points mis en évidence par le pointé de surface dans un ou plusieurs autres volumes correspondant à une ou plusieurs autres valeurs pour un quatrième paramètre (en l'occurrence X3). Dans chacun de ces autres volumes, il est possible de mettre en oeuvre à nouveau des pointés de surfaces classiques, et ainsi de contourner la faille pour à partir d'un point mis en évidence par un de ces pointés de surface poursuivi dans un de ces autres volumes, revenir dans le premier volume par une projection inverse (point et ). On donne ci-après quelques exemples de guides de pointés de surfaces possibles. Exemple 1. On considère des sismiques tridimensionnelles corrigées au déport, Xi et X2 étant des dimensions correspondant au paramètre de position en surface, X3 correspondant au paramètre de déport, X4 correspondant au temps corrigé. Un guide de pointé de surface dans la direction X3 peut être : 1-13 (j) = Adj (i4) où Adj est la fonction d'ajustement au signal sur la trace courante d'un point de la surface et où J4 est une valeur échantillonnée selon la direction X4. Exemple 2. On considère le volume sismique représenté selon deux directions X1, X2 correspondant au paramètre de position au sol, issu d'une acquisition dont la distance source/récepteur est échantillonnée régulièrement avec un pas A3. X3 représente la distance de déport, X4 le temps, échantillonné avec un pas A4 On dispose pour le réflecteur à cartographier d'une carte de vitesse V (h, i2), vitesse de propagation de l'onde au point (ii, i2, i4) du sous-sol. Si on note Adj la fonction d'ajustement du signal sur la trace courante d'un point de la surface, l'utilisation de la loi de Dix permet de définir, pour un point P (i1, J2, is, i4) de la surface, un guide de propagation dans la direction X3 des offsets par la fonction (représentée à la figure 6) r13 (i 2- (i3) 2) X (A3) 2 n,). . L) (A. -'. Ici, le guide de pointé de surfaces dépend de la vitesse moyenne dans le sous-sol, qui est un paramètre externe au volume que l'on cherche à imager. Ce paramètre apporte une information de description dans la direction X3. Exemple 3. Dans le cadre de la mise à jour d'un modèle de vitesse par la méthode de balayage ou scan de vitesse, on fait varier la vitesse initiale d'un facteur de perturbation a. Après les différentes migrations profondeur, on obtient n volumes sismiques échantillonnés suivant la direction X3 correspondant au facteur de perturbation is=a. La perturbation de vitesse étant linéaire, on peut considérer en première approximation que pour une position X1, X2, la profondeur X4 de l'échantillon associé à un réflecteur évolue de manière linéaire d'un volume sismique à l'autre, dans la direction X3. Le guide de pointé de surface est donc défini comme une fonction linéaire (représentée à la figure 7) du type no) =-+jxA, p (j) = Adj'4 + j-i3 j X''3 où As est le pas d'échantillonnage selon la direction X3. Au fur et à mesure que l'on projette les pointés d'un volume sismique à l'autre et que l'on valide les pointés obtenus, on corrige le guide de pointé pour qu'il soit cohérents avec les pointés obtenus. Par exemple, chaque nouveau pointé dans un volume tridimensionnel constitue un nouveau point du guide n. Par interpolation ou par régression, on peut modifier le guide de pointé n pour qu'il corresponde aux points obtenus. Ceci peut permettre dans certains cas, de créer de façon incrémentale le guide de pointé ou d'en déduire certaines variables externes au volume V quantifiables analytiquement et apportant une information descriptive des relations entre les points de surface dans une ou plusieurs des directions. Comme on l'aura compris, la technique proposée permet de prendre en compte des méthodes de propagation 2D et 3D existantes et de les étendre à des volumes de dimensions 4 ou supérieure. On augmente ainsi le nombre de chemins d'extension possibles notamment pour les pointés de surface. On notera que les relations de projection ou guides de pointés peuvent dépendre aussi bien d'informations décrivant la dynamique ou la cinématique de l'environnement sismique que d'informations empiriques obtenues à partir des échantillons de l'acquisition ou d'informations externes. REVENDICATIONS 1. Procédé de traitement de données sismiques, caractérisé en ce qu'en vue de leur analyse en fonction d'au moins quatre dimensions, on met en oeuvre un pointé et/ou une propagation de pointé sur au moins un sousensemble de ces données qui correspond à une collection de ces données échantillonnées selon trois ou deux dimensions, la ou les grandeurs qui correspondent à la ou aux autres dimensions étant alors de valeurs fixes, puis on applique à au moins un des points mis en évidence par ce pointé et/ou cette propagation de pointés une loi de projection donnée afin de déterminer un point qui lui correspond dans un autre sous-ensemble correspondant à une collection de données échantillonnées selon les mêmes dimensions, pour une valeur différente d'au moins une grandeur qui correspond à au moins une autre dimension. 2. Procédé selon la revendication 1, caractérisé en ce qu'on met en oeuvre une propagation de pointé à partir d'au moins un point projeté dans un tel autre sous-ensemble. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce qu'on met en oeuvre une succession d'alternance de propagations de pointés et de projections. 4. Procédé selon l'une des revendications 2 ou 3, caractérisé en ce qu'on applique à au moins un des points déterminé par propagation de pointé dans un tel autre sous-ensemble une loi de projection inverse afin de déterminer un point qui lui correspond dans au moins un sous-ensemble à partir duquel on a mis en oeuvre une projection. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'on met en oeuvre sur un point projeté un test de cohérence pour le valider. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'on affiche une cartographie des points pointés et/ou projetés et/ou déterminés par propagation de pointés dans les différents sous-ensembles considérés. 7. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'un paramètre modifié d'un sous-ensemble à un autre est un déport source/récepteur et/ou un angle d'incidence et/ou un temps d'acquisition et/ou un type d'ondes sismiques et/ou une vitesse de migration et/ou un paramètre prenant en compte au moins l'un de ces paramètres. 8. Procédé selon l'une des revendications précédentes, caractérisé en ce que le nombre de dimensions prises en compte est supérieur à quatre. 9. Procédé de détermination cartographique sur des données sismiques, caractérisé en ce qu'on met en oeuvre le procédé selon l'une des revendications précédentes pour déterminer une cartographie selon au, moins quatre dimensions. 10. Procédé selon la revendication 9, caractérisé en ce qu'on détermine à partir d'un point projeté et/ou pointé et/ou déterminé par propagation de pointés un voisinage de points comportant en plus des points dans un voisinage prédéterminé dudit point dans le sous-ensemble de dimensions deux ou trois auquel ledit point appartient, un ou des points projetés à partir dudit point dans d'autres sous-ensembles et en ce qu'on détermine sur l'ensemble des points de ce voisinage celui de ces points qui optimise un critère donné. 11. Procédé selon l'une des revendications précédentes, caractérisé en ce que dans le cas de données sismiques tridimensionnelles corrigées au déport, les dimensions prises en compte comprennent deux dimensions correspondant aux paramètres de position en surface (X1, X2), une dimension (X3) correspondant au paramètre de déport, une dimension (X4) correspondant au temps corrigé et en ce qu'une loi de propagation selon la direction la direction correspondant au paramètre de déport (X3) est : ij3 (j) = Adj (i4) n-) =) où Adj est la fonction d'ajustement au signal sur la trace courante d'un point de la surface et où i4 est une valeur échantillonnée selon la dimension (X4) correspondant au temps corrigé. 12. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que les dimensions prises en compte comprennent deux dimensions correspondant aux paramètres de position en surface (Xi, X2), une dimension (X3) correspondant au paramètre de déport, une dimension (X4) correspondant au temps, et en ce que les données étant issues d'une acquisition avec une distance source/récepteur échantillonnée régulièrement avec un pas A3, et un temps échantillonné avec un pas A4, la loi de projection dans la direction du déport est r13 (i4) 2) 2 + (i 2- (i3) 2) x (, à3) 2 nO). . (A. - OU V (i1, i2) est une vitesse de propagation de l'onde dont on dispose au point (i1, J2, i4) du sous-sol, où Adj est la fonction d'ajustement du signal sur la trace courante d'un point de la surface. 13. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que pour mettre à jour un modèle de vitesse par un balayage de vitesse, on fait varier la vitesse initiale d'un facteur de perturbation (a) et en ce que après les différentes migrations profondeur, on obtient n volumes sismiques échantillonnés suivant une dimension (X3) correspondant au facteur de perturbation, en ce que d'autres dimensions prises en compte comprennent deux dimensions correspondant aux paramètres de position en surface (X1, X2) et une dimension (X4) correspondant à la profondeur de l'échantillon, et en ce qu'une loi de propagation selon la dimension correspondant au facteur de propagation est du type : r13 j-i, j j n ') =A ±- xA, où As est le pas d'échantillonnage selon cette dimension et où Adj est la fonction d'ajustement du signal sur la trace courante d'un point de la surface.