Contribution of long-term hydrothermal experiments for understanding the smectite-to-chlorite conversion in geological environments

The smectite-to-chlorite conversion is investigated through long-duration experiments (up to 9 years) conducted at 300 °C. The starting products were the Wyoming bentonite MX80 (79 % smectite), metallic iron and magnetite in contact with a Na–Ca chloride solution. The predominant minerals in the run products were an iron-rich chlorite (chamosite like) and interstratified clays interpreted to be chlorite/smectite and/or corrensite/smectite, accompanied by euhedral crystals of quartz, albite and zeolite. The formation of pure corrensite was not observed in the long-duration experiments. The conversion of smectite into chlorite over time appears to take place in several steps and through several successive mechanisms: a solid-state transformation, significant dissolution of the smectite and direct precipitation from the solution, which is over-saturated with respect to chlorite, allowing the formation of a chamosite-like mineral. The reaction mechanisms are confirmed by X-ray patterns and data obtained on the experimental solutions (pH, contents of Si, Mg, Na and Ca). Because of the availability of some nutrients in the solution, total dissolution of the starting smectite does not lead to 100 % crystallization of chlorite but to a mixture of two dominant clays: chamosite and interstratified chlorite/smectite and/or corrensite/smectite poor in smectite. The role of Fe/(Fe + Mg) in the experimental medium is highlighted by chemical data obtained on newly formed clay particles alongside previously published data. The newly formed iron-rich chlorite has the same composition as that predicted by the geothermometer for diagenetic to low-grade metamorphic conditions, and the quartz + Fe-chlorite + albite experimental assemblage in the 9-year experiment is close to that fixed by water–rock equilibrium.     

Fast wavelet-based stochastic simulation using training images

Spatial uncertainty modelling is a complex and challenging job for orebody modelling in mining, reservoir characterization in petroleum, and contamination modelling in air and water. Stochastic simulation algorithms are popular methods for such modelling. In this paper, discrete wavelet transformation (DWT)-based multiple point simulation algorithm for continuous variable is proposed that handles multi-scale spatial characteristics in datasets and training images. The DWT of a training image provides multi-scale high-frequency wavelet images and one low-frequency scaling image at the coarsest scale. The simulation of the proposed approach is performed on the frequency (wavelet) domain where the scaling image and wavelet images across the scale are simulated jointly. The inverse DWT reconstructs simulated realizations of an attribute of interest in the space domain. An automatic scale-selection algorithm using dominant mode difference is applied for the selection of the optimal scale of wavelet decomposition. The proposed algorithm reduces the computational time required for simulating large domain as compared to spatial domain multi-point simulation algorithm. The algorithm is tested with an exhaustive dataset using conditional and unconditional simulation in two- and three-dimensional fluvial reservoir and mining blasted rock data. The realizations generated by the proposed algorithm perform well and reproduce the statistics of the training image. The study conducted comparing the spatial domain filtersim multiple-point simulation algorithm suggests that the proposed algorithm generates equally good realizations at lower computational cost.     

A coupled model for heating and hydratation in unsaturated clays

Knowledge of transport processes of heat and moisture in soils of arid zones is vital to understanding the environmental and economic impacts of many activities: agriculture, waste disposal, geoenvironmental practices and earth sciences. Through extensive review and study on the different aspects of coupled transfer processes in swelling porous media, a general mathematical model for coupled heat, moisture, air flow and deformation problems in clayey soils is proposed in a consistent and unified manner. The model is characterized by the presence of a deformable solid matrix filled with two fluid phases (liquid water and air). In the proposed model, both pore water and air transfers are assumed to be governed by the generalized Darcy’s law. Fully coupled, non-linear partial differential equations are established and then solved by using a Galerkin weighted residual approach in space domain and an implicit integrating scheme in time domain. The obtained model has been finally validated by means of some case tests for the prediction of the thermo-hydro-mechanical behaviour of unsaturated swelling soils. The calculated relative errors between experimental and numerical results are 3% for temperature and 7% for stresses. Consequently, the developed numerical model predicts satisfactory results, compared to experimental test measures. The model is applicable to two-dimensional problems with various initial and boundary conditions; non-linear soil parameters can be easily included in this model.     

Sexual reproduction of Hippospongia communis (Lamarck, 1814) (Dictyoceratida,Demospongiae): comparison of two populations living under contrasting environmental conditions

