抽水引致的地面沉降的参数敏感性分析

目前对地面沉降的研究和计算中,由于地质情况复杂、土体的变异性不容易掌握,对于土体参数的选取多数是根据经验值,具有不确定性。特别是在大面积的地面沉降计算中,这种不确定性更为严重。采用数值分析法,以Biot固结理论为基础,应用有限元软件Adina,建立了单井抽水的地面沉降模型,分析各参数在抽水过程中对地面沉降的影响,从而确定敏感参数,以期对地面沉降研究中的参数选取提供有价值的指导。     

华北平原鲁北地区地下水超采导致地面沉降区域特征及演化趋势预测

鲁北地区作为华北平原地面沉降的重要组成部分,其地面沉降问题日趋严重。以滨州博兴县为工程背景,基于研究区详细水文地质与工程地质资料以及历年地面沉降监测数据,系统分析该地区地下水动态分布及地面沉降分布演化特征。以Biot多孔介质固结理论为基础,建立博兴县地面沉降三维流 固耦合数值模型,还原地面沉降发展过程并预测分析不同地下水开采方案下的沉降演化规律。研究结果表明：博兴县浅层地下水位降幅呈现南大北小特点,深层地下水形成了以县城区为中心的椭圆形地下水区域降落漏斗;地面逐渐形成了分别以博兴县城区、湖滨镇和店子镇为沉降中心的三个小型沉降区,且有相互关联扩展的趋势;地面沉降三维流固耦合模型较为理想还原了研究区地面沉降发展过程,预测在现状地下水开采方案下未来10年内地面沉降仍以较大速率继续发展,累计沉降量超500mm的区域面积不断扩大,当减小20%现状地下水开采量时是较为合理有效的开采方案。     

面向对象方法在Biot固结有限元程序中的应用

面向对象程序设计方法应用于有限元程序编制带来了相当大的优越性。但迄今为止，这种方法在岩土工程中的应用实例很少。至于用面向对象方法进行土力学中常见的Biot固结有限元分析软件的开发，尚未见报道。现将面向对象有限元方法引入Biot固结有限元程序编制中，用C＋＋语言编制中一个Biot固结有限元程序实例。和结构化程序设计方法相比较，采用面向对象方法进行Biot固结有限元程序的编制，使程序的可维护性、可再生性、可扩充性和程序开发的效率都得到了很大的提高。指出将面向对象程序设计方法引入了Biot固结有限元程序的编制等岩土工程的数值计算领域，是今后的发展方向。     

西安市地面沉降机制

在西安市地面沉降基本特征分析的基础上,初步探讨地面沉降带下伏土层的固结性能,进而分析西安市地面沉降机制,初步研究表明,由于土层的超固结性能,西安市地面沉降明显不同于沿海城市（如上海,天津等）的地面沉降。     

Biot固结理论中连续性方程形式的讨论

对于Biot固结方程中的连续性方程形式,不同的研究者存在分歧。在详细分析了Biot推导固结方程的过程后,发现其在推导固结方程时用到的弹性应变能密度函数存在错误。为了得到正确的连续性方程,采用两种思路进行了分析：①通过对土体单元变形的分析,得到了考虑土体颗粒和流体可压缩性的连续性方程;②按照Biot推导固结方程的思路,通过修正推导过程中用到的弹性应变能密度函数得到了连续性方程。上述两种思路得到的连续性方程是等价的,同时将上述连续性方程和其他研究者通过质量守恒得到的连续方程进行了对比,结果是一致的,从而确定了连续性方程的正确形式,并澄清了目前在Biot固结方程中的连续性方程上存在的分歧,所得结论可为固结方程的解析或数值计算提供依据。     

地下水开采—地面沉降模型研究

本文建立三维流支－－一维非线性固结地面沉降模型,在混合井流,降雨滞后补给,初始水头形成,人为边界刻画,水流－固结耦合及软土层固结滞后于地下水开采层水头变化等方面有所改进,所建模型用于苏州市,模拟出若干地面沉降重要特征,地面沉降中心偏离地下水漏斗中心,地面沉降动态滞后于地下水的水头动态,软土层渗透系数随固结过程的变化及头在软土层中的传递特征。     

横观各向同性圆柱土样轴对称Biot固结的解析解

