Influences of loading direction and intermediate principal stress ratio on the initiation of strain localization in cross-anisotropic sand

Since cross-anisotropic sand behaves differently when the loading direction or the stress state changes, the influences of the loading direction and the intermediate principal stress ratio (b = (σ ₂ ? σ ₃)/(σ ₁ ? σ ₃)) on the initiation of strain localization need study. According to the loading angle (angle between the major principal stress direction and the normal of bedding plane), a 3D non-coaxial non-associated elasto-plasticity hardening model was proposed by modifying Lode angle formulation of the Mohr–Coulomb yield function and the stress–dilatancy function. By using bifurcation analysis, the model was used to predict the initiation of strain localization under plane strain and true triaxial conditions. The predictions of the plane strain tests show that the major principal strain at the bifurcation points increases with the loading angle, while the stress ratio decreases with the loading angle. According to the loading angle and the intermediate principal stress ratio, the true triaxial tests were analyzed in three sectors. The stress–strain behavior and the volumetric strain in each sector can be well captured by the proposed model. Strain localization occurs in most b value conditions in all three sectors except for those which are close to triaxial compression condition (b = 0). The difference between the peak shear strength corresponding to the strain localization and the ultimate shear strength corresponding to plastic limit becomes obvious when the b value is near 0.4. The influence of bifurcation on the shear strength becomes weak when the loading direction changes from perpendicular to the bedding plane to parallel. The bifurcation analysis based on the proposed model gives out major principal strain and peak shear strength at the initiation of strain localization; the given results are consistent with experiments.     

DEM simulations of sandstone under true triaxial compressive tests

Numerically simulated true triaxial compression tests (σ ₁ ≥ σ ₂ ≥ σ ₃) are conducted in this study to elucidate the failure mechanism of sandstone using 3D discrete element method (DEM), in particular the effect of the intermediate principal stress (σ ₂). Eight series of tests (σ ₃ = 0, 10, 20, 30, 40, 50, 70, and 100 MPa) are conducted. Within each series, σ ₂ is varied from σ ₂ = σ ₃ to σ ₂ = σ ₁ from test to test. For each test, σ ₁ is raised monotonically to failure while keeping σ ₂ and σ ₃ constant. The DEM simulations reveal the effect of σ ₂ on the variations of peak stress, Young’s modulus, failure plane angles, the brittle–ductile transition, and the evolution of failure modes, the effect beyond the well-understood effect of σ ₃. The simulation is in qualitative agreement with the results obtained experimentally. Detailed analyses performed on the particle-scale responses further the understanding of the microscopic mechanisms. The distribution of contact force becomes more homogeneous with the increase of σ ₃, which leads to the resulting damage being more localized rather than diffused. The interaction between contact force distribution and coalescence of cracks determines the processes and patterns of fracturing in the sample scale. σ ₂ is found to affect the microscopic stress distribution as well as structure evolution, and this effect weakens with the increase of σ ₃.     

Detailed comparison of nine intact rock failure criteria using polyaxial intact coal strength data obtained through PFC<Superscript>3D</Superscript> simulations

Study of intact rock failure criteria is an important topic in rock mechanics. In this study, applicability of nine different intact rock failure criteria is investigated for intact coal strength data. PFC^3D modeling was used to simulate the laboratory polyaxial tests for cubic intact coal blocks of side dimension 110 mm under different confining stress combinations. A modified grid search procedure is proposed and used to find the best-fitting parameter values and to calculate the coefficient of determination (R ²) values for each criterion. Detailed comparisons of the nine criteria are made using the following aspects: R ² values, σ ₁ ? σ ₂ plots for different σ ₃, shapes on the deviatoric plane, linearity or nonlinearity on the meridian planes. Through the comparisons of R ² values, σ ₁ ? σ ₂ plots and meridian lines, the modified Wiebols–Cook and modified Lade criteria were found to fit the intact coal strength data best. The nine failure criteria are categorized into three types based on the appearances on the deviatoric plane.     

Permeability of granular soil employing flexible wall permeameter

