波动理论基础知识及其在地震工程中的初步应用（续）

波动理论基础知识及其在地震工程中的初步应用（续）廖振鹏，杨柏坡（国家地震局工程力学研究所，哈尔滨１５００００）·ＴｏｐｉｃａｌＬｅｃｔｕｒｅ·（１０）ＥＳＳＥＮＴＩＡＬＫＮＯＷＬＥＤＧＥＯＮＷＡＶＥＴＨＥＯＲＹＡＮＤＩＴＳＰＲＥＬＩＭＩＮＡＲＹＡＰＰ...     

带阻尼的一维地震反射波成像

本文是文献［８］的继续，文中提出了地震勘探的带阻尼波动方程反演的新方法，应用此反演算法可以通过微机识别地层的波阻抗，这里把地层分成１００层，地震数据含有１０％的误差，且反演计算时间仅仅是正演计算时间的三倍．文献［１］－［４］的计算方法是基于波动方程的混合问题，而文献［８］所提出的反演算法是基于波动方程的哥西问题，许多例子充分说明后者更有效．     

地震波初至走时的计算方法综述

在地震波场中,初至波到时信息由于初至震相可追踪、易识别性,在地震学领域占有重要的位置,广泛地应用于叠前偏移、叠前速度分析、地震走时层析成像及地震定位等.本文主要介绍了四类具有代表性的计算初至波走时的方法:(1)基于高频近似射线理论方法,如最短路径方法(SPM),及修正后的最短路径方法(MSPM);(2)基于程函方程的数值解方法,如有限差分方法(FD)、快速推进法(FMM)和快速扫描法(FSM);(3)基于惠更斯原理的波前构建法(WFC);(4)基于频率域波动方程数值解法(FWQ).最短路径方法计算精度较高,稳定性较好,但其需要采用更多的网格节点,因此计算效率低;程函方程数值解法无需计算射线路径,具有计算效率高、稳定性较好、易于实现等优势,但其计算精度较低,可以通过引入高阶差分格式得到提高;波前构建法计算精度高,稳定性好,但其需要在射线网格和规则网格之间做网格转换,因此计算效率较低;频率域波动方程方法能适应任意复杂介质,但其计算精度和计算效率较低.     

北半球春季平流层最后增温过程及其年际和年代际变化特征

利用1979～2010年NCEP-DOE 2逐日再分析资料,以北半球春季平流层极夜急流核心纬带(65°～75°N)纬向平均纬向风最后一次转为东风的日期定义为春季平流层最后增温事件(SFW)的爆发日期,研究发现,SFW事件平均在4月中下旬发生,且由平流层高层向低层依次滞后,10 hPa的SFW爆发平均超前50 hPa约13天;爆发当日伴随纬向风场时间变率和行星波辐合的最大值,平流层环流实现由冬向夏的季节转换;过去32年以来SFW的爆发早晚具有显著的年际变化,最早的SFW事件发生在3月中旬,最晚的SFW事件在5月下旬才出现.合成分析表明,SFW爆发偏早(晚)年的春季,纬向风场由西风向东风的转变更为快速(缓慢),爆发前5天至爆发后5天,30 hPa纬向风减小约20 m s-1(5 m s-1),伴随的平流层行星波活动也相对较强(弱);表现在环流异常场上,SFW爆发前后平流层极区环流异常呈反(同)位相分布,表明发生较早的SFW事件主要受波强迫驱动而伴随爆发性增温,而发生较晚的SFW事件则更反映了极涡的季节变化特征.无论SFW偏早还是偏晚年,爆发后极区平流层与对流层温度异常之间均呈反位相关系,反映了SFW爆发事件中的平流层-对流层动力耦合特征.另外,在20世纪90年代中期前后,SFW爆发日期还存在明显的年代际转折,90年代中期之前SFW平均发生日期较之后约偏早11天;与之相联系的是冬末、春初行星波活动在90年代中期之前偏强,而在90年代中期之后有偏弱趋势.     

非对称型强飓风中的准平衡流特征分析Ⅱ:波动结构

