Lattice-Boltzmann simulations of the capillary pressure–saturation–interfacial area relationship for porous media

Hysteresis in the relationship between capillary pressure (P_c), wetting phase saturation (S_w) and nonwetting–wetting interfacial area per volume (a_nw) is investigated using multiphase lattice-Boltzmann simulations of drainage and imbibition in a glass bead porous system. In order to validate the simulations, the P_c–S_w and a_nw–S_w main hysteresis loops were compared to experimental data reported by Culligan et al. [Culligan KA, Wildenschild D, Christensen BS, Gray WG, Rivers ML, Tompson AB. Interfacial area measurements for unsaturated flow through porous media. Water Resour Res 2004;40:W12413]. In general, the comparison shows that the simulations are reliable and capture the important physical processes in the experimental system. P_c–S_w curves, a_nw–S_w curves and phase distributions (within the pores) show good agreement during drainage, but less satisfactory agreement during imbibition. Drainage and imbibition scanning curves were simulated in order to construct P_c–S_w–a_nw surfaces. The root mean squared error (RMSE) and mean absolute error (MAE) between drainage and imbibition surfaces was 0.10 mm⁻¹ and 0.03 mm⁻¹, respectively. This small difference indicates that hysteresis is virtually nonexistent in the P_c–S_w–a_nw relationship for the multiphase system studied here. Additionally, a surface was fit to the main loop (excluding scanning curves) of the drainage and imbibition P_c–S_w–a_nw data and compared to the surface fit to all of the data. The differences between these two surfaces were small (RMSE = 0.05 mm⁻¹ and MAE = 0.01 mm⁻¹) indicating that the P_c–S_w–a_nw surface is adequately represented without the need for the scanning curve data, which greatly reduces the amount of data required to construct the non-hysteretic P_c–S_w–a_nw surface for this data.     

Novel insights on intensity and typology of direct human-nature interactions in protected areas through passive crowdsourcing

Recent advances in geotagging, sharing and automatically analyzing online content from Social Networking Sites (SNS) offer unprecedented opportunities for the analysis of human-nature interactions. Previous studies in this field, however, offer limited insights regarding the benefits of automated content analysis especially at large scales, biases arising from the selection of SNS sources, and the predictive power of visitation models based on SNS data. We explore quantitative and qualitative aspects related to intensity, interests and sentiments associated with on-site experiences in 568 protected areas in Israel and the Palestinian Authority. We analyze counts and content of >100,000 photographs and tweets from four different SNSs, calibrate visitation models and predict visitation in unmonitored sites, cluster sites based on the typology of human-nature interactions reflected in online photographs, and characterize the polarity of sentiments associated with experiences in individual sites and clusters thereof. We find benefits in combining data from multiple sources and controlling for biases related to sites’ photogenicity and type of human-nature interactions. Our results suggest that current best estimates of visitation in unmonitored sites underestimate by 39% the actual number of visits. We discuss how the techniques and findings in this study are applicable in the broader context of the management and conservation of sites of environmental or cultural interest.     

关于波浪作用下粉土质海床中共振现象的解释

本文详细讨论了粉土质床在波浪作用下所呈现共振现象的力学机制。利用有限深海床下波浪载荷对土层响应的研究，系统地比较了理论计算结果和实验结果。提出在波浪作用下粉土质土层中存在着一硬土层，从而引发Ｂｒａｇｇ共振。     

Seismic Wave Propagation and Excitation in Multi-layered Media with Irregular Interfaces Part (Ⅰ): SH Case

Study of seismic wave excitation and propagation in laterally heterogeneous media was an active and important subject in seismology in the past two decades, numerous analytical and numerical efforts have been made in this research field. In this article, I have, first, made a brief review on those developments and then introduced and summarized a unified and efficient method, global generalized reflection-transmission (abbreviated to R/T thereafter) matrices method, for synthetic seismograms in multi-layered media with irregular interfaces developed by the author [24~26]. As demonstrated in this article, this method could be regarded as an extension of the generalized R/T coefficients method for the horizontally layered case [2,5] to the layered media with irregularly shaped interfaces by incorporating the T matrices technique [27,28]. Because of the use of a recursive scheme in computing the global generalized R/T matrices, this method is efficient, particularly for the case with a large number of irre     

考虑土体固液二相性质的地下管线地震反应研究

本文将管线周围的土体视为固-液二相介质，采用考虑地震过程中土体内孔隙水压力的增长和消散的有效应力分析方法，并以非线性的土体本构模型模拟土体动力特性的改变以及在孔隙水压作用下的变形，以求接近土体的真实地震反应。研究表明随着土体软化程度的不同，受到土体约束的管线应力亦发生相应变化。文中详细研究了地震动峰值、管线直径、行波波速和入射角度等对管线应力的影响。     

A semi-empirical approach to model pressure dependence of elastic moduli in granular media accounting for variations of coordination-number and Poisson-ratio

