基于MPI并行算法的电阻率法多种装置数据的三维联合反演

由于地表电阻率法受到浅地表局部异常体的干扰,反演精度受到影响,井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰,但是由于钻井位置的局限性,数据量得不到保障,也会导致反演精度降低.为此,本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先,利用有限差分法实现正演模拟,采用非线性共轭梯度法（NLCG）恢复电阻率结构;其次,调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算,当开辟12个进程时,反演程序获得最大加速比4.51;最后,通过合成数据和实测数据算例证明该反演算法的有效性,也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.     

垂直气柱中大气微量成分总含量的反演

提出了一种利用地面中等光谱分辨率的太阳红外光谱反演整层大气中某几种微量成分总含量的方法,用逐线积分（ Ｌ Ｂ Ｌ） 法计算整层大气分子吸收,采用差分吸收法测量太阳红外光谱得到实际大气中微量成分的吸收,用最大似然法反演出微量成分的总量。详细分析了这种方法的可行性和反演精度,数值模拟分析和实际测量表明该方法有较高的精度,已把该方法用于监测实际大气中甲烷和水汽总含量的变化。     

HITRAN数据库对卫星1.6 μm通道CO2辐射传输特性计算的影响

卫星遥感器中的CO21.6μm弱吸收带通道测量信号可以反映CO2的近地层浓度分布,是温室气体卫星反演的重要通道之一。HITRAN数据库是建立卫星遥感CO2浓度算法依赖的重要分子光谱参数数据库,目前已经更新到了2012版,不同版本中大气分子谱线参数存在差异。本文利用逐线积分模式LBLRTM,研究了最近3个版本HITRAN数据库(04、08、12版)在CO2的弱吸收带通道上大气光学厚度、透过率的差异,发现04版计算的气体光学厚度普遍偏高,可对CO2造成约38 ppm的低估;08版本得到的气体光学厚度与12版本接近,反演相差2 ppm以内。在此基础上,分析了不同HITRAN数据库对整层CO2变化和近地层CO2变化的敏感性,结果表明:04版对整层和近地面大气的变化敏感性最强,并且放大了近地层信号;08版与12版对整层或近地层CO2的敏感度接近,两者经过卷积后得到的信号无差异。     

黑河地区沙漠气溶胶浓度和谱分布特征

本根据黑河地区的观测结果．着重分析了沙漠气溶胶浓度及其谱分布特征。沙漠气溶胶浓度较低，其中细粒子占绝对优势。扬尘时气溶胶浓度显增加，其中粗粒子增加幅度最大．比平时高20倍以上．直径(D）在2.0-10．0μm间的粗粒子浓度比重从晴天的1.9％上升到9．6％。中等降雨对D>4．0μm的沙漠气溶胶粒子的清除效率非常高(>90％)，但对总浓度的清除效果不很明显(42．5％)，尤其是对其中的D<0．4μm的小粒子(32．8％)，因此0．5相似文献   

三维复杂地壳结构非线性走时反演

中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

基于经验模式给定热层温度对热层大气密度反演影响的评估

利用加速度计数据反演热层大气密度算法一般需由经验模式给定热层大气温度,进而计算大气阻尼系数C_D.本文基于CHAMP卫星加速度计数据反演得到大气密度,以2008年为例,利用反演得到的热层大气密度循环迭代修正大气阻尼系数C_D,通过对比修正前后密度偏差,评估经验模式给定热层温度对热层大气密度反演造成的影响.结果表明,经验模式热层温度计算偏差对大气密度反演造成的影响小于5%,而且考虑大气成分的改变则进一步降低了这种影响.     

基于热点参考系的板块绝对运动模型

本文在最新一代板块相对运动模型(MORVEL模型)的基础上,通过最小二乘反演观测热点方向数据建立板块绝对运动模型.经过对前人用于约束板块绝对运动模型的所有热点数据与MORVEL模型的一致性作系统的统计学检验,发现存在离群数据.相应地,提出两种系统剔除离群数据的新方法:逐一剔除法和全局搜索法.结合对总残差的卡方检验和对残差频率分布的正态检验,逐步筛选最优模型,最终得到一个基于热点参考系的板块绝对运动新模型(T87模型).该模型能够合理拟合全球分布的87个热点方向数据,而模型预测的板块绝对运动速率比观测热点火山的迁移速率系统偏小,最大偏差达到4 cm·a~(-1).这样的偏差可由观测热点速率的系统误差或地幔返回流造成的热点间运动解释.不论是哪种解释成立,本文结果表明热点方向数据能够独立、有效地定义全球热点参考系.这样定义的热点参考系可方便地应用于板块绝对运动、地幔对流及真极移的研究.     

食物浓度对模糊秀体溞生长和繁殖的影响

模糊秀体潘是热带水体中枝角类的主要优势种,它在竞争及对热带水体的适应方面的优势尚不清楚。本文以模糊秀体潘为研究对象,观察不同食物浓度下其生长和繁殖的特点,探讨模糊秀体潘成热带水体优势种的原因。以蛋白核小球藻为食物,浓度梯度设为0.05mg／L、0.1mg／L、0.3mg／L、0.6mg／L和1.0mg／L小球藻浓度为0.05mg／L的培养液中,模糊秀体潘存活率很低且未观察到繁殖现象其它4个浓度小球藻喂养的模糊秀体潘体长生长率和最大体长均随食物浓度的升高而增加,最大体长分别达到1.087mm、1.125mm、1.225mm、1.313mm;食物浓度对模糊秀体溢的胚胎发育时间,平均寿命,个体累计产幼数无显著影响：随着食物浓度的升高,模糊秀体滏初次怀卵体长显著增加,繁殖前期的历时相应缩短,平均每胎产幼数增加：内禀增长率随食物浓度升高而增大,达到0.6mg／L后比较稳定。浓度为0.6mg／L时,模糊秀体潘的净生殖率也达到最大值。1.0mg／L组的世代时间最短,内禀增长率接近最大值。0.6mg／L和1.0mg／L为模糊秀体潘生长比较适宜的食物浓度。热带自然水体中,模糊秀体潘相对其它浮游动物个体小和怀卵率低的现象并非食物不足所致;食物充足条件下的最大成体的个体大小与野外样品接近,也说明自然水体中捕食并没有显著影响成体的个体大小。模糊秀体潘体型较小、身体透明、行动敏捷,容易逃避鱼类的捕食,可能是模糊秀体潘能在食物丰富的热带水体中,特别是透明度较低的富营养化湖泊中能够成为枝角类中优势种的主要原因。     

羌塘盆地中部地区地壳S波速度结构及构造意义

羌塘盆地是我国最大的海相盆地.本文根据在羌塘盆地内布设的27个宽频带地震观测台站记录的远震波形数据,利用非线性复谱比反演算法得到各台站下方100 km深度范围内S波速度结构.结果表明.羌塘地区Moho深度较为平缓,平均深度为61 km;北羌塘地壳内低速层广泛分布;北羌塘具有两个较大的沉积盆地,龙尾错和白滩湖坳陷,沉积厚度分别有10 km和15 km.尽管北羌塘下地壳受到强烈的新生代火山岩作用改造,但是这种深部岩浆热作用会加速烃源岩中有机质的热演化历程.北羌塘两个盆地具有很好的油气前景.与北羌塘低速层分布相比,南羌塘下低速层更深,可能与班公怒江洋于中生代的俯冲消减及拉萨地体北向俯冲有关.