基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

基于匹配追踪和遗传算法的大地电磁噪声压制

针对匹配追踪计算量大、大地电磁数据处理效率低的问题,提出基于匹配追踪和遗传算法的大地电磁噪声压制方法.首先,利用Gabor原子构建过完备原子库,并对过完备原子库集合进行划分.然后,借助遗传算法的自适应性,快速搜寻最优匹配原子及所在位置.最后,运用最优匹配原子对待处理信号进行稀疏分解,重构有用信号.通过对计算机模拟的典型强干扰和矿集区实测大地电磁数据进行分析处理,实验结果表明,相对于匹配追踪和正交匹配追踪,文中所提方法能从过完备原子库中快速、自适应地选取最优匹配原子与不同噪声干扰类型高精度的匹配,极大地提升了计算效率;大地电磁时间域序列中的大尺度强干扰被有效剔除,视电阻率曲线更为光滑、连续,低频段的数据质量得到明显改善.     

几种不同矩阵算法的遥感图像几何精纠正效果比较

分别采用Gauss、Householder及Givens三种变换形式进行遥感图像几何精纠正试验，并将试验结果与商业遥感软件PCI的纠正结果进行比较。结果表明，Givens算法在一至五次多项式纠正中，其总体精度优于其它算法。     

Coupling snow accumulation and melt rate modules of monthly water balance models with the Jazim monthly water balance model

Several models for simulation of water balance processes in semi-arid mountainous basins were developed by coupling different modules of existing water balance models (WBM). Snow accumulation and snowmelt rate relationships extracted from the McCabe-Markstrom, Guo, Rao-Al Wagdany and WASMOD-M WBMs, originally developed for basins with humid climate, were coupled with the Jazim WBM, primarily developed for arid basins. Karaj Basin, central Iran, with snowy autumn–winter and dry summer periods, was selected to assess model performance. The model parameters were optimized using a genetic algorithm (GA). All coupled models performed better than the non-modified (original) WBMs in the study basin. The coupled Jazim–McCabe-Markstrom model provided the best performance in simulating low and high monthly flows. It estimated the snowmelt runoff values more accurately than other proposed coupled models because the linear relationships used in the snow module of the McCabe-Markstrom model are more compatible with snow variations in the Karaj Basin.     

场地效应研究的传递台站线性反演法——以芦山地震为例

采用线性反演法研究了地震动的场地效应,并针对线性反演法的局限性,提出了"传递台站"的方法来扩展其应用.通过传递台站,地震事件不再局限于必须被参考台记录到,台站也不再局限于须与参考台站同时获取到地震记录,改进后可大大增加数据样本,使得结果更趋稳定.同时,对于原来难以参与分析的台站,也可得到其场地反应.研究中,还采用了基于遗传算法的非参考台方法,进一步分析了场地反应和验证所提方法的合理性.将所提方法应用于2013年芦山7.0级地震周边区域,分析了这一地区的强震台站的场地反应和地震波衰减因子Q_s值.研究结果显示,当采用一个和二个传递台站进行分析时,所能利用的地震事件和所能求解场地反应的台站数目逐步增加,随着数据样本的增加,分析结果趋于稳定.研究区域的Q_s值拟合结果为Q_s(f)=334.16f~(0.73).遗传算法结果表明,线性反演法中所选取的参考台站有其自身的场地反应.当利用遗传算法获得的参考场的场地反应去校正线性反演法的结果时,对于多数场地,两者的结果符合得很好,说明了遗传算法的结果可以合理地认为是场地较为真实的反应.     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。     

西藏尼洋河水生生物群落时空动态及与环境因子的关系:1.浮游植物

于2008-2009年按照季节调查了西藏尼洋河浮游植物群落的组成、丰度和多样性,并运用多元统计方法定量分析了浮游植物的空间和季节变化特征及其与主要环境因子之间的关系.结果显示,尼洋河浮游植物共计7门29科48属,其中硅藻为优势浮游藻类.浮游植物Shannon-Wiener多样性指数(Pielou均匀度指数)在尼洋河中游(尼洋河下游)最高,其他河段呈下降趋势,符合中间高度膨胀假说.尼洋河沿程浮游植物的总丰度、物种丰富度、Shannon-Wiener多样性指数以及Pielou均匀度指数不存在显著差异,夏季的浮游植物物种丰富度与其他季节的存在显著差异,夏季总丰度与秋、冬季的存在显著差异,冬季的浮游植物Shannon-Wiener多样性指数与春季的存在显著差异,冬季的均匀度指数与其他3个季节的存在显著差异.尼洋河浮游植物季节演替依赖于外源性水源补充,沿程演替则与河道底质有着很大关系.典范对应分析(CCA)表明,尼洋河硅藻门舟形藻科的藻类与理化因子铵态氮、表层pH、表层水温相关,部分蓝藻以及绿藻与水质理化因子也存在着关联.分类回归树(CART)模型预测了尼洋河着生藻类时空分布与主要环境因子之间的定量关系,尼洋河浮游植物群落总丰度和均匀度指数受pH值影响较大,pH值低于8.0的水域浮游植物均匀度指数比pH值高于8.0的水域大,尼洋河浮游植物Shannon-Wiener多样性指数受到河道底质影响较大,底质为黏土的水域浮游植物Shannon-Wiener多样性指数较底质为砂石的大.这些关键环境因子对尼洋河水域浮游植物的时空变化有着重要的指示作用,建议加强对浮游植物及这些环境因子的关注,保障尼洋河水域生态环境的可持续发展.     

GA辅助NLS的GNSS-IR土壤湿度反演方法北大核心CSCD

采用具有阻尼因子的函数模型,使用遗传算法(genetic algorithm,GA)辅助非线性最小二乘(nonlinear least squares,NLS)方法对相位参数进行求解。结果表明:1)相较于标准余弦函数模型,该方法的反演相位与土壤湿度的相关系数有较为明显的提升,反演结果也更加稳定,在5°~15°、5°~20°、5°~25°三个高度角范围内的相关系数均大于0.68,不同高度角之间的相关系数差值小于0.07;2)反演精度有不同程度提高,R 2提高5.72%~76.06%,RMSE减小6.12%~24.24%,MAE减小2.7%~28.3%,将该方案所求相位用于多星线性回归模型后平均RMSE减小10%。