Packed Glass Beads Research Articles

Study of Hydrodynamics and Residence Time Distribution in a Trickle Bed Reactor

This study focuses on investigating the hydrodynamics using computational fluid dynamics (CFD) simulation and residence time distribution (RTD) within a trickle bed reactor (TBR). CFD simulation was performed on a packed column filled with structured packing of glass beads, while the experiments for RTD studies were conducted in a TBR filled with random packing of Raschig rings. The geometry of TBR and size of packing during CFD simulation was kept close to that of the experimental setup. CFD simulation results were obtained which predicted the flow of fluid through the TBR at different vertical planes along the column. The fluctuations of liquid velocity at various positions inside the reactor were predicted by the CFD model. RTD studies were performed for both pulse input and step input using a non-reactive tracer such as sodium hydroxide. The flowrate of water was varied from 0.6 LPH to 6 LPH for pulse input, while for step input the ratio of water to tracer flowrates was kept at 1:1, 1:2 and 2:1. This study allows us to understand how hydrodynamics and residence time distribution may have an effect on the performance of TBR.

Systematic Study of Wettability Alteration of Glass Surfaces by Dichlorooctamethyltetrasiloxane Silanization-A Guide for Contact Angle Modification.

To investigate the effects of wettability on multiphase flow in porous media, glass bead packs or micromodels are commonly used. Their wettability can be altered by the surface treatment method-silanization. Although silanization is widely used for glass wettability modification, comparable systematic approaches over a large range of geometries, treatment conditions, and measurement systems are scarce. In this work, dichlorooctamethyltetrasiloxane (Surfasil) treatment was systematically investigated, resulting in a guide for achieving a wide range of contact angles. Initially, the influence of the Surfasil solvent, treatment time, and Surfasil-to-solvent ratio was investigated on glass plates using the sessile drop method. By varying these variables, it was possible to achieve a wide range of comparable, repeatable, and stable contact angles, from approximately 20-95° for air-water systems. Due to the linear increase of contact angle with larger Surfasil exposure, either due to the time or concentration, contact angle tuning is possible until the critical point. Beyond the critical point of exposure, a system-specific plateau value is reached, independent of the approach. After establishing a clear relationship between the parameters and contact angles, the same treatment parameters were applied to single beads, micromodels, and beadpacks with heptane as the chosen solvent. Optical image analysis was used for the microchips, and micro CT data analysis was used for the bead packs. The treatment appeared to be transferable to all geometries, resulting in similar wetting conditions within the limitations of the measurements. It is concluded that a glass plate can be used as an analogue for obtaining the contact angle alteration trends for more complex porous media with similar compositions. Data analysis methods and surface roughness could have an effect on the obtained contact angle spread.

Prospective use of a cell‐free supernatant from the biosurfactant‐producing marine Planococcus maritimus SAMP MCC 3013 for oil removal and enhanced oil recovery

AbstractThe cell‐free supernatant (CFS) of the biosurfactant (BS) producing pure strain of Planococcus maritimus SAMP MCC 3013 was used for enhanced oil recovery (EOR). An indigenous marine culture of P. maritimus was grown in mineral salt medium (MSM) supplemented with glucose (1%) (w/v) (30°C, 120 rpm, 7 days). CFS was used to demonstrate its application for oil recovery through three assays: pumice stone, saturated gravel, and the glass beads‐packed column. The recovery of light paraffin oil was ∼40% in the pumice stone assay and ∼43% in the oil‐saturated gravel assay, respectively. The glass bead column demonstrated ∼42% of additional 2T engine oil recovery by CFS. Our study showed that BS containing CFS from P. maritimus SAMP MCC 3013 has the ability to recover oil in its crude form without using any sophisticated purification techniques. All three methods operated by us are harmonizing with each other with respect to the findings. Those methods can be employed in large scale operations for the bio‐recovery of oils using an environmentally friendly approach through the functional properties of BS. This is most likely the first report describing the use of a glass bead pack column for EOR applications particularly from Planococcus sp. These results are symbolic of prospective applications of BS containing CFS from indigenous P. maritimus SAMP MCC 3013 under ex‐situ conditions.

Pore size distribution measurement with magnetic resonance [formula omitted] distributions outside the fast diffusion regime

