三维微观孔隙结构与岩石变形机制及输运特性

正岩石是一种有着复杂内部结构的非均匀材料。在较小的尺度上,岩石含有各种不同的矿物颗粒、胶结物、孔隙和微裂纹等微观特征。孔隙是岩石微观结构的重要组分之一,是岩石微观非均匀性的重要体现。孔隙对岩石的力学性质、输运性质和其他岩石物理性质有重要的影响。因此:应该在岩石力学和岩石物理学中描述孔隙的几何特性,如孔隙大小、形状、弯曲     

龙马溪组页岩微观结构、地震岩石物理特征与建模

龙马溪组页岩是目前国内页岩气勘探的主要层位之一.由于岩石物理实验结果具有区域性,龙马溪组页岩的岩石特征与其地震弹性性质的响应规律需要开展相关的实验和理论研究工作予以明确.本研究基于系统的微观结构观察(扫描电镜和CT成像技术)和岩石物理实验来分析龙马溪组页岩样品地震弹性性质的变化规律,并依据微观结构特征建立相应的地震岩石物理表征模型.研究结果表明,石英含量对龙马溪组页岩的孔隙度以及有机碳(TOC)含量具有一定的控制作用,TOC和黄铁矿主要赋存于孔隙中;岩石骨架组成亦受控于石英或粘土含量,在石英含量大于40%(对应粘土含量小于30%)时,以石英、粘土共同作为岩石骨架,而粘土含量大于30%时,则以粘土作为岩石的骨架.因此,岩石骨架组成矿物、TOC含量、孔隙度共同制约龙马溪组页岩的地震弹性性质,富有机质储层岩石通常表现出低泊松比、低阻抗和低杨氏模量的特征,但由于支撑矿物的转换,某些富有机质页岩亦可表现为高阻抗特征.粘土矿物的定向排列仍然是造成页岩样品表现出各向异性的主要原因,各向异性参数与粘土含量具有指数关系.基于龙马溪组页岩的岩性特征及微观结构特征,可以利用自洽模型(SCA)、微分等效模量模型(DEM)和Backus平均模型的有效组合较为准确地建立龙马溪组页岩的地震岩石物理模型,实验结果和测井数据验证了模型的准确性.研究结果可为龙马溪组页岩气储层的测井解释和地震"甜点"预测提供依据.     

地壳岩石矿物电导率实验研究进展

矿物和岩石电导率的实验研究一直都是国内外专家学者所关注的重点内容,目前电导率的实验测量主要采用交流阻抗谱的方法.从实验结果来看,地壳矿物和岩石电导率的影响因素主要分为外部环境和内部性质两部分,包括矿物和岩石所受到的温度和压力大小、内部赋存的流体和熔体、岩石颗粒边缘的碳膜、岩石的颗粒大小和面理方向以及矿物晶格赋存的结构水与晶格方向等.结合地球物理探测结果,前人对下地壳高导层成因提出了各自不同的看法,同时电导率的实验数据对于解释地下岩层结构、动力学特征也带来了一定的帮助.通过总结可以发现对于断层带的电导率性质目前仍然知之甚少,在今后的实验研究中则需要重点关注.     

温度对岩石导电性的影响及其与地震的关系

组成地壳表部岩石的导电性在自然埋藏状态下主要取决于岩石的结构、构造和矿物成分以及岩石所处环境的温度、压力、孔隙溶液的含量、浓度及其化学成分。不同含湿量的各种岩石在外力作用下视电阻率的变化特征国内外已积累了较丰富的实验资料。但是,温度对岩     

鄂尔多斯盆地上三叠统延长组致密砂岩储层绿泥石对润湿性的影响

