FR 3054705 A1 20180202 FR 1657370 A 20160729 1 OUTIL DE GESTION DE MULTIPLES RESSOURCES EN EAU [1] L'invention se rapporte à un système et un procédé permettant d'optimiser l'allocation de plusieurs ressources en eau pour satisfaire la demande en eau produite par au moins deux entités de production d'eau interconnectées, en réduisant la consommation d'énergie et en renforçant l'autonomie en eau du système. Domaine technique [2] De nos jours, des tensions hydriques existent, avec une intensité variable selon les zones géographiques, et notamment selon les ressources hydriques locales. Il confirme le besoin d'outils pour la gestion quantitative de ressources en eau. [3] Ainsi, la présente invention se rapporte à un outil de gestion quantitative, ou outil d'aide à l'exploitation, dédié aux systèmes de production d'eau interconnectés, c'est-à-dire aux systèmes comprenant : - plusieurs entités de production, chacune s'alimentant dans une ou plusieurs ressources, - des éléments demandeurs en eau, et - des ressources en eau. [04] Ces entités de production sont interconnectées dans le but de transférer de l'eau brute ou de l'eau traitée entre-elles. Cet outil d'aide à l'exploitation permet d'assurer la continuité du service de production d'eau en période de stress hydrique ou non, en limitant notamment la consommation d'énergie liée à l'exploitation des systèmes de production d'eau, en maintenant au mieux l'équilibre ressource/demande, et en renforçant l'autonomie du système en utilisant autant que faire se peut l'eau du système. 3054705 2 [5] Il est à noter qu'au sens de la présente demande, on entend par eau traitée toute eau ayant subie un traitement avant sa distribution, et notamment avant sa transmission aux éléments demandeurs d'eau. Cette eau traitée peut être potable ou non, et dédiée à des applications agricoles, 5 industrielles ou relatives à des collectivités. [6] Il est à noter qu'au sens de la présente demande, on entend par eau brute toute eau ayant été prélevée d'une ressource en eau, que celle-ci soit d'ordre naturel ou artificiel. Les eaux brutes ou ressources en eau peuvent être de tous types tels que des cours d'eau, des plans d'eau, des eaux 10 souterraines, etc... [7] Au sens de la présente demande, on entend par stress hydrique un niveau de la ressource en eau inférieur au niveau statistique moyen. [8] Au sens de la présente demande, on entend par autonomie un système capable de satisfaire à une partie au moins de la demande en eau 15 sans importer de l'eau provenant d'un élément, dispositif, ou système extérieur au système selon l'invention. État de la technique 2 0 [09] De nos jours, de multiples outils informatiques dédiés à la gestion quantitative de ressources en eau sont utilisés et peuvent être répertoriés en quatre grandes catégories : les Outils d'Aide à la Décision (OAD), les plateformes génériques d'interfaçage et de gestion de données, les environnements de développement avec des librairies spécifiques 25 associées, et les logiciels de modélisation autonomes. [10] Toutefois, aucun des systèmes actuellement utilisés ne permet de répondre à l'ensemble des attentes. En effet, aucun d'eux ne permet notamment de simuler des périodes de stress hydrique prospectives ou de nouvelles interconnexions entre les systèmes de distribution permettant 30 d'allouer de l'eau à l'ensemble des éléments demandeurs. Par système de distribution, on entend au sens de la présente invention l'ensemble 3054705 3 ressources en eau, entités de production et éléments demandeurs reliés à l'une au moins des entités de production du système de distribution. Il est à noter également que les systèmes actuellement utilisés ne déterminent les priorités de transfert et d'allocation d'eau que selon le critère de 5 consommation d'énergie liée à la production d'eau (optimisation mono- objectif), et non selon les critères de consommation d'énergie liée à la production d'eau et d'équilibre entre les ressources en eau et la demande (optimisation multi-objectifs). Par ailleurs, actuellement, les eaux souterraines (ou aquifères) ne peuvent être des sources à considérer par 10 les systèmes actuellement utilisés qu'à condition qu'un modèle de simulation complexe soit mis en place (modélisation 2D ou 3D des aquifères). [11] Il existe donc un réel besoin de fournir un système de gestion quantitative des ressources en eau palliant ces défauts, inconvénients et 15 obstacles de l'art antérieur. En particulier, il existe un besoin de fournir un système permettant de simuler des stratégies de production d'eau, ces stratégies étant optimisées selon l'équilibre entre les ressources en eau et la demande, en optimisant la consommation d'énergie et en rationalisant les transferts d'eau brute et/ou d'eau traitée entre les éléments interconnectés 20 du système, le but étant de réduire les consommations d'énergie liées aux transferts des volumes d'eau des entités de production aux éléments demandeurs, de façon à fournir l'eau la moins consommatrice d'énergie possible. [12] Par stratégie, on entend au sens de la présente invention, un ensemble 25 d'états définissant chacun des éléments du système et en particulier chacune des entités de production et ressources en eau du système. Description de l'invention 30 3054705 4 [13] Pour résoudre un ou plusieurs des inconvénients cités précédemment, l'invention a pour objet un système de gestion quantitative d'eau comprenant les éléments suivants : - au moins deux entités de production d'eau interconnectées, 5 directement ou indirectement, chaque entité de production étant définie par un débit maximal de production, - au moins une ressource en eau alimentant au moins une entité de production, chaque ressource en eau étant définie par un débit maximal de prélèvement, 10 - au moins un élément demandeur d'eau produite définie par une courbe temporelle préétablie de demande en eau produite en fonction du temps, chaque liaison entre entités de production, ressources en eau et éléments demandeurs étant assurée par un ouvrage de transfert ayant un débit de 15 transfert maximal et/ou minimal prédéterminés, ces ouvrages de transfert pouvant être interconnectés, et chaque entité de production et chaque ressource en eau étant en outre associée à une fonction de pondération P. [14] Selon l'invention, les ressources en eau peuvent être de tous types tels que des cours d'eau, des plans d'eau, des eaux souterraines, etc... 20 [15] Il est à noter que les entités de production d'eau sont dites interconnectées, directement ou indirectement, car celles-ci peuvent être directement reliées entre-elles ou peuvent être reliées entre-elles indirectement en étant connectées à un élément commun. [16] Avec cet objectif en vue, le système selon l'invention, par ailleurs 25 conforme au préambule cité ci-avant, est caractérisé en ce qu'il comprend un premier calculateur adapté pour minimiser la fonction de pondération globale Pg du système, somme de toutes les fonctions de pondération des différents éléments du système, tout en garantissant le respect de la courbe temporelle préétablie de demande en eau produite de chaque élément 30 demandeur et sous contrainte de respect des débits maximaux et/ou des débits minimaux des différents éléments du système. 3054705 5 [17] On entend par fonction de pondération P d'un élément du système, au sens de la présente invention, une fonction de pondération mathématique définie comme étant un coefficient multiplicateur du flux d'eau circulant au travers de l'élément en question. Ce coefficient multiplicateur est défini 5 comme étant une valeur qui peut être fixée arbitrairement, ou à partir de la consommation énergétique associée à cet élément. La valeur de ce coefficient multiplicateur dépend de valeurs constantes et de valeurs variables. [18] Par valeurs constantes, au sens de la présente invention, on entend 10 des valeurs arbitraires, définies par l'utilisateur, priorisant l'allocation d'eau dans les éléments du système. Les valeurs constantes représentent des contraintes qui contraignent le système et limitent les transferts d'eau entre les différents éléments du système. Ces valeurs constantes sont caractéristiques de chacun ou d'un ensemble d'éléments du système. 15 [19] Par contraintes, au sens de la présente invention, on entend les capacités minimales et/ou maximales de transfert entre les éléments du système, les capacités de production d'eau des entités de production et des ressources en eau, les capacités de stockage des entités de production et des ressources en eau, ou encore les débits de restitution en aval des 20 ressources en eau. Il est à noter que ces contraintes peuvent être d'ordre technique ou réglementaire. [20] Par valeurs variables, au sens de la présente invention, on entend des valeurs qui sont fonction des allocations de flux d'eau dans le système. [21] On entend par ouvrage de transfert, au sens de la présente invention, 25 toute construction apte à transférer de l'eau entre deux points relativement éloignés. [22] Selon l'invention, une fonction de pondération P est associée à chaque entité de production et à chaque ressource en eau du système. Pour une stratégie donnée, à chaque ressource en eau et chaque entité de production 30 du système correspond une fonction de pondération. 3054705 6 [23] Au sens de la présente invention, on entend par fonction de pondération globale Pg, la somme des différentes fonctions de pondération P des différentes entités de production et ressources en eau du système. [24] Selon l'invention, la fonction de pondération globale Pg doit être 5 minimisée, tout en tenant compte de la disponibilité des ressources en eau. [25] Ainsi, pour minimiser les fonctions de pondération P de chaque ressource en eau et de chaque entité de production du système, et donc de minimiser la fonction de pondération globale Pg, il est nécessaire de minimiser les valeurs variables qui définissent chacune des fonctions de 10 pondération. [26] Le premier calculateur a pour objectif de générer une stratégie de production déterminée selon l'objectif de minimisation des fonctions de pondération P de chaque ressource en eau et de chaque entité de production pour une satisfaction de la demande de 100%. Le premier 15 calculateur réalise une optimisation mono objectif contraignant l'allocation en eau et la répartition des flux d'eau dans le système pour satisfaire la demande en eau produite. L'algorithme d'optimisation mono objectif génère des résultats définissant la stratégie de production. [27] Par algorithme d'optimisation mono objectif, au sens de la présente 20 invention, on entend une suite de règles opératoires appliquées dans un ordre déterminé, à un nombre fini de données, sur les éléments du système afin de déterminer, en un nombre fini d'étapes, la stratégie optimale répondant à l'objectif de satisfaction de la demande en eau produite, tout en tenant compte de l'hydrologie des ressources mobilisées par le système et 25 en respectant les contraintes du système. [28] Avantageusement, la satisfaction de la courbe temporelle préétablie de demande en eau produite de chaque élément demandeur d'eau produite est en outre sous contrainte de respect des volumes d'eau prélevables dans les ressources en eau, des volumes d'eau produite par les entités de 30 production et des contraintes du système. 3054705 7 [29] Avantageusement, les différents ouvrages de transfert comprennent des vannes commandées par le premier calculateur de sorte que le système délivre aux éléments demandeurs, en temps réel, de l'eau produite selon la courbe temporelle préétablie en minimisant la fonction de pondération 5 globale Pg, tout en respectant les contraintes du système. [30] Avantageusement, l'entité de production associée à la fonction de pondération P la plus faible transmet de l'eau produite en priorité à l'entité de production d'eau associée à la fonction de pondération P la plus élevée et/ou transmet de l'eau produite en priorité à l'élément demandeur associé 10 à la fonction de pondération P la plus élevée. [31] Avantageusement, une entité de production reliée à plusieurs ressources en eau prélève en priorité un volume d'eau de la source associée à la plus faible fonction de pondération P. [32] Avantageusement, le système selon l'invention comprend en outre un 15 second calculateur adapté pour définir des stratégies de production pour satisfaire la demande en minimisant la consommation énergétique, chaque stratégie correspondant à un pourcentage de la courbe temporelle préétablie de l'élément demandeur, et ledit second calculateur étant connecté au premier calculateur de sorte que ledit premier calculateur 20 génère une valeur minimale de la fonction de pondération globale Pg pour chaque stratégie de production, ledit second calculateur étant en outre adapté pour déterminer la meilleure stratégie définie comme étant la stratégie fournissant le pourcentage maximal de satisfaction de la demande pour un ratio variation de la fonction de pondération globale Pg sur variation 25 de la demande minimisé. Le second calculateur réalise une optimisation multi objectifs contraignant l'allocation en eau et la répartition des flux d'eau dans le système pour satisfaire la demande en eau produite et pour minimiser la consommation énergétique. L'algorithme d'optimisation multi objectifs génère des résultats définissant les stratégies de production. 30 [33] Par algorithme d'optimisation multi objectifs, au sens de la présente invention, on entend une suite de règles opératoires appliquées dans un 3054705 8 ordre déterminé, à un nombre fini de données, sur les éléments du système afin de déterminer, en un nombre fini d'étapes, l'ensemble des stratégies optimales répondant à l'objectif de satisfaction de la demande et à l'objectif de minimisation de la consommation énergétique du système, tout en tenant 5 compte de l'hydrologie des ressources mobilisées par le système et en respectant les contraintes du système. [34] Ce second calculateur a pour objectif de générer différents scénarios définissant différents pourcentages de satisfaction de la demande. Ce second calculateur, en collaboration avec le premier calculateur, a 10 également pour but d'associer à chaque pourcentage de satisfaction de la demande, une fonction de pondération globale Pg. Cette fonction de pondération globale Pg correspond à la plus faible fonction de pondération globale Pg satisfaisant au scénario sélectionné. Cette plus faible fonction de pondération globale Pg est définie par la somme des fonctions de 15 pondération P des entités de production et des ressources en eau, et cette fonction de pondération globale Pg représente alors la meilleure stratégie qui satisfait à un scénario particulier. [35] Avantageusement, le système selon l'invention comprend en outre un troisième calculateur adapté pour définir l'hydrologie des ressources et la 20 courbe temporelle préétablie de demande en eau produite sur des éléments demandeurs, en fonction de données historiques analysées par des algorithmes statistiques. [36] Ce troisième calculateur a pour objectif de générer plusieurs scénarios caractérisant les états des ressources du système et/ou caractérisant l'état 25 de la demande en eau des éléments demandeurs. Ces scénarios peuvent être générés statistiquement et peuvent être fonction de données historiques brutes introduites par l'utilisateur. [37] Le troisième calculateur, dont l'usage est facultatif, comprend un ensemble d'algorithmes statistiques qui génèrent les données d'entrée 30 utilisées par les deux premiers calculateurs. Ces données d'entrées sont 3054705 9 relatives soit aux scénarios hydrologiques sur les ressources en eau, soit à la courbe temporelle préétablie de demande en eau produite. [38] Par scénario, on entend au sens de la présente demande, le déroulement préétabli et prévisionnel d'une situation en fonction de 5 l'hydrologie des ressources en eau et/ou de la demande en eau des éléments demandeurs simulées. Les scénarios servent de données d'entrée aux algorithmes d'optimisation mono ou multi objectifs du premier et du second calculateur. [39] Ce troisième calculateur permet donc, par exemple, d'importer 10 directement, par l'utilisateur, des données sur l'hydrologie des ressources et/ou une courbe temporelle préétablie de demande en eau, avec ou sans modélisation statistique. [40] Il est à noter que les deux premiers calculateurs ou les trois calculateurs précités peuvent ne faire l'objet que d'un seul ensemble 15 calculateur, qui serait adapté pour porter et exercer chacune des fonctionnalités des deux premiers ou des trois calculateurs précités. [41] Dans un mode de réalisation particulier de l'invention, au moins une ressource en eau est une nappe aquifère dont le volume d'eau prélevable est déterminé par le niveau piézométrique. 20 [42] L'invention a également pour objet un procédé de gestion quantitative d'eau d'un système comprenant les éléments suivants : - au moins deux entités de production d'eau interconnectées, directement ou indirectement, chaque entité de production étant définie par un débit maximal de production, 25 - au moins une ressource en eau reliée à l'une au moins des entités de production, chaque ressource en eau étant définie par un débit maximal de prélèvement, - au moins un élément demandeur d'eau produite définie par une courbe temporelle préétablie de demande en eau produite en fonction 30 du temps, 3054705 10 chaque liaison entre entités de production, ressources en eau et éléments demandeurs étant assurée par un ouvrage de transfert ayant un débit maximal de transfert prédéterminé, ces ouvrages de transfert pouvant être interconnectés, 5 chaque entité de production et chaque ressource en eau étant en outre associée à une fonction de pondération P. [43] Avec cet objectif en vue, le procédé selon l'invention, par ailleurs conforme au préambule cité ci-avant, est essentiellement caractérisé en ce qu'il comprend les étapes suivantes : 10 - détermination de l'ensemble des fonctions de pondération globales Pg, et - sélection de la fonction de pondération globale Pg la plus basse. [44] Il est à noter que toutes les fonctions de pondération globales Pg déterminées respectent les contraintes de débit des entités de production et 15 des ressources en eau. [45] Dans un mode de réalisation particulier de l'invention, le procédé selon l'invention comprend en outre une étape de détermination par le niveau piézométrique d'un volume disponible et/ou prélevable d'une ressource en eau du système. 20 [46] Cette étape de détermination par le niveau piézométrique d'un volume disponible et/ou prélevable d'une ressource en eau du système peut être appliquée à tout type de système indépendamment du système de l'invention. Elle peut être utilisée en dehors du cadre de l'invention afin d'améliorer les pratiques d'exploitation des captages d'eau souterraine, 25 grâce à la définition d'un niveau critique d'exploitation et à la détermination de volumes d'eau souterraine prélevables. [47] De préférence, l'étape de détermination par le niveau piézométrique d'un volume d'eau prélevable d'une ressource en eau souterraine comprend les étapes suivantes : 3054705 11 - détermination des chroniques temporelles historiques de la nappe aquifère au point de prélèvement en mesurant la cote piézométrique pseudo-statique de la ressource en eau sur une durée prédéfinie ; - détermination d'une relation entre la cote piézométrique pseudo- 5 statique et la transmissivité au point de prélèvement ; - détermination du niveau critique d'exploitation au point de prélèvement ; - sélection d'une cote piézométrique régionale de référence ; - détermination d'une relation empirique entre la cote piézométrique 10 régionale moyenne mensuelle de référence et la cote pseudo- statique moyenne mensuelle au point de prélèvement. [48] Avantageusement, l'étape de détermination par le niveau piézométrique d'un volume d'eau souterraine prélevable comprend en outre les étapes : 15 - transformation de la cote piézométrique régionale mensuelle moyenne de référence mesurée en cote pseudo-statique mensuelle moyenne au point de prélèvement en appliquant la relation empirique ; - calcul du rabattement maximum admissible à chaque pas de 20 temps par interpolation linéaire ; - calcul du volume maximum d'eau prélevable. Brève description de la figure [49] L'invention sera mieux comprise à la lecture de la description qui suit, 25 faite uniquement à titre d'exemple, et en référence à la figure en annexe qui représente un exemple de système selon l'invention. 30 3054705 12 Description détaillée Introduction [50] Les objectifs de l'invention sont, d'une part, de minimiser la fonction de 5 pondération globale Pg d'un système de gestion quantitative d'eau dans le but d'obtenir un système le plus autonome possible, utilisant la plus faible consommation d'énergie possible, tout en respectant ses contraintes. Particulièrement, la présente invention décrit notamment un système de gestion quantitative d'eau comprenant les éléments suivants : 10 a) au moins deux entités de production U d'eau interconnectées, chaque entité de production U étant définie par un débit maximal de production, b) au moins une ressource en eau S reliée à l'une au moins des entités de production U, chaque ressource en eau S étant définie 15 par un débit maximal de prélèvement, c) au moins un élément demandeur D d'eau produite défini par une courbe temporelle préétablie de demande en eau produite en fonction du temps. 20 Les entités de production d'eau interconnectées [51] Au sens de la présente invention, on entend par entité de production U d'eau, toute entité adaptée pour produire de l'eau, potable ou non, ainsi que toute entité adaptée pour transférer de l'eau à un élément quelconque auquel elle est reliée, sans que cette eau transférée ne soit nécessairement 25 traitée. [52] Ainsi, par eau produite, eau traitée, eau transférée, on entend au sens de la présente invention, tous les volumes d'eau en provenance de l'une des entités de production U du système selon l'invention. [53] Il est à noter que, dans l'invention, les entités de production U d'eau 30 ont pour fonction de produire de l'eau sans pour autant qu'il y ait de traitement et/ou de traiter un volume d'eau introduit dans l'entité de 3054705 13 production U, ou de transférer un volume d'eau vers une autre entité du système à laquelle elle est reliée. Les ressources en eau 5 [54] Les ressources en eau S du système peuvent être des ressources en eau brute souterraines (sources, nappes aquifères, karts... ) équipées de captages. Les ressources en eau S peuvent également être des eaux de surface équipées de captages (de prises) d'eau (lacs, cours d'eau...). [55] Selon une variante de réalisation de l'invention, au moins une 10 ressource en eau S du système à laquelle est reliée au moins une entité de production U est une nappe aquifère dont le volume d'eau prélevable est déterminé par le niveau piézométrique. [56] Il est à noter qu'une ressource en eau S peut être reliée à une ou plusieurs entités de production U du système selon l'invention. 15 Les éléments demandeurs d'eau produite [57] Les éléments demandeurs D d'eau produite peuvent être des collectivités, des industriels, des entités agricoles, ou tout autre élément demandeur D désirant être alimenté en eau traitée ou non. 20 [58] Un élément demandeur D est en général relié à une seule entité de production U. Toutefois, selon l'invention, un élément demandeur D peut être relié à différentes entités de production U. Dans cette configuration, selon l'invention, l'élément demandeur D sera alimenté par l'eau en provenance de l'entité de production U associée à la fonction de pondération 25 P la plus faible. Il est à noter qu'un élément de type château d'eau peut également se situer entre une ou plusieurs entités de production P et un ou plusieurs éléments demandeurs D. [59] II est à noter qu'un élément demandeur D peut être relié à une ou plusieurs entités de production U du système selon l'invention. 30 3054705 14 La courbe temporelle préétablie de demande en eau [60] La demande en eau par les éléments demandeurs D est instantanée et peut être suivie sur une courbe temporelle. Cependant, il est également possible de réaliser une courbe temporelle de demande en eau par 5 simulation afin d'anticiper la demande, ou d'observer des différences en fonction des demandes simulées afin de planifier sur une durée prédéterminée les volumes d'eau à prélever et à produire en fonction de la demande en eau. Cette courbe de demande en eau préétablie peut donc être obtenue par simulation, à l'aide du troisième calculateur (voir ci-après), 10 suite à l'utilisation de données antérieures liées au système et relatives à des situations passées, ou à l'aide de données permettant d'anticiper la demande future. Les liaisons entre chaque entité 15 [61] Pour que les différents éléments du système puissent échanger de l'eau entre eux, il est nécessaire que chaque liaison entre entités de production U, ressources en eau S et éléments demandeurs D soit assurée par un ouvrage de transfert C ayant un débit maximal de transfert prédéterminé limitant ainsi les débits d'échange entre les différents éléments 20 du système. [62] Il est à noter que selon l'invention, les ouvrages de transfert C du système reliant chacun des éléments peuvent comprendre des vannes V afin de réguler les échanges d'eau entre les différents éléments du système selon l'invention. 25 [63] Par exemple, selon l'invention, les ouvrages de transfert C pouvant être utilisés sont les aqueducs, les galeries, les canaux, les canalisations, les dispositifs aériens ou souterrains de transfert d'eau. Fonction de pondération 30 [64] Selon l'invention, chaque entité de production U et chaque ressource en eau S du système est associée à une fonction de pondération P. Cette 3054705 15 fonction de pondération P tient compte de nombreux critères techniques variables et constants et est définie comme étant une somme de valeurs variables et de valeurs constantes. [65] Par exemple, les valeurs constantes peuvent être des valeurs 5 arbitraires, définies par l'utilisateur, priorisant l'allocation d'eau dans les éléments du système. [66] Les valeurs variables, quant à elles, peuvent par exemple être fonction : - des priorités définies par l'utilisateur sur les éléments demandeurs 10 D qui régissent l'allocation de flux dans le système ; ainsi, les éléments demandeurs D à plus forte priorité seront desservis les premiers ; et peuvent être également fonction : - des ouvrages de transfert C du système, chaque ouvrage de 15 transfert C du système étant caractérisé par un coefficient. Plus ce coefficient est élevé, plus l'ouvrage de transfert C en question consomme de l'énergie pour réaliser un transfert d'eau. [67] Ces coefficients sont intégrés aux fonctions de pondération P de chaque ressource en eau S et de chaque entité de production U pour 20 satisfaire la courbe temporelle préétablie de demande en eau produite. De cette façon, par exemple, la fonction de pondération P d'une ressource en eau S du système reliée à une entité de production U du système par un ouvrage de transfert C consommant une importante quantité d'énergie est plus élevée que celle considérant un ouvrage de transfert C consommant 25 une faible quantité d'énergie. [68] Selon l'invention, plus la fonction de pondération P d'un élément du système est élevée, plus l'eau qu'elle délivre a demandé une consommation d'énergie élevée et/ou plus les contraintes hydrologiques / hydrauliques / réglementaires s'appliquant à l'élément du système en question sont 30 restrictives. Il est à noter que l'objectif est de délivrer l'eau la moins consommatrice d'énergie possible tout en garantissant au mieux 3054705 16 l'autonomie du système selon l'invention, et tout en respectant les contraintes en vigueur relatives aux ressources en eau S et aux entités de production U, à au moins un élément demandeur D, en fonction également de ses contraintes. Donc, de préférence, une entité de production U reliée à 5 plusieurs ressources en eau S prélève en priorité un volume d'eau de la ressource associée à la plus faible fonction de pondération P de façon à minimiser sa fonction de pondération. [69] Pour une stratégie particulière, c'est-à-dire pour un état caractérisant les ressources du système S et/ou caractérisant la demande en eau des 10 éléments demandeurs D, chaque entité de production U d'eau et chaque ressource en eau S du système est associée à plusieurs fonctions de pondération respectant l'allocation d'eau aux éléments demandeurs et les transferts d'eau satisfaisant à la demande en eau. Ainsi, une stratégie particulière est déterminée pour une période donnée selon l'hydrologie des 15 ressources en eau S (débits des cours d'eau, niveaux piézométriques des aquifères) et de la demande en eau des éléments demandeurs D. [70] Il est à noter que plusieurs stratégies peuvent être le résultat d'un scénario particulier. En effet, le système, pour un scénario particulier, est défini par un ensemble de fonction de pondération P attribué à chaque entité 20 de production U et ressource en eau S du système. Toutefois, pour un scénario donné, différentes stratégies sont possibles, donc différents ensembles de fonctions de pondération P sont possibles. [71] Une fonction de pondération globale Pg définit donc chaque stratégie et correspond à la somme de toutes les fonctions de pondération P de 25 chacun des éléments du système pour un scénario particulier. [72] Il est à noter que la meilleure stratégie pour un scénario particulier correspond à la fonction de pondération globale Pg la plus faible. Le premier calculateur 30 [73] Le système selon l'invention comprend un premier calculateur adapté, dans un premier temps, pour déterminer et calculer l'ensemble des fonctions 3054705 17 de pondération P de chacun des éléments du système pour une stratégie particulière qui satisfait à 100% de la demande en eau des éléments demandeurs D. Dans un second temps, ce premier calculateur est adapté pour calculer et déterminer la fonction de pondération globale Pg, somme 5 des fonctions de pondération de chacun des éléments du système, la plus faible qui permet de satisfaire à 100 % de la demande en eau par les éléments demandeurs D. [74] Il est à noter que selon un mode préférentiel de l'invention, les vannes V des différents ouvrages de transfert C sont commandées par ce premier 10 calculateur de sorte que le système délivre aux éléments demandeurs D, en temps réel, de l'eau produite selon la courbe temporelle préétablie en minimisant la fonction de pondération globale Pg. [75] Il est également à noter que, selon l'invention, plusieurs entités de production U sont reliées entre elles. Ainsi, selon un mode préférentiel de 15 l'invention, l'entité de production U associée à la fonction de pondération P la plus faible transmet de l'eau produite en priorité à l'entité de production d'eau U associée à la fonction de pondération P la plus élevée. De cette façon, si un élément demandeur D n'est relié qu'à une seule entité de production U, la fonction de pondération P de cette entité de production U 20 peut être ajustée suite à des échanges d'eau avec une ou plusieurs autres entités de production U du système auxquelles elle est reliée afin de permettre la transmission d'eau la moins consommatrice d'énergie possible à l'élément demandeur D. Ces échanges d'eau sont commandés par le premier calculateur qui actionne ou non les vannes V des ouvrages de 25 transfert C pour ainsi permettre des échanges d'eau entre les différentes entités de production U interconnectées du système. [76] Aussi, dans le cas où plusieurs entités de production U d'eau sont connectées entre-elles, il est possible selon l'invention que l'entité de production U définie par la plus grande fonction de pondération P transmette 30 de l'eau à celle définie par la plus faible fonction de pondération. Par succession, les entités de production U définies par des fonctions de 3054705 18 pondération élevées peuvent transmettre de l'eau à celles définies par des fonctions de pondération plus faibles. De cette façon, les fonctions de pondération des éléments du système tendent à s'harmoniser pour, autant que faire se peut transmettre l'eau la moins consommatrice d'énergie aux 5 éléments demandeurs D. [77] Dans cette configuration, les fonctions de pondération P et l'équilibre entre les ressources en eau S et la demande en eau sont pris en compte pour réaliser une optimisation mono-objectif : satisfaire 100% de la demande. 10 Second calculateur [78] De préférence, un second calculateur peut également être utilisé et adapté pour définir différentes stratégies de satisfaction de la demande. Ce second calculateur est adapté pour générer, en fonction des scénarios, des 15 stratégies de production répondant à différents pourcentages de satisfaction de la demande. Ces différents pourcentages de satisfaction peuvent être, par exemple, équivalents à 70%, 80%, 8,0,/o, 92% de la demande en eau. Pour chaque pourcentage de satisfaction, ce second calculateur, en collaboration avec le premier calculateur, est apte à lui associer une fonction 20 de pondération globale Pg selon une stratégie de production. Ce second calculateur est en outre adapté pour mettre en avant la fonction de pondération globale Pg la plus faible du système pour un pourcentage de satisfaction particulier de la demande. [79] Ainsi, il est possible de simuler, en fonction de différents paliers de 25 satisfaction de la demande, les échanges d'eau entre les différents éléments du système et également les volumes d'eau à prélever aux ressources en eau S. [80] II est à noter également que le second calculateur est en outre adapté pour déterminer les stratégies optimales qui forment le front de Pareto pour 30 une entité de production U en particulier. Ces stratégies optimales sont définies comme étant les stratégies fournissant un pourcentage maximal 3054705 19 pour un ratio minimisé défini comme la variation de la consommation d'énergie, dont dépend la fonction de pondération P de l'entité de production particulière, sur la variation de la demande. Ainsi, on peut réaliser, par simulation, une courbe de pourcentage de satisfaction de la demande (en 5 %) en fonction de la consommation d'énergie d'une entité de production U particulière du système selon l'invention. À l'aide de cette courbe, on peut observer que la fonction de pondération globale Pg de l'entité de production U augmente en fonction du taux de satisfaction de la demande, jusqu'à un taux de satisfaction, par exemple proche de 99%. On peut alors observer 10 une rupture de pente (la pente diminue) ce qui signifie, qu'après la rupture de pente, la consommation d'énergie va augmenter plus vite que la satisfaction de la demande. Une seconde diminution de pente peut par exemple être observée aux environs de 99,6% de satisfaction de la demande. 15 [81] En conséquence, le résultat le plus intéressant, c'est-à-dire pour lequel la consommation d'énergie est la plus faible, au vu de la courbe pour de telles demandes, correspond au premier changement de pente. En effet, lorsque le taux de satisfaction de la demande est de 99% pour une fonction de pondération X, cette fonction de pondération X augmente plus 20 rapidement au-dessus d'une satisfaction de la demande supérieure à 99%. Ainsi, à l'aide de cette courbe générée par le second calculateur, il est remarquable que pour satisfaire 100% de la demande, il serait plus judicieux d'importer de l'eau provenant d'une autre entité de production U extérieure au système. 25 [82] En d'autres termes, à partir d'un taux de satisfaction de la demande de 99%, la consommation d'énergie de l'entité de production U en question augmente. En conséquence, selon l'invention, le premier calculateur peut ouvrir une vanne V qui relie cette entité de production U à une autre pour que cette entité, lorsqu'elle atteint, dans notre exemple, 99 % de satisfaction, 30 puisse recevoir de l'eau d'une autre entité dans le but de ne pas augmenter 3054705 20 sa consommation d'énergie et d'ainsi ne pas augmenter la fonction de pondération globale Pg du système. [83] Ainsi, à l'aide de ce second calculateur, il est possible d'apprécier à quel moment il est judicieux, par exemple, d'échanger de l'eau entre 5 éléments du système. Dans cette configuration, les fonctions de pondération P et l'équilibre entre les ressources en eau S et les entités demandeurs en eau D sont pris en compte pour réaliser une optimisation multi-objectifs : satisfaire au mieux la demande pour une fonction de pondération globale Pg la plus faible possible, tout en minimisant la consommation d'énergie liée à 10 la production de cette eau. Il est ainsi possible de proposer des stratégies de production d'eau produite déterminées selon deux objectifs antagonistes, la satisfaction de la demande définie par la courbe temporelle préétablie et la consommation en énergie associée. 15 Troisième calculateur [84] Dans un mode de réalisation préférentiel de l'invention, un troisième calculateur peut également être utilisé pour définir les données d'entrée du premier calculateur et du deuxième calculateur : les scénarios de l'hydrologie des ressources S (débits des cours d'eau, niveaux 20 piézométriques des aquifères), ainsi que les scénarios de courbes temporelles de demande en eau produite sur les entités de demande en eau D. Ainsi, ce troisième calculateur simule les scénarios de l'hydrologie des ressources S et de demande en eau sur les éléments demandeurs D en se basant sur des données historiques introduites par un utilisateur. 25 [85] Toutefois, les données d'entrée ne sont pas nécessairement simulées par le troisième calculateur, et peuvent être créées par un utilisateur de à partir de données brutes. [86] Ce troisième calculateur est adapté, par exemple, pour générer des séries de données prospectives sur différents horizons de temps (1 mois, 3 30 mois, 6 mois), en se basant sur l'analyse statistique fréquentielle des séries de données temporelles historiques de chaque ressource S et élément 3054705 21 demandeur D, quantifiant le temps de retour probabiliste des événements observés. [87] Concernant l'état de la ressource en eau S (débits, niveaux piézométriques), ce troisième calculateur peut par exemple proposer trois 5 scénarios types : année humide (périodes de retour de 2,5 ans ; 5 ans ; 10 ans), année médiane, année sèche (période de retour de 2,5 ans, 5 ans, 10 ans). [88] Et concernant la demande en eau, ce troisième calculateur peut proposer les scénarios types suivants : demande minimale, demande 10 médiane, demande maximale, quantile fixé par l'utilisateur. Liens entre les trois calculateurs [89] Il est à noter que les trois calculateurs précédemment mentionnés peuvent ne faire l'objet que d'un seul calculateur (ou deux) adapté pour 15 répondre aux différentes attentes, et apporter les mêmes fonctionnalités de ces trois calculateurs. [90] L'utilisation des trois calculateurs permet, par exemple, de simuler l'hydrologie d'une ou plusieurs ressources en eau S, sur une durée donnée et prédéterminée, et de simuler des valeurs de volume d'eau souhaité par 20 un ou plusieurs éléments demandeurs D afin d'anticiper au mieux les quantités d'eau qu'il serait nécessaire de prélever aux ressources en eau S et les volumes d'eau à produire par les entités de production U du système pour satisfaire les demandes en eau. Cette anticipation peut permettre de minimiser la fonction de pondération globale Pg du système dans la mesure 25 où si la demande en eau est trop élevée au vu du volume d'eau qui peut être produit et/ou prélevé, il est possible d'anticiper les échanges d'eau entre les différentes entités de production U interconnectées du système. [91] Ou encore, les trois calculateurs peuvent être utilisés pour planifier sur une durée prédéterminée les volumes d'eau à prélever dans au moins une 30 ressource en eau S ainsi que les volumes d'eau à produire en fonction de la 3054705 22 courbe temporelle de demande préétablie définie sur au moins un élément demandeur D. [92] Cet ensemble de trois calculateurs est en outre adapté pour rejouer des situations passées afin d'améliorer et/ou perfectionner le système en 5 minimisant les fonctions de pondération P des éléments du système, pour mieux appréhender des situations futures. [93] Cet ensemble de trois calculateurs permet également de simuler les effets de modifications du système selon l'invention : suppression ou ajout d'une entité de production U, suppression ou ajout d'une ressource en eau 10 S, suppression ou ajout d'un ouvrage de transfert C reliant une entité de production U à un quelconque élément du système, ajout, modification, ou suppression d'une contrainte. Ainsi, on peut simuler des aménagements du système et contrôler les variations des fonctions de pondération P des éléments du système, et notamment celles des fonctions de pondérations 15 globales Pg en fonction des scénarios appliqués. [94] Cet ensemble de trois calculateurs peut être utilisé en mode « Gestion opérationnelle » pour une planification hebdomadaire et mensuelle des prélèvements et de la production, ou en mode « Étude et perfectionnement » pour rejouer des situations passées afin de mieux gérer des situations 20 futures, ou pour évaluer l'effet de nouveaux aménagements. Module Volume Prélevable [95] Selon un mode préférentiel de l'invention, au moins une ressource en eau S du système à laquelle est reliée au moins une entité de production U 25 est une nappe aquifère dont le volume d'eau prélevable est déterminé par le niveau piézométrique. Dans cette configuration, selon l'invention, la méthode 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 sur l'état quantitatif 30 de l'aquifère. Cette méthode est simple, robuste et applicable à des aquifères continus (non karstiques, non fracturés), pour de multiples 3054705 23 configurations de captages (forages, puits à drains, champs captants) à l'exception des captages de sources. [96] Il est à noter que l'on peut se trouver dans la configuration où les volumes d'eau prélevables sont déterminés par le débit de la ressource en 5 eau S et non par la valeur de niveau piézométrique. [97] Il est à noter que cette méthode peut être appliquée à tout type de système indépendamment du système décrit dans la présente invention. [98] Cette méthode permet d'évaluer les volumes d'eau prélevables d'une nappe aquifère, et ainsi de déterminer des données essentielles à 10 l'évaluation de la fonction de pondération P d'une telle ressource en eau S. Ces données sont alors transmises au premier calculateur pour qu'il puisse calculer les valeurs variables qui définissent une telle ressource en eau S, puis la fonction de pondération P de cette source S pour calculer la fonction de pondération globale Pg du système en temps réel, et la minimiser. 15 [99] Il est à noter que la gestion quantitative des eaux souterraines et la définition de volumes prélevables au sein d'un aquifère ou d'un système aquifère nécessite de définir au préalable les limites spatiales de l'ensemble considéré. Ces limites sont définies en fonction des caractéristiques géologiques et hydrogéologiques locales. 20 [100] La définition de volumes prélevables d'eau souterraine au sein d'un aquifère ou d'un système aquifère nécessite également de disposer d'une bonne connaissance de l'état des ressources en eau S, 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 25 intrinsèque. Ces limites sont définies en fonction des caractéristiques générales du captage d'exploitation. [101]La gestion quantitative des eaux souterraines nécessite donc une connaissance de l'état des ressources en eau S, des prélèvements et des besoins. 30 [102]Ainsi, selon l'invention, les données d'exploitation transmises au premier calculateur comprennent donc des données relatives : 3054705 24 - aux caractéristiques générales du captage d'exploitation : par exemple, la coupe géologique qui est une représentation schématique de la succession des formations géologiques sur un profil vertical, avec indication des cotes associées, ou encore la 5 coupe technique qui est une représentation schématique des caractéristiques et des cotes des équipements d'un ouvrage souterrain sur un profil vertical [exemples : trou nu, caractéristiques du pré-tubage et du tubage, massif filtrant, cimentation, protection en tête]), 10 - aux données géologiques et hydrogéologiques : par exemple la piézométrie régionale, la transmissivité ou encore le coefficient d'emmagasinement, et - aux données d'exploitation : par exemple les mesures de niveau et cote du repère, les volumes prélevés, les débits nominaux d'une 15 ou plusieurs des pompes. [103] D'autres données d'exploitation comme le nombre d'heures journalières de pompage et le débit instantané moyen journalier sont facultatives. [104] Par ailleurs, cette méthodologie s'appuie sur l'utilisation de la 20 solution analytique de Cooper-Jacob (Cooper, H.H. and C.E. Jacob, 1946. A generalized graphical method for evaluating formation constants and summarizing well field history, Am. Geophys. Union Trans., vol. 27, pp. 526534) couramment utilisée en hydrogéologie quantitative. Elle fournit une estimation du volume maximum d'eau souterraine prélevable en fonction du 25 rabattement, lui-même dépendant du niveau piézométrique statique de la nappe. [105]Cette méthodologie comprend, d'une part, les étapes suivantes de calibration : - détermination des chroniques temporelles historique de la nappe 30 aquifère au point de prélèvement en mesurant la cote 3054705 25 piézométrique pseudo-statique de la ressource en eau S sur une durée prédéfinie ; - détermination d'une relation entre la cote piézométrique pseudostatique et la transmissivité au point de prélèvement ; 5 - détermination du niveau critique d'exploitation au point de prélèvement ; - sélection d'une cote piézométrique régionale de référence ; - détermination d'une relation empirique entre la cote piézométrique mensuelle moyenne de référence et la cote pseudo-statique 10 mensuelle moyenne au point de prélèvement. [106] Par exemple, les chroniques temporelles peuvent être analysées sur un minimum de deux ans. [107] La relation entre la cote piézométrique pseudo-statique et la transmissivité au point de prélèvement peut être déterminée, à l'aide des 15 données d'exploitation de la nappe aquifère, par la relation analytique de Cooper-Jacob. [108]11 est à noter que le niveau critique d'exploitation d'un puits ou d'un forage peut être conditionné par : - la cote altimétrique supérieure de la partie crépinée du tubage ; 20 - la cote d'exploitation locale limite de la nappe, qui peut être : - une cote de maintien de captivité de la nappe, - une cote de non-invasion de biseau salé, - une cote de non-dénoyage d'une zone productive, - une cote réglementaire (piézométrie objectif, seuil d'alerte... ) 25 - la cote de la crépine d'aspiration de la pompe ou la cote de sécurité déclenchant l'arrêt de la pompe. [109] La cote la plus élevée, donc la plus défavorable, sera retenue comme niveau critique d'exploitation de l'ouvrage, notée zNc. La détermination du niveau critique nécessite donc de disposer de la coupe technique et de la 30 coupe géologique de l'ouvrage, et de l'existence éventuelle de cotes réglementaires de gestion de la nappe. 3054 705 26 [110] Dans le cas d'un champ captant, le niveau critique le plus défavorable sera retenu et appliqué au captage conceptuel. Il est à noter que la sélection d'un piézomètre régional de référence nécessite de lister l'ensemble des piézomètres captant la masse d'eau 5 souterraine étudiée. Ces piézomètres peuvent être identifiés facilement en consultant la base de données ADES, (www.ades.eaufrance.fr), pour la France par exemple, 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 (idéalement supérieur à 10 10 ans) seront considérés. [111] Le piézomètre régional de référence sera 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 mensuelle de la cote piézométrique statique mesurée au piézomètre sur la 15 durée d'observation commune. [112] 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 moyenne des cotes mensuelles statiques au piézomètre 20 hspt(t) = hsp(t) + (h' - lisp) [113] Dans le cas de plusieurs piézomètres avec une RMSE proche, celui qui dispose de l'historique le plus long sera retenu. [114]Une relation empirique, le plus souvent linéaire, ou constituée de 25 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 3054705 27 plusieurs régressions linéaires sur un graphique de type scatter-plot, ou d'autres fonctions de corrélation. [115] La méthodologie peut comprendre, d'autre part, les étapes suivantes de simulation : 5 - transformation de la cote piézométriques régionale mensuelle moyenne de référence mesurée en un temps donné en cote pseudo-statiques mensuelle moyenne au point de prélèvement en appliquant la relation empirique ; - calcul du rabattement maximum admissible à chaque pas de 10 temps par interpolation linéaire ; - calcul du volume maximum d'eau prélevable, à chaque pas de temps, en fonction du rabattement maximum admissible et de la relation entre la cote piézométrique pseudo-statique et la transmissivité au captage, à l'aide de la relation de Cooper-Jacob. 15 [116]11 est à noter que cette méthode permet de calculer un volume maximum prélevable théorique sur la base d'une cote piézométrique pseudo-statique 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. 2 0 [117] La sous-étape de transformation de la cote statique au piézomètre de référence en cote piézométrique pseudo-statique à la ressource en eau S sera réalisée uniquement dans le cas d'une simulation prédictive. Une chronique piézométrique mensuelle prédictive de niveaux statiques sera générée sur le piézomètre de référence à l'aide du troisième calculateur. 25 [118] La relation ou l'ensemble de relations linéaires empiriques, obtenues lors de la précédente étape, permettant de calculer 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, seront utilisées. 30 [119] La longueur de l'historique disponible sur les piézomètres régionaux de suivi, et la fiabilité des données mesurées, rend plus fiable l'utilisation du 3054705 28 troisième calculateur sur ce point plutôt que directement le captage, d'où l'utilisation d'une méthode de transposition. [120] 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- 5 statique hp, et le niveau critique d'exploitation, noté zNc. Smax(t) = hp,(t) - zii, [121] 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 10 la relation T = f(hps) obtenue lors de la précédente étape. Il dépend du rabattement maximum admissible Smax, - Vmax Qmax- texp 4TrT(hps)Smax texp 2,303 log (2.25T(hps)texp) r2S texp étant la durée maximum d'exploitation. [122] Dans le calcul du volume maximum prélevable, la valeur de texp est 15 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 qu'il contient. [123]Ainsi, cette méthodologie prend en compte indirectement la situation hydrologique, les effets de recharge/décharge de la nappe et les 20 phénomènes naturels, qui influencent les cotes pseudo-statiques mesurées au captage. Elle prend également en compte indirectement les effets des prélèvements globaux sur la nappe et les activités anthropiques, qui influencent les cotes pseudo-statiques mesurées au captage. [124]Ainsi, selon cette méthode préférentielle selon l'invention, il est 25 possible de réaliser des prévisions de disponibilité future de la ressource en eau S basées sur : - les tendances hydriques connues au jour le jour, 3054705 29 - l'intégration comme variable du changement climatique dans les scénarios à moyen et long terme, [125] Le suivi quantitatif des eaux souterraines est surveillé, selon le mode préférentiel de l'invention, en surveillant les variations de niveau 5 piézométrique des nappes. Toutefois ce suivi quantitatif peut également être surveillé en mesurant le débit des émergences (sources) selon le type d'aquifère. La gestion quantitative des eaux souterraines revient dans la plupart du temps à fixer des niveaux piézométriques objectifs, ce qui revient à fixer des volumes prélevables. 10 [126] L'invention a été illustrée et décrite en détail dans les dessins et la description précédente. Celle-ci doit être considérée comme illustrative et donnée à titre d'exemple et non comme limitant l'invention à cette seule description. De nombreuses variantes de réalisation sont possibles. 15 Exemple de système selon l'invention [127]Selon l'invention, on peut avoir un système de gestion quantitative d'eau comprenant les éléments suivants (voir figure) : - deux entités de production U1 et U2 : une première entité de production U1 et une seconde entité de production U2 d'eau sont 20 interconnectées par deux ouvrages de transfert Cl et C2 comprenant chacun une vanne V1 et V2 ; - la première entité de production U1 est reliée à deux ressources en eau S1 et S1' par deux ouvrages de transfert C3 e C4 et la seconde entité de production U2 est reliée à une seule ressource en eau 25 S2 par un ouvrage de transfert C5 ; - trois éléments demandeurs D1, D1' et Dl" sont reliés à la première entité de production U1 par des ouvrages de transfert C6, C7 et C8, et un seul élément demandeur D2 est relié à la seconde entité de production U2 par un ouvrage de transfert C9, chaque élément 30 demandeur étant caractérisé par une courbe temporelle préétablie de demande en eau produite. 3054705 30 [128] Les éléments demandeurs sont caractérisés par une courbe temporelle préétablie de demande en eau produite à laquelle il faut satisfaire à 100%. Toutefois, il se peut qu'il ne soit pas possible de satisfaire à 100% la demande. Cependant, à l'aide des calculateurs, il est possible de 5 maximiser la satisfaction de la demande tout en minimisant la fonction de pondération globale Pg' du système. [129] Dans cette configuration, chaque entité du système est associée à une fonction de pondération P. En d'autres termes, la première entité de production U1 est associée à la fonction de pondération PU1, la seconde 10 entité de production U2 est associée à la fonction de pondération PU2, les ressources en eau S1, Si' et S2 sont respectivement associées aux fonctions de pondérations PS1, PS1' et PS2. Le but étant de minimiser la fonction de pondération globale Pg', somme de toutes les fonctions de pondération des éléments du système. 15 [130] Pour ce faire, un algorithme d'optimisation détermine l'ensemble des fonctions de pondération globales Pg', permettant de satisfaire à la courbe temporelle préétablie demande en eau, en fonction des différentes fonctions de pondération P des éléments du système. Et finalement, l'algorithme sélectionne la fonction de pondération globale Pg' la plus basse. 20 [131] Selon l'invention, le premier calculateur, qui peut récolter les données de l'algorithme, commande les vannes V1 et V2 pour que les entités de production U1 et U2 s'échangent de l'eau entres elles afin de réduire au maximum la fonction de pondération globale Pg' du système dans le but que cette fonction de pondération globale Pg' corresponde à celle sélectionnée 25 par l'algorithme. [132] Par exemple, on peut supposer que la fonction de pondération PU2 est plus élevée que celle de PU1. Ceci peut être dû, par exemple, au fait que l'exploitation de la ressource S2 soit plus consommatrice d'énergie que celle de S1 et/ou celle de Si', en d'autres termes que la fonction de pondération 30 PS2 soit plus élevée que PS1 ou PS1'. Dans cette situation, il est préférable que la seconde entité de production U2 transmette de l'eau en provenance 3054705 31 de la première entité de production U1 à l'élément demandeur D2. En effet, la première entité de production U1 produit de l'eau à moindre pondération par rapport à l'eau produite par la seconde entité U2. Ainsi, il est judicieux de stopper la production d'eau de la seconde entité de production U2, ou 5 tout du moins de diminuer le prélèvement d'eau à la ressource S2. [133]À l'aide de cet échange entre les entités de production U1 et U2, la fonction de pondération globale Pg' du système peut être minimisée dans la mesure où la fonction de pondération de la ressource en eau S2 n'est plus totalement à considérer. 10 [134]Ainsi, la seconde entité de production U2 transmet de l'eau à moindre pondération à l'élément demandeur D2 et la première entité de production U1 transmet quant à elle de l'eau aux éléments demandeurs D1, D1' et Dl" qui proviennent des ressources S1 et S1'. Il est à noter que si la fonction de pondération PS1 est plus élevée que celle de PS1', l'entité de production U1 15 prélèvera de l'eau en priorité de la ressource S1'. [135] Les indications ci-dessus prévalent lorsque l'on se trouve dans une configuration instantanée, avec commande de vannes, mais aussi en configuration de simulations. [136] Dans les revendications, le mot « comprenant » n'exclut pas d'autres 20 éléments et l'article indéfini « un/une » n'exclut pas une pluralité. REVENDICATIONS1. Système de gestion quantitative d'eau comprenant les éléments suivants : - au moins deux entités de production (U) d'eau interconnectées, directement ou indirectement, chaque entité de production (U) étant définie par un débit maximal de production, - au moins une ressource en eau (S) alimentant au moins une entité de production (U), chaque ressource en eau (S) étant définie par un débit maximal de prélèvement, - au moins un élément demandeur (D) d'eau produite définie par une courbe temporelle préétablie de demande en eau produite en fonction du temps, chaque liaison entre entités de production (U), ressources en eau (S) et éléments demandeurs (D) étant assurée par un ouvrage de transfert (C) ayant un débit maximal prédéterminé, ces ouvrages de transfert (C) pouvant être interconnectés, chaque entité de production (U) et chaque ressource en eau (S) étant en outre associée à une fonction de pondération P, ledit système étant caractérisé en ce qu'il comprend un premier calculateur adapté pour minimiser la fonction de pondération globale Pg du système, somme de toutes les fonctions de pondération P des différents éléments du système, tout en garantissant le respect de la courbe temporelle préétablie de demande en eau produite de chaque élément demandeur (D) et sous contrainte de respect des débits maximaux et/ou des débits minimaux des différents éléments du système. 2. Système selon la revendication 1, caractérisé en ce que la satisfaction de la courbe temporelle préétablie de demande en eau produite de chaque élément demandeur (D) d'eau produite est en outre sous contrainte de respect des volumes d'eau prélevables dans les ressources en eau (S), des 3054705 33 volumes d'eau produite par les entités de production (U) et des contraintes du système. 3. Système selon les revendications 1 ou 2, caractérisé en ce que les 5 différents ouvrages de transfert (C) comprennent des vannes (V) commandées par le premier calculateur de sorte à ce que le système délivre aux éléments demandeurs (D), en temps réel, de l'eau produite selon la courbe temporelle préétablie en minimisant la fonction de pondération globale Pg, tout en respectant les contraintes du système. 10 4. Système selon l'une quelconque des revendications 1 à 3, caractérisé en ce que l'entité de production (U) associée à la fonction de pondération P la plus faible transmet de l'eau produite en priorité à l'entité de production (U) d'eau associée à la fonction de pondération P la plus élevée et/ou 15 transmet de l'eau produite en priorité à l'élément demandeur (D) associé à la fonction de pondération P la plus élevée. 5. Système selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'une entité de production (U) reliée à plusieurs ressources en eau 20 (S) prélève en priorité un volume d'eau de la source associée à la plus faible fonction de pondération P. 6. Système selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il comprend en outre un second calculateur adapté pour définir des 25 stratégies de production pour satisfaire la demande en minimisant la consommation énergétique, chaque stratégie correspondant à un pourcentage de la courbe temporelle préétablie de l'élément demandeur (D), et ledit second calculateur étant connecté au premier calculateur de sorte que ledit premier calculateur génère une valeur minimale de la fonction de 30 pondération globale Pg pour chaque stratégie de production, ledit second calculateur étant en outre adapté pour déterminer la stratégie optimale 3054705 34 définie comme étant la stratégie fournissant le pourcentage maximal de satisfaction de la demande pour un ratio variation de la fonction de pondération globale Pg sur variation de la demande minimisé. 5 7. Système selon l'une quelconque des revendications 1 à 6, caractérisé en ce qu'il comprend en outre un troisième calculateur adapté pour définir l'hydrologie des ressources (S) et la courbe temporelle préétablie de demande en eau produite sur les éléments demandeurs (D), en fonction de données historiques analysées par des algorithmes statistiques. 10 8. Système selon l'une quelconque des revendications 1 à 7, caractérisé en ce qu'au moins une ressource en eau (S) est une nappe aquifère dont le volume d'eau prélevable est déterminé par le niveau piézométrique. 15 9. Procédé de gestion quantitative d'eau d'un système comprenant les éléments suivants : - au moins deux entités de production (U) d'eau interconnectées, directement ou indirectement, chaque entité de production (U) étant définie par un débit maximal de production, 20 - au moins une ressource en eau (S) reliée à l'une au moins des entités de production (U), chaque ressource en eau (S) étant définie par un débit maximal de prélèvement, - au moins un élément demandeur (D) d'eau produite définie par une courbe temporelle préétablie de demande en eau produite en fonction 25 du temps, chaque liaison entre entités de production (U), ressources en eau (S) et éléments demandeurs (D) étant assurée par un ouvrage de transfert (C) ayant un débit maximal de transfert prédéterminé, ces ouvrages de transfert (C) pouvant être interconnectés, 30 chaque entité de production (U) et chaque ressource en eau (S) étant en outre associée à une fonction de pondération P, 3054705 ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes : - détermination de l'ensemble des fonctions de pondération globales Pg, et - sélection de la fonction de pondération globale Pg la plus basse. 5 10. Procédé selon la revendication 9, caractérisé en ce qu'il comprend en outre une étape de détermination par le niveau piézométrique d'un volume disponible et/ou prélevable d'une ressource en eau (S) du système. 10 11. Procédé selon la revendication 10, caractérisé en ce que l'étape de détermination par le niveau piézométrique d'un volume d'eau prélevable d'une ressource en eau (S) comprend les étapes suivantes : - détermination des chroniques temporelles historiques de la nappe aquifère au point de prélèvement en mesurant la cote piézométrique 15 pseudo-statique de la ressource en eau (S) sur une durée prédéfinie ; - détermination d'une relation entre la cote piézométrique pseudostatique et la transmissivité au point de prélèvement ; - détermination du niveau critique d'exploitation au point de prélèvement ; 20 - sélection d'une cote piézométrique régionale de référence ; - détermination d'une relation empirique entre la cote piézométrique mensuelle moyenne de référence et la cote pseudo-statique mensuelle moyenne au point de prélèvement. 25 12. Procédé selon la revendication 11, caractérisé en ce que l'étape de détermination par le niveau piézométrique d'un volume d'eau prélevable comprend en outre les étapes : - transformation de la cote piézométrique régionale mensuelle moyenne de référence mesurée en un temps donné en cote pseudo- 30 statique mensuelle moyenne au point de prélèvement en appliquant la relation empirique ; 3054705 36 - calcul du rabattement maximum admissible à chaque pas de temps par interpolation linéaire ; - calcul du volume maximum d'eau prélevable.