The study of the reproductive processes of benthic invertebrates is essential to the understanding of their population dynamics and is also important in formulating conservation plans, especially for exploited species. The sexual reproduction of Hippospongia communis, the ‘honeycomb’ bath sponge, was studied at two locations in the Mediterranean Sea: the Kerkennah Islands (Tunisia, South Mediterranean), where the mean annual seawater temperature is 19 °C, and Marseille (France, Northwestern Mediterranean Sea), where the mean annual water temperature is 16 °C. The aim of this comparative study was to determine whether different environmental conditions could affect reproduction patterns. At both locations, H. communis was found to contain sexual reproductive elements year‐round. Oogenesis and embryogenesis occurred throughout the year, whereas spermatogenesis occurred during shorter periods between October and November, in both populations. While gametogenesis seemed to be synchronized, indicating that fertilization could occur at the same time at both locations, spawning was observed in late summer in Marseille, whereas it started in late spring for the Kerkennah population. Larval development of H. communis seems to take longer for sponges living at cooler locations such as Marseille. Reproductive effort calculated for both sexes showed significantly higher values for specimens from Kerkennah Islands. By comparing sexual reproductive patterns of populations living in two contrasted environments, we suggest that a change of thermal regime can affect H. communis phenology.     

Active faulting in the western Pyrénées (France): Paleoseismic evidence for late Holocene ruptures

We have identified a 50-km-long active fault scarp, called herewith the Lourdes Fault, between the city of Lourdes and Arette village in the French Pyrénées. This region was affected by large and moderate earthquakes in 1660 (Io = VIII–IX, MSK 64,), in 1750 (Io = VIII, MSK 64) and in 1967 (M_d = 5.3, Io = VIII, MSK 64). Most earthquakes in this area are shallow and the few available focal mechanism solutions do not indicate a consistent pattern of active deformation. Field investigations in active tectonics indicate an East–West trending and up to 50-m-high fault scarp, in average, made of 3 contiguous linear fault sub-segments. To the north, the fault controls Quaternary basins and shows uplifted and tilted alluvial terraces. Deviated and abandoned stream channels of the southern block are likely due to the successive uplift of the northern block of the fault. Paleoseismic investigations coupled with geomorphic studies, georadar prospecting and trenching along the fault scarp illustrate the cumulative fault movements during the late Holocene. Trenches exhibit shear contacts with flexural slip faulting and thrust ruptures showing deformed alluvial units in buried channels. ¹⁴C dating of alluvial and colluvial units indicates a consistent age bracket from two different trenches and shows that the most recent fault movements occurred between 4221 BC and 2918 BC. Fault parameters and paleoseismic results imply that the Lourdes Fault and related sub-segments may produce a M_W 6.5 to 7.1 earthquake. Fault parameters imply that the Lourdes Fault segment corresponds to a major seismic source in the western Pyrénées that may generate earthquakes possibly larger than the 1660 historical event.     

Analysis of subdiffusion in disordered and fractured media using a Grünwald-Letnikov fractional calculus model

The increasing applications of fractional calculus in simulating the anomalous transport behavior in disordered and fractured heterogeneous porous media has grown rapidly over the past decade. In the present study, a temporal fractional flux relationship is employed as a constitutive equation to relate the volumetric flow rate to the gradient of the pore pressure. The novelty of this paper entails interpreting the time fractional derivative operator in the flux relationship by the Grünwald-Letnikov (G-L) definition as opposed to the Caputo interpretation which has been widely considered. Subsequently, a numerical scheme based on the block-centered finite-difference discretization is formulated to handle the resulting non-linear fractional diffusion model. In addition, a linear stability analysis is successfully performed to establish the stability criterion of the developed numerical scheme. An expression for the modified incremental material balance index was derived to assess the effectiveness of the numerical discretization process. Finally, numerical experiments were performed to provide qualitative insights into the nature of pressure evolution in a hydrocarbon reservoir under the influence subdiffusion. In summary, the results establish that subdiffusion regime results in the development of higher pressure drop in the reservoir. This paper will provide a strong foundation for researchers interested in investigating anomalous diffusion phenomena in porous media.     

RETRACTED ARTICLE: Numerical analysis of multiphase flow in porous material

Recent developments in the application of x-ray micro-tomography in laboratory geomechanics have allowed all the individual grains of sand in a test sample to be seen and identified uniquely in 3D. Combining such imaging capabilities with experiments carried out “in situ” within an imaging set-up has led to the possibility of directly observing the mechanisms of deformation as they happen. The challenge has thus become extracting pertinent, quantified information from these rich time-lapse 3D images to elucidate the mechanics at play. This paper presents a new approach (ID-Track) for the quantification of individual grain kinematics (displacements and rotations) of large quantities of sand grains (tens of thousands) in a test sample undergoing loading. With ID-Track, grains are tracked between images based on some geometrical feature(s) that allow their unique identification and matching between images. This differs from Digital Image Correlation (DIC), which makes measurements by recognising patterns between images. Since ID-Track does not use the image of a grain for tracking, it is significantly faster than DIC. The technique is detailed in the paper, and is shown to be fast and simple, giving good measurements of displacements, but suffering in the measurement of rotations when compared with Discrete DIC. Subsequently, results are presented from successful applications of ID-track to triaxial tests on two quite different sands: the angular Hostun sand and the rounded Caicos Ooids. This reveals details on the performance of the technique for different grain shapes and insight into the differences in the grain-scale mechanisms occurring in these two sands as they exhibit strain localisation under triaxial loading.