Understanding the changes in permeability of soil, when soil is subjected to high confining pressure and flow pressure, which may alter the textural and geomechanical characteristics of soil, is of great importance to many geo-engineering activities such as, construction of high-rise buildings near the coast or the water bodies, earthen dams, pavement subgrades, reservoir, and shallow repositories. It is now possible to evaluate the changes in permeability of soil samples under varying conditions of confining pressure and flow pressure using flexible wall permeameter (FWP). In the present study, investigation was carried out on a cylindrical sample of granular soil employing FWP under varied conditions of confining pressure (σ₃)—50–300 kPa, which can simulate the stress conditions equivalent to depth of about 20 m under the earth’s crust, and a flow pressure (f_p)—20–120 kPa, which is mainly present near the small earthen embankment dams, landfill liners, and slurry walls near the soft granular soil with high groundwater table. The obtained results indicate a linear relationship between hydraulic conductivity (k) with effective confining pressure (σ_eff.), k, decreasing linearly with an incremental change in σ_eff.. Further, k increases significantly with an increase in f_p corresponding to each σ_eff., and q increases significantly with increase in the f_p corresponding to each (σ₃). It was also observed that corresponding to the low f_p of 20 kPa, the reduction in k is nonlinear with σ_3. The percentage reduction in k is observed to be 9, 13, and 27% corresponding to σ₃ of 50–100, 100–200, and 200-300 kPa, respectively.     

Distribution rule of the in situ stress state and its influence on the permeability of a coal reservoir in the southern Qinshui Basin,China

The recent development of the coalbed methane (CBM) industry has a significant role in advancing hydraulic fracturing theory and technology. However, further development requires a better understanding of how fractures influence reservoir permeability. In situ stress data from 54 CBM wells in the southern Qinshui Basin, China, were obtained by the injection/falloff test method to analyse the effect of in situ stress on the permeability of the CBM reservoir. The types of in situ stress states were classified, and the coal reservoir permeability under different in situ stress states was analysed. The results indicate that the maximum horizontal principal stress (σ_H), minimum horizontal principal stress (σ_h) and vertical principal stress (σ_v) all have positive linear relationships with the coal seam burial depth. Three in situ stress states were observed from the shallow to deep regions of the CBM reservoir in the study area: σ_H?>?σ_h?>?σ_v, σ_H?>?σ_v?>?σ_h and σ_v?>?σ_H?>?σ_h, which account for 9, 76 and 15% of the test wells, respectively. Coal reservoir permeability decreases with increasing horizontal principal stress, whereas it first decreases with increasing σ_v, then increases and finally decreases. The variation in permeability with σ_v is due to the conversion of the in situ stress states. Coal reservoir permeability has obvious differences under different in situ stress states. The permeability is the largest when σ_v?>?σ_H?>?σ_h, followed by σ_H?>?σ_h?>?σ_v and smallest when σ_H?>?σ_v?>?σ_h. The permeability differences are caused by the fracture propagation shape of the rock strata under different in situ stress states.     

Dewatering Behaviour of a Uranium Ore Slurry Containing Clays

The main objective of this paper was to investigate the dewatering behaviour of a clayey uranium ore slurry. The slurry (containing 28% clay size) exhibited moderate water adsorption (w _l  = 83% and w _p  = 30%). Primarily composed of muscovite (46%) and quartz (30%), the clay minerals included illite (8%), chlorite (5%) and kaolinite (2%) alongside a CEC of 41 (cmol(+)/kg) with Ca²⁺ and Mg²⁺ as the dominant cations. Likewise, the high EC (17,600 μS/cm) and ionic strength (1.15 mol/L) indicated a flocculated microstructure due to the presence of SO₄ ^2? (22,600 mg/L) and Mg²⁺ (1340 mg/L) in the slurry water. Settling included sedimentation and consolidation at low initial solids condition (25–35%) whereas only consolidation was observed at high initial solids contents (40–50%). The average k reduced from 1.2 × 10^?6 m/s (initial s = 25%) to 5.3 × 10^?8 m/s (initial s = 50%) along with a void ratio reduction from 7.4 to 2.6. Due to thixotropic strength, volume compressibility during consolidation showed apparent pre-consolidation at low effective stress (0.3–2 kPa) with a reduction in void ratio from 2.6 to 2.5. The e _s was found to be 2.46 at σ′ = 2 kPa and was followed by a steeper slope with the void ratio reducing to 2.1 at σ′ = 31 kPa. Likewise, the hydraulic conductivity during consolidation decreased from 2.6 × 10^?9 m/s (at e = 2.6) to 2.0 × 10^?10 m/s (at e = 2.1).     

Impact of pore water conductivity and porosity on the electrical conductivity of kaolinite

