弹性介质分向插值无限元

无限元是一种有效的人工边界,可用于处理弹性波的传播问题。在传统动力无限元的基础上,提出了一种采用分向插值技术的新型动力无限元,详细地推导了这种无限元的形函数,建立了完全解析形式的刚度矩阵,以提高计算效率,采用该无限元边界,计算了弹性介质中的线源Lamb问题,通过对比解析解答的地基表面位移,验证了该无限元的有效性。算例分析表明,采用此类无限元时,有限元单元边长建议取不超过1/8剪切波波长,网格边界到激励源点的距离宜取5倍剪切波波长。无限单元中的幅值衰减系数对计算结果影响甚微,建议取较小值。     

Uniform asymptotic Green's functions for efficient modeling of cracks in elastic layers with relative shear deformation controlled by linear springs

We present a uniform asymptotic solution (UAS) for a displacement discontinuity (DD) that lies within the middle layer of a three‐layer elastic medium in which relative shear deformation between parallel interfaces is controlled by linear springs. The DD is assumed to be normal to the two interfaces between the elastic media. Using the Fourier transform method we construct a leading term in the asymptotic expansion for the spectral coefficient functions for a DD in a three‐layer‐spring medium. Although a closed‐form solution will require a solution in terms of an infinite series, we demonstrate how this UAS can be used to construct highly efficient and accurate solutions even in the case in which the DD actually touches the interface. We compare the results using the Green's function UAS solution for a crack crossing a soft interface with results obtained using a multi‐layer boundary element method. We also present results from an implementation of the UAS Green's function approach in a pseudo‐3D hydraulic fracturing simulator to analyze the effect of interface shear deformation on the fracture propagation process. These results are compared with field measurements. Copyright © 2008 John Wiley & Sons, Ltd.     

瑞利波特性刚度矩阵分析方法

沿水平方向从半无限体分层系统截取出有限厚度分层固体系统,截取边界设置为刚性边界。将固体系统分层细化,由刚度矩阵方法可得到表面波在该系统频散特性,利用瑞利波位移在底层衰减特性从表面振型波可筛选出瑞利振型波,这些瑞利振型波与某一频率范围半无限体分层系统瑞利振型波对应,其频率下限与截取厚度有关,截取厚度越大,频率下限越小。同时还分析了分层细化厚度、截取厚度与分析频率关系。     

利用时间分裂的错格伪谱法模拟地震波在基于改进BISQ模型的双相介质中的传播

改进的BISQ(Biot-Squirt)模型中各参数具有明确的物理意义和可实现性,在不引入特征喷流长度的情况下可将Biot流动和喷射流动两种力学机制有机地结合起来;而高精度的地震波场数值模拟技术是研究双相介质地震波传播规律的重要手段。本文从本构方程、动力学方程和动力学达西定律出发,推导了基于改进BISQ模型的双相各向同性介质的一阶速度--应力方程组;采用时间分裂错格伪谱法求该方程组的数值解,模拟半空间及层状双相介质中的地震波场。数值模拟结果表明:①与传统方法相比,时间分裂错格伪谱法波场数值模拟的精度更高,压制网格频散效果更好;②在非黏滞相界情况下,慢纵波呈传播性,而在黏滞相界情况下,慢纵波呈扩散性,以静态模式出现在震源位置;③双相介质分界面处,各类波型复杂的反射透射规律可由数值模拟结果清晰展现。     

Source signature and elastic waves in a half‐space under a momentary shear line impulse

