Two-dimensional simulations of surface deformation caused by slab detachment

Detachment of the deeper part of subducted lithosphere causes changes in a subduction zone system which may be observed on the Earth's surface. Constraints on the expected magnitudes of these surface effects can aid in the interpretation of geological observations near convergent plate margins where detachment is expected. In this study, we quantify surface deformation caused by detachment of subducted lithosphere. We determine the range of displacement magnitudes which can be associated with slab detachment using numerical models. The lithospheric plates in our models have an effective elastic thickness, which provides an upper bound for rapid processes, like slab detachment, to the surface deformation of lithosphere with a more realistic rheology. The surface topography which develops during subduction is compared with the topography shortly after detachment is imposed. Subduction with a non-migrating trench system followed by detachment leads to a maximum surface uplift of 2–6 km, while this may be higher for the case of roll-back preceding detachment. In the latter situation, the back-arc basin may experience a phase of compression after detachment. Within the context of our elastic model, the surface uplift resulting from slab detachment is sensitive to the depth of detachment, a change in friction on the subduction fault during detachment and viscous stresses generated by sinking of the detached part of the slab. Overall, surface uplift of these magnitudes is not diagnostic of slab detachment since variations during ongoing subduction may result in similar vertical surface displacements.     

Seismic landslide simulations in discontinuous deformation analysis

Discontinuous deformation analysis (DDA) is a numerical approach used to simulate the post-failure behavior of a blocky assembly. Three available algorithms incorporate seismic impacts into DDA simulations for earthquake-induced slope failure. The following methods are used: directly applying time-dependent accelerations to falling/sliding blocks (Method 1); adding time-dependent accelerations to base block (Method 2); and time-dependently constraining seismic displacements of the base block (Method 3). However, incorrect absolute movements of falling/sliding and base blocks were obtained using Method 1. Additionally, relative movements between falling/sliding blocks and the base block are opposite to those simulated by the other two algorithms—Methods 2 and 3. Since locating an earthquake-induced landslide before an earthquake is extremely difficult, the seismic movements of base rock are recorded. Method 1 applies recorded seismic data to sliding blocks in conflict with d’Alembert’s principle of mechanics. Additionally, in Method 2, when the computation time step must be longer than the time in seismic data, computational results reveal abnormal base block displacements due to the non-zero velocity recorded at the end time of seismic data in seismic DDA. In this study, a novel algorithm to diminish the velocity of the base rock in the seismic analysis is utilized to modify Method 2. Furthermore, this work confirms that DDA with the modified Methods 2 and 3 is a practical approach for earthquake-induced landslide simulations.     

砂土液化流动变形的简化方法

已有的液化砂土流动特性试验结果表明,砂土在液化流动状态下是剪切稀化非牛顿流体,可以用幂函数表示其剪应力-剪应变率的关系,从而建立了砂土液化流动的本构方程。基于FLAC3D程序的二次开发平台,将液化流动本构方程开发到FLAC3D中,建立了液化流动变形的简化分析方法。通过倾斜场地的液化流动变形分析,发现倾斜场地的液化变形曲线可以用正弦函数曲线描述,这与Towhata的理论分析成果一致,验证了本方法的合理性。分析了液化层坡度、稠度系数、流动指数以及弹性参数等变量对液化变形的影响。计算结果表明,液化变形随液化层坡度的增大而逐渐增大,液化砂土的稠度系数和流动指数对液化流动变形有重要的影响,而弹性参数对变形基本无影响,因此,在实际工程分析中,需要对流动模型参数进行深入研究。     

Estimating inelastic sediment deformation from local site response simulations

Significant insight into the dynamic local site response of a horizontally layered sediment deposit to seismic excitation can be gained from numerical simulations. In this paper we use a nonlinear local site response analysis code SPECTRA to estimate the coseismic sediment deformation at a seismically active site in Lotung, Taiwan. We address some basic issues relevant for interpreting the simulation results, including the impact of noise and baseline offsets present in the input ground motion. We also consider the sensitivity of the predicted deformation responses to statistical variations of sediment constitutive properties. Finally, we apply a suite of hypothetical strong ground motions to the base of the sediment deposit to better understand the pattern of inelastic deformation likely to result from strong seismic shaking.     

Rotation of porphyroblasts in non-coaxial deformation: insights from computer simulations

