A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.     

Direct forecasting of reservoir performance using production data without history matching

The conventional paradigm for predicting future reservoir performance from existing production data involves the construction of reservoir models that match the historical data through iterative history matching. This is generally an expensive and difficult task and often results in models that do not accurately assess the uncertainty of the forecast. We propose an alternative re-formulation of the problem, in which the role of the reservoir model is reconsidered. Instead of using the model to match the historical production, and then forecasting, the model is used in combination with Monte Carlo sampling to establish a statistical relationship between the historical and forecast variables. The estimated relationship is then used in conjunction with the actual production data to produce a statistical forecast. This allows quantifying posterior uncertainty on the forecast variable without explicit inversion or history matching. The main rationale behind this is that the reservoir model is highly complex and even so, still remains a simplified representation of the actual subsurface. As statistical relationships can generally only be constructed in low dimensions, compression and dimension reduction of the reservoir models themselves would result in further oversimplification. Conversely, production data and forecast variables are time series data, which are simpler and much more applicable for dimension reduction techniques. We present a dimension reduction approach based on functional data analysis (FDA), and mixed principal component analysis (mixed PCA), followed by canonical correlation analysis (CCA) to maximize the linear correlation between the forecast and production variables. Using these transformed variables, it is then possible to apply linear Gaussian regression and estimate the statistical relationship between the forecast and historical variables. This relationship is used in combination with the actual observed historical data to estimate the posterior distribution of the forecast variable. Sampling from this posterior and reconstructing the corresponding forecast time series, allows assessing uncertainty on the forecast. This workflow will be demonstrated on a case based on a Libyan reservoir and compared with traditional history matching.     

A multi-domain decomposition-based Fourier finite element method for the simulation of 3D marine CSEM measurements

We introduce a multi-domain decomposition Fourier finite element (MDDFFE) method for the simulation of three-dimensional (3D) marine controlled source electromagnetic measurement (CSEM). The method combines a 2D finite element (FE) method in two spatial dimensions with a hybrid discretization based on a Fourier FE method along the third dimension. The method employs a secondary field formulation rather than the total field formulation. We apply the MDDFFE method to several synthetic marine CSEM examples exhibiting bathymetry and/or multiple 3D subdomains. Numerical results show that the use of the MDDFFE method reduces the problem size by as much as 87 % in terms of the number of unknowns, without any sacrifice in accuracy.     

Dewatering Behaviour of a Uranium Ore Slurry Containing Clays

The main objective of this paper was to investigate the dewatering behaviour of a clayey uranium ore slurry. The slurry (containing 28% clay size) exhibited moderate water adsorption (w _l  = 83% and w _p  = 30%). Primarily composed of muscovite (46%) and quartz (30%), the clay minerals included illite (8%), chlorite (5%) and kaolinite (2%) alongside a CEC of 41 (cmol(+)/kg) with Ca²⁺ and Mg²⁺ as the dominant cations. Likewise, the high EC (17,600 μS/cm) and ionic strength (1.15 mol/L) indicated a flocculated microstructure due to the presence of SO₄ ^2? (22,600 mg/L) and Mg²⁺ (1340 mg/L) in the slurry water. Settling included sedimentation and consolidation at low initial solids condition (25–35%) whereas only consolidation was observed at high initial solids contents (40–50%). The average k reduced from 1.2 × 10^?6 m/s (initial s = 25%) to 5.3 × 10^?8 m/s (initial s = 50%) along with a void ratio reduction from 7.4 to 2.6. Due to thixotropic strength, volume compressibility during consolidation showed apparent pre-consolidation at low effective stress (0.3–2 kPa) with a reduction in void ratio from 2.6 to 2.5. The e _s was found to be 2.46 at σ′ = 2 kPa and was followed by a steeper slope with the void ratio reducing to 2.1 at σ′ = 31 kPa. Likewise, the hydraulic conductivity during consolidation decreased from 2.6 × 10^?9 m/s (at e = 2.6) to 2.0 × 10^?10 m/s (at e = 2.1).     

