Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Weakness and mechanical anisotropy of phyllosilicate-rich cataclasites developed after mylonites of a low-angle normal fault (Simplon Line,Western Alps)

The Simplon Fault Zone is a late-collisional low-angle normal fault (LANF) of the Western Alps. The hanging wall shows evidence of brittle deformation only, while the footwall is characterized by a c. 1 km-thick shear zone (the Simplon Fault Zone), which continuously evolved, during exhumation and cooling, from amphibolite facies conditions to brittle-cataclastic deformations. Due to progressive localization of the active section of the shear zone, the thermal-rheological evolution of the footwall resulted in a layered structure, with higher temperature mylonites preserved at the periphery of the shear zone, and cataclasites occurring at the core (indicated as the Simplon Line). In order to investigate the weakness of the Simplon Line, we studied the evolution of brittle/cataclastic fault rocks, from nucleation to the most mature ones. Cataclasites are superposed on greenschist facies mylonites, and their nucleation can be studied at the periphery of the brittle fault zone. This is characterized by fractures, micro-faults and foliated ultracataclasite seams that develop along the mylonitic SCC′ fabric, exploiting the weak phases mainly represented by muscovite and chlorite. Approaching the fault core, both the thickness and frequency of cataclasite horizons increase, and, as their thickness increases, they become less and less foliated. The fault core itself is represented by a thicker non-foliated cataclasite horizon. No Andersonian faults or fractures can be found in the footwall damage zone and core zone, whilst they are present in the hanging wall and in the footwall further from the fault. Applying a stress model based on slip tendency, we have been able to calculate that the friction coefficient of the Simplon Line cataclasites was <0.25, hence this fault zone is absolutely weak. In contrast with other fault zones, the weakening effect of fluids was of secondary importance, since they accessed the fault zone only after an interconnected fracture network developed exploiting the cataclasite network.     

Effects of surface flux parameterization on the numerically simulated intensity and structure of Typhoon Morakot (2009)

The effects of surface flux parameterizations on tropical cyclone(TC) intensity and structure are investigated using the Advanced Research Weather Research and Forecasting(WRF-ARW) modeling system with high-resolution simulations of Typhoon Morakot(2009).Numerical experiments are designed to simulate Typhoon Morakot(2009) with different formulations of surface exchange coefficients for enthalpy(C_K) and momentum(C_D) transfers,including those from recent observational studies based on in situ aircraft data collected in Atlantic hurricanes.The results show that the simulated intensity and structure are sensitive to C_K and C_D,but the simulated track is not.Consistent with previous studies,the simulated storm intensity is found to be more sensitive to the ratio of C_K/C_D than to C_K or C_D alone.The pressure-wind relationship is also found to be influenced by the exchange coefficients,consistent with recent numerical studies.This paper emphasizes the importance of C_D and C_K on TC structure simulations.The results suggest that C_D and C_K have a large impact on surface wind and flux distributions,boundary layer heights,the warm core,and precipitation.Compared to available observations,the experiment with observed C_D and C_K generally simulated better intensity and structure than the other experiments,especially over the ocean.The reasons for the structural differences among the experiments with different C_D and C_K setups are discussed in the context of TC dynamics and thermodynamics.     

天津市地热资源开发利用特征及存在问题分析

天津地区的地热地质条件良好,蕴藏着丰富的地热资源。在地热资源的开发利用过程中,每年可产生可观的经济、社会和环境效益,但由于开发方式的不尽合理,带来诸多问题。文章阐述了天津开发利用地热资源的特征,分析了存在的问题。认为不同热储层间水力联系的存在和地热井间距的过小,造成开采层间相互影响较大;地热井布局不合理、开采层位相对集中、回灌程度不高导致局部热储层压力大幅度下降,持续的大幅度压力下降将引发吊泵和地面沉降;超标准排放地热尾水会造成资源浪费和热污染;采、灌井间距的过小和相对位置的不合理,存在冷突破风险。针对不同问题,给出了针对性解决问题的建议。     

河岸带潜流层温度示踪流速计算方法

为揭示低温水影响下的河岸带潜流层的温度场和流场分布特性,利用野外水温水位实时监测试验,研究河岸带潜流层温度场在不同季节、不同空间位置上的分布特性,并利用水温资料计算获得地下水流速。结果表明:河岸带潜流层温度场在夏季和冬季分别呈现出"上暖下冷"和 "上冷下暖"的温度分层现象;通过对温度示踪方法的4种计算方法进行分析比较,得到Hatch相位法计算的地下水流速具有较高的准确性,在2014年12月15—31日时段内流速大小为1.03×10-4~7.96×10-4m/s,在空间上,断面深度增加,地下水流速降低,且不同深度流速曲线接近平行。     

