Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling     

Mark Ashworth  Florian Doster 《国际地质力学数值与分析法杂志》2020,44(17):2304-2328

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.  相似文献   

剪切波分裂揭示的青藏高原上地壳地震各向异性基本特征          下载免费PDF全文

郭铁龙  高原 《地球物理学报》2020,63(3):1085-1103

青藏高原整体隆升,构造运动与介质变形强烈,然而由于地震观测数据不足,青藏高原内部上地壳各向异性研究一直是一个空白.本研究使用西藏地区的地震台网(2009年5月—2017年5月)的观测资料,利用剪切波分裂研究青藏高原上地壳地震各向异性特征.由于青藏高原固定地震台站分布稀疏,可用于进行剪切波分裂研究的近场地震事件记录稀少,本研究采用地震事件的单台定位技术,对公开的地震目录里没有震源深度数据的地震事件进行震源位置约束,并引入微震模板匹配定位方法,对连续地震波形进行检索,识别出地震目录里遗漏的新的微震(小地震)事件波形.微震识别获得的新地震事件记录是地震目录里报告的地震事件记录的大约6倍,用于补充研究区的剪切波分裂数据分析.通过数据分析,对比快波偏振方向,证实微震识别获得的数据极大地增加了有效数据的数量,提高了结果的可靠性.研究结果表明,雅鲁藏布江缝合带与班公—怒江缝合带之间的拉萨地块东部地区,台站的快剪切波(快波)偏振方向主要受区域应力场影响,快波偏振方向主要是NS或NNE方向,表明了区域最大主压应力方向;但个别地震台站(当雄台)快波偏振方向受原地主压应力影响,其快波偏振方向既不平行于断裂走向也不平行于区域主压应力方向,揭示出地壳介质的局部变形导致的局部应力方向不同于青藏块体里的其他地区.研究区西部的改则、普兰和研究区北部的双湖,快波偏振方向显示与断裂等构造走向一致的特点.研究区东部的昌都和察隅,快波偏振方向除了与断裂走向(或构造线)一致,还与地表运动的方向相同,揭示了青藏块体东部的深部物质可能的运移方向.这个现象虽然还需更多的研究证实,但这个发现的重要启示是,地震各向异性结合地表变形可用于探讨地壳深部物质的运动.  相似文献   

A simple hypoplastic model with loading surface accounting for viscous and fabric effects of clays     

Merita Tafili  Theodoros Triantafyllidis 《国际地质力学数值与分析法杂志》2020,44(16):2189-2215

  相似文献   

A plasticity constitutive model for unsaturated,anisotropic, nonexpansive soils     

Panagiotis Sitarenios  Michael Kavvadas 《国际地质力学数值与分析法杂志》2020,44(4):435-454

The paper describes and evaluates an incremental plasticity constitutive model for unsaturated, anisotropic, nonexpansive soils (CMUA). It is based on the modified Cam-Clay (MCC) model for saturated soils and enhances it by introducing anisotropy (via rotation of the MCC yield surface) and an unsaturated compressibility framework describing a double dependence of compressibility on suction and on the degree of saturation of macroporosity. As the anisotropic and unsaturated features can be activated independently, the model is downwards compatible with the MCC model. The CMUA model can simulate effectively: the dependence of compressibility on the level of developed anisotropy, uniqueness of critical state independent of the initial anisotropy, an evolving compressibility during constant suction compression, and a maximum of collapse. The model uses Bishop's average skeleton stress as its first constitutive variable, favouring its numerical implementation in commercial numerical analysis codes (eg, finite element codes) and a unified treatment of saturated and unsaturated material states.  相似文献   

乌海地区剪切波分裂研究     

魏建民  张晖  赵星  张浩鸣  郝亮  阿那尔 《震灾防御技术》2020,15(4):777-787

根据乌海地区构造环境，采用SAM方法研究乌海地区地壳各向异性特征，使用乌海地震台2014年1月至2020年6月数字地震波形进行分析。根据65个有效地震记录，得到乌海地区剪切波分裂参数，其中快剪切波平均优势偏振方向为NE63.1°±46.4°，慢剪切波平均时间延迟为（1.13±0.66）ms/km。乌海地震台快剪切波偏振显示出4个优势偏振方向，分别为NE、EW、NNE、NNW向。将得到的各向异性结果与研究区应力场和地质构造进行分析，认为研究区周边复杂的剪切波分裂变化是主压应力场、原地主压应力、断裂带分布共同作用的结果。  相似文献   

