存在流体介质夹层时瑞利波频散曲线的完整求解

应用介质分界面上位移应力连续性的变化关系,采用δ矩阵法和固体层状介质中的瑞利波频散方程,推导了含流体介质夹层情况下瑞利波的频散方程。利用推导的频散方程对理论模型进行模拟计算,并与以往文献计算方法进行比较,模拟结果吻合较好,验证了公式的正确性;同时,模拟计算结果也说明了本文公式的普遍适用性和计算结果的完整性。根据文中对4个理论模型的模拟结果来看,与固体层状介质的情况相比,当存在流体夹层时,瑞利波频散曲线模式中多出了一个文中定义的0阶模式,在今后的多道瞬态瑞利波法应用中,应该可以将其看成是流体介质夹层存在的标志性曲线。     

夹层厚度对瑞利面波频散曲线的影响

文章以弹性波横波速度为参数,利用瑞利面波频散曲线的快速标量传递算法进行正演模拟计算,系统研究软弱层厚度对瑞利面波基阶模频散曲线的影响,分析和总结了频率与相速度（f-c）、频率与波数（f-k）、深度与相速度（h-c）等瑞利面波基阶模曲线随软弱层厚度变化的规律。     

软弱夹层速度对瑞利面波频散曲线的影响

文章以弹性波横波速度为参数,利用瑞利面波频散曲线的快速标量传递算法进行正演模拟计算,系统研究软弱层速度对瑞利面波基阶模频散曲线的影响,分析和总结了频率与相速度（f-c）、频率与波数（f-k）、深度与相速度（h—c）等瑞利面波基阶模曲线随软弱层速度变化的规律。     

Matlab环境下瑞利波有限差分正演与曲线绘制

以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。      

多层地基非轴对称Biot固结的理论解

提出了一种新方法来推导多层地基非轴对称Biot固结问题的传递矩阵,相应的理论推导工作量较少。基于柱坐标系下非轴对称Biot固结的基本方程,通过引入中间变量,并对坐标 进行Fourier级数展开,对时间 和坐标 进行Laplace-Hankel变换,得到了6×6阶和2×2阶的两组常微分方程;然后,两组常微分方程进行关于 的Laplace变换和逆变换,得到了单层地基非轴对称Biot固结问题的传递矩阵;结合边界条件和连续条件,运用传递矩阵法得到多层地基非轴对称Biot固结问题在Laplace-Hankel变换域内的解,通过Laplace-Hankel逆变换得到了该问题物理域的真实解答。编制了计算程序,并进行了数值计算与分析。该方法具有计算速度较快的特点,方便工程应用。     

软弱夹层底层速度对瑞利面波频散曲线的影响

文章利用快速标量传递算法进行瑞利面波频散曲线的正演计算,系统研究软弱夹层底层速度对瑞利面波频散曲线的影响,分析和总结了频率与相速度(f-c)、频率与波数(f-k)、深度与相速度(h-c)等瑞利面波基阶模曲线随软弱夹层底层速度变化的规律.     

锚索抗滑桩内力计算的传递矩阵法

为了改进现有锚索抗滑桩计算中存在的不足,将传递矩阵算法应用于锚索抗滑桩结构内力计算。推导了适合抗滑桩计算的场矩阵和点矩阵,能够完成普通抗滑桩、桩顶作用集中荷载以及桩身有多个受力分段的抗滑桩内力计算,再结合具体的不同点矩阵形式,可以完成预应力锚索抗滑桩内力计算。根据预应力锚索抗滑桩的工作原理,给出了不同受力阶段的传递矩阵计算式。最后,通过具体的计算实例对预应力锚索抗滑桩不同受力阶段的内力和位移进行了计算分析,结果表明,在锚索预应力作用下,受荷段出现较大的剪力和弯矩,位移较小;进入变形协调阶段后,滑面处和锚固段的剪力明显增加,锚固段出现较大弯矩,桩顶出现向前的位移。     

基于传递-反射矩阵法的层状公路结构动力响应研究

利用传递-反射矩阵方法和叠加原理推导了双轮胎振动荷载作用下弹性层状公路结构的动力响应解,并利用Hankel数值逆变换技术进行了数值求解。弹性层状结构的传递-反射矩阵同时具有单个矩阵尺寸小、矩阵元素只含有负指数项2个优点,计算结果不受土层单元厚度与荷载频率大小的影响,可以大幅提高计算效率和精度。通过将数值计算与室内模型试验结果进行比较,验证了方法的正确性。在此基础上,分析了面层厚度及模量对层状公路结构动应力、动位移的影响程度和规律,得出了一些有意义的结论。研究结果可为公路路面、路基结构设计提供参考。     

最优分割算法改进与程序实现

本文论述对“有序量最优分割法”中,类直径矩阵 D 的计算公式推导、分解与交换过程。采用“递推”紧凑算法。充分利用中间性数据的传通,巧妙设计计算流程,比原算法大大减少计算量、节省计算机内存容量,缩短计算机运行时间4倍以上,实现 PC—1500、PB—700 等袖珍型计算机应用程序软件。同样可移植到一切能使用BASIC 语言的微、小型计算机上实现。     