从藏南陆-陆碰撞带深部结构构造演化探讨斑岩铜矿的成岩成矿问题

本文以INDEPTH项目对印度大陆与欧亚大陆碰撞带深部成像结果为基础,从构造演化角度探讨藏南陆-陆碰撞带冈底斯斑岩铜矿带的成矿作用问题。深部探测给出的碰撞带深部结构与侯增谦等地质学家提出的深部结构有较大的异同,如何协调起来以深化对藏南陆-陆碰撞条件下成矿作用的认识,这是本文讨论的中心。藏南碰撞带成矿实际上是在新特提斯大洋岩石圈俯冲形成的冈底斯岩浆弧成矿作用的基础上,再经过陆-陆碰撞挤压强烈改造后的再成矿。碰撞带的深部结构构造演化的特点是:(1)新特提斯大洋岩石圈板块向北连续俯冲了约120 Ma,形成的冈底斯陆缘火山岩浆弧带,这导致了陆缘带地壳增厚并含有大量的地幔岩浆流体物质(如南美安第斯成矿带那样);(2)在印度大陆与冈底斯陆缘弧接近碰撞时,在对挤中新特提斯大洋洋壳与大洋岩石圈地幔发生向上挤出与向下拆沉,并使部分洋壳残片和大洋岩石圈物质保存在中上地壳内;(3)两大陆岩石圈碰撞对接后,印度岩石圈地幔加深达70~80 km并沿地壳底部向北推进,并将加厚地壳内大量的成矿物质、钙碱性岩浆,洋壳及新生的下地壳,以及部分地幔物质从地壳底部将其围限起来,成为后期再成矿的物质基础;(4)查明了碰撞带深部壳/幔间产生了一层中间速度层(相当于MASH层),在中上地壳部位出现一层巨大的部分熔融层;(5)在碰撞挤压下冈底斯带内产生多组断裂构造,大型逆冲断裂系与背冲断裂,并引发了含矿岩浆的再活动,并在浮力(下地壳内)和挤压力作用下多次活动上升生成斑岩型铜矿床;(6)成矿后地表遭受过强烈的风化剥蚀作用,使矿床出露地表。     

日照市河谷平原区成井方法对比研究

在分析日照河谷平原区水文地质条件的基础上,通过大量抽水试验数据,对比分析研究了当地民井与本次施工水井成井方法的差异,指出了当地民井在成井方法上的不足和缺陷,并提出了适宜该区的成井方法建议。     

城市生活垃圾卫生填埋场沉降变形分析

城市生活垃圾填埋场地基沉降和垃圾堆体沉降受到的影响因素多,现行规范的计算方法直接采用与实际相差较大,为减少不均匀沉降对垃圾防渗层的影响、准确确定垃圾储量和保证垃圾堆体的稳定,开展垃圾填埋场沉降研究具有重要意义。垃圾填埋场沉降变形包括地基沉降和垃圾堆体沉降,垃圾堆体沉降又分为主沉降和次沉降,主沉降完成时间短于次沉降。通过对陕西省宝鸡市垃圾填埋场地基土进行沉降计算,最大沉降差超过100 mm,对垃圾防渗层不利影响较大;而对垃圾堆体沉降计算,修正后的计算值与沉降观测值非常接近。考虑不同地区地层差异和垃圾成分差异,积累沉降观测资料具有重要意义。     

Analysis of interaction between ship bottom air cavity and boundary layer

The successful designs of hulls for ships employing drag reduction by air bottom cavitation have been based on solutions of inverse problems of the theory of ideal incompressible fluid. However, prediction of the drag reduction ratio, the air demand by ventilated cavities and the cavity impact on the hull–propeller interaction is impossible in the framework of this theory because all mentioned characteristics depend on interaction of air cavities with the ship boundary layers. Because the known CFD tools are not fitted to ventilated cavitation at low Froude numbers, an analysis of this interaction requires a novel flow model. This model includes the incompressible air flow in the ventilated cavity, the compressible flow of a water–air mixture in the boundary layer on cavities and downstream of them and the curl-free incompressible outer water flow. The provided 2D computations employing this model allows for explanations of the earlier observed effects and for prediction of the air demand by ventilated cavities. The computed velocity profiles downstream of cavities are in the accordance with the available experimental data.