润湿性是储层的重要性质,对提高采收率和油气运移研究至关重要.岩石矿物组成影响储层的润湿性,但对其具体作用仍没有明确认识.文章以鄂尔多斯盆地延长组储层为例,开展了储层岩石学特征和润湿性特征的综合研究,分析了绿泥石对砂岩储层润湿性的影响作用.结果表明,延长组储层内绿泥石非常发育,主要包括碎屑状绿泥石、绿泥石膜和黑云母转化的绿泥石;其中绿泥石膜通常覆盖在孔隙周围,构成孔隙表面的主要矿物.储层受到原油的充注及改造作用后,绿泥石在孔隙尺度具有水润湿和油润湿两类.绿泥石是组成储层油润湿孔壁的主要矿物,对岩石油润湿性的形成起重要作用.在岩芯尺度,受到岩石非均质性和含油性差异等影响,孔隙周围绿泥石的含量与Amott-Harvey指数相关性较差.     

灰岩微观孔隙的三维定量描述及其对变形的影响

岩石的力学行为是由其微观的岩石物理性质和微观结构所最终决定的.因此在岩石物理学中描述孔隙结构及其几何性质如孔隙大小、形状、弯曲度和连通性有着极其重要的意义.研究孔隙结构的传统方法是光学显微镜和扫描电镜.微米CT (microCT)作为一种新的三维成像技术为研究岩石孔隙结构的复杂性提供了新的视角,这对理解岩石的孔隙结构及其对岩石的物理性质的影响,以及岩石破坏的演化及其与破裂微观力学的关系有极大的帮助.     

四川盆地侏罗系致密砂岩弹性特征及岩石物理建模

近年来围绕四川盆地侏罗系陆相致密砂岩已取得了勘探突破,其中川中—川西过渡带具备形成大气田的地质条件,但对该套致密砂岩弹性性质变化规律的研究还较少,致使利用地震方法进行"甜点"储层预测的精度不高.本文利用四川盆地侏罗系沙溪庙组32块样品开展了系统的声学测量,在此基础上,分析了样品弹性性质的变化规律.结合X射线衍射矿物组分分析、扫描电镜、铸体薄片和岩石薄片特征确定了不同成岩作用对岩石储集性能的影响.研究结果表明,研究区致密砂岩储层表现为孔隙型储层,受差异性成岩作用影响,黏土含量、钙质含量和硅质含量的差异以及它们分布特征之间的差异对岩石弹性性质造成了很大的影响.在研究区对岩石物性及弹性性质有明显影响的成岩作用包括早期的钙质胶结作用、压实作用和溶蚀作用,因此针对不同时期的成岩作用对岩石弹性及物性的影响,利用接触-胶结模型、微分等效模量模型和临界孔隙度校正的Hashin-Shtrikman上限模型建立了研究区致密砂岩的岩石物理模型.     

个旧芦塘坝黄铜矿P-A分形模型研究

利用分形方法对微观尺度下矿物微观组构及变化进行研究是目前的热点.本文以个旧芦塘坝矿段不同高度的硫化物矿石中广泛发育的黄铜矿为研究对象.建立起微观条件下的黄铜矿矿物的面积-周长(P-A)分形模型,并研究黄铜矿矿物的大小、不均匀性及空间分布特征.探讨黄铜矿这一标志性矿物在不同高度中的变化情况.结果显示随着深度从1758 m...     

缅甸掸邦地块二叠纪-三叠纪灰岩的岩石磁学初步研究结果

对缅甸东部掸邦地块二叠纪-三叠纪灰岩样品进行了饱和等温剩磁获得曲线及反向场退磁曲线、三轴等温剩磁的系统热退磁、磁化率-温度曲线等岩石磁学测试分析,试图通过岩石磁学方法识别缅甸掸邦地块二叠纪和三叠纪灰岩样品中的磁性矿物种类、含量等信息,为研究样品中主要磁性矿物所携带的特征剩磁分量的磁化年龄等性质提高可靠的基础,从而为进一步研究掸邦地块及滇缅泰地块晚古生代-早中生代的古地理位置等提供有力的证据.结果表明缅甸掸邦地块二叠纪和三叠纪灰岩样品中的主要磁性矿物为磁铁矿,或磁铁矿和针铁矿的组合,少部分三叠纪灰岩样品显示有赤铁矿、磁赤铁矿和针铁矿的组合.     