The purpose of this study is to quantify the magnitudes of surface conduction and pore water conduction from the measured electrical conductivity of kaolinite, with the ultimate goal of estimating the electrical conductivity of kaolinite with a wide range of pore water conductivities (σ _w = 0.013–3.356 S/m) and porosities (n = 0.368–1.0). Therefore, the theoretical background of the electrical conductivity in soils was reviewed, and electrical conductivity measurements on kaolinite were performed using both slurry and consolidation tests in this study. The results of this study demonstrate that the variations of measured electrical conductivity (σ _mix) with n are debatable according to the values of σ _w, because a decrease in n results in both an increase in surface conduction (K _s) and a decrease in pore water conduction (K _w); this causes the relative magnitude of K _s compared to that of K _w to vary with σ _w and n. Consequently, this study develops the relation between the porosity-normalized K _s/K _w and 1/σ _w. Additionally, the surface conductivity of the tested kaolinite is back-calculated and compared with the previous relationship between K _s and zeta potential of kaolinite. The measured and estimated σ _mix values are compared with the varying pore water conductivity and porosity values.     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: <Emphasis Type="Italic">b</Emphasis>-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.     

Sequence signal processing of 2012 Ahar–Varzaghan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> 6.4) of NW Iran

Before starting seismic cycle of Ahar–Varzaghan 2012 event, a partial gap in the form of a pre-seismic calm sequence (seismicity rate, r = 0.46 event/year, b = 1.4) with duration of 303 days spatially has dominated over the entire seismogenic area. From April 17, 2012, to May 31, 2012, r significantly increased to 2.16, indicating strong foreshock sequence, and b value changed to 1.9, remarkably. In the last two months before the mainshock, foreshocks have partially migrated toward the earthquake fault (with a decrease in size, b = 2.0). Significantly, high rate of seismicity and low V _P /V _S (1.64) in the foreshocks sequence and also very high seismicity rate (17.3) and high V _P /V _S (1.76) in the aftershocks sequence make substantial differences between the seismic cycle and the background seismicity. Moreover, a significant E–W migration of the microseismicity was confirmed in the study area.     

High-pressure displacive phase transition of a natural Mg-rich pigeonite

A high-pressure single-crystal X-ray diffraction study has been carried out on a P2₁/c natural Mg-rich pigeonite sample with composition ca. Wo₆En₇₆Fs₁₈ using a diamond anvil-cell. The unit-cell parameters were determined at 14 different pressures to 7.14 GPa. The sudden disappearance of the b-type reflections (h + k = odd) and a strong discontinuity (about 2.8%) in the unit-cell volume indicated a first-order P2₁/c–C2/c phase transition between 4.66 and 4.88 GPa. The P(V) data of the P2₁/c phase were fitted to 4.66 GPa by a third-order Birch–Murnaghan equation of state (BM3 EoS), whereas the limited number of experimental data collected within the C2/c phase between 4.88 and 7.14 GPa were fitted using the same equation of state but with K′ constrained to the value obtained for the P2₁/c fitting. The equation of state coefficients are V ₀ = 424.66(6) ų, K _T0 = 104(2) GPa and K′ = 8(1) for the P2₁/c phase, and V ₀ = 423.6(1) ų, K _T0 = 112.4(8) GPa, and K′ fixed to 8(1) for the C2/c phase. The axial moduli for a, b, and c for the P2₁/c phase were obtained using also a BM3-EoS, while for the C2/c phase only a linear calculation could be performed, and therefore the same approach was applied for comparison also to the P2₁/c phase. In general the C2/c phase exhibits axial compressibilities (β _c > β _a >> β _b) lower than those of the P2₁/c phase (β _b > β _c ≈ β _a; similar to those found in previous studies in clinopyroxenes and orthopyroxenes). The lower compressibility of the C2/c phase compared with that of the P2₁/c could be ascribed to the greater stiffness along the b direction. A previously published relationship between P _c and M2 average cation radius (i.r.) has been updated using all the literature data on P2₁/c clinopyroxene containing large cations at M2 site and our new data. The following weighted regression was obtained: P _c (GPa) = 26(4) ? 28(5) ×  i.r (Å), R ² = 0.97. This improved equation can be used to predict the critical pressure of natural P2₁/c clinopyroxene samples just knowing the composition at M2 site.     

Stability at high-pressure,elastic behaviour and pressure-induced structural evolution of CsAlSi<Subscript>5</Subscript>O<Subscript>12</Subscript>, a potential host for nuclear waste

