简介:Tesseral 2D模型正演软件功能与特点的演示PPT,面向地球物理勘探学习者、处理解释人员及需要开展观测系统设计与方法验证的技术人员。内容围绕地震勘探多解性问题,说明模型正演在采集参数设计、波场传播认识、处理流程测试和解释结果验证中的作用,并逐一介绍Model builder建模、Computation Engine合成记录与波场快照、Viewer可视化、Processing Block速度分析及叠前叠后偏移等功能。PPT结合有限差分正演、P波与SV波计算、梯度/声学/弹性及各向异性介质、吸收衰减和AVO效应等知识点,辅以溶洞、砂泥互层模型示例,便于读者理解复杂地质体的地震响应。资源包共1个PPT文件,约2MB,结构紧凑,适合作为软件功能概览和正演流程入门参考;已有226人学习,可帮助快速建立Tesseral 2D整体认识并对照项目需求查阅。
1. 复杂构造多解性绕不开:Tesseral-2D 正演软件整套资料拆解
解释成果和实钻井位对不上,很多队伍第一反应是偏移参数没调好,回头把处理流程翻来覆去改,改到最后一版和第一版几乎没差。真正的根子往往在前面:观测系统选得对不对、复杂介质下常规处理的假设还成不成立,这些在上机处理之前就该用正演验证。这份 Tesseral-2D 正演软件的资料包,讲的就是怎么在计算机里先造一个已知答案的地质模型,把波传一遍,看记录长什么样,再回头去定采集参数、处理流程和解释方案。它适合做观测系统设计的地球物理人员、要验证 AVO 和多分量响应的方法研究人员,以及刚接触正演、想把有限差分这套东西跑通的新手。
2. 四个功能包怎么串成一条链:从 LAS 曲线到合成炮集
Tesseral-2D 的结构不是一堆并列菜单,而是一条单向数据链:Model builder 造模型,Computation Engine 算波场,Viewer 看结果,Processing Block 把合成记录当成真实资料再处理一遍。理解这条链的输入输出格式,比记住界面按钮重要得多,因为实际项目里卡住人的地方,八成都在两个包的衔接处。
2.1 Model builder:密度-速度模型的构建逻辑与 LAS 曲线接入
Model builder 的产物是一个二维密度-速度剖面,纵横向离散成网格。地震正演的精度和成本都由这个网格决定:网格太粗,薄层和绕射体直接被抹平;网格太细,计算量按网格数乘时间步数增长,一台普通工作站跑一个 2000×1000 的弹性模型能耗掉大半天。
常见做法是先用层状框架搭骨架,再往里填速度。层状框架给出每套地层的顶底界面和横向起伏,速度场可以选择常速、垂向渐变、横向渐变,或者直接由 LAS 测井曲线插值出来。LAS 文件接入这一环要注意:测井曲线给的是深度域的 P 波声波时差,软件读进来后要先换算成速度,再按层位标定把深度对应到模型网格上。如果井曲线和地震层位在深度上有系统偏移,插出来的速度场会整体错位,后面算出来的反射同相轴对不上地质界面。
| 模型要素 | 常见取值/做法 | 影响 |
|---|---|---|
| 网格横向采样 | 2~5 m | 决定横向分辨率与绕射体保真度 |
| 网格纵向采样 | 1~2 m | 薄层调谐效应能否被刻画 |
| 速度来源 | LAS 曲线插值 / 分层常速 / 梯度 | 决定波场真实程度 |
| 密度来源 | 测井密度曲线或 Gardner 关系换算 | 决定反射系数与 AVO 响应 |
| 介质类型 | 声学 / 弹性 / 各向异性 / 吸收 | 决定计算量与波型完整性 |
介质类型的选择是最容易过度设计的一步。只关心纵波构造形态,用声学介质就够了,计算快、波场干净;要看转换波、要研究极化,必须上弹性;薄互层那种非单调的强各向异性响应,得上各向异性介质,否则横波分裂和时差特征根本出不来。我的习惯是先跑声学版本确认几何关系对不对,再逐级加复杂度,每加一级都留一版结果做对比,这样出问题时能立刻定位是哪一层假设引入的偏差。
2.2 Computation Engine:有限差分求解 P/SV 波的参数设定
Computation Engine 用有限差分在时间-空间域推进波场,核心是三个参数:空间差分阶数、时间步长、边界吸收条件。空间差分阶数常见 4 阶到 8 阶,阶数越高,对高速层的频散压制越好,但每个网格点的计算量也线性上升。时间步长受稳定性条件约束,步长取大了波场直接发散,取小了白白烧时间。
# 典型的二维弹性正演参数配置(示意,字段名以软件界面为准) GRID_NX = 2000 # 横向网格数,覆盖剖面长度 GRID_NZ = 1000 # 纵向网格数,覆盖目的层深度 DX = 2.0 # 横向网格间距(m) DZ = 1.0 # 纵向网格间距(m) DT = 0.0005 # 时间步长(s),需满足稳定性条件 T_MAX = 6.0 # 记录长度(s) FD_ORDER = 8 # 空间差分阶数 BOUNDARY = "PML" # 吸收边界,抑制人工反射 SOURCE_TYPE = "P" # 震源类型:纵波震源 SOURCE_COORD = (100, 2) # 震源位置(网格坐标) RECEIVER_SPACING = 10 # 检波器道间距(网格数) MEDIUM = "ELASTIC" # 介质类型:声学/弹性/各向异性/吸收这段参数里最需要盯住的是 DT 和 DX、DZ 的比例关系。稳定性条件近似要求 DT 小于网格间距除以最大速度再乘一个与差分阶数有关的系数,工程师的习惯做法是先按最大层速度算一个上限,再取上限的 0.6~0.8 倍做实际步长,留出安全余量。频散是另一个隐性问题:高速层里波长被网格采样得太稀,波前面会出现一圈圈波纹,看起来像真实的高频干扰,其实是数值频散,解决办法是提高差分阶数或加密网格,绝不是去加滤波。
波场快照是这一环最有价值的输出。资料里给的声波模型例子把波前在 t1=0.350 s、t2=0.510 s、t3=0.810 s、t4=1.110 s、t5=1.420 s 这几个时刻定格下来,首波传播到哪、在哪个界面发生透射和反射、绕射体产生的双曲线顶点落在哪,一眼就能看出来。做观测系统设计时,我会先看快照找出目的层反射波和干扰波在时间上的分离窗口,再回头调最大炮检距和记录长度,而不是凭经验拍一个数。
2.3 Viewer:波场快照与振幅属性的对比查看方式
Viewer 的作用是把计算结果变成能判读的图件。它支持对比不同波场分量、提取信号振幅、读取相关属性值,也能把数据导出到外部系统继续处理。实际用起来,最有价值的是同一模型不同参数下的并排对比:改一个差分阶数、改一个吸收边界,其他都不动,两幅快照叠着看,差异立刻显形。
幅值提取这一步要留神显示增益。软件为了让人看清弱反射,默认会做自动增益,同一张图里强反射被压、弱反射被抬,看着都对,但幅值关系已经不能用来做 AVO 判断了。要看真实振幅,必须关掉自动增益,用统一的色标范围出图。我见过有人拿自动增益后的快照去量振幅比值,得出的 AVO 趋势和理论完全相反,白折腾两天。
2.4 Processing Block:合成炮集的速度分析、CMP 叠加与偏移
Processing Block 让这套软件从"看波场"升级成"预演处理流程"。它能把合成炮集当真实资料走一遍:速度分析、CMP 抽道、叠加、叠前叠后时间偏移、叠前叠后深度偏移。这一环的意义在于,你事先知道模型的真实速度,可以拿处理出来的速度谱和真值对比,看看在某个复杂构造下,常规速度分析会不会把速度估偏,偏多少,叠出来的成像在哪个位置失真。
流程上建议按这个顺序走:先抽 CMP 道集看道集是否平直,再用真速度做一次叠加当参照,然后用常规速度分析重做一遍,两版叠加剖面叠在一起看差异。偏移环节重点看绕射收敛和断面归位,模型里埋的溶洞、砂泥互层、横向变速薄层沉积,都是检验偏移算法边界的现成考题。资料里那套横向参数变化的薄层沉积模型,用层爆炸方式模拟时薄层表现为多相位反射,这个特征在叠加剖面上会变成一个复波,偏移之后能不能收敛回真实层位,直接反映处理流程在该构造下的适用性。
3. 介质模型与 AVO 响应的预演:把解释假设放到已知答案里检验
正演真正的杀伤力不在出一张漂亮图,而在于它给了一个已知答案的考试卷。你可以在模型里埋一个油气藏,改变它的厚度、速度、密度、各向异性强度,看合成记录上的 AVO 曲线怎么变;你也可以把同一套解释假设分别套在正演数据和实测数据上,看偏差出现在哪一步。这一章把介质建模和 AVO 预演拆开讲,因为它们对应两类完全不同的使用人群。
3.1 梯度、声学、弹性与各向异性四类介质的建模范式
四类介质不是难度递进关系,而是对应四类地质问题。声学介质只保留纵波,适合构造形态验证和观测系统几何设计,计算量最小。弹性介质同时算 P 波和 SV 波,用于转换波研究和多分量采集设计,成本大概是声学的好几倍。梯度介质描述速度随深度或横向连续变化,垂向渐变多见于压实作用明显的碎屑岩剖面,横向渐变则对应古地貌、盐丘侧翼这类构造。各向异性介质的参数最多,除了速度密度还要给各向异性系数,薄互层引起的等效各向异性是典型场景。
建模时的血泪经验是:别一上来就堆最复杂的介质。先把目标问题想清楚——如果目的是检验某一层顶面在叠加剖面上会不会被强反射掩盖,声学模型足以回答;如果是判断转换波能不能把某套薄储层分辨出来,就必须上弹性。用错介质类型的代价不是算得慢,是算出来的结论不适用于你要回答的问题,那种错误在报告里最不容易被发现。
3.2 吸收衰减与初至到达时间场的设置
吸收衰减用品质因子 Q 描述,Q 越小衰减越强。真实地层里,浅层疏松介质的 Q 值低,高频成分被吃掉得快,深层致密地层 Q 值高,频带保持得好。正演时如果不加吸收,合成记录的高频成分比实测数据丰富得多,分辨率看着很诱人,实际对比时会发现处理参数按正演结果调,用到实测资料上就偏了。常见做法是给每套地层赋一个经验 Q 值,先做一版无吸收、一版有吸收,对比频带差异,再决定后续处理给多大的反 Q 补偿量。
初至到达时间场在观测系统设计里是硬指标。软件会输出初至时间场,把每个检波点位置的最早到达时间给出来。这个场用来核算最大炮检距是否够用、折射干扰会不会盖住目的层反射、静校正量级大概多大。设计三维工区之前先在二维剖面上把初至时间场看一眼,能省掉很多后期返工。
3.3 薄层、各向异性与吸收三种效应下的 AVO 正演对比
AVO 正演是这套软件的重点场景。做法是固定一个目的层,分别改厚度、各向异性强度、吸收强度,各算一套合成道集,然后抽近、中、远道看振幅随炮检距的变化。三类效应的表现很容易混淆:
- 薄层调谐会让 AVO 曲线出现周期性起伏,厚度变化时截距和梯度同时摆动,看起来像岩性变化,其实是调谐。
- 各向异性会改变远道振幅的衰减速率,尤其是薄互层引起的等效各向异性,远道信息被系统性抬升或压制。
- 吸收衰减让高频在远道损失更多,表现为频谱随炮检距变窄,如果只盯着振幅不看频谱,容易误判成 AVO 异常。
这三种效应叠在一起时,单靠实测道集基本分不开。正演的价值就是先造出单独效应的道集当标尺,再逐步叠加,看组合效应长什么样,最后拿实测去比对,判断观测到的异常最可能由哪种机制主导。资料里提到的"模拟复杂介质的 AVO 效应",落地点就在这套对照实验上,而不是算一个道集就完事。
3.4 溶洞、砂泥互层与薄层沉积三类模型的复现步骤
三类典型模型各有各的考法。溶洞模型考的是绕射和散射:在均匀速度背景里放一个低速度异常体,尺度从几个网格到几十个网格变化,看绕射波的能量、双曲线顶点的位置和横向分辨率。砂泥互层模型考的是薄层调谐:按实际砂泥厚度比建层,砂层速度密度给实测值,看反射波形随砂层厚度变化的调谐曲线。横向变速薄层沉积模型考的是成像:横向速度变化会让常规的共中心点叠加假设失效,正演出来能看到反射同相轴在横向上错断。
复现步骤一般是固定的四步:
- 在 Model builder 里搭背景速度场,先只放一个目的体。
- 设置震源、检波器排列和记录长度,先跑一版看波场快照确认传播正常。
- 换介质类型或改参数,生成多版合成炮集。
- 进 Processing Block 走叠加与偏移,比对成像结果与模型真值。
每一步都保存工程文件,别嫌麻烦。中途改了一个参数发现结果不对,想退回上一版,没有工程存档的话只能从头再搭,那种感觉相当难受。
4. 上机避坑:五条踩过的坑和对应解法
这套软件参数多、耦合强,新手最容易在几个固定位置翻车。下面五条都是从实际跑模型的过程中攒下来的,每条按现象、原因、解决写清楚。
4.1 波场快照出现异常波纹
现象:快照里波前面外围有一圈圈规则的高频波纹,越靠近高速层越密。
原因:数值频散。高速层中的波长被网格采样点数不足,差分近似失效,不是真实物理现象。
解决:先把差分阶数从 4 阶提到 8 阶看是否减轻;仍然明显就加密网格间距,同时按稳定性条件重算时间步长。不要用滤波手段掩盖,那会把真实的高频信息一起损伤。
4.2 模拟中途波场幅值发散
现象:跑到一段时间后幅值指数增长,快照一片乱,最终溢出中断。
原因:时间步长超过稳定性上限,或者吸收边界参数设置不当导致边界反射叠加。
解决:按最大层速度重算时间步长上限,取 0.6~0.8 倍;检查吸收边界层数是否与网格间距匹配,层数太少吸收不干净。发散一旦出现必须停机重跑,不要指望用截断幅值抢救。
4.3 LAS 曲线接入后反射界面整体错位
现象:合成记录上的反射时间与模型设定的层位深度对应不上,整体偏移一个固定量。
原因:测井曲线深度基准与模型顶深基准不一致,或者声波时差换算速度时单位没统一。
解决:接入前先确认曲线深度起点、采样间隔和单位(微秒每米还是微秒每英尺),换算后在模型里放一个已知深度的标志层核对。这个核对步骤只用几分钟,能省掉后面几小时的返工。
4.4 处理流程预演结果和实测规律不一致
现象:用 Processing Block 跑出来的速度谱、叠加效果,跟实际资料的表现差得远。
原因:模型过于理想化,缺少噪声、吸收、各向异性等实际因素,或者合成炮集的观测系统与实测不一致。
解决:先核对观测系统参数(道间距、最大炮检距、最小炮检距)是否与实测工区一致,再逐项加入吸收和随机噪声。正演的作用是给趋势和边界,不是复刻实测,比较时要看规律方向而不是绝对数值。
4.5 大模型跑不完或内存告急
现象:网格数一加大,计算时间从几小时涨到十几个小时,甚至直接内存不足退出。
原因:计算量与网格数、时间步数、介质复杂度的乘积成正比,弹性各向异性模型尤其吃资源。
解决:先用粗网格和短记录长度跑通流程,确认参数无误后再加密;把大剖面切成若干段分别计算;声学能回答的问题不要用弹性去跑。资源规划应在建模之前做,而不是算到一半才发现跑不动。
5. 把一个正演模型变成可复用模板:参数归档与结果核验的固定动作
这套软件用久了会发现,真正省时间的不是算得快,而是每次都能从上一版模型里继承大部分设置。我的做法是给每个正演项目建一个参数档案,把网格参数、介质类型、震源检波器配置、介质分层速度密度表、Q 值表、差分阶数、时间步长全部落成一张表,模型文件另存一版带日期后缀。下次遇到类似构造,拷一份档案改几个字段就能起跑,比在界面里翻菜单重设快得多。
核验环节有一个固定动作值得坚持:每次算完,先不看合成记录,先看三个东西——波场快照里波前是否连续无频散、初至时间场是否随炮检距单调增大、最深层反射的到达时间是否与按层速度手算的双程时一致。这三项任何一项不对,后面的道集和偏移结果都不用看了,直接回建模环节查。手算双程时那一步最容易被跳过,但它是唯一一个不依赖软件自身逻辑的独立校验,出了系统性错误时只有它能兜住。
| 归档项 | 记录内容 | 用途 |
|---|---|---|
| 网格参数 | NX、NZ、DX、DZ、DT | 换工区时快速缩放 |
| 介质配置 | 类型、分层速度密度、Q 值 | 复现同类地质问题 |
| 观测系统 | 道间距、炮检距范围、记录长度 | 与实测工区对齐 |
| 求解参数 | 差分阶数、边界条件、震源类型 | 排查频散与稳定性 |
| 核验结果 | 手算双程时、初至时间场范围 | 独立校验计算正确性 |
进阶一点的用法是把正演结果反过来约束处理。比如在某个横向变速区,先用正演算出常规速度分析在该区的理论误差量级,然后拿实测资料的速度谱去比对,若误差量级吻合,说明实测速度偏差主要来自方法假设而非人为拾取,这时候该做的是换成像方法,而不是反复重拾速度。这个判断链条一旦建立起来,正演就从"出图工具"变成了处理决策的依据。
还有一个习惯:任何一版准备写进报告的合成记录,我都会把对应的参数档案和核验记录一起附上。有一次隔了两个月回头复现一个溶洞模型的结论,靠的就是当时那份档案,不然光凭一张图根本还原不出当时的差分阶数和边界设置。从那以后我每次跑正演都强制走一遍参数归档加三项核验,希望帮到你。
本文还有配套的精品资源,点击获取