Modelagem de movimento horizontal de água no solo utilizando a teoria dos fractais

A matemática fractal tem sido utilizada para caracterizar o transporte de água e solutos em meios porosos e também para simular características físicas e geométricas de meios porosos. O objetivo deste trabalho foi avaliar a correlação entre os parâmetros de infiltração de água sortividade e expoente de tempo (n) e os parâmetros dimensão fractal (D) e expoente de Hurst (H). Para isso, dez colunas horizontais foram simuladas em computador, sendo preenchidas com material de textura franca ou argilosa, puros ou em combinações de camadas alternadas dos dois materiais, seguindo a distribuição de um Conjunto de Cantor determinístico com dimensão fractal 0,63. As simulações de movimento horizontal de água foram realizadas no modelo numérico Hydrus 2D. Os parâmetros sortividade (S) e expoente de tempo (n) da equação de Philip foram estimados para cada simulação, utilizando o módulo de regressão não linear do software SAS®. A sortividade aumentou com o aumento de material de textura franca nas colunas, o que foi atribuído à dependência do parâmetro S no raio capilar. O parâmetro expoente de tempo estimado por regressão não linear mostrou valor menor que o valor empírico de 0,5 sugerido na literatura. A dimensão fractal estimada utilizando o expoente de Hurst foi 17,5 % menor do que a dimensão fractal do Conjunto de Cantor utilizado para gerar as colunas.Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns

Fractal Characterization of Fracture Networks: An Improved Box-counting Technique

Box counting is widely used for characterizing fracture networks as fractals and estimating their fractal dimensions (D). If this analysis yields a power law distribution given by N $\propto$ r−D, where N is the number of boxes containing one or more fractures and r is the box size, then the network is considered to be fractal. However, researchers are divided in their opinion about which is the best box‐counting algorithm to use, or whether fracture networks are indeed fractals. A synthetic fractal fracture network with a known D value was used to develop a new algorithm for the box‐counting method that returns improved estimates of D. The method is based on identifying the lower limit of fractal behavior (rcutoff) using the condition ds/dr → 0, where s is the standard deviation from a linear regression equation fitted to log(N) versus log(r) with data for r \u3c rcutoff sequentially excluded. A set of 7 nested fracture maps from the Hornelen Basin, Norway was used to test the improved method and demonstrate its accuracy for natural patterns. We also reanalyzed a suite of 17 fracture trace maps that had previously been evaluated for their fractal nature. The improved estimates of D for these maps ranged from 1.56 ± 0.02 to 1.79 ± 0.02, and were much greater than the original estimates. These higher D values imply a greater degree of fracture connectivity and thus increased propensity for fracture flow and the transport of miscible or immiscible chemicals

Impact of the capillary pressure-saturation pore-size distribution parameter on geological carbon sequestration estimates

Cost estimates for geologic carbon sequestration (GCS) are vital for policy and decision makers evaluating carbon capture and storage strategies. Numerical models are often used in feasibility studies for the different stages of carbon injection and redistribution. Knowledge of the capillary pressure-saturation function for a selected storage rock unit is essential in applications used for simulating multiphase fluid flow and transport. However, the parameters describing these functions (e.g. the van Genuchten m pore size distribution parameter) are often not measured or neglected compared to other physical properties such as porosity and intrinsic permeability. In addition, the use of average instead of point estimates of m for numerical simulations of flow and transport can result in significant errors, especially in the case of coarse-grained sediments and fractured rocks. Such erroneous predictions can pose great risks and challenges to decision-making. We present a comparison of numerical simulation results based on average and point estimates of the van Genuchten m parameter for different porous media. Forward numerical simulations using the STOMP code were employed to illustrate the magnitudes of the differences in carbon sequestration predictions resulting from the use of height-averaged instead of point parameters. The model predictions were converted into cost estimates and the results indicate that varying m values in GCS modeling can cause cost differences of up to hundreds of millions dollars

Solute and Bacterial Transport through Partially-Saturated Intact Soil Blocks

Steady-state transport of water, chloride and bacteria was measured through intact blocks of Maury and Cecil soils, under partially saturated conditions. Major objectives were to determine if transport occurs uniformly or via preferential flow paths, and if soil physical properties could be used to predict breakthrough. The blocks were instrumented with TDR probes and mounted on a vacuum chamber containing 100 cells that collected eflluent. After each experiment the blocks were sampled for soil physical properties. The fluxes showed no spatial autocorrelation and the eflluent variance was not statistically different between soils. Less than 3% of the influent bacteria appeared in the effluent. Maximum bacterial breakthrough occurred after 0.25 water-filled pore volumes had been leached, and was greater for Cecil soil than for Maury soil. The chloride breakthrough curves were fitted to the convection dispersion equation. The best predictor of dispersivity was volumetric water content (R2 = 0.28, P \u3c 0.01), with dispersivity increasing with decreasing water content. Lower water contents lead to more tortuous flow paths and thus, a broadening of the velocity distribution. Soil structural controls on solute dispersion under partially saturated conditions are likely to be indirect, and related to differences in water content at given flux produced by differences in pore-size distribution

