珊瑚颗粒孔隙结构及渗流特性分析

珊瑚砂的孔隙结构对珊瑚砂的力学性质和渗透特性影响显著,采用高分辨率计算机断层扫描仪（computed tomography,简称CT）对珊瑚砂颗粒进行了扫描,建立了珊瑚砂颗粒的三维孔隙结构模型,并对孔隙的直径、体积、紧密度、球度进行了定量分析,采用数值模拟方法探究了松散堆积珊瑚砂及单颗粒珊瑚砂内孔隙结构的渗透特性。试验结果表明：孔隙数量随着孔隙直径的增加显著减少,且微孔数量较多,大孔数量较少;孔隙的直径与紧密度之间为幂函数关系,与球度之间则表现为线性关系。渗流分析结果显示,松散堆积珊瑚砂的渗透能力显示出各向异性,Z方向的绝对渗透率小于X、Y方向。珊瑚砂颗粒内的孔隙渗流速度较慢,不同颗粒之间渗流能力差距较大,部分颗粒不能发生渗流,且颗粒内孔隙的绝对渗透率仅占总绝对渗透率的1.26%。     

构造置换及其控矿规律——以吉林板石沟铁矿为例

强烈的塑性变形使华北地台东北部太古宙鞍山群中的吉林板石沟铁矿发生强烈的构造置换;造成原始仅有二三层的铁矿褶皱重复,在X(包络线)、Y(枢纽线)方向均被拉断,形成透镜状的复式褶皱勾状体。现有的19个矿组均为这种复式褶皱的转折端,并多呈"Z"型不对称形式。根据以上控矿规律本文提出两个找矿方向,一是包络线方向,另一是枢纽线方向,对1、3矿组具体地做了勘探设计。目前本文的勘探设计已得到勘探验证,新增铁矿储量数千万吨。     

Plagioclase preferred orientation by TOF neutron diffraction and SEM-EBSD

The lattice preferred orientation (LPO) of an anorthosite (composed of andesine) sampled from a highly deformed anorthositic mylonite (Grenville Province, Quebec) was measured by TOF neutron diffraction and SEM-EBSD. The quantitative texture analysis of neutron data was accomplished by using the Rietveld texture analysis with the WIMV algorithm, implemented in the program package Materials Analysis Using Diffraction (MAUD). The texture calculations of the EBSD data were performed by using the program BEARTEX. Analyses from neutron and electron diffraction data gave similar results if EBSD data are smoothed to account for grain statistics. The principal pole figures show (010) roughly parallel to the rock foliation, (001) poles exhibiting a low angle (25°) to the pole to foliation, and (100) poles close to the Y-direction (perpendicular to the lineation and foliation pole). The [100] crystallographic direction shows a maximum in the lineation direction, [010] directions concentrate near the foliation pole. The geological deformation conditions and the constructed pole figure patterns indicate that the preferred orientation could be attributed to intracrystalline slip dominantly on (010) with [100] as slip direction. Elastic properties, calculated by averaging, document weak anisotropy that has implications for the seismic structure of the lower crust.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

Petrogenesis and tectonic setting of the peralkaline Pine Canyon caldera, Trans-Pecos Texas, USA

The Pine Canyon caldera is a small (6–7 km diameter) ash-flow caldera that erupted peralkaline quartz trachyte, rhyolite, and high-silica rhyolite lavas and ash-flow tuffs about 33–32 Ma. The Pine Canyon caldera is located in Big Bend National Park, Texas, USA, in the southern part of the Trans-Pecos Magmatic Province (TPMP). The eruptive products of the Pine Canyon caldera are assigned to the South Rim Formation, which represents the silicic end member of a bimodal suite (with a “Daly Gap” between 57 and 62 wt.% SiO₂); the mafic end member consists primarily of alkali basalt to mugearite lavas of the 34–30 Ma Bee Mountain Basalt. Approximately 60–70% crystallization of plagioclase, clinopyroxene, olivine, magnetite, and apatite from alkali basalt coupled with assimilation of shale wall rock (M_a/M_c = 0.3–0.4) produced the quartz trachyte magma. Variation within the quartz trachyte–rhyolite suite was the result of 70% fractional crystallization of an assemblage dominated by alkali feldspar with subordinate clinopyroxene, fayalite, ilmenite, and apatite. High-silica rhyolite is not cogenetic with the quartz trachyte–rhyolite suite, and can be best explained as the result of  5% partial melting of a mafic granulite in the deep crust under the fluxing influence of fluorine. Variation within the high-silica rhyolite is most likely due to fractional crystallization of alkali feldspar, quartz, magnetite, biotite, and monazite. Lavas and tuffs of the South Rim Formation form A-type rhyolite suites, and are broadly similar to rock series described in anorogenic settings both in terms of petrology and petrogenesis. The Pine Canyon caldera is interpreted to have developed in a post-orogenic tectonic setting, or an early stage of continental rifting, and represents the earliest evidence for continental extension in the TPMP.     

