HyD2002 - Les objets du modèle

There is no translation available.

HyD2002 - Objectifs - Interface graphique - Les objets du modèle - Fichiers d'entrée-sorties - Tutoriel

hydcomp.gif

 


 

Description

Les systèmes considérés sont constitués d'éléments de stockage tant superficiels que souterrains, d'éléments de captage, de transport tant naturels qu'artificiels, d'utilisation et de consommation d'eau.

Les modèles proposés permettent la simulation continue du fonctionnement des systèmes basé sur des règles de gestion très souples. Ces consignes font référence à des priorités de distribution de demandes et une hiérarchisation de la mobilisation des ressources. La mobilisation s'effectue en respectant les courbes de remplissage objectif des réservoirs.

L'utilisation du modèle consiste à modifier, au fur et à mesure des simulations, les consignes de gestion pour améliorer la satisfaction des objectifs.

Formulation

Nous reprenons, ici et au niveau de chacun des composants, la formulation mathématique définie dans le programme WARSYP [WARSYP 1998]. Cette formulation est axée vers l'optimisation de la gestion stratégique, mais elle permet de rendre explicite les caractéristiques des objets et les contraintes de gestion tactique.

Le problème consiste à minimiser une fonction objectif composée de termes pénalisation et de couts sur un horizon décisionnel, sachant les conditions formulées dans les équations (1)-(30) relatives aux divers composants ; les termes de pénalisation sont relatifs aux déficits ou aux déviations par rapport à des niveaux objectifs ; et les termes de couts sont relatifs aux pompages de surface ou dans les aquifères. La fonction objectif peut être exprimée comme suit :

{tex}\small \sum _{t \in T}\left[\sum _{i \in NE}\sum _{m \in M+it}\lambda^+_{imt.} V^+_{imt}+\sum _{i \in NE}\sum _{m \in M-it}\lambda^-_{imt.} V^-_{imt} + \\ \sum _{(e,f) \in NAN} \alpha_{eft.}Q^-_{eft}+ \sum _{n \in ND \cup NC}\sum _{f \in Fn} \delta_{nft.}D^-__{nft.}+ \sum _{n \in NS Ynt.}S^-_{nt} \\ +\sum _{k \in NR} \sigma_{kt}.R^-_{kt}+ \sum _{k \in NP} \sum _{n \in Nk} \mu_ {knt.} Y_{knt}+\sum _{(e,f) \in NA} \phi_{eft.}Q_{eft}+\ \\ \sum _{i \in NC} \omega_i (\sum _{f \in Pi}Q_{ift}- (d_{it}- \sum _{l \in Fi}D^-_{ilt})) \right]{/tex}

Le stockage dans le réservoir et les aquifères et la limitation des déficits cumulés des demandes nécessitent un couplage multipériodes, qui conduit à l'adoption d'un modèle complexe de résolution d'équations linéaires.

Définition

Propriétés structurales

# Systeme en-tête avec le nom du type d'objet
nom label de l'objet
description texte de complément d'information


Propriétés de gestion

simulation index de simulation (1..99). Les fichiers de simulation stockent les propriétés de gestion et les résultats relatifs à un même système d'eau.
mode processus de simulation - 1 à 3 simulation de la gestion tactique en mois, décades, jours respectivement -
periodes année début ; mois début (1..12) ; nombre de mois de l'horizon décisionnel pour l'optimisation, 1 sinon ; nombre de mois de simulation
ecoulement répertoire des données d'apports - unité : hm3
lame ecoulee répertoire des données d'apports exprimées en lames - unité : m
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Systeme
nom Amenagements_A2
description amenagements "1995" - Bras de Sable - AEP 50%
simulation 2
mode 2.decade
periodes 1962 1 1 348
ecoulement C:\ActiveHyd\guadelou\flows\hm3\
lame ecoulee C:\ActiveHyd\guadelou\flows\m\
sorties global an mois decade jour

hydram.gif

 


Figure B2.1. - Système d'eau - graphique associé

 


Ensembles principaux

{tex}\small NT{/tex} pas de temps successifs de l'horizon décisionnel
{tex}\small NE{/tex} ensemble des réservoirs, soit des noeuds de surface avec une capacité de stockage
{tex}\small ND{/tex} ensemble des noeuds demandes
{tex}\small NC{/tex} ensemble des noeuds demandes non consommatrices. {tex}\small NC = NCF \cup NCE{/tex}, où
{tex}\small NCF{/tex} ensemble des noeuds demandes non consommatrices sans capacité de stockage associé
{tex}\small NCE{/tex} ensemble des noeuds demandes non consommatrices associés à des réservoirs, {tex}\small NCE \subset NE{/tex}
{tex}\small NS{/tex} ensemble des noeuds demandes qui requièrent un volume de stockage réservé dans des réservoirs amont, {tex}\small NS \subseteq ND\cup NC{/tex}
{tex}\small NJ{/tex} ensemble des noeuds jonctions
{tex}\small NU{/tex} ensemble des noeuds pompages
{tex}\small NN{/tex} ensemble des noeuds de surface, {tex}\small NN=NE\cup ND\cup NCF\cup NJ\cup NU{/tex}
{tex}\small NK{/tex} ensemble des aquifères, noeuds souterrains avec une capacité de stockage
{tex}\small NA{/tex} ensemble des arcs de surface
{tex}\small NAN{/tex} ensemble des arcs naturels, {tex}\small NAN\subset NA {/tex}
{tex}\small NAB{/tex} ensemble des arcs artificiels, {tex}\small NAB\subset NA {/tex}
{tex}\small \Gamma {/tex} ensemble des pas de temps successifs utilisés pour controler le déficit cumulé de demande en eau sur ces périodes pour chaque demande et chaque période

Ensembles secondaires

