深部软岩巷道锚注支护数值模拟研究

采用大型有限元数值模拟软件ANSYS,对深部软岩巷道锚注支护前后围岩变形破坏规律进行了数值模拟,对锚注支护前后围岩的应力、位移及塑性区的变化情况等进行了系统分析,结果表明,锚注支护显著提高了围岩的强度和承载能力,有效地控制了深部软岩巷道的损伤变形。     

坡面薄层流水动力学特性的实验研究

采用变坡实验水槽研究了坡面薄层水流水动力学特性——流态、流速、水深及阻力系数随流量和坡度的变化规律,实验坡度5°～25°、单宽流量为0.625×10-3～12.5×10-3m3/(s·m).研究结果表明,坡面薄层水流流态与水深密切相关,坡面流的流态基本上呈过渡流和紊流;薄层水流平均流速、水深和Darcy-Weisbach阻力系数主要受流量控制,坡度的影响并不显著,其关系均可用流量和坡度的幂函数形式模拟,相关系数(R2)分别为0.98、0.97和0.95,去掉坡度后相关系数分别下降7%、16%和3%.随着雷诺数的增大,Darcy-Weisbach阻力系数下降.当单宽流量小于0.008m3/(s·m)时,坡度对阻力系数的影响较为显著;当流量大于0.008m3/(s·m)时,阻力系数基本受流量控制.随着流量增大,阻力系数呈幂函数形式递减.     

五矩二流太阳风等离子体特性的数值研究

本文数值求解了各向同性二流太阳风的五矩方程组,得出了1 R_s-2AU区域内太阳风密度、速度、电子和质子温度、它们的热流通量密度q以及非麦克斯韦分布尾部粒子过剩量ξ随日心距离的变化关系.文中比较了二流太阳风五矩模型、四矩模型（ξ=0）和低阶矩模型（不包括q和ξ二个矩方程）的等离子体特性,着重讨论了量ξ对质子温度及其热流通量的影响.结果表明,包括言的五矩方程可改善T_e/T_p和q_e/q_p的计算值与观测值的符合程度.     

典型深街谷内树木空间配置对行人呼吸高度处气流的影响

为揭示树木的不同空间配植方案对行人呼吸高度气流的影响,本文将树木视为均匀多孔介质,通过附加源项法从空气动力学角度用CFD模拟了H/W=2的典型深街谷几何内4种树木配植情景,实验表明,不同空间配置下树木对街谷内行人呼吸高度处局地气流的影响强弱在空间分布模式上差异悬殊：① 均匀种植的树木对街谷内行人呼吸高度的气流起到阻碍作用,不均匀种植则有效提升街谷的整体流速。4种空间配植方案下树木对气流的影响程度不同,阻碍作用从大到小的顺序为均匀间距8 m（Spa8m）>均匀间距6 m（Spa6m）>均匀间距20 m（Spa20m）>不均匀配植（Non-uniform）;对应的平均气流增强指标顺序为$\bar{D}_{spa8}$（-19.31%）<$\bar{D}_{spa6}$（-16.14%）<$\bar{D}_{spa20}$（-10.73%）<$\bar{D}_{non-uniform}$（1.25%）。② 对比不均匀和均匀的种植方案,不均匀植树的街谷内部行人呼吸高度的气流流速比其对照案例（均匀植树Spa8m方案）整体增强了106.49%。街谷中部不种树,在街谷两端配置树木并预留足够的自由空间的不均匀植树方案,能够让角涡渗入街谷中部,促使街谷内部的垂直漩涡和两端的水平角涡运动,增强湍流和垂直交换,有效减少了街道两端“风口效应”和街道中部“风影效应”的区域,改善了整个街谷行人呼吸平面的风环境。④ 合理空间配置的树木能够改善街谷内部的行人风环境。街谷内行人呼吸高度处的气流对局地条件很敏感,树木的局部配置（空间簇集、密度）将引起强烈的空间变化。在既有城市建筑布局条件下,如何通过谨慎的景观设计,利用树木等城市绿化措施有效地改善城市的行人风环境,缓解污染扩散、疾病传播等问题,本文的方法可提供一定的参考。     

北淮阳及其邻接区地壳稳定性研究

运用多学科、多技术相结合的研究方法,从现今构造应力场特征、现代地壳形变特征、主干断裂构造以及区域断裂活动性分析、新构造运动、地震活动、数值模拟等方面入手,对北淮阳及其邻接区地壳稳定性进行了全面而系统的研究,指出本区地震活动受走滑活动断裂控制,中、强地震集中于区内东南部和西北部NW向与NE向断裂交汇处;第四纪无大规模差异升降运动;数值模拟等综合分析表明,信阳地区的区域构造稳定性高,现代地壳形变量小,应力平缓,地震活动少且弱。     

