A multiscale work-analysis approach for geotechnical structures

This paper presents a second-order work analysis in application to geotechnical problems by using a novel effective multiscale approach. To abandon complicated equations involved in conventional phenomenological models, this multiscale approach employs a micromechanically-based formulation, in which only four parameters are involved. The multiscale approach makes it possible a coupling of the finite element method (FEM) and the micromechanically-based model. The FEM is used to solve the boundary value problem (BVP) while the micromechanically-based model is utilized at the Gauss point of the FEM. Then, the multiscale approach is used to simulate a three-dimensional triaxial test and a plain-strain footing. On the basis of the simulations, material instabilities are analyzed at both mesoscale and global scale. The second-order work criterion is then used to analyze the numerical results. It opens a road to interpret and understand the micromechanisms hiding behind the occurrence of failure in geotechnical issues.     

地面地震层析技术在金属矿勘查中的试验研究

由于金属矿区地质构造的复杂性和复杂的地表条件, 适合于探测层状介质的地震反射方法技术在金属矿勘查中的应用受到了限制。为提高地震记录的横向分辨率, 探测复杂的地质构造及隐伏岩体等, 需研究相应的地震方法技术, 地面地震层析技术在金属矿勘查中的试验研究就是为解决上述问题而开展的。本文就非线性地面地震层析技术的基本原理进行了阐述, 并就该方法在金属矿勘查中的试验效果进行了讨论。     

利用钻孔超声成像的图像特征进行岩石风化程度分类

钻孔超声成像方法能精确地反映井壁的光滑情况, 进而判断岩石的风化发育程度, 确定风化层的厚度。阐述了钻孔超声成像的走时和振幅图像特征, 给出了1个在香港地区的应用实例。     

2-D crustal structure of Changbaishan-Tianchi volcanic region determined by seismic traveltime inversion

IntroductionTrial-and-error forward modeling of wide-angle seismic reflection/refraction traveltimes for 2-D velocity structure is extremely time-consuming, even for experienced data interpreters. For wide-angle seismic reflection/refraction experiments that consist of numerous shots along a single line, it is quite difficult through repeated trial-and-error forward modeling to construct a 2-D model that fits the data within acceptable limits (Cerveny, et al, 1977; ZHANG, et al, 200 . In ad…     

Multiparametric Traveltime Inversion

In conventional seismic processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macrovelocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with normal moveout (NMO) velocities, as derived from a previous velocity analysis conducted on common midpoint (CMP) data. This work presents a modified version of the Hubral and Krey algorithm that is designed to extend the original version in two ways, namely (a) it makes an advantageous use of previously obtained common-reflection-surface (CRS) attributes as its input and (b) it also allows for gradient layer velocities in depth. A new strategy to recover interfaces as optimized cubic splines is also proposed. Some synthetic examples are provided to illustrate and explain the implementation of the method.     

基于矩形网格的有限差分走时计算方法

对于大多数速度场，地震波沿射线传播的初至波走时，可以用有限差分外推的方法在二维或三维数值网格上计算出来. 在保证精度的条件下，为提高计算效率和适应性，本文推导了基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法. 另外，该方法对首波和散射波做了合适的处理，而且不会碰到传统射线法存在的阴影区和焦散区等问题. 简单模型和复杂的Marmousi模型试算的结果表明，该方法精度较高并适用于强纵、横向变速的复杂介质. 基于该方法的Kirchhoff叠前深度偏移, 在主要构造和目的层位置的成像效果上基本达到了波动方程法叠前深度偏移的位置成像效果. 由于未考虑续至波等有效能量，在成像的保幅性上不如波动方程法叠前深度偏移的效果，但其计算效率则明显高于全格林函数法和波动方程法.      

Resolution Sensitivity and Scaling of Large-Eddy Simulations of the Stable Boundary Layer

Large-eddy simulations (LES) of the continuously turbulent quasi-equilibrium stable boundary layer (SBL) are conducted with grid lengths in the range of 12.5 m to 2 m, in order to explore resolution sensitivity, and determine at what point grid convergence occurs. The structure of the mean potential temperature, winds, and turbulent fluxes varies significantly over this resolution range. The highest resolution simulations show a significant degree of convergence. The dimensionless momentum diffusivity asymptotes to a value of 0.06, corresponding to a limiting flux Richardson number of 0.15.Using the converged simulations, some scaling hypotheses underpinning first-order and second-order closure models are revisited. The effective Richardson number stability functions of the LES are compared with the forms often used in numerical weather prediction (NWP). The mixing implied by the LES is less than that used in NWP. The commonly used similarity profiles for heat and momentum fluxes, and the scalings for dissipation and pressure covariances are compared with the LES. This information could provide guidance for the next generation of SBL parametrization schemes.