水质大涡模拟数学模型研究

建立了基于弱可压缩流基础上的水质大涡模拟数学模型及污染物质扩散系数与平均运动的关系。数值计算采用有限体积法和预测校正法,固壁边界条件采用"部分滑移条件"。并应用该模型对重庆嘉陵江、长江两江交汇段的流场,污染物质浓度分布进行了计算,计算结果与实测值比较吻合。研究表明,横向扩散与流动状态有关,当支流流量相对较小时,横向扩散较弱,形成岸边污染带,当支流流量相对较大时,横向扩散较强,污染带向江心移动。     

基于GIS的马尾马江片区公共厕所布局分析与优化

公共服务配套设施的合理布局和质量提升是新型城镇化的内涵之一,将GIS技术应用于规划领域可使规划成果具更理性与可实施性。以福州市马尾区马江片区为例,基于出行人群密度、区域繁华程度、公共厕所服务半径三者间的对应关系,对马江片区公共厕所的空间布局进行分析,认为该片区现状公共厕所存在数量与质量上的双重不足。针对存在的不足,提出优化算法与优化布局方案,并基于GIS技术实现可视化模拟。优化后的公共厕所共有53座,较现状的27座增加26座,符合《城市公共厕所设计标准》规定;公共厕所服务区覆盖片区总面积的91.5%,较现状的38.8%增加52.7%,满足城市规划与未来发展的需求。     

基于构造导向的叠后处理技术——以渤海湾盆地A油田为例

对地震资料进行叠后解释性处理时,往往只是单纯的数学运算很少受到地质条件的约束。当地下岩层为水平地层时,这种方法会获得较好的效果。但地下岩层具有一定的倾角,这时就可能会发生分辨率降低、信噪比下降或资料失真等现象。针对这个问题,引入构造导向体来约束叠后资料的处理过程。首先根据原始地震资料生成包含倾角和方位角信息的构造导向体;其次在构造导向体的控制下以相似性为判断条件对原始地震资料进行构造导向的中值滤波或扩散滤波;最后通过原始地震资料和相干体等对处理前、后的资料品质进行了综合的分析对比。方向信息的引入使运算沿着构造的走向进行,消除了地层倾斜的影响,最大程度保证资料真实性的同时使资料的品质得到极大地提高,有效指导了A油田井位的设计和实施。     

江西德兴斑岩铜矿床某些元素的扩散和渗滤作用

以江西德兴斑岩铜矿床之细矿脉为研究对象,文章介绍了细矿脉剖面采集及研究的工作方法,着重讨论了自矿脉至围岩元素的空间展布特征。成矿元素Cu、Ag、S和控矿元素Fe、Rb、K都属确定性分布型。对照元素的扩散和渗滤理想模式,得到如下认识:在成矿成晕的地球化学过程中,成矿元素Cu、Ag、S以及Fe以扩散作用为主,控矿元素K和Rb以渗滤作用为主。     

最优化的矿物扩散地质速率计:北大别黄土岭麻粒岩缓慢冷却的证据

对北大别黄土岭麻粒岩中石榴子石进行了详细的电子探针工作，并根据一颗石榴子石的成分环带，建立了扩散模型并进行了数值模拟，结合最优化的方法，得到了极其缓慢的冷却速率0.1℃/a。从扩散动力学数值模拟的角度对大别山造山带的冷却速率进行了探索性研究。该结果表明，黄土岭麻粒岩在受麻粒岩相改造后经历的是一个极其缓慢、持续时间很长的冷却过程。     

软煤和硬煤的甲烷吸附扩散特性对比

研究煤的吸附和放散特性,对于查明煤与瓦斯突出中甲烷的作用,提高煤层气采收率或矿井抽采率具有重要作用。通过自行研制的实验装置,对河南新安矿和平煤一矿软煤和硬煤甲烷吸附和扩散的性质、相同吸附压力下的吸附量和Langmuir吸附常数的差异进行分析,阐明煤表面能的控制因素,查明煤在初始解吸时刻的扩散系数。结果表明,同一煤层中软煤比共生的硬煤具有更大的Langmuir吸附体积,在同一吸附压力下同阶煤的吸附量与煤的软硬和吸附压力有关,软煤和硬煤的表面能与煤阶有关,短时间解吸时软煤扩散系数大于硬煤,长时间解吸后硬煤的扩散系数大于软煤,且软硬煤的扩散系数均随解吸时间延长和吸附平衡压力的增加而减小。     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Numerical estimation of effective diffusion coefficients for charged porous materials based on micro-scale analyses

