一阶弹性波方程数值模拟中的混合吸收边界条件

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

三角网格谱元法地震波场数值模拟

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

上边界条件对多年冻土地温场数值模拟结果的影响分析

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

上覆硬厚岩层破断运动及断顶控制

硬厚覆岩的破断与结构将发生显著的变化。针对汝箕沟煤矿二2煤层上覆厚层石英砂岩条件,采用现场实测、数值模拟、理论分析及现场试验等方法,研究了硬厚覆岩裂隙发育特征、破断运动、矿压显现特征及断顶控制技术,为硬厚覆岩工作面灾害防治提供了依据。研究表明,硬厚石英砂岩的微观结构致密完整,呈现大面积悬空、大步距破断运动,引起强烈的支架动载,甚至导致工作面风流逆转、沟通上部采空区隐形火区;硬厚岩层破断运移后产生明显的离层空间,并在开切眼上部和工作面上部及采空区中部形成覆岩主裂隙带,且发育高度大于经验计算数值,甚至与上方采空区沟通诱发灾害;实施开切眼深孔断顶爆破及降低工作面采高,有效缩短了硬厚岩层的初次破断步距,降低了支架动载及主裂隙带高度。     

日本海域地壳结构基本特征及其地质意义

利用日本海域内完成的多条地震测深剖面及其它地球物理资料．综合分析该域地壳结构特点以及Benioff带特点．日本海域地壳结构具有陆壳、洋壳区别和两分、三分两类;Benioff带在不同方向上具有不同倾角,延伸深度差别较大;日本海开裂的动力源自太平洋板块西向俯冲所产生的多期挤压和引张应力场．     

基于混合边界条件的有限单元法GPR正演模拟

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

浅析用开氏法测定煤样中氮元素含量的主要影响因素

概述了开氏法测量煤中氮含量的重要性及其开氏法测定的原理,对测定过程中几方面主要因素的影响进行了较详细的论述分析.指出在测定过程中,必须认真、细致的考虑到各方面的影响因素,才能保证测定结果的准确性.     

基于GPU并行加速的叠前逆时偏移方法

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