Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied. First, the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier expanding with respect to azimuth, and the state equation is established by Hankel integral transform method, furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation. Secondly, by the transfer matrixes, the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of ground as well as the contact conditions. Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundaryvalue problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily. At the end of this paper, a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.

Computation of one—dimensional consolidation of double layered ground using differential quadrature method

The authors give the solution to the problem of one-dimensional consolidation of double-layered ground with the use of the differential quadrature meahod.Case studies showed that the computational results for pore-water pressure in soil layer agreed with those of analytical solution;and that in the computationsl re-sults for the interface of soil lwyer also agreed with those of the analytical solution except for the small discrep-ancies during shortly after the start of computation.The advantages of the solution presented in this paper are that compared with the analytical solution,it avoids the cumbersome work in solving the transcendental equation for eigenvalues,and in the case of the Laplace transform solution,it can resolve the precision prob-lem in the numerical solution of long time inverse Laplace transform.Because of the matrix form of the solution in this paper,it is convenient for formulation computational program for engineering practice.The formulas for calculating double-layered ground consolidation may be easily extended to the case of multi-layered soils.

Mode-based equivalent multi-degree-of-freedom system for one-dimensional viscoelastic response analysis of layered soil deposit

Discrete models such as the lumped parameter model and the finite element model are widely used in the solution of soil amplification of earthquakes. However, neither of the models will accurately estimate the natural frequencies of soil deposit, nor simulate a damping of frequency independence. This research develops a new discrete model for onedimensional viscoelastic response analysis of layered soil deposit based on the mode equivalence method. The new discrete model is a one-dimensional equivalent multi-degree-of-freedom（MDOF） system characterized by a series of concentrated masses, springs and dashpots with a special configuration. The dynamic response of the equivalent MDOF system is analytically derived and the physical parameters are formulated in terms of modal properties. The equivalent MDOF system is verified through a comparison of amplification functions with the available theoretical solutions. The appropriate number of degrees of freedom（DOFs） in the equivalent MDOF system is estimated. A comparative study of the equivalent MDOF system with the existing discrete models is performed. It is shown that the proposed equivalent MDOF system can exactly present the natural frequencies and the hysteretic damping of soil deposits and provide more accurate results with fewer DOFs.

Exploring Soil Layers and Water Tables with Ground-Penetrating Radar

Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore its application for mapping in subsurface agricultural soils to a depth of several meters. For a loamy sand and a clayey site on the North China Plain, clay inclusions in the sand were detected; the thickness, inclination, and continuity of the confining clay and silt layers was assessed; and a local water table was mapped. Direct sampling (soil coring and profiling) in the top meter and independent measurement of the water table were utilized to confirm the findings. Also, effective estimates of the dielectric number for the site with the dielectric number of moist clayey soils depending strongly on frequency were obtained. Thus, important properties of soils, like the arrangement and type of layers and in particular their continuity and inclination, could be explored with moderate efforts for rather large areas to help find optimal locations for the time-consuming and expensive measurements which would be necessary to detail a model of the subsurface.

Numerical Simulation of Heat Transfer Characteristics of Horizontal Ground Heat Exchanger in Frozen Soil Layer

A simplified numerical model of heat transfer characteristics of horizontal ground heat exchanger (GHE) in the frozen soil layer is presented and the steady-state distribution of temperature field is simulated. Numerical results show that the frozen depth mainly depends on the soil′s moisture content and ambient temperature. The heat transfer loss of horizontal GHE tends to grow with the increase of the soil′s moisture content and the decrease of ambient temperature. Backfilled materials with optimal thermal conductivity can reduce the thermal loss effectively in the frozen soil. The applicability of the Chinese national standard “Technical Code for Ground Source Heat Pump (GB 50366-2005)” is verified. For a ground source heat pump project, the feasible layout of horizontal GHE should be determined based on the integration of the soil′s structure, backfilled materials, weather data, and economic analysis.

Layer-Constrained Triangulated Irregular Network Algorithm Based on Ground Penetrating Radar Data and Its Application

