A HORIZONTAL 2-D HYDRAULIC NUMERICAL MODEL AND IT’S APPLICATIONS TO FLOOD FORECAST

1 INTRODUCTION There were many flood disasters in China in recent years. When the water level in a river is very high, weak parts of its dike may be destroyed resulted in the submersion of the protected land and properties. It is of significance for decision-makers to exactly predict the processes of flood propagation during flood control. There are many modes of dike bursting, such as seepage destroying by overflow on top of dike caused by dike body sinking induced by piping and soil fl…     

华北地台北缘地球物理场特征与金属矿床预测

对华北地台北缘的地球物理场特征进行了探讨,根据重、磁异常数据反演计算了该地区的莫霍界面、居里界面、磁性界面的起伏.利用地球物理场资料和反演计算结果对该地区的构造格架和断裂进行了推断,同时预测了9个成矿远景区:(1)集宁一呼和浩特金成矿区;(2)张家口一赤城金、银多金属成矿区;(3)密云一高岭金、铜成矿区;(4)青龙一马兰庄金成矿区;(5)秦皇岛、金多金属成矿区;(6)郝家营多金属成矿区;(7)承德地区金、多金属成矿区;(8)赤峰一喀喇沁旗金成矿区;(9)宁城东金、多金属成矿区.     

Hybrid vibration experiments with a bridge foundation system model

In order to improve seismic design technology of bridges, it is necessary to evaluate the vibration characteristics of a bridge–soil system that consists of soil, foundation structure, pier and superstructure. However, there have been few experimental studies on seismic behavior of bridge–soil system. In this paper, we conducted the hybrid vibration experiment on seismic behavior of bridge–soil system, and examined the applicability of hybrid vibration experiment to study seismic response of bridge–soil system. Based on the experiment results, seismic response of bridge was quantitatively studied.     

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.     

地形变场与现代地壳垂直运动及地震活动关系初探(Ⅰ)——地形变场的计算方法和物理意义

地震与地壳形变有着密切的关系 ,但台站地倾斜资料空间连续性不好 .为此在场分析的基础上 ,利用逐步订正法将地倾斜观测资料内插到 1°× 1°的经纬度网格点 ,引入曲面总曲率的概念计算格点上的Gauss曲率 ,并推广到格点化的有限区域上得到表征曲面形态和变形强度的地表形变场 .结果表明 :该方法可以合理地弥补资料不连续的缺陷 ,并在地表形变场上捕捉到一些关键的地表形变特征 ,如地凸区、地凹区、鞍形场等 .为地震分析和预报中如何拓展和利用大范围时空信息提供了一种有益的思路 .     

阿拉尔灌区棉田蒸散量计算模型

根据2000年阿克苏水平衡站有底测坑试验资料，分析了土壤水分有效性函数与土壤相对有效含水率，作物生物学特性函数，与群体叶面积指数的关系，结果表明：（1）土壤水分有效性函数与土壤相对有效含水率呈直线函数关系；（2）作物生物学特性函数与群体叶面积指数呈指数函数关系。选用20cm蒸发器水面蒸发量、作物生物学特性函数和土壤水分有效性函数，应用数理统计方法建立了阿拉尔灌区棉田蒸散量计算模型。该模式仅需常规气象与土壤湿度资料，计算简便，精度较高，便于在缺乏实测资料的地区使用。     

209国道许家湾滑坡的成因及治理对策

结合许家湾滑坡灾害点的地质背景资料和现场勘察结果，从滑坡区的地貌、区域地质构造、地层岩性、土壤特性、水文地质特征以及人类工程活动等方面分析滑坡成因机制。指出同向多层次结构的边坡类型、山高坡陡及地表水汇集走廊的地貌特征、发育多组密集节理的区域地质构造、软硬互层的地层岩性、土体的残余强度特性等是滑坡形成的内因；而人为削坡、强降雨和破坏植被营造水耕梯田等人类活动所造成的流水浸蚀，是滑坡形成的诱发空子。基于滑坡稳定性的定量计算，确定了采用抗滑挡土墙、格构式地梁和垫梁组成抗滑系统，并布置了由树枝笼渗沟构成的排水系统等工程治理措施。现场监测结果表明滑坡治理措施卓有成效。     

二维有限元网格的自适应剖分及程序实现

有限单元法在计算岩土力学领域已得到广泛应用。在有限元的分析计算中，前处理工作即准备计算数据耗时费力。如何简化前处理工作、实现网络自动生成对于推广有限元分析至关重要。作者在综合考虑多种算法原理的基础上，对非规则几何图形、复连通域的二维模型提出了一套新的自适应剖分算法。该算法可以同时生成三角形单元和四边形单元，还能对单元形状进行规则性调整。给出了详细的实施步骤，并编写了相应的应用程序，应用程序同时提供了对单元节点进行优化排序的功能。实际应用表明，能够高效快速地给出良好的剖分网格。     

蒸发皿系数Kp计算方法研究

应用指标回归法和定性（风速和相对湿度）定量（吹程）资料，建立了蒸发皿系数Kp计算方程，并用黄河下游引黄灌区48个观测站实测资料进行了验证计算。计算结果表明，应用该方程，可大大改善Kp值的计算精度。     

改进的最小二乘法在水文分析计算中的应用

在水分析计算中，经常涉及到变量之间的线性或非线性拟合，而在拟合各种特性曲线时，通常应用以实测资料与拟合曲线间的误差平方和最小作为目标函数的方法——最小二乘法，但这种方法忽视了所有实测点应与拟合曲线间的相对误差尽量不超过某一百分比的原则，为了达到上述要求，提出了非线性的加权最小二乘法及线性相关方程的最小距离平方和法，探讨改进了传统的最小二乘法达到优化的效果。最小距离平方和法与常用的图解法相比，本法所得成果较为客观；与传统的单方向(x或y方向)最小二乘回归法相比，所求线性方程不会因坐标系的选取而改变。最后应用算例进行了初步讨论。