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射线劳埃衍射技术对单晶叶片的晶体缺陷进行无损检测,但是这种检测方法主要依赖人工识别,效率低,结果可重复性差,不适合批量化检测。本文结合工程需要,提出对劳埃衍射斑点进行自动识别的算法,主要包括衍射图样的预处理、轮廓检测、轮廓形态筛选及轮廓符合检测等。该算法能够自动检测出衍射图样上的衍射斑点,并最终给出斑点的位置坐标数据及其误差。根据衍射斑点的位置,通过衍射分析算法,计算出单晶叶片上的晶体取向,并最终对叶片的晶体缺陷给出综合评价。      

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

反演瑞雷波频散曲线能有效获取地层横波速度和厚度.但由于其高度的非线性、多参数、多极值等特点,传统的全局搜索方法易出现收敛速度慢、早熟收敛及搜索精度低的问题.鉴于此,本文提出并测试了基于萤火虫优化算法（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变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

基于稀疏表示的地面磁共振信号提取方法

在多孔隙含水层中地面磁共振（surface nuclear magnetic resonance,SNMR）信号呈现多弛豫衰减特性,常规盲源分离方法和单指数拟合方法引起信号严重失真和信息缺失等问题.本文提出了基于稀疏表示的随机噪声背景下多弛豫SNMR信号的提取方法.根据SNMR信号的衰减特征,设计了精确刻画SNMR信号且与随机噪声不相关的离散衰减余弦冗余字典.其次,针对多弛豫SNMR信号稀疏度未知的问题,通过设置合理的残差比阈值控制迭代次数,改进了广义正交匹配追踪（generalized orthogonal matching pursuit,gOMP）算法,使得该方法应用于SNMR信号的提取时,具有更好的自适应性和普适性.再次,鉴于SNMR测量数据为多次独立重复采集的结果,提出了基于数据流的SNMR信号提取策略,在提高算法鲁棒性的同时,保证了信号提取结果的唯一性.最后,通过仿真和实测数据证明了基于gOMP算法的稀疏表示方法可以显著地提升多弛豫SNMR信号的提取质量,降低随机噪声对含水层反演结果的影响,提高SNMR探测能力.     

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

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

汶川地震科学钻探3号井孔附近微震震源参数确定及其意义

基于微震监测仪器频带及波形记录的时频特征,我们筛选出汶川地震科学钻探3号井孔周围微震台阵2012年记录的218个t_s-t_p1 s的微震.通过盖戈法与和达法相结合确定微震震源的几何参数,发现这些微震分布大体呈NE-SW展布,与龙门山断裂带的走向基本一致.无论是采用《地震台站观测规范》中的量规函数(量规函数GF)还是李学政等(2003)的量规函数(量规函数LXZ),近震震级均与矩震级呈现出较好的线性关系,向震级小的一端延伸时都表现为M_LM_W,但采用李学政等的量规函数时该趋势更加明显.同时采用Brune和Boatwright震源谱衰减模型对观测震源谱的拟合表明拐角频率具有模型依赖性:基于Brune模型拟合得到的拐角频率大于基于Boatwright模型得到的拐角频率.无论基于哪一种模型,矩震级与拐角频率、破裂半径的对数线性关系均较弱,与应力降、视应力的对数则存在较好的线性关系.这些关系不支持应力降、折合能量、视应力为常数的观点,表明微震与大地震的震源物理过程存在差异.视应力与应力降成比例,比例系数小于0.5,表明破裂动力学模式符合SavageWood模式.近震震级与矩震级拟合关系Mw=a+bM_L中b的大小与应力降和地震矩的关系有关,△σ∝M_0~γ,则b=1/(1+γ),因此从b的大小可以粗略地判断应力降与地震矩的关系.本文对应于γ≈1的情况,与基于量规函数LXZ得到的近震震级与矩震级关系中系数b=0.53吻合,这说明仅从辐射能量的角度考虑量规函数LXZ较量规函数GF准确.     

基于动态图像变形的PP与PS波层位直接匹配

多分量地震资料的矢量偏移、多波地震资料的联合解释与反演均需要估算纵横波的速度比,实现纵波与转换横波在时间或深度域的匹配.基于DTW,本文实现了一种适用于PP与PS波直接匹配的动态图像变形算法.该算法分为三个部分:首先,使用二阶对称动态规划算法逐样点递归计算PP与PS波走时或深度的误差累积和;其次,在以误差累积和为目标函数的回溯阶段设定变形窗,并在纵横波速比约束的变形窗内递归回溯搜索匹配路径;最后,根据最大相关系数判定准则在匹配路径中确定最佳匹配路径,获得使PP与PS波匹配的拉伸或压缩时移量.利用所获得的拉伸压缩时移量计算纵横波速度比就可以实现PP与PS波之间的匹配.模型与实际陆上多分量地震资料测试结果表明:该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.     

基于卷积完全匹配层的非规则网格时域有限元探地雷达数值模拟

复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟.     

水平层状非均质TI地层中仪器偏心情况下三维感应测井响应高效数值模拟与响应特征分析

本文针对中国石油集团测井有限公司(CPL)三维阵列感应测井仪器(Three Dimensional Induction Tool,TDIT)的真实结构与参数,利用2.5维数值模式匹配技术研究建立水平层状非均质TI地层中仪器偏心情况下三维感应测井响应的高效数值模拟算法.首先,利用2.5维数值模式匹配技术给出水平层状非均质TI模型中偏心磁流源并矢Green函数的半解析解.在此基础上,结合三维阵列感应测井仪器参数以及测量过程中出现的仪器偏心与旋转等实际情况,通过叠加原理推导出该仪器响应的有效计算公式.最后,通过数值模拟结果分析考察仪器偏心、仪器自旋角等参数对测井响应的影响以及变化特征,为多分量感应资料处理与解释提供理论基础.