Hydraulic parameters and solute velocities in triple-porosity karstic-fissured-porous carbonate aquifers: case studies in southern Poland

 Rock and flow parameters of three karstic-fissured-porous aquifers in the Krakow-Silesian Triassic formations were measured using various methods and compared. Though cavern and fissure porosities are shown to be very low (cavern porosity below 0.5% and fracture porosity below 0.2%), they contribute dominantly to the hydraulic conductivity (from about 1.3×10^–6 to about 11×10^–6 m/s). Matrix porosity (2–11%) is shown to be the main water reservoir for solute transport and the main or significant contributor to the specific yield (<2%). Though the matrix porosity is shown to be much larger than the sum of the cavern and fissure porosities, its contribution to the total hydraulic conductivity is practically negligible (hydraulic conductivity of the matrix is from about 5×10^–11 m/s to about 2×10^–8 m/s). On the other hand, the matrix porosity (for neglected cavern and fissure porosities) when combined with tracer ages (or mean travel times) is shown to yield proper values of the hydraulic conductivity (K) by applying the following formula: K≅(matrix porosity×mean travel distance)/(mean hydraulic gradient×mean tracer age). Confirming earlier findings of the authors, this equation is shown to be of great practical importance because matrix porosity is easily measured in the laboratory on rock samples, whereas cavern and fracture porosities usually remain unmeasurable. Received: 21 February 1997 · Accepted: 13 May 1997     

压力溶解对颗粒聚集岩体中热-水-应力耦合影响的数值模拟

在自主研制的孔隙介质热-水-应力耦合有限元程序中引入Taron等提出的颗粒聚集体的压力溶解模型,针对一个假设的实验室尺度且位于非饱和石英颗粒聚集岩体中的高放废物地质处置模型,拟定两种计算工况：（1）孔隙率和渗透系数是压力溶解的函数;（2）孔隙率和渗透系数均为常数,进行4 a处置时段的数值模拟,考察了岩体中的温度、颗粒界面水膜及孔隙中的溶质浓度、迁移和沉淀质量、孔隙率及渗透系数、孔隙水压力、地下水流速和应力的变化、分布情况。研究结果表明：工况1计算终了时,压力溶解使得孔隙率和渗透系数分别下降到初始值的43%～54%、4.4%～9.1%。在核废料释热温度场的作用下,工况1、2中的负孔隙水压力分别为初始值的1.00～1.25倍、1.00～1.10倍,前者表现了压力溶解的明显影响;两种工况的岩体中的应力量值及分布基本相同。     

Evolution of pore space in the Poza Rica trend (Mid-Cretaceous), Mexico

