江苏水文站网优化调整有关问题的思考

在概述全省站网现状基础上，分析了站网优化调整的主要原因，提出了如何优化调整的总体思路和基本构想，进一步明确强调在水文站网优化调整中需要坚持和把握的原则和方法。     

Centrifuge Modeling of Rock Slopes Susceptible to Block Toppling

Summary This paper presents the results of centrifuge tests that were aimed at validating the Goodman-Bray method for rock slope toppling analysis. The Goodman-Bray method was extended by the authors to accommodate non-persistent basal planes of rock columns. Two gypsum column models, with and without anchors were used to represent the failure modes. Measured critical centrifuge accelerations were in agreement with the results obtained from numerical modeling. A background of the toppling slope failures associated with a large hydropower project in China instigated the need for the centrifuge study. The centrifuge model tests used an artificial rock. The observed failure mode did not follow a straight failure plane as proposed by Goodman and Bray. The failures revealed a bi-planar slip surface with a deep-seated portion near the toe of the slope. The outcomes of the centrifuge tests illustrated the need to search for the critical failure surface when performing a toppling analysis. The search technique is similar to that usually performed in a conventional sliding analysis.     

遥测降雨抗差分析方法应用研究

遥测系统观测的雨量资料常会包含粗差，根据降雨量的特征，结合抗差理论，包为民等曾提出三步抗差修正方法。将此法应用到龙颈上水库等10个水库流域的实测资料当中，验证其有效性，并分析比较各步抗差修正的效果及降雨特征与抗差效果之间的关系。结果表明，抗差修正方法可行，且第一步和第三步结合的修正效果最好；降雨中心较明确，多年降雨特征相似，强度变幅稳定的洪水，不容易出现误抗。     

Influence of initial density of cohesionless soil on evolution of passive earth pressure

This paper focuses on the influence of the initial void ratio on the evolution of the passive earth pressure and the formation of shear zones in a dry sand body behind a retaining wall. For the numerical simulation a rigid and very rough retaining wall undergoing a horizontal translation against the backfill is considered. The essential mechanical properties of cohesionless granular soil are described with a micro-polar hypoplastic model which takes into account stresses and couple stresses, pressure dependent limit void ratios and the mean grain size as a characteristic length. Numerical investigations are carried out with an initially medium dense and initially loose sand using a homogeneous and random distribution of the initial void ratio. The geometry of calculated shear zones is discussed and compared with a corresponding laboratory model test.     

坝基集中渗漏通道的温度场反分析与试验研究

渗漏是严重影响堤坝安全的主要形式,确定渗漏通道的位置是进行堤坝渗漏整治的前提和基础。温度场探测堤坝渗漏的研究是其中一种较便捷的方法,但是在理论上还需要进一步完善。在库水和地下水存在着温度差异的基础上,考虑堤坝发生集中渗漏后,渗漏通道中的地下水与周围土体进行热量交换,造成渗流通道附近地层温度发生改变,建立起温度场对渗漏通道的反演模型。在前人研究的基础上,考虑地温梯度的影响,可以建立有集中渗漏通道影响下的二维热传导方程,进而以现场测井的温度数据作为边界条件,对泛定方程进行求解,最终得到渗漏通道的位置。计算结果和实测结果的对比表明,通过温度场反推渗漏通道参数的方法是可行的。     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

Slowness-backazimuth weighted migration: A new array approach to a high-resolution image

We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.     

NUMERICAL STUDY ON THE MIXING OF UNSORTED SEDIMENT PARTICLES DISCHARGED INTO A CROSS-FLOW BY MULTIPHASE PARTICLE-IN-CELL （MP-PIC） METHOD

The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.     

Gaussian beams in inhomogeneous media: A review

The paper outlines the most important results of the paraxial complex geometrical optics (CGO) in respect to Gaussian beams diffraction in the smooth inhomogeneous media and discusses interrelations between CGO and other asymptotic methods, which reduce the problem of Gaussian beam diffraction to the solution of ordinary differential equations, namely: (i) Babich’s method, which deals with the abridged parabolic equation and describes diffraction of the Gaussian beams; (ii) complex form of the dynamic ray tracing method, which generalizes paraxial ray approximation on Gaussian beams and (iii) paraxial WKB approximation by Pereverzev, which gives the results, quite close to those of Babich’s method. For Gaussian beams all the methods under consideration lead to the similar ordinary differential equations, which are complex-valued nonlinear Riccati equation and related system of complex-valued linear equations of paraxial ray approximation. It is pointed out that Babich’s method provides diffraction substantiation both for the paraxial CGO and for complex-valued dynamic ray tracing method. It is emphasized also that the latter two methods are conceptually equivalent to each other, operate with the equivalent equations and in fact are twins, though they differ by names. The paper illustrates abilities of the paraxial CGO method by two available analytical solutions: Gaussian beam diffraction in the homogeneous and in the lens-like media, and by the numerical example: Gaussian beam reflection from a plane-layered medium.