Vermiculite bio-barriers for Cu and Zn remediation: an eco-friendly approach for freshwater and sediments protection

The increase in heavy metal contamination in freshwater systems causes serious environmental problems in most industrialized countries, and the effort to find eco-friendly techniques for reducing water and sediment contamination is fundamental for environmental protection. Permeable barriers made of natural clays can be used as low-cost and eco-friendly materials for adsorbing heavy metals from water solution and thus reducing the sediment contamination. This study discusses the application of permeable barriers made of vermiculite clay for heavy metals remediation at the interface between water and sediments and investigates the possibility to increase their efficiency by loading the vermiculite surface with a microbial biofilm of Pseudomonas putida, which is well known to be a heavy metal accumulator. Some batch assays were performed to verify the uptake capacity of two systems and their adsorption kinetics, and the results indicated that the vermiculite bio-barrier system had a higher removal capacity than the vermiculite barrier (+34.4 and 22.8 % for Cu and Zn, respectively). Moreover, the presence of P. putida biofilm strongly contributed to fasten the kinetics of metals adsorption onto vermiculite sheets. In open-system conditions, the presence of a vermiculite barrier at the interface between water and sediment could reduce the sediment contamination up to 20 and 23 % for Cu and Zn, respectively, highlighting the efficiency of these eco-friendly materials for environmental applications. Nevertheless, the contribution of microbial biofilm in open-system setup should be optimized, and some important considerations about biofilm attachment in a continuous-flow system have been discussed.     

A novel methodological approach for land subsidence prediction through data assimilation techniques

Computational Geosciences - Anthropogenic land subsidence can be evaluated and predicted by numerical models, which are often built over deterministic analyses. However, uncertainties and...     

Finite element analysis of land subsidence above depleted reservoirs with pore pressure gradient and total stress formulations

The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd.     

Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

Ensemble smoothing of land subsidence measurements for reservoir geomechanical characterization

Geomechanical models are often used to predict the impact on land surface of fluid withdrawal from deep reservoirs, as well as investigating measures for mitigation. The ability to accurately simulate surface displacements, however, is often impaired by limited information on the geomechanical parameters characterizing the geological formations of interest. In this study, we employ an ensemble smoother, a data assimilation algorithm, to provide improved estimates of reservoir parameters through assimilation of measurements of both horizontal and vertical surface displacement into geomechanical model results. The method leverages the demonstrated potential of remote sensing techniques developed in the last decade to provide accurate displacement data for large areas of the land surface. For evaluation purposes, the methodology is applied to the case of a disk‐shaped reservoir embedded in a homogeneous, isotropic, and linearly elastic half space, subject to a uniform change in fluid pressure. Multiple sources of uncertainty are investigated, including the radius, R, the thickness, h, and the depth, c, of the reservoir; the pore pressure change, Δp; porous medium's vertical uniaxial compressibility, c_M, and Poisson's ratio, ν, and the ratio, s, between the compressibilities of the medium during loading and unloading cycles. Results from all simulations show that the ensemble smoother has the capability to effectively reduce the uncertainty associated with those parameters to which the variability and the spatial distribution of land surface displacements are most sensitive, namely, R, c, c_M, and s. These analyses demonstrate that the estimation of these parameters values depends on the number of measurements assimilated and the error assigned to the measurement values. Copyright © 2014 John Wiley & Sons, Ltd.     

The role of preconditioning in the solution to FE coupled consolidation equations by Krylov subspace methods

The repeated solution in time of the linear system arising from the finite element integration of coupled consolidation equations is a major computational effort. This system can be written in either a symmetric or an unsymmetric form, thus calling for the implementation of different preconditioners and Krylov subspace solvers. The present paper aims at investigating when either a symmetric or an unsymmetric approach should be better used. The results from a number of representative numerical experiments indicate that a major role in selecting either form is played by the preconditioner rather than by the Krylov subspace method itself. Two other important issues addressed are the size of the time integration step and the possible lumping of the flow capacity matrix. It appears that ad hoc block constrained preconditioners provide the most robust algorithm independently of the time step size, lumping, and symmetry. Copyright © 2008 John Wiley & Sons, Ltd.     

Direct,partitioned and projected solution to finite element consolidation models

Direct, partitioned, and projected (conjugate gradient‐like) solution approaches are compared on unsymmetric indefinite systems arising from the finite element integration of coupled consolidation equations. The direct method is used in its most recent and computationally efficient implementations of the Harwell Software Library. The partitioned approach designed for coupled problems is especially attractive as it addresses two separate positive definite problems of a smaller size that can be solved by symmetric conjugate gradients. However, it may stagnate and when converging it does not prove competitive with a global projection method such as Bi‐CGSTAB, which may take full advantage of its flexibility in working on scaled and reordered equations, and thus may greatly improve its computational performance in terms of both robustness and convergence rate. The Bi‐CGSTAB superiority to the other approaches is discussed and demonstrated with a few representative examples in two‐dimensional (2‐D) and three‐dimensional (3‐D) coupled consolidation problems. Copyright © 2002 John Wiley & Sons, Ltd.     

Parallel inexact constraint preconditioning for ill-conditioned consolidation problems

Constraint preconditioners have proved very efficient for the solution of ill-conditioned finite element (FE) coupled consolidation problems in a sequential computing environment. Their implementation on parallel computers, however, is not straightforward because of their inherent sequentiality. The present paper describes a novel parallel inexact constraint preconditioner (ParICP) for the efficient solution of linear algebraic systems arising from the FE discretization of the coupled poro-elasticity equations. The ParICP implementation is based on the use of the block factorized sparse approximate inverse incomplete Cholesky preconditioner, which is a very recent and effective development for the parallel preconditioning of symmetric positive definite matrices. The ParICP performance is experimented with in real 3D coupled consolidation problems, proving a scalable and efficient implementation of the constraint preconditioning for high-performance computing. ParICP appears to be a very robust algorithm for solving ill-conditioned large-size coupled models in a parallel computing environment.     

Numerical performance of projection methods in finite element consolidation models

Projection, or conjugate gradient like, methods are becoming increasingly popular for the efficient solution of large sparse sets of unsymmetric indefinite equations arising from the numerical integration of (initial) boundary value problems. One such problem is soil consolidation coupling a flow and a structural model, typically solved by finite elements (FE) in space and a marching scheme in time (e.g. the Crank–Nicolson scheme). The attraction of a projection method stems from a number of factors, including the ease of implementation, the requirement of limited core memory and the low computational cost if a cheap and effective matrix preconditioner is available. In the present paper, biconjugate gradient stabilized (Bi‐ CGSTAB) is used to solve FE consolidation equations in 2‐D and 3‐D settings with variable time integration steps. Three different nodal orderings are selected along with the preconditioner ILUT based on incomplete triangular factorization and variable fill‐in. The overall cost of the solver is made up of the preconditioning cost plus the cost to converge which is in turn related to the number of iterations and the elementary operations required by each iteration. The results show that nodal ordering affects the perfor mance of Bi‐CGSTAB. For normally conditioned consolidation problems Bi‐CGSTAB with the best ILUT preconditioner may converge in a number of iterations up to two order of magnitude smaller than the size of the FE model and proves an accurate, cost‐effective and robust alternative to direct methods. Copyright © 2001 John Wiley & Sons, Ltd.     

Correction to: A novel methodological approach for land subsidence prediction through data assimilation techniques

Computational Geosciences - A Correction to this paper has been published: https://doi.org/10.1007/s10596-021-10079-6