大巴山前陆西段叠加构造变形分析及其构造应力场特征

大巴山前陆西段发育典型的横跨叠加褶皱,清楚地记录了自身构造演化历史,为查明大巴山前陆构造体制转换提供了非常有利的条件。本文通过野外褶皱形态分析,结合层面和断层面滑动矢量的观测,获得大巴山前陆西段中生代主要的两期挤压构造应力场,即早期(T3—J1)的构造应力场以近南北向挤压为主,大巴山地区受其控制,发生近东西向展布褶皱变形,这期构造变形与印支运动有关;晚期,即J3—K1期间,大巴山构造体制发生重大转变,本区构造应力场由近南北挤压转变为近东西向强烈挤压,大巴山发生强烈构造变形,大巴山推覆构造带以断裂活动为主,镇巴—城口—房县断裂带强烈向大巴山前陆逆掩。而大巴山前陆则以强烈的近南北向延伸的褶皱变形为主,早期近东西向褶皱遭受强烈改造,形成大巴山前陆西段较为典型的横跨叠加褶皱构造。这期构造变形是燕山运动早期活动的结果。大巴山叠加变形的研究不仅为研究中国盆—山系统变形特征提供了具体例证,而且对大巴山前陆油气勘探起重要指导作用。     

机械光栅式测缝计应用于无砟轨道板上拱监测的可行性分析

高速铁路轨道的稳定性和平顺性是高速铁路正常运营的关键,因此对运营期高速铁路的轨下结构持续地开展变形监测是十分必要的。当前我国对轨道板变形的检测主要是依靠人工肉眼观察式的现场检查和常规水准测量方式进行监测,效率低下,难以在有限的天窗时间内完成辖区内轨道结构的全覆盖检测。基于此,从理论分析和实验测试两方面探讨机械光栅式测缝计应用于无砟轨道板上拱自动化监测的可行性,结果表明机械光栅式测缝计能够抵抗高温、淋水和振动等不良条件的影响,满足轨道板上拱监测的精度要求。     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Geodynamic modelling of the Century deposit,Mt Isa Province,Queensland

This paper is concerned with an example of quantitative modelling of orebody formation as a guide to reducing the risk for future mineral exploration. Specifically, the paper presents a detailed 3–D numerical model for the formation of the Century zinc deposit in northern Queensland. The model couples fluid flow with deformation, thermal transport and chemical reactions. The emphasis of the study is a systems approach where the holistic mineralising system is considered rather than concentrating solely on the mineral deposit. In so doing the complete plumbing system for mineralisation is considered with a view to specifying the critical conditions responsible for the ore deposit occurring where it does and having the size and metal grades that are observed. The numerical model is based on detailed geological, tectonic, isotopic and mineralogical data collected over the past 20 years. The conclusions are that the Century zinc deposit is located where it is because of the following factors: (i) a thermal anomaly is associated with the Termite Range Fault due to advection of heat from depth by fluid flow up the Termite Range Fault; (ii) bedding‐plane fissility in the shale rocks hosting the Century zinc deposit has controlled the wavelength and nature of D₁ folding in the vicinity of the deposit and has also controlled increases in permeability due to hydrofracture of the shales; such hydrofracture is also associated with the production of hydrocarbons as these shales passed through the ‘oil‐window’; (iii) Pb–Zn leached from crustal rocks in the stratigraphic column migrated up along faults normal to the Termite Range Fault driven by topographic relief associated with inversion at the end of the Isan Orogeny; these fluids mixed with H₂S derived at depth moving up the Termite Range Fault to mix with the crustal fluids to precipitate Pb–Zn in a plume downstream from the point of mixing. Critical factors to be used as exploration guides are high temperatures, carbonaceous fissile shales now folded into relatively tight D₁ folds, fault‐controlled plumbing systems that enable fluid mixing, depletion of metals upstream of the deposit and,in particular,a very wide Fe‐depletion halo upstream of the deposit.     

Micromechanical viscoelasto‐plastic models and finite element implementation for rate‐independent and rate‐dependent permanent deformation of stone‐based materials

This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd.     

Constraining sequence stratigraphy in north Australian basins: SHRIMP U–Pb zircon geochronology between Mt Isa and McArthur River

