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.     

Wrench tectonics of Abu Gharadig Basin,Western Desert,Egypt: a structural analysis for hydrocarbon prospects

Mapping based on the interpreted seismic data covering the Abu Gharadig Basin in the northern Western Desert has revealed that the deposition of the Upper Cretaceous succession was controlled by dextral wrench tectonics. This dextral shear accompanied NW movement of the African Plate relative to Laurasian Plate. Structural depth maps of the Cenomanian Bahariya Formation and the Turonian-Coniacian D and A members of Abu Roash Formation display a clear NE-SW anticline dissected by NW-SE normal faults. This anticline represents one of the en echelon folds characterizing the wrench compressional component. The interpreted normal faults reflect the extensional T-fractures associated with the wrenching tectonics. The interaction between the aforementioned NE-SW anticline with the NW-SE extensional faults further confirms the effect of the Upper Cretaceous dextral wrench tectonic. However, the influence of this wrench tectonics was gradually diminishing from the Cenomanian up to the Coniacian times. The NW-SE compressional stress of the dextral wrench compressional component during the Cenomanian up to Coniacian age was greater in NW direction than the SE direction. Three mapped structural closures which are predicted to be potential hydrocarbon traps belonging to the Bahariya Formation and Abu Roash D Member, and are recommended to be drilled in the study area, with potential reservoirs. The regularity of the en echelon array of both anticlines and normal faults within the wrench zones suggests additional closures may be located elsewhere beside the study area.     

A simulation study of capillary transport,preferential retention and distribution of salts in historic sandstone buildings

Salt crystallisation is a major problem of deterioration in historic stone buildings, monuments and sculptures. The capillary rise of soil water is one of the primary sources of salts in stone structures, which evaporates leaving the salts behind. It has been noted that the spatial distribution profile of different species of salts crystallised in historic stone buildings is not homogeneous, i.e. different salts crystallise at different locations. The capillary transport and inhomogeneous spatial distribution of different salts in the porous building materials has been considered to be a result of solubility-dependent crystallisation; however, the factors responsible for this phenomenon are not clearly known. This paper aims to investigate the factors influencing the differential distribution of salts during capillary rise of soil water. In this study, the capillary transport of salts was simulated on two different sandstones—Locharbriggs, a Permo–Triassic, red sandstone and Stoke Hall, a Carboniferous, buff sandstone. The experiments were carried out under controlled environmental conditions to eliminate the possibility of evaporation-driven crystallisation of salts depending on their solubilities. The results indicate that fractionation or differential distribution of salts takes place even in the absence of evaporation and crystallisation. The sandstones exhibit properties like an ion exchange column, and ionic species present in the salt solution show differential distribution within the porous network of sandstone.     

广西田阳县那赖旧石器遗址发现的石制品和玻璃陨石

对广西田阳县那赖旧石器遗址进行实地调查,在百色盆地右江第Ⅳ级阶地的网纹红土中发现了一件玻璃陨石和大量的旧石器（包括12件手斧）。玻璃陨石的特征表明其为原地埋藏,未经后期流水作用搬运。玻璃陨石和石器之间的关系进一步证实,百色盆地第Ⅳ级阶地中保存的玻璃陨石可以代表手斧的埋藏年代。那赖遗址的石制品组合与百色盆地其他遗址的类似。通过与世界上其他阿舍利手斧在形态上和技艺上对比表明,该遗址手斧代表了早中更新世之交东亚的阿舍利技术特点。

     

Markov Chain Random Fields for Estimation of Categorical Variables

