Formation of bed‐scale spatial patterns in dolomite abundance during early dolomitization: Part I. Mechanisms and feedbacks revealed by reaction–transport modelling

Dolomites of varied ages exhibit metre‐scale nested patterns of lateral periodic variation in permeability and porosity and, by inference, dolomite abundance as most examples are 100% dolomite. Two‐dimensional reaction–transport modelling simulations of bed‐scale dolomitization were used to assess whether those patterns in dolomite abundance could form during near‐surface replacement dolomitization. Simulations used a 2 m high and 18 m long model domain, a low‐Mg calcite grainstone precursor and an evaporated Mississippian seawater brine (430 parts per thousand salinity) as the dolomitizing fluid. The domain was initially populated with random variations in porosity and/or grain size. Results reveal that spatial patterns in dolomite abundance emerge when there is as little as 1% dolomite formed, with similarities between the modelled patterns and outcrop‐documented patterns. The nested patterns include a near‐random component that constitutes ≤40% of the total variance, short‐range correlation ranging from 1·5 to 3·3 m and a longer‐range periodic trend with a wavelength up to 6·5 m. The emergence of pattern in dolomite abundance is the result of an autogenic self‐organizing phenomenon. It is triggered by variation in initial calcite reactive surface area that occurs due to the random heterogeneities in initial porosity and/or grain sizes. The pattern develops due to a combination of kinetic disequilibrium reactions (dolomite precipitation and calcite dissolution) and positive feedbacks between dolomite growth, calcite dissolution and fluid flow. Flow is around loci of higher dolomite, lower porosity and higher reactive surface areas, but through loci of lower dolomite, higher porosity and lower reactive surface areas. The resulting less porous/more dolomite and more porous/less dolomite structures at the metre‐scale arise from those localized interactions. This self‐organizing mechanism for pattern formation constitutes a new model for geochemical self‐organization during dolomitization and is the only self‐organization model that is proven applicable to the formation of metre‐scale patterns during early, near‐surface dolomitization.     

Reaction mechanism for the replacement of calcite by dolomite and siderite: implications for geochemistry, microstructure and porosity evolution during hydrothermal mineralisation

Carbonate reactions are common in mineral deposits due to CO₂-rich mineralising fluids. This study presents the first in-depth, integrated analysis of microstructure and microchemistry of fluid-mediated carbonate reaction textures at hydrothermal conditions. In doing so, we describe the mechanisms by which carbonate phases replace one another, and the implications for the evolution of geochemistry, rock microstructures and porosity. The sample from the 1.95 Moz Junction gold deposit, Western Australia, contains calcite derived from carbonation of a metamorphic amphibole—plagioclase assemblage that has further altered to siderite and dolomite. The calcite is porous and contains iron-rich calcite blebs interpreted to have resulted from fluid-mediated replacement of compositionally heterogeneous amphiboles. The siderite is polycrystalline but nucleates topotactically on the calcite. As a result, the boundaries between adjacent grains are low-angle boundaries (<10°), which are geometrically similar to those formed by crystal–plastic deformation and recovery. Growth zoning within individual siderite grains shows that the low-angle boundaries are growth features and not due to deformation. Low-angle boundaries develop due to the propagation of defects at grain faces and zone boundaries and by impingement of grains that nucleated with small misorientations relative to each other during grain growth. The cores of siderite grains are aligned with the twin planes in the parent calcite crystal showing that the reactant Fe entered the crystal along the twin boundaries. Dolomite grains, many of which appear to in-fill space generated by the siderite replacement, also show alignment of cores along the calcite twin planes, suggesting that they did not grow into space but replaced the calcite. Where dolomite is seen directly replacing calcite, it nucleates on the Fe-rich calcite due to the increased compatibility of the Fe-bearing calcite lattice relative to the pure calcite. Both reactions are interpreted as fluid-mediated replacement reactions which use the crystallography and elemental chemistry of the calcite. Experiments of fluid-mediated replacement reactions show that they proceed much faster than diffusion-based reactions. This is important when considering the rates of reactions relative to fluid flow in mineralising systems.     