Based on Brownstein-Tarr (BT) theory, a new method is proposed to determine the pore size distribution which is an important petrophysical property in studying fluid storage and fluid transport in reservoir rocks. For reservoir rocks with large pores, or high surface relaxivity, direct conversion from the magnetic resonance (MR) T2 distribution with the fast diffusion assumption may be unreliable. In this study, the ground mode lifetime T20 distribution was proposed to be utilized to estimate the pore size distribution. The T20 lifetime of large pores within a porous system can be determined directly through experiment. However, the T20 lifetimes of smaller pores may overlap with the first nonground mode T21 lifetimes of the larger pores, which collectively contribute to the MR signal with a shorter lifetime. To accurately determine the T20 lifetime of smaller pores, it is essential to remove the contribution to the distribution of the T21 lifetimes of the larger pores. The T21 distribution can be both experimentally observed and calculated through the proposed method. The proposed pore size distribution method is verified by calculation in the present work. Water-saturated glass bead packs were employed for initial experiments. The calculated pore size, and its distribution, match the geometric pore size. The proposed method was applied to determine the pore size distribution of a Berea sandstone. The results are in good agreement with SEM measurements.

Stick‐Slip Nucleation and Failure in Uniform Glass Beads Detected by Acoustic Emissions in Ring‐Shear Experiments: Implications for Identifying the Acoustic Emissions of Earthquake Foreshocks

AbstractStress accumulation and release reflected by acoustic emissions (AEs) during shearing of granular materials provide important information on failure mechanisms in seismic faults and landslides controlled by stick‐slip. Among many characteristics (amplitude, energy, counts, and frequency) of AE signals generated by stick‐slip, stress changes corresponding to various frequency AEs in different stages of the stick‐slip process are not clear, which limits our knowledge of the characteristics of precursory signals before stick‐slip failure. To better understand the physical mechanisms of granular stick‐slip, we monitored the mechanical and AE signals using high‐frequency (2 MHz) synchronous acquisition during constant‐speed shear of packs of uniform glass beads with different sizes at different normal stresses. The release rate of AE energy was found to accelerate with the dilatation of the sample volume, and the stress drop of stick‐slip was augmented with the increase of normal stress and particle size. Three characteristic events of single cycle stick‐slip were observed in this study: main slip, minor slip, and microslip. We analyzed the AE frequency spectra of these three event types. Both main slip and minor slip corresponded to stress drop and generated high‐frequency AEs (about several hundred kHz), while the AE frequencies generated by microslip were lower (about tens of kHz) and exhibited stress strengthening, which were not apparent in previous studies due to the low frequency of acquisition. We propose that the microslip is mainly due to sliding on grain contacts, while the main slip and minor slip resulted from breakage and reforming of force chains. Low‐frequency AEs from microslip may suggest a crucial precursor of seismic faults and landslides.

Stress sensitivity of elastic moduli in high-porosity cemented sandstone — Heuristic models and experimental data

Quantitative description of the reservoir rock stress sensitivity is critical for seismic modeling and interpretation. Laboratory studies have indicated that the wave velocity of high porosity cemented sandstone is asymmetrically much more sensitive to stress unloading than stress loading. However, there is a lack of a simple rock physics model that is capable of describing the asymmetric stress sensitivity of reservoir rock. A new rock physics model is developed by extending an existing rock physics model based on contact theory combined with elastic bounds. The new model relates the stress softening of rock to cement cracking and crumbling induced by stress release. The diluting factor is introduced to describe the weakening of effective cement and analyze the increase of stress sensitivity when the stress is gradually removed. The combination of the new model and its base model forms a rock physics modeling workflow that can accurately describe the evolution of velocities measured in samples undergoing stress loading and unloading. The model performance is compared with samples representing synthetic weakly cemented glass bead pack and sandstone manufactured with different cements at different forming stress levels. The modeling results are in accordance with the measured stress sensitivities of wave velocities. The sequential model calibration using a simple constraint optimization approach yields important calibration parameters that are indicative of the elastic stress sensitivity and damage behavior in the studied rocks. The model can be particularly interesting for time-lapse monitoring of fluid injection and seismic interpretation of overpressured reservoir rock.

