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

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

鄂尔多斯地块地壳内S波衰减特征的研究

采用分布于鄂尔多斯地块及其周缘的32个地震台站记录到的19次地震事件共487条数字化波形记录, 利用遗传算法反演得到了整个鄂尔多斯地块地壳内S波的三段几何扩散参数、S波非弹性衰减Q值随频率变化的关系式、32个台站的场地响应以及19次地震的震源参数. 计算结果表明, 鄂尔多斯地块地壳内S波非弹性衰减在1 Hz处的Q值比地质构造活跃地区的结果要大得多. 32个地震台站的场地响应在高频部分除个别台站外基本没有表现出明显的放大效应,但在低频区间部分台站之间有一定的差异. 19次地震事件的地震矩对数值与震级ML有很好的线性关系.      

弹性波逆散射微扰论

本文将普遍声逆散射微扰论应用于弹性波层析成像问题,在Born变换下推出了以旋转角为补偿参数的各阶微扰重建公式,实现了对非均匀各向同性散射体内3个参数（质量密度ρ和两个Lamé系数λ,μ）的同时重建. 对于层析成像问题,在弹性波的传播过程中P波与SV波有耦合,但它们不会和SH波发生耦合,于是可以得到3个形式相对简单的标量方程. 在Born变换下,在散射波中引入微扰参数,将散射体的3个参数分别按该微扰参数展开,然后利用二维自由空间的Green函数分别得到散射的P波、SV波和SH波的积分表示. 最后,经一维傅氏变换后,得到Born变换下散射体3个参数的各阶微扰重建公式.     

伪谱法地震波正演模拟的多线程并行计算

地震波的正演模拟,尤其是3D正演模拟,往往涉及大规模的数据存储和计算,问题的规模往往超出计算机的物理内存,或者计算时间让问题的求解者难以忍受,即使采用目前存储和计算能力很强的计算机,其计算费用仍然是十分昂贵的.本文提出一种基于多线程协同使用多CPU和计算域分割的正演模拟并行计算技术,使得问题的求解过程得以加快,大大地缩短了用户等待的时间.为了检验我们的并行算法的可行性,文中以傅利叶正演模拟技术为例,给出了声波和3D各向异性弹性波模拟的例子,并对不同版本（串行、并行）运算效率进行了比较,证实了方法的有效性.     

应用混合变量最小势能原理求解有一个悬空角点弹性矩形薄板的稳定

x将混合变量的最小势能原理推广到求解弹性矩形薄板的稳定问题中,求解了有一个悬空角点弹性矩形薄板的稳定问题,并给出了相应问题确定临界载荷的特征方程及计算结果,为工程中薄板的设计计算提供了有效的参考,尤其是对现代建筑和桥梁中的受压构件的稳定设计和计算提供了一个简捷有效的计算方法。通过文中计算表明,混合变量的最小势能原理适用于有一个悬空角点矩形板稳定问题的求解,从而验证了它用于薄板稳定问题的一般性和通用性。     

非保守弹性动力系统两类变量的广义拟变分原理及其应用

按照广义力和广义位移之间的对应关系,将弹性动力学的各基本方程分别乘上相应的虚量,然后在相应的体积域和面积域上积分,将积分式代数相加,再将代数和在时间域上积分,代入本构关系,并考虑到体积力和面积力均为伴生力,进而建立了非保守弹性动力系统的第1类两类变量的广义拟变分原理;再应用类似的方法,通过代入另一类本构关系,建立了非保守弹性动力系统的第2类两类变量广义拟变分原理.应用第1类两类变量广义拟余能原理给出同时求解1个典型的非保守弹性动力系统的固有频率、变形和内力的计算方法.最后,讨论了有关问题.     

井筒声场研究回顾与展望

声波测井是测井学科的三大方法之一,在储层评价中起到了重要作用．随着陆上油气勘探目标趋向孔隙结构复杂、非均质性强的油气储集层,基于地层各向同性介质的测井理论显露弊端;因此研究复杂地质条件下的井筒声场分布特征、为储层解释和发展测井新方法提供理论依据成为必然．本文对井筒声场研究进行了回顾,并就当前井筒声场研究现状进行了总结,对其发展趋势和应用前景进行了分析。     

弹性波CT技术在大足石刻岩体破碎带探测中的应用

大足石刻是中国现存的重要石质文物之一,具有重要的历史和文化价值,其中有的石刻造像曾被列为国家石质文物保护“一号工程”;最近,文物保护部门又启动了对危岩体的加固和防渗工程。为确定大足石刻宝顶山大佛湾岩体中破碎带发育情况,应用跨孔弹性波CT技术对指定位置的工程岩体进行了探测。CT探测得到的岩体弹性波速成像清楚地揭示了两个钻孔之间岩体中破碎带的规模和产状,与岩体所在的区域地质构造应力场作用协调一致,证明了探测成像结果的可靠性,为大佛湾岩体的渗流和稳定性评估提供了依据。     

地震波衰减的物理机制研究

对近年来有关地震波衰减的物理机制进行评述和初步研究认为,对于地壳介质,造成地震波衰减的主要原因是地壳内存在大量裂隙,裂纹中饱含水或部分含水,地震波传播时引起裂隙中的流体运动,从而造成地震波的衰减,对于上地幔,扩散控制的位错阻尼机制可能是造成地幔地震波衰减的一个重要原因,而上地幔软流层部分熔融的存在也是造成地震波衰减不可忽视的原因．     

玛尼地震的InSAR形变研究与模拟

采用欧洲空间局ERS-2的星载干涉雷达数据,选取1997年11月8日MW7.6级玛尼地震作为研究对象,采用了差分干涉方法,在通过对覆盖同一地区的SAR数据进行差分干涉处理,得到了玛尼地震的视线向同震形变场。经研究发现:该地震形变场呈长轴近北东东向不规则椭圆形分布,地表破裂带长度约为130km,发震断层走向约为78°,断裂为左行走滑特征。断层以南为隆起区,在发震断层附近最大视线向隆起位移量为113.6cm,断层以北为沉降区,最大视线向沉降位移量为170.4cm。基于Okada模型实现了具有复杂结构的4段断层段参数的InSAR形变场数据模拟,获得断层的最大走滑为6m,估计出玛尼地震的标量地震矩M0为2.69×10^20Nm,计算得到的矩震级MW为7.6。证明了研究方法的正确性和研究结论的可靠性。