Effect of 2D spatial variability on slope reliability: A simplified FORM analysis

To meet the high demand for reliability based design of slopes, we present in this paper a simplified HLRF(Hasofere Linde Rackwitze Fiessler) iterative algorithm for first-order reliability method(FORM). It is simply formulated in x-space and requires neither transformation of correlated random variables nor optimization tools. The solution can be easily improved by iteratively adjusting the step length. The algorithm is particularly useful to practicing engineers for geotechnical reliability analysis where standalone(deterministic) numerical packages are used. Based on the proposed algorithm and through direct perturbation analysis of random variables, we conducted a case study of earth slope reliability with complete consideration of soil uncertainty and spatial variability.     

Probabilistic stability analyses of undrained slopes by 3D random fields and finite element methods

A long slope consisting of spatially random soils is a common geographical feature. This paper examined the necessity of three-dimensional(3 D) analysis when dealing with slope with full randomness in soil properties. Although 3 D random finite element analysis can well reflect the spatial variability of soil properties, it is often time-consuming for probabilistic stability analysis. For this reason, we also examined the least advantageous(or most pessimistic) cross-section of the studied slope. The concept of"most pessimistic" refers to the minimal cross-sectional average of undrained shear strength. The selection of the most pessimistic section is achievable by simulating the undrained shear strength as a 3 D random field. Random finite element analysis results suggest that two-dimensional(2 D) plane strain analysis based the most pessimistic cross-section generally provides a more conservative result than the corresponding full 3 D analysis. The level of conservativeness is around 15% on average. This result may have engineering implications for slope design where computationally tractable 2 D analyses based on the procedure proposed in this study could ensure conservative results.     

Origin of natural sulfur-metal chimney in the Tangyin hydrothermal field,Okinawa Trough: constraints from rare earth element and sulfur isotopic compositions

For the first time, we present the rare earth element (REE) and sulfur isotopic composition of hydrothermal precipitates recovered from the Tangyin hydrothermal field (THF), Okinawa Trough at a water depth of 1206 m. The natural sulfur samples exhibit the lowest ΣREE concentrations (ΣREE= 0.65×10^–6–4.580×10^–6) followed by metal sulfides (ΣREE=1.71×10^–6–11.63×10^–6). By contrast, the natural sulfur-sediment samples have maximum ΣREE concentrations (ΣREE=11.54×10^–6–33.06×10^–6), significantly lower than those of the volcanic and sediment samples. Nevertheless, the δEu, δCe, (La/Yb)_N, La/Sm, (Gd/Yb)_N and normalized patterns of the natural sulfur and metal sulfide show the most similarity to the sediment. Most hydrothermal precipitate samples are characterized by enrichments of LREE (LREE/HREE=10.09–24.53) and slightly negative Eu anomalies or no anomaly (δEu=0.48–0.99), which are different from the hydrothermal fluid from sediment-free mid-oceanic ridges and back-arc basins, but identical to the sulfides from the Jade hydrothermal field. The lower temperature and more oxidizing conditions produced by the mixing between seawater and hydrothermal fluids further attenuate the leaching ability of hydrothermal fluid, inducing lower REE concentrations for natural sulfur compared with metal sulfide; meanwhile, the negative Eu anomaly is also weakened or almost absent. The sulfur isotopic compositions of the natural sulfur (δ³⁴S=3.20‰–5.01‰, mean 4.23‰) and metal sulfide samples (δ³⁴S=0.82‰–0.89‰, mean 0.85‰) reveal that the sulfur of the chimney is sourced from magmatic degassing.     

黔东地区航磁特征与岩性构造填图

为了弄清黔东深覆盖地区的成矿地质背景,以2014年该区1:5万比例尺的直升机航磁资料为基础,结合地质、重力等资料,研究了该区不同岩浆岩、断裂构造的航磁特征,在综合分析的基础上,编制了该区的岩性构造图。通过分析对比,本次新圈定具有一定规模的隐伏岩体7处、蚀变岩15处,新圈定不同规模的断裂17条 。上述成果对黔东地区基础地质研究和找矿工作具有重要意义。     

重力方法在祁连山地区构造特征与岩体解释中的应用

依据地质、航磁等资料,结合实测岩石物性资料,以1∶250 000重力资料为基础对祁连山地区(张掖、刚察地区)进行综合研究。研究地质构造与地球物理特征,划分了大地构造单元,探讨了隐伏与半隐伏岩体分布。该区域经多年的重力工作,勘探面积已近3万多km²,根据重力资料研究布格重力场特征并进行分区; 对地质构造单元重新划分,对主要断裂进行探讨; 圈定隐伏与半隐伏岩体,推断其空间展布特征,并对其解释。结果表明: 祁连山地区地壳密度结构横向分区、纵向分层明显,布格重力异常场的分区特征与地质构造分区特征基本一致。结合2.5D重力剖面反演印证了新的推测结果,为综合地质、地球物理找矿提供参考。     

