A Magyar Hidrológiai Társaság IV. Országos Vándorgyűlése III. kötet, Vízminőség-védelem - Környezetvédelem (Győr, 1983. június 29-30.)
III./1. - Dr. Székely Ferenc: Hidrogeológiai védőidomok meghatározásának modellezési kérdései
lanatszerűnek/ tekintjük, akkor a szennyeződés frontszerüen, de a szivárgási frontnál lassabban mozog. Ez azt jelenti, hogy az adszorbeálódó szennyeződésre vonatkozó elérési idő megnő: ahol d = C q/G o - az u.n. egyensúlyi koncentrációeloszlási együttható, G q - az a feltételezett koncentráció, amelyet a C q folyadékkoncentráció mellett adszorbeálódott szennyezőanyag okozna a folyadékfázisban feloldva. mint a t izochron felület által határolt térfogat. Az /l/ és /2/ képletek gyakorlati alkalmazásának megkönnyítése céljából PORTRAH-IV nyelven irt számítógépi programot készítettünk. A program alkalmazását a 2. ábrán vázolt példán keresztül mutatjuk be. Határozzuk meg a trajektóriákat és a rétegirányu t n , valamint a függőleges szivárgás idejével megm ^ növelt t'~ szivárgási idő izochron felületeit egy C=2ooo m oldalhosszúságú, vizzáró peremekkel határolt, négyzetalaku vizvezető rétegben, amely felülről egy féligáteresztő rétegen keresztül járulékos utánpótlódást kap a zérus depresszióju táplálórétegből /leaky aquifer/. A vizvezető réteg transzmisszibilitása T=loo nf/nap, szivárgási tényezője K=loo n/nap, porozitása n =o,l, a féligáteresztő réteg vastagsága m=lo m, * —3 függőleges irányú szivárgási tényezője K^=lo m/nap, porozitása n Q=o,l. A feladatot rácsmódszerrel oldottuk meg. A réteget lefedő blokkok középpontjaiban felvett rácspontok x és y tengely menti távolsága egyformán loo m volt /ekvidisztans rács/, a kut a 2ox2o=4oo blokkból álló modell /1,1/ koordinátájú bal felső blokkjában üzemelt Q^ ^=looo rr/nap hozammal. A permanens depresszióeloszlást numerikus szimulációval hatáTü / roztuk meg. Az ábrán néhány trajektóriát, valamint a t és m—1 t ~ izochrónokat tüntettük fel a táplálóréteg állandó és vizszintes helyzetű piezometrikus nyomásszintjének feltételezésével. Ekkor KHV a réteg egész területére kiterjed, BHV viszont a t„ izochron görbén belüli terület. /6/ alapján meghatározott védőidom kisebb lesz 53
