断层附近地应力扰动带宽度及其影响因素

在建立正断层模型的基础上,利用三维有限元数值模拟方法研究了断层附近地层中的地应力变化规律。其结果表明,由于断层活动,在断层附近普遍发育应力扰动带;在应力扰动带范围内,地应力方向和大小发生明显的变化,断层中部附近应力值普遍较低,而断层端部的应力值通常异常增大。应力扰动带的分布范围主要受断层规模的控制,与断裂带的岩石力学性质、断层走向、断层面形态和边界应力条件等因素也密切相关。随断层长度和断距逐步增大,应力扰动带的宽度相应增加。断裂带的岩石越破碎,其岩石弹性模量越低,断层对地应力的影响宽度越大。断层走向与区域最大水平主应力方向越接近垂直、断距与断层长度的比值越大,区域内的差应力越大,则扰动带宽度与断层规模的比值也越大。选择渤海湾盆地BZ-X油田进行验证,在建立油田实际三维地质模型基础上,根据边界应力条件,利用三维有限元方法对沙河街组二段的地应力分布进行了数值模拟计算。根据断层周围的地应力变化,确定了应力扰动带的分布范围,断层附近应力扰动带宽度的分布规律与正断层模型分析结果相一致。      

Stream tube integration method for analysing subsurface fluid flow

A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd.     

GBW10010a大米标准物质复(研)制及数据特征

随着我国对生态文明建设的重视,自然资源综合调查势在必行,对生物标准物质亦提出了新的需求。当前相关调研工作已经大面积开展,自然资源综合调查、农产品与食品安全评价都需要对生物样品元素组成进行准确测试,需要以生物标准物质作为生物成分测试量值比对和溯源的基础,因此对生物基体标准物质的需求量大幅增加。大米作为主要粮食之一,其食品安全日益受到重视,对大米中的化学成分进行准确的分析测试具有重要的现实意义,因而对大米标准物质的需求量尤为突出,但目前大米成分分析标准物质已供不应求。本文严格按照《标准物质定值的通用原则及统计学原理》（JJF 1343—2012）和《地质分析标准物质的研制》（JJF 1646—2017）等相关规范要求,开展了GBW10010a大米成分分析标准物质的复（研）制工作,包括样品采集、加工制备、均匀性检验、稳定性检验、多家实验室协作定值测试及不确定度评定等关键环节。结果表明：本次复（研）制的大米标准物质定值成分多样、量值准确可靠,符合国家一级标准物质的要求。GBW10010a共定值54项主微量元素,包括Ag、Al、As、B、Ba、Be、Bi、Ca、Cd、Ce、Co、Cr、Cs、Cu、Dy、Er、Eu、Fe、Gd、Ge、Hg、Ho、K、La、Li、Ho、Mg、Mn、Mo、N、Na、Nb、Nd、Ni、P、Pb、Pr、Rb、S、Sb、Sc、Se、Si、Sm、Sr、Tb、Th、Tl、Tm、U、V、Y、Yb、Zn,其中的39项元素给出了标准值及不确定度,包括Ag、Al、As、B、Ba、Ca、Cd、Ce、Co、Cs、Cu、Dy、Er、Fe、Hg、K、Li、Mg、Mn、Mo、N、Na、Nd、Ni、P、Pb、Pr、Rb、S、Sb、Se、Si、Sm、Sr、Tb、Tl、Y、Yb、Zn;15项元素提供参考值,包括Be、Bi、Cr、Eu、Gd、Ge、Ho、Ho、La、Nb、Sc、Th、Tm、U、V。与原有GBW10010大米标准物质相比较,GBW10010a中As、Cd、Co、Cr、Cu、Hg、Mn、Mo、Ni、Zn等重金属元素含量显著下降,其中Cd、Cu、Zn降幅较大,分别下降约39%、43%、38.7%,一定程度上反映了农田生态环境的改善。本批标准物质定值元素总数量增加了6项,新增定值元素Ag、Nb（Nb给出参考值）,并且各项元素不确定度范围整体上有所缩小,如Al、Cd、Cu、Fe、K、Mg、Mo、Na、P、Pb、Se、Zn等对生物易有影响的重要元素,表明了地质分析测试方法技术的进步及定值水平的提高。本批标准物质定值元素涵盖了具有生物效应的大部分主微量元素,适用于农业生态环境地球化学调查与评价、生物样品测试、农产品质量与食品安全评价样品测试时的分析仪器校正、分析方法评价和分析质量监控等多个领域。     

A Discrete Wavenumber Boundary Element Method for study of the 3-D response 2-D scatterers

A mathematical formulation of the 2·5D elastodynamic scattering problem is presented and validated. The formulation is a straightforward extension of the Discrete Wave number Boundary Integral Equation Method (DWBIEM) originally proposed by Kawase¹ for 2D scattering problems and subsequently extended to the 3D problem by Kim and Papageorgiou.² It is demonstrated that the Green's function which is appropriate for a boundary formulation of the 2·5D elastodynamic scattering problem is the one corresponding to a unit force moving on a straight line with constant velocity. Such a Green's function is derived in the present study. The formulation may be used to study the wavefields in models of sedimentary deposits (e.g. valleys) or topography (e.g. canyons or ridges) with a 2D variation in structure but obliquely incident plane waves. The advantage of a 2·5D formulation is that it provides the means for calculations of 3D wavefields in scattering problems by requiring a storage comparable to that of the corresponding 2D calculations. © 1998 John Wiley & Sons, Ltd.     

榍石LA-ICP-MS U-Pb定年中元素分馏的影响及校正研究

