Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

克孜尔水库大坝的面波检测

1999年 3月 15日克孜尔地震之后对克孜尔水库坝体出现的裂缝进行了面波勘探。首先介绍了面波勘探法的基本原理 ,然后通过典型剖面检测与资料处理找出异常区、带并进行解释。试验表明面波勘探法适用于克孜尔水库大坝监测     

基于DEM的分布式水文模型构建方法

基于 DEM的分布式水文模型是现代水文学同高科技 (如计算机技术、 3S技术等 )相结合的产物 ,是研究变化环境下水文循环与水资源演化规律的理想工具 ,代表了水文模型的最新发展方向。从 DEM的特性出发 ,本文探讨并总结了分布式水文模型的特点、建模思路和模型基本结构框图。在流域离散化方面 ,重点介绍了分布式水文模型常用的三种单元划分方法 ;最后 ,针对分布式水文模型构建问题 ,从“输入模块”、“单元水文模型”、“河网汇流模型”三方面 ,阐述了分布式水文模型微结构构建方法。     

不同损失条件下的泥石流预报模型

泥石流预报是泥石流减灾的重要手段之一，然而泥石流形成的复杂性使泥石流预报准确度低，误报和漏报率较高。泥石流误报和漏报都会造成损失，但二者造成的损失有很大的差别。为了减少泥石流误报或漏报造成的损失，应当考虑两种错报造成损失的不同。根据使总平均损失达到最小的原则，建立了不同损失条件下的泥石流预报模型，并将该模型应用到云南东川蒋家沟。     

泥石流多发干旱河谷区植被恢复研究

云南东北部东川小江流域的干旱河谷土地退化十分严重，土地表层的砾石含量已超过60％，有机质含量、全氮、有效氮、有效磷、有效钾的含量相对于残存燥红土明显下降。草坡是当地的主要植被类型，并未起到防治土地退化的应有作用。以合欢（Leucaena glauca)、马桑（Coriaria sinica)种植为主的恢复方式对土地的恢复有明显的改善作用。高密度的乔木树种合欢的种植使林下的草本物种的数量和组成发生了极大变化，较耐荫的植物占据了林分下层，且物种数量较草坡减少，群落多样性指数降低。三种马桑种植方式下，不但原有的草坡群落的物种仍有大量的生存空间，而且为一些新物种创造了生存环境，从而使物种的种类较原有草坡地增加，提高了群落的多样性指数。合欢种植的生物量最高，对解决当地薪柴严重短缺意义重大，马桑种植在增加灌木层生物量的同时，并未改变草坡的物种组成，还可兼顾解决当地的饲草和放牧，应大力推广。     

A New Model for Heat Flow in Extensional Basins: Estimating Radiogenic Heat Production

Radiogenic heat production (RHP) represents a significant fraction of surface heat flow, both on cratons and in sedimentary basins. RHP within continental crust—especially the upper crust—is high. RHP at any depth within the crust can be estimated as a function of crustal age. Mantle RHP, in contrast, is always low, contributing at most 1 to 2 mW/m² to total heat flow. Radiogenic heat from any noncrystalline basement that may be present also contributes to total heat flow. RHP from metamorphic rocks is similar to or slightly lower than that from their precursor sedimentary rocks. When extension of the lithosphere occurs—as for example during rifting—the radiogenic contribution of each layer of the lithosphere and noncrystalline basement diminishes in direct proportion to the degree of extension of that layer. Lithospheric RHP today is somewhat less than in the distant past, as a result of radioactive decay. In modeling, RHP can be varied through time by considering the half lives of uranium, thorium, and potassium, and the proportional contribution of each of those elements to total RHP from basement. RHP from sedimentary rocks ranges from low for most evaporites to high for some shales, especially those rich in organic matter. The contribution to total heat flow of radiogenic heat from sediments depends strongly on total sediment thickness, and thus differs through time as subsidence and basin filling occur. RHP can be high for thick clastic sections. RHP in sediments can be calculated using ordinary or spectral gamma-ray logs, or it can be estimated from the lithology.     

粘性泥石流残留层的形成及减阻作用--以云南蒋家沟泥石流为例

通过对云南蒋家沟粘性泥石流的野外观测 ,讨论了粘性泥石流残留层的定义特征 ,在大量实测数据基础上 ,分析了残留层易形成时段的特性 ;分别对一次泥石流过程、若干阵泥石流和一阵泥石流残留层的变化进行了相关分析 ;最后用观测数据分析了泥石流残留层的减阻作用     

三峡库岸城区滑坡治理与沿江公路建设一体化模式--以重庆万州清泉路滑坡为例

三峡库岸城区滑坡90%属于由第四纪松散土体组成的散体滑坡,滑坡区属于城区重要组成部分,降雨及三峡库水位降落是此类滑坡孕发的动力因素。本文以重庆万州清泉路滑坡为例,基于库岸城区滑坡防治原则,构建了实用的滑坡治理与沿江路建设一体化模式,推导了抗滑桩的内力计算方法;详细分析了一体化模式的结构组成及受力特性。研究成果对于三峡库区城区滑坡的治理与开发利用具有较强的指导借鉴作用。