条带矿柱平均垂直应力的宏观分布

本文选用层状介质理论及其相应的开采沉陷预计分析系统，通过反演给出了条带开采时条带矿柱平均垂直应力宏观分布的图式及其简化计算方法。     

两种计算瑞利波理论频散曲线的传递法模拟对比

对目前计算层状介质中瑞利导波理论频散曲线的2类传递矩阵方法进行总结,分别利用2类传递矩阵方法对假定的层状介质模型进行模拟计算,比较了2类方法的特点。通过模拟计算结果可知,2类传递矩阵算法是一致的,但快速标量传递算法的计算速度比Abo-zena,Menke和Bixing Zhang等推导的传递矩阵算法快1倍多,并且其计算的稳定性和计算精度要优于后者,但在避免高频数值溢出方面,Bixing Zhang等推导的传递矩阵算法比快速标量传递算法要好。     

基于球体层状介质模型的大地电磁正演

目前大地电磁测深法的一维正演理论,是基于平面波垂直入射水平层状介质的假设模型。但由于地球是一个球体,因此有必要研究基于球状介质模型的大地电磁正演理论。这里详细推导了基于球体层状介质模型的大地电磁正演公式,计算了若干理论模型。通过同基于水平层状介质模型的大地电磁正演结果对比,验证了正演公式的正确性。同时,指出当探测周期增加到上万秒时,阻抗相位会增大;而当探测周期增加到数十万秒时,视电阻率会减小。     

利用连分式定义瞬变电磁法全区视电阻率研究

介绍了一种新的全区视电阻率计算方法。该方法利用连分式迭代求解非线性方程技术,直接对均匀半空间电偶源瞬变电磁法观测的垂直磁场与电阻率的非线性方程直接求解。利用该数值方法对层状介质进行了理论计算。数值计算结果表明,采用连分式方法可直接求出反映地电断面结构的瞬变电磁法全区视电阻率。      

磁偶源频率测深视电阻率曲线的正演计算

磁偶源频率测深法是一种有效的地球物理勘探方法。多年来,人们对这种方法的正演理论和反演解释进行了研究。文中利用快速汉克尔数值滤波算法计算了层状介质上直立磁偶极子的电磁场响应以及频率测深视电阻率的振幅响应,这种计算方法与常用的数值积分算法相比,具有计算速度快、精度高等特点。同时详细地讨论了磁偶源频率测深的特点及视电阻率理论曲线的特征。     

磁偶源2.5维瞬变电磁场全空间FDTD数值模拟

为解决模拟瞬变电磁场在3维介质中传播出现的计算量大、对计算机要求高等问题,使用了3维磁偶极源来模拟2维均匀介质和层状介质中方形和薄板低阻体的全空间响应特征。模拟中使用了时域有限差分法,采用非均匀网格和时间步长,分析了瞬变电磁场在均匀介质和层状介质中的传播规律,以及巷道和交界面对场的影响特征,为解释全空间均匀介质和层状介质TEM异常提供参考。     

层状介质背景下地-井方位激电异常分析与解释

为了更好地探究层状介质背景下井旁隐伏矿体的地-井方位激电异常特征,首先研究了地-井三维激电的有限元正演模拟算法,并编制了程序;然后构建三维地电模型并进行正演模拟,讨论层状介质背景下地-井方位激电异常,并分析其异常特征规律;最后通过比值法归一化压制层状介质背景,突出局部异常。理论和算例表明,比值法归一化处理是地-井方位激电观测中消除层状介质背景异常影响的有效方法。     

不均匀介质中定深度的方法技术

利用定深度电法理论对水平层状介质找矿方法，水平方向不均匀介质中的找矿方法技术及用定深度装置研究矿体产状和形态进行了理论和方法的研究，从电阻网络模拟成果，野外方法的实践等对该方法作了进一步的验证和补充。     

三套国外TEM软件正反演结果对比

在对三套国外瞬变电磁软件Maxwell、EM vision和Emigma研究分析的基础上,设计均匀半空间、层状介质和均匀半空间加载板状体及棱柱体的模型,逐一比较分析三套软件的异同,并进行了均匀半空间和层状介质的1D反演,通过模拟计算和结果分析,比较不同条件下各软件的特点和结算结果的准确性。结果表明,在均匀半空间和层状介质条件下,三套软件的正演计算结果能完全吻合,但瞬变电磁3D正演问题还未从根本上得以解决。     

海洋沉积物的合成地震记录以及孔隙度和渗透率的确定

