Determination of multi-component content and construction of digital cores based on CT grey thresholds of altered igneous rocks

Altered igneous reservoirs have low porosity and permeability,compact structure and certain heterogeneity.A simple digital core with certain generality and multi-parameter constraints can be con-structed to characterize the microscopic pore structure and mineral composition.In this paper,based on core X-ray,CT images and whole-rock mineral analysis,threshold segmentation of mass content and grayscale distribution of various minerals in different lithologies of igneous rocks in the buried hill of Huizhou depression is carried out to construct digital core of altered igneous rocks.The results show that after converting the mineral mass content into volume content,the minerals of altered igneous rocks in Huizhou depression can be classified into components.According to the range of grayscale value,components can be divided into six parts.Due to the difference of the content of components in different lithologies of igneous rocks,differentiated grayscale threshold segmentation is needed to obtain the digital core for a single lithology.The final digital core generation process includes two steps:building a single component digital core,and stacking and combining.This kind of universal digital core model can support the subsequent pore scale numerical simulation and comprehensive rock physics research.

Elastic Characteristics of Digital Cores from Longmaxi Shale Using Lattice Spring Models

Effective medium methods for the attribution of micro-structures to macro elastic properties of shales are important for the prediction of sweet spots in the shale-gas production.With X-ray micro-computed tomography(XMCT),the micro-structures of shale core samples from Longmaxi Formation are visualized and characterized by 3D digital images.As an efficient alternative to conventional effective medium methods for estimating elastic properties,we propose a consistent workflow of lattice spring modeling(LSM)to emulate the digital cores using three types of lattices.Particular attention is paid to investigate the effective Young’s moduli,Poisson’s ratios,and preferred orientations,by uniaxial compression tests along two directions.Within elastic deformation,the impact of lattice arrangements on the anisotropy is even more than those of stress disturbances and micro-structural features.Compared with analytical approximations and theoretical predictions,the LSM numerical scheme shows general applicability for heterogeneous porous rocks.

Poroelastic finite-difference modeling for ultrasonic waves in digital porous cores

Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with smallscale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference （FD） method of Biot＇s poroelastic equations, with an arbitrary even-order （2L） accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer （CPML） absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultra- sonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thick- nesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Qsc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.

Relationship between rock uniaxial compressive strength and digital core drilling parameters and its forecast method

The rock uniaxial compressive strength(UCS)is the basic parameter for support designs in underground engineering.In particular,the rock UCS should be obtained rapidly for underground engineering with complex geological conditions,such as soft rock,fracture areas,and high stress,to adjust the excavation and support plan and ensure construction safety.To solve the problem of obtaining real-time rock UCS at engineering sites,a rock UCS forecast idea is proposed using digital core drilling.The digital core drilling tests and uniaxial compression tests are performed based on the developed rock mass digital drilling system.The results indicate that the drilling parameters are highly responsive to the rock UCS.Based on the cutting and fracture characteristics of the rock digital core drilling,the mechanical analysis of rock cutting provides the digital core drilling strength,and a quantitative relationship model(CDP-UCS model)for the digital core drilling parameters and rock UCS is established.Thus,the digital core drilling-based rock UCS forecast method is proposed to provide a theoretical basis for continuous and quick testing of the surrounding rock UCS.

Study of the numerical simulation of tight sandstone gas molecular diffusion based on digital core technology

Diffusion is an important mass transfer mode of tight sandstone gas. Since nano-pores are extensively developed in the interior of tight sandstone, a considerable body of research indicates that the type of diffusion is mainly molecular diffusion based on Fick＇s law. However, accurate modeling and understanding the physics of gas transport phenomena in nanoporous media is still a challenge for researchers and traditional investigation（analytical and experimental methods） have many limitations in studying the generic behavior. In this paper, we used Nano-CT to observe the pore structures of samples of the tight sandstone of western of Sichuan. Combined with advanced image processing technology, threedimensional distributions of the nanometer-sized pores were reconstructed and a tight sandstone digital core model was built, as well the pore structure parameters were analyzed quantitatively. Based on the digital core model, the diffusion process of methane molecules from a higher concentration area to a lower concentration area was simulated by a finite volume method. Finally, the reservoir＇s concentration evolution was visualized and the intrinsic molecular diffusivity tensor which reflects the diffusion capabilities of this rock was calculated. Through comparisons, we found that our calculated result was in good agreement with other empirical results. This study provides a new research method for tight sandstone digital rock physics. It is a foundation for future tight sandstone gas percolation theory and numerical simulation research.

