巨厚含水层注浆浆液扩散半径及影响因素研究——以峰峰矿区为例

总结了饱水裂隙岩体中浆液的球形扩散模型和柱状扩散模型。通过对不同注浆模型中各参数,如浆液黏度、岩石孔隙率、注浆孔半径、岩体渗透系数及启动压力等对注浆扩散半径的影响分析,认为注浆孔半径和岩石孔隙率对宾汉姆流型浆液扩散影响不明显,而影响浆液扩散半径的主要因素是浆液黏度、岩体渗透系数及启动压力梯度。最后针对峰峰矿区,通过计算提出了包括注浆压力和注浆结束压力的注浆参数理论值。     

Helium measurements of pore fluids obtained from the San Andreas Fault Observatory at Depth (SAFOD, USA) drill cores

⁴He accumulated in fluids is a well established geochemical tracer used to study crustal fluid dynamics. Direct fluid samples are not always collectable; therefore, a method to extract rare gases from matrix fluids of whole rocks by diffusion has been adapted. Helium was measured on matrix fluids extracted from sandstones and mudstones recovered during the San Andreas Fault Observatory at Depth (SAFOD) drilling in California, USA. Samples were typically collected as subcores or from drillcore fragments. Helium concentration and isotope ratios were measured 4?C6 times on each sample, and indicate a bulk ⁴He diffusion coefficient of 3.5?±?1.3?×?10^?C8 cm²?s^?C1 at 21°C, compared to previously published diffusion coefficients of 1.2?×?10^?C18 cm²?s^?C1 (21°C) to 3.0?×?10^?C15 cm²?s^?C1 (150°C) in the sands and clays. Correcting the diffusion coefficient of ⁴He_water for matrix porosity (??3%) and tortuosity (??6?C13) produces effective diffusion coefficients of 1?×?10^?C8 cm^2?s^?C1 (21°C) and 1?×?10^?C7 (120°C), effectively isolating pore fluid ⁴He from the ⁴He contained in the rock matrix. Model calculations indicate that <6% of helium initially dissolved in pore fluids was lost during the sampling process. Complete and quantitative extraction of the pore fluids provide minimum in situ porosity values for sandstones 2.8?±?0.4% (SD, n?=?4) and mudstones 3.1?±?0.8% (SD, n?=?4).     

A numerical study of mineral alteration and self-sealing efficiency of a caprock for CO2 geological storage

Geochemical interactions of brine–rock–gas have a significant impact on the stability and integrity of the caprock for long-term CO₂ geological storage. Invasion of CO₂ into the caprock from the storage reservoir by (1) molecular diffusion of dissolved CO₂, (2) CO₂-water two-phase flow after capillary breakthrough, and (3) CO₂ flow through existing open fractures may alter the mineralogy, porosity, and mechanical strength of the caprock due to the mineral dissolution or precipitation. This determines the self-enhancement or self-sealing efficiency of the caprock. In this paper, two types of caprock, a clay-rich shale and a mudstone, are considered for the modeling analyses of the self-sealing and self-enhancement phenomena. The clay-rich shale taken from the Jianghan Basin of China is used as the base-case model. The results are compared with a mudstone caprock which is compositionally very different than the clay-rich shale. We focus on mineral alterations induced by the invasion of CO₂, feedback on medium properties such as porosity, and the self-sealing efficiency of the caprock. A number of sensitivity simulations are performed using the multiphase reactive transport code TOUGHREACT to identify the major minerals that have an impact on the caprock’s self-sealing efficiency. Our model results indicate that under the same hydrogeological conditions, the mudstone is more suitable to be used as a caprock. The sealing distances are barely different in the two types of caprock, both being about 0.6 m far from the interface between the reservoir and caprock. However, the times of occurrence of sealing are considerably different. For the mudstone model, the self-sealing occurs at the beginning of simulation, while for the clay-rich shale model, the porosity begins to decline only after 100 years. At the bottom of the clay-rich shale column, the porosity declines to 0.034, while that of mudstone declines to 0.02. The sensitive minerals in the clay-rich shale model are calcite, magnesite, and smectite-Ca. Anhydrite and illite provide Ca²⁺ and Mg²⁺ to the sensitive minerals for their precipitation. The mudstone model simulation is divided into three stages. There are different governing minerals in different stages, and the effect of the reservoir formation water on the alteration of sensitive minerals is significant.     