一种新型有限差分网格剖分方法在大地电磁一维正演中的应用

介绍了有限差分算法在大地电磁测深法中的应用,推导了一维层状介质和连续介质模型下的有限差分算法,提出了一种新的网格剖分方法,并通过两层和三层介质模型的正演计算,验证了算法的正确性以及网格剖分方法的合理性。设计了一种连续介质模型,与两种不同程度近似的多层层状介质模型进行了正演结果的对比,指出了研究连续介质模型正演的必要性。对于连续介质模型,计算中剖分得越详细,其结果越接近真实值,而随着误差的逐渐减小,精确度的提高不再明显,但资源与时间的消耗将大幅增加。因此,无限制地减小网格步长是不必要的。同时,为了兼顾误差大小与网格步长的关系,采用了一种将一次网格剖分进行二次差分计算的方法得到新的结果,使计算结果明显得到了改善。     

计算层状介质中轴对称柱面瑞利面波频散函数的δ矩阵法

笔者将δ矩阵法应用于计算层状介质中轴对称柱面瑞利面波的频散函数,得到了六阶δ矩阵法、五阶δ矩阵法、快速δ矩阵法3种方法,很好地解决了高频数值精度丢失问题以及高频数值溢出问题,并提高了计算速度,数值计算及工程应用验证了上述方法的有效性,且表明了这些方法也完全适用于平面瑞利面波频散曲线的求取问题。     

基阶和高阶瑞雷波频散曲线反演对比分析

通过瑞雷波频散曲线反演,可以建立近地表横波速度模型。瑞雷波频散曲线具有多个模态,一般情况下,基阶模态是具有最大能量的模态。目前对瑞雷波频散曲线的反演大多数是基于基阶模态,对于高阶模态的反演比较少。通过对基阶和高阶瑞雷波频散曲线反演对比分析,认为高阶模态和基阶模态频散曲线反演有比较好的一致性趋势,相对而言高阶模态频散曲线反演的横波速度比基阶要高;高阶和基阶模态能量加权反演,不仅保留了频散曲线基阶模态能量强,勘探深度深的特点,而且对浅层低速层也具有更好的响应,可有效提高反演的精度。     

软夹层介质瑞利波合并化频散曲线模拟分析

针对含有软夹层介质的地质模型,研究其瑞利波合并化频散曲线随夹层厚度和埋深的变化规律,分析产生“之”字形现象的根本原因。利用传递矩阵方法计算瑞利波多模式频散曲线,以及它们在自由表面的垂直位移强度,分别按能量最大法和相对位移权值法获得合并化的频散曲线,从而获得软夹层厚度或埋深变化的合并化频散曲线。由计算结果分析可知,两种方法所获得的频散曲线的整体变化规律基本相同,“之”字形随软弱夹层厚度的增加或埋深的变浅而增多;按能量最大法合成的频散曲线,“之”字形现象的出现是由于模式之间的跳跃所引起;按位移权值法计算合成频散曲线时,产生“之”形的原因是由于从单个模式的位移占优跳至多模式的位移共同占优所致。      

三层层状介质中的多导波模式及其频散和位移特征

采用传递矩阵法求解三层均匀各向同性的层状弹性固体介质中瑞利波频散曲线及自由表面的垂直位移强度，对不同情形下多导波模式的频散曲线和位移曲线特征进行了理论和数值模拟研究。声源中心频率变化时的计算结果表明，位移的激发强度与声波的频率存在密切的关系；数值模拟结果表明，在不同模型情况下接收到的导波模式是不同的。因此，在对野外勘探中，震源频率的选择和资料解释应引起注意。     

基于附加层法的非规则剖面的瑞利波频散特性研究

对于某些非规则水平层状地质剖面,采用传递矩阵法获得的瑞利波频散曲线可能会出现截止或截断现象,但采用附加层的方法可避免这些现象的发生。通过采用附加层法前后相速度与激发强度大小的对比,来论证附加层法的正确性。对3个典型的非规则剖面进行模拟计算,利用传递矩阵法算得添加附加层前后瑞利波频散曲线,以及某些特征点处位移分量随深度的变化情况。通过对计算结果的对比分析可知,附加层法不仅可计算出地质剖面的实数导波,还可算出其泄漏模式波;在固定的频率点上,附加层的埋深越大,则添加附加层后算得的激发强度大小与原模型的误差越小;当附加层的埋深固定时,添加附加层后计算出的位移曲线与原位移曲线之间的误差随着对比点频率的增大而减小。可见,今后根据实际工作的精度要求,可用附加层法计算某些非规则剖面的瑞利波频散曲线。     

求取道路结构型地层瑞利波频散曲线的方法

求取道路结构型地层的瑞利波频散曲线,过去一般是采用刚度矩阵法、有限单元法等方法在复数域进行求解.作者在吸取前人某些研究思路的基础上,提出了添加附加层的方法,使原复数域的求解问题转化到实数域进行,然后用二分法进行求解.附加层法提高了正演计算的速度,且模拟计算结果表明了新方法的可行性.