Experimental modelling of shallow magma emplacement: Application to saucer-shaped intrusions

Saucer-shaped dolerite and sandstone intrusions are common in sedimentary basins world-wide. We have conducted a series of scaled experiments simulating the process of magma emplacement in sedimentary basins, with particular attention on the formation of saucer-shaped sills. The model materials were (1) cohesive fine-grained silica flour, representing brittle crust; and (2) molten low-viscosity oil, representing magma. The experiments were performed in both homogeneous and layered models. In all the experiments, oil injection resulted in doming of the surface. In the homogeneous models, the injected oil formed cone sheets and sub-vertical dykes. Cone sheets formed for shallow injection (1–3 cm), and vertical dykes formed for deeper injection (4–5 cm). In layered models, the injected oil always formed saucer-shaped intrusions. Our experimental results show that (1) sill intrusion results in the formation of a dome, with melt erupting at the rim; (2) layering controls the formation of sills and saucer-shaped sills; (3) saucer-shaped sills are fed from the bottom and the fluid flows upward and outward; and (4) the diameter of saucer-shaped sills increase with increasing emplacement depth. The systematic relation between domes and sills and the depth-dependence of sill diameters show that saucer-shaped intrusions result from the interaction between a growing flat-lying shallow sill and doming of the free surface. We conclude that saucer-shaped intrusions represent fundamental geometries formed by shallow magma intrusion in stratified basins.     

Early deformation fabrics of melange in the Mesozoic Waipapa Terrane, North Island, New Zealand

Abstract The Waipapa Terrane in northern North Island, New Zealand, is a Mesozoic accretion-ary complex formed along the Gondwana margin. It contains abundant melange rocks with distinctive characteristics. Precise analyses of their mesoscopic fabrics in Waiheke Island near Auckland have revealed the following sequence of deformation. The earliest phase of deformation of the sandstone/mudstone association, which is the main constituent of this terrane, originated by chaotic mixing of sand and mud due to liquidization of water-saturated, poorly consolidated sediments. The second phase was characterized by hydrofracturing and subsequent forceful injection of ductile mud into rather brittle sand. Local intrusions of sand forming dykes and sills followed these events, as well as intrusions of pelagic/hemipelagic green argillite originally underlying the sandstone/mudstone association. An abundant occurrence of these mixing and multi-stage injection/intrusion fabrics strongly suggest that the Waipapa Terrane around the study area was a site of high pore-fluid pressure. Scaly-foliated melange fabrics with monoclinic symmetry, originating from layer-parallel shearing, were then locally superimposed on the pre-existing melange fabrics. Similar scaly-foliated fabrics also developed in the chert beds originally intercalated between the green argillite and the uppermost part of the oceanic crust. These scaly fabrics might have been related to the regional stacking and juxtaposition of the accreted sediments. The sequence and variation in style of deformation forming the melange fabrics presumably reflected changes in porosity and state of compaction of accreted sediments in a shallow tectonic level.     

Emplacement mechanisms of sill complexes: Information from the geochemical architecture of the Golden Valley Sill Complex,South Africa

The well-exposed Golden Valley Sill Complex, Karoo basin, South Africa, consists of four large sills (ca. 100 m thick; long axes: 13–24 km), one small sill (55–80 m thick; long axis: 4 km; forming an appendix to one of the large sills), and two large dykes (15–20 m thick; 25 and 70 km long), plus some minor intrusions. Field mapping shows physical connections between the small sill and one of the large sills, but no other connection between the large sills, or connections between the sills and the large dykes.     

Depositional processes and gas pore pressure in pyroclastic flows: an experimental perspective