Modeling of Cs+ diffusion and retention in the DI-A2 experiment (Mont Terri). Uncertainties in sorption and diffusion parameters

In the DI-A2 experiment several non-reactive and reactive tracers were injected as a pulse in a packed-off borehole in the Opalinus Clay. Unlike the previous DI-A1 test, the design of the Teflon filter in the injection borehole forced the water to flow through the filter and the open space between the filter and the borehole wall (the filter itself did not act as a diffusion barrier between the circulating solution and the rock). The decrease in tracer concentration in the liquid phase was monitored during a period of a year. Afterwards, the borehole section was overcored and the tracer profiles in the rock were analyzed. A main interest of this experiment was to understand the chemical behavior of sorbing tracers: Cs⁺ (stable), ⁸⁵Sr²⁺, ⁶⁰Co²⁺ and Eu³⁺ (stable). The complete dataset (except for Eu³⁺ because of strong sorption to experimental equipment) was analyzed in a previous study with a 2D diffusion–reaction model and the derived diffusion and sorption parameters were compared with laboratory data. As in DI-A1, a difference by a factor of about 2 for sorption (magnitude of the Freundlich isotherm) was obtained between in situ and laboratory batch sorption experiments.Recent experimental and modeling studies have shown equivalent Cs⁺ sorption on intact and disaggregated Opalinus Clay samples. In view of these developments, new modeling of Cs⁺ diffusion and retention in the DI-A2 experiment has been performed using CrunchFlow. The calculations include transport by diffusion and a multisite cation exchange model to account for the retention of Cs⁺. The new results show that upscaling of Cs⁺ sorption from laboratory to field is no longer required. However, a difference in sorption by a factor of about 2 is still explained by the use of different versions of the same cation exchange model (a small difference in the selectivity coefficient for one type of site). This uncertainty in sorption leads to an uncertainty in the effective diffusion coefficient (D_e) for Cs⁺, also by a factor of 2 (2–4 × 10⁻¹⁰ m²/s). Clearly, the values of D_e obtained are correlated with the strength of sorption in the model, with stronger sorption leading to larger D_e values. Discrimination between the two versions of the exchange model is not possible when using only the results of the in situ test. Additionally, during early times (t < 10 days) the drop in Cs⁺ concentration in the circulation system is slower than expected. Due to the experimental setup, this slow decrease in concentration cannot be caused by the filter in the contact between borehole and rock. Poor mixing in the circulation system could explain this effect.     

巴东组典型红层软岩单轴压缩试验细观组构特征

红色粉砂质泥岩和泥质粉砂岩是巴东组特殊性岩土的代表,具有遇水易崩解特性。根据岩样单轴压缩试验,基于PFC_2D程序,建立中硬岩、软岩单轴压缩试验数值模型;模型选用平行黏结颗粒接触模型。根据应力-应变曲线特征,采用分段定量参数标定法对软岩试样的细观参数进行标定;且通过调整法向和切向黏结强度比值σ_c/τ_c,控制试样单轴压缩破坏模式;分别对不同风化程度的粉砂质泥岩和泥质粉砂岩在饱和、天然状态下单轴压缩试验过程进行模拟;并分析了试验过程中颗粒法向接触力、切向接触力、配位数、孔隙度等细观组构参数的分布特征和演化规律。研究结果表明:该数值方法能够很好地模拟中硬岩、软岩的单轴压缩试验过程,除中风化粉砂质泥岩外,风化程度和浸水条件仅能影响试样颗粒接触力在统计角度范围的大小,不能影响其分布形式;同时,风化程度和浸水条件延长了配位数显著下降的加载时间,配位数随裂缝数量的增加而减小;浸水条件和风化程度对泥质粉砂岩的孔隙度演化影响较大;风化程度对粉砂质泥岩孔隙度演化影响较大,而浸水条件影响较小。      

Active and passive measurements of radon diffusion coefficient from building construction materials