{tex}\small M_i {/tex} ensemble des zones correspondant au découpage du réservoir i pour chaque période
{tex}\small F_n {/tex} ensemble des niveaux correspondant au découpage de la demande en eau du noeud n, {tex}\small n\in ND \cup NC {/tex}
{tex}\small P_e \subset NN {/tex} ensemble des noeuds directement à l'aval du noeud e, {tex}\small e \in NN-ND {/tex}. Si {tex}\small f \in P_e {/tex}, {tex}\small (e,f) \in NA{/tex}
{tex}\small P_e^{\prime} \subset NN {/tex} ensemble des noeuds à l'aval du noeud e, {tex}\small e \in NN-ND{/tex}
{tex}\small Q_f \subset NN-ND {/tex} ensemble des noeuds directement à l'amont du noeud f, {tex}\small f \in NN{/tex}. Si {tex}\small Q_f {/tex}, {tex}\small (e,f) \in NA{/tex}
{tex}\small Q_f^{\prime}\subset NN-ND{/tex} ensemble des noeuds à l'amont du noeud f, {tex}\small f \in NN{/tex}
{tex}\small NR{/tex} ensemble des aquifères avec recharge artificielle
{tex}\small N^k \subseteq NJ{/tex} ensemble des noeuds jonctions qui permettent la recharge artificielle de l'aquifère k, {tex}\small k \in NR{/tex}
{tex}\small N^k \subseteq Np{/tex} ensemble des aquifères avec pompage des demandes en eau
{tex}\small N_k \subseteq ND {/tex} ensemble des noeuds demandes avec pompage possible dans l'aquifère k, {tex}\small k \in NP{/tex}
{tex}\small J_k \subset NE \cup NJ{/tex} ensemble des réservoirs ou des noeuds jonctions avec échange avec l'aquifère k, {tex}\small k \in NK{/tex}
{tex}\small k(i){/tex} aquifère en relation avec le réservoir ou la jonction i, {tex}\small i \in NE \cup NJ{/tex}
{tex}\small j(n){/tex} noeud jonction vers lequel l'eau du noeud demande n est partiellement restituée au système, {tex}\small n \in ND{/tex}

Les paramètres du modèle

Paramètres d'état

{tex}\small b_{it}{/tex} apport extérieur au réservoir ou au noeud jonction i au pas de temps t, {tex}\small i \in NE \cup NJ{/tex}, {tex}\small t \in NT{/tex}. Si le noeud i est en tête du bassin, c'est à dire s'il n'a pas de noeud en amont, bit représente l'apport naturel. Mais s'il y a des noeuds plus en amont, bit représente l'écoulement net, correspondant à l'apport naturel du bassin intermédiaire à l'aval des noeuds amont.
{tex}\small r_{kt}{/tex} recharge extérieure de l'aquifère k pour le pas de temps t, {tex}\small k \in NK{/tex},
{tex}\small t \in NT{/tex}
{tex}\small V_{i0}{/tex} valeur du stockage du réservoir i au début de l'horizon décisionnel, {tex}\small i \in NE{/tex}
{tex}\small Z_{k0}{/tex} valeur du stockage de l'aquifère k au début de l'horizon décisionnel, {tex}\small k \in NK{/tex}
{tex}\small V^{m+}it V^{m-}it{/tex}
capacités de stockage maximale et minimale du réservoir i au pas de temps t, {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}
{tex}\small V^{m+}imt {/tex} capacité maximale de stockage de la zone m du réservoir i au pas de temps t, {tex}\small m \in M_i{/tex}, {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}tel que {tex}\small V^{m+}it = V^{m+}iMt{/tex}
{tex}\small Z^{m+}k Z^{m-}k{/tex} capacité de stockage maximale et minimale de l'aquifère k, {tex}\small k \in NK{/tex}
{tex}\small f^eit{/tex} coefficient de perte par évaporation par rapport au volume stocké dans le réservoir i pour le pas de temps t, iÎNE, tÎNT
{tex}\small f^1it{/tex} coefficient de perte par infiltration profonde par rapport au volume stocké dans le réservoir i pour le pas de temps t, {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}
{tex}\small f^1jt{/tex} coefficient de perte en ligne non récupérable par rapport au volume transité dans la jonction j pour le pas de temps t, {tex}\small j \in NJ{/tex}, {tex}\small t \in NT{/tex}
{tex}\small f^iit{/tex} coefficient de perte par infiltration dans l'aquifère k(i) par rapport au volume stocké dans le réservoir i pour le pas de temps t, {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}
{tex}\small f^ijt{/tex} coefficient de perte par infiltration dans l'aquifère k(j) par rapport au volume transité dans la jonction j pour le pas de temps t, {tex}\small j \in NJ{/tex}, {tex}\small t \in NT{/tex}
{tex}\small f^fkit{/tex} coefficient de réalimentation du réservoir ou du noeud jonction j par rapport au volume stocké dans l'aquifère k pour le pas de temps t, {tex}\small k \in NK{/tex}, {tex}\small i \in J_k{/tex}, {tex}\small t \in NT{/tex}
{tex}\small q^{m+}eft{/tex} limite supérieure de transit dans l'arc de surface (e, f), pour {tex}\small (e,f) \in NA-NAB{/tex}
{tex}\small \rho _n{/tex}
coefficient de restitution de la demande en eau du noeud demande n vers la jonction {tex}\small j(n){/tex}, {tex}\small n \in ND{/tex}
{tex}\small \xi^\tau n{/tex} fraction de demande en eau cumulée sur t pas de temps, dont le déficit est très fortement pénalisé, {tex}\small n \in ND \cup NC{/tex}, {tex} \small \tau \in \Gamma{/tex}
{tex}\small P^{m+}ut{/tex} limite supérieure de pompage du noeud u pour le pas de temps t, {tex}\small u \in NU{/tex}, {tex}\small t \in NT{/tex}
{tex}\small P^{m+}uet{/tex} limite supérieure de transit dans l'arc artificiel (u, e) pour le pas de temps t, {tex} \small (u,e) \in NAB {\tex}, {tex}\small u \in NU{/tex}, {tex}\small t \in NT{/tex}
{tex}\small R^{m+}kt{/tex} limite supérieure de recharge artificielle de l'aquifère k pour le pas de temps t, {tex}\small k \in NU{/tex}, {tex}\small t \in NT{/tex}
{tex}\small X^{m+}jkt{/tex} limite supérieure de recharge artificielle de l'aquifère k à partir du noeud jonction j pour le pas de temps t, {tex}\small k \in NU{/tex}, {tex}\small j \in N^k{/tex}
, {tex}\small t \in NT{/tex}
{tex}\small m^{m+}kt{/tex} limite supérieure de pompage dans l'aquifère k pour le pas de temps t, {tex}\small k \in NP{/tex}, {tex}\small t \in NT{/tex}
{tex}\small y^{m+}knt{/tex} limite supérieure de pompage de la demande n dans l'aquifère k pour le pas de temps t, {tex}\small k \in NU{/tex}, {tex}\small N \in n_K{/tex}, {tex}\small t \in NT{/tex}

 

Paramètres de niveaux objectifs