The depositional processes and gas pore pressure in pyroclastic flows are investigated through scaled experiments on transient, initially fluidized granular flows. The flow structure consists of a sliding head whose basal velocity decreases backwards from the front velocity (U _f) until onset of deposition occurs, which marks transition to the flow body where the basal deposit grows continuously. The flows propagate in a fluid-inertial regime despite formation of the deposit. Their head generates underpressure proportional to U _f ² whereas their body generates overpressure whose values suggest that pore pressure diffuses during emplacement. Complementary experiments on defluidizing static columns prove that the concept of pore pressure diffusion is relevant for gas-particle mixtures and allow characterization of the diffusion timescale (t _d) as a function of the material properties. Initial material expansion increases the diffusion time compared with the nonexpanded state, suggesting that pore pressure is self-generated during compaction. Application to pyroclastic flows gives minimum diffusion timescales of seconds to tens of minutes, depending principally on the flow height and permeability. This study also helps to reconcile the concepts of en masse and progressive deposition of pyroclastic flow units or discrete pulses. Onset of deposition, whose causes deserve further investigation, is the most critical parameter for determining the structure of the deposits. Even if sedimentation is fundamentally continuous, it is proposed that late onset of deposition and rapid aggradation in relatively thin flows can generate deposits that are almost snapshots of the flow structure. In this context, deposition can be considered as occurring en masse, though not strictly instantaneously.     

Quality of pTRM acquisition in pyrrhotite bearing contact-metamorphic limestones: possibility of a continuous record of Earth magnetic field variations

The rockmagnetic and palaeomagnetic signal in pyrrhotite bearing limestones of different contact-metamorphic settings were investigated related to intrusions ranging from small sills to large magmatic complexes. Magnetic susceptibility, the pyrrhotite/magnetite ratio and thermal modelling serve as an investigative tool to define three metamorphic zonations: a contact zone of a mixed magnetic assemblage and low susceptibilities, a pyrrhotite bearing transitional zone, where full thermoremanent magnetisations (TRMs) are acquired due to temperatures above the Curie temperature of pyrrhotite (T_c,po), and a marginal zone containing pyrrhotite and magnetite generated at temperatures below T_c,po. The fact that TRMs can consist of independent pTRMs is successfully tested by modified Thellier experiments. It is shown that a metamorphic environment with low fluid circulation provides a scenario for the recording of independent pTRMs. Multicomponent behaviour of the NRM residing in samples from the transitional zone can therefore be interpreted as a continuous record of Earth magnetic field variations.     

Role of fluids in surface deformation caused by the 1999 Chi‐Chi earthquake in Taiwan

The 1999 Chi‐Chi earthquake significantly altered the landscape of central Taiwan. Surface deformation produced by the earthquake along the trace of the Chelungpu thrust can be classified into two styles: (1) uplift without significant surface rupture, and (2) uplift accompanied by surface rupture. Here we examine areas that exhibited the first style of deformation (e.g. Wufeng). Seismic stress at the time of the main shock may have been relieved by high pore‐fluid pressure in a 300‐m‐thick sand and gravel aquifer. Along the thrust fault, frictional heating of these sediments resulted in thermal expansion and an increase in pore‐fluid pressure. High pore‐fluid pressure damped seismic‐wave energy and enhanced intergranular slips of unconsolidated sandy and gravel sediments, which were possibly assisted by sulphuric acid corrosion, leading to a high sulphate content in the groundwater (c. 70 mg L^?1). These changes permitted surface folding and terrace‐style uplifting to occur without significant rupture. In contrast, other areas in which the second style of deformation is dominant (e.g. Fengyuen‐Shihkang) have thin (0–10 m) sand and gravel deposits and lower concentrations of sulphate (c. 30 mg L^?1) in groundwater. In these areas, sediments were heated but not sufficiently to produce significant thermal expansion and increase in pore‐fluid pressure; accumulation of stress in these locations led to rupture at the ground surface, with the formation of steep fault scarps. The areas exhibiting the first deformation style are characterized by the presence of high pore‐fluid pressure, frictional heat conduction, and possibly chemical corrosion related to sulphuric acid attack and formation of sulphate, in contrast to those involving significant uplift and surface rupture. The areal distribution of these two surface deformation styles suggests that the aforementioned fluid‐related subsurface processes may have altered the characteristics of sediments and caused diverse responses to the quake. Copyright © 2002 John Wiley & Sons, Ltd.     

Determining magma flow in sills,dykes and laccoliths and their implications for sill emplacement mechanisms

