Theoretical Analysis of Simulating the Locked-In Stress in Rock Pore by Thermal Expansion of Hard Rubber

Rocks are composed of mineral particles and micropores between mineral which has a great influence on the mechanical properties of rocks. In this paper, based on the theory of locked-in stress developed by academician Chen Zongji, the locked-in stress problem in underground rock is simulated by the thermal expansion of hard rubber particles. The pore inclusion in rock is assumed to be uniformly distributed spherical cavities. Using the thermal stress theory, the stress of rock with a spherical pore inclusion is equivalent to the thermal stress generated by the spherical hard rubber inclusion. The elastic theory formula of the temperature increment and the equivalent pore pressure of the spherical hard rubber inclusion is derived. The numerical simulation of the rock mass model with a spherical hard rubber inclusion is carried out and compared to the theoretical calculation results; the results show that they are consistent. The method proposed by this paper for simulating stress distribution in rock by thermal stress is reasonable and feasible; it has a positive meaning for further study of mechanic phenomenon of rock with micropore inclusion.

A novel box-counting method for quantitative fractal analysis of threedimensional pore characteristics in sandstone

Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.

Computation of fractal dimension of rock pores based on gray CT images

The characterization of pore structure in rocks is relevant in determining their various mechanical behaviors. Digital image processing methods integrated with fractal theory were applied to analyze images of rock slices obtained from industry CT, elucidating the characteristics of rock pore structure and the relationship between porosity and fractal dimensions. The gray values of pixels in CT images of rocks provide comprehensive results with respect to the attenuation coefficients of various materials in corresponding rock elements, and these values also reflect the effect of rock porosity at various scales. A segmentation threshold can be determined by inverse analysis based on the pore ratios that are measured experimentally, and subsequently binary images of rock pores can be obtained to study their topological structures. The fractal dimension of rock pore structure increases with an increase in rock pore ratio, and fractal dimensions might differ even if pore ratios are the same. The more complex the structure of a rock, the larger the fractal dimension becomes. The experimental studies have validated that fractal dimension calculated directly from gray CT images of rocks can give an effective complementary parameter to use alongside pore ratios and they can suitably represent the fractal characteristics of rock pores.

Static dissolution-induced 3D pore network modification and its impact on critical pore attributes of carbonate rocks

To determine the effect of dissolution on pore network development in carbonate rocks, dissolution experiments, X-Ray microtomography, and thin section analysis were conducted on argillaceous limestone and grain limestone samples at different temperatures and constant pH, HCl concentration. The relationship between Ca^(2+) concentration and time was revealed through the experiments; pore size distribution before and after dissolution indicate that there is no correlation between the temperature and pore size variation, but pore size variation in grain limestone is more significant, indicating that the variation is mainly controlled by the heterogeneity of the rock itself(initial porosity and permeability) and the abundance of unstable minerals(related to crystal shape, size and mineral type). At different temperatures, the two kinds of carbonate rocks had very small variation in pore throat radius from 0.003 mm to 0.040 mm, which is 1.3 to 3.5 times more, 1.7 on average of the original pore throat radius. Their pore throat length varied from 0.05 mm to 0.35 mm. The minor changes in the pore throat radius, length and connectivity brought big changes to permeability of up to 1 000×10^(-3) μm^2.

Comparison of the Pore Structure of Ultralow-Permeability Reservoirs Between the East and West Subsags of the Lishui Sag Using Constant-Rate Mercury Injection