The gaseous state and chemical inert behavior of radon make it important tracer for the radon transport study through the building materials. The radon resistant property of building construction materials is important parameter to control the indoor radon levels in living and workplaces. The materials with higher radium content can be made less severe by the use of some building materials of low diffusion coefficient and diffusion length. This makes the study of radon diffusion through building material more important along with the study of exhalation and radioactivity content. Keeping this in mind the radon diffusion study was carried out through different building construction materials used for wall and floor by active and passive techniques. The diffusion coefficient from these building materials measured by passive methods varied from (0.9 ± 0.5) × 10^?7 to (22.95 ± 13.19) × 10^?6 m²s^?1 and radon diffusion length varied from 0.21 to 3.31 m for cement, soil, sand, wall putty and plaster of Paris (POP) etc. The radon diffusion coefficient measured by active technique varied from 1.93 × 10^?10 to 1.36 × 10^?7 m²s^?1 for samples with definite geometry like paper, polyethylene, marble, granite etc. The radon diffusion coefficient and diffusion length depend upon the porosity and density of materials for powder samples.     

钻井液流态对孔壁稳定性的影响

沉积地层中孔壁的稳定性与钻井液流态有着密切关系。结合河南叶舞凹陷盆地2 386 m盐矿钻探工程,采用不同的泵量和上返流速在不同孔段进行了钻探试验。试验和理论计算表明:在水敏性地层（黏土、泥岩、砂质泥岩等）浅层岩石孔隙率和塑性指数较高时,在钻井液层流流态情况下,主要以缩径为主,最大缩径率达26.6%;中部和深部同等地层,在紊流流态情况下,缩径和超径同时存在。其中,随着孔深的增加,缩径率呈减小趋势,超径率则有增大趋势,泥岩地层最大超径率达157.4%,盐岩地层超径率高达225.7%。      

Anion exclusion effects in compacted bentonites: Towards a better understanding of anion diffusion

Diffusion of ³⁶Cl⁻ in compacted bentonite was studied using through-diffusion, out-diffusion and profile analysis techniques. Both the bulk dry density of the bentonite and the composition of the external solution were varied. Increasing the bulk dry density of the bentonite resulted in a decrease of both the effective diffusion coefficient and the Cl-accessible porosity. Increasing the ionic strength of the external solutions resulted in an increase of both the effective diffusion coefficient and the Cl-accessible porosity. This can be explained by anion exclusion effects (Donnan exclusion). At high ionic strength values (I ? 1 M NaCl) the Cl-accessible porosity approaches the interparticle porosity. This interparticle porosity is the difference between the total and interlayer porosity of the bentonite. The interlayer porosity was found to depend on the degree of compaction. Up to a bulk dry density of 1300 kg m⁻³ the interlayer is built up of 3 water layers. Between 1300 and 1800 kg m⁻³ the interlayer water is reduced from 3 to 2 layers of water. Above 1800 kg m⁻³ evidence for a further decrease to 1 layer of water was found. These findings are in agreement with X-ray data found in the literature showing a decrease of the basal spacing of montmorillonite (the main clay mineral in bentonite) with increasing degree of compaction. The relationship between the effective diffusion coefficient of Cl⁻ and the diffusion-accessible porosity can be described by an empirical relationship analogous to Archie’s law. To predict the effective diffusion coefficient of Cl⁻ in compacted bentonite, the diffusion coefficient of Cl⁻ in water, the bulk dry density and the ionic strength of the pore water have to be known.     

A double medium model for diffusion in fluid-bearing rock

The concept of a double porosity medium to model fluid flow in fractured rock has been applied to model diffusion in rock containing a small amount of a continuous fluid phase that surrounds small volume elements of the solid matrix. The model quantifies the relative role of diffusion in the fluid and solid phases of the rock. The fluid is the fast diffusion path, but the solid contains the volumetrically significant amount of the diffusing species. The double medium model consists of two coupled differential equations. One equation is the diffusion equation for the fluid concentration; it contains a source term for change in the average concentration of the diffusing species in the solid matrix. The second equation represents the assumption that the change in average concentration in a solid element is proportional to the difference between the average concentration in the solid and the concentration in the fluid times the solid-fluid partition coefficient. The double medium model is shown to apply to laboratory data on iron diffusion in fluid-bearing dunite and to measured oxygen isotope ratios at marble-metagranite contacts. In both examples, concentration profiles are calculated for diffusion taking place at constant temperature, where a boundary value changes suddenly and is subsequently held constant. Knowledge of solid diffusivities can set a lower bound to the length of time over which diffusion occurs, but only the product of effective fluid diffusivity and time is constrained for times longer than the characteristic solid diffusion time. The double medium results approach a local, grain-scale equilibrium model for times that are large relative to the time constant for solid diffusion.     