The elastic and structural behaviour of the synthetic zeolite CsAlSi₅O₁₂ (a = 16.753(4), b = 13.797(3) and c = 5.0235(17) Å, space group Ama2, Z = 2) were investigated up to 8.5 GPa by in situ single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions. No phase-transition occurs within the P-range investigated. Fitting the volume data with a third-order Birch–Murnaghan equation-of-state gives: V ₀ = 1,155(4) ų, K _T0 = 20(1) GPa and K′ = 6.5(7). The “axial moduli” were calculated with a third-order “linearized” BM-EoS, substituting the cube of the individual lattice parameter (a ³, b ³, c ³) for the volume. The refined axial-EoS parameters are: a ₀ = 16.701(44) Å, K _T0a = 14(2) GPa (β_a = 0.024(3) GPa^?1), K′ _a = 6.2(8) for the a-axis; b ₀ = 13.778(20) Å, K _T0b = 21(3) GPa (β_b = 0.016(2) GPa^?1), K′ _b = 10(2) for the b-axis; c ₀ = 5.018(7) Å, K _T0c = 33(3) GPa (β_c = 0.010(1) GPa^?1), K′ _c = 3.2(8) for the c-axis (K _T0a:K _T0b:K _T0c = 1:1.50:2.36). The HP-crystal structure evolution was studied on the basis of several structural refinements at different pressures: 0.0001 GPa (with crystal in DAC without any pressure medium), 1.58(3), 1.75(4), 1.94(6), 3.25(4), 4.69(5), 7.36(6), 8.45(5) and 0.0001 GPa (after decompression). The main deformation mechanisms at high-pressure are basically driven by tetrahedral tilting, the tetrahedra behaving as rigid-units. A change in the compressional mechanisms was observed at P ≤ 2 GPa. The P-induced structural rearrangement up to 8.5 GPa is completely reversible. The high thermo-elastic stability of CsAlSi₅O_12, the immobility of Cs at HT/HP-conditions, the preservation of crystallinity at least up to 8.5 GPa and 1,000°C in elastic regime and the extremely low leaching rate of Cs from CsAlSi₅O₁₂ allow to consider this open-framework silicate as functional material potentially usable for fixation and deposition of Cs radioisotopes.     

Accuracy assessment of near-shore bathymetry information retrieved from Landsat-8 imagery

The improvement in the capabilities of Landsat-8 imagery to retrieve bathymetric information in shallow coastal waters was examined. Landsat-8 images have an additional band named coastal/aerosol, Band 1: 435–451 nm in comparison with former generation of Landsat imagery. The selected Landsat-8 operational land image (OLI) was of Chabahar Bay, located in the southern part of Iran (acquired on February 22, 2014 in calm weather and relatively low turbidity). Accurate and high resolution bathymetric data from the study area, produced by field surveys using a single beam echo-sounder, were selected for calibrating the models and validating the results. Three methods, including traditional linear and ratio transform techniques, as well as a novel proposed integrated method, were used to determine depth values. All possible combinations of the three bands [coastal/aerosol (CB), blue (B), and green (G)] have been considered (11 options) using the traditional linear and ratio transform techniques, together with five model options for the integrated method. The accuracy of each model was assessed by comparing the determined bathymetric information with field measured values. The standard error of the estimates, correlation coefficients (R ² ) for both calibration and validation points, and root mean square errors (RMSE) were calculated for all cases. When compared with the ratio transform method, the method employing linear transformation with a combination of CB, B, and G bands yielded more accurate results (standard error = 1.712 m, R ² _calibration = 0.594, R ² _validation = 0.551, and RMSE =1.80 m). Adding the CB band to the ratio transform methodology also dramatically increased the accuracy of the estimated depths, whereas this increment was not statistically significant when using the linear transform methodology. The integrated transform method in form of Depth = b ₀  + b ₁ X _CB  + b ₂ X _B  + b ₅ ln(R _CB )/ln(R _G ) + b ₆ ln(R _B )/ln(R _G ) yielded the highest accuracy (standard error = 1.634 m, R ² _calibration = 0.634, R ² _validation = 0.595, and RMSE = 1.71 m), where R _i (i = CB, B, or G) refers to atmospherically corrected reflectance values in the i ^th band [X _i  = ln(R _i -R _{deep water})].     

Fracture evolution and energy mechanism of deep-buried carbonaceous slate

