北大别穹隆变形组构一期或多期形成:来自韧性变形组构叠加数值模拟的证据

矿物拉伸线理和面理是研究岩石韧性变形运动学的基本要素。然而自然界许多应变带内拉伸线理和面理产状分布往往并不单一。如何正确认识这类变形组构,并利用这些组构开展构造变形分析至今仍备受争议。本文基于Eshelby理论,通过数值模拟,研究了岩石经历两期不同韧性变形后面理和拉伸线理的产状特征。模拟结果显示岩石经历两期构造变形时,岩石中的较软矿物仅能保留下最后一期构造变形运动学信息。这一认识得到了北大别穹隆内变形组构的支持。北大别穹隆内产状复杂的拉伸线理只记录了造山后伸展期变形。本次数值模拟结果对野外构造观察和分析以及确定构造变形时间具有指导意义。     

一种模拟韧性剪切带中椭球状标志体运动、变形的新方法

韧性剪切带中赋存的椭球状应变标志体(砾石、碎斑等)是研究剪切带运动学、动力学的重要应变标志体。传统研究中椭球状标志体通过对野外露头的观测或与实验岩石学结合,判断剪切带运动方向、探讨运动学和流变学特征。随着数值模拟技术在韧性剪切带中的引入和推广,国内外许多学者试图恢复椭球状标志体的运动轨迹和变形特征,并取得了显著的成果。然而,国内文献对于模拟韧性剪切带椭球状标志体的定量及模拟研究甚少,研究方法也鲜为介绍。基于此,针对韧性剪切带中椭球状标志体变形的最新研究进展,详细介绍建立在Jeffery理论和Eshelby理论之上的数值模拟思路和方法,并利用Mathcad软件模拟了给定条件下的椭球状标志体的运动轨迹、变形特征。     

渗流场中岩石流变特性的数值模拟

渗流对岩石的流变特性有显著的影响,但限于试验条件,相关的试验成果较少。而数值计算技术已成为进行岩石力学研究的重要方法,因此,采用数值模拟方法对渗流场中岩石的流变特性进行了研究。对已有的岩石室内流变试验进行模拟来验证数值方法的有效性,在压缩流变试验的数值模拟中对试件分别施加单向、三向渗流场来模拟不同渗流条件对岩石流变特性的影响。模拟结果表明,在渗流方向上渗流产生的动水压力会导致试件的流变变形增大。此外,不同渗流场对流变的影响也有所差异,三向渗流场中试件受到的轴向渗透力小于单向渗流场,而侧向渗透力则大于单向渗流情况,这也导致三向渗流场中试件的轴向应变小于单向渗流场,而侧向应变大于单向渗流场。     

岩石破坏的机理再认识

岩石的各种强度准则一直是解释岩石或岩体破坏机理的主要理论。此外,从流变现象出发,通过数值模拟、实验对照、理论分析、进一步对岩石破坏机理的解释进行了探讨,认为受力岩石或岩体内差异流动变形也是导致岩石破坏的原因。     

白银厂矿区劈理成因的磁组构分析

将磁性组构方法用于甘肃白银厂矿区的构造分析.结果表明该区岩石的磁性组构与岩石组构具有同一分布规律.岩石的磁线理与岩石的拉长线理L~a近于一致,磁面理与第一期劈理S_1基本平行,并垂直于岩石的磁化率椭球的最小主轴,也垂直于应变椭球的最短主轴方位,故推断第一期劈理S_1的成因主要属压扁成因.同时受到后期的多次变形与褶皱迭加影响.S_1轴面的优选面状分布发生变位不大,表明该区多次构造变形主要来自早期的同一主应力场作用所致.矿区岩石的磁面理与应变椭球拉长轴面或第一期劈理面S_1存在某些的角度差,表明该区岩石的劈理面由于受后期多次构造变形影响,可能发生剪切运动或旋转机制.     

海底隧道风化花岗岩流变试验研究

采用全自动三轴流变试验机对厦门东通道海底隧道风化槽地段岩石进行三轴压缩流变试验,研究了岩石在不同围压和不同应力水平作用下轴向应变随时间的变化规律。试验表明,强风化花岗岩时效特性较全风化花岗岩更加明显;围压对岩石蠕变变形存在很大影响,围压越大,蠕变变形量越小。通过分析岩石压缩流变过程中应力-应变关系可知,蠕变变形在弹性阶段不对岩石整体构成明显损伤。而进入塑性阶段后,黏塑性变形对岩石破坏影响较大。建议通过支护增加围压,以提高隧道围岩屈服强度和减小流变变形,防止隧道失稳破坏。     

粒度对韧性剪切带岩石变形的影响