Investigation of Various Parameters Effect on Cerchar Abrasivity Index with PFC3D Modeling

The Cerchar test is one of the appropriate and routine tests for determining the rock abrasion; but as for the costs and pin wear measurement errors in laboratory procedures and lack of access to laboratory equipment, using of numerical modeling can lead to use of greater number of samples required during the course of mechanized excavation and reduce the costs and errors in the laboratory test. In this study, the Cerchar abrasivity test was modeled using PFC3D (Particle Flow Code in 3 Dimensions) software. In order to verify the simulation results, Cerchar laboratory test results obtained by Rostami (Rock Mech Rock Eng 47(5):1905–1919, 2014) were compared with the numerical modeling results. In modeling studies, the effects of some parameters such as apply load, test speed, pin hardness and scratching distance on pin wear were investigated. As a conclusion of the study, good agreement between modeling and experimental results was obtained for a given condition. As in the experiment with various loads in both laboratory tests and modeling, with increasing applied load the Cerchar abrasivity index also increased in the experiment with Rockwell hardness HRC (An abbreviation for Rockwell Hardness measured on the C scale. The Rockwell test determines the hardness by measuring the depth of penetration of an indenter under a large load compared to the penetration made by a preload, that on the C scale use from a indenter with 120° cone and 150 kgf load) 42.     

Seasonal residence and counterurbanization: the role of second homes in population redistribution in Finland

Since the late twentieth century, many developed countries have experienced population deconcentration, labelled as counterurbanization. There has been an academic discussion on the meaning, validity and universality of this concept, drivers of counterurbanization and its impacts on rural areas. To date, research on counterurbanization mostly apply static and discrete definitions of residence, migration and population, which is an increasingly simplistic view in the contemporary reality of a growing multitude of forms of mobility, often related to dual residence. Particularly large-scale quantitative studies on counterurbanization are confined by existing statistical practices. This paper attempts to overcome this obstacle and describe the transformation of the settlement system in Finland acknowledging the role of second home mobility. To achieve this goal, it introduces two alternative measures of population, seasonal and average population, and analyses their spatial dynamics between the years 1990 and 2010 based on georeferenced grid statistical data. The study finds that although registered population has been concentrating during the period in analysis, seasonal population has been increasingly dispersed due to the growing number of second homes. It shows that the counterurbanization process, though not noticed by conventional statistics, does occur in Finland, manifested by seasonal rather than permanent moves. The article concludes that various forms of mobility should be taken into account when analysing the urban–rural population dynamics and transformations of settlement systems as well as in rural development planning.     

岩石内应力的X射线-中子散射测量方法

岩石内应力的储存和释放对深地资源和能源开发具有重要约束。矿物晶格间距的动态变化是揭示岩石内应力演化的重要指示。X射线衍射法是当前获取矿物晶体结构的主要手段,可以准确测定矿物的特征衍射峰。但是,对于具有复杂结构和矿物组成的岩石,X射线在穿透过程中携带的统计信息,难以反映矿物晶格间距的细节特征。中子射线的能量高、穿透深度大,且与原子核反应,因而能更准确地刻画矿物元素位置和结构细节,但中子的波长较长、校准困难,结果具有不确定性,导致中子射线在岩石内应力研究中未能发挥应有作用。在分析X射线、中子技术优势和固有局限的基础上,提出了岩石内应力的X射线-中子衍射测量方法,重点剖析了该方法的基本原理,提出了该方法的关键技术与研究前景。得出了如下结论:利用X射线在确定岩石矿物晶格参数的优势,结合加热处理等技术,可实现岩石矿物无应力条件下晶格间距的标定; 采用中子衍射技术可精确测量岩石矿物的衍射偏移峰,实现岩石内应力深度轮廓的精准刻画; 结合X射线与中子的优势,可实现岩石内应力绝对值的精确测量。提出的技术手段有利于揭示岩石内应力的微观储存和释放过程,以及长期存留的物质条件和物理力学机制,并有望为岩石矿物组成和微观结构研究提供一种新的技术方案。     