本文提出了海洋沉积物垂直地震剖面（VSP）的数值模拟，在多孔弹性水平层状介质中，用垂直入射波计算所得的理论地震记录包括有用毕奥理论预测的频散和衰减效应，根据毕奥理论激发出快纵波和慢纵波而且在界面上存在波型转换，由于在整个慢纵波传播过程中存在能量损失机制，所以在反射和透射系数的计算中，我们考虑了这种效应，最后，我们探讨了从合成地震记录资料中求取孔隙度和渗透率的谱比率法，根据弱框架近似法，求出了速度和单位衰减的解析表达式，一旦计算出与频率有关的速度和单位衰减，就可以用提供的弱框架近似法计算出海洋沉积物的孔隙度和渗透率，关于合成模拟的谱比率计算表明；与VSP法相联系的毕奥理论能够用用来确定海洋沉积物的孔隙度和渗透率。     

存在流体介质夹层时瑞利波频散曲线的完整求解

应用介质分界面上位移应力连续性的变化关系,采用δ矩阵法和固体层状介质中的瑞利波频散方程,推导了含流体介质夹层情况下瑞利波的频散方程。利用推导的频散方程对理论模型进行模拟计算,并与以往文献计算方法进行比较,模拟结果吻合较好,验证了公式的正确性;同时,模拟计算结果也说明了本文公式的普遍适用性和计算结果的完整性。根据文中对4个理论模型的模拟结果来看,与固体层状介质的情况相比,当存在流体夹层时,瑞利波频散曲线模式中多出了一个文中定义的0阶模式,在今后的多道瞬态瑞利波法应用中,应该可以将其看成是流体介质夹层存在的标志性曲线。     

计算层状介质中轴对称柱面瑞利面波频散函数的δ矩阵法

笔者将δ矩阵法应用于计算层状介质中轴对称柱面瑞利面波的频散函数,得到了六阶δ矩阵法、五阶δ矩阵法、快速δ矩阵法3种方法,很好地解决了高频数值精度丢失问题以及高频数值溢出问题,并提高了计算速度,数值计算及工程应用验证了上述方法的有效性,且表明了这些方法也完全适用于平面瑞利面波频散曲线的求取问题。     

Analytical layer element solutions for deformations of transversely isotropic multilayered elastic media under nonaxisymmetric loading

This paper presents the analytical layer element solutions for deformations of transversely isotropic elastic media subjected to nonaxisymmetric loading at an arbitrary depth. The state vectors for the nonaxisymmetric problem are deduced through the substitution of the Hu Hai‐chang solutions into the basic equations for the transversely isotropic elastic media. From the state vectors, the analytical layer element of a single layer is obtained in the Hankel transformed domain. The analytical layer element is an exact and symmetric stiffness matrix whose elements are without positive exponential functions, which can not only simplify the calculation but also improve the stability of computation. On the basis of the continuity conditions between adjacent layers, the global stiffness matrix is obtained by assembling the interrelated layer elements. The solutions for the multilayered elastic media in the transformed domain are obtained by solving the algebraic equation of the global stiffness matrix, which satisfies the boundary conditions. The actual solutions in the physical domain are further obtained by inverting the Hankel transform. Finally, some cases are analyzed to verify the solutions and evaluate the influences of the transversely isotropic character and stratified character of the media on the load–displacement responses. The numerical results show that the variations of the elastic properties between layers have a great effect on the displacements of the multilayered media. Copyright © 2014 John Wiley & Sons, Ltd.     

Transient dynamic response of multilayered saturated media subjected to impulsive loadings

This paper presents a stable and efficient method for calculating the transient solution of layered saturated media subjected to impulsive loadings by means of the analytical layer element method. Starting with the field equations based on Biot's linear theory for porous, fluid‐saturated media, and the seepage continuity equation, an analytical layer element for a single layer is established by applying Laplace‐Hankel integral transform. The global stiffness matrix in the transform domain for a layered saturated half‐space subjected to a transient circular patch loading is obtained by assembling the layer elements of each layer. The displacements in the time domain are derived by Laplace‐Hankel inverse transform of the global stiffness matrix. Numerical examples are conducted to verify the accuracy of the method and to demonstrate the influences of type of transient loading, buried depth of loading, permeability, and stratification of materials on the transient response of the multilayered saturated poroelastic media.     

A macroscopic elastoplastic constitutive law for multilayered media: Application to reinforced earth material

This paper refers to the formulation of a macroconstitutive law for elastoplastic multilayered media, using both thermodynamical approach and homogenization procedure. The latter shows that both the microstress and strain tensors are constant in each constituent. This result permits the formulation of a macroconstitutive law involving the microplastic strains as internal parameters. Both the macrofree energy and dual dissipation potential are computed and used to derive explicit expressions for the macrocompliance tensor, yield surface and hardening rule. Finally, this law is used to analyse the behavior of a reinforced earth material under various loading paths.     

考虑参数空间变异性多层土坡系统可靠度分析