The transient deformation of an elastic half‐space under a line‐concentrated impulsive vector shear load applied momentarily is disclosed in this paper. While in an earlier work, the author gave an analytical–numerical method for the solution to this transient boundary‐value problem, here, the resultant response of the half‐space is presented and interpreted. In particular, a probe is set up for the kinematics of the source signature and wave fronts, both explicitly revealed in the strained half‐space by the solution method. The source signature is the imprint of the spatiotemporal configuration of the excitation source in the resultant response. Fourteen wave fronts exist behind the precursor shear wave S: four concentric cylindrical, eight plane, and two relativistic cylindrical initiated at propagating centres that are located on the stationary boundaries of the solution domain. A snapshot of the stressed half‐space reveals that none of the 14 wave fronts fully extend laterally. Instead, each is enclosed within point bounds. These wave arresting points and the two propagating centres of the relativistic waves constitute the source signature. The obtained 14 wave fronts are further combined into 11 disparate wave fronts that are grouped into four categories: an axis of symmetry wave—so named here by reason of being a wave front that is contiguous to the axis of symmetry, three body waves, five surface waves and two inhibitor waves—so named here by reason that beyond them the material motion dies out. Of the three body waves, the first is an unloading shear wave, the second is a diffracted wave and the third is a reflected longitudinal two‐branch wave. Of the two inhibitor waves, the first is a two‐joint relativistic wave, while the second is a two‐branch wave. The wave system, however, is not the same for all the dependent variables; a wave front that appears in the behaviour of one dependent variable may not exist in the behaviour of another. It is evident from this work that Saint–Venant's principle for wave propagation problems cannot be formulated. Therefore, the above results are valid for the particular proposed model for the momentary line‐concentrated shear load. The formulation of the source signature, the wave system, and their role in the half‐space transient deformation are presented here. Copyright © 2003 John Wiley & Sons, Ltd.     

The physical Properties of marine sediments as inferred from Wavelet Correlation Analysis

Seismic waves traveling in the water/sediment or sub-bottom sediment interface have been the subject of considerable interest in underwater acoustics in recent years. Some progress has been made in understanding the propagation and attenuation characteristics of interface waves in different geological environments. However, the generating mechanisms are poorly understood. In particular, what is the acoustic-seismic energy conversion process? As seismic waves involve both time and space parameters it should be able to relate directly the propagation characteristics of the ocean bottom interface waves to the shear properties of the sediments over the propagation area. To address these problems we have applied the wavelet correlation method (WCR) to examine the variations of bottom characteristics and their role in coupling waterborne sound into the sea bottom. To confirm the validity of the developed modeling technique, we applied the wavelet correlation analysis for synthetic seismograms and field data. In this discussion paper we present images of the first and second shear modes and the interface wave component as a function of arrival time and frequency. We also discuss the possibility of inverting phase and group velocity information directly from the wavelet cross-correlation function and propose the way of using the WCR method to predict physical dynamic parameters of marine sediments.     

Theory of Volumetric Moving Dislocation in Poroelastic Halfspace and Characterization of Magma Intrusion Events

The undrained change in pore fluid pressure that accompanies dike intrusion may be conveniently represented as a moving volumetric dislocation. The concept of a dilation center was developed to represent the field of undrained pressure change in a saturated linear elastic medium. Since instantaneous pore fluid pressures can be developed to a considerable distance from the dislocation, monitoring the rate of pressure generation and subsequent pressure dissipation in a fully coupled manner enables certain characteristics of the resulting dislocation to be defined. The principal focus of this study is the application of dislocation-based methods to analyze the behavior of the fluid pressure response induced by intrusive dislocations in a semi-infinite space, such as dike intrusion, hydraulic fracturing and piezometar insertion. Partially drained pore pressures result from the isothermal introduction of volumetric moving pencil-like dislocations described as analogs to moving point dislocation within a semi-infinite saturated elastic medium. To represent behavior within the halfspace, an image dislocation is positioned under the moving coordinate frame fixed to the front of the primary moving dislocation, to yield an approximate solution for pore pressure for constant fluid pressure conditions. Induced pore pressures are concisely described under a minimum set of dimensionless parameter groupings representing propagation velocity, and relative geometry. Charts defining induced pore fluid pressure at a static measuring point provide a meaningful tool for determining unknown parameters in data reduction. Two intrusive events at Krafla, Iceland are examined using the type curve matching techniques. Predicted parameters agree favorably with field data.     

