3D numeričkim modelom provedena je numerička analiza dinamike pronosa slane vode u koritu rijeke Neretve na dionici 0+000.0 m - 25+000.0 m. Osim postojećeg stanja korita rijeke analizirane su i hipotetske situacije sa izvedbom jedne brane u profilu 2+000.0 m (velika dubina i mala širina) ili 4+000.0 m (mala dubina i velika širina), te izvedbom 4 brane u profilima 2/4/6/8+000.0. Regulacija podizanja i spuštanja brane vezana je uz srednji izmjereni protok na mjernoj postaji Metković u prethodnom razdoblju od 24 h (QSR-24) ili 6 h (QSR-6). U situacijama QSR-24 > 300 i 350 m3/s, odnosno QSR-6 > 300 m3/s, brana je potpuno spuštena, dok se za manje srednje vrijednosti protoka brana podiže do određene razine (od -2 m do -5 m) kroz jedan sat. Korišten je 3D numerički model Mike 3fm, zasnovan na metodi konačnih volumena. Za kalibraciju numeričkog modela korišteni su izmjereni protoci, razine i salinitet na monitoring postajama u 2018. godini, dok su za validaciju modelskih rezultata korišteni podaci iz 2017. godine. Podaci o temperaturi i salinitetu mora u akvatoriju ušća Neretve preuzeti su iz nacionalnih programa praćenja morskog okoliša. Modelska dinamika uzvodno-nizvodne translacije slanog klina u nestacionarnim uvjetima strujanja u suglasju je s mjerenjima. Modelska debljina razdjelne zone između slatke i slane vode veća je od izmjerene, kao posljedica numeričke difuzije i usvojene proračunske diskretizacija u vertikalnom smjeru od 1 m. Mjerene vrijednosti i rezultati modeliranja upućuju na zaključak da pri nestacionarnim protocima većim od 700 m3/s dolazi do potpunog istiskivanja slanog klina iz riječnog korita, dok pri protoku manjem od 150 m3/s slani klin napreduje do uzvodne stacionaže 20+909.0 m (Metković) i dalje uzvodnije. Ukoliko se želi spriječiti intruzija slanog klina uzvodno od stacionaže brane (2+000.0 m ili 4+000.0 m) potrebno je izvesti branu s podizanjem do dubine -2 m, uz uvjet podizanja brane QSR-24 300 and 350 m3/s and QSR-6 > 300 m3/s, the dam was fully lowered whereas at smaller mean discharge values the dam is lifted up to a certain level (from -2 m to -5 m) during one hour. The 3D numerical model Mike 3fm, based on the finite volume method, was implemented. The discharges, levels and salinity measured at monitoring stations in 2018 were used for the model calibration and the 2017 data for the model results validation. Data on the sea temperature and salinity in the waters of the Neretva river mouth were taken from the national marine environment monitoring programmes. The model dynamics of the upstream-downstream translation of the salt “wedge” in non-stationary circulation conditions is in line with the measurements. The model thickness of the dividing zone between the freshwater and saltwater is larger than the one measured as a result of the numerical diffusion and adopted discretization of the calculation in the vertical direction of 1 m. The measured values and modeling results point to the conclusion that a complete extrusion of the salt “wedge” from the riverbed occurs at non-stationary discharges larger than 700 m3/s whereas the salt “wedge” advances upstream to 20+909.0 m (Metković) and further at discharges lower than 150 m3/s.
If the intrusion of the salt “wedge” upstream of the dam (2+000.0 m or 4+000.0 m) is to be prevented, a dam with a lifting to the depth of -2 m, under the condition of dam lifting at QSR-24 300 und 350 m3/s, bzw. QSR-6 > 300 m3/s, ist die Schleuse geschlossen, während im Falle von niedrigeren mittleren Durchflusswerten die Schleuse bis zu einer bestimmten Höhe (von -2 m bis -5 m) auf eine Stunde geöffnet wird. Das numerische 3D-Modell Mike 3fm, das auf der Methode der finiten Elemente gründet, wurde eingesetzt. Zur Kalibrierung des numerischen Modells wurden die in 2018 an Überwachungsstationen gemessenen Werte von Durchflüssen, Wasserständen und Salinität genutzt, während zur Validierung der Modellergebnisse die Angaben aus dem Jahr 2017 genutzt wurden. Die Angaben zur Temperatur und Salinität an der Mündung der Neretva wurden aus den nationalen Programmen zur Überwachung der Meeresumwelt übernommen. Die Modelldynamik des Vordringens des Salzkeils flussaufwärts und flussabwärts bei der nichtstationären Strömung stimmt mit den Messungen überein. Die Modelldichte der Übergangszone zwischen Süßwasser und Salzwasser ist höher als die gemessene Dichte, was als Folge der numerischen Diffusion und angenommenen Diskretisierung in der vertikalen Richtung von einem Meter gesehen wird. Die gemessenen Werte und die Ergebnisse der Modellierung weisen auf die Schlussfolgerung hin, dass bei nichtstationären Durchflüssen, die großer als 700 m3/s sind, zu einer völligen Verdrängung des Salzkeils aus dem Flussbett kommt, während bei dem Durchfluss unter 150 m3/s der Salzkeil bis zur Stationierung 20+909.0 m (Metković) und weiter flussaufwärts vordringt. Wenn man die Intrusion des Salzkeils flussaufwärts von der Stationierungsschleuse (2+000.0 m oder 4+000.0 m) verhindern will, sollte eine Schleuse gebaut werden, die bis zur Tiefe -2 m unter der Bedingung der Schleusenöffnung QSR-24 < 300 m3/s geöffnet wird