Dolomite characteristics and diagenetic model of the Calcari Grigi Group (Asiago Plateau,Southern Alps – Italy): an example of multiphase dolomitization

A multidisciplinary study, conducted over the carbonate platform deposits of the Liassic Calcari Grigi Group (Southern Alps), highlighted how the use of outcrop analogues can contribute to better define the distribution of dolomitic bodies related to fault networks, to characterize the petrophysical properties of the dolomitic sequence and unravel a complex diagenetic history. This study was carried out in the Asiago Plateau (southernmost part of the eastern Southern Alps, northern Italy) which provides excellent outcrops of the Jurassic Calcari Grigi Group. The dolomitization of the Jurassic sequence is variable in terms of stratigraphic extension and geographic distribution. In the studied localities the dolomitization is generally limited to the Mount Zugna Formation and is characterized by an undulatory front, with ‘sub‐vertical dolomitic chimneys’ along the major faults. Within this unit, and often associated with faults, stacked high‐porosity and permeability bed‐parallel dolomitic bodies are developed that show excellent petrophysical properties. The dolomitic intervals are characterized by pervasive unimodal and patchy polymodal dolomite crystals. Thin section, cathodoluminescence, isotopic and fluid inclusion analyses were used to constrain the paragenetic evolution of the sequence which is similar in all the studied localities. The first dolomitization stage is marked by zoned dolomite crystals with a dull luminescent core. The porosity is thought to have increased after this stage, with dark blue luminescent dolomite accompanied by the corrosion of older crystals. The appearance of saddle dolomite marks the onset of the porosity reduction stage, ending with the infilling of vugs and the remaining open pores with calcite cement. The diagenetic evolution locally stopped at the saddle dolomite stage with the complete occlusion of the remaining pores. Paragenetic and fluid‐inclusion data suggest an evolutionary trend of increasing temperatures and decreasing salinity toward brackish fluids responsible for dolomite and calcite precipitation. The integration of the available data seem to indicate that the diagenetic evolution of the study area is related to: (i) the interplay between evolving fluids (from marine to brackish); (ii) the burial of the sequence (increasing temperature); and (iii) the evolution of the hydrogeological system (fault and fracture network, fluid mixing). This complex paragenetic evolution is strongly linked to the evolution of the porosity framework that evolved from a good, widespread network in the early stages of the burial history to a confined system in the later stages due to reduction of porosity by the deposition of late calcite and dolomite cements.     

渤中凹陷北部陡坡带热液活动及其对湖相碳酸盐岩储层的影响

渤海海域渤中凹陷北陡坡带沙一段发育湖相碳酸盐岩及相关的混积岩类油藏,探明储量超亿方,产量高,埋深超过3 200~3 800 m.分析表明:湖相碳酸盐岩段及混积岩段发生了强烈的白云岩化作用,溶蚀作用强,储层物性好,孔隙度可超过25%.早期研究认为生屑的溶蚀是储层物性改善的重要因素.进一步研究表明,除了上述因素外,储层还受到了深部热流体的影响.通过包裹体激光拉曼分析证实,储层中存在大量含CO₂、N₂、H₂S、H₂的包裹体,CO₂最高含量超过70%;研究区多个构造富含CO₂,碳同位素及氦同位素显示其为幔源成因;对储层中大量的自生黄铁矿分析表明,大部分黄铁矿Co/Ni比值大于1,显示了热液成因;在缝洞及大孔隙内沉淀典型微晶鞍形白云石,具弧形晶面,波状消光,扫描电镜下鞍形白云石内见密集的生长纹.综合分析表明,富CO₂深部幔源热液侵位时间为3 Ma以来,侵位后造成了白云石晶间残余方解石的大量选择性溶解;酸性流体对粒内长石也产生了强烈的溶蚀作用,形成了大量砾石及岩屑的铸模孔并沉淀高岭石,流体上述溶解作用对储层段物性改善起到了决定性影响.随着方解石及长石的溶解,流体PH值的升高,盐度的增加,储层出现明显的白云岩化作用,晚期白云石的较高的均一温度及鞍形白云石证实了深部热液对白云岩化的贡献.幔源流体主要来自深大断裂及早期火山后期的排气、排液作用.      