Quantification of Hydrological, Geochemical, and Mineralogical Processes Governing the Fate and Transport of Uranium over Multiple Scales in Hanford Sediments

A long-term measure of the DOE Environmental Remediation Sciences Division is to provide sufficient scientific understanding to allow a significant fraction of DOE sites to incorporate coupled biological, chemical, and physical processes into decision making for environmental remediation and long-term stewardship by 2015. Our research targets two related, major obstacles to understanding and predicting contaminant transport at DOE sites: the heterogeneity of subsurface geologic media, and the scale dependence of experimental and modeled results

Investigating predictors of superior face recognition ability in police super-recognisers

There are large individual differences in the ability to recognise faces. Super-recognisers are exceptionally good at face memory tasks. In London, a small specialist pool of police officers (also labelled ‘super-recognisers’ by the Metropolitan Police Service) annually makes 1000’s of suspect identifications from closed-circuit television footage. Some suspects are disguised, have not been encountered recently or are depicted in poor quality images. Across tests measuring familiar face recognition, unfamiliar face memory and unfamiliar face matching, the accuracy of members of this specialist police pool was approximately equal to a group of non-police super-recognisers. Both groups were more accurate than matched control members of the public. No reliable relationships were found between the face processing tests and object recognition. Within each group, however, there were large performance variations across tests, and this research has implications for the deployment of police worldwide in operations requiring officers with superior face processing ability

Implications of Fecal Bacteria Input from Latrine-Polluted Ponds for Wells in Sandy Aquifers

Ponds receiving latrine effluents may serve as sources of fecal contamination to shallow aquifers tapped by millions of tube-wells in Bangladesh. To test this hypothesis, transects of monitoring wells radiating away from four ponds were installed in a shallow sandy aquifer underlying a densely populated village and monitored for 14 months. Two of the ponds extended to medium sand. Another pond was sited within silty sand and the last in silt. The fecal indicator bacterium E. coli was rarely detected along the transects during the dry season and was only detected near the ponds extending to medium sand up to 7 m away during the monsoon. A log-linear decline in E. coli and Bacteroidales concentrations with distance along the transects in the early monsoon indicates that ponds excavated in medium sand were the likely source of contamination. Spatial removal rates ranged from 0.5-1.3 log10/m. After the ponds were artificially filled with groundwater to simulate the impact of a rain storm, E. coli levels increased near a pond recently excavated in medium sand, but no others. These observations show that adjacent sediment grain-size and how recently a pond was excavated influence how much fecal contamination ponds receiving latrine effluents contribute to neighboring groundwater

Modeling water movement in horizontal columns using fractal theory Modelagem de movimento horizontal de água no solo utilizando a teoria dos fractais

Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns.A matemática fractal tem sido utilizada para caracterizar o transporte de água e solutos em meios porosos e também para simular características físicas e geométricas de meios porosos. O objetivo deste trabalho foi avaliar a correlação entre os parâmetros de infiltração de água sortividade e expoente de tempo (n) e os parâmetros dimensão fractal (D) e expoente de Hurst (H). Para isso, dez colunas horizontais foram simuladas em computador, sendo preenchidas com material de textura franca ou argilosa, puros ou em combinações de camadas alternadas dos dois materiais, seguindo a distribuição de um Conjunto de Cantor determinístico com dimensão fractal 0,63. As simulações de movimento horizontal de água foram realizadas no modelo numérico Hydrus 2D. Os parâmetros sortividade (S) e expoente de tempo (n) da equação de Philip foram estimados para cada simulação, utilizando o módulo de regressão não linear do software SAS®. A sortividade aumentou com o aumento de material de textura franca nas colunas, o que foi atribuído à dependência do parâmetro S no raio capilar. O parâmetro expoente de tempo estimado por regressão não linear mostrou valor menor que o valor empírico de 0,5 sugerido na literatura. A dimensão fractal estimada utilizando o expoente de Hurst foi 17,5 % menor do que a dimensão fractal do Conjunto de Cantor utilizado para gerar as colunas

Impact of the capillary pressure-saturation pore-size distribution parameter on geological carbon sequestration estimates

Cost estimates for geologic carbon sequestration (GCS) are vital for policy and decision makers evaluating carbon capture and storage strategies. Numerical models are often used in feasibility studies for the different stages of carbon injection and redistribution. Knowledge of the capillary pressure-saturation function for a selected storage rock unit is essential in applications used for simulating multiphase fluid flow and transport. However, the parameters describing these functions (e.g. the van Genuchten m pore size distribution parameter) are often not measured or neglected compared to other physical properties such as porosity and intrinsic permeability. In addition, the use of average instead of point estimates of m for numerical simulations of flow and transport can result in significant errors, especially in the case of coarse-grained sediments and fractured rocks. Such erroneous predictions can pose great risks and challenges to decision-making. We present a comparison of numerical simulation results based on average and point estimates of the van Genuchten m parameter for different porous media. Forward numerical simulations using the STOMP code were employed to illustrate the magnitudes of the differences in carbon sequestration predictions resulting from the use of height-averaged instead of point parameters. The model predictions were converted into cost estimates and the results indicate that varying m values in GCS modeling can cause cost differences of up to hundreds of millions dollars