Investigation of the resonance of nonwetting droplets in constricted capillary tubes

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. ID1-ID17 ◽  
Author(s):  
Chao Zeng ◽  
Wen Deng ◽  
Jing Fan ◽  
Kelly H. Liu

Fluids in porous media could resonate due to passing seismic waves when meeting certain conditions of excitation frequency and pore geometry. This study of the resonance of fluids can contribute to applications in hydrocarbon microtremor analysis, enhanced oil recovery methods, and environmental remediation of nonaqueous phase liquid contamination. In this study, the analytical expression of the frequency response function and temporal response function is derived to characterize the resonance of two-phase immiscible fluids in the constricted tube. The computational fluid dynamics modeling is used to validate this resonance theory by achieving a good agreement in the prediction of output/input amplification ratio in the frequency and time domains.

SPE Journal ◽  
2013 ◽  
Vol 18 (06) ◽  
pp. 1092-1101 ◽  
Author(s):  
Amgad Salama ◽  
Shuyu Sun ◽  
M.F.. F. El-Amin

Summary The flow of two or more immiscible fluids in porous media is widespread, particularly in the oil industry. This includes secondary and tertiary oil recovery and carbon dioxide (CO2) sequestration. Accurate predictions of the development of these processes are important in estimating the benefits and consequences of the use of certain technologies. However, this accurate prediction depends—to a large extent—on two things. The first is related to our ability to correctly characterize the reservoir with all its complexities; the second depends on our ability to develop robust techniques that solve the governing equations efficiently and accurately. In this work, we introduce a new robust and efficient numerical technique for solving the conservation laws that govern the movement of two immiscible fluids in the subsurface. As an example, this work is applied to the problem of CO2 sequestration in deep saline aquifers; however, it can also be extended to incorporate more scenarios. The traditional solution algorithms to this problem are modeled after discretizing the governing laws on a generic cell and then proceed to the other cells within loops. Therefore, it is expected that calling and iterating these loops multiple times can take a significant amount of computer time. Furthermore, if this process is performed with programming languages that require repeated interpretation each time a loop is called, such as Matlab, Python, and others, much longer time is expected, particularly for larger systems. In this new algorithm, the solution is performed for all the nodes at once and not within loops. The solution methodology involves manipulating all the variables as column vectors. By use of shifting matrices, these vectors are shifted in such a way that subtracting relevant vectors produces the corresponding difference algorithm. It has been found that this technique significantly reduces the amount of central-processing-unit (CPU) time compared with a traditional technique implemented within the framework of Matlab.


2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Tatiana I. Becker ◽  
Yuriy L. Raikher ◽  
Oleg V. Stolbov ◽  
Valter Böhm ◽  
Klaus Zimmermann

Abstract Magnetoactive elastomers (MAEs) are a special type of smart materials consisting of an elastic matrix with embedded microsized particles that are made of ferromagnetic materials with high or low coercivity. Due to their composition, such elastomers possess unique magnetic field-dependent material properties. The present paper compiles the results of investigations on MAEs towards an approach of their potential application as vibrating sensor elements with adaptable sensitivity. Starting with the model-based and experimental studies of the free vibrational behavior displayed by cantilevers made of MAEs, it is shown that the first bending eigenfrequency of the cantilevers depends strongly on the strength of an applied uniform magnetic field. The investigations of the forced vibration response of MAE beams subjected to in-plane kinematic excitation confirm the possibility of active magnetic control of the amplitude-frequency characteristics. With change of the uniform field strength, the MAE beam reveals different steady-state responses for the same excitation, and the resonance may occur at various ranges of the excitation frequency. Nonlinear dependencies of the amplification ratio on the excitation frequency are obtained for different magnitudes of the applied field. Furthermore, it is shown that the steady-state vibrations of MAE beams can be detected based on the magnetic field distortion. The field difference, which is measured simultaneously on the sides of a vibrating MAE beam, provides a signal with the same frequency as the excitation and an amplitude proportional to the amplitude of resulting vibrations. The presented prototype of the MAE-based vibrating unit with the field-controlled “configuration” can be implemented for realization of acceleration sensor systems with adaptable sensitivity. The ongoing research on MAEs is oriented to the use of other geometrical forms along with beams, e.g. two-dimensional structures such as membranes.


1984 ◽  
Vol 24 (06) ◽  
pp. 606-616 ◽  
Author(s):  
Charles P. Thomas ◽  
Paul D. Fleming ◽  
William K. Winter

