Vízügyi Közlemények, 1988 (70. évfolyam)
1. füzet - Székely Ferenc: Szivárgási és advekciós transzportfolyamatok numerikus modellezése rétegzett hidrogeológiai rendszerben
30 Székely Ferenc Numerical modeling of seepage and advective transport prjeesses in stratified hvdrogeological systems by Dr. Ferenc SZÉKELY, M. E. hydrogeologist In case of si rali tied la \ ors 1 lie non-steady slate flow of subsurface waters has been calculated by the author by use of vertical leakage theory the changes in the concentration of solutions - on the other hand - are handled with due consideration on the parameters of advection, of sorption and of the first-order (e.g. radioactive) decay processes defining the FSRADA processes or models. Under certain circumstances a substantial role is played by the phenomenon of vertical mixing in the aquifers described by another model called FSRADAM. The mathematical model for the simulation of the FSRADAM processes consisted of Eqs. 1 and 2 in layers with low permeability and identified as к = 1 . 3 ni- I . and of Eqs. 3 to 9 in the aquifers identified as к = 2, 4,..., ni. The investigated flow field was subdivided by the author by use of vertical planes drawn parallel to the x and y axes into rectangle grid-elements where the top and bottom of elements follow the middle of thickness of aquitards. The half-points of the vertical edges of the grid-elements of the water bearing layers are called "nods" or "grid-points". Another important geometric element is the nodal-block which is composed of maximum 4 quartes of grid elements around the nods (Fig. 2). The H-function appeared as a second-order (parabolic hyperboloid) distribution in the gridelements of the water conducting layers. The equipotential curves (10) consisted of hyperbolas with low permeability the H-function was linear in the vertical direction. Formally, the distribution of pressure could be regarded as trilinear. The simulation of concentrations, in other words, the calculation of concentrations in a given block has been executed by the author by use of the so-called piston flow model set up for the description of local advection. For practical application a demonstrative example is presented by the author (Fig. 3). Results are shown in Tables I and II. # # # La modélisation numérique des processus d'infiltration et de transport advectif dans le système hydrogéologique stratifié par Dr. SZÉKELY F., hydrogéologue Pour calculer les caractéristiques de l'écoulement non-permanent des eaux souterraines, l'auteur utilise la théorie d'infiltration verticale qui est le plus férquemment appliqée. LA variation de la concentration des solutés est calculée à partir de l'advecti on, de la sorption et de la décomposition primaire (par example: la décomposition radioactive). L'ensemble de ces processus est considéré comme un processus FSRADA ou un modèle. Dans certaines conditions, le mélangement vertical dans les aqufères joue également un rôle important. Le pi ocessus du mélangement est décrit par un modèle dit FSRADAM. Le modèle mathématique de la simulation du processus FSRADAM est constitué par les équations de 1 à 2 cncernant les couches de k= 1,3,... ni1 qui sont faiblement perméables. Les équations de 3 à 9 concernent les couches aquifères de к = 2,4, ...ni (Figure I). Au cours de la simulation, l'auteur décompose le domaine d'infltration étudié par les plans verticaux parallèles aux axes x et y pour avoir des éléments de grille rectangle, ai nsi la surface de dédoublement relie les centres d'épaisseur des couches à faible perméabilité.