Estimation of dolomite formation: Dolomite precipitation and dolomitization

Reactive-transport models are developed here that produce dolomite via two scenarios: primary dolomite (no CaCO₃ dissolution involved) versus secondary dolomite (dolomitization, involving CaCO₃ dissolution). Using the available dolomite precipitation rate kinetics, calculations suggest that tens of meters of thick dolomite deposits cannot form at near room temperature (25-35°C) by inorganic precipitation mechanism, though this mechanism will provide dolomite aggregates that can act as the nuclei for dolomite crystallization during later dolomitization stage. Increase in supersaturation, Mg⁺²/Ca⁺² ratio and CO₃^-2 on the formation of dolomite at near room temperature are subtle except for temperature.This study suggests that microbial mediation is needed for appreciable amount of primary dolomite formation. On the other hand, reactive-transport models depicting dolomitization (temperature range of 40 to 200°C) predicts the formation of two adjacent moving coupled reaction zones (calcite dissolution and dolomite precipitation) with sharp dolomitization front, and generation of >20% of secondary porosity. Due to elevated temperature of formation, dolomitization mechanism is efficient in converting existing calcite into dolomite at a much faster rate compared to primary dolomite formation.     

Permeability generation and resetting of tracers during metamorphic fluid flow: implications for advection-dispersion models

Advection-dispersion fluid flow models implicitly assume that the infiltrating fluid flows through an already fluid-saturated medium. However, whether rocks contain a fluid depends on their reaction history, and whether any initial fluid escapes. The behaviour of different rocks may be illustrated using hypothetical marble compositions. Marbles with diverse chemistries (e.g. calcite + dolomite + quartz) are relatively reactive, and will generally produce a fluid during heating. By contrast, marbles with more restricted chemistries (e.g. calcite + quartz or calcite-only) may not. If the rock is not fluid bearing when fluid infiltration commences, mineralogical reactions may produce a reaction-enhanced permeability in calcite + dolomite + quartz or calcite + quartz, but not in calcite-only marbles. The permeability production controls the pattern of mineralogical, isotopic, and geochemical resetting during fluid flow. Tracers retarded behind the mineralogical fronts will probably be reset as predicted by the advection-dispersion models; however, tracers that are expected to be reset ahead of the mineralogical fronts cannot progress beyond the permeability generating reaction. In the case of very unreactive lithologies (e.g. pure calcite marbles, cherts, and quartzites), the first reaction to affect the rocks may be a metasomatic one ahead of which there is little pervasive resetting of any tracer. Centimetre-scale layering may lead to the formation of self-perpetuating fluid channels in rocks that are not fluid saturated due to the juxtaposition of reactants. Such layered rocks may show patterns of mineralogical resetting that are not predicted by advection-dispersion models. Patterns of mineralogical and isotopic resetting in marbles from a number of terrains, for example: Chillagoe, Marulan South, Reynolds Range (Australia); Adirondack Mountains, Old Woman Mountains, Notch Peak (USA); and Stephen Cross Quarry (Canada) vary as predicted by these models. Received: 3 February 1997 / Accepted: 26 June 1997     

Diagenesis of the Khuff Formation (Permian–Triassic), northern United Arab Emirates