Borehole stress monitoring analysis on advanced abutment pressure induced by Longwall Mining

The paper presented the research on the dynamic advanced abutment stress induced by longwall mining with borehole stress meters on mining side coal mass. Twenty vibrating wire borehole stress meters were installed into the extracting coal mass wall of a first mining roadway of 910 m depth in Zhuji Coal Mine, China, and were used to monitor dynamic changes in vertical and horizontal stresses. Three months of continuous monitoring and further analysis showed that the impacting distance of advanced abutment stress induced by mining in the strike of the working face along its central axis was the farthest, greater than 250 m (the face length is 220 m); it gradually decreased in the radial direction of the face from its central axis outward; the pressure peak was located within 24 m in the front of the mining coal wall; non-synchronous caving of the layered mudstone roof at the stope occurred. Comparison between vertical and horizontal stress increments indicated that the horizontal stress was much smaller than the vertical stress in the coal mass of mining side, while the latter’s magnitude determined the drastic degree of mine pressure manifestation. The study has been applied to determine the advanced support length of the working face and further provide a reliable basis to forecast such dynamic disasters as rock burst, coal and gas outburst, etc., as well as to design the asymmetric supports on both sides of a gateway.     

横观各向同性水敏性地层斜井眼坍塌压力确定

在以泥页岩为目标储层的非常规油气藏水平井开发中, 钻遇井壁的稳定性问题突出, 一方面是由于在设计钻井液密度时简单采用各向同性介质模型和Mohr-Coulomb强度准则, 另一方面则是由于选用的钻井液性能欠佳, 难以保证建井全过程的井壁稳定。为此, 将地层视为横观各向同性介质, 采用Mogi-Coulomb强度准则判别岩石基体和层理面的稳定性, 同时考虑泥页岩水化作用对岩石力学强度参数的影响, 建立了水敏性地层井壁稳定分析模型, 分析了坍塌压力随井斜、方位、钻井液性能、钻井时间等因素的变化规律。研究结果表明:对于具有显著层理面结构的泥页岩地层, 采用横观各向同性模型比各向同性介质模型更能描述强度的各向异性;采用Mogi-Coulomb准则比Mohr-Coulomb强度准则得到的坍塌压力值与实际吻合程度更高, 在保证井壁稳定的前提下有利于降低钻井液密度, 实现钻井提速;提高钻井液封堵能力、减小泥页岩吸水扩散能力, 有利于延长坍塌周期。利用东濮、威远、焦石坝等工区多口页岩气水平井的实钻资料对本文方法进行了验证, 证实本文方法确定的钻井液密度窗口下值更能满足安全钻井需求。     

Emission reduction technology and application research of surface borehole methane drainage in coal mining-influenced region

Windblown methane is an important gas resulting in atmospheric greenhouse effect. Therefore, reduction in windblown methane is one of the important measures to mitigate atmospheric greenhouse effect. In China, weak coal seam of low permeability is common in coal mines, so it is beset with difficulties to decrease the methane emission rate by means of gas drainage from the virgin coal seam, further to decrease the windblown mine gas. Utilizing the pressure relief and permeability and fluidity improvement effect in coal mining an approach to release methane emission through surface borehole was established, for example establishing a comprehensive deformation fracture model of surface borehole in extraction area based on quantitative rule of overlying rock movement in pit and forming a technology to select the surface borehole arrangement site in extraction area on the basis of deformation of bore frame structure and distribution characteristics of extraction flow field. And optimization technology of shape and structure of surface borehole in extraction area on the basis of ultimate stress analysis of surface borehole casing was given. The technology overcomes effectively the problem that surface borehole casing is vulnerable to premature fracture due to impact of strata movement on the surface borehole, and further increases the drainage result of the surface borehole. The technology has been test in China Shanxi Jincheng Sihe Coal Mine, achieving good results, including 12,000 Nm³/d pure methane drainage rate from single borehole, 85% methane concentration and 1.1 million Nm³ accumulative methane drainage, which demonstrate practicability and advanced performance of the technology.     