ABSTRACT The Poza Rica trend of the Tampico embayment, Mexico, will ultimately produce more than 2.3 × 10⁹ barrels of oil from Mid-Cretaceous (Albian-Cenomanian) basin-margin deposits. Bioclastic grainstone, packstone, and wackestone are interbedded with polymictic lime breccia and dolomitized debris; all were deposited by sediment gravity flow. Indigenous sediment was pelagic lime mud. Typical reservoir porosities are about 10%; permeabilities average 2 md and rarely exceed 100 md. Porosity is largely the result of selective dissolution of rudist fragments, which were originally aragonite. Detailed petrographic study, with emphasis on the diagenetic products, allows quantitative assessment of porosity at various diagenetic stages from original sediment to reservoir rock. A relatively simple diagenetic history is reflected by about 90% of the samples studied: primary porosity was reduced through lithification of matrix mud and initial cementation by clear, equant to bladed, non-ferroan calcite. Later dissolution produced extensive skelmoldic and minor vuggy porosity. Subsequently, non-ferroan calcite cement reduced porosity before the emplacement of hydrocarbons. Reconstructed sediment porosities are comparable to, but lower than, modern counterparts. The initial phase of cementation and presumed lithification of mud greatly reduced porosity in all lithologies, but appreciably more porosity persisted in grainstone and packstone than in wackestone or mudstone. Dissolution produced a porosity resurrection, which exceeded that of the initial sediment in some grainstones. Calcite cementation and local multiphase quartz cementation and dolomitization reduced porosity to present average values of 8–12% in grain-supported rocks and 3% in mud-supported rocks. The greater persistence of primary porosity and, therefore, permeability in grain-supported rocks probably accounts for their greater secondary porosity development and ultimate reservoir quality. Geometrically averaged permeabilities range from only 0.17 md in wackestone to 3.85 md in dolomite, but differ significantly with rock type and grain size. Permeability increases with porosity in all lithologies; the rate of increase is greater at higher porosities and with coarser grain sizes. The agent for both early cementation and development of secondary porosity appears to have been meteoric water. Subaerial exposure appears to be ruled out, however, by a basin-margin depositional environment and continued burial beneath Upper Cretaceous pelagic sediments. Early exposure to meteoric water can be explained by a hydrologic head developed during penecontemporaenous exposure that produced cavernous porosity in the adjacent Golden Lane trend. Descending meteoric water likely emerged as submarine springs along the Tamabra trend. Deposition of pelagic limestone during the Turonian blanketed part of the Golden Lane escarpment to enhance development of a large freshwater lens; gaps in the blanket localized springs and influenced flow patterns within the Poza Rica field. Analogous freshwater circulation exists today in northern Florida.     

节理粗糙度对岩块体积的影响规律研究

岩体结构直接控制着岩体力学、水力学特性,岩块体积分布规律可直接反映岩体结构特征。目前研究多将节理面简化为平面,忽略了粗糙度对岩块体积的影响。利用Hurst指数H和高度均方根R_q表征节理粗糙度,研究了节理粗糙度对岩块数量和体积分布规律的影响。结果表明：(1)H、R_q控制节理的粗糙起伏形态,粗糙度随着R_q的增大而增大、随着H的增大而减小;(2)块体的数量总体上与节理粗糙度正相关,即随着粗糙度的增大而增加;(3)R_q和H对岩块体积分布的影响,主要通过改变小体积块体所占比例实现,块体的平均体积和中位体积随着粗糙度的增大而减小;(4)当节理正交且间距相近时,块体数量随粗糙度的变化可以分为3个阶段,即稳定阶段、初始增长阶段和快速增长阶段,R_q和H共同控制着各区间的相对分布范围。最后,基于摄影测量采集的数据建立了岩体节理模型,对大连某公路边坡的岩块体积分布进行了研究,验证了上述结论的准确性。     

Comparative study of tracer diffusion of HTO, Na and Cl in compacted kaolinite, illite and montmorillonite

The through-diffusion of HTO, ²²Na⁺ and ³⁶Cl⁻ in kaolinite, homo-ionic Na-illite and homo-ionic Na-montmorillonite was measured at a high degree of compaction as a function of the salt concentration in the ‘external solution’, i.e. in the solution in contact with the clay sample. The clays were chosen for this study because of their differences in the number and nature of ion exchange sites leading to different proportions of interlayer-, inter-particle and free pore water. It was found that the diffusive mass transfer of Na⁺ in Na-montmorillonite and Na-illite increased with decreasing external salt concentration, while the opposite trend was observed for the diffusion of Cl⁻. These trends are more pronounced in the case of Na-montmorillonite than in Na-illite, while almost no salt effect was observed for kaolinite. Similarly no salt effect was observed for the diffusion of HTO through all of the clays tested. These observations are in agreement with a conceptual model where it is assumed that cations diffuse preferentially in the interlayer or diffuse double-layer porosity, while anions are almost completely excluded from these regions. In the case of Na⁺ diffusion, the salt effects can be explained by an influence on the concentration gradient of diffusing cations, while in the case of Cl⁻ the external salt concentration has an effect on the accessible porosity. Effective diffusion coefficients of Cl⁻ fulfil the same relationship to porosity as those of the uncharged HTO, when using accessible porosities for such a comparison. Furthermore it is shown that pore diffusion coefficients for the three tracers are fairly well correlated with the respective diffusion coefficients in bulk water, if the effective diffusion coefficients for Na⁺ are derived from calculated tracer concentration gradients in the interlayer or diffuse double-layer porosities.     