构造模拟研究表明, (1) 粒度是控制岩石流变性质的重要因素之一, 粒度越粗, 岩石粘度越高; (2) 粒状矿物含量与基质的比值越低, 粒度对岩石流变性质的影响就越小, 反之则越大; (3) 岩石的粒度差异对韧性剪切带的发育程度有着重要影响, 粒度越细, 韧性剪切带中有限应变越大.这些结果与野外所观察的地质现象及应变测量结果完全一致.      

海底泥流的流变试验及强度模型

海底泥流是海底斜坡发生失稳滑动后,经复杂的水土交换作用演变成的稀式海底滑坡体,表现出土体和流体的双重特征,而现有的测试方法难以获得低剪切应变率下连续变化的强度值,不能很好地揭示其综合强度特性。利用新型全流动贯入仪和RST流变仪,对模拟海底泥流在不同剪切应变率下的流变和强度特性进行了多组试验研究,并基于试验结果分析了不排水剪切强度、屈服应力和表观黏度与含水率的相关性。基于剪切稀化理论,提出分阶段拟合模式来描述海底泥流从低剪切应变率到中高剪切应变率的流变关系;通过多种常规流变模式拟合结果的对比分析,显示出新流变模式的适用性和优势。考虑强度软化的影响,建立了全剪切应变率范围内海底泥流的不排水剪切强度模型,为海底泥流运动过程的数值模拟和灾害评价提供科学依据。     

脆性岩石破裂演化过程的三维细胞自动机模拟

在三维条件下定义了实体细胞自动机的基本组件,综合运用弹塑性理论、细胞自动机自组织演化理论、统计原理以及岩石力学理论等,建立了模拟岩石三维破裂过程的张量型细胞自动机模型,并开发了相应的数值模拟软件EPCA3D。利用EPCA3D对脆性岩石进行单轴压缩破裂过程模拟,通过破裂模式、岩样内部破裂情况、全应力-应变曲线和声发射曲线等的分析,揭示岩石破裂的机制。结果表明,EPCA3D能够较好地模拟非均质岩石在载荷作用下微裂纹的萌生、扩展和贯通的全过程。利用该模拟系统,采用应力-应变线性组合的加载控制方式,研究了不同围压对岩石全应力-应变曲线的影响。研究结果表明：对于脆性非均质岩石,低围压下容易表现为II类行为,且强度较低;高围压下容易表现为I类变形行为,且强度较高。     

模拟岩石应变软化变形全过程的统计损伤本构模型研究

建立可以模拟岩石应力-应变全过程的本构模型一直是岩土工程中研究的热点之一.首先,在深入探讨应变软化类岩石受力变形全过程每个阶段特征的基础上,并针对Lemaitre应变等价性理论的不足之处,基于更合理的假定建立了能够反映应变软化类岩石体积变化以及峰后出现残余强度的新型损伤模型;然后,通过引入统计损伤理论并考虑损伤阈值对岩石变形过程的影响,建立了能够模拟应变软化类岩石受力变形全过程的统计损伤本构模型,该模型还可以反映岩石的体积变化以及初始孔隙率对其变形过程的影响;最后,通过理论计算与试验结果及他人模型的对比分析,表明了该模型的合理性与可行性.     

The motion of deformable ellipsoids in power-law viscous materials: Formulation and numerical implementation of a micromechanical approach applicable to flow partitioning and heterogeneous deformation in Earth’s lithosphere

Earth’s lithosphere is heterogeneous in rheology on a wide range of observation scales. When subjected to a tectonic deformation, the incurred flow field can vary significantly from one rheologically distinct element to another and the flow field in an individual element is generally different from the bulk averaged flow field. Kinematic and mechanical models for high-strain zones provide the relations between prescribed tectonic boundary conditions and the resulting bulk flow field. They do not determine how structures and fabrics observed on local and small scales form. To bridge the scale gap between the bulk flow field and minor structures, Eshelby’s formalism extended for general power-law viscous materials is shown to be a powerful means. This paper first gives a complete presentation of Eshelby’s formalism, from the classic elastic inclusion problem, to Newtonian viscous materials, and to the most general case of a power-law viscous inhomogeneity embedded in a general power-law viscous medium. The formulation is then implemented numerically. The implications and potential applications of the approach are discussed. It is concluded that the general Eshelby formalism together with the self-consistent method is a powerful and physically sound means to tackle large plastic deformation of Earth’s lithosphere.     

新世纪构造地质学两大支柱理论: 最大有效力矩准则与变位形分解

传统构造地质学用摩尔-库伦准则和贝克尔的应变椭球体理念分别解释地壳中的脆性断层和塑性变形,将变形局部化的韧性剪切带形成解释为平行应变椭球体的圆切面,却无法解释变形局部化的共轭剪切带稳定夹角~110°面对应缩短方向。变形局部化是独立于脆性和塑性变形外的变形领域,受最大有效力矩准则控制。20世纪末提出的变位形分解理念,摆脱连续介质力学的束缚,合理地说明广泛存在的走滑断层平行俯冲带或逆冲断层带。非均匀变形和非连续介质力学理念的建立,为地质学与力学的结合开辟了新的前景。文章试用上述两理念概略分析中国和邻区中新生代构造格局,以期引发讨论。      

