FR2958041A1 - Method for time domain reflectometry type measuring of water content in natural material e.g. soil, involves determining primary instant as instant in which number of digital samples is increased continuously above average value - Google Patents

Method for time domain reflectometry type measuring of water content in natural material e.g. soil, involves determining primary instant as instant in which number of digital samples is increased continuously above average value Download PDF

Info

Publication number
FR2958041A1
FR2958041A1 FR1052094A FR1052094A FR2958041A1 FR 2958041 A1 FR2958041 A1 FR 2958041A1 FR 1052094 A FR1052094 A FR 1052094A FR 1052094 A FR1052094 A FR 1052094A FR 2958041 A1 FR2958041 A1 FR 2958041A1
Authority
FR
France
Prior art keywords
curve
instant
point
inflection
digital samples
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR1052094A
Other languages
French (fr)
Other versions
FR2958041B1 (en
Inventor
Albert Marchiol
Daniel Meline
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Agence Nationale pour la Gestion des Dechets Radioactifs ANDRA
Original Assignee
Agence Nationale pour la Gestion des Dechets Radioactifs ANDRA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Agence Nationale pour la Gestion des Dechets Radioactifs ANDRA filed Critical Agence Nationale pour la Gestion des Dechets Radioactifs ANDRA
Priority to FR1052094A priority Critical patent/FR2958041B1/en
Publication of FR2958041A1 publication Critical patent/FR2958041A1/en
Application granted granted Critical
Publication of FR2958041B1 publication Critical patent/FR2958041B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • G01N33/246Earth materials for water content
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N22/00Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
    • G01N22/04Investigating moisture content

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The method involves determining a constant average value, which represents a noise level of a digital signal, from digital samples of the signal for determining a departure instant (tdi), where duration separating the departure instant and an end instant (tfi) measures a return-trip propagation time of the signal. The primary instant as an instant in which a number (N) of the digital samples is increased continuously above the average value is determined, where N is an integer number higher than or equal to 2. An independent claim is also included for a device for measuring water content of a medium.

Description

1 PROCEDE DE MESURE DE TENEUR EN EAU D'UN MILIEU ET DISPOSITIF ASSOCIE 1 METHOD FOR MEASURING THE WATER CONTENT OF A MEDIUM AND ASSOCIATED DEVICE

DESCRIPTION Domaine technique et art antérieur L'invention concerne un procédé de mesure de teneur en eau d'un milieu ainsi qu'un dispositif de mesure de teneur en eau qui met en oeuvre le procédé de l'invention. L'invention s'applique avantageusement à la mesure de la teneur en eau dans des matériaux naturels (sols, sables, argiles, argilites, etc.), cimentaires (ciments, mortiers, bétons), ou remaniés. Le procédé de l'invention se situe dans le cadre général des mesures effectuées selon le principe de la réflectométrie dans le domaine temporel plus communément appelées mesures TDR (TDR pour « Time Domain Reflectometry »). Une mesure TDR consiste à mesurer le temps que met une onde électromagnétique pour se propager le long d'électrodes placées dans un milieu. La connaissance de ce temps permet de déterminer la vitesse de propagation de l'onde dans le milieu et, en conséquence, la permittivité relative du milieu. Dans le domaine des très hautes fréquences, la permittivité relative du milieu dépendant fortement de la teneur en eau, il est alors possible de déterminer celle-ci. La figure 1 représente un schéma de principe de dispositif apte à effectuer une mesure TDR. 2 Le dispositif comprend un réflectomètre 1, une antenne ou une sonde S placée dans le milieu M que l'on étudie, un système de transport de signal 2 qui relie le réflectomètre 1 à la sonde S et un dispositif de traitement de signal 3 relié au réflectomètre 1. La sonde S est faite, par exemple, de deux électrodes métalliques a1, a2 qui se terminent en circuit ouvert, l'une des deux électrodes métalliques étant reliée à la masse du dispositif. Le dispositif de traitement 3 délivre une commande Cd qui est appliquée au réflectomètre 1. La commande Cd permet, tout d'abord, de lancer l'acquisition du signal de mesure et, ensuite, de transférer un signal numérique R(t) qui résulte de cette acquisition vers le circuit de traitement 3. Pour l'acquisition du signal de mesure, le réflectomètre 1 émet une impulsion qui se propage dans le système de transport 2. Lorsque l'impulsion arrive sur le connecteur d'entrée CN qui relie le système de transport 2 à la sonde S, une première partie de l'onde est réfléchie et une deuxième partie est transmise le long des électrodes a1, a2. Lorsque la deuxième partie de l'onde arrive sur le circuit ouvert qui termine la sonde, cette deuxième partie de l'onde est totalement réfléchie. Le signal de mesure ainsi constitué est alors numérisé et transféré au dispositif de traitement sous l'action de la commande Cd. La figure 2 illustre le principe de calcul du temps de propagation de l'onde électromagnétique qui parcourt l'antenne, selon l'art antérieur. Sur la figure 2 sont représentées deux courbes Cl, C2 du signal 3 numérique R(t) mesuré par le réflectomètre 1, dans le cas d'une sonde placée dans de l'argile. Les courbes Cl et C2 ont des formes différentes. Cela est dû au fait que les câbles de liaison qui sont utilisés pour l'obtention de ces courbes ont des longueurs différentes : la longueur des câbles du système de transport 2 est plus grande dans le cas de la courbe C2 que dans le cas de la courbe Cl. Chaque courbe Cl, C2 comprend un premier maximum (respectivement maxi, max2) qui représente la première réflexion, une décroissance du signal vers un minimum (respectivement miner mine) et, à partir du minimum, une forte remontée du signal qui présente un point d'inflexion (respectivement I2r I2) et représente la deuxième réflexion. Selon l'art antérieur, l'instant qui représente la première réflexion est déterminé des deux tangentes à la courbe premier maximum et l'instant deuxième réflexion est déterminé par l'intersection encadrent le représente la l'intersection minimum et de par qui qui de la tangente à la courbe au niveau du la tangente à la courbe au niveau du point d'inflexion de la forte remontée du signal. La durée At qui correspond au temps de réflexion aller-retour de l'onde qui se propage le long de l'antenne est alors donnée par la différence entre l'instant qui représente la deuxième réflexion et l'instant qui représente la première réflexion. En référence à la courbe Cl, le temps de propagation de l'onde électromagnétique s'écrit en conséquence : 4 Ott = t22 - t11 De même, en référence à la courbe C2, il vient : At2 = t12 - t21 Il apparaît ainsi clairement que le temps de propagation Ott peut être très différent du temps de propagation At2. A titre d'exemple non limitatif, pour une sonde de 60cm, la durée At2 est égale à 1Ons pour un câble de longueur 40m et la durée Ott est égale à 12ns pour un câble de longueur 20m. TECHNICAL FIELD AND PRIOR ART The invention relates to a method for measuring the water content of a medium and a device for measuring the water content which implements the method of the invention. The invention is advantageously applied to the measurement of the water content in natural materials (soils, sands, clays, argillites, etc.), cementitious materials (cements, mortars, concretes) or remodeled. The method of the invention is situated in the general context of measurements carried out according to the time domain reflectometry principle more commonly known as TDR (Time Domain Reflectometry) measurements. A TDR measurement consists of measuring the time that an electromagnetic wave takes to propagate along electrodes placed in a medium. The knowledge of this time makes it possible to determine the speed of propagation of the wave in the medium and, consequently, the relative permittivity of the medium. In the field of very high frequencies, the relative permittivity of the medium strongly depending on the water content, it is then possible to determine it. Figure 1 shows a block diagram of a device capable of performing a TDR measurement. 2 The device comprises a reflectometer 1, an antenna or a probe S placed in the medium M that is being studied, a signal transport system 2 which connects the reflectometer 1 to the probe S and a connected signal processing device 3 The sensor S is made, for example, of two metal electrodes a1, a2 which terminate in an open circuit, one of the two metal electrodes being connected to the ground of the device. The processing device 3 delivers a command Cd which is applied to the reflectometer 1. The command Cd makes it possible, first of all, to start the acquisition of the measurement signal and then to transfer a digital signal R (t) that results from this acquisition to the processing circuit 3. For the acquisition of the measurement signal, the reflectometer 1 emits a pulse which propagates in the transport system 2. When the pulse arrives on the CN input connector which connects the transport system 2 to the probe S, a first part of the wave is reflected and a second part is transmitted along the electrodes a1, a2. When the second part of the wave arrives on the open circuit that ends the probe, this second part of the wave is totally reflected. The measurement signal thus constituted is then digitized and transferred to the processing device under the action of the control Cd. FIG. 2 illustrates the principle of calculating the propagation time of the electromagnetic wave which travels over the antenna, according to FIG. prior art. In FIG. 2 are represented two curves C1, C2 of the digital signal R (t) measured by the reflectometer 1, in the case of a probe placed in clay. Curves C1 and C2 have different shapes. This is because the connecting cables which are used to obtain these curves have different lengths: the length of the cables of the transport system 2 is greater in the case of the curve C2 than in the case of the curve Cl. Each curve C1, C2 comprises a first maximum (respectively maxi, max2) which represents the first reflection, a decrease of the signal towards a minimum (respectively miner mine) and, from the minimum, a strong rise of the signal which presents a point of inflection (respectively I2r I2) and represents the second reflection. According to the prior art, the instant which represents the first reflection is determined from the two tangents to the maximum first curve and the second reflection moment is determined by the intersection frame represents the minimum intersection and by whom which of the tangent to the curve at the tangent to the curve at the point of inflection of the strong rise of the signal. The duration At which corresponds to the round-trip reflection time of the wave propagating along the antenna is then given by the difference between the instant which represents the second reflection and the moment which represents the first reflection. With reference to curve C1, the propagation time of the electromagnetic wave is written accordingly: 4 Ott = t22 - t11 Similarly, with reference to curve C2, it comes: At2 = t12 - t21 It thus appears clearly that the Ott propagation time can be very different from the At2 propagation time. As a non-limiting example, for a 60cm probe, the duration At2 is equal to 1Ons for a cable length 40m and the duration Ott is equal to 12ns for a cable length 20m.

Cette forte dépendance de la mesure du temps de propagation à la longueur des câbles de liaison est un réel inconvénient. Il en résulte en effet une grande incertitude dans les résultats de calcul de la teneur en eau. C'est tout particulièrement l'instant qui représente la première réflexion qui est défini avec beaucoup d'imprécision. Ceci représente un inconvénient majeur pour déterminer, de façon reproductible, la constante diélectrique relative d'un milieu et, partant, la teneur en eau de ce milieu. This strong dependence of the measurement of the propagation time on the length of the connecting cables is a real disadvantage. This results in a great deal of uncertainty in the results of calculating the water content. It is especially the moment which represents the first reflection which is defined with a lot of inaccuracy. This represents a major disadvantage to reproducibly determine the relative dielectric constant of a medium and hence the water content of this medium.

L'invention ne présente pas cet inconvénient. Exposé de l'invention En effet, l'invention concerne un procédé de mesure de teneur en eau d'un milieu par 25 réflectométrie, le procédé comprenant : - la formation d'un signal numérique représentatif d'un signal réfléchi par une antenne placée dans le milieu, le signal numérique étant constitué d'une succession d'échantillons numériques 30 dont la variation en fonction du temps constitue une courbe, - la détermination, à partir du signal numérique, d'un premier instant et d'un deuxième instant ultérieur au premier instant, la durée qui sépare le deuxième instant du premier instant mesurant 5 un temps de propagation aller-retour d'un signal qui se propage le long de l'antenne, - la détermination, à partir des échantillons numériques qui précèdent le premier instant, d'une valeur moyenne sensiblement constante qui représente un niveau de bruit du signal numérique, et - la détermination du premier instant comme un instant où un nombre N d'échantillons numériques successifs s'élèvent continument au dessus de la valeur moyenne sensiblement constante, N étant un nombre entier supérieur ou égal à 2, conformément à la théorie de Nyquist sur la fréquence d'échantillonnage d'un signal. Selon une caractéristique supplémentaire du procédé de l'invention, le deuxième instant est calculé, par un circuit de traitement, comme un instant où une tangente à la courbe en un point de minimum de la courbe et une tangente à la courbe en un point d'inflexion d'une partie de la courbe qui s'élève continument à partir du point de minimum s'interceptent. Selon une autre caractéristique supplémentaire de l'invention, le point d'inflexion est déterminé par les étapes suivantes : 6 The invention does not have this disadvantage. The invention relates to a method for measuring the water content of a medium by reflectometry, the method comprising: the formation of a digital signal representative of a signal reflected by a placed antenna in the medium, the digital signal consisting of a succession of digital samples 30 whose variation as a function of time constitutes a curve, - the determination, from the digital signal, of a first instant and a second instant subsequent to the first instant, the time which separates the second instant from the first instant measuring a round trip delay of a signal propagating along the antenna, - the determination, from the digital samples which precede the first instant, a substantially constant average value which represents a noise level of the digital signal, and - the determination of the first instant as a moment when a number N of samples successive numerals rise continuously above the substantially constant average value, N being an integer greater than or equal to 2, in accordance with Nyquist's theory of the sampling frequency of a signal. According to a further characteristic of the method of the invention, the second instant is calculated, by a processing circuit, as a moment when a tangent to the curve at a minimum point of the curve and a tangent to the curve at a point d inflection of a part of the curve which rises continuously from the point of minimum are interposed. According to another additional characteristic of the invention, the point of inflection is determined by the following steps:

- filtrage, à partir du point de minimum de la courbe, de la partie de la courbe qui s'élève à partir du point de minimum, - détermination d'une courbe croissante au sens strict à partir de la partie de courbe filtrée obtenue à l'étape de filtrage, - détermination du point d'inflexion comme étant le point d'inflexion de la courbe croissante au sens strict. filtering, from the minimum point of the curve, the part of the curve that rises from the minimum point, determining a strictly increasing curve from the filtered curve portion obtained at the filtering step, - determining the point of inflection as being the point of inflection of the increasing curve in the strict sense.

Selon encore une autre caractéristique supplémentaire de l'invention, après la détermination de la courbe croissante au sens strict et avant la détermination du point d'inflexion, une étape d'interpolation par spline cubique élimine des points aberrants de la courbe. Selon encore une autre caractéristique supplémentaire de l'invention, le filtrage est effectué par érosion morphologique. L'invention concerne également un dispositif de mesure de teneur en eau d'un milieu par réflectométrie, le dispositif comprenant : - des moyens de formation d'un signal numérique représentatif d'un signal analogique réfléchi par une antenne placée dans le milieu, le signal numérique étant constitué d'une succession d'échantillons numériques dont la variation en fonction du temps constitue une courbe, et - des moyens de détermination, à partir du signal numérique, d'un premier instant et d'un deuxième instant ultérieur au premier instant, la durée qui sépare le deuxième instant du premier instant mesurant 7 un temps de propagation aller-retour d'un signal qui se propage le long de l'antenne. Les moyens de détermination du premier instant comprennent une unité de calcul apte à déterminer . - à partir des échantillons numériques qui précèdent le premier instant, une valeur moyenne sensiblement constante qui représente un niveau de bruit du signal numérique, et - le premier instant comme un instant où un nombre N d'échantillons numériques successifs s'élèvent continument au dessus de la valeur moyenne sensiblement constante, N étant un nombre entier supérieur ou égal à 2. According to yet another additional characteristic of the invention, after the determination of the curve increasing in the strict sense and before the determination of the point of inflection, a cubic spline interpolation step eliminates outliers of the curve. According to yet another additional characteristic of the invention, the filtering is carried out by morphological erosion. The invention also relates to a device for measuring the water content of a medium by reflectometry, the device comprising: - means for forming a digital signal representative of an analog signal reflected by an antenna placed in the medium, the digital signal consisting of a succession of digital samples whose variation as a function of time constitutes a curve, and - means for determining, from the digital signal, a first instant and a second instant subsequent to the first moment, the time between the second instant and the first instant measuring a round trip delay of a signal propagating along the antenna. The means for determining the first instant comprise a calculation unit that can be determined. from the digital samples which precede the first instant, a substantially constant mean value which represents a noise level of the digital signal, and the first moment as a moment when a number N of successive digital samples rise continuously above of the substantially constant average value, N being an integer greater than or equal to 2.

De façon particulièrement avantageuse, le procédé de l'invention permet d'obtenir des mesures de teneur en eau beaucoup plus reproductibles que le procédé de l'art antérieur lorsque les conditions de mesure changent (modification de longueur des câbles de liaison). Brève description des figures D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel fait en référence aux figures jointes, parmi lesquelles : - La figure 1, déjà décrite, représente un schéma de principe de dispositif apte à effectuer une mesure TDR ; - La figure 2, déjà décrite, illustre le principe de calcul du temps de propagation de l'onde 8 Particularly advantageously, the method of the invention makes it possible to obtain much more reproducible water content measurements than the method of the prior art when the measurement conditions change (modification of the length of the connecting cables). BRIEF DESCRIPTION OF THE FIGURES Other features and advantages of the invention will appear on reading a preferred embodiment with reference to the appended figures, among which: FIG. 1, already described, represents a device block diagram able to perform a TDR measurement; FIG. 2, already described, illustrates the principle of calculating the propagation time of the wave 8

électromagnétique qui parcourt les antennes, selon l'art antérieur ; - La figure 3 représente un schéma de principe de dispositif apte à effectuer une mesure TDR selon l'invention ; - La figure 4 représente une unité de calcul apte à mettre en oeuvre le procédé de l'invention ; - La figure 5 représente une vue de détail de l'unité de calcul de la figure 4 ; - La figure 6 illustre le principe de calcul du temps de propagation de l'onde électromagnétique qui parcourt les antennes, selon l'invention. Sur toutes les figures, les mêmes références désignent les mêmes éléments. electromagnetic which travels the antennas, according to the prior art; FIG. 3 represents a block diagram of a device able to perform a TDR measurement according to the invention; FIG. 4 represents a calculation unit able to implement the method of the invention; FIG. 5 represents a detailed view of the computing unit of FIG. 4; FIG. 6 illustrates the principle of calculating the propagation time of the electromagnetic wave traveling through the antennas, according to the invention. In all the figures, the same references designate the same elements.

Exposé détaillé d'un mode de réalisation préférentiel de l'invention. La figure 3 représente un schéma de principe de dispositif apte à effectuer une mesure TDR selon l'invention. Le dispositif de mesure de l'invention diffère du dispositif de l'art antérieur par le dispositif de traitement UT. Le dispositif de traitement UT reçoit les données numériques qui constituent le signal R(t). Le dispositif UT comprend (cf. figures 4 et 5) une unité de calcul U1 qui calcule l'instant de départ tdi de la durée At, une unité de calcul U2 qui calcule l'instant de fin tfi de la durée At et une unité de calcul U3 qui calcule la durée At par différence entre l'instant de fin et l'instant de départ. L'unité de traitement UT est, par exemple, un microcontrôleur, un ordinateur, 9 une carte de traitement, etc. Le dispositif UT comprend également des moyens (non représentés sur la figure 4) pour délivrer la commande Cd qui pilote le réflectomètre. DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT OF THE INVENTION FIG. 3 represents a block diagram of a device able to perform a TDR measurement according to the invention. The measuring device of the invention differs from the device of the prior art by the processing device UT. The processing device UT receives the digital data constituting the signal R (t). The device UT comprises (see FIGS. 4 and 5) a calculation unit U1 which calculates the start time tdi of the duration At, a calculation unit U2 which calculates the end time tfi of the duration At and a unit calculation U3 which calculates the duration At by difference between the end time and the start time. The processing unit UT is, for example, a microcontroller, a computer, a processing card, etc. The device UT also comprises means (not shown in FIG. 4) for delivering the command Cd which drives the reflectometer.

La figure 6 illustre le principe de calcul du temps de propagation de l'onde qui parcourt l'antenne selon l'invention. Sur la figure 6, les courbes C1 et C2 sont identiques aux courbes représentées sur la figure 2 et représentent, en conséquence, les mêmes mesures. Selon l'invention, pour chaque courbe Ci (i=1, 2), le point de départ de la durée At de propagation de l'onde est calculé en un point Pi(tdi, Xi) où l'instant tdi est l'instant où la courbe s'élève d'un nombre N de points successifs (par exemple N=5) pour prendre une valeur Xi supérieure à une valeur moyenne Xoi qui représente la partie horizontale de la courbe qui précède l'arrivée de la première réflexion (niveau de bruit). FIG. 6 illustrates the principle of calculating the propagation time of the wave traveling through the antenna according to the invention. In FIG. 6, the curves C1 and C2 are identical to the curves represented in FIG. 2 and represent, consequently, the same measurements. According to the invention, for each curve Ci (i = 1, 2), the starting point of the propagation time At of the wave is calculated at a point Pi (tdi, Xi) where the instant tdi is the moment when the curve rises by a number N of successive points (for example N = 5) to take a value Xi greater than an average value Xoi which represents the horizontal part of the curve which precedes the arrival of the first reflection (noise level).

L'unité de calcul U1 calcule ainsi, tout d'abord, la valeur moyenne Xoi à partir des données relatives à la courbe Ci. La valeur moyenne Xoi est donnée par la tangente à la courbe Ci, calculée en un point quelconque de la partie horizontale de la courbe Ci qui représente le niveau de bruit du signal. L'instant tdi est alors calculé comme l'instant auquel le signal numérique s'élève de N échantillons numériques au dessus de la tangente de valeur moyenne . A titre d'exemple non limitatif, N est égal à 5. The calculation unit U1 thus calculates, first of all, the average value Xoi from the data relating to the curve Ci. The mean value Xoi is given by the tangent to the curve Ci, calculated at any point of the part horizontal curve Ci which represents the noise level of the signal. The instant tdi is then calculated as the instant at which the digital signal rises from N digital samples above the mean value tangent. By way of non-limiting example, N is equal to 5.

L'instant tfi de fin de la durée At est calculé, par exemple, par la méthode des tangentes 10 (unité de calcul U2). L'unité de calcul U3 calcule ensuite la durée At comme la différence des instants tfi et tai Un mode de réalisation préférentiel de l'unité de calcul U2 est détaillé en référence à la figure 5. L'unité de calcul U2 comprend les unités de calcul élémentaires suivantes : une unité de filtrage Ua qui filtre, à partir du point de minimum de la courbe, les points situés sur la partie de la courbe qui présentent une forte remontée du signal ; - une unité de détermination Ub d'une courbe croissante au sens strict à partir de la courbe filtrée délivrée par l'unité de filtrage; - une unité de détermination Uc du point d'inflexion situé sur la courbe croissante au sens strict obtenue précédemment ; et - une unité de détermination Ud de l'instant tfi à l'intersection de la tangente à la courbe au niveau du minimum et de la tangente à la courbe au point d'inflexion. Dans un mode de réalisation particulier de l'invention, l'unité de filtrage Ua met en oeuvre la méthode de traitement communément appelée « érosion morphologique » qui est connue, à ce jour, dans le domaine du traitement d'image. L'érosion morphologique, appliquée de façon avantageuse aux échantillons numériques du signal, permet ici d'obtenir une courbe dépourvue de toute partie décroissante, la pente de la courbe obtenue étant strictement positive ou nulle. Dans un perfectionnement de l'invention, après 11 détermination d'une courbe croissante au sens strict et avant détermination du point d'inflexion, le procédé comprend une étape d'interpolation par spline cubique qui élimine les points aberrants de la courbe.5 The time tfi end of the duration At is calculated, for example, by the method of tangents 10 (computing unit U2). The calculation unit U3 then calculates the duration At as the difference of the instants tfi and tai. A preferred embodiment of the calculation unit U2 is detailed with reference to FIG. 5. The calculation unit U2 comprises the following elementary calculations: a filter unit Ua which filters, from the minimum point of the curve, the points situated on the part of the curve which show a strong rise of the signal; a determination unit Ub of a curve increasing in the strict sense from the filtered curve delivered by the filtering unit; a determination unit Uc of the inflection point situated on the strictly increasing curve obtained previously; and a determination unit Ud of the instant tfi at the intersection of the tangent to the curve at the minimum and the tangent to the curve at the point of inflection. In a particular embodiment of the invention, the filter unit Ua implements the treatment method commonly known as "morphological erosion" which is known, to date, in the field of image processing. Morphological erosion, advantageously applied to the digital samples of the signal, makes it possible here to obtain a curve devoid of any decreasing part, the slope of the curve obtained being strictly positive or zero. In an improvement of the invention, after determining a strictly increasing curve and before determining the point of inflection, the method comprises a cubic spline interpolation step which eliminates the outliers of the curve.

Claims (10)

REVENDICATIONS1. Procédé de mesure de teneur en eau d'un milieu (M) par réflectométrie, le procédé comprenant : - la formation d'un signal numérique (R (t) ) représentatif d'un signal réfléchi par une antenne placée dans le milieu (M), le signal numérique étant constitué d'une succession d'échantillons numériques dont la variation en fonction du temps constitue une courbe (Ci), - la détermination, à partir du signal numérique, d'un premier instant (tdi) et d'un deuxième instant (tfi) ultérieur au premier instant, la durée (At) qui sépare le deuxième instant du premier instant mesurant un temps de propagation aller-retour d'un signal qui se propage le long de l'antenne, caractérisé en ce que la détermination du premier instant comprend : - une étape de détermination, à partir des échantillons numériques qui précèdent le premier instant (tdi), d'une valeur moyenne sensiblement constante (Xoi) qui représente un niveau de bruit du signal numérique, et - une étape de détermination du premier instant (tdi) comme un instant où un nombre N d'échantillons numériques successifs s'élèvent continument au dessus de la valeur moyenne sensiblement constante, N étant un nombre entier supérieur ou égal à2.30 13 REVENDICATIONS1. A method for measuring the water content of a medium (M) by reflectometry, the method comprising: - forming a digital signal (R (t)) representative of a signal reflected by an antenna placed in the medium (M ), the digital signal consisting of a succession of digital samples whose variation as a function of time constitutes a curve (Ci), - the determination, from the digital signal, of a first instant (tdi) and of a second instant (tfi) subsequent to the first instant, the duration (At) which separates the second instant from the first instant measuring a round trip delay of a signal that propagates along the antenna, characterized in that the determination of the first instant comprises: - a step of determining, from the digital samples which precede the first instant (tdi), a substantially constant mean value (Xoi) which represents a noise level of the digital signal, and - a dice step terminating the first instant (tdi) as a time when a number N of successive digital samples rise continuously above the substantially constant mean value, N being an integer greater than or equal to 2.30 2. Procédé selon la revendication 1, dans lequel le deuxième instant (tdi) est calculé, par un circuit de traitement (U3), comme un instant où une tangente à la courbe (Ci) en un point de minimum de la courbe (mini) et une tangente à la courbe en un point d'inflexion (Pi) d'une partie de la courbe qui s'élève continument à partir du point de minimum s'interceptent. 2. Method according to claim 1, wherein the second instant (tdi) is calculated, by a processing circuit (U3), as a moment when a tangent to the curve (Ci) at a minimum point of the curve (min. ) and a tangent to the curve at a point of inflection (Pi) of a portion of the curve that rises continuously from the point of minimum intersect. 3. Procédé selon la revendication 2, dans lequel le point d'inflexion est déterminé par les étapes suivantes . - filtrage (Ua), à partir du point de minimum de la courbe, de la partie de la courbe qui s'élève à partir du point de minimum, détermination (Ub) d'une courbe croissante au sens strict à partir de la partie de courbe filtrée obtenue à l'étape de filtrage, - détermination (Uc) du point d'inflexion comme étant le point d'inflexion de la courbe croissante au sens strict. 3. Method according to claim 2, wherein the point of inflection is determined by the following steps. filtering (Ua), starting from the minimum point of the curve, the part of the curve that rises from the point of minimum, determining (Ub) of a curve increasing in the strict sense from the part filtered curve obtained in the filtering step, - determining (Uc) of the inflection point as the point of inflection of the curve increasing in the strict sense. 4. Procédé selon la revendication 3, dans lequel, après la détermination de la courbe croissante au sens strict et avant la détermination du point d'inflexion, une étape d'interpolation par spline cubique est effectuée pour éliminer des points aberrants de la courbe. 14 4. The method of claim 3, wherein, after determining the curve increasing in the strict sense and before the determination of the point of inflection, a cubic spline interpolation step is performed to remove outliers of the curve. 14 5. Procédé selon la revendication 3 ou 4, dans lequel le filtrage est effectué par érosion morphologique. 5. The method of claim 3 or 4, wherein the filtering is performed by morphological erosion. 6. Dispositif de mesure de teneur en eau d'un milieu (M) par réflectométrie, le dispositif comprenant : - des moyens de formation d'un signal numérique (R(t)) représentatif d'un signal réfléchi par une antenne placée dans le milieu (M), le signal numérique étant constitué d'une succession d'échantillons numériques dont la variation en fonction du temps constitue une courbe (Ci), - des moyens de détermination, à partir du signal numérique, d'un premier instant (tdi) et d'un deuxième instant (tfi) ultérieur au premier instant, la durée (At) qui sépare le deuxième instant du premier instant mesurant un temps de propagation aller-retour d'un signal qui se propage le long de l'antenne, caractérisé en ce que les moyens de détermination du premier instant comprennent une unité de calcul (U2) apte à déterminer : - à partir des échantillons numériques qui précèdent le premier instant, une valeur moyenne sensiblement constante (Xoi) qui représente un niveau de bruit du signal numérique, et - le premier instant (tdi) comme un instant où un nombre N d'échantillons numériques successifs s'élèvent continument au dessus de la valeur moyenne sensiblement constante, N étant un nombre entier supérieur ou égal à 2. 15 6. Device for measuring the water content of a medium (M) by reflectometry, the device comprising: - means for forming a digital signal (R (t)) representative of a signal reflected by an antenna placed in a the medium (M), the digital signal consisting of a succession of digital samples whose variation as a function of time constitutes a curve (Ci), - means for determining, from the digital signal, a first moment (tdi) and a second time (tfi) subsequent to the first instant, the duration (At) between the second instant of the first instant measuring a round trip delay of a signal that propagates along the antenna, characterized in that the means for determining the first instant comprise a calculation unit (U2) able to determine: - from the digital samples which precede the first instant, a substantially constant average value (Xoi) which represents a level of bru it of the digital signal, and the first instant (tdi) as a moment when a number N of successive digital samples rise continuously above the substantially constant average value, N being an integer greater than or equal to 2. 15 7. Dispositif selon la revendication 6 et qui comprend, en outre, un circuit de traitement (U3) apte à déterminer le deuxième instant (tfi) comme un instant où une tangente à la courbe (Ci) en un point de minimum de la courbe et une tangente à la courbe en un point d'inflexion d'une partie de la courbe qui s'élève continument à partir du point de minimum s'interceptent. 7. Device according to claim 6 and which further comprises a processing circuit (U3) capable of determining the second instant (tfi) as a moment when a tangent to the curve (Ci) at a minimum point of the curve. and a tangent to the curve at a point of inflection of a portion of the curve that rises continuously from the point of minimum intersect. 8. Dispositif selon la revendication 7, dans lequel le circuit de traitement comprend : - une unité de filtrage (Ua) qui filtre, à partir du point de minimum de la courbe, la partie de la courbe qui s'élève à partir du point de minimum, une unité de détermination (Ub) d'une courbe croissante au sens strict à partir de la partie de courbe filtrée délivrée par l'unité de filtrage, et une unité de détermination (Uc) du point d'inflexion comme étant le point d'inflexion de la courbe croissante au sens strict. 8. Device according to claim 7, wherein the processing circuit comprises: - a filter unit (Ua) which filters, from the point of minimum curve, the part of the curve which rises from the point at least a determination unit (Ub) of a strictly increasing curve from the filtered curve portion delivered by the filter unit, and a determination unit (Uc) of the inflection point as the point of inflection of the increasing curve in the strict sense. 9. Dispositif selon la revendication 8, dans lequel une unité de calcul met en oeuvre une étape d'interpolation par spline cubique qui élimine des points aberrants de la courbe. 9. Device according to claim 8, wherein a calculation unit implements a cubic spline interpolation step which eliminates aberrant points of the curve. 10. Dispositif selon la revendication 8 ou 9, dans lequel l'unité de filtrage effectue un filtrage par érosion morphologique. The device of claim 8 or 9, wherein the filtering unit performs morphological erosion filtering.
FR1052094A 2010-03-23 2010-03-23 METHOD FOR MEASURING THE WATER CONTENT OF A MEDIUM AND ASSOCIATED DEVICE Active FR2958041B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
FR1052094A FR2958041B1 (en) 2010-03-23 2010-03-23 METHOD FOR MEASURING THE WATER CONTENT OF A MEDIUM AND ASSOCIATED DEVICE

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR1052094A FR2958041B1 (en) 2010-03-23 2010-03-23 METHOD FOR MEASURING THE WATER CONTENT OF A MEDIUM AND ASSOCIATED DEVICE

Publications (2)

Publication Number Publication Date
FR2958041A1 true FR2958041A1 (en) 2011-09-30
FR2958041B1 FR2958041B1 (en) 2013-03-15

Family

ID=42935443

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1052094A Active FR2958041B1 (en) 2010-03-23 2010-03-23 METHOD FOR MEASURING THE WATER CONTENT OF A MEDIUM AND ASSOCIATED DEVICE

Country Status (1)

Country Link
FR (1) FR2958041B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9784093B2 (en) 2014-05-08 2017-10-10 WellGauge, Inc. Well water depth monitor
FR3112613A1 (en) * 2020-07-20 2022-01-21 Universite Grenoble Alpes INSTRUMENT AND METHOD FOR ANALYZING A COMPLEX MEDIUM TO DETERMINE ITS PHYSICO-CHEMICAL PROPERTIES

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030042916A1 (en) * 2001-09-04 2003-03-06 Anderson Scott Knudson Absolute-reading soil moisture and conductivity sensor

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030042916A1 (en) * 2001-09-04 2003-03-06 Anderson Scott Knudson Absolute-reading soil moisture and conductivity sensor

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
ANONYMOUS: "INSTRUCTION MANUAL TDR100, Revision 2/10", 17 February 2010 (2010-02-17), pages I-III;28 - 32, XP002606967, Retrieved from the Internet <URL:http://www.campbellsci.com/documents/manuals/tdr100.pdf> [retrieved on 20101018] *
NEZAFAT M ET AL: "A new enhanced morphological filter and signal recovery", EUROCON'2001, TRENDS IN COMMUNICATIONS, INTERNATIONAL CONFERENCE ON. JUL. 4-7, 2001, PISCATAWAY, NJ, USA,IEEE, vol. 1, 4 July 2001 (2001-07-04), pages 91 - 94, XP010554275, ISBN: 978-0-7803-6490-5 *
ROBINSON D A ET AL: "Considerations for improving the accuracy of permittivity measurement using time domain reflectometry: Air-water calibration, effects of cable length", SOIL SCIENCE SOCIETY OF AMERICA JOURNAL JANUARY/FEBRUARY 2003 SOIL SCIENCE SOCIETY OF AMERICA US, vol. 67, no. 1, January 2003 (2003-01-01), pages 62 - 70, XP002606968 *
SOUTO F J ET AL: "Design and air-water calibration of a waveguide connector for TDR measurements of soil electric permittivity in stony soils", BIOSYSTEMS ENGINEERING, vol. 101, no. 4, 1 December 2008 (2008-12-01), ACADEMIC PRESS, UK, pages 463 - 471, XP025691779, ISSN: 1537-5110, [retrieved on 20081108], DOI: 10.1016/J.BIOSYSTEMSENG.2008.09.015 *
STEVEN R. EVETT: "THE TACQ COMPUTER PROGRAM FOR AUTOMATIC TIME DOMAIN REFLECTOMETRY MEASUREMENTS:II. WAVEFORM INTERPRETATION METHODS1", TRANSACTIONS OF ASAE, vol. 1, no. 6, 2000, pages 1 - 12, XP002606966, Retrieved from the Internet <URL:http://ddr.nal.usda.gov/bitstream/10113/1943/1/IND23220843.pdf> [retrieved on 20101018] *
TIMLIN D J ET AL: "Comparison of three methods to obtain the apparent dielectric constant from time domain reflectometry wave traces", SOIL SCIENCE SOCIETY OF AMERICA JOURNAL 1996 JUL-AUG SOIL SCIENCE SOC OF AMERICA, vol. 60, no. 4, July 1996 (1996-07-01), pages 970 - 977, XP002606969 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9784093B2 (en) 2014-05-08 2017-10-10 WellGauge, Inc. Well water depth monitor
US10145235B2 (en) 2014-05-08 2018-12-04 WellGauge, Inc. Well water depth monitor
FR3112613A1 (en) * 2020-07-20 2022-01-21 Universite Grenoble Alpes INSTRUMENT AND METHOD FOR ANALYZING A COMPLEX MEDIUM TO DETERMINE ITS PHYSICO-CHEMICAL PROPERTIES
WO2022018377A1 (en) * 2020-07-20 2022-01-27 Universite Grenoble Alpes Instrument and method for analysing a complex medium in order to determine its physicochemical properties

Also Published As

Publication number Publication date
FR2958041B1 (en) 2013-03-15

Similar Documents

Publication Publication Date Title
EP1097354B1 (en) Crossed measurements of flowmeter sound signals
FR2632734A1 (en) METHOD OF ESTABLISHING A STRATIGRAPHIC MODEL OF THE BASEMENT FROM AN ACOUSTIC IMPEDANCE PROFILE AND A SEISMIC SECTION
WO2019091705A1 (en) Structural health monitoring for an industrial structure
EP1615049B1 (en) Fast coherent treatment for codes with periodical line spectra
EP0247908A1 (en) Liquid level pick-up head using elastically guided waves
EP1751494B1 (en) Device and method for the dimensional characterisation of a cylindrical object
EP0195792A1 (en) Method and device for detecting by using foucault currents typical defects on a product.
FR2687227A1 (en) METHOD FOR DETERMINING A FIELD - SIGNATURE SIGNATURE OF A PLURALITY OF SEISMIC SOURCES.
EP0550333A1 (en) Method and apparatus for measuring gas flow rate
FR2958041A1 (en) Method for time domain reflectometry type measuring of water content in natural material e.g. soil, involves determining primary instant as instant in which number of digital samples is increased continuously above average value
EP0490740A1 (en) Method and apparatus for pitch period determination of the speech signal in very low bitrate vocoders
EP2849994B1 (en) Method for predicting at least one movement of a ship under the effect of the waves
FR2930344A1 (en) METHOD AND DEVICE FOR NON-DESTRUCTIVE ULTRASONIC CONTROL OF THE POROSITY OF A COMPOSITE MATERIAL PART
CA2278408C (en) Procedure and device for measuring the flow rate of a liquid stream
EP2941653A1 (en) Method for analysing a cable by compensating for the dispersion experienced by a signal when it is propagated within said cable
EP1360524B1 (en) Method for determining ocean current and associated device
Wong et al. One-dimensional signal processing techniques for airborne laser bathymetry
EP4155724A1 (en) Method for detecting a fault in a structure of a device by vibro-acoustic modulation
FR3086064A1 (en) METHOD FOR DETERMINING A DEPTH, OR A BATHYMMETRIC PROFILE, BASED ON AN AVERAGE CELERY PROFILE OF SOUND, METHOD FOR DETERMINING SUCH A CELERY PROFILE, AND ASSOCIATED SONAR SYSTEM
CA1270937A (en) Method and device for determining the long range emission signature of a set of seismic emissions
CA2134856A1 (en) Process for determining the distance ambiguity rank in radar echoes
FR2512860A1 (en) Digital surface controller for ground probing tests - contains microprocessor performing calculations and issuing commands resulting from pressure and water level sensor inputs
EP3213112B1 (en) Method for processing a radar signal in land/sea detection mode; processing system and associated computer program product
EP3646045B1 (en) Reflectometry system for analysing faults in a transmission line
FR2626666A1 (en) Method for measuring the thickness of road covering layers using pulsed radar

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 7

PLFP Fee payment

Year of fee payment: 8

PLFP Fee payment

Year of fee payment: 9

PLFP Fee payment

Year of fee payment: 11

PLFP Fee payment

Year of fee payment: 12

PLFP Fee payment

Year of fee payment: 13

PLFP Fee payment

Year of fee payment: 14

PLFP Fee payment

Year of fee payment: 15