塑性扩容梯度对孔隙介质岩体中T-H-M耦合作用的二维有限元分析

由双重孔隙-裂隙介质热-水-应力耦合模型退化为单一孔隙介质模型,将其与岩体扩容梯度引入笔者所研制的二维有限元程序中,使用Mohr-Coulomb准则,计入塑性扩容对岩体孔隙率及渗透系数的影响,针对一个假设的实验室尺度且位于饱和孔隙介质岩体中的高放废物地质处置库模型,拟定不同扩容梯度值的5种工况,进行4年处置时段的数值模拟,考察了岩体中的温度、正应力、塑性区、孔隙率及渗透系数、孔隙水压力和地下水流速的变化、分布情况。结果主要显示,相比于不考虑扩容梯度的工况,考虑扩容梯度工况的正应力、孔隙率及渗透系数、孔隙水压力和地下水流速等的分布与塑性区的分布有明显的对应关系,呈现了某种"剪切带效应";正应力量值、塑性区面积、孔隙率及渗透系数、孔隙水压力、地下水流速等均随所取扩容梯度值的变大而增加。     

Strontium and oxygen isotope profiles across marble-silicate contacts,Lizzies Basin,East Humboldt Range,Nevada: constraints on metamorphic permeability contrasts and fluid flow

 Sr isotope profiles across marble-silicate rock contacts are used in conjunction with previously published oxygen isotope profiles to constrain fluid movement, porosity and permeability contrasts in migmatitic rocks from Lizzies Basin in the East Humboldt Range, Nevada. The ¹⁸O/¹⁶O systematics in the high-grade sequence have been interpreted to reflect infiltration of ∼2×10² m³/m² of a relatively low ¹⁸O hydrous fluid through the sequence, but with preservation of δ¹⁸O anomalies in thin marble bands due to a 30-fold lower porosity in the marble compared with silicate rocks (Wickham and Peters 1992). The Sr isotope profiles confirm that tracer exchange between marble and silicate rock was primarily by diffusion, and in one case, indicate that porosities differed by less than a factor of four in the ∼10 cm boundary layer which exhibits diffusive modification of ⁸⁷Sr/⁸⁶Sr ratios. This contrasts with modelling of the oxygen isotope profiles which imply porosity contrasts >10 for one marble band and >50 for a second marble band. Either strontium and oxygen isotope diffusion reflect different events (possible if fluid Sr contents varied with time) or porosity varied substantially with the silicate rocks. Oxygen isotope profiles in the deeper part of the metamorphic section in which δ¹⁸O values of silicate rocks have been homogenised and lowered, indicate similar diffusion distances (and thus porosity-time evolution) to oxygen isotopic profiles higher in the section. Comparison of strontium and oxygen isotope diffusion distances constrains fluid Sr contents to between ∼50 and ∼500 ppm deep in the section, but less than ∼10 ppm higher in the section. The difference is related to release of relatively saline, Sr-rich fluids, by the abundant leucogranites and associated skarns deep in the section (cf. Peters and Wickham 1995). Received: 9 December 1994/Accepted: 13 April 1995     

Microporosity in Alberta Plains coals

The nature of porosity in 11 Alberta Plains coals is scanned using different characterization techniques, namely, gas adsorption, helium and mercury displacement, and mercury porosimetry. Surface areas calculated from carbon dioxide adsorption at 25°C (varying between 75 and 506 m²/g) are consistently higher than those calculated from nitrogen at −196°C (varying between 1 and 9 m²/g). The samples from Paskapoo deposits have slightly higher surface areas as compared to samples from Horseshoe Canyon deposits. The carbon content in the tested coals (73.4 to 76.8%) is too narrow of range to determine the influence of this parameter on surface area and porosity.Total pore volumes are measured in the range 1.2 to 3000 nm and then divided into macropores (30 to 3000 nm), transitional pores (1.2 to 30 nm) and micropores (less than 1.2 nm). Almost all the surface area in coals is located in the micropores. Coals vary widely in their total pore volumes and porosities range between 9.5 and 22.3%. Percent pore volume contained in micropores varies between 29.5 and 75.9, in transitional pores between 1.1 and 20.5 and in macropores between 18.5 and 50.0.     

