巴里坤地震台SQ-70DSi水平摆观测数据干扰分析

基于2018—2021年巴里坤地震台水平摆观测资料,选取典型受干扰曲线,分析其变化趋势,总结仪器运行中存在的各类干扰因素。研究发现,干扰主要包括自然干扰(降雨、气压、雷击)、人为干扰(进洞巡检、调零标定、仪器维修)、仪器自身故障(数采故障、光源灯故障)等。     

基于GPS的青藏块体东北缘近期地壳水平运动与变形研究

以"中国地壳运动观测网络"工程1999-2001年GPS观测数据以及"大陆强震机理与预测"国家基础研究项目青藏块体东北缘GPS网1999,2001,2003年的观测资料为基础,应用GAM IT/GLOBK软件对数据进行精密解算,获得了青藏块体东北缘1999-2001年相对于鼎新基准站以及欧亚板块的水平运动速度场。结合板块构造以及大地动力学理论,对青藏块体进行了块体的划分,分析了块体间的相互作用,并进一步分析了该地区的地壳水平运动特征,得出了以下结论:(1)GPS技术用于地壳运动的监测是可行的;(2)1999-2001年期间,祁连-海原断裂带GPS运动位移与地质结果十分吻合;(3)印度板块与欧亚板块的碰撞以碰撞后印度大陆继续以50 mm/a的速度向北推进是青藏块体东北缘的运动动力;(4)甘青块体西部子块体以北东向挤压运动为主;东部子块体主要表现为以海原断裂带为北边界的顺时针旋转为主;(5)鄂尔多斯块体存在逆时针旋转运动,其旋转速度在块体西南侧的六盘山断裂带一带最大,向北至块体西北部明显减小。     

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.     

Seasonal evaluation of rainfall estimation by four cumulus parameterization schemes and their sensitivity analysis

Sustainable water resources management require scientifically sound information on precipitation, as it plays a key role in hydrological responses in a catchment. In recent years, mesoscale weather models in conjunction with hydrological models have gained great attention as they can provide high‐resolution downscaled weather variables. Many cumulus parameterization schemes (CPSs) have been developed and incorporated into three‐dimensional Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) mesoscale model 5 (MM5). This study has performed a comprehensive evaluation of four CPSs (the Anthes–Kuo, Grell, Betts–Miller and Kain–Fritsch93 schemes) to identify how their inclusion influences the mesoscale model's precipitation estimation capabilities. The study has also compared these four CPSs in terms of variability in rainfall estimation at various horizontal and vertical levels. For this purpose, the MM5 was nested down to resolution of 81 km for Domain 1 (domain span 21 × 81 km) and 3 km for Domain 4 (domain span 16 × 3 km), respectively, with vertical resolutions at 23, 40 and 53 vertical levels. The study was carried out at the Brue catchment in Southwest England using both the ERA‐40 reanalysis data and the land‐based observation data. The performances of four CPs were evaluated in terms of their ability to simulate the amount of cumulative rainfall in 4 months in 1995 representing the four seasonal months, namely, January (winter), March (spring), July (summer) and October (autumn). It is observed that the Anthes–Kuo scheme has produced inferior precipitation values during spring and autumn seasons while simulations during winter and summer were consistently good. The Betts–Miller scheme has produced some reasonable results, particularly at the small‐scale domain (3 km grid size) during winter and summer. The KF2 scheme was the best scheme for the larger‐scale (81 km grid size) domain during winter season at both 23 and 53 vertical levels. This scheme tended to underestimate rainfall for other seasons including the small‐scale domain (3 km grid size) in the mesoscale. The Grell scheme was the best scheme in simulating rainfall rates, and was found to be superior to other three schemes with consistently better results in all four seasons and in different domain scales. 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.     

Nonlinear analysis of laterally loaded rigid piles in cohesive soil

This article presents a method for the nonlinear analysis of laterally loaded rigid piles in cohesive soil. The method considers the force and the moment equilibrium to derive the system equations for a rigid pile under a lateral eccentric load. The system equations are then solved using an iteration scheme to obtain the response of the pile. The method considers the nonlinear variation of the ultimate lateral soil resistance with depth and uses a new closed‐form expression proposed in this article to determine the lateral bearing factor. The method also considers the horizontal shear resistance at the pile base, and a bilinear relationship between the shear resistance and the displacement is used. For simplicity, the modulus of horizontal subgrade reaction is assumed to be constant with depth, which is applicable to piles in overconsolidated clay. The nonlinearity of the modulus of horizontal subgrade reaction with pile displacement at ground surface is also considered. The validity of the developed method is demonstrated by comparing its results with those of 3D finite element analysis. The applications of the developed method to analyze five field test piles also show good agreement between the predictions and the experimental results. The developed method offers an alternative approach for simple and effective analysis of laterally loaded rigid piles in cohesive soil. Copyright © 2011 John Wiley & Sons, Ltd.     

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.