Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

岩质高边坡岩体变形参数及松弛带厚度研究

预测岩质高边坡开挖后岩体变形模量的变化及松弛带厚度，是分析岩质高边坡在开挖后变形（位移）和作好防护设计的重要资料，运用波动力学关于平均应力与体积模量、岩体纵波速度与弹性模量、变形模量间的关系，通过部分实测资料及边坡应力场有限元分析的资料，分别建立了纵波速度与岩体变形模、岩体应力间的关系，研究了开挖边坡岩体变形模量的变化，预测了岩体松弛带的厚度。     

模拟试验法确定桩基承载力

通过室内中型剪切摩擦试验确定桩与地层间的C、Φ值,依据莫尔—库仑准则确定桩侧摩阻力。桩侧法向应力采用桩土上部荷载引起的作用于桩身的侧向应力,桩端承载力通过室内三轴试验确定,两者之和即为单桩承载力。该方法用于肇源松花江大桥桩基承载力的测试,得出的桩承载力与实测值吻合较好,取得了预期的效果。     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

A Modeling Study of the Effects of Anomalous Snow Cover over the Tibetan Plateau upon the South Asian Summer Monsoon

The effect of anomalous snow cover over the Tibetan Plateau upon the South Asian summer monsoon is investigated by numerical simulations using the NCAR regional climate model (RegCM2) into which gravity wave drag has been introduced. The simulations adopt relatively realistic snow mass forcings based on Scanning Multi-channel Microwave Radiometer (SNINIR) pentad snow depth data. The physical mechanism and spatial structure of the sensitivity of the South Asian early summer monsoon to snow cover anomaly over the Tibetan Plateau are revealed. The main results are summarized as follows. The heavier than normal snow cover over the Plateau can obviously reduce the shortwave radiation absorbed by surface through the albedo effect, which is compensated by weaker upward sensible heat flux associated with colder surface temperature, whereas the effects of snow melting and evaporation are relatively smaller.The anomalies of surface heat fluxes can last until June and become unobvious in July. The decrease of the Plateau surface temperature caused by heavier snow cover reaches its maximum value from late April to early May. The atmospheric cooling in the mid-upper troposphere over the Plateau and its surrounding areas is most obvious in May and can keep a fairly strong intensity in June. In contrast, there is warming to the south of the Plateau in the mid-lower troposphere from April to June with a maximum value in May.The heavier snow cover over the Plateau can reduce the intensity of the South Asian summer monsoon and rainfall to some extent, but this influence is only obvious in early summer and almost disappears in later stages.     

剔除高密度电阻率法三极装置浅部不均匀体效应影响的方法及效果

基于总结出的高密度电阻率法三极测量装置下的二维正演拟断面中浅部不均匀体畸变规律和特征,设计了一套能有效剔除三极测量装置浅部不均匀体效应影响的数据处理方法,成功应用于野外实测数据的处理。     

巫山县污水处理厂高填方地基湿化变形试验研究

结合巫山县污水处理厂高填方工程，选取压实粘性土样，在室内进行了接近实际应力路径下的“单线法”湿化变形试验，研究了土样在不同的压实度、固结压力和湿化历史下的湿化机理、湿化前后的应力．应变关系及强度变化规律。取得的湿化变形试验成果表明，湿化使土样的变形模量和强度均有所降低；湿化变形不仅包括湿化体应变，而且包括湿化偏应变；随着压实度的增加，湿化变形减小。试验成果为巫山县污水处理厂高填方地基的填筑设计、控制湿化变形的措施和附加压密措施提供了依据。     

多孔介质非饱和导水率预测的分形模型

多孔介质非饱和导水率是地下水污染预测与评价的重要参数。根据分形几何的基本原理和方法,推导出了与Campbell经验公式在形式上完全一致的多孔介质非饱和导水率的预测公式。公式中的幂指数为介质孔隙分维和随机行走分维的函数,分别体现了多孔介质的静态性质与动态性质对其中水分运动的影响,但静态性质的影响是主要的,即导水率主要受多孔介质的结构控制。根据文献中报道的大量数据,利用笔者推导的预测公式计算得到的幂指数的统计值与试验测定的幂指数的统计值基本一致,说明推导的理论公式预测多孔介质非饱和导水率是较为可靠的。