FR2599508A1 - Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches - Google Patents

Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches Download PDF

Info

Publication number
FR2599508A1
FR2599508A1 FR8607964A FR8607964A FR2599508A1 FR 2599508 A1 FR2599508 A1 FR 2599508A1 FR 8607964 A FR8607964 A FR 8607964A FR 8607964 A FR8607964 A FR 8607964A FR 2599508 A1 FR2599508 A1 FR 2599508A1
Authority
FR
France
Prior art keywords
medium
wave
detector
source
propagation
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
Application number
FR8607964A
Other languages
English (en)
Other versions
FR2599508B1 (fr
Inventor
Charles Naville
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CGG SA
Original Assignee
Compagnie Generale de Geophysique SA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Compagnie Generale de Geophysique SA filed Critical Compagnie Generale de Geophysique SA
Priority to FR8607964A priority Critical patent/FR2599508B1/fr
Priority to US07/057,279 priority patent/US4789969A/en
Publication of FR2599508A1 publication Critical patent/FR2599508A1/fr
Application granted granted Critical
Publication of FR2599508B1 publication Critical patent/FR2599508B1/fr
Expired legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H5/00Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/04Measuring characteristics of vibrations in solids by using direct conduction to the detector of vibrations which are transverse to direction of propagation
    • G01H1/08Amplitude

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

L'INVENTION PERMET LA MESURE DE L'ANISOTROPIE DE PROPAGATION D'UNE ONDE TRANSVERSE ENTRE DEUX POINTS DE REFERENCE DONNES D'UN MILIEU M NON ISOTROPE, NOTAMMENT UNE ONDE DE CISAILLEMENT DANS UNE ROCHE STRATIFIEE OU FISSUREE. SELON L'INVENTION : ON DISPOSE LE LONG D'UN MEME RAYON VECTEUR UNE SOURCE S ET DEUX DETECTEURS G, G, LES DETECTEURS ETANT PLACES A L'ENDROIT DES POINTS DE REFERENCE RESPECTIFS; ON EXCITE LA SOURCE PAR UN SIGNAL D'EXCITATION PRODUISANT UNE ONDE TRANSVERSE DANS LE MILIEU; ON RECUEILLE LES SIGNAUX DE MESURE RESPECTIFS RESULTANTS PRODUITS PAR CHACUN DES DETECTEURS; A PARTIR DU SIGNAL D'EXCITATION ET DES SIGNAUX DE MESURE, ON DETERMINE LA FONCTION DE TRANSFERT DU MILIEU LE LONG DE CHACUN DES DEUX TRAJETS SOURCE-DETECTEUR RESPECTIFS; ON DEDUIT LA FONCTION DE TRANSFERT DIFFERENTIELLE DU MILIEU ENTRE LES DEUX POINTS DE REFERENCE. ON PEUT EGALEMENT PREVOIR D'AUTRES CONFIGURATIONS, NOTAMMENT DEUX SOURCES ET UN DETECTEUR, OU APPLIQUER LE PROCEDE A LA SISMIQUE REFLEXION AVEC DES ONDES A INCIDENCE OBLIQUE OU MEME A INCIDENCE NULLE.

Description

La présente invention concerne un procédé de mesure de l'anisotropie de
propagation d'une onde transverse entre deux points de référence donnés d'un
milieu non isotrope à étudier.
L'invention peut être également appliquée à la mesure de l'anisotropie de réflexion d'une onde transverse en un point de référence donné, propre à créer une réflexion non isotrope d'onde transverse, dans
un milieu sensiblement isotrope.
On décrira plus particulièrement l'application de ce procédé à la prospection géophysique par mesure de l'anisotropie de propagation ou de réflexion des ondes de cisaillement dans les roches; néanmoins, cette application n'est en aucune façon limitative, 15 et le procédé de l'invention peut être appliqué à de nombreux domaines, très différents, dans la mesure ou il peut être employé chaque fois que l'on est en présence d'une onde progressive comportant des composantes transverses: tel est le cas notamment de l'onde électro20 magnétique, ou d'une onde acoustique se propageant dans un milieu solide (ou assimilable à un solide), ce qui permet d'appliquer l'invention aux techniques radar,
sonar et d'imagerie médicale.
En d'autres termes, le procédé de l'invention est susceptible de s'appliquer, de façon générale, chaque fois que l'on est en présence d'un milieu qui propage de façon non isotrope les ondes de cisaillement ou des ondes analogues, ou même d'un milieu isotrope, lorsque
c'est une réflexion qui produit l'anisotropie.
On va tout d'abord expliquer le phénomène d'anisotropie de propagation d'une onde transverse en prenant comme exemple le cas de la propagation d'une onde de cisaillement dans une roche, en référence aux
figures 1 et 2.
Le phénomène de l'anisotropie de propagation d'une onde de cisaillement souvent désignée "onde S" (pour "shear") - est typique des matériaux fissurés ou stratifiés, et l'étude de la propagation des ondes de cisaillement dans ces matériaux donne des indications
précieuses sur la structure du milieu traversé.
Sur la figure 1, on a illustré une onde incidente de cisaillement SH/SV ayant par exemple une orientation à 45 par rapport à l'horizontale, donc présentant des composantes à la fois dans la direction verticale V et dans la direction horizontale H. Cette onde incidente pénètre dans un milieu stratifié M, anisotrope, qui va créer un phénomène de discrimination 15 des ondes de cisaillement ("shear-wave splitting") le long du trajet de propagation: les composantes verticale et horizontale vont se propager avec des polarisations, des vitesses et des atténuations différentes dans le milieu M, produisant ainsi, à partir d'une même onde 20 incidente SV/SH, deux ondes résultantes SV et SH, polarisées respectivement verticalement et horizontalement, présentant un décalage temporel T. Ainsi, la forme d'onde, qui est la somme vectorielle des deux composantes SV
et SH, va graduellement se modifier au fur et à mesure 25 de sa propagation au travers du milieu anisotrope.
De la même façon, la figure 2 illustre la propagation d'une onde de cisaillement SN/SE se propageant le long d'un rayon vecteur sensiblement vertical, et ayant des composantes à la fois dans une direction 30 Nord N et Est E. Cette onde incidente SN/SE pénètre dans un milieu M anisotrope présentant une fissuration orientée par exemple dans le sens Nord-Sud. Corsme dans le cas précédent, les deux composantes Nord SN et Est SE de l'onde incidente vont se propager à des vitesses différentes en raison de l'anisotropie, modifiant graduellement la forme d'onde et créant un décalage temporel T. La composante se propageant le plus 5 rapidement est celle qui est polarisée parallèlement aux faces des fissures (c'est-àdire dans le sens Nord-Sud, dans l'exemple choisi); en outre, cette composante ne sera que faiblement atténuée. En revanche,
la composante orientée Ouest-Est va se propager plus 10 lentement et subir une atténuation plus forte.
L'étude de l'anisotropie du milieu comporte donc l'étude des paramètres suivants: la recherche des orientations correspondant respectivement à la plus faible et à la plus grande vitesse 15 de propagation (recherche des directions principales d'anisotropie), ce qui permet d'évaluer l'orientation des fissures, qui seront perpendiculaires à la direction de compression horizontale minimale; le décalage temporel T entre les composantes corres20 pondant aux deux directions principales d'anisotropie; la mesure de l'atténuation des différentes composantes, notamment leur atténuation différentielle
dans les deux directions principales d'anisotropie.
Ceci fournira des indications sur la nature du milieu traversé et permettra de quantifier l'importance de l'anisotropie, fournissant ainsi une indication de la densité de fracturation ou de la
porosité du milieu traverse.
La présente invention a précisément pour objet de proposer une méthode permettant de déterminer ces différents paramètres, tout en utilisant un appareillage classique (générateur d'onde de cisaillement, géophones omnidirectionnels deux ou trois composantes,...) et par application de méthodes sismiques
habituelles (sismique réflexion ou sismique de puits).
L'invention se base essentiellement sur l'effet différentiel de propagation des ondes de cisaillement le long de deux trajets de mesure le long d'un même rayon vecteur, c'est-à-dire le long de deux
trajets dont l'un recouvre complètement l'autre.
L'application du procédé de l'invention 10 suppose que, hormis les réflexions, la propagation s'effectue de façon sensiblement linéaire, ce qui est généralement le cas rencontré en sismique dans un
milieu rocheux.
Plus précisément, selon l'invention, pour 15 mesurer l'anisotropie de propagation de l'onde transverse entre deux points de référence donnés du milieu non isotrope étudier, selon un premier mode de mise en oeuvre: on dispose le long d'un même rayon vecteur une source 20 et deux détecteurs, les détecteurs étant placés à l'endroit des points de référence respectifs, on excite la source par un signal d'excitation produisant une onde transverse dans le milieu, on recueille les signaux de mesure respectifs résul25 tants produits par chacun des détecteurs, à partir du signal d'excitation et des signaux de mesure, on détermine la fonction de transfert du milieu le long de chacun des deux trajets source-détecteur respectifs, 30. on déduit la fonction de transfert différentielle du
milieu entre les deux points de référence.
Cette hypothèse s'applique particulièrement à la sismique de puits, la direction du rayon vecteur étant la direction de forage du puits (technique dite
"VSP" ou "PSV"; Profil Sismique Vertical).
Dans un second mode de mise en oeuvre: on dispose le long d'un même rayon vecteur deux sources et un détecteur, les sources étant placées à l'endroit des points de référence respectifs, on excite les sources par les signaux d'excitation respectifs produisant des ondes transverses dans le milieu, 10.on recueille les signaux de mesure respectifs résul- tants produits par le détecteur, à partir des signaux d'excitation et des signaux de mesure, on détermine la fonction de transfert du milieu le long de chacun des deux trajets source15 détecteur respectifs, on déduit la fonction de transfert différentielle du
milieu entre les deux points de référence.
Dans ce dernier cas, les sources peuvent être aussi bien des sources réelles (générateurs d'ondes de cisaillement) que des sources virtuelles
situées chacune au point de réflexion d'une onde incidente produite par une source réelle distante.
Les sources réelles distantes peuvent notamment être des sources d'ondes de compression, de sorte que l'onde incidente est une onde longitudinale prcueisant des composantes d'onde transverse après réflexion. Ceci permet de s'affranchir des effets d'anisotropie de propagation d'onde transverse du
milieu traversé entre les sources réelles et les sources 30 virtuelles.
Cette dernière hypothèse s'applique tout particulièrement à la sismique réflexion, avec une configuration dans laquelle le détecteur est situé à
distance des sources, les ondes incidentes et réfléchies se propageant obliquement par rapport à la verticale.
Dans un troisième mode de mise en oeuvre: 5. on dispose, en alignement de propagation avec les deux points de référence, une source et un détecteur, on excite la source par un signal d'excitation produisant une onde transverse dans le milieu, de manière que cette onde se propage le long d'un même 10 rayon vecteur dans un premier sens jusqu'aux deux points de référence, cette onde se réfléchissant ensuite à chacun des points de référence et se propageant, en retour, le long dudit rayon vecteur en sens inverse jusqu'au détecteur, 15. on recueille les signaux de mesure respectifs résultants produits par le détecteur, on détermine la fonction de transfert du milieu le long de chacun des deux trajets source-point de référence-détecteur respectifs, 20. on en déduit la fonction de transfert différentielle
du milieu entre les deux points de référence.
Cette hypothèse s'applique particulièrement au cas de la sismique réflexion de puits avec source et détecteur situés sur une même verticale (en un même point ou non), configuration dans laquelle les ondes se propagent à la verticale et subissent des réflexions en
des endroits situés à la verticale de ce point.
Comme indiqué plus haut, l'invention s'applique également à la mesure de l'anisotropie de réflexion 30 d'une onde transverse en un point de référence donné,
propre à créer une réflexion non isotrope d'onde transverse, dans un milieu sensiblement isotrope.
Dans ce cas, selon l'invention:
2599.08
on dispose, en alignement de propagation avec le point de référence, une source et un détecteur, on excite la source par un signal d'excitation produisant une onde transverse dans le milieu, de ma5 nière que cette onde se propage le long d'un rayon vecteur dans un premier sens jusqu'au point de référence, o elle se réfléchit et se propage, en retour, le long dudit rayon vecteur en sens inverse jusqu'au détecteur, 10. on recueille les signaux de mesure respectifs résultants produits par le détecteur, on détermine la fonction de transfert du réflecteur
au point de référence.
D'autres caractéristiques et avantages de
l'invention seront exposés dans la description détaillée
qui va suivre des différents modes de mise en oeuvre, en référence aux dessins annexés, sur lesquels: - les figures 1 et 2, déjà commentées, illustrent le phénomène d'anisotropie de propagation d'une onde 20 transverse dans un milieu non isotrope, - la figure 3 illustre le premier mode de mise en oeuvre de l'invention, - la figure 4 illustre le second mode de mise en oeuvre de l'invention, dans leauel les sources sont 25 des sources réelles, - la figure 5 illustre une variante du mode de mise en oeuvre précédent, dans lequel les sources d'ondes de cisaillement sont des sources virtuelles, - la figure 6 illustre le troisième mode de mise 30 en oeuvre de l'invention, - la figure 7 illustre les mesures successives à opérer pour obtenir une image tridimensionnelle de
l'anisotropie du milieu.
259950O
Sur la figure 3, la référence S désigne une source d'onde de cisaillement, et G1 et G2 deux géophones, les organes S, G1 et G2 étant alignés le long
d'un trajet de propagation sensiblement rectiligne.
Ce trajet est par exemple vertical, l'hypothèse d'étude étant alors celle de la figure 2. Dans cette configuration, l'axe Z est alors la verticale (ou la direction de forage du puits), les organes S, G1 et G2 étant, par exemple,placés à différents niveaux dans 10 un même puits, et le milieu M à étudier étant celui
situé entre les deux géophones G et G2.
Par ailleurs, S et G1 peuvent éventuellement
être confondus.
L'onde produite par la source S sera une 15 onde de cisaillement ayant ses composantes SX et Sy dans le plan horizontal (X, Y) normal à la direction de propagation Z. Le signal exprimé sous forme matricielle, reçu par le géophone G1 sera:
(I () = F(S,G1) (
F(S,G1) étant la fonction de transfert (opérateur r.atriciel de dimensions (2,2)) entre S et G1. De la même façon, le signal capté par le géophone G2 sera:
(II) | = -F(G1,G2) * F(S,G1) (
Y2/S P(G1,G2) étant la fonction de transfert du milieu M entre G1 et G2, et * étant l'opérateur de convolution, d'o: (IIbis) y2
= F (G1,G2)
X1 Y1 Il est ainsi possible de déterminer la fonction de transfert F(G1, G2) 2 partir des mesures
(X1,Y1) et (X2,Y2) relevées aux.points G1 et G2.
Les signaux de source et de géophones sont
bien entendu exprimés dans le même repère géographique 10 fixe (X, Y).
La fonction de transfert F(G1,G2), qui est la fonction à déterminer pour évaluer l'anisotropie du milieu M entre G1 et G2, est de la forme suivante (en transformée de Fourier):
(III) F(G1,G2) = K P
1 exp (ioT1 O o Q2 exp(iT2 Q2 exp (iw-r2) p qui peut également se mettre sous la forme: (IIIbis) F(G1,G2) = K' P-1 1 0o -Q2 exp i (2- 1 p P étant la matrice de passage du repère fixe (X,Y) au repère (orthogonal ou non) des deux ondes principales
de cisaillement.
(T2-T1) étant le décalage temporel entre les composantes 5 correspondant aux deux directions principales d'anisotropie entre les points G1 et G2, Q2/Q1 étant l'atténuation relative des composantes selon les deux directions principales d'anisotropie,
K et K' étant des filtres isotropes de transmission.
Les deux directions principales d'anisotropie sont déterminées en utilisant pour G et G2 des géophones omnidirectionnels à réponse selon trois axes (ou, tout au moins, selon deux axes, si l'on peut orienter de façon précise le géophone par rapport à la direction 15 du rayon vecteur), en recherchant la direction selon laquelle le signal arrive le premier et subit la plus
faible atténuation.
Ayant ainsi déterminé P, il est possible de tirer de la fonction de transfert déduite des mesures 20 relevées aux points G1 et G2 les paramètres T 2-T1 et Q2/Q1 permettant de quantifier l'anisotropie du milieu M.. Dans le cas d'une conversion P/S en transmission (conversion d'une onde de pression (longitudinale) 25 en une onde de cisaillement (transversale) à l'interface entre deux milieux différents), une source virtuelle peut être assimilée à G1, là o la conversion se produit. Commela direction de polarisation de l'onde P avant
conversion est connue, il est possible de calculer le 30 tenseur F(G1,G2) à partir de l'expression (IIbis).
La figure 4 illustre un second mode de mise en oeuvre du procédé: dans ce cas, au lieu de disposer une source et deux géophones le long d'un même rayon vecteur, on dispose deux sources S2 et S1 et un seul géophone G, le but étant d'étudier l'anisotropie du
milieu M situé entre les deux sources S1 et s2.
Les autres hypothèses sont les mêmes que dans le cas de la figure 3. Le géophone G va tout d'abord recevoir une première onde de cisaillement en provenance de la source S1 (dont on connaît l'orientation et l'amplitude des composantes), qui va produire la réponse 10 suivante: Xl sx
(IV) ( = F(S1,G)
Egalité qui permet de calculer F(S1,G) si la source n'est pas polarisée linéairement dans une seule direction. De la même façon, le géophone va recevoir 15 l'onde émise par la source S2, qui aura traversé le milieu M entre S2 et Sl, ainsi que le milieu entre S et G. Le signal recueilli pourra être exprimé par: X2 sx S2 (V) = F(S1,G) * F(S2,S1) s 2/ Sy 2 Compte tenu du fait que F(S1,G) est connu, cette expression peut également se mettre sous la 20 forme: X2 s X2 (Vbis) F (S1G) F(S2,S1)
Y2 %SY2
On peut en déduire l'opérateur recherché
F(S2,S1), qui est lui aussi de la forme (III) en première approximation.
La configuration correspondant à ce mode de mise en oeuvre (c'est-à-dire deux sources et un géophone) s'applique de façon très avantageuse à la sismique réflexion, dans laquelle on étudie des profils bidimensionnels en ondes converties réfléchies P/S, c'est-à-dire o l'on émet des ondes de pression (ondes 10 longitudinales), sans composantes de cisaillement, qui se transforment après - et seulement près - réflexion
en ondes comportant des composantes de cisaillement.
On utilisera de préférence la technique dite "WLP" ("Wide Line Profiling" ou "profils à extension latérale") ou-les études de sismique à trois dimensions
pour émettre un signal virtuel de source en onde S, polarisé linéairement dans plusieurs directions successives.
Cette configuration a été illustrée fiqure 5, 20 les traits fins représentant des ondes de pression (ondes longitudinales), et les traits épais des ondes
de cisaillement (ondes transverses).
Les deux ondes de pression créées par les sources P1 et P2 se propagent parallèlement jusqu'aux 25 points S'1 et SI 2 o se produit une réflexion engendrant des composantes d'onde de cisaillement. On. se ramène ainsi au cas de la figure 4, les sources réelles S1 et S2 étant remplacées par deux sources virtuelles S' et S'2, et le milieu à étudier étant le milieu M i 2' entre les deux points de réflexion, c'est-à-dire en fait entre deux interfaces de stratification de la roche, d'o l'intérêt de cette configuration. Le fait que, entre les sources réelles et les sources virtuelles, les ondes ne se propagent que selon un mode longitudinal, permet d'éliminer toute influence de l'anisotropie du milieu entre P1 et S' d'une part et P2 et S' d'autre part, anisotropie d'autant plus gênante que les deux ondes ne se propagent pas selon le même rayon vecteur; en revanche, l'influence de l'anisotropie du milieu sera déterminante le long du trajet entre S2' S' et G, et on a vu plus 15 haut comment il est possible de déterminer et d'éliminer l'influence de l'anisotropie du milieu entre S' et G. La figure 6 illustre un troisième mode de mise en oeuvre de l'invention, dans lequel l'onde de 20 cisaillement se propage à incidence (verticale ou non) normale aux réflecteurs, de manière à étudier les réflexions se produisant en des points R1 et R2 situés sur le même rayon double o sont situés à la fois la source S et le géophone G. La source S et le géophone G peuvent être, comme illustré, situés en un même point; ils peuvent tout aussi bien se trouver en des points différents, mais toujours situés en alignement avec R et R2 -(par
exemple à deux niveaux différents d'un puits foré en 30 direction de R1 et R2).
Comme il n'est plus possible, du fait de l'incidence normale, d'obtenir une conversion onde de pression/onde de cisaillement, la source S doit être une source d'onde de cisaillement. Les deux trajets aller-retour seront les trajets S R1G d'une part, et
S R2 G, d'autre part.
Connaissant le signal (Sx, Sy) émis par la 5 source et le signal recueilli par l'unique géophone,
il s'agit de déterminer les caractéristiques d'anisotropie du milieu M situé entre les deux points de réflexion Ri et R2.
Le signal recueilli en G à la suite du premier trajet S R1 G de l'onde de cisaillement est: xi SX
(VI) = F(R1,S) * F(S,R1) * K
Y1 De la même façon, le signal recueilli à la suite du second trajet S R2 G est: x 2 sx
(VII) =F(R1,S)*F(R2,R1)*F(R1,R2)*F(S,R1)*K2
Y; s Y2 Syt K1 et K2 sont des opérateurs à une dimension représentatifs de la stratification, c'est-à-dire de paramètres sans incidence sur la polarisation de l'onde: réflexions (supposées isotropes) en R1 et R2 et transmissions (également supposées isotropes) de S à R1, et de R1 à G.
2599508 15
L'opérateur recherché est le produit de convolution:
(VIII)
F(R2,R1) * F (R1,R2)
Cette expression étant de la forme: /Q2 exp (2i-1) o (IX) K P 1 o Q2 exp (2iïw2) p forme analogue à l'expression (III) évaluée dans les 5 cas précédents, dès lors que l'on aura pu déterminer
au préalable la fonction de transfert F(S,R1).
Si, entre R1 et S, les deux directions principales d'anisotropie sont à angle droit dans chacun des milieux traversés, on aura alors:
(X) F(S,R1) = TF (R1,S)
Ce qui permet de calculer facilement F(S,R1) connaissant le produit:
(XI) F(R1,S) * TF (S,R1)
déterminé à partir de deux mesures avec la source S dans deux directions différentes. Cette détermination
peut par exemple être obtenue par la méthode de factori15 sation de Cayley.
Dans le cas contraire, il faut procéder pas à pas et déterminer au préalable l'anisotropie dans chacun des n milieux traversés entre S et R1, de façon à établir le produit: n Qlkexp(iwTlk) o
(XII) F(S,R) _ K P-1 -P
k = 1 k kk k O Q2 kexp (iWT2k) afin de pouvoir déduire F(R2,R1) dans l'expression IX. Dans la configuration de la figure 6 avec R1 et R2 confondus (c'est-à-dire un trajet aller-retour le long d'un même rayon), il est possible d'appliquer l'invention à la détermination de l'anisotropie due à la 10 réflexion au point R1/R2, la propagation des ondes de cisaillement étant alors supposée sensiblement isotrope
dans le milieu traversé entre S et R1/R2.
Si le milieu est anisotrope, on peut également déterminer au préalable son anisotropie grâce à des réflexions fixes isotropes connues (méthode utilisable
en particulier en technique radar).
Pour terminer, on peut formuler les remarques suivantes: il est possible de multiplier les directions pour 20 avoir une image tridimensionnelle de l'anisotropie, comme illustré figure 7: en se plaçant par exemple dans la configuration avec deux sources et un géophone, on fait tourner le rayon vecteur de manière à faire propager l'onde de cisaillement dans plu25 sieurs directions dans le milieu M. dans les exemples décrits plus haut, on a envisagé des sources et des géophones places dans un puits ou bien en surface. Cette différence est sans aucune incidence sur la mise en oeuvre du procédé, et il est possible d'envisager toutes les combinaisons possibles dans une même configuration (sources en
puits et géophones en surface, etc.).
une calibration préalable des géophones et des sources peut être nécessaire, pour déterminer à la fois l'anisotropie propre à l'instrument et l'aniso'O tropie de couplage instrument-roche (notamment dans le cas d'un couplage par un outil de puits appliquant la source ou le géophone contre un point de la
paroi cylindrique du puits).

Claims (7)

REVENDICATIONS
1. Un procédé de mesure de l'anisotropie de propagation d'une onde transverse entre deux points de référence donnés d'un milieu (M) non isotrope à étudier, notamment un procédé de prospection géophysique par mesure de l'anisotropie de propagation des ondes de cisaillement dans les roches, caractérisé en ce que: on dispose le long d'un même rayon vecteur une source 10 (S) et deux détecteurs (G1, G2), les détecteurs étant placés à l'endroit des points de référence respectifs, on excite la source par un signal d'excitation produisant une onde transverse dans le milieu, on recueille les signaux de mesure respectifs résul15 tants produits par chacun des détecteurs, à partir du signal d'excitation et des signaux de mesure, on détermine la fonction de transfert du milieu le long de chacun des deux trajets source-détecteur respectifs, 20. on déduit la fonction de transfert différentielle du
milieu entre les deux points de référence.
2. Un procédé de mesure de l'anisotropie de propagation d'une onde transverse entre deux points de référence donnés d'un milieu fM) non isoro-e à 9tu25 dGer notamment un procédé de prospection géophysique par mesure de l'anisotropie de propagation des ondes de cisaillement dans les roches, caractérisé en ce que: 30. on dispose le long d'un même rayon vecteur deux sources (S1, S2; S'1, S'2) et un détecteur (G), les
sources étant placées à l'endroit des points de réfé-
rence respectifs, À on excite les sources par les signaux d'excitation respectifs produisant des ondes transverses dans le milieu, 5. on recueille les signaux de mesure respectifs résultants produits par le détecteur, à partir des signaux d'excitation et des signaux de mesure, on détermine la fonction de transfert du milieu le long de chacun des deux trajets source-détec10 teur respectifs, on déduit la fonction de transfert différentielle du
milieu entre les deux points de référence.
3. Un procédé selon la revendication 2, dans lequel les sources (S1, S2) sont des sources réelles. 15
4. Un procédé selon la revendication 2, ran lequel le sources (S1, S'i2) sont des sources virtuelles situées chacune au point de réflexion d'une onde
incidente produite par une source réelle (P'1, P2) distante.
5. Un procédé selon la revendication 4, dans
lequel l'onde incidente est une onde longitudinale produisant des composantes d'onde transversale après réflexion.
6. Un procédé de mesure de l'anisotropie de pro25 pagation d'une onde transverse entre deux points de référence donnés, propres à créer des réflexions sensiblement isotropes d'onde transverse, d'un. milieu (M) non isotrope à étudier, notamment un procédé de prospection géophysique 30 par mesure de l'anisotropie de propagation des ondes de cisaillement dans les roches, caractérisé en ce que: on dispose, en alignement de propagation avec les deux points de référence (R1, R2), une source (S) et un détecteur (G), on excite la source par un signal d'excitation produisant une onde transverse dans le milieu, de manière que cette onde se propage le long d'un même rayon vecteur dans un premier sens jusqu'aux deux points de référence, cette onde se réfléchissant ensuite à chacun des points de référence et se propageant, en retour, le long dudit rayon vecteur en sens inverse 10 jusqu'au détecteur, on recueille les signaux de mesure respectifs résultants produits par le détecteur, on détermine la fonction de transfert du milieu le long de chacun des deux trajets source-point de réfé15 rence-détecteur respectifs, on en déduit la fonction de transfert différentielle
du milieu entre les deux points de référence.
7. Un procédé de mesure de l'anisotropie de 20 réflexion d'une onde transverse en un point de réfé20 rence donné, propre à créer une réflexion non isotrope d'onde transverse, dans un milieu sensiblement isotrope, notamment un procédé de prospection géophysique par mesure de l'anisotropie de réflexion des ondes de cisaillement dans les roches, 25 caractérisé en ce que: on dispose,en alignement de propagation avec le point de référence, une source (S) et un détecteur (G), on excite la source par un signal d'excitation pro30 duisant une onde transverse dans le milieu, de manière que cette onde se propage le long d'un rayon vecteur dans un premier sens jusvu'au point de référence, o. elle se réfléchit et se propage, en retour, le long dudit rayon vecteur en sens inverse jusqu'au détecteur, on recueille les signaux de mesure respectifs résul5 tants produits par le détecteur, on détermine la fonction de transfert du réflecteur
au point de référence.
FR8607964A 1986-06-03 1986-06-03 Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches Expired FR2599508B1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
FR8607964A FR2599508B1 (fr) 1986-06-03 1986-06-03 Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches
US07/057,279 US4789969A (en) 1986-06-03 1987-06-02 Method of measuring the anisotropy of propagation or reflection of a transverse wave, particularly a method of geophysical prospecting by measurement of the anisotropy of propagation or of reflection of shear waves in rocks

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR8607964A FR2599508B1 (fr) 1986-06-03 1986-06-03 Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches

Publications (2)

Publication Number Publication Date
FR2599508A1 true FR2599508A1 (fr) 1987-12-04
FR2599508B1 FR2599508B1 (fr) 1989-08-04

Family

ID=9335938

Family Applications (1)

Application Number Title Priority Date Filing Date
FR8607964A Expired FR2599508B1 (fr) 1986-06-03 1986-06-03 Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches

Country Status (2)

Country Link
US (1) US4789969A (fr)
FR (1) FR2599508B1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2648562A1 (fr) * 1989-06-16 1990-12-21 Geophysique Cie Gle Procede et systeme de determination de l'anisotropie azimutale dans un domaine d'un milieu materiel

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4903244A (en) * 1984-07-20 1990-02-20 Standard Oil Company (Now Amoco Corporation) Multisource multireceiver method and system for geophysical exploration
US4933913A (en) * 1986-10-30 1990-06-12 Amoco Corporation Method of seismic surveying for resolving the effects of formation anisotropy in shear wave reflection seismic data
US4912979A (en) * 1987-10-14 1990-04-03 Amoco Corporation Method and apparatus for detecting and measuring elastic anisotropy
US5027332A (en) * 1987-10-14 1991-06-25 Amoco Corporation Method for geophysical exploration
US4907206A (en) * 1988-07-29 1990-03-06 Amoco Corporation Method for estimating formation lithology
US4995008A (en) * 1989-12-27 1991-02-19 Exxon Production Research Company Method of using a circularly-polarized source to characterize seismic anisotropy
US5540093A (en) * 1994-11-15 1996-07-30 Western Atlas International Method for optimizing the alignment of a horizontal borehole relative to the strike of rock-layer stress planes
FR2759172B1 (fr) * 1997-02-05 1999-03-05 Inst Francais Du Petrole Methode de traitement de donnees sismiques de puits multi-composantes orientees
US5982706A (en) * 1997-03-04 1999-11-09 Atlantic Richfield Company Method and system for determining normal moveout parameters for long offset seismic survey signals
US5996415A (en) * 1997-04-30 1999-12-07 Sensys Instruments Corporation Apparatus and method for characterizing semiconductor wafers during processing
US6019000A (en) * 1997-11-20 2000-02-01 Sensys Instruments Corporation In-situ measurement of deposition on reactor chamber members
FR2831962B1 (fr) * 2001-11-08 2004-06-25 Geophysique Cie Gle Procede de traitement sismique, notamment pour la compensation de birefringence sur des traces sismiques
US6842400B2 (en) * 2001-12-18 2005-01-11 Halliburton Energy Services, Inc. Acoustic logging apparatus and method
US6772067B2 (en) * 2001-12-19 2004-08-03 Halliburton Energy Services, Inc Acoustic logging apparatus and method for anisotropic earth formations
US20040244490A1 (en) * 2003-01-22 2004-12-09 Steven Turner Method and apparaus for measuring acoustic properties of materials to determine other physical propertiesfield of the invention
US7175599B2 (en) * 2003-04-17 2007-02-13 Brigham And Women's Hospital, Inc. Shear mode diagnostic ultrasound
CN103697999B (zh) * 2013-12-30 2015-10-14 中国科学院武汉岩土力学研究所 一种高应力硬岩tbm施工隧道微震波速实时获取方法
CN113295776B (zh) * 2021-05-28 2022-04-08 东南大学 一种瑞利波信号的高效检测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2295633A1 (fr) * 1974-12-16 1976-07-16 Electronique Appliquee Dispositif operateur d'ondes acoustiques transversales
US4080836A (en) * 1977-03-07 1978-03-28 Rockwell International Corporation Method of measuring stress in a material
EP0169076A2 (fr) * 1984-07-20 1986-01-22 Amoco Corporation Procédé de mesure des formations terrestres souterraines

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3003577A (en) * 1957-06-05 1961-10-10 Texaco Inc Apparatus for seismic exploration
US3622965A (en) * 1970-06-18 1971-11-23 Western Geophysical Co Method of geophysical prospecting by measuring the attenuation of seismic waves in the earth
US3946598A (en) * 1974-07-11 1976-03-30 Robin M. Towne And Associates, Inc. Method and apparatus for determining the dynamic parameters of soil in situ
US4138894A (en) * 1974-12-16 1979-02-13 L'electronique Appliquee Transversal acoustic wave operating device
FR2524645A1 (fr) * 1982-04-06 1983-10-07 Thomson Csf Polarimetre acoustique
US4648039A (en) * 1983-12-30 1987-03-03 Schlumberger Technology Corporation Compressional/shear wave separation in vertical seismic profiling
US4713968A (en) * 1986-05-22 1987-12-22 Mobil Oil Corporation Method and apparatus for measuring the mechanical anisotropy of a material

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2295633A1 (fr) * 1974-12-16 1976-07-16 Electronique Appliquee Dispositif operateur d'ondes acoustiques transversales
US4080836A (en) * 1977-03-07 1978-03-28 Rockwell International Corporation Method of measuring stress in a material
EP0169076A2 (fr) * 1984-07-20 1986-01-22 Amoco Corporation Procédé de mesure des formations terrestres souterraines

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
1982 ULTRASONICS SYMPOSIUM PROCEEDINGS, éditor B.R.McAvoy, Sponsored by the IEEE Group on Sonics & Ultrasonics, 27-29 octobre 1982, vol. 2, San Diego, CA., pages 861-863, IEEE; J.ROUGE: "Acoustical polarimeter" *
JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, vol. 79, no. 1, janvier 1986, pages 84-90, Acoustical Society of America, Woodbury, New York, US; H.W.JONES et al.: "The detection of liquids and viscoelastic substances trapped under solid surfaces" *
ULTRASONICS, vol. 23, no. 5, septembre 1985, pages 215-222, Butterworth & Co. (Publishers) Ltd., Guildford, Surrey, GB; D.R.ALLEN et al.: "Ultrasonic SH wave velocity in textured aluminium plates" *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2648562A1 (fr) * 1989-06-16 1990-12-21 Geophysique Cie Gle Procede et systeme de determination de l'anisotropie azimutale dans un domaine d'un milieu materiel

Also Published As

Publication number Publication date
FR2599508B1 (fr) 1989-08-04
US4789969A (en) 1988-12-06

Similar Documents

Publication Publication Date Title
FR2599508A1 (fr) Procede de mesure de l'anisotropie de propagation ou de reflexion d'une onde transverse, notamment procede de prospection geophysique par mesure de l'anisotropie de propagation ou de reflexion des ondes de cisaillement dans les roches
EP0889331B1 (fr) Méthode pour modéliser en 3D l'impédance d'un milieu hétérogène
EP2369369B1 (fr) Procédé de surveillance d'un site de stockage géologique de gaz par inversion stratigraphique de données sismiques
US20040230389A1 (en) Determination of waveguide parameters
FR2831961A1 (fr) Methode de traitement de donnees sismiques de puits en amplitude preservee absolue
EP1669778B1 (fr) Méthode pour déterminer de l'information spéculaire dans l'imagerie sismique avant sommation
CA2786411A1 (fr) Procede et dispositif d'acquisition de donnees sismiques marines
FR2647554A1 (fr) Procede et dispositif de mesure de l'anisotropie azimutale d'une formation
Taweesintananon et al. Distributed acoustic sensing for near-surface imaging using submarine telecommunication cable: A case study in the Trondheimsfjord, Norway
EP0124442B1 (fr) Procédé et dispositif d'holographie acoustique utilisant un faisceau ultrasonore limité dans l'espace
EP0797780B1 (fr) Methode de traitement de traces sismiques reflexion enregistrees pour des deports variables
CA2226389C (fr) Methode de traitement de donnees sismiques de puits multi-composantes orientees
FR2743896A1 (fr) Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement
EP0752114A1 (fr) Methode d'exploration d'un milieu a tectonique complexe
EP0944847B1 (fr) Procede de traitement sismique et notamment procede de prospection sismique 3d mettant en oeuvre une migration des donnees sismiques
EP1226456B1 (fr) Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties
EP0610123B1 (fr) Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques
EP0909961B1 (fr) Procédé de séparation d'ondes en sismique de puits pour des acquisitions de type à déport croissant
FR2800170A1 (fr) Perfectionnement aux procedes de traitement sismique
Poletto* et al. Dual seismic fields from distributed acoustic sensors
EP0458947B1 (fr) Procede et dipsositif d'acquisition de donnees sismiques de puits selon deux directions opposees
FR2675590A1 (fr) Procede de traitement pour l'obtention d'une section somme a offsets nuls.
Nasralla Distributed Acoustic Sensing (DAS) As a New Tool For Subsurface Imaging
Boëlle et al. Deep target Pre Stack Time Migration and Migration Velocity Analysis using converted PS data
CA2365336A1 (fr) Methode et algorithme d'utilisation des ondes de surface

Legal Events

Date Code Title Description
ST Notification of lapse