基于正交贴体网格的VTI介质地震波数值模拟     

刘志强  李钢柱  黄磊  许磊  李文宝 《地球科学与环境学报》2023,(2):168-178

为了解决规则网格有限差分法在处理起伏地表自由边界条件时需进行复杂坐标转换和插值运算的问题,以及对起伏地表进行阶梯状离散近似所产生的虚假散射波问题,将计算流体力学中的正交贴体网格生成技术引入到起伏地表下VTI介质的网格剖分中,采用基于同位网格的DRP/opt MacCormack有限差分法模拟起伏地表下VTI介质中的地震波场。网格的正交性使实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明,本文提出的网格剖分法能够有效消除阶梯状网格产生的虚假散射波,从而提高起伏地表下VTI介质地震波数值模拟精度。  相似文献   

黏弹性介质中瑞雷波有限差分数值模拟与波场分析          下载免费PDF全文

王常波  田坤  刘立彬  李国磊  张学涛 《CT理论与应用研究》2019,28(2):153-166

瑞雷面波经常被用来反演地表浅层横波速度,受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质,而实际地层更接近黏弹性介质,对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质,并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟,再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型,并给出广义标准线性固体的具体算例,实施自由表面条件时采用声学-弹性边界近似法,通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面,采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算,数值解与解析解的对比验证了本文方法的准确性与有效性,正演结果的对比表明波场尤其是面波频散会受黏弹性影响,因此有必要在面波勘探中考虑黏弹性因素。   相似文献   

结构混合试验边界协调技术与应用     

张钰 《地震科学进展》2019,(10):73-73

  相似文献   

被动桩水平位移荷载施加位置的探讨     

上官士青  杨敏  李卫超 《岩土力学》2015,36(10):2934-2938

被动桩是指一种由于土体水平位移而发生挠曲变形的桩。目前不论是室内试验还是数值计算,常在模型边界对土体施加水平位移荷载（位移边界条件）形成土体的位移场,用来研究被动桩的变形特征,但不同研究者采用的水平位移施加位置（边界位置）与被动桩的距离存在较大差异。结合一个工程案例,采用有限单元法对水平位移加载位置进行了单变量参数分析。研究表明,随着水平位移加载位置与被动桩距离的增加,桩身变形显著减小。这表明了施加位移边界后地基中存在明显的应力扩散。当该距离小于5.5倍桩径时,在模型边界施加同等大小的水平位移所需应力显著增加。同时讨论了水平荷载加载位置与被动桩距离的合理范围,认为应在5～8倍桩径范围之内。  相似文献   

一阶弹性波方程数值模拟中的混合吸收边界条件   总被引：2，自引：1，他引：1       下载免费PDF全文

任志明  刘洋 《地球物理学报》2014,57(2):595-606

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.  相似文献   

三角网格谱元法地震波场数值模拟   总被引：5，自引：0，他引：5  

刘有山  滕吉文  徐涛  刘少林  司芗  马学英 《地球物理学进展》2014,(4)

谱元法结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度,是近年来发展的一种重要的地震波场数值模拟方法.经典的谱元法采用四边形(六面体)网格,利用一维Gauss-Legendre-Lobatto(GLL)积分的张量积得到对角的质量矩阵,以大大提高计算效率,但是四边形(六面体)网格不能够灵活地刻画复杂的几何模型的弯曲界面.为此,在谱元法中引入三角形(四面体)网格到二维(三维)是十分必要的.不同于经典的谱元法,在非结构化网格中不能使用GLL积分的张量积,使得非结构化网格的谱元法的实现存在着诸多的困难.目前,比较流行的三角网格谱元法,通过使用KoornwinderDubiner(KD)正交多项式,并正交化这些KD多项式构建基函数,同时利用重合的插值节点和积分节点以获取对角的质量矩阵;它所使用的积分点为优化的点集——Fekete点,且这些积分点能与四边形网格完全耦合.相比于四边形,三角网格谱元法能显著提高复杂模型的描述能力,对起伏地表模型有很大优势.本文引入高效的最佳匹配层(PML)吸收边界条件,并通过数值试验将三角网格谱元法与经典的谱元法进行对比研究.相比于经典的谱元法,三角网格谱元法显著缺点为较低的计算精度.对于7阶谱元,为了能够精确地模拟面波,三角网格谱元法需要在每个最短的面波波长内至少有11个采样点,然而经典的谱元法仅需4个采样点,并且前者所需的内存量约为后者的5.5倍.  相似文献   