This diagenetic study (including fieldwork, petrographic, fluid inclusion, and stable isotope investigations) deals with the outcrop of Upper Permian–Lower Triassic carbonate rocks, which are equivalent to the Khuff Formation. The studied succession, which outcrops in the Ras Al Khaimah region, northern United Arab Emirates, comprises three formations, including the Bih, the Hagil, and the Ghail formations. The study focuses on unraveling the conditions and fluid compositions encountered during diagenesis of the succession. Emphasize is also made on linking diagenesis to major stratigraphic surfaces and to highlight reservoir property evolution and heterogeneity of the studied rocks. The evolution of fluids and related diagenetic products can be summarized as follows: (1) formation of near-surface to shallow burial, fine-crystalline dolomite (dolomite matrix) through pervasive dolomitization of carbonate sediments by modified marine pore waters; (2) formation of coarse-crystalline dolomite cement by highly evolved marine pore waters (13–23 wt.% NaCl eq.) at elevated temperatures (120–208°C), and (3) calcite cementation by highly saline fluid (20–23 wt.% NaCl eq.) at high temperature (170–212°C). A final calcite cement generation has been formed by the percolation of meteoric fluids during uplift. Fracture- and vug-filling diagenetic minerals are mainly restricted to the mid-Bih breccia marker level, suggesting preferential focused fluid flow through specific stratigraphic surfaces as well as along tectonic-related structures. Reservoir properties have been evolved as result of the interplay of the original sedimentary texture and the diagenetic evolution. Porosity is higher in the Bih Formation, which is dominated by dolomitized packstones and grainstones, than in the Hagil and Ghail formations, consisting mainly of dolomitized mudstones and wackestones. Image analyses were used to quantify the visual porosity in thin sections. The highest porosity values were measured in the Bih Formation, which is characterized by significant amounts of vug- and fracture-filling cements. This feature is attributed to the increase of porosity owing to substantial dissolution of abundant intergranular and vug-filling cements. In contrast, the Hagil and Ghail formations, which consist of finer-grained rock than the Bih Formation, were less cemented, and thus, the porosity enhancement by cement dissolution was insignificant.     

雪峰山西侧重点白云岩储层成岩、储集特征 及质量影响因素

雪峰山西侧地区下古生界中—上寒武统和震旦系灯影组白云岩是潜在的重要油气储集层。作为主要目的层段,本区的白云岩类型主要包括结晶白云岩、残余颗粒白云岩、灰质白云岩和藻白云岩四大类。白云岩储层总体上均属于低—特低孔、低—特低渗储层。影响储层发育的成岩事件包括重结晶作用、白云岩化作用、溶蚀作用、破裂作用等。研究表明,白云岩结晶程度可以影响其初始孔隙的发育程度,具粉晶—细晶结晶程度的白云岩具有最好的孔渗性。白云岩化作用的强弱对储层物性具有较为直接的影响。最佳白云岩储层中主要储集空间包括晶间孔隙(包括剩余晶间孔和晶间溶蚀孔)、次生溶孔、洞和裂缝。总体上,较强的白云岩化作用、一定程度的重结晶、较发育的微裂缝、较好的区域渗透性、晚期压溶缝合线发育和晚期较大构造缝不发育等因素有利于白云岩储层的发育。     

Reaction-induced fluid flow in synthetic quartz-bearing marbles

The reaction kinetics and fluid expulsion during the decarbonation reaction of calcite+quartz=wollastonite+CO₂ in water-absent conditions were experimentally investigated using a Paterson-type gas apparatus. Starting materials consisted of synthetic calcite/quartz rock powders with variable fractions of quartz (10, 20, and 30 wt%) and grain sizes of 10 µm (calcite) and 10 and 30 µm (quartz). Prior to reaction, samples were HIPed at 700 °C and 300 MPa confining pressure and varying pore pressures. Initial porosity was low at 2.7–6.3%, depending on pore pressure during HIP and the amount and grain size of quartz particles. Samples were annealed at reaction temperatures of 900 and 950 °C at 150 and 300 MPa confining pressures, well within the wollastonite stability field. Run durations were between 10 min and 20 h. SEM micrographs of quenched samples show growth of wollastonite rims on quartz grains and CO₂-filled pores between rims and calcite grains and along calcite grain boundaries. Measured widths of wollastonite rims vs. time indicate a parabolic growth law. The reaction is diffusion-controlled and reaction progress and CO₂ production are continuous. Porosity increases rapidly at initial stages of the reaction and attains about 10–12% after a few hours. Permeability at high reaction temperatures is below the detection limit of 10^–21 m² and not affected by increased porosity. This makes persistent pore connectivity improbable, in agreement with observed fluid inclusion trails in form of unconnected pores in SEM micrographs. Release of CO₂ from the sample was measured in a downstream reservoir. The most striking observation is that fluid release is not continuous but occurs episodic and in pulses. Ongoing continuous reaction produces increase in pore pressure, which is, once having attained a critical value (P_crit), spontaneously released. Connectivity of the pore space is short-lived and transient. The resulting cycle includes pore pressure build-up, formation of a local crack network, pore pressure release and crack closure. Using existing models for plastic stretching and decrepitation of pores along with critical stress intensity factors for the calcite matrix and measured pore widths, it results that P_crit is about 20 MPa. Patterns of fluid flow based on mineralogical and stable isotope evidence are commonly predicted using the simplifying assumption of a continuous and constant porosity and permeability during decarbonation of the rock. However, simple flow models, which assume constant pore pressure, constant fluid filled porosity, and constant permeability may not commonly apply. Properties are often transient and it is most likely that fluid flow in a specific reacting rock volume is a short-lived episodic process.Editorial responsibility: J. Hoefs     