Geochemistry of groundwater from a rhyolite aquifer,Northwest Iran

Chemical data on groundwater composition in rhyolitic hard rock aquifers with limited global occurrence are rarely found. In this research geochemistry of Mahabad Rhyolite Aquifer, NW Iran, was studied considering major ions, silica and trace elements measured in wet and dry seasons. Based on the results, the mean silica content was 18 mg l^?1, less than the average of the rhyolitic waters. However, the relatively higher electrical conductivity (EC) of 418 µS cm^?1 was measured. Based on a PHREEQCI model, the weathering of the silicate minerals and dissolution of carbonated intercalations turns groundwater dominantly into Ca–HCO₃ type, enhancing EC, pH and silica concentration along the flow path. Trace elements of Sr, Ba and Pb were measured at highest concentrations, the later with an average value of 83 ppb exceeds the drinking guidelines. Cluster analysis confirms biotite weathering and barite dissolution as the main sources of the trace elements in the groundwater. The results signify geochemical features of rhyolitic groundwater which can be a useful tracer of mixing in flow systems containing variety of aquifers including rhyolites.     

High-resolution vertical profiles of groundwater electrical conductivity (EC) and chloride from direct-push EC logs

Elevated groundwater salinity associated with produced water, leaching from landfills or secondary salinity can degrade arable soils and potable water resources. Direct-push electrical conductivity (EC) profiling enables rapid, relatively inexpensive, high-resolution in-situ measurements of subsurface salinity, without requiring core collection or installation of groundwater wells. However, because the direct-push tool measures the bulk EC of both solid and liquid phases (EC_a), incorporation of EC_a data into regional or historical groundwater data sets requires the prediction of pore water EC (EC_w) or chloride (Cl^?) concentrations from measured EC_a. Statistical linear regression and physically based models for predicting EC_w and Cl^? from EC_a profiles were tested on a brine plume in central Saskatchewan, Canada. A linear relationship between EC_a/EC_w and porosity was more accurate for predicting EC_w and Cl^? concentrations than a power-law relationship (Archie’s Law). Despite clay contents of up to 96%, the addition of terms to account for electrical conductance in the solid phase did not improve model predictions. In the absence of porosity data, statistical linear regression models adequately predicted EC_w and Cl^? concentrations from direct-push EC_a profiles (EC_w = 5.48 EC_a + 0.78, R ² = 0.87; Cl^? = 1,978 EC_a – 1,398, R ² = 0.73). These statistical models can be used to predict EC_w in the absence of lithologic data and will be particularly useful for initial site assessments. The more accurate linear physically based model can be used to predict EC_w and Cl^? as porosity data become available and the site-specific EC_w–Cl^? relationship is determined.     

Dynamic process of the Wenjiagou rock landslide in Sichuan Province,China

A rock avalanche is a geological event that is always sudden, rapid and with a long run-out, and can result in large loss of lives and property. The Wenjiagou rock avalanche was a high-speed rock landslide caused by a strong earthquake, in Mianzhu, Sichuan Province, southwest China. In this study, we reproduce the movement and deposition processes of the sliding mass by numerical simulation. We analyze the effects of the friction coefficient of each slip surface and the strength of the parallel bonds and contact stiffness between particles on the dynamic process and deposit features using three-dimensional particle flow code (PFC3D). The simulation results agree with the field measurements when the friction coefficient is 0.2, parallel bond strength is 2 MPa, and contact stiffness is 2?×?10⁸ kN/m. The landslide lasted about 115 s from the initial movement to the final deposition at the exit of the valley. The maximum velocity of the sliding mass was 114 m/s.     

Diffusion of municipal waste contaminants in compacted lateritic soil treated with bagasse ash

