有限元法与伪谱法混合求解弹性波动方程   总被引：6，自引：0，他引：6  

马德堂  朱光明 《地球科学与环境学报》2004,26(1):61-64

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.  相似文献   

虚谱法求解波动方程的算法改进     

马德堂  朱光明 《地球科学与环境学报》2000,22(2):63-66

虚谱法利用付里叶变换和有限差分求解波动方程 ,是地震波场正演模拟的重要方法之一 ,但其常规算法运算量大 ,运算时间长。为此 ,我们在常规算法的基础上作了些改进 ,主要是提出了反付里叶变换的快速算法 ,并运用正反付里叶变换的快速算法和避免重复运算等手段 ,使运算量大大减少 ,可节省运算时间三倍多  相似文献   

双重网格井间地震层析成像技术     

马德堂  朱光明  张文波 《地球科学与环境学报》2005,27(4):83-86

提出了一种不同于常规的井间地震层析成像技术，即在较细的网格上进行射线追踪，以提高射线路径和旅行时的计算精度，在较粗的网格上进行层析反演成像，使网格像素上射线覆盖的最低次数达到一定要求，以提高成像质量。分析了采用这种双重网格地震层析成像技术的必要性，模型实验证实了此技术实用性和有效性。  相似文献   

井间地震初至波分析   总被引：1，自引：0，他引：1  

张文波  朱光明  马德堂  杜有成 《地球科学与环境学报》2006,28(1):70-74

在模拟对比不同模型条件下初至波的基础上,分析了造成初至不清晰两种可能的原因,一种是初至波为续至波所干扰,另一种初至波在垂直方向上的投影很小甚至为零。并指出可以综合利用水平分量和垂直分量记录来提高初至时间拾取的精度。  相似文献   

基于双线性插值的三维横向各向同性介质初至波射线追踪   总被引：1，自引：0，他引：1  

马德堂  杨春雨  李绪宣 《地球物理学进展》2014,(3)

给出了基于旅行时双线性插值(LTI)的三维横向各向同性(VTI)介质中的地震初至波旅行时计算及射线追踪方法,通过三维均匀及水平层状VTI介质模型的试算,证实了该方法的准确性及适应性.该方法可用于三维VTI介质中的深度偏移及层析成像.  相似文献   

一种改进的二阶弹性波动方程的最佳匹配层吸收边界条件   总被引：3，自引：0，他引：3  

刘有山  刘少林  张美根  马德堂 《地球物理学进展》2012,27(5)

在分析前人有关二阶波动方程的最佳匹配层(PML)吸收边界条件的构建方式的基础上,讨论了这些PML吸收边界条件在计算效率和计算精度方面的不足,并给出了一种新的改进PML吸收边界条件.本文的核心思想是在频率域引入一个中间变量,避免了褶积运算,对该变量做傅里叶逆变换得到时间域的精确解.通过模型试算,把本文的改进算法与前人的算法的计算精度和计算时间进行比较.对比分析表明,本文提出的PML吸收边界条件计算量小且精度较高,并且是一种稳定的算法.  相似文献