In this study, the differences in reservoir parameters, such as pore radius, throat radius, and pore-throat ratio, between the east and west subsags of the Lishui Sag are analyzed by using data obtained from a constant-rate mercury injection experiment. Furthermore, the quality of the reservoirs in the two subsags is systematically evaluated. Results show that the throat radius of the Lishui west subsag is larger than that of the east subsag, and this parameter has a positive correlation with reservoir quality. However, the pore-throat ratio of the east subsag is larger than that of the west subsag, which has an inverse relationship with reservoir quality. The main reasons for this reservoir difference can be attributed to sedimentation and diagenesis. The sedimentary facies types of the Lishui east subsag are the fan delta, shore lake, shallow lake, and shore shallow lake;their sandstone composition maturity is low;the clay mineral content is high;and the rock has undergone strong diagenesis. Therefore, the physical conditions of the reservoir are poor. However, the sandstones in the Lishui west subsag have weak cementation and compaction, mainly with an intergranular pore structure type, which leads to good connectivity between pores. Therefore, the storage performance and seepage capacity of the Lishui west subsag are better than those of the east subsag;the west subsag is the main area of oil and gas accumulation, as confirmed in the process of exploration and development.

Pore-scale investigation of residual oil displacement in surfactant–polymer flooding using nuclear magnetic resonance experiments

Research on the Gangxi III area in the Dagang Oilfield shows that there was still a significant amount of oil remaining in oil reservoirs after many years of polymer flooding.This is a potential target for enhanced oil recovery（EOR）.Surfactant–polymer（SP） flooding is an effective chemical EOR method for mobilizing residual oil and improving displacement efficiency macroscopically,but the microscopic oil displacement efficiency in pores of different sizes is unclear.Nuclear magnetic resonance（NMR） is an efficient method for quantifying oil saturation in the rock matrix and analyzing pore structures.In this paper,the threshold values of different pore sizes were established from the relationship between mercury injection curves and NMR T2 spectrums.The distribution and migration of residual oil in different flooding processes was evaluated by quantitatively analyzing the change of the relaxation time.The oil displaced from pores of different sizes after the water flood,polymer flood,and the SP flood was calculated,respectively.Experimental results indicate that（1） the residual oil in medium pores contributed the most to the incremental oil recovery for the SP flood,ranging from 40 % to 49 %,and small pores usually contributed ／30 %;（2） the residual oil after the SP flood was mainly distributed in small and medium pores;the residual oil in medium pores accounted for 47.3 %–54.7 %,while that trapped in small pores was 25.7 %–42.5 %.The residual oil in small and medium pores was the main target for EOR after the SP flood in oilfields.

A Study on Organic Inclusions in Clastic Reservoir Rocks and Their Application to the Assessment of Oil and Gas Accumulations

Terrigenous clastic reservoir rocks are widespread in China, and nearly all the industrial oil and gas accumulations in eastern China occur in the clastic rocks. The study shows that organic inclusions are mostly distributed in the secondary fissures and pores which were formed in the process of oil-rock interaction, rather than in the cements or secondary enlargements.The organic inclusions are dominantly organic gas-rich or are composed of pure hydrocarbons.Homogenization temperatures range mainly from 120℃ to 130℃, which shows a relatively high maturity of organic matter. Vertical and horizontal temperature changes provide the grounds for the investigation of hasin evolution and thermal fluid-kinetics model. Fluorescence spectral characteristics of the organic inclusions indicate that oils and gases in the area studied probably have experienced two-stage or two-time migration. Micro-fluorescence research is one of the effective approaches to oil/source correlation and oil migration-stage determination. The abundance and occurrence of organic inclusions is one of the indicators of oil and gas abundance and accumulation in rock layers. With the help of other information, organic inclusions can provide the basis for the prospective assessment of oil and gas in clastic reservoir rocks.

Thermal effects on clay rocks for deep disposal of high-level radioactive waste

Thermal effects on the Callovo-Oxfordian and Opalinus clay rocks for hosting high-level radioactive waste were comprehensively investigated with laboratory and in situ experiments under repository relevant conditions:(1) stresses covering the range from the initial lithostatic state to redistributed levels after excavation,(2) hydraulic drained and undrained boundaries, and(3) heating from ambient temperature up to 90℃-120℃ and a subsequent cooling phase. The laboratory experiments were performed on normal-sized and large hollow cylindrical samples in various respects of thermal expansion and contraction, thermally-induced pore water pressure, temperature influences on deformation and strength, thermal impacts on swelling, fracture sealing and permeability. The laboratory results obtained from the samples are consistent with the in situ observations during heating experiments in the underground research laboratories at Bure and Mont-Terri. Even though the claystones showed significant responses to thermal loading, no negative effects on their favorable barrier properties were observed.

