基于多孔介质理论的地基土变形模量估算方法

土的变形模量是基础沉降弹性分析理论所必需的基本参数,而目前我国的岩土工程勘察报告一般并不提供土的变形模量。通过对前人研究结果的总结分析,基于多孔介质理论,并考虑土的泊松比、孔隙比及土体扰动等影响因素,笔者提出了一种根据土的压缩模量估算变形模量的方法,对于应用弹性理论计算基础沉降和充分利用已有研究资料都具有实际的意义。     

Consolidation analysis of a cross-anisotropic homogeneous elastic soil using a finite layer numerical method

This paper presents a finite layer procedure for Biot's consolidation analysis of layered soils using a cross-anisotropic elastic constitutive model. The program is first verified using published results. Then, using this program, the influences of cross-anisotropy on the immediate settlement, the final settlement, and the consolidation behaviour are investigated by changing one model parameter at one time. The results obtained using the cross-anisotropic elastic model are compared with results using an isotropic elastic model. It is found that the cross-anisotropy has very large influences on the immediate settlement, the final settlement, and the consolidation behaviour. Curves or tables of the immediate settlement coefficients, the final settlement coefficients, and the average degree of consolidation are obtained and presented in the paper. These curves or tables can be easily used to estimate the immediate settlement, the final settlement, and the consolidation settlement of a cross-anisotropic soil. Copyright © 2004 John Wiley & Sons, Ltd.     

准各向同性裂隙岩体中有效动弹性参数与弹性波速关系的研究

从弹性波在介质中传播的特点和弹性波法在岩体工程中的实际应用情况方面,分析了弹性波在含裂隙岩体介质中的传播特征。将岩体中裂隙呈随机分布的裂隙岩体似为准各向同性裂隙岩体,根据能量平衡原理、岩石断裂力学理论和卡氏最小功能定理(Castiglano’s theorem),确定了准各向同性裂隙岩体的有效动弹性参数( 、 、 )与裂隙密度参数(Xv)的数学关系。在此基础上,根据各向同性介质中传播的弹性波理论,建立了准各向同性裂隙岩体中弹性波速( 、 )和 / 与裂隙密度参数(Xv)之间的关系。理论研究表明,在此裂隙岩体中,有效动弹性参数和弹性波速与裂纹密度参数之间相互的关系均近似成倒数函数的非线性关系;裂隙密度参数增加,有效动弹性模量( )也减少,弹性波速( , )也减小;在裂隙密度参数较小时,此裂隙岩体的有效弹性模量和弹性波速随裂隙密度参数的变化较大,说明裂隙体的有效弹性参数和弹性波速对裂隙体非常敏感。     

Application of extended Kalman filter to back analysis of the natural stress state accounting for measuring uncertainties

In this work, we present a numerical procedure for determining the nature stress state in the rock mass around a tunnel. A finite element method is applied for analyzing the direct problems of tunneling during the back analysis of parameter estimation, in which a no‐tension elastic–plastic model is used to simulate the elastic–tensile and elastic–plastic‐tensile failure states which often occur in the cases of underground excavation in heavily jointed rock masses. By considering the natural stress state as random parameters of the tunneling system, the Kalman filter method is employed for feedback analysis to modify the parameter values in a statistical context, which uses the prior information in the process of estimation and employs a set of displacements obtained from field measurements. To verify the effectiveness of the proposed method of inverse analysis, the developed numerical procedure is applied to a synthetic example of deep tunnels in yielding rock masses. The relative importance of the a priori and updating information is investigated, as is the importance of their uncertainty. The results show great potential of the proposed approach. Copyright © 2010 John Wiley & Sons, Ltd.     

PROPAGATION OF TORSIONAL SURFACE WAVES IN VISCOELASTIC MEDIUM

It is well known that an elastic homogeneous half-space does not allow torsional surface waves to propagate. The present paper attempts to find out the possibility of propagation of such waves in a viscoelastic half-space. The study reveals that although the homogeneous elastic half-space does not allow torsional surface waves to propagate, a viscoelastic half-space does so. The wave is damped due to the viscoelastic parameter. It has also been found that as the viscoelastic parameter decreases, the medium becomes elastic and the torsional surface waves ceases to propagate.     

