A micro–macro investigation of the capillary strengthening effect in wet granular materials

Wet granular materials are three-dimensionally simulated by the discrete element method with water bridges incorporated between particles. The water bridges are simplified as toroidal shapes, and the matric suction is constantly maintained in the material. A comparison with experimental tests in the literature indicates that the toroidal shape approximation may be one of the best choices with high practicability and decent accuracy. Mechanical behaviours of wet granular materials are studied by triaxial tests. Effects of particle size distributions and void ratios are investigated systematically in this study. The hydraulic limit of the pendular state is also discussed. It gives the capillary cohesion function which is not only determined by the degree of saturation but also positively correlated to relative density and particle size polydispersity and inversely proportional to mean particle size. Furthermore, the capillary strengthening effect is also analysed microscopically in aid of the Stress–Force–Fabric relationship, mainly in fabric anisotropy, coordination number and stress transmission pattern, which revealed the micro-mechanisms of the additional effective stress induced by capillary effect.     

An adjoint-based optimization method for gas production in shale reservoirs

In this work, we consider a new model for flow in a multiporosity shale gas reservoir constructed within the framework of an upscaling procedure where hydraulic fractures are treated as (\(n-1\)) interfaces (\(n=2,3\)). Within this framework, the hydrodynamics is governed by a new pressure equation in the shale matrix which is treated as a homogenized porous medium composed of organic matter (kerogen aggregates with nanopores) and inorganic impermeable solid (clay, calcite, quartz) separated from each other by a network of interparticle pores of micrometer size. The solution of the pressure equation is strongly influenced by the constitutive response of the retardation parameter and effective hydraulic conductivity where the former incorporates gas adsorption/desorption in the nanopores of the kerogen. By focusing our analyses on this nonlinear diffusion equation in the domain occupied by the shale matrix, an optimization strategy seated on the adjoint sensitivity method is developed to minimize a cost functional related to gas production and net present value in a single hydraulic fracture. The gradient of the objective functional computed with the adjoint formulation is explored to update the controlled pressure drop aiming to optimize production in a given window of time. The combination of the direct approach and gradient-based optimization using the adjoint formulation leads to the construction of optimal production scenarios under controlled pressure decline in the well. Numerical simulations illustrate the potential of the methodology proposed herein in optimizing gas production.     

Oil geochemistry derived from the Qinjiatun–Qikeshu oilfields: insight from light hydrocarbons

A suite of 27 oils from the Qinjiatun–Qikeshu oilfields in the Lishu Fault Depression of the Songliao Basin was analyzed using whole oil gas chromatography. In combination with the relative distribution of C₂₇, C₂₈, and C₂₉ regular steranes, detailed geochemical analyses of light hydrocarbons in oil samples revealed crude oils characterized by the dual input of lower aquatic organisms and higher terrestrial plants. Several light hydrocarbon indicators suggest that the liquid hydrocarbons have maturities equivalent to vitrinite reflectances of around 0.78%–0.93%. This is consistent with the maturity determination of steranes C₂₉ 20S/(20S + 20R) and C₂₉ ααβ/(ααα + αββ). Crude oils derived from the two distinct oilfields likely both have source rocks deposited in a lacustrine environment based on light hydrocarbon parameters and on higher molecular weight hydrocarbon parameters. The results show that light hydrocarbon data in crude oils can provide important information for understanding the geochemical characteristics of the Qinjiatun–Qikeshu oils during geologic evolution.     

Distribution and geochemical significance of trace elements in shale rocks and their residual kerogens

There is a dearth of information about the distribution of trace elements in kerogen from shale rocks despite several reports on trace element composition in many shale samples. In this study, trace elements in shale rocks and their residual kerogens were determined by inductively coupled plasma–mass spectrometry. The results from this study show redox-sensitive elements relatively concentrated in the kerogens as compared to the shales. This may be primarily due to the adsorption and complexation ability of kerogen, which enables enrichment in Ni, Co, Cu, and Zn. For the rare earth elements (REEs), distinct distribution characteristics were observed for shales dominated by terrigenous minerals and their kerogen counterparts. However, shales with less input of terrigenous minerals showed similar REE distribution patterns to their residual kerogen. It is speculated that the distribution patterns of the REEs in shales and kerogens may be source-related.     