Compacted clay can minimize infiltration of liquid into waste or control the release of contaminated liquids to the surrounding soils and groundwater. Compacted lateritic soil treated with up to 12 % bagasse ash and municipal solid waste (MSW) leachate sourced from a domestic waste land fill were used in diffusion test studies to access the diffusion characteristics of some inorganic species present in the municipal solid waste leachate. Diffusion set-up were prepared containing 0, 4, 8 and 12 % bagasse ash—soil mixes compacted at 2 % wet of optimum using the modified proctor effort. The set up was saturated with water for 30 days before the introduction of MSW leachate and initiation of diffusion test for another 90 days. After diffusion testing, water content within the soil column showed a decrease with depth. Diffusion test results generally showed that diffusion is an active means of transport of chemical species even at very low flow rates in the compacted soil-bagasse ash mixes, and the effective diffusion coefficient is affected by bagasse ash. The pore fluid concentration profile for the various chemical species tested showed that the compacted soil-bagasse ash mix has the capacity to attenuate Ca²⁺, Pb²⁺ and Cr³⁺ ions.     

Availability of coastal groundwater discharge as an alternative water resource in a large-scale reclaimed land,Korea

This study was conducted to identify the availability of coastal groundwater discharge (CGD), subsurface fluids flowing from inland through the coastal area to sea, as an alternative water resource for a large-scale reclaimed land. The behaviors of stable isotopes indicated that groundwater originated from inland precipitation and traveled as CGD along the coast line. Most of the groundwater samples collected from domestic wells installed along the old coast line were considered to be relatively fresh from the correlation analysis among chemical constituents. The average electrical conductivity (EC) values of the samples were identified as averaging 1,125–1,297 μS cm^?1, corresponding to appropriate crop growth. A weathered-rock layer in a small catchment within the reclaimed land was proved to be a main CGD pathway, with electrical resistivity anomalies ranging from 7 to 14 Ω m. Five monitoring wells were placed in this catchment to delineate the occurrence of CGD. Long-term vertical EC profiling results for the monitoring wells indicated that CGD occurs within a depth of 30 m below the ground surface. Annual monitoring data for groundwater level and EC demonstrated that the water quality of CGD was improved by introducing fresh terrestrial groundwater. A remarkable improvement in water quality (EC decrease of 900–1,600 μS cm^?1) of CGD was observed during the saline water pumping test that explains how CGD could be an alternative water resource for the reclaimed land.     

Effect of confining pressure on the diffusion of HTO, 36Cl− and 125I− in a layered argillaceous rock (Opalinus Clay): diffusion perpendicular to the fabric

The through- and out-diffusion of HTO, ³⁶Cl⁻ and ¹²⁵I⁻ in Opalinus Clay, an argillaceous rock from the northern part of Switzerland, was studied under different confining pressures between 4 and 15 MPa. The direction of diffusion and the confining pressure were perpendicular to the bedding. Confining pressure had only a small effect on diffusion. An increase in pressure from 4 to 15 MPa resulted in a decrease of the effective diffusion coefficient of ∼20%. Diffusion accessible porosities were not measurably affected. The values of the effective diffusion coefficients, D_e, ranged between (5.6±0.4)×10⁻¹² and (6.7±0.4)×10⁻¹² m² s⁻¹ for HTO, (7.1±0.5)×10⁻¹³ and (9.1±0.6)×10⁻¹³ m² s⁻¹ for ³⁶Cl⁻ and (4.5±0.3)×10⁻¹³ and (6.6±0.4)×10⁻¹³ m² s⁻¹ for ¹²⁵I⁻. The rock capacity factors, α, measured were circa 0.14 for HTO, 0.040 for ³⁶Cl⁻ and 0.080 for ¹²⁵I⁻. Because of anion exclusion effects, anions diffuse slower and exhibit smaller diffusion accessible porosities than the uncharged HTO. Unlike ³⁶Cl⁻, ¹²⁵I⁻ sorbs weakly on Opalinus Clay resulting in a larger rock capacity factor. The sorption coefficient, K_d, for ¹²⁵I⁻ is of the order of 1–2×10⁻⁵ m³ kg⁻¹. The effective diffusion coefficient for HTO is in good agreement with values measured in other sedimentary rocks and can be related to the porosity using Archie's Law with exponent m=2.5.     

Recovery of groundwater in the Sanriku region contaminated by the tsunami inundation from the 2011 Tohoku earthquake

