scholarly journals Quantification of Uncertainty and Best Practice in Computing Interfacial Curvature from Complex Pore Space Images

Materials ◽  
2019 ◽  
Vol 12 (13) ◽  
pp. 2138 ◽  
Author(s):  
Takashi Akai ◽  
Qingyang Lin ◽  
Abdulla Alhosani ◽  
Branko Bijeljic ◽  
Martin Blunt

Recent advances in high-resolution three-dimensional X-ray CT imaging have made it possible to visualize fluid configurations during multiphase displacement at the pore-scale. However, there is an inherited difficulty in image-based curvature measurements: the use of voxelized image data may introduce significant error, which has not—to date—been quantified. To find the best method to compute curvature from micro-CT images and quantify the likely error, we performed drainage and imbibition direct numerical simulations for an oil/water system on a bead pack and a Bentheimer sandstone. From the simulations, local fluid configurations and fluid pressures were obtained. We then investigated methods to compute curvature on the oil/water interface. The interface was defined in two ways; in one case the simulated interface with a sub-resolution smoothness was used, while the other was a smoothed interface extracted from synthetic segmented data based on the simulated phase distribution. The curvature computed on these surfaces was compared with that obtained from the simulated capillary pressure, which does not depend on the explicit consideration of the shape of the interface. As distinguished from previous studies which compared an average or peak curvature with the value derived from the measured macroscopic capillary pressure, our approach can also be used to study the pore-by-pore variation. This paper suggests the best method to compute curvature on images with a quantification of likely errors: local capillary pressures for each pore can be estimated to within 30% if the average radius of curvature is more than 6 times the image resolution, while the average capillary pressure can also be estimated to within 11% if the average radius of curvature is more than 10 times the image resolution.

SPE Journal ◽  
2017 ◽  
Vol 22 (03) ◽  
pp. 940-949 ◽  
Author(s):  
Edo S. Boek ◽  
Ioannis Zacharoudiou ◽  
Farrel Gray ◽  
Saurabh M. Shah ◽  
John P. Crawshaw ◽  
...  

Summary We describe the recent development of lattice Boltzmann (LB) and particle-tracing computer simulations to study flow and reactive transport in porous media. First, we measure both flow and solute transport directly on pore-space images obtained from micro-computed-tomography (CT) scanning. We consider rocks with increasing degree of heterogeneity: a bead pack, Bentheimer sandstone, and Portland carbonate. We predict probability distributions for molecular displacements and find excellent agreement with pulsed-field-gradient (PFG) -nuclear-magnetic-resonance (NMR) experiments. Second, we validate our LB model for multiphase flow by calculating capillary filling and capillary pressure in model porous media. Then, we extend our models to realistic 3D pore-space images and observe the calculated capillary pressure curve in Bentheimer sandstone to be in agreement with the experiment. A process-based algorithm is introduced to determine the distribution of wetting and nonwetting phases in the pore space, as a starting point for relative permeability calculations. The Bentheimer relative permeability curves for both drainage and imbibition are found to be in good agreement with experimental data. Third, we show the speedup of a graphics-processing-unit (GPU) algorithm for large-scale LB calculations, offering greatly enhanced computing performance in comparison with central-processing-unit (CPU) calculations. Finally, we propose a hybrid method to calculate reactive transport on pore-space images by use of the GPU code. We calculate the dissolution of a porous medium and observe agreement with the experiment. The LB method is a powerful tool for calculating flow and reactive transport directly on pore-space images of rock.


2021 ◽  
pp. 105501
Author(s):  
W.H. Wu ◽  
D.G. Eskin ◽  
A. Priyadarshi ◽  
T. Subroto ◽  
I. Tzanakis ◽  
...  

1960 ◽  
Vol 26 (12) ◽  
pp. 1167-1170
Author(s):  
Tuneo IKEUTI ◽  
Wataru SIMIDU
Keyword(s):  