Experimental calibration of sapphirine–spinel Fe–Mg exchange thermometer: Implication for constraints on P–T condition of Howard Hills, Napier Complex, East Antarctica

The Fe²⁺–Mg distribution coefficients between sapphirine and spinel:

were experimentally determined at pressures of 9–13 kbar and temperatures of 950–1150 °C using a natural ultrahigh-temperature (UHT) granulite with paragenesis of these minerals from the Napier Complex in East Antarctica [X_Mg = Mg / (Fe + Mg); X_Fe = Fe / (Fe + Mg)]. A new sapphirine–spinel geothermometer has been obtained as:

We applied the exchange thermometer to UHT or high-grade metamorphic rocks that were reported from various complexes in the world. If the K_D values of 2.63–4.34 obtained from low-Cr mineral pairs such as X_Cr^Spr < 0.016 and X_Cr^Spl < 0.047 were substituted into the equation, their temperature conditions would be estimated as 806–1050 °C at 11 kbar. The X_Cr means Cr / (Al + Cr(+ Fe³⁺)). These temperatures are reasonable retrograde or near peak metamorphic condition.     

伊拉克艾哈代布油田中-下白垩统碳酸盐岩储层孔隙结构及控制因素

伊拉克艾哈代布油田中-下白垩统碳酸盐岩储层孔隙结构类型划分及其差异性的受控因素和成因的认识是制约该类储层分类评价的一个关键问题。本文综合运用了常规薄片、铸体薄片、扫描电镜、电子探针、常规物性、压汞分析、阴极发光、稳定同位素分析等手段及统计分析方法对储层的主要孔隙类型、喉道类型、喉道分布类型进行了研究,并在此基础上确定了高孔中渗细喉型、中高孔低渗细喉型、中高孔特低渗微喉型、低孔特低渗极微喉型四种孔隙结构类型及其划分依据和方案,同时指出了在该类储层中好的孔隙结构具有网络状、偏细型、分选较差的孔喉特征。通过分析明确了孔隙结构差异性的主要控制因素为岩石组构和溶蚀作用,认为该套储层的孔喉空间及结构是以沉积形成的原生孔隙受同生期或准同生期的溶蚀改造为基础,并在后期表生岩溶作用下的改造而成;而先期受岩石组构控制下的溶蚀作用控制了储层孔隙结构的差异性。通过本文研究将对中东地区该类储层的形成及分类评价研究具有较好的指导意义。     

页岩基质微观孔隙结构分析新方法

页岩基质孔隙主要包含有机孔隙和无机孔隙,页岩油气在有机孔隙和无机孔隙中的渗流机理不同,对页岩中有机孔隙和无机孔隙的微观结构进行定量表征具有重要意义.首先通过扫描电子显微镜(scanning electron microscope,简称SEM)实验分别获取具有代表性的页岩无机孔隙和有机孔隙扫描电镜图像,其中,无机孔隙相对较大,其图像的分辨率较低,有机孔隙相对较小,其图像的分辨率较高;然后,通过图像处理和马尔可夫链蒙特卡洛(Markov chain Monte Carlo,简称MCMC)法重构出相应的无机孔隙数字岩心和有机孔隙数字岩心,并提出局部叠加法构建同时包含无机孔隙和有机孔隙的页岩基质孔隙数字岩心;最后对无机孔隙数字岩心、有机孔隙数字岩心和基质孔隙数字岩心的结构特征进行了对比分析.结果表明,局部叠加法构建的页岩基质孔隙数字岩心能够同时描述页岩中的无机孔隙和有机孔隙结构特征,无机孔隙本身连通性较差,有机孔隙本身连通性较好,有机孔隙的局部孔隙度和局部渗透率较高,对页岩中的流体渗流有着重要作用.该方法为页岩中不同的孔隙结构特征描述和油气在纳米尺度孔隙中的传输模拟提供了一个可靠的研究平台.      

