立井冻结法凿井工程中的热工计算

本在分析有关实测资料的基础上，应用近似方法求解了二类边界条件下有相变的平面轴对称的热传导方程，给出了人工冻结凿井设计中用以确定冻结壁交圈时间的相邻两冻结管周围冻结土体轴向相交接的时间，冻结壁形成变化规律，冻结壁平均温度和温度和内外冻结壁厚度关系。将计算结果与实测资料比较表明，中所给计算方法和公式有较高的精度。     

基于有限体积法的冻土水热耦合程序开发及验证

土体冻结和融化时的水分迁移、相变与传热是一个相互影响的耦合过程。采用基于有限体积法的开源软件OpenFOAM,编制描述土体冻融过程的水热耦合计算程序。首先,基于土体水分和热量迁移基本方程、水分相变与温度的平衡方程,同时考虑相变对水分特征参数和热特性参数的影响以及相变潜热对传热过程的影响,建立冻土水热耦合数学模型。然后,采用基于多面体网格的有限体积方法对水热耦合控制方程进行空间离散,采用全隐式向后差分方法对方程进行时间离散,由此编制冻土水热耦合计算程序。该程序具有良好的几何适应性、质量和能量守恒性,具备面向复杂问题的并行计算功能。最后,采用该程序对两组不同温度边界条件的室内土体冻结试验进行数值模拟,并与试验结果进行对比,结果表明该程序可以较为准确地模拟土体冻结过程中温度场和水分场的演化特征。     

基于有限体积法的冻土水热耦合程序开发及验证

土体冻结和融化时的水分迁移、相变与传热是一个相互影响的耦合过程。本文采用基于有限体积法的开源软件OpenFOAM,编制描述土体冻融过程的水热耦合计算程序。首先,基于土体水分和热量迁移基本方程、水分相变与温度的平衡方程,同时考虑相变对水分特征参数和热特性参数的影响,以及相变潜热对传热过程的影响,建立冻土水热耦合数学模型。然后,采用基于多面体网格的有限体积方法对水热耦合控制方程进行空间离散,采用全隐式向后差分方法对方程进行时间离散,由此编制冻土水热耦合计算程序。该程序具有良好的几何适应性、质量和能量守恒性,具备面向复杂问题的并行计算功能。最后,采用该程序对两组不同温度边界条件的室内土体冻结试验进行数值模拟,并与实验结果进行对比,结果表明该程序可以较为准确地模拟土体冻结过程中温度场和水分场的演化特征。     

浅埋集输油管线拟稳态温度场及热工计算

从二维相变热传导方程出发，应用保角变换将研究区域规整化，在热力等效原理基础上导出冬季冻结期输油管拟稳态温度公式，并由此给出保温隔热层厚度，管线散热量，管道轴向温降，管道上部最大冻深，冻结速率，系统热阻热工公式。     

饱和正冻土中的水、热、力场耦合模型

土体在冻结过程中 ,不论是融土区、过渡区 (指正在形成的分凝冰及冻结缘区 )还是冻土区所涉及的问题不仅是温度场、水分场问题 ,而且涉及到力学场问题 ,并且力学场对土体的变形过程 (冻胀、压密 )及分凝冰的形成起着重要的作用 .依据连续介质力学、热力学原理 ,提出了土体冻结过程中的三场耦合方程 .并指出 ,土体冻结过程中的体积变化与应力状态、补 -排水条件、冻结条件密切相关 ;体变包含弹性部分、孔隙水变化引起的固结或膨胀部分、相变引起的体变及冻土区粘塑变形部分 ;水分驱动力控制着水分迁移量 ,它是由土基质势、溶质势、冰基质势、孔隙水压力及重力势组成 ;冻结缘现象是由于该区各点处的水分驱动力与相变温度同时满足相变条件而形成 ,冻结缘区含冰量的多少取决于冷能供给情况 ;分凝冰的形成及发展是引起较大冻胀的主要来源 ;分凝冰的形成条件由温度场及力学场统一控制 .即温度要达到相平衡条件 ;孔隙水应力除用于抵抗土骨架有效应力外 ,还需克服土颗粒间的粘聚力     

寒区圆形截面隧道温度场的解析解

根据冻土地区的实际情况地圆形隧道正冻区和未冻区的热传导方程进行简化,应用无量纲量和摄动技术对简化方程进行求解,给出了圆形隧道冻结过程的近似解析解。通过和有限元数值解结果比较,发现该近似解具有较高的精度,能满足工程实际的精度要求,该解既可用于工程估算 用来检验数值计算的计算程序及结果。     

饱和软黏土中埋地管道冻结模型试验研究

为研究软土地区埋地管道在土体冻结过程中的管道受力机理, 开展了饱和软黏土中地埋管道冻结模型试验。通过人工冻结技术, 近似还原了管土受冻过程, 研究了人工冻结过程中土体温度场、 水分场、 位移场分布情况, 以及管道的力学特性。结果表明： 在冻结过程中, 土体温度场的变化直接影响着土中水分场的分布; 在冻结锋面前缘存在着剧烈的水分迁移现象, 大量的水分向冻结锋面迁移, 使得土体产生线性冻胀; 冻胀发展速率受外部荷载的直接影响; 当冻结发展到管道处时, 位于冻胀和非冻胀过渡段位置处的管身出现应力最大值。研究结果对于正冻土中管道的安全评估具有重要的意义。     