Comparison of AMSR-2 wind speed and sea surface temperature with moored buoy observations over the Northern Indian Ocean

The Northern Indian Ocean (NIO) is unique due to seasonal reversal of wind patterns, the formation of vortices and eddies which make satellite observations arduous. The veracity of sea surface wind (SSW) and sea surface temperature (SST) products of sun-synchronous AMSR-2 satellite are compared with high-temporal moored buoy observations over the NIO. The two year-long (2013–2014) comparisons reveal that the root-mean-square-error (RMSE) of AMSR-2 SST and SSW is \(<0.4{^{\circ }}\hbox {C}\) and \(<1.5\hbox { ms}^{-1}\), respectively, which are within the error range prescribed for the AMSR-2 satellite (\(\pm 0.8{^{\circ }}\hbox {C}\), \(\pm 1.5\hbox { ms}^{-1})\). The SST–wind relation is analyzed using data both from the buoy and satellite. As a result, the low-SST is associated with low-wind condition (positive slope) in the northern part of the Bay of Bengal (BoB), while low SST values are associated with high wind conditions (negative slope) over the southern BoB. Moreover, the AMSR-2 displayed larger slope for SST–wind relation and could be mainly due to overestimation of SST and underestimation of wind as compared to the buoy. The AMSR-2 SSW exhibited higher error during post-monsoon followed by monsoon season and could be attributed to the high wind conditions associated with intense oceanic vortices. The study suggests that the AMSR-2 products are reliable and can be used in tropical air–sea interactions, meso-scale features, and weather and climate studies.     

Early Permian transgressive–regressive cycles: Sequence stratigraphic reappraisal of the coal-bearing Barakar Formation,Raniganj Basin,India

The present research is an attempt to assess the Barakar Formation of the Raniganj Gondwana Basin, India, in the frame of fluvio-marine (estuarine) depositional systems using sequence stratigraphic elements. Analysis of predominant facies associations signify deposition in three sub-environments: (i) a river-dominated bay-head delta zone in the inner estuary, with transition from braided fluvial channels (FA-B1) to tide-affected meandering fluvial channels and flood plains (FA-B2) in the basal part of the succession; (ii) a mixed energy central basin zone, which consists of transitional fluvio-tidal channels (FA-B2), tidal flats, associated with tidal channels and bars (FA-B3) in the middle-upper part of the succession; and (iii) a wave-dominated outer estuary (coastal) zone (FA-B4 with FA-B3) in the upper part of the succession. Stacked progradational (P1, P2)–retrogradational (R1, R2) successions attest to one major base level fluctuation, leading to distinct transgressive–regressive (T–R) cycles with development of initial falling stage systems tract (FSST), followed by lowstand systems tract (LST) and successive transgressive systems tracts (TST-1 and TST-2). Shift in the depositional regime from regressive to transgressive estuarine system in the early Permian Barakar Formation is attributed to change in accommodation space caused by mutual interactions of (i) base level fluctuations in response to climatic amelioration and (ii) basinal tectonisms (exhumation/sagging) related to post-glacial isostatic adjustments in the riftogenic Gondwana basins.     

Rootless calc-silicate folds in granite: An implication towards syn- to post-plutonic emplacement

Deformation of the Champaner Group of rocks that form a part of Southern Aravalli Mountain Belt, western India, occurred during the Grenville orogeny (ca. 1400–935 Ma). Two phases of deformation are recorded: \(\hbox {D}_{1}\), persistent throughout the group and characterised by westerly plunging tight isoclinal folds and \(\hbox {D}_{2}\), a localized phase of deformation associated with shortening of the earlier folds from the eastern margin. Both the phases of deformation are in association with the syn-tectonically emplaced Godhra granite. The present work records rootless calc-silicate folds in granite belonging to the older formation, located at the eastern fringe of the Champaner Group. Field evidences suggest superimposition of Type 2 interference pattern trending NE–SW over rootless Type 0 of varying trends from NW–SE to N–S. The superposed pattern obtained from the field study differs in terms of structural trends with the neighbouring Precambrian stratigraphic units. These stratigraphic units include the Champaner Group to which the study area belongs, the Kadana Formation of the Lunavada Group and Pre-Chamapaner Gneissic Complex. Rootless character of folds found within the study area imply syn-post plutonic emplacement of Godhra granite.     