给出了横观各向同性圆柱土样轴对称Biot固结的封闭级数解,该解由满足非齐次边界条件的弹性力学解与满足齐次边界条件的渗流拟动态解的叠加构成。通过算例分析,探讨了横观各向土体的Biot固结规律和Mandel-Cryer效应。     

土体自重固结压缩对地面沉降影响研究

以南通市为例,建立了土体自重固结压缩、地下水开采与地面沉降三维全耦合模型.在对模型进行识别、验证的基础上,模拟预测了南通市在地下水停止开采,仅在土体自重固结压缩影响下,自2010年12月31日-2025年12月31日各含水层组地下水流场变化和地面沉降的发展趋势.结果表明,在此期间,南通市由土体自重固结压缩引起的最大地面沉降量为2.42 mm,最大地面沉降速率为0.16 mm/a,土体白重固结压缩对南通市地面沉降的影响极为有限.     

多层地基非轴对称Biot固结的理论解

提出了一种新方法来推导多层地基非轴对称Biot固结问题的传递矩阵,相应的理论推导工作量较少。基于柱坐标系下非轴对称Biot固结的基本方程,通过引入中间变量,并对坐标 进行Fourier级数展开,对时间 和坐标 进行Laplace-Hankel变换,得到了6×6阶和2×2阶的两组常微分方程;然后,两组常微分方程进行关于 的Laplace变换和逆变换,得到了单层地基非轴对称Biot固结问题的传递矩阵;结合边界条件和连续条件,运用传递矩阵法得到多层地基非轴对称Biot固结问题在Laplace-Hankel变换域内的解,通过Laplace-Hankel逆变换得到了该问题物理域的真实解答。编制了计算程序,并进行了数值计算与分析。该方法具有计算速度较快的特点,方便工程应用。     

超固结土层垂向形变机理及临界水位控沉分析——以西安南郊钻孔土样研究为例

超固结地层发生地面沉降需要一个大于超固结应力的起始附加应力来驱动。依据有效应力原则,临界水位与超固结应力之间存在线性相关关系,则超固结地层垂向形变机理亦可通过临界水位描述为：当人为开采地下水引起的起始水位降深大于临界水位时,地面沉降将明显发生。西安南郊地质技校钻孔埋深100-300m内土样试验资料分析及计算结果表明：该地层普遍呈现超固结特性;第一、二承压含水组的临界水位值分别为63.44m及73.76m;将“临界水位”作为控制地面沉降与最大限度开采地下水资源的主控指标来考虑,是地面沉降的有效防治措施。     

Geochemistry of Thermal Waters of Continental Margin of Far East of Russia

Studied waters belong to warm(T=30-50℃),alkaline(pH=8.9-9.3),low mineralized(TDS235 mg/1)Na-HCO_3 or Na-SO_4-HCO_3 thermal waters with high content of SiO_2(up to 81 mg/l)and F(up to 3.9 mg/1),occur on modern volcano-tectonic rejuvenated areas of Eastern Sikhote-Alin orogenic belt.Low~3He concentration as well as N_2/O_2 and N_2/Ar ratios exclude influence of deep mantle fluid.New rare earth element data constrain our understandmg of water-rock interaction occurring in the water source region.Meteoric origin of waters is proved by stable isotope values varying from-71‰to-136.1‰and from-10.8‰to-18.8‰forδ~2U andδ~(18)O respectively.REE patterns reflect high pH,resultfing from water-rock interaction and oxidative conditions.Calculations of deep aquifer temperature using Na-K and quartz geothermometers show 116.8-131.1°C and 82.2-125.8℃respectively.Presence of deep faults both with abnormal thermal gradient(~45-50 K/km)define unique geochemical shape of thermal waters of Sikhote-Alin,area,where no present volcanic activity is registered.     

黄土路基沉降曲线特征分析

在分析郑西客运专线试验段路基沉降实测资料的基础上,提出沉降变形常规预测方法受观测时间及单次观测结果的影响比较大,不利于用来预测最终沉降。通过分析现场沉降变形与时间关系曲线,提出路基沉降的发展主要受固结的影响,其发展是一个从沉降增加到趋缓的过程,可以明显分为固结变形及流变变形,并用非饱和土等效固结理论加以解释,提出了一种利用s- 曲线特征预测最终沉降量的方法,经检验对比发现其预测结果受观测时间及单次观测结果的影响比较小,优于其它常规预测方法     

沉管隧道基槽爆破施工对既有堤岸稳定性影响的数值仿真分析