龙马溪页岩弹性各向异性与矿物分布之间的关系探讨

岩石弹性各向异性特征是普遍存在的,但导致岩石弹性各向异性的原因复杂且仍然存在一定争议.本研究以龙马溪页岩为例,试图建立页岩弹性各向异性和矿物分布之间关系.研究使用无损超声波探测获取岩石弹性各向异性参数,并使用背散射技术获取岩石矿物分布特征.研究通过引进变异系数来描述矿物或孔裂隙在不同方向的差异性,并通过2个正交方向的变异系数建立微观异质性指标,用于表征微观尺度上矿物或孔裂隙在不同方向的差异程度.微观背散射图像显示,龙马溪组页岩主要由石英和伊利石构成,且这两种矿物和孔裂隙在X和Y面上都有定向排列特征.相应地,它们的变异系数在X面和Y面上均表现出随角度增大而降低的特征;而在Z面,石英、伊利石和孔裂隙变异系数变化不明显,这与Z面上矿物和孔裂隙无明显方向性的特点一致.无损超声波探测结果显示,波速在X和Y面上随角度增加而减小,这与主要矿物和孔裂隙变异系数变化趋势相同;而在Z面,波速变化不大,与主要矿物和孔裂隙变异系数变化不明显的特征一致.以上观测结果说明,宏观波速与矿物的微观变异系数明显相关,暗示岩石弹性各向异性与矿物分布直接相关.     

Multifractality in Seismicity Spatial Distributions: Significance and Possible Precursory Applications as Found for Two Cases in Different Tectonic Environments

We explore fractal properties of two observed seismicity distributions prior to the 2003 M _w 7.4 Colima, Mexico and 1992 M _w 7.3 Landers, USA earthquakes, together with several mathematical fractal distributions and two non-fractal ones, in order to estimate minimum reliable sample sizes, determine whether fractality for observed seismicity is essentially different from random uniform distributions, and explore the possibility of extracting premonitory information from fractal characteristics of seismicity before large earthquakes. Sample sizes above 800 events for whole catalogs appear to be sufficient to maintain ordered multifractality and to yield dimension estimates that vary smoothly and reliably. Fractal estimates appear to be best for whole catalogs that include aftershocks. The fractal characteristics of spatial distributions of seismicity are essentially different from those of the uniform random distribution, which is the null hypothesis of a non-fractal distribution with minimum information. The fractal dimensions and afractality measures of seismicity distributions change with time and show distinctive behaviors associated with foreshocks and main events, although these behaviors are different for each example. Results suggest the possibility of a priori identification of foreshocks to large earthquakes. A combination of fractal dimension and afractality measures over time may be helpful in large earthquake premonitory studies.     

Fractal dimension of fault systems in Japan: Fractal structure in rock fracture geometry at various scales

Based on fault maps, whether or not the fracture geometry of rocks is self-similar, was examined by using a box-counting algorithm. The statistical self-similarity (fractal structure) of the fault fracture systems holds well at the scale of about 2 to 20 km. The fractal dimension in Japan varied from 1.05 to 1.60. The fractal dimension is about 1.5–1.6 at the central part of the Japan Arc, and decreases with distance from the center. At a smaller scale, the fractal structure also holds well in the rock fracture geometry. The fractal dimension of the North Izu Peninsula fault system (branching faults) is 1.49 at the scale of 0.625 to 10 km, the fractal dimension of rock fracture geometry at the scale order of 10^–1 to 10^–2 meters is about 1.49–1.61. The upper limit of the fractal dimension of rock fracture geometry is about 1.6, judging from the estimation of fractal dimension on actual fracture geometry of rocks. This value may impose a restraint on modeling of faulting and the fracture process of rocks.     