Three-dimensional seismic data from the Faeroe-Shetland Basin provides detailed information on the relationships between sills, dykes, laccoliths and contemporaneous volcanic activity. The data shows that sills are predominantly concave upwards, being complete or partial versions of radially or bilaterally symmetrical forms that possess flat inner saucers connected to a flat outer rim by a steeply inclined sheet. Such morphologies are only partially modified by pre-existing faults. Sills can be sourced from dykes or the steep climbing portions of deeper sills. Both sills and dykes can provide magma to overlying volcanic fissures and sills can be shown to feed shallow laccoliths. Magma flow patterns, as revealed by opacity rendering, suggest that sills propagate upwards and outwards away from the magma feeder. As an individual sill can consist of several leaves emplaced at different stratigraphic levels, and as a sill or dyke can provide magma to volcanic fissures, other sills and laccoliths, the data suggests that neutral buoyancy concepts may not provide a complete explanation for the mechanism and level of sill emplacement. Instead, the data suggests that the presence of lithological contrasts, particularly ductile horizons such as overpressured shales may permit sill formation at any level below the neutrally buoyant level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Ken Thomson–deceased, April 2007     

Early Miocene to Quaternary evolution of volcanism and the basin formation in western Anatolia: a review

Western Anatolia, largely affected by extensional tectonics, witnessed widespread volcanic activity since the Early Miocene. The volcanic vents of the region are represented by epicontinental calderas, stratovolcanoes and monogenetic vents which are associated with small-scale intrusions as sills and dykes. The volcanic activity began with an explosive character producing a large ignimbritic plateau all over the region, indicating the initiation of the crustal extension event. These rhyolitic magmas are nearly contemporaneous with granitic intrusions in western Anatolia. The ignimbrites, emplaced approximately contemporaneous with alluvial fan and braided river deposits, flowed over the basement rocks prior to extensional basin formation. The lacustrine deposits overlie the ignimbrites. The potassic and ultrapotassic lavas with lamprophyric affinities were emplaced during the Late Miocene–Pliocene. The volcanic activities have continued with alkali basalts during the Quaternary.     

Forward induced seismic hazard assessment: application to a synthetic seismicity catalogue from hydraulic stimulation modelling

The M _w 3.2-induced seismic event in 2006 due to fluid injection at the Basel geothermal site in Switzerland was the starting point for an ongoing discussion in Europe on the potential risk of hydraulic stimulation in general. In particular, further development of mitigation strategies of induced seismic events of economic concern became a hot topic in geosciences and geoengineering. Here, we present a workflow to assess the hazard of induced seismicity in terms of occurrence rate of induced seismic events. The workflow is called Forward Induced Seismic Hazard Assessment (FISHA) as it combines the results of forward hydromechanical-numerical models with methods of time-dependent probabilistic seismic hazard assessment. To exemplify FISHA, we use simulations of four different fluid injection types with various injection parameters, i.e. injection rate, duration and style of injection. The hydromechanical-numerical model applied in this study represents a geothermal reservoir with preexisting fractures where a routine of viscous fluid flow in porous media is implemented from which flow and pressure driven failures of rock matrix and preexisting fractures are simulated, and corresponding seismic moment magnitudes are computed. The resulting synthetic catalogues of induced seismicity, including event location, occurrence time and magnitude, are used to calibrate the magnitude completeness M _c and the parameters a and b of the frequency-magnitude relation. These are used to estimate the time-dependent occurrence rate of induced seismic events for each fluid injection scenario. In contrast to other mitigation strategies that rely on real-time data or already obtained catalogues, we can perform various synthetic experiments with the same initial conditions. Thus, the advantage of FISHA is that it can quantify hazard from numerical experiments and recommend a priori a stimulation type that lowers the occurrence rate of induced seismic events. The FISHA workflow is rather general and not limited to the hydromechanical-numerical model used in this study and can therefore be applied to other fluid injection models.     

Natural hydraulic cracking: numerical model and sensitivity study