利用非对称波分量分解和小波分析的方法,对准平衡动力模型下非对称强飓风内中尺度波动的空间结构和时间序列特征进行分析.结果表明,平衡流场中1波扰动占主要地位且具有涡旋Rossby波的典型结构特征,准平衡流各波数下扰动的空间分布反映了中尺度波动的混合性质;模式大气和准平衡垂直运动的全局功率谱中,超过信度检验的强波动信号中不仅包含分别表征重力波和涡旋Rossby波的高频和低频波动信号,还存在表征具有物理性质不可分特性的混合涡旋Rossby-重力波的中频波动.混合波的出现建立了不同频段波动之间的能量交换通道,其信号的变化对飓风系统的强弱变化具有一定的指示作用.非平衡垂直运动的波动强信号则主要集中在高频和低频区域,反映了在飓风强度变化情况下,与高频重力波有关的快波调整过程所引起的垂直扰动的振荡和传播.强垂直风切变对飓风内中尺度波动的切向和径向传播具有重要影响,当环境垂直风切变很强时,准平衡1波扰动在径向和切向方向上均呈"驻波"形态,随着环境垂直风切变的减弱,1波扰动以混合波波速逆基本气流传播.     

基于平面波照明的偏移成像补偿

受地下复杂构造和地震数据采集系统的影响,地震波对地下目标的照明出现不均匀,在地震数据的偏移成像中出现成像阴影.根据地震数据最小二乘偏移/反演理论,和把地震波场照明结果作为最小二乘偏移/反演中的Hessian矩阵的近似对偏移成像进行补偿的原理,提出一种应用平面波照明结果对平面波偏移成像结果进行补偿以消除偏移成像阴影的方法.这种基于平面波照明的偏移成像补偿方法相对于局部角度域的照明偏移成像补偿方法具有计算效率上的优势.     

基于波动理论的结构损伤识别研究进展

目前人们对于结构的使用安全越来越重视,结构在日常使用或灾后的损伤识别检测也变得尤为重要。近年来国内外对于波在结构中的传播理论进行了深入研究,基于波动理论的结构损伤识别方法也取得了一定进展。文章首先介绍波在介质中的传播以及在各种类型结构中的传播规律和传播特性,其次从基于波传播理论的结构损伤识别、基于Lamb波的结构损伤识别、波动理论和神经网络相结合、波动理论与其他技术或算法的融合4个方面对国内外基于波动理论对结构损伤识别方法的研究成果进行综述。     

各向异性弹性波动方程多分量联合叠后逆时偏移

为了更好地描述地下介质的各向异性特性,实现准确的波场成像,开展了各向异性介质逆时偏移研究。采用高阶交错网格有限差分法推导出各向异性弹性波动方程叠后逆时深度偏移方程,证明各向异性介质向各向同性介质转化的条件,并研究Thom son参数对弹性波场的影响,合成了各向同性介质和各向异性介质多分量零偏移距剖面,实现了多分量模拟记录联合叠后逆时偏移。计算结果表明,能够准确地实现多分量波场叠后逆时成像,使地质层位中的断层、断点等目标成像清晰准确,且偏移成像精度较高,因此各向异性逆时偏移法是一种高精度的深度域偏移成像方法,可以指导实际天然地震资料和人工地震资料的数据处理。     

一维有耗散波动方程的奇性分析与小波

对系数有强奇性（间断）的波动方程,用通常的线性化简的方法时往往会将数值小但奇性强的项略去,导致结果严重失真。利用小波变换这一工具,可以在化简时保留奇性的主要部分,使所得的结果从奇性分析的观点看来更为精确。此方法曾被用来处理系数有间断的一维波动方程,得到了与精确解的奇性主部完全一致的解．在本文中,我们改进了用小波变换作奇性化简的方法,对系数有间断的一维有耗散波动方程求得了与精确解奇性主部完全一致的解。这说明利用小波分析作奇性化简的方法对高频近似及奇性分析问题是普遍适用的．     

声波测试中振子衰减的线列阵模型

从声波在岩体中传播产生衰减的基本原理出发，建立了声波衰减振子的线列阵模型。该模型是将传播声波路径的岩体分为若干个振子，从波动概念出发，以振动原理作依据建立波动方程．为探讨岩体声波衰减机制提供了一种新的研究方法，把振动和波动直观地联系起来．单一方向上有限点的声波测试技术有可能被多方向的声测技术所代替．该模型的建立，使模态分析和参数识别方法应用到岩体声波测试技术中成为可能。