Visualisation of fluid flow mechanisms through a viscous-porous rock-analogue medium – experiment and model results

Author(s):  
Reinier van Noort ◽  
Lawrence Hongliang Wang ◽  
Viktoriya Yarushina

<p>Understanding fluid flow patterns in the shallow and deep earth is one of the major challenges of modern earth sciences. Fluid flow may be slow and pervasive, or fast and focused. In the deep earth, focused fluid flow may result in, for example, dikes, veins, volcanic diatremes and gas venting systems. In the shallow Earth, focused fluid flow can be found in the form of fluid escape pipes and gas conducting chimneys, mud volcanoes, sand injectites, pockmarks, hydrothermal vent complexes, etc.</p><p>Focused fluid flow has been reproduced in visco-plastic models of flow through porous materials. However, the mechanisms that cause fluid flow to focus along such relatively narrow channels, with transiently elevated permeability, have not been investigated thoroughly in experiments. We have carried out experiments in a transparent Hele-Shaw cell. In our experiments, a hydrous fluid is injected into an aggregate of viscous grains, and the mechanisms by which this injected fluid flows are recorded using a digital camera. Our experiments demonstrate a dependence of fluid flow mechanisms on the injection rate. At low injection rate, we observe the formation of a slowly-rising diapir. As this diapir slowly rises through the porous medium, it is fed by transient, focused fluid flow following the path of the rising diapir. Once the diapir escapes through the surface of our aggregate, continued fluid flow through the porous aggregate is focused and transient. At high injection rate, instead of a diapir fed by focused fluid flow, an open channel forms as a result of local fluidization of the granular material.</p><p>Our experimental observations are interpreted through visco-plastic models simulating the experimental conditions. These numerical models can reproduce the diapirs observed in our experiments at low flow rate by assuming flow through a layered porous aggregate, with a layer with relatively high bulk viscosity overlying a layer with relatively low bulk viscosity. For low injection rates, such a model reproduces focused fluid flow in the low-viscosity layer, that feeds into a slowly rising diapir in the high-viscosity layer. This model observation thus suggests that the passage of the rising diapir in our experiments leaves a trail, where the aggregate bulk viscosity is lowered and along which ongoing fluid flow can focus transiently.</p>

2020 ◽  
Author(s):  
Philipp Eichheimer ◽  
Marcel Thielmann ◽  
Wakana Fujita ◽  
Gregor J. Golabek ◽  
Michihiko Nakamura ◽  
...  

Abstract. Fluid flow on different scales is of interest for several Earth science disciplines like petrophysics, hydrogeology and volcanology. To parameterize fluid flow in large-scale numerical simulations (e.g. groundwater and volcanic systems), flow properties on the microscale need to be considered. For this purpose experimental and numerical investigations of flow through porous media over a wide range of porosities are necessary. In the present study we sinter glass bead media with various porosities. The microstructure, namely effective porosity and effective specific surface, is investigated using image processing. We determine flow properties like hydraulic tortuosity and permeability using both experimental measurements and numerical simulations. By fitting microstructural and flow properties to porosity, we obtain a modified Kozeny-Carman equation for isotropic low-porosity media, that can be used to simulate permeability in large-scale numerical models. To verify the modified Kozeny-Carman equation we compare it to the computed and measured permeability values.


2020 ◽  
Author(s):  
Hongliang Wang ◽  
Viktoriya Yarushina ◽  
Yury Podladchikov

<p>Fluid flow instability in deforming porous rock, commonly known as porosity waves, has been used to explain formation of seismic chimneys, one of the most important expressions for the localized fluid flow in the subsurface. Experiments show that volumetric deformation of rocks is strongly coupled with shear deformation, leading to shear-induced decompaction at low confining pressure and shear-enhanced compaction at higher confining pressure. Previous studies introduce a weakening factor of R for bulk viscosity in the viscous deforming regime. While it has successfully reproduced the channelized fluid flow in numerical models, it cannot investigate the effect of shear deformation. More controversially, negative effective pressure (P<sub>t</sub>-P<sub>f</sub>) is required for the channel formation. Here, we develop a viscoplastic rheology that takes into account effects of shear stress and plastic failure on the volumetric deformation, consistent with experimental data. A dilation pressure is naturally introduced through viscoplastic strain-rate when plastic failure occurs under high fluid pressure and shear stress condition. Our model results show that this new rheology can produce channelized fluid flow without negative effective pressure in the model.</p><p>In order to apply our models into real geological setting, we test the effects of reservoir properties, geological layering, transport properties of the layers and faults.  Our results show that fluid channel initiates at local topography highs in the reservoir and a high-permeability fault can also trigger the initiation of fluid channels. Fluid channels can have different length and time scales in different layers, depending on bulk viscosity and permeability of the layers.</p>


Author(s):  
Jizu Lv ◽  
Zhenxian Zhang ◽  
Chengzhi Hu ◽  
Minli Bai

