Seepage System Research Articles

Development and application of a thermo-hydro-mechanical coupling test system for a single rock fracture

A thermo-hydro-mechanical coupling test system was developed to study the physical and mechanical properties of a single rock fracture under the interaction of multi-physics field. The test system mainly contained a loading system, a shear-flow box and sealing system, a high-temperature environment box and a seepage system. The test system can conduct the normal compression and shear-flow coupling test under complex stress paths within 300 °C. Two load cells were set at both ends of the shear direction to solve problem of inaccurate shear force monitoring. The unique shear-flow box was designed to solve the dynamic sealing problem. To verify the reliability of the test system, the seepage tests during heating–cooling process and a shear-flow coupling test for a single granite fracture under various boundary conditions were carried out. The test system provides research equipment for understanding the interaction mechanism of multi-physical field in deep underground engineering.

Inversion Model for Permeability Coefficient Based on Random Forest–Secretary Bird Optimization Algorithm: Case Study of Lower Reservoir of C-Pumped Storage Power Station

The geological complexity of the karst regions presents significant challenges, with the permeability coefficient being a critical parameter for accurately analyzing seepage behavior in hydraulic engineering projects. To overcome the limitations of traditional inversion methods, which often exhibit low computational efficiency, poor accuracy, and instability, this study utilizes a finite-element forward model and orthogonal experimental design to establish a sample set for permeability-coefficient inversion. A surrogate model for seepage calculation based on the Random Forest (RF) algorithm is subsequently developed. Furthermore, the Secretary Bird Optimization Algorithm (SBOA) is incorporated to propose an intelligent RF–SBOA inversion method for permeability-coefficient estimation, which is validated through a case study of the C-pumped storage power station. The results demonstrate that the RF model’s predictions for water levels at four boreholes closely align with the measured data, outperforming models such as CART, BP, and SVR. The SBOA effectively identifies the optimal geological permeability coefficient, with the borehole water-level inversion achieving a maximum relative error of only 0.128%, which meets the accuracy requirements for engineering applications. Additionally, the computed distribution of the natural seepage field is consistent with the typical distribution patterns observed in mountain seepage systems. During the normal water-storage phase, both the calculated seepage flow and gradient comply with engineering standards, while the seepage-field distribution aligns with empirical observations. This inversion model provides a rapid and accurate method for estimating the permeability coefficient of strata in the project area, with potential applicability to permeability inversion in other engineering geology contexts, thus demonstrating considerable practical value for engineering applications.

GIS-Based Mapping of the Water Quality and Geochemical Assessment of the Ionic Behavior in the Groundwater Aquifers of Middle Ganga Basin, Patna, India

The study implemented Geographic Information System (GIS) techniques and multivariate hydrogeochemical analysis to evaluate the spatial-temporal and seasonal variation in the groundwater quality of Patna, India. For this purpose, sixty groundwater samples were collected and analyzed for major anions and cations during the pre-monsoon, monsoon, and post-monsoon seasons of 2019-2020. The physicochemical parameters such as pH, EC (Electrical Conductivity), TDS (Total Dissolved Solids), TH (Total Hardness), Ca2+, Mg2+, Na+, K+, HCO3-, Cl-, SO42- were considered to evaluate the water quality index. The result revealed degradation in groundwater quality from pre-monsoon (49.21) to post-monsoon (74.48). EC, TDS, TH, Mg2+, Na+, Ca2+, K+, and HCO3- ions were found accountable for high WQI values at various sampling sites during different seasons. Spatial maps showed that 45 % of the sampling stations exhibited poor quality in all three seasons, where the eastern part of the studied region was revealed to be the most affected area. The application of multivariate statistical methods and hydrogeochemical investigation has clearly defined the dominant role of the weathering process, and reverse ion exchange mechanism in controlling the aquifer’s ionic chemistry. Moreover, poor seepage system, and waste leachate from the surface have been found as the main cause of high levels of Na+, K+, and Cl- in the eastern part of Patna.

Pore-scale insights into CO2-water two-phase flow and implications for benefits of geological carbon storage