红层软岩遇水作用的孔隙结构多重分形特征

红层软岩内部孔隙具有随机、多样化的分布特点,孔隙结构的变化是影响其宏观力学性能的关键所在。由SEM扫描电镜获取岩样不同饱水时间下的细观结构图像,根据盒维数法计算出孔隙的分形维数,发现随着饱水时间的加长,孔隙的分形维数呈现增大趋势。同时对孔隙的数量、大小进行统计分析,得出不同饱水时间下岩样内部孔隙的分布特征。基于多重分形理论,采用统计矩的方法对孔隙结构进行定量表征。结果表明,孔隙结构的分配函数与q阶次趋于线性关系,验证了该结构的自相似性与无标度性,由广义分形维数D(0)>D(1)>D(2)说明了孔隙具有多重分形特征,由多重分形谱参数分析了孔隙结构的不规则性与复杂程度,更好地表征孔隙大小各异的分布情况。结合孔隙结构的多重分形特征与岩样抗压强度,建立起孔隙结构变化与其力学性能的关联性,对研究红层软岩的损伤过程具有一定的指导意义。     

The influence of petrophysical properties on the salt weathering of porous building rocks

The influence of pore structure, water transport properties and rock strength on salt weathering is evaluated by means of a thorough rock characterisation and a statistical analysis. The pore structure was described in terms of its porosity, pore size distribution (quantified by mean pore radius) and specific surface area, density and water transport was characterised by means of water permeability (saturated flow) and capillary imbibition (unsaturated flow); whilst the rock strength test was carried out using uniaxial compressive strength, compressional and shear wave velocities, dynamic elastic constants and waveform energy and attenuation were obtained from the digital analysis of the transmitted signal. A principal component analysis and a stepwise multiple regression model was carried out in order to examine the direct relationships between salt weathering and petrophysical properties. From the principal component analysis, two main components were obtained and assigned a petrophysical meaning. The first component is mostly linked to mechanical properties, porosity and density whereas the second component is associated with the water transport and pore structure. Salt weathering, quantified by the percentage of weight loss after salt crystallisation, was included in both principal components, showing its dependence on their petrophysical properties. The stepwise multiple regression analysis found that rock strength has a predominant statistical weight in the prediction of salt weathering, with a minor contribution of water transport and pore structure parameters.     

济南岩溶含水介质孔隙结构特征

以济南泉域为例,基于2004-2019年济南岩溶泉域大气降水资料、地下水水位动态监测资料和岩芯的压汞实验结果,采用分形理论中孔隙分形维数法和数理统计的方法,探究含水岩层孔隙结构对泉水位动态的影响。结果表明:（1）炒米店组和三山子组2号样的孔隙发育均匀,孔喉分布集中,岩层储水性能好、储水空间较大;马家沟组北庵庄段的孔隙发育均匀,孔喉分布集中,岩层储水性能差,储水空间小;三山子组1号样和马家沟组东黄山段的孔隙发育不均匀,孔喉分布分散,岩石内部储水空间大;（2）岩层渗透率主要受孔喉分布的影响,济南泉域碳酸盐岩含水层的孔隙储水空间大、储水性能好,济南泉域独特的孔隙结构在维持趵突泉水位动态稳定中起到重要作用。      

Pre-eruptive conditions revealed by mega- and pheno-cryst compositions from the Quaternary Erzincan Volcanics, Eastern Turkey: Insights into the magma processes