Order of diagenetic events controls evolution of porosity and permeability in carbonates

The properties of carbonate rocks are often the result of multiple, diagenetic events that involve phases of cementation (porosity occlusion) and dissolution (porosity enhancement). This study tests the hypothesis that the order of these events is a major control on final porosity and permeability. A three-dimensional synthetic model of grainstone is used to quantify trends that show the effect of early cementation, non-fabric selective dissolution, and then a second-generation of (post-dissolution) cement. Models are 3 mm³ with a resolution of 10 μm. Six simple paragenetic sequences are modelled from an identical starting sediment (without accounting for compaction) where the same diagenetic events are placed in different sequences, allowing for quantification of relative changes in the resultant porosity and permeability for each diagenetic event, the trajectory through time, as well as for each final rock. All modelled paragenetic sequences result in reductions in porosity and permeability, but the order of diagenetic events controls the trajectory and final rock properties. Differences in the order of early cement precipitation alone produce variable final values, but all follow the porosity–permeability relationship as expressed by the Kozeny-Carman equation. However, final values for the sequences which include a phase of dissolution fall on a new curve, which departs from that predicted by the Kozeny-Carman relationship. This allows an alternative form of porosity–permeability relationship to be proposed: κ = 2280ϕ–30,400, where ϕ is porosity (%) and κ is permeability (mD). Hence while the Kozeny-Carman relationship predicts porosity–permeability changes that occur with cementation, it is unable to capture accurately changes within the pore network as a result of dissolution. Although the results may be dependent on the properties of the initial carbonate sediment and simplified diagenetic scenarios, it is suggested that this new porosity–permeability relationship may capture some generalized behaviour, which can be tested by modelling further sediment types and diagenetic scenarios.     

Fluid and enthalpy production during regional metamorphism