An anisotropic elastoplastic model for soft clays based on logarithmic contractancy

A new constitutive model for soft structured clays is developed based on an existing model called S‐CLAY1S, which is a Cam clay type model that accounts for anisotropy and destructuration. The new model (E‐SCLAY1S) uses the framework of logarithmic contractancy to introduce a new parameter that controls the shape of the yield surface as well as the plastic potential (as an assumed associated flow rule is applied). This new parameter can be used to fit the coefficient of earth pressure at rest, the undrained shear strength or the stiffness under shearing stress paths predicted by the model. The improvement to previous constitutive models that account for soil fabric and bonding is formulated within the contractancy framework such that the model predicts the uniqueness of the critical state line and its slope is independent of the contractancy parameter. Good agreement has been found between the model predictions and published laboratory results for triaxial compression tests. An important finding is that the contractancy parameter, and consequently the shape of the yield surface, seems to change with the degree of anisotropy; however, further study is required to investigate this response. From published data, the yield surface for isotropically consolidated clays seems ‘bullet’ or ‘almond’ shaped, similar to that of the Cam clay model; while for anisotropically consolidated clays, the yield surface is more elliptical, like a rotated and distorted modified Cam clay yield surface. © 2015 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons Ltd.     

Charge – transfer energy in computer modeling of structure and properties of minerals

As is known, a variable parameter for partial ionicity is often involved in the techniques of modeling structure and properties of inorganic crystals by energy-minimization methods. In such a case, to make a calculated energy comparable to experimental estimate (atomization energy) there is a need to involve an intraatomic energy term, the so-called charge-transfer energy (Urusov and Eremin 1995). For this purpose a set of charge-transfer energies for valence states of many elements is compiled. The procedure based on aspiration of equalizing the electronic levels electronegativities is proposed to estimate the total contribution of the charge-transfer energy in the cohesive energy of a crystal. The other prescribed quantities are repulsive potential parameters which are related to the energetic propertis of valence shells of bonded atoms. The final potential function consists of an effective Coulomb attraction, repulsion and covalent Morse type term, properly weighted. The energy-minimization procedure is performed for various values of the ionicity degree parameter of the crystal under study. Thereafter the charge-transfer energy correction is introduced in order for the cohesive energy calculations to be made. The calculations of structure, energy and elastic properties were performed using the METAPOCS code with ionicity variation for halides (NaCl, CaF₂), oxides (MgO, Al₂O₃, SiO₂ quartz and SnO₂), sulphide ZnS and silicate CaSnSiO₅. The resulting cohesive energies usually agree with experiment much better than those calculated in pure ionic or covalent approximations. As to elastic properties of crystals under consideration, they are generally functions of the ionicity parameter too and in most instances the predicted properties at intermediate values of the ionicity degree are in a satisfactory accordance with experimental data.     

一种土的非线性弹性本构模型参数的反演方法

旨在提出一种土的非线性弹性本构模型参数反演的方法。以现今普遍实行的地基载荷试验为基础,依据遗传算法的组合优化理论,采用正演计算和遗传算法优化相结合的方式,建立了土层非线性弹性本构模型参数反演的方法;并依据某黄土场地地基载荷试验数据,实施了黄土土层非线性弹性本构模型参数反演的全过程。计算结果表明,所建立的方法可以实现土层非线性弹性本构模型中相互关联的多个参数的组合优化,并在对初始值要求较低的情况下,可以获得良好的参数反演精度。从而为土的变形特性分析和土与其中及相邻结构的共同作用分析,提供了较好的土体本构模型参数的确定方法。     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。     

Numerical analysis of axially loaded vertical piles and pile groups

A numerical method, based on a simplified elastic continuum boundary element method, is presented for the settlement analysis of axially loaded vertical piles and pile groups. The soil flexibility coefficients are evaluated using the analytical solutions for a layered elastic half space. Results are presented and compared with existing published solutions for the following cases: (i) piles in homogeneous soil, (ii) piles in finite soil layer, (iii) piles end-bearing on stiffer layer, (iv) piles socketted into stiffer bearing layer, and (v) piles in Gibson soil. Reasonably good agreement is obtained between the present solutions and existing published solutions.     

