收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

基于模式识别算法的高频瑞雷波频散曲线非线性反演研究

宋先海  
【摘要】: 利用瑞雷波推断地下介质结构和岩土力学参数,已在研究地球内部结构、地壳及地幔物质组成、大地构造和地震灾害预测等区域及全球地震学领域、近表面地球物理工程和超声无损检测等领域中获得了广泛应用。特别是在近表面地球物理研究中,由于高频瑞雷波勘探具有快速、轻便、无损、衰减小、抗干扰能力强、浅层分辨率高且不受各地层速度关系的影响等优点,已在实际生产和科研中得到了广泛应用。由于瑞雷波理论的复杂性及重要性,近十年来,参与瑞雷波勘探技术研究的人越来越多。尤其是最近几年,国内外已掀起了瑞雷波勘探技术的研究热潮。可以预料,瑞雷波勘探将会成为本世纪浅层或超浅层地球物理勘察和工程质量无损检测的重要手段之一。 瑞雷波频散曲线的反演解释是瑞雷波勘探技术的关键步骤之一,通过反演频散曲线可直接进行地质分层和获取重要的岩土力学参数——横波速度。和大多数其它地球物理反演问题一样,瑞雷波频散曲线的反演是一典型的高度非线性、多参数、多极值的地球物理反演问题和最优化控制问题。目前现有的局部线性化方法,其反演结果的可靠性严重依赖于初始模型的选取,只有当初始模型接近真实模型时才能获得较好的反演结果,否则将很容易陷入局部极小值中得出错误的地质解释;同时该方法还需要计算偏导数信息,雅可比矩阵的求取精度也将直接影响反演结果的质量。现有的非线性全局优化算法,如遗传算法(GA)和模拟退火(SA),用于瑞雷波反演是相当有用的,但由于它们本质上是对可行解空间进行随机取样,执行的是非确定性的随机搜索过程,全局最优解的获取需要多次成功的独立试验后(例如20次)加以平均,且并不一定保证每次反演总能搜索到全局最优解。所以计算工作量大、收敛速度慢、计算时间长,在实际应用中受到了限制。 鉴于此,本文围绕瑞雷波反演技术中存在的主要问题,结合代数拓扑学、应用数学和最优化理论等前沿交叉学科的最新研究成果,研究了基于模式识别算法的高频瑞雷波频散曲线非线性反演方法。本文研究的模式识别算法是一类确定性的全局优化算法,其反演机理可以概括为:(1)在每一次迭代时,算法通过当前模式搜索足够的方向以确保一个全局最优下降方向最终能被识别。(2)算法采用合理的反向追踪机制可对试验步长进行有效控制,避免了模式识别过程中过长或过短的试验步长,从而更加有利于算法向最有前景的全局最优点运动。(3)算法特有的隐含弹性网格结构,使其具有隧道效应特征,从而可使算法有效地逃逸局部极小值的吸引。因此,模式识别算法可有效地克服上述瑞雷波反演问题。 本文以“基于模式识别算法的高频瑞雷波频散曲线非线性反演研究”这一策略为核心问题进行了深入地研究与探讨,取得了以下主要研究成果: 1、调研、总结了瑞雷波勘探技术和模式识别算法的研究现状;指出了当前用于瑞雷波反演的局部线性化方法和非线性全局优化算法所存在的不足;分析了模式识别算法用于高频瑞雷波频散曲线非线性反演的可行性和优越性。 2、深入研究了一种新的确定性全局优化算法——模式识别算法。深入研究并总结了模式识别算法的基本原理、反演机制、算法实施及反演的一般流程;针对其在实际反演中的具体情况,探讨了算法实施过程中的一些关键技术问题。如初始网格步长和模式的选取、扩展因子和收缩因子的选取、完全预测和完全识别策略的执行,以及算法终止准则等问题。 3、执行了有效的数值模拟测试。本文首先通过一个高度非线性、多参数、多极值的数学函数“Rastrigin”(有100个极大值和90个极小值)具体测试、分析了模式识别算法的执行效率,并通过反复试验获得了最优的反演参数。然后用另外18个著名数学函数进一步检验了反演参数选取的合理性、算法的有效性和程序的正确性。这些测试函数不仅具有高度非线性、多参数、多极值的特点,而且目标函数形态各异,因此能有效地模拟高频瑞雷波频散曲线非线性反演问题。测试结果表明:模式识别算法是一种不依赖于初始模型的确定性全局优化算法,不仅计算精度高,而且由于算法执行的是确定性的模式识别过程,最优解一次计算便可完成,而不需要多次独立试验后取平均值,所以计算时间短。 本文数值模拟测试后建议:在最优化反演过程中,初始网格步长取为1、扩展因子取为2、收缩因子取为0.5、选取最大正基模式、同时采用完全预测和完全识别策略,则在多次模式识别后最终迭代网格步长将会近似收敛到零,此时一般将会获得较好的反演结果。 4、实现了快速稳定的多模式瑞雷波频散曲线正演模拟算法。由于瑞雷波反演的大部分时间都将用在正演模拟计算上,因此正演是反演的前提和瓶颈,成为提高反演速度和反演精度的关键因素之一。而以往的瑞雷波频散曲线正演模拟算法(如Haskell算法),不仅计算速度慢,而且易出现高频数值精度丢失和高频数值溢出等不稳定现象。为此,本文分析并实现了快速矢量传递算法,该方法是建立在快速δ矩阵算法基础之上的,不仅表达形式简单、易于编程实现,而且能有效地克服上述问题。例如,利用本文给出的模型测试结果表明:同等条件下快速矢量传递算法比Haskell方法计算速度约快40%左右;Haskell算法在频率4000Hz左右出现了高频数值精度丢失问题;Menke算法在频率24kHz左右即出现了高频数值溢出问题:而快速矢量传递算法在频率高达40kHz仍未出现高频数值精度丢失和高频数值溢出等不稳定现象。该算法的高频数值稳定性不仅对近表面研究具有重要的意义,而且对超声表面波在无损检测中的实际应用(如混凝土质量检测、裂缝探测、超声探伤等)也具有重要的研究价值。 5、研究并探讨了多模式表面波叠加耦合机理。多层介质正演模拟结果表明瑞雷波往往表现为多个模式,在实际工作中利用表面波谱分析方法提取的频散曲线往往是一条,而不是理论计算的多条,且在某些公路型等特殊地质结构中频散曲线有时还会出现“之”字形回折现象。为此,本文分析并实现了多模式瑞雷波频散曲线叠加计算方法,并通过三个常见的典型地质模型(递增模型、软夹层模型和硬夹层模型)进行了叠加频散曲线的正演模拟计算,探讨了多模式表面波的叠加耦合机理。其重要的理论与实际意义在于: (1)能够在理论上对多模式表面波共同叠加作用的耦合机理有一个深入地了解和认识,并能对实测频散曲线中出现的“之”字形回折现象的成因予以解释。利用本文模型研究结果表明:在大多数情况下,基阶波与其它各模式波相比,其能量仍占主导地位,叠加频散曲线实际上就是基阶波频散曲线。而对于某些地下含有明显的软夹层或硬夹层等特殊地质结构,基阶波的能量则可能在某些频段(例如高频段40~100Hz或中频段13~23Hz)处于次要地位,此时叠加频散曲线是基阶波与各高模式波在相应频段相互叠加共同耦合作用的结果,即会出现“之”字形回折现象,从而为前人在实践中总结的利用频散曲线的拐点进行地质分层提供了重要的理论依据。 (2)在已知地质信息较少的情况下,通过叠加频散曲线正演模拟计算能够有利于指导我们对实测频散曲线进行合理的地质分层,尤其是对于“之”字形频散曲线。 (3)为多模式表面波的反演提供了另一个有潜力的、可能有效的反演途径。使用叠加频散曲线进行反演可能会避免多模式表面波反演中易犯的正常模式误识别的错误和需要指定各模式数据权重的困扰。 6、成功进行了大量理论模型试算。(1)为了进一步检验模式识别算法程序的正确性和该算法反演瑞雷波频散曲线的有效性(计算精度和计算时间),本文针对瑞雷波在实际近表面应用中经常遇到的典型地质结构,设计了大量有代表性的地质模型,并用模式识别算法进行了理论模型试算,测试了该算法反演瑞雷波频散曲线的可靠性。(2)由于实测数据不可避免地会含有噪声,为了验证算法的稳定性,作者还对其抗噪能力进行了检验。(3)考虑到实际表面波相速度反演还常受许多因素干扰,本文进一步探讨了频带范围、纵波速度和密度中的误差、频点数和初始模型、层数和层厚度对模式识别算法执行效率的影响。(4)为了进一步演示模式识别算法的优越性,本文还将其与两种常用的非线性全局优化算法(GA和SA)进行了对比测试。本文理论模型试算结果表明:模式识别算法是一种不依赖于初始模型的确定性全局优化算法,不仅计算精度高、具有较强的抗干扰能力(稳定性好),而且由于算法执行的是确定性的模式识别过程,最优解一次计算便可完成,所以计算速度快。例如,在同等条件下,模式识别算法计算速度一般比GA快90倍左右,比SA快80倍左右,能有效地用于高频瑞雷波频散曲线非线性反演。 7、成功基于模式识别算法执行了多模式表面波联合反演研究。考虑到在某些含低速夹层等特殊地质结构中,高阶波较基阶波在某些频段(如高频段)有时可能会更加发育,此时利用多模式表面波联合反演是必要的。为此,本文首先通过一个典型的地质模型,分析了多模式瑞雷波的敏感性;然后通过模式识别算法进行了多模式表面波联合反演研究,校验了理论模拟结果的正确性,并探讨了充分利用多阶表面波联合反演的优点。研究结果表明:(1)对于本文给出的地质结构,深层介质基阶波敏感性一般仅分布在11Hz左右较窄的频带范围内,且在深层(例如第三层和第四层)还会出现敏感性重叠现象,从而降低了基阶波深层敏感性和分辨率,易使反演出现多解性。(2)高阶波敏感性则分布在较宽的频带范围内,且基本不会出现深层敏感性重叠现象。(3)随着模数的增高,表面波峰值敏感性会向高频方向移动,且敏感性分布得更宽、相互分离得更好。(4)高阶波较基阶波对地质结构具有更强的敏感性,充分利用高阶波不仅可使反演过程稳定、提高横波速度的反演精度和计算结果的可靠性,而且还能够提高模型纵向分辨率。(5)利用的高阶波越多,计算结果的可靠性就越高,分辨能力就越强(高阶模式数据比低阶模式数据具有更强的分辨能力)。 8、执行了有效的典型实例分析研究。通过一个来自某高速公路路基上的典型实例反演研究,进一步检验了模式识别算法的有效性和实用性。 9、实现了与本课题相关的重要程序源代码。具体包括:多模式瑞雷波正演模拟算法(如Haskell算法、Schwab-Knopoff算法、Menke算法和快速矢量传递算法):多模式表面波叠加频散曲线正演模拟算法;广义模式识别算法(GPS),遗传算法(GA),模拟退火(SA),人工神经网络(SWIANN,Surface Wave Inversion by Artificial Neural Networks);粒子群优化算法(PSO);Occam算法;f-k变换频散曲线提取算法;基于相移法的多模式表面波频散曲线成像工具(Dispersion Curves Imaging Tools,DCIT)及其它相关程序,并编写了相应的显示程序模块。 本文的主要创新点为: 1、在国内外首次将基于最大正基模式的广义模式识别算法(GPS)应用于高频瑞雷波频散曲线非线性反演研究,并获得了成功。本文大量典型的数值试验、理论模型试算、抗噪能力检验以及与其它全局优化算法对比分析等均证明了GPS算法具有明显的优势,不仅反演精度高,而且计算速度快。本文的研究同时也为地球物理非线性反演和其它全局优化控制领域提供了另一个重要的研究方向——开发研究确定性的全局优化算法,以提高算法的执行效率。 2、在国内外首次提出了基于GPS算法的多阶表面波联合反演策略,并通过广义模式识别算法证明和演示了多模式瑞雷波联合反演不仅能使反演过程稳定、提高横波速度的反演精度和计算结果的可靠性,而且还能够提高模型纵向分辨率;同时利用的高阶波越多,计算结果的可靠性就越高,分辨能力就越强。 3、深入分析了多模式表面波的相互叠加耦合机理,从而为前人在实践中总结的利用频散曲线的拐点进行地质分层提供了重要的理论依据,并为多模式表面波联合反演提供了另一个有潜力的发展方向。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 宋先海,肖柏勋,余才盛,顾汉明,张学强,刘江平;瞬态瑞雷波反演横波的SVD算法及其应用[J];地质与勘探;2004年01期
2 肖翔,徐果明,朱良保;面波频散曲线快速追踪算法[J];地震地磁观测与研究;2004年03期
3 苏玉杰,薛桂玉;瑞雷波法在滑坡滑动面勘探中的应用[J];中国农村水利水电;2005年03期
4 王延坤;刘江平;王万合;甘文兵;;对缺失频段的瑞雷波频散曲线反演的误差分析[J];工程地球物理学报;2006年01期
5 马志杰;董连成;;Love波相速度频散曲线敏感性比较[J];山西建筑;2008年12期
6 徐华全;顾勤平;刘远志;;多道瞬态面波法在覆盖层调查中的应用[J];内蒙古石油化工;2008年17期
7 董月昌;鲍文刚;张维冈;刘喜武;;超软地层的声波特征研究[J];内蒙古石油化工;2008年24期
8 禹建功;吴斌;;正交各向异性空心圆柱体的周向SH波[J];声学技术;2009年03期
9 孙广开;焦阳;李光海;赵玉博;朱英奎;;超声导波管道缺陷检测数值模拟[J];河北工业科技;2010年01期
10 胡明顺;潘冬明;李娟娟;陈涛;夏暖;;基于频散曲线合成面波地震记录的方法[J];煤田地质与勘探;2010年02期
11 金淼;王水强;;瑞雷面波在探测驳岸基础中的应用[J];上海地质;2010年01期
12 杨小慧;李德春;于鹏飞;;煤层中瑞利型槽波的频散特性[J];物探与化探;2010年06期
13 魏文涛;;多道瞬态瑞雷波法在铁路路基检测中的应用[J];西铁科技;2010年04期
14 常永红;;面波在水库坝基检测中的运用实例[J];新疆有色金属;2011年03期
15 李长征,肖柏勋;含层状液体介质瑞雷面波频散曲线正演分析[J];工程地球物理学报;2005年03期
16 吴加和;沈骥千;陈红萍;徐永清;倪金才;李德华;杨素霞;;地裂缝发育区面波勘探技术研究[J];能源技术与管理;2006年06期
17 朱翔鹏;张学强;严哲;;瑞雷波在探测地下防空洞的应用研究[J];工程地球物理学报;2007年01期
18 邵广周;李庆春;梁志强;;液体表层层状介质导波频散曲线研究[J];地球物理学报;2007年03期
19 周波;马志杰;;Love波相速度频散曲线计算[J];山西建筑;2009年17期
20 秦臻;姚姚;张才;宋建勇;;频率域速度谱中提取面波的多模频散曲线方法[J];人民长江;2009年11期
中国重要会议论文全文数据库 前10条
1 阎石;何彬彬;赵乃志;;管结构导波频散曲线绘制与试验验证[A];第20届全国结构工程学术会议论文集(第Ⅱ册)[C];2011年
2 张海澜;王秀明;林伟军;;井孔频散曲线的奇异特征[A];2008年全国声学学术会议论文集[C];2008年
3 姚华建;肖翔;徐果明;;基于图像处理的相速度频散曲线自动追踪方法[A];2001年中国地球物理学会年刊——中国地球物理学会第十七届年会论文集[C];2001年
4 鲁来玉;张碧星;鲍光淑;;瑞利波“之”字形频散曲线分析[A];2000年中国地球物理学会年刊——中国地球物理学会第十六届年会论文集[C];2000年
5 董连成;陶夏新;;从地脉动台阵三分量观测中推算工程场地勒夫波频散曲线[A];第八届全国振动理论及应用学术会议论文集摘要[C];2003年
6 王永华;;瞬态瑞雷波勘探技术在西藏公路工程勘察中的应用[A];第八次全国岩石力学与工程学术大会论文集[C];2004年
7 魏继祖;单娜琳;梁芳敏;;面波多阶模态频散曲线的BP神经网络反演问题研究[A];中国地球物理·2009[C];2009年
8 周阳;黄真萍;叶琪;张远帆;;瞬态瑞雷面波法研究及其在边坡工程中的应用[A];2010年全国工程地质学术年会暨“工程地质与海西建设”学术大会论文集[C];2010年
9 鲁来玉;张碧星;汪承灏;;瑞利波高阶模式的实验和反演研究[A];中国声学学会2006年全国声学学术会议论文集[C];2006年
10 祁生文;孙进忠;万志清;;瞬态瑞雷波勘探方法的一点改进[A];第七届全国岩石动力学学术会议文集[C];2000年
中国博士学位论文全文数据库 前10条
1 刘雪峰;层状介质中瑞利面波多模式性质及其在正反演中的应用[D];哈尔滨工业大学;2011年
2 卢建旗;多道面波分析方法及其应用研究[D];中国地震局工程力学研究所;2013年
3 汪利民;起伏地表三维高频瑞雷面波传播特性研究[D];中国地质大学;2013年
4 张逍越;面向结构健康监测的Lamb波力学建模研究[D];南京航空航天大学;2012年
5 邵广周;多阶模式瑞利波频散特征与反演研究[D];长安大学;2009年
6 杨天春;瑞利波“之”字形频散与道路结构频散曲线的正演研究[D];中南大学;2004年
7 白朝旭;地震波场中瑞雷波研究及在工程建设中应用[D];吉林大学;2008年
8 肖柏勋;高模式瑞雷面波及其正反演研究[D];中南大学;2000年
9 李晶;面波在地震波场中的特性研究及其应用[D];成都理工大学;2006年
10 宋先海;基于模式识别算法的高频瑞雷波频散曲线非线性反演研究[D];中国地质大学;2008年
中国硕士学位论文全文数据库 前10条
1 万远收;瞬态瑞利波法频散曲线的仿真研究[D];华中科技大学;2013年
2 邓乐翔;瑞雷波场正演模拟及频散曲线的提取[D];长安大学;2010年
3 陈艳艳;瑞雷波频散曲线提取及快速反演分析应用研究[D];哈尔滨工业大学;2010年
4 周文宗;溶洞地质条件下瑞雷波频散曲线正反演研究[D];中南大学;2013年
5 徐凯;超声Lamb波信号传播特性研究[D];南京信息工程大学;2012年
6 潘世旻;模态分离在瑞雷波频散曲线提取中的应用研究[D];哈尔滨工业大学;2011年
7 毛承英;基于改进遗传算法的瑞雷波频散曲线反演[D];中南大学;2010年
8 秦波;基于Rayleigh波理论的地基土动参数反演[D];北京交通大学;2011年
9 刘亚亮;超声表面波(LSAWs)法表征low-k薄膜杨氏模量的快速实现[D];天津大学;2010年
10 陈强;基于Lamb波的薄板钢结构损伤识别研究[D];哈尔滨工程大学;2013年
中国重要报纸全文数据库 前2条
1 白文起 刘兰宁 赵海燕;为了千年古镇的地质安全[N];地质勘查导报;2008年
2 白文起 刘兰宁 赵海燕;为了移民的居住安全[N];中国矿业报;2008年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978