结合爆破数值模拟中爆破荷载的特点,分析了各种荷载施加方法的特点及爆破震动场模拟的要求,基于圣维南荷载等效原理,提出了爆破荷载的等效施加方法,即将作用在炮孔壁上的爆破荷载等效后,施加在同排炮孔中心线或者面上来模拟爆破荷载的作用。结果表明,等效施加方法在近区存在一定差异,但在中远区吻合得较好。利用该等效荷载施加方法,采用三维实体模型计算分析了某沉管隧道基槽爆破施工对既有堤岸稳定性影响,计算了周围不同位置、不同深度位置处的振动响应,并在此基础上,分析研究了减小爆破振动对既有结构物影响的减振措施,可为类似工程提供借鉴。     

Structure of sediments of kaolinite

In an assembly of clay particles placed in a fluid, each particle is typically subjected to: (1) double-layer repulsive forces; (2) van der Waals attractive forces; and (3) contact mechanical forces. The study presented here concerns an approximate, quantitative analysis of clay suspensions, with considerations to the first two - the physico-chemical forces. Using recent theories to calculate the physico-chemical forces between two clay particles in an approximate model of an assembly, the equilibrium void ratio of a clay suspension is computed. The mechanical forces are ignored in the analysis. The results serve to verify the validity of physico-chemical theories employed and help interpret experimental data more fundamentally in terms of the system variables.     

岩石坚硬程度对核电厂地基-基础地震响应特征的影响分析

考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明：随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期（高频段）方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。     

Strains of scattering of near-field of a point source

Three dimensional scattering of near-field is studied for dilatation and rotation in the time domain. The perturbation method is applied to solve the equation of motion for the first order scattering from a weak inhomogeneity in an otherwise homogeneous medium. The inhomogeneity is assumed close enough to the point source so that the near-field intermediate wave is dominating over the far-field sphericalP andS pulses. The integral expressions are derived to relate dilatation and rotation of scattering to the radial fluctuations of velocities and density in the inhomogeneity. These integrals are solved to calculate the strains of scattering from (a part of) an inhomogeneous spherical shell of arbitrary curvature. Variable curvature may allow the shape of inhomogeneity volume element to change uniformly from spherical to rectangular. Rotation of scattering from a spherical shell is independent ofP wave velocity inhomogeneity. Dilatation of scattering does not involveS wave velocity inhomogeneity but its gradient. The back scattering results are obtained as a special case. Strains are computed numerically, for hypothetical models to study the effects of various parameters viz., velocity inhomogeneity, distance of source from inhomogeneity and from receiver, and thickness of inhomogeneity. The curvature of the spherical shell is varied to study the effects of the shape of inhomogeneous volume element on scattering.     

不同形式诱导缝的等效强度承载性能的数值分析

将考虑尺寸效应规律的混凝土等效断裂模型应用到不同形式诱导缝等效强度公式中,得到了修正后的不同形式诱导缝的等效强度公式。同时,得到了诱导缝等效强度的变化曲线,并比较了单双向间隔诱导缝的优劣,得到双向间隔诱导缝优于单向间隔诱导缝的结论,为高拱坝中诱导缝的设计提出了建议。     

土石坝心墙水力劈裂分析

简要分析了水力劈裂的发生条件、力学机理。基于水压楔劈机理,利用有限元方法建立了一种水力劈裂的发生判定方法。该方法假定心墙预先存在局部渗透弱面（裂缝）,通过将裂缝位置的单元材料改为裂缝软材料,考虑库水进入裂缝后对裂缝周围土体的作用,建立水力劈裂分析的平面有限元模型,确定裂缝端部垂直于裂缝面的正应力,进而依据该正应力判断水力劈裂发生的可能性,该方法同时可模拟水力劈裂的发展过程。     

土的可恢复剪胀的一种解释

从土的各向异怀角度对土的可恢复剪胀现象进行了解释。基于各向异性情况下的土体弹性本构关系理论分析，认为土的可恢复剪胀现象可部分归因于土的各向异性引起的弹性剪胀。借助有关土体弹性参数实验结果，研究了应力诱导各迥异性对土体弹性剪胀的影响，结果表明：随土体应力诱导各向异性的增大，土体的弹性剪胀也增大。从土体弹性剪胀角度研究了土的卸荷体缩条件，认为土体卸荷体缩取决于加载应力路径的应力增量比，给出了土体出现卸荷体缩的区域。