In this paper,a layer-constrained triangulated irregular network（ LC-TIN） algorithm is proposed for three-dimensional（ 3 D） modelling,and applied to construct a 3 D model for geological disease information based on ground penetrating radar（ GPR） data. Compared with the traditional TIN algorithm,the LCTIN algorithm introduced a layer constraint to the discrete data points during the 3 D modelling process,and it can dynamically construct networks from layer to layer and implement 3 D modelling for arbitrary shapes with high precision. The experimental results validated this method,the proposed algorithm not only can maintain the rationality of triangulation network,but also can obtain a good generation speed. In addition,the algorithm is also introduced to our self-developed 3 D visualization platform,which utilized GPR data to model geological diseases. Therefore the feasibility of the algorithm is verified in the practical application.

Stress dependence of elastic wave dispersion and attenuation in fluid-saturated porous layered media

The fluid-saturated porous layered(FSPL)media widely exist in the Earth's subsurface and their overall mechanical properties,microscopic pore structure and wave propagation characteristics are highly relevant to the in-situ stress.However,the effect of in-situ stress on wave propagation in FSPL media cannot be well explained with the existing theories.To fill this gap,we propose the dynamic equations for FSPL media under the effect of in-situ stress based on the theories of poroacoustoelasticity and anisotropic elasticity.Biot loss mechanism is considered to account for the stress-dependent wave dispersion and attenuation induced by global wave-induced fluid flow.Thomsen's elastic anisotropy parameters are used to represent the anisotropy of the skeleton.A plane-wave analysis is implemented on dynamic equations yields the analytic solutions for fast and slow P waves and two S waves.Modelling results show that the elastic anisotropy parameters significantly determine the stress dependence of wave velocities.Vertical tortuosity and permeability have remarkable effects on fast and slow P-wave velocity curves and the corresponding attenuation peaks but have little effect on S-wave velocity.The difference in velocities of two S waves occurs when the FSPL medium is subjected to horizontal uniaxial stress,and the S wave along the stress direction has a larger velocity,which implies that the additional anisotropy other than that induced by the beddings appears due to horizontal stress.Besides,the predicted velocity results have the reasonable agreement with laboratory measurements.Our equations and results are relevant to a better understanding of wave propagation in deep strata,which provide some new theoretical insights in the rock physics,hydrocarbon exploration and stress detection in deep-strata shale reservoirs.

Measurement of a thin layer's thickness using independent component analysis of ground penetrating radar data

To detect overlapped echoes due to the thin pavement layers,we present a thickness measurement approach for the very thin layer of pavement structures.The term "thin" is relative to the incident wavelength or pulse.By means of independent component analysis of noisy signals received by a single radar sensor,the overlapped echoes can be successfully separated.Once the echoes from the top and bottom side of a thin layer have been separated,the time delay and the layer thickness determination follow immediately.Results of the simulation and real data verify the feasibility of the presented method.

Finite layer and triangular prism element method to subsidence prediction and stress analysis in underground mining

The application of the finite layer & triangular prism element method to the 3D ground subsidence and stress analysis caused by mining is presented. The layer elements and the triangular prism elements have been alternatively used in the numerical simulation system, the displacement pattern, strain matrix, elastic matrix, stiffness matrix, load matrix and the stress matrix of the layer element and triangular prism element have been presented. By means of the Fortran90 programming language, a numerical simulation system based on finite layer & triangular prism element have been built up, and this system is suitable for subsidence prediction and stress analysis of all mining condition and mining methods. Comparing with the infinite element method, this approach dramatically reduces the size of the set of equations that need to be solved, and greatly reduces the amount of data preparation required. It not only saves the internal storage, and the computation time, but also decreases the cost.

我国黄土高原地区地震动衰减关系研究的若干进展

