scholarly journals Numerical simulations of sand production in oil wells using the CFD-DEM-IBM approach

2020 ◽  
Author(s):  
Aigerim Rakhimzhanova ◽  
Colin Thornton ◽  
Yerlan Amanbek ◽  
Yong Zhao

Three-dimensional CFD-DEM-IBM simulations of sand production in a sandstone formation, using periodic boundaries, were performed using 10000 frictional elastic spheres bonded together and compressed at 1 MPa of overburden pressure. Sand production simulation geometry and procedure are proposed, in which the cone penetration test (CPT) has been used to investigate a physical perforation penetration of the cemented sandstone material with the real-world grain size distribution from the Ustyurt-Buzachi Sedimentary Basin. The Immersed Boundary Method (IBM) was adapted for the sand production simulation geometry to simulate the fluid flow near the well casing. Oil with low viscosity and density was used as an injection fluid (reservoir fluid). Erosion near the perforation tunnel due to the pressure drawdown was examined, where the production of sand particles was initiated during the first flow due to the drag force that lifted the sand particles from the perforation damage zone. At the beginning of the simulation a sand arch was captured around the perforation tunnel and due to the fluid flow it collapsed and perforation cavity became larger. The amount and mass of produced sand particles were calculated.

2000 ◽  
Vol 3 (02) ◽  
pp. 127-131 ◽  
Author(s):  
J.J. Zhang ◽  
C.S. Rai ◽  
C.H. Sondergeld

Summary The mechanical strength of a reservoir formation is the most crucial information required for predicting sand production and recommending sand control completion. So far, the only reliable technique to obtain the formation strength quantitatively is to perform laboratory tests on core samples. The laboratory tests require substantial volumes of cores, which in most cases are not available. In this study, we present a new method to avoid this restriction. The significance of this approach is in its simplicity and efficiency in constructing a reliable mechanical failure envelope. The key results of this study, based on measurements on a variety of sandstones, are the following.A single normalized failure envelope characterizes sandstone formations. This universal curve makes it possible to construct the failure envelope for a sandstone formation from the knowledge of critical pressure.There exists a correlation between the critical pressure and the compressional wave velocity (at equivalent depths of burial).The failure envelope for a sandstone formation can then be constructed simply from compressional wave velocities. These velocities are generally accessible from conventional logging data. Introduction A major problem encountered during hydrocarbon production is the influx of sand, or sand production.1 It can cause severe damage to both production equipment and the producing formation. Furthermore, remediation processes after sanding are extremely difficult, costly or often impossible. Although significant research has been conducted on sand production, we are still at an embryonic stage in predicting sand influx. Several factors determine sand production. The most critical factors are (1) formation strength; (2) in-situ stress; and (3) production rate. The hydrocarbon production process is associated with reservoir depletion, which results in a decrease of reservoir pore pressure. Consequently, the effective overburden pressure, defined as total overburden pressure minus pore pressure, increases. Formation collapse is most likely if the effective stress exceeds the formation strength. In addition, production rate increase, which is associated with large fluid pressure gradients near the borehole, tends to draw the sand into the wellbore. Generally, one can estimate the in-situ stress. For example, the horizontal minimum stress can be measured from hydraulic fracture testing,2-4 and the overburden pressure from overburden density data. The production rate is a controllable parameter. The parameter of concern is the formation strength, which is the focus of this study. The most reliable technique for obtaining mechanical strength data is triaxial testing of core samples in the laboratory. With appropriate arrangements of applied stresses one can determine a failure envelope in stress space. Such a failure envelope quantifies the stress conditions under which the material fails. Although the laboratory test can provide dependable mechanical strength data, it is not followed routinely simply because it is time-consuming and costly. Moreover, in most cases, a sufficient amount of core is not available. Traditionally the mechanical strength, or Mohr-Coulomb failure criterion, is estimated from P- and S-wave velocities and density log data5-8 based on the correlation of Deere and Miller.9 This approach estimates uniaxial compressive strength, and assumes a constant frictional angle. In this study, we seek an alternative method of estimating overall nonlinear mechanical strength in a three-dimensional stress space. Basic Concept of Failure Envelope Formation collapse is an indication that the in-situ stress is beyond the failure limit of the formation material. This failure stress limit is a quantitative parameter that defines the formation mechanical strength. For a one-dimensional state of stress, the mechanical strength can be simply quantified with a single parameter: the uniaxial compressive strength. However, because the in-situ formation stress is three dimensional and anisotropic, a more complicated mathematical expression involving all the stresses is required to quantify the mechanical strength. This quantitative expression of mechanical strength is known as the failure envelope or failure criterion.


Author(s):  
Fred E. Hossler