中国中西部盆地海相白云岩主要形成机制与模式

塔里木、四川及鄂尔多斯盆地是中国中西部海相白云岩的主要发育地区。三大盆地重点层系海相白云岩新近的成因研究表明,大规模准同生白云岩和埋藏成因白云岩的发育均与蒸发台地密切相关。蒸发台地中由海水浓缩形成的富Mg~(2+)卤水一方面在准同生期,通过蒸发泵和下渗机制交代碳酸钙沉积物而形成与蒸发岩共生的准同生白云岩,另一方面作为富含Mg~(2+)的地层孔隙水,在准同生-浅埋藏期乃至中、深埋藏期,通过侧向渗透、侧向与垂向压实排挤和垂向热对流机制与粗结构的碳酸钙沉积物发生交代反应,在蒸发岩系周边和上下形成广泛分布的埋藏成因白云岩。与热流体作用有关的白云石化主要依靠构造断裂、裂缝、不同级次的层序界面、孔洞层等输导体系发生,分布较局限。热流体云化常表现为对先期白云岩进行叠加改造而形成热水改造白云岩。热流体性质不一,可以是深埋藏混合热水、深部循环水、地幔深部的岩浆热液等。白云石(岩)的生物成因不仅表现为微生物作用导致白云石直接沉淀,还表现为生物的存在与活动为白云石化作用提供Mg~(2+)和云化流体通道。由微生物和宏观藻释放出Mg~(2+),在埋藏期对方解石进行交代是各种富含藻类的灰岩中部分白云石的重要形成机制。生物扰动可明显改善岩石的孔渗性,从而显著促进白云石化作用的发生。     

渗透回流模式白云石化作用对碳酸盐岩储层的影响

白云石化作用对含油气的碳酸盐岩储层有明显的改造作用,而前人的研究大部分是针对白云石的成因模式,未对其形成过程给出具体描述,更未对白云石化作用对储层的改造程度做出定量分析。针对白云石化中的“渗透回流模式”,建立了白云石化作用的一维模型,利用多组分、多相流溶质运移模拟软件TOUGHREACT,研究渗透回流模式对储层物性影响的深度范围。结果表明,当高Mg/Ca比的浓缩海水进入含水层后,导致方解石溶解,白云石沉淀,储层孔隙度与渗透率均比初始值大,模拟5万年时,依靠重力作用,渗透回流白云石化作用能够影响到50 m深的储层范围,产生一定规模的白云岩储层。     

碳酸盐岩溶解-沉淀热力学模型及其在塔北地区的应用

 碳酸盐岩作为一种化学岩类其物性受成岩作用的影响十分明显,溶解作用是否发育直接关系到碳酸盐岩储层质量的好坏。碳酸盐矿物在埋藏环境中的溶解度与地层水的化学成分、环境的温度、pH等有着密切的关系。利用地层水化学资料,根据化学热力学的理论,通过方解石及白云石与地层水作用的平衡反应,计算吉布斯自由能变化ΔG,以此来判断方解石及白云石的溶解趋势,并建立了封闭系统中这两种矿物与水作用的热力学模型。根据这个理论模型,对塔北地区159口井奥陶系及6口井的寒武系地层水进行了ΔG的计算,结果显示,塔北地区碳酸盐岩与地层水反应的ΔG基本为负值,表明该地区总体处于碳酸盐岩溶解的环境,寒武系白云岩反应的自由能小于奥陶系灰岩,即在这种深埋条件下白云岩比灰岩更易溶解,这一点在实际岩芯资料上及前人的实验模拟中均有显示。塔北于奇地区自由能的计算结果也为负值,说明该地区强烈的充填作用并非发生在现今条件下,于奇地区现今的埋藏环境是有利于碳酸盐矿物溶解的,这一结论也在该地区奥陶系及寒武系岩芯上得到了验证。
      

