Simulations of solute transport in fractured porous media using 2D percolation networks with uncorrelated hydraulic conductivity fields

2004 ◽  
Vol 12 (6) ◽  
pp. 613-627 ◽  
Author(s):  
Christine Rivard ◽  
Fr�d�rick Delay
2007 ◽  
Vol 7 (1) ◽  
pp. 59-66 ◽  
Author(s):  
D.E. Rolston

The science of soil-water physics and contaminant transport in porous media began a little more than a century ago. The first equation to quantify the flow of water is attributed to Darcy. The next major development for unsaturated media was made by Buckingham in 1907. Buckingham quantified the energy state of soil water based on the thermodynamic potential energy. Buckingham then introduced the concept of unsaturated hydraulic conductivity, a function of water content. The water flux as the product of the unsaturated hydraulic conductivity and the total potential gradient has become the accepted Buckingham-Darcy law. Two decades later, Richards applied the continuity equation to Buckingham's equation and obtained a general partial differential equation describing water flow in unsaturated soils. For combined water and solute transport, it had been recognized since the latter half of the 19th century that salts and water do not move uniformly. It wasn't until the middle of the 20th century that scientists began to understand the complex processes of diffusion, dispersion, and convection and to develop mathematical formulations for solute transport. Knowledge on water flow and solute transport processes has expanded greatly since the early part of the 20th century to the present.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-13 ◽  
Author(s):  
Chi Yao ◽  
Chen He ◽  
Jianhua Yang ◽  
Qinghui Jiang ◽  
Jinsong Huang ◽  
...  

An original 3D numerical approach for fluid flow in fractured porous media is proposed. The whole research domain is discretized by the Delaunay tetrahedron based on the concept of node saturation. Tetrahedral blocks are impermeable, and fluid only flows through the interconnected interfaces between blocks. Fractures and the porous matrix are replaced by the triangular interface network, which is the so-called equivalent matrix-fracture network (EMFN). In this way, the three-dimensional seepage problem becomes a two-dimensional problem. The finite element method is used to solve the steady-state flow problem. The big finding is that the ratio of the macroconductivity of the whole interface network to the local conductivity of an interface is linearly related to the cubic root of the number of nodes used for mesh generation. A formula is presented to describe this relationship. With this formula, we can make sure that the EMFN produces the same macroscopic hydraulic conductivity as the intact rock. The approach is applied in a series of numerical tests to demonstrate its efficiency. Effects of the hydraulic aperture of fracture and connectivity of the fracture network on the effective hydraulic conductivity of fractured rock masses are systematically investigated.


Sign in / Sign up

Export Citation Format

Share Document