基于量测位移的隧洞围岩弹性抗力系数反演方法

在水工隧洞工程中,岩石弹性抗力系数K反映了岩石弹性抗力的大小,是设计环节中需要确定的一个极为重要的基础计算参数。水工隧洞工程中需用水压法试验测定弹性抗力系数,由于水压法测定抗力系数不仅耗时、费工、且试验工艺复杂。本文通过对比分析弹性状态下深埋圆形隧洞开挖支护与水压法试验力学效应关系,提出一种基于量测隧洞位移的围岩抗力系数的反演方法,该法可根据量测隧洞位移直接反演得到弹性抗力系数。结合北京市南水北调配套工程南干渠工程试验段工程实践,对浅埋水工隧洞位移进行监测,利用提出的反演方法,得到了隧洞围岩的弹性抗力系数。提出的反演方法为复杂水工隧洞不同围岩应力状态问题的弹性抗力系数反演提供参考和借鉴。     

卸荷模量数取值的研究

以三峡二期围堰计算为例,说明Duncan-Chang非线性弹性E-μ模型中的卸荷模量数kur对计算精度有一定的影响。通过试验资料得出加荷模量数k取决于土的成分和密度,并且得到加荷模量数与卸荷加荷模量数之比的经验关系,因而卸荷模量数可直接根据加荷模量数的大小来确定,由此所确定的卸荷模量数比较符合实际。     

修正剑桥模型参数对计算结果的影响

修正剑桥模型是土工计算中广泛应用的弹塑性模型,模型参数的取值是计算分析的重要内容之一,运用实测位移资料采用反分析方法可以有效地确定模型参数。通过有限元计算分析了模型中各参数的变化引起的钻孔灌注桩码头前沿的水平位移的改变,计算结果表明,影响水平位移的主要参数是与弹性模量有关的无量纲参数,与弹性模量有关的指数,泊松比和原始各向等压加荷试验常数,其中尤以与弹性模量有关的无量纲参数对计算结果的影响最为显著,这为反分析参数的选择提供了依据。     

应用CH4/N2指标估算塔里木盆地天然气热成熟度

通过对塔里木盆地满加尔凹陷低成熟煤岩在开放在线程序升温体系的热模拟实验,获得煤岩在不同温度点的热分解产物CH4、CO和N2,并对模拟后残留产物进行了镜质体反射率测试,建立了镜质体反射率与对应受热温度之间的关系。煤岩镜质体反射率与热温度之间的二阶关系表达为:Ro=0.0014×T+0.109,r=0.9931(Ro<0.6%);Ro=0.0067×T-1.5855,r=0.9996(Ro>0.6%)。通过上述方程建立CH4/N2值与煤岩镜质体反射率之间的对应关系,并利用CH4/N2值对塔里木盆地库车坳陷天然气热成熟度进行了预测。预测结果与实际地质分析结果相吻合,说明CH4/N2值可以作为塔里木盆地煤成气热成熟度预测指标应用于油气勘探中。     

Effect of fabric on shear wave velocity in granular soils

Acta Geotechnica - The small-strain elastic shear wave velocity ($$V_S$$) is a basic mechanical property of soils and is an important parameter in geotechnical engineering. Recently, $$V_S$$ has...     

On hyperelasticity for clays

Under purely elastic conditions a constitutive model for clays should be conservative. This can be achieved by using a hyperelastic stress–strain relation derived from an elastic potential. To circumvent problems connected with formulation of a suitable potential a simple concept of back-stress elasticity is proposed. Several elastic models are critically reviewed, two of which had been originally incorrectly classified as hyperelastic. A novel elastic potential which can be used to generate various hyperelastic models is proposed. Its salient feature is that the resulting stiffness is a homogeneous function of order 1 with respect to stress. An example of Mathematica calculation with tensorial differentiation is enclosed (Wolfram S. Mathematica. A system for doing mathematics by computer. Addison-Wesley Publishing, 1988). ©     