Models for regional metamorphism have been constructed to determine the thermal effects of reaction enthalpy and the amount of fluid generated by dehydration metamorphism. The model continental crust contains an average of 2.9 wt % water and dehydrates by a series of reactions between temperatures of 300 and 750° C. Large scale metamorphism is induced by instantaneous collision belt thickening events which double the crustal thickness to 70 km. After a 20 Ma time lag, erosion due to isostatic rebound restores the crust to its original thickness in 100 Ma. At crustal depths greater than 10 km, where most metamorphism takes place, fluid pressure is unlikely to deviate significantly from lithostatic pressure. This implies that lower crustal porosity can only be maintained if rock pores are filled by fluid. Therefore, porosities are primarily dependent on the rate of metamorphic fluid production or consumption and the crustal permeability. In the models, permeability is taken as a function of porosity; this permits estimation of both fluid fluxes and porosities during metamorphism. Metamorphic activity, as measured by net reaction enthalpy, can be categorized as endothermic or exothermic depending on whether prograde dehydration or retrograde hydration reactions predominate. The endothermic stage begins almost immediately after thickening, peaks at about 20 Ma, and ends after 40 to 55 Ma. During this period the maximum and average heat consumption by reactions are on the order 11.2·10^–14 W/cm³ and 5.9·10^–14 W/ cm³, respectively. The maximum rates of prograde isograd advance decrease from 2.4·10^–8 cm/s, for low grade reactions at 7 Ma, to 7·10^–10 cm/s, for the highest grade reaction between 45 and 58 Ma. Endothermic cooling reduces the temperature variation in the metamorphic models by less than 7% (40 K); in comparison, the retrograde exothermic heating effect is negligible. Dehydration reactions are generally poor thermal buffers, but under certain conditions reactions may control temperature over depth and time intervals on the order of 1 km and 3 Ma. The model metamorphic events reduce the hydrate water content of the crust to values between 1.0 and 0.4 wt % and produce anhydrous lower crustal granulites up to 15 km in thickness. In the first 60 Ma of metamorphism, steady state fluid fluxes in the rocks overlying prograde reaction fronts are on the order of 5·10^–11 g/cm²-s. These fluid fluxes can be accommodated by low porosities (<0.6%) and are thus essentially determined by the rate of devolitalization. The quantity of fluid which passes through the metamorphic column varies from 25000 g/cm², within 10 km of the base of the crust, to amounts as large as 240000 g/cm², in rocks initially at a depth of 30 km. Measured petrologic volumetric fluid-rock ratios generated by this fluid could be as high as 500 in a 1 m thick horizontal layer, but would decrease in inverse proportion of the thickness of the rock layer. Fluid advection causes local heating at rates of about 5.9·10^–14 W/cm³ during prograde metamorphism and does not result in significant heating. The amount of silica which can be transported by the fluids is very sensitive to both the absolute temperature and the change in the geothermal gradient with depth. However, even under optimal conditions, the amount of silica precipitated by metamorphic fluids is small (<0.1 vol %) and inadequate to explain the quartz veining observed in nature. These results are based on equilibrium models for fluid and heat transport that exclude the possibility of convective fluid recirculation. Such a model is likely to apply at depths greater than 10 km; therefore, it is concluded that large scale heat and silica transport by fluids is not extensive in the lower crust, despite large time-integrated fluid fluxes.     

Estimation of Block Sizes for Rock Masses with Non-persistent Joints

Summary  Discontinuities or joints in the rock mass have various shapes and sizes. Along with the joint orientation and spacing, the joint persistence, or the relative size of the joint, is one of the most important factors in determining the block sizes of jointed rock masses. Although the importance of joint persistence on the overall rock mass strength has long been identified, the impact of persistence on rock strength is in most current rock mass classification systems underrepresented. If joints are assumed to be persistent, as is the case in most designs, the sizes of the rock blocks tend to be underestimated. This can lead to more removable blocks than actually exist in-situ. In addition, a poor understanding of the rock bridge strength may lead to lower rock mass strengths, and consequently, to excessive expenditure on rock support. In this study, we suggest and verify a method for the determination of the block sizes considering joint persistence. The idea emerges from a quantitative approach to apply the GSI system for rock mass classification, in which the accurate block size is required. There is a need to statistically analyze how the distribution of rock bridges according to the combination of joint orientation, spacing, and persistence will affect the actual size of each individual block. For this purpose, we generate various combinations of joints with different geometric conditions by the orthogonal arrays using the distinct element analysis tools of UDEC and 3DEC. Equivalent block sizes (areas in 2D and volumes in 3D) and their distributions are obtained from the numerical simulation. Correlation analysis is then performed to relate the block sizes predicted by the empirical equation to those obtained from the numerical model simulation. The results support the concept of equivalent block size proposed by Cai et al. (2004, Int. J. Rock Mech. Min. Sci., 41(1), 3–19).     

Welded tuff porosity characterization using mercury intrusion,nitrogen and ethylene glycol monoethyl ether sorption and epifluorescence microscopy