广州某地铁人工冻结法施工热力分析

以广州某地铁采用人工冻结法施工为例,考虑相变和衬砌水泥水化热的生成,根据带相变瞬态温度场的热量控制微分方程,采用伽辽金法推导出的有限元计算公式计算分析了此地铁人工冻结施工过程中的温度场分布规律.并根据温度场的分布特点,确定了开挖时的人工冻结壁厚度.最后,对冻结壁和其周围土体分别应用不同的力学本构模型进行施工过程的力学模拟.结果表明,采用人工冻结法作为一种辅助的施工方法,不但能满足施工过程中结构的强度要求,而且能满足环境沉降量的要求.     

基于SPH方法的非饱和多孔介质相变耦合研究

在考虑相变的热能平衡方程和非饱和水分迁移质量控制方程的基础上,建立温度场-水分场的耦合模型,并采用一种无网格粒子算法（SPH）进行数值求解。其中,耦合方程中考虑了水流传热以及温度势对水流的直接驱动,在不考虑相变的情况下,该耦合模型可退化为常温下的水-热耦合模型,故可用于模拟冻融循环的相关问题。从求解热能平衡方程中的含冰量出发,实现解耦并对半无限单向冻结条件下介质内非稳态温度场和体积含水率分布场进行模拟,将耦合作用下的温度场与不耦合的解析解进行对比,反映出水分迁移对温度场存在较大影响。最后,求解了路基边坡在季节性周期温度边界下,温度场、水分场分布的演变规律,并评估了边坡阴阳面受热不均对水热两场分布的影响。计算结果基本能反映土冻结相变的实际物理过程,光滑粒子算法可以用于尝试解决冻土领域的其他相关问题。     

双圈管冻结壁温度场形成特性及影响因素

为研究双圈管冻结壁温度场变化规律，以淮南某矿副井为研究对象，根据其相关地质参数，利用FLAC3D软件数值模拟双圈管冻结壁温度场形成过程以及不同因素对冻结壁平均温度的影响。研究结果表明:双圈管内土体温度最低，两侧温度逐渐升高；冻结孔间距越小，交圈时间越早，内外圈管交圈之后形成封闭的未冻承压水仓，对冻结壁不利，冻结锋面向内侧扩展速度大于向外侧扩展速度；双圈管冻结壁平均温度与冻结时间呈对数关系下降，有效厚度在内外圈管交圈后增长十分明显，且与冻结时间呈对数关系上升；双圈管主、界面温度场曲线随冻结时间近似由马鞍形分布逐渐转变为梯形分布，界面温度场扩展速度大于主面温度场；土体初始温度、盐水温度及导热系数对冻结壁平均温度影响均较大，土体初始温度和盐水温度越低、导热系数越大，冻结壁平均温度越低。研究成果为相关冻结工程的设计和施工提供参考。      

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

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

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

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

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.     

土石坝心墙水力劈裂分析

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

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.     

土钉支护失稳分析

某基坑采用复合土钉墙支护,开挖至坑底后,基坑北侧发生失稳事故。本文就其事故发生原因作一探讨。     

地震次生地质灾害风险评估集对态势分析

5．12四川汶川8级强震引发了大量次生地质灾害,加剧了灾情并严重影响抗震救灾。地震次生地质灾害受多重因素制约,其孕育与致灾过程具有不确定性。利用集对分析原理与方法,对地震次生地质灾害进行同异反态势分析,拟定了风险评估的集对分析同一度、差异度、对立度等指标体系的构建原则与赋值标准,对承灾体系统的不确定性及其作用作了刻画与分析,建立了不同风险分区代表性的集对分析联系度表达式,为地震灾害研究提供了可资借鉴的新思路和技术方法。     

下扬子区三叠纪古地理演化

下扬子区在三叠纪期间接受了早-中三叠世海相碳酸盐沉积、中-晚三叠世海陆交互相和陆相湖沼沉积.沉积相带的空间分布和古地理格局的变迁明显地受到区域构造的控制.区内盆地发育与演化是在华南板块与扬子板块、扬子板块与华北板块相互作用的背景下进行的.早-中三叠世末期的印支运动是下扬子区构造-古地理格局改变的决定性因素,使下扬子海盆闭合,沉积类型由海相沉积变为陆相沉积.下扬子区三叠纪古地理的演化特征也为板块碰撞提供了沉积证据.     

多年冻土融沉性分类研究

在冻土工程中,冻土的融沉性评价是工程地质勘察的主要内容之一,融沉性分类是冻土地基基础设计施工的重要依据．根据345个冻土原状样品融沉压缩试验数据,提出了细砾、砂土、粉土、黏性土、泥炭化黏性土和泥炭质土等6类土的融沉系数一含水量或融沉系数一超塑含水量线性回归方程式,得到与各融沉性分级相应的界限含水量或界限超塑含水量．最后...