wolfhece.eva.hydrogramme_mono

Hydrogramme Synthétique MonoFréquance - HSMF Mono Frequency Synthetic Hydrograph - MFSH

@author : Pierre Archambeau - ULiege - HECE @date : 2023

Module Contents

wolfhece.eva.hydrogramme_mono._objectif(x, hyd, nbsteps, dt, istart, iend, vobj)

Fonction objectif à minimiser : écart entre volume attendu (vobj) et calculé (maximum sur base d’une convolution sur l’ensemble de l’hydrogramme) variable : débit à l’heure courante (x)

On ajoute tout de même la convolution de la fin d’intervalle dans la fonction obj car si le vrai maximum est indépedant de x,

l’algo d’optimisation pourrait ne pas converger correctement, autrement cela ne change pas la position de l’optimum

class wolfhece.eva.hydrogramme_mono.Hydro_HSMF(Q, durees, temps_montee, dt, label='')

_init_HSMF()

plot()

_plot_Q(ax: matplotlib.pyplot.Axes, decal=0, label=None, lw=2, xup_bound=None)

opt_hydro(opti_method='brent')

Calcul de l’hydrogramme par optimisation progressive

Hypothèse : les débits sont donnés avec une résolution maximale horaire

Méthode :

• la phase d’initialisation crée une montée en crue linéaire sur base du temps de montée et du premier débit horaire

• on recherche ensuite le débit suivant:

  • suppoosé constant sur une durée correspondant à l’intervalle de discrétisation de durée

  • par résolution d’un problème de minimisation de l’écart de volume vis-à-vis du volume théorique

  • le volume peut apparaître n’importe où dans l’hydrogramme entre 0 et iend (recherche par np.max du produit de convolution)

• une fois le problème résolu, on incrémente l’espace par identification de l’index du maximum

Constat :

• il n’y a aucune garantie que l’hydrogramme total respecte le temps de base

• le temps de base réel peut en effet être plus court si au moins un volume n’est pas totalement situé sur un inervalle d’intégration à droite du pic

• plus le temps de montée est long, plus le point précédent sera rencontré

wolfhece.eva.hydrogramme_mono.dt = 1