Small-scale ductile shear zones as transposed rheologically weak domains: A numerical modeling investigation and practical application

We test the hypothesis that small ductile shear zones are developed from initial rheologically weak domains. We regard weak domains as ellipsoidal inhomogeneities and apply Eshelby's formalism extended for power-law viscous materials to investigate the kinematics and finite strain evolution of the partitioned flow fields in weak domains. We show that, under an imposed bulk flow field, weak domains, regardless of their initial shapes and orientations, deform into zone-like features at relatively small bulk strains and the deformation paths inside weak domains have all characteristics expected in ductile shear zones. We apply our model to the Cap de Creus area Spain, where abundant small ductile shear zones exist. To relate the fabrics inside these shear zones with the regional deformation, we take a multi-scale approach. We assume that the area contained many weak domains which were randomly shaped and oriented initially and were deformed into shear zones eventually in response to the regional flow. We constrain the regional flow field by the fabric patterns across the area and compute numerically the partitioned flow fields in individual weak domains. The latter are related to fabrics inside shear zones. Our model reproduces first-order features of field-observed fabrics. Although the deformation path of each shear zone was close to simple shearing, the deformation of the whole belt was not. Our approach also resolves the strain compatibility problem for a finite-sized shear zone embedded in a far less deformed country rock.     

A fundamental problem with the kinematic interpretation of geological structures

Non-steady flows are ubiquitous in nature. Not only do imposed boundary conditions vary with time and rock rheology change during the course of deformation but also deformation is generally heterogeneous and all of these conditions lead to non-steady flow histories.In modern kinematic analysis, flow apophyses, the instantaneous stretching axes and the vorticity vector, collectively referred to as the `eigen directions' of deformation, are used in place of Bruno Sander's kinematic axes. For deformation with a steady flow history, this practice is well justified and has led to great advances in structural geology. But for non-steady flow histories, the geometrical relationships among eigen directions (flow pattern) vary with time. This makes it inappropriate to correlate an association of structures and fabrics with certain `time-invariant' flow patterns and hence raises the question: How should we interpret geological structures and fabrics kinematically, without the assumption of homogeneous and steady deformation? We suggest that the answer lies primarily in forward-modeling of deformation, based on a knowledge of rock properties.     

Continuum hydrodynamics of dry granular flows employing multiplicative elastoplasticity

We present a Lagrangian formulation for simulating the continuum hydrodynamics of dry granular flows based on multiplicative elastoplasticity theory for finite deformation calculations. The formulation is implemented within the smoothed particle hydrodynamics (SPH) method along with a variant of the usual dynamic boundary condition. Three benchmark simulations on dry sands are presented to validate the model: (a) a set of plane strain collapse tests, (b) a set of 3D collapse tests, and (c) a plane strain simulation of the impact force generated by granular flow on a rigid wall. Comparison with experimental results suggests that the formulation is sufficiently robust and accurate to model the continuum hydrodynamics of dry granular flows in a laboratory setting. Results of the simulations suggest the potential of the formulation for modeling more complex, field-scale scenarios characterized by more elaborate geometry and multi-physical processes. To the authors’ knowledge, this is the first time the multiplicative plasticity approach has been applied to granular flows in the context of the SPH method.     

Pure shear to simple shear-dominated transpression during the Eburnean major D2 deformation,Daléma sedimentary basin,eastern Senegal

Field studies in the Palaeoproterozoïc Daléma basin, Kédougou-Kéniéba Inlier, reveal that the main tectonic feature comprises alternating large shear zones relatively well-separated by weakly deformed surrounding rock domains. Analysis of the various structures in relation to this major D₂ phase of Eburnean deformation indicates partitioning of sinistral transpressive deformation between domains of dominant transcurrent and dominant compressive deformation. Foliation is mostly oblique to subvertical and trending 0–30° N, but locally is subhorizontal in some thrust-motion shear zones. Foliation planes of shear zones contain a superimposed subhorizontal stretching lineation which in places cross-cuts a steeply plunging stretching lineation which is clearly expressed in the metasedimentary rocks of weakly deformed surrounding domains. In the weakly deformed domains, the subhorizontal lineation is absent, whereas the oblique to subvertical lineation is more fully developed. Finite strain analyses of samples from surrounding both weakly deformed and shearing domains, using finite strain ratio and the Fry method, indicate flattened ellipsoid fabrics. However, the orientation of the long axis (X) of the finite strain ellipsoid is horizontal in the shear zones and oblique within the weakly deformed domains. Exceptionally, samples from some thrust zones indicate a finite strain ellipsoid in triaxial constriction fabrics with a subhorizontal long axis (X). In addition, the analysis of the strain orientation starting from semi-ductile and brittle structures indicates that a WNE–ESE (130° N to 110° N) orientation of strain shortening axis occurred during the Eburnean D₂ deformation.     