改进的无网格法及其在冻土区桩基温度场中的应用

采用带权正交基函数对传统的无网格法中的基函数进行了改进,避免了计算过程中可能出现的矩阵不可逆情况,而且编程容易实现且计算效率高。并将其应用到冻土区桩基包括了热传导和相变潜热的温度场中,同时考虑了混凝土水化热释放对桩周冻土的影响,将计算结果和有限元计算结果及现场实测数据进行了比较分析,从趋势和最大数值看都反映了实测曲线的趋势,又把不同深度处桩侧和不同桩径处温度随时间的变化规律计算结果和有限元计算结果做了对比,都验证了该方法的可行性和优越性。     

辽宁赛马碱性岩体异性石化学成分特征及其蚀变组合对碱性岩浆-热液演化的指示意义

辽宁赛马碱性岩体早年因产铀矿而闻名,该岩体主要由响岩、霞石正长岩和异霞正长岩组成,其中铀、锆和稀土等元素矿化主要集中于异霞正长岩岩浆阶段。异性石是异霞正长岩中特征的锆-稀土矿物,主要分为两期,晚期异性石表现出更加富集Nb、REE等高场强元素的特点。早期异性石经历了一系列的热液蚀变,根据蚀变强弱程度,蚀变矿物组合可分为:(1)异性石+钠锆石+霓石±钠沸石;(2)异性石+钠锆石+锆石+钠沸石±霓石;(3)异性石假晶,假晶主要由残余异性石+钠锆石+锆石+钠沸石+霓石+钾长石+铈硅磷灰石组成。相比于岩浆锆石,蚀变组合中次生锆石具有富Ca、Al、Fe的特点,与异性石本身化学成分和流体性质密切相关。通过对异性石及其蚀变组合的精细矿物学研究,我们得知假晶的形成可能是异性石"溶解-再沉淀"的结果,致使假晶形成的流体至少包括:(1)占主导的富Na(±K)、Al、F的自交代流体;(2)少量晚期富Ca流体。假晶中次生锆石和铈硅磷灰石的结晶说明了Zr和REE等高场强元素的热液活动性,自交代碱性流体和富Ca流体在此过程中起到"搬运"和"提纯"的作用,这对认识碱性岩稀有、稀土成矿机制具有重要的指示意义。     

基于重力及重力梯度数据的位场偏移成像

为了验证位场偏移快速成像方法的有效性以及对比基于重力和重力梯度数据的偏移成像的效果,笔者利用MATLAB编写了重力及重力梯度数据位场偏移成像的算法,并进行了模型正演计算,发现基于梯度数据的偏移成像的横向分辨率、抗噪性和抗背景场能力好于重力数据的偏移成像,并且能够较好的确定模型的中心埋深。将此方法应用于葫芦岛市黑鱼沟的某铅锌矿洞的深度计算中,验证了该方法的有效性。     

基于改进K-SVD字典学习方法的地震数据去噪

为实现更好的地震数据去噪技术,笔者引入一种新的算法:快速迭代收缩阀值法(FISTA),通过FISTA和K-奇异值分解(K-SVD)不断迭代更新K-SVD字典,利用更新得到的K-SVD字典对地震数据进行稀疏表示,去除稀疏系数中较小的数值,使数据中的随机噪声得到压制。对层状模型合成地震记录,Marmousi模型合成地震记录以及实际地震数据进行对比实验,得出FISTA算法较OMP算法能更好地提高地震数据的信噪比,同时有效地保护了反射信号。     

2017年九寨沟7.0级地震序列震源机制解和构造应力场特征

九寨沟余震序列的震源机制和构造应力场有助于认识本次地震的发震构造和孕震机理。本文基于四川区域地震台网的波形资料, 采用波形拟合(CAP)方法和P波初动+振幅比(HASH)方法反演得到2017年8月8日九寨沟7.0级地震序列中59次ML≥3.0地震的震源机制解, 并基于该结果采用阻尼线性逆推法(DRSSI), 计算研究区域的平均构造应力场, 给出该区域的应力场特征。结果显示, 利用CAP方法反演得到的本次主震的最佳双力偶机制解节面I: 走向248°/倾角86°/滑动角–169°, 节面II: 走向157°/倾角79°/滑动角–4°, 矩震级为Mw6.31, 矩心深度5 km, 属走滑型地震事件; 大部分余震的震源机制解错动类型与主震一致, 矩心深度集中在3~10 km; 应力场反演结果显示, 该区域周边的应力性质为走滑型, 最大主应力方向呈NWW–SEE向, 与该区域的应力场方向一致, 表明本次地震主要受区域应力的控制。结合该区域的地震地质构造等已有研究成果, 分析认为此次地震的发震断层为走向NW–SE、倾向SW的左旋走滑断裂——树正断裂, 巴颜喀拉块体向E-SE向的水平运动受到华南块体的强烈阻挡导致此次地震的发生, 汶川地震的发生对本次地震具有一定的促进作用。