Characteristics of parametamorphic rock reservoirs in Pingxi area, Qaidam Basin, NW China

Based on core, thin-section, scanning electron microscopy(SEM) and well logging data, the characteristics of the parametamorphic rock reservoirs in the Pingxi area were analyzed by means of whole rock X-ray diffraction and micron CT scanning. The parametamorphic rock reservoirs mainly had three types of rocks: slate, crystalline limestone and calc-schist; the original rocks were Ordovician-Silurian marine clastic and carbonate rocks. The three types of parametamorphic rock reservoirs developed three types and six sub-types of reservoir space. The first type of reservoir space was fractures, including structural, weathered and dissolution fractures; the second type was dissolved porosities, including dissolved pores and caves; the third type was nano-sized intercrystalline porosities. The three types of parametamorphic rock reservoirs were different widely in the quantity, volume and radius of pore-throats, and were strongly affected by the type and development degree of fractures. The parametamorphic rock reservoirs were formed by metamorphism, weathering, structural fragmentation and dissolution. Metamorphism reformed the parametamorphic rock reservoirs significantly, breaking the traditional constraint of finding weathering crust at top. The parametamorphic rock reservoirs experienced five formation stages, and their distribution was controlled by rock type, metamorphic degree, ancient geomorphology, and weathering intensity.

Reservoir Rock Facies──Case Presentation of Sandstone Type

The reservoir rock facies is presented by its principal indieators to show the aspects of sedimentary facies,diagenetic change and oil-gas indication. These indicators used by the author are 1) the petrologic──of which the sedimentary facies and vitrinite reflectance are included, 2) the physical ── pore/throat di-ameter ratio and coordination number of throat connecting pore, and 3) the geochemical──photochemical parameters of individual organic inclusion. Based on the above mentioned quantitative indicators,the defined reservoir rock facies may not only be used for itself evaluation, but also may be put in facies column or facies-palaeogeographic map to predict or to trace oil-gas reservoir. Microscope photometry,micro-FT-IR and Laser Raman methods were used for studying all the aforesaid parameters by means of thinned polished sections from core or cemented cuttings, except the 3-D others, from parallel and vertical to bedding or some duplicate core samples, that the rose fluorescein preparation must be soaked in for the convenience of studying pore throat structure.

Integrated Rock Typing of the Grainstone Facies in a Sequence Framework：a Case from the Jurassic Arab Formation in the Persian Gulf

The late Jurassic Arab Formation, a significant carbonate-evaporite reservoir rock in the Persian Gulf, is characterized by frequent grainstone facies. For rock type identification and reservoir characterization, core description, petrographic studies and pore system evaluation are integrated for Balal oil field in the Persian Gulf. The grainstone facies are developed into three shoal subenvironments on a carbonate ramp platform： leeward, central and seaward. Compaction, dissolution, cementation, anhydrite mineralization and dolomitization are the main diagenetic processes affecting the depositional pore system. Considering depositional and diagenetic features and pore types, the grainstones are classified into six rock types （RT 1 to RT6）. Rock types 1, 2 and 5 have large pore throat sizes with inter- granular and touching vug pore types. In rock type 3, moldic pores lead to high porosity and low permeability. Rock types 4 and 6 are cemented by anhydrite, calcite and dolomite. Generally, RTs 1, 2, 3 and 5 are related to late Transgressive systems tract （TST） and early Highstand systems tracts （HST） and show fair to good reservoir quality. In contrast, RTs 4 and 6 of late HST system tract show lower poroperm values, due to evaporite mineralization. Characterization of the grainstone facies provides a comprehensive understanding of the reservoir zones of the Arab Formation.

