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.     

滇西动态重力网的分形特征及空间分辨力

介绍了用分形学描述地球物理观测网非均匀分布和联合应用分形特征与Ｓｈａｎｎｏｎ定理确定最佳网格距的方法。研究了滇西动态重力网的分形特征，它的分形维Ｄｆ＝１．２７，最佳网格距ｒ０＝２８ｋｍ。讨论了１９８６年至１９９４年间网区发生的８次Ｍｓ＞５．０级地震前重力场变化的可信性，指出这些变化的波长均大于１００ｋｍ，至少２倍于网的最小空间分辨力，因而是可信的     

重力插值方法研究

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

地球化学景观吸引子的发现及其意义

勘查地球化学对异常评价的最大不确定性源于对地球化学景观特性的认识不够。我们应用构造相空间的方法分析了浙江儒岙幅（1:5万）水系沉积物测量的Ag含量数据，揭示了地球化学景观是具有低维吸引子（DZ≈2．9）的混饨系统。这一发现将对勘查地球化学中采样密度的选择、地球化学场模拟、异常评价等根本性问题产生重大影响。     

On the practice of estimating fractal dimension

Coastlines epitomize deterministic fractals and fractal (Hausdorff-Besicovitch) dimensions; a divider [compass] method can be used to calculate fractal dimensions for these features. Noise models are used to develop another notion of fractals, a stochastic one. Spectral and variogram methods are used to estimate fractal dimensions for stochastic fractals. When estimating fractal dimension, the objective of the analysis must be consistent with the method chosen for fractal dimension calculation. Spectal and variogram methods yield fractal dimensions which indicate the similarity of the feature under study to noise (e.g., Brownian noise). A divider measurement method yields a fractal dimension which is a measure of complexity of shape.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

多波束数据的海底数字地形模型构建

提出大批量多格式原始多波束数据的DTM构建方法，以满足大区域、大比例尺海底地形制图的需要。在多种格式原始多波束数据接口和系统内部标准数据结构的基础上，通过对数据文件、数据种类和数据运算量的有效组织管理，实施边读入、边权重配赋的网格插值，分析了高斯、指数和距离平方反比权重函数的适用性，及最小值、最大值和平均值的实用性。在权重配赋网格插值基础上，提出分形fBM和张力样条配合使用的方程式插值方法，保证DTM数据的有效外延和地形分辨。整套算法效率高，并能有效保证DTM的精度，对存储在外部介质的数据遍历一次即可完成网格插值。     

富钴结壳分布的分形特征

利用我国在西太平洋国际海底两个海山获得的富钴结壳资料进行分形研究后发现，富钴结壳的分布具有变维分形分布的特点，结壳丰度呈二阶累计和分形分布，结壳厚度呈多阶不定累计和分形分布，通过计算得到了它们各自的分维数，进而对两个海山的结壳资源进行了初步评价。