Brevet FR3155909A1 — publié le 20250530

PROCÉDÉ PLUVIOMÉTRIQUE D’ESTIMATION D’UN VOLUME D’EAU LIQUIDE PRÉCIPITANTE AU-DESSUS D’UNE POSITION GÉOGRAPHIQUE L’invention concerne un procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique mis en œuvre par un dispositif de réception et de traitement de signaux GNSS comprenant : une antenne apte à capter des signaux GNSS et positionnée au niveau de ladite position géographique; et une unité de calcul apte à recevoir et à traiter les signaux GNSS captés par ladite antenne ; le procédé comprenant les étapes suivantes : une étape de réception d’au moins un signal GNSS, dit signal reçu, dont le spectre est composé d’au moins deux fréquences distinctes ; une étape de détermination d’un retard total au zénith dudit signal reçu au niveau de ladite position géographique ; une étape de détermination d’un volume d’eau liquide précipitante au-dessus de ladite position géographique à partir uniquement dudit retard total au zénith. Fig. pour l’abrégé : Fig. 1 Procédé pluviométrique d’estimation du volume d’eau précipitante et dispositif météorologique correspondant Domaine technique de l’invention L’invention se rapporte au domaine technique de la pluviométrie satellitaire et des objets connectés pluviométriques. Elle concerne plus particulièrement un procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique prédéterminée par analyse de signaux GNSS. L’invention a également pour objet un dispositif météorologique mettant en œuvre un tel procédé et un système d’estimation pluviométrique mettant en œuvre un tel procédé d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique. Arrière-plan technologique L’estimation de quantités pluviométriques, et notamment l’estimation du taux de précipitation, est une composante essentielle dans la prédiction météorologique. L’estimation du taux de précipitation est cruciale dans la prédiction de scénarios météorologiques extrêmes tels que les crues (en particulier les crues dites « éclair »), les débordements de cours d’eau, ainsi que des phénomènes associés tels que les ruptures de barrages ou encore les inondations de zones habitées. Dès lors, la gestion réussie de la sécurité civile face à ces scénarios est conditionnée par un lancement d’alerte précoce, la plus exacte possible, permettant de prévenir les populations concernées le plus rapidement possible et, le cas échéant, d’exécuter des plans de sauvegarde et d’évacuation des zones potentiellement affectées afin de mettre la population qui y demeure en sécurité. Il existe aujourd’hui un ensemble riche de techniques pluviométriques permettant d’estimer le taux de précipitation ; ces techniques sont basées sur l’un ou plusieurs des éléments suivants : réseaux de stations météorologiques, radars météorologiques ; des satellites permettant d’observent les nuages et leur évolution. Bien que capables de détecter les précipitations et d’indiquer leur intensité et leur localisation, les principaux inconvénients des techniques précitées résident dans leur faible résolution spatiale due notamment à une absence, totale ou partielle, de couverture du réseau. Or, une meilleure résolution spatiale et/ou une densité plus importante, notamment de l’ordre du kilomètre, sont nécessaires pour pouvoir estimer avec précision le taux de précipitation afin de pouvoir réaliser des prédictions suffisamment pertinentes, notamment aux abords des villages en bord de rivière. Par ailleurs, les technologies susmentionnées sont difficilement accessibles par des entités locales comme des mairies ou des associations civiles. Ceci est dû, d’une part, à leur coût et technicités élevés ; et d’autre part, au manque d’adaptabilité et de déploiement desdites technologies sur une petite échelle. Il existe donc un besoin pour un procédé permettant d’estimer avec une grande précision des quantités pluviométriques, telles que le taux de précipitation, sur une zone géographique donnée, en particulier pour des zones dont la surface est de l’ordre du kilomètre carré. Il existe également un besoin pour un système d’estimation pluviométrique robuste, fiable et à faible coût permettant de mettre en œuvre le procédé d’estimation pluviométrique recherché ci-dessus. De plus, ledit système doit pouvoir être déployé simplement et rapidement tout en étant résistant aux intempéries et facilement réparable. Objectifs de l’invention La présente invention a pour objectif de remédier à tout ou partie des inconvénients de l’art antérieur, notamment ceux exposés ci-avant, en proposant une solution permettant de déterminer le volume d’eau liquide précipitante au-dessus d’une position géographique à partir de la seule mesure d’un retard total au zénith d’un signal GNSS. L’invention vise à fournir un procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique. L’invention vise en particulier à fournir, dans au moins un mode de réalisation, un dispositif météorologique apte à réceptionner et traiter des signaux GNSS et à mettre en œuvre le procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique. L’invention vise aussi à fournir, dans au moins un mode de réalisation de réalisation, un système d’alerte hydrologique dans une zone géographique comprenant des cours d’eau et capable d’alerter sur des risques d’inondation et/ou de crue. Pour ce faire, l’invention concerne un procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique donnée. Le procédé pluviométrique comprend les étapes suivantes : une étape de réception, par une antenne, dite antenne de réception, agencée à ladite position géographique prédéterminée et apte à capter des signaux GNSS, d’au moins un signal GNSS (21), dit signal reçu, dont le spectre est composé d’au moins deux fréquences distinctes ; une étape de détermination par une unité de calcul apte à recevoir et à traiter des signaux GNSS, d’un retard total au zénith de chaque signal reçu par ladite antenne de réception à ladite position géographique ; une étape de détermination, par ladite unité de calcul, d’un volume d’eau liquide précipitante au-dessus de ladite position géographique prédéterminée à partir uniquement dudit retard total au zénith déterminé à l’étape précédente et de données atmosphériques relatives à ladite position géographique prédéterminée. Un procédé pluviométrique selon l’invention permet donc d’estimer le volume d’eau liquide précipitante contenu dans la colonne d’atmosphère au-dessus de la position géographique donnée à partir du seul retard total au zénith du signal GNSS reçu. Si on le souhaite, une pluralité de signaux GNSS peuvent être reçus par l’antenne de réception en vue d’améliorer la précision de l’estimation du volume d’eau liquide précipitante au-dessus de ladite colonne de ladite position géographique. La pluralité de signaux GNSS peut notamment provenir d’une pluralité de satellites distincts, notamment des satellites appartenant à des constellations de satellites distinctes telles que GPS, Galileo, GLONASS, Bei-Dou. Ces différentes constellations émettent des signaux dans des bandes de fréquences bien définies. L’invention peut employer des signaux dont le spectre comprend des fréquences dans l’une quelconques desdites bandes, notamment les bandes GPS L1, L2 et L5 ; et les bandes Galileo E1, E6, E5a, E5b. Il est d’usage d’employer L1 = 1575, 42 MHz et L2 = 1277, 60 MHz comme deux fréquences porteuses utilisées pour transmettre une information codée sous forme binaire. La précision dans l’estimation du volume d’eau liquide précipitante contenu dans la colonne de fluide au-dessus de la position géographique donnée est fonction de différents paramètres, l’un d’entre eux est la précision dans la détermination du retard total au zénith qui, à son tour, dépend de la connaissance précise de la localisation géographique de l’antenne de réception. Dans un mode de réalisation de l’invention, l’unité de calcul déterminant un volume d’eau liquide précipitante au-dessus de ladite position géographique prédéterminée à partir uniquement dudit retard total au zénith est une unité de calcul distante, notamment faisant partie d’un serveur informatique distant. Cette caractéristique permet de configurer l’unité de calcul selon les besoins précis en termes de puissance de calcul. Notamment afin de réaliser les calculs mathématiques nécessaires pour la détermination dudit volume d’eau liquide précipitante. Dans un mode de réalisation de l’invention, le procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique donnée comprend une étape préliminaire de positionnement de l’antenne de réception par une méthode de positionnement absolu. En procédant ainsi, la précision de la détermination du retard total au zénith, et par conséquent de celle de l’estimation du volume d’eau liquide précipitante, sont améliorées. Dans un mode de réalisation cumulatif de l’invention, la méthode de positionnement absolu est une méthode de positionnement ponctuel précis. La méthode de positionnement ponctuel précis permet un positionnement précis à l’aide d’un seul récepteur et ne dépend pas de la disponibilité de données provenant de stations de référence à proximité de l’utilisateur. De plus, cette méthode permet d’assurer une précision centimétrique de positionnement de l’antenne de réception. En variante, la méthode positionnement absolue est une méthode par système de renforcement par satellite SBAS ou encore une méthode de positionnement ponctuel précis PPP-RTK. Dans un mode de réalisation de l’invention, l’étape de réception d’au moins un signal GNSS comprend en outre les sous-étapes de prétraitement suivantes : une sous-étape d’horodatage de réception, par l’unité de calcul, d’au moins un signal reçu ; une sous-étape de décodage, par l’unité de calcul, de l’horodatage d’émission dudit signal reçu ; une sous-étape de détermination et de correction des retards et/ou des perturbations engendrés par la traversée de l’ionosphère dudit signal reçu à partir desdites au moins deux fréquences distinctes composant le spectre dudit signal reçu ; une sous-étape de détermination, par l’unité de calcul, d’au moins un retard total oblique du au moins à un signal reçu à partir desdits horodatages d’émission et de réception dudit au moins un signal reçu. Ainsi, les sous-étapes de décodage d’horodatage de l’émission et d’horodatage de la réception du signal permettent de déterminer le temps de parcours du signal entre le satellite au moment de l’émission dudit signal et l’antenne de réception. L’ionosphère étant un milieu dispersif, celle-ci propage les différentes fréquences composant le signal émis par le satellite à des vitesses différentes, engendrant ainsi un retard du signal, dit retard ionosphérique. Ce retard peut être corrigé au moyen de modèles mathématiques tels que les modèles « ionosphere free », « de Klobuchar » et « de NeQuick-G ». Dans le cas du modèle « ionosphere free », les effets ionosphériques de premier ordre sur les mesures du code et de la phase de la porteuse dépendent de l'inverse du carré de la fréquence du signal. Par conséquent, les récepteurs bi-fréquence peuvent éliminer les effets ionosphériques par une combinaison linéaire des mesures du code ou de la porteuse. Le passage du signal GNSS à travers de la troposphère engendre un retard du signal, dit retard troposphérique. Ce retard résulte d’un changement de l’indice de réfraction du milieu de propagation dû à l’humidité de l’air et aux variations de pression et de température. Le retard troposphérique est ainsi fortement corrélé avec la saison, l’heure locale et la latitude. La troposphère est composée de gaz neutres et ne constitue pas un milieu dispersif pour les signaux électromagnétiques, ce qui implique que le retard atmosphérique affecte les fréquences GNSS de la même façon. Les espèces principales de l’atmosphère proche, l’azote et l’oxygène, correspondent à la partie sèche de la troposphère. Ces effets contribuent pour 90 % au délai total et sont à dynamique lente, soit environ 1 % de variation en quelques heures, et sont donc assez bien modélisés. La partie humide de la troposphère, c’est-à-dire la vapeur d’eau, est responsable des 10 % du délai restant ; la raison principale de ce retard est dû au moment dipolaire des molécules de vapeur d’eau. Ainsi, le délai troposphérique s’écrit comme la somme de ces deux contributions. Dans le contexte de la présente invention, on entend par « retard total oblique » le retard total que subit le signal radio GNSS en raison de l'atmosphère neutre le long du trajet entre un satellite et une antenne réceptrice au sol. Dans un mode de réalisation de l’invention, la valeur d’un retard total au zénith du signal reçu au niveau de la position géographique est déterminée à partir dudit retard total oblique du signal reçu précédemment obtenu. Selon cette caractéristique, le retard total oblique est transformé en un retard total au zénith permettant d’effectuer la détermination du volume d’eau liquide précipitante sur la colonne verticale au-dessus de la position géographique d’intérêt. Avantageusement, cette transformation est effectuée à l’aide d’un fonction de rabattement décrivant la dépendance du retard par rapport à l’angle d’élévation ce qui permet de projeter un délai oblique en un délai au zénith en fonction de l’élévation locale de la ligne de visée. Des exemples de telles fonctions de rabattement sont celles connues sous le nom de fonctions de rabattement de « Mariani » , de « Davis » et de « Herring ». La fonction de rabattement de « Mariani » est donnée par la fraction continue suivante : De par la nature du retard total au zénith du signal reçu, celui-ci fournit une information sur la masse d'eau sous forme liquide et gazeuse dans la colonne au-dessus de la position géographique. Dans un mode de réalisation de l’invention, l’étape de détermination d’un volume d’eau liquide précipitante au-dessus de ladite position géographique comporte les sous-étapes suivantes : une sous-étape de détermination d’une composante sèche dudit retard total au zénith du signal reçu ; une sous-étape de détermination du volume d’eau liquide précipitante au-dessus de la position géographique à partir dudit retard total au zénith du signal reçu et de la composante sèche dudit retard total au zénith du signal reçu; En procédant ainsi, l’invention cherche à différencier les contributions dites « sèche » et « humide » constituant ledit retard total au zénith du signal reçu. Si on le souhaite, la sous-étape de détermination d’une composante sèche du retard total au zénith du signal reçu, notée ZHD, peut être obtenue par un modèle physique hydrostatique, des exemples non limitatifs de tels modèles sont le modèle « Saastamoinen », le modèle de « Hopfield » ou encore le modèle connu sous le nom de « MOPS ». Les expressions mathématiques des modèles de « Hopfield » et « Saastamoinen » sont données respectivement par les équations suivantes : Dans les expressions mathématiques précédentes, P est la surface de pression, T la température, h la hauteur ellipsoïdale et la latitude. Dans un mode de réalisation de l’invention, la sous-étape de détermination du volume d’eau liquide précipitante au-dessus de la position géographique à partir dudit retard total au zénith est effectuée par un algorithme d’apprentissage automatique. Ce mode de réalisation avantageux permet de disposer d’un outil pour déterminer rapidement le volume d’eau liquide précipitante à partir de la seule valeur du retard total au zénith. Cette caractéristique est particulièrement intéressante dans la prise de décisions associées au volume des précipitations, des tels cas sont rencontrés par exemple dans l’évaluation des risques d’inondation et/ou des risques de débordement de cours d’eau. Par ailleurs, le pouvoir de généralisation des algorithmes d’apprentissage automatique permet de traiter l’ensemble de cas possibles selon les combinaisons des paramètres météorologiques. Dans un mode de réalisation de l’invention, l’algorithme d’apprentissage automatique est un prédicteur scalaire prenant en entrée une valeur de retard total au zénith d’un signal GNSS ; ledit prédicteur scalaire étant entraîné sur un ensemble de données relatives à des valeurs de retard total au zénith de signaux GNSS, de la position géographique et de données météorologiques. En procédant ainsi, l’invention permet d’établir les liens de dépendance entre les différents paramètres météorologiques, la position géographique et les valeurs de retard total au zénith. Ainsi, il est possible d’estimer le volume d’eau liquide précipitante à partir de la seule valeur du retard total au zénith. Selon un mode de réalisation avantageux, ledit prédicteur scalaire est un réseau de neurones à graphes dont l’un des sommets est constitué de l’architecture du réseau « MetNet », cette partie permet de prendre en compte les différentes conditions météorologiques. L’invention utilise la technique d’apprentissage par transfert en complétant le graphe par : un sommet d’entrée constitué d’un réseau de neurones de type perceptron à couches denses prenant en compte le retard total au zénith ; et par un sommet de sortie constitué d’un seul neurone et d’une fonction d’activation renvoyant une valeur réelle correspondant audit volume d’eau liquide précipitante recherché. Dans un mode de réalisation alternatif de l’invention, la sous-étape de détermination du volume d’eau liquide précipitante au-dessus de la position géographique à partir dudit retard total au zénith est effectuée par un modèle physique de la partie d’eau gazeuse de la troposphère. Cette caractéristique avantageuse permet alternativement de déduire le volume d’eau liquide précipitante en soustrayant au retard total au zénith la partie sèche et la partie humide sous forme de vapeur. L’invention a également pour objet un dispositif météorologique capable de réceptionner et de traiter des signaux GNSS et apte à mettre en œuvre le procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une position géographique selon l’invention. Le dispositif météorologique permet de mesurer le volume d’eau précipitante au-dessus de ladite position géographique. Dans un mode de réalisation de l’invention, ledit dispositif météorologique comprend des moyens de communications avec un serveur météorologique distant, comprenant notamment des données météorologiques telles que la température, l'humidité, la vitesse et la direction du vent, les précipitations, la pression atmosphérique, ainsi que des prévisions à court et à long terme d’orages et d’autres phénomènes météorologiques. Avantageusement, les données météorologiques sont régulièrement mises à jour et permettent au dispositif météorologique de disposer des informations météorologiques précises et actualisées. De plus, le dispositif météorologique peut également recevoir des alertes météorologiques en temps réel en cas de conditions météorologiques extrêmes. Avantageusement les dispositifs météorologiques peuvent être des dispositifs météorologiques de surveillance de cours d’eau lorsque ceux-ci sont positionnées au niveau d’un cours d’eau tels que décrits dans décrite dans le document de brevet WO2023174922 du demandeur décrivant un dispositif dédié à mesurer la hauteur de la surface de l'eau et la vitesse de surface des systèmes hydrologiques par télédétection. Selon cette caractéristique avantageuse, la détermination du volume d’eau précipitante au-dessus de l’emplacement dudit dispositif pluviométrique permet de déterminer la contribution fluviale directement intégré sur le cours d’eau au niveau dudit emplacement géographique. L’invention a également pour objet un système météorologique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une zone géographique comprenant une pluralité de dispositifs météorologiques, selon l’invention, distribués au sein de ladite zone géographique. La pluralité de dispositifs permet de couvrir une zone plus large qu’une seule position géographique et d’estimer par pondération et/ou extrapolation un volume d’eau liquide précipitante dans ladite zone géographique. Dans un mode de réalisation de l’invention, un dispositif météorologique parmi la pluralité de dispositifs météorologiques est un dispositif maître chargé de centraliser les mesures effectuées par l’ensemble des dispositifs météorologiques distribués dans ladite zone géographique. Selon cette caractéristique avantageuse, il est possible de faire coopérer l’ensemble des dispositifs météorologiques afin de centraliser les informations des volumes d’eau liquide précipitante calculés par chacun desdits dispositifs. Dans un mode de réalisation de l’invention, chaque dispositif météorologique peut être désigné, à n’importe quel moment souhaité, maître par rapport à l’ensemble complémentaire de dispositifs météorologiques. Selon cette caractéristique avantageuse, il est possible de faire fonctionner le système même lorsqu’une panne est détectée dans l’un des dispositifs météorologiques et notamment le dispositif maître. Avantageusement encore, des tels dispositifs météorologiques de surveillance des cours d’eau peuvent comporter des moyens de mesure du niveau de l’eau s’écoulant dans un cours d’eau ainsi que des moyens de mesure du débit dudit cours d’eau, notamment au travers de scanners LIDAR, tels que décrit dans le brevet du demandeur. L’invention a également pour objet un système d’alerte hydrologique dans une zone géographique présentant des cours d’eau, ledit système d’alerte comprenant : une pluralité de dispositifs météorologiques de surveillance de cours d’eau déployées au niveau de chaque cours d’eau de ladite zone géographique, chacun desdits dispositifs météorologiques de surveillance comprenant un dispositif météorologique selon l’invention et un dispositif de mesure du niveau desdits cours d’eau; un serveur distant de surveillance météorologique apte à communiquer avec la pluralité de dispositifs météorologiques de surveillance ; caractérisé en ce que ledit système d’alerte est apte à émettre une alerte de crue, pour l’un quelconque desdits cours d’eau, à partir du niveau desdits cours d’eau et/ou d’une estimation d’un volume d’eau liquide précipitante dans ladite zone géographique obtenue par le procédé selon l’invention et d’un seuil prédéterminé, dit seuil d’alerte. Selon cette caractéristique avantageuse, la présence d'une pluralité de dispositifs météorologiques de surveillance de cours d'eau déployés au niveau de chaque cours d'eau permet une surveillance continue et en temps réel des conditions météorologiques et hydrologiques. Ces dispositifs météorologiques de surveillance fournissent des données précises et actualisées sur les conditions météorologiques et le niveau des cours d'eau, ce qui permet d'évaluer de manière fiable le risque d'inondation. De plus, l'utilisation d'un serveur distant de surveillance météorologique permet une communication efficace et rapide entre le serveur et les dispositifs météorologiques de surveillance. Cette communication en temps réel facilite la collecte, l'analyse et l'interprétation des données météorologiques et hydrologiques, permettant ainsi au système d'alerte d'émettre des alertes de crue précises et pertinentes. L'alerte de crue est déclenchée non seulement à partir du niveau des cours d'eau, mais aussi à partir d'une estimation du volume d'eau liquide précipitante dans la zone géographique, ce qui renforce la fiabilité et la précision du système d'alerte. Le seuil d’alerte prend typiquement la forme d’un volume d’eau en mètres cubes, mais il peut également être défini en termes de hauteur en mètres, en termes d’un débit en mètres cubes par seconde, etc. Dans un mode de réalisation de l’invention, l’unité de calcul est une unité de calcul distante, notamment une unité de calcul d’un ordinateur distant exécuté dans un serveur informatique. Liste des Figures D'autres buts, caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante donnée à titre uniquement non limitatif et qui se réfère aux figures annexées dans lesquelles : FIG. 1 est un schéma-bloc d’un procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante selon un mode de réalisation de l’invention. FIG. 2 est une vue schématique d’un dispositif de réception et de traitement de signaux GNSS apte à mettre en œuvre le procédé de la FIG. 1 selon un mode de réalisation de l’invention. FIG. 3 est une vue schématique des configurations géométriques définies par la propagation de signaux GNSS depuis des satellites vers un dispositif météorologique de surveillance apte à mettre en œuvre le procédé de la FIG. 1 selon un mode de réalisation de l’invention. FIG. 4 est une vue schématique d’un dispositif météorologique de surveillance selon un mode de réalisation de l’invention. FIG. 5 est une vue schématique d’un système d’alerte hydrologique dans une zone géographique comprenant une pluralité de dispositifs météorologiques de surveillance selon un mode de réalisation de l’invention. Description détaillée d’un mode de réalisation de l’invention Sur les figures, les échelles et les proportions ne sont pas strictement respectées et ce, à des fins d’illustration et de clarté. En outre, les éléments identiques, similaires ou analogues sont désignés par les mêmes références dans toutes les figures. La FIG. 1 représente schématiquement et partiellement un schéma-bloc d’un procédé pluviométrique d’estimation d’un volume d’eau liquide précipitante V au-dessus d’une position géographique donnée G. Le procédé commence par une étape préliminaire E0 de positionnement par une méthode de positionnement ponctuel précis d’une antenne de réception 11, dite antenne de réception, apte à capter des signaux GNSS 21 agencée à ladite position géographique G prédéterminée. Consécutivement à l’étape préliminaire E0 de positionnement de l’antenne de réception 11, le procédé comporte une étape E1 de réception, par ladite antenne 11, d’au moins un signal GNSS 21, dit signal reçu, dont le spectre est composé d’au moins deux fréquences distinctes, en l’occurrence des fréquences L1 = 1575, 42 MHz et L2 = 1277, 60 MHz. L’étape E1 de réception du au moins un signal GNSS 21 comprend en outre les sous-étapes de prétraitement suivantes : une sous-étape E1.1 d’horodatage de réception, par une unité de calcul 12 apte à recevoir et à traiter des signaux GNSS de chaque signal reçu 21 par ladite antenne de réception 11 à ladite position géographique G ; une sous-étape E1.2 de décodage, par l’unité de calcul 12, de l’horodatage d’émission dudit au moins un signal reçu . Les sous-étapes E1.1 et E1.2 permettent de déterminer le temps de parcours du signal 21 entre le satellite 20 au moment de l’émission dudit signal et l’antenne de réception. 11 L’étape E1 de réception du au moins un signal GNSS 21 comprend en outre les sous-étapes de prétraitement suivantes : une sous-étape E1.3 de détermination et de correction, par la méthode de « Klobuchar », des retards et/ou des perturbations engendrés par la traversée de l’ionosphère dudit signal reçu 21 à partir desdites au moins deux fréquences distinctes composant le spectre dudit signal reçu 21 ; une sous-étape E1.4 de détermination, par l’unité de calcul 12, d’au moins un retard total oblique du au moins un signal reçu 21 à partir desdits horodatages d’émission et de réception dudit au moins un signal reçu 21. Consécutivement, le procédé comporte une étape E2 de détermination par l’unité de calcul 12, d’un retard total au zénith τ dudit signal reçu 21 au niveau de ladite position géographique G. Le retard total au zénith τ est déterminé à partir du retard total oblique du signal reçu 21 à l’aide d’une fonction de rabattement de « Mariani ». Enfin, le procédé termine par une étape E3 de détermination, par l’unité de calcul 12, d’un volume V d’eau liquide précipitante au-dessus de ladite position géographique G à partir uniquement dudit retard total au zénith τ déterminé à l’étape E2 et de données atmosphériques relatives à ladite position géographique prédéterminée G. Pour ce faire l’unité de calcul exécute un algorithme d’apprentissage automatique sous forme d’un prédicteur scalaire prenant en entrée la valeur de retard total au zénith τ d’un signal GNSS 21; ledit prédicteur scalaire étant entraîné sur un ensemble de données relatives à des valeurs de retard total au zénith τ de signaux GNSS 21, de la position géographique G et de données météorologiques historiques. La FIG. 2 représente schématiquement un dispositif 10 de réception et de traitement de signaux GNSS 21 émis par un satellite 20 ; le dispositif 10 est apte à mettre en œuvre le procédé de la FIG. 1 selon un mode de réalisation de l’invention. Le dispositif 10 comprend une antenne 11 de réception de signaux GNSS 21 et une unité de calcul 12. L’antenne 11 est positionnée au niveau de la position géographique G à l’aide d’une méthode de positionnement ponctuel précis. La FIG. 3 représente une vue schématique des configurations géométriques définies par la propagation de signaux GNSS 21 depuis des satellites 20 vers un dispositif météorologique de surveillance positionné au point G et apte à mettre en œuvre le procédé de la FIG. 1 selon un mode de réalisation de l’invention. Sont également représentées l’ionosphère 100 et la troposphère 200. Le centre de la Terre est représenté par le point O et l’angle α d’élévation du satellite 20 est l’angle formé par l’espace tangent T à la Terre au point G. Est également représenté la direction du zénith Z en G. La FIG. 4 représente une vue schématique d’un dispositif 10 météorologique de surveillance, réception et traitement de signaux GNSS selon un mode de réalisation de l’invention. Le dispositif 10 comprend une antenne 11 de réception de signaux GNSS, une unité de calcul 12 apte à analyser lesdits signaux GNSS. Le dispositif 10 comprend en outre une mémoire de stockage 15 apte à enregistrer les mesures de volume d’eau liquide précipitante au cours du temps, ainsi que des moyens de communication 17 aptes à communiquer avec des dispositifs externes tels que des serveurs distants. Le dispositif 10 est également doté d’une source de puissance électrique rechargeable (non représentée) ainsi que des moyens d’alimentation électrique et de recharge électrique de ladite source de puissance électrique sous forme de panneaux solaires 16. En fin, le dispositif 10 est muni d’une caméra permettant l’enregistrement et la détermination de la vitesse de surface d’un cours d’eau ; et d’un capteur 13 du niveau d’un cours d’eau, sous forme de scanner LIDAR 13. On a représenté en FIG. 5 une vue schématique d’un système d’alerte hydrologique dans une zone géographique comprenant une pluralité de dispositifs météorologiques de surveillance 10 selon un mode de réalisation de l’invention La pluralité de dispositifs météorologiques de surveillance 10 sont positionnées au niveau d’un cours d’eau 50 et mesurent régulièrement le niveau du cours d’eau ainsi que déterminent la valeur du volume d’eau liquide précipitante au-dessus de l’emplacement de chacun desdits dispositifs par le procédé de la FIG. 1 . Les dispositifs 10 maintiennent des communications avec un serveur distant 60 contenant des données météorologiques régulièrement mises à jour. Ledit serveur distant 60 comprenant une unité de calcul (non représentée) et une unité de communication (non représentée) permettant audit serveur de communiquer avec la pluralité de dispositifs météorologiques de surveillance 10 ; l’unité de calcul du serveur distant 60 est apte à exécuter le procédé de détermination d’un volume d’eau liquide précipitante au-dessus de chaque position géographique donnée par l’emplacement géographique desdits dispositifs d’alerte météorologique. Les dispositifs 10 sont aptes à émettre une alerte à un dispositif d’émission d‘alertes 40 d’inondation et/ou de crues au niveau de la zone géographique couverte par la pluralité des dispositifs 10 à partir d’une estimation d’un volume V d’eau liquide précipitante dans ladite zone géographique et en comparent ledit volume à un seuil prédéterminé. L’invention ne se limite pas aux modes de réalisation décrits. En particulier on pourra prévoir d’utiliser plusieurs dispositifs météorologiques de surveillance sur différents cours d’eau. On pourra également prévoir l’utilisation de nouvelles constellations de satellites et/ou des nouvelles bandes de fréquences. Procédé pluviométrique d’estimation d’un volume (V) d’eau liquide précipitante au-dessus d’une position géographique (G) prédéterminée comprenant les étapes suivantes : une étape (E1) de réception, par une antenne (11), dite antenne de réception, agencée à ladite position géographique (G) prédéterminée et apte à capter des signaux GNSS (21), d’au moins un signal GNSS (21), dit signal reçu, dont le spectre est composé d’au moins deux fréquences distinctes ; une étape (E2) de détermination, par une unité de calcul (12) apte à recevoir et à traiter des signaux GNSS (21), d’un retard total au zénith (τ) de chaque signal reçu (21) par ladite antenne de réception (11) à ladite position géographique (G) ; une étape (E3) de détermination, par ladite unité de calcul (12), d’un volume (V) d’eau liquide précipitante au-dessus de ladite position géographique prédéterminée (G) à partir uniquement dudit retard total au zénith (τ) déterminé à l’étape précédente et de données atmosphériques relatives à ladite position géographique prédéterminée (G). Procédé selon la revendication 1, caractérisé en ce qu’il comprend une étape préliminaire (E0) de positionnement de l’antenne de réception (11) par une méthode de positionnement absolu. Procédé selon la revendication 2, caractérisé en ce que la méthode de positionnement absolu est une méthode de positionnement ponctuel précis. Procédé selon l’une des revendications 1 à 3, caractérisé en ce que l’étape (E1) de réception d’au moins un signal GNSS comprend en outre les sous-étapes de prétraitement suivantes : une sous-étape (E1.1) d’horodatage de réception, par l’unité de calcul (12), de chaque signal reçu (21) ; une sous-étape (E1.2) de décodage, par l’unité de calcul (12), de l’horodatage d’émission de chaque signal reçu (21) ; une sous-étape (E1.3) de détermination et de correction des retards et/ou des perturbations engendrés par la traversée de l’ionosphère (100) de chaque signal reçu (21) à partir desdites fréquences distinctes composant le spectre de chaque signal reçu (21) ; une sous-étape (E1.4) de détermination, par l’unité de calcul (12), d’au moins un retard total oblique de chaque signal reçu (21) à partir desdits horodatages d’émission et de réception de ce signal reçu (21). Procédé selon la revendication 4, caractérisé en ce que l’étape (E2) de détermination d’un retard total au zénith (τ) de chaque signal reçu (21) à la position géographique prédéterminée (G) s’effectue à partir dudit retard total oblique du signal reçu (21). Procédé selon l’une des revendications 1 à 5, caractérisé en ce que l’étape (E3) de détermination d’un volume (V) d’eau liquide précipitante au-dessus de ladite position géographique (G) comporte les sous-étapes suivantes : une sous-étape (E3.1) de détermination d’une composante sèche dudit retard total au zénith (τ) du signal reçu (21) ; une sous-étape (E3.2) de détermination du volume (V) d’eau liquide précipitante au-dessus de la position géographique (G) à partir dudit retard total au zénith (τ) du signal reçu (21) et de la composante sèche dudit retard total au zénith (τ) du signal reçu (21). Procédé selon la revendication 6, caractérisé en ce que la sous-étape (E3.2) de détermination du volume d’eau liquide précipitante au-dessus de la position géographique (G) à partir dudit retard total au zénith (τ) est effectuée par un algorithme d’apprentissage automatique. Procédé selon la revendication 7, caractérisé en ce que l’algorithme d’apprentissage automatique est un prédicteur scalaire prenant en entrée une valeur de retard total au zénith (τ) d’un signal GNSS (21) ; ledit prédicteur scalaire étant entraîné sur un ensemble de données relatives à des valeurs de retard total au zénith (τ) de signaux GNSS (21), de la position géographique (G) et de données météorologiques. Dispositif météorologique (10) comprenant une antenne de réception (11) apte à capter des signaux GNSS (21) agencée au niveau d’une position géographique (G) prédéterminée et une unité de calcul (12) apte à traiter les signaux GNSS (21) reçus par ladite antenne (11), ledit dispositif météorologique (10) étant caractérisé en ce qu’il met en œuvre le procédé pluviométrique d’estimation d’un volume (V) d’eau liquide précipitante au-dessus d’une position géographique (G) selon l’une des revendications 1 à 8. Système météorologique d’estimation d’un volume d’eau liquide précipitante au-dessus d’une zone géographique comprenant une pluralité de dispositifs météorologiques (10) selon la revendication 9 distribuées au sein de ladite zone géographique. Système d’alerte hydrologique dans une zone géographique comprenant des cours d’eau (50) , ledit système d’alerte comprenant : une pluralité de dispositifs météorologiques (10) de surveillance de cours d’eau (50) déployées au niveau de chaque cours d’eau (50) de ladite zone géographique selon la revendication 9 ; un serveur distant (60) de surveillance météorologique apte à communiquer avec la pluralité de dispositifs météorologiques de surveillance (10); caractérisé en ce que ledit système d’alerte est apte à émettre une alerte de crue, pour l’un quelconque desdits cours d’eau, à partir d’une estimation d’un volume (V) d’eau liquide précipitante dans ladite zone géographique obtenue par le procédé selon l’une des revendications 1 à 8 , lorsque ledit volume d’eau liquide précipitante estimé est supérieur à un seuil prédéterminé.