WO2004086557A2 - Procede pour determiner des reponses impulsionnelles d’un milieu vis-à-vis de la transmission d’ondes entre differents points - Google Patents

Procede pour determiner des reponses impulsionnelles d’un milieu vis-à-vis de la transmission d’ondes entre differents points Download PDF

Info

Publication number
WO2004086557A2
WO2004086557A2 PCT/FR2004/000684 FR2004000684W WO2004086557A2 WO 2004086557 A2 WO2004086557 A2 WO 2004086557A2 FR 2004000684 W FR2004000684 W FR 2004000684W WO 2004086557 A2 WO2004086557 A2 WO 2004086557A2
Authority
WO
WIPO (PCT)
Prior art keywords
medium
points
waves
signals
equal
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.)
Ceased
Application number
PCT/FR2004/000684
Other languages
English (en)
Other versions
WO2004086557A3 (fr
Inventor
Matthias Fink
Julien De La Gorgue De Rosny
Claire Julia
Thomas Folegot
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.)
Centre National de la Recherche Scientifique CNRS
Universite Paris Diderot Paris 7
Original Assignee
Centre National de la Recherche Scientifique CNRS
Universite Paris Diderot Paris 7
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 Centre National de la Recherche Scientifique CNRS, Universite Paris Diderot Paris 7 filed Critical Centre National de la Recherche Scientifique CNRS
Priority to US10/550,429 priority Critical patent/US7656748B2/en
Priority to DE602004025690T priority patent/DE602004025690D1/de
Priority to AT04742298T priority patent/ATE459014T1/de
Priority to EP04742298A priority patent/EP1620749B1/fr
Publication of WO2004086557A2 publication Critical patent/WO2004086557A2/fr
Publication of WO2004086557A3 publication Critical patent/WO2004086557A3/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/895—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum
    • G01S15/8952—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum using discrete, multiple frequencies
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046—Techniques for image enhancement involving transmitter or receiver