Tectonic significance of Triassic mafic rocks in the June Complex,Sanandaj–Sirjan zone,Iran

This study concentrates on the petrological and geochemical investigation of mafic rocks embedded within the voluminous Triassic June Complex of the central Sanandaj–Sirjan zone (Iran), which are crucial to reconstruct the geodynamics of the Neotethyan passive margin. The Triassic mafic rocks are alkaline to sub-alkaline basalts, containing 43.36–49.09 wt% SiO₂, 5.19–20.61 wt% MgO and 0.66–4.59 wt% total alkalis. Based on MgO concentrations, the mafic rocks fall into two groups: cumulates (Mg# = 51.61–58.94) and isotropic basaltic liquids (Mg# = 24.54–42.66). In all samples, the chondrite-normalized REE patterns show enrichment of light REEs with variable (La/Yb)_N ratios ranging from 2.48 to 9.00, which confirm their amalgamated OIB-like and E-MORB-like signatures. Enrichment in large-ion lithophile elements and depletion in high field strength elements (HFSE) relative to the primitive mantle further support this interpretation. No samples point to crustal contamination, all having undergone fractionation of olivine + clinopyroxene + plagioclase. Nevertheless, elemental data suggest that the substantial variations in (La/Sm)_PM and Zr/Nb ratios can be explained by variable degrees of partial melting rather than fractional crystallization from a common parental magma. The high (Nb/Yb)_PM ratio in the alkaline mafic rocks points to the mixing of magmas from enriched and depleted mantle sources. Abundant OIB alkaline basalts and rare E-MORB appear to be linked to the drifting stage on the northern passive margin of the Neotethys Ocean.     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.     

Flooding a landscape: impact of Holocene transgression on coastal sedimentology and underwater archaeology in Kiladha Bay (Greece)

Franchthi Cave, bordering Kiladha Bay, in Greece, is a key archaeological site, due to its long occupation time, from?~?40,000 to?~?5000 year BP. To date, no clear evidence of Neolithic human dwellings in the cave was found, supporting the assumption that Neolithic people may have built a village where there is now Kiladha Bay. During the Neolithic period/Early Holocene, wide areas of the bay were indeed emerged above sea level. Bathymetric and seismic data identified a terrace incised by a valley in?~?1 to 2 m sediment depth. Eight sediment cores, up to 6.3-m-long, were retrieved and analysed using petrophysical, sedimentological, geochemical, and chronostratigraphic methods. The longest core extends into the exposure surface, consisting of a layer of carbonate rubble in a finer matrix, representing weathering processes. Dated organic remains place this unit at?~?8500 cal year BP. It is overlain by stiff silty mud representing an estuarine environment. This mud is capped by reduced sediments with roots marking an exposure surface. A shell-layer, dated to?~?6300 cal year BP, overlies this terrestrial sequence, reflecting the marine transgression. This layer occurs at 10.8 mbsl, 7.7 m deeper than the global sea level at that time, suggesting tectonic subsidence in the area. It is overlain by finer-grained marine carbonate-rich sediments. The top of the core shows traces of eutrophication, pebbles and marine shells, all likely a result of modern anthropogenic processes. These results are interpreted in the context of human occupation: the exposed surface contains pottery sherds, one dating to the Early to Middle Neolithic period, indicating that Neolithic people were present in this dynamic landscape interacting with a migrating coastline. Even if the artefacts are isolated, future investigations of the submerged landscape off Franchthi Cave might lead to the discovery of a Neolithic village, which eventually became buried under marine sediments.