Abstract A mathematical model describing one-dimensional (1D), isothermal flow of a ternary, two-phase surfactant system in isotropic porous media is presented along with numerical solutions of special cases. These solutions exhibit oil recovery profiles similar to those observed in laboratory tests of oil displacement by surfactant systems in cores. The model includes the effects of surfactant transfer between aqueous and hydrocarbon phases and both reversible and irreversible surfactant adsorption by the porous medium. The effects of capillary pressure and diffusion are ignored, however. The model is based on relative permeability concepts and employs a family of relative permeability curves that incorporate the effects of surfactant concentration on interfacial tension (IFT), the viscosity of the phases, and the volumetric flow rate. A numerical procedure was developed that results in two finite difference equations that are accurate to second order in the timestep size and first order in the spacestep size and allows explicit calculation of phase saturations and surfactant concentrations as a function of space and time variables. Numerical dispersion (truncation error) present in the two equations tends to mimic the neglected present in the two equations tends to mimic the neglected effects of capillary pressure and diffusion. The effective diffusion constants associated with this effect are proportional to the spacestep size. proportional to the spacestep size. Introduction In a previous paper we presented a system of differential equations that can be used to model oil recovery by chemical flooding. The general system allows for an arbitrary number of components as well as an arbitrary number of phases in an isothermal system. For a binary, two-phase system, the equations reduced to those of the Buckley-Leverett theory under the usual assumptions of incompressibility and each phase containing only a single component, as well as in the more general case where both phases have significant concentrations of both components, but the phases are incompressible and the concentration in one phase is a very weak function of the pressure of the other phase at a given temperature. pressure of the other phase at a given temperature. For a ternary, two-phase system a set of three differential equations was obtained. These equations are applicable to chemical flooding with surfactant, polymer, etc. In this paper, we present a numerical solution to these equations paper, we present a numerical solution to these equations for I D flow in the absence of gravity. Our purpose is to develop a model that includes the physical phenomena influencing oil displacement by surfactant systems and bridges the gap between laboratory displacement tests and reservoir simulation. It also should be of value in defining experiments to elucidate the mechanisms involved in oil displacement by surfactant systems and ultimately reduce the number of experiments necessary to optimize a given surfactant system.


2021 ◽  
Author(s):  
Omar Chaabi ◽  
Emad W. Al-Shalabi ◽  
Waleed Alameri

Abstract Low salinity polymer (LSP) flooding is getting more attention due to its potential of enhancing both displacement and sweep efficiencies. Modeling LSP flooding is challenging due to the complicated physical processes and the sensitivity of polymers to brine salinity. In this study, a coupled numerical model has been implemented to allow investigating the polymer-brine-rock geochemical interactions associated with LSP flooding along with the flow dynamics. MRST was coupled with the geochemical software IPhreeqc. The effects of polymer were captured by considering Todd-Longstaff mixing model, inaccessible pore volume, permeability reduction, polymer adsorption as well as salinity and shear rate effects on polymer viscosity. Regarding geochemistry, the presence of polymer in the aqueous phase was considered by adding a new solution specie and related chemical reactions to PHREEQC database files. Thus, allowing for modeling the geochemical interactions related to the presence of polymer. Coupling the two simulators was successfully performed, verified, and validated through several case studies. The coupled MRST-IPhreeqc simulator allows for modeling a wide variety of geochemical reactions including aqueous, mineral precipitation/dissolution, and ion exchange reactions. Capturing these reactions allows for real time tracking of the aqueous phase salinity and its effect on polymer rheological properties. The coupled simulator was verified against PHREEQC for a realistic reactive transport scenario. Furthermore, the coupled simulator was validated through history matching a single-phase LSP coreflood from the literature. This paper provides an insight into the geochemical interactions between partially hydrolyzed polyacrylamide (HPAM) and aqueous solution chemistry (salinity and hardness), and their related effect on polymer viscosity. This work is also considered as a base for future two-phase polymer solution and oil interactions, and their related effect on oil recovery.


Author(s):  
Zhongxian Liu ◽  
Jiaqiao Liu ◽  
Sibo Meng ◽  
Xiaojian Sun