多层土坡在岩土工程实际中十分常见,不仅土体参数存在一定的空间变异性,而且土体框架呈现明显的层状分布特征,然而目前对考虑土体参数空间变异性的多层土坡稳定可靠度研究的远远不够。提出了基于多重响应面边坡系统可靠度分析的蒙特卡洛模拟（MCS）方法,给出了计算流程图,系统地研究了考虑土体参数空间变异性的多层土坡系统可靠度问题。结果表明,提出方法能够有效地分析考虑参数空间变异性低失效概率水平的多层土坡系统可靠度问题,并且具有较高的参数敏感性分析计算效率。     

Response of multilayered transversely isotropic medium due to axisymmetric loads

A novel procedure associated with the precise integration method (PIM) and the technique of dual vector is proposed to effectively calculate the magnitude and distribution of deformations in a homogeneous multilayered transversely isotropic medium. The planes of transverse isotropy are assumed to be parallel to the horizontal surface of the soil system. The linearly elastic medium is subjected to four types of vertically acting axisymmetric loads prescribed either at the external surface or in the interior of the soil medium. There are no limits for the thicknesses and number of soil layers to be considered. By virtue of the governing equations of motion and the constitutive equations of the transversely isotropic elastic body, and based on the Hankel integral transform and a dual vector formulation in a cylindrical coordinate system, the partial differential motion equations can be converted into first‐order ordinary differential matrix equations. Applying the approach of PIM, it is convenient to obtain the solutions of ordinary differential matrix equations for the continuously homogeneous multilayered transversely isotropic elastic soil in the transformed domain. The PIM is a highly accurate algorithm to solve the sets of first‐order ordinary differential equations, which can ensure to achieve any desired accuracy of the solutions. What is more, all calculations are based on the standard method with the corresponding algebraic operations. Computational efforts can be reduced to a great extent. Finally, numerical examples are provided to illustrate the accuracy and effectiveness of the proposed approach. Some more cases are analyzed to evaluate the influences of the elastic parameters of the transversely isotropic media on the load‐displacement responses. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of maximum saturated depth in two-layered drainage blankets over the barrier in landfill cover system

For accurately calculating the maximum saturated depth over the landfill barrier and analyzing the variation of liquid head with inflow time, a new method is presented that considers the unsteady-state condition. Based on the principle of water balance and the extended Dupuit assumption, the solution of the new method can be calculated. Also, this method is capable of revealing the effect of the multilayered drainage media. One of the advantages of the method presented in this paper is that it can be used to draw a phreatic surface in layered drainage media under unsteady-state condition and estimate the location of the maximum liquid depth, for both homogeneous conditions and two-layered drainage conditions. Based on the developed analysis method, the effects of a parametric variation on the maximum saturated depth and the variation of saturated depth with time are presented. The calculated results show that with the increase of slope gradients of cover system, the saturated depth of cover system is easy to reach a stable value. Moreover, when the drainage layer consists of two layers with different thickness and hydraulic conductivity, the increase of hydraulic conductivity of the upper drainage layer will make it easier for the saturated depth of cover system to reach a stable and lower value.     

Axisymmetric thermal consolidation of multilayered porous thermoelastic media due to a heat source

This paper presents an analytical layer element solution to axisymmetric thermal consolidation of multilayered porous thermoelastic media containing a deep buried heat source. By applying the Laplace–Hankel transform to the state variables involved in the basic governing equations of porous thermoelasticity, the analytical layer elements that describe the relationship between the transformed generalized stresses and displacements of a finite layer and a half‐space are derived. The global stiffness matrix equation is obtained by assembling the interrelated layer elements, and the real solutions in the physical domain are achieved by numerical inversion of the Laplace–Hankel transform after obtaining the solutions in the transformed domain. Finally, numerical calculations are performed to demonstrate the accuracy of this method and to investigate the influence of heat source's types, layering, and the porous thermoelastic material parameters on thermal consolidation behavior. Copyright © 2015 John Wiley & Sons, Ltd.     

地埋点热源作用下层状半空间的热-力耦合响应

采用解析层元法对存在地下点热源的岩土工程问题进行解答。首先从热弹性力学三维问题的基本控制方程出发,利用拉普拉斯-傅里叶积分变换推导出其在变换域内单层介质及下卧半空间的解析层元;然后结合有限单元法原理组装得到总刚度矩阵,结合边界条件,得到其在变换域内的解答,最后应用拉普拉斯-傅里叶积分逆变换技术,得到物理域内的解。编制了相应的计算程序,算例分析表明：该结果与已有文献吻合较好,该方法求解层状半空间的热-力耦合响应问题具有较好的适用性和较高的精度;层状岩土介质体系中,热扩散系数对温度及地表隆起的变化过程影响显著,但对其初始值和最终稳定值影响不明显;分层特性对岩土介质温度分布及地表位移变化过程均有显著影响。