S'abonner à un flux RSS
 

Coefficient de ruissellement (HU) : Différence entre versions

De Wikhydro
(21 révisions intermédiaires par un utilisateur sont masquées)
Ligne 1 : Ligne 1 :
 
''<u>Traduction anglaise</u> : Runoff coefficient''
 
''<u>Traduction anglaise</u> : Runoff coefficient''
  
<u>Dernière mise à jour</u> : 3/5/2020
+
<u>Dernière mise à jour</u> : 6/5/2020
  
Coefficient traduisant le rapport entre l’une des caractéristiques de [[Pluie brute (HU)|pluie brute]] tombée sur un [[Bassin versant (HU)|bassin versant]] (en général la hauteur totale ou l'intensité moyenne) et l’une des caractéristiques de [[Pluie nette (HU)|pluie nette]] écoulée à son [[Exutoire (HU)|exutoire]] pendant une durée donnée (le plus souvent la hauteur totale ou l'intensité maximum).  
+
Coefficient traduisant le rapport entre l’une des caractéristiques de la [[Pluie brute (HU)|pluie brute]] précipitée sur un [[Bassin versant (HU)|bassin versant]] et l’une des caractéristiques de la [[Pluie nette (HU)|pluie nette]] écoulée à son [[Exutoire (HU)|exutoire]] pendant une durée donnée.
  
==Résumé==
+
Le concept de coefficient de ruissellement, en particulier en hydrologie urbaine, consiste à supposer qu’à l’échelle d’un élément de surface réceptrice, voire de la totalité du bassin, la pluie nette ou le ruissellement peuvent s’exprimer sous la forme d’une fraction <math>C</math> de la pluie brute. Cette fraction peut être relative :
 