The effective medium theory based on the Hertz–Mindlin contact law is the most popular theory to relate dynamic elastic moduli (or elastic velocities) and confining pressure in dry granular media. However, many experimental results proved that the effective medium theory predicts pressure trends lower than experimental ones and over-predicts the shear modulus. To mitigate these mispredictions, several evolutions of the effective medium theory have been presented in the literature. Among these, the model named modified grain contact theory is an empirical approach in which three parametric curves are included in the effective medium theory model. Fitting the parameters of these curves permits to adjust the pressure trends of the Poisson ratio and the bulk modulus. In this paper, we present two variations of the modified grain contact theory model. First, we propose a minor modification in the fitting function for the porosity dependence of the calibration parameters that accounts for non-linearity in the vicinity of the critical porosity. Second, we propose a major modification that reduces the three-step modified grain contact theory model to a two-step model, by skipping the calibration parameter–porosity fit in the model and directly modelling the calibration parameter–pressure relation. In addition to an increased simplicity (the fitting parameters are reduced from 10 to 6), avoiding the porosity fit permits us to apply the model to laboratory data that are not provided with accurate porosity measurements. For this second model, we also estimate the uncertainty of the fitting parameters and the elastic velocities. We tested this model on dry core measurements from literature and we verified that it returns elastic velocity trends as good as the original modified grain contact theory model with a reduced number of fitting parameters. Possible developments of the new model to add predictive power are also discussed.     

Isolation of plane shear wave using water saturated trench barrier

The article presents a theoretical investigation of the propagation of normally and obliquely incident plane shear waves past a rectangular trench filled up with water saturated soil sandwiched between anisotropic elastic media. The motivation for this work is due to the effectiveness of the isolation of waves by the in-filled rectangular trench. Unlike the most of the previous researchers, this model considers the soil–structure interaction effects and directly determines the influence of barrier in the form of in-filled trench on the mode of wave propagation. It is of interest to determine the reflection and transmission coefficients, and the energy partition distribution of shear waves in the in-filled rectangular trench showing the influence of barrier on the propagation of waves. An extensive parametric study through numerical computation is carried out to investigate the influence of the material properties of the in-filled trench and the amplitude ratios on shear waves. The in-filled trench barrier directly declines the intensity of waves significantly in such a way that the waves do not create any hazards to the nearby structures, if exists at all.     

Structure of residual oil as a function of wettability using pore-network modelling

In the water flooding of mixed-wet porous media, oil may drain down to relatively low residual oil saturations (S_or). Various studies have indicated that such low saturations can only be reached when oil layers in pore corners are included in the pore-scale modelling. These processes within a macroscopic porous medium can be modelled at the pore-scale by incorporating the fundamental physics of capillary dominated displacement within idealised pore network models. Recently, the authors have developed thermodynamic criteria for oil layer existence in pores with non-uniform wettability which takes as input geometrically and topologically representative networks, to calculate realistic S_or values for mixed-wet and oil-wet sandstones [16, 21]. This previous work is developed in this paper to include (i) the visualisation of the 3D structure of this residual oil, and (ii) a statistical analysis of this “residual/remaining” oil. Both the visualisation and the statistical analysis are done under a wide range of wettability conditions, which is reported for the first time in this paper.The structure of residual oil for strongly water wet systems is well known (where residual = remaining oil) and our model agrees with this but this structure changes radically for mixed wet systems (where residual ≠ remaining) and this has not yet been visualised experimentally. We find that for more water-wet systems high final residual oil saturations are reached at relatively small amounts of water injected and this oil is present in the pores as bulk oil. On the other hand, for more oil-wet systems we find a slow decrease of the amount of remaining oil with increasing amounts of injected water. During the process, the remaining connectivity of the oil phase is increasingly provided by oil layers only, hence the slow drainage. The final residual oil saturation, only reached in the theoretical limit of an infinite amount of injected water, is almost entirely contained in large number of (relatively low volume) oil layers, which are present in pores of most radius sizes.     

The effect of microporosity on transport properties in porous media

Sizeable amounts of connected microporosity with various origins can have a profound effect on important petrophysical properties of a porous medium such as (absolute/relative) permeability and capillary pressure relationships. We construct pore-throat networks that incorporate both intergranular porosity and microporosity. The latter originates from two separate mechanisms: partial dissolution of grains and pore fillings (e.g. clay). We then use the reconstructed network models to estimate the medium flow properties. In this work, we develop unique network construction algorithms and simulate capillary pressure–saturation and relative permeability–saturation curves for cases with inhomogeneous distributions of pores and micropores. Furthermore, we provide a modeling framework for variable amounts of cement and connectivity of the intergranular porosity and quantifying the conditions under which microporosity dominates transport properties. In the extreme case of a disconnected inter-granular network due to cementation a range of saturations within which neither fluid phase is capable of flowing emerges. To our knowledge, this is the first flexible pore scale model, from first principles, to successfully approach this behavior observed in tight reservoirs.     

Rayleigh波的频散方程高频近似分解和多模式激发数目

对于在Rayleigh波勘探中经常提取到的多模式频散曲线,如何在理论上对其特征进行描述是一个很有意义的问题. 本文根据Rayeigh波频散方程中传递矩阵具有的高频近似特征,提出了传递矩阵的高频近似分解公式,据此导出了频散方程的高频近似分解公式,并在此基础上,定义了Rayleigh波的4种基本模式(R模、S模、R型周期模和S型周期模). 根据频散方程具有的周期性特征,给出了周期模频散曲线之间的平均间隔近似公式,以及在任意频率和相速度段内频散方程根的数目（即多模式的激发数目）的预测公式. 在这些结果的基础上,本文最后提出了频散方程数值搜根的一种新方法.