The overall benefits of geological carbon storage (GCS) depend primarily on CO2 storability and injectability, expressed as saturation and relative permeability, respectively. The effects of GCS schemes on these two properties, the macroscopic response indicators of a two-phase seepage system, are closely related to pore-scale two-phase behaviors. However, the comprehensive effects of capillary number (Ca) and wettability (θ) on saturation and relative permeability are poorly understood. Here we proposed a digital rock physics (DRP) technique workflow for the phase field method and systematically investigated how these effects control two-phase seepage at pore scale through the high-resolution visualization results obtained. We created a Ca-θ phase diagram identified by four pore-scale displacement mechanisms, including finger-like invasion, burst, cooperative filling and coexistence of concave and convex interfaces, to illustrate the comprehensive effects of Ca and θ. We found that the relative permeability of the defending phase (water in this work) is determined by the net effect of the direct driving and viscous coupling effects. We organized comprehensive Ca-θ diagrams and revealed the favorable conditions for CO2 injectability and storability. Our results demonstrate that GCS schemes, mainly about capillary number and wettability, can significantly influence CO2 storage performance via the two-phase flow at pore scale, which should be considered carefully. This work provides valuable insights into the selection of an optimal GCS scheme and contributes to an in-depth understanding of multiphase seepage at pore scale.

Nonlinear evolution characteristics and seepage mechanical model of fluids in broken rock mass based on the bifurcation theory

With the deep extension of coal mining in China, fault water inrush has become one of the major disasters threatening the safety production of coal mine. Based on the control equations of steady state and non-Darcy seepage in fractured rock mass, the multi-parameter nonlinear dynamic seepage equations of fractured rock mass are established in this paper. Based on the nonlinear dynamics theory, the function of the state variable in the system is derived, and the influence of the gradual change of non-Darcy flow factors on the structural stability of seepage system is studied. The research achievements show that there are three branches in the equilibrium state of the seepage system. Specifically, the stability of the equilibrium state changes abruptly near the limit parameter. The seepage dynamic system of fractured rock mass has the delayed bifurcation, and the coal mine disaster such as fault water inrush occurs easily at the bifurcation point. The research results are of great significance to enrich the theory of fault water inrush in coal mine, and to reveal the disastrous mechanism of fault water inrush and guide its prevention and control technology in coal mine, which can provide the theoretical reference for predicting the water seepage stability in fractured rock mass.

Dynamic Pressure Analysis of Shale Gas Wells Considering Three-Dimensional Distribution and Properties of the Hydraulic Fracture Network

A major challenge in transient pressure analysis for shale gas wells is their complex transient flow behavior and fracturing parameters. While numerical simulations offer high accuracy, analytical models are attractive for transient pressure analysis due to their high computational efficiency and broad applicability. However, traditional analytical models are often oversimplified, making it difficult to capture the complex seepage system, and three-dimensional fracture characteristics are seldom considered. To address these limitations, this study presents a comprehensive hybrid model that characterizes the transient flow behavior and analyzes the pressure response of a fractured shale gas well with a three-dimensional discrete fracture. To achieve this, the hydraulic fracture is discretized into several panels, and the transient flow equation is numerically solved using the finite difference method. Based on the Langmuir adsorption isotherm and the pseudo-steady diffusion in matrix and Darcy flow in the network of micro-fractures, a reservoir model is established, and the Laplace transformation is adopted to solve the model analytically. The transient responses are obtained by dynamically coupling the flow in the reservoir and the discrete fracture. The precision of the proposed model is validated using the commercial numerical simulator, Eclipse. A series of transient pressure dynamic curves are drawn to make a precise observation of different flow regimes, and the effects of several parameters on transient pressure response are also examined. The results show that the shale gas well testing interpretation curves comprise nine flow stages. The pressure drop of shale gas reservoirs is lower than that of conventional gas reservoirs due to the replenishment of desorbed gas. The artificial fracture flow capacity, fracture length, and height are the main engineering factors affecting the pressure responses of shale gas wells. Maximizing the degree and scope of reconstruction can enhance the gas well production capacity during fracturing construction. The research results also indicate that our model is a reliable semi-analytical model for well test interpretations in real case studies.

Local Seismicity and Sediment Deformation in the West Svalbard Margin: Implications of Neotectonics for Seafloor Seepage

