FR3032799B1 - Procede pour la determination de proprietes d'un milieu et dispositif pour la determination de proprietes d'un milieu - Google Patents
Procede pour la determination de proprietes d'un milieu et dispositif pour la determination de proprietes d'un milieu Download PDFInfo
- Publication number
- FR3032799B1 FR3032799B1 FR1650833A FR1650833A FR3032799B1 FR 3032799 B1 FR3032799 B1 FR 3032799B1 FR 1650833 A FR1650833 A FR 1650833A FR 1650833 A FR1650833 A FR 1650833A FR 3032799 B1 FR3032799 B1 FR 3032799B1
- Authority
- FR
- France
- Prior art keywords
- propagation
- medium
- signals
- reception
- acoustic wave
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000005540 biological transmission Effects 0.000 claims abstract description 109
- 239000000126 substance Substances 0.000 claims abstract description 22
- 230000010363 phase shift Effects 0.000 claims abstract description 19
- 230000000644 propagated effect Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims description 40
- 238000004364 calculation method Methods 0.000 claims description 23
- 238000010897 surface acoustic wave method Methods 0.000 claims description 22
- 230000001902 propagating effect Effects 0.000 claims description 16
- 230000001934 delay Effects 0.000 claims description 5
- 230000003534 oscillatory effect Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 2
- 239000000758 substrate Substances 0.000 description 39
- 238000005259 measurement Methods 0.000 description 13
- 230000003993 interaction Effects 0.000 description 11
- 238000005070 sampling Methods 0.000 description 9
- 238000011144 upstream manufacturing Methods 0.000 description 8
- 235000019687 Lamb Nutrition 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000002787 reinforcement Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004026 adhesive bonding Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 235000011837 pasties Nutrition 0.000 description 1
- 239000007779 soft material Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
- G01N29/024—Analysing fluids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/022—Liquids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02818—Density, viscosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0423—Surface waves, e.g. Rayleigh waves, Love waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0426—Bulk waves, e.g. quartz crystal microbalance, torsional waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0427—Flexural waves, plate waves, e.g. Lamb waves, tuning fork, cantilever
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/105—Number of transducers two or more emitters, two or more receivers
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
Abstract
La présente invention concerne en particulier un procédé pour la détermination de propriétés physiques et/ou chimiques d'un milieu (M) sur la base d'au moins une première et une deuxième onde acoustique, lesquelles se sont propagées au moins partiellement à travers le milieu (M) à partir d'au moins un émetteur (SE1, SE2; SE3, SE4) jusqu'à au moins un récepteur (SE2, SE1; SE4, SE3), dans lequel on détermine, à partir de signaux de réception (yr 0, yr 1; ys 0, ys 1) générés au niveau d'au moins deux récepteurs, une différence de temps de propagation (Δt) des ondes acoustiques et/ou un temps de propagation absolu (t0) d'une onde acoustique, et au moyen d'une différence de temps de propagation (Δt) déterminée et/ou d'un temps de propagation absolu (t0) déterminé, des propriétés physiques et/ou chimiques du milieu (M) étant déterminées, telles par exemple une vitesse moyenne d'écoulement (vm). Selon l'invention, - on utilise au moins deux signaux d'émission (x0, x1) différents, sensiblement de même fréquence de base (f), lesquels présentent un déphasage prescrit entre eux et grâce auxquels on génère respectivement une première et une deuxième onde acoustique de telle sorte que, pour chaque direction de propagation, deux signaux de réception (yr 0, yr 1; ys 0, ys 1) sont générés, lesquels relèvent de signaux d'émission (x0, x1) différents, et - une différence de temps de propagation (Δt) et/ou un temps de propagation absolu (t0) étant déterminés grâce aux quatre signaux de réception (yr 0, yr 1; ys 0, ys 1). On met en outre à disposition un dispositif pour la détermination de propriétés physiques et/chimiques d'un milieu (M), lequel est en particulier réalisé et prévu pour la mise en œuvre d'un procédé selon l'invention.
Description
Procédé pour la détermination de propriétés d’un milieu et dispositif pour la détermination de propriétés d’un milieu
Description
La présente invention concerne un procédé pour la détermination de propriétés physiques et /ou chimiques d'un milieu à l’aide d’au moins une première et une deuxième onde acoustique, lesquelles se sont respectivement propagées au moins partiellement à travers le milieu à partir d'un émetteur vers un récepteur selon la présente invention, ainsi qu'un dispositif pour la détermination de propriétés physiques et/ou chimiques d’un milieu selon la présente invention.
Le milieu dont les propriétés physiques et/ou chimiques doivent être déterminées par un procédé du genre est de préférence un liquide ou un matériau tendre, en particulier un milieu hautement visqueux, semblable à une pâte ou pâteux. Concernant les ondes acoustiques utilisées pour la détermination des propriétés, il s'agit par exemple d’ondes ultrasonores qui sont générées par un émetteur correspondant grâce à un signal d’émission.
Dans un procédé du genre, on génère généralement deux ondes acoustiques grâce à un signal d’émission, lesquelles se propagent au moins partiellement à travers le milieu le long de directions de propagation identiques ou différentes avant d’être respectivement réceptionnées au niveau d’un récepteur situé dans la direction de propagation respective. Des ondes acoustiques sont par exemple générées dans un milieu s’écoulant, d’une part dans une première direction de propagation en direction d’écoulement du milieu, et d’autre part dans une deuxième direction de propagation à l'encontre de la direction d’écoulement du milieu. A partir des signaux de réception générés au niveau de récepteurs respectifs, il est alors possible de déterminer une différence de temps de propagation et d’en déduire par exemple la vitesse d’écoulement (moyenne) du milieu. Si l’on détermine encore, en variante ou additionnellement, des temps de propagation absolus d'une onde acoustique d’un émetteur vers un récepteur à l’aide des signaux de réception, il est possible d'obtenir d’autres conclusions sur des propriétés physiques et/ou chimiques du milieu, telles par exemple la densité, la température ou la composition.
Le document WO 2008/034878 A2 a fait connaître un dispositif du genre, dans lequel des ondes acoustiques de surface sont générées, lesquelles couplent des ondes sonores volumiques dans le milieu respectif dans un guide d’ondes. Grâce au découplage répété d’ondes de surface à ceux des endroits où l’onde sonore volumique touche une paroi délimitant le milieu, des ondes acoustiques de surface sont de nouveau réceptionnées au niveau d’un récepteur, dont les temps de propagation et différences de temps de propagation sont caractéristiques pour le milieu ainsi que les propriétés physiques et/ou chimiques de celui-ci.
Dans un dispositif décrit dans le document WO 2008/034878 A2, ainsi que le procédé mis en œuvre avec, le traitement des signaux de réception générés au niveau des récepteurs respectifs pour une onde acoustique réceptionnée - ici, une onde acoustique de surface - revêt ainsi une importance décisive. C’est ainsi que la détermination d’une différence de temps de propagation ou d’un temps de propagation absolu à partir des signaux de réception générés au niveau du récepteur n’est en aucun cas triviale, et s’accompagne le cas échéant d’efforts importants de calcul. Selon les informations à extraire des signaux de réception, on utilise des procédés les plus variés pour le traitement des signaux. On connaît par exemple le fait d’utiliser des signaux d’émission modulés afin de, grâce aux signaux de réception obtenus, pouvoir en déduire de manière plus fiable les propriétés du milieu. A cet effet, on utilise par exemple une modulation d’amplitude en quadrature, en bref : modulation IQ, avec laquelle, par rapport à des signaux d’émission non modulés, on peut aussi régulièrement parvenir à une amélioration nette de la résolution pouvant être obtenue. D’autres procédés connus reposent généralement sur des signaux de réception échantillonnés et soumis à une conversion analogique vers numérique, des fonctions de corrélation, tout comme également des transformées de Hilbert et Wavelett, étant formées.
Dans le document DE 102 06 134 A1, on révèle par ailleurs un procédé pour la détermination d’une vitesse d’écoulement d’un milieu au moyen d’ondes ultrasonores, lequel doit êtrS-fésistant par rapport à des perturbations, telles un bruit, lors du traitement du signal. Pour ce faire, des signaux de réception analogiques générés sont numérisés, une corrélation croisée est effectuée, et une conversion de Hilbert est réalisée afin de calculer un rapport de phases entre des ondes acoustiques réceptionnées et calculer une différence de temps à partir de cela. Un dispositif adapté pour cela afin de traiter les signaux de réception (dispositif de traitement de signaux) est par conséquent comparativement complexe et s’accompagne de coûts non négligeables.
Sur ces bases, le but de la présente invention consiste à mettre à disposition un procédé amélioré pour la détermination de propriétés physiques et/ou chimiques d'un milieu au moyen d’ondes acoustiques, lequel est amélioré concernant l’effort de calcul lors du traitement de signaux des signaux de réception, et lequel permet en particulier une détermination particulièrement fiable et robuste ainsi que moins onéreuse d’une vitesse d’écoulement d’un milieu au moyen d’ondes acoustiques. Un dispositif pour la mise en œuvre de ce procédé doit en outre être mis à disposition.
Ce but est atteint avec le procédé selon la présente invention, tout comme avec le dispositif selon l’invention.
Les procédés présentés ci-dessous interviennent dans ce cas pour une détermination d’une différence de temps de propagation entre des ondes acoustiques réceptionnées s’étant propagées le long de directions de propagation identiques et/ou différentes et une détermination d’un temps de propagation absolu d'une onde acoustique d'au moins un émetteur vers au moins un récepteur, dans lequel une onde acoustique s’est toujours au moins partiellement propagée à travers le milieu vers un récepteur. Dans ce cas, les deux procédés selon l’invention peuvent bien entendu également être combinés entre eux sans aucun problème, de telle sorte qu'après la détermination d’un temps de propagation absolu par exemple et d'une distance connue entre l’émetteur et le récepteur, il est possible d’en déduire la vitesse du son des ondes acoustiques dans le milieu et qu’un facteur de proportionnalité est ainsi mis à disposition, lequel peut être utilisé, ensemble avec une différence de temps de propagation, pour le calcul d’une vitesse (moyenne) d’écoulement du milieu en écoulement.
Les deux procédés partent respectivement du principe - qu'au moins une première et une deuxième onde acoustique se sont propagées au moins partiellement d’un émetteur vers un récepteur à travers le milieu, et - que les ondes acoustiques sont respectivement générées par un signal d’émission, - que la première onde acoustique se propage à travers le milieu le long d’une première direction.de propagation et que la deuxième onde acoustique se propage à travers le milieu le long d’une deuxième direction de propagation de préférence différente, par exemple opposée, à la première direction de propagation, et - que, pour une première onde acoustique réceptionnée après sa propagation à travers le milieu le long de la première direction de propagation, un premier signal de réception est généré, et que pour une deuxième onde acoustique réceptionnée après sa propagation à travers le milieu le long de la deuxième direction de propagation, un deuxième signal de réception est généré.
On prévoit selon l’invention que, pour la détermination de la différence de temps de propagation des ondes acoustiques et/ou d’un temps de propagation absolu d’une onde acoustique d’un émetteur vers un récepteur, on utilise au total au moins quatre signaux de réception - deux premiers et deux deuxièmes signaux de réception, lesquels relèvent d’au moins deux signaux d’émission différents. Ces signaux d’émission présentent une fréquence de base sensiblement identique et un déphasage prescrit entre eux. C’est ainsi qu’en utilisant le signal d’émission identique - également le même dans une variante - on génère ainsi respectivement une première et une deuxième onde acoustique dans chaque direction de propagation, de sorte que pour chaque direction de propagation, il y a deux signaux de réception. C’est ainsi qu’un signal d’émission est par exemple utilisé pour générer simultanément ou successivement deux ondes acoustiques se propageant dans le milieu le long de la première et de la deuxième direction de propagation, dans lequel, pour une onde acoustique se propageant le long de la première direction de propagation - par exemple en aval d’une direction d’écoulement d’un milieu s’écoulant - on génère un premier signal de réception, et pour une onde acoustique se propageant le long de la deuxième direction de propagation - vers l’amont - un deuxième signal de réception est généré. Un autre signal d’émission à décalage de phase avec une fréquence de base identique est alors utilisé selon l’invention pour générer deux autres ondes acoustiques se propageant dans le milieu le long de la première et de la deuxième direction de propagation, dans lequel, dans ce cas également, pour une onde acoustique se propageant le long de la première direction de propagation - vers l’aval - un autre premier signal de réception qui en diffère est généré, et pour une onde acoustique se propageant le long de la deuxième direction de propagation - vers l’amont - un autre deuxième signal de réception qui en diffère est généré.
Dans ce cas, pour la génération des deux ondes acoustiques, un signal d’émission peut être émis par deux émetteurs différents espacés spatialement l’un de l’autre. En variante, il est en principe également envisageable qu’un seul émetteur situé entre deux récepteurs déclenche tout d’abord, avec un signal d’émission ou avec deux signaux d’émission identiques, deux ondes acoustiques en direction des deux récepteurs et déclenche ensuite, avec un signal d’émission à fréquence de base identique avec un décalage de phase par rapport au signal d’émission précédemment envoyé ou avec deux signaux d’émission à décalage de phase par rapport au signal d’émission précédemment envoyé, deux autres ondes acoustiques en direction des récepteurs. La seule chose qui importe est de générer une paire de signaux de réception dans chaque direction de propagation, lesquels relèvent de deux signaux d’émission à décalage de phase l’un par rapport à l’autre avec une fréquence de base identique.
On entend ici par signaux d’émission à fréquence de base sensiblement identique, des signaux d’émission qui divergent l’un de l’autre concernant leur fréquence de base d’un maximum de 1/1000 par rapport à une fréquence centrale. L’idée de base de la présente invention consiste dans ce cas, à partir des quatre signaux de réception dans leur ensemble, lesquels relèvent de deux signaux d’émission différents à fréquence de base identique et déphasage prescrit entre eux, grâce à une utilisation minimale en matériel lors du traitement des signaux, de manière rapide et avec une incertitude de mesure comparativement faible, à déterminer une différence de temps de propagation et/ou un temps de propagation absolu, et à déterminer après cela des propriétés physiques et/ou chimiques du milieu. C’est ainsi qu’avec la présente invention, grâce à l’utilisation selon l’invention de deux signaux d’émission différents, aucune génération onéreuse de signaux de différence ou de synchronisation des signaux d’émission et de réception n’est nécessaire. Il n’est pas non plus nécessaire de prévoir un traitement numérique onéreux des signaux, étant donné qu’en particulier pour une numérisation éventuelle de signaux de réception ou d’une numérisation éventuelle de signaux qui en résultent, il est possible de renoncer à la transformation de Fourier qui est sinon régulièrement avantageuse.
De manière avantageuse, les quatre signaux de réception analogiques ou numérisés sont ici directement pris en compte dans un calcul, à savoir, en particulier en tant que paramètres dans une équation de calcul mise en œuvre par un dispositif de traitement de signaux, afin de calculer une différence de temps de propagation et/ou un temps de propagation absolu. Par conséquent, on ne réalise par exemple pas une simple moyenne de différences de temps de propagation déjà calculées avec une paire de signaux de réception. Dans une variante de réalisation respectivement préférée, chacun des signaux de réception entre plutôt directement dans une règle de calcul pour une différence de temps de propagation et/ou une règle de calcul pour un temps de propagation absolu, de sorte qu’un calcul peut seulement avoir lieu en cas de présence de l’ensemble des quatre signaux de réception.
Pour le calcul d’une différence de temps de propagation, on prévoit par exemple que des signaux de réception pour des ondes acoustiques réceptionnées, lesquelles ont été générées par des signaux d'émission différents, soient multipliés entre eux, de préférence de manière croisée.
Selon la configuration d’un dispositif de traitement de signaux utilisé, il est par exemple possible de prévoir ici une multiplication de signaux de réception analogiques, avec un calcul qui s’en suit d’une différence de temps de propagation à partir de signaux numérisés, c’est-à-dire résultant d’un échantillonnage (anglais : sampling). Une variante en guise d’alternative pour un calcul d’une différence de temps de propagation consisterait à prévoir, par rapport à cela, d’abord une multiplication de signaux de réception déjà numérisés.
Indépendamment de cela, pour la détermination d’une différence de temps de propagation At à partir de signaux de réception γθ(1), yir(t); y®(t) et yis(t), on fait appel à une formule de la forme qui suit :
[Equation 1]
Dans l’équation 1, les variables munies d’un indice « 0 » représentent ceux des signaux qui relèvent d’un premier signal d’émission x0, alors que les signaux désignés par un indice « 1 » relèvent d’un deuxième signal d’émission xt qui en est décalé en phase avec une fréquence de base identique. Un exposant « r » désigne en outre ceux des signaux qui relèvent d’une première onde acoustique réceptionnée dans une première direction de propagation, par exemple en aval d’un milieu en écoulement. Les signaux désignés par un exposant « s » relèvent quant à eux d’une deuxième onde acoustique réceptionnée, de préférence dans une deuxième direction de propagation opposée, par exemple en amont.
Avec les valeurs pour un numérateur Z(t) et un dénominateur N(t), on détermine ensuite directement une différence de temps de propagation At dans une variante de réalisation préférée, d’après la formule suivante :
[Equation 2]
Ici, f représente une fréquence de base pour les deux signaux d’émission.
En cas de fréquences de base f0 et fi divergeant peu des signaux d’émission différents x0, Xi, cette fréquence de base f peut également être déterminée par exemple à partir de la moyenne arithmétique des deux fréquences de base f0 et h :
[Equation 2.1] C’est ainsi que les signaux de réception individuels y'(t), y^t); y®(t), y®(t) n’ont pas besoin d’être convertis grâce à un traitement ultérieur fastidieux en un signal analytique, mais on prévoit déjà du côté d’émission, ou en tant que référence, au moins deux signaux différents x0 (t) et Xi (t) avec une même fréquence de base f et un déphasage prescrit, de préférence de 90°.
Les signaux d’émission utilisés présentent alors de préférence la forme suivante :
[Equations 3]
Dans lesquelles φ désigne un angle de départ pouvant être librement sélectionné et w(t) représente une fonction de fenêtre, laquelle affiche les signaux porteurs à décalage de phase de préférence de 90° des signaux d’émission x0 et Xi, et les masque de nouveau après une durée déterminée. Une telle fonction de fenêtre peut par exemple également comprendre une séquence de modulation codée. Une plage particulièrement grande concernant la durée totale d’un signal d’émission réalisé en tant que rafale d’émission pour l’évaluation de la différence de phase ou de temps de propagation est obtenue ici lors du choix d’une fenêtre de Blackman-Nuttal.
Un signal d’émission x0, χΊ comprend ainsi de préférence respectivement une impulsion oscillatoire ou rafale impulsionnelle, dans lequel les enveloppantes pour les impulsions oscillatoires individuelles des signaux d’émission à décalage de phase entre eux sont sensiblement identiques entre elles. C'est ainsi que les signaux de réception obtenus à partir des signaux d’émission pour les ondes acoustiques se propageant au moins partiellement à travers le milieu peuvent être recalculés directement de manière particulièrement simple sur une base temporelle commune - également décalée de manière définie - en signaux d’amplitude et de phase ainsi que des temps de propagation et différences de temps de propagation, sans références supplémentaires comme cela est par exemple nécessaire pour une modulation IQ et sans divisions par des valeurs proches de 0, comme on le montrera encore plus en détails dans ce qui suit.
La détermination d’une différence de phase ou de temps de propagation pour les différentes premières et deuxièmes ondes acoustiques repose, lors d’une détermination d’une vitesse (moyenne) d’écoulement vm d’un milieu en écoulement, sur les différents temps de propagation dans le milieu en fonction de la direction.
On se sert ici du fait que, dans des directions de propagation différentes, par exemple vers l’aval et par conséquent en direction d’écoulement d’une part, et vers l’amont et par conséquent à l’encontre de la direction d’écoulement d’autre part, différentes vitesses du son apparaissent, lesquelles dépendent de la vitesse d’écoulement. C’est ainsi qu’une vitesse du son en direction d’écoulement est de cr = c0 + vm, et une vitesse du son à l’encontre de la direction d’écoulement est de cs = c0 - vm, où c0 est respectivement la vitesse du son dans le milieu au repos.
Si l’émetteur et le récepteur forment alors une distance l0 de manière absolue et sont éloignés l’un de l’autre d’une distance I dans la projection sur la direction d’écoulement du milieu, le long de cette distance, on obtient les temps de propagation du son correspondants pour les ondes acoustiques individuelles
[Equation 4.1] en direction d’écoulement ou
[Equation 4.2] à l’encontre de la direction d’écoulement.
Pour une différence de temps de propagation At des première et deuxième ondes acoustiques, on obtient par conséquent
[Equation 5]
C’est ainsi que la différence de temps de propagation At est quasi proportionnelle à la vitesse d’écoulement moyenne vm.
Le facteur correspondant de proportionnalité l/(c0)2 correspondant peut être prescrit, par exemple grâce à la détermination (unique) d’un temps de propagation absolu d’une onde acoustique d’un émetteur vers un récepteur. Pour ce faire, on peut utiliser un procédé selon l’invention pour la détermination d’un temps de propagation absolu. En cas de distance connue entre l'émetteur et le récepteur d’un agencement déterminé d’émetteur-récepteur et d’une vitesse du son connue dans un milieu déterminé, dont la vitesse d’écoulement doit être mesurée et/ou surveillée, le facteur de proportionnalité peut également être réglé de manière fixe dans le dispositif de traitement de signaux.
Au moyen des deux signaux d’émission x0, Xi à décalage de phase, on génère ainsi les signaux de réception
[Equations 6]
Une émission simultanée de l’ensemble des quatre signaux d’émission, c’est-à-dire une émission d’un signal d’émission Xo en direction d’écoulement et à l’encontre de la direction d’écoulement ainsi qu’une émission d’un signal d’émission Xt en direction d'écoulement et à l’encontre de la direction d’écoulement, est considérée comme avantageuse ici afin de multiplier par exemple entre eux de manière croisée les signaux de réception analogiques γθ (t), yj (t); y® (t), y® (t) conformément à l’équation 1, et ensuite procéder à un échantillonnage avec un faible taux d’échantillonnage pour la numérisation des valeurs Z(t) et N(t) obtenues de ce fait.
En variante, il est également possible d’émettre les signaux d’émission individuels les uns après les autres et de générer ainsi des ondes acoustiques successivement. Les signaux d'émission individuels sont donc générés séquentiellement par le ou les émetteur(s) respectif(s), et les signaux de réception sont également générés
séquentiellement au niveau du récepteur ou des récepteurs. Dans ce cas, on travaille de préférence avec un taux d’échantillonnage élevé.
Avec quatre signaux de réception pouvant être différenciés les uns des autres, lesquels relèvent de deux signaux d’émission à décalage de phase entre eux, comme déjà mentionné en préambule, non seulement une détermination d’une différence de temps de propagation, mais aussi la détermination d’un temps de propagation absolu sont possibles. Dans ce cas, il ne serait par ailleurs pas obligatoire que les ondes acoustiques se propagent à travers le milieu le long d’au moins deux première et deuxième directions de propagation différentes l’une de l’autre. Egalement en cas de directions de propagation identiques, sur la base de la solution selon l’invention, il serait possible de déterminer un temps de propagation absolu au moyen de quatre signaux de réception différenciables les uns des autres.
Dans ce cas, on travaille de préférence avec une valeur de départ estimée et une valeur de différence calculée à partir des quatre signaux de réception. On additionne en l’occurrence une valeur de différence calculée avec précision à partir des quatre signaux de réception à la valeur de départ estimée (grossièrement). Ce qui importe dans ce cas, c’est la définition d’un critère solide de temps de propagation, par exemple la position du maximum de corrélation croisée dans la plage de temps de décalage. On s’efforce globalement à ce que l’incertitude de mesure de la valeur de départ estimée pour le temps de propagation soit inférieure à la moitié d’une durée de période 1/f du signal porteur utilisé pour les signaux d’émission (par exemple xoet Xt conformément aux équations 3 ci-dessus).
Dans ce cas, on présuppose généralement un calcul général pour un temps de propagation absolu t0 à partir d’une valeur de départ t0 et d’une valeur de différence At’ de manière à ce que s’applique :
[Equation 7]
Ce faisant, pour le calcul de la valeur de différence At’, il s’est avéré par exemple avantageux d’utiliser au moins l’un des termes suivants :
[Equations 8]
Ici, f est une valeur estimée pour la fréquence centrale des signaux de réception γθ, y,; y®, y® et tk une constante de temps prescrite avec tk s 0.
La fonction « antan2 » représente dans ce cas une fonction d’arc tangente bien connue en programmation et à déterminer avec un faible effort de calcul, laquelle peut par exemple être définie via les propriétés suivantes :
[Equations 9]
Une valeur de différence At’, et par conséquent un temps de propagation absolu t0 peuvent alors être définis par calcul, par exemple conformément aux équations 7 et 8 ci-dessus, à partir du rapport de formule suivant, de manière comparativement simple sur le plan du calcul grâce à un dispositif de traitement de signaux sans transformée de Fourier:
[Equation 10]
Une formule du genre qui suit constitue une possibilité en variante :
[Equation 11]
Un calcul selon la méthode suivante constitue une autre variante :
[Equations 12]
Une règle de calcul de la forme suivante s’est par ailleurs également avérée avantageuse:
[Equation 13.1] avec
[Equations 13.2]
Une erreur lors de l’estimation de la fréquence centrale f peut être corrigée en ce que les temps de propagation des groupes de signaux considérés sont approximés Λ par une droite. La valeur de fonction de ces droites avec la valeur de départ t0 estimée est alors la meilleure approximation pour le temps de propagation absolu t0, et est alors largement indépendante de la fréquence centrale f mise en équation pour le calcul.
Un ajustage (ou étalonnage du zéro) peut par exemple avoir lieu par réglage de l’angle de départ φ de signaux d’émission x0 et Xi ou l’addition de la constante de temps tk au terme de compensation chrono proportionnel. Λ
Une évaluation pour la valeur de départ t0 lors du calcul du temps de propagation absolu t0 peut également avoir lieu à l’aide des quatre signaux de réception γθ(ί), y' (t)i yo(t)> yî(t) générés lors de la réception des ondes acoustiques. Dans ce cas, une valeur de départ peut par exemple être prescrite pour le calcul du temps de propagation au moyen d’au moins une enveloppe de signal yhuei qui est donnée par la formule suivante:
[Equation 14.1] ou plus précisément
[Equations 14.2]
Par conséquent, grâce par exemple à l’enveloppe de signal, il est tout d’abord A. possible d’estimer une valeur de départ t0, et ensuite, via une valeur de différence
At’ (calculée dans ce qui suit), de déterminer un temps de propagation absolu t0. Avec le temps de propagation absolu t0, il est par exemple possible de calculer le facteur de proportionnalité l/(c0)2, et après une détermination d’une différence de
temps de propagation At, de fixer une vitesse d’écoulement du milieu en écoulement.
Dans une variante de réalisation, on peut en outre prévoir qu'en cas d’une détermination de plusieurs valeurs pour des différences de temps de propagation et/ou pour des temps de propagation absolus d’une pluralité d’ondes acoustiques, le calcul d'une moyenne et/ou une intégration des différences de temps de propagation déterminées et/ou des temps de propagation soient effectués pour minimiser d’éventuelles incertitudes de mesure. D’une manière connue par ailleurs, il est bien entendu possible, au moyen d’une différence de temps de propagation déterminée et/ou au moyen d’un temps de propagation absolu déterminé des première et deuxième ondes acoustiques générées, de déterminer de la même manière une vitesse du son, une concentration, une consistance et/ou une température du milieu respectif, telle un niveau de remplissage du milieu dans un espace intérieur, une épaisseur d’une paroi jouxtant le milieu et/ou une distance de deux parties de paroi jouxtant le milieu, comme cela est également expliqué par exemple dans le document WO 2008/034878 A2 en rapport avec des ondes acoustiques de surface.
Un autre aspect de la présente invention est la mise à disposition d’un dispositif tel que défini ci-dessous, lequel est en particulier réalisé et prévu pour la mise en oeuvre du procédé selon l’invention.
Un tel dispositif pour la détermination de propriétés physiques et/ou chimiques d’un milieu présente généralement ce qui suit : - un guide d’ondes acoustique, lequel présente un espace intérieur à remplir avec le milieu, - au moins deux émetteurs, lesquels sont respectivement réalisés et prévus pour la génération d’ondes acoustiques dans le guide d'ondes grâce à un signal d'émission de l’émetteur de telle sorte qu’une onde acoustique générée se propage le long du guide d’ondes et au moins partiellement à travers le milieu, - au moins deux récepteurs, lesquels sont réalisés et prévus de manière espacée spatialement l’un par rapport à l’autre le long du guide d’ondes et respectivement pour la génération d’un premier ou deuxième signal de réception pour des première ou deuxième ondes acoustiques atteignant le récepteur respectif, dans lequel une première onde acoustique se propage respectivement le long d’une première direction de propagation, et un premier signal de réception étant généré à cet effet au niveau d’un récepteur, et une deuxième onde acoustique se propageant respectivement le long d’une deuxième direction de propagation identique à la première direction de propagation ou différente de celle-ci, et un deuxième signal de réception étant généré à cet effet au niveau d’un récepteur, et - un dispositif de traitement de signaux, au moyen duquel on détermine, à partir des signaux de réception, • une différence de temps de propagation entre des premières et deuxièmes ondes acoustiques réceptionnées lorsque la deuxième direction de propagation est différente de la première direction de propagation, et/ou, • un temps de propagation absolu d’une onde acoustique d’un émetteur vers un récepteur lorsque la deuxième direction de propagation est différente de la première direction de propagation ou est identique à celle-ci, et lequel est réalisé et prévu pour déterminer, avec une différence de temps de propagation déterminée et/ou un temps de propagation absolu déterminé, des propriétés physiques et/ou chimiques du milieu.
Par analogie au procédé selon l’invention, les émetteurs du dispositif sont dans ce cas également réalisés et prévus pour générer, avec deux signaux d’émission différents à fréquence de base identique et déphasage prescrit entre eux, respectivement des premières et deuxièmes ondes acoustiques de telle sorte que dans la première direction de propagation, deux premiers signaux de réception pouvant être différenciés l’un de l’autre, et dans la deuxième direction de propagation, deux deuxièmes signaux de réception pouvant être différenciés l’un de l’autre sont générés par les récepteurs du dispositif. Le dispositif de traitement de signaux est alors réalisé et prévu pour déterminer une différence de temps de propagation et/ou un temps de propagation absolu grâce aux quatre signaux de réception générés, dans lequel on utilise de préférence à cet effet respectivement les équations susmentionnées et déposées dans le dispositif de traitement de signaux.
Une utilisation de transducteurs d’ultrasons en tant qu’émetteur ou récepteur s’est avérée particulièrement avantageuse, lesquels sont réalisés et prévus pour la génération d’ondes acoustiques de surface dans un guide d’ondes du dispositif. Par analogie avec un dispositif connu de par le document WO 2008/034878 A2, un tel dispositif se caractérise alors entre autres par un guide d’ondes acoustique, lequel comprend au moins deux éléments de guidage situés en vis-à-vis l’un de l’autre, lesquels délimitent un espace intérieur à remplir avec un milieu et lesquels forment, lors du remplissage de l’espace intérieur avec un milieu, respectivement avec une surface intérieure, une surface de délimitation avec le milieu. Au niveau de la surface intérieure respective, au moins une partie de l’onde acoustique de surface générée par un émetteur est convertie en ondes sonores volumiques du milieu, et au moins une partie des ondes sonores volumiques est à son tour convertie en ondes acoustiques de surface du guide d’ondes, lesquelles peuvent alors être réceptionnées au niveau du récepteur respectif. C’est justement avec un tel dispositif qu’il est particulièrement efficace de mettre en œuvre le procédé selon l’invention étant donné qu’ici, les quatre signaux de réception différenciables sont formés pour un périmètre de bilan comparativement petit et que d’excellents résultats sont obtenus pour le traitement des signaux avec un investissement nettement réduit lesquels, comme on a pu le confirmer grâce à des simulations et expériences correspondantes, présentent une incertitude de mesure extrêmement faible. C’est ainsi qu’avec un effort de calcul nettement réduit, on peut même obtenir une précision de mesure plus élevée le cas échéant lors de la détermination d’une vitesse d’écoulement d’un milieu qui s’écoule.
Dans une variante de réalisation possible d’un dispositif selon l’invention, le dispositif de traitement de signaux présente des moyens pour la multiplication de signaux de réception analogiques et des moyens pour le calcul qui s’en suit d'une différence de temps de propagation à partir de signaux numérisés, c’est-à-dire résultant d’un échantillonnage (anglais : « sampling »). Les moyens pour la multiplication des signaux de réception analogiques comprennent alors, par exemple des multiplicateurs, additionneurs et/ou inverseurs. Pour la numérisation, on utilise par exemple au moins un convertisseur analogique/numérique du dispositif.
En variante, le dispositif de traitement de signaux peut présenter des moyens pour la numérisation des signaux de réception et des moyens pour la multiplication qui s’en suit de signaux de réception numérisés pour le calcul d’une différence de temps de propagation. Dans ce cas, une variante de réalisation prévoit en guise d’exemple que les signaux de réception soient amenés vers un convertisseur analogique/numérique via un multiplexeur et que seulement les signaux numérisés provenant du convertisseur analogique/numérique soient soumis à un traitement ultérieur, par exemple conformément aux équations ci-dessus, afin de calculer à partir de cela une différence de temps de propagation At et/ou un temps de propagation absolu t0.
En principe, il est en outre possible de prévoir que la fréquence de base des deux signaux d’émission différents, à décalage de phase l’un par rapport à l’autre, tout comme leur déphasage, soit variable, et ainsi réglable. Conformément à l’équation 3 ci-dessus, on préfère que le déphasage des deux signaux d’émission (x0, x^ soit compris entre 45° et 135° et soit par exemple sensiblement de 90° ou de π/2.
Pour la génération des quatre ondes acoustiques par les deux signaux d’émission différents et la réception des quatre signaux de réception associés, il est bien entendu possible de prévoir des émetteurs ou récepteurs individuels dans le dispositif de mesure. On préfère cependant que le dispositif présente au moins une unité d’émission et de réception pouvant fonctionner au choix en tant qu’émetteur ou récepteur, de sorte que via une unité d’émission et de réception, dans un mode d’émission, il est possible de générer un signal d’émission pour une onde acoustique dans l’une des directions de propagation tout comme également de réceptionner via celle-ci, dans un mode de réception, une onde acoustique qui s’est propagée dans l’autre direction de propagation vers l’unité d’émission et de réception, et de pouvoir générer un signal de réception correspondant.
Dans un développement possible, concernant l’unité d’émission et de réception, il s’agit d’un transducteur d’ultrasons, lequel convient de préférence également pour un procédé impulsion-écho afin de générer des signaux d’émission et réceptionner des ondes acoustiques sur la base desquelles le transducteur d’ultrasons génère des signaux de réception. D’autres variantes possibles de la présente invention sont également données par les sous-revendications. D’autres avantages et caractéristiques de la présente invention apparaîtrons par ailleurs clairement dans la description qui suit d’exemples de réalisation à l’aide des figures.
Sur les figures, on peut voir : sur la figure 1A, schématiquement, une variante de réalisation d’un dispositif selon l’invention qui est réalisé et prévu pour exécuter le procédé selon l’invention et en particulier pour déterminer une vitesse d’écoulement d’un milieu en écoulement ; sur la figure 1B, schématiquement, un dispositif de traitement de signaux du dispositif de la figure 1A, grâce auquel des signaux de réception analogiques sont en particulier multipliés entre eux et additionnés avant d’être numérisés et soumis à un traitement ultérieur, afin de calculer une différence de temps de propagation; sur la figure 2A, un diagramme montrant la variation d'amplitude normalisée des signaux de réception pour un dispositif selon les figures 1A et 1B, lorsque celui-ci utilise des signaux d’émission conformément à l’équation 3 indiquée plus haut pour la génération d’ondes acoustiques; sur la figure 2B, un diagramme montrant la variation dans le temps de signaux Z(t) et N(t) obtenus à partir des signaux de réception de la figure 2A conformément à l’équation 1 ci-dessus; sur la figure 2C, des différences de temps de propagation superposées à partir des signaux des figures 2A et 2B avec 1000 variations de signaux dans le temps simulées, bruitées et décalées différemment, la différence réelle de temps de propagation ayant été de 10ns; sur la figure 3, schématiquement, une autre variante de réalisation d’un dispositif selon l’invention, où un dispositif de traitement de signaux présente un démultiplexeur et un multiplexeur ainsi qu’un convertisseur analogique/numérique afin de numériser des signaux de réception générés avant leur traitement ultérieur et en particulier leur multiplication et addition; sur la figure 4, un dispositif connu de l’état de la technique pour la détermination d’une vitesse d’écoulement d’un milieu en écoulement au moyen d’ondes acoustiques de surface.
Sur la vue tronquée de la figure 4, un dispositif (de mesure) connu par ailleurs est montré partiellement, lequel est réalisé et prévu pour la détermination de propriétés physiques et/ou chimiques d’un milieu M, en particulier pour la détermination ou mesure d’une vitesse d’écoulement du milieu M s’écoulant. Un guide d’ondes acoustique avec deux substrats 1, 2 en tant qu’éléments de guidage du guide d’ondes est une partie du dispositif de mesure, dans lequel le milieu M s’écoule à travers le guide d’ondes. Les substrats 1, 2 qui sont situés l’un en face de l’autre et dont les surfaces 11, 21 (intérieures) orientées l’une vers l’autre s’étendent parallèlement entre elles le long d’une direction d’extension principale du guide d’ondes, sont fabriqués en un matériau non piézoélectrique.
Ces substrats 1, 2 sont situés l’un en face de l’autre à une distance a, et sont bordés ici par deux plaques situées l’une en face de l’autre à cette distance a ou par des parties de paroi situées l’une en face de l’autre, d’un tube, lesquelles forment un espace intérieur 3 (en forme de canal) du guide d’ondes. Dans l’espace intérieur 3, le milieu M à mesurer, qui est représenté schématiquement par des lignes ondulées, est rempli, le milieu M pouvant s’écouler à travers l’espace intérieur 3. La direction d’écoulement du milieu M liquide ou fluide à travers l’espace intérieur 3 est en principe au choix. Ici, il y a un écoulement à partir d’une ouverture d’entrée vers une ouverture de sortie le long de la direction d’étendue principale du guide d’ondes et parallèlement aux surfaces intérieures 11, 21. La direction d’écoulement est indiquée sur la figure 4 par des flèches à l’ouverture d’entrée ou de sortie.
On associe aux deux (premier et deuxième) substrats 1 et 2 du guide d’ondes du dispositif une (première) unité d’émission et de réception SE1 ou une (deuxième) unité d’émission et de réception SE2. Dans ce oas, chaque unité d’émission et de réception SE1, SE2 peut être exploitée dans au moins deux modes de fonctionnement différents, d’une part en tant qu’émetteur et d’autre part en tant que récepteur, pour exciter ou réceptionner des ondes acoustiques. C’est ainsi que la première unité d’émission et de réception SE1 du deuxième substrat 2 peut par exemple (tout d’abord) être exploitée en tant qu’émetteur, alors que la deuxième unité d’émission et de réception SE2 du premier substrat 1 est exploitée en tant que récepteur.
Les unités d’émission et de réception SE1, SE2 sont respectivement disposées au niveau d’une surface extérieure 22 ou 12 du substrat 2 ou 1 respectif, laquelle est située en vis-à-vis d’une surface intérieure 21 ou 11 respectivement tournée vers l’espace intérieur 5 avec le milieu M. Concernant les deux unités d’émission et de réception SE1, SE2, il s’agit de préférence de convertisseurs piézoélectriques avec des électrodes interdigitées. La fixation d’une unité d'émission et de réception SE1, SE2 sur le substrat 2, 1 respectif s’effectue de préférence par collage, de sorte que celle-ci peut être montée de manière simple et aisée. En variante, on peut également prévoir d’autres types de fixation.
La première unité d’émission et de réception SE1 de la figure 4 se trouve ici dans la zone d’une première extrémité du guide d’ondes, alors que la deuxième unité d’émission et de réception SE2 est disposée dans la zone d’une autre, deuxième, extrémité du guide d’ondes et que le guide d’ondes s’étend entre ces deux extrémités le long d’une direction d’extension principale dans la vue en coupe transversale représentée.
Grâce à une unité d’émission et de réception, par exemple SE1, travaillant en tant qu’émetteur, des ondes acoustiques de surface OW2 sont générées dans le substrat 2 au moyen d’un signal d’émission prescrit, de préférence impulsionnel. Une partie de l’énergie de ces ondes acoustiques de surface OW2 générées est couplée dans le milieu M au niveau de la surface de délimitation de la surface intérieure 21 en tant qu’onde sonore volumique VW1. Sur la figure 4, la propagation d’une’ onde sonore volumique VW1 et la direction de propagation de l’onde sonore volumique VW1 est dans ce cas représentée schématiquement par une ligne en pointillés et par une flèche à côté de cette ligne en pointillés. Grâce à respectivement deux flèches en sens contraire à côté des lignes en pointillés, on exprime le fait que les ondes sonores volumiques se propagent le long de parties de chemin PA1 à PA7 représentées par les lignes en pointillés dans une direction dans . un mode de fonctionnement du dispositif, et dans l’autre direction dans un autre mode de fonctionnement du dispositif.
Les deux substrats 1, 2, lesquels forment les surfaces intérieures 12, 21, se composent de préférence d’un matériau non piézoélectrique et présentent une épaisseur d, laquelle est définie en tant que distance des surfaces intérieures 11,21 de la surface extérieure 12 ou 22 respectivement associée. Ici, l’épaisseur d est inférieure ou égale à la longueur d’onde des ondes acoustiques de surface respectivement générées. De ce fait, des ondes acoustiques de surface qui se propagent à l’intérieur des substrats 1, 2 peuvent présenter des propriétés d’onde de telle sorte qu’elles se propagent le long des surfaces intérieures 11, 21 tout comme également le long des surfaces extérieures 12, 22 des substrats 1, 2 en forme de plaques. On excite ainsi des ondes de Lamb ou des ondes dans la plage de transition entre des ondes de Lamb et des ondes de Rayleigh. En fonction de l’épaisseur d des plaques des substrats 1, 2, il y aura dans ce cas des ondes acoustiques de surface essentiellement sous la forme d’ondes de Lamb (d inférieur à la longueur d’onde des ondes acoustiques de surface) ou sous la forme d’ondes dans la plage de transition entre des ondes de Lamb et des ondes de Rayleigh (d égal à la longueur d’onde des ondes acoustiques de surface). Dans tous les cas, les ondes acoustiques de surface se propagent le long des deux surfaces 11, 12 ou 21, 22 des substrats 1 et 2.
Comme montré sur la figure 4, les ondes acoustiques de surface OW2 s’étendent ainsi par exemple en partant de l’unité d’émission et de réception SE1 travaillant en tant qu’émetteur, le long de la direction d’extension du deuxième substrat 2 et en particulier le long de la surface intérieure 21 de ce substrat 2. Une partie de l’énergie d’onde sonore des ondes acoustiques de surface OW2 se propageant le long de la surface intérieure 21 du substrat 2 est couplée dans le milieu M se trouvant à l’intérieur de l’espace intérieur 3, de sorte que les ondes sonores volumiques VW1 sont générées à l’intérieur du milieu M. Dans ce cas, une direction de propagation de ces ondes sonores volumiques WV1 couplées dedans est inclinée de l’ordre d’un angle caractéristique Δ, lequel n’est pas représenté, par rapport à la normale sur la surface 21 plate.
Les ondes sonores volumiques VW1 se propagent dans le milieu M respectivement le long d'un chemin P1. Ce chemin peut être subdivisé en différentes parties PA1, PA2, PA3, PA4, PA5, PA6, PA7 représentées par des lignes en pointillés. Chacune de ces parties de chemin s’étend entre un (deuxième) substrat 2 et l’autre (premier) substrat 1. Dès que l’onde sonore volumique VW1 est arrivée à un emplacement d’interaction au niveau de la surface intérieure 11 du substrat 1 situé en vis-à-vis, une partie de son énergie est couplée dans le substrat 1 de telle sorte que des ondes acoustiques de surface OW1 y sont générées, par exemple sous la forme d’ondes de Lamb ou d’ondes de surface dans la plage de transition entre des ondes de Lamb et des ondes de Rayleigh, lesquelles se propagent le long de ce substrat 1.
En outre, à chaque moment où l’onde sonore volumique VW1 atteint la surface intérieure 11 ou 21 d’un des substrats 1, 2, il y a une interaction de cette onde sonore avec le substrat 1,2 correspondant. Dans ce cas, il y a en règle générale un échange d’énergie, à base d’énergie acoustique entre le substrat 1, 2, en particulier l’onde de surface OW1, OW2 du substrat respectif 1, 2 et l’onde sonore volumique VW1. L’onde sonore volumique VW1 est au moins partiellement réfléchie et modifie ce faisant sa direction de propagation. Si l’interaction consiste en un couplage d’énergie à partir de l’onde sonore volumique VW1 dans l’onde de surface OW1 concernée, alors l’amplitude de l’onde de surface OW1 augmente du fait de ce couplage, et l’amplitude de l’onde sonore volumique VW1 diminue. En variante, en fonction des propriétés du substrat et du milieu M ainsi que de l’onde, il peut également y avoir un couplage d’énergie de l’onde de surface OW1 dans l’onde sonore volumique VW1.
Grâce à l’interaction de l’onde sonore volumique VW1 avec les substrats 1,2 le long de son chemin P1, plusieurs emplacements d’interaction sont ainsi définis. Au niveau de ces emplacements d’interaction, l’onde sonore volumique VW1 entre respectivement en interaction avec un substrat 1, 2 et les ondes de surface OW1, OW2 apparaissant dans le substrat 1, 2. C’est ainsi que des (premiers) trains d’ondes comprenant des ondes sonores volumiques VW1 se propagent globalement sur un chemin de propagation P1 sensiblement en forme de zigzag dans le milieu entre la première unité d’émission et de réception SE1 et la deuxième unité d’émission et de réception SE2 le long de la direction d’extension principale du guide d’ondes. En raison de l’interaction des ondes sonores volumiques VW1 avec le premier substrat 1 situé en vis-à-vis du deuxième substrat 2 au niveau de la surface intérieure 11 de celui-ci, il y a une excitation des ondes de surface OW1 qui se propagent sur le substrat 1 et peuvent finalement être réceptionnées au niveau de l’unité d’émission et de réception SE2 travaillant en tant que récepteur. Entre les emplacements d’interaction du deuxième substrat 2, à savoir ceux des emplacements où les ondes sonores volumiques VW1 interagissent avec le deuxième substrat 2, les ondes de surface OW1 se propagent sans renforcement, mais subissent un autre renforcement à l’emplacement d’interaction suivant (le cas échéant). Grâce à la mesure des trains d’ondes arrivant au niveau du récepteur E, en particulier des ondes de surface OW1, lesquelles ont été excitées par interaction avec les ondes sonores volumiques VW1, il est possible de déterminer le temps de propagation de trains d’ondes entre les première et deuxième unités d’émission et de réception SE1, SE2. C’est ainsi qu’à partir d’ondes acoustiques de surface OW1 (ou de groupes d’ondes de surface OW1), lesquelles arrivent les unes après les autres à l’unité d’émission et de réception SE2 travaillant en tant que récepteur, il est possible d’en déduire la vitesse du son à l’intérieur du milieu M, en particulier lorsque le temps de propagation des trains d’ondes entre les unités d’émission et de réception SE1, SE2 est déterminé. Etant donné que les temps de propagation mesurés des ondes acoustiques de surface OW1 couplées via les ondes sonores volumiques VW1 au niveau des emplacements d’interaction respectifs sont influencés par les propriétés du milieu M, il est ainsi possible de déterminer, grâce à une unité d’évaluation électronique (non représentée ici) à laquelle les signaux d’une unité d’émission et de réception SE2 sont retransmis, des propriétés physiques et/ou chimiques du milieu M à mesurer. Dans ce cas, on génère respectivement, à partir d’une onde acoustique réceptionnée à l’unité d’émission et de réception SE1 ou SE2, laquelle s’est propagée en tant qu’onde sonore volumique VW1 (ou VW2) au moins en partie à travers le milieu M à partir de l’autre unité d’émission et de réception SE1, SE2 jusqu’à cette unité d’émission et de réception SE2, SE1, un signal de réception.
Dans ce cas, il faut tenir compte du fait qu’au niveau d’un récepteur ou d’une unité d’émission et de réception exploitée en mode de réception SE2, ce sont généralement exclusivement des ondes de surface OW1 qui sont réceptionnées, ces ondes de surface OW1 relevant des ondes sonores volumiques VW1 touchant le premier substrat 1. Les différences de temps déterminées entre la réception d’une ou plusieurs onde(s) de surface sont également utilisées de manière correspondante en tant que base pour la détermination de la vitesse d’écoulement, comme cela va encore être expliqué plus en détails dans ce qui suit.
Pour ce faire, le dispositif de la figure 4 peut être exploité via un multiplexeur dans deux modes de fonctionnement différents, c’est-à-dire qu’également la deuxième unité d’émission et de réception SE2 précédemment exploitée dans un mode de réception est exploitée dans un mode d’émission, et que la première unité d’émission et de réception SE1 précédemment exploitée en mode d’émission est exploitée en mode de réception. Ce faisant, grâce à la deuxième unité d’émission et de réception SE2 exploitée en mode d’émission, une onde acoustique de surface OW1 est alors excitée dans le premier substrat 1. Au moins une partie de l’énergie de cette onde acoustique de surface OW1 est convertie comme précédemment en énergie d’une onde sonore volumique VW2 se propageant dans le milieu M, laquelle se propage alors sur un chemin de propagation P2 à partir de la deuxième unité d’émission et de réception SE2 jusqu’à la première unité d’émission et de réception SE1 à travers le milieu M. L’onde de surface OW1 excitée par la deuxième unité d’émission et de réception SE2 d’une deuxième onde acoustique part ainsi de cette unité d’émission et de réception SE2, de sorte que l’onde sonore volumique VW2 ainsi excitée se propage sensiblement à l’opposé de l’onde sonore volumique VW1 précédente et en forme de zigzag à travers le milieu M en direction de la première unité d’émission et de réception SE1.
Grâce à la commutation entre le mode d’émission et de réception des unités d’émission et de réception SE1, SE2, la direction de propagation des ondes sonores volumiques VW1, VW2 est par conséquent inversée à l’intérieur du milieu M le long de la direction de propagation principale du guide d’ondes. Il en résulte qu’en fonction du mode de fonctionnement, des (premiers ou deuxièmes) trains d’ondes, ou premières ou deuxièmes ondes acoustiques, se propagent d’une part entre les deux unités d’émission et de réception SE1 et SE2 de manière à ce qu’ils présentent des parties de chemin PA1 à PA7 s’étendant dans le milieu, le long desquels ils ont un vecteur de propagation avec une composante de vitesse de propagation vectorielle en direction de l’écoulement du milieu M (premiers trains d’ondes), et se propagent d’autre part de manière à ce qu’ils présentent des parties de chemin PA7 à PA1 s’étendant dans le milieu, le long desquelles ils ont un vecteur de vitesse de propagation avec une composante de vitesse de propagation vectorielle dans le sens contraire à l’écoulement du milieu M (deuxièmes trains d’ondes). Il en résulte que les temps de propagation des ondes de surface réceptionnées à la première unité d’émission et de réception SE1 d’un deuxième train d’ondes divergent des temps de propagation pour des ondes de surface réceptionnées à la deuxième unité d’émission et de réception SE2 d’un premier train d’ondes en raison de l’écoulement du milieu M. Grâce à une mesure des temps de propagation (absolus) des deux trains d’ondes et/ou à la détermination d’une différence des temps de propagation des trains d’ondes opposés l’un par rapport à l’autre, il est ainsi en principe possible de déterminer la vitesse d’écoulement du milieu M. En outre, grâce à la propagation des trains d’ondes le long du guide d’ondes, il est possible d’en déduire des informations sur la densité ou la concentration de substances dans le milieu M.
Le fait de parvenir, à l’aide de signaux réceptionnés pour des ondes acoustiques arrivant (signaux de réception), à une conclusion quant à une différence de temps de propagation ou même au temps de propagation absolu d’une onde acoustique entre les unités d’émission et de réception SE1, SE2 n’est cependant pas trivial et nécessite le cas échéant un énorme effort de calcul. Cela vaut par ailleurs également pour d’autres dispositifs de mesure avec lesquels, à l’aide d’ondes acoustiques, des propriétés physiques et/ou chimiques d’un milieu sont déterminées. C’est là qu’intervient alors un procédé selon l’invention, lequel prévoit entre autres l’émission de deux signaux d’émission à décalage de phase entre eux pour la génération d’ondes acoustiques et avec lequel il est possible de réaliser plus efficacement l’évaluation de signaux réceptionnés en particulier dans le cas d’un dispositif avec un guide d’ondes selon la figure 4, et en particulier de déterminer beaucoup plus facilement la détermination d’une différence de temps de propagation d’ondes acoustiques se propageant le long de directions différentes à l’aide de signaux réceptionnés. C’est ainsi qu’un dispositif de mesure selon la figure 4 connu par l’état de la technique est amélioré de telle sorte que des émetteurs du dispositif génèrent, avec deux signaux d’émission différents à fréquence de base identique et déphasage prescrit entre eux, respectivement des premières et deuxièmes ondes acoustiques le long des différentes directions de propagation, de préférence en direction d’écoulement et en opposition à la direction d’écoulement, de telle sorte que dans chaque direction de propagation, respectivement deux signaux de réception yô< yî y®, y®différents sont générés par les récepteurs du dispositif. On prévoit en outre un dispositif de traitement de signaux, par exemple SVA ou SVB, grâce auquel on détermine une différence de temps de propagation At et/ou un temps de propagation absolu t0 grâce aux quatre signaux de réception γθ, y’, y® et y® générés.
Dans une variante de réalisation de la figure 1A, on prévoit alors un générateur de signaux 4 couplé avec plusieurs unités d’émission et de réception SE1, SE2, SE3 et SE4 ou intégré dedans, de sorte que via deux unités d’émission et de réception SE1 et SE2 ainsi que SE3 et SE4 espacées spatialement l’une de l’autre le long de l’étendue longitudinale du guide d’ondes l’une par rapport à l’autre, il est respectivement possible d’utiliser un signal d’émission χΊ ou un signal d’émission x0 à décalage de phase par rapport à celui-ci, de même fréquence de base, pour la génération d’un onde acoustique (de surface).
Dans ce cas, de manière correspondante à l’équation 3 déjà mentionnée plus haut, les signaux d’émission x0, χΊ individuels présentent de préférence la forme suivante :
Ici également, les unités d’émission et de réception SE1 à SE4 individuelles sont réalisées de telle sorte qu’elles puissent générer des ondes acoustiques tout comme également réceptionner des ondes acoustiques. C’est ainsi qu’une unité d’émission et de réception SE1 par exemple, génère une onde acoustique se propageant en direction d’écoulement (première direction de propagation), laquelle est réceptionnée à l’unité d’émission et de réception SE2 située en aval et y génère un signal de réception y*· Cette unité d’émission et de réception SE2 génère à son tour une deuxième onde acoustique se propageant en amont (dans la deuxième direction de propagation opposée), laquelle est réceptionnée à l’unité d’émission et de réception SE1 située en amont et conduit à la génération d’un autre signal de
réception y®. Dans ce cas, les deux signaux de réception yir et y® relèvent d’un signal d’émission x-i.
On génère de la même manière, via le deuxième agencement d’émission et de réception SE3-SE4, avec des signaux d’émission x0, des signaux de réception ys0 et yro·
Grâce à un dispositif de traitement de signaux SVA selon la figure 1 B, on soumet les signaux de réceptions γθ, yir; y®, y® analogiques individuels à un traitement ultérieur. Dans ce cas, conformément à l’équation 6 ci-dessus, ce qui suit s’applique pour les signaux de réception individuels :
Grâce à l’expression du dispositif de traitement de signaux SVA montrée sur la figure 1 B, qui comprend ici un total de quatre multiplicateurs 5.1,5.2, 5.3 et 5.4 ainsi que deux additionneurs 7.1 et 7.2 et deux convertisseurs analogiques/numériques 8.1 et 8.2, les signaux de réception analogiques γθ, yi; y®, y®sont soumis à un traitement ultérieur et numérisés conformément à l’équation 1 ci-dessus pour donner des signaux formant une enveloppe Z(t) et N(t).
Via une unité de calcul, représentée ici en tant que moyen de traitement de signaux 9, on obtient alors à partir de ces signaux numérisés formant une enveloppe Z(t) et N(t), de manière correspondante à l’équation 2, des valeurs pour une différence de temps de propagation At entre les ondes acoustiques se propageant en amont et en aval. A partir de cela, en cas d’une distance l0 connue entre les unités d'émission et de réception SE1 et SE2 ou SE3 et SE4 individuelles ainsi que d’une vitesse du son c0 connue dans le milieu M qui se trouve dans l’espace intermédiaire 3, il est possible de parvenir à la conclusion d’une vitesse d’écoulement (moyenne) du milieu M conformément aux équations 4.1,4.2 et 5 présentées plus haut.
Il est par ailleurs bien entendu également possible d’obtenir un facteur de proportionnalité (l/(c0)2) avec un dispositif selon la figure 1A à partir d'un temps de
propagation absolu t0 à l’aide des quatre signaux de réception γθ, yj; y®, y® en ce que le dispositif de traitement de signaux SVA présente par exemple des moyens pour la mise en œuvre de l’équation 7 mentionnée plus haut.
La détermination de la différence de temps de propagation At avec un dispositif de traitement de signaux SVA frappe de par sa simplicité car elle permet, en utilisant des composants de technique de commutation simples, de déduire de manière fiable la différence de temps de propagation des ondes acoustiques dans les différentes directions de propagation, en particulier sans avoir besoin de recourir à une transformation de Fourier lors du traitement des signaux. A l’aide des figures 2A, 2B et 2C, des variations dans le temps des signaux, déterminées par simulation, doivent être représentées et en particulier la précision de mesure atteinte lors de la détermination d’une différence de temps de propagation At doit être montrée à l’aide de la figure 2C.
La figure 2A montre dans ce cas une variation d’amplitude normalisée pour les signaux de réception yo. yi; yo> y*individuels. Sur la figure 2B on représente de nouveau les signaux Z(t) et N(t) (voir équation 1) qui en résultent.
Sur la figure 2C, on représente alors les différences de temps de propagation At qui en sont déduites et déterminées en fonction du temps d’échantillonnage TAbt. tors de l’échantillonnage via le dispositif de traitement de signaux SVA. Un total de 1000 différentes variations dans le temps de signaux bruitées et décalées ont été simulées, et les résultats ont été superposés sur la figure 2C. Dans ce cas, on peut voir qu’avec la mise en œuvre du procédé selon l’invention, la valeur prescrite de différence de temps de propagation de 10ns est toujours bien atteinte au barycentre du signal, et ce avec la configuration comparativement simple et ainsi d’un coût avantageux du dispositif de traitement de signaux SVA.
Alors qu’avec la variante de réalisation selon la figure 1A, on prévoit une émission simultanée de l’ensemble- des quatre signaux d’émission ou une commande simultanée de quatre unités différentes d’émission et de réception SE1 à SE4 sous la forme de transducteurs d’ultrasons - de préférence intégrés dans un seul capteur à ultrasons d’un dispositif (de mesure) selon l’invention - avec deux signaux d’émission x0 et X! différents, il est bien entendu également possible de prévoir une émission décalée dans.Je temps de signaux d’émission x0 et χυ La figure 3 montre alors schématiquement une variante de réalisation possible.
Ici, une variante de réalisation d’un dispositif de mesure selon l’invention en particulier pour la détermination d’une vitesse d’écoulement (moyenne) vm d’un milieu M se trouvant dans un espace intérieur 3 est dotée d’un dispositif de traitement de signaux SVB en variante, lequel comprend entre autres un multiplexeur 40 et un démultiplexeur 90. Le générateur de signaux 4 prescrivant les signaux d’émission x0 et χΊ, lequel fait de nouveau partie ici d’un agencement avec deux unités d’émission et de réception SE1 et SE2, est couplé avec un multiplexeur 40 afin d’amener de manière décalée dans le temps aux unités d’émission et de réception ÔE1 et SE2 les signaux d’émission x0 et Xt individuels. C’est ainsi que chaque unité d’émission et de réception SE1 ou SE2 peut générer, avec différents signaux d’émission x0, χΊ, des premières ou deuxièmes ondes acoustiques dans la direction d’écoulement ou à l’encontre de la direction d’écoulement.
Les unités d’émission et de réception SE1 et SE2 amènent alors successivement au démultiplexeur 90 les signaux de réception yo> yr yo- yf sous f°rme analogique, et ceux-ci sont alors numérisés via un convertisseur analogique/numérique 8 couplé en aval. Ici, les signaux numérisés sont ensuite utilisés pour le calcul de la différence de temps de propagation At et/ou du temps de propagation absolu t0 dans le moyen de traitement de signaux 9, par exemple conformément aux équations 7, 8 et 10 représentées plus haut.
Avec les dispositifs selon l’invention et le procédé pouvant être mis en œuvre par ce biais, il est possible de déterminer les variations en fonction du temps du montant et de la phase de signaux ultrasonores transmis ou réfléchis et les temps de propagation de signal ou différences de temps de propagation qui en sont déduits, d’une manière particulièrement simple et efficace. L’utilisation du matériel peut en outre être limitée à un minimum. Le procédé mis en œuvre pour la détermination d’une différence de temps de propagation At et/ou un procédé pour la détermination d’un temps de propagation absolu tO est dans ce cas extrêmement robuste et s’accompagne d’une faible incertitude de mesure lorsque, comme représenté, on fait appel à quatre signaux de réception yo- yr yo- yî différents, lesquels relèvent de deux signaux d’émission x0, Xi différents de même fréquence de base f et d’un décalage de phase prescrit entre eux.
Pour compléter, il convient encore de remarquer que la séquence fréquentielle ou la variation dans le temps de la fréquence des signaux d’émission x0, Xi, est en effet de préférence identique pour chaque mesure effectuée, mais qu’il est possible de faire varier la fréquence de base f pour différentes mesures, c’est-à-dire en particulier différents agencements de mesure et/ou milieux.
On préfère par ailleurs en principe des signaux d’émission impulsionnels résultant en des signaux de réception impulsionnels yo (t)> yi 0); yo (t), yT CO manière correspondante à la figure 2A. En principe, une excitation continue d’ondes acoustiques au niveau des émetteurs respectifs ou unités d’émission et de réception SE1 à SE4 respectives serait cependant également envisageable.
Il est par ailleurs évident qu’avec les dispositifs représentés et les procédés mis en œuvre avec, il est possible de déterminer non seulement une vitesse (moyenne) d’écoulement vm d’un milieu M qui s’écoule, mais aussi - en variante ou de manière complémentaire - de réaliser une mesure de la concentration, de l’épaisseur, de la distance, de la température et/ou du niveau de remplissage, comme cela est déjà bien connu pour des dispositifs comparables qui travaillent avec des ondes acoustiques, en particulier des ondes ultrasonores. La même chose vaut pour la détermination de variations en fonction du temps de l’amplitude momentanée ou de la phase momentanée des ondes acoustiques réceptionnées et des valeurs qui en sont déduites pour la densité et la viscosité du milieu M.
Liste des numéros de référence 1 (premier) substrat 11 surface intérieure 12 surface extérieure 2 (deuxième) substrat 21 surface intérieure 22 surface extérieure 3 espace intérieur 4 générateur de signaux 40 multiplexeur 5.1 - 5.4 multiplicateur 6 inverseur 7.1,7.2 additionneur 8 convertisseur analogique/numérique (convertisseur A/N) 8.1, 8.2 convertisseur analogique/numérique (convertisseur A/N) 9 moyen de traitement de signaux 90 démultiplexeur a distance d épaisseur f fréquence de base M milieu OW1, OW2 onde sonore de surface P1, P2 chemin PA1 - PA7 partie de chemin SE1, SE2, SE3, SE4 unité d’émission et de réception SVA, SVb dispositif de traitement de signaux t0 temps de propagation vm vitesse moyenne d'écoulement VW1,VW2 onde sonore volumique w fonction de fenêtre x signal d’émission yr premier signal de réception (en aval) ys deuxième signal de réception (en amont) φ angle de départ ^t différence de temps de propagation
Claims (25)
- Revendications1. Procédé pour la détermination de propriétés physiques et/ou chimiques d’un milieu sur la base d’au moins une première et une deuxième onde acoustique, lesquelles se sont respectivement propagées au moins partiellement à travers le milieu à partir d’au moins un émetteur vers au moins un récepteur, dans lequel - les ondes acoustiques sont respectivement générées par un signal d’émission (x0, x^, - la première onde acoustique se propage à travers le milieu (M) le long d’une première direction de propagation, et la deuxième onde acoustique se propage à travers le milieu (M) le long d’une deuxième direction de propagation différente de la première direction de propagation, et - pour une première onde acoustique réceptionnée après sa propagation à travers le milieu (M) le long de la première direction de propagation, un premier signal de réception (y', y J) est généré, et pour une deuxième onde acoustique réceptionnée après sa propagation à travers le milieu (M) le long de la deuxième direction de propagation, un deuxième signal de réception (y* , y®) est généré, dans lequel on détermine, à partir des signaux de réception (y', yir; y®, y®) , une différence de temps de propagation (Δΐ) des ondes acoustiques, et au moyen d’une différence de temps de propagation (Δί) déterminée, des propriétés physiques et/ou chimiques du milieu (M) étant déterminées, caractérisé en ce que - au moins deux signaux d’émission (x0, Xi) différents sensiblement de même fréquence de base (f) sont utilisés, lesquels présentent un déphasage prescrit entre eux et par lesquels sont générés respectivement la première et la deuxième onde acoustique de telle sorte que pour chaque direction de propagation, deux signaux de réception ( y', y' ; y®, y® ) sont générés, lesquels relèvent de signaux d’émission (x0, Xi) différents, et - une différence de temps de propagation (At) des ondes acoustiques se propageant le long des deux directions de propagation différentes est déterminée grâce aux quatre signaux de réception
- 2. Procédé selon la revendication 1, caractérisé en ce que les quatre signaux de réceptionentrent directement dans un calcul pour calculer une différence de temps de propagation (At), et une différence de temps de propagation (At) étant seulement calculée en cas de présence de l’ensemble des quatre signaux de réception i
- 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que, pour la détermination d’une différence de temps de propagation (At), des signaux de réceptioni pour des ondes acoustiques sont multipliés entre eux, lesquels ont été générés par des signaux d’émission (x0, Xi) différents.
- 4. Procédé selon la revendication 3, caractérisé en ce que l’on prévoit une multiplication de signaux de réception (γθ, yj; y®, y®) analogiques, avec un calcul qui s’ensuit d’une différence de temps de propagation (At) à partir de signaux numérisés, ou bien en ce que l’on prévoit, pour un calcul d’une différence de temps de propagation (At), une multiplication de signaux de réception (γθ, y'; y®, y®) numérisés.
- 5. Procédé selon l’une des revendications 1 à 4, caractérisé en ce que l’on utilise, pour la détermination d’une différence de temps de propagation (At), une formule de la forme suivante:où t est le temps et Yo étant le premier signal de réception pour l’onde acoustique générée par un premier signal d’émission (x0), Yi étant le premier signal de réception pour l’onde acoustique générée par un deuxième signal d’émission (x^, Yoétant le deuxième signal de réception pour l’onde acoustique générée dans l’autre direction de propagation par un premier signal d’émission (x0) yi et étant le deuxième signal de réception pour l’onde acoustique générée dans l’autre direction de propagation par un deuxième signal d’émission (Xi)·
- 6. Procédé selon la revendication 5, caractérisé en ce que l’on détermine une différence de temps de propagation At d’après la formule suivante :où f est ici une fréquence de base pour les deux signaux d’émission (x0, Xi).
- 7. Procédé selon la revendication 6, caractérisé en ce que les deux signaux d’émission (x0, Xi) différents présentent des fréquences de base (f0, fi) divergeant peu l’une de l’autre, et en ce qu’une fréquence de base f pour la détermination de la différence de temps de propagation At est déterminée à partir de la moyenne arithmétique des deux fréquences de base (f0, fi).
- 8. Procédé selon l’une des revendications précédentes, caractérisé en ce que les signaux d’émission (x0, Xi) comprennent respectivement une impulsion oscillatoire et en ce que des enveloppes pour les impulsions oscillatoires individuelles des signaux d’émission (x0, Xi) à décalage de phase entre eux sont sensiblement identiques entre elles.
- 9. Procédé selon l’une des revendications précédentes, caractérisé en ce que des ondes acoustiques sont générées simultanément ou les unes après les autres par les signaux d’émission (x0, Xi).
- 10. Procédé selon l’une des revendications précédentes, caractérisé en ce que l’on détermine, au moyen de la différence de temps de propagation (At), une vitesse d’écoulement (vm) d’un milieu (M) en écoulement, et en ce que les ondes acoustiques sont générées par les signaux d’émission (x0, Xi) d’une manière telle, que des premières ondes acoustiques se propagent à travers le milieu (M) le long d’une première direction de propagation dans la direction d’écoulementdu milieu (M), et que des deuxièmes ondes acoustiques se propagent à travers le milieu (M) le long d’une deuxième direction de propagation à l’encontre de la direction d’écoulement du milieu (M).
- 11. Procédé selon la revendication 10, caractérisé en ce que l’on calcule une vitesse d’écoulement grâce à la différence de temps de propagation (Δί) et à un facteur de proportionnalité (l/(c0)2) connu.
- 12. Procédé pour la détermination de propriétés physiques et/ou chimiques d’un milieu sur la base d’au moins une première et une deuxième onde acoustique, lesquelles se sont respectivement propagées au moins partiellement à travers le milieu à partir d’au moins un émetteur vers au moins un récepteur, dans lequel - les ondes acoustiques sont respectivement générées par un signal d’émission (x0, Xi), - la première onde acoustique se propage à travers le milieu (M) le long d’une première direction de propagation, et la deuxième onde acoustique se propage à travers le milieu (M) le long d’une deuxième direction de propagation identique à la première direction de propagation ou différente de celle-ci, et - pour une première onde acoustique réceptionnée après sa propagation à travers le milieu (M) le long de la première direction de propagation, un premier signal de réception (γθ, y') est généré, et pour une deuxième onde acoustique réceptionnée après sa propagation à travers le milieu (M) le long de la deuxième direction de propagation, un deuxième signal de réception (y® , y®) est généré, dans lequel on détermine, à partir des signaux de réception (γθ, y J; y®, y®), un temps de propagation absolu (t0) d’une onde acoustique d’un émetteur (SE1, SE2, SE3, SE4) vers un récepteur (SE2, SE1, SE4, SE3), et au moyen d’un temps de propagation absolu (t0) déterminé, des propriétés physiques et/ou chimiques du milieu (M) étant déterminées, caractérisé en ce que - au moins deux signaux d’émission (x0, Xi) différents sensiblement de même fréquence de base (f) sont utilisés, lesquels présentent un déphasage prescrit entre eux et avec lesquels on génère respectivement une première et une deuxième onde acoustique de telle sorte que, pour chaque direction de propagation, deux signaux de réceptionsont générés, lesquels relèvent de signaux d’émission (x0, Xi) différents, et - un temps de propagation absolu (t0) d’une onde acoustique entre un émetteur (SE1, SE2, SE3, SE4) et un récepteur (SE2, SE1, SE4, SE3) est déterminé grâce aux quatre signaux de réception
- 13. Procédé selon la revendication 12, caractérisé en ce que l’on détermine un temps de propagation absolu (t0) à partir d’une valeur de départ estimée ( t0 ) et d’une valeur de différence (At’) calculée à partir des quatre signaux de réception
- 14. Procédé selon la revendication 13, caractérisé en ce que, pour le calcul de la valeur de différence (At’), on utilise au moins l’un des termes suivants:f étant ici une valeur estimée pour la fréquence centrale des signaux de réception (γθ, yj; y®, yf) et tk une constante de temps prescrite avec tk > 0 et yS étant le premier signal de réception pour l’onde acoustique générée par un premier signal d’émission (x0), y® étant le premier signal de réception pour l’onde acoustique généré par un deuxième signal d’émission (xO, Yo étant le deuxième signal de réception pour l’onde acoustique générée dans l’autre direction de propagation par un premier signal d’émission (x0) YÏ etétant le deuxième signal de réception pour l’onde acoustique générée dans l’autre direction de propagation par un deuxième signal d’émission (xi)·
- 15. Procédé selon la revendication 13 ou 14, caractérisé en ce que l’on prescrit une valeur de départ ( t0 ) pour le calcul du temps de propagation au moyen d’au moins une enveloppe de signal d’une paire de signaux de réception (γθ, yj; γθ, yS ) , dans lequel les signaux de réception ide la paire de signaux de réception relèvent de signaux d’émission (x0, χΊ) différents, mais d’ondes acoustiques s’étant propagées dans la même direction de propagation.
- 16. Procédé selon la revendication 15, caractérisé en ce qu’une enveloppe de signal yhuei est donnée grâce à la formule suivante oùy0 est un premier signal de réception pour l’onde acoustique générée par un premier signal d’émission (x0) et y! est un deuxième signal de réception pour l’onde acoustique générée par un deuxième signal d’émission (x^.
- 17. Procédé selon l’une des revendications 11 à 16, caractérisé en ce qu’au moyen du temps de propagation absolu (t0) déterminé, on détermine une vitesse d’écoulement (vm) d’un milieu (M) s’écoulant et en ce que les ondes acoustiques sont générées par les signaux d’émission (x0, x^ de telle sorte que des premières ondes acoustiques se propagent à travers le milieu (M) le long d’une première direction de propagation dans la direction d’écoulement du milieu (M) et que des deuxièmes ondes acoustiques se propagent à travers le milieu (M) le long d’une deuxième direction de propagation à l’encontre de la direction d’écoulement du milieu (M), dans lequel on calcule, à l’aide du temps de propagation absolu (t0) déterminé, un facteur de proportionnalité (l/(c0)2), et avec le facteur de proportionnalité (l/(c0)2) et une différence de temps de propagation (At) mesurée, en particulier avec une différence de temps de propagation (At) déterminée selon un procédé des revendications 1 à 11 entre au moins une première et une deuxième onde acoustique, lesquelles se sont propagées aumoins partiellement à travers le milieu (M) en écoulement le long de différentes directions de propagation, une vitesse d’écoulement (vm) étant déterminée.
- 18. Procédé selon l’une des revendications 1 à 17 précédentes, caractérisé en ce que, lors d’une détermination de plusieurs valeurs pour des différences de temps de propagation (Δί) et/ou pour des temps de propagation absolus (t0) d’une pluralité d’onde acoustiques, on effectue le calcul de la moyenne et/ou une intégration des différences de temps de propagation (Δΐ) et/ou des temps de propagation (t0) déterminés.
- 19. Procédé selon l'une des revendications 1 à 18 précédentes, caractérisé en ce que, lors du procédé respectif, on détermine au moyen d’une différence de temps de propagation (Δί) déterminée et/ou ou moyen d’un temps de propagation absolu (t0) déterminé, une vitesse du son, une concentration, unedensité, une température, un niveau de remplissage du milieu dans un espace intérieur, une épaisseur d’une paroi jouxtant le milieu et/ou une distance de deux parties de paroi jouxtant le milieu.
- 20. Dispositif pour la détermination de propriétés physiques et/ou chimiques d’un milieu, en particulier pour la mise en œuvre d’un procédé selon l’une des revendications 1 à 19, avec - un guide d’ondes acoustique, lequel présente un espace intérieur (3) à remplir avec le milieu (M), - au moins deux émetteurs (SE1, SE2; SE3, SE4), lesquels sont respectivement réalisés et prévus pour la génération d’ondes acoustiques dans le guide d’ondes grâce à un signal d’émission (x0, Xi) de l’émetteur (SE1, SE2; SE3, SE4) de telle sorte qu’une onde acoustique générée se propage le long du guide d’ondes et au moins partiellement à travers le milieu (M), - au moins deux récepteurs (SE2, SE1; SE4, SE3), lesquels sont réalisés et prévus de manière espacée spatialement l’un de l’autre le long du guide d’ondes et respectivement pour la génération d’un premier ou deuxième signal de réception (γθ, y'; y®, y®) pour une première ou deuxième onde acoustique atteignant le récepteur (SE2, SE1 ; SE4, SE3) respectif, dans lequel une première onde acoustique se propage respectivement le long d’une première direction de propagation et pour ce faire, un premier signal de réception (γθ, yi) étant généré au niveau d’un récepteur (SE2, SE1 ; SE4, SE3), et une deuxième onde acoustique se propageant respectivement le long d’une deuxième direction de propagation identique à la première direction de propagation ou différente de celle-ci et pour ce faire, un deuxième signal de réception (y®, y®) étant généré au niveau d’un récepteur (SE2, SE1; SE4, SE3), et - un dispositif de traitement de signaux (SVA, SVB), à l’aide duquel on détermine, à partir des signaux de réception i• une différence de temps de propagation (At) entre des premières et deuxièmes ondes acoustiques réceptionnées lorsque la deuxième direction de propagation est différente de la première direction de propagation, • un temps de propagation absolu (t0) d’une onde acoustique d’un émetteur (SE1, SE2; SE3, SE4) vers un récepteur (SE 2, SE1; SE4, SE3) lorsque la deuxième direction de propagation est différente de la première direction de propagation ou identique à celle-ci, et lequel est réalisé et prévu pour déterminer, avec une différence de temps de propagation (At) déterminée et/ou un temps de propagation absolu (t0) déterminé, des propriétés physiques et/ou chimiques du milieu (M), caractérisé en ce que - les émetteurs (SE1, SE2; SE3, SE4) du dispositif sont réalisés et prévus pour générer, avec deux signaux d’émission (x0, Xi) différents sensiblement de même fréquence de base (f) et à décalage de phase prescrit entre eux, respectivement des premières et deuxièmes ondes acoustiques de telle sorte que dans la première direction de propagation, deux premiers signaux de réception différenciables, et dans la deuxième direction de propagation, deux deuxièmes signaux de réception ( γθ, y'; y®, y®) différenciables sont générés par les récepteurs (SE2, SE1, SE4, SE3) du dispositif, et- le dispositif de traitement de signaux (SVA, SVB) est réalisé et prévu pour déterminer une différence de temps de propagation (At) et/ou un temps de propagation absolu (t0) grâce aux quatre signaux de réception (γθ, y^; y®, y®).
- 21. Dispositif selon la revendication 20, caractérisé en ce que les émetteurs (SE1, SE2; SE3, SE4) sont réalisés et prévus pour générer des ondes acoustiques de surface (OW1, OW2) dans le guide d’ondes.
- 22. Dispositif selon la revendication 20 ou 21, caractérisé en ce que le dispositif est réalisé et prévu pour la détermination d’une vitesse d’écoulement (vm) d’un milieu (M) s’écoulant à travers l’espace intérieur (3) du guide d’ondes.
- 23. Dispositif selon l’une des revendications 20 à 22, caractérisé en ce que le dispositif de traitement de signaux (SVA, SVB) présente des moyens (5.1, 5.2, 5.3, 5.4) pour la multiplication de signaux de réception analogiques (γθ, y'; y®, y®) et des moyens (9) pour le calcul qui s’ensuit d’une différence de temps de propagation (At) à partir de signaux numérisés, ou des moyens (8) pour la numérisation des signaux de réception (γθ, y'; y®, y®) et des moyens (9) pour la multiplication qui s’ensuit de signaux de réception numérisés pour le calcul d’une différence de temps de propagation (At).
- 24. Dispositif selon l'une des revendications 20 à 23, caractérisé en ce que la fréquence de base (f) des deux signaux d’émission (x0, x^ est variable.
- 25. Dispositif selon l’une des revendications 20 à 22, caractérisé en ce que le déphasage des deux signaux d’émission (x0, x^ est variable entre 45° et 135°, et est en particulier de 90°.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE1020151022001 | 2015-02-16 | ||
DE102015102200.1A DE102015102200B4 (de) | 2015-02-16 | 2015-02-16 | Verfahren zur Bestimmung von Eigenschaften eines Mediums und Vorrichtung zur Bestimmung von Eigenschaften eines Mediums |
Publications (2)
Publication Number | Publication Date |
---|---|
FR3032799A1 FR3032799A1 (fr) | 2016-08-19 |
FR3032799B1 true FR3032799B1 (fr) | 2019-10-25 |
Family
ID=56552196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1650833A Active FR3032799B1 (fr) | 2015-02-16 | 2016-02-02 | Procede pour la determination de proprietes d'un milieu et dispositif pour la determination de proprietes d'un milieu |
Country Status (4)
Country | Link |
---|---|
US (1) | US10605779B2 (fr) |
CN (1) | CN105891326B (fr) |
DE (1) | DE102015102200B4 (fr) |
FR (1) | FR3032799B1 (fr) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102017214314A1 (de) * | 2017-08-17 | 2019-02-21 | Robert Bosch Gmbh | Verfahren zum Erkennen einer Position eines mobilen Roboterfahrzeugs sowie Roboterfahrzeug und System |
DE102019218387B4 (de) | 2019-11-27 | 2021-12-09 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Messanordnung und Verfahren zur Bestimmung von Eigenschaften eines zu extrudierenden Materials während der Durchführung eines Schnecken-Extrusionsprozesses |
CN113375737B (zh) * | 2020-06-15 | 2024-01-12 | 郑州大学 | 时差式超声气体流量计的流速计量方法 |
CN112161557B (zh) * | 2020-08-31 | 2023-04-11 | 武汉轻工工程技术有限公司 | 一种水平位移测量方法 |
Family Cites Families (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3575050A (en) | 1968-12-04 | 1971-04-13 | Panametrics | Fluid flowmeter |
US3935735A (en) | 1974-09-03 | 1976-02-03 | Badger Meter, Inc. | Ultrasonic flow meter |
JPS5829854B2 (ja) | 1977-07-26 | 1983-06-25 | 富士電機株式会社 | 超音波式測定装置 |
DE3046254A1 (de) | 1980-12-08 | 1982-07-15 | Siemens AG, 1000 Berlin und 8000 München | Ultraschallstroemungsmesser |
DE3347420C2 (de) | 1983-12-29 | 1987-03-19 | Nukem Gmbh, 6450 Hanau | Verfahren zur Messung des zeitlichen Abstandes von zwei elektrischen Impulsen und Vorrichtung zur Durchführung des Verfahrens |
CH672187A5 (fr) | 1986-10-23 | 1989-10-31 | Landis & Gyr Gmbh | |
DE3823177A1 (de) | 1988-07-08 | 1990-01-11 | Flowtec Ag | Ultraschall-durchflussmessung unter auswertung von phasenverschiebungen |
DE58905910D1 (de) | 1988-07-08 | 1993-11-18 | Flowtec Ag | Verfahren und anordnung zur durchflussmessung mittels ultraschallwellen. |
JPH09311126A (ja) * | 1996-05-24 | 1997-12-02 | Toshiba Corp | 撹拌機能付き濃度計装置 |
US6390999B1 (en) | 1996-06-28 | 2002-05-21 | Rocky Mountain Research, Inc. | Method and apparatus for flow measurement with temperature and density compensation |
US6062091A (en) | 1997-04-22 | 2000-05-16 | Baumoel; Joseph | Method and apparatus for determining ultrasonic pulse arrival in fluid using phase correlation |
FR2781048B1 (fr) | 1998-07-10 | 2000-09-15 | Faure Herman | Mesures croisees des signaux acoustiques d'un debitmetre |
DE19843806A1 (de) | 1998-09-24 | 2000-03-30 | Gerd Stange | Verfahren und Vorrichtung zur Durchflußmessung mit schallinduzierter Frequenzmischung in aktiven Meßaufnehmern |
DE19957905C5 (de) | 1999-10-01 | 2011-12-01 | Ge Inspection Technologies Gmbh | Verfahren und Vorrichtung zur zerstörungsfreien Prüfung von Schweißnähten mittels Ultraschall |
US6435037B1 (en) | 2000-01-06 | 2002-08-20 | Data Sciences International, Inc. | Multiplexed phase detector |
KR100349504B1 (ko) | 2000-04-24 | 2002-08-21 | 주식회사 창민테크 | 초음파 유속 측정장치 |
US20020062681A1 (en) * | 2000-11-30 | 2002-05-30 | Livingston Richard A. | Oxygen sensor and flow meter device |
JP4135056B2 (ja) | 2001-02-15 | 2008-08-20 | 横河電機株式会社 | 超音波流量計 |
US7007546B2 (en) * | 2003-12-22 | 2006-03-07 | Biode Inc. | Measurement, compensation and control of equivalent shear rate in acoustic wave sensors |
US8214168B2 (en) | 2004-09-07 | 2012-07-03 | Transonic Systems, Inc. | Noninvasive testing of a material intermediate spaced walls |
DE102005015456A1 (de) | 2005-04-04 | 2006-10-05 | Viasys Healthcare Gmbh | Verfahren zur Bestimmung der zeitlichen Lage eines Wellenpakets sowie Flussmessgerät |
WO2008034878A2 (fr) | 2006-09-20 | 2008-03-27 | Fachhochschule Coburg | procédé et dispositif de détermination des caractéristiques d'un milieu sous forme de liquide ou de matériau mou |
US7270015B1 (en) * | 2006-11-29 | 2007-09-18 | Murray F Feller | Thermal pulsed ultrasonic flow sensor |
CN101373951A (zh) | 2007-08-23 | 2009-02-25 | 李代甫 | 一种数字鉴相的方法及其装置 |
WO2010034713A2 (fr) | 2008-09-26 | 2010-04-01 | Hochschule Für Angewandte Wissenschaften Fachhochschule Coburg | Procédé et dispositif pour déterminer des caractéristiques d'un milieu |
DE102008043956A1 (de) | 2008-11-21 | 2010-05-27 | Robert Bosch Gmbh | Verfahren zur Bestimmung einer Strömungsgeschwindigkeit eines fluiden Mediums |
US8085156B2 (en) * | 2009-04-08 | 2011-12-27 | Rosemount Inc. | RF cavity-based process fluid sensor |
DE102009003020A1 (de) | 2009-05-12 | 2010-11-18 | Robert Bosch Gmbh | Laufzeitmessungskorrektur bei einem Strömungssensor |
EP2336732A1 (fr) * | 2009-12-15 | 2011-06-22 | Kamstrup A/S | Compteur de consommation et partie d'écoulement et boîtier formés par une structure polymère monolithique |
DE102009049067A1 (de) | 2009-10-12 | 2011-04-14 | Robert Bosch Gmbh | Verfahren und Vorrichtung zur verbesserten Ultraschall- Laufzeitdifferenzmessung |
KR101844098B1 (ko) * | 2010-09-03 | 2018-03-30 | 로스 알라모스 내셔널 씨큐어리티 엘엘씨 | 파이프 내 유체의 음향 특성을 비침투적으로 결정하는 방법 |
CN201828295U (zh) | 2010-10-28 | 2011-05-11 | 利尔达科技有限公司 | 一种相位差式超声波流量计量装置 |
DE102010063937A1 (de) | 2010-12-22 | 2012-06-28 | Sonotec Ultraschallsensorik Halle Gmbh | Verfahren und Vorrichtung zur Messung von Laufzeiten von Ultraschallimpulsen |
FR2972801B1 (fr) | 2011-03-14 | 2017-04-14 | Inst Polytechnique Grenoble | Systeme d'inspection par emission-reception d'ondes acoustiques |
US8700344B2 (en) * | 2011-04-20 | 2014-04-15 | Neptune Technology Group Inc. | Ultrasonic flow meter |
US9689726B2 (en) | 2012-07-09 | 2017-06-27 | Texas Instruments Incorporated | Flow meter |
DE102012212901A1 (de) | 2012-07-24 | 2014-01-30 | Continental Automotive Gmbh | Verfahren und Vorrichtung zum Betreiben einer akustischen Messvorrichtung |
CN104236647A (zh) | 2013-10-30 | 2014-12-24 | 安徽恩米特微电子股份有限公司 | 基于单片机解决方案的超声波流量计 |
US9424674B2 (en) * | 2014-04-01 | 2016-08-23 | Saudi Arabian Oil Company | Tomographic imaging of multiphase flows |
-
2015
- 2015-02-16 DE DE102015102200.1A patent/DE102015102200B4/de active Active
-
2016
- 2016-02-02 FR FR1650833A patent/FR3032799B1/fr active Active
- 2016-02-12 US US15/043,428 patent/US10605779B2/en active Active
- 2016-02-16 CN CN201610087772.2A patent/CN105891326B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
FR3032799A1 (fr) | 2016-08-19 |
US20160238570A1 (en) | 2016-08-18 |
CN105891326B (zh) | 2020-11-06 |
CN105891326A (zh) | 2016-08-24 |
DE102015102200A1 (de) | 2016-08-18 |
US10605779B2 (en) | 2020-03-31 |
DE102015102200B4 (de) | 2022-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR3032799B1 (fr) | Procede pour la determination de proprietes d'un milieu et dispositif pour la determination de proprietes d'un milieu | |
EP0261736B1 (fr) | Dispositif d'exploration par echographie ultrasonore d'un milieu en mouvement | |
EP1097354B1 (fr) | Mesures croisees des signaux acoustiques d'un debitmetre | |
EP0225667B1 (fr) | Appareil d'exploration par échographie ultrasonore de milieux en mouvement et notamment d'écoulements sanguins ou d'organes tels que le coeur | |
AU2016387414B2 (en) | Improved beam shaping acoustic signal travel time flow meter | |
FR2472755A1 (fr) | Appareil de mesure de la vitesse d'un courant de fluide par ultrasons | |
CA2153192C (fr) | Dispositif et methode de mesure de profil de vitesse dans un fluide polyphasique | |
US9645003B2 (en) | Efficient dispersion correction for FMCW-radar in a pipe | |
FR2701111A1 (fr) | Débitmètre de fluide. | |
EP0561690B1 (fr) | Procédé et dispositif de détermination du passage à une distance présélectionnée d'un point réflecteur à l'aide du temps de propagation d'une onde continue | |
FR2599503A1 (fr) | Procede de mesure des parametres de la phase solide d'une suspension et dispositif mettant en oeuvre ce procede | |
JP2836989B2 (ja) | ドップラ流速計 | |
JP2019070647A (ja) | 連続波レーダ測定によって媒体の充填レベルを特定するための方法および充填レベル測定装置 | |
FR2784182A1 (fr) | Debitmetre doppler a frequences multiples | |
FR2790554A1 (fr) | Methode et dispositif de mesure du temps de propagation d'un signal, en particulier d'un signal ultrasonore | |
CN108474766B (zh) | 声学确定介质特性的方法和借助反射元件声学确定介质特性的设备 | |
CN108496075B (zh) | 用于确定介质特性的方法和用于确定介质特性的设备 | |
CA2278408C (fr) | Procede et dispositif pour mesurer la vitesse d'ecoulement d'une veine fluide | |
EP0769152A1 (fr) | Procede et dispositif pour mesurer la vitesse d'un mobile | |
CH433781A (fr) | Dispositif de mesure de l'épaisseur d'un échantillon | |
EP0162987B1 (fr) | Procédé et installation pour mesurer le debit d'un fluide par ultrasons | |
FR3050828B1 (fr) | Procede de determination d'une vitesse d'ecoulement d'un fluide s'ecoulant dans un troncon de conduite et dispositif associe | |
JP2005337937A (ja) | 気泡センサ | |
US12123758B2 (en) | Flow sensor performing multi level down sampling and method thereof | |
US20220251948A1 (en) | Flow sensor performing multi level down sampling and method thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLFP | Fee payment |
Year of fee payment: 3 |
|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLFP | Fee payment |
Year of fee payment: 5 |
|
PLFP | Fee payment |
Year of fee payment: 6 |
|
PLFP | Fee payment |
Year of fee payment: 7 |
|
CD | Change of name or company name |
Owner name: ENDRESS+HAUSER FLOW DEUTSCHLAND AG, DE Effective date: 20220119 |
|
PLFP | Fee payment |
Year of fee payment: 8 |
|
PLFP | Fee payment |
Year of fee payment: 9 |