基于电成像的孔洞型碳酸盐岩储层连通性评价

碳酸盐岩储层非均质性强,孔隙结构复杂,相较于常规碎屑岩储层,其孔隙度、渗透率关系的规律及特征不明显。碳酸盐岩储层的连通性与孔洞特征具有强相关性,且连通性是影响储层渗透率的重要因素,因此,对碳酸盐岩孔洞连通性的评价是孔洞型碳酸盐岩储层评价的关键问题。以伊拉克东南部B油田FQ-1井M组MB21亚油组孔洞型碳酸盐岩储层为例,基于地层微电阻率扫描成像(FMI)测井的高分辨率动态图像,利用最大类间方差算法(OTSU)进行二次迭代分割,识别电成像图像孔洞;根据8连通域标记法标记连通孔洞;提取每个连通孔洞的特征孔径、内部最长通道距离、内部最长通道与水平方向夹角、平均通道宽度这4个几何参数,最终构建孔洞型碳酸盐岩储层连通指数。连通指数在伊拉克东南部B油田FQ-1井M组MB21亚油组孔洞型碳酸盐岩储层的应用表明,利用此方法建立的连通指数与FQ-1井碳酸盐岩储层岩心渗透率相关系数达0.8437,并且连通指数曲线的趋势符合岩心渗透率在储层不同深度的分布情况,证明了连通指数可以指示储层孔隙结构较好、溶蚀孔洞发育的深度段,并可以划分储层中连通性较好的深度区间,该参数对于寻找碳酸盐岩储层渗透层具有重要意义。

流固耦合作用下深部岩石动态力学响应研究进展

深部岩石处于高地应力、高渗透压和强动态扰动的复杂地质环境之中,3者作用下岩石体更加容易发生损伤破裂,诱发突涌水、渗漏、井喷等工程地质灾害,因而探究流固耦合作用下岩石的动态力学响应是开展岩石工程建设的前提之一。近年来,国内外众多学者在考虑水和不同应力状态下的岩石动态力学实验研究方面取得了显著的成果。为给工程建设提供更加全面的指导并为后续研究奠定基础,从实验装置、测试结果以及围压与水的作用机理层面,对上述工作进行了回顾与总结。首先介绍了分离式霍普金森压杆测试装置的基本原理,以及用于模拟深部岩石赋存环境所进行的装置改进,包括围压分离式霍普金森压杆实验系统和孔压(渗透压)耦合的分离式霍普金森压杆实验系统,简要分析了各类装置在研究流固耦合作用下岩石动力学问题时的优势和不足。其次,总结了考虑不同应力状态(单向加载、三向围压加载)的流固耦合作用下岩石的动态力学响应特性。详细介绍了固定预设孔压、渗透压耦合作用下深部岩石的动态力学响应及其随孔隙水压、渗透压变化的规律。随后,概述了围压对岩石动力学性质的影响机理,分析了不同围压条件下的影响规律;总结了水对岩石动态力学性质的强化、弱化微观机制和定量表达。最后,对流固耦合作用下深部岩石的动态力学响应进行了概括总结,并对未来实验研究工作和深部赋存条件下岩石动态力学的研究方向进行了展望。

注热化学流体吞吐解除岩石储层反凝析伤害评价