AbstractIn the Fram Strait, mid‐ocean ridge spreading is represented by the ultra‐slow system of the Molloy Ridge, the Molloy Transform Fault and the Knipovich Ridge. Sediments on oceanic and continental crust are gas charged and there are several locations with documented seafloor seepage. Sedimentary faulting shows recent stress release in the sub‐surface, but the drivers of stress change and its influence on fluid flow are not entirely understood. We present here the results of an 11‐month‐long ocean bottom seismometer survey conducted over the highly faulted sediment drift northwards from the Knipovich Ridge to monitor seismicity and infer the regional state of stress. We obtain a detailed earthquake catalog that improves the spatial resolution of mid‐ocean ridge seismicity compared with published data. Seismicity at the Molloy Transform Fault is occurring southwards from the bathymetric imprint of the fault, as supported by a seismic profile. Earthquakes in the northern termination of the Knipovich Ridge extend eastwards from the ridge valley, which together with syn‐rift faulting identified in seismic reflection data, suggests that a portion of the currently active spreading center is buried under sediments away from the bathymetric expression of the rift valley. This hints at the direct link between crustal rifting processes and faulting in shallow sediments. Two earthquakes occur close to the seepage system of the Vestnesa Ridge further north from the network. We suggest that deeper rift structures, reactivated by gravity and/or post‐glacial subsidence, may lead to accommodation of stress through shallow extensional faults, therefore impacting seepage dynamics.

Characterizing Håkon Mosby Mud Volcano (Barents Sea) cold seep systems by combining ROV-based acoustic data and underwater photogrammetry

Cold-seeps have a unique geo-ecological significance in the deep-sea environment. They impact the variability of present-day submarine sedimentary environments, affecting the evolution of the landscape over time and creating a variety of submarine landforms, one of which is Mud Volcanoes (MVs). MVs form due to the extrusion of mud, fluids, and gas, mainly methane, from deeper sedimentary layers. These natural gas seepage systems could significantly affect climate change and the global carbon cycle. We present a comprehensive method that combines ROV-based multibeam mapping and underwater photogrammetry to enhance the understanding of the physical relationships between geomorphic units characterizing the Håkon Mosby Mud Volcano (HMMV) and the distribution of associated habitats. HMMV is indeed characterized by high thermal and geochemical gradients from its center to the margins resulting in a clear zonation of chemosynthetic communities. Our approach integrates multi-resolutions and multi-sources data acquired using a work-class ROV. The ROV-based microbathymetry data helped to identify the different types of fine-scale submarine landforms in the central part of HMMV. This revealed three distinct geomorphic units, with the central hummocky region being the most complex. To further study this area, ROV images were analyzed using a defined Structure from Motion workflow producing millimetric resolution 2D and 3D models. Object-Based Image Analysis (OBIA), applied on orthomosaics, allowed us to obtain a fine classification of main benthic communities covering a total area of 940m2, including the active seepage area of the hummocky rim. Four major substrate types were identified in these regions: uncovered mud, bacterial mats high-density, bacterial mats low-density, sediments and tubeworms. Their relationship with terrain morphology and seepage activity were investigated at different spatial scales, contributing to a deeper understanding the ecological functioning of cold seep ecosystems in MVs. The applied workflow is proposed as an innovative processing technique for future studies on cold-seep systems. Geomorphic and ecological processes in extreme environments are inherently linked and marked by spatial patterns typifying associated habitats and sedimentary environments. This is poorly investigated in previous studies, leaving a substantial gap in the geomorphological drivers responsible for habitat distribution and extent in cold seep systems.

Fractured Zone Detection Using Petrophysical Logs by Rescaled Windowed R/S Analyses and Discrete Fracture Network (DFN) of Yanchang Formation in Mahuangshan-Hongliugou Area, Ordos Basin, China