Preparation of replicas of the complex arrangement of blood vessels in various organs and tissues has been accomplished by infusing low viscosity resins into the vasculature. Subsequent removal of the surrounding tissue by maceration leaves a model of the intricate three-dimensional anatomy of the blood vessels of the tissue not obtainable by any other procedure. When applied with care, the vascular corrosion casting technique can reveal fine details of the microvasculature including endothelial nuclear orientation and distribution (Fig. 1), locations of arteriolar sphincters (Fig. 2), venous valve anatomy (Fig. 3), and vessel size, density, and branching patterns. Because casts faithfully replicate tissue vasculature, they can be used for quantitative measurements of that vasculature. The purpose of this report is to summarize and highlight some quantitative applications of vascular corrosion casting. In each example, casts were prepared by infusing Mercox, a methyl-methacrylate resin, and macerating the tissue with 20% KOH. Casts were either mounted for conventional scanning electron microscopy, or sliced for viewing with a confocal laser microscope.


2021 ◽  
pp. 105678952110286
Author(s):  
H Zhang ◽  
J Woody Ju ◽  
WL Zhu ◽  
KY Yuan

In a recent companion paper, a three-dimensional isotropic elastic micromechanical framework was developed to predict the mechanical behaviors of the innovative asphalt patching materials reinforced with a high-toughness, low-viscosity nanomolecular resin, dicyclopentadiene (DCPD), under the splitting tension test (ASTM D6931). By taking advantage of the previously proposed isotropic elastic-damage framework and considering the plastic behaviors of asphalt mastic, a class of elasto-damage-plastic model, based on a continuum thermodynamic framework, is proposed within an initial elastic strain energy-based formulation to predict the behaviors of the innovative materials more accurately. Specifically, the governing damage evolution is characterized through the effective stress concept in conjunction with the hypothesis of strain equivalence; the plastic flow is introduced by means of an additive split of the stress tensor. Corresponding computational algorithms are implemented into three-dimensional finite elements numerical simulations, and the outcomes are systemically compared with suitably designed experimental results.


2021 ◽  
Vol 233 ◽  
pp. 109189
Author(s):  
Bin Yan ◽  
Wei Bai ◽  
Sheng-Chao Jiang ◽  
Peiwen Cong ◽  
Dezhi Ning ◽  
...  

1999 ◽  
Vol 7 ◽  
pp. 408-417 ◽  
Author(s):  
J. H. Strickland ◽  
L. A. Gritzo ◽  
R. S. Baty ◽  
G. F. Homicz ◽  
S. P. Burns

Author(s):  
Mohammad Mehdi Tavakol ◽  
Mohammad Eslami

Fluid flow around single or multiple bluff bodies mounted on a surface has great significance in science and engineering. Understanding the characteristics of different vortices formed around wall-mounted bodies is quite necessary for different applications. Although the case of a single surface mounted cube has been studied extensively, only little attention has been paid to the flow around two or more rectangular blocks in array. Therefore, a CFD code is developed to calculate three dimensional steady state laminar fluid flow around two cuboids of arbitrary size and configuration mounted on a surface in free stream conditions. The employed numerical scheme is finite volume and SIMPLE algorithm is used to treat pressure and velocity coupling. Results are presented for two rectangular blocks of the different size mounted on a surface in various inline arrangements. Streamlines are plotted for blocks of different size ratio. Velocity and pressure distributions are also plotted in the wake region behind the obstacles. It is shown that how the behavior of flow field and vortical structures depend on the respective size and location of the larger block in comparison with the case of two inline wall mounted cubes of the same size.


1999 ◽  
Vol 28 ◽  
pp. 83-89 ◽  
Author(s):  
A. J. Payne ◽  
D.J. Baldwin

AbstractThis work attempts to explain the fan-like landform assemblages observed in satellite images of the area covered by the former Scandinavian ice sheet (SIS). These assemblages have been interpreted as evidence of large ice streams within the SIS. If this interpretation is correct, then it calls into doubt current theories on the formation of ice streams. These theories regard soft sediment and topographic troughs as being the key determinants of ice-stream location. Neither can be used to explain the existence of ice streams on the flat, hard-rock area of the Baltic Shield. Initial results from a three-dimensional, thermomechanical ice-sheet model indicate that interactions between ice flow, form and temperature can create patterns similar to those mentioned above. The model uses a realistic, 20 km resolution gridded topography and a simple parameterization of accumulation and ablation. It produces patterns of maximum ice-sheet extent, which are similar to those reconstructed from the area’s glacial geomorphology. Flow in the maximum, equilibrium ice sheet is dominated by wedges of warm, low-viscosity, fast-flowing ice. These are separated by areas of cold, slow-flowing ice. This patterning appears to develop spontaneously as the modelled ice sheet grows.


Sign in / Sign up

Export Citation Format

Share Document