Summary An indirect boundary element method (IBEM) is developed to model the two-dimensional (2D) diffraction of seismic waves by a fluid-filled crack in a fluid-saturated poroelastic half-space, using Green's functions computed considering the distributed loads, flow, and fluid characteristics. The influence of the fluid-filled crack on the diffraction characteristics is investigated by analyzing key parameters, such as the excitation frequency, incident angle, crack width and depth, and medium porosity. The results for the fluid-filled crack model are compared to those for the fluid-free crack model under the same conditions. The numerical results demonstrate that the fluid-filled crack has a significant amplification effect on the surface displacements, and that the effect of the depth of the fluid-filled crack is more complex compared to the influence of other parameters. The resonance diffraction generates an amplification effect in the case of normally incident P waves. Furthermore, the horizontal and vertical displacement amplitudes reach 4.2 and 14.1, respectively. In the corresponding case of the fluid-free crack, the vertical displacement amplitude is only equal to 4.1, indicating the amplification effect of the fluid in the crack. Conversely, for normally incident SV waves at certain resonance frequencies, the displacement amplitudes above a fluid-filled crack may be lower than the displacement amplitudes observed in the corresponding case of a fluid-free crack.


Author(s):  
Christian Rohde ◽  
Lars von Wolff

We consider the incompressible flow of two immiscible fluids in the presence of a solid phase that undergoes changes in time due to precipitation and dissolution effects. Based on a seminal sharp interface model a phase-field approach is suggested that couples the Navier–Stokes equations and the solid’s ion concentration transport equation with the Cahn–Hilliard evolution for the phase fields. The model is shown to preserve the fundamental conservation constraints and to obey the second law of thermodynamics for a novel free energy formulation. An extended analysis for vanishing interfacial width reveals that in this limit the sharp interface model is recovered, including all relevant transmission conditions. Notably, the new phase-field model is able to realize Navier-slip conditions for solid–fluid interfaces in the limit.


2020 ◽  
Vol 21 (2) ◽  
pp. 339
Author(s):  
I. Carneiro ◽  
M. Borges ◽  
S. Malta

In this work,we present three-dimensional numerical simulations of water-oil flow in porous media in order to analyze the influence of the heterogeneities in the porosity and permeability fields and, mainly, their relationships upon the phenomenon known in the literature as viscous fingering. For this, typical scenarios of heterogeneous reservoirs submitted to water injection (secondary recovery method) are considered. The results show that the porosity heterogeneities have a markable influence in the flow behavior when the permeability is closely related with porosity, for example, by the Kozeny-Carman (KC) relation.This kind of positive relation leads to a larger oil recovery, as the areas of high permeability(higher flow velocities) are associated with areas of high porosity (higher volume of pores), causing a delay in the breakthrough time. On the other hand, when both fields (porosity and permeability) are heterogeneous but independent of each other the influence of the porosity heterogeneities is smaller and may be negligible.


Geofluids ◽  
2017 ◽  
Vol 2017 ◽  
pp. 1-11 ◽  
Author(s):  
Yunfeng Dai ◽  
Zhifang Zhou ◽  
Jin Lin ◽  
Jiangbo Han

To describe accurately the flow characteristic of fracture scale displacements of immiscible fluids, an incompressible two-phase (crude oil and water) flow model incorporating interfacial forces and nonzero contact angles is developed. The roughness of the two-dimensional synthetic rough-walled fractures is controlled with different fractal dimension parameters. Described by the Navier–Stokes equations, the moving interface between crude oil and water is tracked using level set method. The method accounts for differences in densities and viscosities of crude oil and water and includes the effect of interfacial force. The wettability of the rough fracture wall is taken into account by defining the contact angle and slip length. The curve of the invasion pressure-water volume fraction is generated by modeling two-phase flow during a sudden drainage. The volume fraction of water restricted in the rough-walled fracture is calculated by integrating the water volume and dividing by the total cavity volume of the fracture while the two-phase flow is quasistatic. The effect of invasion pressure of crude oil, roughness of fracture wall, and wettability of the wall on two-phase flow in rough-walled fracture is evaluated.


2018 ◽  
Vol 96 (9) ◽  
pp. 1016-1028 ◽  
Author(s):  
Pramod Kumar Yadav ◽  
Sneha Jaiswal

The present problem is concerned with two-phase fluid flow through a horizontal porous channel in the presence of uniform inclined magnetic field. The micropolar fluid or Eringen fluid and Newtonian viscous fluid are flowing in the upper and lower regions of the horizontal porous channel, respectively. In this paper, the permeability of each region of the horizontal porous channel has been taken to be different. The effects of various physical parameters like angles of inclination of magnetic field, viscosity ratio, micropolarity parameter, etc., on the velocities, micro-rotational velocity of two immiscible fluids in horizontal porous channel, wall-shear stress, and flow rate have been discussed. The result obtained for immiscible micropolar–Newtonian fluids are compared with the results of two immiscible Newtonian fluids. The obtained result may be used in production of oil from oil reservoirs, purification of contaminated ground water, etc.


Sign in / Sign up

Export Citation Format

Share Document