基于超完备字典学习的缺失地震数据重构方法

地震勘探目标区域环境的复杂多变性导致采集的地震数据存在不完整或者不规则等问题,针对这一问题,本文在压缩感知相关理论的支撑下,提出了基于超完备字典学习的缺失地震数据重构方法.首先利用K-SVD字典学习技术对地震样本数据进行训练,建立超完备字典对地震数据进行稀疏表示,然后引入高斯随机采样矩阵作为测量矩阵对地震数据进行采样;在数据重构阶段采用分段正交匹配追踪算法实现缺失地震数据的重构.通过与传统的地震数据重构方法对比,本文算法的重构效果在峰值信噪比、信噪比等指标上均优于对比算法,证明了超完备字典学习方法能更好的根据地震数据特征进行稀疏表示,从而获得较好的重构效果.     

2017年9月临城震群遗漏地震检测及发震构造初探

以河北测震台网记录的连续波形数据为基础,采用模板匹配滤波方法,检测2017年9月4日—6日河北临城震群遗漏地震。利用波形互相关扫描,检测到19次被测震台网常规分析遗漏的地震,震级范围为M_(L )0.0—1.2,并标定新检测地震事件的P波和S波到时,估计其震中位置及震级。通过地震精定位和最大地震震源机制解,认为此次临城地震序列的主要发震构造应为近EW向活动断层。     

地震波模拟中多轴复频移近似完全匹配层的吸收效果及稳定性

在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他非分裂完全匹配层,其还具有不改变方程的形式、易于实现等优势.但是当入射波角度较大,边界吸收效果变弱,且残留在边界中的能量使近似完全匹配层变得极其不稳定.多轴复频移近似完全匹配层的提出就是为了改善对大角度入射波的吸收并且提高边界的稳定性.通过实验模拟和矩阵特征值灵敏度来研究多轴复频移近似完全匹配层的吸收效果及稳定性.结果表明该方法不仅能够吸收掠入波,而且对常规入射波的吸收也得到提升,同时拥有更好的稳定性.     

高分辨率遥感影像中建筑物轮廓信息矢量化

针对高分辨率遥感影像的特点,提出了基于多类分割与模板匹配的建筑物轮廓矢量化方法:首先对影像进行多尺度SVM分割获取建筑物候选点;然后在Radon变换结合主轴分析获取建筑物主方向的基础上,引入多类分割思想,构建用于分割的能量函数,基于α-扩展算法解算能量函数,将轮廓线分割为3类边缘线段;接着构建形状先验的边缘模板,进行精确定位以获取边缘线的实际位置,去除锯齿状变形,降低提取结果的影响;最后相邻边缘线相正交得到拐角点,依次连接每一个拐角点,得到规则化的建筑物轮廓。相比于同类其他方法,此方法考虑了边缘点的方向信息和相邻边缘点趋于同一类的先验知识,可得到近似全局最优的边缘线段分割结果,避免了规则化过程中选择初始点和处理顺序的麻烦和不利影响,同时充分利用了影像特征,对边缘线段进行精确定位,减弱建筑物提取结果误差的影响。对不同影像的实验结果证明此方法可得到规则化的建筑物外轮廓线,不考虑提取过程中遗漏的建筑物,矢量化结果平均准确度为89%、完整度98%、几何形状相似性87%、整体质量85%。     

利用全波形匹配方法确定水力压裂诱发地震震源机制

准确确定水力压裂诱发地震的震源机制对于描述裂缝破裂类型以及工区地应力状态十分重要.本文采用全波形匹配的方法确定诱发地震震源机制解,在拟合波形的同时对实际数据P波初动极性和纵横波振幅比也进行匹配,并在此基础上发展了一种基于邻域算法分级优化确定震源机制解的新方法.新方法的优势在于它对解具有更强的约束,并且能够解决现有方法中存在的权值系数选择问题.合成数据测试结果表明新方法能够得到可靠的震源机制解.本文将该方法应用于国内某页岩气田水力压裂诱发地震监测数据,反演结果表明诱发地震的震源类型以走滑破裂为主,表明该地区的最大水平主应力大于垂直主应力.     

一种附加钟速约束的GNSS TDCP测速改进方法