Porosity of welded tuff from Snowshoe Mountain, Colorado, was characterized by mercury intrusion porosimetry (MIP), nitrogen sorption porosimetry, ethylene glycol monoethyl ether (EGME) gas phase sorption and epifluorescence optical microscopy. Crushed tuff of two particle-size fractions (1-0.3 mm and less than 0.212 mm), sawed sections of whole rock and crushed tuff that had been reacted with 0.1 N hydrochloric acid were examined. Average MIP pore diameter values were in the range of 0.01–0.02μm. Intrusion volume was greatest for tuff reacted with 0.1 N hydrochloric acid and least for sawed tuff. Cut rock had the smallest porosity (4.72%) and crushed tuff reacted in hydrochloric acid had the largest porosity (6.56%). Mean pore diameters from nitrogen sorption measurements were 0.0075–0.0187 μm. Nitrogen adsorption pore volumes (from 0.005 to 0.013 cm³/g) and porosity values (from 1.34 to 3.21%) were less than the corresponding values obtained by MIP. More than half of the total tuff pore volume was associated with pore diameters < 0.05μm. Vapor sorption of EGME demonstrated that tuff pores contain a clay-like material. Epifluorescence microscopy indicated that connected porosity is heterogeneously distributed within the tuff matix; mineral grains had little porosity. Tuff porosity may have important consequences for contaminant disposal in this host rock.     

密度测井资料在确定岩层孔隙度上的应用

孔隙度是储层评价的重要参数之一.密度测井是一种重要的孔隙度测井方法,在石油勘探领域常常被用作确定岩层孔隙度.通过研究密度测井的基本原理、仪器构成、刻度方法等,推导出不同情况下孔隙度的计算公式.并利用该公式进行Matlab编程,对实际测井数据进行处理,计算了多个井段的岩层孔隙度,证明了密度测井在孔隙度计算方面具有很高的应用价值.     

A method for estimating crack-initiation stress of rock materials by porosity

Crack-initiation stress of a rock under compression is the stress level that marks the initiation of the rock microfracturing process or in other words, the onset of new damage to the rock. This paper proposed a simple methodology with justifications to explore the feasibility of using total and effective porosities as estimators of crack-initiation stress of brittle crystalline rock materials under uniaxial compression. The validity/applicability of the proposed method was examined by an experimental study of granitic materials from Malanjkhand, Madhya Pradesh. It was found that effective porosity depicts better correlation with crack-initiation stress than with uniaxial compressive strength of the granitic materials. On the other hand, total porosity does not show any perceptible correlation with uniaxial compressive strength and crack-initiation stress. Plausible reasons for the nature of the obtained results were also explained in view of rock failure process under compression. It is concluded that following the proposed method, effective porosity can be used as a physical index to obtain a quick estimate of crack-initiation stress of the investigated rocks empirically.     

基于图像分割的煤岩割理CT图像各向异性特征

利用工业CT对自然煤岩样进行断层扫描观测。针对煤岩裂隙系统的多尺度、各向异性特征,应用Canny算子图像分割与方向性边缘检测技术,提取煤岩CT图像割理的总体特征、水平方向和垂直方向特征;根据特征图像计算了煤岩样的总体分形维数、孔隙度,各向异性分形维数与孔隙度及其在三维空间中的分布;讨论了分形维数与孔隙度、渗透率之间的关系,并根据煤岩样的分形维数、孔隙度对实际工程岩体的孔隙度和渗透率进行了外推计算。研究表明,煤岩样不同扫描断面的分形维数和孔隙度不同,同一煤样同一断面不同方向的分形维数与孔隙度亦不相同。利用图像分割与边缘检测对工业CT图像进行分析,可以对煤岩的各向异性分形维数与孔隙度在2D与3D空间进行精细描述。      

Mineralogical alteration and associated permeability changes induced by a high-pH plume: Modeling of a granite core infiltration experiment