Porphyroblast inclusion trails have the potential to provide critical information about tectonometamorphic events. Recently, however, traditional interpretations of inclusion trails have been called into question by the suggestions that porphyroblasts do not rotate during non-coaxial deformation and that apparent spiral inclusion trails can be generated in coaxial deformation. We present a new computer model that simulates inclusion trail development. Model results suggest: (1) that the extent of porphyroblast rotation is controlled by conditions at the porphyroblast-matrix boundary; (2) that curved inclusion trails may develop in unrotated porphyroblasts; (3) that classic "snowball" inclusion trails are most simply explained by rotational growth histories; and (4) that some of the observations used to support the view that porphyroblasts do not rotate (e.g. weakly sigmoidal inclusion trails, apparent truncations of inclusion trails) can be accounted for by variations in the growth rate of rotating porphyroblasts.     

Discrete element method simulations of bio-cemented sands

Microbially induced calcite precipitation (MICP) has emerged as a novel soil improvement method. In this paper, 3-D discrete element method (DEM) simulations are used to explore the behavior of MICP-cemented sands. Comparisons of the macro-scale response of numerical and physical specimens are made. Microstructure analyses indicate a shear band formed in the numerical specimens, consistent with physical experiments. The bond breakage pattern in numerical specimens is explored and compared to observed measurements from physical specimens. The relationship between dilatancy and stress-strain behavior is evaluated. The results indicate DEM is an effective technique to capture the mechanical behavior of MICP-cemented sand.     

应用RMi法估算岩体变形模量

RMR法与Q法是国内外八、九十年代岩体质量和参数估算的常用方法,然而RMR法和Q法对质量较差的岩体不太实用,而RMi法是一种既适用于软岩又适用于硬岩的岩体质量评价和参数估算方法。本文在介绍RMi法基本原理的基础上,对贵州省鱼简河水库坝基岩体变形模量进行确定。     

Innovative passive reinforcements for the gradual stabilization of a landslide according with the observational method

A large number of landslides occur in North-Eastern Italy during every rainy period due to the particular hydrogeological conditions of this area. Even if there are no casualties, the economic losses are often significant, and municipalities frequently do not have sufficient financial resources to repair the damage and stabilize all the unstable slopes. In this regard, the research for more economically sustainable solutions is a crucial challenge. Floating composite anchors are an innovative and low-cost technique set up for slope stabilization: it consists in the use of passive sub-horizontal reinforcements, obtained by coupling a traditional self-drilling bar with some tendons cemented inside it. This work concerns the application of this technique according to the observational method described within the Italian and European technical codes and mainly recommended for the design of geotechnical works, especially when performed in highly uncertain site conditions. The observational method prescribes designing an intervention and, at the same time, using a monitoring system in order to correct and adapt the project during realization of the works on the basis of new data acquired while on site. The case study is the landslide of Cischele, a medium landslide which occurred in 2010 after an exceptional heavy rainy period. In 2015, some floating composite anchors were installed to slow down the movement, even if, due to a limited budget, they were not enough to ensure the complete stabilization of the slope. Thanks to a monitoring system installed in the meantime, it is now possible to have a comparison between the site conditions before and after the intervention. This allows the evaluation of benefits achieved with the reinforcements and, at the same time, the assessment of additional improvements. Two stabilization scenarios are studied through an FE model: the first includes the stabilization system built in 2015, while the second evaluates a new solution proposed to further increase the slope stability.

     

HyMap-C机载高光谱仪定标方法与结果评价

HyMap-C机载高光谱仪能够在400~2 500 nm波长范围内获取144个精细光谱通道,对于地球表面地物特征有极强的探测能力.为保证高光谱遥感数据的可靠性与定量化分析的精度,作者分别采用单色准直光法、漫反射板定标法和内部参考标准法对HyMap-C机载高光谱仪进行了实验室光谱定标、绝对辐射定标以及飞行过程中航线上相对辐射定标.通过航空飞行测量与地面实际测量的光谱反射率数据进行对比验证,结果表明两者整体谱型与局部吸收特征非常一致,定标结果可靠,数据可信.文中阐述的定标原理与方法为我国自主研制机载高光谱仪的定标提供了依据与思路,定标后获取到的高光谱分辨率数据可广泛应用于地质调查领域特别是矿物填图中.     

Virtual element method for geomechanical simulations of reservoir models

