地震勘探原理教学中引入数值模拟方法探索.docVIP

地震勘探原理教学中引入数值模拟方法探索 　　[摘要] 研究中针对地震波传播理论教学中公式多而繁杂、内容抽象等特点，提出了将数值模拟方法引入到“地震勘探原理”课程教学中来，将地震勘探的理论教学与数值模拟技术紧密结合在一起。通过采用数值模拟，可将一些抽象的地震波动力学和运动学传播理论进行直观、形象、动态的展示，激发学生的学习兴趣，调动学生的求知欲望，从而收到事半功倍的教学效果。 　　[关键词] 地震勘探；课程教学；数值模拟；地震波；传播理论；教学效果 　　[中图分类号] G424 [文献标识码] A [文章编号] 1674-893X（2014）05?0054?03 　　地震勘探课程在应用地球物理学中是一门非常重要的课程。对于地震勘探课程的学习，一般可分为地震勘探原理和地震资料数据处理与解释两部分。地震勘探原理这门课程主要是学习地震勘探的基本理论、基本原理、野外数据采集工作方法等，为学习后续课程-地震资料处理与解释打下基础。地震勘探原理课程教学大纲要求学生必须掌握地震波动力学和运动学的基本概念、基本原理与分析方法，但是有关地震波动力学和运动学方面的内容公式较多，推导较繁，而且比较抽象，因此学生在学习的过程中普遍反映不容易理解，学习效果差。因此，如何提高地震勘探原理课程的学习效果是任课教师必须认真思考的问题[1]。本文针对“地震勘探原理”课程公式多、内容抽象等特点，结合我校应用地球物理的专业特色，探讨了在该课程课堂教学中使用数值模拟的方法来提高学生的学习兴趣，收到了较好的教学效果[2-4]。 　　一、地震波场数值模拟简介 　　地震波模拟是根据给定地下介质的结构模型和相应物理参数来模拟地震波的传播过程，从而研究地震波在地下介质中的传播规律。由于模拟过程中可直观、形象、动态地显示地震波动力学和运动学传播特征，非常容易调动学生的学习兴趣和求知欲望，可以收到事半功倍的教学效果。一般来说，地震波模拟可分为物理模拟和数值模拟两种方法。物理模拟是在实验室内将野外的地质构造和地质体按照一定的比例制作成物理模型，然后利用超声波或激光超声波等方法对野外地震勘探方法进行模拟；数值模拟就是利用有限差分、有限元等数值方法求解波动方程，从而获得已知模型的地震波传播。采用物理模拟方法存在费用高、选材困难等缺点，并且不适合于课堂理论教学；而数值模拟只需要一台较高计算速度的计算机就可全部解决问题，既简单方便，成本又低，且非常适合于课堂理论教学。 　　地震波场数值模拟方法主要有几何射线法和波动方程法。几何射线法也称为射线追踪法，其主要理论基础是，在高频近似条件下，地震波的主能量沿射线轨迹传播，即根据地震波的传播规律确定地震波在实际地层中传播的射线途径，并运用惠更斯原理和费马原理来重建射线路径，利用程函方程等计算射线的旅行时间。射线法的主要优点是概念明确、显示直观、运算方便、适应性强，缺点是旅行时的计算在一定程度上是近似的，特别是对复杂构造进行三维射线追踪时繁琐且误差较大。波动方程数值模拟方法是以弹性（粘弹性）理论及牛顿力学为基础的，求解双典型偏微分方程-波动方程为手段的一种数值模拟方法。这种方法不仅能保持地震波的运动学特征，而且还能保持地震波的动力学特征。根据地震波场数值模拟所采用的这两种不同方法的特点，我们在地震勘探原理课程讲授过程中采用波动方程弹性波数值模拟方法，这种方法的模拟结果可以使学生更好地了解地震波在传播过程中的运动学和动力学特征。目前用于波动方程数值模拟的方法主要有限差分法、有限元法、虚谱法等，在这些方法中有限差分法在计算精度、计算效率上都占有较大的优势，因此我们在数值模拟时采用了有限差分模拟方法[5]。 　　二、应用实例 　　地震波数值模拟在地震勘探原理课程教学中的应用比较广泛，在大部分的知识点讲解中都可采用数值模拟来对一些方法原理进行动态演示。在本文中主要通过两个知识点的应用来说明数值模拟方法在地震勘探原理课程教学中的效果。 　　（一）地震波的传播及界面处的反射、折射、透射和波形转换 　　在讲解地震波场的基本知识一节中，主要讲解地震波的传播特点，其中一些名词如波前、波后、球面波、平面波等，学生理解起来比较抽象，一些概念只能靠死记硬背。在讲解地震纵、横波的传播特点时，对于波在界面处发生的反射、透射、折射以及波形转换等，由于学生对波场的概念没有直观的认识，所以老师讲解起来也非常抽象和费力。在这种情况下，我们将数值模拟技术引入到课堂教学中，通过开发地震波场实时模拟软件，只需设置好模型和参数，软件将以动画的形式动态展示波场传播的全过程。 　　下面以两层介质为例进行说明：模型大小为600m×300m，网格大小为 ，时间采样间隔为 。模型上层纵波波速为3 000m/s，横波波速为1 732m/s，密度为1.8g/cm3；下层纵波波速为4