Within the framework of the HPF project (Hyperalkaline Plume in Fractured Rock) at the Grimsel Test Site (Switzerland), a small scale core infiltration experiment was performed at the University of Bern. A high-pH solution was continuously injected, under a constant pressure gradient, into a cylindrical core of granite containing a fracture. This high-pH solution was a synthetic version of solutions characteristic of early stages in the degradation of cement. The interaction between the rock and the solutions was reflected by significant changes in the composition of the injected solution and a decrease in the permeability of the rock. Changes in the mineralogy and porosity of the fault gouge filling the fracture were only minor. One-dimensional reactive transport modeling, using a modified version of the GIMRT code, was used to interpret the results of the experiment. Dispersive and advective solute transport, mineral reaction kinetics and a coupling between porosity and permeability changes were taken into account. In order to obtain a reasonable agreement between models and experimental results, reactive surface areas of the order of 10⁵ m²/m³ rock had to be used. These values are much smaller than the values measured for the fault gouge filling the fracture, which are in the order of 10⁶–10⁷ m²/m³ rock. However, the results could be improved by adding a small fraction of fine grained mineral, which could explain the high initial peaks in Al and Si concentration. With the inclusion of this fine grained fraction, the initial surface areas in the model were within the range of the measured specific surface areas of the fault gouge. The fact that the decrease in permeability was significant despite the minor changes in mineralogy, suggests that permeability may be controlled by changes in the structure of the rock (pore geometries) rather than by only the bulk volumetric (porosity) changes.     

On the theory of consolidation with double porosity—III A finite element formulation

We employ a finite element method to consider the numerical solution of Aifantis' equations of double porosity.^1–5 The physical and mathematical foundations of the theory of double porosity were considered in two previous papers, ^6,7 where analytical solutions of the relevant equations are also given. For practical purposes and for the completeness of the presentation, we re-derive the basic equations without resorting to the more rigorous scheme of multiporosity theory.^1–6 Instead, we follow the more familiar approach of Biot^8,9 which we now extend from the case of single porosity to the case of double porosity. Then a general finite element formulation of the relevant eqations is thoroughly discussed and explicitly described for the two-dimensional case where four-node rectangular non-conforming isoparametric elements are employed. Three numerical examples are presented in detail to illustrate the method and assess the differences between single and double porosity models in consolidation.     

Modelling approaches for anion-exclusion in compacted Na-bentonite

The accessible porosity for Cl⁻ in bentonite is smaller than the the total porosity due to anion repulsion (exclusion) by the surface of montmorillonite, the main mineral in bentonite. The accessible porosity is a function of the bentonite density and the salt concentration. Anion exclusion data were gathered from the literature, reprocessed in a coherent data set, and modelled using four different models. Very simple models, with or without anion access to the interlayer space, are successful in reproducing trends in anion exclusion in bentonite as a function of ionic strength in the external solution and montmorillonite bulk dry densities in the bentonite. However, a model that considers clay microstructure changes as a function of bentonite compaction and ionic strength is necessary to reproduce observed trends in the data for all experimental conditions within a single model. Our predictive model excludes anions from the interlayer space and relates the interlayer porosity to the ionic strength and the montmorillonite bulk dry density. This presentation offers a good fit for measured anion accessible porosities in bentonites over a wide range of conditions and is also in agreement with microscopic observations.     

砂岩储层孔隙保存的定量预测研究

对北方地区若干盆地的砂岩孔隙发育特征的研究表明,地温场、地质年代和盆地沉降方式对砂岩孔隙的演化和保存有制约作用。地温梯度每增加1℃,砂岩孔隙度平均减小约7%;在地温梯度2-4℃/100m范围内,有效储层的保存深度差异可达2500-3000m.地质年代每增加1Ma,砂岩孔隙度降低约0.018%-0.009%.地层超压可最大保存5%-7%的孔隙度。盆地沉降方式不同引起的孔隙保存量的差异为2%-5%,相应的有效孔隙保存深度的差值约1000m。     

西南天山中新生代盆地成矿流体运移规律

乌拉根铅锌矿床和萨热克铜矿床是西南天山中新生代盆地最有代表性的两个层控砂砾岩型矿床,乌拉根铅锌矿床产于下白垩统克孜勒苏群第五岩性段(K1 kz5)的粗砂质细砾岩中,后期经历了弱的构造改造作用;萨热克铜矿床产于上侏罗统库孜贡苏组上段(J3 k2)杂砾岩中,其北矿段后期经历了强烈的构造改造作用,南矿段可见岩浆热液蚀变作用后...