基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

利用改进的粒子群算法反演视震源时间函数

视震源时间函数的提取是研究震源参数的重要途径. 本文提出了利用改进的粒子群(PSO)算法反演视震源时间函数的方法, 以水平线方法得到的结果作为PSO算法的初值, 并对PSO算法的惯性因子和学习因子进行改进, 提高计算效率. 采用改进的PSO算法对模拟数据进行了反演计算, 并与映射Landweber反褶积(PLD)方法和遗传算法(GA)进行了对比分析. 结果表明, 相对于PLD方法, 改进的PSO算法反演结果与真实结果误差更小; 相对于遗传算法, 改进的PSO算法计算效率提高了5倍以上. 最后, 利用改进的算法对2005年10月8日巴基斯坦克什米尔M_W7.6地震的P波视震源时间函数进行了提取, 结果表明此次地震P波视震源时间函数在25 s之内, 震源沿西北向破裂. 该结果与张勇等的结果一致.     

利用PLSA技术进行图像分割

引用概率潜语义分析(PLSA)技术作了航空影像图像分割的尝试,介绍了PLSA图像分割的原理方法和分割的具体算法。通过三幅航空影像分割实验,并与生物地理学图像分割(BBOIS)方法和基于吸引子算法的结果作了对比,结果表明,本文方法具有一定的优越性,是一个很有潜力的图像分割方法。     

一种基于分割的机载LiDAR点云数据滤波

针对当前滤波算法在处理地形不连续区域或存在复杂建筑物区域时容易过分“腐蚀”地形并难以去除一些低矮植被的不足,提出了一种基于分割的机载LiDAR点云滤波算法。首先,对原始点云基于地表连续性进行分割;然后,在移除点数目较小的粗差点集之后采用对分割点集建立缓冲区的方法,区分地面和非地面点集;在较大地物经过迭代分割基本移除之后,使用约束平面的方法移除高度较小的地表附着物以实现滤波。实验结果表明,与经典滤波算法相比,该算法提高了地面点的分类精度,在滤除地物信息的同时能有效地保留地形特征。     

GPS单频机电离层延迟改正新算法

探讨了几种新的电离层延迟改正算法,通过算例检验了新方案的效率和可行性,对不同精度用户选取电离层延迟改正方案给出了建议。     

卫星重力梯度测量数据的粗差探测

分别讨论了阈值法、Dixon检验法和小波算法应用于卫星重力梯度数据粗差探测的效果。为了克服这3种单一粗差探测方法的不足,提出了卫星重力梯度数据的粗差探测组合方案,即联合Dixon检验和小波算法,以及联合Dixon检验、阈值法和小波算法的组合方案。模拟结果表明,两种组合方案均能有效改善粗差探测的效果。     

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

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

     

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

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

     

基于萤火虫和蝙蝠群智能算法的瑞雷波频散曲线反演

反演瑞雷波频散曲线能有效获取地层横波速度和厚度.但由于其高度的非线性、多参数、多极值等特点,传统的全局搜索方法易出现收敛速度慢、早熟收敛及搜索精度低的问题.鉴于此,本文提出并测试了基于萤火虫优化算法（FA）和带惯性权重的蝙蝠优化算法（WBA）的新的瑞雷波频散曲线反演策略.在瑞雷波频散曲线反演中,FA全局搜索能力强,但后期搜索精度低,而WBA局部搜索能力强,搜索精度高,但易出现早熟收敛.故本文将二者结合,提出了一种新的优化策略,称其为WFBA,即在反演前期使用FA,后期使用WBA,很好地解决了FA后期搜索精度低及WBA早熟收敛的问题.本文首先反演了三个典型理论模型的无噪声、含噪声的数据,验证了WFBA对瑞雷波数据反演的有效性与稳定性.然后将WFBA与WBA、FA单独反演以及不含惯性权重的FBA和粒子群优化算法（PSO）反演的结果进行了对比,说明了WFBA相对于WBA、FA、FBA和PSO具有更稳定、收敛速度更快、求解精度更高等优点.最后,反演了来自美国怀俄明地区的实测资料,检验了WFBA对瑞雷波数据反演的实用性.理论模型试算和实测资料分析表明,WFBA很适用于瑞雷波频散曲线的定量解释,具有很高的实用性价值.

     

四川芦山7.0级强震及其余震序列重定位

本研究采用双差定位法对2013年4月20日发生在龙门山断裂带上的四川芦山7.0级强震及主震后48小时内504次余震序列进行重新定位,最终得到328个精定位地震事件.结果表明,余震在水平方向上主要沿龙门山断裂带的山前断裂西南段的大川—双石断裂发生,扩展模式以西南向为主(约23 km)兼有弱东北向(约12 km),并非简单的单侧扩展.在深度方向上,余震主要以铲状结构分布在18～22 km之间.通过拟合精定位后的小震空间分布特征,显示本次地震断层面倾角在15～25 km深度范围内由主震处约44°逐渐向西南向扩展增大至约73°,可能表明断层往西南破裂过程中走滑分量逐渐增强,与2008年汶川地震引起中央断裂倾角由西南向东北变化相类似.