Fifty‐five new SHRIMP U–Pb zircon ages from samples of northern Australian ‘basement’ and its overlying Proterozoic successions are used to refine and, in places, significantly change previous lithostratigraphic correlations. In conjunction with sequence‐stratigraphic studies, the 1800–1580 Ma rock record between Mt Isa and the Roper River is now classified into three superbasin phases—the Leichhardt, Calvert and Isa. These three major depositional episodes are separated by ～20 million years gaps. The Isa Superbasin can be further subdivided into seven supersequences each 10–15 million years in duration. Gaps in the geological record between these supersequences are variable; they approach several million years in basin‐margin positions, but are much smaller in the depocentres. Arguments based on field setting, petrography, zircon morphology, and U–Pb systematics are used to interpret these U–Pb zircon ages and in most cases to demonstrate that the ages obtained are depositional. In some instances, zircon crystals are reworked and give maximum depositional ages. These give useful provenance information as they fingerprint the source(s) of basin fill. Six new ‘Barramundi’ basement ages (around 1850 Ma) were obtained from crystalline units in the Murphy Inlier (Nicholson Granite and Cliffdale Volcanics), the Urapunga Tectonic Ridge (‘Mt Reid Volcanics’ and ‘Urapunga Granite’), and the central McArthur Basin (Scrutton Volcanics). New ages were also obtained from units assigned to the Calvert Superbasin (ca 1740–1690 Ma). SHRIMP results show that the Wollogorang Formation is not one continuous unit, but two different sequences, one deposited around 1730 Ma and a younger unit deposited around 1722 Ma. Further documentation is given of a regional 1725 Ma felsic event adjacent to the Murphy Inlier (Peters Creek Volcanics and Packsaddle Microgranite) and in the Carrara Range. A younger ca 1710 Ma felsic event is indicated in the southwestern McArthur Basin (Tanumbirini Rhyolite and overlying Nyanantu Formation). Four of the seven supersequences in the Isa Superbasin (ca 1670–1580 Ma) are reasonably well‐constrained by the new SHRIMP results: the Gun Supersequence (ca 1670–1655 Ma) by Paradise Creek Formation, Moondarra Siltstone, Breakaway Shale and Urquhart Shale ages grouped between 1668 and 1652 Ma; the Loretta Supersequence (ca 1655–1645 Ma) by results from the Lady Loretta Formation, Walford Dolomite, the upper part of the Mallapunyah Formation and the Tatoola Sandstone between ca 1653 and 1647 Ma; the River Supersequence (ca 1645–1630 Ma) by ages from the Teena Dolomite, Mt Les and Riversleigh Siltstones, and Barney Creek, Lynott, St Vidgeon and Nagi Formations clustering around 1640 Ma; and the Term Supersequence (ca 1630–1615 Ma) by ages from the Stretton Sandstone, lower Doomadgee Formation and lower part of the Lawn Hill Formation, mostly around 1630–1620 Ma. The next two younger supersequences are less well‐constrained geochronologically, but comprise the Lawn Supersequence (ca 1615–1600 Ma) with ages from the lower Balbirini Dolomite, and lower Doomadgee, Amos and middle Lawn Hill Formations, clustered around 1615–1610 Ma; and the Wide Supersequence (ca 1600–1585 Ma) with only two ages around 1590 Ma, one from the upper Balbirini Dolomite and the other from the upper Lawn Hill Formation. The Doom Supersequence (<1585 Ma) at the top of the Isa Superbasin is essentially unconstrained. The integration of high‐precision SHRIMP dating from continuously analysed stratigraphic sections, within a sequence stratigraphic context, provides an enhanced chronostratigraphic framework leading to more reliable interpretations of basin architecture and evolution.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.     

Late-Hercynian shearing during crystallization of granitoid magmas (Sila massif,southern Italy): regional implications

Shearing of regional extent, involving granitoids and underlying mid-crustal rocks of the Sila massif (Calabria, Italy), is analysed in this paper. The deformed granitoids are affected by a wide NNW-SSE oriented deformation zone, stretching for about 60 km, from the neighbourhood of Cecita Lake to Cropani village. Meso- and micro-structures in granitoids, close to the boundary with underlying migmatitic paragneiss, indicate that deformation developed from melt-present to solid-state conditions. Simultaneous tectonics and magmatism activated a plutonic accretionary process at mid-crustal levels. This took place at about 300 Ma and involved hybrid magmas with a dominat contribution from a mantle source. The deformation regime remained steady for a long time during magma crystallization and cooling in subsolidus conditions. The regional top-to-the-W sense of shear in the present geographic coordinates, recorded in the deformed granitoids, seems geometrically consistent with the coeval direction of maximum extension found in another sector of the southern Hercynian belt, suggesting the original position of the Sila basement in this context. Magmatic ativity ended with the intrusion of mafic and felsic magams affected by a very weak deformation, ongoing during the final strain increments of the late-Hercynian stage.