Quaternary Erzincan Volcanics (QEVs) from the Erzincan Basin consist of mega- and pheno-cryst-bearing high-K calc-alkaline dome lavas. Fourteen nearly phenocrystic domes, with a range of basaltic-andesite, andesite, dacite and rhyolite compositions, were emplaced in the North Anatolian Fault Zone. The emplacement ages yielded by the unspiked K–Ar technique range from 102 to 140 ka. The andesitic domes (each less than 3 km in diameter) contain amphibole megacrysts. Amphibole compositions show a linear variation from ferro-edenite, edenite to pargasite from rhyolite to andesite. Pargasitic amphibole megacrysts scattered into the groundmass are very similar in composition to the microlites. All plagioclases are 53 mol%. Oscillation types are An₃₂₋₅₀ whose variations range from 10 to 16 mol% An and have 10–150 μm in thickness. Pre-eruptive conditions, calculated from mega- and pheno-cryst composition, using pyroxene and two oxide thermometers and the Al-in-hornblende barometer, ranged from 918 to 837 °C and 6.6 to 4.3 kbar for andesitic magma, 824–755 °C and 4.6–4.2 kbar for dacitic magma to 803–692 °C and 4.3–3.9 kbar for rhyolitic magma, which correspond to a depth of >10 km for storage region of the crust. The fO₂ values vary from −14.25 to −15.35 log units which are plotted just below nickel–nickel oxide (NNO) buffers. The systematic decrease in thermobarometric results from andesite to rhyolite is consistent with a single magma reservoir moving upward through the crust followed by fractional crystallization. Textural and compositional relationships of mega- and pheno-crystic phases suggest that magma mixing, fluid input to the reservoir and fractional crystallization processes, with a small amount crustal contamination play key role in evolution of the QEVs.     

Cretaceous isochron ages from K–Ar and Ar / Ar dating of eclogitic rocks in the Tso Morari Complex, western Himalaya, India

The Tso Morari Complex, which is thought to be originally the margin of the Indian continent, is composed of pelitic gneisses and schists including mafic rock lenses (eclogites and basic schists). Eclogites studied here have the mineral assemblage Grt + Omp + Ca-Amp + Zo + Phn + Pg + Qtz + Rt. They also have coesite pseudomorph in garnet and quartz rods in omphacite, suggesting a record of ultrahigh-pressure metamorphism. They occur only in the cores of meter-scale mafic rock lenses intercalated with the pelitic schists. Small mafic lenses and the rim parts of large lenses have been strongly deformed to form the foliation parallel to that of the pelitic schists and show the mineral assemblages of upper greenschist to amphibolite facies metamorphism. The garnet–omphacite thermometry and the univariant reaction relations for jadeite formation give 13–21 kbar at 600 °C and 16–18 kbar at 750 °C for the eclogite formation using the jadeite content of clinopyroxene (X_Jd = 0.48).

Phengites in pelitic schists show variable Si / Al and Na / K ratios among grains as well as within single grains, and give K–Ar ages of 50–87 Ma. The pelitic schist with paragonite and phengite yielded K–Ar ages of 83.5 Ma (K = 4.9 wt.%) for paragonite–phengite mixture and 85.3 Ma (K = 7.8 wt.%) for phengite and an isochron age of 91 ± 13 Ma from the two dataset. The eclogite gives a plateau age of 132 Ma in Ar/Ar step-heating analyses using single phengite grain and an inverse isochron age of 130 ± 39 Ma with an initial ⁴⁰Ar / ³⁶Ar ratio of 434 ± 90 in Ar/Ar spot analyses of phengites and paragonites. The Cretaceous isochron ages are interpreted to represent the timing of early stage of exhumation of the eclogitic rocks assuming revised high closure temperature (500 °C) for phengite K–Ar system. The phengites in pelitic schists have experienced retrograde reaction which modified their chemistry during intense deformation associated with the exhumation of these rocks with the release of significant radiogenic ⁴⁰Ar from the crystals. The argon release took place in the schists that experienced the retrogression to upper greenschist facies metamorphisms from the eclogite facies conditions.     

考虑渗流和剪胀的圆形巷道围岩广义SMP准则解

