Réservoir linéaire (modèle du) (HU) : Différence entre versions

Version du 8 janvier 2020 à 15:34

Traduction anglaise : Linear reservoir

Dernière mise à jour : 8/1/2020

Fonction de transfert très simple largement utilisée en hydrologie urbaine pour représenter la transformation d'un hydrogramme de pluie nette en hydrogramme à l'exutoire.

Formulation mathématique

Il s’agit d’un modèle « à contrôle aval » qui repose sur l’équation de continuité (1) et sur une équation de stockage (2) :

  Avec :

• $ Vs(t) $ : stockage à l’instant t ;
• $ Qs(t) $ : débit à l’exutoire à l’instant t ;
• $ Qe(t) $ : débit entrant ;
• $ K $ : paramètre du modèle, homogène à un temps.

La combinaison de ces deux équations conduit à une équation différentielle du premier ordre, dont la solution pour $ K $ constant est :

Avec

• $ Q_s(0) $ : débit à $ t = 0 $ résultant par exemple d’un écoulement précédent ;
• $ Q_b $ : débit de base permanent éventuel.

D’un point de vue théorique, le paramètre $ K $ correspond au décalage temporel des centres de gravité de $ Q_e(t) $ et $ Q_s(t) $. D’après l’équation (3), l’hydrogramme unitaire instantané du modèle est :

Le maximum de l’opérateur $ O(t) $ se produit à $ t=0 $ (Voir figure 3) et l’on peut donc s’attendre à ce que le modèle réponde plus rapidement que la réalité. D’un point de vue physique, l’équation de l’hydrogramme unitaire instantané est analogue à la décharge transitoire d’un condensateur électrique. C’est la raison pour laquelle le modèle agit comme un filtre « passe-bas » des pluies entrées du modèle en amortissant les composantes de hautes fréquences de ces pluies.

Figure 3 : Représentation de la fonction $ O(t) $.

Utilisation pratique

La fonction $ Q_e(t) $ n'est généralement pas connue sous une forme analytique ; elle doit donc être approchée par une fonction simple pour permettre l'intégration du modèle. La démarche la plus courante consiste à la représenter par une fonction en escalier, discrétisée sur un pas de temps constant. La figure 4 présente une telle décomposition.

Figure 4 : Discrétisation du débit à l'entrée.

Sur le plan analytique, la fonction $ Q_e(t) $ peut alors être écrite comme une succession de $ n $ segments horizontaux de droites :

$ Q_e(t) $ = $ Q_{e1} $ si $ 0 ≤ t < Δt $

$ Q_e(t) $ = $ Q_{e2} $ si $ Δt ≤ t < 2.Δt $

...

$ Q_e(t) $ = $ Q_{ei} $ si $ (i-1).Δt ≤ t < i.Δt $

...

$ Q_e(t) $ = $ Q_{en} $ si $ (n-1)Δt ≤ t < n.Δt $

Il est alors facile d'intégrer $ Q_s(t) $ :

soit finalement :

En posant $ Q_{si} = Q_{s(i.Δt)} $, on obtient ainsi une formule de récurrence facile à utiliser :

Valeur du paramètre K

De nombreuses relations empiriques ont été proposées pour déterminer le paramètre K pour des bassins versants non jaugés. En France, on utilise généralement divers ajustements proposés par Desbordes (1974) comme, par exemple :

$ K = 5,93.A^{0,441} $

$ K = 0,494.A^{-0,0076}.IMP^{-0,512}.I^{-0,401}.L^{0,608} $

$ K = 3,55.A^{0,27}.(1 + IMP)^{-1,9}.I^{-0,36}.L^{0,15}.D_p^{0,21}.H_p^{-0,07} $

avec :

• $ A $ : surface du bassin versant (hectares) ;
• $ D_p $ : durée de la période de "pluie critique" du bassin (de l'ordre de grandeur du temps de réponse du bassin)(mn) ;
• $ H_p $ : hauteur de pluie pendant cette durée (mm) ;
• $ I $ : pente du plus long parcours (%) ;
• $ IMP $ : coefficient d'imperméabilisation (%) ;
• $ L $ : longueur du plus long parcours de l'eau (collecteur principal) (m).

Les équations de prédétermination de K peuvent être utilisées au stade des projets ou des simulations de réseaux existants à condition que le bassin ne comporte pas d’ouvrages importants de stockage des eaux pluviales.

Conclusions générales sur le modèle

Le modèle du réservoir linéaire, sous sa forme la plus simple, apparaît bien adapté à des situations assez variables allant de petits bassins versants fortement urbanisés à des bassins versants de plusieurs centaines d'hectares. La surprenante aptitude d'un modèle linéaire très simple à représenter des phénomènes complexes et fortement non linéaires semble avoir des origines multiples.

• On peut par exemple avancer que la combinaison de transformations élémentaires non linéaires, spatialement distribuées, et isolées par des singularités, peut conduire à une transformation globale approximativement linéaire. Dans ce cas, le modèle pourrait être alors physiquement adéquat du moins dans certaines conditions.
• On peut aussi penser que les bonnes performances du modèle peuvent être dues à la complexité même des systèmes urbains de collecte et d'évacuation des eaux pluviales. Les contraintes urbanistiques, techniques, financières, voire administratives, dans la réalisation des réseaux semblent en effet conduire à une certaine régularité de la géométrie des surfaces urbaines et des réseaux de collecte [Thibault, 1987].

En tout état de cause, l'utilisation de ce modèle pour simuler des bassins versants urbanisés, équipés de réseaux d'assainissement relativement développés, et dont la surface est comprise entre quelques hectares et quelques centaines d'hectares paraît totalement justifiée (O'Loughlin et al, 1996)

Bibliographie :

• Desbordes M. (1974) : "Réflexions sur les méthodes de calcul des réseaux urbains d'assainissement" ; thèse DI ; Université des Sciences et Techniques du Languedoc ; Montpellier ; 171p. ; 1974.
• O'Loughlin, G., Huber, W. and Chocat, B. (1996) : Rainfall-Runoff Processes And Modelling; Journal of Hydraulic Research, Vol. 34, 1996, No. 6, pp. 733-751
• Thibault S. (1979) : Eléments pour une phénoménologie en hydrologie urbaine ; thèse Docteur ingénieur ; INSA Lyon -France ; 361 p. ; 1979.