Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

中国兴凯湖北岸平原晚全新世花粉记录及泥炭沼泽形成

通过兴凯湖北岸平原泥炭剖面高分辩率花粉分析研究，对晚全新世花粉划分4个组合带。XKH-4组合时期(1857～1746aB．P．)为沼泽发育前期，这一时期花粉浓度小，陆生草本植物占优势，气候干冷。XKH-3组合时期(1746～1287aB．P．)为沼泽发育早期，这一时期花粉浓度较大，且水生植物花粉含量为剖面最高，喜温落叶阔叶植物大发展时期，气候温和湿润。XKH-2组合时期(1287～602aB．P．)为沼泽发育中期，这一时期花粉浓度最小，但以陆生草本植物为主，木本植物为辅，水生植物急剧减少，针叶植物出现两次高峰，气候向冷干方向发展。XKH-1组合时期(602aB．P．至今)为沼泽发育盛期，这一时期花粉浓度最大，陆生草本植物大发展时期，气候波动较大。     

常州市数字化测绘工程及GIS建库概述

数字化地形图和地籍图是城市规划、建设和管理及国土资源管理的重要基础资料．也是开展数字城市和信息化建设的最基础平台。通过对系统总体结构的设计和应用成果介绍，分析了目前基础地理数据生产、管理和应用的现状，介绍了基于GIS的常州市数字化地形图、地籍图测绘工程总体目标及基础地理信息系统总体设计．并概述了整个数字化测绘工程的实施过程和特点。     

开放数据库互联与地理信息系统

论述了开放数据库互联的基本原理 ,从技术角度探讨了在地理信息系统中应用开放数据库互联技术的实际价值 ,并建立了一个在客户 /服务器模式下两者集成的可行性模型。     

基于PSInSAR的成都地区地表沉降空间分布研究

针对成都地区土地规划利用的需要,采用永久散射体合成孔径雷达干涉测量技术(Permanent Scatterer InSAR,PSInSAR)对成都地区地表沉降空间分布情况进行了实验研究。结果表明,在该研究时段内研究区域最大沉降量达到24 mm,最大抬升量达到18 mm。研究区西北至东南方向地表在研究时段内发生了抬升,东北至西南方向出现了地表沉降现象,并对可能造成地表沉降的主要原因进行了分析。     

富互联网应用程序技术在地理信息系统开发中的应用研究

结合富互联网应用程序技术的最新进展,以Adobe公司的Flex技术为例,研究了基于四叉树的地图金字塔数据模型与瓦片双缓存技术、Flex与地理信息服务的三种后台交互技术,以及Flex的组件及其皮肤属性、动画、滤镜和FXG图形格式等各种用户体验技术在GIS开发中的应用.     

3维数字地球快速缓冲区分析算法

提出一种应用在3维数字地球中的通过图形处理器(GPU)快速实现矢量数据缓冲区分析的算法。使用一张4通道的纹理图作为容器将地理实体的矢量数据传入GPU,利用GPU的高效并行特性,将目标缓冲区纹理中的每个像素所对应的矢量坐标与原实体进行距离量算,在一次渲染中得到缓冲区纹理,最后提取出缓冲区纹理的边界。选择中国的流域和湖泊矢量数据,将本文算法与两种传统的CPU算法进行了缓冲区分析计算、测试和对比。结果显示,本文算法相对于传统矢量算法效率提高了9—16倍,相对于传统栅格算法效率提高11—20倍。实验证明,该算法计算简单,效果明显,特别是随着数据量增大,缓冲区计算速度显著优于传统算法,并能有效解决传统矢量法缓冲区分析中的数据自相交问题。     

区域规划空间信息系统建设与应用研究——以江苏省城乡规划信息系统为例

区域性城乡规划信息系统建设开展较晚,与城市级规划信息系统不同,其应用需求、数据构成、技术架构具有其独特性。结合新一轮江苏省城镇体系规划修编工作,开展全省城乡规划信息数据建库和应用平台的开发工作,系统根据区域规划的编制和研究应用需求,制定了面向未来区域规划工作需求的数据标准,集成了全省多源、多尺度、多时态的空间数据,实现了跨专业、跨平台空间信息的整合与共享。为江苏省城镇体系规划以及其他相关编制提供数据分析、规划决策支撑等服务,并预留未来与其他各类地理信息系统的应用接口。