In order to study the influence of confining pressure and water content on the mechanical properties, fracture evolution and energy damage mechanism of deep-buried carbonaceous slate, uniaxial and triaxial compression tests were carried out under natural and saturated states and acoustic emission monitored. The deep-buried carbonaceous slate samples were obtained at a depth of 1020 m from the Lanjiayan tunnel in Sichuan province, China, where the maximum in situ stress has been measured at 44.2 MPa. The results suggest that water has a significant softening effect on the strength and deformation characteristics of carbonaceous slate, but the effect decreases with an increase in the confining pressure. When both the confining pressure and water content are increased, the acoustic emission events and dissipated energy gradually increase at the pre-peak and post-peak stages. Thus, the AE evolution type seen in the natural state under low confining pressure usually presents as a main shock-type event, and it changes to a foreshock–main shock–after shock event when saturated and at high confining pressures. Based on the S-shaped energy evolution law, the damage evolution process of carbonaceous slate was analyzed. The damage stress thresholds σ _ea and σ _eb were obtained, which can be considered as the thresholds of the rock entering the energy-hardening and energy-softening stages. Finally, a new brittleness energy index BDE is proposed to describe the influence of confining pressure and water content on the damage mechanism of deep-buried carbonaceous slate.     

High-pressure phase transitions of CaRhO<Subscript>3</Subscript> perovskite

High-pressure phase transitions of CaRhO₃ perovskite were examined at pressures of 6–27 GPa and temperatures of 1,000–1,930°C, using a multi-anvil apparatus. The results indicate that CaRhO₃ perovskite successively transforms to two new high-pressure phases with increasing pressure. Rietveld analysis of powder X-ray diffraction data indicated that, in the two new phases, the phase stable at higher pressure possesses the CaIrO₃-type post-perovskite structure (space group Cmcm) with lattice parameters: a = 3.1013(1) Å, b = 9.8555(2) Å, c = 7.2643(1) Å, V _m  = 33.43(1) cm³/mol. The Rietveld analysis also indicated that CaRhO₃ perovskite has the GdFeO₃-type structure (space group Pnma) with lattice parameters: a = 5.5631(1) Å, b = 7.6308(1) Å, c = 5.3267(1) Å, V _m  = 34.04(1) cm³/mol. The third phase stable in the intermediate P, T conditions between perovskite and post-perovskite has monoclinic symmetry with the cell parameters: a = 12.490(3) Å, b = 3.1233(3) Å, c = 8.8630(7) Å, β = 103.96(1)°, V _m  = 33.66(1) cm³/mol (Z = 6). Molar volume changes from perovskite to the intermediate phase and from the intermediate phase to post-perovskite are –1.1 and –0.7%, respectively. The equilibrium phase relations determined indicate that the boundary slopes are large positive values: 29 ± 2 MPa/K for the perovskite—intermediate phase transition and 62 ± 6 MPa/K for the intermediate phase—post-perovskite transition. The structural features of the CaRhO₃ intermediate phase suggest that the phase has edge-sharing RhO₆ octahedra and may have an intermediate structure between perovskite and post-perovskite.     

Properties of the aftershocks sequences of the 2000 and 2002 earthquakes in Aksehir-Afyon graben,west-central Anatolia,Turkey

The b-value of the Gutenberg–Richter’s frequency–magnitude relation and the p-value of the modified Omori law, which describes the decay rate of aftershock activity, were investigated for more than 500 aftershocks in the Aksehir-Afyon graben (AAG) following the 15 December 2000 Sultandagi–Aksehir and the 3 February 2002 Çay–Eber and Çobanlar earthquakes. We used the Kandilli Observatory’s catalog, which contains records of aftershocks with magnitudes ≥2.5. For the Çobanlar earthquake, the estimated b-values for three aftershock sequences are in the range 0.34 ≤  b ≤ 2.85, with the exception of the one that occurred during the first hour (4.77), while the obtained p-values are in the range 0.44 ≤ p ≤ 1.77. The aftershocks of the Sultandagi earthquake have a high p-value, indicating fast decay of the aftershock activity. A regular increase of b can be observed, with b < 1.0 after 0.208 days for the Çay–Eber earthquake. A systematic and similar increase and decrease pattern exists for the b- and p-values of the Çobanlar earthquakes during the first 5 days.     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

The stress regime in a Rotliegend reservoir of the Northeast German Basin