Definitions

  • the present invention relates to methods for determining impulse responses of a medium with respect to the transmission of waves between different points.
  • the invention relates to a method for determining impulse responses of a medium with respect to the transmission of waves between different points, method comprising:
  • Responses' pulse can be determined to wave focusing purposes in the medium, for example for the purposes of imaging medium or selective communication between different points of the medium, particularly in the context of the implementation of the "DORT" method (Decomposition of the Operator
  • the process of acquiring impulse responses between the points in question in the medium is however relatively slow, insofar as it is carried out sequentially, for the different emission points i one after the other, as described for example in document OA-02/32 316 within the framework of a method of focusing acoustic waves other than the "DORT" method above.
  • This slowness lengthens the implementation time of the process, and when the medium is changing, it does not allow the impulse responses to be obtained quickly enough to be able to be validly used in this medium.
  • a process of the kind in question is characterized in that during the step
  • the signals ei (t) are transmitted simultaneously by the said N emission points, these signals ei (t) having a duration T and each being a sum of n substantially monochromatic elementary signals, of the same amplitude and of frequencies respective f 0 , i + k. ⁇ f, where f 0 , i is a predetermined frequency specific to point i, k is an integer between 0 and n, n is an integer at least equal to 2 and ⁇ f is a step predetermined frequency, the respective frequencies f 0 , i specific to the different points i being distinct from each other and included in a frequency band of width ⁇ f, and in that during step (c), each impulse response hij (t) from a correlation signal between the signal ei (t) emitted at point i and the signal rj (t) picked up at point j.
  • the impulse responses are calculated simultaneously between N transmission points and M reception points (possibly confused with the emission points), which accelerates the process of acquiring impulse responses.
  • the waves transmitted in the medium between the emission points and the reception points are acoustic waves; during step (a), the medium in which the waves are emitted is reverberant; the frequency step ⁇ f is less than or equal to 1 / ⁇ , where ⁇ is the temporal dispersion of the medium (that is to say the duration of the signal picked up after emission of a pulse signal in the medium); the frequency step ⁇ f is substantially equal to 1 / ⁇ , where ⁇ is the time dispersion of the medium; - the duration T is at least equal to N / ⁇ f; the duration T is at least equal to N. ⁇ , where ⁇ is the time dispersion of the medium; the elementary signals have random phases; the waves are transmitted with a certain bandwidth, the frequencies fOi include a minimum frequency fO and the number n is determined so that the frequency band between fO and f0 + [(n + 1). ⁇ f] substantially covers said bandwidth; the reception points are confused with the emission points.
  • Figure 1 is a schematic view of a device suitable for implementing the present invention, in one embodiment.
  • the method according to one embodiment of the invention is implemented by a device 1 which is adapted to transmit and receive waves in a medium 2 which is homogeneous or not, generally reverberant f which can be solid, liquid or gaseous.
  • the waves in question are preferably acoustic and the description below will only refer to such acoustic waves, but said waves could possibly be of another nature, for example of electromagnetic nature.
  • the device 1 comprises for example a computer 3, or any other similar control device, which controls the emissions of acoustic waves and their listening in the medium 2.
  • the computer 3 can for example control a dedicated electronic device 4 which can in particular comprise an electronic central unit CPU associated with a central memory MEM, this central unit CPU controlling itself N channels (N being an integer at least equal to 2) each connected to a transducer Tl-TN making it possible to transmit and pick up signals acoustics in environment 2. All channels are synchronized with each other.
  • Each measurement channel of the device 4 can for example comprise an El-EN sampler connected to the corresponding transducer T1-TN, and each sampler can itself be connected to a memory associated with a central unit, jointly referenced Cl-CN.
  • the device 1 can be used in particular for the purpose of imaging the medium 2, or even to establish selective communication with a point of said medium, or the like, for example using the "D.O.R.T.” method.
  • the computer 3 controls the device 1 so that it successively follows the following steps:
  • all the signals ei (t) are transmitted simultaneously during step (a) by the N transducers, but these signals transmitted are orthogonal to one another so that the information specific to each emission point can then be separated from each other others in the signals received rj (t).
  • each signal emitted ei (t) is equal to the sum of n substantially monochromatic elementary signals, of the same amplitude and of respective frequencies fo, i + k. ⁇ f, where: f 0 / i is a predetermined frequency specific to the transducer i , k is an integer between 0 and n, n is an integer at least equal to 2, - and ⁇ f is a predetermined frequency step, advantageously less than or equal to 1 / ⁇ and preferably substantially equal to 1 / ⁇ , where ⁇ is an average value of the temporal dispersion of the medium 2, that is to say the duration of the signal picked up after emission of a pulse signal (Dirac function) in the medium 2.
  • the respective frequencies f 0 , i specific to the different points i are distinct from each other and included in a frequency band of width ⁇ f: advantageously, these frequencies fo, i are separated two by two by a difference ⁇ f / N.
  • BP is the passband with which the acoustic waves are emitted and received by the transducers
  • fO is called the lowest of the frequencies fOi
  • fO and n can advantageously be determined so that the frequency band of the signals elementary (band between fO and f0 + [(n + 1). ⁇ f]) substantially covers said passband.
  • the elementary signals composing the signals ei (t) have random phases. All the signals ei (t) have the same duration T, which is advantageously at least equal to N / ⁇ f, and preferably at least equal to N. ⁇ .
  • each impulse response hij (t) is calculated from a correlation signal between the signal ei (t) emitted at point i and the signal rj (t) received at point j.
  • this correlation signal is window by means of a door function ⁇ (t) of width 1 / ⁇ f
  • impulse responses hij (t) can advantageously be determined by the formula:
  • reception points j at which the acoustic wave is picked up are merged with the emission points and constituted by the Tl-TN transducers, but the reception points in question may if necessary be constituted by a second array of transducers (not shown), distinct from the array of Tl-TN transducers.
  • the acoustic waves which are generated by the emission of the signals ei (t) are therefore picked up at M reception points, M being a non-zero integer, which M reception points can be either confused with the N reception points.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Reduction Or Emphasis Of Bandwidth Of Signals (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

Procédé pour déterminer des réponses impulsionnelles d'un milieu (2) vis-à-vis de la transmission d'ondes entre différents points (T1-TN), on émet des ondes dans le milieu en générant en chaque point i des signaux ei(t) constitués chacun de n signaux élémentaires de fréquences respectives espacées deux à deux d'un pas deltaf et distinctes des fréquences des signaux élémentaires correspondant aux autres points, puis on capte des signaux rj (t) en des points j après transmission desdites ondes dans le milieu, et on calcule chaque réponse impulsionnelle hij(t) à partir d'un signal de corrélation entre le signal ei(t) émis au point i et le signal rj(t) capté au point j.

Description

Procédé pour déterminer des réponses impulsionnelles d' n milieu vis-à-vis de la transmission d'ondes entre différents points .
La présente invention est relative aux procédés pour déterminer des réponses impulsionnelles d'un milieu vis-à-vis de la transmission d' ondes entre différents points .
Plus particulièrement, l'invention concerne un procédé pour déterminer des réponses impulsionnelles d'un milieu vis-à-vis de la transmission d'ondes entre différents points, procédé comprenant :
(a) au moins une étape d'émission au cours de laquelle on émet des ondes dans le milieu en générant des signaux ei(t) à partir d'un nombre N de points d'émission appartenant au milieu, où N est un nombre entier au moins égal à 2 et i est un indice compris entre 1 et N qui désigne l'un desdits N points d'émission,
(b) au moins une étape de réception au cours de laquelle on capte des signaux rj(t) à partir desdites ondes après transmission dans ledit milieu, en un nombre M de points de réception appartenant au milieu (ces points de réception peuvent être confondus avec les points d'émission), où M est un entier naturel non nul et j est un indice compris entre 1 et M qui désigne l'un desdits M points de réception,
(c) et au moins une étape de détermination desdites réponses impulsionnelles ij (t) entre chaque point d'émission i et chaque point de réception j à partir des signaux émis ei(t) et captés rj (t) .
Les réponses ' impulsionnelles peuvent être déterminées aux fins de focalisation d'ondes dans le milieu, par exemple aux fins d'imagerie du milieu ou de communication sélective entre différents points du milieu, notamment dans le cadre de la mise en œuvre de la méthode "D.O.R.T." (Décomposition de l'Opérateur
Retournement Temporel) décrite notamment par Prada et al.
(C. Prada and M. Fink, "Eigenmodes of the time reversai operator: a solution to sélective focusing in multiple- target média." Wave Motion , 20, pp 151-163 (1994)).
Le processus d'acquisition des réponses impulsionnelles entre les points en question du milieu est toutefois relativement lent, dans la mesure où il est effectué séquentiellement, pour les différents points d'émission i les uns à la suite des autres, comme décrit par exemple dans le document O-A-02/32 316 dans le cadre d'un procédé de focalisation d'ondes acoustiques autre que la méthode "D.O.R.T." susmentionnée.
Cette lenteur allonge le temps de mise en œuvre du procédé, et lorsque le milieu est changeant, elle ne permet pas d'obtenir les réponses impulsionnelles assez vite pour pouvoir être valablement utilisées ensuite dans ce milieu.
La présente invention a notamment pour but de pallier ces inconvénients. A cet effet, selon l'invention, un procédé du genre en question est caractérisé en ce qu'au cours de l'étape
(a), on fait émettre simultanément les signaux ei(t) par lesdits N points d'émission, ces signaux ei(t) ayant une durée T et étant chacun une somme de n signaux élémentaires sensiblement monochromatiques, de même amplitude et de fréquences respectives f0,i+k.δf , où f0,i est une fréquence prédéterminée propre au point i, k est un nombre entier compris entre 0 et n, n est un nombre entier au moins égal à 2 et δf est un pas fréquentiel prédéterminé, les fréquences f0,i respectives propres aux différents points i étant distinctes les unes des autres et comprises dans une bande de fréquences de largeur δf, et en ce qu'au cours de l'étape (c) , on calcule chaque réponse impulsionnelle hij (t) à partir d'un signal de corrélation entre le signal ei(t) émis au point i et le signal rj (t) capté au point j.
Grâce à ces dispositions, on calcule simultanément les réponses impulsionnelles entre N points d' émission et M points de réception (éventuellement confondus avec les points d'émission), ce qui accélère le processus d' acquisition des réponses impulsionnelles .
Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : - les fréquences f0/i respectives propres aux différents points i sont séparées deux à deux par un écart δf/N ; au cours de l'étape (c) , ledit signal de corrélation est fenêtre au moyen d'une fonction porte π(t) de largeur 1/δf ; au cours de l'étape (c) , les réponses impulsionnelles hij (t) sont déterminées par la formule :
Figure imgf000005_0001
les ondes transmises dans le milieu entre les points d'émission et les points de réception sont des ondes acoustiques ; au cours de l'étape (a), le milieu où sont émises les ondes est réverbérant ; le pas fréquentiel δf est inférieur ou égal à 1/τ, où τ est la dispersion temporelle du milieu (c'est-à- dire la durée du signal capté après émission d'un signal impulsionnel dans le milieu) ; le pas fréquentiel δf est sensiblement égal à 1/τ, où τ est la dispersion temporelle du milieu ; - la durée T est au moins égale à N/δf ; la durée T est au moins égale à N.τ, où τ est la dispersion temporelle du milieu ; les signaux élémentaires présentent des phases aléatoires ; les ondes sont émises avec une certaine bande passante, les fréquences fOi comprennent une fréquence minimale fO et le nombre n est déterminé pour que la bande de fréquence comprise entre fO et f0+ [ (n+1) . δf] recouvre sensiblement ladite bande passante ; les points de réception sont confondus avec les points d'émission.
D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante de deux de ses modes de réalisation, donnés à titre d'exemples non limitatifs, en regard des dessins joints .
Sur le dessin, la figure 1 est une vue schématique d'un dispositif adapté pour mettre en œuvre la présente invention, dans un mode de réalisation.
Dans l'exemple non limitatif représenté sur le dessin, le procédé selon un mode de réalisation de l'invention est mis en œuvre par un dispositif 1 qui est adapté pour émettre et recevoir des ondes dans un milieu 2 homogène ou non, généralement réverbérant f qui peut être solide, liquide ou gazeux. Les ondes en question sont de préférence acoustiques et la description ci-après ne fera référence qu'à de telles ondes acoustiques, mais lesdites ondes pourraient éventuellement être d'une autre nature, par exemple de nature électromagnétique.
Le dispositif 1 comporte par exemple un ordinateur 3, ou tout autre dispositif de commande similaire, qui commande les émissions d'ondes acoustiques et leur écoute dans le milieu 2. A cet effet, l'ordinateur 3 peut par exemple commander un dispositif électronique dédié 4 qui peut notamment comporter une unité centrale électronique CPU associée à une mémoire centrale MEM, cette unité centrale CPU commandant elle-même N voies (N étant un entier au moins égal à 2) reliées chacune à un transducteur Tl-TN permettant d'émettre et de capter des signaux acoustiques dans le milieu 2. Toutes les voies sont synchronisées entre elles.
Chaque voie de mesure du dispositif 4 peut comporter par exemple un échantillonneur El-EN relié au transducteur correspondant Tl-TN, et chaque échantillonneur peut être lui-même relié à une mémoire associée à une unité centrale, référencées conjointement Cl-CN.
Le dispositif 1 peut être utilisé notamment aux fins d'imagerie du milieu 2, ou encore pour établir une communication sélective avec un point dudit milieu, ou autre, par exemple en utilisant la méthode "D.O.R.T."
(Décomposition de l'Opérateur Retournement Temporel) décrite notamment par Prada et al. (C. Prada and M. Fink,
"Eigenmodes of the time reversai operator: a solution to sélective focusing in multiple-target média." Wave Motion ,
20, pp 151-163 (1994) ) .
Dans tous les cas, il est nécessaire de déterminer les réponses impulsionnelles entre les transducteurs Tl-TN.
A cet effet, l'ordinateur 3 commande le dispositif 1 de façon qu'il suive successivement les étapes suivantes :
(a) au moins une étape d'émission au cours de laquelle on émet des ondes dans le milieu en générant des signaux ei(t) à partir des N transducteurs Tl-TN (ces N transducteurs peuvent le cas échéant ne constituer qu'une partie des transducteurs reliés au dispositif 4, auquel cas les réponses impulsionnelles relatives à l'ensemble des transducteurs seraient acquises en plusieurs passes successives, ces passes restant toutefois en nombre très inférieur au nombre total de transducteurs), i étant un indice compris entre 1 et N qui désigne l'un desdits N transducteurs ;
(b) au moins une étape de réception au cours de laquelle les N transducteurs Tl-TN captent des signaux rj (t) à partir desdites ondes après transmission dans ledit milieu, où j est un indice compris entre 1 et N désignant l'un desdits transducteurs ;
(c) et au moins une étape de détermination desdites réponses impulsionnelles hij (t) entre chaque point d'émission i et chaque point de réception j à partir des signaux émis ei(t) et captés rj (t) .
Pour permettre d' acquérir en une seule fois tous les signaux nécessaires à la détermination des réponses impulsionnelles concernant les N transducteurs susmentionnés (c'est-à-dire tout ou partie des transducteurs reliés au dispositif 4, comme expliqué ci- dessus), tous les signaux ei(t) sont émis simultanément au cours de l'étape (a) par les N transducteurs, mais ces signaux émis sont orthogonaux entre eux de façon que les informations propres à chaque point d' émission puissent ensuite être séparées les unes des autres dans les signaux captés rj (t) .
Ainsi, chaque signal émis ei(t) est égal à la somme de n signaux élémentaires sensiblement monochromatiques, de même amplitude et de fréquences respectives fo,i+k.δf , où : f0/i est une fréquence prédéterminée propre au transducteur i, k est un nombre entier compris entre 0 et n, n est un nombre entier au moins égal à 2, - et δf est un pas fréquentiel prédéterminé, avantageusement inférieur ou égal à 1/τ et de préférence sensiblement égal à 1/τ, où τ est une valeur moyenne de la dispersion temporelle du milieu 2, c'est-à-dire la durée du signal capté après émission d'un signal impulsionnel (fonction de Dirac) dans le milieu 2.
Les fréquences f0,i respectives propres aux différents points i sont distinctes les unes des autres et comprises dans une bande de fréquences de largeur δf : avantageusement, ces fréquences fo,i sont séparées deux à deux par un écart δf/N. De plus, si BP est la bande passante avec laquelle les ondes acoustiques sont émises et reçues par les transducteurs, et si on nomme fO la plus faible des fréquences fOi, alors fO et n peuvent avantageusement être déterminés pour que la bande des fréquences des signaux élémentaires (bande comprise entre fO et f0+ [ (n+1) . δf] ) recouvre sensiblement ladite bande passante.
Avantageusement, les signaux élémentaires composant les signaux ei(t) présentent des phases aléatoires. Tous les signaux ei(t) ont une même durée T, qui est avantageusement au moins égale à N/δf, et de préférence au moins égale à N.τ.
Au cours de l'étape (c) , on calcule chaque réponse impulsionnelle hij (t) à partir d'un signal de corrélation entre le signal ei(t) émis au point i et le signal rj (t) capté au point j .
Dans l f exemple considéré ici, ce signai de corrélation vaut: cij(t) = - θ)dθ (cette intégrale est faite par
Figure imgf000009_0001
exemple sur une durée au moins égale à T, dont le début coïncide avec le début du signal rj (t) ) .
Avantageusement, ce signal de corrélation est fenêtre au moyen d'une fonction porte π(t) de largeur 1/δf
(valant 1 pendant une période de 1/δf et 0 en dehors de cette période) , auquel cas les réponses impulsionnelles hij (t) peuvent avantageusement être déterminées par la formule :
Figure imgf000009_0002
On notera que, dans l'exemple particulier représenté sur le dessin et décrit ci-dessus, les points de réception j auxquels on capte l'onde acoustique sont confondus avec les points d'émission et constitués par les transducteurs Tl-TN, mais les points de réception en question peuvent le cas échéant être constitués par un second réseau de transducteurs (non représenté) , distinct du réseau de transducteurs Tl-TN.
Dans le cas. général, les ondes acoustiques qui sont générées par l'émission des signaux ei(t) sont donc captées en M points de réception, M étant un nombre entier non nul, lesquels M points de réception peuvent être soit confondus avec les N points d'émission (auquel cas M = N) , soit non totalement confondus avec les points d'émission, soit entièrement distincts des points d'émission.