{tex}\small \eta _{it}{/tex}
niveau de stockage objectif du réservoir i pour le pas de temps t, {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}
{tex}\small M^+ _{it}{/tex} = { {tex}\small m \in M_i / V^{m+}_{it} > \eta_{it}{/tex}, pour {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}}
{tex}\small M^- _{it}{/tex} = {{tex}\small m \in M_i / V^{m+}_{it} \leq \eta_{it}{/tex}, pour {tex}\small i \in NE{/tex}, {tex}\small t \in NT{/tex}}
{tex}\small d_{nt}{/tex} demande en eau du noeud demande n et pour le pas de temps t, {tex}\small n \in ND \cup NC{/tex}, {tex}\small t \in NT{/tex}
{tex}\small s_{nt}{/tex} volume de stockage réservé objectif dans les réservoirs amont pour satisfaire de futur besoins potentiels au noeud demande n pour le pas de temps t, {tex}\small n \in NS{/tex}, {tex}\small t \in NT{/tex}
{tex}\small q^{m-}_{eft}{/tex} débit réservé associé à l'arc naturel (e, f) pour le pas de temps t, {tex}\small (e,f) \in NAN{/tex} {tex}\small t \in NT{/tex}
{tex}\small r^{m-}_{kt}{/tex} volume de recharge artificielle objectif de l'aquifère k pour le pas de temps t, {tex}\small k \in NR{/tex}, {tex}\small t \in NT{/tex}

 

Paramètres de pénalisation

  

{tex}\small \lambda^+imt, \lambda^-imt{/tex}
coefficient de pénalités sur l'écart positif et négatif par rapport au stockage objectif du niveau de stockage dans la zone m du réservoir i pour le pas de temps t, {tex}\small t \in NT{/tex}, {tex}\small t \in NT{/tex}, et {tex}\small m \in M^+it{/tex}et {tex}\small m \in M^-it{/tex}respectivement
{tex}\small \delta_{nft}{/tex}
coefficient de pénalités sur le déficit de la demande de niveau f dans le noeud demande n pour le pas de temps t, {tex}\small n \in ND \cup NC{/tex}, {tex}\small f \in F_n{/tex}, {tex}\small t \in NT{/tex}.
Si le niveau 1 correspond au niveau le plus proche de la demande objectif, on a {tex}\small \delta_{nlt} \leq \delta_{n2t} ...\delta_{nFnt}{/tex}
{tex}\small \gamma_{nt}{/tex} coefficient de pénalités sur le déficit de stockage réservé dans les réservoirs amont pour satisfaire de futur besoins potentiels au noeud demande n pour le pas de temps t, {tex}\small n \in NS{/tex}, {tex}\small t \in NT{/tex}
{tex}\small \alpha_{eft}{/tex}
coefficient de pénalités sur le déficit de débit réservé relatif à l'arc naturel (e, f) pour le pas de temps t, {tex}\small (e,f) \in NA{/tex}, {tex}\small t \in NT{/tex}
{tex}\small \sigma_{kt}{/tex} coefficient de pénalités sur le déficit de recharge artificielle de l'aquifère k pour le pas de temps t, {tex}\small k \in NR{/tex}, {tex}\small t \in NT{/tex}
{tex}\small \phi_{eft}{/tex} coefficient de cout de transit dans l'arc de surface (e, f) pour le pas de temps t, {tex}\small (e,f) \in NA{/tex}, {tex}\small t \in NT{/tex}
{tex}\small \mu_{knt}{/tex} coefficient de cout de pompage de la demande n dans l'aquifère k pour le pas de temps t, {tex}\small k \in NP{/tex}, {tex}\small n \in NK{/tex}, {tex}\small t \in NT{/tex}
{tex}\small \omega_{it}{/tex}
coefficient de pénalités de déversement au niveau d'une demande non consommatrice i pour le pas de temps t, {tex}\small i \in NC{/tex}, {tex}\small t \in NT{/tex}

  Les variables du modèle

{tex}\small V_{it}{/tex} volume stocké dans le réservoir i pour le pas de temps t, {tex}\small i\in NE {/tex}, {tex}\small t\in NT {/tex}
{tex}\small V^+_{imt}{/tex} volume de stockage en excès par rapport au stockage objectif dans la zone m du réservoir i pour le pas de temps t, {tex}\small m\in M^+_{it}{/tex}, {tex}\small i\in NE {/tex}, {tex}\small t\in NT {/tex}
{tex}\small V^-_{imt}{/tex} mÎvolume de stockage en déficit par rapport au stockage objectif dans la zone m du réservoir i pour le pas de temps t, {tex}\small m\in M^-_{it}{/tex}, {tex}\small i\in NE {/tex}, {tex}\small t\in NT {/tex}
{tex}\small S_{int}{/tex} volume de stockage réservé dans le réservoir i pour satisfaire des besoins futurs potentiels de la demande n pour le pas de temps t, {tex}\small n\in NS {/tex}, {tex}\small i\in NE {/tex}, {tex}\small t\in NT {/tex}
{tex}\small Z_{kt}{/tex} volume stocké dans le réservoir k pour le pas de temps t, {tex}\small k\in NK {/tex}, {tex}\small t\in NT {/tex}
{tex}\small Q_{eft}{/tex} flux d'eau dans l'arc de surface (e, f) pour le pas de temps t, {tex}\small (e,f)\in NA{/tex}, {tex}\small t\in NT {/tex}
{tex}\small X_{jkt} {/tex} volume de recharge artificielle de l'aquifère k à partir du noeud jonction j pour le pas de temps t, {tex}\small k\in NR {/tex}, {tex}\small j\in N^k{/tex}, {tex}\small t\in NT {/tex}
{tex}\small Y_{knt} {/tex} volume pompé de la demande n dans l'aquifère k pour le pas de temps t, {tex}\small k\in NP {/tex}, {tex}\small n\in N_k{/tex}, {tex}\small t\in NT {/tex}
{tex}\small D^-\{nft} {/tex} déficit de demande en eau dans le niveau f du noeud demande n pour le pas de temps t, {tex}\small n\in ND \cup NC{/tex}, {tex}\small f\in F_n{/tex} , {tex}\small t\in NT {/tex}
{tex}\small S^-_{nt}{/tex} déficit sur le stockage réservé pour satisfaire de futurs besoins potentiels de la demande n pour le pas de temps t, {tex}\small n\in NS{/tex}, {tex}\small t\in NT {/tex}
{tex}\small Q^-_{eft}{/tex} déficit sur le débit réservé associé à l'arc naturel (e, f) pour le pas de temps t, {tex}\small (e,f)\in NAN{/tex}, {tex}\small t\in NT {/tex}
{tex}\small R^-_{kt}{/tex} déficit sur la recharge artificielle de l'aquifère k pour le pas de temps t, {tex}\small n\in NS{/tex}, {tex}\small t\in NT {/tex}

 


Réservoir

Description

Un réservoir est considéré comme un noeud de surface possédant une capacité de stockage.

A chaque pas de temps, le volume d'eau stocké dans un réservoir est diminué des pertes par évaporation et infiltration et augmenté par les dessertes des arcs amont, par les apports nets naturels, et éventuellement par réalimentation souterraine. Le volume net est en partie déversé en aval et en partie stocké pour le pas de temps suivant.