白云石化作用及白云岩储层研究进展

白云石化作用和白云岩储层一直都是碳酸盐岩研究中的重要领域。近年来,随着实验分析技术的进步以及油气勘探的深入,对白云岩的研究也取得了诸多进展:1数值模拟技术逐步应用到白云岩研究中,实现了白云石化模式研究由定性到定量的转变;2对微生物白云石化的研究不断加强,识别出了微生物相关白云石的特殊显微形貌特征并对其生物矿化机制进行了分析;3对已有白云石化模式的重新审视:包括对混合水白云石化的修正、对回流白云石化的扩展以及对构造—热液白云石化模式的丰富;4在白云石化与孔隙相关关系的研究中,突破了白云石化增孔的传统认识,逐步强调白云石化在孔隙保存方面的作用;5注重研究白云岩结构、成岩环境、原始相带以及白云岩形成之后的溶蚀改造等因素对白云岩储层发育的控制作用。在未来的研究中,应加强3个方面内容,一是对白云岩结构演化规律的定量研究,二是在成岩流体示踪方面要加强对新技术手段(如二元同位素、Mg同位素)的使用,同时注意借鉴成矿流体研究中的成熟技术和方法,三是要加强对深部白云岩储层形成和保存机制的探索。     

库车凹陷克拉2气田深层优质储层成因及成岩作用模式

基于成岩作用的影响因素及库车盆地的实际,深入研究了克拉2气田次生孔隙发育段储层岩石成分特征及早期古气候条件下胶结物对于成岩作用及次生孔隙的影响,分析了地层压力系统中超高压带、构造裂缝、自生矿物分带的成因、作用、特点及其与储层物性、成岩环境的相关关系,评价了储层中长石、方解石的溶解,白云石对方解石的交代,压实、胶结等主要成岩作用类型对次生孔隙的具体影响和贡献,认为该区早期碱性成岩环境下方解石胶结物的溶解是形成该区次生孔隙发育带的主导因素,并根据该区自生矿物流体包裹体及埋藏史特征,建立了该区成岩作用及孔隙演化的模式。     

川西中泥盆统观雾山组白云岩储层储集空间类型、成因及演化

川西地区观雾山组白云岩储层储集空间主要为孔洞和裂缝两类.为弄清观雾山组白云岩储层孔洞成因、孔洞充填期次及演化,首先,通过对孔洞型白云岩储层发育规律与沉积相、层序关系的分析,结合第一期白云石胶结物形成与围岩白云石化的先后顺序,认为川西地区观雾山组白云岩储层孔洞为相控准同生岩溶形成;针对孔洞内不同期次白云石和方解石胶结物的包裹体均一温度、碳氧同位素、激光原位U-Pb同位素定年、锶同位素、稀土元素等分析,认为孔洞内胶结物形成于封闭的成岩环境,成岩流体为受下伏碎屑岩地层水加入改造的中泥盆世海水.观雾山组白云岩储层储集空间经历了三个演化阶段:沉积期-白云石化之前的孔洞及裂缝形成阶段、白云石化期间的围岩白云石化及第一期白云石胶结物形成阶段和中—深埋藏成岩期的孔隙定型阶段,其中中—深埋藏阶段是孔隙减少的主要阶段,造成约250的孔隙损失.     

甘肃酒泉盆地青西油田裂缝特征及成因分析