Claims

REVENDICATIONS
1. Procédé pour déterminer des réponses impulsionnelles d'un milieu (2) vis-à-vis de la transmission d'ondes entre différents points (Tl-TN), procédé comprenant :
(a) au moins une étape d'émission au cours de laquelle on émet des ondes dans le milieu (2) en générant des signaux ei(t) à partir d'un nombre N de points d'émission (Tl-TN) appartenant au milieu, où N est un nombre entier au moins égal à 2 et i est un indice compris entre 1 et N qui désigne l'un desdits N points d'émission,
(b) au moins une étape de réception au cours de laquelle on capte des signaux rj (t) à partir desdites ondes après transmission dans ledit milieu, en un nombre M de points de réception (Tl-TN) appartenant au milieu, où M est un entier naturel non nul et j est un indice compris entre 1 et M qui désigne l'un desdits M points de réception,
(c) et au moins une étape de détermination desdites réponses impulsionnelles hij (t) entre chaque point d'émission i et chaque point de réception j à partir des signaux émis ei(t) et captés rj (t) , caractérisé en ce qu'au cours de l'étape (a), on fait émettre simultanément les signaux ei(t) par lesdits N points d'émission (Tl-TN), ces signaux ei(t) ayant une durée T et étant chacun une somme de n signaux élémentaires sensiblement monochromatiques, de même amplitude et de fréquences respectives f0,ι+k.δf , où f0,i est une fréquence prédéterminée propre au point i, k est un nombre entier compris entre 0 et n, n est un nombre entier au moins égal à 2 et δf est un pas fréquentiel prédéterminé, les fréquences f0,i respectives propres aux différents points i étant distinctes et comprises dans une bande de fréquences de largeur δf, et en ce qu'au cours de l'étape (c) , on calcule chaque réponse impulsionnelle hij (t) à partir d'un signal de corrélation entre le signal ei (t) émis au point i et le signal rj (t) capté au point j .
2. Procédé selon la revendication 1, dans lequel les fréquences f0,i respectives propres aux différents points i sont séparées deux à deux par un écart δf/N.
3. Procédé selon la revendication 1 ou la revendication 2, dans lequel au cours de l'étape (c) , ledit signal de corrélation est fenêtre au moyen d'une fonction porte π(t) de largeur 1/δf.
4. Procédé selon la revendication 3, dans lequel au cours de l'étape (c) , les réponses impulsionnelles hij (t) sont déterminées par la formule : hij(t) ≈ TI(t).jei(θ -t).rj(θ)dθ .
5. Procédé selon l'une quelconque des revendications précédentes, dans lequel les ondes transmises dans le milieu entre les points d'émission et les points de réception sont des ondes acoustiques .
6. Procédé selon l'une quelconque des revendications précédentes, dans lequel au cours de l'étape
(a), le milieu où sont émises les ondes est réverbérant.
7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le pas fréquentiel δf est inférieur ou égal à 1/τ, où τ est la dispersion temporelle du milieu.
8. Procédé selon la revendication 7 , dans lequel le pas fréquentiel δf est sensiblement égal à 1/τ, où τ est la dispersion temporelle du milieu.
9. Procédé selon l'une quelconque des revendications précédentes, dans lequel la durée T est au moins égale à N/δf.
10. Procédé selon l'une quelconque des revendications précédentes, dans lequel la durée T est au moins égale à N.τ, où τ est la dispersion temporelle du milieu.
11. Procédé selon l'une quelconque des revendications précédentes, dans lequel les signaux élémentaires présentent des phases aléatoires.
12. Procédé selon l'une quelconque des revendications précédentes, dans lequel les ondes sont émises avec une certaine bande passante, les fréquences fOi comprennent une fréquence minimale fO et le nombre n est déterminé pour que la bande de fréquence comprise entre fO et f0+[ (n+1) . δf] recouvre sensiblement ladite bande passante.
13. Procédé selon l'une quelconque des revendications précédentes, dans lequel les points de réception sont confondus avec les points d'émission (Tl- TN) .
PCT/FR2004/000684 2003-03-24 2004-03-19 Procede pour determiner des reponses impulsionnelles d’un milieu vis-à-vis de la transmission d’ondes entre differents points Ceased WO2004086557A2 (fr)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US10/550,429 US7656748B2 (en) 2003-03-24 2004-03-19 Method of determining impulse responses from a medium in relation to the transmission of waves between different points
DE602004025690T DE602004025690D1 (de) 2003-03-24 2004-03-19 Verfahren zur bestimmung von impulsantworten von einem medium in bezug auf die übertragung von wellen zwischen verschiedenen punkten
AT04742298T ATE459014T1 (de) 2003-03-24 2004-03-19 Verfahren zur bestimmung von impulsantworten von einem medium in bezug auf die übertragung von wellen zwischen verschiedenen punkten
EP04742298A EP1620749B1 (fr) 2003-03-24 2004-03-19 Procédé pour déterminer des réponses impulsionnelles d'un milieu vis-à-vis de la transmission d'ondes entre différents points

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR03/03545 2003-03-24
FR0303545A FR2853075B1 (fr) 2003-03-24 2003-03-24 Procede pour determiner des reponses impulsionnelles d'un milieu visa-vis de la transmission d'ondes entre differents points

