13.07.2015 Views

Assistance au calage de modèles numériques en hydraulique ... - TEL

Assistance au calage de modèles numériques en hydraulique ... - TEL

Assistance au calage de modèles numériques en hydraulique ... - TEL

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

3.2 MÉTHODES D’AJUSTEMENT DES PARAMÈTRESfonctions coût ont <strong>de</strong> fait été testées dans la littérature hydrologique pour comparer<strong>de</strong>ux hydrogrammes. Diskin et Simon (1977) et surtout Gre<strong>en</strong> et Steph<strong>en</strong>son (1986)prés<strong>en</strong>t<strong>en</strong>t un inv<strong>en</strong>taire critique <strong>de</strong> ces différ<strong>en</strong>ts critères. Parmi ceux-ci, on peut releverl’effici<strong>en</strong>ce du modèle proposé par Nash et Sutcliffe (1970) d’après une analogiestatistique avec le coeffici<strong>en</strong>t <strong>de</strong> détermination (Q o représ<strong>en</strong>te la moy<strong>en</strong>ne <strong>de</strong>s débitsobservés) :R 2 =( p∑i=1)( ) 2 (2 ∑ pQ i o − Q o −( p∑i=1i=1)( ) 2 2Q i o − Q i c(Q i o − Q o) 2) 2(3.14)La société HR WALLINGFORD a lancé il y a 10 ans <strong>de</strong> cela une étu<strong>de</strong> sur l’optimisation<strong>en</strong> hydr<strong>au</strong>lique fluviale, et notamm<strong>en</strong>t sur les fonctions coût à utiliser. L’utilisationd’une fonction alambiquée – que nous ne reproduirons pas ici – a ainsi été préconiséepour t<strong>en</strong>ir compte du dé<strong>calage</strong> temporaire év<strong>en</strong>tuel <strong>de</strong>s hydrogrammes (Morriset Anastasiadou-Parth<strong>en</strong>iou, 1994; Lavedrine et Anastasiadou-Parth<strong>en</strong>iou, 1995).Anastasiadou-Parth<strong>en</strong>iou et Samuels (1998) ont finalem<strong>en</strong>t mis <strong>en</strong> œuvre la fonctionsuivante :où :J = 1p − 2p−1∑i=2[1 −[ min(0, Ki ) ] 4 ] ( Q i o − Q i c) 2(3.15)K i =[1 +Q i+1o− 2 Q i o + Q i−1o(∆t) 2( Qi+1o− Q i−1o2 ∆t)] 3/2Comparer <strong>de</strong>s lignes d’e<strong>au</strong>Le <strong>de</strong>uxième cas correspond à l’utilisation <strong>de</strong> données internes <strong>au</strong> domaine modélisé: soit <strong>de</strong>s nive<strong>au</strong>x d’e<strong>au</strong> mesurés le long du cours d’e<strong>au</strong> à <strong>de</strong>s instants précis, soit <strong>de</strong>slaisses <strong>de</strong> crue. Des fonctions coût différ<strong>en</strong>tes <strong>de</strong> celles utilisées <strong>en</strong> hydrologie sont doncà déterminer. Deux fonctions transparaiss<strong>en</strong>t <strong>de</strong> la littérature : le critère <strong>de</strong>s moindrescarrés et le critère minimax.Le critère <strong>de</strong>s moindres carrés est le plus couramm<strong>en</strong>t utilisé dans ce cas <strong>de</strong> figure :J =p∑i=1(h i c − h i o) 2(3.16)Ce critère représ<strong>en</strong>te un écart moy<strong>en</strong> <strong>en</strong>tre valeurs observées et calculées. À ce titre, cecritère n’est pas très s<strong>en</strong>sible <strong>au</strong>x points aberrants ou horsains dès lors que le nombretotal <strong>de</strong> points est suffisant. De nombreuses variantes <strong>de</strong> ce critère ont été employées <strong>en</strong>hydr<strong>au</strong>lique fluviale :– normalisé, soit par les valeurs observées (Ramesh et al., 2000), soit par leurmoy<strong>en</strong>ne temporelle (B<strong>en</strong>nett, 1975) ;72

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!