Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Ice-Core Study of the Link between Sea-Salt Aerosol,Sea-Ice Cover and Climate in the Antarctic Peninsula Area

Three ice cores and a set of snow pit samples collected on James Ross Island, Antarctic Peninsula, in 1979, 1981 and 1991 have been analyzed for water stable isotope content D or ¹⁸O (isotopic temperature) and major chemical species. A reliable and detailed chronological scale has been established first for the upper 24.5 m of water equivalent (1990–1943) where various data sets can be compared, then extended down to 59.5 m of water equivalent (1847) with the aid of seasonal variations and the sulphate peak reflecting the 1883 Krakatoa volcanic eruption. At James Ross Island, sea-salt aerosol is generally produced by ice-free marine surfaces during the summer months, although some winter sea-salt events have been observed. For the upper part of the core (1990–1943), correlations (positive or negative) were calculated between isotopic temperature, chloride content (a sea-salt indicator), sea-ice extent, regional atmospheric temperature changes and atmospheric circulation. The D and chloride content correlation was then extended back to 1847, making it possible to estimate decadal sea-ice cover fluctuations over the study period. Our findings suggest that ice-core records from James Ross Island reflect the recent warming and sea-ice decrease trends observed in the Antarctic Peninsula area from the mid-1940s.     

Petrography of alluvial sands as a past and present environmental indicator: Case of the Loire River (France)

Four sample sets of the Upper and Middle Loire river sands were analyzed in order to study the impact of natural and anthropogenic factors on their petrographic composition in space (on an 800 km stretch) and time. Composition was determined by modal analysis of three sand-size fractions using a polarizing optical microscope and calculated for each sample (“standard sand” = Sst). The watershed is composed mainly of endogenic (Massif Central) and sedimentary (southern Parisian Basin) rocks. B-set sands collected in channels for different water flows in 1996 show that Sst compositions vary by only 5 %. Present-day sands in the Upper Loire and Middle Loire have very high petrographic immaturity comparing to others worldwide fluvial sands, although bio-climatic conditions favor sand maturation by source-rock weathering in the watershed. This shows the strong impact of the Massif Central on sediment yield due to relief rejuvenation as a consequence of the formation of the Alps during the Quaternary. Fluvial sands stored during the Weichselian and the Holocene in the Middle Loire floodplain, although partly weathered since their deposition, show higher inputs from the endogenic rocks of the Massif Central than present-day deposits. This can be explained by Weichselian periglacial conditions and the development of crop farming since the Neolithic, which favored mechanical erosion, particularly in the Massif Central which is characterized by a cold, humid climate and steep slopes. The upstream-downstream change in the composition of presently deposited sand is low in the diked area. It shows however that basalt and some heavy mineral grains are vulnerable to abrasion during transport and indicates a marked sediment yield from ancient sediment stored in the floodplain. This is in line with the high incision of the river bed over the last 150 years partly due to dam construction and aggregate mining.     

云南澜沧-耿马Ms7.6地震的完全库仑破裂应力变化与后续地震的动态、静态应力触发

计算和研究了1988年云南澜沧-耿马Ms7.6地震产生的完全库仑破裂应力变化的时空演化图像, 对继澜沧地震13 min后发生的耿马M7.2地震的应力触发问题进行了分析, 同时也对澜沧-耿马地震后24 d内发生的Ms5.0~6.9后续强余震的应力触发问题进行了探讨. 结果显示, 澜沧地震断层破裂产生的完全库仑破裂应力变化空间分布图像具有很强的非对称性, 正值的动态和静态库仑破裂应力变化区域均具有与强余震分布位置吻合较好的现象. 耿马Ms7.2地震受到了澜沧Ms7.6地震产生动态和静态库仑破裂应力的触发作用; 绝大多数后续强余震受到了动、静态库仑破裂应力的综合触发作用.     

Using slowness and azimuth fluctuations as new observables for four‐dimensional reservoir seismic monitoring

Four‐dimensional imaging using geophysical data is of increasing interest in the oil and gas industries. While travel‐time and amplitude variations are commonly used to monitor reservoir properties at depth, their interpretation can suffer from a lack of information to decipher the parts played by different parameters. In this context, this study focuses on the slowness and azimuth angle measured at the surface using source and receiver arrays as complementary observables. In the first step, array processing techniques are used to extract both azimuth and incidence angles at the source side (departure angles) and at the receiver side (arrival angles). In the second step, the slowness and angle variations are monitored in a laboratory environment. These new observables are compared with traditional arrival‐time variations when the propagation medium is subject to temperature fluctuations. Finally, field data from a heavy‐oil permanent reservoir monitoring system installed onshore and facing steam injection and temperature variations are investigated. The slowness variations are computed over a period of 152 days. In agreement with Fermat's principle, strong correlations between the slowness and arrival‐time variations are highlighted, as well as good consistency with other techniques and field pressure measurements. Although the temporal variations of slowness and arrival time show the same features, there are still differences that can be considered for further characterization of the physical changes at depth.     

The OSIRIS‐REx target asteroid (101955) Bennu: Constraints on its physical,geological, and dynamical nature from astronomical observations

We review the results of an extensive campaign to determine the physical, geological, and dynamical properties of asteroid (101955) Bennu. This investigation provides information on the orbit, shape, mass, rotation state, radar response, photometric, spectroscopic, thermal, regolith, and environmental properties of Bennu. We combine these data with cosmochemical and dynamical models to develop a hypothetical timeline for Bennu's formation and evolution. We infer that Bennu is an ancient object that has witnessed over 4.5 Gyr of solar system history. Its chemistry and mineralogy were established within the first 10 Myr of the solar system. It likely originated as a discrete asteroid in the inner Main Belt approximately 0.7–2 Gyr ago as a fragment from the catastrophic disruption of a large (approximately 100‐km), carbonaceous asteroid. It was delivered to near‐Earth space via a combination of Yarkovsky‐induced drift and interaction with giant‐planet resonances. During its journey, YORP processes and planetary close encounters modified Bennu's spin state, potentially reshaping and resurfacing the asteroid. We also review work on Bennu's future dynamical evolution and constrain its ultimate fate. It is one of the most Potentially Hazardous Asteroids with an approximately 1‐in‐2700 chance of impacting the Earth in the late 22nd century. It will most likely end its dynamical life by falling into the Sun. The highest probability for a planetary impact is with Venus, followed by the Earth. There is a chance that Bennu will be ejected from the inner solar system after a close encounter with Jupiter. OSIRIS‐REx will return samples from the surface of this intriguing asteroid in September 2023.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.