甘肃酒泉盆地青西油田泥砾岩相和泥云岩相储层中发育两种裂缝系统。采用多种方法,系统研究了不同岩相区裂缝的产状特征、电性响应特征和物性特征,并重点分析了控制裂缝形成的成岩演化成因。研究结果表明:泥砾岩相中发育有砾缘缝、粒内溶蚀缝和构造裂缝;泥云岩相中发育层面溶蚀缝、泄水裂缝及构造裂缝。裂缝系统不但对储层渗流能力有较大程度的改善,而且也具有一定程度的储集能力。母岩矿物的选择性（差异）溶蚀是控制本区裂缝和次生孔隙发育的主要因素：泥砾岩成分复杂,各种碎屑物及胶结物之间的交代、溶蚀以及不稳定碎屑颗粒的选择性（差异）溶蚀是砾缘缝、成岩溶蚀缝及大量次生孔隙在泥砾岩储层中形成的主要原因;泥云岩中存在的二元和三元纹层结构从根本上决定了裂缝的发育程度,碳酸钙的白云岩化和沸石化是层面溶蚀缝形成的主要成岩作用。     

鲁西豫东地区奥陶系顶部岩溶储集层特征及有利控制因素*

基于野外露头、岩心、薄片和测录井等地质及地球物理资料分析,并结合区域地质背景,深入研究了鲁西豫东(东濮)地区奥陶系顶部峰峰组和上马家沟组岩溶储集层发育特征、演化过程,并进一步分析优质储集层形成的主控因素,预测了有利储集层分布区。结果表明,奥陶系顶部储集层主要岩性为颗粒石灰岩、云质石灰岩、泥晶石灰岩、粉细晶白云岩、含膏白云岩和藻粘结白云岩等,裂缝和溶洞是主要储集空间类型。白云岩类储集层平均孔隙度为2.5%,平均渗透率为7.43×10^-3μm²;石灰岩类储集层平均孔隙度为2.2%,平均渗透率为2.72×10^-3μm²。储集层经历了沉积—准同生期成孔(寒武纪—中奥陶世)→风化淋滤(晚奥陶世—早石炭世)→矿物充填(晚石炭世—白垩纪末期)→溶蚀改造(古近纪)4个演化阶段。沉积相、成岩作用及构造应力控制了研究区优质储集层的形成。颗粒滩及潮坪相带、断裂发育带、白云石化作用和准同生岩溶、表生岩溶、埋藏岩溶作用6项叠合区储集层质量最优;颗粒滩及潮坪相带、白云石化作用和准同生岩溶作用、表生岩溶作用4项叠合区,以及断裂发育带、表生岩溶作用、埋藏岩溶作用3项叠合区储集层质量次之;颗粒滩及潮坪相带、断裂发育带、白云石化作用和准同生岩溶作用、埋藏岩溶作用5项叠合区储集层质量一般。     

Reactive transport modeling of the interaction between a high-pH plume and a fractured marl: the case of Wellenberg

In the context of the proposed low- and intermediate-level radioactive waste repository at Wellenberg (Switzerland), calculations simulating the interaction between hyperalkaline solutions and a fractured marl, at 25 °C, have been performed. The aim of these calculations is to evaluate the possible effects of mineral dissolution and precipitation on porosity and permeability changes in such a fractured marl, and their impact on repository performance. Solute transport and chemical reaction are considered in both a high-permeability zone (fracture), where advection is important, and the wall rock, where diffusion is the dominant transport mechanism. The mineral reactions are promoted by the interaction between hyperalkaline solutions derived from the degradation of cement (a major component of the engineered barrier system in the repository) and the host rock. Both diffusive/dispersive and advective solute transport are taken into account in the calculations. Mineral reactions are described by kinetic rate laws. The fluid flow system under consideration is a two-dimensional porous medium (marl, 1% porosity), with a high-permeability zone simulating a fracture (10% porosity) crossing the domain. The dimensions of the domain are 6 m per 1 m, and the fracture width is 10 cm. The fluid flow field is updated during the course of the simulations. Permeabilities are updated according to Kozeny's equation. The composition of the solutions entering the domain is derived from modeling studies of the degradation of cement under the conditions at the proposed underground repository at Wellenberg. Two different cases have been considered in the calculations. These 2 cases are representative of 2 different stages in the process of degradation of cement (pH 13.5 and pH 12.5). In both cases, the flow velocity in the fracture diminishes with time, due to a decrease in porosity. This decrease in porosity is caused by the precipitation of calcite (replacement of dolomite by calcite) and other secondary minerals (brucite, sepiolite, analcime, natrolite, tobermorite). However, the decrease in porosity and flow velocity is much more pronounced in the lower pH case. The extent of the zone of mineral alteration along the fracture is also much more limited in the lower pH case. The reduction of porosity in the fractures would be highly beneficial for repository performance, since it would mean that the solutions coming from the repository and potentially carrying radionuclides in solution would have to flow through low-conductive rock before they would be able to get to higher-conductive features. The biggest uncertainty in the reaction rates used in the calculations arises from the surface areas of the primary minerals. Additional calculations making use of smaller surface areas have also been performed. The results show that the smaller surface areas (and therefore smaller reaction rates) result in a smaller reactivity of the system and smaller porosity changes.     

