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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 14
- 238000002310 reflectometry Methods 0.000 title claims description 7
- 239000002689 soil Substances 0.000 title description 2
- 239000005445 natural material Substances 0.000 title 1
- 238000012545 processing Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 9
- 230000003628 erosive effect Effects 0.000 claims description 5
- 230000000877 morphologic effect Effects 0.000 claims description 5
- 230000001594 aberrant effect Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 15
- 239000000523 sample Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000001902 propagating effect Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000004567 concrete Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000004570 mortar (masonry) Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/246—Earth materials for water content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
- G01N22/04—Investigating 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
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)
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)
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)
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 |
-
2010
- 2010-03-23 FR FR1052094A patent/FR2958041B1/en active Active
Patent Citations (1)
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)
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)
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 |