为长期有效解除凝析气藏反凝析带来的储层伤害,提出了一种注热化学流体吞吐解除凝析油堵塞孔喉的新方法。针对致密砂岩、碳酸盐岩和页岩3类岩石,开展了注热化学流体吞吐岩心驱油实验,评价了热化学流体解除凝析油堵塞的可行性,并基于核磁共振技术分析了吞吐前后3类岩石微观孔隙结构的变化特征,明确了热化学流体解除反凝析伤害机理。实验结果表明,热化学流体(NH_(4)Cl和NaNO_(2))在乙酸溶液催化下能够迅速释放大量气体和蒸汽,提高岩心压力和温度,降低原油黏度,使凝析油由液态转变为气态,减小渗流阻力。致密砂岩、碳酸盐岩和页岩岩心注热化学流体吞吐累积凝析油采出程度分别为65.7%,73.9%和46.3%。其中,致密砂岩和碳酸盐岩岩心仅需2轮吞吐即可有效清除55.5%和67.6%的凝析油堵塞。而页岩岩心需要延长焖井时间及增加吞吐次数方能提高凝析油的采出程度。热化学流体吞吐后,致密砂岩、碳酸盐岩和页岩岩心平均孔径分别由0.37,1.04和0.0021μm增大至0.84,2.04和0.0058μm。热化学流体吞吐效果是注CH_(4)吞吐和注甲醇吞吐效果的1.85和1.32倍。

基于储层砂岩孔隙结构的地震岩石物理模型

在实验室内对新疆车排子地区春光油田沙湾组储层砂岩样品进行了系统地震岩石物理实验,在确定储层砂岩中不同孔隙类型组成条件下,基于非接触模型综合考虑储层砂岩孔隙度以及孔隙形状对于干燥岩石骨架弹性性质的影响,并利用实验结果对模型参数进行标定进而建立表征目标储层岩石骨架地震弹性性质的岩石物理模型,模型结果与实验结果吻合良好。依据实验及模型结果,利用喷射流作用模型替代传统Gassmann方程进行流体替换,同样能够较好地对超声频率实验结果进行解释。

深层超深层油气藏高应力下数字岩心构建方法

深层超深层油气藏由于埋藏深,其地应力达200 MPa,会显著改变储层岩石孔隙的微观结构。数字岩心是孔隙尺度数值模拟的重要载体,但是现有数字岩心重构方法是基于常温常压下岩心的扫描图像重构,不能反映高应力下的孔隙结构。为此,提出了一种基于离散元法考虑高应力影响的数字岩心重构方法。首先,采用分水岭算法分割CT图像,利用球面谐波分析方法建立轮廓数据库,并在PFC^(3D)中建立Clump(团簇)模板库;然后,根据孔隙度和粒径分布使用模板库中的Clump建立离散元模型,并用两点相关和线性路径相关函数曲线评价模型的准确性;随后,标定颗粒间微观力学参数,并加载应力模拟得到不同应力下的数字岩心;最后,分析了不同应力下数字岩心的孔隙几何拓扑结构,计算孔隙度和渗透率。以Bentheim砂岩为例,构建了其不同应力下的数字岩心,研究结果表明,应力增大,导致孔隙和喉道半径缩小、喉道伸长、连通性变差、孔隙度和渗透率减小。研究结果为深层超深层油气藏孔隙尺度模拟提供了技术途径。

基于数字岩心的碳酸盐岩复杂孔隙特征研究——以普光气田飞仙关组储层为例

川东北普光气田飞仙关组为一套深层海相碳酸盐岩储层,储集空间类型、组合方式多样,孔隙结构复杂,储层非均质性强,利用常规测井资料难以有效表征其储层的孔隙结构及渗透力。文中应用数字岩心技术精细刻画孔隙结构,即采用全直径岩心CT扫描、微米级X-ray CT扫描、大面积扫描电镜拼图成像技术等,建立纳米—厘米级别的多维度多尺度的岩心数字化表征。对大数据图像特征参数提取和分析,建立了不同类型储层的孔渗关系,将整个储层类型分为三区一带:微裂缝区、粒(晶)间孔主控区、铸模孔区、粒间孔加铸模孔过渡带。研究明确了粒(晶)间孔结构是渗透率的主控因素,更好地明确了该类储层孔隙结构对渗流能力的影响,揭示了礁滩相碳酸盐岩地层的微观渗流特性。

基于岩石物理模型的碳酸盐岩物性参数替换方法

