2025年大学《地球物理学》专业题库- 地震波传播模拟技术研究

2025年大学《地球物理学》专业题库——地震波传播模拟技术研究考试时间：______分钟总分：______分姓名：______一、选择题（每小题2分，共20分。下列选项中，只有一项是符合题目要求的。）1.在地震波传播模拟中，选择中心差分格式对一维波动方程进行离散，其主要目的是为了达到：A.简化计算过程B.便于处理边界条件C.获得二阶空间和时间上的精度D.减少对计算机内存的需求2.当地震波从速度较低的介质入射到速度较高的介质界面时，反射波和透射波的振幅通常表现为：A.反射波增强，透射波减弱B.反射波减弱，透射波增强C.反射波和透射波均增强D.反射波和透射波均减弱3.在使用有限差分法模拟二维地震波传播时，为了保证数值解的稳定性，时间步长Δt的选择必须满足特定的条件，该条件通常被称为：A.网格间距条件B.介质参数条件C.Courant-Friedrichs-Lewy（CFL）条件D.边界条件匹配条件4.PML（完美匹配层）边界条件在地震波传播模拟中的主要作用是：A.模拟自由表面B.模拟吸收边界，减少波的全反射和泄漏C.模拟层间反射D.模拟介质不均匀性5.对于各向异性介质，地震波的传播速度不仅取决于位置，还与波的传播方向有关，这种特性最适合采用哪种数值方法进行模拟？A.有限元法（FEM）B.谱元法（SEM）C.常用的各向同性介质有限差分法D.无需特殊数值方法，差分法同样适用6.在地震波传播模拟结果中，绕射现象通常出现在：A.介质边界处B.介质内部存在不连续体（如断层、盐丘）的尖顶或边缘处C.波源附近D.吸收边界附近7.地震走时曲线是地震学中重要的成果之一，它在地震波传播模拟中主要用于：A.判断介质是否存在断裂B.确定地下介质的速度结构C.计算震源机制解D.规划地震观测台站8.将有限差分法、有限元法等数值模拟方法得到的地震波形与实际观测到的地震记录进行对比分析，其主要目的是：A.验证模拟软件的功能B.评估模拟结果的精度和可靠性C.学习数值方法的基本原理D.研究震源特性9.地震波传播模拟技术在油气勘探中的主要应用不包括：A.模拟地震反射剖面B.计算油气藏的储层参数C.进行地震偏移成像D.预测地震地质灾害10.随着计算技术的发展，现代地震波传播模拟技术正朝着哪个方向发展？A.更简单的数值模型B.更低精度的离散化C.更高精度、更大规模（更高维度、更复杂模型）和更高效的计算D.仅限于理论研究，不再进行实际应用模拟二、填空题（每空2分，共20分。）1.地震波在均匀介质中传播时，其振幅随距离的增加而呈指数衰减，这种现象称为__________。2.有限差分法（FDTD）是一种基于将偏微分方程离散化为__________方程的数值方法。3.在地震波传播模拟中，为了模拟地表或自由界面，常采用__________边界条件。4.共中心点（CDP）道集是地震勘探中常用的数据组织方式，它对应于以震源和检波器连线为轴线的__________剖面。5.地震波在介质界面处发生的反射和折射现象，其强度由__________决定。6.偏移成像技术本质上是通过模拟地震波从像点出发的__________来重建地下结构。7.地震波传播模拟需要设置震源模型，常见的震源类型包括__________震源和点震源。8.为了提高数值模拟的精度，除了采用高阶差分格式外，还可以考虑使用__________等数值方法。9.地震波传播模拟结果的可视化通常包括波场图和__________的绘制。10.将机器学习算法与地震波传播模拟相结合，旨在提高模拟效率或实现__________的预测。三、简答题（每题5分，共25分。）1.简述一维波动方程的推导过程及其物理意义。2.解释什么是数值模拟中的稳定性问题，并简述有限差分法中CFL条件的作用。3.简述FDTD方法模拟二维水平层状介质中P波传播的基本步骤。4.说明吸收边界条件（如PML）在地震波传播模拟中是如何工作的，并解释其作用原理。5.简述地震波传播模拟技术在工程地震学中的一个具体应用实例。四、计算题（每题7分，共14分。）1.一维波动方程为∂²u/∂t²=v²*∂²u/∂x²，其中v为波速。试用中心差分格式（空间二阶，时间二阶）对上述方程进行离散，并推导出时间步长Δt与空间步长Δx之间的关系（CFL条件）。2.假设一束P波从速度为v₁的介质入射到速度为v₂的介质界面，入射角为θ₁。试用斯涅尔定律（Snell'sLaw）计算反射角θᵣ和透射角θₜ。假设v₁=3000m/s，v₂=1500m/s，入射角θ₁=30°，请计算反射角和透射角的具体数值（要求写出公式和计算过程）。五、论述题（每题11分，共22分。）1.论述有限差分法（FDTD）在模拟地震波传播中的主要优缺点，并与其他两种常见的数值方法（如FEM或SEM）进行比较。2.结合地震勘探的实际应用，论述地震波传播模拟结果的分析与解释在油气勘探中具有哪些重要意义。试卷答案一、选择题1.C2.B3.C4.B5.A6.B7.B8.B9.D10.C二、填空题1.衰减2.代数3.吸收4.共偏移距5.反射系数（或斯涅尔定律/振幅透射系数与反射系数的关系）6.反射路径7.水平8.有限元法/谱元法/高阶格式9.走时曲线10.勘探目标属性（如流体饱和度/孔隙度）三、简答题1.解析思路：从质点运动方程或牛顿第二定律出发，考虑小扰动假设，推导出连续介质中质点振动满足的偏微分方程。物理意义：描述了弹性介质中微小扰动（地震波）在介质中传播的规律，其中第一项表示波的传播，第二项表示介质对扰动的恢复力。2.解析思路：数值模拟稳定性是指当模拟时间步长或空间步长选择不当，导致数值解在迭代过程中出现震荡、发散，无法收敛到真实解的现象。CFL条件是有限差分法中保证数值解稳定性的一个必要条件，它限制了时间步长与空间步长及波速之间的比例关系，确保数值格式是稳定的。满足CFL条件是进行有效FDTD模拟的前提。3.解析思路：首先建立二维水平层状介质模型，定义速度结构；然后设定震源位置和类型，以及检波器排列；接着，在模拟区域内划分网格；然后，根据FDTD更新方程，逐时间步长迭代计算每个网格点上的P波位移（或压力）分量；最后，收集震源位置周围或特定测线上的计算结果，生成波场图。4.解析思路：吸收边界条件设计成对离开计算域的波具有吸收作用的边界，使其能量被“吸收”而不是反射回计算域内部。PML通过在计算域外层额外设置一层具有一定“虚”介电常数或磁导率的介质，使得波在PML层中因能量耗散而衰减。其作用原理是基于匹配阻抗的思想，使计算域边界的阻抗与PML层的阻抗逐渐过渡，从而最大限度地减少波在边界处的反射。5.解析思路：例如，在核工程中，利用地震波传播模拟技术研究核电站所在区域的地质结构、断层分布和地震活动性，评估地震风险，为选址、设计抗震结构提供依据。具体可模拟地震波在地下结构中的传播路径和强度衰减，预测可能发生的地面运动，分析其对核电站关键设施的影响。四、计算题1.解析思路：对一维波动方程进行中心差分离散，得到关于时间n+1和空间m+1/2（或m-1/2）时刻波场值u^n+1_m与前一时间n和空间m（或m+1/2，m-1/2）时刻波场值u^n_m及前一时间n+1和空间m（或m+1/2，m-1/2）时刻波场值u^n-1_m的关系式。然后，为了获得二阶精度，对空间二阶差分和時間二阶差分进行泰勒展开，并整理，最终得到更新公式。最后，将时间步长Δt用空间步长Δx和波速v表示，得到CFL条件。推导过程（简要）：∂²u/∂t²≈(u^(n+1)_m-2u^n_m+u^n-1_m)/(Δt)²∂²u/∂x²≈(u^n_m+1-2u^n_m+u^n_m-1)/(Δx)²代入波动方程，要求空间和时间二阶精度：(u^(n+1)_m-2u^n_m+u^n-1_m)/(Δt)²=v²*(u^n_m+1-2u^n_m+u^n_m-1)/(Δx)²整理得：u^(n+1)_m=2u^n_m-u^n-1_m+(vΔt/Δx)²*(u^n_m+1-2u^n_m+u^n_m-1)令r=vΔt/Δx，则：u^(n+1)_m=2u^n_m-u^n-1_m+r²(u^n_m+1-2u^n_m+u^n_m-1)CFL条件：为了保证稳定性，要求该差分格式是无条件稳定或有条件稳定。对于此二阶中心差分格式，是条件稳定的，要求r≤1，即vΔt≤Δx。这就是CFL条件。2.解析思路：应用斯涅尔定律，即入射角θ₁的正弦值与反射角θᵣ的正弦值相等，且等于透射角θₜ的正弦值与两个介质波速之比，即sinθ₁=sinθᵣ和sinθ₁=(v₂/v₁)sinθₜ。根据这两个关系式，可以解出反射角θᵣ和透射角θₜ。将已知数值代入公式进行计算即可。计算过程：sinθᵣ=sinθ₁=sin30°=0.5θᵣ=arcsin(0.5)=30°sinθₜ=sinθ₁*(v₁/v₂)=sin30°*(3000/1500)=0.5*2=1θₜ=arcsin(1)=90°反射角θᵣ=30°，透射角θₜ=90°。五、论述题1.解析思路：首先分别阐述FDTD方法的优点和缺点。优点可包括：原理简单直观，物理图像清晰，易于编程实现（尤其一维和二维），能够自然地处理复杂的边界条件和介质不连续性，是发展最快和应用最广的方法之一。缺点可包括：对于复杂三维模型计算量巨大，计算效率相对较低（尤其未优化时），在处理强吸收或低波速介质时可能存在数值反射，精度受网格尺寸限制等。然后，与其他方法比较，FEM优点是适应复杂几何形状能力强，能提供位移和应力场的分布，缺点是推导复杂，编程实现相对困难。SEM优点是精度高，尤其适用于波动方程，效率高，缺点是推导和编程比FDTD更复杂，主要用于特定问题（如波动问题）。总结时强调FDTD的普及性和优势，但也需认识到其局限性以及FEM、SEM等在特定场景下的价值。2.解析思路：首先说明地震波传播模拟是获取地下结构信息的重要手段，模拟结果可以提供比理论计算更丰富的地下信息。然后，从油气勘探角度论述其意义：*储层识别与评价：通过模拟波在含油气储层和围岩中的传播差异（如振幅、频率变化），帮助识别潜在的储层位置和评价其物性。*构造解释：模拟断层、盐丘等复杂地质构造对地震波的反射、折射和绕射影响，辅助解释实际地震剖面，确定构造形态和幅度。*岩性预测：结合岩石物理模型，模拟不同岩性对波传播的影响，预测地下岩性分布。