Natural hydrofracturing caused by overpressure plays an important role in geopressure evolution and hydrocarbon migration in petroliferous basins. Its mechanism is quite well understood in the case of artificial hydraulic fracturing triggered by high-pressure fluid injection in a well. This is not so for natural hydraulic fracturing which is assumed to initiate as micro-cracks with large influence on the permeability of the medium. The mechanism of natural hydraulic cracking, triggered by increasing pore pressure during geological periods, is studied using a fracturing model coupled to the physical processes occurring during basin evolution. In this model, the hydraulic cracking threshold is assumed to lie between the classical failure limit and the beginning of dilatancy. Fluid pressure evolution is calculated iteratively in order to allow dynamic adjustment of permeability so that the fracturing limit is always preserved. The increase of permeability is interpreted on the basis of equivalent fractures. It is found that fracturing is very efficient to keep a stress level at the rock’s hydraulic cracking limit: a fracture permeability one order of magnitude larger than the intrinsic permeability of the rock would be enough. Observations reported from actual basins and model results strongly suggest that natural hydraulic cracking occurs continuously to keep the pressure at the fracturing limit under relaxed stress conditions.     

Impact of Poroelastic Response of Sandstones on Geothermal Power Production

During geothermal power production using a borehole doublet consisting of a production and injection well, the reservoir conditions such as permeability k, porosity φ and Skempton coefficient B at the geothermal research site Gross Schoenebeck/Germany will change. Besides a temperature decrease at the injection well and a change of the chemical equilibrium, also the pore pressure p _p will vary in a range of approximately 44 MPa ± 10 MPa in our reservoir at ?3850 to ?4258 m depth. This leads to a poroelastic response of the reservoir rocks depending on effective pressure p _eff (difference between mean stress and pore pressure), resulting in a change in permeability k, porosity φ and the poroelastic parameter Skempton coefficient B. Hence, we investigated the effective pressure dependency of Flechtinger sandstone, an outcropping equivalent of the reservoir rock via laboratory experiments. The permeability decreased by 21% at an effective pressure range from 3 to 30 MPa, the porosity decreased by 11% (p _eff = 6 to 65 MPa) and the Skempton coefficient decreased by 24% (p _eff = 4 to 25 MPa). We will show which mechanisms lead to the change of the mentioned hydraulic and poroelastic parameters and the influence of these changes on the productivity of the reservoir. The most significant changes occur at low effective pressures until 15 to 20 MPa. For our in situ reservoir conditions p _eff = 43 MPa a change of 10 MPa effective pressure will result in a change in matrix permeability of less than 4% and in matrix porosity of less than 2%. Besides natural fracture systems, fault zones and induced hydraulic fractures, the rock matrix its only one part of geothermal systems. All components can be influenced by pressure, temperature and chemical reactions. Therefore, the determined small poroelastic response of rock matrix does not significantly influence the sustainability of the geothermal reservoir.     

Fluid Pressure Variation in a Sedimentary Geothermal Reservoir in the North German Basin: Case Study Groß Schönebeck

The Rotliegend of the North German basin is the target reservoir of an interdisciplinary investigation program to develop a technology for the generation of geothermal electricity from low-enthalpy reservoirs. An in situ downhole laboratory was established in the 4.3 km deep well Groβ Schönebeck with the purpose of developing appropriate stimulation methods to increase permeability of deep aquifers by enhancing or creating secondary porosity and flow paths. The goal is to learn how to enhance the inflow performance of a well from a variety of rock types in low permeable geothermal reservoirs. A change in effective stress due to fluid pressure was observed to be one of the key parameters influencing flow properties both downhole and in laboratory experiments on reservoir rocks. Fluid pressure variation was induced using proppant-gel-frac techniques as well as waterfrac techniques in several different new experiments in the borehole. A pressure step test indicates generation and extension of multiple fractures with closure pressures between 6 and 8.4 MPa above formation pressure. In a 24-hour production test 859 m³ water was produced from depth indicating an increase of productivity in comparison with former tests. Different depth sections and transmissibility values were observed in the borehole depending on fluid pressure. In addition, laboratory experiments were performed on core samples from the sandstone reservoir under uniaxial strain conditions, i.e., no lateral strain, constant axial load. The experiments on the borehole and the laboratory scale were realized on the same rock types under comparable stress conditions with similar pore pressure variations. Nevertheless, stress dependences of permeability are not easy to compare from scale to scale. Laboratory investigations reflect permeability variations due to microstructural heterogeneities and the behavior in the borehole is dominated by the generation of connections to large-scale structural patterns.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Small scale arcuate intrusions on Saint Helena,South Atlantic

Sixteen arcuate intrusions have been emplaced at extremely high levels into the basaltic shield volcanoes of Saint Helena. These intrusions are of special interest because of their small size and modes of emplacement. The arcuate masses are of three distinctive types:

1. Irregular, steeply inward-dipping, basic sheets with diameters of 150 m to 450 m infill tensional fractures originating at depths of about 500 m beneath the volcano surface.
2. Strongly curved sheets from 25 m to 750 m in diameter are cross-sections of inclined intrusions which in three dimensions resemble single sticks of celery. These intrusions, varying from basalt to trachyte in composition, are infilled tensional fractures originating at « point » pressure sources, inclined to the horizontal, at depths of about 500 m below the surface.
3. Salic intrusions with near-vertical sides and gently inclined roof-infillings have outer diameters of 350 m to 1070 m. Pressure exerted onto the flanks of the volcano by domed, convex upwards, areas of a magma chamber roof, at a depth of about 2 km, caused near-vertical ring fractures to form. Formation of a sub-horizontal cross fracture and subsequent intrusion of magma produced the « roof-infillings » by updoming the overlying basalts or sinking of the enclosed block, or combinations of the two processes. Two intrusions of this third type are multiple.

     

门源M6.9地震诱发地质灾害特征研究

2022年1月8日发生的门源M6.9地震诱发了崩塌、滑坡、砂土液化、地裂缝等多种同震地质灾害。通过对门源M6.9地震地质灾害进行现场调查,得出了地质灾害的分布特征和各类型地质灾害的主要特点,分析了地震地质灾害不发育的原因,并对地震地质灾害的长期效应进行了分析预测。研究结果表明：门源地震诱发地质灾害主要分布在震中附近;崩塌、落石总体规模较小,滑坡多为岩质滑坡,且以冰碛物和表层岩土体的溜滑为主。受表层土体冻结和孔隙水压力消散的影响,饱和砂土液化沿较窄的地裂缝呈串珠状分布,喷出物多为粉细砂。地震形成了4条左旋左阶斜列的地表破裂带,并在极震区内形成了大量的地裂缝。断层破碎带对地震动的阻隔作用、覆盖层薄、地表土冻结可能是造成本次地震地质灾害总体不发育的主要原因;地震产生的大量地裂缝导致斜坡和堆积体的稳定性减弱,在耦合集中降雨、冻融作用等因素后可能诱发滑坡灾害,松散堆积于沟床处的崩滑物作为物源,可能会增加地震影响区泥石流灾害的风险。     

Secondary welding of submarine,pumice-lithic breccia at Mount Chalmers,Queensland, Australia

Very thick units of massive pumice and lithic clast-rich breccia in the Early Permian Berserker beds at Mount Chalmers, Queensland, are deposits from cold, water-supported, volcaniclastic mass flows emplaced in a below-wave base submarine setting. Adjacent to syn-volcanic andesitic and rhyolitic sills and dykes, the pumice-lithic breccia shows a well-developed eutaxitic texture. The eutaxitic foliation is parallel to intrusive contacts and extends as far as a few metres away from the contact. At these sites, pumice clasts are strongly flattened and tube vesicles within the pumice clasts are compacted and aligned parallel to the direction of flattening. Some lenticular pumice clasts contain small (2 mm), round, quartz-filled amygdales and spherulites. Further away from the sills and dykes, the pumice clasts have randomly oriented, delicate tube vesicle structure and are blocky or lensoid in shape. Round amygdales were generated by re-vesiculation of the glass and the spherulites indicate devitrification of the glass at relatively high temperatures. The eutaxitic texture is therefore attributed to re-heating and welding compaction of glassy pumice-lithic breccia close to contacts with intrusions. In cases involving sills, secondary welding along the contacts formed extensive, conformable, eutaxitic zones in the pumice-lithic breccia that could be mistaken for primary welding compaction in a hot, primary pyroclastic deposit.     

Computed tomography scan imaging of natural open fractures in a porous rock; geometry and fluid flow

