FR2799552A1 - Traitement de mesures de formes d'ondes acoustiques provenant d'outils de diagraphie de sondage en matrice - Google Patents
Traitement de mesures de formes d'ondes acoustiques provenant d'outils de diagraphie de sondage en matrice Download PDFInfo
- Publication number
- FR2799552A1 FR2799552A1 FR0012733A FR0012733A FR2799552A1 FR 2799552 A1 FR2799552 A1 FR 2799552A1 FR 0012733 A FR0012733 A FR 0012733A FR 0012733 A FR0012733 A FR 0012733A FR 2799552 A1 FR2799552 A1 FR 2799552A1
- Authority
- FR
- France
- Prior art keywords
- matrix
- receivers
- nsb
- rsb
- rref
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Des procédés pour traiter des enregistrements de diagraphie acoustique obtenus à partir d'un outil de diagraphie comportant une matrice de R récepteurs, comportent la sélection d'une dimension de matrice secondaire de Rsb récepteurs, avec Rsb < R; le calcul du nombre de matrices secondaires Nsb possibles à partir de la matrice de R récepteurs; l'identification d'un récepteur Rref dans la matrice pour servir de récepteur de référence pour chaque matrice secondaire; le traitement de données dans le plan lenteur-temps pour chaque matrice secondaire par rapport au récepteur de référence Rref ; et l'empilement des données traitées dans le plan lenteur-temps pour toutes les matrices secondaires.
Description
TRAITEMENT <B>DE MESURES DE</B> FORMES <B>D'ONDES ACOUSTIQUES</B> PROVENANT <B>D'OUTILS DE</B> DIAGRAPHIE <B>DE SONDAGE EN</B> MATRICE <U>Domaine technique</U> La présente invention concerne des procédés pour traiter des mesures de formes d'ondes acoustiques, particulier des mesures de formes d'ondes acoustiques effectuées dans le but de caractériser les propriétes de formations souterraines. L'invention concerne en particulier des procédés de traitement de mesures acoustiques réalisées avec des outils en matrice.
Etat <U>de la technique antérieure</U> Les outils de diagraphie acoustique, destinés être utilisés pour caractériser des formations souterraines entourant un sondage en mesurant effets de la formation sur un signal acoustique se propageant à travers celui-ci, sont bien connus. Ces outils comportent couramment une matrice de transducteurs, habituellement des récepteurs, destinés à améliorer la sensibilité et la précision des mesures effectuées. Un exemple d'un outil en matrice destiné à être utilisé dans un sondage pour des mesures acoustiques est l'outil DSI de Schlumberger, qui est représenté schématiquement sur la figure 1. L'outil DSI comprend une partie émettrice 10 possédant une paire de sources de dipôles (supérieure et inférieure) 12, agencées orthogonalement dans le plan radial et une source monopôle 14. Un joint d'isolation acoustique 16 raccorde la partie émettrice 10 à une partie réceptrice 18 qui contient une matrice de huit stations réceptrices séparées, contenant chacune deux paires d'hydrophones, l'une orientée en ligne avec l'une des sources de dipôles, l'autre avec la source orthogonale. Une cartouche de circuits électroniques 20 est connectée au sommet de la partie réceptrice 18 et permet une communication entre l'outil et une unité de commande 22 située à la surface, par 'intermédiaire d'un câble électrique 24. Avec un tel outil, il est possible d'effectuer des mesures à la fois de monopoles et de dipôles. L'outil DSI possède plusieurs modes de fonctionnement d'acquisition de données tous pouvant être combinés pour acquérir (numériser) des formes d'ondes. Les modes sont : mode de dipôles supérieur et inférieur (UDP, =)-formes d'ondes enregistrées à partir de paires de récepteurs alignées avec la source de dipôles respective utilisée pour générer le signal ; mode à dipôles croisés-formes d'ondes enregistrées à partir de chaque paire de récepteurs pour amorcer la source de dipôles en ligne et croisés ; modes Stoneley- formes d'ondes de monopoles provenant l'amorçage à basse fréquence de la source monopôle ; modes P et S (P & S)-formes d'ondes de monopoles provenant de l'amorçage à haute fréquence de l'émetteur monopôle ; et premier mode de mouvement-seuil de monopôle croisant des données provenant de l'amorçage à haute fréquence de la source monopôle.
Une façon de déterminer les lenteurs de cisaillement et de Stoneley en compression à partir de ces mesures consiste à utiliser un traitement de cohérence lenteur-temps (STC). Le traitement STC est une technique d'analyse de formes d'ondes complètes qui a pour but de trouver toutes les ondes se propageant dans la forme d'onde composite. Le traitement adopte un algorithme de similitude pour détecter des arrivées qui sont cohérentes sur la matrice de récepteurs et estime leur lenteur. L'algorithme de base fait avancer une fenêtre temporelle de longueur fixe sur longueurs d'ondes, par petits échelons se recouvrant, sur une plage de temps d'arrivée potentiels. Pour chaque position dans le temps, la position de fenêtre est déplacée linéairement dans le temps, sur les formes d'ondes de la matrice de récepteurs, en commençant par un déplacement correspondant à l'onde la plus rapide attendue et en passant à l'onde la plus lente attendue. Pour chaque déplacement, une fonction de cohérence est calculée pour mesurer la similitude des ondes à l'intérieur de la fenêtre. Lorsque la fenêtre temporelle et le déplacement correspondent au temps d'arrivée et à la lenteur d'un composant particulier, les formes d'ondes à l'intérieur de la fenêtre sont presque identiques, fournissant une valeur de cohérence élevée. De cette manière, l'ensemble de formes d'ondes provenant de la matrice est examiné sur une certaine plage de temps d'arrivée possibles et de lenteurs pour les composantes de la forme d'ondes. Le traitement STC produit des tirs de contours de cohérence (similitude) dans le plan lenteur/temps d'arrivée. Les régions de cohérence forte correspondent à des arrivées particulières dans les formes d'ondes. La lenteur et le temps d'arrivée à chaque crête de cohérence sont comparés aux caractéristiques de propagation attendues des arrivées recherchées et celles qui correspondent le mieux à ces caractéristiques sont conservées. La classification des arrivées de cette manière produit un enregistrement continu de lenteur en fonction de la profondeur. Pour des ondes dispersives, le traitement STC est modifié pour prendre en compte l'effet de la fréquence. Puisque la sortie du traitement STC est un tracé de cohérence, la cohérence de chaque arrivée peut etre utilisée comme indicateur de qualité, des valeurs supérieures entraînant un caractère répétitif de la mesure plus important. Lors du traitement de formes 'ondes de dipôles, l'une des crêtes de cohérence correspond au mode en flexion, mais avec une lenteur qui est toujours plus grande (plus lente) que la lenteur de cisaillement réelle. Une correction calculée à l'avance est utilisée pour supprimer ce biais. D'autres détails concernant le traitement STC peuvent trouver dans Kimball C. V.<I>and</I> Marzetta T.L., "Sernblance processing <I>of</I> borehole acoustic array data" Geophysics, <I>Vol. 49, No. 3 (mars 1984) ; pages 274-28I.</I>
Pour compenser les variations des mesures dues forage plutôt que celles dues à la formation, une série mesures sont effectuées sur un intervalle dans lequel les propriétés de la formation sont attendues peu varier, s'il y a lieu. Dans sa forme la plus simple, l'intervalle correspond à l'importance de matrice de récepteurs et aux formes d'ondes au niveau de chaque station réceptrice, mesurées pour un amorçage donné d'une source ("matrice de récepteurs" ou "mode récepteur" ou "Rec."). Dans le traitement STC simple toutes les stations réceptrices sont considérées. Dans traitement STC à tirs multiples (MSTC), des matrices secondaires de stations réceptrices à l'intérieur de la matrice de récepteurs sont considérées, par exemple, matrice secondaire de cinq stations réceptrices dans une matrice de récepteurs de huit stations receptrices (d'autres nombres de stations réceptrices dans la matrice secondaire peuvent être utilisés en fonction des exigences). Dans ce cas, le même intervalle de formation correspondant à l'étendue d'une matrice secondaire de cinq stations réceptrices peut être mesuré plusieurs fois tandis que l'outil sonde le forage, les cinq stations constituant la matrice secondaire étant choisies au niveau de chaque amorçage de source pour mesurer le même intervalle de formation. Une autre approche, appelée "mode émetteur" ou "matrice pseudo-émetteur" ("Tra."), prend les formes d'ondes provenant d'amorçage de sources séquentielles ' mesure l'émetteur passe le long de l'intervalle à mesurer.
Pour compenser le mouvement de l'outil entre les mesures, une station réceptrice effectivement immobile ou une matrice secondaire doit être utilisée. Ceci peut être obtenu en modifiant la station réceptrice considérée de façon que sa position dans le forage soit effectivement immobile lorsque l'émetteur se déplace à travers l'intervalle. Une compensation de forage ("BHC") peut être réalisée pour les résultats modes P et S, en traitant les formes d'ondes de la matrice de récepteurs et de la matrice de pseudo-émetteurs et en moyennant les résultats. Un exemple schématique d'un outil de matrice de huit récepteurs effectuant des mesures sur un certain intervalle de profondeur d'une longueur de cinq récepteurs pour des amorçages successifs d'un émetteur, est représenté sur la figure 2. En utilisant la matrice secondaire de cinq récepteurs, quatre mesures sont effectuées même intervalle de profondeur, avec l'émetteur dans les positions i, i+1, i+2 et i+3. Le traitement de chaque matrice secondaire produit une série de tracés temporels de ralentissement comme représenté la figure 3(a). S'il est désiré combiner les informations obtenues à partir de cette série d'amorçages d'émetteurs, il est nécessaire de prendre en compte le décalage temporel entre chaque matrice secondaire dû au mouvement de la source. Ceci a été obtenu précédemment en projetant chaque tracé temporel de lenteur sur l'axe de lenteur (figure 3(b)), suivi par l'empilement des quatre tracés en un (figure 3(c)). D'autres détails concernant cette approche peuvent se trouver dans HSU K. and<I>Chang S. K.,</I> "Multiple-shot processing <I>of</I> array sonic waveforms", Geophysics, <I>Vol. 52, No. 10 (octobre</I> <I>1987) ; pages 1376-1390.</I> Il faut comprendre que l'effondrement des tracés ST séparés sur l'axe de lenteur fournit des informations utiles dérivées de l'axe des temps de chaque tracé. US 4 543 648 décrit également des procédés de traitement d'un tir à un autre de formes d'ondes acoustiques.
<U>Exposé de l'invention</U> Cette invention fournit des procédés pour traiter enregistrements de diagraphie acoustique obtenus à partir d'un outil de diagraphie comportant une matrice R récepteurs, comprenant la sélection d'une dimension de matrice secondaire Rsb récepteurs, avec Rsb < R ; le calcul du nombre de matrices secondaires possibles à partir de la matrice de R récepteurs ; l'identification d'un récepteur Rref dans matrice pour servir de récepteur de référence pour chaque matrice secondaire ; le traitement de données dans le plan lenteur- temps pour chaque matrice secondaire par rapport récepteur de référence Rref ; et l'empilement des données traitées dans le plan lenteur-temps pour toutes les matrices secondaires. L'empilement peut être réalisé en utilisant une sommation algébrique ou géométrique, une moyenne algébrique étant préférée. Le traitement dans le plan lenteur temps peut comprendre l'empilement de cohérences lenteur-temps pour chaque matrice secondaire ou traitement de données pour obtenir des piles comprenant l'énergie de cohérence et l'énergie totale dans plan lenteur temps pour les matrices secondaires, puis l'obtention du rapport des deux pour déterminer données de cohérence lenteur-temps empilées.
Il est préféré utiliser le nombre maximum possible de matrices secondaires pour une matrice de récepteurs donnée pour maximiser la quantité d'informations obtenues à partir de la formation. Pour une matrice de R recepteurs (fournie parfois ci-dessous par Nrcvr) , le nombre maximum de matrices secondaires Nsb de Rsb récepteurs en longueur, est fourni par Nsb = R _ Rsb + 1 (1) Pour un outil à matrice de huit récepteurs, une dimension de matrice secondaire de cinq récepteurs permet d'enregistrer quatre ensembles de formes d'ondes pour une profondeur donnée.
L'identification des traces (formes d'ondes) à enregistrer peut être obtenue en calculant une matrice géométrique G de dimension Nsb sur Rsb, qui fournit les récepteurs à considérer pour chaque tir à une profondeur donnée. L'identification des traces à partir de matrice complète pour chaque tir peut être obtenue en calculant une matrice de tirs multiples M de dimension Nsb. Nrcvr + Nsb - 1 (2) Les traces qui ne sont pas considérées à partir d'un tir donné reçoivent des valeurs de 0 dans une telle matrice. Le récepteur de référence rref pour un ensemble de matrices secondaires est identifié d'apres la relation (1), où Nsb identifie également la position du récepteur de référence rref dans la matrice.
<U>Brève description des dessins</U> La figure 1 représente un outil de diagraphie acoustique de l'art antérieur ; la figure 2 représente un dessin schématique d'une mesure de tirs multiples d'une région de forage ; les figures 3(a) à 3(c) représentent un tracé illustrant un aménagement de traitement de l'art antérieur ; et la figure 4 représente un organigramme avec les principales étapes de traitement d'un procédé selon l'invention.
<U>Exposé d'un mode de réalisation</U> L'invention va être décrite au moyen d'un exemple, en référence à une matrice de huit récepteurs (outil du type DSI) utilisant des matrices secondaires de cinq récepteurs en longueur. Il faut comprendre que le nombre de récepteurs dans la matrice peut être différent en fonction du type d'outil utilisé. Dans cette description, le terme "récepteur est utilisé pour se référer à des détecteurs uniques (par exemple, hydrophones) ou à une station avec un certain nombre de détecteurs disposés autour de la circonférence de l'outil à un certain nombre de niveaux (par exemple, quatre hydrophones au niveau de chaque station réceptrice). Le nombre de récepteurs dans chaque matrice secondaire est également choisi par l'utilisateur. La matrice secondaire la plus petite comporte trois récepteurs car le traitement d'un plan lenteur-temps nécessite l'utilisation de paires de mesures de temps d'arrivée, typiquement, à partir de paires adjacentes de récepteurs. La dimension de la matrice secondaire détermine la résolution de la détermination de la lenteur effectuée finalement, plus la matrice est petite, plus le lit pouvant etre mesuré est mince. Pour un outil avec un espacement de récepteurs de 6", ceci signifie qu'un outil à huit récepteurs peut mesurer des intervalles de profondeur entre l' et 3,5'.
Des descriptions générales des techniques d'acquisition et de traitement pouvant être appliquées à cette invention peuvent se trouver dans la demande de brevet UK No. 9907620.0 et dans US 4 8 236, US 5 278 805 et US 5 594 706.
Le traitement de la présente invention tente d'améliorer la résolution des calculs matrice complète STC (DSTC) dispersive en traitant les résultats provenant de matrices secondaires et en suréchantillonnant l'intervalle intéressant. Bien que ceci réduise la quantité de données disponibles provenant de chaque matrice secondaire, fait que l'outil ait des occasions multiples d'échantillonner la région avec des matrices secondaires différentes de la même taille signifie que ceci peut être compensé avec la redondance des données obtenues en empilant les plans ST.
La méthodologie de traitement générale de la présente invention est la suivante Dans le cas de tirs multiples à profondeur A1 est la forme d'onde appliquée a l'entrée du traitement STC. i est le nombre de matrices secondaires (i = , N) .
Fi est le plan ST calculé pour la matrice secondaire i.
Une trame est représentée comme une fonction deux paramètres Fi = Fi (t, S) , où t est le temps et S la lenteur. Pour N matrices secondaires, N = 1, 5 par exemple, cinq trames F = {F1, F2, F3, F4, F5} sont obtenues, où F est la famille de trames à profondeur Z.
Le de la technique à tirs multiples consiste a calculer
Mais puisqu'il existe un décalage de temps entre deux trames consécutives, il est nécessaire de prendre en compte ce décalage de temps avant d'empiler les plans Entre les trames Fi et Fi-,, l'outil se déplace d'une distance connue d = Oz (voir la figure 2). Cette distance d est connue. Le mouvement de l'outil entre des tirs consécutifs signifie qu'il existe un décalage de temps entre les différentes matrices secondaires. Pour trame Fi(t,S) et un point A(ti, Si) situé dans ce plan ST, pour décaler ce point dans le temps, il est nécessaire de le décaler de la valeur tti = ti.d S (ti) (4) Pour l'ensemble de la trame, après décalage de temps, ceci donne F'i rS) = Fi(t-(i-1)d.S(t),S) (5) Cette équation, appelée équation de décalage de temps, fournit la relation entre une trame de données avant après correction de temps. La connaissance de la façon de corriger ce décalage de temps permet de définir une équation générale de la technique de traitement à tirs multiples.
Mais puisqu'il existe un décalage de temps entre deux trames consécutives, il est nécessaire de prendre en compte ce décalage de temps avant d'empiler les plans Entre les trames Fi et Fi-,, l'outil se déplace d'une distance connue d = Oz (voir la figure 2). Cette distance d est connue. Le mouvement de l'outil entre des tirs consécutifs signifie qu'il existe un décalage de temps entre les différentes matrices secondaires. Pour trame Fi(t,S) et un point A(ti, Si) situé dans ce plan ST, pour décaler ce point dans le temps, il est nécessaire de le décaler de la valeur tti = ti.d S (ti) (4) Pour l'ensemble de la trame, après décalage de temps, ceci donne F'i rS) = Fi(t-(i-1)d.S(t),S) (5) Cette équation, appelée équation de décalage de temps, fournit la relation entre une trame de données avant après correction de temps. La connaissance de la façon de corriger ce décalage de temps permet de définir une équation générale de la technique de traitement à tirs multiples.
commençant par l'équation (5) ci-dessus, il est possible de la réécrire sous la forme d'une equation générale F (t,s) = Fi(t,s) *8(t-(i-1)d.s(t) ,s) (6) * est l'opérateur de convolution sur la variable t.
Pour rééchantillonner la trame pour l'empiler avec les autres, la fonction shah est utilisée, définie par
Pour un signal s(t), il est possible de réécrire ce signal sous la forme (t) = s (t) shah(-) (8) Ainsi, le signal s (t) et discrétisé avec taux d'échantillonnage La réécriture de l'équation sur une trame décalée dans le temps et après rééchantillonnage donne
réagençant l'équation (9), on obtient
Ceci définit la transformation à appliquer sur chaque trame avant empilement.
Pour un signal s(t), il est possible de réécrire ce signal sous la forme (t) = s (t) shah(-) (8) Ainsi, le signal s (t) et discrétisé avec taux d'échantillonnage La réécriture de l'équation sur une trame décalée dans le temps et après rééchantillonnage donne
réagençant l'équation (9), on obtient
Ceci définit la transformation à appliquer sur chaque trame avant empilement.
Connaissant la transformation qu'il est nécessaire d'appliquer sur une trame i, il est possible d'écrire l'équation générale de traitement de tirs multiples. Pour un niveau Z donné, celle-ci est
nT) (11) La partie 8(t-(i-1)d.S(t)) est 'opération de décalage dans le temps. Puisque celle- ' est la partie la plus coûteuse en temps CPU, elle est optimisée en utilisant l'opérateur de décalage de temps optimum déterminé de la manière suivante En général, l'opérateur de décalage dans le temps peut être décomposé en deux parties Za = Zn + Z# où n= [a] (: partie entière) et a=a-n(@(x j < _1: partie réelle ou fractionnaire).
nT) (11) La partie 8(t-(i-1)d.S(t)) est 'opération de décalage dans le temps. Puisque celle- ' est la partie la plus coûteuse en temps CPU, elle est optimisée en utilisant l'opérateur de décalage de temps optimum déterminé de la manière suivante En général, l'opérateur de décalage dans le temps peut être décomposé en deux parties Za = Zn + Z# où n= [a] (: partie entière) et a=a-n(@(x j < _1: partie réelle ou fractionnaire).
Zn peut être obtenu par simple rééchantillonnage. Il existe un grand nombre de manières pour décaler les données de a échantillons. Une opération du déphasage dans le domaine fréquentiel est une possibilité. Toutefois, lorsque les données sont dans le domaine temporel, il peut être coûteux de transformer les données du domaine temporel à fréquentiel et d'y revenir à la fin. Le processus propose ici une convolution standard utilisant un opérateur tronqué. Une quelconque série dans le temps peut être exprimée par
En conséquence, un décalage appliqué à la série dans le temps est équivalent à l'application d'un décalage avec un changement de signe sur fonction $.
En conséquence, un décalage appliqué à la série dans le temps est équivalent à l'application d'un décalage avec un changement de signe sur fonction $.
Pour une série dans le temps de quatre échantillons, par exemple, les opérateurs pour la transformation de Fourier rapide, le décalage dans le temps et la transformation inverse, peuvent être exprimés sous la forme matricielle i) Transformation de Fourier rapide
Opérateurs de décalage dans le temps dans le domaine fréquentiel
ou A = na/2 Opérateur de transformation de Fourier inverse
'application de ces opérateurs la fonction delta donne (n1 - 7, `
Opérateurs de décalage dans le temps dans le domaine fréquentiel
ou A = na/2 Opérateur de transformation de Fourier inverse
'application de ces opérateurs la fonction delta donne (n1 - 7, `
_ <SEP> (sin'8=sin9)Z-'+(cos'-9+cose)+
<tb> 2 <SEP> (sin' <SEP> 6 <SEP> - <SEP> sin <SEP> 8)Z <SEP> + <SEP> (cos' <SEP> <B>0</B> <SEP> - <SEP> cos <SEP> 6)Z' L'équation (11) est l'équation générale de tirs multiples qui ne fait aucune hypothèse sur l'outil utilisé mais qui suppose un milieu homogène. Pour éviter le problème d'un temps calcul important pour le rééchantillonnage dans le temps et le décalage dans le temps pendant l'exécution procédé général, la présente invention tire avantage du fait que l'espacement des récepteurs et mouvement de l'outil d'un incrément de profondeur sont les mêmes pour calculer le décalage dans temps pendant le traitement d'un plan ST pour chaque matrice secondaire.
<tb> 2 <SEP> (sin' <SEP> 6 <SEP> - <SEP> sin <SEP> 8)Z <SEP> + <SEP> (cos' <SEP> <B>0</B> <SEP> - <SEP> cos <SEP> 6)Z' L'équation (11) est l'équation générale de tirs multiples qui ne fait aucune hypothèse sur l'outil utilisé mais qui suppose un milieu homogène. Pour éviter le problème d'un temps calcul important pour le rééchantillonnage dans le temps et le décalage dans le temps pendant l'exécution procédé général, la présente invention tire avantage du fait que l'espacement des récepteurs et mouvement de l'outil d'un incrément de profondeur sont les mêmes pour calculer le décalage dans temps pendant le traitement d'un plan ST pour chaque matrice secondaire.
Dans le cas présent, le récepteur de référence de chaque matrice secondaire est remplacé avant de calculer la similitude de cette matrice secondaire pour prendre en compte le mouvement de l'outil et le décalage dans le temps. L'approche préférée consiste à utiliser comme référence le premier récepteur de la première matrice secondaire, deuxième récepteur de la deuxième matrice secondaire et ainsi de suite. À chaque fois, la similitude pour chacune de ces matrices secondaires est calculée. Ceci prend en compte le décalage de temps entre différentes matrices secondaires, de sorte qu'il est possible d'empiler les similitudes immédiatement après calcul, le décalage dans le temps entre les différentes matrices secondaires étant pris en compte. I1 faut comprendre que l'effet de remplacement du récepteur de référence de cette manière a pour effet d'utiliser la sortie d'un seul récepteur comme référence pour chaque tir (car ce récepteur change de position dans la matrice secondaire à mesure que l'outil progresse dans le forage). Dans l'exemple de la figure 2, le récepteur 4 est le récepteur de référence. Après avoir effectué la correction décalage de temps et la similitude pour chaque matrice secondaire calculée, il est nécessaire d'empiler plans ST. Ceci peut être réalisé en utilisant une moyenne algébrique ou une moyenne géométrique.
La moyenne algébrique peut être calculée comme suit .
où Fi est le plan ST calculé pour une matrice secondaire et N est le nombre de matrices secondaires utilisées pour calculer le résultat de tirs multiples final .
où Fi est le plan ST calculé pour une matrice secondaire et N est le nombre de matrices secondaires utilisées pour calculer le résultat de tirs multiples final .
moyenne géométrique peut être calculée comme suit .
Puisque l'effet d'un mauvais plan ST est moindre si la moyenne algébrique est utilisée plutôt que la moyenne géométrique, la moyenne algébrique est la forme préférée d'empilement.
Puisque l'effet d'un mauvais plan ST est moindre si la moyenne algébrique est utilisée plutôt que la moyenne géométrique, la moyenne algébrique est la forme préférée d'empilement.
I1 existe deux autres moyens pour calculer le résultat de tirs multiples final. Une carte de similitude de plans ST est définie comme le rapport entre l'énergie de cohérence Ci(S,t) l'énergie totale Ti(S,t), en conséquence
Il est possible de calculer le plan final en calculant le plan ST de tirs multiples final selon l'équation 12 ou en empilant d'abord l'énergie de cohérence et l'énergie totale séparément puis, en calculant les tirs multiples finaux comme le rapport entre l'empilement de l'énergie de cohérence et de l'énergie totale. Dans la seconde technique, ce qui suit est calculé
La comparaison des résultats entre les deux procédés ne présente aucune différence significative dans la précision du résultat et le premier procédé préféré, car il est plus facile à mettre en oeuvre avec les procédés STC actuellement utilisés.
Il est possible de calculer le plan final en calculant le plan ST de tirs multiples final selon l'équation 12 ou en empilant d'abord l'énergie de cohérence et l'énergie totale séparément puis, en calculant les tirs multiples finaux comme le rapport entre l'empilement de l'énergie de cohérence et de l'énergie totale. Dans la seconde technique, ce qui suit est calculé
La comparaison des résultats entre les deux procédés ne présente aucune différence significative dans la précision du résultat et le premier procédé préféré, car il est plus facile à mettre en oeuvre avec les procédés STC actuellement utilisés.
Un traitement efficace de données à tirs multiples nécessite que ce qui suit soit déterminé 1. Le nombre de matrices secondaires utilisees pour calculer les tirs multiples.
2. L'identification des traces d'entrée requises pour chaque matrice secondaire comprenant le résultat de tirs multiples à une profondeur donnée.
3. La sélection automatique du récepteur de référence pour chaque matrice secondaire selon configuration utilisée.
Pour un outil avec R récepteurs, Rsb est le nombre de récepteurs dans une matrice secondaire. Ce nombre est choisi par l'utilisateur et est lié à la résolution à obtenir. Par définition, sa limite inférieure est Le nombre de change pas le traitement réel mais seule l'interprétation physique du résultat du traitement. est le nombre maximum de matrices secondaires utilisées pour calculer les tirs multiples. Cette valeur est égale au nombre de profondeurs de mesure dont l'outil déplacé pour effectuer les mesures tirs multiples et peut être définie par Nsb = R - Rsb + 1<B>(16)</B> Pour savoir quels sont les récepteurs nécessaires a utiliser pour calculer les tirs multiples, une "matrice géométrique" est mise en oeuvre. En calculant des tirs multiples, il est nécessaire de modifier le recepteur de référence pour chaque matrice secondaire lors du calcul du plan ST pour prendre en compte la correction de décalage de temps. Pour Nsb matrices secondaires de Rsb récepteurs, la matrice géométrique a pour dimensions G=G (Nsb, Rsb) et est définie par
Le premier récepteur de la première matrice secondaire rn est égal à la valeur de la matrice secondaire nécessaire pour calculer les tirs multiples.
Le premier récepteur de la première matrice secondaire rn est égal à la valeur de la matrice secondaire nécessaire pour calculer les tirs multiples.
conséquence, rn = Nsb, le deuxième récepteur de cette matrice secondaire est r,+, et ainsi de suite. Chaque colonne de la matrice peut être calculée par ri = rn + i , i= 0, Rsb où i est le numéro du récepteur à l'intérieur de matrice secondaire. Le récepteur de référence est toujours le même que le premier récepteur de la première matrice secondaire. Ainsi, si cette valeur est calculée au début du calcul, tous les paramètres pour le traitement STC sont fixés pour ce calcul et il n'est pas nécessaire de redéfinir le récepteur de référence pour chaque matrice secondaire. De plus, la matrice géométrique permet une identification facile au niveau de chaque trame des récepteurs qui sont activés et désactives à cette profondeur. En conséquence, il est simplement nécessaire de multiplier l'état du récepteur par le vecteur courant pour savoir quels récepteurs doivent etre utilisés, ce qui simplifie la gestion des données. Il est simplement nécessaire savoir si un récepteur est actif ou non pour une trame donnée, car les données de tous les récepteurs sont calculées au début et les résultats utilisés permutation circulaire de la colonne de la matrice pour chaque nouvelle trame. Dans l'exemple de la figure 2, la matrice géométrique est comme suit
Le récepteur numéro 4 est le récepteur de référence pour chaque matrice secondaire.
Le récepteur numéro 4 est le récepteur de référence pour chaque matrice secondaire.
La matrice de tirs multiples M prolonge ce concept pour identifier la contribution de chaque récepteur dans la matrice d'outils aux résultats de tirs multiples. Cette matrice a pour dimensions M (Nsb, Nrcvr+Nsb-1) , comme représentait dans 1 'équation 19 ci-dessous
Chaque colonne de cette matrice représente une position de l'outil en fonction de la profondeur et rmin est le récepteur le plus proche de l'émetteur rmax est le plus éloigné. Chaque élément de la matrice devient 1 lorsqu'il est nécessaire d'utiliser la forme d'onde enregistrée par le récepteur approprié et 0 dans le cas contraire. La connaissance de cette matrice permet une identification facile des traces récepteurs requises et ainsi, facilite la gestion des données. Il est également possible de savoir d'après cette matrice quelle matrice secondaire dépend de chaque trace de récepteur.
Chaque colonne de cette matrice représente une position de l'outil en fonction de la profondeur et rmin est le récepteur le plus proche de l'émetteur rmax est le plus éloigné. Chaque élément de la matrice devient 1 lorsqu'il est nécessaire d'utiliser la forme d'onde enregistrée par le récepteur approprié et 0 dans le cas contraire. La connaissance de cette matrice permet une identification facile des traces récepteurs requises et ainsi, facilite la gestion des données. Il est également possible de savoir d'après cette matrice quelle matrice secondaire dépend de chaque trace de récepteur.
La matrice géométrique et la matrice de tirs multiples sont calculées une fois au début du traitement et pendant le traitement, les données d'entrée sont décalées avant d'exécuter chaque calcul.
La figure 4 représente un organigramme des étapes de base d'un procédé selon l'invention comme décrit ci- dessus.
La présente invention trouve des applications dans le domaine de la caractérisation de formations souterraines entourant des forages, par exemple dans l'industrie pétrolière et du gaz.
Claims (8)
1. Procédé pour traiter des enregistrements de diagraphie acoustique obtenus à partir d'un outil de diagraphie comportant une matrice de R récepteurs, comprenant (i) la sélection d'une dimension de matrice secondaire de Rsb récepteurs, avec Rsb < R (ii) le calcul du nombre de matrices secondaires Nsb possibles à partir de la matrice de R récepteurs ; (iii) l'identification d'un récepteur Rref dans la matrice pour servir de récepteur de référence pour chaque matrice secondaire ; (iv) le traitement de données dans le plan lenteur-temps pour chaque matrice secondaire par rapport au récepteur de référence Rref ; et (v) 'empilement des données traitées dans le plan lenteur-temps pour toutes les matrices secondaires.
2. Procédé selon la revendication 1 dans lequel l'empilement est réalisé en utilisant une moyenne algébrique ou une moyenne géométrique.
3. Procédé selon la revendication ou 2, dans lequel l'étape d'empilement comprend le calcul de cohérences lenteur-temps pour chaque matrice secondaire suivi l'empilement ou le calcul de l'énergie de cohérence et de l'énergie totale pour les matrices secondaires utilisées, puis l'obtention rapport pour déterminer les données de cohérence lenteur-temps empilées
4. Procédé selon la revendication 1, 2 ou , dans lequel Nb comprend le nombre maximum possible de matrices secondaires à partir de la matrice de R récepteurs selon Nsb = R - Rsb + 1.
5. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre le calcul d'une matrice géométrique de dimension par Rsb, qui identifie les récepteurs particuliers à considérer à chaque profondeur.
6. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre le calcul d'une matrice de tirs multiples de dimensions Nsb par R + Nsb-1, qui identifie la contribution de chaque récepteur dans l'outil pour une matrice secondaire donnée.
7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le récepteur de référence Rref est identifié par la relation Rref = R - Rsb + 1 où Rref est le numéro du récepteur de réference dans la matrice.
8. Procédé selon l'une quelconque des revendications précédentes, dans lequel le nombre de récepteurs dans chaque matrice secondaire est choisi pour déterminer la résolution de la mesure.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GBGB9923461.9A GB9923461D0 (en) | 1999-04-01 | 1999-10-06 | Processing sonic waveform measurements from array borehole logging tools |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2799552A1 true FR2799552A1 (fr) | 2001-04-13 |
FR2799552B1 FR2799552B1 (fr) | 2006-12-08 |
Family
ID=10862106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0012733A Expired - Fee Related FR2799552B1 (fr) | 1999-10-06 | 2000-10-05 | Traitement de mesures de formes d'ondes acoustiques provenant d'outils de diagraphie de sondage en matrice |
Country Status (5)
Country | Link |
---|---|
US (1) | US6459993B1 (fr) |
CA (1) | CA2322128C (fr) |
FR (1) | FR2799552B1 (fr) |
GB (1) | GB2357841B (fr) |
NO (1) | NO333705B1 (fr) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2372322B (en) * | 2000-10-16 | 2003-04-16 | Schlumberger Holdings | Method for determining formation slowness particularly adapted for measurement while drilling applications |
US6845325B2 (en) * | 2001-11-08 | 2005-01-18 | Schlumberger Technology Corporation | Global classification of sonic logs |
EP3064242A1 (fr) | 2003-04-28 | 2016-09-07 | Advanced Circulatory Systems Inc. | Ventilateur et procédés de traitement de traumatisme crânien et d'hypotension |
US7075215B2 (en) * | 2003-07-03 | 2006-07-11 | Pathfinder Energy Services, Inc. | Matching layer assembly for a downhole acoustic sensor |
US6995500B2 (en) * | 2003-07-03 | 2006-02-07 | Pathfinder Energy Services, Inc. | Composite backing layer for a downhole acoustic sensor |
US7036363B2 (en) * | 2003-07-03 | 2006-05-02 | Pathfinder Energy Services, Inc. | Acoustic sensor for downhole measurement tool |
US7513147B2 (en) * | 2003-07-03 | 2009-04-07 | Pathfinder Energy Services, Inc. | Piezocomposite transducer for a downhole measurement tool |
US6957572B1 (en) | 2004-06-21 | 2005-10-25 | Schlumberger Technology Corporation | Apparatus and methods for measuring mud slowness in a borehole |
US7587936B2 (en) | 2007-02-01 | 2009-09-15 | Smith International Inc. | Apparatus and method for determining drilling fluid acoustic properties |
US7672193B2 (en) * | 2007-05-21 | 2010-03-02 | Schlumberger Technology Corporation | Methods and systems for processing acoustic waveform data |
US8117907B2 (en) | 2008-12-19 | 2012-02-21 | Pathfinder Energy Services, Inc. | Caliper logging using circumferentially spaced and/or angled transducer elements |
US12032110B2 (en) | 2016-05-25 | 2024-07-09 | Schlumberger Technology Corporation | Elastic parameter estimation |
WO2018080450A1 (fr) * | 2016-10-25 | 2018-05-03 | Halliburton Energy Services, Inc. | Détermination de lenteur d'onde de volume de formation rocheuse à résolution améliorée à partir d'ondes guidées par le trou de forage |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4543648A (en) * | 1983-12-29 | 1985-09-24 | Schlumberger Technology Corporation | Shot to shot processing for measuring a characteristic of earth formations from inside a borehole |
EP0660136A2 (fr) * | 1993-12-20 | 1995-06-28 | Schlumberger Limited | Traitement en fond du puits de l'information de form d'onde sonore |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU509996B2 (en) * | 1976-05-17 | 1980-06-05 | Schlumberger Technology B.V. | Compensating fcr atilt in boreholes |
AU730715B2 (en) * | 1996-10-09 | 2001-03-15 | Baker Hughes Incorporated | Method of obtaining improved geophysical information about earth formations |
US6453240B1 (en) | 1999-04-12 | 2002-09-17 | Joakim O. Blanch | Processing for sonic waveforms |
-
2000
- 2000-09-22 GB GB0023310A patent/GB2357841B/en not_active Expired - Fee Related
- 2000-10-03 US US09/678,454 patent/US6459993B1/en not_active Expired - Lifetime
- 2000-10-04 CA CA002322128A patent/CA2322128C/fr not_active Expired - Fee Related
- 2000-10-05 NO NO20005014A patent/NO333705B1/no not_active IP Right Cessation
- 2000-10-05 FR FR0012733A patent/FR2799552B1/fr not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4543648A (en) * | 1983-12-29 | 1985-09-24 | Schlumberger Technology Corporation | Shot to shot processing for measuring a characteristic of earth formations from inside a borehole |
EP0660136A2 (fr) * | 1993-12-20 | 1995-06-28 | Schlumberger Limited | Traitement en fond du puits de l'information de form d'onde sonore |
Non-Patent Citations (1)
Title |
---|
KAI HSU AND SHU-KONG CHANG: "Multiple-shot processing of array sonic waveforms", GEOPHYSICS, vol. 52, no. 10, 10 October 1987 (1987-10-10), pages 1376 - 1390, XP002348113 * |
Also Published As
Publication number | Publication date |
---|---|
GB0023310D0 (en) | 2000-11-08 |
NO20005014D0 (no) | 2000-10-05 |
FR2799552B1 (fr) | 2006-12-08 |
US6459993B1 (en) | 2002-10-01 |
CA2322128A1 (fr) | 2001-04-06 |
CA2322128C (fr) | 2008-09-30 |
GB2357841A (en) | 2001-07-04 |
NO333705B1 (no) | 2013-08-26 |
GB2357841B (en) | 2001-12-12 |
NO20005014L (no) | 2001-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2534696A1 (fr) | Procede et dispositif pour l'etude indirecte des ondes de cisaillement | |
EP0176408B1 (fr) | Procédé et dispositif de diagraphie pour l'inspection acoustique d'un sondage muni d'un tubage | |
EP1173781B1 (fr) | Methode et dispositif de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques | |
FR3068479A1 (fr) | Compensation de la reponse angulaire pour vsp dsp | |
EP0354112B1 (fr) | Méthode pour obtenir un modèle représentatif d'un milieu hétérogène et notamment du sous-sol | |
FR2646513A1 (fr) | Procede et dispositif de diagraphie pour l'inspection acoustique d'un sondage muni d'un tubage | |
EP2253970B1 (fr) | Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway | |
FR2782389A1 (fr) | Procede et dispositif pour corriger les effets des deplacements du navire dans les mesures de sismique maritime | |
CA1325840C (fr) | Procede d'acquisition de donnees sismiques et dispositif de mise en oeuvre du procede | |
EP3129809B1 (fr) | Focalisation sismique adaptative | |
EP0555148B1 (fr) | Procédé pour déterminer une signature en champ lointain d'une pluralité de sources sismiques | |
FR2532058A1 (fr) | Appareil et procede de diagraphie acoustique et procede de reduction du bruit du aux ondes de compression et de stoneley | |
WO2007021857A2 (fr) | Procede de representation d'etirement d'ondelettes dans des donnees sismiques | |
FR2647554A1 (fr) | Procede et dispositif de mesure de l'anisotropie azimutale d'une formation | |
FR2799552A1 (fr) | Traitement de mesures de formes d'ondes acoustiques provenant d'outils de diagraphie de sondage en matrice | |
GB2125550A (en) | Underwater seismic testing | |
FR2478324A1 (fr) | Procede de prospection et d'exploration sismiques | |
FR2792731A1 (fr) | Diagraphie de puits par induction electromagnetique et des moyens de mise en oeuvre | |
GB2440276A (en) | Processing seismic data acquired using three or more over/under streamers and/or sources | |
FR2575299A1 (fr) | Procede d'affichage de donnees acoustiques de diagraphies de sondages par production d'empilages de temps de parcours | |
FR2743896A1 (fr) | Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement | |
EP0053525A1 (fr) | Procédé d'exploration sismique par la technique du profil sismique vertical et installation pour sa mise en oeuvre | |
EP0786671B1 (fr) | Procédé de prospection sismique marine au moyen d'un couple de capteurs hydrophone et géophone | |
CA2182111C (fr) | Methode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu | |
EP0944847B1 (fr) | Procede de traitement sismique et notamment procede de prospection sismique 3d mettant en oeuvre une migration des donnees sismiques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 16 |
|
ST | Notification of lapse |
Effective date: 20170630 |