Structural fractures generally develop in the upper crust strata and are usually distributed in a convergent pattern, forming structural fracture zones with a specific strike. Fracture zones control the reservoir seepage system and seriously affect the migration and accumulation of oil and gas in tight sandstone reservoirs. Therefore, characterizing the characteristics of the fracture zones for tight oil exploration and development effectively is essential. In this paper, the variable scale fractal method is introduced to calculate the petrophysical log, and a new curve H is built. An intensity log is to characterize the intensity of structural fracture development. The H curve is in a good linear relationship with the intensity curve after the comparison of the H curve and intensity curve in 32 wells. A quantitative relationship between H and the intensity curve is established. Based on the parameters obtained from the core and image logs, the discrete fracture network model was established using H curves from more than 300 wells, and the structural fracture zone was analyzed. The model shows that the fracture zones formed by structural fractures are in S-N and NW-SE directions in the study area. The orientation of the structural fracture zone is consistent with that of the fractured fault zone and fault, and the development of the fractured zone is consistent with the regional tectonic evolution characteristics. The characteristics of the fracture zone explain the distribution law of oil accumulation and groundwater salinity in the study area.

Applications of remote sensing and near-surface geophysical methods for the investigations of a morainic ridge susceptible to landslide. A case study from the high lands of the Himalayas, Pakistan

Numerous factors are necessary to comprehensively comprehend the complex phenomenon of landslides and, ultimately, facilitate effective measures for their protection and management. This paper investigates the case history of the Maicher landslide in the Nagar District (Northern Pakistan) by using a multidisciplinary combined technique approach based on conventional geophysical measurements and modern remote sensing techniques. The first technique is to identify the rate of deformations with real-time landslide monitoring sensors and satellite images. Care is devoted to data acquisition and processing using the ad-hoc methodology to minimize errors and ensure high reliability. For the second technique (electrical resistivity tomography, surface waves measurements, and ground-penetrating radar), convenient methods for landslide characterization are used and herein presented. The combination of information taken from the different techniques allows the results to be validated by conventional and modern approaches and enhances the accuracy of each method. The results suggest that shear surface planes at a depth of 35–40 m have developed over time. The analysis also provided information on controlling factors and possible expansion of the different slope sections. This work highlighted that geophysical parameters (shear wave velocity Vs, dielectric constant ε, and resistivity value ρ) could provide reliable qualitative and quantitative information about the shallow layer, bedrock contacts, possible sliding surface and subsurface water saturation, and seepage system of a landslide. The geophysical results reflect an overview of the geological model and are effective in terms of the slope stabilization process.

Authigenic Carbonate Precipitation at Yam Seep Controlled by Continuous Fracturing and Uplifting of Four‐Way Closure Ridge Offshore SW Taiwan

AbstractHydrocarbon seeps are common manifestations of gas leakage from the seafloor. However, the fate of methane seepage within the gas hydrate stability zone at active margins is poorly constrained. This study presents a 40‐thousand‐year record of hydrocarbon seepage archived by a ∼5‐m long core composed of authigenic carbonate from the Yam Seep area, Four‐Way Closure Ridge off SW Taiwan. Different carbonate microfacies could be distinguished: Consolidated microcrystalline aragonite representing lithified host sediments intercalated by pure aragonite present in 10–50 cm thick intervals in the core. These aragonite intervals are interpreted as having precipitated within former fractures in the host rock. High resolution U‐Th dating of these aragonites is interpreted to record the minimum age of the opening of these fractures. The chronology of aragonite precipitation throughout the core suggests a record of continuous seepage from ∼41 to 2 ka that fluctuated in intensity over this time period. The chronology of putative fracturing events and observed carbonate precipitation suggest (a) an active period of fracturing and seepage from ∼37 to 27 ka, (b) a more quiescent period from ∼27 to 16 ka, followed by (c) another active period from ∼16 to 12 ka. A schematic model illustrates the evolution of carbonate formation within the core influenced by faulting, fracturing, erosion, gas hydrate accumulation, and aragonite precipitation and provides a unique 40,000‐year‐old record of methane seepage and crucial insights into the dynamics of long‐term seepage systems at active margins.

Study on the Difference in CO2 and Air Injection in the Fracture/Matrix Dual Medium of Continental Shale Reservoirs.