2.5D Simulation of basin-edge effects on the ground motion characteristics

The effects of basin-edge and soil velocity on the ground motion characteristics have been simulated using 2.5D modeling. One of the most significant advantages of the 2.5D simulation is that 3D radiation pattern can be generated in a 2D numerical grid using double-couple shear dislocation source. Further, 2.5D numerical modeling avoids the extensive computational cost of 3D modeling. The responses of basin-edge model using different soil velocities revealed that surface waves were generated near the edge of the basin and propagated normal to the edge, towards the basin. Further, the results depict increase of amplification, duration and surface wave generation with the decrease in soil velocity.     

线电流源声频大地电磁测深的二维正演计算及响应特点

论述了如何在大地电磁测深二维正演计算的基础上，结合线电流源声频大地电磁测深（ＣＳＡＭＴ）的边界条件，利用三角形６节点有限元法计算ＴＥ极化方式ＣＳＡＭＴ二维响应的方法；给出了４个模型的计算结果     

海底地震仪(OBS)主动源转换横波的应用

横波（S波）为偏振波,具有不同于纵波（P波）的特性,对于地震探测具有特殊的意义.在被动源地震探测中已得到广泛的应用,如接收函数、S波分裂等.在主动源（气枪）海底地震（OBS）探测中,震源在水中,S波为地层转换波,其应用还不多.本文在介绍转换S波的产生、模式、处理和识别的基础上,以实例为切入点简述其具体应用.这些应用主要是基于1D/2D转换S波,用于揭示海底岩石类型、推断地壳性质、共轭陆缘问题、判定地幔蛇纹石化、估算天然气水合物的饱和度和预测流体等.目前在南海已获得了大量的2D和3D的OBS转换S波数据,可将转换S波的研究逐步从2D发展到3D研究,同时结合其他地球物理资料进行共同分析.利用转换S波的研究,有利于揭示南海扩张停止后形成的海山下不同地层的岩性和判定上地幔低速的性质等.      

On the strength of transversely isotropic rocks

Accurate prediction of strength in rocks with distinct bedding planes requires knowledge of the bedding plane orientation relative to the load direction. Thermal softening adds complexity to the problem since it is known to have significant influence on the strength and strain localization properties of rocks. In this paper, we use a recently proposed thermoplastic constitutive model appropriate for rocks exhibiting transverse isotropy in both the elastic and plastic responses to predict their strength and strain localization properties. Recognizing that laboratory‐derived strengths can be influenced by material and geometric inhomogeneities of the rock samples, we consider both stress‐point and boundary‐value problem simulations of rock strength behavior. Both plane strain and 3D loading conditions are considered. Results of the simulations of the strength of a natural Tournemire shale and a synthetic transversely isotropic rock suggest that the mechanical model can reproduce the general U‐shaped variation of rock strength with bedding plane orientation quite well. We show that this variation could depend on many factors, including the stress loading condition (plane strain versus 3D), degree of anisotropy, temperature, shear‐induced dilation versus shear‐induced compaction, specimen imperfections, and boundary restraints.     

SILENT BOUNDARY CONDITIONS FOR WAVE PROPAGATION IN SATURATED POROUS MEDIA

Wave propagation both in one- and in two-dimensional saturated elastic porous media is analysed by means of a two-field finite element model with silent boundaries. An extension of the elastic ‘multidirectional’ transmitting boundary to two-phase media is developed to simulate the silent boundary condition. The theoretical assessment and the numerical formulation of the first-order silent boundary technique is presented in detail. Some examples are used to demonstrate the reliability of the first-order method, especially for problems with plane and axisymmetric waves having various angles of incidence. Finally, the wave propagation along a pile shaft is presented, to simulate a common non-destructive dynamic pile test.     

煤矿巷道内地震勘探波场特征