Le réservoir peut être découpé en zones où l'on pénalisera plus ou moins des déficits et des surplus de stockage par rapport à une courbe objectif. Il est par exemple ainsi possible de définir une tranche de stockage d'eau de moindre qualité qu'il faut au maximum éviter de solliciter ou de constituer un creux pour l'amortissement des crues.

Formulation

Bilan en eau dans les réservoirs {tex}\small \forall i \in NE{/tex}, {tex}\small t \in NT{/tex}

{tex}\small V ^prime_{it} = V_{it-1}+ \sum _{e \in Q_i}Q_{eit}+ b_{it}{/tex} (1a)
{tex}\small V^\prime{it}- f^e_{it.}V_{it}-f^1_{it}V_{it}+f^\prime_{it.}-i^\prime_{it.} = \sum _{j \in Pi} Q_{ijt} + V_{it}{/tex} (1b)
{tex}\small f^ \prime _{it} = f^f_{kit.}Z_{kt-1}{/tex} si {tex}\small\exists k / i \in J_k {/tex}
        0, sinon
(1c)
et
{tex}\small i^ \prime = f^i_{it.}V_{it}.V_{it}{/tex} si {tex}\small k(i) \not=0{/tex}

        0, sinon
(1d)
et {tex}\small V^{m-}it \leq V_{it} \leq V ^{m+}it{/tex} (2)

Déficits et surplus de stockage par rapport à une courbe objectif

  {tex}\small V_it = \eta _ {it}+ \sum _{m \in M+it} V^+ _{imt} - \sum _{m \in M-it} V^- _{imt} {/tex}
(3)
{tex}\small V^+ _{imt}\leq V^{m++}_{imt} = V^{m+}_{im} - \eta it {/tex}
                               {tex}\small V^{m+}_{imt} - V^{m+}_{im-1t}{/tex}  ,  si {tex}\small \eta _{it} < V^{m+}_{im-1t}{/tex}
{tex}\small \forall m \in M^+_{it}{/tex} (4)
et {tex}\small V^-_{imt} \leq V^{m+-}_{imt} = \eta _{it}-V^{m+}_{imt}{/tex}
                               {tex}\small V^{m+}_{im+1t}-V^{m+}_{imt}{/tex}  ,  si {tex}\small V^{m+}_{im+1t}\leq \eta _{it}{/tex}
{tex}\small \forall m \in M^-_{it}{/tex} (5)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales
# Reservoir en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier du noeud de surface
localisation coordonnées (x, y) du noeud de surface
      Courbe de remplissage
hauteurs nombre de points n (2..12) de la courbe de remplissage
hauteur n valeurs de niveau - unité : m
surface n valeurs de surface du plan d'eau - unité : ha
stockage n valeurs de stockage - unité : hm3
   
apport naturel nom de la chronique d'apport, coefficient correcteur
L evaporation variation mensuelle interannuelle de la lame d'eau évaporée (L=V/S.Cev) - unité : m
pluie directe nom de la chronique d'apport sous forme de lame, coefficient correcteur
Propriétés de gestion
C perte variation mensuelle interannuelle des coefficients de pertes par infiltration profonde - voir équations (1a)-(1b)
infiltration 0 si pas d'infiltration, index a de l'aquifère sinon et dans ce cas nous avons la définition suivante
C infiltration variation mensuelle interannuelle des coefficients de pertes par infiltration dans l'aquifère a
stockage objectif variation mensuelle interannuelle du stockage objectif - unité : hm3
      Découpage en zones
zones nombre de zones considérées M (1..5). Pour chaque zone m, nous avons les 3 définitions suivantes
volume m variation mensuelle interannuelle de la limite supérieure de stockage de la zone m (1..M) - unité : hm3
penalites positives m variation mensuelle interannuelle des pénalités associées au surplus de stockage par rapport à l'objectif de la zone m (1..M)
penalites negatives m variation mensuelle interannuelle des pénalités associées au déficit de stockage par rapport à l'objectif de la zone m (1..M)
   
stockage initial niveau de stockage initial - unité : hm3
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Reservoir
nom Letaye
index noeud 1
localisation 55.30 43.40
hauteurs 9
hauteur 18 19 20 21 22 23 24 25 26
surface 0 1.44 3.53 5.99 8.7 11.63 14.74 18.01 21.43
stockage 0 0.006 0.031 0.078 0.151 0.253 0.385 0.548 0.745
apport naturel LETAYE 1
L evaporation 0.123 0.131 0.162 0.167 0.169 0.17 0.175 0.166 0.142 0.126 0.107 0.109
pluie directe RETENUE1 1
C perte 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04
infiltration 0
stockage objectif 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55
zones 2
volume 0 0 0 0 0 0 0 0 0 0 0 0 0
volume 1 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05
volume 2 0.745 0.745 0.745 0.745 0.745 0.745 0.745 0.745 0.745 0.745 0.745 0.745
penalites positives 1 0 0 0 0 0 0 0 0 0 0 0 0
penalites negatives 1 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
penalites positives 2 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
penalites negatives 2 2.95 2.95 2.95 2.95 2.95 2.95 2.95 2.95 2.95 2.95 2.95 2.95
stockage initial 0.5
sorties global an mois

hgresgc5.jpg
Figure B6.1. - Réservoir - graphique associé


 

Jonction

Description

Une jonction est un noeud de surface où l'écoulement connaît une modification : confluence, prise en rivière, apport intermédiaire, etc.

Pour chaque pas de temps, le volume d'eau entrant à un noeud jonction correspond à la desserte par les arcs amont et éventuellement aux apports nets naturels, à la réalimentation par un aquifère et à la restitution d'une partie d'une demande. Le volume sortant correspond à l'alimentation des arcs aval et éventuellement aux pertes par infiltration et à la recharge artificielle d'un aquifère.

Formulation

Bilan en eau dans les jonctions {tex} \small\forall j{\in}NJ, t{\in}NT {/tex}