In-situ stresses have significant impact, either positive or negative, on the short and long term behaviour of fractured reservoirs. The knowledge of the stress conditions are therefore important for planning and utilization of man-made geothermal reservoirs. The geothermal field Groß Schönebeck (40 km north of Berlin/Germany) belongs to the key sites in the northeastern German Basin. We present a stress state determination for this Lower Permian (Rotliegend) reservoir by an integrated approach of 3D structural modelling, 3D fault mapping, stress ratio definition based on frictional constraints, and slip-tendency analysis. The results indicate stress ratios of the minimum horizontal stress S _hmin being equal or increasing 0.55 times the amount of the vertical stress S _V (S _hmin ≥ 0.55S _V ) and of the maximum horizontal stress S _Hmax ≤ 0.78–1.00S _V in stress regimes from normal to strike slip faulting. Thus, acting stresses in the 4,100-m deep reservoir are S _V  = 100 MPa, S _hmin = 55 MPa and S _Hmax = 78?100 MPa. Values from hydraulic fracturing support these results. Various fault sets of the reservoir are characterized in terms of their potential to conduct geothermal fluids based on their slip and dilatation tendency. This combined approach can be adopted to any other geothermal site investigation.     

High-pressure behavior of natural single-crystal epidote and clinozoisite up to 40 GPa

The comparative compressibility and high-pressure stability of a natural epidote (0.79 Fe-total per formula unit, Fe_tot pfu) and clinozoisite (0.40 Fe_tot pfu) were investigated by single-crystal X-ray diffraction and Raman spectroscopy. The lattice parameters of both phases exhibit continuous compression behavior up to 30 GPa without evidence of phase transformation. Pressure–volume data for both phases were fitted to a third-order Birch–Murnaghan equation of state with V ₀ = 461.1(1) ų, K ₀ = 115(2) GPa, and \(K_{0}^{'}\) = 3.7(2) for epidote and V ₀ = 457.8(1) ų, K ₀ = 142(3) GPa, and \(K_{0}^{'}\) = 5.2(4) for clinozoisite. In both epidote and clinozoisite, the b-axis is the stiffest direction, and the ratios of axial compressibility are 1.19:1.00:1.15 for epidote and 1.82:1.00:1.19 for clinozoisite. Whereas the compressibility of the a-axis is nearly the same for both phases, the b- and c-axes of the epidote are about 1.5 times more compressible than in clinozoisite, consistent with epidote having a lower bulk modulus. Raman spectra collected up to 40.4 GPa also show no indication of phase transformation and were used to obtain mode Grüneisen parameters (γ _i) for Si–O vibrations, which were found to be 0.5–0.8, typical for hydrous silicate minerals. The average pressure coefficient of Raman frequency shifts for M–O modes in epidote, 2.61(6) cm^?1/GPa, is larger than found for clinozoisite, 2.40(6) cm^?1/GPa, mainly due to the different compressibility of FeO₆ and AlO₆ octahedra in M3 sites. Epidote and clinozoisite contain about 2 wt% H₂O are thus potentially important carriers of water in subducted slabs.     

Clinopyroxenes from mantle-related xenocrysts in alkaline basalts from Hannuoba (China): augite–pigeonite exsolutions and their thermal significance

Optically homogeneous augite xenocrysts, closely associated with spinel–peridotite nodules, occur in alkali basalts from Hannuoba (Hebei province, China). They were studied by electron and X-ray diffraction to define the occurrence and significance of pigeonite exsolution microtextures. Sub-calcic augite (Wo₃₄) exsolved into En_62–62Fs_25–21Wo_13–17 pigeonite and En_46–45Fs_14–14Wo_40–42 augite, as revealed by TEM through diffuse coarser (001) lamellae (100–300 Å) and only incipient (100) thinner ones (<70 Å). C2/c augite and P2₁/c pigeonite lattices, measured by CCD-XRD, relate through a(Aug)?a(Pgt), b(Aug)?b(Pgt), c(Aug)≠c(Pgt) [5.278(1) vs 5.189(1)Å] and β(Aug)≠β(Pgt) [106.55(1) vs 108.55(2)°]. Cell and site volumes strongly support the hypothesis that the augite xenocrysts crystallised at mantle depth from alkaline melts. After the augite xenocrysts entered the magma, (001) lamellae first formed by spinodal decomposition at a T_min of about 1,100 °C, and coarsened during very rapid transport to the surface; in a later phase, possibly on cooling, incipient (100) lamellae then formed.