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

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

层状侵入体韵律层理形成过程的定量模拟——以四川攀枝花层状岩体为例

层状侵入体的韵律层理是一种特殊的岩石构造,长期以来其成因倍受关注,成因假说很多。本文在综合这些观点关于岩浆物理,化学性质的研究,尤其是在岩浆房地对流特征研究成果的基础上,提出了新的成因解释,认为固化带前锋处矿物的结晶,熔浆组分的扩散及固化带前锋本身的推进等３个因素的相互耦合是形成韵律层理的主要机制,这个过程是发生在一定过冷却条件下的,在这种定性认识基础上,笔者提出了旨在对韵律层理形成过程进行模拟计     

青藏高原气象学的研究进展和问题

分高原天气学、高原气候学、高原及邻近地区的大气环流、以及高原数值预报和模拟四方面简要回顾了新中国成立以来我国（也兼及国外）青藏高原气象学的主要进展，也提出今后研究中应注意的有关问题。     

大埋深条件下降雨入渗补给的初步分析

运用SWAP模型对太行山山前平原--河北省栾城县田间土壤水分运动进行了数值模拟.从土水势及土壤含水量等方面,初步分析了降雨灌溉入渗补给条件下土壤水分运动的规律.通过分析和讨论,探求了大埋深条件下降雨入渗过程,初步揭示了降雨入渗对地下水位变化的影响.     

小浪底水库淤积形态数值模拟

水库淤积形态是影响库容分布、水库排沙的一项重要因素。小浪底水库近坝段淤积泥沙粒径极细,具有流动性,针对其细颗粒泥沙淤积特点,揭示了水库细颗粒淤积物的流变特性与流型特征;通过引入水、淤积物、床面之间的界面受力分析,构建了细颗粒淤积物失稳流动描述模式,并与水沙输移模型相耦合,建立了考虑细颗粒淤积物流动特性的水库淤积形态模拟方法,在此基础上对小浪底水库淤积形态进行了验证分析。研究结果表明:低密度细颗粒淤积物为宾汉型流体,淤积平衡坡降较小,当其密度大于1.25 g/cm 3后,流动性快速减弱;考虑细颗粒淤积物流动特征的水库淤积形态模拟结果与实测结果吻合较好。研究成果可为水库淤积形态形成机理及其对水沙调控的响应研究提供技术支撑。     

甘孜州色达县典型牵引式滑坡地质灾害形成机理分析--以甘孜州甲学乡甲热滑坡为例

甲热滑坡位于甘孜州色达县甲学乡境内,属于典型的牵引式滑坡,滑坡体平均坡度约35°,主要由碎石土组成,滑坡长约195 m,宽约280 m,平均厚约15 m,滑坡体总方量约60万m3,属中型滑坡。滑坡坡体中上部变形现象较明显,若发生滑动,将直接威胁当地41户160人居民的生命财产安全。本文在进行野外现场基础地质调查后,对滑坡体的变形破坏机制进行了定性分析,并利用三位有限元数值模拟进行了定量分析,对以后研究同类滑坡的稳定性具有重要意义。     

论地震数值预报——关于我国地震预报研究发展战略的思考

随着中国工业化和社会现代化进程的加快 ,中国城市化进程必将进一步加速。如何减轻地震灾害的问题正变得日益严峻。尽管地震预测是一个举世公认的国际性科学难题 ,但在强化各种减轻地震灾害措施的同时 ,仍须大力推进地震预测研究。为此 ,需要打破长期徘徊在以地震前兆异常监测为基础的经验性预测局面 ,把注意力尽快转向研究以动力学为基础的数值预报。以GPS为代表的空间对地观测技术 ,岩石圈巨型高分辨率宽频带流动地震台阵观测技术以及数值模拟技术已经为地震数值预报研究提供了前所未有的技术基础。以地震数值预报为目标的GPS阵列地壳形变连续观测 ,高分辨率地壳上地幔结构探测 ,地壳动力学 ,地震孕育和破裂过程的理论、模拟试验和实际观测 ,数据同化和计算软件的开发应成为今后研究的重点。现在的问题是 ,需要积极借助数值天气预报的经验 ,强化多学科 ,多部门的组织协调 ,尽早在有条件的地区开展地震动力学数值预报的科学试验。地震数值预报研究必将极大地促进中国地震科学基础研究和地震预报的进一步发展。     

深厚覆盖层高压旋喷凝结体的宏观力学参数研究

基于现有工程试验资料,通过将高压旋喷凝结体分成成芯凝结体和不成芯凝结体两组材料组成,提出了深厚覆盖层高压旋喷凝结体的统计力学模型。该模型除能反映成芯凝结体的小尺度力学性质外,还可以反映不成芯凝结体的力学性质,也可考虑试验钻孔中不同芯样获得率的影响。应用随机模拟理论,给出了大体积高喷体内成芯凝结体与不成芯凝结体的空间位置及各小尺度力学参数的抽样方法。并通过结合弹塑性有限元和随机模拟两方法,建立了高喷体宏观力学性质的计算机模拟方法。应用所提三维弹塑性有限元随机模拟方法,分析研究了向家坝水电站深厚覆盖层高压旋喷凝结体的宏观力学参数。     

东昆仑断裂粘滑错动对青藏铁路变形效应的数值模拟

东昆仑断裂是青藏高原北部现今仍在强烈活动的地震断裂之一,该断裂的未来地震活动及其突发性粘滑错动是青藏铁路面临的重大工程地质问题。基于东昆仑断裂的运动学特征,通过分别加入8 m和3 m的水平左旋位移,模拟了东昆仑断裂未来地震活动时震中位于铁路线附近和远离铁路2种情形下的铁路变形效应。模拟结果表明：震中位于铁路线附近时,断裂南侧基岩和第四系均发生8 m的左旋走滑位移,而铁路附近的第四系水平位移明显减小,铁轨和道床没有明显的断错,表现为4~5 m的连续左旋弯曲变形;铁路东、西两侧形成NE向的张裂陷和NW向的地震鼓包,而道床和铁轨的垂直位移幅度较小。震中远离铁路时的变形效应与震中位于铁路线附近时的变形相似,但位移幅度减小,铁轨和道床形成1~2 m的连续左旋弯曲变形。因此,东昆仑断裂未来再次发生7~8级强烈地震时,无论地震震中远离铁路还是在铁路附近,其断裂的突发性粘滑错动都将导致青藏铁路的大变形和破坏。     

喜马拉雅山珠峰绒布冰川流域径流模拟

基于2009年5-10月喜马拉雅山北坡珠峰绒布冰川流域实测水文气象数据、 50 m分辨率DEM和中国第一次冰川编目资料, 在HYCYMODEL水文模型中加入冰川消融子模块, 模拟了绒布冰川流域径流过程.冰川消融子模块以海拔5 180 m基站的实测日气温、 日降水作为模型输入, 把气温、 降水插值到该流域40个高程带中, 分别计算各高程带的冰川消融和裸地蒸发, 并考虑液态降水对冰面的加热作用.野外气象观测表明: 2009年5-10月流域海拔5 180~5 750 m内, 月气温递减率在0.63~0.73 ℃·(100m)^-1之间, 均值为0.70 ℃·(100m)^-1; 同期降水观测显示, 海拔5 180 m以下降水梯度为-7.3 mm·(100m)^-1, 该高度之上降水梯度为22 mm·(100m)^-1. HYCYMODEL水文模型的敏感性检验表明, 该流域径流变化主要受气温影响, 降水变化引起的径流变化较小, 气温和降水变化对流域径流的影响是非线性的.