简介面向地震勘探学习者、处理解释人员与高校相关专业师生的PPT围绕Tesseral-2D模型正演软件的功能及特点展开重点说明模型正演在缓解复杂介质资料多解性、辅助观测系统与采集参数设计、验证处理流程和解释结果等方面的作用。资源包为1个PPT文件约2MB单文件结构便于集中浏览和课堂交流。目前已有226人学习。内容按软件结构与能力组织介绍Model builder、Computation Engine、Viewer和Processing Block四个功能包涉及LAS测井曲线建模、有限差分合成记录、波场快照、速度分析与CMP叠加、叠前叠后时间/深度偏移等环节同时讨论梯度、声学、弹性及各向异性介质能量衰减、初至时间场和复杂介质AVO效应并结合薄层、溶洞、砂泥互层等模型说明地震响应特征便于读者建立从模型设计到结果查看的完整认识。1. 多解性从哪来正演为什么是解释结论的后悔药一份解释成果能对得上井不代表它是对的。同一套反射同相轴可能同时是砂泥薄互层、含层间高速夹层的溶洞体、带各向异性的页岩层序三者共同产生的它们的走时差落在观测误差以内振幅谱也相近。这就是地震勘探绕不开的多解性地下数据只有一套能拟合它的地质模型有无数个。2D 模型正演软件存在的意义不是去造地震记录而是把一条解释结论丢进波动方程里跑一遍看它在物理上究竟能不能产生你观测到的那套波场。它是解释结论的后悔药也是把玄学式的经验判断换成可证伪实验的那道工序。搞懂它的底层机制本质上是搞懂地震记录里每一个相位、每一段振幅到底是被介质的哪一部分加工出来的。2. Model builder一条 LAS 曲线怎么变成可计算的 2D 密度-速度网格正演里最容易被低估的一步是建模。有限差分算得快不快、结果像不像八成在模型建好的那一刻就已经定了。Model builder 做的事情说起来朴素把一段二维地质剖面的密度和速度填进规则网格让后面的波动方程有东西可算。但填这个动作里藏着三组决策——填什么量、填多细、界面怎么过渡每一个都会在最后的记录上留下指纹而且往往是事后最难反推的指纹。2.1 模型的最小数据单元为什么 Vp、Vs、ρ 缺一不可网格里每个节点至少有三个数纵波速度 Vp、横波速度 Vs、密度 ρ。声学近似下只需要 Vp 和 ρ因为流体不传播剪切波一旦切到弹性介质Vs 就必须存在它决定了横波和转换波的速度也决定了泊松比进而决定了 AVO 的梯度项。密度则直接进波阻抗界面的法向反射系数是 R (ρ₂V₂ − ρ₁V₁)/(ρ₂V₂ ρ₁V₁)。三个数缺一个后面的计算要么降级成声学要么用经验关系凑——而凑出来的 Vs往往是 AVO 分析里最大的误差源。各向异性介质还要再加三个 Thomsen 参数ε 描述水平与垂直 P 波速度的相对差异δ 控制近垂直方向的波前曲率γ 描述 S 波分裂程度。VTI 介质里 P 波相速度可以写成 v(θ) Vp₀(1 δ·sin²θcos²θ ε·sin⁴θ)θ 是与对称轴的夹角。这个式子在正演内部并不是每步解析求值而是被离散进刚度张量的 c₁₁、c₁₃、c₃₃、c₄₄、c₆₆ 五个分量通过本构关系把应变换成应力。也就是说建模界面上填的 ε、δ、γ最后都变成了矩阵里的数字填错一个波前形状就歪而且歪得很隐蔽——近偏移距几乎看不出来远偏移距才暴露。2.2 LAS 曲线的重采样与横向外推从一维井到二维剖面LAS 格式的测井曲线是正演最扎实的输入因为它来自真实测量比任何经验公式都可靠。声波时差曲线给出 Vp密度曲线给出 ρ横波曲线有则更好。这里第一个坑就是单位声波时差常用 μs/ft而模型网格用 m/s两者之间差一个 0.3048 的换算因子。一旦搞错速度整体偏 3.28 倍合成记录的走时会离谱到无法解释。第二个坑是采样率测井采样常在 0.125 m 量级而模型网格的水平与垂直步长通常是 210 m直接按最近点抽稀会把薄层抹掉正确做法是先做抗混叠低通再抽点。第三件事是横向外推。一口井只给出一条一维曲线2D 剖面要靠地震解释层位作为骨架把井上的速度按层序外推到整条剖面层内插值、断层处断开、尖灭处收敛。这一步的精度直接决定了构造形态能不能被正确模拟。下面是一段把 LAS 转成模型网格的示意代码import numpy as np # 1) 声波时差转纵波速度AC 常见单位为 us/ft乘 0.3048 换算到 m def ac_to_vp(ac_us_per_ft): return 0.3048 * 1e6 / ac_us_per_ft # 结果单位 m/s # 2) 缺少横波曲线时的经验估计碎屑岩泥岩线近似仅作初值 def vp_to_vs_initial(vp_mps): return (vp_mps - 1360.0) / 1.16 # 单位 m/s非碎屑岩需替换关系式 # 3) 深度域测井重采样到模型网格先低通抗混叠再线性插值 def resample_to_grid(depth, value, grid_z, dz, filter_len9): kernel np.hanning(filter_len) kernel / kernel.sum() smoothed np.convolve(value, kernel, modesame) # 抗混叠防止薄层被抽稀抹平 return np.interp(grid_z, depth, smoothed) # grid_z 为等间距深度序列 # 4) 组装二维网格每一列来自一口井或一个外推控制点 def build_grid(grid_z, well_curves, lateral_positions): nx, nz len(lateral_positions), len(grid_z) vp np.zeros((nz, nx)); vs np.zeros((nz, nx)); rho np.zeros((nz, nx)) for j, curve in enumerate(well_curves): vp[:, j] resample_to_grid(curve[depth], ac_to_vp(curve[ac]), grid_z, 5.0) vs[:, j] resample_to_grid(curve[depth], curve[vs], grid_z, 5.0) rho[:, j] resample_to_grid(curve[depth], curve[den], grid_z, 5.0) return np.array([vp, vs, rho])函数里最关键的是filter_len与dz的匹配滤波器长度要覆盖至少两个网格步长否则低通等于没做。vp_to_vs_initial只适用于以石英为主的碎屑岩序列遇到碳酸盐岩或含气层段必须换成对应关系否则泊松比会算得离谱。build_grid返回的是三维数组切片索引顺序是深度优先——这一点在导出成计算引擎输入时要特别小心行列颠倒的模型跑出来的剖面看起来哪儿都对就是不对。2.3 介质类型的开关声学、弹性、各向异性在数值上差在哪里这套工具提供了四种介质模型梯度介质、声学介质、弹性介质、弹性各向异性介质。梯度介质指速度随深度连续变化线性或指数都可以它带来的直观效果是射线弯曲——波前不再是同心圆弧而是被连续折射成弧线这对模拟大套泥岩的连续变速很有用。声学介质是成本最低的一档只解压力场适合验证构造形态和首波路径。弹性介质要同时解两个速度分量和三个应力分量计算量和内存大约翻三到四倍。各向异性再往上叠加三个参数体积和耗时会继续涨。这四档不是精度越高越好的线性关系而是匹配问题的选择。要验证观测系统覆盖次数、看绕射波在偏移后是否收敛声学就够了要做 AVO 或者多分量极化分析必须上弹性要解释为什么远偏移距的时距曲线系统性偏离双曲线、为什么 NMO 校正后仍有剩余时差那就只能上各向异性。用错档位不是算不出来而是算出一个看起来很合理的错误答案这种错误最难排查。2.4 平滑还是台阶速度界面在网格上的处理真实界面是间断的规则网格是离散的两者天然冲突。界面正好落在两个节点之间时反射系数会被摊到一到两个网格上等效于给界面做了一次空间低通结果是反射同相轴后面拖出一串高频尾巴。常见做法是在界面处做一到两个网格的平滑过渡或者用等效介质把薄互层合并成一个带各向异性参数的整体块。后一种做法是权宜之计值得多说一句。把薄层组等效成 VTI 介质能在近偏移距把时距曲线对得很准但薄层调谐效应被藏进了 ε、δ 里到了远偏移距和频率域就会失真。这是设计上的取舍而不是缺陷——当你只需要验证构造和采集参数时等效模型的计算量比逐层建模低一个量级值这个价钱当你要研究薄层振幅随偏移距的变化时就必须老老实实把每一层画出来。心里清楚自己在用哪一种比用哪一种更重要。3. Computation Engine有限差分是怎么把波一步步推过网格的模型建好之后计算引擎要做的事可以概括成一句话在每一个时间步用当前时刻的应力去更新速度再用更新后的速度去更新应力如此往复。这句话听起来简单但它背后是一整套关于离散化、稳定性、边界处理和物理近似的工程决策。理解这些决策才能解释为什么同一个模型换个参数结果会从完美直接翻车到完全不能用。3.1 交错网格与 CFL 条件为什么 dt 一错整个波场就炸弹性波一阶速度-应力方程组在空间上通常用交错网格离散速度定义在网格面上应力定义在网格中心两者在空间上错开半个网格、在时间上错开半个步长。这样做的好处是差分模板更紧凑、界面处不连续量的处理更自然代价是实现上索引容易搞混。下面是二维情况的一步更新示意# 二维弹性交错网格单步更新示意差分算子 D 内部展开为 2N 阶模板 def step(vx, vz, sxx, szz, sxz, dt, rho, c11, c13, c33, c44): # 先由速度求应变率再更新应力时间上应力与速度错开半步 dvx_dx D(vx, axis1); dvz_dz D(vz, axis0) dvx_dz D(vx, axis0); dvz_dx D(vz, axis1) sxx dt * (c11 * dvx_dx c13 * dvz_dz) szz dt * (c13 * dvx_dx c33 * dvz_dz) sxz dt * c44 * (dvx_dz dvz_dx) # 速度用刚更新完的应力求梯度保证能量守恒性质不被破坏 vx dt / rho * (D(sxx, axis1) D(sxz, axis0)) vz dt / rho * (D(sxz, axis1) D(szz, axis0)) return vx, vz, sxx, szz, sxzc11到c44就是建模阶段填进去的刚度张量分量各向同性时 c11 c33 ρVp²、c44 ρVs²各向异性时按 Thomsen 参数展开。rho是密度数组逐点不同。D是沿指定方向的高阶差分算子阶数越高同样网格精度越好但边界处需要的辅助节点也越多。这段代码里最容易被忽视的是更新顺序先应力后速度、且速度用的是新应力这个顺序保证了离散系统在长时间积分下不发散。顺序写反短时间看不出问题跑上几千步就会缓慢长出一个虚假的直流偏移。时间步长 dt 不能随便给它受 CFL 条件约束。对二阶时间、2N 阶空间差分稳定性要求大致是 dt ≤ h / (Vmax · √(Σ|aᵢ|))其中 aᵢ 是差分系数h 是最小网格步长Vmax 是全模型的最大速度。这里有个血泪经验Vmax 必须取整条剖面里的最大值包括你建模时随手在底部加的那个高速基底。只按地表速度算 dt前几千步风平浪静波一碰到高速层立刻数值爆炸屏幕上全是 NaN。3.2 P 波与 SV 波的分离弹性张量里的独立分量2D 情况下P 波和 SV 波在平面内耦合传播SH 波是解耦的单独一套方程。这也是为什么 2D 正演特别适合研究极化地震勘探SV 分量和 P 分量的耦合关系被完整保留在 c13 这一项里它直接控制着转换波的强度。当界面两侧的 Vs 对比强烈比如砂岩进入含气层Vs 骤降c13 的变化会让一部分 P 波能量在界面上转成 SV转换波出现在时间剖面上就是那组不该出现的同相轴。能精确算出这个转换是弹性模型相对声学模型最大的价值。声学模型没有 Vs也就没有转换波用它做多分量采集设计会漏掉一整个波型。反过来说如果只是要看构造形态、验证 CMP 覆盖声学模型跑得快、迭代次数多反而更有效率。这种按问题选物理近似的思路在这套工具里贯穿始终。3.3 各向异性与吸收Thomsen 参数和 Q 值怎么进方程各向异性进入方程的方式前文说过是通过刚度张量。它在结果上的表现有三种一是远偏移距的时距曲线不再是双曲线走时差随偏移距呈更高阶增长二是同一反射层的 AVO 曲线在小角度范围内就被抬升或压低三是速度分析的 NMO 速度会明显小于水平层速度两者之差恰好与 δ 成正比。如果解释人员不知道模型是各向异性的很容易把它当成处理参数没调好反复折腾速度谱这属于典型的方向性错误。吸收衰减进入方程的方式则更麻烦一些。弹性模型是完全弹性的波传播不损失能量要模拟吸收要么引入粘弹性机制在刚度里加复数模量或者引入记忆变量要么用 Q 值构造松弛机制。原料里提到考虑波的能量衰减在这个语境下它影响的不只是振幅随时间的衰减还有相位——不同频率衰减不同子波在传播过程中会被重塑主频降低、相位延迟。做 AVO 时这一点尤其关键因为 AVO 的梯度项对振幅随偏移距的变化非常敏感吸收造成的相位漂移会让梯度估计算出系统性偏差。3.4 首波、初至时间场与波场快照Viewer 里的四个时刻在说什么波场快照是这套工具最有辨识度的输出。以声波模型为例快照分别在 t₁ 0.350 s、t₂ 0.510 s、t₃ 0.810 s、t₄ 1.110 s、t₅ 1.420 s 四个时刻实际是五个冻结波场可以直接看到波前在介质中的位置和形状。这不是为了好看而是为了把时间剖面上一团叠在一起的相位还原成地下正在发生什么。在这类快照里能看清三样东西。第一是直达波与首波的分离当存在高速层时波以临界角入射后沿界面滑行一边滑行一边向上辐射能量形成折射波首波。快照上能直接看到首波的波前是一条与界面平行、斜向上抬起的直线段它在时间剖面上就是那个初至之后的稳定同相轴。利用首波的斜率差反推高速层速度是声波模型正演最经典的用法之一。第二是绕射波溶洞、断点这类尺度小于波长的异常体不产生规则反射而是向各个方向辐射能量快照上表现为一个点源向四周扩散的圆弧这正是溶洞模型正演里最需要观察的特征。第三是薄层内的多次反射砂泥互层模型在快照上能看到能量在薄层之间来回弹最终在时间剖面上表现为多相位反射层——一组长得不像一波的复合相位。这套工具还有一项输出叫初至到达时间场即对每个网格点记录波第一次到达的时刻形成一个 t(x, z) 场。它等价于把波场快照里的波前位置全部记录下来对做层析反演、检验偏移速度模型特别有用——因为它是逐点的、可比较的而剖面只能定性地看。4. Processing Block合成炮集为什么还要走一遍速度分析和偏移算出合成炮集就完事了吗不是。正演真正的价值是它提供了一个答案已知的数据集可以拿去做处理流程的对照实验。你在地面上永远不知道真实的地下长什么样所以你无法判断一次速度分析做砸了到底是参数问题还是资料问题但在合成数据上这个判断是可以做的因为模型就是你手里的标准答案。Processing Block 提供的就是这层能力速度分析、CMP 抽道集与叠加、叠前叠后时间偏移、叠前叠后深度偏移全部可以在这套合成数据上跑一遍。4.1 CMP 道集、速度谱与叠加把正演结果当成野外资料对待合成炮集要按野外数据的流程走先定义观测系统炮点位置、检波点位置、偏移距范围、道间距、最大最小偏移距再抽 CMP 道集然后做速度分析。速度分析通常是相似系数谱对每个 CMP 位置在一系列试探速度上做 NMO 校正然后计算道集叠加后的相似系数峰值位置就是最佳叠加速度。合成数据的好处是峰值应该出现在哪里你是知道的——层速度可以直接从模型里读出来。这一步最容易暴露的问题恰恰是峰值不在预期位置。如果模型里存在明显的横向变速常规的单一速度谱会给出一个被平均过的折中值叠加剖面上表现为同相轴中途抖动甚至断掉。把合成数据的叠加结果和真实数据的叠加结果放在一起对比如果两者呈现同类型的抖动就基本可以判断问题出在地质条件而不是处理流程这个判断在真实资料上是很难下的。4.2 时间偏移与深度偏移速度模型错了构造会把构造推到哪时间偏移用的是均方根速度场深度偏移用的是层速度模型。两者在正演数据上的表现差异是理解偏移算法最直观的方式。合成一条已知的倾斜界面用一个偏高的偏移速度去做时间偏移同相轴会被推得比真实位置更陡甚至出现过偏移的钩状用一个偏低的偏移速度同相轴会被压平倾角变小。深度偏移对速度误差更敏感因为它直接换算深度速度偏 5%成像深度就偏 5%而且是随深度累积的。这个实验的意义在于它把偏移速度模型不准这个抽象说法变成了一个可以量化的误差传递关系。做实际资料处理时看到成像道集上残余的曲率能大致判断速度该往哪个方向调、调多少。这种手感靠真实数据是练不出来的因为真实数据没有标准答案。4.3 AVO 响应各向异性、吸收、薄层如何改写截距-梯度关系AVO 是这套工具用得最深的场景之一。常规 AVO 用 Shuey 近似反射系数写成 R(θ) ≈ A B·sin²θA 是截距近零偏移距反射B 是梯度振幅随偏移距的变化率。把不同偏移距上的振幅拟合成一条直线A、B 就成了岩性和流体识别的特征量。但在四类复杂情况下这条直线会弯曲甚至整体平移。各向异性会引入一个与 δ 成正比的附加项表现为小角度范围内梯度被系统性抬升吸收会改变相位谱使远偏移距的振幅出现额外的衰减趋势拟合出的 B 偏小薄层会让顶底反射在时间上相互干涉振幅随偏移距的变化不再是单界面的响应而是顶底两个界面干涉的结果表现为 A、B 交会图上的点沿着一条异常的线分布这条线在常规解释里对应的岩性是完全错误的。这正是原料里强调模拟复杂介质的 AVO 效应各向同性、各向异性、吸收、梯度、薄层的原因——这四种机制如果不分开建模是没法在实测数据上归因的。4.4 处理流程的闭环为什么参数测试要在已知答案上做把正演、处理、对比串起来就构成了一个闭环建模提出一个地质假设→ 正演计算这个假设下的波场→ 处理用各种参数去恢复它→ 对比看恢复出来的结果与模型差多少。这个闭环里每一步的误差来源是可以隔离的。举个具体例子假设某个工区远偏移距的能量很弱看起来像是采集不够。在合成数据上把同一套观测系统套用到模型上跑一遍如果合成数据在远偏移距也弱那是介质本身的速度和界面结构导致的跟采集无关如果合成数据远偏移距很强那问题可能出在实测资料的处理链路上比如 NMO 拉伸切除过多、或者吸收补偿过度。没有合成数据这一步这两种可能性靠经验是分不开的——这就是指导解释、验证解释结果的实际含义。5. 避坑与排查正演结果和实测对不上时先查这五处正演跑出结果不难跑出能用来支撑解释的结果才是真正的门槛。下面五条是这类工具在实际使用中最高频的翻车点每一条都对应一个可复现的现象和一条明确的排查路径。第一条同相轴后面拖着一串涟漪现象是同相轴主瓣后面跟着一系列高频、衰减很快的小振荡看起来像信号被加了回响。原因是网格频散差分算子只在低频段逼近真实导数当波长接近网格尺寸时相速度开始偏离真实速度低频走得快、高频走得慢波包被拉长成尾巴。解决方向有三个按性价比排序提高空间差分阶数从 4 阶到 8 阶效果立竿见影、减小网格步长代价是计算量按平方涨、降低震源子波主频最省事但会损失分辨率。判断是否已经满足要求看每个最短波长被多少个网格点采样一般要保证 810 个点。第二条波场快照的边缘出现从外面回来的波现象是快照边缘出现与内部波场不连续、方向指向模型内部的同相轴时间剖面上表现为深部出现本不该存在的倾斜能量。原因是人工边界把能量全反射回来了PML 或吸收层没有生效。常见成因有两类吸收层太薄对大角度入射的吸收效果差或者模型本身贴边界面直接落在网格边界上吸收层没有空间。解法是把吸收层加厚到 1020 个网格并检查模型四周是否留了缓冲带。这个坑的隐蔽之处在于反射回来的能量常常看起来挺像真的反射只有对比不同吸收层厚度下的结果才能确认。第三条跑到几千步之后波场数值飞掉现象是快照上出现密集的、幅度极大的条纹或者直接变成 NaN。原因基本只有一个dt 超过了 CFL 允许的上限。而 dt 之所以会超多半是因为算上限时用的 Vmax 是地表速度而模型底部塞了一个高速基底。解法是遍历整个模型包括吸收层和边界填充区取最大速度重新计算 dt。另外当模型里有极薄的高速夹层时最小网格步长和最大速度可能同时出现在同一位置这种组合最危险建议把这种网格在建模阶段就做等效处理。第四条合成记录整体比实测深或浅了一截现象是同一套反射的同相轴合成记录比实测记录在时间轴上整体偏移了几十毫秒构造形态倒是对的。原因是深度域与时间域混用测井是深度域的模型网格是深度域的而合成记录经时间域输出中间有一次时深转换而输入层位如果本身就是时间域解释成果就会多换算一次。另一个高频成因是声波时差单位没换算μs/ft 当成 μs/m 用速度整体差 3.28 倍这种情况连构造形态都会跟着变形。排查顺序建议先核单位、再核域、最后核基准面。第五条AVO 梯度反号或者交会图整体平移现象是合成数据拟合出的梯度项与实测资料符号相反或者截距-梯度交会图上所有点整齐地偏到另一侧。原因通常是约定不一致震源子波极性SEG 正极性约定、密度单位g/cm³ 与 kg/m³ 混用、反射系数的符号定义任一处不一致都会导致整体反号。这类问题不会让剖面看起来有任何异常只会让最终的岩性解释完全反转。最经济的排查办法是在零偏移距上先对齐把合成记录与实测记录的近偏移距道做互相关看时移和极性确认之后再谈 AVO。6. 设计取舍为什么有限差分 交互建模是这类工具的合理答案正演算法有好几类每一种都有它舒服的区间。方法复杂介质适应性波场快照高频/大模型计算成本射线追踪弱高频近似难处理绕射、薄层不能直接输出强低反射率法中层状介质横向变化受限部分可输出强低有限差分强任意横向变速、各向异性、吸收完整输出弱受频散和网格限制高表格里的最后一列是关键。有限差分的成本高换来的是什么介质都能算、什么波型都能看的通用性。对于以验证解释为目标的工作通用性比效率重要得多你事先并不知道异常波是绕射、是转换波还是薄层调谐只有保留完整波场才能在快照上把它认出来。射线追踪算得快但它只给走时和振幅遇到溶洞这类尺度接近波长的目标射线假设直接失效。这就是这套工具选择有限差分作为唯一计算引擎的原因——它是一个不会预先替你排除可能性的方法。交互式建模是第二个取舍。把建模、计算、查看、处理做成四个可以来回切换的功能包代价是每个包单独拿出来都不算最深好处是迭代速度快。一次正演的完整循环——改模型、算波场、看快照、抽道集、做偏移——如果在四五个工具之间倒数据中间任何一个格式转换出错都会让人怀疑人生。集成在一套界面里模型改了立刻能重算重算完立刻能在同一个窗口里对比波场分量、提取振幅、算属性这种反馈速度对学习型使用场景是决定性的。第三层取舍在 2D 本身。2D 正演假设介质在垂直剖面方向不变这在构造走向稳定的工区成立在三维构造上不成立侧向绕射和离面能量全部丢失。但 2D 的计算量比 3D 低一个量级以上迭代次数可以多一个量级对于理解波场机制、设计观测系统、验证处理流程这几件事2D 的性价比远高于 3D。真正需要 3D 的场景是面元设计、方位各向异性、以及绕射体精确成像那些问题 2D 确实答不了。我从那以后每次拿到一套正演结果都会先做三件事核一遍 CFL 和每波长网格点数看一眼快照边缘有没有回波再把零偏移距道和实测对齐一次极性。这三步花不了十分钟但能挡掉绝大多数结论看着对、其实是数值假象的坑。这套资料把建模、计算、查看、处理四个环节的参数和界面都摊开摆在那里对着自己的工区跑一遍比读十篇原理文章都快——资源就在手边拿来当标准答案用比自己拍脑袋靠谱得多。希望帮到你。本文还有配套的精品资源点击获取