Deformation partitioning during transpression in response to Early Devonian oblique convergence,northern Appalachian orogen,USA

Transpressive deformation was distributed heterogeneously within the Central Maine belt shear zone system, which formed in response to Early Devonian oblique convergence during the Acadian orogeny in the northern Appalachians. ‘Straight’ belts are characterized by tight folds, S>L fabrics and sub-parallel form lines, and asymmetric structures that together indicate dextral–SE-side-up kinematics. In contrast, intervening zones between ‘straight’ belts are characterized by open folds and L≫S fabrics. Within both types of zone, metasedimentary rocks have fabrics defined by the same minerals at the same metamorphic grade, including a penetrative, moderately to steeply NE-plunging mineral lineation. Thus, we interpret accumulation of plastic deformation and regional metamorphic (re-) crystallization to have been synchronous across the Central Maine belt shear zone system. Discordance between inclusion trails in regionally developed porphyroblasts of garnet and staurolite and matrix fabrics in ‘straight’ belt rocks records shortening by tightening of folds and greater reorientation of matrix fabrics with respect to porphyroblasts. Kinematic partitioning of flow was responsible for the contrasting states of finite deformation recorded in the Central Maine belt shear zone system. Perturbations in the flow were caused by serially developed thrust-ramp anticlines in the stratigraphic succession immediately above the Avalon-like basement, at which décollement of the shear zone system was initially rooted. General shear deformation at the ramps involved strain softening with an enhanced component of noncoaxial flow. In contrast, deformation during extrusion in the intervening zones involved strain hardening with a greater component of coaxial flow. Part of the thickening stratigraphic succession exceeded T_solidus, reflected by the occurrence of migmatites and granites. The latter were partly sourced from the underlying Avalon-like basement that was involved in the deformation and melting.     

超高压——高压剪切带

超高压－高压剪切带是大陆碰撞造山过程中变质－变形分解作用形成的一类特殊的剪切带。剪切带内的构造以含柯石英、金刚石等超高压－高压矿物的榴辉岩、面理化超基性－基性岩、片麻岩、大理岩为特征，具有复杂的变形组构和变形历史。最终形成布丁－基质或残斑－基质结构。深入研究超高压－高压剪切带的几何学、运动学和流变学特征，对分析超高压－高压变质岩石的形成－折返过程及碰撞造山动力学，具有重要的实际和理论意义。     

大别—苏鲁区UHP变质岩构造学及流变学演化

在大别—苏鲁区的30个关键位置,对UHP/HP变质岩进行详细构造解析、大比例尺(1∶10000)制图并在区域尺度上进行观察和对比,以便揭示它们的构造几何学、变形条件和流变学演化。初步的研究结果指出,广泛出露的UHP/HP榴辉岩相岩石形成一个巨大的UHP/HP变质带,提供了一个观察中朝与扬子克拉通之间三叠纪大陆深俯冲-碰撞带过程的窗口。观察的显微构造及组构指出,UHP/HP变质带内岩石变形机制,无论是在榴辉岩相阶段还是在榴辉岩相后阶段,都是以塑性流变为主,其力学行为和组构特征都受组成矿物的强度、强度差等流变学特征,以及变形物理环境如压力、温度、应变速率、差异应力和流体含量等的制约。在俯冲/碰撞带内的变形分解作用于岩石圈不同层次及不同的构造阶段都曾发生,而且,在不同尺度上,应变局部化形成具高应变的剪切带网络,且一般显示典型的布丁-基质或碎斑-基质构造及流变学型式。根据构造、岩石、变质作用及地质年代学资料,借助于岩石圈流变学基本原理,提出一个大别—苏鲁区UHP/HP变质岩石流变学演化的工作模式,它涉及早期扬子与中朝克拉通间三叠纪(~250~230Ma)大陆深俯冲/碰撞、UHP/HP变质岩形成,相继深埋岩石的多期折返。特别强调UHP/HP岩石向地壳表层的折返,主要是构造过程,地面侵蚀作用是次要的。     

Analysis of three-dimensional strain modified uniform distributions: andalusite fabrics from a granite aureole

A method of analysing finite strain from an initial random distribution of lines or planes is presented. The shape of the strain ellipsoid is mapped to the eigenvalue ratios of an orientation matrix. The method is used to interpret the andalusite fabrics from the aureole of the Ardara granite, County Donegal, Ireland. The high flattening strains determined in the aureole are accompanied by oblate strains in the granite margins and support the concept of deformation produced by expansion of the pluton.