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

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

富氟花岗岩中萤石岩浆成因的新证据

通过富氟花岗岩-HF-H2O体系的熔化-结晶实验,在620～700℃、1×108Pa条件下获得了石英+碱性长石+黄玉+云母+萤石的矿物组合.在该实验体系中首次发现了八面体萤石.实验结果表明,在岩浆条件下可以形成萤石,从而为富氟花岗岩中萤石的岩浆成因提供了新证据.     

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

辽东地区五龙金矿是华北克拉通东部一个典型的石英脉型金矿床,其成矿过程尚存争议.多种证据表明五龙金矿成矿流体属于岩浆热液,成矿时代又与距离较近的三股流岩体的侵位年龄接近,因而有一种观点认为二者密切相关.目前尚未有充足的地质证据表明三股流岩体释放了五龙金矿成矿所需要的含金岩浆热液,那么三股流岩体是否为其提供必要的能量呢？利用现有的地球化学和地球物理数据,建立三股流岩体仅作为热源驱动大气降水和释放岩浆热液的对流数值模型,模拟三股流岩体周边的流场和温度场.数值模拟结果表明,三股流岩体自侵位后可在其周边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级强烈地震时,无论地震震中远离铁路还是在铁路附近,其断裂的突发性粘滑错动都将导致青藏铁路的大变形和破坏。     

丰城三溪峡——东乡石马岭推覆构造是扬子板块与华南板块的直接分界线

扬子板块与华南板以丰城三溪峡－东乡石马岭推覆构造分野。该推覆构造发育在萍乡－广丰深断裂与遂川－德兴深断裂的复合部位，具有可行－逆冲推覆性质。原地系统为扬子地层分区中元古界牛头岭组、高桥组和将军岭组，外来系统是华南地层分区新元古界源里组，分别相当于蓟县系修水组和青白口系上施组同期沉积。根据推覆变形形式、构造岩特征和变形程度，可将推覆体划分为包括糜棱岩亚带、糜棱岩化亚带在内的变形主带和劈理化带。韧性剪切变形与推覆变形同期产生，但影响范围更广泛。     

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

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

一种新的、快速的碳、氮、硫同位素测定手段──EA－IRMS连线分析技术

一种新的、快速的碳、氮、硫同位素测定手段──ＥＡ－ＩＲＭＳ连线分析技术储雪蕾（中国科学院地质研究所，北京１０００２９）关键词元素分析仪－同位素比值质谱仪，连线分析碳，氮，硫同位素碳、氮、硫三种元素的稳定同位素在地学、生态学、环境科学等领域受到了重视，...     

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

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

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

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