Many wells in the Sanriku region used as sources for water supply systems were heavily contaminated by the tsunami of the 2011 great Tohoku earthquake on March 11 in 2011. To better understand the nature of the groundwater contamination by the tsunami inundation and to clarify the recovery process of contaminated groundwater at the study wells, groundwater monitoring has been conducted once or twice yearly since early summer in 2011. High and abnormal values of electric conductivity (EC), chloride ion concentration (CIC), Na⁺, Ca⁺, heavy metal ions, and heavier isotopes of the contaminated groundwater were also obtained in April and June 2011. The chemical elements have rapidly and exponentially decreased as a result of effective pumping of the contaminated groundwater from the study wells and because of abundant rainfall in 2011. In April 2015 (about 4 years after the tsunami inundation), the CIC and EC of the contaminated groundwater of two study wells in Minamisanriku town had reached pre-inundation values. The estimated residence times of groundwater of the two study wells were 105–118 days in the full-day pumping stage and 910–1000 days in the daytime-only pumping stage.     

Paragenesis of Cretaceous to Eocene carbonate reservoirs in the Ionian fold and thrust belt (Albania): relation between tectonism and fluid flow

ABSTRACT This paper examines the diagenetic history of dual (i.e. matrix and fracture) porosity reservoir lithologies in Cretaceous to Eocene carbonate turbidites of the Ionian fold and thrust belt, close to the oil‐producing centre of Fier–Ballsh (central Albania). The first major diagenetic event controlling reservoir quality was early cementation by isopachous and syntaxial low‐Mg calcite. These cements formed primarily around crinoid and rudist fragments, which acted as nucleation sites. In sediments in which these bioclasts are the major rock constituent, this cement can make up 30% of the rock volume, resulting in low effective porosity. In strata in which these bioclasts are mixed with reworkedmicrite, isopachous/syntaxial cements stabilized the framework, and matrixporosity is around 15%. The volumetric importance of these cements, their optical and luminescence character (distribution and dull orange luminescence) and stable isotopic signal (δ¹⁸O and δ¹³C averaging respectively; ?0·5‰ VPDB and +2‰ VPDB) all support a marine phreatic origin. Within these turbidites and debris flows, several generations of fractures alternated with episodes of cementation. A detailed reconstruction of this history was based on cross‐cutting relationships of fractures and compactional and layer‐parallel shortening (LPS) stylolites. The prefolding calcite veins possess orange cathodoluminescence similar to that of the host rock. Their stable isotope signatures (δ¹⁸O of ?3·86 to ?0·85‰ VPDB and δ¹³C of – 0·14 to + 2·98‰ VPDB) support a closed diagenetic rock‐buffered system. A similar closed system accounts for the selectively reopened and subsequently calcite‐cemented LPS stylolites (δ¹⁸O of ?1·81 to ?1·14‰ VPDB and δ¹³C of +1·52 to +2·56‰ VPDB). Within the prefolding veins, brecciated host rock fragments and complex textures such as crack and seal features resulted from hydraulic fracturing. They reflect expulsion of overpressured fluids within the footwall of the frontal thrusts. After folding and thrust sheet emplacement, some calcite veins are still rock buffered (δ¹⁸O of ?0·96 to +0·2‰ VPDB and δ¹³C of +0·79 to +1·37‰ VPDB), whereas others reflect external (i.e. extraformational) and thus large‐scale fluid fluxes. Some of these veins are linked to basement‐derived fluid circulation or originated from fluid flow along evaporitic décollement horizons (δ¹⁸O around +3·0‰ VPDB and δ¹³C around +1·5‰ VPDB). Others are related to the maturation of hydrocarbons in the system (δ¹⁸O around ?7·1‰ VPDB and δ¹³C around +9·3‰ VPDB). An open joint system reflecting an extensional stress regime developed during or after the final folding stage. This joint system enhanced vertical connectivity. This open joint network can be explained by the high palaeotopographical position and the folding of the reservoir analogue within the deformational front. The joint system is pre‐Burdigalian in age based upon a dated karstified discordance contact. Sediment‐filled karst cavity development is linked to meteoric water infiltration during emergence of some of the structures. Despite its sediment fill, the karst network is locally an important contributor to reservoir matrix porosity in otherwise tight lithologies. Development of secondary porosity along bed‐parallel and bed‐perpendicular (i.e. layer‐parallel shortening) stylolites is interpreted as a late‐stage diagenetic event associated with migration of acidic fluids during hydrocarbon maturation. Development of porosity along the LPS system enhanced the vertical reservoir connectivity.