同位素地质年代学是探索地质体时空演化及地球动力学等问题的基础学科,应用最为广泛的当属含铀副矿物的U-Pb定年技术。榍石具有相对较低的U-Pb体系封闭温度,并广泛发育于岩浆岩、各类变质岩、热液成因岩石以及少量沉积岩中,是一种理想的中高温地质事件定年矿物。利用电感耦合等离子体质谱法(ICP-MS)测定榍石U-Pb年龄时,不可避免地要解决高普通铅以及元素分馏效应对测试的影响。本文对榍石LA-ICP-MS实验过程中的元素分馏行为进行研究,采用相同基体的标准样品与未知样品对比,发现了榍石不同颗粒之间元素分馏行为不一致的现象;同时采用不同的元素分馏校正方法,分别应用于锆石、独居石和榍石进行对比研究,认为分馏行为一致的副矿物定年可以采用"指数法"和"均值法"对数据进行校正,但是对于榍石这种分馏行为不一致的副矿物,定年时只有采用"截距法"对数据进行校正才可以获得正确的年龄。进而将此结论应用于秦岭造山带老牛山地区岩浆成因榍石样品,得到的结果与锆石年龄一致,表明"截距法"可以避免分馏行为不一致导致的校正不准确的问题。本研究成果为榍石LA-ICP-MS U-Pb定年方法的完善提供了一种思路。     

A continuum model of layered rock masses with non-associative joint plasticity

Layered rock masses can be modelled either as standard, orthotropic continua if the layer bending can be neglected or as Cosserat continua if the influence of layer bending is essential. This paper presents a finite element smeared joint model based on the Cosserat theory. The layers are assumed to be elastic with equal thickness and equal mechanical properties. All the cosserat parameters are expressed through the elastic properties of layers, layer thickness and joint stiffness. Plastic-slip as well as tensile-opening of layer interface (joint) are accounted for in a manner similar to the conventional non-associative plasticity theory. As an application, the behaviour of an excavation in a layered rock mass is examined. The displacement and stress fields given by smeared joint models based on the Cosserat continuum and the conventional anisotropic continuum approaches are compared with those obtained from the discrete joint model. The conventional anisotropic continuum model is found to break-down completely when the effective shear modulus in the direction parallel to layering is low in comparison to the shear modulus of the intact layer, whereas the Cosserat model is found to be capable of accurately reproducing complex load–deflection patterns irrespective of the differences in shear moduli. © 1998 John Wiley & Sons, Ltd.     

The complex one-sided integrals of Cauchy and Hadamard and application to boundary element method

The definitions of complex integrals of Cauchy and Hadamard with the singular point coinciding with the end point of the integration curve are proposed. It is shown that the new integrals satisfy most of the properties of the regular ones, including the change of variables. It is also shown that the Cauchy principal value (CPV) and Hadamard finite-part (HFP) integrals can be considered as a sum of the new type integrals. The application to numerical solution by the boundary element method (BEM) and the complex hypersingular integral equation (CHSIE) for the multiregions of interacting elastic bodies and bodies with cracks and holes is discussed. The different ways to place the collocation points are considered. The numerical results for the problems of circular hole and circular elastic inclusion in infinite plate indicated that the appropriate choice of the approximating functions leads to a high accuracy of the calculation. Applications of the new technique to geomechanics problems are discussed. © 1998 John Wiley & Sons, Ltd.     

数字地图图像的基本形态学变换和结构元素的分解

本文是作者关于“数学形态学与数字地图图像的分析和识别”研究成果的一部分。文章提出了地图图形的集合概念和地图图像形态分析的基元确定原则,证明了结构元素的基元最小分解定理;同时,介绍了基本的形态学变换,统一并扩展了带约束条件的数学形态学变换模型;从而,奠定了数字地图图像分析和识别的数学基础。     

Analysis of settlement and lateral deformation of soft clay foundation beneath two full-scale embankments

In 1986, the Malaysian Highway Authority constructed a series of trial embankments on the Muar Plain (soft marine clay) with the aim of evaluating the effectiveness of various ground improvement techniques. This study investigates the effect of two such ground improvement schemes: (a) preloading of foundation with surface geogrids and synthetic vertical drains and (b) sand compaction piles. The paper is focused on the finite element analysis of settlements and lateral displacements of the soft foundation. In scheme (a), the numerical predictions are compared with the field measurements. In scheme (b), only the numerical analysis is presented and discussed in the absence of reliable measurements due to the malfunctioning of the electronic extensometer and inclinometer system during embankment construction. The current analysis employs critical state soil mechanics, and the deformations are predicted on the basis of the fully coupled (Biot) consolidation model. The vertical drain pattern is converted to equivalent drain walls to enable plane strain modelling, and the geogrids are simulated by linear interface slip elements. The effect of sand compaction piles is investigated considering both ideal drains and non-ideal drains, as well as varying the pile stiffness. © 1997 by John Wiley & Sons, Ltd.     

Elastoplastic deformation around underground openings under biaxial initial stress field

Estimation of elastoplastic deformation around an underground opening induced by the excavation of it, especially displacement and strain field in plastic region, is presented in this paper, as well as the formulation for calculating the displacement and strain in the plastic region around the underground opening by the coupled Boundary Element Method - Characteristics Method (BEM-CM). In this method, the non-associated flow rule is adopted to calculate the displacement and strain field in the plastic region, which is determined by the integration of the displacement along characteristics lines under the boundary condition of the elastic displacement on an elastoplastic interface analysed. It is shown that this method is one of the accurate and effective methods for estimating not only the shape and extent of the plastic region but also the state of the displacement and strain in the plastic region around the underground opening, comparing the theoretical solution with numerical results by this method for a circular opening under hydrostatic initial stress condition. Furthermore, this method is applied to rectangular and horse-shoe shaped openings and the characteristics of the strain field in the plastic region are discussed. © 1997 by John Wiley & Sons, Ltd.