基于微结构-尺度双分解的页岩随机介质模拟和非均质特征分析

页岩储层矿物颗粒、孔/裂隙、干酪根等微观结构呈现明显的尺度化分布特征,常规的单结构单尺度随机介质模拟方法难以完整描述和重构微观尺度的页岩储层介质,本文提出了一种微结构-尺度双分解的随机介质模拟方法.基于龙马溪组页岩数字岩心,将岩心切片按照占比分解为脆性矿物、孔隙、干酪根及背景介质四种类型,对脆性矿物、孔隙和干酪根三种微结构进行尺度分解,通过优化随机介质模型参数,实现精确模拟不同尺度的微结构组分,再按占比进行微结构-尺度双合成.结果表明,微结构-尺度双分解随机介质模拟大幅度提高强非均质页岩储层介质的建模精度.     

Ultra‐high temperature metamorphism recorded in Fe‐rich olivine‐bearing migmatite from the Khondalite Belt,North China Craton

We document the first occurrence of Fe‐rich olivine‐bearing migmatitic metapelite in the Khondalite Belt, North China Craton. Petrological analyses revealed two exotic assemblages of orthopyroxene+spinel+olivine and orthopyroxene+spinel+cordierite. Phase relation modelling suggests that these assemblages are diagnostic of ultra‐high temperature (UHT) metamorphism in the Fe‐rich system, with temperatures from 1,000 to 1,050°C at 0.6 GPa. U–Th–Pb SIMS analyses on zircon reveal a similar age of c. 1.92 Ga for the olivine‐bearing migmatite and an adjacent gabbronoritic intrusion that is therefore identified as the heat source for the UHT metamorphism. These results, coupled with additional analysis of the famous Tuguiwula sapphirine‐bearing granulite, lead to a re‐appraisal of the P–T path shape and heat source for the UHT metamorphism. We suggest that UHT metamorphism, dated between 1.92 and 1.88 Ga, across the whole Khondalite belt, proceeded from a clockwise P–T evolution with an initial near‐isobaric heating path at ~0.6–0.8 GPa, and a maximum temperature of 1,050°C followed by a cooling path with minor decompression to ~0.5 GPa. Considering our results and previous works, we propose that the orogenic crust underwent partial melting at temperature reaching 850°C and depth of ~20 to ~30 km during a period of c. 30 Ma, between 1.93 and 1.90 Ga. During this time span, the partially molten crust was continuously intruded by mafic magma pulses responsible for local greater heat supply and UHT metamorphism above 1,000°C. We propose that the UHT metamorphism in the Khondalite belt is not related to an extensional post‐collisional event, but is rather syn‐orogenic and associated with mafic magma supplies.     

A non-hydrostatic algorithm for free-surface ocean modelling

