Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

SAMPLING REMOTELY SENSED IMAGERY FOR STORAGE,RETRIEVAL, AND RECONSTRUCTION*

A major difficulty in remote sensing is handling the many data from sensors aboard aircraft and satellites. In this paper we identify an optimal procedure for sampling remotely sensed data before their storage or on their retrieval. The procedure depends on spatial correlation in the scene and uses kriging to estimate values that have been lost. An example in which data from an airborne multispectral scanner could be diminished to only about one tenth without serious loss of precision illustrates the method.     

重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. Copyright © 2003 John Wiley & Sons, Ltd.     

中国开辟区不利开采地形类型及其分布

利用1998年"大洋一号"调查船近海底作业所获的深拖海底电视和照相系统资料,对中国开辟区不利开采地形类型及其分布特征进行了研究,并对海底地形的坡度进行了计算,结果表明:(1)中国开辟区东、西两区海底的东向坡和西向坡是均衡分布的。东、西两区海底地形坡度的分布特征明显不同,东区海底地形坡度以0～15°为主,而西区坡度大于3°的地形则很少。坡度大于5°的地形主要分布于坡脚处,特别是坡度超过15°的地形总是出现在水深变化的转折点处。(2)研究区不利开采地形可分为坡度为5～15°和大于15°两种类型,东区地形坡度大于5°的不利开采地形比西区多14%,东区坡度大于15°的地形为3.6%,比西区(仅为1.4%)多,东、西两区都有少量坡度近90°的地形出现,并且这种地形在东区所占比率比西区大,说明东区断崖、陡坎和断裂沟槽比西区多,地形也更为复杂多变。     

地形坡度对多金属结核分布的控制作用

选取我国东太平洋多金属结核开辟区内的一个区域作为研究区,利用人工神经网络中应用最为广泛的BP网络,建立控制多金属结核分布的地质因素与多金属结核分布之间的映射关系,探讨地形坡度对多金属结核分布的控制作用.结果表明,在经度、纬度、水深、坡度四个因素中,坡度对结核分布影响程度最大;多金属结核主要分布于坡度小于5°的地方;当坡度小于5°时,丰度与地形坡度呈正相关,品位与地形坡度呈负相关,丰度与品位呈负相关;品位与坡度似呈指数关系,坡度增大,品位降低.     

应用数字地形模型方法分析漳州港海底冲淤变化

利用漳州港海域无围堰填海造地前后的水深测量数据，结合海域地质钻探、表层和柱状沉积物资料，用数字地形模型方法计算了该区海底的冲淤量，分析了冲淤变化及其影响因素，结果显示：（1）淤积主要发生在近岸海域，其主要原因是近岸3-10m厚的海底淤泥层在回填土的挤压和推移下发生形变，产生淤浅假象；另一个原因是回填时部分泥沙流失入海，淤积在近岸海域，致使淤积强度增大。（2）由近岸往外，淤积强度逐渐减弱，并过渡为产中刷状态。     

Triple diagram method for the prediction of wave height and period

Many formulations have been developed so far to predict the wave height and period from fetch length and wind blowing duration for a constant wind speed. This study aimed to predict wave parameters from fetch length and meteorological factors by using triple diagram methodology based on Kriging principles. Proposed model results were compared with Joint North Sea Wave Project (JONSWAP) model which is used so commonly in the ocean and coastal engineering studies. For the implementation of the methodology hourly wave and wind data were obtained from a buoy located in Lake Ontario. Numerical and graphical comparisons demonstrated that the proposed method outperforms the classical formulation.     

Separation of Waves by Interpolation Method

1 .ＩｎｔｒｏｄｕｃｔｉｏｎＴｏｍｅａｓｕｒｅｔｈｅｒｅｆｌｅｃｔｉｏｎｏｆｉｎｃｉｄｅｎｔｗａｖｅｓｐｒｏｄｕｃｅｄｂｙａｐｈｙｓｉｃａｌｍｏｄｅｌｉｎａｗａｖｅｆｌｕｍｅｉｓａｃｏｍｍｏｎｐｒｏｂｌｅｍｗｉｔｈｐｈｙｓｉｃａｌｔｅｓｔｓｏｆｗａｖｅａｃｔｉｏｎｏｎｃｏａｓｔａｌｓｔｒｕｃｔｕｒｅｓ .Ｗａｖｅｓｇｅｎｅｒａｔｅｄｂｙｔｈｅｗａｖｅ ｍａｋｅｒｐｒｏｐａｇａｔｅｆｏｒｗａｒｄｉｎｔｈｅｗａｖｅｆｌｕｍｅａｎｄａｒｅｒｅｆｌｅｃｔｅｄｂｙｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌ,ａｎｄｔｈｅ…