To clarify the impact of different displacement media on the enhanced oil recovery of continental shale and realize the efficient and reasonable development of shale reservoirs, this paper takes the continental shale of the Lucaogou Formation in the Jimusar Sag in the Junggar Basin (China, Xinjiang) as the research object and uses real cores to build the fracture/matrix dual-medium model. Computerized tomography (CT) scanning is used to visually compare and analyze the influence of fracture/matrix dual-medium and single-matrix medium seepage systems on oil production characteristics and clarify the difference between air and CO2 in enhancing the oil recovery of continental shale reservoirs. Through a comprehensive analysis of the production parameters, the whole oil displacement process can be divided into three stages: the oil-rich and gas-poor stage, oil and gas coproduction stage, and gas-rich and oil-poor stage. Shale oil production follows the rule of fractures first and matrix second. However, for CO2 injection, after the crude oil in the fractures is recovered, the oil in the matrix migrates to the fractures under the action of CO2 dissolution and extraction. Overall, the oil displacement effect of CO2 is better than that of air, resulting in a 5.42% higher final recovery factor. Additionally, fractures can increase the permeability of the reservoir, which can greatly enhance oil recovery in the early oil displacement process. However, as the amount of injected gas increases, its impact gradually decreases, and ultimately, it is consistent with the recovery of nonfractured shale, which can achieve nearly the same development effect.

Magnetic properties of gas hydrate-bearing sediments and their association with iron geochemistry in the Sea of Marmara, Turkey

The anaerobic oxidation of methane, a key geochemical process that is involved in the cycling of sulfate and iron (oxyhydr)oxides in marine sediments, results in the formation of iron sulfides. Although ferrimagnetic iron sulfides have been identified in seepage systems, the link between iron migration and sediment magnetic properties remains poorly understood. Here, we investigate two cores from the Sea of Marmara to evaluate biogeochemical iron cycling and iron sulfide mineralogy in gas hydrate-bearing sediments. Magnetic analyses indicate the presence of greigite and pyrrhotite in a core from a hydrate-rich site with a high hydrocarbon flux, which contrasts with a lack of these minerals in a core characterized by only mild seepage. This is supported by the results of rock magnetic and scanning electron microscope analyses of the sediments. The presence of authigenic greigite is critical for assessing local redox records and together with the occurrence of monoclinic pyrrhotite may suggest specific diagenetic processes in gas hydrate environments. Our analysis demonstrates the usefulness of these ferrimagnetic minerals, with a high saturation isothermal remanent magnetization to magnetic susceptibility ratio (SIRM/χ > 15 kAm−1) and a high index of hysteresis parameters (DJH > 0.2) indicative of magnetic mineralogy changes, for evaluating variability in the intensity of seepage fluxes and for estimating gas hydrate distributions.

CAGE18-5 Cruise report : Remotely-operated vehicle (ROV) investigations of active gas seepage sites in the Barents Sea

Cruise CAGE18-5 was the first research expedition under the helm of UiT The Arctic University of Norway in Tromsø with the new ice-going research vessel R/V Kronprins Håkon. This new vessel provides new opportunities to collect cross-disciplinary data for addressing the objectives of the Norwegian Centre of Excellence for Arctic Gas Hydrate, Environment and Climate, CAGE. CAGE investigates Arctic gas hydrate and methane seepage systems in order to better understand the effects they may have on our oceans, ecosystems and global climate. The new research vessel and its facilities allows CAGE access to a state-of-art remotely-operated vehicle, the ROV Ægir 6000, opening a new domain for experimental work and acquisition of sample material from the seafloor. The overall goal of cruise CAGE18-5 therefore was to utilize the ROV in order to provide guided video imagery and to study active gas seepage systems at the gas-hydrate pingo (GHP) site located in the outer Storfjord Trough, and at Storbanken in the northeastern part of the Barents Sea. The cruise may be known as: CAGE18_5

Impact of Gas Saturation and Gas Column Height at the Base of the Gas Hydrate Stability Zone on Fracturing and Seepage at Vestnesa Ridge, West-Svalbard Margin