Pore structure of selected Chinese coals with heating and pressurization treatments

Pore structure of Chinese coals with heating and pressurization treatments was studied using small angle X-ray scattering(SAXS),N2adsorption/desorption isotherms and scanning electron microscope(SEM).SAXS was performed for some samples after heat treatment at seven elevated temperatures from 25 to 250°C at 0 MPa and for other samples with hydrostatic pressure treatment at 0,5,10,15 and 20 MPa at the room temperature.The results show that N2adsorption isotherm together with SAXS could be a comprehensive method to evaluate the pore shape and the pore size distribution:the pore shapes are generally spherical for low rank coal and they are mainly ellipsoidal for high rank coal.All these measurements were then interpreted using the fractal theory to reveal relationship between surface fractals and coal rank,and the evolution of surface fractals under heating and pressurization treatments.The results show that surface fractal dimension(Ds)changes with different treating temperature and pressure and maximum vitrinite reflectance(Ro,m).Especially in the bituminous stage,Ds shows an increasing trend with Ro,m under varied temperatures.Moreover,Ds shows an increasing trend with increasing temperature before 200°C,and a decreasing trend after 200°C.Furthermore,the results show that Ds has a more complex relationship with Ro,m under varied treating temperature than that under varied treating pressure.     

Damage and Changes in Mechanical Properties of a Gabbro Thermally Loaded up to 1,000°C

Thermal loading of rocks at high temperatures induces changes in their mechanical properties. In this study, a hard gabbro was tested in the laboratory. Specimens were slowly heated to a maximum temperature of 1,000°C. Subsequent to the thermal loading, specimens were subjected to uniaxial compression. A drastic decrease of both unconfined compressive strength and elastic moduli was observed. The thermal damage of the rock was also highlighted by measuring elastic wave velocities and by monitoring acoustic emissions during testing. The micromechanisms of rock degradation were investigated by analysis of thin sections after each stage of thermal loading. It was found that there is a critical temperature above which drastic changes in mechanical properties occur. Indeed, below a temperature of 600°C, microcracks start developing due to a difference in the thermal expansion coefficients of the crystals. At higher temperatures (above 600°C), oxidation of Fe²⁺ and Mg²⁺, as well as bursting of fluid inclusions, are the principal causes of damage. Such mechanical degradation may have dramatic consequences for many geoengineering structures.     

Comparative study on CO2 sources in soil developed on carbonate rock and non-carbonate rock in Central Guizhou

In this paper, by using concentration and carbon stable isotope the CO₂ sources of soil profiles developed on limestone, dolostone and claystone basements in Central Guizhou, China are comparatively studied. The results show that CO₂ concentration of soil profiles developed on different basements is different, having the following sequence: limestone ; dolostone;claystone. Below the soil depth of 20 cm from the surface the δ₁₃ value of CO₂ in soil profile developed on limestone ranges from -12.811%. - -13.492%.(PDB), that in soil profile developed on dolostone varys from -13.212%. - -14.271%.(PDB) and that in soil profile developed on claystone is about-20.234%. - -21.485%.(PDB). Taking the carbon isotope of soil organic matter and carbonate rock as two isotopic endmembers, the proportion of soil CO₂ generated by dissolution of carbonate rock is calculated, about 21%–25% for soil profile developed on limestone basement, 19%–21% for soil profile developed on dolostone basement. There is almost no influx of CO₂ generated by the dissolution of carbonate rock in soil profile developed on claystone basement.     

Clustering of mining-induced seismic events in equivalent dimension spaces