GNSS(global navigation satellite system)载波相位历元间差分（time-differenced carrier phase,TDCP）测速方法因其高精度、无收敛等优势而被广泛应用于动态导航定位中,但受限于观测环境引起的多维粗差与数据异常,TDCP的实际测速精度和稳定性通常难以保证。本文提出了一种附加钟速约束的TDCP测速改进方法,利用GNSS接收机钟速低频变化的特性,通过多项式拟合进行接收机钟速参数的短期预测,当钟速估计值偏离预测时,将拟合值作为约束条件回代到观测模型中重新进行速度估计,从而克服粗差的影响,提升TDCP测速的精度与稳定性。通过Ublox-F9P接收机和Huawei Mate40手机在城市复杂环境下的车载实测数据进行了算法验证,结果表明：附加钟速约束后,Ublox-F9P和Huawei Mate40的E、N、U方向测速精度RMS分别为0.11 m/s、0.10 m/s、0.14 m/s和0.20 m/s、0.19 m/s、0.18m/s,水平测速精度分别提升了84%和45%,高程测速精度分别提升了33%和10%。     

非稳态地震记录时变子波估计

地震子波估计是地震资料处理和解释中的一个关键问题,子波估计的可靠性会直接影响反褶积和反演的准确度.现有的子波估计方法分为确定型和统计型两种类型,本文通过结合这两类方法,利用确定型的谱分析法和统计型的偏度最大化方法,分别提取时变子波的振幅和相位信息,得到估计的时变子波.这种方法不需要对子波进行任何时不变或相位等的假设,具有对时变相位的估计能力.进而利用估计时变子波进行非稳态反褶积,提高地震记录的保真度,为精细储层预测和描述提供高质量的剖面.理论模型试算验证了方法的可行性,通过实际地震资料的处理应用,表明该方法能有效地提取出子波时变信息.     

基于反射界面自动拾取的AVA反演方法研究

本文将逐点匹配滤波思想应用于叠前AVA反演,提出了一种反射界面自动拾取方法,从而将叠前非线性反演问题转化为两步线性反演进行求解.首先基于逐点计算的方式获取每个采样点对应的残差值,然后在残差曲线上自动拾取波谷以定位反射界面,最后在频率域采用最小二乘反演对叠前参数进行更准确的估计.相比于传统的同时估计反射界面位置和幅值的非线性反演方法,该方法不需要关于参数的个数、分布等先验假设,并且具有运算速度快、求解稳定的优点.模型试验以及实际地震资料处理验证了该方法的有效性.     

弱形式时域完美匹配层

应用高精度人工边界条件可有效提升近场波动数值模拟计算效率.完美匹配层是吸收层形式高精度人工边界条件,匹配层内场方程和界面条件通常分别采用复坐标延伸技术变换强形式无限域内波动方程和界面条件得到,亦曾将无限域界面条件当作匹配层界面条件.场方程和界面条件构建过程相互独立,可能出现匹配不合理而引发数值失稳、计算精度低下等问题.本文提出采用复坐标延伸技术变换弱形式无限域波动方程以构建完美匹配层的方法.弱形式波动方程耦合了波动方程及界面条件,进而规避了变换后所得场方程与界面条件之间的匹配不合理问题.新方法可直接建立弱形式匹配层,在此基础上亦可给出强形式匹配层.弱形式便于有限元离散,强形式便于有限差分离散.基于弱形式完美匹配层,结合勒让德谱元建立了弹性介质近场波动谱元模拟方案.利用算例验证了新方案的精度及数值稳定性.本文工作可直接推广至多相耦合介质近场波动数值模拟.     

改进的高精度匹配追踪方法研究及应用

本文在常规算法的基础上,提出了一种改进的匹配追踪方法.该方法通过引入新的完备库构建策略,并以雷克子波作为原子,利用原始地震数据最大相关性估计完备库的原子的位置和能量,使得分解精度进一步提高,极大地提高了算法的适应性;为避免原子间隔过小问题,引进了最小原子间距,使得分解效率和分解质量进一步提高.模型试算和实际资料应用表明,本文算法不仅提高了信号稀疏表示的质量,加快了收敛速度,而且算法的适应性和分解精度也得到提高;该方法能够较好地挖掘地震有效信息,提高地震解释精度.