Investigating microscopic seepage characteristics and fracture effectiveness of tight sandstones:a digital core approach

Microscopic seepage characteristics are critical for the evaluation of tight sandstone reservoirs.In this study,a digital core approach integrating microscopic seepage simulation and CT scanning was developed to characterize microscopic seepage and fracture effectiveness(the ratio of micro-fractures that contributes to fluid flow)of tight sandstones.Numerical simulations were carried out for characterizations of tight sandstones.The results show that the axial permeability of the investigated cylindrical tight sandstone from Junggar Basin in China is 0.460μm~2,while the radial permeability is 0.3723μm~2,and the axial and radial effective fracture ratios are 0.4387 and 0.4806,respectively,indicating that cracks are not fully developed and the connectivity between micro-cracks is poor.Directional permeability that is difficult to measure by laboratory experiments can be obtained readily using the proposed method in this paper.The results provide important information for improving the exploration and development of tight sandstone reservoirs.

Investigation of image segmentation effect on the accuracy of reconstructed digital core models of coquina carbonate

Digital core models reconstructed using X-ray tomography(X-CT)enable the quantitative characterization of the pore structure in three dimensions(3D)and the numerical simulation of petrophysics.When the X-CT images accurately reflect the micro structures of core samples,the greyscale threshold in the image segmentation determines the accuracy of digital cores and the simulated petrophysical properties.Therefore,it is vital to investigate the comparison parameter for determining the key greyscale threshold and the criterion to describe the accuracy of the segmentation.Representative coquina digital core models from X-CT are used in this work to study the impact of grayscale threshold on the porosity,pore percolation,connectivity and electrical resistivity of the pore scale model and these simulations are calculated by Minkowski functions,component labeling and fi nite element method,respectively,to quantify the pore structure and simulate electrical resistivity.Results showed that the simulated physical properties of the digital cores,varied with the gradual increase of the greyscale threshold.Among the four parameters related to the threshold,the porosity was most sensitive and chose as the comparison parameter to judge the accuracy of the greyscale threshold.The variations of the threshold change the micro pore structures,and then the electrical resistivity.When the porosity of the digital core model is close to the experimental porosity,the simulated porosity exponent matches the experimental porosity exponents well.The good agreement proved that the porosity is the critical comparison parameter to describe the accuracy of image segmentation.The criterion is that the porosity of the digital core after segmentation should be close to the experimental porosity.

Reconstructing the 3D digital core with a fully convolutional neural network

In this paper, the complete process of constructing 3D digital core by fullconvolutional neural network is described carefully. A large number of sandstone computedtomography (CT) images are used as training input for a fully convolutional neural networkmodel. This model is used to reconstruct the three-dimensional (3D) digital core of Bereasandstone based on a small number of CT images. The Hamming distance together with theMinkowski functions for porosity, average volume specifi c surface area, average curvature,and connectivity of both the real core and the digital reconstruction are used to evaluate theaccuracy of the proposed method. The results show that the reconstruction achieved relativeerrors of 6.26%, 1.40%, 6.06%, and 4.91% for the four Minkowski functions and a Hammingdistance of 0.04479. This demonstrates that the proposed method can not only reconstructthe physical properties of real sandstone but can also restore the real characteristics of poredistribution in sandstone, is the ability to which is a new way to characterize the internalmicrostructure of rocks.

3D characterization of porosity and minerals of low-permeability uranium-bearing sandstone based on multi-resolution image fusion

In the process of in situ leaching of uranium,the microstructure controls and influences the flow distribution,percolation characteristics,and reaction mechanism of lixivium in the pores of reservoir rocks and directly affects the leaching of useful components.In this study,the pore throat,pore size distribution,and mineral composition of low-permeability uranium-bearing sandstone were quantitatively analyzed by high pressure mercury injection,nuclear magnetic resonance,X-ray diffraction,and wavelength-dispersive X-ray fluorescence.The distribution characteristics of pores and minerals in the samples were qualitatively analyzed using energy-dispersive scanning electron microscopy and multi-resolution CT images.Image registration with the landmarks algorithm provided by FEI Avizo was used to accurately match the CT images with different resolutions.The multi-scale and multi-mineral digital core model of low-permeability uranium-bearing sandstone is reconstructed through pore segmentation and mineral segmentation of fusion core scanning images.The results show that the pore structure of low-permeability uranium-bearing sandstone is complex and has multi-scale and multi-crossing characteristics.The intergranular pores determine the main seepage channel in the pore space,and the secondary pores have poor connectivity with other pores.Pyrite and coffinite are isolated from the connected pores and surrounded by a large number of clay minerals and ankerite cements,which increases the difficulty of uranium leaching.Clays and a large amount of ankerite cement are filled in the primary and secondary pores and pore throats of the low-permeability uraniumbearing sandstone,which significantly reduces the porosity of the movable fluid and results in low overall permeability of the cores.The multi-scale and multi-mineral digital core proposed in this study provides a basis for characterizing macroscopic and microscopic pore-throat structures and mineral distributions of low-permeability uranium-bearing sandstone and can better understand the seepage characteristics.