In this paper, we study the use of virtual element method (VEM) for geomechanics. Our emphasis is on applications to reservoir simulations. The physical processes behind the formation of the reservoirs, such as sedimentation, erosion, and faulting, lead to complex geometrical structures. A minimal representation, with respect to the physical parameters of the system, then naturally leads to general polyhedral grids. Numerical methods which can directly handle this representation will be highly favorable, in particular in the setting of advanced work-flows. The virtual element method is a promising candidate to solve the linear elasticity equations on such models. In this paper, we investigate some of the limits of the VEM method when used on reservoir models. First, we demonstrate that care must be taken to make the method robust for highly elongated cells, which is common in these applications, and show the importance of calculating forces in terms of traction on the boundary of the elements for elongated distorted cells. Second, we study the effect of triangulations on the surfaces of curved faces, which also naturally occur in subsurface models. We also demonstrate how a more stable stabilization term for reservoir application can be derived.     

冻结黏土三轴试验微观变形机理的离散元分析

基于三维颗粒离散单元法,赋予颗粒相应的细观参数,并采用黏结发生在接触颗粒间有限范围内的模型来考虑冻土颗粒中冰的胶结作用,建立了冻结黏土三维离散元数值模型.在相同围压、不同温度和相同温度、不同围压下对冻结黏土的室内三轴试验进行数值模拟,对比了数值试验与室内测试的应力-应变曲线,两者吻合较好.数值模拟结果表明:围压增大会使得接触黏结逐渐失效,在剪切带中胶结冰的破坏区域将增大,而温度的降低则会产生相反结果,这些微观变化都将对冻结黏土的宏观力学变形产生较大影响,同时,细观参数对温度的依赖性也很明显.冻结黏土三轴试验微观变形离散元模拟思路及方法可为今后运用离散单元法研究冻土力学行为提供一定的参考.     

盾构法施工隧道地层变形时空统一预测方法研究

从经典的Mindlin弹性理论公式出发,基于盾构施工工作面推力为圆形均布荷载、盾构外壁与土层的摩擦力均匀分布等假定的基础上,考虑盾构施工推进过程中盾构机位置不断变化的实际情况,建立了盾构施工推进过程的力学计算模型,由此推导了盾构推进工作面附加压力以及盾构机外壁与土层摩擦力作用下隧道周边任意位置地层变形的计算公式。基于随机介质理论,充分考虑隧道开挖施工引起的地层变形的时间和空间发展过程,推导了盾构推进地层损失引起的地层变形计算公式。通过工程实例计算,对盾构隧道施工地层变形的时空变化过程进行了研究,计算结果同实测结果具有良好的一致性。     

土石料湿化变形计算方法探讨

土石料的湿化是一种普遍现象,虽然可以利用非饱和理论计算湿化变形,但人们更愿意采用常规的初应变法进行计算,湿化试验则常用单线法或双线法进行。对湿化试验和相关计算方法进行分析,提出考虑湿化变形及其影响的新思路,（1）在某一应力状态下的湿化变形是由两部分组成的,一部分是由于材料本身刚度变“软”引起的湿化变形,另一部分变形是湿化引起的附加变形;（2）在双线法试验中直接将相同应力状态下对应干土和湿土的应变相减计算是不合适的,应当采用相同应力状态下各自的切线模量计算;（3）基于广义位势理论,根据应力和应变的对偶特性,可以把湿化附加变形变换为一种湿化等效应力,虽然其是一种近似处理方法,但它比较简单,容易直观地量化湿化的作用。文中思路对于颗粒破碎、干湿循环引起的变形、流变变形等也具有借鉴意义。     

边坡大变形损伤破坏的高效无网格模拟

首先讨论了Lagrangian和Eulerian无网格近似的联系和区别,然后基于稳定节点积分和增量本构理论,建立了分析边坡静动力破坏的高效大变形无网格法,并给出了详细的计算流程。该方法采用弹塑性损伤耦合本构关系来模拟岩土类材料的破坏演化过程,其中屈服函数采用Drucker-Prager准则,损伤准则为基于应变的各项同性损伤函数。由于无网格近似和稳定节点积分具有非局部近似的特性,在保证空间离散稳定性和提高计算效率的同时,也可准确有效地模拟应变集中所形成的剪切带的发生与扩展,通过数值算例验证了方法的有效性。     

沉降曲线预测方法-变形过程指数法

