基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.     

STA/LTA—AIC算法对地震P波震相拾取稳定性影响

选取区域地震台网记录的地震波形数据,使用STA/LTA算法与STA/LTA—AIC算法,进行地震P波震相初至到时自动拾取,对地方震及震中距较大的震相进行P波震相拾取效果分析,发现:STA/LTA算法对于地方震P波震相识别精度较高,与STA/LTA—AIC算法拾取的P波震相初至到时相差不大;震中距变大后,STA/LTA算法对P波拾取位置相对于最佳位置向后延迟,STA/LTA—AIC算法有效矫正了STA/LTA算法拾取位置的延迟问题,与人工拾取位置差别可忽略不计。     

单晶高温合金X射线衍射斑点自动识别算法研究

单晶高温合金展现出优越的抗疲劳性能和高温蠕变性能,广泛应用于航空发动机和燃气轮机的热端部件。但是,其制备过程中会产生晶体取向偏离、杂晶等缺陷。目前国际上已经普遍使用X射线劳埃衍射技术对单晶叶片的晶体缺陷进行无损检测,但是这种检测方法主要依赖人工识别,效率低,结果可重复性差,不适合批量化检测。本文结合工程需要,提出对劳埃衍射斑点进行自动识别的算法,主要包括衍射图样的预处理、轮廓检测、轮廓形态筛选及轮廓符合检测等。该算法能够自动检测出衍射图样上的衍射斑点,并最终给出斑点的位置坐标数据及其误差。根据衍射斑点的位置,通过衍射分析算法,计算出单晶叶片上的晶体取向,并最终对叶片的晶体缺陷给出综合评价。      

Crustal velocity structure of central Gansu Province from regional seismic waveform inversion using firework algorithm

The firework algorithm (FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude (M _W) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region, inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15, 15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s, respectively.     

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

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

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

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用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).遗传算法结果表明,线性反演法中所选取的参考台站有其自身的场地反应.当利用遗传算法获得的参考场的场地反应去校正线性反演法的结果时,对于多数场地,两者的结果符合得很好,说明了遗传算法的结果可以合理地认为是场地较为真实的反应.