Simplified scaling relations for the depth-dependence and IFT reduction of fluid imbibition in gas-saturated reservoir rocks

Imbibition of water into fractured unconventional reservoirs is a significant factor contributing to the large volumes of permanently lost water by hydraulic fracturing. The lost water partitions at the fracture-matrix interface thereby restricting gas flow into the fractures. Thus, minimizing the large water footprint during the fracturing process is very attractive for water conservation and gas production. A key step towards achieving this goal is quantifying uncertainties needed to reliably predict water imbibition and to evaluate the success of efforts aimed at minimizing it. In this contribution, we present new scaling relations for water imbibition in gas-saturated rocks that account for the depth-dependence of fluid properties and chemical alteration of IFT through the use of a surfactant. We then propose a new method for quantifying uncertainties relating to rock properties, two-phase fluid properties, and the driving force for imbibition. We applied the proposed approach to analyze glass-bead pack and shale experiments involving water imbibition into air-filled pores, and validated our model predictions against published values from X-ray CT images and transient permeability measurements. Our work provides a methodology for quantifying the uncertainties pertaining to water imbibition, thereby enabling optimization of post-fracturing interventions aimed at minimizing water loss in deep unconventional shale gas reservoirs.

Modeling the Transition between Localized and Extended Deposition in Flow Networks through Packings of Glass Beads

We use a theoretical model to explore how fluid dynamics, in particular, the pressure gradient and wall shear stress in a channel, affect the deposition of particles flowing in a microfluidic network. Experiments on transport of colloidal particles in pressure-driven systems of packed beads have shown that at lower pressure drop, particles deposit locally at the inlet, while at higher pressure drop, they deposit uniformly along the direction of flow. We develop a mathematical model and use agent-based simulations to capture these essential qualitative features observed in experiments. We explore the deposition profile over a two-dimensional phase diagram defined in terms of the pressure and shear stress threshold, and show that two distinct phases exist. We explain this apparent phase transition by drawing an analogy to simple one-dimensional mass-aggregation models in which the phase transition is calculated analytically.

Investigation of ultrasound propagation through dense glass bead packings immersed in water

Conventional ultrasound imaging generally relies on a single scattering assumption and a constant sound hypothesis. However, in dense granular sediments where the glass beads diameter is comparable to the wavelength (typically 500 μm), both hypothesizes are no longer valid, resulting in a loss of resolution and contrast, up to the situation where a target to image could totally vanish in the fog. To better characterize this strongly scattering medium, we have conducted several ultrasonic experiments in reflection and in transmission, either on the diffused wave or on the coherent one. Both the transport and scattering mean free paths are estimated in the imaging bandwidth (1–5 MHz), alongside with the phase velocity and extinction length of the coherent wave. These results are compared with an analytical model, i.e. Generalized Coherent Potential Approximation*, which considers a coated sphere buried in an effective medium. Finally, we show a dispersive beamformer which takes into account the measured phase and group velocity. This allows us to obtain a more resolved and contrasted image in such a scattering medium. * Jing X, Sheng P, Zhou M. Phys Rev Lett. 1991

Acoustic Localization of an Intruder in a Strongly Scattering Medium

Localizing an intruder submerged in a strongly scattering medium, such as a dense granular suspension, is a practical challenge. Here we extract the coherent ultrasonic echo from a steel ball submerged in a dense glass-bead packing saturated by water, by using a standard single-element ultrasonic transducer thanks to a configuration averaging process. Different configurations of the granular packing are created by the nonaffine motion of the beads with a mixing blade, akin to the Brownian motion, in the vicinity of the intruder. We investigate the efficiency of this process to reduce the so-called material noise from multiply scattered ultrasound, as a function of the configuration number, the shear rate and the blade-intruder distance. Nonaffine motions of the beads in the shear zone induced by the blade are then analyzed on the basis of a split-bottom rotating shear cell. This method helps to develop not only ultrasonic imaging tools of buried objects in turbid marine sediments, but also the local rheology based on a ball falling monitored by ultrasonic tracking.

Application of critical path analysis to streaming potential coupling coefficient in porous media

