Brevet FR2841344A1 — publié le 20031226

L'invention est relative à un procédé de lissage d'une propriété de subsurface dans une structure géologique représentée par des mesures sismiques. L'invention est également relative à un dispositif de lissage d'une propriété de subsurface dans une structure géologique représentée par des mesures sismiques. L'invention est enfin relative à un produit-programme d'ordinateur permettant le fonctionnement d'un dispositif programmable de lissage d'une propriété de subsurface dans une structure géologique représentée par des mesures sismiques. Le document WO 01/63323 Al décrit un procédé de traitement de données sismiques comprenant les étapes consistant à : un volume de données sismiques couvrant un volume prédéterminé de terre ; déterminer pour chaque voxel du volume de données sismiques l'orientation locale des données sismiques ; déterminer pour chaque voxel s'il existe un bord dans le voisinage, et exécuter une opération de lissage sur chaque voxel dans le volume de données sismiques, dans lequel la direction de l'opération de lissage est l'orientation locale des données, et dans lequel l'opération de lissage ne dépasse pas le bord, de manière à obtenir un volume de données sismiques traitées, où la quantité associée à chaque voxel dans le volume de donnés traitées est le résultat obtenu par l'exécution de l'opération de lissage dans le voxel dans le volume de données sismiques. Le document WO 02/13139 Al décrit un procédé de traitement d'images sismiques comprenant les étapes : obtenir un ensemble de données d'une image initiale bidimensionnelle ou tridimensionnelle, dans lequel chaque élément de l'ensemble de données est l'intensité initiale d'image du point de l'image ; calculer pour chaque point les dérivées partielles de l'élément dans n directions pour obtenir un ensemble de valeurs dérivé des dérivées partielles ; pour chaque point une matrice carrée structurale symétrique à partir des valeurs des dérivées partielles; effectuer une itération point à point pondérée par une variable voisine de zéro lorsque le point est voisin d'un bord et voisine de 1 lorsqu'il est très éloigné d'un bord ; et répéter ces étapes un nombre de fois pour obtenir l'image traitée. Un premier but de l'invention est d'améliorer le lissage des propriétés de subsurface dans une structure géologique représentée par des mesures sismiques. Un deuxième but de l'invention est de permettre un lissage simple et rapide des propriétés de subsurface, tout en ne lissant pas les discontinuités. L'invention a pour objet un procédé de lissage d'une propriété de subsurface dans une structure géologique représentée par mesures sismiques, dans lequel on construit une fonction continue Sij,k(t) par interpolation ou approximation des traces sismiques discrètes d'une matrice sismique multidimensionnelle, ladite fonction étant désignée comme "trace sismique locale continue", comportant les étapes suivantes : a). utiliser comme décalage optimal hij,pq,k de deux traces sismiques locales continues voisines Sij,k(t) et Spq,k(t), la valeur de décalage rendant maximale leur fonction de corrélation ; b) retenir comme voisinage conditionnel d'une trace sismique locale continue "centrale" Sij,k(t) le sous-voisinage consistant en des traces voisines Spq,k(t) correspondant à des décalages optimaux hij,pq,k associés à des corrélations Rij,pq,k(hij,pq,k) supérieures à un seuil prédéterminé compris entre 0 et 1 ; c). choisir une propriété de subsurface à lisser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique locale continue "centrale" de référence ; d). décaler les propriétés de subsurface du voisinage conditionnel en translatant la variable courante de la valeur du décalage optimal hij,pq,k; e). prendre comme valeur lissée au point (i, j, k) une moyenne des propriétés de subsurface décalées à l'étape d). Selon d'autres caractéristiques de l'invention : - la moyenne des propriétés de surface de l'étape e) est une moyenne qui peut être pondérée, par exemple par la valeur de corrélation maximale correspondant au décalage optimal, - la moyenne des propriétés de subsurface de l'étape e) est choisie dans l'ensemble suivant : moyenne arithmétique, moyenne géométrique, moyenne harmonique, - on choisit comme propriété de subsurface à lisser l'amplitude réfléchie et captée par des géophones, - si nécessaire, l'ensemble du processus de lissage est répété un certain nombre de fois pour améliorer le lissage, - on visualise des coupes de matrice multidimensionnelle de propriété lissée sur un écran de visualisation. L'invention a également pour objet un dispositif pour la mise en #uvre du procédé selon l'invention, comportant des moyens pour utiliser comme décalage optimal de deux traces sismiques locales continues voisines la valeur de décalage rendant maximale leur fonction de corrélation, des moyens pour décaler les propriétés de subsurface du voisinage conditionnel en translatant la variable courante de la valeur du décalage optimal hij,pq,k, des moyens pour choisir une propriété de subsurface à lisser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique locale continue "centrale" de référence, des moyens pour décaler les propriétés de subsurface du voisinage conditionnel en translatant la variable courante de la valeur du décalage optimal hij,pq,k et des moyens pour prendre comme valeur lissée au point (i, j, k) une moyenne des propriétés de subsurface décalées à l'étape d). Selon d'autres caractéristiques de l'invention : -le dispositif comporte des moyens de mémorisation et des moyens de visualisation de paramètres sismiques déterminés à l'aide du procédé selon l'invention. L'invention a enfin pour objet un produit-programme d'ordinateur, comportant des éléments de code de programme pour exécuter les étapes d'un procédé selon l'invention, lorsque ledit programme est exécuté par un ordinateur. L'invention sera mieux comprise grâce à la description qui va suivre donnée à titre d'exemple non limitatif en référence aux dessins annexés dans lesquels : - La figure 1 représente schématiquement une matrice sismique tridimensionnelle et une trace sismique locale continue. - La figure 2 représente schématiquement, un exemple de voisinage local Nij,k consistant en un ensemble de trace sismique locale continues Spq,k(t) voisines de la trace centrale de référence Sij,k(t) elle même située a la verticale du géophone Gij et interpolant les amplitudes sismiques autour de t=k. - La figure 3 représente schématiquement, une section verticale de voisinage local conditionnel Nij,k(r) d'une trace sismique locale continue de référence Sij,k(t). - La figure 4 représente schématiquement, une section verticale d'une matrice sismique tridimensionnelle avec un décalage optimal hij,pq,k et un voisinage conditionnel Nij,k(r). - La figure 5 représente schématiquement, une section verticale d'une matrice sismique tridimensionnelle analogue à la figure 4 avec présence d'une faille ou discontinuité . - La figure 6 représente schématiquement, un organigramme fonctionnel d'un procédé selon l'invention. En référence aux figures 1 à 5, les éléments identiques ou fonctionnellement équivalents sont désignés ou repérés de manière identique. Sur la figure 1, une matrice sismique tridimensionnelle est obtenue par relevé de mesures enregistrées par des géophones Gij disposés sur un réseau x, y en des points de coordonnées i, j. Le relevé des mesures échantillonnées dans le temps est représenté suivant un axe t descendant représentatif de la profondeur ou d'une verticale descendante à partir de la surface de la terre ou de la mer. Les mesures sont caractérisées par leur amplitude, par exemple une amplitude relevée par le géophone Gij au temps ou à la profondeur d'échantillonnage tk. La mesure discrète effectuée par le géophone Gij au temps ou à la profondeur tk est appelée amplitude sismique Sij,k. L'ensemble des amplitudes sismiques correspondant à un géophone Gij de coordonnées i, j est une matrice unidimensionnelle (5;1, Sy2, Sijk .... e SijN) appelées trace sismique discrète, car cette matrice unidimensionnelle correspond à la trace selon le point de coordonnées horizontales i, j de la matrice sismique tridimensionnelle obtenue par mesures sismiques. L'axe vertical t orienté selon une verticale descendante désigne usuellement le temps, mais peut également être traité pour représenter une profondeur à partir de la surface. L'invention concerne aussi bien l'application à une troisième coordonnée t représentative du temps, qu'à une troisième coordonnée t représentative de la profondeur. A partir de la trace sismique discrète située à la verticale d'un géophone Gij on définit, par interpolation ou approximation des valeurs discrètes autour de t=tk=k, une fonction continue Sij,k(t) qui est désignée comme "trace sismique locale continue". Les méthodes d'approximation ou d'interpolation de valeurs discrètes pour engendrer une fonction continue sont nombreuses, et comprennent notamment les interpolations ou approximations polynomiales, ainsi que les interpolations ou approximations trigonométriques polynomiales. Toute autre variante d'interpolation ou d'approximation fournissant une fonction continue peut également être appliquée à la présente invention pour fournir une "trace sismique locale continue". Sur la figure 2, plusieurs traces sismiques locales continues définissent un exemple de voisinage d'une trace sismique locale continue "centrale" Sij,k(t) de référence. Le voisinage d'une trace sismique locale continue Sij,k(t) est défini comme l'ensemble des traces sismiques locales continues dont les indices spatiaux horizontaux p, q sont voisins des indices spatiaux horizontaux i, j de la trace sismique locale continue de référence. A titre d'exemple, les coordonnées spatiales horizontales p, q correspondant au géophone Gpq sont voisines des coordonnées spatiales horizontales i, j correspondant au géophone Gij si les valeurs absolues des différences i-p et j-q sont inférieures à des entiers donnés, par exemple à 2. Dans ce cas, comme représenté à la figure 2, la trace sismique locale continue Sij,k(t) est associée à huit traces sismiques locales continues voisines entourant la trace sismique locale continue "centrale" Sij,k(t). Dans le cas de traces sismiques locales continues produites par des mesures sismiques, le profil des horizons géologiques introduit des décalages verticaux entre des traces sismiques locales continues voisines. Dans le but de déterminer des relations entre deux traces sismiques continues locales voisines Sij,k(t) et Spq,k(t) centrées sur la même coordonnée verticale d'échantillonnage t=tk=k et correspondant à des coordonnées spatiales i, j et p, q différentes, on calcule la fonction de corrélation Rij,pq,k(h) des deux traces sismiques locales continues. La fonction de corrélation de deux traces sismiques locales continues voisines est obtenue par la formule suivante où le numérateur correspond à la fonction de covariance de Sij,k(t) et Spq,k(t) tu +e obtenue par l'expression suivante Cl,,k(h) = jS1Jk(t). S,k (t + h) . dt . tek-à Dans cette intégrale définissant Cij,pq,k(h) , le paramètre A définit une fenêtre verticale d'investigation autour de t=tk=k. Par exemple, si les traces sismiques locales continues Sij,k(t) et Spq,k(t) sont des polynômes trigonométriques de la forme suivante interpolant les donnés sismiques Cij.pq.k(h) = A-.cos(sû)h)+B.sm(sû)h) s=1 Un résultat mathématique classique et connu est que la fonction de corrélation lorsque cette fonction de corrélation se rapproche de 1. L'étude des fonctions de corrélation des traces sismiques locales continues Sij,k(t) et Spq,k(t) permet de définir comme décalage optimal la valeur hij,pq,k de h correspondant au maximum de la fonction de corrélation Rij,pq,k(h) , c'est-à-dire au maximum de la fonction de corrélation le plus proche de 1. Contrairement à l'art existant basé sur une formulation discrète de la fonction de corrélation Rij,pq,k(h), l'utilisation d'une formulation continue des traces sismiques et donc de la fonction de corrélation permet d'obtenir un décalage optimum qui n'est pas contraint à être un multiple entier du pas d'échantillonnage suivant l'axe vertical correspondant à la variable t. La définition du décalage optimal hij,pq,k de deux traces sismiques locales continues voisines permet d'obtenir une première approximation des horizons passant par cette trace sismique locale continue de référence Sij,k(t), comme représenté aux figures 3 à 5. Contrairement à l'art existant, le fait d'utiliser des décalages optimaux qui ne sont pas contraint à être un multiple entier du pas d'échantillonnage suivant l'axe vertical permet d'éviter des erreurs connues sous le nom d'aliasing . Sur ces figures 3 à 5, on définit un voisinage conditionnel Nij,k(r) de la trace sismique continue locale de référence Sij,k(t) qui est un sous-voisinage du voisinage initial des traces sismiques locales continues retenues pour le calcul de corrélation et de définition des décalages optimaux. Le voisinage conditionnel Nij,k(r) est choisi de manière que, pour toute trace sismique continue locale Spq,k(t) appartenant à Nij,k(r), le décalage optimal hij,pq,k correspondant est tel que la corrélation Rij,pq,k(hij,pq,k) entre Sij,k(t) et Spq,k(t) est supérieure à un seuil prédéterminé r compris entre 0 et 1. Comme on le voit sur la figure 5, la fixation du seuil r permet également d'englober dans un voisinage conditionnel une discontinuité ou une faille géologique, ce qui constitue un avantage important par rapport à l'art antérieur. Les dispositions précitées fournissent ainsi une modélisation continue permettant la mise en #uvre de l'invention. En référence à la figure 6, un procédé de lissage ou propriété de subsurface dans une structure géologique représenté par des mesures sismiques comporte une première étape 100 de modélisation numérique, comme décrit en référence aux figures 1 à 5, pour définir des traces sismiques locales continues d'une matrice sismique multidimensionnelle, calculer les décalages optimaux de traces sismiques locales continues voisines et définir les voisinages conditionnels des traces sismiques locales continues "centrales" de référence. Après cette première étape 100 de modélisation numérique, une étape 101 permet de choisir la propriété à lisser sur un voisinage conditionnel. Cette propriété à lisser sur un voisinage conditionnel peut consister en toute fonction multidimensionnelle de la subsurface à analyser. En particulier, la propriété à lisser peut simplement consister en l'amplitude du signal sismique réfléchi par les horizons de la subsurface en direction des géophones de mesure. A l'étape 102, on remplace les propriétés du voisinage conditionnel par les propriétés décalées de la valeur du décalage optimal hij,pq,k, de manière à constituer un ensemble de propriétés relatives au même horizon. Ainsi, dans le cas d'une trace locale sismique continue centrale de référence Sij,k(t), on remplace toutes les valeurs Spq,k(t) du voisinage conditionnel par les valeurs Spq,k(t+hij,pq,k) . On effectue ensuite à l'étape 103 une moyenne des valeurs décalées Spq,k(t+hij,pq,k) établies à l'étape 102 et on choisit comme valeur de la propriété au point de coordonnées spatiales i,j et de coordonnée temporelle ou de profondeur t=tk=k , la moyenne ainsi calculée. L'invention s'applique à plusieurs variantes d'établissement de moyennes : on peut par exemple utiliser une moyenne pondérée ou non pondérée ; ou, alternativement, on peut également utiliser une moyenne arithmétique, une moyenne géométrique, une moyenne harmonique , pondérée ou non pondérée. Dans le cas d'une moyenne pondérée, ont peut avantageusement utiliser comme facteur de pondération la valeur maximale de la fonction de corrélation Rij,pq,k (hij,pq,k) , en reflétant ainsi le degré de confiance accordable à chaque décalage optimal correspondant à un maximum d'une fonction de corrélation. Le procédé de lissage ou de filtrage d'une propriété selon la présente invention permet ainsi d'améliorer le contraste et la précision des propriétés de subsurface dans une matrice sismique multidimensionnelle. Si nécessaire, afin d'améliorer le lissage, l'ensemble du processus décrit par la presente invention peut être répété sur le lissage obtenu lors d'une iteration précédente. L'étape 104 constitue une étape itérative dans laquelle l'opérateur peut définir le nombre d'itérations adequates à la résolution de son problème. L'invention est de préférence mise en #uvre sur un dispositif comportant des moyens pour accomplir successivement les étapes du procédé décrit en référence à la figure 6. En particulier, ce dispositif comporte des moyens de visualisation des paramètres sismiques ou propriétés obtenues à l'étape 105 et des moyens de mémorisation de programme et de calculs intermédiaires. De préférence, ce dispositif est un dispositif programmable commandé par un produit-programme d'ordinateur, réalisé pour mettre en #uvre un procédé selon l'invention. L'invention décrite en référence à plusieurs objets particuliers n'y est nullement limitée mais couvre au contraire toute modification de forme et toute variante de réalisation dans le cadre et l'esprit de l'invention, l'essentiel étant d'utiliser des valeurs de propriétés recalées par rapport aux décalages optimaux déterminés comme valeurs correspondant à la maximisation des fonctions de corrélation de traces sismiques locales continues ou même discrètes. REVENDICATIONS 1. Procédé de lissage d'une propriété de subsurface dans une structure géologique représentée par mesures sismiques, dans lequel on construit une fonction continue Sij,k(t) par interpolation ou approximation des traces sismiques discrètes d'une matrice sismique multidimensionnelle, ladite fonction Sij,k(t) étant désignée comme "trace sismique locale continue", comportant les étapes suivantes : a). utiliser comme décalage optimal de deux traces sismiques locales continues voisines Sij,k(t) et Spq,k(t), la valeur de décalage rendant maximale leur fonction de corrélation, ce décalage optimal n'étant pas contraint à être un multiple entier du pas d'échantillonnage vertical ; b). retenir comme voisinage conditionnel d'une trace sismique locale continue "centrale" Sij,k(t) le sous-voisinage consistant en des traces voisines Spq,k(t) correspondant à des décalages optimaux associés à des corrélations Rij,pq,k(h) supérieures à un seuil prédéterminé compris entre 0 et 1 ; c). choisir une propriété de subsurface à lisser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique locale continue "centrale" de référence ; d). décaler les propriétés de subsurface du voisinage conditionnel en translatant la variable courante de la valeur du décalage optimal (hij,pq,k) ; e). prendre comme valeur lissée au point (i, j, k) une moyenne des propriétés de subsurface décalées à l'étape d). 2. Procédé selon la revendication 1, dans lequel la moyenne des propriétés de surface de l'étape e) est une moyenne pondérée, par exemple par la valeur de corrélation maximale correspondant au décalage optimal. 3. Procédé selon la revendication 1 ou 2, dans lequel la moyenne des propriétés de subsurface de l'étape e) est choisie dans l'ensemble suivant : moyenne arithmétique, moyenne géométrique, moyenne harmonique pondérée ou non. 4. Procédé selon la revendication 1, dans lequel on choisit comme propriété de subsurface à lisser, l'amplitude réfléchie et captée par des géophones. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel on itère le procéssus de lissage en l'appliquant sur le résultat d'un lissage antérieur. 6. Procédé selon la revendication 1 ou la revendication 5, dans lequel on visualise une matrice multidimensionnelle de propriété lissée sur un écran de visualisation. 7. Dispositif pour la mise en #uvre du procédé selon l'une quelconque des revendications 1 à 4, comportant des moyens pour utiliser comme décalage optimal de deux traces sismiques locales continues voisines la valeur hij,pq,k de décalage rendant maximale leur fonction de corrélation, des moyens pour décaler les propriétés de subsurface du voisinage conditionnel en translatant la variable courante de la valeur du décalage optimal hij,pq,k, des moyens pour choisir une propriété de subsurface à lisser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique locale continue "centrale" de référence, des moyens pour décaler les propriétés de subsurface du voisinage conditionnel en translatant la variable courante de la valeur du décalage optimal hij,pq,k et des moyens pour prendre comme valeur lissée au point (i, j, k) une moyenne des propriétés de subsurface décalées à l'étape d). 8. Dispositif selon la revendication 7, comportant des moyens de mémorisation et des moyens de visualisation de paramètres sismiques déterminés à l'aide du procédé selon l'une quelconque des revendications 1 à 6. 9. Produit-programme d'ordinateur, comportant des éléments de code de programme pour exécuter les étapes du procédé selon l'une quelconque des revendications 1 à 5, lorsque ledit programme est exécuté par un ordinateur. 10. Produit-programme d'ordinateur, comportant des éléments de code de programme pour exécuter les étapes du procédé selon la revendication 6, lorsque ledit programme est exécuté par un ordinateur.