碳酸盐岩的物性参数在油气勘探和开发的过程中起着重要的作用,且碳酸盐岩为具有多重孔隙类型的岩石,复杂的孔隙类型使得孔隙度与弹性参数之间的关系非常离散。本文基于岩石物理模型,提出一种碳酸盐岩物性参数替换的方法,首先对碳酸盐岩储层进行岩石物理建模,对模型中的等效孔隙纵横比进行反演,在进行替换时,保持其他参数不变,只改变孔隙度、方解石含量、含水饱和度和孔隙形状的体积分数中的一项。再结合AVO理论,进行正演模拟,正演模拟揭示孔隙度与孔隙形状的变化对地震响应的影响更为强烈,远远大于方解石含量、含水饱和度变化的影响,方解石含量和含水饱和度的变化对于地震响应的影响较微弱。实际资料应用表明,文章提出的碳酸盐岩物性参数的替换方法可以有效地分析物性参数及孔隙形状变化的影响,表征岩石的物理性质并判断岩石的孔隙类型。

基于三维数字岩心的油页岩裂解孔隙重构及表征研究

油页岩孔隙结构作为油气渗流的主要通道,高温热解过程中将产生复杂变化,然而其变化规律不明严重制约储集层分类评价和高效开发。因此,开展不同温度及不同有机质含量条件下油页岩热解对孔隙结构影响的研究对于油页岩的有效开发具有重要意义。针对热解前后油页岩孔缝结构变化精细表征难的问题,本文通过编程方法提出了一种基于三维数字岩心重构的裂缝及孔隙结构表征新方法,定量分析了热解过程中样品裂缝宽度、裂缝方向、分形维数及孔隙度等参数的变化规律。研究结果表明:(1)本文所提出的裂缝性数字岩心精细化处理方法结合等位精细化中面提取算法,在准确计算岩石物性参数的基础上可实现对缝宽及裂缝方向分布的精确计算,可达到对裂缝孔隙精细表征的目的;(2)随着热解温度的升高,油页岩样品孔隙度由18.3%升至20.9%,且裂缝宽度不断增大,油页岩样品物性得到持续改善,但当热解进行到一定程度后油页岩孔隙结构变化趋势放缓;(3)相同温度条件下随着有机质含量的升高,热解后油页岩样品孔隙度由12.1%升至29.3%,但4种样品分形维数未随着有机质含量升高呈现出单调递增趋势,表明热解裂缝复杂性除受有机质含量影响外,也受到样品中有机质初始分布的影响。本文可为油页岩热解后孔隙结构精细表征提供新思路。

高功率激光水平辐照花岗岩致裂成孔特性及成孔能效评价

针对当前硬岩隧道钻爆法施工钻孔速度慢、钻头损耗大的工程难题,使用高功率激光水平辐照硬岩快速成孔替代传统PDC钻头钻孔,开展高功率激光水平辐照花岗岩成孔的可行性试验,分析高功率激光水平辐照花岗岩的致裂成孔特性;通过单轴压缩试验、扫描电子显微镜(SEM-EDS)和X射线衍射仪(XRD)衍射试验分析激光作用前后岩石的力学性能和微观物相演变;综合对比激光射孔与PDC钻头钻孔的破岩速度和比能。研究结果表明:岩石表面在激光持续辐照后形成“类爆炸坑”,产生的熔融物在重力作用下堆积在孔底,吸收了孔底处大部分的激光能量,抑制了孔底岩石的热裂反应;岩石的纵向破坏受重力影响导致上下孔壁熔融物分布不均匀,随着激光持续辐照,岩孔中心将不断向下偏离光斑入射点;经功率为4.5 kW的激光水平辐照10 s后,岩样的单轴抗压强度从63.31 MPa降低至42.36 MPa,其原因是钠长石分解引起岩样中黑云母的质量分数增加;功率为4.5 kW激光射孔的破岩速度较PDC钻头钻孔的破岩速度提高95.03%,比能增加15.3%,证明了高功率激光水平辐照硬岩成孔的有效性。