The Vestnesa Ridge, located off the west Svalbard margin, is a >60 km long ridge consisting of fine-grained sediments that host a deep-marine gas hydrate and associated seepage system. Geological and geophysical observations indicate the predominance of vertical fluid expulsion through fractures with pockmarks expressed on the seafloor along the entire ridge. However, despite the apparent evidence for an extended free gas zone (FGZ) below the base of the gas hydrate stability zone (BGHSZ), present-day seafloor seepage has been confirmed only on the eastern half of the sedimentary ridge. In this study, we combine the relationships between aqueous phase pressure, capillary pressure, sediment clay fraction, porosity, and total stress to simulate how much gas is required to open preexisting fractures from the BGHSZ towards the seafloor. Data from four specific sites with different lithology and pressure regime along the ridge are used to constrain the simulations. Results demonstrate that fracturing is favored from the FGZ (with gas saturations < 0.1 and gas column heights < 15 m) towards the seafloor. Neglecting the capillary pressure overpredicts the size of the gas column by up to 10 times, leading to erroneous maximum gas vent volume predictions and associated ocean biosphere consequences. Further parametric analyses indicate that variations in the regional stress regime have the potential to modify the fracture criterion, thus driving the differences in venting across the ridge. Our results are in line with independent geophysical observations and petroleum system modeling in the study area, adding confidence to the proposed approach and highlighting the importance of the capillary pressure influence on gas pressure.

Damage constitutive model of gas-bearing coal using industrial CT scanning technology

Determining the deformation characteristics and mechanical behavior of gas-bearing coal is beneficial to understanding the phenomena and rules of dynamic disasters during coal mining. In this study, a computed tomography (CT) scanner and a triaxial loading seepage system were used to perform CT scanning tests on gas-bearing coal under triaxial compressive conditions. We processed the CT images in different deformation stages of the gas-bearing coal and analyzed them with a digital image processing technology to depict dynamic evolution damage trend. In addition, we employed a fracture area projection method to propose a new damage variable and theoretically constructed a new damage constitutive model. The results showed that the coal specimens became prone to increased instability and damage as gas pressure increased. The growth and expansion of fissures in gas-bearing coal were heterogeneous and anisotropic and mainly experienced three processes, that is, closure in the compaction stage and elasticity stage, expansion in the strain-hardening stage, and sharp propagation and interconnection in the post-peak stage. With the progression of fracture propagation, the degree of damage gradually decreased then increased. The proposed new damage variable was more suitable than the traditional damage variable for depicting the damage degree and damage evolution process. Furthermore, the proposed new damage constitutive model well mirrored the deformation behavior of gas-bearing coal under triaxial compressive loads.

Nonlinear Seepage Behaviors of Pore-Fracture Sandstone under Hydro-Mechanical Coupling

This work focused on the nonlinear seepage behaviors of flow in pore-fracture media. Natural sandstones were selected to prefabricate single-fracture specimens with different inclinations (0–90°). Seepage tests of combined media were performed under different confining pressures (8–10 MPa) and different water pressures (3–7 MPa) in a triaxial pressure chamber. The fitting analysis of experimental data showed that Forchheimer’s law described the nonlinear characteristics of flow in the pore-fracture media. Linear term coefficient a and nonlinear term coefficient b of the sandstone samples with different inclinations changed more obviously with the increased inclination. When the fracture inclination was greater than 30°, a and b values had a sudden jump. The nonlinear inertial-parameter equation of fluid flow in pore-fracture media was proposed based on non-Darcy flow coefficient β and inherent permeability k. The applicability of the following methods to evaluate Darcy’s law was discussed, including normalized hydraulic conductivity, pressure gradient ratio, and discharge ratio. The three methods were able to determine critical parameters and distinguish linear and nonlinear flow. Furthermore, it was specified for the first time that when β was negative, critical nonlinear effect E was −0.1, and Forchheimer’s coefficient F0 was −0.091. In the −∇P-Q relationship, the fitting curve was convex to the −∇P axis, and the increase of Q was higher than the linear increase, presenting the nonlinearity of overflow. On the one hand, the fractures and pores were compressed under the confining pressure due to the prefabricated fractures of different shapes and different inclinations. A higher seepage water pressure was needed to stabilize the seepage system with the excessive flow rate. On the other hand, the barrier effect of the fluid inside the rock was completely lost because the fluid expanded the seepage channel. Its permeability was changed, leading to seepage instability.

Origin and Periodic Behavior of Short Duration Signals Recorded by Seismometers at Vestnesa Ridge, an Active Seepage Site on the West-Svalbard Continental Margin