Multi-dimensional Markov chain conditional simulation (or interpolation) models have potential for predicting and simulating categorical variables more accurately from sample data because they can incorporate interclass relationships. This paper introduces a Markov chain random field (MCRF) theory for building one to multi-dimensional Markov chain models for conditional simulation (or interpolation). A MCRF is defined as a single spatial Markov chain that moves (or jumps) in a space, with its conditional probability distribution at each location entirely depending on its nearest known neighbors in different directions. A general solution for conditional probability distribution of a random variable in a MCRF is derived explicitly based on the Bayes’ theorem and conditional independence assumption. One to multi-dimensional Markov chain models for prediction and conditional simulation of categorical variables can be drawn from the general solution and MCRF-based multi-dimensional Markov chain models are nonlinear.     

一种改进的基于非高斯性最大化的预测反褶积算法

The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.     

Turbulent kinetic energy budget in a gravel-bed channel flow

The present experimental investigation focuses on the characteristics of near bed turbulence in a fully rough, uniform open-channel flow over a gravel-type bed. Due to bed topography small scale heterogeneity, the flow is not uniform locally in the near bed region and a double averaging methodology is applied over a length scale much larger than the gravel size. The double-averaged Turbulent Kinetic Energy (TKE) budget derived in the context of the present flow over a gravel bed differs from the TKE budget written for flow over a vegetation canopy. The non-constant shape of the roughness function measured in our gravel bed leads to an additional bed-induced production term which is null for vertical roughness elements, such as simplified vegetation elements. The experimental estimation of the terms of the TKE budget reveals that the maximum turbulent activity takes place away from the reference plane, near the roughness crests. However, within the interface sublayer the work of the bed induced velocity fluctuations against the Reynolds stress is of the same magnitude as the main turbulence production term. Consequently, the characteristics of the TKE budget have similarities with uniform flows over canopies and strongly differ from uniform flows over smooth and transitionally rough flows over sedimentlike beds.     

Validity of the Bouchet’s complementary relationship at 102 observatories across China

Bouchet in 1963 hypothesized that for large homogeneous land surface with minimum advection of heat and moisture, there exists a 1:1 complementary relationship of potential and actual evaporation coupled through land-atmosphere feedbacks. The complementary relationship has been widely used to estimate regional actual evaporation and explain the pan evaporation paradox. We examine the standardized potential evaporation (potential evaporation divided by wet environment evaporation) at 102 observatories at different elevations across China. Generally, the relationship is appropriate at the low elevations (<1000 m). With the increase of elevation, vapor transfer power becomes much less than radiation energy budget because of lower vapor pressure deficit and stronger global solar radiation. As a result, at the high elevations (over 1000 m), the excess energy resulted by limited moisture availability is not enough to be converted into drying power of the air. This result suggests that the complementary relationship is asymmetric at the high elevations. Supported by the Presidential Special Award Foundation, the Chinese Academy of Sciences (Grant No. O7R70020SD) and the National Key Technology R & D Program (Grant No. 2006BAC08B0408)     

地震区多层石结构抗震抗剪能力的模糊可靠性分析

本文通过考虑地震烈度的随机性和模糊性,以及石结构强度破坏等级界限的模糊性,研究多层石结构抗震抗剪能力的可靠性问题.将地震地面运动模拟为含有模糊烈度参数的平稳过滤有色噪声的随机过程,引入抗剪强度破坏指数,建立了石结构抗震抗剪能力的模糊安全准则,为研究多层石结构的抗震抗剪能力可靠性提供了一个合理的评定方法。     

大跨度桥梁强震动数字监测系统的研制

主跨1 196 m的龙江特大桥是我国第一座位于高烈度地震区的大跨度悬索桥,频繁的地震活动对大桥结构安全构成严重威胁。大桥建设时建立了健康监测系统,其中包含多个结构地震反应观测点,随后布设了场地强震动观测系统,这是对大桥地震响应观测的进一步补充和完善。本文详细介绍场地强震动观测系统的设计方案、系统构成、通信方式等。该系统获取的观测结果,可为龙江大桥在地震作用下结构健康诊断提供依据,为大桥结构地震响应分析提供可靠的地震动输入,有助于了解大桥结构在地震动作用下的反应特征, 提高大跨桥梁的抗震设计水平。