Kinetic-simulating experiment combined with GC-IRMS analysis: application to identification and assessment of coal-derived methane from Zhongba Gas Field (Sichuan Basin, China)

The main purpose of this study is to model the δ¹³C values of methane derived from coal by combining kinetic-simulating experiment with the gas chromatography-isotope ratio mass spectrum (GC-IRMS) analysis. The stable carbon isotopic variation of methane in pyrolysates with heating temperature indicates that the assumptions for both a constant kinetic isotope effect (α) and a uniform initial isotopic composition (δ¹³C_o) are impractical for explaining the carbon isotopic fractionation during coalification. For purposes of simplification, two approaches are used in this paper to deal with the heterogeneity of terrestrial organic matter. One is that, assuming a uniform initial isotopic composition (i.e., δ¹³C_{i, o}=δ¹³C_o) for all methane-generating precursors in coal, the isotopic variation of methane is fitted by adjusting ΔE_{a, i} (E_{a¹³C, i}−E_{a¹²C, i}) for each hypothetical reaction. The other is that, assuming a constant kinetic isotope effect during the whole gas formation, that is all ΔE_{a, i} values are identical, the modeling of methane isotopic composition is achieved by changing the ¹³CH₄ generation potential of each reaction (f_{i, ¹³C}), namely, by adjusting the initial δ¹³C value (δ¹³C_{i, o}) for each methane-generating precursor. Results of the kinetic calculation shows that the two simulating methods can yield a similar result at a geological heating rate of 2 °C/My, which further demonstrates that those natural gases with methane δ¹³C value being approximately −36‰ are possibly sourced from the upper Triassic coal measure strata in the Northwestern Sichuan Basin.     

On the long-term seismic hazard analysis in the Zhangjiakou–Penglai seismotectonic zone, China

The Zhangjiakou–Penglai seismotectonic zone (ZPSZ) lies in the northern part of North China and extends along the Zhangjiakou–Beijing–Tianjin–Bohai Bay–Penglai–Yellow Sea. It is about 900 km long and some 250 km wide in a northwest direction. The great Sanhe-Pinggu (M_S=8.0) earthquake occurred on September 1679 and the Tangshan (M_S=7.8) earthquake on July 1976 caused serious economic and life losses. According to some differences in crust structure and regional tectonic stress field, the ZPSZ is divided into western and eastern segment by the 117°E line for study on long-term seismic hazard analysis. An analysis of Gutenberg–Richter's empirical relation of earthquake-frequency and time process of historic and recent earthquakes along the eastern and western segments shows that the earthquake activity obeys a Poisson process, and these calculations indicate that the earthquake occurrence probability of M_S=6.0–6.9 is 0.77–0.83 in the eastern segment and the earthquake occurrence probability of M_S=7.0–7.9 is 0.78–0.80 in the western segment of the ZPSZ during a period from 2005 to 2015.     

A Physical and Numerical Investigation on the Stability of Shallow Tunnels in Strain Softening Media

Summary The behaviour up to failure of shallow underground openings is discussed on the basis of some laboratory, small-scale model tests and of finite element simulation. The experimental results are first illustrated. They were obtained from two-dimensional (plane strain) and three-dimensional tunnel models tested under standard gravity conditions. Then, the phenomenon of strain localisation that characterizes the medium surrounding the model tunnels is discussed, recalling two alternative approaches for its numerical interpretation. On this basis, a finite element procedure for strain softening analyses is outlined and applied to the simulation of the tests in both two- and three-dimensional conditions. The comparison between experimental and numerical results leads to some conclusions on the influence of strain localisation on the overall behaviour of shallow tunnels and on the stability of their headings.     

A time-space based approach for mapping rainfall-induced shallow landslide hazard

A rainfall-induced shallow landslide is a major hazard in mountainous terrain, but a time-space based approach is still an unsettled issue for mapping rainfall-induced shallow landslide hazards. Rain induces a rise of the groundwater level and an increase in pore water pressure that results in slope failures. In this study, an integrated infinite slope analysis model has been developed to evaluate the influence of infiltration on surficial stability of slopes by the limit equilibrium method. Based on this new integrated infinite slope analysis model, a time-space based approach has been implemented to map the distributed landslide hazard in a GIS (Geographic Information Systems) and to evaluate the shallow slope failure induced by a particular rainfall event that accounts for the rainfall intensity and duration. The case study results in a comprehensive time-space landslide hazard map that illustrates the change of the safety factor and the depth of the wetting front over time.     

Petrology and Geochemistry of Early Cretaceous Bimodal Continental Flood Volcanism of the NW Etendeka, Namibia. Part 1: Introduction, Mafic Lavas and Re-evaluation of Mantle Source Components

