Brevet FR3054689A1 — publié le 20180202

PROCEDE DE DETERMINATION TEMPORELLE D'UN VOLUME D'EAU PRELEVABLE MAXIMUM ADMISSIBLE D'UNE SOURCE D'EAU SOUTERRAINE [1] L'invention se rapporte à un procédé de détermination temporelle d'un volume d'eau prélevable maximum admissible d'une source d'eau souterraine, ledit volume d'eau étant prélevé en un point de prélèvement et l'état hydrogéologique de la source d'eau souterraine étant qualifié par des mesures de cotes piézométriques sur un piézomètre de référence. Domaine technique [2] De nos jours, les arrêtés de limitation des usages de l'eau traduisent des tensions récurrentes sur les ressources en eau souterraines, et notamment sur les nappes, en période d'étiage. Ces tensions sont exacerbées par la saisonnalité des besoins en eau liés à l'irrigation agricole ou à l'afflux touristique pour l'alimentation en eau potable. Pour un système de production d'eau, ces tensions peuvent engendrer des coûts d'exploitation supplémentaires liés à des achats d'eau, des coûts de pompage ou des coûts de traitement. Le contexte du changement climatique pourrait amplifier ces tensions à moyen terme. [3] Le volume d'eau prélevable d'une ressource en eau souterraine est souvent difficile à estimer en raison de la complexité des hydrosystèmes, la ressource d'eau souterraine pouvant par exemple être une nappe contenue dans un aquifère. État de la technique [04] De nombreux procédés de détermination d'un volume prélevable existent à l'échelle d'une ressource en eau souterraine. Mais la difficulté 3054689 2 consiste à évaluer le volume d'eau réellement disponible sur un ouvrage de prélèvement de la source en eau souterraine, c'est-à-dire au niveau d'un système de production d'eau constitué d'un captage d'eau souterraine (point de prélèvement) ou d'un ensemble de captages (points de prélèvements) 5 formant un champ captant. L'estimation de ce volume prélevable dépend des caractéristiques du captage lui-même et des caractéristiques hydrodynamiques locales de la source d'eau souterraine. [5] L'estimation d'un volume prélevable à l'échelle d'un système de production d'eau mobilisant un ou plusieurs captages d'eau souterraine 10 dépend à la fois du volume d'eau prélevable dans la ressource en eau souterraine et de la configuration des captages eux-mêmes. En théorie, seul un modèle spatialisé 2D ou 3D permettrait d'effectuer ce genre de simulations et de simuler une série de données piézométriques prospectives en un point donné de la source d'eau souterraine. Ce type d'outil est 15 complexe à mettre en oeuvre, mais il existe également des procédés simplifiés graphiques, statistiques et/ou analytiques permettant d'estimer un volume prélevable pour un forage donné. [6] Des procédés actuellement utilisés distinguent plusieurs termes pour définir les limites d'une ressource ou d'un système de production : limite 20 hydrogéologique, limite potentielle, limite de la ressource/système, et « Deployable Output ». [7] Il est à noter que le « Deployable Output » correspond au volume d'eau brute pouvant être mis en distribution après traitement, à partir d'une ou plusieurs ressources, tenant compte, pour une demande en eau déterminée: 25 - des limites hydrogéologiques, - des contraintes physiques et opérationnelles (dues aux capacités physiques des infrastructures de pompage, transport et traitement, aux contraintes hydrauliques en termes de hauteur d'eau à respecter ou de pression minimum/maximum), 3054689 3 - des contraintes de qualité d'eau (périodes associées à une dégradation de la qualité d'eau, ou niveau cible à ne pas dénoyer pour des aquifères, etc.), - des contraintes réglementaires (autorisations de prélèvement). 5 Il est à noter également que chacun de ces paramètres peut constituer un facteur limitant pour le « Deployable Output ». Une fois que le « Deployable Output » a été déterminé, les pertes en eau et exports sont soustraits et les imports ajoutés, afin d'obtenir les volumes disponibles en distribution. Toutefois, cette approche n'est pas satisfaisante car elle ne fournit qu'une 10 valeur unique et très conservative du volume maximum prélevable. [08] D'autres procédés se basent sur une estimation de la limite exploitable à long-terme d'un forage. Par exemple, pour un aquifère captif, des procédés analytiques permettent de calculer une limite d'exploitation sécuritaire théorique. Cette limite représente généralement le débit pouvant 15 être maintenu pendant une longue durée, sans franchir le rabattement disponible. Il est à noter que le rabattement disponible correspond à la différence entre le niveau d'eau en l'absence de pompage et la cote du toit de l'aquifère captif. Cette limite est alors considérée comme un indicateur plutôt que comme une valeur limite fiable. 20 [09] Toutefois, les procédés actuellement utilisés et dédiés à l'estimation du volume maximum prélevable sur les points de prélèvement d'eau souterraine sur une longue durée fournissent des estimations très sécuritaires, car les captages ne sont généralement pas exploités en continu et les aquifères se rechargent en parallèle. L'application de ces procédés 25 sur des durées aussi longues, par exemple vingt ans, est peu réaliste, car les conditions d'applicabilité sur de telles durées ne sont plus satisfaites. De plus, ces procédés font l'hypothèse d'un prélèvement à débit constant tout au long de cette longue durée, ce qui est également peu réaliste. En d'autres termes, les procédés actuellement utilisés ne répondent pas aux objectifs 30 recherchés, car ils fournissent une valeur unique de volume maximum prélevable, non actualisée en fonction de la situation hydrogéologique. 3054689 4 [10] Il existe donc un réel besoin de fournir un procédé palliant ces défauts, inconvénients et obstacles de l'art antérieur. En particulier, il existe un besoin de fournir un procédé permettant de répondre à l'ensemble des attentes, notamment : 5 - l'évaluation d'un volume d'eau souterraine prélevable réaliste et pas simplement sécuritaire ; - l'évaluation d'un volume d'eau souterraine prélevable variable dans le temps selon les différentes contraintes ; - une utilisation simple, applicable à un grand nombre de captages. 10 Description de l'invention [11] Pour résoudre un ou plusieurs des inconvénients cités précédemment, 15 l'invention a pour objet un procédé de détermination temporelle d'un volume d'eau prélevable maximum admissible d'une source d'eau souterraine, ledit volume d'eau étant prélevé en un point de prélèvement et l'état hydrogéologique de la source d'eau souterraine étant qualifié par des mesures de cotes piézométriques sur un piézomètre de référence, 20 ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes : a) mesure en continu par un premier capteur de niveau de cotes piézométriques sur le point de prélèvement, ledit capteur disposant d'un premier historique de données disponible sur une durée prédéterminée révolue ; 25 b) mesure en continu par un deuxième capteur de niveau de cotes piézométriques sur le piézomètre de référence, ledit deuxième capteur disposant d'un deuxième historique de données disponible sur la durée prédéterminée révolue ; ledit procédé comprenant en outre les étapes suivantes mises en oeuvre par 30 une machine de calcul : c) traitement des données du premier capteur de niveau pour créer une chronique de cotes piézométriques pseudo-statiques et une 3054689 5 chronique de cotes piézométriques dynamiques sur la durée prédéterminée révolue ; d) détermination d'une chronique de rabattements au point de prélèvement sur ladite durée prédéterminée révolue, un rabattement 5 étant défini comme étant la différence, à un instant donné, entre la cote piézométrique pseudo-statique et la cote piézométrique dynamique ; e) détermination de valeurs de transmissivité apparentes au point de prélèvement, permettant de reproduire au mieux la chronique de 10 rabattements déterminés sur ladite durée prédéterminée révolue, en utilisant la relation analytique de Cooper-Jacob afin d'associer une valeur de transmissivité apparente à différentes classes de cotes piézométriques pseudo-statiques, l'association constituant une première relation ; 15 f) détermination d'un niveau critique d'exploitation audit point de prélèvement ; g) sélection d'un piézomètre de référence disposant de mesures de cotes piézométriques sur ladite durée prédéterminée révolue ; h) calcul de moyennes temporelles desdites cotes piézométriques 20 pseudo-statiques et des cotes piézométriques mesurées sur le piézomètre de référence sur ladite durée prédéterminée révolue ; i) détermination d'une deuxième relation entre lesdites moyennes temporelles desdites cotes piézométriques pseudo-statiques et les moyennes temporelles desdites cotes piézométriques mesurées sur 25 le piézomètre de référence sur ladite durée prédéterminée révolue ; j) détermination d'un rabattement maximum admissible pour chaque valeur de cote piézométrique pseudo-statique, un rabattement maximum admissible étant défini comme étant la différence entre une cote piézométrique pseudo-statique et ledit niveau critique 30 d'exploitation ; 3054689 6 k) détermination du volume prélevable maximum admissible de ladite source d'eau souterraine prélevé audit point de prélèvement en utilisant la relation de Cooper-Jacob et lesdites première et deuxième relations. 5 [12] La ressource en eau peut être une ressource en eau brute souterraine, comme des sources, des nappes, des karsts, exploitée par un ou plusieurs captages. De préférence, la ressource en eau est une nappe dont le volume d'eau prélevable est déterminé par le niveau piézométrique. Il est à noter que cette ressource en eau peut être reliée à une ou plusieurs entités de 10 production d'eau telles que des usines de traitement. [13] Au sens de la présente invention, on entend par niveau critique d'exploitation, le niveau en dessous duquel le captage d'eau altère la bonne régénérescence en eau de la ressource et/ou altère le matériel de captage. [14] Avantageusement, les moyennes temporelles sont des moyennes 15 mensuelles. [15] L'invention a également pour objet un programme informatique comprenant des instructions adaptées à la mise en oeuvre de chacune des étapes du procédé ci-avant décrit lorsque ledit programme est exécuté sur une machine de calcul. 2 0 [16] L'invention a encore pour objet un système comprenant des moyens adaptés à la mise en oeuvre de chacune des étapes du procédé décrit ci-avant. [17] L'invention sera mieux comprise à la lecture de la description qui suit, faite uniquement à titre d'exemple et fournie à titre indicatif et nullement 25 limitatif. Description détaillée [18] Selon un mode de réalisation de l'invention, un procédé de 30 détermination temporelle d'un volume d'eau prélevable maximum admissible est appliqué à une nappe. Le volume d'eau est alors prélevé en 3054689 7 un point de prélèvement qui peut par exemple être équipé d'un dispositif de captage quelconque utilisé et connu de l'état de l'art. L'état hydrogéologique de la nappe est alors qualifié par des mesures de cotes piézométriques sur un piézomètre de référence. 5 [19] La ressource en eau étant une nappe, le volume d'eau prélevable qu'elle contient est alors déterminé par niveau piézométrique. Dans cette configuration, le procédé de détermination par le niveau piézométrique permet de prendre indirectement en compte la situation hydrogéologique de la nappe, et les phénomènes naturels et anthropiques ayant un effet 10 sur l'état quantitatif de l'aquifère. Cette méthode est simple, robuste et applicable à des aquifères continus c'est-à-dire non karstiques et non fracturés, et pour de multiples configurations de dispositifs de captage comme les forages, les puits à drains, et les champs captants, par exemple, à l'exception des captages de sources. 15 [20] Il est à noter que la gestion quantitative des eaux contenues dans la nappe et la définition de volumes prélevables au sein de cette nappe nécessite de définir au préalable des limites spatiales de l'ensemble considéré. Ces limites sont définies en fonction des caractéristiques géologiques et hydrogéologiques locales. 20 [21] Par exemple, on dispose généralement de données géologiques et hydrogéologiques de la piézométrie régionale, de la transmissivité ou encore du coefficient d'emmagasinement, et de données d'exploitation (par exemple des mesures de niveau d'eau et une cote du repère utilisé, des volumes prélevés, des débits nominaux d'une ou plusieurs des pompes). 25 [22] La définition de volumes prélevables de la nappe nécessite également de disposer d'une bonne connaissance de son état hydrogéologique, d'un bilan hydraulique précis de toutes les entrées et sorties d'eau du système (naturelles ou anthropiques), et d'évaluer sa capacité de stockage intrinsèque. 30 [23] Par exemple, on dispose comme caractéristiques générales du captage d'exploitation, des coupes géologiques qui sont des 3054689 8 représentations schématiques de la succession des formations géologiques sur un profil vertical, avec indication des cotes associées, ou encore des coupes techniques qui sont des représentations schématiques des caractéristiques et des cotes des équipements d'un ouvrage souterrain, par 5 exemple, sur un profil vertical. On peut citer par exemple, des trous nus, les caractéristiques du pré-tubage et du tubage, massif filtrant, cimentation, protection en tête. [24] La gestion quantitative d'une nappe nécessite donc une connaissance de son état, de sa recharge, des prélèvements et des besoins. 10 [25] D'autres données d'exploitation comme le nombre d'heures journalières de pompage et le débit instantané moyen journalier peuvent facultativement être également prises en compte. [26] Le procédé selon l'invention comprend une étape a) consistant à mesurer en continu, par un premier capteur de niveau les cotes 15 piézométriques sur le point de prélèvement, ledit capteur disposant d'un premier historique de données disponibles sur une durée prédéterminée révolue. Il est nécessaire de disposer de valeurs de cotes piézométriques de la nappe s'étendant sur une durée prédéterminée, qui peut par exemple être d'au moins deux ans. Ce premier capteur peut fonctionner avec d'autres 20 capteurs afin de permettre ultérieurement de mettre en avant des cotes pseudo-statiques, définissant l'état de la nappe sans l'influence du prélèvement au point de prélèvement, ainsi que des cotes dynamiques définissant l'état de la nappe en phase de pompage au point de prélèvement. Ces cotes sont des valeurs d'entrée nécessaires à la mise en 25 oeuvre du procédé selon l'invention. [27] Le procédé comprend également une étape b), réalisée après l'étape a), consistant à mesurer en continu par un deuxième capteur de niveau des cotes piézométriques sur le piézomètre de référence, ledit deuxième capteur disposant d'un deuxième historique de données disponible sur la durée 30 prédéterminée révolue. Cette étape permettra de sélectionner ultérieurement le piézomètre de référence adéquat associé à la nappe dans 3054689 9 laquelle le volume d'eau est prélevé. En d'autres termes, cette étape permet de lier le comportement de la nappe au point de prélèvement avec le comportement de la nappe au piézomètre de référence. Cette liaison est possible, comme on le verra dans la suite de cette description, en comparant 5 les moyennes temporelles des cotes piézométriques pseudo-statiques au point de prélèvement et les moyennes temporelles des cotes piézométriques au piézomètre de référence sur la même durée prédéterminée révolue. [28] Après l'étape b) , le procédé comprend une étape c) mise en oeuvre 10 par une machine de calcul et consistant à traiter des données du premier capteur de niveau pour créer une chronique de cotes piézométriques pseudo-statiques et une chronique de cotes piézométriques dynamiques sur la durée prédéterminée révolue. [29] Pour ce faire, des données d'exploitation journalières sur le captage 15 peuvent être récupérées sur une période d'au moins deux ans. Ces données d'exploitation peuvent, par exemple, renseigner sur : - la cote pseudo-statique (valeur de profondeur journalière maximum, en m) de la nappe étudiée ; - la cote dynamique (valeur de cote journalière minimum, en m) ; 20 - la durée journalière de pompage (en h) ; - le débit instantané moyen journalier (en m3/h) ; - le volume prélevé journalier (en m3/j). [30] Ensuite, une étape d) est réalisée et consiste à déterminer une chronique de rabattements au point de prélèvement sur ladite durée 25 prédéterminée révolue, un rabattement étant défini comme étant la différence, à un instant donné, entre une ladite cote piézométrique pseudostatique et ladite cote piézométrique dynamique. [31] Après l'étape d), une étape e) est appliquée et consiste à déterminer des valeurs de transmissivité apparentes au point de prélèvement, 30 permettant de reproduire au mieux la chronique de rabattements sur ladite durée prédéterminée révolue, en utilisant la relation analytique de Cooper- 3054689 10 Jacob afin d'associer une valeur de transmissivité apparente à différentes classes de cotes piézométriques pseudo-statiques, l'association constituant une première relation. [32] Cette première relation peut être déterminée par la relation analytique r2 S . 2.303Q (2.25Tt 5 de Cooper-Jacob suivante : s(r,t) = 4-irT log ) Celle-ci est couramment utilisée en hydrogéologie quantitative. Elle permet notamment de fournir une estimation du volume maximum d'eau souterraine prélevable en fonction du rabattement, lui-même dépendant du niveau piézométrique statique de la nappe. 10 [33] En particulier, cette relation de Cooper-Jacob est utilisée afin de calculer un rabattement théorique s(r,t), en fixant les paramètres nécessaire à l'application de cette relation, c'est-à-dire la transmissivité (T), le coefficient d'emmagasinement (S), la distance radiale au point de prélèvement (r), la durée du prélèvement journalier (t), et le débit instantané 15 moyen journalier (Q). [34] On peut donc ensuite représenter des couples de points rabattement / cote piézométriques journalières pseudo-statique sur un graphique pour des valeurs proches de débit journalier moyen. Le graphique permet de segmenter visuellement les valeurs de rabattement observées selon 20 plusieurs classes de cotes piézométriques journalières pseudo-statiques. Pour chaque classe, les valeurs de transmissivité apparente peuvent être ajustées afin que la valeur de l'indicateur d'écart (communément désigné par RMSE) entre rabattement observé et rabattement simulé avec la relation de Cooper-Jacob soit minimum. 25 [35] Dans une étape f) suivante, on détermine un niveau critique d'exploitation audit point de prélèvement. Il est à noter que le niveau critique d'exploitation d'un captage tel qu'un puits ou un forage peut être conditionné par : - la cote altimétrique supérieure de la partie crépinée du tubage ; 30 - la cote d'exploitation locale limite de la nappe, qui peut être : o une cote de maintien de captivité de la nappe, 3054689 11 o une cote de non-invasion de biseau salé, o une cote de non-dénoyage d'une zone productive, o une cote réglementaire (piézométrie objectif, seuil d'alerte...), - la cote de la crépine d'aspiration de la pompe ou la cote de sécurité 5 déclenchant l'arrêt de la pompe. [36] La cote piézométrique la plus élevée, donc la plus défavorable, puisqu'elle minimisera le rabattement maximum admissible et donc le volume prélevable, sera retenue comme niveau critique d'exploitation de l'ouvrage, notée zNc. La détermination du niveau critique nécessite donc de 10 disposer de la coupe technique et de la coupe géologique de l'ouvrage, et de l'existence éventuelle de cotes réglementaires de gestion de la nappe. Il est à noter que dans le cas d'un champ captant, le niveau critique le plus défavorable sera retenu et appliqué à un captage conceptuel qui représentera le champ captant par un unique point de prélèvement. Il est à 15 noter que le niveau critique d'exploitation peut être fixé antérieurement et représenter un seuil d'alerte. [37] Ensuite, on procède à une étape g) de sélection du piézomètre de référence. Il est à noter que la sélection du piézomètre de référence, qui est généralement un piézomètre régional de référence, nécessite de lister 20 l'ensemble des piézomètres captant la masse d'eau souterraine de la nappe étudiée. Ces piézomètres peuvent être identifiés facilement en consultant des bases de données, et en effectuant une recherche des stations de suivi du niveau d'eau par masse d'eau. Parmi ces piézomètres, seuls ceux qui sont en activité et qui disposent d'un historique suffisant sur une durée 25 prédéterminée révolue (idéalement supérieur à 10 ans) seront considérés. Cet historique est connu suite aux mesures acquises par le deuxième capteur. [38] Ensuite, dans une étape h), on calcule des moyennes temporelles desdites cotes piézométriques pseudo-statiques et des cotes 30 piézométriques mesurées sur les piézomètres de référence listés sur ladite durée prédéterminée révolue. 3054689 12 [39] De préférence, les moyennes temporelles peuvent être des moyennes mensuelles. Le piézomètre de référence est retenu de la manière suivante : minimisation de la RMSE entre moyenne mensuelle de la cote piézométrique pseudo-statique au captage et la translation de la moyenne 5 mensuelle de la cote piézométrique statique mesurée au piézomètre sur la durée d'observation commune. [40] Cette translation (notée hspt) est obtenue en ajoutant à chaque valeur de la chronique mensuelle statique (lisp) du piézomètre l'écart relatif entre la moyenne des cotes mensuelles pseudo-statiques au captage (fisc) et la 10 moyenne des cotes mensuelles statiques au piézomètre (lisp) comme indiqué ci-après : hspt(t) = hsp(t) + (h' - lisp) [41] Dans le cas de plusieurs piézomètres avec une RMSE proche, celui 15 qui dispose de l'historique le plus long sera idéalement retenu. [42] Ensuite, dans une autre étape i), on détermine une deuxième relation entre lesdites moyennes temporelles desdites cotes piézométriques pseudo-statiques et lesdites cotes piézométriques mesurées sur le piézomètre de référence sur ladite durée prédéterminée révolue. Cette 20 étape permet ensuite, dans le cas de scénarios prospectifs sur l'état hydrogéologique de la nappe, de déterminer une moyenne temporelle de la cote piézométrique pseudo-statique au captage à partir d'une moyenne temporelle de la cote piézométrique au piézomètre de référence en utilisant ladite deuxième relation. 25 [43] La deuxième relation, le plus souvent linéaire, ou constituée de plusieurs segments linéaires, sera retenue et permettra d'exprimer la cote piézométrique pseudo-statique mensuelle moyenne théorique au captage en fonction de la cote piézométrique mensuelle moyenne au piézomètre de référence. Cette relation empirique pourra être obtenue en utilisant une ou 3054689 13 plusieurs régressions linéaires sur un graphique de type scatter-plot, ou d'autres fonctions de corrélation. [44] Après cela, une étape j) est réalisée et consiste à déterminer le rabattement maximum admissible pour chaque valeur de cote 5 piézométrique pseudo-statique, un rabattement maximum admissible étant défini comme étant la différence entre une cote piézométrique pseudostatique et ledit niveau critique d'exploitation. Il est à noter que le rabattement maximum admissible est considéré comme étant la différence entre la cote piézométrique pseudo-statique et le niveau critique 10 d'exploitation, c'est-à-dire le seuil d'alerte. Un rabattement maximum admissible plus élevé autorisera un volume prélevable plus grand. Le rabattement maximum admissible, Smax, variable au cours du temps t, est défini comme la différence entre le niveau piézométrique pseudo-statique hp, et le niveau critique d'exploitation, noté zNc. 15 Smax(t) = hp,(t) - z' [45] Le rabattement maximum admissible tient compte indirectement des prélèvements et des phénomènes de recharge et décharge naturels de la nappe. Ces derniers influencent le niveau pseudo-statique. [46] Puis finalement, une étape k) est réalisée et consiste à déterminer le 2 0 volume prélevable maximum admissible de ladite source d'eau souterraine audit point de prélèvement en utilisant la relation de Cooper-Jacob et lesdites première et deuxième relations. [47] Le volume maximum prélevable Vmax est calculé, à chaque pas de temps de la simulation, à l'aide de la relation de Cooper-Jacob (1946) et de 25 la relation T = f(hps) obtenue lors de la précédente étape. Il dépend du rabattement maximum admissible sm', 3054689 14 Vmax Qmax- texp 4TtT(hps)smax texp 2,303 log (2.25T(hps)texp) r2S texp étant la durée maximum d'exploitation. [48] Dans le calcul du volume maximum prélevable, la valeur de texp est fixée par défaut à 20 heures/jour. Dans le cas d'un champ captant, le volume maximum prélevable sera le volume prélevable par l'ensemble des captages 5 qu'il contient. [49] Ainsi, ce procédé prend en compte indirectement la situation hydrogéologique, les effets de recharge/décharge de la nappe et les phénomènes naturels, qui influencent les cotes pseudo-statiques mesurées au captage. Elle prend également en compte indirectement les effets des 10 prélèvements globaux sur la nappe et les activités anthropiques, qui influencent les cotes pseudo-statiques mesurées au captage. [50] Ainsi, selon cette méthode préférentielle selon l'invention, il est possible de réaliser des prévisions de disponibilité future de la ressource en eau, en particulier de la nappe, basées sur : 15 - les tendances hydriques connues au jour le jour, - l'intégration comme variable du changement climatique dans les scénarios à moyen et long terme. [51] Le suivi quantitatif des eaux souterraines est surveillé, selon le mode préférentiel de l'invention, en surveillant les variations de niveau 20 piézométrique des nappes. Toutefois ce suivi quantitatif peut également être surveillé en mesurant le débit des émergences (sources) selon la nature du captage. [52] Il est à noter que ce procédé permet de calculer un volume maximum prélevable théorique sur la base d'une cote piézométrique pseudo-statique 25 et d'un niveau critique d'exploitation sur un captage. Ce volume maximum prélevable théorique n'est pas nécessairement atteignable en conditions d'exploitation. 3054689 15 Programme informatique [53] Il est à noter qu'un programme d'ordinateur comprenant des instructions adaptées à la mise en oeuvre de chacune des étapes du 5 procédé ci-dessus décrites peut être élaboré. Ainsi, une machine de calcul est apte à exécuter ce programme d'ordinateur afin de gagner en efficacité et en rapidité. Système 10 [54] Par ailleurs, un système comprenant des moyens adaptés à la mise en oeuvre de chacune des étapes précitées peut être réalisé. [55] L'invention a été illustrée et décrite en détail dans la description précédente. Celle-ci doit être considérée comme illustrative et donnée à titre 15 d'exemple et non comme limitant l'invention à cette seule description. De nombreuses variantes de réalisation sont possibles. [56] Dans les revendications, le mot « comprenant » n'exclut pas d'autres éléments et l'article indéfini « un/une » n'exclut pas une pluralité. REVENDICATIONS1. Procédé de détermination temporelle d'un volume d'eau prélevable maximum admissible d'une source d'eau souterraine, ledit volume d'eau étant prélevé en un point de prélèvement et l'état hydrogéologique de la source d'eau souterraine étant qualifié par des mesures de cotes piézométriques sur un piézomètre de référence, ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes : a) mesure en continu par un premier capteur de niveau de cotes piézométriques sur le point de prélèvement, ledit capteur disposant d'un premier historique de données disponible sur une durée prédéterminée révolue ; b) mesure en continu par un deuxième capteur de niveau de cotes piézométriques sur le piézomètre de référence, ledit deuxième capteur disposant d'un deuxième historique de données disponible sur la durée prédéterminée révolue ; ledit procédé comprenant en outre les étapes suivantes mises en oeuvre par une machine de calcul : c) traitement des données du premier capteur de niveau pour créer une chronique de cotes piézométriques pseudo-statiques et une chronique de cotes piézométriques dynamiques sur la durée prédéterminée révolue ; d) détermination d'une chronique de rabattements au point de prélèvement sur ladite durée prédéterminée révolue, un rabattement étant défini comme étant la différence, à un instant donné, entre la cote piézométrique pseudo-statique et la cote piézométrique dynamique ; e) détermination de valeurs de transmissivité apparentes au point de prélèvement, permettant de reproduire au mieux la chronique de rabattements déterminés sur ladite durée prédéterminée révolue, en 3054689 17 utilisant la relation analytique de Cooper-Jacob afin d'associer une valeur de transmissivité apparente à différentes classes desdites cotes piézométriques pseudo-statiques, l'association constituant une première relation ; 5 f) détermination d'un niveau critique d'exploitation audit point de prélèvement ; g) sélection d'un piézomètre de référence disposant de mesures de cotes piézométriques sur ladite durée prédéterminée révolue ; h) calcul de moyennes temporelles desdites cotes piézométriques 10 pseudo-statiques et des cotes piézométriques mesurées sur le piézomètre de référence sur ladite durée prédéterminée révolue ; i) détermination d'une deuxième relation entre lesdites moyennes temporelles desdites cotes piézométriques pseudo-statiques et les moyennes temporelles desdites cotes piézométriques mesurées sur 15 le piézomètre de référence sur ladite durée prédéterminée révolue ; j) détermination d'un rabattement maximum admissible pour chaque valeur de cote piézométrique pseudo-statique, un rabattement maximum admissible étant défini comme étant la différence entre une cote piézométrique pseudo-statique et ledit niveau critique 20 d'exploitation ; k) détermination du volume prélevable maximum admissible de ladite source d'eau souterraine prélevé audit point de prélèvement en utilisant la relation de Cooper-Jacob et lesdites première et deuxième relations. 25 2. Procédé selon la revendication 1, selon lequel lesdites moyennes temporelles sont des moyennes mensuelles. 3. Programme informatique comprenant des instructions adaptées à la 30 mise en oeuvre de chacune des étapes du procédé selon les revendications 1 ou 2 lorsque ledit programme est exécuté sur une machine de calcul. 3054689 18 4. Système comprenant des moyens adaptés à la mise en oeuvre de chacune des étapes du procédé selon les revendications 1 ou 2.