In order to describe diffusive transport of solutes through a porous material, estimation of effective diffusion coefficients is required. It has been shown theoretically that in the case of uncharged porous materials the effective diffusion coefficient of solutes is a function of the pore morphology of the material and can be described by the tortuosity (tensor) (Bear, 1988 [1]). Given detailed information of the pore geometry at the micro-scale the tortuosity of different materials can be accurately estimated using homogenization procedures. However, many engineering materials (e.g., clays and shales) are characterized by electrical surface charges on particles of the porous material which strongly affect the (diffusive) transport properties of ions. For these type of materials, estimation of effective diffusion coefficients have been mostly based on phenomenological equations with no link to underlying micro-scale properties of these charged materials although a few recent studies have used alternative methods to obtain the diffusion parameters (Jougnot et al., 2009; Pivonka et al., 2009; Revil and Linde, 2006 2, 3 and 4). In this paper we employ a recently proposed up-scaled Poisson–Nernst–Planck type of equation (PNP) and its micro-scale counterpart to estimate effective ion diffusion coefficients in thin charged membranes. We investigate a variety of different pore geometries together with different surface charges on particles. Here, we show that independent of the charges on particles, a (generalized) tortuosity factor can be identified as function of the pore morphology only using the new PNP model. On the other hand, all electro-static interactions of ions and charges on particles can consistently be captured by the ratio of average concentration to effective intrinsic concentration in the macroscopic PNP equations. Using this formulation allows to consistently take into account electrochemical interactions of ions and charges on particles and so excludes any ambiguity generally encountered in phenomenological equations.     

Rapid cooling and exhumation of eclogitic rocks from the Great Caucasus, Russia

Diffusion modelling of growth-zoned garnet is used in combination with standard geothermometric and geobarometric techniques to estimate cooling and denudation rates from the mafic eclogites of the Red Cliff area, Great Caucasus, Russia. Euhedral garnet porphyroblasts exhibit different degrees of prograde growth zoning depending on the size of the grain (100 μm to several mm in diameter). Zoning patterns are mainly expressed in terms of Fe–Mg exchange, with 100*Mg/(Mg+Fe) increasing from 18–20 to 33–37 from core to rim. Geothermobarometry yields conditions of 680±40 °C and a minimum of 1.6±0.2 GPa and of 660±40 °C and 0.8±0.2 GPa for the high-pressure and retrograde stages of equilibration, respectively. A temperature of 600±40 °C has been recorded for the late-stage metamorphic overprint in the mica schists surrounding the eclogites. Relaxation of garnet zoning profiles was modelled for three different hypothetical P–T –t trajectories, all with an initial temperature of 680 °C and a pressure change of 0.8 GPa. The first two trajectories involve decompression associated with regular cooling down to 660 °C (near isothermal) and 600 °C. The third path is a two-step trajectory comprising near-isobaric cooling down to 600 °C followed by isothermal decompression to 0.8 GPa. These P–T trajectories cover as wide a range of pressure and temperature changes endured by the rocks as possible, thus representing extreme cases for calculating cooling and exhumation rates. Calculations indicate that the zoning pattern of the smallest garnet (i.e. garnet for which the zoning is most easily eliminated during post-growth processes) along the different paths can be preserved for the following average exhumation and cooling rates: path 1, 143 mm a^?1 and 102 °C Ma^?1; path 2, 60 mm a^?1 and 171 °C Ma^?1; path 3, 11–30 mm a^?1 and 200–400 °C Ma^?1. These results are discussed in light of theoretical P–T–t paths extracted from thermal models of regions of thickened crust, and from analogue models of accretionary wedge and continental lithosphere subduction.     

Multiscale anisotropic diffusion for ringing artifact suppression in geophysical deconvolution data

Ringing artifact degradations always appear in the deconvolution of geophysical data. To address this problem, we propose a postprocessing approach to suppress ringing artifacts that uses a novel anisotropic diffusion based on a stationary wavelet transform (SWT) algorithm. In this paper, we discuss the ringing artifact suppression problem and analyze the characteristics of the deconvolution ringing artifact. The deconvolution data containing ringing artifacts are decomposed into different SWT subbands for analysis, and a new multiscale adaptive anisotropic filter is developed to suppress these degradations. Finally, we demonstrate the performance of the proposed method and describe the experiments in detail.