上边界条件对多年冻土地温场数值模拟结果的影响分析   总被引：1，自引：1，他引：0  

曹元兵  盛煜  吴吉春  李静  宁作君  胡晓莹  冯子亮  王生廷 《冰川冻土》2014,36(4):802-810

以玛多地区多年冻土为背景,建立多年冻土地温场的数值计算模型,以不同的方式考虑近 60 a 来的气温变化构成不同的上边界条件,通过模型计算分析不同上边界条件下的不同时期温度场、未来冻土退化特征.结果表明：在上边界条件中采用气象站实测近60 a波动温度值和采用近60 a平均恒定值时,浅层冻土地温差异明显,且越浅层地温与越近时间的上边界条件相关.预测未来100 a冻土地温变化趋势发现,相同升温速率和升温初始温度条件下,上边界采用实测60 a波动温度值对冻土退化过程影响较小;升温初始温度值提高到与趋势线衔接后,冻土退化起始时间从约第45年提前到约第20年;60 a实测温度和升温初始温度值均提高到与其初始温度场上边界条件衔接后,冻土退化起始时间从约第20年提前到约第15年;冻土退化从开始到完全退化经历时间为25 a左右.  相似文献   

MT有限元模拟中截断边界的影响     

汤井田  薛帅 《吉林大学学报(地球科学版)》2013,43(1):267-274

利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。  相似文献   

基于混合边界条件的有限单元法GPR正演模拟   总被引：6，自引：6，他引：0       下载免费PDF全文

冯德山  陈承申  王洪华 《地球物理学报》2012,55(11):3774-3785

从Maxwell方程组出发,推导了探地雷达(GPR)有限元波动方程.阐述了透射边界条件和Sarma边界条件的原理,推导了这两种边界条件的理论公式;通过在衰减层内加入过渡带优化了Sarma边界条件的加载方法,压制了介质区和衰减层交界面处的人为反射.考虑到透射边界条件与Sarma边界条件不同的理论机制,提出了一种结合透射边界条件和Sarma边界条件的混合边界条件,它利用Sarma边界条件对到达边界区域的GPR波能量衰减功能和透射边界对GPR波能量的透射功能,使GPR波经过Sarma边界条件的衰减吸收后,再通过透射边界条件将剩余能量透射出去,集成了二者的优势.并以二维均匀模型中的中心脉冲激励源方式为例,通过Matlab程序实现,以GPR的全波场快照的直观方式,对比了有、无边界条件及不同边界条件对人工截断边界的处理效果,说明了该混合边界条件对到达截断边界处的GPR波的处理优于单一边界条件.最后,以基于混合边界条件的有限单元法对两个典型的GPR地电模型进行了正演模拟,指导了GPR数据处理与工程实践.  相似文献   

基于GPU并行加速的叠前逆时偏移方法   总被引：1，自引：1，他引：0  

石颖  陆加敏  柯璇  田东升  王菲 《东北石油大学学报》2012,36(4):111-115,11

为了提高复杂地下介质的成像精度和偏移算法的计算效率,提出可高效对地下复杂构造进行准确成像的GPU加速叠前逆时偏移方法.该方法采用双程声波方程进行波场延拓,突破倾角限制,借助于高阶有限差分方法实现叠前逆时偏移成像;利用GPU(Graphic Processing Unit)并行加速技术对波场延拓和成像进行计算,相比于传统算法,其计算效率有较大提高,可以解决叠前逆时偏移算法计算量过大问题;在获取波场信息过程中,也采用随机边界条件,实施以计算换存储策略,解决逆时偏移计算中的海量存储问题.模型测试结果表明,该方法能够高效和高精度地对地下复杂地质体成像.  相似文献