本文基于广义SMP准则,综合考虑中间主应力、孔隙水压力和围岩的剪胀特性,建立了巷道围岩的理想弹塑性模型。根据弹塑性理论,推导了渗流作用下的围岩应力场、位移场和塑性区半径的统一解析解。结合具体算例对SMP准则计算得到的围岩应力分布和塑性区半径与Mohr-Coulomb（M-C）准则进行了对比,并对孔隙水压力和围岩剪胀角等影响因素进行了分析。研究结果表明：SMP准则计算得到的塑性区半径小于M-C准则,说明M-C准则相较于SMP准则更为保守;孔隙水压力对巷道围岩的位移场具有很大的影响,表现为孔隙水压力越大,巷道洞壁附近围岩位移量越大,且围岩塑性区半径和弹塑性交界面处的峰值切向应力与孔隙水压力成正比;采用相关联流动法则会低估围岩的强度,不考虑剪胀则会低估围岩的实际变形。     

基于离散裂隙网络模型的节理岩体渗透张量及特性分析

节理岩体几何结构非常复杂,研究其渗流特性对于指导含水岩层稳定性分析具有重要价值。应用离散裂隙网络模型DFN方法,基于VC++6.0软件平台,建立了平面渗流分析方法,分析了节理岩体不同几何分布情况下的渗透率张量特征,通过定义渗流定向性系数对岩体渗流的定向性特征进行了定量分析。结果表明：单组节理岩体渗流具有明显的各向异性特征,渗流定向性随着节理角度变化显著;节理随着节理贯通性增加,节理渗透率呈现对数增加趋势;两组节理情况下,各向异性特征随着节理组间夹角变化;两组节理岩体渗流特征研究中,正交分布下,岩体仍存在各向异性,但渗流定向性系数较低;当节理倾角服从正态分布时,随着节理倾角标准差增大,渗透率增加;两组节理夹角不同时,节理渗透主方向倾角随着夹角增大而相应增大,基本沿两组节理夹角方向的角平分线方向。     

Anisotropy of rock elasticity behavior and of gas migration in a Variscan Carboniferous rock mass in the South Limburg, The Netherlands

Investigations of the engineering geological, physical and rock mechanical properties have been carried out on Carboniferous core samples of the research wells Geverik, Kastanjelaan and Heugem in South Limburg, The Netherlands over a total length of 2687 m. On the basis of a mathematical modeling of the deformation and anisotropy of the Variscan-folded Carboniferous, the following parameters were determined: density, porosity, water resumption and specific permeability for gas.

