三维复杂地壳结构非线性走时反演

中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景.     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

淮南顾北矿F104断层两侧岩溶水化学形成机制及导隔水性评价

利用断层两盘岩溶水化学特征判断其导隔水性对于我国华北型煤田水害防治具有十分重要的实践意义。以顾北矿F₁₀₄断层两侧太原组岩溶地下水为研究对象,在分析断层两盘水文地质条件基础上,采用Piper三线图、离子组合比和主成分法,分析了主要组分来源及水-岩作用差异性,并采用PHREEQC软件对岩溶地下水进行反向水化学模拟。结果表明:南北区岩溶水均存在方解石和白云石的溶解和沉淀现象,南区阳离子交换吸附和脱硫酸作用程度强于北区,而北区黄铁矿氧化和岩盐溶解作用较南区明显,南北两区水化学环境及水-岩作用存在显著差异,进而推断F₁₀₄断层具有较好的阻水性,且影响了其两侧的氧化-还原环境及温度差异,控制着地下水径流方向和水-岩作用程度。      

浅议带有内边界的多边形裁剪算法的追踪过程

任意多边形的裁剪算法在地图制图过程中和GIS空间数据处理中具有重要意义。本文依照Weiler-Atherton算法,在不带内边界的无拓扑关系任意多边形裁剪算法基础上,论述了带内边界的无拓扑关系任意多边形裁剪算法的步骤和裁剪子区域多边形的追踪过程,为研究更为复杂的任意多边形裁剪算法提供了思路,使Weiler-Atherton算法应用更为广泛。     

基于倾斜摄影模型地质界线半自动追踪研究

传统二维地质填图生成的地形地质图是地形曲面上对三维地质模型的降维表达,是三维地质模型在地表窗口的三维展现。本文以高分辨率倾斜模型数据(.OSGB格式文件)为源数据,研究三维地质填图的地层界线采集与数据组织方法。研究了以人工干预的方式半自动追踪地层界线的方法,既减少了操作员的工作量,又提高了工作效率,同时还避免了全自动识别精度低、需要攻克的技术难度大等问题。     

基于压缩感知理论的地震数据降噪方法

针对地震数据在采集处理过程中存在的随机噪声,本文从压缩感知的角度,给出了一种地震数据降噪方法.其基本思路是:首先对含有随机噪声地震数据通过离散余弦变换进行稀疏表示,然后选取随机高斯矩阵为测量矩阵,并计算出传感矩阵,在地震数据重构阶段,采用正交匹配追踪算法对地震数据进行重构;通过实验方法对比,本文方法的降噪效果在峰值信噪比、信噪比、均方误差指标上均优于对比方法,证明了本文方法对地震数据中的随机非平稳噪声有较好的压制效果,提高了地震数据的信噪比.     

水力压裂微地震事件定位方法综述

随着国内外页岩气勘探开发的不断深入,水力压裂微地震监测技术重要性越来越大,而该技术的核心是微地震事件定位.近些年前人对微地震事件的定位方法做了比较深入的研究,形成了众多定位方法.为此,依据目标函数构造及其求解方式对微地震事件定位方法进行了分类,主要介绍了两大类具有代表性的微地震事件定位方法:(1)基于走时的射线追踪定位方法;(2)基于波形的偏移定位方法.然后对其算法思路流程、优缺点等方面进行了归纳总结,分析了各类算法的研究现状及两种典型算法的应用实例,最后对微地震事件定位方法的发展趋势进行了展望.     

随机波浪下泰勒离散系数的时域解

利用Wolk提出的粒子追踪方程,通过等分频率法划分不规则波谱,利用MATLAB做粒子运动模拟计算,得到无因次化泰勒离散系数K/D随时间t变化的曲线;通过与Huang等得到的P-M谱的泰勒离散系数K/D计算结果比较证明了本计算方法的可靠性。采用该方法研究了不规则波条件下,波序列(同一谱型不同波面序列)和谱型(谱峰周期、有效波高、谱峰升高因子)对波浪离散系数的影响;计算结果表明:同一谱型不同波序列对泰勒纵向离散系数稳定值和稳定时间无影响;不规则波谱峰周期越大,纵向离散系数K/D越小,稳定时间越短;有效波高越大,纵向离散系数K/D越大,稳定时间越长;谱峰升高因子越大,泰勒离散系数K/D越大,稳定时间越长;与规则波相比,不规则波的泰勒离散系数K/D的值略小10%~30%。     

地球自由振荡的保结构模拟

大地震等诸多激励均能激发全球自由振荡现象,通常表现为驻波形式的全球整体振荡.现有的地震波数值模拟方法多为非保结构方法,无法压制长时程计算中的积累误差.本文采用优化的三阶辛格式谱元法,对地球自由振荡及全球尺度的地震波传播进行了长时程模拟.通过与传统的基于Newmark算法的谱元方法结果对比分析,明确验证了本文所得优化的三阶辛格式谱元法在模拟地球自由振荡等大规模长时程问题上的优越性和准确性.上述进展在方法论层面为今后探测、刻画全球尺度地球非均匀结构的驻波数值方法奠定了部分基础,并为相关研究领域提供了新的选择.     

2D各向异性TTI介质中地震波场的射线法和波动方程数值模拟对比

尽管射线理论是基于波动方程高频假设下的一种近似理论,两者本应相互佐证共同发展完善,然而实际上确是平行独立的.鉴于此,本文提出了一种可以相互验证波场数值模拟或射线追踪数值模拟方法正确性和有效性的佐证方法.为了实现上述目的,我们首先将用于各向同性介质中多次波预测的分区多步不规则最短路径算法推广至各向异性TTI介质中,使其具备追踪各向异性TTI介质中多次波的能力.然后,通过伪谱法波场数值模拟和分区多步不规则最短路径算法预测得到的波场快照图、单炮地震记录剖面、以及合成理论地震图的分析对比,不但可以相互验证两种算法的正确性和有效性,同时对地震波在各向异性TTI介质中的传播规律和特点有了深刻的理解和认识.结果表明:本文提出的各向异性介质中多次波射线追踪算法是正确有效的.