利用远震研究四川盆地及其周缘瑞利面波相速度和方位各向异性          下载免费PDF全文

王小龙  雷兴林  马胜利  张烨  郭志  闻学泽 《地球物理学报》2020,63(2):445-459

基于四川盆地及周边的245个宽频带台站2010年9月—2014年9月期间的远震记录,提取双台路径瑞利面波相速度频散资料,反演得到四川盆地20～120s的高分辨率瑞利面波相速度及各向异性空间分布.在丰富区域地球物理基础数据的同时,结合已有研究成果对地壳上地幔变形耦合进行探讨,结果表明短周期(20～30s)的相速度分布与四川盆地的地质构造特征相吻合,作为川滇地块、松潘—甘孜地块和四川盆地之间的边界——龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用;松潘—甘孜地块,特别是川滇地块中下地壳普遍表现为明显的低速异常,表明中下地壳相对软弱;而四川盆地的中下地壳整体呈现相对高速,表明四川盆地具有相对坚硬的中下地壳.研究区域东南角接近北扬子地块与南扬子地块的缝合部位,呈现高速异常.四川盆地南部和东南邻区不同周期均具有较强的各向异性,且快波方向较为一致,反映这些地区不同深度变形耦合较好.四川盆地西部、北部及东北部邻区,不同周期的各向异性快波方向变化较大,不同深度变形耦合较差.这些特征与绕喜马拉雅东构造结的物质流动被扬子地块的高速地壳阻挡的宏观认识基本一致.  相似文献   

Seismic amplitude inversion for the transversely isotropic media with vertical axis of symmetry     

Feng Zhang  Tuo Zhang  Xiang-Yang Li 《Geophysical Prospecting》2019,67(9):2368-2385

Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data.  相似文献   

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads     

Juan Wang  Hai‐Sui Yu 《国际地质力学数值与分析法杂志》2014,38(4):331-348

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

岩体结构面各向异性变异系数的尺寸效应研究     

洪陈杰  黄曼  夏才初  罗战友  杜时贵 《岩土力学》2020,41(6):2098-2109

岩体结构面的各向异性与尺寸效应特征对其力学性质影响较大，综合考虑两者之间的相互关系对工程岩体稳定性评价具有重要意义。全面分析不同尺度结构面的各向异性分布，提出考虑正交方向三维形貌参数的各向异性变异系数AVC3D。通过渐进覆盖法统计4组天然岩体结构面10种不同采样尺寸的各向异性变异系数及其尺寸效应规律。结果表明：各向异性变异系数均随结构面尺寸增大而减小至定值，其与结构面尺寸具有较好的负对数函数关系；通过归一化处理，各向异性变异系数与结构面尺寸呈现较好的线性函数关系，表明各向异性变异系数具有分形结构，分维数D可以实现其尺寸效应的规律统计。该方法揭示了影响结构面形貌特征方向性变化的最大有效倾角 和粗糙系数C的增长幅度会随采样尺寸增大趋于稳定的变化机制，体现了当结构面达到各向异性尺寸效应的阈值后，会呈现出稳定的各向异性规律。  相似文献   

内蒙古阿巴嘎地区地壳方位各向异性研究          下载免费PDF全文

强正阳  吴庆举  何静  李永华 《地球物理学报》2019,62(8):2946-2958

本文利用架设在内蒙古阿巴嘎地区38个宽频带地震台站记录到的远震数据，通过拟合P波接收函数径向Pms转换波到时和叠加不同方位切向分量，确定了地壳各向异性参数.结果表明，大部分台站Pms延迟时间在0.35 s左右；而少数台站时差较大，推测可能受到索伦缝合带附近地壳残留倾斜界面影响.各向异性快波方向变化范围在N95°E到N180°E之间，平均为N130.6°E±19.1°，推测中下地壳矿物在ENE-WSW向区域主压应力作用下发生晶格定向排列可能是导致地壳各向异性的主要成因.研究区壳幔变形特征和机制不同，属于解耦变形.  相似文献