An original implementation of a non-hydrostatic, free-surface algorithm based on a pressure correction method is proposed for ocean modelling. The free surface is implemented through an explicit scheme combined with a mode-spitting method but the depth-averaged velocity and the position of the free surface are updated at each non-hydrostatic iteration. The vertical momentum equation is also integrated up to the surface enabling a natural and accurate treatment of the surface layer. The consistent specification of the numerical schemes provides balanced transfers of potential and kinetic energy. This algorithm is well-suited for implementation as a non-hydrostatic kernel on originally hydrostatic free-surface ocean models such as Symphonie (http://poc.obs-mip.fr/pages/research_topics/modelling/symphonie/symphonie.htm) for which it has originally been developed.Energy balances associated with the propagation of short surface waves and solitary waves are presented for two dedicated well-documented configurations over closed domains. The buoyancy flux, the work rate of the pressure force together with the power of the advective terms are evaluated and discussed for the generation and the propagation of these two types of waves. The dissipation rate is in particular shown to be several orders of magnitude smaller than the work rates of the hydrostatic and non-hydrostatic pressure forces confirming the necessity for the exchanges of energy to be numerically balanced. The algorithm is subsequently applied to the complex generation of non-linear solitary internal waves by surface tides over Georges Bank, in the Gulf of Maine. The generation and the propagation of the observed non-linear and non-hydrostatic features in this region are correctly reproduced.     

Robe’s circular restricted three-body problem with a Roche ellipsoid-triaxial versus oblate system

This paper investigates Robe’s circular restricted three-body problem for two cases: with a Roche ellipsoid-triaxial system and with a Roche ellipsoid-oblate system. Without ignoring any component in both problems, a full treatment is given of the buoyancy force. The relevant equations of motion are established, and the special case where the density of the fluid and that of the infinitesimal mass are equal (D=0) is discussed. The location of the libration point and its stability when the infinitesimal mass is denser than the medium (D>0) are studied and it is found that the point (0,0,0) is the only libration point and this point is stable.     

A fast computation of the diurnal secondary ion production in the ionosphere of Titan

We propose an analytic model that allows rapid computation of the secondary ion production due to electron impact from the primary photo-production in the ionosphere of Titan. The model parameters are given for each of the 5 major ion productions (N⁺₂, CH⁺₄, N⁺, CH⁺₃, N⁺⁺₂) as well as for the electron production.     

PAMPRE: A dusty plasma experiment for Titan's tholins production and study

Organic aerosols play a significant role in the properties and evolution of Titan's atmosphere. But our knowledge of them and their physico-chemical mechanisms of formation and evolution are currently limited to a few data obtained by Titan observations from the Earth or from space probes. For this reason, laboratory experiments are developed to simulate the atmospheric chemistry and produce analogues of these aerosols in order to understand better their properties and how they are formed. The plasma discharges are the most efficient devices for the production of such analogues. However, the existing plasmas simulations introduce experimental biases compared with the conditions of aerosols production in Titan's atmosphere: chemistry is induced by electrons instead of photons; the solid analogues are produced and deposited on solid surfaces; direct analysis of the particles inside the reactive chamber is not easy. In order to avoid some of these experimental problems, we have developed another method of production of Titan's aerosols analogues. It is based on a capacitively coupled radio-frequency (RF) cold plasma system at low pressure in a N₂-CH₄ gaseous mixture. In this plasma, solid particles produced from the gas phase are in levitation, thus preventing any wall effect on their production, and allowing the study of the formation and growth of the particles directly in the plasma. Moreover, the electron energy distribution of this plasma can be compared with the solar spectrum. This article describes the RF plasma experiment and presents the first results obtained with an initial N₂-CH₄ (90:10) gaseous mixture which produced our first studied analogues of Titan's aerosols.     

Change in carbon flux (1960–2015) of the Red River (Vietnam)

Global riverine carbon concentrations and fluxes have been impacted by climate and human-induced changes for many decades. This paper aims to reconstruct the longterm carbon concentrations and carbon fluxes of the Red River, a system under the coupled pressures of environmental change and human activity. Based on (1) the relationships between particulate and dissolved organic carbon (POC, DOC) or dissolved inorganic carbon (DIC), and suspended sediments (TSS) or river water discharge and on (2) the available detailed historical records of river discharge and TSS concentration, the variations of the Red River carbon concentration and flux were estimated for the period 1960–2015. The results show that total carbon flux of the Red River averaged 2555?±?639 kton C year^?1. DIC fluxes dominated total carbon fluxes, representing 64% of total, reflecting a strong weathering process from carbonate rocks in the upstream basin. Total carbon fluxes significantly decreased from 2816 kton C year^?1 during the 1960s to 1372 kton C year^?1 during the 2010s and showed clear seasonal and spatial variations. Organic carbon flux decreased in both quantity and proportion of the total carbon flux from 40.9% in 1960s to 14.9% in 2010s, reflecting the important impact of dam impoundment. DIC flux was also reduced over this period potentially as a consequence of carbonate precipitation in the irrigated, agricultural land and the reduction of the Red River water discharge toward the sea. These decreases in TSS and carbon fluxes are probably partially responsible for different negatives impacts observed in the coastal zone.     

A compound Ca‐, Al‐rich inclusion from CV3 chondrite Northwest Africa 3118: Implications for understanding processes during CAI formation

A calcium‐aluminum‐rich inclusion 3N from the Northwest Africa (NWA) 3118 CV3 carbonaceous chondrite is a unique cm‐sized compound object, primarily a forsterite‐bearing type B (FoB) CAI, that encloses at least 26 smaller CAIs of different types, including compact type A (CTA), B, C, and an ultra‐refractory inclusion. Relative to typical type A and B CAIs found elsewhere, the bulk compositions of the types A and B CAIs within 3N more closely match the bulk compositions predicted by equilibrium condensation of a gas of solar composition. Being trapped within the FoB melt may have protected them from melt evaporation that affected most “stand‐alone” CAIs. 3N originated either as an aggregate of many smaller (mostly types A, B, C) CAIs plus accreted Fo‐bearing material (like an amoeboid olivine aggregate) which experienced partial melting of the whole, or else as a FoB melt droplet that collided with and trapped many smaller solid CAIs. In the former case, 3N recorded the earliest accretion of pebble‐sized bodies known. In the latter case, the presence of a large number of individual refractory inclusions within 3N suggests a very high local density of refractory solids in the immediate region of the host CAI during the brief time while it was melted. Collisions would have occurred on time scales of hours at most, assuming a melt solidification interval for the host CAI of 300–400 °C (maximum) and a cooling rate of ~10 °C/h.     

The Piànico tephra: an early Middle Pleistocene record of intraplate volcanism in the Mediterranean

ABSTRACT Distal tephra are a valuable record of the volcano-tectonic evolution of an area under study. Here, we document the case of the Early Middle Pleistocene rhyolitic tephra of Piànico, discovered in the Southern Italian Alps. The geochemical characteristics of Piànico are unique among the distal tephra outcropping in Italy and indicate an intraplate volcanism at the source. The alkali composition and trace elements show a striking resemblance with the rhyolitic complex of the Euganean Hills, located 170 km from Piànico. However, these rhyolites are much older (Oligocene). Alternatively, the source of this intraplate volcanic episode could be located in the southern Tyrrhenian Sea.