Na depletion in modern adakites via melt / rock reaction within the sub-arc mantle

Interaction between slab-derived melt and mantle peridotite and the role of slab melt as a metasomatizing agent in the sub-arc mantle is being increasingly recognized. Adakite, the slab melt erupted on the surface, usually exhibits anomalously high MgO, CaO, Cr and Ni contents that indicate interaction with mantle peridotitite. Here we note that Cenozoic adakites have Na₂O contents below 5.8 wt.% with ∼95% samples lower than 5.0 wt.%, and are generally depleted in this component relative to experimental basalt partial melts (mostly beyond 5.0 wt.% and up to 9.0 wt.% Na₂O) produced under 1.5-3.0 GPa conditions that are most relevant to adakite production. We interpret the adakite Na depletion to be also a consequence of the melt / rock reaction that takes place within the hot mantle wedge. During ascent and reaction with mantle peridotite, primary adakite melts gain mantle components MgO, CaO, Cr and Ni but lose Na₂O, SiO₂ and perhaps K₂O to the mantle, leading to Na-rich mantle metasomatism. Selective assimilation of predominately mantle clinopyroxene, some spinel and minor olivine at high T/P has been considered to be an important process in producing high-Mg adakites from primary low-Mg slab melts [Killian, R., Stern, C. R., 2002. Constraints on the interaction between slab melts and the mantle wedge from adakitic glass in peridotite xenoliths. Eur. J. Mineral. 14, 25-36]. In such a process, Na depletion in the assimilated melt is the result of dilution due to the increase in melt mass. Phase relationships in the reaction system siliceous melt + peridotite and quantitative calculation suggest that assimilation of mantle clinopyroxene, olivine and spinel and fractional crystallization of sodic amphibole and orthopyroxene, under conditions of moderate T/P and increasing melt mass, is also an important process that modifies the composition of adakites and causes the Na depletion.     

交通荷载下加筋土桥台工作性能试验研究

加筋土桥台因其具有良好的复合体特性,被广泛应用于各项工程建设,研究它在实际工程条件下的工作性能对其设计和推广具有重要意义。加筋土桥台在服役过程中主要承受交通荷载作用,通过室内大比例缩尺模型的加速应力试验,研究交通荷载条件下加筋土桥台的工作性能,考虑加筋间距与筋材刚度的影响,综合分析了不同交通荷载等级和循环次数下桥台变形与筋材应变的变化规律。研究结果表明,加筋土桥台在交通荷载作用下整体变形随循环次数趋于收敛,其顶部沉降、面板水平位移与筋材应变均远小于规范给出的阈值;各变形值随加载时间呈阶梯状发展,不同条件下变形收敛趋势有所差异;增大筋材刚度和减小加筋层间距有利于控制桥台的变形,筋材刚度对变形的控制作用要优于加筋间距;交通荷载作用下面板水平应变与顶部竖向应变并不满足2倍关系,二者之间的定量关系有待进一步研究。     

陆源碎屑岩中的古土壤在地层分析中的应用

对各种古土壤的分类标准进行比较之后,重点介绍了Mack的古土壤分类及其土壤的演变序列。Mack的古土壤分类主要依据稳定矿物和保存在土壤中的古地貌属性,将古土壤划分为9种。在地层序列中古土壤的识别标志主要包括颜色、粒度、生物成因构造（植物生长痕迹）、原生沉积构造退化特征、层次化和边界特征、宏观结构特征、矿物及地球化学特征等。可以通过野外观察、室内显微镜观察、实验室化验及地震、测井资料来进行鉴别。陆相盆地具有很强的分隔性,盆地相变又非常复杂,寻找区域性较稳定的层序界面比较困难。越来越多的研究证实用古土壤划分陆相成因地层层序、解决地区及全球对比具有重要的实用价值。