Fluid–rock interactions associated with regional tectonics and basin evolution

An integrated approach consisting of fracture analysis, petrography, carbon, oxygen and strontium‐isotope analyses, as well as fluid‐inclusion micro‐thermometry, led to a better understanding of the evolution of fluid–rock interactions and diagenesis of the Upper Permian to Upper Triassic carbonates of the United Arab Emirates. The deposited carbonates were first marked by extensive early dolomitization. During progressive burial, the carbonates were affected by dolomite recrystallization as well as precipitation of vug and fracture‐filling dolomite, quartz and calcite cements. After considerable burial during the Middle Cretaceous, sub‐vertical north–south oriented fractures (F1) were cemented by dolomite derived from mesosaline to hypersaline fluids. Upon the Late Cretaceous maximum burial and ophiolite obduction, sub‐vertical east–west fractures (F2) were cemented by dolomite (Dc2) and saddle dolomite (Ds) derived from hot, highly saline fluids. Then, minor quartz cement has precipitated in fractures from hydrothermal brines. Fluid‐inclusion analyses of the various diagenetic phases imply the involvement of increasingly hot (200°C) saline brines (20 to 23% NaCl eq.). Through one‐dimensional burial history numerical modelling, the maximum temperatures reached by the studied rocks are estimated to be in the range of 160 to 200°C. Tectonically‐driven flux of hot fluids and associated diagenetic products are interpreted to have initiated during the Late Cretaceous maximum burial and lasted until the Oligocene–Miocene compressional tectonics and related uplift. The circulation of such hydrothermal brines led to partial dissolution of dolomites (Dc2 and Ds) and to precipitation of hydrothermal calcite C1 in new (mainly oriented north–south; F3) and pre‐existing, reactivated fractures. The integration of the obtained data confirms that the diagenetic evolution was controlled primarily by the interplay of the burial thermal evolution of the basin and the regional tectonic history. Hence, this contribution highlights the impacts of regional tectonics and basin history on diagenetic processes, which may subsequently affect reservoir properties.     

深层—超深层白云岩储集层: 机遇与挑战*

深层—超深层白云岩储集层是未来油气勘探的重要领域和目标,是潜在的优质油气储集体或储集层。查明白云岩的成因、孔隙发育与保存机理以及流体性质与演化路径,建立符合实际地质状况的储集层地质模型,对于预测深层—超深层白云岩储集层的分布规律是非常重要的。疏理了白云岩成因的几个核心问题的现状与发展趋势,指出了白云石化作用模式的演变过程和发展方向,提出了建立深层—超深层白云岩多阶段动态演化模式应考虑的问题。中国三大克拉通盆地(四川盆地、塔里木盆地和鄂尔多斯盆地)深层—超深层赋存丰富的白云岩,并且已经有一些重要的油气发现,为开展深层—超深层白云岩油气勘探提供了好的机遇,但深埋藏、多期构造改造、叠加和复杂的流体演化历史也是此项研究重要的挑战。