地震动的衰减关系是指地震动随震级、距离和场地条件等变化的经验关系,近年来在我国亦被称为地震动预测方程,是估计地震动及其影响场的主要方法之一,在地震区划和重大工程场地地震安全性评价中被广泛应用。黄土是一种特殊土,在我国广泛分布。我国黄土高原地区地质构造复杂,新构造活动强烈,中强地震频发,地震动的衰减关系有其特殊性,总结我国黄土高原地区地震动衰减关系的研究成果对促进黄土高原地区抗震研究有重要意义。在简要介绍国内外地震动衰减关系研究的基础上,全面系统地总结我国学者在黄土高原地区的地震烈度衰减关系,基岩和土层场地地震动峰值加速度、峰值速度、峰值位移以及反应谱衰减关系方面的研究成果;评述和讨论在黄土覆盖地区地震动衰减研究领域存在的问题和今后的研究方向。文章的研究工作对从事黄土高原地区地震工程研究的科技人员有重要参考价值。

探地雷达在白洋淀湖底地层结构探测中的技术攻关与实践

白洋淀湖底地层调查对于白洋淀生态地质调查、湖底生态清淤及湿地生态保护具有非常重要的现实意义。为系统探明白洋淀湖底地层结构,从区域尺度对湖底地层实现完整刻画,为钻探调查提供精确的靶区数据支撑,创新了探地雷达应用范围,采用船载方式,透过水体开展了湖底地层勘查。剖析白洋淀湖底结构探地雷达勘探技术难点,分析了含底泥和不含底泥2种情形下水上勘探电磁波的传播特征和湖底原生地层的地球物理响应条件,建立了淀水+淤泥层+湖底原生地层的湖底简化结构模型。通过对比国内外主流探地雷达系统不同硬件设备的勘探能力与适应条件,优选出双皮划艇搭载50 MHz低频组合天线的工作模式,结合适用于水下电磁波弱信号提取的数据处理技术,获得了白洋淀湖底15 m深度范围的地层结构数据。结果表明,此种工作模式适用于湖区水深小于5 m的水域,原生地层电磁波整体呈现弱反射特征,湖底15 m深度范围内存在2~3组砂层反射界面。淤泥层会大大削弱湖底原生地层的电磁波反射强度。存在厚砂层的水域,电磁波会形成强反射波组,砂层展布形态刻画清晰,砂层厚度变化在0~3 m之间。本项研究实现了湖底地层结构建模分析、天线勘探模式优选、数据处理过程优化等技术集成。该技术成果为白洋淀淀区地下水与地表水三维地质建模提供了基础依据,为湖底地层钻探勘查提供了靶区支撑,有力推动了白洋淀湿地水-陆一体化探测技术方法体系的构建。

地表基质调查分层及分层测试指标体系设计与构建

地表基质是地球表层物质与能量循环的主要载体,地表基质调查是一项创新性工作,明确地表基质调查分层及分层测试指标是开展地表基质调查的基础。本文基于地表基质调查的内涵和目标定位,深入分析了多目标区域地球化学调查、全国第三次土壤普查和全国农用地土壤污染详查等调查深度和测试指标;利用近三年自然资源综合调查指挥中心在吉林省梨树县、辽宁省丹东市、黑龙江省哈尔滨市三个地表基质调查试点项目实测数据,运用变动系数法分析了不同地区地表基质层的地球化学元素垂向变化特征,结合表生地质作用(风化为主)最大深度、地表植被根系最大深度、浅层地下水位波动下限、基岩顶面深度等因素,将地表基质在垂向上大致分为三层:表层(0~2m)、中层(2~10m)和深层(10~20m),基岩埋深浅于20m则以揭露到基岩为准;表层(0~2m),称为生产层,包括支撑农业生产的耕作层,主要服务国家粮食安全,以收集资料为主;中层(2~10m),称为生态层,主要支撑植被生长和群落演替,服务国土空间开发保护与格局优化,是地表基质调查的重点;深层(10~20m),称为沉积层,主要研究地表基质从基岩—深层—中层—表层演替规律和特征,以部署少量工程控制为宜。基于地表基质分层调查意义,结合元素在垂向上的变化规律,提出4+N型地表基质分层测试指标体系,即4类必测指标和针对特殊工作需要的N类选测指标,以期加快推进地表基质调查标准化、规范化。

地中电偶极源极低频电磁场的分布与传播特征