Streaming potential is a passive hydrogeophysical method that can be used to monitor water flow in porous media. However, a quantitative use of this method requires a good understanding of the streaming potential generation in water saturated porous media. In this work, we used the Critical Path Analysis (CPA) method to propose a new model to estimate the streaming potential coupling coefficient (SPCC) in porous media. By comparing the proposed model with 109 experimental data points for both unconsolidated samples such as uniform glass bead and sand packs and consolidated samples such as sandstones, carbonates and with 36 simulated data points, we demonstrated that the CPA is applicable to model the SPCC for porous media with both narrow and broad pore size distributions. Additionally, from the obtained result for the uniform grain packings, we found the rc–d relationship between the critical pore radius rc and grain diameter d that was validated by using 70 experimental data points available in literature. The obtained rc–d relationship may be helpful for applications of CPA-based models for uniform grain packings. We believe that the CPA-based model is an useful alternative approach for modeling the streaming potential coupling coefficient in porous media along with two other models available in literature.

Permeability in the thin section

We compute the porosity and specific surface area of a two-dimensional (2D) pore-scale segmented image of a natural rock sample, add existing relations linking tortuosity to porosity, and then use the Kozeny-Carman equation to obtain estimates of the absolute permeability of the three-dimensional (3D) host sample. These estimates are very close to the true permeability, measured or computed, of the host sample in rocks with well-connected pore space, including high-porosity sand, glass bead packs, low-to-medium porosity sandstones, and some carbonates. They are higher than the true permeability in carbonates with relatively large pores connected by constricted conduits. They fail in impermeable synthetically constructed samples with disconnected pores. Still, this extremely simple method proves to be viable for a large class of natural rocks. Moreover, it can be used to construct a permeability-porosity trend from a microscopic 2D image by selecting subsections of this image and computing their permeability-porosity pairs. Because the measured or computed permeability-porosity pair of the host 3D sample falls upon such a trend, we argue that this trend may be applicable at a spatial scale much coarser that of the pore-scale image.

Wettability-dependent wave velocities and attenuation in granular porous media

Understanding seismic wave propagation in granular porous media is important for subsurface characterization. The presence of fluids, their distribution, and the prevailing wettability condition result in additional complexities. Although it is known that wave propagation in dry granular porous media is dominated by the presence of force chains, their influence in (partially) saturated granular porous media with different wettability conditions remains largely unexplored. To make progress in this direction, we have designed laboratory experiments by combining core flooding and ultrasonic measurements in glass bead packings that are chemically treated to alternate the wettability. The P- and S-wave velocity-saturation relation and attenuation-saturation relation are retrieved from the waveforms for water- and gas-wetting samples. The results demonstrate that there is a transition from an attenuating but stable P-wave pulse at low and moderate saturation to a set of incoherently scattered waves at high saturation. The incoherent scattering in the gas-wetting case is negligibly small, whereas it is more pronounced in the water-wetting case. We interpret these observations in terms of the wettability-dependent ability for water to penetrate into grain contacts. In the water-wetting case, liquid bridges are thought to locally reinforce the force chains and to increase their characteristic length scale. This leads to an increase in P-wave velocity and promotes incoherent scattering because the ratio of dominant wavelength to characteristic length scale decreases. In the gas-wetting case, however, the presence of gas prevents the water from direct contact with the glass beads and therefore stops the formation and growth of the liquid bridges within the force-chain network.

Experimental study of hydraulic properties in grain packs: Effects of particle shape and size distribution

Uniform and multi-dispersed grain packs have been frequently used to conceptually study flow in porous media. Numerical simulations were previously used to address the effect of particle shape on characteristics, such as pore space fractal dimension, moisture characteristic curve (MCC) and saturated hydraulic conductivity (SHC) of grain packs. However, experimental observations are still required since fractal-based approaches have been extensively proposed to model various properties in porous media. In this study, 16 angular sand and 16 spherical glass bead samples with different particle size distributions (PSDs) from well- to poorly-sorted were packed. The MCC was measured using the combination of sandbox and pressure plates methods. The pore space fractal dimension (DMCC), calculated from the measured MCC, ranged from 0.80 to 2.86 in sand and from ‐0.18 to 2.81 in glass bead packs, which indicated that DMCC may be negative in homogeneous media (e.g., glass bead packs) consistent with several studies in the literature. Results showed greater DMCC for the sand packs than the glass bead packs with the same geometric mean diameter values and PSDs. This clearly demonstrated the effect of particle shape on DMCC in the studied packs. The critical path analysis (CPA) approach was used to estimate the SHC measured using the constant-head method. We found that the CPA estimated the SHC accurately, within a factor of four of the measurements on average. Although the CPA is theoretically known to be accurate in media with broad pore size distributions, we experimentally found that it estimated the SHC in various types of grain packs reasonably well.