Publications (2)

Publication Number Publication Date
WO2004086557A2 true WO2004086557A2 (fr) 2004-10-07
WO2004086557A3 WO2004086557A3 (fr) 2004-10-28

Family

ID=32947109

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2004/000684 Ceased WO2004086557A2 (fr) 2003-03-24 2004-03-19 Procede pour determiner des reponses impulsionnelles d’un milieu vis-à-vis de la transmission d’ondes entre differents points

Country Status (7)

Country Link
US (1) US7656748B2 (fr)
EP (1) EP1620749B1 (fr)
AT (1) ATE459014T1 (fr)
DE (1) DE602004025690D1 (fr)
ES (1) ES2341005T3 (fr)
FR (1) FR2853075B1 (fr)
WO (1) WO2004086557A2 (fr)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2932339A1 (fr) * 2008-06-09 2009-12-11 Centre Nat Rech Scient Procede et dispositif de sondage par propagation d'ondes
WO2014009424A1 (fr) * 2012-07-12 2014-01-16 Universite Pierre Et Marie Curie (Paris 6) Procédé et système d'inspection d'une zone
WO2016181054A1 (fr) 2015-05-12 2016-11-17 Electricite De France Procédé et dispositif de sondage par propagation d'ondes

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2903827B1 (fr) * 2006-07-11 2009-01-23 Centre Nat Rech Scient Procede et dispositif de transmission d'ondes.
FR2916054B1 (fr) * 2007-05-07 2010-09-10 Altran Technologies Altran Procede et dispositif de detection acoustique en milieu liquide.
GB0721694D0 (en) * 2007-11-05 2007-12-12 Univ Bristol Methods and apparatus for measuring the contents of a search volume

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1423366A (en) * 1972-07-21 1976-02-04 Seiscom Ltd Broad line seismic profiling using simultaneously radiating sources
US4167879A (en) * 1978-06-06 1979-09-18 Panametrics, Inc. Method and apparatus for examining a solid
FR2492109B1 (fr) * 1980-10-10 1985-07-05 Thomson Csf Systeme d'imagerie a emissions multiples et simultanees
US4715020A (en) * 1986-10-29 1987-12-22 Western Atlas International, Inc. Simultaneous performance of multiple seismic vibratory surveys
US4823326A (en) * 1986-07-21 1989-04-18 The Standard Oil Company Seismic data acquisition technique having superposed signals
US6005827A (en) * 1995-03-02 1999-12-21 Acuson Corporation Ultrasonic harmonic imaging system and method
US5608690A (en) * 1995-03-02 1997-03-04 Acuson Corporation Transmit beamformer with frequency dependent focus
US5721710A (en) * 1995-09-29 1998-02-24 Atlantic Richfield Company High fidelity vibratory source seismic method with source separation
US6193659B1 (en) * 1997-07-15 2001-02-27 Acuson Corporation Medical ultrasonic diagnostic imaging method and apparatus
US6138512A (en) * 1997-07-30 2000-10-31 Iowa State University Research Foundation, Inc. Method and apparatus for determining source location of energy carried in the form of propagating waves through a conducting medium
CA2277119C (fr) * 1997-11-14 2003-06-03 Western Atlas International, Inc. Acquisition et traitement de donnees sismiques par utilisation d'une distorsion non lineaire dans un signal de force de sol
FR2815717B1 (fr) * 2000-10-20 2003-01-10 Centre Nat Rech Scient Procede et dispositif non invasif de focalisation d'ondes acoustiques
AUPS051102A0 (en) * 2002-02-13 2002-03-07 Curtin University Of Technology Method of producing continuous, orthogonal signals and method of their use for examining and for detecting changes in a body

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2932339A1 (fr) * 2008-06-09 2009-12-11 Centre Nat Rech Scient Procede et dispositif de sondage par propagation d'ondes
WO2010001027A1 (fr) * 2008-06-09 2010-01-07 Centre National De La Recherche Scientifique - Cnrs - Procede et dispositif de sondage par propagation d'ondes.
US8814795B2 (en) 2008-06-09 2014-08-26 Centre National de la Recherche Scientifique —CNRS— Sounding method and device using wave propagation
WO2014009424A1 (fr) * 2012-07-12 2014-01-16 Universite Pierre Et Marie Curie (Paris 6) Procédé et système d'inspection d'une zone
FR2993368A1 (fr) * 2012-07-12 2014-01-17 Univ Paris Curie Procede et systeme d'inspection d'une zone.
US9507036B2 (en) 2012-07-12 2016-11-29 Universite Pierre Et Marie Curie (Paris 6) Method and system for inspecting a zone
WO2016181054A1 (fr) 2015-05-12 2016-11-17 Electricite De France Procédé et dispositif de sondage par propagation d'ondes
US10267914B2 (en) 2015-05-12 2019-04-23 Electricite De France Method and device for probing by wave propagation

