一种适用于各向异性土体的破坏准则

自然界的土体通常具有各向异性的特点,而传统的破坏准则大多只适用于各向同性的土体。结合应力张量和反映材料各向异性状态的组构张量,定义了修正偏应力及其不变量,提出了适用于各向异性土体材料的破坏准则。给出了共轴条件下正交各向异性和横向各向异性材料在一般应力空间的破坏曲线以及不同应力区中主应力系数b与摩擦角的关系曲线,并分析了它们的特性以及与各向同性材料相应曲线的区别。通过与真三轴试验数据的比较,表明该准则能很好地描述各向异性土体材料的强度特点。     

Study on elastic response of structures to near-fault ground motions through record decomposition

Accelerograms recorded near active faults have some important characteristics that make them different from those recorded in far-fault regions. High-frequency components in acceleration records and long-period velocity pulses are among notable specifications of such ground motions. In this paper, a moving average filtering with appropriate cut-off frequency has been used to decompose the near-fault ground motions into two components having different frequency contents: first, Pulse-Type Record (PTR) that possesses long-period pulses; second, the relatively high-frequency BackGround Record (BGR), which does not include large velocity pulses. Comparing the results with those extracted through wavelet analysis shows that moving average filter is an appropriate and efficient tool for near-fault records decomposition. The method is applied to decompose a suite of 91 selected near-fault records and the elastic response of structures is examined through their response to the decomposed parts. The results emphasizes that in contrast with ordinary far-fault earthquake records, response spectra of near-fault ground motions typically have two distinct local peaks, which are representatives of the high- and low-frequency components, i.e., BGR and PTR, respectively. Moreover, a threshold period is identified below which the response of structures is dominated by BGR while PTR controls the response of structures with periods longer than this period.     

地形对地震烈度衰减的影响

本文采用阶跃截面固体的一维弹性波传播模型,研究了地形对地震烈度衰减的影响。采用固体入射纵波来模拟无地形影响的地震烈度衰减;采用固体透射波来模拟有地形影响的地震烈度衰减。利用阶跃截面处的放松连续条件,得出了烈度与高度差的近似定量关系。同时,以新疆地区逆冲走滑型地震为例,讨论了烈度与高度差的近似定量关系描述地形对地震烈度衰减影响的有效性。     

A new series solution method for two-dimensional elastic wave scattering along a canyon in half-space

This paper presents a semi-analytical method for studying the two-dimensional problem of elastic wave scattering by surface irregularities in a half-space. The new method makes use of the member of a c-completeness family of wave functions to construct the scattering fields, and then applies equal but opposite tractions to those of the foregoing constructed scattering fields on the horizontal surface of the half-space to produce additional scattering fields. These additional scattering fields are a series of Lamb's solutions. Thus the whole scattering field constructed in the series automatically satisfies the Navier equations, the condition of zero traction on the half-space surface, and the radiation boundary conditions at infinity. Using the traction-free conditions along the canyon surface, the coefficients of the series solutions are determined via a least-squares method. For incident P, SV, and Rayleigh waves, the numerical results are presented for the scattering displacements in the vicinity of a semi-circular canyon in the half-space.     

Shear-tensile crack as a tool for reliable estimates of the non-double-couple mechanism: West Bohemia-Vogtland earthquake 1997 swarm

Shear-tensile crack is a model for an earthquake mechanism that is more constrained than the moment tensor but that can still describe a non-shear focus. As such, the shear-tensile crack model is more robust than the moment tensor model and yields more reliable estimates for the earthquake mechanism. Such an advantage verifies the credibility of the non-double-couple component found for some events of the 1997 West Bohemia-Vogtland earthquake swarm. As expected, in several cases, a significantly resolved non-double-couple component was obtained where the moment tensor approach failed. Additionally, for non-shear sources, the shear-tensile crack model offers optimization of the Poisson number within the focus, concurrently with retrieval of the mechanism. However, results obtained for the joint inversion of the 1997 swarm indicate that resolution is low. A series of synthetic experiments indicated that limited observations during 1997 were not the cause. Rather, hypothetical experiments of both very good and extremely poor network configurations similarly yielded a low resolution for the Poisson number. Applying this method to data for recent swarms is irrelevant because the small non-double-couple components detected within the inversion are spurious and, thus, the events are pure double-couple phenomena.     

Source parameters of small and moderate earthquakes in the area of the 2009 L’Aquila earthquake sequence (central Italy)

The main goal of this study is to provide moment tensor solutions for small and moderate earthquakes of the 2009 L’Aquila seismic sequence (central Italy). The analysis was performed by using data coming from the permanent Italian seismic network run by the Istituto Nazionale di Geofisica e Vulcanologia (INGV) and the “Cut And Paste” (CAP) method based on broadband waveform inversion. Focal mechanisms, source depths and moment magnitudes are determined through a grid search technique. By allowing time shifts between synthetics and observed data the CAP method reduces dependence of the solution on the assumed velocity model and on earthquake location. We computed seismic moment tensors for 312 earthquakes having local magnitude in the range between 2.7 and 5.9. The CAP method has made possible to considerably expand the database of focal mechanisms from waveform analysis in the lowest magnitude range (i.e. in the neighborhood of magnitude 3) without overlooking the reliability of results. The obtained focal mechanisms generally show NW–SE striking focal planes in agreement with mapped faults in the region. Comparisons with the already published solutions and with seismological and geological information available allowed us to proper interpret the moment tensor solutions in the frame of the seismic sequence evolution and also to furnish additional information about less energetic seismic phases. Focal data were inverted to obtain the seismogenic stress in the study area. Results are compatible with the major tectonic domain. We also obtained a relation between moment and local magnitude suitable for the area and for the available magnitude range.     

三维/三维构造下大地电磁相位张量数值模拟

当地表存在三维非均匀电导率分布时,区域大地电磁响应发生畸变. 以往对这种畸变研究多假设近地表为三维,区域构造为一维或二维. 对于更一般的三维/三维构造,为了分析并消除这种畸变影响,真实反映地下三维区域构造信息,本文实现了三维大地电磁相位张量积分方程数值算法,并研究在不同地质模型下相位张量响应. 结果表明,相位张量不仅可以反映一般三维构造信息,亦可有效反映复杂近地表构造下三维区域构造信息,而无须假设区域构造为一维或二维,证明相位张量具有较强抗近地表局部非均匀构造干扰能力,能够保持更为一般的三维区域构造信息. 为了加快正演计算,同时保持一定精度,算法采用了积分方程多网格法.     

结构保持的图像序列自适应超分辨率重建

基于概率运动估计的超分辨算法能克服传统超分辨率方法过分依赖于准确运动估计的缺点,为了改善其重建效果,通过分析图像的局部结构模式,给出了局部空域活跃度的定义表征图像的局部结构特征,设计了形状、尺度自适应的非局部空时邻域窗口函数,进而提出了一种结构保持的超分辨率自适应重建算法。实验结果显示,本文方法具有更好的重建性能。     

2009年11月至2010年11月27次中强地震的快速矩张量解

经过半个世纪的发展,地震震源的矩张量表示及其反演已日臻成熟(Dziewonski,Woodhouse,1983).矩张量解不但能够较好地反映地震震源的特性,还能为其它研究提供有用的参考.一些科研机构和组织,如美国地质调查局(USGS)和全球矩心矩张量研究组(GCMT)已经能够在震后5小时内准实时地发布矩     

大地电磁测深数据处理——4元素和6元素阻抗张量元素的比较研究(英文)

阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期(超过100秒),逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。