Premium
An indirect assessment on the impact of connectivity of conductivity classes upon longitudinal asymptotic macrodispersivity
Author(s) -
Fiori Aldo,
Boso Francesca,
de Barros Felipe P. J.,
De Bartolo Samuele,
Frampton Andrew,
Severino Gerardo,
Suweis Samir,
Dagan Gedeon
Publication year - 2010
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2009wr008590
Subject(s) - mathematics , univariate , probability density function , function (biology) , space (punctuation) , domain (mathematical analysis) , statistical physics , statistics , physics , multivariate statistics , mathematical analysis , evolutionary biology , biology , linguistics , philosophy
Solute transport takes place in heterogeneous porous formations, with the log conductivity, Y = ln K , modeled as a stationary random space function of given univariate normal probability density function (pdf) with mean 〈 Y 〉, variance σ Y 2 , and integral scale I Y . For weak heterogeneity, the above mentioned quantities completely define the first‐order approximation of the longitudinal macrodispersivity α L = σ Y 2 I Y . However, in highly heterogeneous formations, nonlinear effects which depend on the multipoint joint pdf of Y , impact α L . Most of the past numerical simulations assumed a multivariate normal distribution (MVN) of Y values. The main aim of this study is to investigate the impact of deviations from the MVN structure upon α L . This is achieved by using the concept of spatial correlations of different Y classes, the latter being defined as the space domain where Y falls in the generic interval [ Y , Y + Δ Y ]. The latter is characterized by a length scale λ ( Y ), reflecting the degree of connectivity of the domain (the concept is similar to the indicator variograms). We consider both “symmetrical” and “non‐symmetrical” structures, for which λ ( Y ′) = λ (− Y ′) (similar to the MVN), and λ ( Y ′) ≠ λ (− Y ′), respectively, where Y ′ = Y − 〈 Y 〉. For example, large Y zones may have high spatial correlation, while low Y zones are poorly correlated, or vice versa. The impact of λ ( Y ) on α L is investigated by adopting a structure model which has been used in the past in order to investigate flow and transport in highly heterogeneous media. It is found that the increased correlation in the low conductive zones with respect to the high ones generally leads to a significant increase in α L , for the same global I Y . The finding is explained by the solute retention occurring in low Y zones, which has a larger effect on solute spreading than high Y zones. Conversely, α L decreases when the high conductivity zones are more correlated than the low Y ones. Dispersivity is less affected by the shape of λ ( Y ) for symmetrical distributions. It is found that the range of validity of the first‐order dispersivity, i.e., α L = I Y σ Y 2 , narrows down for non‐symmetrical structures.