FR 3112613 A1 20220121 FR 2007582 A 20200720 La présente invention a pour objet un instrument et un procédé de détermination du profil basé sur la mesure, à différentes fréquences, de la répartition, selon l’axe de propagation, de l’amplitude d’une onde stationnaire de champ électrique créée dans une ligne de transmission insérée dans le milieu stratifié à étudier. Figure pour l’abrégé : figure 1. INSTRUMENT ET PROCEDE D’ANALYSE D’UN MILIEU COMPLEXE POUR EN DETERMINER SES PROPRIETES PHYSICO-CHIMIQUES La présente invention concerne de manière générale un instrument de type sonde et un procédé pour l’analyse en profondeur d’un milieu complexe pour en déterminer les propriétés physico-chimiques. Par milieu complexe, on entend, au sens de la présente invention, des milieux tels que des sols, notamment agricoles, des murs, ou le manteau neigeux. Il s’agit essentiellement de milieux naturels qui peuvent être stratifiés. Pour l’analyse en profondeur de tels milieux, l’homme de l’art connait des méthodes classiques telles que la mesure des propriétés diélectriques du milieu à analyser en fonction de la fréquence du générateur sinusoïdal. Elle permet de détecter finement la nature des matériaux présents à l’endroit même de la mesure. L’homme de l’art connait également la méthode de réflectométrie temporelle (usuellement désignée par l’acronyme en anglais TDR pour « Time Do main Reflectometry »). Cette méthode consiste à envoyer un échelon de tension dans le milieu à analyser et à enregistrer la tension réfléchie. Elle permet de détecter la variation de la permittivité électrique du milieu en fonction de la profondeur (changements de milieux, couches). Ces méthodes classiques présentent toutefois l’inconvénient de ne pas pouvoir analyser complètement en profondeur un milieu complexe : elles ne permettent d’analyser qu’un seul paramètre tel que la capacité électrique pour la première, à l’endroit même de la mesure, tandis que dans la méthode TDR, une analyse en profondeur est possible, mais elle ne donne que l’évolution en profondeur de la permittivité électrique lorsque la conductivité du milieu n’est pas trop élevée. Afin de pouvoir fournir une image en profondeur de milieux complexes en une seule mesure invasive à l’aide d’un instrument unique donnant l’évolution de plusieurs paramètres physiques, le Demandeur a mis au point un instrument et un procédé de détermination du profil basé sur la mesure, à différentes fréquences, de la répartition, selon l’axe de propagation, de l’amplitude d’une onde stationnaire de champ électrique créée dans une ligne de transmission insérée dans le milieu stratifié à étudier. La vitesse de propagation de l’onde et son atténuation sont déterminées pour différentes couches du milieu à analyser. Cette vitesse et cette atténuation dépendent directement de la constante diélectrique ε r et de la conductivité électrique σ du milieu. L’analyse de la variation fréquentielle, et éventuellement temporelle, de ces deux paramètres et la mesure en parallèle du profil de température permet d’en déduire pour chaque couche des paramètres physiques secondaires qui dépendent de l’application visée : par exemple l’humidité ou la teneur en eau liquide, la salinité, la concentration en intrants, le taux de matière organique, ou l’équivalent en eau (pour la neige). Plus particulièrement, la présente invention a pour objet un instrument pour l’analyse d’un milieu complexe, en particulier un milieu naturel, pouvant être stratifié, afin d’en déterminer les propriétés physico-chimiques, cet instrument comprenant : - une ligne de transmission dite bifilaire, de longueur L, comprenant deux conducteurs métalliques parallèles et disposés en vis-à-vis symétriquement par rapport à un axe central x et étant terminée par un circuit ouvert permettant une réflexion totale de l’onde, - un générateur sinusoïdal RadioFréquence (usuellement désigné par l’acronyme RF) pour alimenter les conducteurs métalliques, - une pluralité de N détecteurs RF disposés entre les conducteurs métalliques en étant régulièrement espacés les uns des autres, chaque détecteur RF étant muni d’une antenne imprimée, chaque détecteur RF étant apte à convertir la puissance du signal capté par l’antenne qui lui est associée en une tension continue, - une carte superviseur pour contrôler le générateur sinusoïdal RF et la pluralité de détecteurs RF et leurs antennes imprimées associées. Le pas spatial de mesure est discret : en d’autres termes, les détecteurs RF et leurs antennes associées sont régulièrement répartis le long de la ligne de transmission et distants les uns des autres. De manière avantageuse, l’instrument selon l’invention peut également comprendre au moins un microcontrôleur comprenant un convertisseur analogique numérique A/N, le convertisseur analogique numérique A/N étant destiné à convertir chaque tension mesurée par les détecteurs RF en une valeur numérique destinée à être transmise à la carte superviseur. De manière avantageuse, les antennes imprimées peuvent être disposées entre lesdits conducteurs métalliques de manière décentrée par rapport à l’axe central x de la ligne de transmission. Ce décentrage des antennes imprimées a pour effet d’augmenter le rapport signal (capté par une antenne)/bruit. De manière avantageuse, les détecteurs RF peuvent être regroupés en modules comprenant chacun 8 détecteurs RF chainés en série et un microcontrôleur. La présente invention a également pour objet un procédé pour analyser, selon la profondeur x, un milieu complexe comprenant au moins une couche de matériau solide et/ou liquide de manière à déterminer des propriétés physiques de ladite couche de matériau solide et/ou liquide. Si le milieu complexe comporte une pluralité de couches, chaque couche doit être considérée comme un milieu continu (donc sans discontinuité brutale de ses propriétés physique au sein de la couche). Le procédé selon l’invention comprend les étapes suivantes : - introduction dans le milieu complexe à analyser, d’un instrument selon l’invention, la surface dudit milieu complexe définissant l’origine de l’analyse x=0 ; - alimentation, à l’aide du générateur sinusoïdal, de la ligne de transmission par un signal sinusoïdal de fréquence f variant entre 2 MHz et 2 GHz, de manière à générer un champ électrique E, induisant un courant dans chacune des antennes imprimées dont la puissance est convertie en une tension continue par le détecteur RF qui lui est associé, la propagation dudit champ électrique E le long de la ligne de transmission et entre les conducteurs métalliques dudit instrument se traduisant par l’apparition d’au moins une onde stationnaire de longueur d’onde λ et d’amplitude V 20 (z) en fonction de l’abscisse z avec z=L-x ; - conversion par le convertisseur analogique numérique A/N (usuellement désigné par l’acronyme en anglais ADC pour « Analog to Digital Converter ») des différentes tensions continues ainsi obtenues par les détecteurs RF en des valeurs numériques ; - transmission des valeurs numériques ainsi obtenues à la carte superviseur, qui est programmée pour en réaliser le post-traitement et les transformer en une courbe représentant l’évolution de l’amplitude V 20 (z) de l'onde stationnaire dans la couche de matériau solide et/ou liquide, selon l’abscisse z, avec z=L-x (en d’autres termes, les N détecteurs RF donnent indirectement une mesure spatialement échantillonnée de V 20 (z) ; - détermination par interpolation des valeurs des profondeurs x où apparaissent des tensions minimales et des tensions maximales et de l’amplitude des valeurs maximales de V 20 (z) ; - calcul de la permittivité électrique complexe ε*= ε r – j.σ/(2.π.f.ε 0 ) de ladite couche de matériau solide et/ou liquide en fonction de la profondeur x et pour chaque fréquence f de mesure ; avec ε 0 désignant la constante diélectrique du vide, ε r (partie réelle de la permittivité complexe ε*) désignant la constante diélectrique dans ladite couche de matériau solide et/ou liquide, et σ désignant la conductivité électrique σ de ladite couche de matériau solide et/ou liquide, f désignant la fréquence du signal sinusoïdal, et j désignant l’opérateur mathématique tel que j 2 =-1 le calcul de la permittivité électrique complexe ε* comprenant deux étapes : - une première étape comprenant la détermination de la demi-longueur d’onde λ/2 de l'onde stationnaire correspondant à la distance entre deux minima successifs de la courbe d’évolution de l’amplitude V20(z) de l'onde stationnaire, le calcul de la vitesse c de l'onde stationnaire et de la constante diélectrique ε r dans ladite couche de matériau solide et/ou liquide ; et - une deuxième étape de calcul de l’atténuation exponentielle α de l'onde stationnaire entre deux maxima successifs de la courbe d’évolution de l’amplitude V20(z) de l'onde stationnaire, cette atténuation dépendant directement de la partie imaginaire de la permittivité électrique complexe ε* de ladite couche de matériau solide et/ou liquide, ladite partie imaginaire dépendant elle-même directement de la conductivité électrique σ de ladite couche de matériau solide et/ou liquide. La superposition de l’onde progressive et de l’onde régressive se traduit par l’apparition d’une onde stationnaire (son amplitude est fixe à une abscisse z donnée mais elle reste bien sûr sinusoïdale en temps). Pour une ligne de transmission avec pertes, l’origine de l’axe z=0 étant prise sur la charge, on démontre que l’expression générale complexe de cette tension stationnaire s’écrit selon l’équation (1) : V(z) = Vi.(e γ z + Γ L .e- γ z ) (1) avec Vi définissant la tension imposée par le générateur de tension sinusoïdal en entrée de la ligne de transmission, considérée comme la référence des phases, Γ L étant le coefficient de réflexion sur la charge, qui est sensiblement égal à 1 car la charge est un circuit ouvert, z définissant l’abscisse selon la relation z=L-x L’allure générale de l’amplitude de V(z) est alors donnée par la , qui montre une courbe obtenue par la discrétisation sur N points de mesure de la formule suivante développée de l’équation (1) : (2) avec Vi, α (atténuation en Np/m), Γ L (de l’ordre de 1), β (constante de phase en rad/m) étant des paramètres qui sont fonction du milieu. La constante diélectrique effective (tenant compte de la géométrie de la ligne) peut être déduite directement de β, et la conductivité électrique déduite de α via des formules des lignes de propagation. Notons que δ exprime le fait que le dernier capteur n’est pas à l’abscisse z=0. En pratique, la période de l’ondulation observée étant égale à λ/2, on en déduit directement ε r effectif pour chaque demi-période spatiale détectée et, à la suite, la constante diélectrique ε r connaissant la loi de correspondance entre ε r effectif et ε r après l’étalonnage de l’instrument dans quatre milieux de constante diélectrique connue (air, sable sec, billes de verre et alcool éthylique). L’enveloppe de la courbe permet également de déduire α donc de remonter à la conductivité σ. Pour cela, nous mesurons les distances entre minima d’une part et la décroissance des maxima de l’amplitude d’autre part. De manière avantageuse, le procédé selon l’invention peut en outre comprendre une étape expérimentale de mesure de la température dudit milieu stratifié à analyser, pour augmenter la précision d’analyse. Cette étape pourra être par exemple réalisée à l’aide d’un profileur de température annexe. Selon un premier mode de réalisation du procédé selon l’invention, le milieu à analyser est un sol humide de type de type sol agricole, ou sol contenant des déchets granulaires ou une tourbière. Dans ce mode de réalisation, le procédé selon l’invention comprend en outre une étape de calcul avec des équations analytiques du taux d’humidité à partir de la constante diélectrique ε r dans la ou les couches de matériau solide suivi de l’évolution temporelle de ce taux en faisant une requête dans une base de données alimentées par les mesures de l’instrument. De préférence, le procédé selon ce mode de de réalisation, adapté pour l’analyse d’un sol humide, peut en outre comporter les étapes suivantes : - mesure du profil de conductivité σ du milieu à analyser calculé à partir de la partie imaginaire de la permittivité complexe dans la ou les couches de matériau solide ou liquide suivi de l’évolution temporelle de cette conductivité en fonction de la profondeur x en faisant une requête dans une base de données alimentées par les mesures de l’instrument, - à partir des mesures réalisées aux étapes précédentes, on en déduit en analysant la dérive temporelle de l’humidité et de la conductivité σ pour chaque fréquence de mesure et en la comparant avec des équations analytiques provenant d’une modélisation physico-chimique du milieu, le taux de matière organique (dans la ou les couches du milieu à analyser, et le taux d’intrants (N, P, K) et le taux de salinité. Selon un deuxième mode de réalisation du procédé selon l’invention, le milieu à analyser est un milieu naturel stratifié de type manteau neigeux. Dans ce mode de réalisation, le procédé selon l’invention comprend en outre les étapes suivantes : - mesures des pertes et de la constante diélectrique du manteau neigeux autour d’une fréquence de 2 GHz par détermination respectivement des maxima et des minima de l’onde stationnaire en hautes fréquences dans la ou lesdites couches de matériau solide ; - détermination de l’amplitude du champ électrique détecté à 2 MHz, pour en déduire la constante diélectrique du milieu dans la ou les couches de matériau solide ; - calcul à partir des différentes valeurs déterminées aux étapes précédentes, de la proportion volumique de glace, d’eau et d’air dans la ou lesdites couches de matériau solide, pour en déduire directement la hauteur de neige et des valeurs de Teneur en Eau Liquide (TEL) et d’Equivalent en Eau de la Neige (EEN) dans la ou les couches de matériau solide. La présente invention a encore pour objet l’utilisation du procédé selon l’invention déterminer si le milieu analysé (en particulier s’il s’agit d’un milieu naturel) est stratifié et comprend au moins deux couches de matériau solide et/ou liquide différentes l’une de l’autre par nature : cette détermination de l‘éventuelle stratification du milieu analyser est alors réalisée par observation de l’apparition d’un changement/ rupture au niveau de l’évolution de V 20 (z) , de la longueur d’onde λ et de l’atténuation α de l’onde stationnaire qui varient en fonction de la nature des couches de matériau solide et/ou liquide du milieu analysé. A partir d’un nombre suffisant de capteurs (dont le nombre dépend de la précision voulue sur la constante diélectrique), il est possible de déterminer si la ligne est plongée dans un milieu stratifié car au moins deux longueurs d’onde différentes seront détectées. D’autres avantages et particularités de la présente invention résulteront de la description qui va suivre, donnée à titre d’exemple non limitatif et faite en référence aux figures annexées et aux exemples : est une représentation schématique en coupe d’un exemple d’instrument de type sonde selon l’invention ; est une représentation schématique en coupe de l’instrument de type sonde illustré sur mis en place dans un sol complexe de type sol humide ; est une représentation typique de l’évolution spatiale de l’amplitude de l’onde de tension stationnaire V20(z) dans la ligne en présence d’une conductivité non nulle ; est un exemple de mesure réalisée à 1206 MHz montrant l’évolution de l’amplitude V20(z) de l'onde stationnaire dans une parcelle de terrain agricole (les points expérimentaux étant des petits carrés et la courbe modèle étant une ligne continue ) en fonction de la profondeur x ; est un exemple de mesure réalisée à 585 MHz montrant l’évolution de l’amplitude V 20 (z) de l'onde stationnaire dans de l’air, puis dans un sol de type remblais sec, conformément à la disposition illustrée par ; montre l’évolution de l’amplitude de l'onde stationnaire en microvolts en fonction de la profondeur x, à la fréquence f de 427 MHz, pour un sol de type manteau neigeux (exemple 3), montre l’évolution de la constante diélectrique à chaque abscisse, à f=3MHz, pour le sol de type manteau neigeux de l’exemple 3 ; montre l’évolution de l’amplitude de la tension réfléchie au niveau du sol dans un sable sec et humide, à l’aide d’une sonde TDR (exemple 4). Les figures et montrent un exemple d’instrument de mesure de type sonde selon la présente invention. est commentée dans la partie descriptive qui précède. Les figures et sont commentées dans les exemples 1 et 2 qui suivent. Cet instrument 1 comprend : - une ligne de transmission 3 bifilaire comprenant deux conducteurs métalliques 30, 31 parallèles et disposés en vis-à-vis symétriquement par rapport à un axe central x et étant terminée par un circuit ouvert permettant une réflexion totale de l’onde, - un générateur sinusoïdal RF 4 pour alimenter les conducteurs métalliques 30, 31, - une pluralité de N détecteurs RF 51, 52, 53, 54, 55, 56, 57, 58 disposés entre les conducteurs métalliques 30, 31 en étant régulièrement espacés les uns des autres, - chaque détecteur RF 51, 52, 53, 54, 55, 56, 57, 58 étant muni d’une antenne imprimée 61, 62, 63, 64, 65, 66, 67, 68, ces antennes imprimées 61, 62, 63, 64, 65, 66, 67, 68 étant de préférence disposées entre les conducteurs métalliques 30, 31 de manière décentrée par rapport à l’axe central x de la ligne de transmission 3, - chaque détecteur RF 51, 52, 53, 54, 55, 56, 57, 58 étant apte à convertir la puissance du signal capté par l’antenne 61, 62, 63, 64, 65, 66, 67, 68 qui lui est associée en une tension continue, - dans l’exemple de réalisation montré par les figures et , les détecteurs RF 51, 52, 53, 54, 55, 56, 57, 58 sont regroupés en modules comprenant chacun 8 détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) chainés et un microcontrôleur ADC/PIC 7. Ce microcontrôleur 7 comprend un convertisseur analogique numérique A/N pour convertir chaque tension mesurée par les détecteurs RF 51, 52, 53, 54, 55, 56, 57, 58 en une valeur numérique destinée à être transmise à une carte superviseur 8, qui contrôle le générateur sinusoïdal RF 4 et la pluralité de détecteurs RF 51, 52, 53, 54, 55, 56, 57, 58 et leurs antennes imprimées 61, 62, 63, 64, 65, 66, 67, 68 associées. EXEMPLES EXEMPLE 1 : détermination à l’aide du procédé selon l’invention de l’évolution de l’amplitude V20(z) de l'onde stationnaire dans une parcelle de sable humide agricole. On cherche à déterminer le profil en profondeur d’une parcelle de sable humide. Pour cela, on dispose l’instrument selon l’invention tel que représenté sur et dans la parcelle en creusant des avant-trous avec une tarière et, si nécessaire, avec un perforateur, puis en introduisant verticalement l'instrument dans les trous ainsi formés montre de manière très visible une atténuation exponentielle de l'amplitude de l'onde stationnaire due à la conductivité σ du milieu et une diminution graduelle avec la profondeur de la longueur d'onde guidée due à une augmentation de l'humidité avec x, ce qui est conforme avec ce qui était attendu et mesuré par une sonde capacitive. Les points de mesure sont indiqués par les carrés, le modèle issu de l'équation V20(z) et de la prise en compte du gradient d'humidité dans la variation de la constante diélectrique est indiqué en trait plein. Les points de mesure sont indiqués par les carrés, le modèle issu de l'équation V20(z) et de la prise en compte du gradient d'humidité dans la variation du constant diélectrique est indiqué en trait plein. EXEMPLE 2 : détermination à l’aide du procédé selon l’invention de l’évolution de l’amplitude V20(z) de l'onde lorsque l'instrument est plongé pour moitié dans un milieu de type remblais saturé en eau, l'autre partie restant dans l'air. On cherche à déterminer le profil en profondeur d'un sol de type remblais saturé en eau. Pour cela, on dispose l’instrument selon l’invention tel que représenté sur dans le sol conformément à la disposition montrée par et la procédure décrite dans l'exemple 1. On détermine alors, conformément au procédé selon l’invention, l’évolution de l’amplitude de l'onde stationnaire en microvolts en fonction de la profondeur x, comme illustrée par à la fréquence f de 585 MHz. permet de directement observer à partir des extrema une constante diélectrique effective de 1,6 sur la partie supérieure de l'instrument et de 5 sur la partie introduite dans le milieu. Ceci correspond parfaitement à de l’air sur la partie supérieure et à une humidité volumique de 25 % sur la partie immergée qui était saturée d’eau, valeur quasiment égale à celle mesurée par analyse gravimétrique sur la même parcelle (23%). On voit aussi l’interface à x=16 cm entre air et remblais, indiqué par la flèche. EXEMPLE 3 : détermination à l’aide du procédé selon l’invention de l’évolution de l’amplitude V20(z) de l'onde lorsque l'instrument est plongé dans un manteau neigeux. On cherche à déterminer le profil en profondeur d'un sol de type manteau neigeux. Pour cela, on dispose l’instrument selon l’invention tel que représenté sur dans le sol conformément à la disposition montrée par et la procédure décrite dans l'exemple 1. On détermine alors, conformément au procédé selon l’invention, l’évolution de l’amplitude de l'onde stationnaire en microvolts en fonction de la profondeur x, comme illustrée par à la fréquence f de 427 MHz. A 3 MHz (cf. ]) les détecteurs RF fournissent une amplitude directement proportionnelle au champ électrique qui dépend de la constante diélectrique à chaque abscisse. La courbe de représente cette constante diélectrique en intégrant le signal sur 5 détecteurs afin de réduire l'influence du bruit de mesure. et permettent ainsi de directement calculer la permittivité complexe à ces deux fréquences pour des couches d'épaisseur supérieure ou égale à 10 cm. Pour une de ces couches, on a trouvé ε r = 1,52 à 427 MHz et ε r = 1,6 à 3MHz. On en déduit en utilisant un modèle de loi de mélange des proportions volumiques de glace I = 29, 2 % et d'eau W = 0,5 %, ainsi qu'une densité D de 272 kg/m 3 . Ces résultats sont à comparer à ceux obtenus expérimentalement à l’aide d’une sonde capacitive associée à une méthode de pesage : I = 32,4 %, W = 0,6 %, densité D de 330 kg/m 3 . EXEMPLE 4 : (COMPARATIF) : détermination du profil d'humidité dans un sable sec et humide à l’aide la méthode de l’art antérieur de réflectométrie temporelle (TDR). On cherche à déterminer le profil en profondeur d'un sol sableux. Pour cela, on dispose l’instrument de TDR (typiquement une sonde constituée d'une ligne de transmission en trident métallique alimentée par un générateur d'impulsion) dans le sol. On détermine alors l’évolution de l’amplitude de la tension réfléchie au niveau du sol comme illustrée par . montre clairement que plus l'humidité du milieu augmente (donnée par le paramètre VMC (acronyme désignant en anglais « Volumetric Moisture Content » et noté ici Θ en %), plus l'impulsion réfléchie est déformée par les pertes que subit l'onde progressive et régressive tout au long du trajet sur la ligne trident, ce qui ne permet pas de déterminer avec une bonne précision le profil d'humidité et de conductivité du milieu mesuré, surtout si celui-ci est stratifié. Instrument (1) pour l’analyse d’un milieu complexe (2) pour en déterminer les propriétés physico-chimiques, ledit instrument comprenant : une ligne de transmission (3) de longueur L comprenant deux conducteurs métalliques (30, 31) parallèles et disposés en vis-à-vis symétriquement par rapport à un axe central (x) et étant terminée par un circuit ouvert permettant une réflexion totale de l’onde, un générateur sinusoïdal RF (4) pour alimenter lesdits conducteurs métalliques (30, 31), - une pluralité de détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) disposés entre lesdits conducteurs métalliques (30, 31) en étant régulièrement espacés les uns des autres, chaque détecteur RF (51, 52, 53, 54, 55, 56, 57, 58) étant muni d’une antenne imprimée (61, 62, 63, 64, 65, 66, 67, 68), chaque détecteur RF (51, 52, 53, 54, 55, 56, 57, 58) étant apte à convertir la puissance du signal capté par l’antenne (61, 62, 63, 64, 65, 66, 67, 68) qui lui est associée en une tension continu, une carte superviseur (8) pour contrôler ledit générateur sinusoïdal RF (4) et la pluralité de détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) et leurs antennes imprimées (61, 62, 63, 64, 65, 66, 67, 68) associées. Instrument (1) selon la revendication 1, comprenant en outre au moins un microcontrôleur (7) comprenant un convertisseur analogique numérique A/N, ledit convertisseur analogique numérique A/N étant destiné à convertir chaque tension mesurée par les détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) en une valeur numérique destinée à être transmise à la carte superviseur (8). Instrument (1) selon l’une quelconque des revendications 1 et 2, dans lequel les antennes imprimées (61, 62, 63, 64, 65, 66, 67, 68) sont disposées entre lesdits conducteurs métalliques (30, 31) de manière décentrée par rapport à l’axe central (x) de la ligne de transmission (3). Instrument (1) selon l’une quelconque des revendications 1 à 3, dans lequel les détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) sont regroupés en modules comprenant chacun 8 détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) chainés et un microcontrôleur. Procédé pour analyser, selon la profondeur x, un milieu complexe (2) comprenant au moins une couche de matériau solide et/ou liquide (20), de manière à déterminer des propriétés physiques de ladite couche de matériau solide et/ou liquide (20), ledit procédé comprenant les étapes suivantes : (100) introduction dans ledit milieu complexe (2) à analyser, d’un instrument (1) tel que défini selon l’une quelconque des revendications 1 à 3, la surface dudit milieu complexe (2) définissant l’origine de l’analyse x=0 ; (110) alimenter, à l’aide dudit générateur sinusoïdal (4), ladite ligne de transmission (2) par un signal sinusoïdal de fréquence f variant entre 2 MHz et 2 GHz, de manière à générer un champ électrique E, induisant un courant dans chacune des antennes imprimées (61, 62, 63, 64, 65, 66, 67, 68) dont la puissance est convertie en une tension continue par le détecteur RF (51, 52, 53, 54, 55, 56, 57, 58) associé, la propagation dudit champ électrique E le long de la ligne de transmission (3) et entre les conducteurs métalliques (30, 31) dudit instrument (1) se traduisant par l’apparition d’au moins une onde stationnaire de longueur d’onde λ et d’amplitude V20(z) en fonction de l’abscisse z avec z=L-x : (120) conversion par le convertisseur analogique numérique A/N des différentes tensions continues ainsi obtenues par les détecteurs RF (51, 52, 53, 54, 55, 56, 57, 58) en des valeurs numériques ; (130) transmission desdites valeurs numériques ainsi obtenues à ladite carte superviseur (8) programmée pour en réaliser le post-traitement et les transformer en une courbe représentant l’évolution de l’amplitude V20(z) de l'onde stationnaire dans ladite couche de matériau solide et/ou liquide (20), selon l’abscisse z, avec z=L-x ; (150) détermination par interpolation des valeurs des profondeurs x où apparaissent des tensions minimales et des tensions maximales et de l’amplitude des valeurs maximales de V20(z) ; (170) calcul de la permittivité électrique complexe ε* = ε r – j.σ/(2.π.f.ε 0 ) de ladite couche de matériau solide et/ou liquide (20) en fonction de la profondeur x et pour chaque fréquence f de mesure, avec ε 0 désignant la constante diélectrique du vide, ε r (partie réelle de la permittivité complexe ε*) désignant la constante diélectrique dans ladite couche de matériau solide et/ou liquide (20), et σ désignant la conductivité électrique de ladite couche de matériau solide et/ou liquide (20), j désignant l’opérateur mathématique tel que j 2 =-1, ledit calcul comprenant deux étapes : (1701) une première étape comprenant la détermination de la demi-longueur d’onde λ/2 de l'onde stationnaire correspondant à la distance entre deux minima successifs de la courbe d’évolution de l’amplitude V20(z) de l'onde stationnaire, le calcul de la vitesse c de l'onde stationnaire et de la constante diélectrique ε r dans ladite couche de matériau solide et/ou liquide (20) ; (1702) une deuxième étape de calcul de l’atténuation exponentielle α de l'onde stationnaire entre deux maxima successifs de la courbe d’évolution de l’amplitude V20(z) de l'onde stationnaire, cette atténuation dépendant directement de la partie imaginaire de la permittivité électrique complexe ε* de ladite couche de matériau solide et/ou liquide (20), ladite partie imaginaire dépendant elle-même directement de la conductivité électrique σ de ladite couche de matériau solide et/ou liquide (20). Procédé selon la revendication 5, comprenant en outre une étape (160) de mesure de la température dudit milieu stratifié (1) à analyser. Procédé selon l’une quelconque des revendications 5 et 6, selon lequel le milieu (1) à analyser est un sol humide, ledit procédé comprenant en outre l’étape suivante : (180) calcul avec des équations analytiques du taux d’humidité à partir de de la constante diélectrique ε r dans la ou lesdites couches de matériau solide et suivi de l’évolution temporelle de ce taux en faisant une requête dans une base de données alimentée par les mesures de l’instrument. Procédé selon la revendication 7, comprenant en outre les étapes suivantes : (190) mesure du profil de conductivité σ du milieu (1) à analyser calculé à partir de la partie imaginaire de la permittivité complexe dans la ou lesdites couches de matériau solide ou liquide et possibilité de suivre l’évolution temporelle de cette conductivité σ en fonction de la profondeur x en faisant une requête dans une base de données alimentée par les mesures de l’instrument , (200) à partir des mesures réalisées aux étapes (180) et (190), on en déduit en analysant la dérive temporelle de l’humidité et de la conductivité σ pour chaque fréquence de mesure et en la comparant avec des équations analytiques provenant d’une modélisation physico-chimique du milieu, e taux de matière organique (MO) dans la ou les couches (20) du milieu à analyser, le taux d’intrants (N, P, K) et de salinité. Procédé selon l’une quelconque des revendications 5 et 6, selon lequel le milieu (1) à analyser est un milieu naturel stratifié de type manteau neigeux, ledit procédé comprenant en outre l’étape suivante : (200) mesures des pertes et de la constante diélectrique du manteau neigeux autour d’une fréquence de 2 GHz par détermination respectivement des maxima et des minima de l’onde stationnaire en hautes fréquences (HF) dans la ou lesdites couches de matériau solide ; (210) détermination de l’amplitude du champ électrique détecté à 2 MHz permettant d’en déduire la constante diélectrique du milieu dans la ou lesdites couches de matériau solide ; (220) Calcul à partir des différentes valeurs déterminées aux étapes (200) et (210) de la proportion volumique de glace, d’eau et d’air dans la ou lesdites couches de matériau solide (10,11), permettant d’en déduire directement la hauteur de neige et des valeurs de Teneur en Eau Liquide (TEL) et d’Equivalent en Eau de la Neige (EEN) dans la ou lesdites couches de matériau solide (10,11). Utilisation du procédé selon l’une quelconque des revendications 5 ou 6, pour déterminer si le milieu (1) analysé est stratifié et comprend au moins deux couches de matériau solide et/ou liquide (20, 21) différentes l’une de l’autre par nature, par observation de l’apparition d’un changement au niveau de l’évolution de V 20 (z), de la longueur d’onde λ et de l’atténuation de l’onde stationnaire qui varient en fonction de la nature des couches de matériau solide et/ou liquide (20, 21) dudit milieu (1) analysé.