声波加pml边界条件(共10页).docx

雷克子波因计算机资源有限，不便在太大的区域做PML的计算，一般只在边界上取有限宽度的区域作为PML计算区域。根据相关文献的研究，PML边界区域最少长度应为半个波长6。本文综合考虑了效果与开销等因素，选取了边界上50层作为PML的计算区域。常规计算区域与PML边界区域的如图2-4所示。衰减系数式中为PML层的厚度，为层内的点距PML与非PML的边界的距离，为纵波速度，那么在PML边界区域内，对于式（2-13）即为理论反射系数，一般取0.001较为合适，为方向的空间步长。，可看作为在常规的计算方程基础上，减去一项进行PML的阻尼修正项。因本文中只考虑各项同性介质中的地震波传播规律，故可做假设。在此利用一下三个假设：因为以上三个近似精度均为时间方向上的近似，且时间精度均为二阶精度，因交错网格技术的时间精度为二阶，故以上近似不影响本式的计算精度。故可得：（2-18a）同理：（2-18b）