An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore‐water pressure

The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd.     

Rockburst research at Falconbridge's Strathcona Mine,Sudbury, Canada

Microseismic systems at five mines in the Sudbury Basin provide the basic data for Falconbridge Limited's rockburst research. Daily and long-term analysis of this data as well as underground observations have confirmed the fault-slip mechanism at three mines. A detailed analysis of the complete history of Falconbridge Mine is being conducted and Distinct Element numerical models are being used to simulate both the stick-slip behaviour of faults and the dynamic effects of the induced vibrations on rock and backfill. University research includes acoustic tomographic imaging of the rock mass based on seismic wave propagation and collection of full microseismic waveforms to allow application of advanced seismic and statistical analysis techniques.Formerly with Mines Technical Services, Falconbridge Limited, Sudbury Operations, Falconbridge , Ont., P0M 1S0.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal, Canada, August 30, 1987.     

Laboratory and field investigations of rockburst phenomena using concurrent geotomographic imaging and acoustic emission/microseismic techniques

The progress made on three phases of a research project, started in 1986 to investigate mining induced seismicity/rockburst phenomena using concurrent geotomographic imaging and microseismic monitoring techniques, is described. Phase I is the geotomographic software development and laboratory calibration trials. Phase II is the enhancement of traditional microseismic monitoring instrumentation with a waveform acquisition system, so that source mechanism studies can be carried out on mining induced seismic events. Phase III is the field trials of the hybrid technique which will be used to monitor changing rock mass physical properties, in response to mining. Preliminary results from all three phases are given, together with an outline of current and future research planned.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal 1987.     

Complete synthetic seismograms for high-frequency multimodeSH-waves

We present an efficient scheme to compute high-frequency seismograms (up to 10 Hz) forSH-waves in a horizontally stratified medium with the mode summation method. The formalism which permits the computation of eigenvalues, eigenfunctions and related integral quantities is discussed in detail. Anelasticity is included in the model by using the variational method. Phase velocity, group velocity, energy integral and attenuation spectra of a structure enable the computation of complete strong motion seismograms, which are the basic tool for the interpretation of near-source broad-band data.Different examples computed for continental structures are discussed, where one example is the comparison between the observed transversal displacement recorded at station IVC for the November 4, Brawley 1976 earthquake and synthetic signals. In the case of a magnitudeM _L =5.7 earthquake in the Friuli seismic area we apply the mode summation method to infer from waveform modeling of all three components of motion of observed data some characteristics of the source.On leave from Institute of Geophysics, ETH, Zürich, Switzerland.     

一步彩虹全息术的成象位置分析

本文从标量衍射理论出发,结合一定的物理意义,直接导出了一步彩虹全息术的成象位置公式;并根据可观察再现象的条件,从理论上导出了较为完整的成象位置关系表;并从中选出一部分成象位置关系进行了实验验证。     

频谱成像技术在稠油热采地震监测中的应用

频谱成像技术可有效的描述地质反射层厚度的非连续性和岩性的非均质性,其在理论上主要是依据薄层反射的调谐原理,过去通常采用以离散傅里叶变换为基础的算法,但是,该方法存在着明显的局限性,因为估算的地震振幅谱的重要特征是所选时窗长度的函数.如果所选时窗过短,振幅谱会与变换窗函数褶积,失去频率的局部化特征,而且过短的时窗会使子波的旁瓣呈现为单一反射的假象.增加时窗长度,会改善频率的分辨率.但如果所选时窗过长,时窗内的多个反射会使振幅谱以槽痕为特征,很难分清单个反射的振幅谱特征.由于在实际运用中,以傅里叶变换相关的算法的时窗问题,难以选择好时窗长度,而且无法定量分析时窗长度产生的偏差,因而会使振幅谱的估算产生偏差.以小波变换为基础的时频分析技术成了非平稳性信号的重要分析工具,在很多实际应用中已取代了傅里叶变换的分析方法.以小波变换为基础的瞬时谱分析技术能得到精确的时频分析结果,同时避免了时窗问题.它反映出了储层在纵向上时间及厚度上变化情况和横向上的地质不连续性的信息,因此能使解释人员快速而有效地描述储层特征的空间变化.     

ACTIVE system for monitoring volcanic activity: A case study of the Izu-Oshima Volcano,Central Japan

A system is proposed for the monitoring of changes in the underground structure of an active volcano over time by applying a transient electromagnetic method. The monitoring system is named ACTIVE, which stands for Array of Controlled Transient-electromagnetics for Imaging Volcano Edifice. The system consists of a transmitter dipole used to generate a controlled transient electromagnetic (EM) field and an array of receivers used to measure the vertical component of the transient magnetic field at various distances, with automatic operation of both units. In order to verify the performance of the proposed system, numerical and field experiments were carried out by application of the system to the Izu-Oshima volcano, where a remarkable change in the apparent DC resistivity over time had been detected in association with the eruption in 1986.     

玄武岩内非均质性对地震资料的影响:印度德干圈闭的一个研究实例(英文)

油气勘探需要重新开发玄武岩下目的层的成像技术。我们今天遇到的最重要的问题之一是玄武岩下地震成像。近年来,由于利用长炮检距这个问题部分已经得到解决。然而,由于地表以及玄武岩的内部非均质性引起的复杂波形,近炮检距的资料仍然不能充分地被利用。影响玄武岩下成像的近垂直入射资料对于了解玄武岩层内部结构是十分有用的。对比较均匀的玄武岩目标层可选用转换波。这里我们利用几个合成模型例子重点描述了更接近现实的非均质玄武岩流引起的实际困难。模拟计算了全波地震记录以有助于了解玄武岩内沉积物对地震资料的影响。本文介绍了印度德干圈闭的一个研究实例。首先探讨了夹层沉积物对整个地震成像的影响。其次利用该区声测井资料以反射系数法计算全波场响应与实际地震资料对比, 通过一系列速度-深度剖面探讨了利用模式转换波(顶底玄武岩界面上的P波转化为S波或相反)对印度库奇地区玄武岩下成像的可行性。通过与野外资料相比较证明玄武岩中多个薄层的效应严重恶化我们所要解释和利用的图象的质量。