高速铁路桥梁沉降控制是保证施工质量的重要措施.通过对京沪高速铁路桥梁大量沉降观测资料的积累,对软土地基沉降过程具有一定实际应用价值的沉降变形规律进行曲线回归,以预测其沉降发展规律.为了对桥梁基础沉降曲线全过程进行拟合和预测,通过各种理论计算和大量沉降观测资料的分析、统计归纳,首次提出了变形过程指数法.该成果为适用于京沪高速铁路及其他相关工程的桥涵基础的沉降预测提拱参考依据.     

非连续变形分析方法中的阻尼问题研究

定量研究了非连续变形分析(DDA)方法中的黏性阻尼与数值阻尼。首先,基于Newmark直接积分法,推导了块体系统的运动方程。其次,通过动力学中的黏性阻尼理论建立了DDA中动力系数、时间步长与黏性阻尼比的关系式,探讨了DDA中的常加速度积分方案的数值阻尼分区及阻尼比计算方法,进而得到两种阻尼共同作用时的阻尼比表达式,并分析了频域内阻尼比的分布情况。最后,以谐振激励下的块体振动为例,通过对比不同阻尼作用下块体位移的DDA计算值与理论解,验证了本文提出的阻尼比计算公式的正确性。研究表明：黏性阻尼对低频的衰减作用明显,数值阻尼则可以很快地消除高频干扰,而二者共同作用下可降低阻尼的频率相关性。该研究成果为DDA的振动、波动等动力计算的阻尼取值提供了理论依据。     

Natrolite,part II: Determination of deformation electron densities by the X-X method

A single crystal of natrolite, Na₂Al₂Si₃O₁₀ ·2H₂O (space group Fdd2), was studied by X-ray diffraction methods at room temperature. The intensities were measured in a complete sphere of reflection up to sinΘ/ λ=0.903 Å^?1. A refinement of high-order diffraction data yielded residuals of R/(F)=0.9%, R_w(F)=0.8%, GoF=1.40 for 1856 high-angle reflections (0.7≤sinΘ/ λ≤0.903 Å^?1) and R(F)=1.0%, R_w(F)=1.2%, GoF=3.07 for all 3471 independent reflections in the complete sphere of reflection. The X-X method was used to calculate deformation electron densities (DED) in natrolite. Within all tetrahedra, residual electron density-was found in the T-O bond directions indicating a considerable covalent contribution to the chemical bond. The range of the interatomic peak heights was from 0.19 to 0.34 e/ų in the SiO₄ tetrahedra and from 0.11 to 0.23 e/ų in the AlO₄ tetrahedron. The ionic contribution to the chemical bond manifests itself in the displacement of the peaks towards the oxygen atoms. Charge displacement due to interaction of nonframework cations with framework oxygen atoms as well as electron densities attributable to the lone pair orbitals in the water molecule have been observed.     

机械预切槽法开挖软土隧道地层变形研究

软土中的机械预切槽法是超前支护的一种，在预筑拱的保护下可以安全高效地开挖隧道。在隧道力学和弹塑性理论的基础上，利用自行开发的三维程序，分析了机械预切槽法与常规全断面法开挖隧道时引起的地层变形。计算结果显示：机械预切槽法形成的预支护作用在控制地层变形方面较常规全断面法具有明显的优势，预切槽法可以减少隧道施工对地层变形的影响范围和影响程度，可以用于对地表沉降控制较严格地区的隧道施工。     

Numerical simulations for coupled rock deformation and gas leak flow in parallel coal seams

Based on the new viewpoint of interaction mechanics for solid and gas, gas leakage in parallel deformable coal seams can be understood. That is, under the action of varied geophysical fields, the methane gas flow in a double deformable coal seam can be essentially considered to be compressible with time-dependent and mixed permeation and diffusion through a pore-cleat deformable, heterogeneous and anisotropic medium. From this new viewpoint, coupled mathematical models for coal seam deformation and gas leak flow in parallel coal seams were formulated and the numerical simulations for slow gas emission from the parallel coal seams are presented. It is found that coupled models might be close to reality. Meanwhile, a coupled model for solid deformation and gas leak flow can be applied to the problems of gas leak flow including mining engineering, gas drainage engineering and mining safety engineering in particular the prediction of the safe range using protective layer mining where coal and gas outbursts can efficiently be prevented. This revised version was published online in July 2006 with corrections to the Cover Date.