为了研究巷道内地震波的波场特征,识别出波场中的各种波,利用交错网格高阶有限差分正演方法模拟了煤矿巷道内点震源激发及平面波震源激发产生的弹性波场。通过波场分析得出煤矿巷道内的共炮点道集中除了有常规地面地震中的各种波外,还有在巷道内上下界面的多次反射波和多次折射转换波,是巷道内地震勘探的主要干扰波。利用平面波震源在巷道内激发产生的波场中,巷道下方的空洞在平面波震源激发的地震记录上,其多次反射绕射波互相叠合在一起,表现为“糖葫芦”状的明显异常特征。这些都为巷道内地震勘探资料处理和解释提供了重要依据。      

On the mechanics of contact between a flexible transducer diaphragm located at a rigid boundary and an elastic material

In this paper we examine the mechanics of contact between a flexible elastic diaphragm and an isotropic elastic halfspace region where contact is established through the application of a normal stress. It is shown that a convenient solution for the maximum displacement and maximum flexural moment in the flexible diaphragm can be obtained by employing a variational technique. Copyright © 2006 John Wiley & Sons, Ltd.     

含砾滑带土三维颗粒流模型建模方法研究

提出了基于CT扫描技术和PFC三维颗粒流程序建立含砾滑带土三维建模方法。首先,通过CT扫描技术获得含砾滑带土三轴试样的一系列二维切片图像,采用二值化和砾石边界识别方法,得到了砾石边界点云数据;其次,利用砾石边界点云数据,采用逆向重构得到各个砾石的三维模型;最后,将重构的三维砾石模型导入PFC三维颗粒流程序,建立颗粒流数值模型。计算结果表明:三轴剪切过程中,砾石的欧拉角随着试样剪切变形的发展而变化;与无砾石试样相比,砾石周围土体小球的运动方向更加杂乱,说明了砾石对含砾滑带土的剪切破坏有显著的影响。通过与三轴剪切试验过程中砾石的空间运动规律进行对比分析,验证了提出的建模方法是可行的。采用该方法研究砾石的空间运动特征,对于揭示砾石在含砾滑带土剪切破坏过程中的作用具有重要的意义。     

Earthquake activity related to seismic cycles in a model for a heterogeneous strike-slip fault

We investigate the evolution of seismicity within large earthquake cycles in a model of a discrete strike-slip fault in elastic solid. The model dynamics is governed by realistic boundary conditions consisting of constant velocity motion of regions around the fault, static/kinetic friction and dislocation creep along the fault, and 3D elastic stress transfer. The fault consists of brittle parts which fail during earthquakes and undergo small creep deformation between events, and aseismic creep cells which are characterized by high ongoing creep motion. This mixture of brittle and creep cells is found to generate realistic aftershock sequences which follow the modified Omori law and scale with the mainshock size. Furthermore, we find that the distribution of interevent times of the simulated earthquakes is in good agreement with observations. The temporal occurrence, however, is magnitude-dependent; in particular, the small events are clustered in time, whereas the largest earthquakes occur quasiperiodically. Averaging the seismicity before several large earthquakes, we observe an increase of activity and a broadening scaling range of magnitudes when the time of the next mainshock is approached. These results are characteristics of a critical point behavior. The presence of critical point dynamics is further supported by the evolution of the stress field in the model, which is compatible with the observation of accelerating moment release in natural fault systems.     

P波作用下跨断层隧道轴线地震响应分析

基于近场波动有限元方法并结合黏弹性人工边界条件,针对3D断层场地,通过求解等效二维场地地震响应,从而获得3D场地地震动输入的自由场响应,并将自由场响应转化为3D模型边界面上的等效节点力,从而建立含断层3D场地P波入射的倾斜输入方法。自由场算例验证所提方法具有较好精度,进而基于建立的输入方法,开展了跨断层隧道地震响应的数值模拟研究。数值模拟结果表明:在P波作用下,隧道跨断层部位处于拉、压、剪切的复杂受力状态,且断层处的隧道衬砌地震响应明显大于其他部位的地震响应;围岩的力学性质与断层的力学性质相差越大,断层处衬砌的地震响应放大越明显;断层处衬砌地震响应随断层深度的增加而增加。另外,跨断层隧道的地震响应受P波入射角度的影响较大,随P波入射角度的增加,断层处隧道衬砌的轴力、弯矩先增加后减小,而剪力具有逐渐减小的规律。     