New Structural Percolation Transition in Fractional Wet 3D‐Porous Media: A Comparative μ‐CT Study

AbstractNatural porous media possess spatial heterogeneous (fractional) wettability, which controls the multiphase flow behavior, that is, displacement and trapping. We used mixed hydrophilic (contact angle θ = 0°) and hydrophobic (θ = 100°) 1 mm sands and 1 mm glass beads as model systems for fractional wet porous media. Both porous media have the same morphological characteristics (Minkowski functions). A comparative μ‐CT study of such fractional wet systems is lacking in the literature. A previous study (Geistlinger & Zulfiqar, 2020, https://doi.org/10.1029/2019wr026826) showed a percolation transition of the invading fluid from compact to fractal displacement when the wettability changed from imbibition to drainage. This correlates to a transition from no to maximal trapping. For the first time, we observed in fractional wet glass bead packs, a second structural percolation transition with a percolation threshold of approximately 50% hydrophobic glass beads. Below this threshold, the percolating cluster consisted of hydrophilic glass beads, which resulted in compact displacement without any trapping. Above the threshold, the percolating cluster consisted of hydrophobic glass beads, which resulted in fractal displacement with maximal trapping. Trapping occurred only in hydrophobic environments/configurations, that is, when the trapped clusters were surrounded by hydrophobic glass beads. Fractional wet glass beads and sands show an opposite trapping behavior, that is, the trapping efficiency increases for glass beads, whereas it decreases for sands with increasing percentage of hydrophobic grains. Interestingly, the trapping efficiency approaches the same limit value for pure hydrophobic sands and glass beads. This demonstrates that bypass trapping is the dominant trapping mechanism, because it depends strongly on the pore connectivity.

Effect of acoustically-induced elastic softening on sound absorption coefficient of hollow glass beads with inner closed cavities.

This study investigated the effects of sound wave magnitude on the sound absorption characteristics of granular materials. The normal incidence sound absorption coefficients of the hollow glass bead packings were measured with different incident sound pressure levels. The measurement results exhibited spectral peaks of sound absorption coefficients at several frequencies, implying that the first peak was caused by the resonance of the hollow glass beads as a layer. Although the first peak in the sound absorption coefficients did not vary when an incident sound pressure was low, the first peak shifted to a lower frequency when the sound pressure exceeded a certain magnitude. These results indicate that the hollow glass beads are softened at high incident sound pressures. Furthermore, by modelling the velocity-dependent elasticity of the hollow glass beads, the velocity and elasticity in the hollow glass beads were simulated in the direction of the sound propagation. The simulation results indicate that the velocity and elasticity profiles depend on the sound wave magnitude. Finally, the velocities of the hollow glass beads under acoustic excitation were measured by inserting an accelerometer into the beads. The measurement results demonstrated that the velocity profile depends on the sound wave magnitude, which parallels the simulation results.

Inhibited nanobubble transport in a saturated porous medium: Effects of deposited colloidal particles

Nanobubbles have recently attracted much interest for their practical use and potential applications in environmental issues. The pre-existence of deposited colloidal particles in porous media commonly occurs when nanobubbles applied to porous media interact with deposited colloidal particles. However, the current understanding of the effect of the interactions with pre-deposited colloidal particles on nanobubble transport in saturated porous media remains incomplete, and the effects are often overlooked. Therefore, we performed 1D column experiments with sequential injections of colloidal and nanobubble suspensions to study the effect of pre-deposited materials on the retention and release of colloids and nanobubbles in packed glass beads. In this study, we used resonant mass measurements to measure the number concentrations of colloids and nanobubbles during transport experiments for the first time to distinguish between coexisting solid colloidal particles and nanobubbles with different buoyancies in mixed effluent during transport. The nanobubble retention increased because of the pre-existence of deposited colloidal particles, indicating that the deposited colloidal particles act as additional deposition sites and physical obstacles for nanobubbles through physicochemical (including hydrophobic) interactions. This study also provides a future reference for the applicability of resonant mass measurement to cotransport experiments of different buoyant particles, including colloids, nanobubbles, and oil droplets.

Single-projection reconstruction technique for positioning monodisperse spheres in 3D with a divergent x-ray beam

The measurement of the position of single-sized spheres in 3D from a single, divergent, radiographic projection is addressed in the present study with the development of a novel method. Generally speaking, the location of the shadow cast by a single sphere on a detector defines a source-detector ray; the position of the particle along this ray is identified by the strong prior knowledge of its radius and the size of the shadow. For a dense assembly of equal-sized particles whose projections overlap, a novel Fourier transform based technique is introduced to give a first 3D determination of the particle centres. The uncertainty of this measurement is calculated from synthetic data with a known noise distribution. A further refinement of this measurement is performed based on the minimisation of the projection residual. The combined approach is validated both on synthetic data, and on real radiographs of a glass bead packing. The effect of noise on the measurement uncertainty is evaluated. The technique is made available to the community in the open source python package radioSphere.

Evaporation-induced salt crystallization and feedback on hydrological functions in porous media with different grain morphologies

Evaporation from porous media in the presence of dissolved salts is an important topic in relation to several environmental issues. In this study, we investigated the interaction of evaporation, salt crystallization and the internal structure of the porous material. Porous columns with different internal structures were generated through random packing of angular quartz sand and round glass beads. Evaporation-induced salt efflorescence on the surface of porous media was recorded by a camera with micro lens. The internal pore structure before and after salt crystallization was visualized using the x-ray computed tomography and quantified by the extracted pore networks. Hydraulic properties in the absence/presence of salt crystals were simulated numerically using the pore-network approach and partially measured experimentally. The evaporative water loss ratio in the angular quartz sand was greater than that in the glass beads with the same sieving size. The evaporation of saline water in angular quartz sand and round glass beads can be separated into three stages similar to water evaporation. Salt efflorescence was generated a crust dome architecture on the surface of the angular sand, but was in tight patchy contact with the round glass beads. The presence of salt crust interrupted the exchange of vapor between the porous media and atmosphere, hindering evaporation. Salt subflorescence clogged some of the pores, and the pore size distribution shifted to smaller sizes part. The porosity decreased and the geometric tortuosity of the porous materials increased owing to salt subflorescence. The presence of subflorescence salt crystals decreased the intrinsic permeability and relative diffusion coefficient of different porous media. The relationship between the vapor diffusion coefficients and porosity is linear in the presence/absence of salt crystallization, considering the geometric tortuosity. The findings indicate that the feedback of salt precipitation on hydraulic properties and vapor diffusion must be considered to accurately quantify hydrological processes, particularly in arid and semi-arid areas where intensive evaporation is accompanied by soil salinization and water shortage.

Three-Dimensional Imaging of Dissolution-Driven Convection in Simple Porous Media Using X-ray CT Scanning

Motivated by geological carbon dioxide (CO2) storage and the associated process of CO2 convective mixing in brine, we present an experimental method to visualize dissolution-driven convection in a three-dimensional porous medium. Owing to the difficulty of imaging the time-dependent convective process within an opaque medium, experiments so far have largely focused on two-dimensional systems. However, this imposes a strong restriction on the lateral spreading of the plumes. To explore the flow pattern within a three-dimensional medium, we present an experimental procedure employing X-ray computed tomography that enables the visualization of the spatial and temporal evolution of the mixing process non-invasively. To mimic the dissolution of CO2 in brine, we use potassium iodide (KI) in water. Here, we show the results with high spatial and temporal resolution for an experiment with a homogeneous packing of glass beads (Ra=5300). We discuss the evolution of the average solute mass fraction, compare it to numerical results of a purely diffusive process and report a strong enhancement of the dissolution process when convection dominates over diffusive mass transfer effects.