非均匀介质中地震波传播的数值模拟

地震波场数值模拟方法多种多样,各种方法都有各自的特点.这里推出一种全新地震波场模拟方法--基于Forsyte广义正交多项式的迭积微分算子法,该方法将计算数学中的Forsyte多项式,应用到地震波传播的数值模拟中,它同时具有广义正交多项式迭积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型数值模拟结果说明了该方法的可行性.将该方法的计算结果与傅氏变换伪谱法、错格高阶有限差分法相比较,结果说明,该方法在波场模拟方面具有较好的发展潜力,并具有自身独到的优越性.     

辽东三股流岩体是否为五龙金矿成矿提供必要的能量？

辽东地区五龙金矿是华北克拉通东部一个典型的石英脉型金矿床,其成矿过程尚存争议.多种证据表明五龙金矿成矿流体属于岩浆热液,成矿时代又与距离较近的三股流岩体的侵位年龄接近,因而有一种观点认为二者密切相关.目前尚未有充足的地质证据表明三股流岩体释放了五龙金矿成矿所需要的含金岩浆热液,那么三股流岩体是否为其提供必要的能量呢？利用现有的地球化学和地球物理数据,建立三股流岩体仅作为热源驱动大气降水和释放岩浆热液的对流数值模型,模拟三股流岩体周边的流场和温度场.数值模拟结果表明,三股流岩体自侵位后可在其周边1 km范围内形成300℃以上的热场,并维持数十万年,因而有利于形成岩浆热液型矿床.然而,五龙金矿所在的位置受三股流岩体侵位后的热场影响不显著.如果三股流岩体释放了充足的岩浆热液,岩浆热液可在岩体边缘和接触带形成热液蚀变和金矿化,这与现有的地质证据不符.故五龙金矿可能与三股流岩体没有直接成因关系,而金成矿热液来源仍需深入研究.      

TRMM遥感降水数据在伊洛瓦底江流域的精度检验和校正方法研究

在伊洛瓦底江流域利用13个雨量站数据对TRMM日降水数据进行了精度检验,提出了一种基于水量平衡的TRMM数据校正方法。该方法首先以SWAT模型为基础,以不同子流域内的水量平衡为目标,计算各个子流域的TRMM数据校正系数;其次以子流域内的坡度值为自变量,校正系数为因变量,对二者进行回归分析,在没有水文控制站的子流域根据坡度求得对应的校正系数,对各个子流域的降水数据分别进行校正。结果显示,校正前TRMM年降水量和月降水量与实测值有较高相关性,但偏差较大;校正后年降水量与实测值基本一致;月降水量与实测曲线吻合较好;日降水量的预报等级从中等提高到了良好,校正前2个站预报等级差,7个良好,4个中等,校正后6个良好,2个中等,5个优秀。结果表明此方法在伊江流域能在不同时间尺度上显著提高TRMM数据的精度,为流域水文水资源分析提供数据支持。     

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

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

裂隙含水层水文地质参数反演——以黑龙江七台河市应急水源地为例

为了获取裂隙含水层水文地质参数,以黑龙江省七台河市应急水源地为例,在抽水试验和示踪试验基础上,利用数值模拟软件GMS(Groundwater Modeling System)建立地下水流数值模型和溶质运移模型。通过研究溶质运移单域模型和两区模型,发现在裂隙含水层中两区模型能更好的描述穿透曲线的拖尾现象。实例研究表明利用数值模拟方法反演的含水层参数与解析法结果很接近,数值模拟的水均衡误差为0.024%,参数可靠性较高。研究发现该地区裂隙连通性较好,可作为应急水源地。     

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

基于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水文模型的敏感性检验表明, 该流域径流变化主要受气温影响, 降水变化引起的径流变化较小, 气温和降水变化对流域径流的影响是非线性的.     

多组分离子交换吸附反应-运移的土柱试验及模拟

利用室内一维饱和垂直土柱对保守性溶质运移和反应性溶质运移分别进行了试验,并用水文地球化学模拟软件PHREEQC模拟了保守性溶质Cl-的运移以及K+、Ca2+、Na+、Mg2+4种主要阳离子的反应运移过程,分析了土壤中发生的离子交换吸附和盐分运移过程,进一步模拟分析了不同入流溶液种类、次序以及有无置换反应情况下出流溶液的差别。结果表明:这些条件的不同均会对出流溶液的种类和浓度产生明显影响,表明水文地球化学反应对地下水中污染物的迁移会产生很大影响,在相关领域的研究中不宜忽略。     

自适应非线性有限元的多重网格法求解及应用

描述了非线性自适应有限元的多重网格法求解过程。 通过算例说明了方法的有效性以及自适应多重网格有限元用于开挖施工过程模拟分析的优越性。     

离心模型试验与数值模拟相结合研究采空边坡渐进破坏特性

以贵州六盘水地区煤层开采为工程背景,采用离心模型试验与数值模拟相结合的方法研究了采空条件下边坡的渐进破坏过程及破坏模式。基于离心模型试验得到了开采过程中边坡的位移场,并采用FLAC~(3D)对离心模型试验过程进行了数值模拟,数值与试验结果得出的边坡位移时程曲线吻合程度较高,塑性区分布与试验点对分析得到的滑裂面破坏过程与路径基本一致。通过对比、融合离心模型试验观测与数值模拟结果,表明贯通坡脚与采空区的滑裂面由下向上形成,且总体发生剪切错动破坏,贯通坡肩与采空区的滑裂面由上至下发生复杂的破坏过程,初始阶段滑裂面发生张拉破坏,中间阶段同时发生张拉与剪切破坏,末尾阶段发生剧烈的剪切滑移破坏。     

发展新一代矿产勘探技术——澳大利亚玻璃地球计划的启示

在全球矿产勘查日益萎缩的今天，人们极为重视发展新一代大型矿床的勘探技术。澳大利亚最近开展了一项雄心勃勃的计划，即玻璃地球计划，目的是研制三维可视化和地质模拟等技术，使大陆表层一千米“像玻璃一样透明”。该计划以信息技术为核心，注重技术之间的综合。从中得到的启示是，应以提高精度为出发点，加强矿产勘查技术综合示范工作；以信息技术为核心，大力发展三维可视化技术和模拟技术；注重多部门的协作，发挥社会主义集中力量办大事的优越性。