Abstract In addition to the increase of thermal conductivity, heat transferring for nanofluids strengthening mechanism also includes the changes of the flow characteristics, therefore it is needed to take an in-depth research on nanofluids flow characteristics. However previous visualization experiment do not quantitatively analyze the change of flow characteristics after nano-particles is added, do not reveal the mechanism of nanofluids changing the characteristics of the fluid in intense turbulent flow condition. Therefore in this paper, by means of particle image velocimetry, quantitatively study SiO2-water nanofluids flow characteristics in intense turbulent flow condition and analyze the influence of SiO2-water nanofluids on turbulent flow energy by measuring the pressure drop caused by fluid flowing through the channel. Fluid flow through rectangular convex channel (channel composed of continuous staggered rectangular convex platform) to obtain the steady intense turbulent flow. The rectangular convex channel makes the fluid flow through obstacles for several times, so that the flow direction changes for several times, and vortexes are generated in the local scope, which makes turbulence enhance and increases minor loss. In this way, flow can be in the intense turbulent state under a low flow rate, which meets the experiment requirement and is convenient to compare the influence of nano-particles on flow resistance and energy loss. The experiment takes the quantitative PIV experimental research on pure water and the volume fraction of 0.5% SiO2-water nanofluids respectively in the Reynolds number is 2300, 2500, 3000, 4000. Through the experiment, we can obtain nanofluids turbulent flow condition fluctuating velocity, turbulence kinetic, energy loss and so on, and the fluid flow velocity vector, streamline and vorticity graph. Through the quantitative comparison of the spiral numbers, the vorticity distribution, and energy loss, analyze strengthening effect and influence of flow resistance on basic fluid after adding nanoparticles.


2015 ◽  
Author(s):  
◽  
Mijo Nikolić

This thesis deals with the problem of localized failure in rocks, which occurs often in civil engineering practice like in dam failure, foundation collapse, stability of excavations, slopes and tunnels, landslides and rock falls. The risk of localized failure should be better understood in order to be prevented. The localized failure in rocks is usually characterized by a sudden and brittle failure without warning in a sense of larger and visible deformations prior to failure. This happens also under the strong influence of material heterogeneities, pre-existing cracks and other defects. The three novel numerical models, incorporating the localized failure mechanisms, heterogeneity of rock and pre-existing cracks and other defects, are presented in this thesis. First model deals with 2D plane strain two-phase rock composite, where stronger phase represents the intact rock and weaker phase initial defects. Second model represents the extension of the previous model towards the 3D space, where full set of 3D failure mechanisms is considered. Heterogeneous properties are taken here through the random distribution and Gauss statistical variation of material properties. The latter model is also used for the analysis of intact rock core specimens geometrical shape deviations influencing the uniaxial compressive strength. Third model is a 2D, dealing with volumetric fluid-structure interaction and localized failure under the influence of fluid flow through the porous rock medium. The discrete beam lattice approach is chosen for general framework for three models, where Voronoi cells represent the rock grains kept together by Timoshenko beams as cohesive links. The enhanced kinematics characterized for embedded discontinuity approach is added upon standard kinematics of cohesive links. This serves for the macrocrack propagation in all failure modes and their combinations, between the rock grains. The fracture process zone formation followed by micro-cracks coalescence, preceding the localized failure, is considered as well. Fluid flow is governed by a Darcy law, while coupling conditions obey Biot’s theory of poroplasticity. The results of the numerical models were verified by the benchmarks available from literature in 2D case. The 3D model was validated against the experimental results conducted on 90 rock specimens, where even slight geometrical deviations of specimens are considered. Presentation of these models, as well as their implementation aspects are given in full detail. Embedded discontinuity concept and the local nature of enhancements required to capture the displacement discontinuities leads to the very efficient approach with static condensation of enhanced degrees of freedom and technique that can be efficiently incorporated into finite element code architecture.


2020 ◽  
Author(s):  
Reinier van Noort ◽  
Viktoriya Yarushina

<p>Seismic chimneys have been observed in sediments overlying reservoirs containing different fluids, such as water, hydrocarbons, or CO2. Furthermore, such chimneys have been linked to pockmarks and gas seepages on the seafloor. Visco-plastic models show how these chimneys can form by focused fluid flow through viscous, porous materials. However, the mechanisms that cause fluid flow to focus along such relatively narrow pathways with transiently elevated permeability have not been investigated thoroughly in experiments.</p><p>We present analogue experiments carried out in a transparent Hele-Shaw cell, in which a fluid is injected into an aggregate of viscous grains, leading to transient focused fluid flow. Fluid flow is imaged using a digital camera, and our observations are compared to models describing chimney formation.</p>


2020 ◽  
Author(s):  
Philipp Eichheimer ◽  
Marcel Thielmann ◽  
Wakana Fujita ◽  
Gregor J. Golabek ◽  
Michihiko Nakamura ◽  
...  

<div> <div> <div> <p>Fluid flow on different scales is of interest for several Earth science disciplines like petrophysics, hydrogeology and volcanology. To parameterize fluid flow in large-scale numerical simulations (e.g. groundwater and volcanic systems), flow properties on the microscale need to be considered. For this purpose experimental and numerical investigations of flow through porous media over a wide range of porosities are necessary. In the present study we sinter glass bead media with various porosities, representing shallow depth crustal sediments. The microstructure, namely effective porosity and effective specific surface, is investigated using image processing. We furthermore determine flow properties like hydraulic tortuosity and permeability using both experimental measurements and numerical simulations. By fitting microstructural and flow properties to porosity, we obtain a modified Kozeny-Carman equation for isotropic low-porosity media, that can be used to simulate permeability in large-scale numerical models. To verify the modified Kozeny-Carman equation we compare it to the numerically computed and experimentally measured permeability values.</p> </div> </div> </div>


2016 ◽  
Vol 19 (5) ◽  
pp. 391-404 ◽  
Author(s):  
B. M. Shankar ◽  
I. S. Shivakumara ◽  
Chiu-On Ng

Sign in / Sign up

Export Citation Format

Share Document