辅助微分方程完全匹配层在声波方程数值模拟中的应用     

赵建国  史瑞其  陈竞一  赵维俊  王宏斌  潘建国 《吉林大学学报(地球科学版)》2014,44(2):675-682

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。  相似文献   

《化工矿产地质》电子类参考文献著录格式(暂定)     

《化工矿产地质》2005,27(2):102-102

由于国家对电子类参考文献尚未颁布参考文献著录细则,为避免著录格式的混乱,本刊参照华东师范大学学报(自然科学版)等相关期刊关于电子类参考文献著录格式,特暂定出本刊电子类参考文献著录格式：  相似文献   

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法   总被引：2，自引：1，他引：1       下载免费PDF全文

薛昭  董良国  李晓波  刘玉柱 《地球物理学报》2014,57(4):1209-1223

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.  相似文献   

汶川地震断裂带科学深钻井孔附近微震观测问题及微震数据转换     

叶庆东  周晓峰  丁志峰  周乐 《地震地磁观测与研究》2014,35(5)

为了配合汶川地震科学深钻,弄清楚钻孔附近浅层的断层面结构,中国地震局地球物理研究所分别在四川省绵竹市天池乡和四川省绵阳市南坝镇布设15套南非矿山地震研究所(IMS)生产的矿山地震仪.分析矿山地震仪天然微震监测中仪器布设、数据采集和数据处理等方面遇到的问题,介绍与该仪器对应的数据的文件格式,并实现该数据与通用地震数据格式SAC二进制格式的转换.  相似文献   

GUIDE TO AUTHOR Sciences in Cold and Arid Regions has no page charges     

《寒旱区科学》2010,(3):F0003-F0003

  相似文献   

Transition from wavelets to ripples in a laboratory flume with a diverging channel   总被引：3，自引：1，他引：2  

William B. RAUEN  Roger A. FALCONER 《国际泥沙研究》2008,23(1):1-12

An experimental investigation on the initiation and development of bed forms on a bed of fine silica sand was conducted under alluvial flow conditions in a laboratory flume with a diverging channel. The main aims of the study were to assess： i） the steepness of bed forms in the transition stage of development; and ii） the threshold height of wavelets （ηt） that triggered the start of ripple development. Detailed bed profile measurements were carried out using an acoustic Doppler probe, traversed longitudinally over the sediment bed at various experimentation times. The bed form dimensions were extracted from such bed profile records and analysed for the wavelet, transition and equilibrium stages. It was found that the steepness of ripples in the transition and equilibrium stages were similar, confirming predictions of previous mathematical model simulations. A lognormal distribution fitted the wavelet length data. The wavelet threshold height was estimated as ηt ≈ 7 mm, or ηt≈ 80 in wall units. Such a height magnitude suggested that ripple development could be triggered by the wavelets reaching the outer flow zone of a turbulent boundary layer. The ηt value obtained corresponded generally to the intersection point between two predictive equations for bed form dimensions. A formulation was developed to predict ηt as a function of the sediment grain size, which was confirmed for the fine sand used in this study.  相似文献   

一维土层非线性地震反应分析的解析递推格式法     

甘杨  李凡  李大华 《吉林大学学报(地球科学版)》2006,36(4):631-635

引入“结构动力学数值分析解析递推格式法”求解软土地基地震反应分析问题。就一维地基非线性地震反应分析问题,将解析递推格式方法与等效线性化方法在各种场地土的计算结果进行了比较,结果表明：(1) 对软土层地震反应分析时,两种方法计算结果差别较大; （2）入射波幅值较强时,两种方法计算结果差别较大;（3）对于中等硬度土层、入射波幅值较弱时,计算结果相似, 并验证了该方法的实用有效性。  相似文献   

南海海底地震仪异常数据的分析和处理   总被引：4，自引：2，他引：2       下载免费PDF全文

王强  丘学林  赵明辉  黄海波  敖威 《地球物理学报》2016,59(3):1102-1112

海底地震仪(Ocean Bottom Seismometer,OBS)数据处理至关重要,是获取深部地壳结构的基础与前提.2006年实施OBS2006-2测线时,有2台OBS(OBS03,OBS06)数据出现异常,无法使用.由于海上航次花费巨大,采集到的数据弥足珍贵.本文采用数据格式检查、邻近台站对比分析、重采样等方法,成功地对这2台OBS数据进行了解编处理,得到了这两个台站的综合地震记录剖面;利用上述方法对2011年实施的OBS973-3测线中的异常台站OBS03进行了分析处理,同样得到了OBS03台站的综合地震剖面;通过查看两次海上实验班报发现,OBS2006-2测线之OBS06与OBS973-3测线之OBS03内部Sedis编号相同,为同一台记录仪器,再一次验证上述处理方法正确可行;然后对OBS2006-2测线2个台站进行震相识别与走时拾取后,利用前人纵波速度模型开展了射线追踪与走时模拟.此次对异常OBS数据的重新处理工作,不仅为OBS探测提供了宝贵的数据处理经验,而且将提高OBS2006-2测线地壳结构的可靠性和约束性,具有重要的研究意义.  相似文献   

适用于声波方程数值模拟的时间-空间域隐式有限差分算子优化方法          下载免费PDF全文

陈东  梁文全  辛维  杨长春 《地球物理学报》2016,59(4):1491-1497

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.  相似文献   

DGN数据向CASS 6.1格式转换的实现     

秦智慧 《测绘与空间地理信息》2009,32(5):130-131,134

数据格式之间的无损转换是不同平台之间实现数据共享的一个重要问题。本文结合DGN数据向CASS6.1的转换实例,分析了两种数据文件的结构特点,研究出两种数据格式转换的思路,并针对转换中的技术难点,提出了解决的方法。  相似文献