Mechanical properties of granular materials: A variational approach to grain‐scale simulations

The mechanical properties of cohesionless granular materials are evaluated from grain‐scale simulations. A three‐dimensional pack of spherical grains is loaded by incremental displacements of its boundaries. The deformation is described as a sequence of equilibrium configurations. Each configuration is characterized by a minimum of the total potential energy. This minimum is computed using a modification of the conjugate gradient algorithm. Our simulations capture the nonlinear, path‐dependent behavior of granular materials observed in experiments. Micromechanical analysis provides valuable insight into phenomena such as hysteresis, strain hardening and stress‐induced anisotropy. Estimates of the effective bulk modulus, obtained with no adjustment of material parameters, are in agreement with published experimental data. The model is applied to evaluate the effects of hydrate dissociation in marine sediments. Weakening of the sediment is quantified as a reduction in the effective elastic moduli. Copyright © 2008 John Wiley & Sons, Ltd.     

Gravitational potential as a source of earthquake energy

Some degree of tectonic stress within the earth originates from gravity acting upon density structures. The work performed by this “gravitational tectonics stress” must have formerly existed as gravitational potential energy contained in the stress-causing density structure.According to the elastic rebound theory (Reid, 1910), the energy of earthquakes comes from an elastic strain field built up by fairly continuous elastic deformation in the period between events. For earthquakes resulting from gravitational tectonic stress, the elastic rebound theory requires the transfer of energy from the gravitational potential of the density structures into an elastic strain field prior to the event.An alternate theory involves partial gravitational collapse of the stress-causing density structures. The earthquake energy comes directly from a net decrease in gravitational potential energy. The gravitational potential energy released at the time of the earthquake is split between the energy released by the earthquake, including work done in the fault zone and an increase in stored elastic strain energy. The stress associated with this elastic strain field should oppose further fault slip.     

叠前弹性参数反演及FFP分析技术在东海A气田开发中的应用

东海盆地西湖凹陷A气田主力储层为浅水三角洲分流河道沉积砂体,储层横向变化快,非均质性强,实施开发井存在较大风险,且目的层段砂岩与泥岩纵波阻抗值叠置,基于叠后地震的属性分析及纵波阻抗反演技术无法区分岩性和进行烃类检测。在岩石物理分析确定岩性和流体敏感弹性参数的基础上,基于叠前弹性参数反演技术对岩性进行预测,进而通过基于多个弹性参数的FFP(Facies Fluids Probabilities)分析技术对含油气砂岩概率进行定量预测。2项技术在A气田开发中取得了良好应用效果: 通过纵横波速比(v_p/v_s)参数反演,清晰刻画了主力层河道砂体形态,据此指导水平井优化并获得高达95%的高砂岩钻遇率; 通过v_p/v_s和纵波阻抗双参数联合的FFP分析技术预测了含气砂岩空间概率分布,据此预测了潜力气层并指导调整井成功挖潜。这2项技术对于东海油气田砂泥波阻抗叠置区的储层精细描述、油气潜力预测具有良好的推广应用价值。     

Linking dynamic elastic parameters to static state of stress: toward an integrated approach to subsurface stress analysis

Stress is the most important parameter to understand basin dynamics and the evolution of hydrocarbon systems. The state of stress can be quantified by numerical geo-mechanical modelling techniques. These techniques require static elastic parameters of the rocks as input, while tectonic and gravitational forces are given as explicit boundary conditions to compute the local state of stress at different scales. We developed a technique to determine the density and elastic constants at seismic frequencies using full Zoeppritz inversion on angle-dependent seismic reflection data. The dynamic elastic parameters as obtained from seismic data differ from their static equivalents, which are necessary to determine the static state of stress. The dynamic elastic parameters are related to their static equivalents through experimentally obtained relations. In these rock-physics experiments, the static and dynamic elastic parameters are measured simultaneously during different external loading conditions. The experiments used here are all carried out in a tri-axial pressure machine under equal axial stresses. Then pre-stack seismic data analysis in combination with the relation between the static and dynamic elastic parameters, from the rock-physics experiments, provides the input parameters for geo-mechanical modelling.