+
* soit aux hauteurs totales de précipitation ou aux volumes mis en jeu au cours d’un événement pluvieux, on parle alors de coefficient « volumétrique » de ruissellement <math>Cv</math> (on parle également dans ce cas de coefficient d'écoulement) ;
La notion de coefficient de ruissellement résulte d’une simplification majeure des processus de transformation de la pluie (pluie brute) en écoulement (pluie nette) sur un élément de surface, éventuellement étendu à l’échelle d’une unité hydrologique ou d’un bassin versant. Cette transformation associe une équation de continuité, traduisant la conservation de la masse (ici des volumes) et une équation dynamique traduisant la mise en mouvement de l’eau sur l’élément sous l’effet des forces en présence (gravité et frottements dus à la rugosité de la surface et à la viscosité de l’eau).
+
* soit aux flux mis en jeu et l’on parle alors de coefficient de ruissellement « instantané » <math>C(t)</math>.
 
+
Le concept de coefficient de ruissellement, en particulier en hydrologie urbaine, consiste à supposer qu’à l’échelle d’un élément de surface réceptrice, voire de la totalité du bassin, la pluie nette ou le ruissellement peuvent s’exprimer sous la forme d’une fraction <math>C</math> de la pluie brute. Cette fraction est relative :
+
* soit aux volumes mis en jeu au cours d’un événement pluvieux, on parle alors de coefficient « volumétrique » de ruissellement <math>Cv</math> (on parle également dans ce cas de coefficient d'écoulement) ;
+
* soit aux flux mis en jeu et l’on parle, alors, de coefficient de ruissellement « instantané » <math>C(t)</math>.
+
 
+
En pratique, l’utilisation d’un coefficient de ruissellement constant présente un intérêt dans les projets d’assainissement pour lesquels les précipitations de projet sont de fréquence faible et autorisent une simplification de la [[Fonction de production et fonction de transfert (HU)|fonction de production]].
+
  
 
==Concepts de base==
 
==Concepts de base==
Ligne 35 : Ligne 29 :
  
  
* pour les coefficients volumétriques de ruissellement, à l'échelle d'un bassin versant on peut écrire :
+
* pour les coefficients volumétriques de ruissellement (que l'on appelle parfois coefficients de production ou coefficients d'apport), à l'échelle d'un bassin versant on peut écrire :
  
  
<center><math>Cv_A = \frac{V_s}{V_b} = \frac{\int(Q_s(t).dt}{\iiint i_b(x,y,t).dx.dy.dt}\quad(4)</math></center>
+
<center><math>Cv_A = \frac{V_s}{V_b} = \frac{\int(Q_s(t).dt}{\iiint {i_b(x,y,t).dx.dy.dt}}\quad(4)</math></center>
  
 
Avec :  
 
Avec :  
Ligne 51 : Ligne 45 :
 
* <math>V_s</math> : volume total sortant du bassin versant pour la pluie.
 
* <math>V_s</math> : volume total sortant du bassin versant pour la pluie.
  
Les relations entre les coefficients locaux (ou élémentaires) de ruissellement (équations (1) et (3)), qui pourraient avoir un sens physique (et être localement mesurés), et les coefficients globaux, à l'échelle d'une unité hydrologique (équations (2) et (4)), ne sont pas évidentes sans hypothèses particulières. En fait, la formulation des coefficients <math>C_A(t)</math> et <math>Cv_A</math> , à l'échelle d'un bassin versant, est essentiellement conceptuelle, et repose sur l'hypothèse d'une analogie de comportement entre la parcelle élémentaire et la totalité du bassin versant.
+
Les relations entre les coefficients locaux (ou élémentaires) de ruissellement (équations <math>(1)</math> et <math>(3)</math>), qui pourraient avoir un sens physique (et être localement mesurés), et les coefficients globaux, à l'échelle d'une unité hydrologique (équations <math>(2)</math> et <math>(4)</math>), ne sont pas évidentes sans hypothèses particulières. En fait, la formulation des coefficients <math>C_A(t)</math> et <math>Cv_A</math>, à l'échelle d'un bassin versant, est essentiellement conceptuelle, et repose sur l'hypothèse d'une analogie de comportement entre la parcelle élémentaire et la totalité du bassin versant.  
  
==Principes de modélisation==
+
==Modélisation des pertes par un coefficient de ruissellement global constant==
  
===Coefficient de ruissellement constant===
+
===Domaine d'utilisation du modèle===
  
C'est la modélisation la plus simple et la plus classique en hydrologie urbaine. Elle semble être adaptée à des unités fortement urbanisées, homogènes, à surfaces actives à peu près constantes. En réalité, cette constance n'est vraie qu'après satisfaction des pertes initiales PI_A. Dans ces conditions, on peut écrire :
+
Déterminer le coefficient volumétrique de ruissellement pour une pluie observée pour laquelle on dispose de mesures représentatives de la hauteur de pluie précipitée sur le bassin versant et du volume écoulé à l'exutoire est relativement simple. La difficulté consiste à prévoir la valeur à attribuer à ce coefficient pour une pluie non encore observée, en particulier si l'on s'intéresse à une autre grandeur que le volume total produit. Il existe par exemple des modèles qui utilisent cette notion pour calculer le débit maximum sortant du bassin versant ([[Méthode rationnelle (HU)|méthode rationnelle]], [[Méthode de Caquot (HU)|méthode de Caquot]], ...). On peut également utiliser une [[Fonction de production et fonction de transfert (HU)|fonction de production]] très simple pour calculer la pluie nette (ou le débit instantané de pluie nette) en fonction de l'intensité instantanée de pluie brute :
  
  
<center><math>Cv_A = \frac{V_s}{V_b} = \frac{V_s}{A.(Hm_A-PI_A)}\quad(5)</math></center>
+
<center><math>i_n(t) = Cv_A.i_b(t)\quad(5)</math></center>
 +
 
 +
<u>Nota</u> : L'intensité de pluie nette, comme le débit de pluie nette obtenu en multipliant l'intensité de pluie nette par la surface du bassin versant, sont des grandeurs intermédiaires virtuelles utilisées en entrée du modèle de transfert pour calculer le débit de sortie du bassin versant.
 +
 
 +
===Introduction de la notion de pertes initiales===
 +
 
 +
Une amélioration simple du modèle consiste à associer au coefficient de ruissellement constant des [[Perte initiale (HU)|pertes initiales]] <math>PI_A</math> qui représentent la hauteur de pluie nécessaire pour que le ruissellement puisse commencer. Ce nouveau paramètre permet de tenir compte du comportement des bassins versants pour les pluies les plus faibles. On peut alors écrire :
 +
 
 +
 
 +
<center><math>Cv_A = \frac{V_s}{V_b} = \frac{V_s}{A.(H_b-PI_A)}\quad(6)</math></center>
 +
 
 +
Soit :
 +
 
 +
 
 +
<center><math>V_s = Cv_A.A.(H_b-PI_A)    \quad(6bis)</math></center>
  
 
Avec :
 
Avec :
 
* <math>Cv_A</math> : coefficient volumétrique de ruissellement ;  
 
* <math>Cv_A</math> : coefficient volumétrique de ruissellement ;  
* <math>V_s</math> : volume total sortant du bassin versant pour la pluie ; <math>V_s=int(Q_s(t).dt}</math> <math>(m^3</math>) ;
+
* <math>V_s</math> : volume total sortant du bassin versant pour la pluie ; <math>V_s=\int{Q_s(t).dt}</math> (<math>m^3</math>) ;
 
* <math>PI_A</math> : perte initiale : hauteur de pluie minimum nécessaire pour que le ruissellement commence (<math>m</math>) ;
 
* <math>PI_A</math> : perte initiale : hauteur de pluie minimum nécessaire pour que le ruissellement commence (<math>m</math>) ;
* <math>Hm_A</math> : Hauteur d'eau totale moyenne précipitée sur le bassin versant (<math>m</math>).
+
* <math>H_b</math> : Hauteur d'eau totale moyenne précipitée sur le bassin versant (<math>m</math>).
 +
 
 +
<u>Nota</u> : En théorie on devrait calculer <math>H_b</math> par une relation de la forme : <math>H_b = \frac{\iiint {i_b(x,y,t).dx.dy.dt}}{A}</math> ; en pratique on se contente en général de faire la moyenne des hauteurs d'eau mesurées sur les différents pluviomètres présents sur le bassin versant (éventuellement en les pondérant par une surface).
 +
 
 +
===Détermination des paramètres dans le cas où l'on dispose de mesures===
 +
 
 +
Si l'on dispose de mesures simultanées de pluie et de débit à l'exutoire du bassin versant, il est possible de tester l'adéquation statistique du modèle linéaire très simple <math>(6)</math> aux observations. Si le modèle est adéquat, les paramètres <math>a</math> et <math>b’=b/a</math> de l'ajustement (''figure 1'') donneront les valeurs moyennes du coefficient de ruissellement <math>CV_A</math> et des pertes initiales <math>PI_A</math>. Si l'on ne souhaite pas introduire le terme correspondant aux pertes initiales, il suffit d'imposer la valeur <math>b=0</math>.
 +
 
 +
 
 +
[[File:coef_ruis_1.JPG|600px|center|thumb|<center>''<u>Figure 1</u> : Détermination des paramètre <math>a</math> et <math>b’</math>.''</center>]]
 +
 
 +
En pratique, dans la plupart des cas, on observe que la courbe n'est pas vraiment linéaire et que le coefficient global de ruissellement a tendance à augmenter avec la hauteur totale de pluie. Il est ainsi possible de considérer des ajustements différents (valeurs différentes des coefficients <math>a</math> et <math>b</math>) par classe de pluies (voir ''figure 2'').
 +
 
 +
 
 +
[[File:coef_ruis_2.JPG|600px|center|thumb|<center>''<u>Figure 2</u> : Valeurs différentes du coefficient de ruissellement selon la hauteur totale de pluie.''</center>]]
 +
 
 +
 
 +
<u>Nota 1</u> : Dans ce cas, seules les valeurs de <math>a</math> et <math>b</math> associées aux pluies les plus faibles peuvent effectivement être associées à une perte initiale et à un coefficient de ruissellement.
 +
 
 +
<u>Nota 2</u> : Les campagnes de mesures sont généralement courtes (de quelques mois à quelques années). Donc, même si on dispose généralement de suffisamment de données pour ajuster correctement les valeurs correspondant aux pluies faibles, voire moyennes, le nombre de pluies fortes, et encore plus, de pluies très fortes, est presque toujours insuffisant pour réaliser un tel ajustement statistique. Il est alors nécessaire d'utiliser les méthodes décrites dans le § suivant.
 +
 
 +
 
 +
===Évaluation des paramètres dans le cas où l'on ne dispose pas de mesures locales===
 +
 
 +
Dans beaucoup de situations on ne dispose pas de mesures, ou en tous cas on ne dispose pas de suffisamment de mesures, et il est alors nécessaire d'évaluer les paramètres <math>Cv_A</math> et éventuellement <math>PI_A</math> en fonction des informations disponibles.
 +
 
 +
====Évaluation du paramètre PI<sub>A</sub>====
 +
 
 +
Même si ce paramètre n'est pas très important pour les pluies fortes, son évaluation correcte est nécessaire dès lors que l'on s'intéresse à des pluies faibles ou moyennes. Chocat (1978) a proposé de relier les pertes initiales à la pente moyenne du bassin versant et à la nature des surfaces par des relations linéaires de la forme :
 +
 
 +
 
 +
<center><math>PI_A = a + b.(3 - I) \quad(7)</math>)</center>
 +
 
 +
Avec :
 +
* <math>I</math> : pente du bassin versant en %.
 +
 
 +
Le ''tableau 1'' indique des valeurs possibles de <math>a</math> et <math>b</math> selon la pente et la nature du sol.
 +
 
 +
 
 +
[[File:coef_ruis_4.JPG|400px|center|thumb|<center>''<u>Tableau 1</u> : Coefficients pour le calcul des pertes initiales ; <u>Source</u> : Chocat (1977).''</center>]]
 +
 
 +
Si on comprend les pertes initiales comme la hauteur de pluie nécessaire au début du ruissellement, ces valeurs semblent surestimer les pertes initiales pour les petites pluies. Certaines surfaces (toits en pente, surfaces imperméables situées à proximité immédiate des avaloirs par exemple) commencent en effet à produire du ruissellement pour des hauteurs d'eau plus faibles.
 +
 
 +
====Évaluation du paramètre Cv<sub>A</sub>====
 +
 
 +
Ce paramètre est le plus important, particulièrement pour les pluies fortes. C'est donc celui qu'il faut prédéterminer avec le plus d'attention.
 +
 
 +
=====Assimilation du coefficient de ruissellement au coefficient d'imperméabilisation=====
 +
 
 +
La façon la plus simple d'estimer <math>Cv_A</math>, préconisée par l'[[Instruction technique (HU)|instruction technique de 1977]], consiste à considérer qu'il est égal au coefficient d'imperméabilisation <math>C_{imp}</math> :
 +
 
 +
 
 +
<center><math>C_{imp} =\frac{A_{imp}}{A}  \quad(8)</math></center>
 +
 
 +
Avec :
 +
* <math>A_{imp}</math> : Surface imperméabilisée.
 +
 
 +
Dans un aménagement traditionnel, avec une grande partie des surfaces imperméables directement connectées au réseau, ce mode de calcul est acceptable pour les pluies moyennes à fortes (typiquement des hauteurs de pluie en quelques heures comprises entre 30 et 60mm). Il majore cependant le coefficient de ruissellement pour les pluies faibles (même en prenant en compte une perte initiale) et peut le minorer pour les pluies extrêmes pour lesquelles la contribution des surfaces perméables peut devenir importante, en particulier dans les zones péri-urbaines.
 +
 
 +
=====Améliorations possibles=====
 +
 
 +
Différentes améliorations ont ainsi été proposées (Schaake ''et al.'', 1967, Viessman ''t al.'', 1977, Kidd, 1979, Desbordes, 1987), la plupart avec des modèles de la forme :
 +
 
 +
 
 +
<center><math>Cv_A = a.C_{imp}+b.I+c \quad(9)</math></center>
 +
 
 +
avec :
 +
* <math>I</math> : pente moyenne du bassin versant.
 +
 
 +
Les valeurs des paramètres  <math>a</math>, <math>b</math> et <math>c</math> peuvent être déterminées par des techniques d'analyse multivariable pour un ensemble de bassins versants donnés pour lesquels on dispose de mesures et ensuite généralisées à des bassins versants supposés avoir le même comportement hydrologique. Elles peuvent également être définies pour des classes de pluies particulières.
 +
 
 +
=====Proposition memento technique ASTEE=====
 +
 
 +
Par exemple, le memento ASTEE (2017) propose de distinguer le comportement du bassin versant pour les pluies courtes et intenses et son comportement  pour les pluies longues, en utilisant les relations <math>(10)</math> et <math>(10bis)</math> suivantes :
 +
 
 +
* Pour les pluies courtes et intenses :
 +
 
 +
 
 +
<center><math>Cv_A = Cr_i.C_{imp}+Cr_p.(1-C_{imp}) \quad(10)</math></center>
 +
 
 +
 
 +
* Pour les pluies longues :
 +
 
 +
<center><math>Cv_A = Ca_i.C_{imp}+Ca_p.(1-C_{imp}) \quad(10bis)</math></center>
 +
 
 +
Avec :
 +
* <math>Cr_i</math> compris entre 0,5 et 0,7 pour les pluies faibles, entre 0,7 et 1 pour les pluies moyennes à forte et valant 1 pour les pluies très fortes ;
 +
* <math>Cr_p</math> valant 0 pour les pluies faibles et moyenne, compris entre 0,1 et 0,4 pour les pluies fortes et tendant vers 1 pour les pluies très fortes ;
 +
* <math>Ca_i</math>  compris entre 0,7 et 1 pour les pluies faibles à forte et valant 1 pour les pluies très fortes ;
 +
* <math>Ca_p</math> valant 0 pour les pluies faibles et moyenne, compris entre 0,1 et 0,4 pour les pluies fortes et tendant vers 1 pour les pluies très fortes.
 +
 
 +
Cette approche constitue une avancée notable par rapport à la démarche consistant à choisir <math>Cv_A</math> égal à <math>C_imp</math> pour toutes les pluies. Il est cependant dommage que le tableau fourni en synthèse simplifie nettement le texte qui l'accompagne en ne retenant que les valeurs maximum des fourchettes des coefficients.
 +
 
 +
=====Pour aller plus loin=====
 +
 
 +
Malgré son intérêt, la méthode précédente n'est cependant pas totalement satisfaisante car elle repose sur différentes hypothèses, non explicitées et généralement non vérifiées :
 +
* <u>Hypothèse 1</u> : les surfaces dites imperméables sont effectivement imperméables : c'est ce qui justifie le choix d'un coefficient <math>Ca_i</math> égal à 1 (même pour les pluies faibles dans le tableau récapitulatif). Or les voiries, qui représentent 10 à 25 % des surfaces considérées comme imperméables, présentent en réalité une [[Capacité d’infiltration (HU)|capacité d'infiltration]] qui varie, selon leur état, de 0,3 à 3 mm/h, c'est à dire qui est d'un ordre de grandeur voisin des intensités que l'on peut avoir pour une pluie longue et peu intense.
 +
* <u>Hypothèse 2</u> : le sol et le réseau d'assainissement constituent les deux seuls exutoires possibles pour une goutte d'eau qui tombe sur un élément de surface donné. Cette hypothèse est également souvent fausse, une partie de l'eau peut en effet s'écouler soit vers un autre exutoire naturel (ruisseau par exemple) soit vers un espace où elle pourra se stocker puis s'infiltrer; il peut s'agir d'un ouvrage spécifiquement conçu pour cette fonction ou plus simplement d'une pelouse ou d'un parc. Par exemple une bonne partie des eaux précipitées sur les accès aux garages des maisons individuelles s'écoulent simplement vers la pelouse voisine et ne rejoignent le réseau d'assainissement que pour des pluies particulièrement intenses.
 +
* <u>Hypothèse 3</u> : l'élément de surface sait, dès le début de la pluie, quelle sera sa durée et son intensité et se comporte en conséquence! En effet, la seule façon d'obtenir un coefficient de ruissellement égal à 1 pour une pluie exceptionnelle, serait que ce coefficient soit égal à 1 dès le début de la pluie...
 +
 
 +
Le fait que ces trois hypothèses ne soient en général pas vérifiées montrent que la proposition précédente conduit encore souvent à majorer les coefficients de ruissellement, et ceci pour toutes les gammes de pluies. Ce n'est pas forcément dommageable dans le cas du dimensionnement des ouvrages car cette majoration va dans le sens de la sécurité, mais ce peut être dommageable dans le cas du diagnostic d'un réseau existant.
 +
 
 +
Pour aller plus loin il paraît nécessaire d'étudier de beaucoup plus près la façon dont les différentes surfaces urbaines sont connectées ou non au système collectif de drainage et contribuent aux écoulements selon les caractéristiques de la pluie. Ceci nécessite des études détaillées de terrain, menées à l'échelle de la ville, voire du quartier, et un couplage avec des mesures effectuées sur des bassins versants représentatifs des différentes situations rencontrées.
 +
 
 +
<u>Nota</u> : La notion de connexion au réseau va probablement devenir centrale dans les années à venir. On parle beaucoup de désimperméabiliser la ville, mais le principal enjeu pour maîtriser les flux, sera surtout de [[Déconnexion des surfaces imperméables (HU)|déconnecter]] le maximum des surfaces imperméables des réseaux, de préférence en gérant l'eau à la parcelle, au plus près de l'endroit où elle arrive au sol.
 +
 
 +
==Conclusions==
 +
 
 +
La modélisation de la fonction de production par de simples coefficients de ruissellement (ou de perte) moyens, à l'échelle d'un bassin versant est très fréquente en hydrologie urbaine. Elle n'est cependant cohérente, au regard de la fonction de transfert, que si cette dernière traduit un processus linéaire non distribué dans l'espace. Ces conditions ne sont approximativement vérifiées que dans les cas, assez fréquents cependant, de bassins versants urbanisés de petite taille, d'occupation des sols homogènes et de surfaces actives peu variables au cours d'une averse, ou d'une averse à l'autre, soumis à des précipitations relativement importantes. Dans les autres cas, l'utilisation de coefficients de ruissellement n'aura que peu d'intérêt, car elle s'éloignera de la réalité physique des processus de production et de transfert du ruissellement.  
  
<u>Nota</u> : En théorie on devrait calculer Hm_A par une relation de la forme : Hm_A = \frac{/iint ib(x,y,t).dx.dy.dt}{A} ; en pratique on se contente en général de faire la moyenne des hauteurs d'eau mesurées sur les différents pluviomètres présents sur le bassin versant (éventuellement en les pondérant par une surface.
+
<u>Bibliographie</u> :  
 +
* ASTEE (2017) : Mémento technique 2017 ; téléchargeable sur : https://www.astee.org/publications/memento-technique-2017/
 +
* Chocat,B. (1978) : Un modèle de simulation des écoulements dans les réseaux d'assainissement pluvial ; thèse Docteur ingénieur ; INSA Lyon ; 304p.
 +
* Desbordes, M. (1987) : Contribution à l'analyse et à la modélisation des mécanismes hydrologiques en milieu urbain ; thèse d'état ; université de Montpellier ; France ; 242 p.
 +
* Kidd, C.H.R. (1979) : Rainfall runoff process over urban surfaces ; Wallingford Institute of hydrology ; Report 53 ; 84 p.
 +
* Schaake, J.C., Geyer, J.C., Knapa, J.W. (1967) : Experimental evaluation of rational method ; journal of the Hydraulics division ; vol 93 ; HY6 ; pp 353-370.
 +
* Viessman, W.V., Knapp, J.W., Lewis, G.I., Harbaugh, T.E. (1977) : Introduction to hydrology ; Ed. Harper and Row ; New-York ; 1977 ; 704 p.
  
 
[[Catégorie:Dictionnaire_DEHUA]]
 
[[Catégorie:Dictionnaire_DEHUA]]
 +
[[Catégorie:Dimensionnement_des_ouvrages_d'écoulement_(HU)]]
 +
[[Catégorie:Dimensionnement_des_ouvrages_de_stockage_(HU)]]
 +
[[Catégorie:Modélisation_de_la_transformation_pluie-débit_(HU)]]

Version du 24 juin 2020 à 16:58

Traduction anglaise : Runoff coefficient

Dernière mise à jour : 6/5/2020

Coefficient traduisant le rapport entre l’une des caractéristiques de la pluie brute précipitée sur un bassin versant et l’une des caractéristiques de la pluie nette écoulée à son exutoire pendant une durée donnée.

Le concept de coefficient de ruissellement, en particulier en hydrologie urbaine, consiste à supposer qu’à l’échelle d’un élément de surface réceptrice, voire de la totalité du bassin, la pluie nette ou le ruissellement peuvent s’exprimer sous la forme d’une fraction $ C $ de la pluie brute. Cette fraction peut être relative :

  • soit aux hauteurs totales de précipitation ou aux volumes mis en jeu au cours d’un événement pluvieux, on parle alors de coefficient « volumétrique » de ruissellement $ Cv $ (on parle également dans ce cas de coefficient d'écoulement) ;
  • soit aux flux mis en jeu et l’on parle alors de coefficient de ruissellement « instantané » $ C(t) $.

Sommaire

Concepts de base

Il existe différentes façons d'exprimer le coefficient de ruissellement selon l'échelle de temps et l'échelle d'espace que l'on considère.

  • pour les coefficients de ruissellement instantanés, à l'échelle d'un élément de surface on peut écrire :


$ C(x,y,t) = 1 - \frac{p(x,y,t)}{i_b(x,y,t)}=\frac{i_n(x,y,t)}{i_b(x,y,t)} \quad(1) $
  • pour les coefficients de ruissellement instantanés, à l'échelle d'un bassin versant on peut écrire :


$ C_A(t) = \frac{Q_s(t)}{\iint i_b(x,y,t).dx.dy} \quad(2) $
  • pour les coefficients volumétriques de ruissellement, à l'échelle d'un élément de surface on peut écrire :


$ Cv(x,y) = 1 - \frac{H_p(x,y)}{H_b(x,y)}=\frac{H_n(x,y)}{H_b(x,y)} \quad(3) $


  • pour les coefficients volumétriques de ruissellement (que l'on appelle parfois coefficients de production ou coefficients d'apport), à l'échelle d'un bassin versant on peut écrire :


$ Cv_A = \frac{V_s}{V_b} = \frac{\int(Q_s(t).dt}{\iiint {i_b(x,y,t).dx.dy.dt}}\quad(4) $

Avec :

  • $ i_b (x, y, t) $ : intensité de pluie brute à l'instant $ t $ au point de coordonnées $ x, y $ ;
  • $ i_n (x, y, t) $ : intensité de pluie nette à l'instant $ t $ au point de coordonnées $ x, y $ ;
  • $ p (x, y, t) $ : pertes au point de coordonnées $ x, y $ et au temps t ;
  • $ Q_s(t) $ : débit à l'exutoire du bassin versant au temps $ t $ ;
  • $ H_p(x,y) $ : hauteur totale de pertes au point de coordonnées $ x, y $ ;
  • $ H_b(x,y $) : hauteur totale de pluie brute au point de coordonnées $ x, y $ ;
  • $ H_n(x,y) $ : hauteur totale de pluie nette au point de coordonnées $ x, y $ ;
  • $ V_b $ : volume total de pluie précipité sur la bassin versant (volume de pluie brute) ;
  • $ V_s $ : volume total sortant du bassin versant pour la pluie.

Les relations entre les coefficients locaux (ou élémentaires) de ruissellement (équations $ (1) $ et $ (3) $), qui pourraient avoir un sens physique (et être localement mesurés), et les coefficients globaux, à l'échelle d'une unité hydrologique (équations $ (2) $ et $ (4) $), ne sont pas évidentes sans hypothèses particulières. En fait, la formulation des coefficients $ C_A(t) $ et $ Cv_A $, à l'échelle d'un bassin versant, est essentiellement conceptuelle, et repose sur l'hypothèse d'une analogie de comportement entre la parcelle élémentaire et la totalité du bassin versant.

Modélisation des pertes par un coefficient de ruissellement global constant

Domaine d'utilisation du modèle

Déterminer le coefficient volumétrique de ruissellement pour une pluie observée pour laquelle on dispose de mesures représentatives de la hauteur de pluie précipitée sur le bassin versant et du volume écoulé à l'exutoire est relativement simple. La difficulté consiste à prévoir la valeur à attribuer à ce coefficient pour une pluie non encore observée, en particulier si l'on s'intéresse à une autre grandeur que le volume total produit. Il existe par exemple des modèles qui utilisent cette notion pour calculer le débit maximum sortant du bassin versant (méthode rationnelle, méthode de Caquot, ...). On peut également utiliser une fonction de production très simple pour calculer la pluie nette (ou le débit instantané de pluie nette) en fonction de l'intensité instantanée de pluie brute :


$ i_n(t) = Cv_A.i_b(t)\quad(5) $

Nota : L'intensité de pluie nette, comme le débit de pluie nette obtenu en multipliant l'intensité de pluie nette par la surface du bassin versant, sont des grandeurs intermédiaires virtuelles utilisées en entrée du modèle de transfert pour calculer le débit de sortie du bassin versant.

Introduction de la notion de pertes initiales

Une amélioration simple du modèle consiste à associer au coefficient de ruissellement constant des pertes initiales $ PI_A $ qui représentent la hauteur de pluie nécessaire pour que le ruissellement puisse commencer. Ce nouveau paramètre permet de tenir compte du comportement des bassins versants pour les pluies les plus faibles. On peut alors écrire :


$ Cv_A = \frac{V_s}{V_b} = \frac{V_s}{A.(H_b-PI_A)}\quad(6) $

Soit :


$ V_s = Cv_A.A.(H_b-PI_A) \quad(6bis) $

Avec :

  • $ Cv_A $ : coefficient volumétrique de ruissellement ;
  • $ V_s $ : volume total sortant du bassin versant pour la pluie ; $ V_s=\int{Q_s(t).dt} $ ($ m^3 $) ;
  • $ PI_A $ : perte initiale : hauteur de pluie minimum nécessaire pour que le ruissellement commence ($ m $) ;
  • $ H_b $ : Hauteur d'eau totale moyenne précipitée sur le bassin versant ($ m $).

Nota : En théorie on devrait calculer $ H_b $ par une relation de la forme : $ H_b = \frac{\iiint {i_b(x,y,t).dx.dy.dt}}{A} $ ; en pratique on se contente en général de faire la moyenne des hauteurs d'eau mesurées sur les différents pluviomètres présents sur le bassin versant (éventuellement en les pondérant par une surface).

Détermination des paramètres dans le cas où l'on dispose de mesures

Si l'on dispose de mesures simultanées de pluie et de débit à l'exutoire du bassin versant, il est possible de tester l'adéquation statistique du modèle linéaire très simple $ (6) $ aux observations. Si le modèle est adéquat, les paramètres $ a $ et $ b’=b/a $ de l'ajustement (figure 1) donneront les valeurs moyennes du coefficient de ruissellement $ CV_A $ et des pertes initiales $ PI_A $. Si l'on ne souhaite pas introduire le terme correspondant aux pertes initiales, il suffit d'imposer la valeur $ b=0 $.


Figure 1 : Détermination des paramètre $ a $ et $ b’ $.

En pratique, dans la plupart des cas, on observe que la courbe n'est pas vraiment linéaire et que le coefficient global de ruissellement a tendance à augmenter avec la hauteur totale de pluie. Il est ainsi possible de considérer des ajustements différents (valeurs différentes des coefficients $ a $ et $ b $) par classe de pluies (voir figure 2).


Figure 2 : Valeurs différentes du coefficient de ruissellement selon la hauteur totale de pluie.


Nota 1 : Dans ce cas, seules les valeurs de $ a $ et $ b $ associées aux pluies les plus faibles peuvent effectivement être associées à une perte initiale et à un coefficient de ruissellement.

Nota 2 : Les campagnes de mesures sont généralement courtes (de quelques mois à quelques années). Donc, même si on dispose généralement de suffisamment de données pour ajuster correctement les valeurs correspondant aux pluies faibles, voire moyennes, le nombre de pluies fortes, et encore plus, de pluies très fortes, est presque toujours insuffisant pour réaliser un tel ajustement statistique. Il est alors nécessaire d'utiliser les méthodes décrites dans le § suivant.


Évaluation des paramètres dans le cas où l'on ne dispose pas de mesures locales

Dans beaucoup de situations on ne dispose pas de mesures, ou en tous cas on ne dispose pas de suffisamment de mesures, et il est alors nécessaire d'évaluer les paramètres $ Cv_A $ et éventuellement $ PI_A $ en fonction des informations disponibles.

Évaluation du paramètre PIA

Même si ce paramètre n'est pas très important pour les pluies fortes, son évaluation correcte est nécessaire dès lors que l'on s'intéresse à des pluies faibles ou moyennes. Chocat (1978) a proposé de relier les pertes initiales à la pente moyenne du bassin versant et à la nature des surfaces par des relations linéaires de la forme :


$ PI_A = a + b.(3 - I) \quad(7) $)

Avec :

  • $ I $ : pente du bassin versant en %.

Le tableau 1 indique des valeurs possibles de $ a $ et $ b $ selon la pente et la nature du sol.


Tableau 1 : Coefficients pour le calcul des pertes initiales ; Source : Chocat (1977).

Si on comprend les pertes initiales comme la hauteur de pluie nécessaire au début du ruissellement, ces valeurs semblent surestimer les pertes initiales pour les petites pluies. Certaines surfaces (toits en pente, surfaces imperméables situées à proximité immédiate des avaloirs par exemple) commencent en effet à produire du ruissellement pour des hauteurs d'eau plus faibles.

Évaluation du paramètre CvA

Ce paramètre est le plus important, particulièrement pour les pluies fortes. C'est donc celui qu'il faut prédéterminer avec le plus d'attention.

Assimilation du coefficient de ruissellement au coefficient d'imperméabilisation

La façon la plus simple d'estimer $ Cv_A $, préconisée par l'instruction technique de 1977, consiste à considérer qu'il est égal au coefficient d'imperméabilisation $ C_{imp} $ :


$ C_{imp} =\frac{A_{imp}}{A} \quad(8) $

Avec :

  • $ A_{imp} $ : Surface imperméabilisée.

Dans un aménagement traditionnel, avec une grande partie des surfaces imperméables directement connectées au réseau, ce mode de calcul est acceptable pour les pluies moyennes à fortes (typiquement des hauteurs de pluie en quelques heures comprises entre 30 et 60mm). Il majore cependant le coefficient de ruissellement pour les pluies faibles (même en prenant en compte une perte initiale) et peut le minorer pour les pluies extrêmes pour lesquelles la contribution des surfaces perméables peut devenir importante, en particulier dans les zones péri-urbaines.

Améliorations possibles

Différentes améliorations ont ainsi été proposées (Schaake et al., 1967, Viessman t al., 1977, Kidd, 1979, Desbordes, 1987), la plupart avec des modèles de la forme :


$ Cv_A = a.C_{imp}+b.I+c \quad(9) $

avec :

  • $ I $ : pente moyenne du bassin versant.

Les valeurs des paramètres $ a $, $ b $ et $ c $ peuvent être déterminées par des techniques d'analyse multivariable pour un ensemble de bassins versants donnés pour lesquels on dispose de mesures et ensuite généralisées à des bassins versants supposés avoir le même comportement hydrologique. Elles peuvent également être définies pour des classes de pluies particulières.

Proposition memento technique ASTEE

Par exemple, le memento ASTEE (2017) propose de distinguer le comportement du bassin versant pour les pluies courtes et intenses et son comportement pour les pluies longues, en utilisant les relations $ (10) $ et $ (10bis) $ suivantes :

  • Pour les pluies courtes et intenses :


$ Cv_A = Cr_i.C_{imp}+Cr_p.(1-C_{imp}) \quad(10) $


  • Pour les pluies longues :
$ Cv_A = Ca_i.C_{imp}+Ca_p.(1-C_{imp}) \quad(10bis) $

Avec :

  • $ Cr_i $ compris entre 0,5 et 0,7 pour les pluies faibles, entre 0,7 et 1 pour les pluies moyennes à forte et valant 1 pour les pluies très fortes ;
  • $ Cr_p $ valant 0 pour les pluies faibles et moyenne, compris entre 0,1 et 0,4 pour les pluies fortes et tendant vers 1 pour les pluies très fortes ;
  • $ Ca_i $ compris entre 0,7 et 1 pour les pluies faibles à forte et valant 1 pour les pluies très fortes ;
  • $ Ca_p $ valant 0 pour les pluies faibles et moyenne, compris entre 0,1 et 0,4 pour les pluies fortes et tendant vers 1 pour les pluies très fortes.

Cette approche constitue une avancée notable par rapport à la démarche consistant à choisir $ Cv_A $ égal à $ C_imp $ pour toutes les pluies. Il est cependant dommage que le tableau fourni en synthèse simplifie nettement le texte qui l'accompagne en ne retenant que les valeurs maximum des fourchettes des coefficients.

Pour aller plus loin

Malgré son intérêt, la méthode précédente n'est cependant pas totalement satisfaisante car elle repose sur différentes hypothèses, non explicitées et généralement non vérifiées :

  • Hypothèse 1 : les surfaces dites imperméables sont effectivement imperméables : c'est ce qui justifie le choix d'un coefficient $ Ca_i $ égal à 1 (même pour les pluies faibles dans le tableau récapitulatif). Or les voiries, qui représentent 10 à 25 % des surfaces considérées comme imperméables, présentent en réalité une capacité d'infiltration qui varie, selon leur état, de 0,3 à 3 mm/h, c'est à dire qui est d'un ordre de grandeur voisin des intensités que l'on peut avoir pour une pluie longue et peu intense.
  • Hypothèse 2 : le sol et le réseau d'assainissement constituent les deux seuls exutoires possibles pour une goutte d'eau qui tombe sur un élément de surface donné. Cette hypothèse est également souvent fausse, une partie de l'eau peut en effet s'écouler soit vers un autre exutoire naturel (ruisseau par exemple) soit vers un espace où elle pourra se stocker puis s'infiltrer; il peut s'agir d'un ouvrage spécifiquement conçu pour cette fonction ou plus simplement d'une pelouse ou d'un parc. Par exemple une bonne partie des eaux précipitées sur les accès aux garages des maisons individuelles s'écoulent simplement vers la pelouse voisine et ne rejoignent le réseau d'assainissement que pour des pluies particulièrement intenses.
  • Hypothèse 3 : l'élément de surface sait, dès le début de la pluie, quelle sera sa durée et son intensité et se comporte en conséquence! En effet, la seule façon d'obtenir un coefficient de ruissellement égal à 1 pour une pluie exceptionnelle, serait que ce coefficient soit égal à 1 dès le début de la pluie...

Le fait que ces trois hypothèses ne soient en général pas vérifiées montrent que la proposition précédente conduit encore souvent à majorer les coefficients de ruissellement, et ceci pour toutes les gammes de pluies. Ce n'est pas forcément dommageable dans le cas du dimensionnement des ouvrages car cette majoration va dans le sens de la sécurité, mais ce peut être dommageable dans le cas du diagnostic d'un réseau existant.

Pour aller plus loin il paraît nécessaire d'étudier de beaucoup plus près la façon dont les différentes surfaces urbaines sont connectées ou non au système collectif de drainage et contribuent aux écoulements selon les caractéristiques de la pluie. Ceci nécessite des études détaillées de terrain, menées à l'échelle de la ville, voire du quartier, et un couplage avec des mesures effectuées sur des bassins versants représentatifs des différentes situations rencontrées.

Nota : La notion de connexion au réseau va probablement devenir centrale dans les années à venir. On parle beaucoup de désimperméabiliser la ville, mais le principal enjeu pour maîtriser les flux, sera surtout de déconnecter le maximum des surfaces imperméables des réseaux, de préférence en gérant l'eau à la parcelle, au plus près de l'endroit où elle arrive au sol.

Conclusions

La modélisation de la fonction de production par de simples coefficients de ruissellement (ou de perte) moyens, à l'échelle d'un bassin versant est très fréquente en hydrologie urbaine. Elle n'est cependant cohérente, au regard de la fonction de transfert, que si cette dernière traduit un processus linéaire non distribué dans l'espace. Ces conditions ne sont approximativement vérifiées que dans les cas, assez fréquents cependant, de bassins versants urbanisés de petite taille, d'occupation des sols homogènes et de surfaces actives peu variables au cours d'une averse, ou d'une averse à l'autre, soumis à des précipitations relativement importantes. Dans les autres cas, l'utilisation de coefficients de ruissellement n'aura que peu d'intérêt, car elle s'éloignera de la réalité physique des processus de production et de transfert du ruissellement.

Bibliographie :

  • ASTEE (2017) : Mémento technique 2017 ; téléchargeable sur : https://www.astee.org/publications/memento-technique-2017/
  • Chocat,B. (1978) : Un modèle de simulation des écoulements dans les réseaux d'assainissement pluvial ; thèse Docteur ingénieur ; INSA Lyon ; 304p.
  • Desbordes, M. (1987) : Contribution à l'analyse et à la modélisation des mécanismes hydrologiques en milieu urbain ; thèse d'état ; université de Montpellier ; France ; 242 p.
  • Kidd, C.H.R. (1979) : Rainfall runoff process over urban surfaces ; Wallingford Institute of hydrology ; Report 53 ; 84 p.
  • Schaake, J.C., Geyer, J.C., Knapa, J.W. (1967) : Experimental evaluation of rational method ; journal of the Hydraulics division ; vol 93 ; HY6 ; pp 353-370.
  • Viessman, W.V., Knapp, J.W., Lewis, G.I., Harbaugh, T.E. (1977) : Introduction to hydrology ; Ed. Harper and Row ; New-York ; 1977 ; 704 p.
Outils personnels