2021 ◽  
Author(s):  
Tianhua Zhang ◽  
Shiduo Yang ◽  
Chandramani Shrivastava ◽  
Adrian A ◽  
Nadege Bize-Forest

Abstract With the advancement of LWD (Logging While Drilling) hardware and acquisition, the imaging technology becomes not only an indispensable part of the drilling tool string, but also the image resolution increases to map layers and heterogeneity features down to less than 5mm scale. This shortens the geological interpretation turn-around time from wireline logging time (hours to days after drilling) to semi-real time (drilling time or hours after drilling). At the same time, drilling motion is complex. The depth tracking is on the surface referenced to the surface block movement. The imaging sensor located downhole can be thousands of feet away from the surface. Mechanical torque and drag, wellbore friction, wellbore temperature and weight on bit can make the downhole sensor movement motion not synchronized with surface pipe depth. This will cause time- depth conversion step generate image artifacts that either stop real-time interpretation of geological features or mis-interpret features on high resolution images. In this paper, we present several LWD images featuring distortion mechanism during the drilling process using synthetic data. We investigated how heave, depth reset and downhole sensor stick/slip caused image distortions. We provide solutions based on downhole sensor pseudo velocity computation to minimize the image distortion. The best practice in using Savitsky-Golay filter are presented in the discussion sections. Finally, some high-resolution LWD images distorted with drilling-related artifacts and processed ones are shown to demonstrate the importance of image post-processing. With the proper processed images, we can minimize interpretation risks and make drilling decisions with more confidence.


2011 ◽  
Vol 239-242 ◽  
pp. 2650-2654
Author(s):  
Fu Chen ◽  
Jie He ◽  
Ping Guo ◽  
Yuan Xu ◽  
Cheng Zhong

According to the mechanisms of carbon dioxide miscible flooding and previous researchers’ work on synthesis of CO2-soluble surfactant, Citric acid isoamyl ester was synthesized, and it’s oil solubility and the rate of viscosity reduction both in oil-water system and oil were evaluated. And then we found that this compound can solve in oil effectively; the optimum mass of Citric acid isoamyl ester introduced in oil-water system is 0.12g when the mass ratio of oil and water is 7:3 (crude oil 23.4g, formation water 10g) and the experimental temperature is 50°C , the rate of viscosity reduction is 47.2%; during the evaluation of the ability of Citric acid isoamyl ester to decrease oil viscosity, we found that the optimum dosage of this compound in 20g crude oil is 0.2g when the temperature is 40°C, and the rate of viscosity reduction is 7.37% at this point.


Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-15 ◽  
Author(s):  
Wenbin Xu ◽  
Zhihui Liu ◽  
Jie Liu ◽  
Yongfei Yang

Well test analysis requires a preselected model, which relies on the context input and the diagnostic result through the pressure logarithmic derivative curve. Transient pressure outer boundary response heavily impacts on the selection of such a model. Traditional boundary-type curves used for such diagnostic purpose are only suitable for single-phase flow in a homogeneous reservoir, while practical situations are often much more complicated. This is particularly true when transient pressure is derived during the field development phase, for example, from permanent down-hole gauge (PDG), where outer boundary condition such as an active aquifer with a transition zone above it plays a big role in dominating the late time pressure response. In this case, capillary pressure and the total mobility in the transition zone have significant effect on the pressure response. This effect is distinctly different for oil-water system and gas water system, which will result in the pressure logarithmic derivatives remarkably different from the traditional boundary-type curves. This paper presents study results derived through theoretical and numerical well testing approaches to solve this problem. The outcome of this study can help in understanding the reservoir behavior and guiding the management of mature field. According to the theoretical development by Thompson, a new approach was derived according to Darcy’s law, which shows that pressure response in the transition zone is a function of total effective mobility. For oil-water system, the total effective mobility increases with an increase in the radius of transition zone, while for gas-water system, the effect is opposite.


Sign in / Sign up

Export Citation Format

Share Document