The bimodal NW Etendeka province is located at the continentalend of the Tristan plume trace in coastal Namibia. It comprisesa high-Ti (Khumib type) and three low-Ti basalt (Tafelberg,Kuidas and Esmeralda types) suites, with, at stratigraphicallyhigher level, interstratified high-Ti latites (three units)and quartz latites (five units), and one low-Ti quartz latite.Khumib basalts are enriched in high field strength elementsand light rare earth elements relative to low-Ti types and exhibittrace element affinities with Tristan da Cunha lavas. The unradiogenic²⁰⁶Pb/²⁰⁴Pb ratios of Khumib basalts are distinctive, most plottingto the left of the 132 Ma Geochron, together with elevated ²⁰⁷Pb/²⁰⁴Pbratios, and Sr–Nd isotopic compositions plotting in thelower ¹⁴³Nd/¹⁴⁴Nd part of mantle array (EM1-like). The low-Tibasalts have less coherent trace element patterns and variable,radiogenic initial Sr (     

Inferring the mass fraction of floc-deposited mud: application to fine-grained turbidites

Fine sediment deposition in the ocean is complicated by the cohesive nature of muds and their tendency to flocculate. The result is disaggregated inorganic grain size (DIGS) distributions of bottom sediment that are influenced by single‐grain and floc deposition. This study outlines a parametric model that characterizes bottom sediment DIGS distributions. Modelled parameters are then used to infer depositional conditions that account for the regional variation in the grain sizes deposited by turbidity currents on the Laurentian Fan–Sohm Abyssal Plain, offshore south‐eastern Canada. Results indicate that, on the channellized Laurentian Fan, the mass fraction of floc‐deposited mud increases only slightly downslope. The small evolution in this fraction arises because sediment concentration and turbulent energy are associated in turbidity currents. On the Sohm Abyssal Plain, however, the mass fraction of floc‐deposited mud decreases, probably as a result of lower sediment concentration at this source‐distal site. Estimates of the mass fraction of mud deposited as flocs suggest that floc deposition is the dominant mode by which sediment is lost from suspension, although single‐grain deposition contributes more to the depositional flux in proximal areas where high energy breaks flocs and in distal areas where low sediment concentration limits floc formation. It is concluded that, throughout the dispersal system, changes in the fraction of flocculated mud deposited from turbidity currents reflect changes in sediment concentration and energy downslope.     

A Finite Volume Scheme for the Transport of Radionucleides in Porous Media

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.     

Application of AVO Analysis to Seismic Data for Detection of Gas below Methane Hydrate Stability Zone in Nankai Trough Area

Abstract. For the purpose of development of methane hydrate, occurring in the deep marine subsurface, as a resource, the most important issue is to understand the methane hydrate system (generation, migration and accumulation) as well as to delineate the methane hydrate reservoir properties. We have applied the Amplitude Versus Offset (AVO) analysis to the seismic data acquired in the Nankai Trough, offshore Japan, in order to confirm the occurrence of gas just below the methane hydrate-bearing zone, assuming that gas will show a so-called Class-3 AVO response. Knowledge of the amount and occurrence of gas in the sediment below methane hydrate-bearing zone is one of the keys to understand the methane hydrate system.
We have utilized the qualitative analysis of AVO methodology to delineate how gas is located below the BSR, which is thought to be the reflection event from the interface between the methane hydrate-bearing zone and the underlying gas-bearing zone. In the region of MITI Nankai Trough Well PSW-3, we observe two BSRs separated by 25 ms. After AVO modeling using well data, we applied AVO attribute analysis and attribute crossplot analysis to the seismic data. Finally we applied an offset-amplitude analysis to CMP gather data at specific locations to confirm the results of AVO attribute analysis. The AVO analysis shows that there is very little gas located in the underlying sediment below methane hydrate-bearing zone. This result supports the fact that we could not obtain any clear evidence of gas occurrence just below the methane hydrate-bearing zone in the Nankai Trough well drilling.     

Regional moment tensor uncertainty due to mismodeling of the crust

The retrieval of earthquake moment tensor (MT) requires the response of the medium, in which seismic waves travel from the hypocenter to the stations, to be known. In inverting long-period (LP) seismic data (teleseismic and LP regional records), a gross earth model is sufficient; with decreasing periods, a more detailed model is needed. This is the case when waveforms of weak earthquakes at regional distances are to be inverted. Regional moment tensors (RMTs) of mostly Mediterranean earthquakes are determined on a routine basis by the Swiss Seismological Survey (SED) by using averaged models of the earth's crust. By inverting broad-band records of the M_w=4.8 earthquake near Udine, N Italy, on Feb. 14, 2002, we tested the sensitivity of the MT solution with respect to possible errors in the earth model used and in the location of the hypocenter depth. We perturbed the P and S velocities and the thickness in the 1-D earth model in the range from 3% to 30% of the parameter values and constructed estimates of confidence regions of the MT and error bars of the source time function (STF) and scalar moment in three frequency bands. Similarly, these error characteristics were determined assuming a mislocation in the hypocenter depth. We found that, in the band of periods from 25 to 50 s, the mechanism is resolved well (at the confidence level 95% at least) up to an earth model uncertainty of 30%, in the passband 10–25 s up to about 10%, but it is undetermined completely at periods of 5–10 s. An error in hypocenter depth of as much as double the value reported by the location procedure does not destroy the resolution of the mechanism at periods above 10 s. In the RMT catalog of the SED, earthquakes of M_w greater than about 3.5 are processed at periods above 30 s; thus, the solutions for these events are robust with respect to a possible uncertainty in the earth model used. Mechanisms of weaker earthquakes, retrieved from short periods, should be interpreted with caution.     

Curve-Fitting Schemes for Fragmentation Data

Fragmentation measurements in the form of sieve passing and mass fraction data were used to test the capability of three different distributions to fit the observed data over a wide range in fragment size and mass. These distributions were based on Rosin-Rammler, lognormal and simple sigmoidal (S-shaped) functions, having 2 input parameters for the single-component versions and 5 input parameters for the two-component versions. Provided convergence was achieved in the non-linear curve-fitting technique, the two-component versions always provided superior fits to the observed data. However, these versions were very sensitive to variations in the values chosen for the input parameters. In this particular regard, the two-component sigmoidal function was the most robust. The present results also show that the two-component lognormal function provided the best fit to the fragmentation data in a general sense, and the two-component Rosin-Rammler function provided the worst fit. However, there was not a significant difference between any of the three methods.