学位论文 > 优秀研究生学位论文题录展示

隧道地震超前预报中波场分离与反演方法的数值模拟研究

作　者: 王朝令
导　师: 刘争平
学　校: 西南交通大学
专　业: 地球探测与信息技术
关键词: 隧道地震超前预报数值模拟有限元 Radon变换 τ-p变换 F-K滤波波场分离微分算子运动积法基因算法地震反演有限差分交错网格
分类号: U452.11
类　型: 博士论文
年　份: 2012年
下　载: 285次
引　用: 1次
阅　读: 论文下载

内容摘要

隧道属于地下建筑工程,施工中常遇到各种不良地质地段,断层、破碎带等隧道地质灾害时有发生,随着我国隧道工程数量、规模不断扩大,深、大、长隧道增加,施工过程中遇到的工程地质、水文地质条件越来越复杂,隧道地质灾害超前探测的理论和实践研究意义重大,其中以隧道地震超前预报方法应用最为广泛,隧道地震超前预报方法以围岩介质的弹性物理性质差异为基础,是地球物理探测方法在隧道、地下井巷、采场等工作面上的运用,通过观测地下物理场的分布及其变化规律来研究工作面前方是否存在地质异常以及地质体赋存状况及规模。本论文通过弹性力学的本构方程、Cauchy方程和Navier方程三大基本方程,构建了波动方程,通过分析有限元原理,研究了弹性波场的有限元数值模拟方法,分析了弹性动态问题到波场模拟的过程,用ANSYS软件实现了地震波场的2D数值模拟,并将其应用在隧道地震超前预报的数值模拟中开展研究。通过对波场的数值模拟,加深了对其在隧道传播的认识,基于层状地质模型和隧道地震反射法的基本原理,建立多种倾角的模型,采用微分算子、运动积法、τ-p变换、抛物型Radon变换、F-K滤波等波场分离方法对取得数据进行分离；以交错网格差分方法为基础,开展基于基因算法的反演研究。主要认识和创新点如下：(1)实现有限元软件ANSYS模拟隧道地震全波场,并分析了单元类型和网格划分对数值模拟影响、以及粘弹性边界的加载等模拟的关键技术。对于频散效应而言,根据现场岩体速度,当网格划分为1m时,可以很好的消除频散效应；针对边界条件,采用本文所提出的粘弹性边界加载在模型周边,可以有效的消除边界反射效应。(2)应用ANSYS有限元软件数值模拟了隧道地震波场,分析了TSP系统和隧道多道超前预报(TMS)的波场传播特征。数值模拟结果显示,隧道地震波场是包含P波,S波和面波的多波多分量,以及多方向干扰反射波构成的复杂波场。由于传播速度和方向的差异,采用直线排列观测时,来自预报前方的反射波与其他方向干扰反射波的时距曲线特征是不同的。多波多分量和多方向反射波的地震波场时距关系曲线虽然可出现双曲,直线等多种特征,但由于隧道内观测空间有限,排列都比较短,导致观测到的双曲半径很大,曲率非常小,所以隧道内波场的时距曲线都可近似于直线关系,直线关系的斜率由波场的传播速度和方向确定。(3)研究和分析了多种地震波场分离方法的基本原理和算法。数值模型应用研究表明：微分算子方法可以分离出纵横波波场,但由于其计算的特点,只适合对于模拟数据进行研究时采用。运动积方法可以很好分离开纵横波波场,但其要求数据比较纯净,因此对于现场数据不是很适用。(4)线性Radon(τ-p)变换可以很好的分离时距曲线为直线或近似为直线关系的波场,对于隧道前方地层陡倾角变化在90°—45°范围,都可以有效的分离多方向和多波反射事件,提取主要来自掌子面前方的纵波和横波反射剖面。类似的,F-K滤波对于时距曲线为直线或近似为直线关系隧道地震波场分离来说,也是一种比较理想的方法,但由于F-K滤波采用的是扇形滤波区域,其分辨率与数据的频率有关。(5)由于双曲特征干扰波或横波的高视速度影响,以及横波的强能量影响,纵波分离记录中或多或少残留着这类干扰阴影。而横波由于速度很低,在Radon域中与高视速度干扰波分离较开,故横波分离的效果比纵波好。因此,在后续的处理中,应以分离的横波剖面作为主要的基础资料。(6)采用交错网格差分方法进行波动方程的数值求解,讨论引入PML边界的条件,并对其频散特性进行分析,比较了有限元方法和交错网格差分方法进行波场数值模拟的差异。相比于有限元方法,尽管交错网格差分方法得到的波场为单波波场,基于其快速的计算特性,将其作为反演的拟合度函数计算的正演是比较适合的。(7)基于交错网格差分波动方程数值正演方法和具有全局寻优的基因算法,提出了深度(里程)域的一维速度反演算法。分析和讨论了基因反演算法中各个参数对反演结果的影响,取得了最优的反演参数组合。对数值模型进行反演运算,验证了其有效性和准确性。(8)通过在湘桂铁路衡阳至永州段长达7KM的刘家冲隧道的隧道地震超前预报的实际应用,验证了本文研究的波场分离方法和反演算法的实际可行性和可信性,顺利完成了该隧道的地震超前预报工作,并形成了一套可有效应用于现场的隧道地震超前预报的数据采集和处理系统。