Numerical simulation of non-Archie electrophysical property of saturated rock with lattice Boltzmann method

The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.

Calculation of oil saturation in clay-rich shale reservoirs:A case study of Qing 1 Member of Cretaceous Qingshankou Formation in Gulong Sag,Songliao Basin,NE China

The targeted reservoir,which is referred as the first member of Cretaceous Qingshankou Formation in Gulong Sag,Songliao Basin,NE China,is characterized by the enrichment of clay and lamellation fractures.Aiming at the technical challenge of determining oil saturation of such reservoir,nano-pores were accurately described and located through focused ion beam scanning electron microscopy and quantitative evaluation of minerals by scanning electron microscopy based on Simandoux model,to construct a 4D digital core frame.Electrical parameters of the shale reservoir were determined by finite element simulation,and the oil saturation calculation method suitable for shale was proposed.Comparison between the results from this method with that from real core test and 2D nuclear magnetic log shows that the absolute errors meet the requirements of the current reserve specification in China for clay-rich shale reservoir.Comparison analysis of multiple wells shows that the oil saturation values calculated by this method of several points vertically in single wells and multiple wells on the plane are in agreement with the test results of core samples and the regional deposition pattern,proving the accuracy and applicability of the method model.

Wormholes propagation for fractured-vuggy formation:Laboratory tests,numerical simulation and field application

The propagation of wormhole is vital important for matrix acidizing and acid fracturing in carbonate reservoirs.While the formation of acid dissolved wormhole is derived from heterogeneous physical and chemical transportations and reactions.Alveolate dissolved pores,krast caves,and natural fissures are the major reservoir spaces for the Sinian dolomite formation in the Anyue gas field of the Sichuan Basin.There were four categories of formation,which are matrix dominated,inter-breccia dissolved pore dominated,dissolved pore and cave dominated,and fissure and cave dominated,based on the development intensity and connectedness of caves and fissures.The caves and fissures make the wormhole formation and propagation particularly complicated.Firstly,the 3-D topological structure of dissolved pores,vugs,fissures and throats inside cores is quantitatively scanned by CT imaging technology for its feature of vivid and damage-free.Secondly,3-D patterns of wormhole are obtained with CT scanning after core flooding by acid.Additionally,the porethroat network model is reconstructed with digital cores technology.Then,the size and ratio of pore and throat before and after core flooding by acid is analyzed and the absolute permeability of pore scale flow is numerically simulated to understand the fundamental influence of pores and vugs distribution and connectedness on wormhole propagation.Lastly,the wormhole pattern gained by CT scanning and simulating with two-scale model is compared.Meanwhile,the corrected two-scale model is utilized to simulate the wormhole propagation for matrix acidizing and acid fracturing of Sinian fractured-vuggy dolomite in Anyue gas field,Sichuan Basin.The optimized injection rate and volume were in agreement with the characteristic matrix acidizing operating curve,which indicates that the two-scale model was suitable for matrix acidizing optimization design of such formations.In addition,the simulated acid etched fracture length with considering the dynamic wormhole leakoff was consistent with the well testing interpreted result.

A novel pore-fracture dual network modeling method considering dynamic cracking and its applications

Unconventional reservoirs are normally characterized by dual porous media, which has both multi-scalepore and fracture structures, such as low permeability or tight oil reservoirs. The seepage characteristicsof such reservoirs is mainly determined by micro-fractures, but conventional laboratory experimentalmethods are difficult to measure it, which is attribute to the dynamic cracking of these micro-fractures.The emerging digital core technology in recent years can solve this problem by developing an accuratepore network model and a rational simulation approach. In this study, a novel pore-fracture dualnetwork model was established based on percolation theory. Fluid flow in the pore of two scales, microfracture and matrix pore, were considered, also with the impact of micro-fracture opening and closingduring flow. Some seepage characteristic parameters, such as fluid saturations, capillary pressure, relative permeabilities, displacement efficiency in different flow stage, can be predicted by proposedcalculating method. Through these work, seepage characteristics of dual porous media can be achieved.