利用层状大地中偶极源响应的正演算法,计算地中电偶极源激发的极低频地震电磁场并分析其在地球环境下的传播特征.设计了多个水平层状地球模型,分别模拟和展示了深埋地中的电偶极源的响应和空间分布特征.重点分析了含地壳波导+LAI波导模型的高阻大地中,倾斜电偶极源激励响应随观测点的偏移距和垂直位置、激励源深度、地壳波导结构和参数变化时各场量响应的幅频特征.研究表明,用倾斜电偶极源模拟和分析高阻大地中孕震电磁辐射的响应及在大地电磁系统中的传播特性是可行的;模拟的高阻大地中的电磁辐射在地壳波导和LAI波导中均表现有慢衰减或幅值增强特性,但两个波导效应具有不同的频率选择性;高阻大地中的电磁辐射在波导的高阻介质中具有幅值强、衰减慢的特点.建议在高阻地层出露的地表、井中或海底的高阻岩层中以及大气层中布设测站.

基于多层感知机技术的地铁盾构施工参数预测

在地铁工程建设中,盾构法施工技术已经得到了广泛的应用,盾构掘进参数的合理预测对提高施工安全性及降低操作难度具有较大实际意义.以中国杭州机场快线地铁隧道某标段为工程背景,以隧道直径范围内土层摩擦角、黏聚力、压缩模量、重度以及隧道顶部埋深、盾构机预设刀盘转速、推进速度作为输入,以盾构施工时的注浆量、注浆压力、出土量、总推力和刀盘扭矩为输出,建立基于多层感知机的盾构掘进参数预测模型.通过对比不同超参数组合情况下的模型在数据集上的预测表现,挑选出适合于该工程盾构施工参数的预测模型.使用实测数据对模型预测效果进行验证,预测值与实测数据总体变化规律一致,平均误差在20%以内.建立的多层感知机模型预测结果较为合理,具有较好的预测精度,可用于复合地层条件下同类型盾构掘进参数的预测.

基于土壤分层当量物性的地埋管换热器换热模型研究

文章综合考虑了土壤分层与地下水渗流等情况,结合实测数据,通过地埋管外壁温度计算得出地埋管内部流体温度与不同分层土壤的热物性参数,根据土壤当量物性进行分层,据此建立地埋管换热模型。经过对比验证实验发现,分层换热模型比常规均质模型更接近实测值,误差更小;基于土壤分层当量物性的地埋管换热器换热模型模拟出水温度与实测值平均绝对误差为0.21℃,模拟进、出口换热温差与实测值平均绝对误差为0.14℃,具有较高的准确度,可为工程设计提供依据并作为后续地埋管换热优化研究的基础。

水平多层spiral型地埋管数值模拟研究及换热特性分析

本文基于COMSOL软件建立了水平spiral型多层地埋管群传热三维数值模拟模型,并对其换热性能和温度场进行分析研究。结果表明:对比水平spiral多层并联管群来说,远离注入端的地埋管管内温度相对较低。地埋管群周围地层温度有所升高,远离注入端的地层温度变化相对较小。地埋管群的换热量随循环流体的流速增加而增加,但增幅降低。在进行换热器设计时,应合理控制循环水的流速。3层螺旋地埋管群的换热量略高,但层间干扰明显,建议选用2层并联螺旋管群。

基于G-TF/SF方法的分层地表电磁环境研究

本文将基于广义总场散射场(generalized total-field/scattered-field,G-TF/SF)方法的时域有限差分(finite-difference time-domain,FDTD)方法用于分层地表的电磁环境计算中,四周采用卷积完全匹配层(convolutional perfectly matched layer,CPML)吸收边界条件。G-TF/SF方法可有效削减由地面四周截断产生的边缘效应,提高分层地表电磁环境计算精度且G-TF/SF边界只需要加载入射场,无需计算反射场与透射场。G-TF/SF方法将部分边界埋在完全匹配层(perfectly matched layer,PML)区域中,对埋在PML层中G-TF/SF边界的入射波乘以相应的系数,保证将合适的入射波等效到主网格区域进行计算。本文用G-TF/SF方法模拟计算了水平单层地表及不同分层地表情形下的电磁环境,并将水平单层地表结果与用解析法计算的结果进行对比;将不同分层情形与采用有限积分技术得到的结果对比,验证了G-TF/SF方法计算结果的正确性。