VDPS多分量地震资料处理软件培训手册.doc

内部资料 请勿外传 VDPS 多分量地震资料处理软件多分量地震资料处理软件 培训手册 恒泰艾普石油天然气技术服务股份有限公司 20112011 年年 5 5 月月 1010 日日 目录目录 简介简介.1 数据准备数据准备.3 水平分量旋转水平分量旋转 .3 三维面元化三维面元化.5 速度建模速度建模.7 为速度分析创建为速度分析创建 ACP 超道集超道集.8 转换波叠加速度分析转换波叠加速度分析.9 转换波叠加速度解释转换波叠加速度解释.10 转换波正常时差校正转换波正常时差校正.12 纵横波速度比解释纵横波速度比解释.13 转换波偏移速度解释转换波偏移速度解释.18 偏移作业提交以及监控管理偏移作业提交以及监控管理.19 简介简介 用途用途 LD-VDPSTMprocess(以以下下简简称称 VDPS) 是一个多分量地震数据处理软 件包，整合了多分量资料常规处理，转换波参数估计，转换波偏移等多分 量地震数据处理方法，具有完备的二、三维多分量地震资料处理功能。 本软件应用范围广泛，能用于： 复杂地质构造区域成像，如气云下成像，玄武岩下成像，透过泥岩 底辟，泥火山成像； 提高浅层分辨率，进行更详细的浅层结构成像； 纵波反射弱地区的目标成像，由于纵波反射弱的地区，转换波会有 较强反射，能很好成像； 提供弹性参数。能提供纵横波速度比，各向异性参数。可分辨岩性， 如砂岩/页岩，检测浅层水流，岩性构图。可对比解释纵横波处理 资料。 本软件系统拥有众多优势技术，包括多分量资料观测系统定义的图形显 示，多分量资料旋转和极性反转，高分辨率宽带子波约束反褶积，模拟退 火和遗传算法剩余静校正，纵横波速度比分析，高精度转换波叠前 Kichhoff 时间偏移等技术，实现了处理与解释构造形态的迭代、处理解释 速度模型的迭代，有效地保证了处理和解释成果的质量，尤其在复杂地表 低信噪比资料和高分辨率处理方面具有优势。 特点特点 基于转换波多参数理论，使用了多种速度，速度比和各向异性参数，并 将这些参数从多分量数据抽离出来，用之建立准确的各向异性速度模 型。在此基础上发展了转换波参数分析、建模方法。精确的参数建模， 为分别建立适应转换波资料处理的精确参数场提供强有力的技术保障。 各向异性参数模型具有岩性以及油气方面的应用，例如利用纵横波速度 进行岩性分辨，检测浅层水流，进行岩性构图，研究各向异性。 基于层状各向同性(VTI)介质假设，对转换波绕射方程高次项进行优化， 使得波场传播更符合实际情况。 基于转换波绕射曲线方程，利用多种参数进行转换波偏移，显著地提高 了处理效果。科学的偏移算法，使处理结果断层清晰，构造准确，更 易于构造解释。 具有良好的保幅特性，能够进一步完成各种以波形为基础的反演处理。 提供了具有独立特性的纵横波剖面，具有更详细的浅层结构成像。 本软件支持 Window 和 Linux 操作平台； 系统数据接口通过导入、导出完成各种格式和类型数据交换。从而支持 多种标准格式以及其它处理系统规定格式和自由格式的数据。 总控统一、界面风格统一、数据接口统一，界面友好，性能稳定，操作 简单，国际通用方式。 支持网络分布式计算、并行处理（PC-Cluster），并集交互和批量于一 体。 数据准备数据准备 将非 DSK 格式数据转化为 DSK 格式。 对多分量地震数据体作必要的预处理，例如道编辑、静校正等。 将地震数据体分别储存为水平分量数据和垂直分量数据。 炮点和检波点的坐标、高程转换为统一的单位(m)。 利用 Geometry 观测系统定义观测系统。把野外观测 SPS 格式文件或 UKOOA 标准格式文件，解编加载到系统中交互显示关于测线几何排列形态信 息。计算炮点、检波点静校正量的参数和道编辑信息等。根据任务要求的面元 大小，定义面元大小，并根据交互形成的图件检查面元大小是否一致、覆盖次 数是否均匀、道集内各道的方位角和炮检距分布是否均匀，从而检查观测系统 定义的正确性。将以上信息以文件的方式存放在系统几何库中。 利用 setthd 模块将几何排列信息直接存放到地震数据道头中去。 水平分量旋转水平分量旋转 实际采集中，检波器水平分量坐标是按实际坐标 x,y 排列的，极化方向与径 向分量(Radidad,炮检距方向)不一致；从而横交线(cross-line，y)也与法线 (transversem，t)不一致，因此需要将水平分量进行旋转才能得到处理所需的径 向水平分量。 启动模块：启动模块： （1）在 DPSpro 流程流中，该模块界面为： （2）在命令行中，键入：cxacp if 输入文件名前缀 of 输出文件名 前缀 ，就会出现参数卡： 参数说明：参数说明： 参数名类型取值范围功能备注 DSKIChar * 第二个输入文件名 Second File for input DSKOChar* 第二个输出文件名 Second File for input ANGLEfloat 旋转角度 Angle for rotation 运行后得到沿跑检径向方向 R 和切向方向 T 的数据。 三维面元化三维面元化 转换波资料需要选排成共转换点道集，从而使道集内的各道具有相同的转 换点。通过计算转换点坐标将数据划分面元得到共转换点道集。为了简化计算 和提高计算效率，实际工作中只是针对目的层确定统一的纵横波速度比，依据 转换波时差方程计算出目的层转换点，依据该转换点抽取 ACP 道集。 计算地震道 ACP 值、ACP 线号以及相应的 AOFFSET 值，形成 ACP 道集。 启动模块：启动模块： （1）在 DPSpro 流程流中，该模块界面为： （2）在命令行中，键入：cxacp if 输入文件名前缀 of 输出文件名 前缀 ，就会出现参数卡： 参数说明：参数说明： 参数名类型取值范围功能备注 GAMMAdouble 纵横波速度比 Vp/Vs ratio CDPDXdouble CDP 间隔 distance between two cdp trace LDXdouble 测线间隔 distance between two lines X0double 3 维测区起始 X 坐标 x coordinate of origin in 3D surevy Y0double 3 维测区起始 Y 坐标 y coordinate of origin in 3D surevy X1double 3 维测区炮线 X 坐标 x coordinate of shot line in 3D surevy Y1double 3 维测区炮线 Y 坐标 y coordinate of shot line in 3D survey DOFFSETdouble 偏移距间隔 offset distance between trace ORICDPdouble 起始 ACP 编号 original acp number ORILINEdouble 起始 线 编号 original line number 运行模块得到修改道头字后的数据，再运行 SISORT 模块可以重排数据得 到 ACP 道集。 速度建模速度建模 建立初始的速度模板进行速度分析。 启动模块启动模块 在命令行中，键入：cxvelmb 就会出现参数卡： 参数说明：参数说明： 参数名类型取值范围功能备注 VelNameChar * 速度文件名 Velocity file Name FisrtCdpInt 起始 CDP 号 First CDP number in the line LastCdpint 终止 CDP 号 Last CDP number in the line IncCdpint CDP 号增量 Increment of CDP No Timefloat 时间采样点 time samples Velocityfloat 相应速度 velocity corresponding time samples G0float 相应 g0 G0 corresponding time samples Gefloat 相应 ge Ge corresponding time samples Chifloat 相应 ge Chi corresponding time samples 为速度分析创建为速度分析创建ACP超道集超道集 为分析速度，创建 ACP 超道集来提高信噪比。 启动模块：启动模块： （1）在 DPSpro 流程流中，该模块界面为： （2）在命令行中，键入：gsuper if 输入文件名前缀 of 输出文件名 前缀 ，就会出现参数卡： 参数说明：参数说明： 参数名类型取值范围 功能备注 PHdouble CDP：起始数，结束数，增量，每组个数 CDP：First,Last,Increment,Group THdouble 线号：起始数，结束数，增量，每组个数 Line：First,Last,Increment,Group OFFSdouble 偏移距：起始量，结束量，增量 Offsets：First,Last,Increment 转换波叠加速度分析转换波叠加速度分析 利用转换波参数模型进行转换波叠加速度分析，得到转换波速度谱。 启动模块：启动模块： （1）在 DPSpro 流程流中，该模块界面为： （2）在命令行中，键入：vcscan if 输入文件名前缀 of 输出文件名 前缀 ，就会出现参数卡： 参数说明：参数说明： 参数名类型取值范围功能备注 VELFILEChar * 速度文件名 Velocity file Name SCANPAR M Int 扫描参数:起始速度，扫描个数， 扫描增量，平滑长度 Vel Start,Scan Number,Increment, LenSmooth MinOffsetfloat 最小偏移距 MinOffset MaxOffsetfloat 最大偏移距 MaxOffset 转换波叠加速度解释转换波叠加速度解释 由于下行纵波和上行转换横波的不对称路径，转换波时距曲线不再是双曲 线，动校正为非双曲线动校正。利用转换点方程和转换波时差方程，推导了转 换点方程和转换波时差方程。将绕射曲线被分解为两个分量：第一个部分仅仅 与速度有关，用于描述绕射曲线中时差的双曲线表现；另一个分量依赖于速度 和其他参数，用于描述时差的非双曲线表现。通过交互分析得到转换波速度， 纵横波叠加速度比以及各向异性参数。其中转换波速度控制近偏移距双曲时差， 纵横波叠加速度比控制因非对称射线路径产生的非双曲时差，各向异性参数控 制着远偏移距非双曲时差。 启动模块：启动模块： 在命令行键入 Psvnmopk，就可启动 Psvnmopk 转换波速度解释分析工具。 弹出工区选择对话框，选择工区、测线；如上图所示。选择工区后弹出当 前工区下的测线，选择测线弹出当前测下已经准备好的速度分析数据。选择已 经准备好的速度分析数据，选择初始参考速度库(PSV 库)。选择后点击 OK。系 统弹出速度分析系统交互主界面。如下图所示： 用户可以在左侧面板中调整转换波的叠加速度，然后依次调节和， eff 直至右侧面板中 ACP 道集的时矩曲线与同相轴相吻合。 初始时默认为单谱模式拾取，速度谱更新工具按钮和编辑切换工具按 eff 钮自动隐藏；当全谱模式拾取时谱拾取工具按钮、交互拾取工具按钮和 C V eff 各项异性谱拾取工具按钮自动隐藏。当拾取谱工具按钮按下时，开始扫描速 度谱，各个工具按钮不能激活，当扫描完毕时才自动激活。可以采用单一谱模 式解释速度，各个道集逐一解释，实时更新谱、谱和谱，当解释一 C V C V eff 轮后再解释，最后解释谱，逐一迭代依次避近，直到最后所有谱道集 C V eff 实时更新很少变化时为止； 也可以采用全谱模式解释，当在谱面板解释 C V 速度后，同时调整和，实时更新谱，实时解释。 C V eff 更详细操作请参考PSVNMOPK 转换波叠加速度解释工具用户使用手册. doc。 转换波正常时差校正转换波正常时差校正 根据转换波绕射时距曲线公式，采用转换波多参数文件进行转换波正常时 差校正。 启动模块：启动模块： （1）在 DPSpro 流程流中，该模块界面为： （2）在命令行中，键入：cxnmocw if 输入文件名前缀 of 输出文件 名前缀 ，就会出现参数卡： 参数说明：参数说明： 参数名类型取值范围功能备注 ACTIONenumerate 校正方向：正向/反向/偏移速度 NMO Direction Forward/Inverse/MIGVEL VELFILEchar 速度模型文件 The velocity model file obtianed from xva 参数卡输入 运行模块得到正常时差校正后的 ACP 道集 纵横波速度比解释纵横波速度比解释 纵横波速度比（）的作用。一、作为转换波叠前/后时间偏移的必要参数； 0 二、纵横波波组匹配的需要。 转换波动校正对纵横波速度比的变化不敏感。纵横波速度比需要利用纵波 数据和转换波数据对比，交互地解释纵横波速度比曲线； 扫描速度比谱，更加直观方便进行分析； 根据速度比拉伸/压缩转换波叠加剖面，检查纵波剖面和转换波剖面匹配情 况。 模块启动：模块启动： 在命令行键入 GammaPick，就可启动 GammaPick 转换波垂向速度比解释分 析工具。 启动 GammaPick 时，会弹出工区选择对话框，选择工区、测线，如下图所 示。选择工区后弹出当前工区下的测线，选择测线弹出当前测下纵波和转换波 叠后数据。选择叠后数据，也可以选择加载已存在的库(pfn 库)。选择后点击 0 OK。 系统弹出分析系统交互主界面。如下所示。 0 在叠后剖面上点击选择地震道，扫描相应谱，如上图所示。根据谱中 0 0 能量团汇聚情况点击拾取曲线。在叠后剖面中，红线标记的地震道为当前 0 谱所在 CDP 点位置。黄线标记为已经拾取过曲线的 CDP 位置。 0 0 拾取曲线后，可点击按钮，按照拾取的曲线拉伸/压缩 PS 波叠后剖 0 0 面到 P 波旅行时，查看拾取的曲线是否正确。 0 拾取曲线结束后，点击启动剖面窗口，如下图所示。 0 0 剖面窗口剖面窗口 0 根据叠置显示的 P 波叠后剖面，拾取层数据，如下图所示。 层拾取 通过层位控制，对曲线的取值进行解释，如下图所示： 0 解释曲线 0 可点击查看时间切片，检查水平空间展布，如下图所示 0 时间切片 0 编辑和修改完毕后，保存曲线到参数库。 0 更详细操作请参考VDPSgammapick 纵横波速度比分析使用手册。 转换波偏移速度解释转换波偏移速度解释 叠前时间偏移的速度模型与叠加速度模型是不同的,需要通过对共成像点道 集（CIG）剩余时差的分析来完成修正，输入是通过叠前时间偏移得到的 CIG 道集，而非 ACP 道集。参考前述转换波叠加速度分析的说明，通过 Kirchhoff 叠前时间偏移得到一种反 NMO 共成像点道集，该道集只移除了非双曲线表现 的时差，保留了双曲线特征。对于该道集，应用双曲线速度分析，得到的速度 能向最优偏移速度收敛，因此能通过这种策略进行速度更新。 偏移速度分析中只有需要修正。 C V 启动模块：启动模块： 在命令行键入 psvmigpk，就可启动 Psvmigpk 转换波偏移速度解释分析工具。 弹出工区选择对话框，选择工区、测线；如图。选择工区后弹出当前工区 下的测线，选择测线弹出当前测下已经准备好的速度分析数据。选择已经准备 好的速度分析数据，选择初始参考速度库(PSV 库)。选择后点击 OK。系统弹出 速度分析系统交互主界面。如下图所示： 用户可以在左侧面板中调整转换波的偏移速度直至右侧面板中 ACP 道集的 时矩曲线与同相轴相吻合。 当拾取谱工具按钮按下时，开始扫描速度谱，各个工具按钮不能激活，当 扫描完毕时才自动激活。实时更新谱，实时解释。 更详细操作请参考PSVMIGPK 转换波偏移速度解释工具用户使用手册。 偏移作业提交以及监控管理偏移作业提交以及监控管理 在命令行中，键入：psvmig path，就会出现转换波并行叠前时间偏移作业 提交管理界面。其中 path 为要处理工区路径，例如/ht3/d/dpsprj/3d3c/。有关工 区路径、文件夹结构请参见DPS 工区管理使用说明。 界面说明 转换波并行叠前时间偏移作业提交管理界面 Run 提交叠前时间偏移任务，提交后监控该任务。 Monitor 监控当前工区路径下，指定的任务。 Quit 退出转换波并行叠前时间偏移作业提交管理界面 Job Setting 转换波并行叠前时间偏移作业参数卡编辑界面,如上图 1 所示， 主要涉及下列内容： Data for migration 数据选取，如红框中表格所示 Parameters for migration 偏移参数定制，定义偏移所需变量的输入值，如蓝 框中表格所示； Work Directory 临时工作路径设置，如黑框所示； Job 任务名设置，定义本次作业名称，如绿框所示； Save 保存按钮，保存以上四项内容到磁盘文件； Load 载入按钮，载入已保存的任务文件。 图 2 转换波并行叠前时间偏移作业节点选择编辑界面 Node Setting 转换波并行叠前时间偏移作业节点选择编辑界面 选择运行任务的计算节点。为了使用方便，右侧可以直接选择全部节点 （蓝色箭头所示），或者全部不选择（黄色箭头所示）。 转换波并行叠前时间偏移作业监控界面 Monitor 转换波并行叠前时间偏移作业监控界面 RunningStatus 运行状态。监控并行任务在每个节点上的运行状态，包括节 点名，PID 号，CPU 和内存占用情况。如红框所示。 进度条 监控任务运行进度。 View Log 查看任务日志，查看监控任务日志信息。 Kill Current Job 结束当前监控任务。 Clear Memory 结束当前所有的叠前时间偏移任务。 功能以及相关操作介绍 任务编辑 点击 Job Setting 表单选择转换波并行叠前时间偏移作业提交管理界面。 输入地震记录文件。具体操作如下： 添加添加 右键点击作业提交管理界面中红色框体区域，弹出选择表单，如下图所示。 根据 DPS 系统工区文件夹结构，选择 Browse from Database 将从数据库选 择文件方式将直接从 path 路径下 data/文件夹下开始选取地震数据，选择 Browse from others 将从用户根目录出开始选择数据文件，两者均支持多选。 选择数据后，数据将列出在表框中，继续执行上述操作可以追加数据文件。 图 4 选择输入地震记录文件表单。 移除移除 鼠标左键双击单个数据文件名称，可以从表单中移除该文件，如图 5 所示。 右键点击表单弹出选择表单，选择 Clear 可以移除所有选择的数据文件。 图 5 移除单个地震数据 定制偏移参数定制偏移参数 在作业提交管理界面中蓝色框体所示表格中： 指定速度文件。点击 VELFILE 右侧单元格中按钮，选择任务所需的速度文 件，或者在该单元格中文本框直接输入,操作如图 6 所示。 图 6 速度文件指定 指定输出文件。点击 OUTFILE 右侧单元格中按钮，指定任务输出文件的路 径和名称，或者在该单元格中文本框直接输入,操作如图 7 所示。 图 7 速度文件指定 指定偏移参数值。在各参数变量对应的 Value 栏中，填入对应的数值。各 参数意义和描述如下表所示 表 1 偏移参数变量说明 参数名类型取值范围功能备注 cdpminint Minimum cdp number of output section cdpmaxint Maximum cdp number of output section x0Int X coordinate of origin in 3D survey y0Int Y coordinate of origin in 3D survey x1Int X coordinate of shot line in 3D survey y1Int Y coordinate of shot line in 3D survey StackInt0,1Output CIP gather 0: CIP gather 1: setction invnmoint0,1Flag for inverse nmo 0:Forward 1:inverse cdplineInt The cdp line in the 3D survey linedxFloat The distance between cdp lines (m) cdpdxFloat The distance between adjacent cdp bins (m) cdpincInt The increment of cdp in output section seabedFloat The depth of seabed (m) TminFloat Minimum time of the migrated output TmaxFloat Maximum time of the migrated output maxoffsetFloat Maximum offset for CIP gather minoffsetFloat Minimum offset for CIP gather doffsetint Distance between adjacent trace in CIP tapertureFloat Times corresponding to aperture apertureFloat Tap the elipse if dip larger than it 指定临时工作目录指定临时工作目录 在 work Directory 文本框中输入或者点击右侧 Browse 按钮选择临时工作目 录。请在提示下选择。任务运行时将在该临时工作目录下建立以项目名和任务 名联合命名的临时文件夹，用于存放临时文件，在任务结束或者任务终止时删 除该临时文件夹。 选择临时工作目录提示 命名任务命名任务 在 Job 文本框中输入任务名称，用于保存任务和建立临时文件夹。 保存任务保存任务 点击 Save 按钮，将 Job Setting 表单中所有内容保存到 path 下 sjf/文件夹中， 以 job 名命名，格式是后缀为 pjf 的 ASIC II 文件。若该命名文件已存在，将收 到提示，选择“是”覆盖原 pjf 文件，“否”忽略保存要求。 保存任务 载入任务载入任务 点击 Load 按钮，弹出对话框，选择 pjf 文件载入的任务。 节点选择节点选择 点击 Node Setting 表单选择转换波并行叠前时间偏移节点选择界面。 转换波并行叠前时间偏移作业提交管理界面在运行时会自动载入默认节点 设置，通过勾选 Name 栏中框体选择该节点，并调节 CPUs 栏中数据选择使用 CPU 数目。 当载入默认节点设置失败时或者需要改变节点基本节点时，右键点击表单 弹出选项选择节点配置文件，载入节点信息，如下图所示 选择节点文件菜单 提交并