Boundary element analysis of transversely isotropic bi-material halfspaces with inclined planes of isotropy and interfaces

In this paper, a single-region boundary element method (BEM) is presented for the analysis of transversely isotropic bi-material halfspaces with arbitrarily inclined planes of isotropy and material interfaces. The proposed BEM uses the fundamental solution of a transversely isotropic bi-material fullspace and five boundary element techniques. Infinite boundary elements are introduced to consider the far-fields of a transversely isotropic bi-material halfspace. The effective integration methods are proposed for dealing with various integrals in the discretized boundary integral equation. The stresses at internal points are obtained using the coordinate transformation of kernel functions, and the stresses on the boundary surface are calculated using an improved traction recovered method. Numerical verifications of displacements and stresses for a benchmark problem are conducted, and excellent agreement with previously published results is obtained. Numerical examples are presented to illustrate the influence of non-horizontal or horizontal planes of isotropy in bi-material halfspaces on the displacements and stresses induced by the tractions on the horizontal boundary surface. Results reveal that the elastic fields vary clearly with the dip angle of the isotropic plane and the stresses across the bi-material interface are closely related to the ratios of the elastic parameters of the bi-material.     

边坡模糊可靠性分析隶属函数取值界限研究

针对岩土参数取值所具有的随机性、模糊性与区间性特征,建议采用隶属函数表征参数取值,分别构造了岩土参数三角形分布、正态分布与拟正态分布隶属函数,并建立了采用隶属函数取值界限表示的不同截集水平岩土参数区间值确定方法,提出分别采用均值加减2.5倍标准差与3倍标准差确定岩土参数三角形分布与拟正态分布隶属函数取值界限,进而以此为基础采用模糊点估计分析方法对工程算例进行分析,结果表明,在截集水平个数为11时,模糊点估计方法所得边坡稳定模糊可靠性指标与蒙特卡洛法所得随机概率可靠性指标精度相当,截集水平个数为9时,评价结果偏安全,建议在工程分析时,截集水平取0.1⁰.9间的9个,该方法计算过程直接简单,计算代价低,更具工程实用性。     

The method of fundamental solution for 3‐D wave scattering in a fluid‐saturated poroelastic infinite domain

By using a complete set of poroelastodynamic spherical wave potentials (SWPs) representing a fast compressional wave P_I, a slow compressional wave P_II, and a shear wave S with 3 vectorial potentials (not all are independent), a solution scheme based on the method of fundamental solution (MFS) is devised to solve 3‐D wave scattering and dynamic stress concentration problems due to inhomogeneous inclusions and cavities embedded in an infinite poroelastic domain. The method is verified by comparing the result with the elastic analytical solution, which is a degenerated case, as well as with poroelastic solution obtained using other numerical methods. The accuracy and stability of the SWP‐MFS are also demonstrated. The displacement, hoop stress, and fluid pore pressure around spherical cavity and poroelastic inclusion with permeable and impermeable boundary are investigated for incident plane P_I and SV waves. The scattering characteristics are examined for a range of material properties, such as porosity and shear modulus contrast, over a range of frequency. Compared with other boundary‐based numerical strategy, such as the boundary element method and the indirect boundary integral equation method, the current SWP‐MFS is a meshless method that does not need elements to approximate the geometry and is free from the treatment of singularities. The SWP‐MFS is a highly accurate and efficient solution methodology for wave scattering problems of arbitrary geometry, particularly when a part of the domain extends to infinity.