Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression an...Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression and permeability tests to investigate the mechanical and seepage properties of tight sandstone.A digital core of tight sandstone was built using Computed Tomography(CT)scanning,which was divided into matrix and pore phases by a pore equivalent diameter threshold.A fluid-solid coupling model was established to investigate the seepage characteristics at micro-scale.The results showed that increasing the confining pressure decreased porosity,permeability,and flow velocity,with the pore phase becoming the dominant seepage channel.Cracks and large pores closed first under increasing pressure,resulted in a steep drop in permeability.However,permeability slightly decreased under high confining pressure,which followed a first-order exponential function.Flow velocity increased with seepage pressure.And the damage mainly occurred in stress-concentration regions under low seepage pressure.Seepage behavior followed linear Darcy flow,the damage emerged at seepage entrances under high pressure,which decreased rock elastic modulus and significantly increased permeability.展开更多
The pore structures of the Majiagou Formation in the Ordos Basin are complex,featuring micro-and nano-scale intra-crystalline and inter-crystalline pores that significantly impact hydrocarbon storage and flow.Precisel...The pore structures of the Majiagou Formation in the Ordos Basin are complex,featuring micro-and nano-scale intra-crystalline and inter-crystalline pores that significantly impact hydrocarbon storage and flow.Precisely characterizing the rock internal structures is crucial for reservoir exploration and development.However,it is difficult to accurately characterize the pore structure of rock using traditional imaging methods to meet the simulation requirements.In this context,this study focuses on high-resolution 3D digital core reconstruction using the SliceGAN model.Specifically,the Modular Automated Processing System(MAPS)image and Quanti-tative Evaluation of Minerals by Scanning Electron Microscopy(QEMSCAN)image were combined to divide MAPS into three categories:pore,dolomite,and calcite.Then,through the SliceGAN algorithm,the 3D digital core was reconstructed.To evaluate the reconstruction,the auto-correlation function,two-point probability function,porosity,mineral content,and specific surface area were employed.The results show that the SliceGAN can effectively capture the micro-features in the core,and the internal structure of the generated core was consistent with that of the original core.This study provided a new sight for reconstructing cores with complex pore structures and strong heterogeneity and innovatively supports tight carbonate reservoir characterization and evaluation.展开更多
In this paper, we obtained three dimensional digital cores using X-ray CT to describe the rock microstructure and applied the open morphology algorithm to simulate oil and formation water distribution in the pore spac...In this paper, we obtained three dimensional digital cores using X-ray CT to describe the rock microstructure and applied the open morphology algorithm to simulate oil and formation water distribution in the pore space at different water saturations during the oil-displacing water flood process. The resistivity, formation factor, and resistivity index of rocks were calculated using the finite element method （FEM） and we studied the effect of rock wettability on electrical properties. The numerical simulation results indicate that the simulated formation factor and resistivity index of the water wet rock agrees well with experiments over the whole range of water saturation and extends the traditional resistivity experiment. The rock wettablilty has a large influence on the rock resistivity index. The resistivity and saturation exponent of oil wet rock are obviously larger than three of water wet rock.展开更多
The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Lo...The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element （FE） methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity （0.75） and kerogen content （ca. ≤3%） over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity （4.5%） of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.展开更多
This study sheds light on how pore structure characteristics and varying dynamic pressure conditions influence the permeability of tight sandstone reservoirs,with a particular focus on the Paleozoic reservoirs in the ...This study sheds light on how pore structure characteristics and varying dynamic pressure conditions influence the permeability of tight sandstone reservoirs,with a particular focus on the Paleozoic reservoirs in the Qingshimao Gas Field.Using CT scans of natural core samples,a three-dimensional digital core was constructed.The maximum ball method was applied to extract a related pore network model,and the pore structure characteristics of the core samples,such as pore radius,throat radius,pore volume,and coordination number,were quantitatively evaluated.The analysis revealed a normally distributed pore radius,suggesting a high degree of reservoir homogeneity and favorable conditions for a connected pore system.However,it was found that the majority of throat radii measured less than 1μm,which severely restricted fluid flow and diminished permeability.Over 50%of the pores measured under 100μm^(3),further constraining fluid movement.Additionally,30%-50%of the pore network was composed of isolated and blind-end pores,which significantly impaired formation connectivity and reduced permeability.Based on this,the lattice Boltzmann method(LBM)was used for pore-scale flow simulation to investigate the influence mechanism of pore structure characteristics and dynamic-static parameters such as displacement pressure difference on the permeability performance of the considered tight sandstone reservoirs for various pressure gradients(0.1,1,and 10 MPa).The simulations revealed a strong relationship between pressure differential and both the number of streamlines and flow path tortuosity.When the pressure differential increased to 1 MPa,30 streamlines were observed,with a tortuosity factor of 1.5,indicating the opening of additional seepage channels and the creation of increasingly winding flow paths.展开更多
A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method whi...A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.展开更多
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 geologica...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.展开更多
The conventional digital core models are usually small in size and have difficulty in representing the complex structures of heterogeneous rocks;Therefore,the parameters of simulated rock physics are difficult to be r...The conventional digital core models are usually small in size and have difficulty in representing the complex structures of heterogeneous rocks;Therefore,the parameters of simulated rock physics are difficult to be referenced.In this study,we propose a feasible simulation method for obtaining multi-scale and multi-component digital cores based on three types of sandstone samples.In the proposed method,the plug and subplug samples are scanned via micro-computed tomography at different resolutions.Furthermore,the images are precisely registered using the proposed hybrid image registration method.In case of high-resolution images,the traditional segmentation method is used to segment the cores into pores and minerals.Subsequently,we established the relations between the gray values and the porosity/mineral content in case of the low-resolution images based on the registered domains and the relation curves were applied to the segmentation of the low-resolution images.The core images constitute the multi-scale and multi-component digital core models after segmentation.Further,the elastic properties of the three samples were simulated at both fine and coarse scales based on the multi-scale and multi-component digital core models,and four component models were considered.The results show that the multi-scale and multi-component digital core models can overcome the representative limits of the conventional digital core models and accurately characterize pores and minerals at different scales.The numerical results of the elastic modulus are more representative at large scales,and considerably reliable results can be obtained by appropriately considering the minerals.展开更多
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 acc...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.展开更多
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 chara...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.展开更多
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 ...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.展开更多
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...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.展开更多
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...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.展开更多
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 characteri...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.展开更多
Research and development of water lock inhibiting measures is very crucial in verifying the link mechanism between the internal factors of water lock and its extent of damage.Based on conventional water-lock physics e...Research and development of water lock inhibiting measures is very crucial in verifying the link mechanism between the internal factors of water lock and its extent of damage.Based on conventional water-lock physics experiments,however,only the consequence of macro water lock damage can be investigated,while the microscopic mechanism cannot be studied.In this paper,3D digital cores of low-permeability sandstones were prepared by means of high-resolution micro-CT scan,and their equivalent pore network model was built as well.Virtual“imbibition”experiments controlled by capillary force were carried out by using pore-scale flow simulation.Then the link mechanism between the microscopic internal factors(e.g.wettability,water saturation and pore-throat structure parameters)and the water-lock damage degree was discussed.It is shown that the damage degree of water lock reduces gradually as the wettability transits from water wet to gas wet.Therefore,the water lock damage can be reduced effectively and gas-well productivity can be improved so long as the capillary environment is changed from strong water wettability to weak gas wettability.The more different the initial water saturation is from the irreducible water saturation,the more serious the water lock damage is.The damage degree of water lock is in a negative correlation with the coordinate number,but a positive correlation with the poreethroat ratio.Based on the existing research results,water lock tends to form in the formations composed of mediumsized throats.It is concluded that there is a critical throat radius,at which the water lock is the most serious.展开更多
The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the res...The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the research included studying the dependence of rock permeability on the stress-strain state in the vicinity of the wells,and physical modeling of the implementation of the method of increasing the permeability of the wellbore zone-the method of directional unloading of the reservoir.The digital part of the research included computed tomography(CT)-based computer analysis of the internal structure,pore space characteristics,and filtration properties before and after the tests.According to the results of physical modeling of deformation and filtration processes,it is found that the permeability of rocks before fracture depends on the stress-strain state insignificantly,and this influence is reversible.However,when downhole pressure reaches 7-8 MPa,macrocracks in the rock begin to grow,accompanied by irreversible permeability increase.Porosity,geodesic tortuosity and permeability values were obtained based on digital studies and numerical modeling.A weak degree of transversal anisotropy of the filtration properties of rocks was detected.Based on the analysis of pore size distribution,pressure field and flow velocities,high homogeneity and connectivity of the rock pore space is shown.The absence of pronounced changes in pore space characteristics and pore permeability after non-uniform triaxial loading rocks was shown.On the basis of geometrical analysis of pore space,the reasons for weak permeability anisotropy were identified.The filtration-capacitance properties obtained from the digital analysis showed very good agreement with the results of field and laboratory measurements.The physical modeling has confirmed the efficiency of application of the directional unloading method for the reservoir under study.The necessary parameters of its application were calculated:bottomhole geometry,stage of operation,stresses and pressure drawdown value.展开更多
Low-salinity waterflooding,as a promising enhanced oil recovery method,has exhibited exciting results in various experiments conducted at different scales.For carbonate rock,pore-scale understanding of the fluid distr...Low-salinity waterflooding,as a promising enhanced oil recovery method,has exhibited exciting results in various experiments conducted at different scales.For carbonate rock,pore-scale understanding of the fluid distribution and remaining oil after low-salinity waterflooding is essential,especially the geometry and topology analysis of oil clusters.We performed the tertiary low-salinity waterflooding and employed X-ray micro-CT to probe the pore-scale displacement mechanism,fluid configuration,oil recovery,and remaining oil distribution.We found that the core becomes less oil-wet after low-salinity waterflooding.Furthermore,we analyzed the oil-rock and oil-brine interfacial areas to further support the wettability alteration.By comparing images after high-salinity waterflooding and low-salinity waterflooding,it is proven that wettability alteration has a significant impact on the behavior of the two-phase flow.Our research demonstrates that low-salinity waterflooding is an effective tertiary enhanced oil recovery technology in carbonate,which changes the wettability of rock and results in less film and singlet oil.展开更多
During the development of tight reservoirs,as the formation pressure decreases,the effective stress on the rock matrix increases.This leads to rock deformation and,consequently,changes in the physical parameters of th...During the development of tight reservoirs,as the formation pressure decreases,the effective stress on the rock matrix increases.This leads to rock deformation and,consequently,changes in the physical parameters of the reservoirs,affecting reservoir productivity.Therefore,it is crucial to identify the influence of stress variation on the pore structure of tight reservoirs at the pore scale.Based on 3D grayscale images of representative samples under different confining pressures using micrometer computed tomography,this study obtained corresponding 3D digital cores through binary segmentation using the maximum class spacing algorithm.Based on digital cores with the same physical size under different confining pressures,the binary data volume was subtracted using a Boolean algorithm to obtain the pore variation space under different confining pressures.The pore variation and permeability damage rate under different confining pressures were then determined.Pore network models of digital cores under different confining pressures were extracted using the maximal ball algorithm,and the corresponding pore-throat radius distribution,coordination number distribution,length/diameter ratio,and shape factor distribution were calculated.The results show that as the confining pressure increases,the pore-throat distribution curve shifts to the direction of smaller pore-throats,the average coordination number decreases,and the throat length/diameter ratio distribution curve and shape factor distribution curve shift to the direction of larger pore-throats.This is because compaction reduces pores and throats,decreases the overall pore-throat connectivity,narrows the throats,and makes the crosssections of the pore-throats more circular.This research serves as a basic platform for efficient development under stress sensitivity in tight reservoirs,and its findings have important academic significance and practical application value.展开更多
Coral reef limestone(CRL)constitutes a distinctive marine carbonate formation with complex mechanical properties.This study investigates the multiscale damage and fracture mechanisms of CRL through integrated experime...Coral reef limestone(CRL)constitutes a distinctive marine carbonate formation with complex mechanical properties.This study investigates the multiscale damage and fracture mechanisms of CRL through integrated experimental testing,digital core technology,and theoretical modelling.Two CRL types with contrasting mesostructures were characterized across three scales.Macroscopically,CRL-I and CRL-II exhibited mean compressive strengths of 8.46 and 5.17 MPa,respectively.Mesoscopically,CRL-I featured small-scale highly interconnected pores,whilst CRL-II developed larger stratified pores with diminished connectivity.Microscopically,both CRL matrices demonstrated remarkable similarity in mineral composition and mechanical properties.A novel voxel average-based digital core scaling methodology was developed to facilitate numerical simulation of cross-scale damage processes,revealing network-progressive failure in CRL-I versus directional-brittle failure in CRL-II.Furthermore,a damage statistical constitutive model based on digital core technology and mesoscopic homogenisation theory established quantitative relationships between microelement strength distribution and macroscopic mechanical behavior.These findings illuminate the fundamental mechanisms through which mesoscopic structure governs the macroscopic mechanical properties of CRL.展开更多
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...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.展开更多
基金financially supported by the National Natural Science Foundation of China(Nos.42272153 and 42472195)the Research Fund of PetroChina Tarim Oilfield Company(No.671023060003)the Research Fund of China National Petroleum Corporation Limited(No.2023ZZ16YJ04).
文摘Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression and permeability tests to investigate the mechanical and seepage properties of tight sandstone.A digital core of tight sandstone was built using Computed Tomography(CT)scanning,which was divided into matrix and pore phases by a pore equivalent diameter threshold.A fluid-solid coupling model was established to investigate the seepage characteristics at micro-scale.The results showed that increasing the confining pressure decreased porosity,permeability,and flow velocity,with the pore phase becoming the dominant seepage channel.Cracks and large pores closed first under increasing pressure,resulted in a steep drop in permeability.However,permeability slightly decreased under high confining pressure,which followed a first-order exponential function.Flow velocity increased with seepage pressure.And the damage mainly occurred in stress-concentration regions under low seepage pressure.Seepage behavior followed linear Darcy flow,the damage emerged at seepage entrances under high pressure,which decreased rock elastic modulus and significantly increased permeability.
基金financially supported by the open fund of Key Lab-oratory of Exploration Technologies for Oil and Gas Resources(Yangtze University),Ministry of Education,NO PI2023-03the open foundation of the National Engineering Laboratory for Exploration and Develop-ment of Low-Permeability Oil&Gas Fields and the National Natural Science Foundation of China(No.42474159).
文摘The pore structures of the Majiagou Formation in the Ordos Basin are complex,featuring micro-and nano-scale intra-crystalline and inter-crystalline pores that significantly impact hydrocarbon storage and flow.Precisely characterizing the rock internal structures is crucial for reservoir exploration and development.However,it is difficult to accurately characterize the pore structure of rock using traditional imaging methods to meet the simulation requirements.In this context,this study focuses on high-resolution 3D digital core reconstruction using the SliceGAN model.Specifically,the Modular Automated Processing System(MAPS)image and Quanti-tative Evaluation of Minerals by Scanning Electron Microscopy(QEMSCAN)image were combined to divide MAPS into three categories:pore,dolomite,and calcite.Then,through the SliceGAN algorithm,the 3D digital core was reconstructed.To evaluate the reconstruction,the auto-correlation function,two-point probability function,porosity,mineral content,and specific surface area were employed.The results show that the SliceGAN can effectively capture the micro-features in the core,and the internal structure of the generated core was consistent with that of the original core.This study provided a new sight for reconstructing cores with complex pore structures and strong heterogeneity and innovatively supports tight carbonate reservoir characterization and evaluation.
基金sponsored by NSFC(Grant No.40574030)CNPC research project(Grant No.06A30102)
文摘In this paper, we obtained three dimensional digital cores using X-ray CT to describe the rock microstructure and applied the open morphology algorithm to simulate oil and formation water distribution in the pore space at different water saturations during the oil-displacing water flood process. The resistivity, formation factor, and resistivity index of rocks were calculated using the finite element method （FEM） and we studied the effect of rock wettability on electrical properties. The numerical simulation results indicate that the simulated formation factor and resistivity index of the water wet rock agrees well with experiments over the whole range of water saturation and extends the traditional resistivity experiment. The rock wettablilty has a large influence on the rock resistivity index. The resistivity and saturation exponent of oil wet rock are obviously larger than three of water wet rock.
基金supported by the Chinese Academy of Sciences Strategic Leading Science and Technology projects(Grant No.XDB10010400)the China Postdoctoral Science Foundation(Grant No.2015M570142)
文摘The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element （FE） methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity （0.75） and kerogen content （ca. ≤3%） over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity （4.5%） of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.
文摘This study sheds light on how pore structure characteristics and varying dynamic pressure conditions influence the permeability of tight sandstone reservoirs,with a particular focus on the Paleozoic reservoirs in the Qingshimao Gas Field.Using CT scans of natural core samples,a three-dimensional digital core was constructed.The maximum ball method was applied to extract a related pore network model,and the pore structure characteristics of the core samples,such as pore radius,throat radius,pore volume,and coordination number,were quantitatively evaluated.The analysis revealed a normally distributed pore radius,suggesting a high degree of reservoir homogeneity and favorable conditions for a connected pore system.However,it was found that the majority of throat radii measured less than 1μm,which severely restricted fluid flow and diminished permeability.Over 50%of the pores measured under 100μm^(3),further constraining fluid movement.Additionally,30%-50%of the pore network was composed of isolated and blind-end pores,which significantly impaired formation connectivity and reduced permeability.Based on this,the lattice Boltzmann method(LBM)was used for pore-scale flow simulation to investigate the influence mechanism of pore structure characteristics and dynamic-static parameters such as displacement pressure difference on the permeability performance of the considered tight sandstone reservoirs for various pressure gradients(0.1,1,and 10 MPa).The simulations revealed a strong relationship between pressure differential and both the number of streamlines and flow path tortuosity.When the pressure differential increased to 1 MPa,30 streamlines were observed,with a tortuosity factor of 1.5,indicating the opening of additional seepage channels and the creation of increasingly winding flow paths.
基金sponsored by NSFC(Grant No.40574030)CNPC Research Project(Grant No.06A30102)
文摘A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.
基金the Natural Science Foundation of China(Nos.51874188,51927807,41941018 and 51704125)the State Key Laboratory for GeoMechanics and Deep Underground Engineering,China University of Mining&Technology(No.SKLGDUEK1717)+1 种基金the Major Scientific and Technological Innovation Project of Shandong Province,China(No.2019SDZY04)the Project of Shandong Province Higher Educational Youth Innovation Science and Technology Program(No.2019KJG013).
文摘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.
基金supported by the National Natural Science Foundation of China Research(Nos.41574122 and 41374124)National Science and Technology major Project(No.2016ZX05006002-004)。
文摘The conventional digital core models are usually small in size and have difficulty in representing the complex structures of heterogeneous rocks;Therefore,the parameters of simulated rock physics are difficult to be referenced.In this study,we propose a feasible simulation method for obtaining multi-scale and multi-component digital cores based on three types of sandstone samples.In the proposed method,the plug and subplug samples are scanned via micro-computed tomography at different resolutions.Furthermore,the images are precisely registered using the proposed hybrid image registration method.In case of high-resolution images,the traditional segmentation method is used to segment the cores into pores and minerals.Subsequently,we established the relations between the gray values and the porosity/mineral content in case of the low-resolution images based on the registered domains and the relation curves were applied to the segmentation of the low-resolution images.The core images constitute the multi-scale and multi-component digital core models after segmentation.Further,the elastic properties of the three samples were simulated at both fine and coarse scales based on the multi-scale and multi-component digital core models,and four component models were considered.The results show that the multi-scale and multi-component digital core models can overcome the representative limits of the conventional digital core models and accurately characterize pores and minerals at different scales.The numerical results of the elastic modulus are more representative at large scales,and considerably reliable results can be obtained by appropriately considering the minerals.
基金We thank Patrick Corbett of Herriot–Watt University for providing the CT scans of the samples.The investigation is financially supported by the National Science&Technology Major Special Project(No.2016ZX05006-002)China Postdoctoral Science Foundation Funded Project(No.2018M632716)+1 种基金Shandong Province Post Doctor Innovative Project Special Fund,Open Project Fund of the National and Local Joint Engineering Research Center of Shale Gas Exploration and Development(No.YiqKTKFGJDFLHGCYJZX444-201901)Chongqing Basic Research and Frontier Exploration Project(No.cstc2018jcyjax0503).
文摘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.
基金financially supported by the National Natural Science Foundation of China(Grant No.41972138)the Technology Major Project of China(Grant No.ZD2019-183007,2016ZX05002-002)。
文摘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.
基金supported by Open Fund (PLN1506) of State Key Laboratory of Oil and Gas Reservoir Geology and ExploitationChinese National Natural Science Foundation (41502287)+2 种基金Chongqing Basic and Frontier Research Projects (CSTC2015JCYJBX0120)Chongqing City Social Undertakings and Livelihood Protection Science and Technology Innovation Special Project (CSTC2017SHMSA120001)Chongqing Land Bureau Science and Technology Planning Project (CQGT-KJ-2017026,CQGTKJ-2015044,CQGT-KJ-2015018, CQGT-KJ-2014040)
文摘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.
基金the National Natural Science Foundation of China(No.41274129)Chuan Qing Drilling Engineering Company's Scientific Research Project:Seismic detection technology and application of complex carbonate reservoir in Sulige Majiagou Formation and the 2018 Central Supporting Local Co-construction Fund(No.80000-18Z0140504)the Construction and Development of Universities in 2019-Joint Support for Geophysics(Double First-Class center,80000-19Z0204)。
文摘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.
基金supported by the National Natural Science Foundation of China (40925013)the Strategic Leading Science and Technology Programme (Class B) of the Chinese Academy of Sciences (Grant No.XDB10010400)the China National Major Science and Technology Project (2011ZX05023-005004)
文摘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.
基金Supported by Project of the National Natural Science Foundation of China (No. 42072323)
文摘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.
基金Project supported by the National Natural Science Foundation of China for the Youth,“Study on the Effect of Shale Bedding Architecture on Features of Ultrasonic Wave and Its Application”(No.41502287)the S&T Program of Chongqing Land and Resources Bureau“Evaluation on Gas Saturation in Shale based on Digital Core Technology”(No.CQGT-KJ-2015018)Sinopec Key Geophysical Laboratory Open Fund Project“Apparent Characteristics of Shale Reservoir and Digital Rock Physics Laboratory Study”(No.33550006-15-FW2099-0015).
文摘Research and development of water lock inhibiting measures is very crucial in verifying the link mechanism between the internal factors of water lock and its extent of damage.Based on conventional water-lock physics experiments,however,only the consequence of macro water lock damage can be investigated,while the microscopic mechanism cannot be studied.In this paper,3D digital cores of low-permeability sandstones were prepared by means of high-resolution micro-CT scan,and their equivalent pore network model was built as well.Virtual“imbibition”experiments controlled by capillary force were carried out by using pore-scale flow simulation.Then the link mechanism between the microscopic internal factors(e.g.wettability,water saturation and pore-throat structure parameters)and the water-lock damage degree was discussed.It is shown that the damage degree of water lock reduces gradually as the wettability transits from water wet to gas wet.Therefore,the water lock damage can be reduced effectively and gas-well productivity can be improved so long as the capillary environment is changed from strong water wettability to weak gas wettability.The more different the initial water saturation is from the irreducible water saturation,the more serious the water lock damage is.The damage degree of water lock is in a negative correlation with the coordinate number,but a positive correlation with the poreethroat ratio.Based on the existing research results,water lock tends to form in the formations composed of mediumsized throats.It is concluded that there is a critical throat radius,at which the water lock is the most serious.
基金supported by the Russian Science Foundation(Grant No.22-11-00273).
文摘The paper presents the results of comprehensive studies of filtration and capacitance properties of highly porous reservoir rocks of the aquifer of an underground gas storage facility.The geomechanical part of the research included studying the dependence of rock permeability on the stress-strain state in the vicinity of the wells,and physical modeling of the implementation of the method of increasing the permeability of the wellbore zone-the method of directional unloading of the reservoir.The digital part of the research included computed tomography(CT)-based computer analysis of the internal structure,pore space characteristics,and filtration properties before and after the tests.According to the results of physical modeling of deformation and filtration processes,it is found that the permeability of rocks before fracture depends on the stress-strain state insignificantly,and this influence is reversible.However,when downhole pressure reaches 7-8 MPa,macrocracks in the rock begin to grow,accompanied by irreversible permeability increase.Porosity,geodesic tortuosity and permeability values were obtained based on digital studies and numerical modeling.A weak degree of transversal anisotropy of the filtration properties of rocks was detected.Based on the analysis of pore size distribution,pressure field and flow velocities,high homogeneity and connectivity of the rock pore space is shown.The absence of pronounced changes in pore space characteristics and pore permeability after non-uniform triaxial loading rocks was shown.On the basis of geometrical analysis of pore space,the reasons for weak permeability anisotropy were identified.The filtration-capacitance properties obtained from the digital analysis showed very good agreement with the results of field and laboratory measurements.The physical modeling has confirmed the efficiency of application of the directional unloading method for the reservoir under study.The necessary parameters of its application were calculated:bottomhole geometry,stage of operation,stresses and pressure drawdown value.
基金the National Key Research and Development Program of China(2022YFE0203400)the National Natural Science Foundation of China(Nos.U23A20595,52034010,52288101)+2 种基金the Qingdao Natural Science Foundation(No.23-2-1-230-zyyd-jch)the Fundamental Research Funds for the Central Universities(No.23CX10004A)the Program for Changjiang Scholars and Innovative Research Team in University(No.IRT_16R69)。
文摘Low-salinity waterflooding,as a promising enhanced oil recovery method,has exhibited exciting results in various experiments conducted at different scales.For carbonate rock,pore-scale understanding of the fluid distribution and remaining oil after low-salinity waterflooding is essential,especially the geometry and topology analysis of oil clusters.We performed the tertiary low-salinity waterflooding and employed X-ray micro-CT to probe the pore-scale displacement mechanism,fluid configuration,oil recovery,and remaining oil distribution.We found that the core becomes less oil-wet after low-salinity waterflooding.Furthermore,we analyzed the oil-rock and oil-brine interfacial areas to further support the wettability alteration.By comparing images after high-salinity waterflooding and low-salinity waterflooding,it is proven that wettability alteration has a significant impact on the behavior of the two-phase flow.Our research demonstrates that low-salinity waterflooding is an effective tertiary enhanced oil recovery technology in carbonate,which changes the wettability of rock and results in less film and singlet oil.
文摘During the development of tight reservoirs,as the formation pressure decreases,the effective stress on the rock matrix increases.This leads to rock deformation and,consequently,changes in the physical parameters of the reservoirs,affecting reservoir productivity.Therefore,it is crucial to identify the influence of stress variation on the pore structure of tight reservoirs at the pore scale.Based on 3D grayscale images of representative samples under different confining pressures using micrometer computed tomography,this study obtained corresponding 3D digital cores through binary segmentation using the maximum class spacing algorithm.Based on digital cores with the same physical size under different confining pressures,the binary data volume was subtracted using a Boolean algorithm to obtain the pore variation space under different confining pressures.The pore variation and permeability damage rate under different confining pressures were then determined.Pore network models of digital cores under different confining pressures were extracted using the maximal ball algorithm,and the corresponding pore-throat radius distribution,coordination number distribution,length/diameter ratio,and shape factor distribution were calculated.The results show that as the confining pressure increases,the pore-throat distribution curve shifts to the direction of smaller pore-throats,the average coordination number decreases,and the throat length/diameter ratio distribution curve and shape factor distribution curve shift to the direction of larger pore-throats.This is because compaction reduces pores and throats,decreases the overall pore-throat connectivity,narrows the throats,and makes the crosssections of the pore-throats more circular.This research serves as a basic platform for efficient development under stress sensitivity in tight reservoirs,and its findings have important academic significance and practical application value.
基金the National Key Research and Development Program of China(No.2021YFC3100800)the National Natural Science Foundation of China(Nos.42407235 and 42271026)the Project of Sanya Yazhou Bay Science and Technology City(No.SCKJ-JYRC-2023-54).
文摘Coral reef limestone(CRL)constitutes a distinctive marine carbonate formation with complex mechanical properties.This study investigates the multiscale damage and fracture mechanisms of CRL through integrated experimental testing,digital core technology,and theoretical modelling.Two CRL types with contrasting mesostructures were characterized across three scales.Macroscopically,CRL-I and CRL-II exhibited mean compressive strengths of 8.46 and 5.17 MPa,respectively.Mesoscopically,CRL-I featured small-scale highly interconnected pores,whilst CRL-II developed larger stratified pores with diminished connectivity.Microscopically,both CRL matrices demonstrated remarkable similarity in mineral composition and mechanical properties.A novel voxel average-based digital core scaling methodology was developed to facilitate numerical simulation of cross-scale damage processes,revealing network-progressive failure in CRL-I versus directional-brittle failure in CRL-II.Furthermore,a damage statistical constitutive model based on digital core technology and mesoscopic homogenisation theory established quantitative relationships between microelement strength distribution and macroscopic mechanical behavior.These findings illuminate the fundamental mechanisms through which mesoscopic structure governs the macroscopic mechanical properties of CRL.
基金Supported by the PetroChina"Fourteenth Five-Year Plan"Prospective Basic Technology Research Project(2021DJ4002)PetroChina Major Oil and Gas Project(2021ZZ10-01).
文摘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.