Dynamic, static and uniaxial compression tests were performed in order to obtain direction-dependent data like dynamic elasticity modulus, static elasticity modulus (Young's modulus), shear modulus and the Poisson's ratios. Additionally, the stress–strain behaviors were determined by uniaxial strength, tensile strength and shear strength tests. For the evaluation of the anisotropy of the rock mass caused by tectonic movements, ultrasonic measurements were of greatest importance. The results show that the ultrasonic velocities of p- and s-waves are dependent on the rock type, the tectonic configuration and the degree of saturation.

The stress–strain behavior of mud, sand and marly limestones shows an orthogonal anisotropy, while limestones and siliceous limestones are quasi-isotropic. The anisotropy is related and parallel to the main tectonic axis and also to the rock texture.     

Kinematic significance of fold- and fault-related fracture systems in the Zagros mountains, southern Iran

Initiation and formation of folds and the Kazerun high-angle fault zone, in the Zagros fold-and-thrust belt, were related to the continuing SW–NE oriented contraction that probably initiated in the Late Cretaceous, and intensified, starting in Miocene, when the Arabian and Eurasian plates collided. The contraction that led to folding and thrusting of the Phanerozoic sequence in the belt has led to the strike–slip reactivation of basement faults that formed during the Precambrian. Two major systems of fractures have developed, under the same regional state of contraction, during the folding and strike–slip faulting processes. Folding led to the formation of a system of fold-related fractures that comprises four sets of fractures, which include an axial and a cross-axial set that trend parallel and perpendicular to the confining fold axial trace, respectively, and two oblique sets that trend at moderate angles to the axial trace. Slip along high-angle, strike–slip faults formed a system of fractures in the damage zone of the faults (e.g., Kazerun), and deformed folds that existed in the shear zone by rotating their axial plane. This fault-related fracture system is made of five sets of fractures, which include the two sets of Riedel shear fractures (R and R′), P- and Y-shear fractures, and an extensional set.

Remote sensing analysis of both fracture systems, in a GIS environment, reveals a related kinematic history for folding outside of the Kazerun shear zone and faulting and deformation (fracturing and rotation of folds) within the Kazerun fault zone. Rotation of the folds and formation of the five sets of the fault-related fractures in the Kazerun shear zone are consistent with a dextral motion along the fault. The mean trends of the shortening directions, independently calculated for the fold- and fault-related fracture systems, are remarkably close (N53 ± 4°E and N50 ± 5°E, respectively), and are perpendicular to the general NW–SE trend of the Zagros fold-and-thrust belt. Although segments of the Kazerun fault are variably oriented within a narrow range, the angular relationships between sets of fault-related fractures and these segments remain constant.     

Solute distribution in porous rhyolite as evaluated by sequential centrifugation

Pore water in a porous rhyolite, having a porosity of 27% and pore radii ranging from >25 μm to 0.008 μm, was centrifugally extracted stepwise with increasing centrifugal speed to examine the potential variations of the compositions of pore water and their relationships to reaction and transport occurring in the rock. The rock was soaked for from 1 h to 7 days in an aqueous solution prior to centrifugation. To evaluate the effect of adsorption under minimum effect of dissolution, Li⁺ and Br⁻ were added to the solution as tracer ions. As centrifugal speed increased, water was extracted in order of large to small pores and the thickness of residual water film became thinner. The concentrations of ions dissolving from the rock (Na⁺, K⁺, Ca²⁺, etc.) after 7 days of immersion were relatively constant in pores of 1-10 μm radii and exponentially increased by 3-100 fold with decreasing pore radius to 0.1 μm. These ions are dissolved from the rock and transported toward the exterior of the rock by diffusion. The calculation using a reactive-transport equation showed that the observed concentration changes reflect the change in solute distribution profile with pore size. The concentration of Si after 7 days of immersion was approximately constant or slightly decreased with increasing centrifugal speed, which appears to be controlled by the solubility. The concentration of Li⁺ decreased with increasing centrifugal speed after 1 h of immersion but the trend changed after 7 days of reaction. Initial behavior of Li⁺ is explained by adsorption on pore walls, and the change of trend is explained by desorption of that previously adsorbed, slight amounts of dissolution, and inflow from the outside of the rock. The change in concentration of Br⁻ with increasing centrifugal speed was small, probably because Br⁻ was not adsorbed on the surfaces. The sequential centrifugation thus provides information on the solute distribution associated with reaction and transport occurring in rock pores.     

不同成因结构面各向异性特征及其剪切力学特性

岩体结构面具有明显的各向异性,其直接影响岩体的变形特性、力学特性与渗流特性,因此对结构面的各向异性特征开展量化分析尤为重要。针对不同成因(劈裂、剪切)结构面,利用结构面量化参数(节理粗糙度系数JRC、节理平均倾角θ、分形维数D_B)对其各向异性特征进行了分析,研究了各向异性特征对其剪切力学特性的影响。研究结果表明:(1)在劈裂断裂结构面中,平行劈裂方向的JRC与θ值普遍大于垂直方向,且随角度变化波动较小,D_B在对角线方向变化较大,其值与所取剖面线长度有关;而在剪切断裂结构面中,平行剪切方向的JRC和θ值与垂直方向无明显差别,但D_B同样在对角线方向变化较大;(2)在评价结构面各向异性时,采用θ、D_B等参数评价时,劈裂断裂结构面与剪切断裂结构面各向异性系数无明显差别,采用JRC作为评价参数时,其各向异性系数差异较大,能较好反映不同结构面之间的差异特征;(3)剪切断裂结构面的峰值剪切荷载和法向位移均高于劈裂断裂结构面,两种结构面的剪胀角达到峰值时的剪切位移相近,剪切断裂结构面的开度分布较为集中且普遍较大,劈裂断裂结构面开度分布则较为分散。