{tex}\small F_{jt} = \rho_{jt}^{\prime} + \sum_{e\in Q_j}Q_{ejt} + f_{jt}^{\prime} + b_{jt} {/tex} (6a)
{tex}\small p^\prime jt = p_n \sum _{f \in Fn (d_{nft}+ D^-{nft}){/tex}
si {tex}\small \exists i \in ND / j=j(n){/tex}
        0, sinon
(6b)
et {tex}\small f^\prime jt = f^f_{kjt}.Z_{kt-1} {/tex} si {tex}\small \exists k / i \in J_k{/tex}
        0, sinon
(6c)
et {tex}\small F_{jt}-f^ljt.F_{jt} = \sum _{f \in P_jQ_{jft}+ \chi_{jt}+ i^\prime_{jt}{/tex}

(6d)
{tex}\small \chi_{jt} = X_{jkt}{/tex} si {tex}\small \exists k \in NR / j \in N^k{/tex}
        0, sinon
(6e)
{tex}\small i^\prime = f^i_{jt}.F_{jt}{/tex} si {tex}\small k(j) \not= 0{/tex}
        0, sinon
(6f)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales
# Jonction en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier du noeud de surface
localisation coordonnées (x, y) du noeud de surface
apport naturel nom de la chronique d'apport, coefficient correcteur
Propriétés de gestion
C perte variation mensuelle interannuelle des coefficients de pertes par infiltration profonde - voir équations (6a)-(6d)
infiltration 0 si pas d'infiltration, index a de l'aquifère sinon et dans ce cas nous avons la définition suivante
C infiltration variation mensuelle interannuelle des coefficients de pertes par infiltration dans l'aquifère a - voir équations (6a)-(6f) -
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Jonction
nom vernou
index noeud 3
localisation 17.60 27.21
apport naturel TRAVERSE 0.64
C perte 0 0 0 0 0 0 0 0 0 0 0 0
infiltration 0
sorties global an mois

hgjonvc4.jpg
Figure B7.1. - Jonction - graphique associé


 

 

Pompage

Description

Un pompage est un noeud de surface où une certaine quantité d'eau peut être dérivée vers des arcs artificiels.

Pour chaque pas de temps, le volume d'eau entrant correspond à la desserte par les arcs amont. Le volume sortant correspond à l'alimentation des arcs aval, avec une limite supérieure pour la desserte des arcs artificiels.

Formulation

Bilan en eau dans les pompages et les demandes non consommatrices {tex}\small \forall i \in NCF \cup NU {/tex}, {tex}\small t \in NT {/tex}

   {tex}\small \sum _{(u,e) \in NAB}Q_{uet} \leq p^{m+}ut{/tex}    (7)

Contrainte de limite de pompage {tex}\small \forall u \in NU{/tex}, {tex}\small t\in NT {/tex}

 

{tex}\small \sum _{(u,e) \in NAB}Q_{uet} \leq p^{m+}ut{/tex}
(23)
{tex}\small 0 \leq Q_{uet} \leq p^{m+}_{uet}{/tex}
{tex}\small \forall (u,e) \in NAB{/tex} (24)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales
# Pompage en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier du noeud de surface
localisation coordonnées (x, y) du noeud de surface
Propriétés de gestion
limite variation mensuelle interannuelle de la limite supérieure de pompage - unité : hm3
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Pompage
nom Pompa1
index noeud 1
localisation 5.00 5.00
limite 3 3 3 3 3 3 3 3 3 3 3 3
sorties global an mois

hgpom1c4.jpg
Figure B8.1. - Pompage - graphique associé


 

Consommation en eau

Description

Une consommation en eau est un noeud de surface correspondant à une demande où l'eau est consommée totalement ou partiellement : demande urbaine, demande industrielle, demande d'irrigation, etc.

Pour chaque pas de temps, le volume d'eau entrant correspond à la desserte par les arcs amont et éventuellement au pompage dans un aquifère. Le volume sortant correspond à la restitution éventuelle d'une partie de l'eau au système à partir d'un noeud jonction à l'aval de la demande.

La demande peut être découpée en différents niveaux où les déficits de desserte en eau sont pénalisés de manière croissante.

Formulation

Déficit de demande en eau{tex}\small \forall n \in ND{/tex}, {tex}\small t \in NT{/tex}

{tex}\small y_{nt}+\sum_{e \in Q_n}Q_{ent} = d_{nt}-\sum_{f \in F_n}D^-_{nft}{/tex} (15a)
{tex}\small y_{nt}=Y_{knt}{/tex} si {tex}\small \exists k \in NP / n \in N_k{/tex}
        0, sinon
(15b)
et {tex}\small 0 \leq D^-_{nft} \leq d^\prime _{nft}{/tex} {tex}\small \forall f \in F_n{/tex}, {tex}\small n \in ND{/tex}, {tex}\small t \in NT{/tex} (16)
{tex}\small d^\prime _{nft}{/tex} correspond à la demande en eau du niveau f, donc au déficit maximal pour ce niveau, tel que
{tex}\small \sum _{f\in F_n}d^\prime_{nft}= d_{nt}{/tex} (17)

Les contraintes suivantes ne sont utilisées que pour les modèles d'optimisation prenant en compte un horizon prévisionnel.

Limite de déficits cumulés dans des mois consécutifs {tex}\small \forall n \in ND, t=1,...,|NT|- \tau + 1, \tau \in \Gamma{/tex}

{tex}\small \sum_1 \leq 0 \leq \tau \sum_{{f \in F_n}D^-_{nft + 0-1}}{/tex} (26)
et {tex}\small 0 \leq \sum_{{f \in F_n}D^-_{nft} \leq d_{nt}{/tex} (27)

Volume réservé dans des réservoirs {tex}\small \forall i \in NE{/tex},{tex}\small t \in NT{/tex}

{tex}\small \sum _{n \in P^\prime i \cap NS} S_{int} \leq V_{it} - V^{m-}it{/tex}
(28)

Déficit de stockage réservé dans des réservoirs amont pour satisfaire des besoins futurs potentiels {tex}\small \forall i \in NE{/tex}, {tex}\small t \in NT{/tex}

{tex}\small \sum _{i \in Q^\prime n \cap NE} S_{int} = s_{nt}-S^-_{nt}{/tex} (29)
{tex}\small 0 \leq S^-_{nt} \leq s_{nt}{/tex}
(30)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales

# Consommation en eau en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier du noeud de surface
localisation coordonnées (x, y) du noeud de surface
      Demande en eau
demande variation mensuelle interannuelle de la demande en eau - unité : hm3
source adresse et coefficient de la demande en eau de référence. Cette demande peut être une chronique de débits (ecoulement.nom de la chronique) ou une caractéristique (débit transité dans un lien, demande en eau, etc) d'un autre composant (type.nom du composant). Si elle existe, cette demande se substitue à la variation mensuelle interannuelle.
   
retour 1 entier ; 1 réel : index du noeud de retour ; coefficient de restitution, valeur comprise entre 0 et 1
   
aquifere 0 si pas de pompage, index de l'aquifère sinon et dans ce cas nous avons les 2 définitions suivantes
limite de pompage variation mensuelle interannuelle de la limite supérieure de pompage - unité : hm3
cout de pompage variation mensuelle interannuelle du cout fictif de pompage
Propriétés de gestion
      Découpage en zones
niveaux nombre de niveaux de demandes F (1..5) . Pour chaque niveau f, nous avons les 2 définitions suivantes
taux f variation mensuelle interannuelle de participation du niveau f (1..F) à la demande en eau - valeurs comprises entre 0 et 1 - pour chaque mois Somme(rate 1 .. rate F)=1
penalites f variation mensuelle interannuelle des pénalités associées aux déficits de niveau f (1..F)
   
restriction conditionnement de la desserte à un indicateur, les demandes en eau correspondants aux niveaux associés à des pénalités inférieures ou égales à la pénalité limite ne seront pas desservies si la valeur de l'indicateur est inférieure à la valeur limite définie. L'indicateur peut être lu dans le fichier d'apports ou correspondre au stockage d'un réservoir, au volume transité dans un arc, etc.
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Consommation en eau
nom deshaut
index noeud 26
localisation 43.80 37.00
demande 0.62 0.7 0.775 0.75 0.62 0.6 0.62 0.62 0.6 0.62 0.6 0.62
source - 1
retour 0 0
aquifere 0
niveaux 2
taux 1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
taux 2 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9
penalites 1 4 4 4 4 4 4 4 4 4 4 4 4
penalites 2 5 5 5 5 5 5 5 5 5 5 5 5
restriction 0 0 -
sorties global an mois

hgcondc4.jpg
Figure B9.1. - Consommation en eau - graphique associé


 

Description

Une zone irriguée est une aire de consommation d’eau où le besoin est calculé à partir d’un bilan hydrique prenant en compte la réserve en eau du sol pour une culture.

Pour chaque pas de temps, le volume d'eau entrant correspond à la desserte par les arcs amont et éventuellement au pompage dans un aquifère. Le volume sortant correspond à la restitution éventuelle d'une partie de l'eau au système à partir d'un noeud jonction à l'aval de la demande.

La demande peut être découpée en différents niveaux où les déficits de desserte en eau sont pénalisés de manière croissante.

Bilan hydrique

Le modèle classique simple à réservoir a été retenu pour le calcul du bilan hydrique. On considère ainsi que toute la pluie arrivant sur le sol s'infiltre jusqu'à saturation du sol, qui correspond au remplissage du réservoir dont la capacité constitue la réserve utile (RU). Tout apport d'eau supérieur est considéré comme perdu : la partie correspondant à l'eau gravitaire le sera par percolation, le reste par ruissellement de surface.

Le bilan peut donc s'écrire simplement :

Niveau_reserve_final = Niveau_reserve_initial + Pluie + Dose_apport - ETR
où   ETR : évapotranspiration réelle.

L'évapotranspiration

L'évapotranspiration réelle donne la consommation réelle de la culture en fonction du stress hydrique. Elle est difficile à évaluer, on se contente des résultats semi-empiriques visant à donner sa variation en fonction de l'humidité du sol ou plutôt en fonction de l'état de la réserve.

La réserve utile (RU) est décomposée en réserve facilement utilisable (RFU) et en réserve difficilement utilisable (RDU). RU et RFU dépendent de la profondeur du sol. Trois types de fonctions sont disponibles pour exprimer leurs variations dans l'année.

Tant qu'il y a de l'eau dans la RFU, c'est à dire que la réserve reste supérieure à la RDU, on considère que la plante reste à son niveau de consommation optimal, ETR = ETM. Lorsque la réserve R devient inférieure à RDU, on considère la relation ETR / ETM = R / RDU.

L'ETM, évapotranspiration maximale d'une culture, correspond donc à la consommation du périmètre planté quand celui-ci est convenablement approvisionné en eau.

Dans la pratique, ETM est calculé par ETM = Kc . ETP, où Kc, coefficient cultural, dépend du stade de développement du végétal.

On peut ainsi définir la variation du coefficient cultural au cours de l'année. L'ETP -EvapoTranspiration Potentielle - est définie à partir d'une (ou plusieurs) chroniques d'évaporation.

Demande en eau

Pour assurer une récolte optimale, il faut maintenir la réserve en eau à un niveau supérieur à la RDU. Mais la réserve ne doit pas être forcément pleine.

On a donc introduit la notion de niveau d'irrigation. Ce seuil est exprimé par le pourcentage de remplissage de RFU. Par convention il n'y a pas d'irrigation si la réserve en eau du sol est supérieure à ce seuil d'alerte ou si ce seuil est nul. Ce seuil peut être variable sur l'année.

Le volume de demande en eau peut ainsi s'exprimer :

Vol_demande = Deficit x Surface_irriguee / Efficience_apport

où   Deficit : déficit hydrique par rapport à la réserve pleine

soit   Deficit = max (Niveau_RFU - Niveau_reserve_initial + ETM - Pluie, 0)

où   ETM : évapotranspiration maximale ;
       Efficience_apport : rapport entre volume réellement utilisé et volume d'apport

Formulation

Déficit de demande en eau {tex}\small \forall n \in ND{/tex} , {tex}\small t \in NT{/tex}

{tex}\small y_{nt} + \sum_{{e \in Q_n}Q_{ent}}=d_{nt}- \sum_{{f \in F_n}D^-_{nft}}{/tex} (15a)
{tex}\small y_{nt}=Y_{knt}{/tex}si {tex}\small \exists k \in NP / n \in N^k{/tex}
        0, sinon
(15b)
et {tex}\small 0 \leq D^-_{nft} \leq d^\prime _{nft}{/tex} {tex}\small \forall f \in F_n{/tex}, {tex}\small n \in ND {/tex}, {tex}\small t \in NT{/tex} (16)
{tex}\small d^\prime _{nft}{/tex}correspond à la demande en eau du niveau f, donc au déficit maximal pour ce niveau, tel que
{tex}\small \sum _{{f \in F_n}d^\prime_{nft}}= d_{nt}{/tex} (17)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales

# Zone irriguee en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier du noeud de surface
localisation coordonnées (x, y) du noeud de surface
      Demande en eau
surface surface irriguée - unité : ha
pluie adresses et coefficients de participation des stations pluviométriques (max 6)
evaporation variation mensuelle interannuelle de l'ETP EvapoTranspiration Potentielle - unité : m
Kc coefficient cultural
ru réserve utilisable - unité : m
rfu réserve facilement utilisable - unité : m
prof.rac. profondeur racinaire - unité : m
niveau niveau d'irrigation, valeur comprise entre 0 et 1 - si 0 alors pas d'irrigation - 1 correspond au maintien de la RFU pleine
efficience taux d'apport : efficience * demande = deficit * surface
   
retour 1 entier ; 1 réel : index du noeud de retour ; coefficient de restitution, valeur comprise entre 0 et 1
   
aquifere 0 si pas de pompage, index de l'aquifère sinon et dans ce cas nous avons les 2 définitions suivantes
limite de pompage variation mensuelle interannuelle de la limite supérieure de pompage - unité : hm3
cout de pompage variation mensuelle interannuelle du cout fictif de pompage


Propriétés de gestion

      Découpage en zones
niveaux nombre de niveaux de demandes F (1..5) . Pour chaque niveau f, nous avons les 2 définitions suivantes
taux f variation mensuelle interannuelle de participation du niveau f (1..F) à la demande en eau - valeurs comprises entre 0 et 1 - pour chaque mois Somme(rate 1 .. rate F)=1
penalites f variation mensuelle interannuelle des pénalités associées aux déficits de niveau f (1..F)
   
restriction conditionnement de la desserte à un indicateur, les demandes en eau correspondants aux niveaux associés à des pénalités inférieures ou égales à la pénalité limite ne seront pas desservies si la valeur de l'indicateur est inférieure à la valeur limite définie. L'indicateur peut être lu dans le fichier d'apports ou correspondre au stockage d'un réservoir, au volume transité dans un arc, etc.
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Zone irriguee
nom ban_Pcov
index noeud 147
localisation 27.36 21.01
surface 2000
pluie DUMANOIR 0.24 NEUFCHAT 0.23 DOUVIRFA 0.21 NESTY 0.12 MOREAU 0.12 BOISDEBO 0.08
evaporation 0.108 0.107 0.152 0.137 0.162 0.164 0.149 0.139 0.131 0.149 0.101 0.091
Kc 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1
ru 0.036
rfu 0.012
prof.rac. 1 1 1 1 1 1 1 1 1 1 1 1
niveau 1 1 1 1 1 1 1 1 1 1 1 1
efficience 1
retour 0 0
aquifere 0
niveaux 2
taux 1 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
taux 2 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8
penalites 1 7 7 7 7 7 7 7 7 7 7 7 7
penalites 2 10 10 10 10 10 10 10 10 10 10 10 10
restriction 0 0 -
sorties global an mois

hgzirbc4.jpg
Figure B10.1. - Zone irriguée - graphique associé

 


Utilisation d'eau

Description

Une utilisation d'eau correspond à une demande non consommatrice. Elle associe une demande à un noeud de surface dont l'apport en eau est restitué complètement au système : centrale hydro-électrique, définition de droits d'eau, etc. La demande peut être découpée en différents niveaux où les déficits de desserte en eau sont pénalisés de manière croissante.

Formulation

Bilan en eau dans les pompages et les demandes non consommatrices {tex}\forall i \in NC \cup NU{/tex}, {tex}\small t \in NT{/tex}

   {tex}\small \sum _{{e \in Q_i}Q_{eit} = \sum _{{f \in Pi}Q_{ift}}{/tex}    (7)

Déficit dans les demandes non consommatrices {tex} \small \forall i \in NC {/tex}, {tex}\small t \in NT{/tex}

{tex}\small \sum _{{f \in Pi}Q_{ift}} \geq d_{it} - \sum _{{l \in Fi}D^-_{ilt}}{/tex}
(12)
et {tex}\small 0 \leq D^-{ilt} \leq d^\prime_{ilt}{/tex} {tex}\small \forall l \in F_i{/tex} (13)
{tex}\small d^\prime_{ilt}{/tex} est la limite supérieure de déficit pour le niveau l tel que
{tex}\small \sum _{{l \in Fi}d^\prime_{ilt}=d_{it}}{/tex} (14)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales

# Utilisation d'eau en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier du noeud de surface
localisation coordonnées (x, y) du noeud de surface
      Demande en eau
demande variation mensuelle interannuelle de la demande en eau - unité : hm3
source adresse et coefficient de la demande en eau de référence. Cette demande peut être une chronique de débits (ecoulement.nom de la chronique) ou une caractéristique (débit transité dans un lien, demande en eau, etc) d'un autre composant (type.nom du composant). Si elle existe, cette demande se substitue à la variation mensuelle interannuelle.

Propriétés de gestion  

Découpage en zones

niveaux nombre de niveaux F de découpage de la demande en eau. Pour chaque niveau f, nous avons les 2 définitions suivantes
taux f variation mensuelle interannuelle de participation du niveau f (1..F) à la demande en eau - valeurs comprises entre 0 et 1 - pour chaque mois Somme(rate 1 .. rate F)=1
penalites f variation mensuelle interannuelle des pénalités associées aux déficits de niveau f (1..F)
restriction conditionnement de la desserte à un indicateur, les demandes en eau correspondants aux niveaux associés à des pénalités inférieures ou égales à la pénalité limite ne seront pas desservies si la valeur de l'indicateur est inférieure à la valeur limite définie. L'indicateur peut être lu dans le fichier d'apports ou correspondre au stockage d'un réservoir, au volume transité dans un arc, etc.
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Utilisation d'eau
nom AepInTrav
index noeud 39
localisation 19.03 31.64
demande 2.2 2 2.2 2.1 2.2 2.1 2.2 2.2 2.1 2.2 2.1 2.2
source - 0
niveaux 1
taux 1 1 1 1 1 1 1 1 1 1 1 1 1
penalites 1 2 2 2 2 2 2 2 2 2 2 2 2
restriction 0 0 -
sorties global an mois

hgutiac4.jpg
Figure B11.1. - Utilisation d'eau - graphique associé

 


Aquifère

Description

Un aquifère est considéré comme un noeud souterrain possédant une capacité de stockage.

Pour chaque pas de temps, le volume d'eau entrant à un aquifère correspond aux apports extérieurs naturels, à l'infiltration à partir de réservoirs ou de noeuds de jonction et éventuellement à de la recharge artificielle à partir d'un noeud jonction. Le volume sortant correspond à la réalimentation de réservoirs ou de noeuds de jonction et à des pompages si le volume net de l'aquifère les permet.

Formulation

Bilan en eau dans les aquifères {tex}\small \forall k \in NK{/tex}, {tex}\small t \in NT{/tex}

{tex}\small Z_{kt-1}+ \sum _ {{j \in Nk}/k \in NR X_jkt} + \sum _{{j \in NJ / k(j) = k f^i_{jt}F_jt}} + \sum _ {{i \in NE} / k(i)} = k f^i_{it} V ^\prime_{it}+ r_{kt} = \sum _{{n \in N_k}/k \in NP Y_{knt} + f^\prime _{kt}+ Z_{kt}}{/tex} (8a)
{tex}\small f^\prime _{kt}= \sum _{i \in NE \cup NJ / i \in Jk f^fkitZ_{kt-1}{/tex} si {tex}\small J_k \not = \not \bigcirc{/tex}        0, sinon (8b)
et {tex}\small Z^{m-} \leq Z_{kt} \leq Z^{m+}k{/tex} (9)

Déficit de recharge artificielle {tex}\small \forall k \in NR{/tex}, {tex}\small t \in NT{/tex}

{tex}\small r^{m-}kt - R^-kt \leq \sum j \in Nk Xjkt \leq r^{m+}kt{/tex} (18)
{tex}\small 0 \leq Xjkt \leq x^{m+}jkt{/tex} (19)
et {tex}\small 0 \leq R^-kt \leq r^{m-}kt{/tex} (20)

Limite de pompage dans les aquifères {tex}\small \forall k \in NP{/tex}, {tex}\small t \in NT{/tex}

{tex}\small \sum _{n \in Nk Yknt \leq m^{m+}kt{/tex} (21)
{tex}\small Y_{knt} \leq y^{m+}knt{/tex} (22)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales
# Aquifere en-tête avec le nom du type d'objet
nom label de l'objet
index identificateur entier de l'aquifère
localisation coordonnées (x, y) du noeud de surface
stockage mini capacité minimale de stockage - unité : hm3
stockage maxi capacité maximale de stockage - unité : hm3
apport naturel nom de la chronique d'apport, coefficient correcteur
Propriétés de gestion
      Transfert vers la surface
noeuds sortie nombre de noeuds de surface avec échange d'eau - Pour chaque noeud, nous avons les 2 définitions suivantes
noeud sortie index i du noeud de surface
C sortie variation mensuelle interannuelle des coefficients d'échange relatifs au réservoir ou à la jonction i - voir équations (8a)-(8b)
   
pompage variation mensuelle interannuelle des limites supérieures de pompage - unité : hm3
      Recharge artificielle
recharge totale variation mensuelle interannuelle des limites supérieures de recharge artificielle - unité : hm3
noeuds recharge nombre de noeuds de recharge - Pour chaque noeud, nous avons les 2 définitions suivantes
noeud recharge index j du noeud de recharge
limite recharge variation mensuelle interannuelle des limites supérieures de recharge depuis le noeud jonction j - unité : hm3
recharge objectif variation mensuelle interannuelle des volumes objectifs de recharge artificielle - unité : hm3
penalites recharge variation mensuelle interannuelle des pénalités associées aux déficits de recharge artificielle
   
stockage initial niveau de stockage initial - unité : hm3
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Aquifere
nom Aquif1
index noeud 1
localisation 5.00 5.00
stockage mini 0
stockage maxi 0
apport naturel - 1
noeuds sortie 1
noeud sortie 1
C sortie 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
pompage 0 0 0 0 0 0 0 0 0 0 0 0
recharge totale 0 0 0 0 0 0 0 0 0 0 0 0
noeuds recharge 1
noeud recharge 1
limite recharge 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
recharge objectif 0 0 0 0 0 0 0 0 0 0 0 0
penalites recharge 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
stockage initial 0
sorties global an mois

hgaquic4.jpg
Figure B12.1. - Aquifère - graphique associé

 


Arc de surface

Description

Un arc de surface est un lien qui assure le transit d'eau entre 2 noeuds de surface : un noeud amont et un noeud aval.

Formulation

Limite de transit dans les arcs {tex}\small \forall (e,f) \in NA-NAN{/tex}, {tex}\small t \in NT{/tex}

{tex}\small 0 \leq Q_{eft} \leq q^{m+}_{eft}{/tex} (25)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales
# Arc en-tête avec le nom du type d'objet
nom label de l'objet
noeuds index des noeuds de surface amont et aval
Propriétés de gestion
flux max variation mensuelle interannuelle des limites supérieures de flux d'eau - unité : hm3
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Arc
nom c900_lez
noeuds 8 28
flux max 4.018 3.629 4.018 3.888 4.018 3.888 4.018 4.018 3.888 4.018 3.888 4.018
sorties global an mois

hgarclc4.jpg
Figure B13.1. - Arc de surface - graphique associé

 


Arc naturel

Description

Un arc naturel est un lien qui assure l'écoulement entre 2 noeuds du bassin hydrographique. On lui associe généralement un débit réservé qui peut varier au cours de l'année.

Formulation

Déficit sur le débit réservé {tex}\small \forall (e,f) \in NAN{/tex}, {tex}\small t \in NT{/tex}

{tex}\small q^{m-}_{eft}-Q^-_{eft} \leq Q_{eft}{/tex} (10)
et {tex}\small 0 \leq Q^-_{eft} \leq q^{m-}_{eft}{/tex} (11)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales

# Arc naturel en-tête avec le nom du type d'objet
nom label de l'objet
noeuds index des noeuds de surface amont et aval


Propriétés de gestion

flux max variation mensuelle interannuelle des limites supérieures de flux d'eau - unité : hm3
flux min variation mensuelle interannuelle des débits réservés - unité : hm3
penalites variation mensuelle interannuelle des pénalités associées aux déficits sur les débits réservés
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Arc naturel
nom vernou
noeuds 6 7
flux max 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
flux min 0.375 0.339 0.375 0.363 0.375 0.363 0.375 0.375 0.363 0.375 0.363 0.375
penalites 1 1 1 1 1 1 1 1 1 1 1 1
sorties global an mois

hgarnvc4.jpg
Figure B14.1. - Arc naturel - graphique associé

 


Arc artificiel

Description

Un arc artificiel est un lien qui assure le transit d'eau, entre 2 noeuds de surface, sous une certaine limite : canal, conduite, etc.

Les arcs artificiels peuvent être associés à un pompage qui introduit une autre limite de prélèvement.

Formulation

Limite de transit dans les arcs {tex}\small \forall (e,f) \in NA-NAN{/tex}, {tex}\small t \in NT{/tex}

{tex}\small 0 \leq Q_{eft} \leq q^{m+}_{eft}{/tex}
(25)

Note : Les différents termes sont décrits dans les ensembles d'objets , les paramètres et les variables du modèle.

Définition

Propriétés structurales
# Arc artificiel en-tête avec le nom du type d'objet
nom label de l'objet
noeuds index des noeuds de surface amont et aval
Propriétés de gestion
flux max variation mensuelle interannuelle des limites supérieures de flux d'eau - unité : hm3
sorties type des résultats à sauvegarder. Les chroniques annuelles et les cumuls mensuels sont toujours stockés. La sauvegarde des autres chroniques - mois, décades, jours - est optionnelle.

Exemple

# Arc artificiel
nom Arc a1
noeuds 1 2
flux max 2.491 2.25 2.491 2.411 2.491 2.411 2.491 2.491 2.411 2.491 2.411 2.491
sorties global an mois

hgutiac4.jpg

Figure B15.1. - Arc artificiel - graphique associé

 

FaLang translation system by Faboba
APT Logo fra      logo brgm web frlogo inraeLogo Institut Agro Mpl petit

 

Enregistrer

Enregistrer

Enregistrer

Enregistrer

Enregistrer

Enregistrer