Also Published As

Publication number Publication date
FR2853075B1 (fr) 2005-06-10
US7656748B2 (en) 2010-02-02
EP1620749A2 (fr) 2006-02-01
WO2004086557A3 (fr) 2004-10-28
US20060187754A1 (en) 2006-08-24
FR2853075A1 (fr) 2004-10-01
ES2341005T3 (es) 2010-06-14
DE602004025690D1 (de) 2010-04-08
ATE459014T1 (de) 2010-03-15
EP1620749B1 (fr) 2010-02-24

Similar Documents

Publication Publication Date Title
EP2294452B1 (fr) Procede et dispositif de sondage par propagation d'ondes
CA2222890C (fr) Methode de transmission de donnees sismiques compressees
CA1325840C (fr) Procede d'acquisition de donnees sismiques et dispositif de mise en oeuvre du procede
FR2904168A1 (fr) Systeme d'estimation de la qualite de reception d'une transmission numerique.
EP0732803A1 (fr) Procédé et dispositif de démodulation par échantillonnage
EP0676649A1 (fr) Procédé et dispositif radar de mesure de distance
FR2950701A1 (fr) Procede de separation de trains d'impulsions radar imbriques
EP3847450A1 (fr) Procede de detection et de caracterisation par ultrasons de defauts dans un materiau heterogene
EP1620749B1 (fr) Procédé pour déterminer des réponses impulsionnelles d'un milieu vis-à-vis de la transmission d'ondes entre différents points
FR2784182A1 (fr) Debitmetre doppler a frequences multiples
Socheleau et al. Characterisation of time-varying underwater acoustic communication channel with application to channel capacity
FR2890177A1 (fr) Procede et systeme d'acquisition de donnees sismiques
EP0594474A1 (fr) Procédé et dispositif de trajectographie passive d'objets sous-marins
EP2168336B1 (fr) Méthode de détection de présence de signaux étalés spectralement
EP0197582B1 (fr) Procédé et appareil d'exploration de milieux par échographie ultrasonore
EP3671250A1 (fr) Interféromètre numérique à sous-échantillonnage
EP2851703A1 (fr) Procédé réalisant conjointement la synchronisation, l'identification, la mesure, l'estimation du filtre de propagation et la localisation d'émetteurs utiles et interferants
EP0580695B1 (fr) Procede de compensation du mouvement de l'antenne pour un sonar
FR3033067A1 (fr) Procede de determination de la bathymetrie a partir d'images satellitaires optiques quasi synchrones
EP3470871A1 (fr) Procédé de détection de signaux radar
CA2134856A1 (fr) Procede de determination du rang d'ambiguite en distance d'echos radar
EP3674717B1 (fr) Dispositif pour la determination du spectre d'un signal sur une large bande de frequences par echantillonnages multiples et procede associe
FR3131641A3 (fr) Procédé de génération d’un signal de test pour capteur RF
FR3165538A1 (fr) Procédé de traitement d’un signal reçu par un récepteur, le signal reçu étant issu de deux signaux à modulation de phase continue
FR2916054A1 (fr) Procede et dispositif de detection acoustique en milieu liquide.

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): BW GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2006187754

Country of ref document: US

Ref document number: 10550429

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2004742298

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2004742298

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 10550429

Country of ref document: US