A Magyar Hidrológiai Társaság XXX. Országos Vándorgyűlése (Kaposvár, 2012. július 4-6.)
1. szekció: A vízkárelhárítási szakterület időszerű feladatai - Laurinyecz Pál (KÖVIZIG) - Faragó Gábor László (ADUKÖVIZIG): Az időben változó paraméterű kaszkádmodell alkalmazási lehetőségei
legyen, valamint a Nash-Sutcliffe 1 féle modellhatékonysági tényező konvergáljon, 1-hez. A konstans paraméterű kaszkád eredményei tükrözik a beduzzasztás figyelmen kívül hagyását. A 3. ábrán a futtatás eredménye látható, a négyzetes átlaghiba 25 m 3 s-1, az NSC=0,95. Az optimalizálás során a következő értékeket kaptuk: nopt=2, kopt=1,7 nap-1. 1 A Nash-Sutcliffe féle tényező a hidrológiai modellek értékeléséhez használt mérőszám, amit a következőképpen definiálunk (Nash et al., 1970): E'-1(Q0 - Qm)2 Z;=1(Q0 - _ Ahol: Qo-a megfigyelt vízhozam a t időben, Q*m- a modellezett vízhozam t időben, Q 0 - a megfigyelt vízhozamok átlaga. A együttható értéke -x>, és 1 között változhat, Ha a hatékonyság, NSC=1, akkor a modellezett értékek tökéletesen megegyeznek a megfigyeltekével, amennyiben viszont, NSC <0, abban az esetben a megfigyelt értékek átlaga jobb mint a modell által számított, más szóval a rezídiumok varianciája (számláló), nagyobb, mint a mért értékek varianciája (nevező). Nap 150 100-100 0------50-----50 0 200 400 600 800 1000 1200 1400 1600 1800 Qcsongrád[m3s-1] 3. ábra: Konstans paraméterű kaszkád eredményei, a Tisza vízhozam idősora valamint a hibák korreláció vizsgálata a tiszai vízhozamokkal Figure 3: Results of DLCM with constants storage coefficient and the discharges of Tisza. Regression between DLCM errors and Tisza flow Ahogy láthatjuk a 3. ábrán a kaszkád eredményei megfelelőnek mondhatóak a kis- középvizes periódusokban, azonban árhullámok idején az eredmények egyértelműen hibásak. Regressziós vizsgálatokra használtuk föl az adódó hiba idősort ennek eredményeképpen szoros összefüggést kaptunk a Tisza 600 m 3 s-1 fölötti csongrádi szelvényének vízhozamai és a kaszkád hibái között. Vegyük észre, hogy kopt-1=K0 vagyis a hullám átlagos átvonulási ideje a folyószakaszon, ami értelemszerűen csökken mikor a Tiszán levonuló árhullám visszaduzzasztja a Köröst. Ezután az átvonulási idő felújításának összefüggését a következő hatvány alakban kerestük: K(Q(t) - Qcnt) = a(Q (t) - Qcnt)b + K0 [15] ahol: K - hullám átvonulási ideje [T], K0 - hullám átvonulási ideje visszaduzzasztás nélküli esetben [T], a - optimalizálandó paraméter [L 2 T-3], b - optimalizálandó paraméter [-], Q(t) - Tisza vízhozamai, Qcrit - Vízhozam küszöbérték NSC = 1 Q0) 2 5