Short duration events (SDEs) are reported worldwide from ocean-bottom seismometers (OBSs). Due to their high frequency (4–30 Hz) and short duration, they are commonly attributed to aseismic sources, such as fluid migration related processes from cold seeps, biological signals, or noise. We present the results of a passive seismic experiment that deployed an OBS network for 10-month (October 2015–July 2016) at an active seepage site on Vestnesa Ridge, West Svalbard continental margin. We characterize SDEs and their temporal occurrence using the conventional short-time-average over long-time-average approach. Signal periodograms show that SDEs have periodic patterns related to solar and lunar cycles. A monthly correlation between SDE occurrences and modelled tides for the area indicates that tides have a partial control on SDEs recorded over 10 months. The numbers of SDEs increase close to the tidal minima and maxima, although a correlation with tidal highs appears more robust. Large bursts of SDEs are separated by interim quiet cycles. In contrast, the periodicity analysis of tremors shows a different pattern, likely caused by the effect of tidally controlled underwater currents on the instrumentation. We suggest that SDEs at Vestnesa Ridge may be related to the dynamics of the methane seepage system which is characterized by a complex interaction between migration of deep sourced fluids, gas hydrate formation and seafloor gas advection through cracks. Our observation from this investigated area offshore west-Svalbard, is in line with the documentation of SDEs from other continental margins, where micro-seismicity and gas release into the water column are seemingly connected.

Experimental Study of the Influence of Moisture Content on the Pore Structure and Permeability of Anthracite Treated by Liquid Nitrogen Freeze-Thaw.

The liquid nitrogen freeze–thaw (LN2-FT) method has been widely used to improve the coal permeability in the coalbed methane (CBM) production. However, the influence of moisture content on the permeability of coal treated by LN2-FT remains unclear, limiting the broad application of this technique. A novel seepage system was proposed to analyze the permeability evolution of anthracite coal samples treated by LN2-FT. Moreover, variations of the pore structure were analyzed using scanning electron microscopy (SEM), mercury intrusion porosimetry (MIP), and low-field NMR. The results showed that pores and fractures appeared on the coal surface after the LN2-FT treatment. As the moisture content of the coal increased, more pores and fractures tended to be formed during the LN2-FT treatment. The total pore volume, porosity, and average pore diameter of the anthracite coal after the treatment were 1.77, 2.44, and 5.58 times higher, respectively, than that of the raw coal. The change in the specific surface area exhibited three trends as the moisture content of the coal samples increased: a slow descent, a steady increase, and a rapid descent. Moreover, it was found that the LN2-FT treatment increased the connections between pores and fractures, improving gas migration in the coal. Furthermore, the LN2-FT treatment significantly increased the permeability of the anthracite coal samples. The higher the coal moisture, the higher the permeability of the coal samples after the LN2-FT treatment. Hence, the LN2-FT technique can substantially improve the permeability of coal reservoirs, providing essential information for the efficient utilization of CBM.

Chemical and isotopic features of seepage gas from mud volcanoes in southern margin of the Junggar Basin, NW China

The chemical and stable carbon isotopic composition of seepage gases from three mud volcanoes were analyzed and compared with their corresponding reservoir gases in the southern margin of the Junggar Basin. The results showed a clear alteration of the seepage gases with respect to their chemical and isotopic composition. The seepage gases released from vents/craters were relatively enriched in methane and had a decrease in ethane and propane contents when compared to their corresponding reservoir gases. Advective process during gas ascent to the surface was the dominant post-genetic process that induced such molecular fractionation. The diffusive process during gas migration through sandy rocks or liquid horizons and gas-water-rock-microbe interactions also contributed to the molecular fractionation. High gas flux and eruptive behaviors triggered by earthquakes can however restrain this kind of fractionation. The seepage gases had relatively lighter carbon isotopes in CH4, similar carbon isotopes in C2H6 and heavier carbon isotopes in C3H8 when compared to their corresponding reservoir gases. Diffusive process led to the lighter carbon isotopes of the seepage gases due to mass-dependent kinetic isotope effect, whereas gas-water-rock-microbe interactions contributed to the heavier ones. The high variability of δ13CCO2 values might be caused by the differently biodegraded seepage systems or by the high sensitivity of δ13CCO2 values to gas-water-rock interactions. Such molecular and isotopic fractionations caused by geochemical processes provide distinctive indicators for a better understanding of gas migration from reservoirs to the ground surface. It is necessary to take these geochemical processes into consideration when we assess gas origin and source rocks based on seepage gases.