54与80坐标系转换数学模型研究

54北京坐标系与80西安坐标系转换计算,属54系与80系不同参考椭球下高斯—克吕格投影数据转换计算,因全国不同区域重力场的变化而无法用一个固定的参数或公式推算,能否另辟途径实现其精确算法。应用多元统计分析基础理论,研究二者互换随机数学模型,以福建省区域为例,成功实现大批数据坐标转换。     

室内常规土工试验的试样制备问题探讨

笔者依据国内近十年来颁布的各类土工试验规范和多年从事室内土工试验的经验,现就试样的制备问题,综合提出制备试样“十四种要素”以达到抛砖引玉的目的。     

一阶弹性波交错网格时间高阶差分格式及稳定性分析

弹性波模拟或逆时偏移时,对空间偏导数采用高阶差分格式可提高计算精度,但这种算法的稳定性条件过于严格,要求差分离散的时间步长必须足够小以确保算法稳定。在常规空间高阶差分格式的基础上,将速度(应力)对时间的高阶导数转化为不同精度的应力(速度)对空间的差分,得到了一种新的基于交错网格的时间高阶、空间高阶差分格式。通过对交错网格时间高阶差分格式稳定性的分析,认为该算法的稳定性条件较常规算法宽松,在弹性波场的求解过程中可以采用更大的时间步长。     

张裂型盆地地热参数的垂向变化与琼东南盆地热流分布特征

海洋热流数据是开展海洋地球动力学研究和油气资源评价的基础数据.为深入认识琼东南盆地的地热特征,本文首先利用耦合沉积作用与岩石圈张裂过程的数值模型分析了张裂型盆地主要地热参数的垂向变化特征;并通过钻孔资料的详细分析,获得了琼东南盆地44口钻孔的热流数据;结合海底地热探针获取的热流数据,对琼东南盆地地热特征及其主要影响因素进行了简要分析.结果表明:沉积作用的热披覆效应对表层热流有较明显的抑制作用,由于沉积物生热效应与披覆效应的共同作用,同一钻孔处海底表层热流与钻孔深度3000~4000 m处热流或与海底间的平均热流差异很小,可以一起用于分析琼东南盆地的热流分布特征;莺歌海组、乐东组热导率随深度变化小于黄流组及其下地层热导率的变化,钻孔沉积层平均热导率约为1.7 W·(m·K)-1,钻孔地层生热率一般低于2.5 μW·m-3,平均生热率为1.34 μW·m-3,平均地温梯度主要介于30~45 ℃/km,热流介于50~99 mW·m-2,陆架区热流主要集中于60~70 mW·m-2,深水区钻孔具有较高的地温梯度和热流值;从北部陆架与上陆坡区往中央坳陷带,热流值从50~70 mW·m-2,增高为65~85 mW·m-2,并且往东有升高趋势,在盆地东部宝岛凹陷、长昌凹陷与西沙海槽北部斜坡带构成一条热流值高于85 mW·m-2的高热流带.进一步分析认为,琼东南盆地现今热流分布特征是深部热异常、强烈减薄岩石圈的裂后冷却作用、晚期岩浆热事件、地壳与沉积层的生热贡献以及沉积作用的热披覆效应等多种主要因素综合作用的结果.     

吕宋海峡M2内潮的三维数值模拟

A three-dimensional isopycnic-coordinate internal tidal model is employed to investigate the generation,propagation, vertical structure and energy conversion of M2 internal tides in the Luzon Strait(LS) with mooring observations. Simulated results, especially the tidal current amplitudes, agree well with observations,demonstrating the reasonability and accuracy of the model. Results indicate that M2 internal tides mainly propagate into three directions horizontally, i.e., eastward towards the western Pacific Ocean, westward towards the Dongsha Island and southwestward towards the South China Sea Basin. In the horizontal direction, tidal current amplitudes decrease as distance increases away from the LS; in the vertical direction, they show an obvious decreasing tendency with depth. Between the double ridges of the LS, a clockwise gyre of M2 baroclinic energy flux appears, which is caused by reflections of M2 internal tides at supercritical topographies, and resonance of M2 internal tides happens along 19.5° and 21.5°N due to the heights and separation distance of the double ridges. The total energy conversion in the LS is about 14.20 GW.