地下非均匀非饱和带中地下洞室的渗流问题数值模拟--介质参数的灵敏度分析

影响任何一个地下废物处置场所长期行为的一个重要因素是进入地下废物贮置洞室内的渗流量的大小，预测地下洞室中的渗流量是困难的，特别是当地下废物贮置洞室位于非均质非饱和带中，三维数值模型用于研究地下非均质中非饱和流动及流入地下洞室中的渗流量，讨论了非均质模型与均质模型的比较以及数值剖分尺度对计算结果的影响，进入洞室中的渗漏率随着引入系统中入渗量的增加而增大，选用4个参数来衡量渗流场的非均质程度：(1)介质的平均渗透率K0；(2)Van Genuchten参数α；(3)渗透率空间分布相关尺度；(4)渗透率空间分布变化的标准差σ，根据一个随机实现的渗透率分布，通过改变平均渗透率来研究其对流入洞室中渗透量的影响，对一个固定的入渗率而言，流入洞室中的渗漏率将随着Van Genuchten参数α的减小而减小，模拟结果表明流入洞室的渗漏率与介质的平均渗透率相关，即随平均渗透率的增大，流入洞室的渗漏率亦增大，流入洞室的渗漏率还高度依赖于非均质渗透率场的空间分布相关长度与标准差，一个大的相关长度或高的标准差均能导致流入洞室的渗漏率增大。     

Generalized Coarse Graining Procedures for Flow in Porous Media

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.     

沉积盆地流体—岩石相互作用研究的现状

沉积盆地流体－岩石相互作用研究主要是通过实验地球化学和岩石学，地球化学模拟方法开展的，综述了这一领域的研究现状，主要包括有机酸来源和分布及期 对矿物稳定性的影响，地层水成因与演化、烃类与岩石间的氧化还原反应、以及储层润湿性的变化等，并展望了其 发展趋势。     

非均质土坝稳定性的渗流场和应力场耦合分析

綦于多孔介质的渗流特性和土的非线性本构关系，研究了渗流场与应力场的耦合效应，并针对一实际十坝建立了非线性耦合分析模型，探讨了位于软基之上的非均质土坝的渗流与应力耦合规律。分析结果表明，在渗流作用下，坝休和坝基中渗流与应力的相互作用是不应忽视的，其作用结果将会使土坝、特别足防渗心端附近应力水平和渗透坡降均有不同程度的提高，表明考虑渗流场和应力场耦合作用的计算结果对土坝稳定安全性不利。     

Scaling procedures for heterogeneous unconfined aquifers

Complex flows in heterogeneous confined and unconfined aquifers is a phenomenon that continues to present difficulties in flow mapping and modelling in the field, laboratory, and through numerical simulations. It is often the case with complicated phenomena that transformative scaling and reduction of the problem through symmetry is of great efficacy in the formation of predictive models in both the laboratory and computational settings. A detailed a study of the application of a broad class of Lie scaling transformations on a set of equations representing the groundwater flows in heterogeneous confined and unconfined aquifers has produced a set of scaling relationships between the spatial variables, hydrologic variables, and parameters. The set of scaling transformations preserve the structure of the equations in the sense that the scaling transformations leave the initial‐boundary value system representing the invariant groundwater flows. This theoretical approach elucidates not only the scaling relationships but also the properties that hydrologic variables and parameters must satisfy in order for calling to be possible. Validation of the theory developed is carried out through a series of four numerical simulations using the USGS modflow ‐2005 software package. The results of these experiments demonstrate that the derived scaling transformations can effectively form predictive models of large‐scale phenomena at small scales with negligible error in many cases. Comments on the limitations of the approach and directions for future research are made in the closing sections. Copyright © 2016 John Wiley & Sons, Ltd.     

Application of topo-edaphic factors and remotely sensed vegetation indices to enhance biomass estimation in a heterogeneous landscape in the Eastern Arc Mountains of Tanzania

Estimating tropical biomass is critical for establishment of conservation inventories and landscape monitoring. However, monitoring biomass in a complex and dynamic environment using traditional methods is challenging. Recently, biomass estimates based on remotely sensed data and ecological variables have shown great potential. The present study explored the utility of remotely sensed data and topo-edaphic factors to improve biomass estimation in the Eastern Arc Mountains of Tanzania. Twenty-nine vegetation indices were calculated from RapidEye data, while topo-edaphic factors were taken from field measurements. Results showed that using topo-edaphic variables or vegetation indices, biomass could be predicted with an R² of 0.4. A combination of topo-edaphic variables and vegetation indices improved the prediction accuracy to an R² of 0.6. Results further showed a decrease in biomass estimates from 1162 ton ha^?1 in 1980 to 285.38 ton ha^?1 in 2012. This study demonstrates the value of combining remotely sensed data with topo-edaphic variables in biomass estimation.     

Geochemical constraints on the distribution of trace elements and volatiles in fluorapatite from the Panasqueira tin-tungsten deposit (Portugal)

A large, euhedral crystal of fluorapatite (ca. 19.5 × 20.0 mm) from the Panasqueira tin-tungsten deposit (Portugal) was investigated in terms of the distribution of trace elements by using several microanalytical techniques. The studied material represents almost pure fluorapatite with minor amounts of other cations (mainly Sr, Mn, REE and Fe), OH and Cl. Particular interest was given to the distribution of rare earth elements with respect to the crystallographic orientation. A broad range of analytical techniques were used, including optical microscopy coupled with cathodoluminescence imaging, electron probe microanalysis (EPMA), laser ablation – inductively coupled plasma mass spectrometry (LA-ICPMS), Raman microspectroscopy, and simultaneous thermal analysis coupled with mass spectrometry. The investigated crystal consists of the main crystal with a distinct core and rim (Ap_2core and Ap_2rim, respectively), which grew on a previous, euhedral crystal (Ap₁). The fluorapatite demonstrates various types of zoning: regular oscillatory, irregular, and internal sectoring, which is also reflected in trace elements concentrations. The rim Ap_2rim has lower concentrations of Mn, Sr and Fe, and significantly higher concentrations of REE compared to the core Ap_2core and older crystal Ap₁. Furthermore, the rim Ap_2rim is strongly depleted in Th, U and Pb. The entire crystal shows elevated Eu contents, expressed as a strong positive anomaly in chondrite-normalized REE patterns. With regards to the volatiles, F concentrations are constant in Ap₁, Ap_2core and Ap_2rim, whereas Cl is below the EPMA detection limit. The Ap_2rim was the only part of the investigated material containing OH and CO₃, which were observed in the Raman spectra. Furthermore, part of the crystal Ap_2core is extensively altered, likely due to fluid-induced metasomatic processes. LA-ICPMS U-Pb dating yielded highly discordant dates due to common Pb content. A lower intercept age of 297 ± 13 Ma (MSWD = 0.13) indicates the age of the fluorapatite crystallization. The overall analytical data constrain growth and post-growth processes, including crystallization of Ap₁ and Ap_2core, which both have typical hydrothermal Sn-W deposit characteristics, whereas Ap_2rim is related to a carbonate stage of the mineralization in the Panasqueira deposit.     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。     

A new load balancing scheme on heterogeneous database cluster

IntroductionAlthoughthe development of cluster has gainedmore and more advantages of low-cost ,high-per-formance and high-availability computing, effi-cient resource usage such as CPU, memory andI/Ois still a keyto achieve better performanceincluster syst…