High energy release during seismic events induced by mining operation is one of the major dangers perturbing production in underground mines. In this work, temporal changes of seismic event parameters for one of the Rudna Mine (Poland) panels are investigated. The study aim was to find whether the temporal clustering of smaller events in different parameters can be observed before and after the high energy events (M_l?≥?3) in the mining panel. The method chosen for analysis was the study of temporal variation of fractal dimension of the seismic events parameter sets composed from: the interevent epicentral distance (dr), logarithm of seismic energy (lE), and interevent energy coefficient (dlE), which is the absolute difference between logarithms of energy of two consecutive events. Temporal variations study was performed in equivalent dimension (ED) space. The transformation of the seismic source parameters into ED space allowed to estimate and compare the temporal changes of the fractal dimension of different parameter spaces using the same method—correlation fractal dimension, and then easily compare the obtained temporal changes of fractal dimension of different parameter sets. The effect of grouping is expressed by decrease of fractal dimension, which is connected with the similarity of events parameter values. The temporal changes of the fractal dimension of seismicity before the strong induced events would indicate some initiation phase of the process leading to the high energy release. In the case of the studied Rudna Mine panel, the temporal behavior of the fractal dimension values in different parameter spaces before seismic events showed significant changes before three out of four events with CLVD dominant source mechanisms.     

CO2-induced dissolution of low permeability carbonates. Part II: Numerical modeling of experiments

We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO₂ = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO₂ Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO₂ used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8).     

Fractal characterization and frequency properties of near-fault ground motions

This study explores the irregularity and complexity of strong earthquake ground motions from the perspective of fractal geometry, and constructs a relation with the frequency content of the ground motions. The box-counting fractal dimensions and five representative period parameters of near-fault ground motions from the Chi-Chi and Northridge earthquakes are calculated and compared. Numerical results indicate that the acceleration and velocity time histories of ground motions present the statistical fractal property, and the dominant pulses of near-fault ground motions have a significant influence on their box dimensions and periods. Further, the average box dimension of near-fault impulsive ground motions is smaller, and their irregular degree of wave forms is lower. Moreover, the box dimensions of ground motions reflect their frequency properties to a large extent, and can be regarded as an alternative indicator to represent their frequency content. Finally, the box dimension D of the acceleration histories shows a considerably negative correlation with the mean period T. Meanwhile, the box dimension of the velocity histories Dye is negatively correlated with the characteristic period T and improved characteristic period Tgi.     

Triassic mid-oceanic sedimentation in Panthalassa Ocean: Sambosan accretionary complex, Japan

Abstract   The Sambosan accretionary complex of southwest Japan was formed during the uppermost Jurassic to lowermost Cretaceous and consists of basaltic rocks, carbonates and siliceous rocks. The Sambosan oceanic rocks were grouped into four stratigraphic successions: (i) Middle Upper Triassic basaltic rock; (ii) Upper Triassic shallow-water limestone; (iii) limestone breccia; and (iv) Middle Middle Triassic to lower Upper Jurassic siliceous rock successions. The basaltic rocks have a geochemical affinity with oceanic island basalt of a normal hotspot origin. The shallow-water limestone, limestone breccia, and siliceous rock successions are interpreted to be sediments on the seamount-top, upper seamount-flank and surrounding ocean floor, respectively. Deposition of the radiolarian chert of the siliceous rock succession took place on the ocean floor in Late Anisian and continued until Middle Jurassic. Oceanic island basalt was erupted to form a seamount by an intraplate volcanism in Late Carnian. Late Triassic shallow-water carbonate sedimentation occurred at the top of this seamount. Accumulation of the radiolarian chert was temporally replaced by Late Carnian to Early Norian deep-water pelagic carbonate sedimentation. Biotic association and lithologic properties of the pelagic carbonates suggest that an enormous production and accumulation of calcareous planktonic biotas occurred in an open-ocean realm of the Panthalassa Ocean in Late Carnian through Early Norian. Upper Norian ribbon chert of the siliceous rock succession contains thin beds of limestone breccia displaced from the shallow-water buildup resting upon the seamount. The shallow-water limestone and siliceous rock successions are nearly coeval with one another and are laterally linked by displaced carbonates in the siliceous rock succession.