2025年大学《地球物理学》专业题库- 构造地震波速成像技术的研究与应用

2025年大学《地球物理学》专业题库——构造地震波速成像技术的研究与应用考试时间：______分钟总分：______分姓名：______一、名词解释（每题3分，共15分）1.层速度2.地震层析成像3.静校正4.共中心点道集5.波速成像分辨率二、简答题（每题5分，共25分）1.简述影响地震波速度的主要地质因素。2.地震反射法中进行动校正的主要目的是什么？3.与地震反射法相比，地震层析成像方法在探测路径结构方面有何优势？4.简述走时层析成像中，数据依赖性（走时残差）的主要来源。5.简述利用地震波速成像技术进行城市地下结构探测的主要流程。三、计算题（每题10分，共20分）1.在一维水平层状介质中，测得某测线上一炮点与两个检波点之间的直达P波走时分别为t_A=1.0s,t_B=1.2s,t_C=1.5s（假设检波点B位于炮点与检波点C的中点）。已知P波在地表介质中的速度v_surface=1500m/s。若第一层厚度h_1=1000m，第二层顶界面深度为z_2。请利用走时数据，计算第二层底界面深度z_2（假设各层介质均匀、水平、连续，且P波在各层中的速度分别为v1=2500m/s,v2=3500m/s）。2.假设某区域进行二维地震层析成像，采用有限差分法建立正演模型。已知震源位于坐标(0,0)，检波点位于坐标(x_d,0)，地下存在一个速度异常体（速度异常Δv）。请简述利用有限差分法计算检波点处接收到的理论走时（或走时残差）的基本步骤。说明需要设置的参数和需要执行的运算。四、论述题（每题12.5分，共25分）1.论述地震层析成像反演中非线性问题的特点，并比较两种主要的反演算法（如梯度下降法与迭代法）在解决此类问题时的基本思想与区别。2.试分析地震波速成像技术在工程场地勘察中的应用潜力与面临的挑战，并提出一些建议。试卷答案一、名词解释1.层速度：指地震波在垂直于等时面方向上传播的速度，是介质密度和体积波（P波和S波）速度的乘积的平方根，与岩层的孔隙度、流体性质等因素密切相关，是地震资料解释中常用的参数。**解析思路：*定义层速度，并点明其物理意义（垂直传播速度）和与介质性质（密度、波速）的关系，强调其在解释中的重要性。2.地震层析成像：利用地震波（如P波、S波）的传播时间（或其他物理量）数据，通过建立观测数据与地下介质参数（通常是波速）之间的关系，建立数学逆问题，以确定地下介质参数分布的一种成像方法，类似于医学CT成像。**解析思路：*定义地震层析成像，说明其利用的物理量（走时等）、解决的反问题类型（参数分布），并类比医学成像帮助理解。3.静校正：指消除或减弱由近地表不均匀性（如地形起伏、表层速度变化、盐丘、断层等）引起的地震波旅行时差异，使来自同一地质界面的反射波到达同一深度（或时间）位置的处理步骤。**解析思路：*定义静校正，说明其目的（消除近地表影响），指明主要校正的对象（地形、速度变化等）和最终效果（统一反射波位）。4.共中心点道集：在地震反射剖面数据处理中，将同一炮点接收到的、来自不同检波点且位于通过炮点和检波点连线的垂直平面上的所有道集中在一起形成的数据集。**解析思路：*定义共中心点道集，明确其形成方式（同一炮点、不同检波点、通过炮检连线的垂直平面），这是动校正前的关键道集。5.波速成像分辨率：指地震波速成像方法能够区分的最小地下结构尺寸或两个地质体开始分离的最小速度差异。分辨率受到多种因素影响，如信号噪声水平、波长、成像几何、反演算法等。**解析思路：*定义分辨率，说明其衡量的是区分能力（最小尺寸或速度差），并指出影响分辨率的因素。二、简答题1.简述影响地震波速度的主要地质因素。**解析思路：*列举主要因素并简述其影响：岩石类型（岩浆岩、沉积岩、变质岩）、孔隙度与饱和度（孔隙度增加，速度通常降低）、流体性质（密度、压缩性，如盐水波速大于淡水）、温度（温度升高，速度一般降低）、压力（压力增大，速度一般升高）、围压、岩石的胶结程度与孔隙填充物（胶结好、填充物致密，速度高）。2.地震反射法中进行动校正的主要目的是什么？**解析思路：*说明动校正的目的：消除由炮点与检波点之间距离（偏移距）变化引起的、与地下界面倾角相关的旅行时差异，使得反射波同相轴在地形校正后的道集中趋于水平，从而能够将共中心点道集中的同相轴叠加成地震剖面。3.与地震反射法相比，地震层析成像方法在探测路径结构方面有何优势？**解析思路：*指出层析成像的优势：不依赖明确的反射界面，能够探测到连续介质内部的速度异常体，即使在没有明显反射波或反射波弱的区域也能成像；成像信息与射线路径上的介质性质直接相关。4.简述走时层析成像中，数据依赖性（走时残差）的主要来源。**解析思路：*列举数据依赖性的主要来源：射线路径弯曲（非直射路径，如过弯道、绕射）、近地表效应（如地形起伏引起的走时修正未完全补偿）、速度模型不精确（背景模型与实际模型差异）、噪声与随机误差（观测走时存在测量不确定性）。5.简述利用地震波速成像技术进行城市地下结构探测的主要流程。**解析思路：*描述主要流程：场地勘察与数据采集方案设计（选择合适方法如微震、面波或折射/反射，确定测线布置、仪器参数）、现场数据采集（确保数据质量）、数据预处理（去噪、滤波、坐标转换、基线校正等）、速度分析（获取浅层速度结构）、波速成像正演与反演（构建初始模型，进行迭代求解得到地下速度结构图像）、成像结果解释（结合地质资料解释地下结构特征如断层、空洞、覆盖层厚度等）。三、计算题1.在一维水平层状介质中，测得某测线上一炮点与两个检波点之间的直达P波走时分别为t_A=1.0s,t_B=1.2s,t_C=1.5s（假设检波点B位于炮点与检波点C的中点）。已知P波在地表介质中的速度v_surface=1500m/s。若第一层厚度h_1=1000m，第二层顶界面深度为z_2。请利用走时数据，计算第二层底界面深度z_2（假设各层介质均匀、水平、连续，且P波在各层中的速度分别为v1=2500m/s,v2=3500m/s）。**解析思路：**利用t_A和v_surface计算炮点处地表到检波点A的表层速度：x_A=v_surface*t_A。*利用t_B和v_surface计算炮点处地表到检波点B的表层速度：x_B=v_surface*t_B。由于B位于A和C的中点，x_B=x_A+v_surface*(t_C-t_B)/2。验证此关系：x_B=1500*1.2=1800m;x_A=1500*1.0=1500m;v_surface*(t_C-t_B)/2=1500*(1.5-1.2)/2=150m。x_A+150=1650m。此处原题t_B=1.2s,t_C=1.5s对应x_B=1800m，x_C=2700m，与B为中点矛盾。修正：假设题目意图是B点走时1.2s，C点走时1.5s，且B为AC中点。则x_C=v_surface*t_C=1500*1.5=2250m。x_B=x_A+(x_C-x_A)/2=x_A+(2250-x_A)/2。联立x_B=1500*1.2=1800m，解得x_A=1350m。*利用t_A和v_surface计算第一层底界面的走时：t_1=h_1/v1=1000m/2500m/s=0.4s。则第一层顶界面到A点的走时为t_A'=t_A-t_1=1.0s-0.4s=0.6s。A点到第一层顶界面的距离为v1*t_A'=2500m/s*0.6s=1500m。因此，第一层顶界面深度为1500m。*利用x_B和第一层顶界面深度计算第二层顶界面深度：z_2=第一层顶界面深度+x_B=1500m+1800m=3300m。*（可选检查）利用t_C和v_surface计算第二层顶界面到C点的走时：t_2'=x_C/v_surface=2250m/1500m/s=1.5s。第二层顶界面到第二层底界面的走时为t_2''=t_C-t_2'=1.5s-1.5s=0s。这表明C点位于第二层底界面处。*计算第二层厚度：h_2=v2*t_2''=3500m/s*0s=0m。这表明根据给定数据，C点恰好在第二层底界面，即z_2=h_1+h_2=1000m+0m=1000m。这与z_2=3300m矛盾。重新审视题目条件：假设题目意图是使用t_A,t_B,t_C三个走时，且B为AC中点。x_A=1500*1.0=1500m。x_B=1500*1.2=1800m。x_C=1500*1.5=2250m。B为AC中点，x_B=(x_A+x_C)/2，验证x_A=1500m,x_C=2250m,(1500+2250)/2=1875m。与x_B=1800m矛盾。结论：原题数据条件存在矛盾，无法得到唯一解。若假设t_B=1.2s,t_C=1.5s,B为AC中点，则x_A=1500m,x_B=1800m,x_C=2250m。需重新设定题目条件或采用其他已知条件（如直接给出x_B或x_C）。2.假设某区域进行二维地震层析成像，采用有限差分法建立正演模型。已知震源位于坐标(0,0)，检波点位于坐标(x_d,0)，地下存在一个速度异常体（速度异常Δv）。请简述利用有限差分法计算检波点处接收到的理论走时（或走时残差）的基本步骤。说明需要设置的参数和需要执行的运算。**解析思路：**需要设置的参数：地下介质的速度模型（包括正常速度v_0和异常速度v_异常=v_0+Δv的分布），震源位置(x_s,y_s)，检波点位置(x_d,0)，网格步长(dx,dy)，时间步长(dt)，计算总时间T_max，震源子波函数（如Ricker波）。*基本步骤：1.建立速度模型网格：将二维区域离散化为网格点，根据速度模型赋值每个网格点上的速度值（v_0或v_异常）。2.初始化震源和检波点：在震源位置(x_s,y_s)放置震源子波，在检波点(x_d,0)放置检波器（初始时无信号）。3.时间步进迭代：a.从t=0开始，以时间步长dt逐步增加时间。b.在每个时间步dt，根据当前时刻t的速度模型，计算波从震源传播到每个网格点所需的距离，并更新波到达该网格点的时间（例如，利用有限差分近似导数求解波动方程的初值问题）。c.根据波动方程的解，判断在时刻t是否波前到达检波点(x_d,0)。如果是，则将检波点在该时刻的波场值累加到检波器的总接收信号中。4.重复步骤3，直到达到总计算时间T_max。5.获取理论走时/信号：获取检波点处的总接收信号。可以通过分析信号到达的时间点得到理论走时（例如，信号能量峰值或零交叉点对应的时间）。走时残差则是计算得到的理论走时与观测走时之间的差值。*需要执行的运算：在每个时间步，需要进行大量的向量/矩阵运算（如有限差分格式的求解），包括对速度、位置、波场值进行更新，判断波到达等。具体运算依赖于所采用的有限差分格式（如显式、隐式、Leapfrog等）。四、论述题1.论述地震层析成像反演中非线性问题的特点，并比较两种主要的反演算法（如梯度下降法与迭代法）在解决此类问题时的基本思想与区别。**解析思路：**非线性特点：*模型方程的非线性：通常地震层析成像的正演问题（计算走时/波场）是基于非线性波动方程（或其简化形式）建立的，导致模型方程本身是非线性的。*观测数据与模型参数的强非线性关系：地下介质的速度结构变化对观测到的走时/振幅等数据的影响是非线性的，尤其是在存在较大速度异常或复杂路径时。*不适定性：层析成像反演问题通常是病态的（不适定的），即观测数据对模型参数的微小变化非常敏感，导致解对噪声和数据误差极其敏感，存在多解性。*计算复杂性：求解非线性逆问题需要大量的计算资源，尤其是当模型维度较高或采用复杂的正演方法时。*算法比较：*梯度下降法（及其变种）：**基本思想：*基于最速下降原理，通过计算目标函数（如走时残差的平方和）关于模型参数的梯度，迭代地朝着减小目标函数值的方向更新模型参数。**特点：*简单直观，计算量相对较小（每次迭代只需计算梯度），但收敛速度可能较慢，尤其是在接近最优解时可能陷入局部最小值。对初始模型的选择比较敏感。变种如共轭梯度法可以加速收敛。*迭代法（如高斯-牛顿法、Levenberg-Marquardt法、SIRT）：**基本思想：*在每次迭代中，利用当前模型和观测数据的近似关系，构建一个线性化的近似方程组，求解该方程组得到模型参数的更新量。高斯-牛顿法直接使用二阶导数（Hessian矩阵），Levenberg-Marquardt法在高斯-牛顿法基础上加入阻尼项，SIRT（同步/异步迭代重加权）则采用加权最小二乘的形式进行迭代。**特点：*通常收敛速度比梯度下降法快，尤其是在非线性程度较高的问题中。高斯-牛顿法收敛快但可能不稳定。Levenberg-Marquardt法鲁棒性较好，能同时处理非线性问题和重加权问题。SIRT方法概念简单，适用于大型稀疏系统。但迭代法通常需要计算导数（或其近似），计算量可能更大。两种迭代法都可能陷入局部最小值。**区别：*梯度下降法沿负梯度方向搜索，思想简单但可能慢；迭代法通过构建近似线性方程组进行搜索，通常收敛更快但实现更复杂，且对初始值和参数设置更敏感。2.试分析地震波速成像技术在工程场地勘察中的应用潜力与面临的挑战，并提出一些建议。**解析思路：**应用潜力：*探测隐伏断裂构造：精细刻画基岩断裂的位置、延伸方向、活动性，评估其对场地稳定性的影响。*获取浅层地速结构：确定地表以下几十米甚至上百米范围内的平均速度剖面，是进行地震小区划、场地响应分析、工程桩基设计、浅层地热勘探等的基础。*识别不良地质体：如隐伏的溶洞、空洞、软弱夹层、古河道、人工填土等，这些体往往具有不同于周围介质的速度特征，可用于场地地质灾害评估和地基处理方案设计。*提供岩土参数：获取的层速度剖面可以结合密度资料估算剪切波速度、泊松比、场地土层厚度等岩土工程参数。*三维成像：结合多道观测和三维数据处理，可以构建地下三维速度结构模型，提供更全面的场地信息。*面临的挑战：*近地表复杂效应：工程场地近地表通常结构复杂，速度横向变化快，存在地形起伏、覆盖层不均、高阶反射波干扰等问题，给数据采集和处理带来困难。*信号质量：由于场地环境通常较为