Computed tomography scan imaging techniques have been used on core samples to investigate the effect of natural open fractures on reservoir flow in the Snøhvit Gas Condensate Field. Firstly, computed tomography (CT) scanning was used to describe the 3D geometrical properties of the fracture network including orientation and fracture density. Two types of fractures were observed: F1 fractures are short and stylolite related and F2 fractures are longer, cross-cutting the core and without any obvious relationship to stylolites. Secondly, monitoring of single and two phase flow experiments on samples containing these two types of natural open fractures was performed under 10 and 80 bar net confining pressure while using CT scanning. 1-phase miscible flooding experiment shows approximately 3 times higher flooding velocity in an open F2 fracture than in the matrix. 2-phase flooding by gas injection into a 100% water saturated core gave early gas breakthrough due to flow in the fracture and thereafter very little water production. The flow experiments showed that the presence of open fractures has a significant local effect on fluid flow even in a case with relatively high matrix porosity (200–300 mD). The sample containing F1 fractures showed a complex flow pattern influenced both by open fractures and stylolites. The CT scan data enables an exact representation of the fracture network in core scale simulation models and therefore improves the understanding of fracture influence on flow in a fractured porous medium. CT scanning of core samples provides an effective tool for integrating geology and fluid flow properties of a porous fractured medium.     

Are open fractures necessarily aligned with maximum horizontal stress?

Fluid flow in fractured rock is an increasingly central issue in recovering water and hydrocarbon supplies and geothermal energy, in predicting flow of pollutants underground, in engineering structures, and in understanding large-scale crustal behaviour. Conventional wisdom assumes that fluids prefer to flow along fractures oriented parallel or nearly parallel to modern-day maximum horizontal compressive stress, or S_Hmax. The reasoning is that these fractures have the lowest normal stresses across them and therefore provide the least resistance to flow. For example, this view governs how geophysicists design and interpret seismic experiments to probe fracture fluid pathways in the deep subsurface. Contrary to these widely held views, here we use core, stress measurement, and fluid flow data to show that S_Hmax does not necessarily coincide with the direction of open natural fractures in the subsurface (>3 km depth). Consequently, in situ stress direction cannot be considered to predict or control the direction of maximum permeability in rock. Where effective stress is compressive and fractures are expected to be closed, chemical alteration dictates location of open conduits, either preserving or destroying fracture flow pathways no matter their orientation.     

Shear and tension hydraulic fractures in low permeability rocks

Laboratory hydrofracture experiments were performed on triaxially stressed specimens of oil shale and low-permeability granite. The results show that either shear or tension fractures could develop depending on the level of differentials stress, even in specimens containing preexisting fractures. With 1 kb of confining pressure and differential stress greater than 2kb, hydraulic fluid diffusion into the specimens reduced the effective confining pressure until failure occurred by shear fracture. Below 2kb of differential stress, tension fractures occurred. These results suggest that hydraulic fracturing in regions of significant tectonic stress may produce shear rather than tension fractures. In this casein situ stress determinations based on presumed tension fractures would lead to erroneous results.     

Numerical modeling of formation damage by two-phase particulate transport processes during CO2 injection in deep heterogeneous porous media

Prediction of CO₂ injection performance in deep subsurface porous media relies on the ability of the well to maintain high flow rates of carbon dioxide during several decades typically without fracturing the host formation or damaging the well. Dynamics of solid particulate suspensions in permeable media are recognized as one major factor leading to injection well plugging in sandstones. The invading supercritical liquid-like fluid can contain exogenous fine suspensions or endogenous particles generated in situ by physical and chemical interactions or hydrodynamic release mechanisms. Suspended solids can plug the pores possibly leading to formation damage and permeability reduction in the vicinity of the injector. In this study we developed a finite volume simulator to predict the injectivity decline near CO₂ injection wells and also for production wells in the context of enhanced oil recovery. The numerical model solves a system of two coupled sets of finite volume equations corresponding to the pressure-saturation two-phase flow, and a second subsystem of solute and particle convection-diffusion equations. Particle transport equations are subject to mechanistic rate laws of colloidal, hydrodynamic release from pore surfaces, blocking in pore bodies and pore throats, and interphase particle transfer. The model was validated against available laboratory experiments at the core scale. Example results reveal that lower CO₂ residual saturation and formation porosity enhance CO₂-wet particle mobility and clogging around sinks and production wells. We conclude from more realistic simulations with heterogeneous permeability spanning several orders of magnitude that the control mode of mobilization, capture of particles, and permeability reduction processes strongly depends on the type of permeability distribution and connectivity between injection and production wells.