全文目录

摘要  6-8
Abstract  8-13
第一章绪论  13-23
  1.1 研究意义  13-14
  1.2 国内外研究现状  14-20
    1.2.1 隧道超前预报的研究现状  14-15
    1.2.2 隧道地震波场数值模拟的研究现状  15-16
    1.2.3 隧道超前预报波场分离研究现状  16-17
    1.2.4 隧道地震波场的反演研究现状  17-20
  1.3 课题来源及主要研究内容  20-21
    1.3.1 课题来源  20
    1.3.2 主要研究内容  20-21
  1.4 研究方法和技术路线  21
  1.5 课题的主要成果和创新点  21-22
  1.6 章节内容  22-23
第二章隧道地震超前预报  23-33
  2.1 隧道施工的主要地质灾害  23-25
    2.1.1 断层、断层破碎带  23-24
    2.1.2 其它地质灾害  24-25
  2.2 隧道地震预报方法  25-32
    2.2.1 负视速度法  25-26
    2.2.2 水平声波剖面技术(HSP)  26-27
    2.2.3 TSP与TGP系统  27-29
    2.2.4 TRT地震层析成像方法  29-30
    2.2.5 隧道地震CT成像(TST)技术  30-31
    2.2.6 多道隧道地震反射预报方法(TMS)  31-32
  本章小结  32-33
第三章地震波场有限元数值模拟方法  33-69
  3.1 弹性波动方程  33-38
    3.1.1 弹性本构方程  34-35
    3.1.2 Cauchy方程  35-36
    3.1.3 Navier方程  36-38
    3.1.4 波动方程的构建  38
  3.2 弹性波场的有限元数值模拟方法  38-43
    3.2.1 有限元方法  38-39
    3.2.2 弹性平面问题的有限元方法  39-42
    3.2.3 弹性动态问题到波场模拟的过程分析  42-43
  3.3 ANSYS弹性波场的数值模拟  43-59
    3.3.1 ANSYS瞬态动力学方程解法  44
    3.3.2 单元类型选择  44-46
    3.3.3 介质吸收  46-49
    3.3.4 网格剖分的频散效应  49-51
    3.3.5 震源加载方式  51-59
  3.4 弹性波场模拟的边界条件  59-68
    3.4.1 边界条件的建立  59-63
    3.4.2 边界条件在数值模拟中的实现  63-66
    3.4.3 边界条件稳定性分析  66-68
  本章小结  68-69
第四章隧道地震波场数值模拟  69-78
  4.1 边墙激发隧道地震波场  69-72
  4.2 掌子面激发地震波场  72-75
  4.3 隧道地震波场记录的时距关系分析  75-77
    4.3.1 TSP观测系统的时距关系  75-76
    4.3.2 隧道多道超前预报系统(TMS)的时距关系  76-77
  本章小结  77-78
第五章隧道反射地震波场分离数值模拟研究  78-126
  5.1 基于波动方程的纵波和横波波场分离  79-82
    5.1.1 基本原理  79-80
    5.1.2 数值模拟研究  80-82
  5.2 运用运动积法的纵波和横波波场分离  82-84
    5.2.1 基本原理  83
    5.2.2 数值模拟研究研究  83-84
  5.3 RADON变换波场分离  84-110
    5.3.1 直线型Radon(τ-p)变换的基本原理  85-87
    5.3.2 抛物Radon变换的基本原理  87-90
    5.3.3 双曲Radon变换  90-91
    5.3.4 数值模拟研究研究  91-110
  5.4 频率-波数域的二维视速度滤波  110-125
    5.4.1 基本原理  111-112
    5.4.2 数值模拟研究  112-125
  本章小结  125-126
第六章隧道地震反射波波动方程反演  126-154
  6.1 一维波动方程反演的数学模型  126-127
  6.2 交错网格有限差分方法  127-134
    6.2.1 交错网格有限差分方法  128-130
    6.2.2 边界条件处理  130-132
    6.2.3 波动方程交错网格有限差分模拟的Matlab实现  132-133
    6.2.4 交错网格有限差分(FDM)与有限元(FEM)的模拟对比  133-134
  6.3 基因反演算法  134-137
  6.4 反演算法的实现  137-143
    6.4.1 数值模型  138
    6.4.2 倾角时差校正  138-139
    6.4.3 交错网格正演模型与目标函数  139-140
    6.4.4 地震子波提取  140-143
  6.5 基因算法参数对反演影响研究  143-152
    6.5.1 初始种群对基因算法的影响  143-145
    6.5.2 种群数量的设定  145
    6.5.3 繁殖代数设定  145-146
    6.5.4 选择函数  146-148
    6.5.5 变异函数  148-149
    6.5.6 交叉函数  149-152
  6.6 数值模型反演  152-153
  本章小结  153-154
第七章实例  154-160
  7.1 工程概况和现场数据采集  154-155
    7.1.1 工程概况  154-155
    7.1.2 现场数据采集  155
  7.2 数据处理  155-157
  7.3 速度反演  157-159
  本章小结  159-160
结论  160-162
致谢  162-163
参考文献  163-171
攻读博士学位期间发表的论文及科研成果  171

隧道地震超前预报中波场分离与反演方法的数值模拟研究

内容摘要

全文目录

相似论文