双狐-2007版变速成图手册[精华]

双狐-2007版变速成图手册[精华] 2007版 编写人: 舒永兵 保定双狐软件有限公司 二OO七年10月 前 言 双狐变速成图系统是在多年的塔里木变速成图攻关和生产过程中逐渐发展和完善起来的，塔里木盆地及其周边探区，构造运动剧烈，地表及地下地质情况都非常复杂，速度研究和变速成图难度大，双狐变速成图针对不同的地质情况集成了包括量板法、平均速度法、层位控制法和模型层析法在内的四种成图方法。 地震地质资料处理解释数据量大、难度大，速度 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 复杂，成图难度大，因此国际上最先进的设备、技术、数学理论都应用于此。双狐软件的主要开发者是在多年的油田科研生产中脱颖而出的若干青年技术专家，他们熟练掌握了国际最先进的勘探技术、计算机技术，并不断自发地开拓进取，攻破若干国际难题，创造了良好的国际影响，同时他们的软件产品也得到了国内外油田专家的认可。 保定双狐软件有限公司 地址:河北省保定市国家高新技术产业开发区火炬产业园五号楼 邮编:071051 主页: E-Mail:css@doublefox.com 技术服务热线:(0312)3180972 联系电话:(0312)3336700,09 传真:(0312)3183177 概述...................................................... 5 系统简介 .................................................................................................5 软件安装 .................................................................................................6 计算机配置要求.......................................................................................7 第一节 数据管理 ........................................... 8 (一)、数据转换 ...................................................................................8 (二)、数据应用 ................................................................................. 13 (三)、测量成果管理 .......................................................................... 16 第二节 常用工具 .......................................... 21 (一)、图形工具 ................................................................................. 21 (二)、数据处理工具 .......................................................................... 29 (三)、其他工具 ................................................................................. 33 (四)、数据库查询工具....................................................................... 38 第三节 速度分析 .......................................... 42 (一)、速度校正 ................................................................................. 42 (二)、速度场分析.............................................................................. 50 (三)、井速度分析.............................................................................. 59 (四)、其它工具 ................................................................................. 65 第四节 变速空校 .......................................... 70 (一)、时深转换 ................................................................................. 70 (二)、变速成图流程生成 ................................................................... 74 (三)、其它 ........................................................................................ 93 附录一:双狐速度研究软件 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 数据格式 ....................... 112 测量成果...........................................................................................112 解释层位T0数据(用于自动勾绘等值线) .......................................112 断层文件...........................................................................................113 速度谱文件 .......................................................................................114 速度谱与测量成果匹配后的叠加速度场文件(*.dfv) ........................115 用于平均速度校正的VSP测井数据 ...................................................115 做深度校正的钻井分层数据...............................................................116 概述 系统简介 地震速度研究几乎与地震方法本身同样久远，在数字处理技术出现之前的整个时代，地球物理学家利用单次覆盖资料，计算各个反射层的速度，仅是简单的估算法。自从有了多次覆盖数字处理技术以来，地球物理学家和地质学家仅能获得速度信息(平均速度、均方根速度、层速度为一空间数据体的速度场)，利用手工加计算器方式获得最终成果图。随着石油勘探开发程度的不断深入，所面对的是地质构造复杂、地震资料信噪比较低、速度谱资料也相应比较差的领域，利用常速充填获得的构造图已远远不能满足地震勘探开发领域的生产需要，因此速度研究显得越来越重要，尤其是山地的速度研究，速度的精度直接影响构造图的精度。建立精细可靠的速度场，获得精细准确的构造图，是我们进行速度研究的最终目的。 双狐变速成图系统是在塔里木多年的速度攻关研究过程中形成的一套成图系统，包括了量板法、平均速度法、层位控制法、模型层析法四种成图方法及各种速度研究应用模块。 对于没有速度谱资料的地区，采用量板法和附近井资料建立速度场，进行构造成图。 对于地层比较平缓的地区，采取平均速度法，不仅有利于小幅度构造解释，而且有利于提高构造图的精度。即在很大程度上不会因为速度场的问题导致小幅度构造的淹没，并能更准确的反映构造幅度和形态。 对于地层倾角较小，地层成层性较好的地区，可以采用层位控制法，对层速度进行平滑，然后建立速度场进行构造成图。 针对地层倾角大，构造幅度也比较大而且速度谱资料品质相对比较好的的地区，采用模型层析法，消除了地层大顷角的折射效应，准确定位构造高点，获取构造的真实形态。 系统集成了优秀的数学算法、突破了大数据量集成处理、解决了逆断层的一次成图问题、逆断层速度建模等关键技术，是一套来源于生产实践的应用软件，双狐变速成图系统2007版如图(1)所示，根据生产实际中的需要及应用要求，系统分为四个部分:数据管理、常用工具、速度分析和变速空校。 图1 双狐变速成图系统界面 软件安装 下面以在Windows xp中文版下安装双狐变速成图为例，介绍其安装方法，步骤如下: 1、在光驱中插入双狐变速成图安装光盘，点击setup.exe运行。如图2,3所示: 图2 图3 2、点击“下一步” 按钮后出现如图4所示的软件许可协议介绍，详细阅读后选择接受协议,点击“下一步”按钮，出现如图5所示的用户信息界面: 图4 图5 4、在用户信息界面中输入用户名、公司名称。点击“下一步”按钮，出现如图6所示界面。点击“更改”按钮可自行定制软件的安装路径，或直接点击“下一步”按钮，接受软件默认的安装路径。 5、点击“下一步”后出现如图7所示界面，点击 “安装”进行软件安装,然后点击完成即可。 图6 图7 计算机配置要求 操作系统:WindowsNT、Windows2000、WindowsMe、WindowsXP 硬件配置:内存512M以上，硬盘120G以上。 推荐系统: Windows2000、WindowsXP 第一节 数据管理 对于任何一个工区，要开展速度研究工作，基础数据是前提条件，基础数据的搜集整理将直接影响到最终的研究结果。基础数据包括研究区域内各目的层的反射层T0时间、工区测量成果数据、井位数据、钻井分层数据、VSP测井数据、叠加速度谱数据，声波测井数据、浮动基准面数据等相关数据。软件界面如图1-1所示，包括数据转换、数据应用及测量成果管理三部分。 图1-1 基础数据管理 (一)、数据转换 数据格式转换: 在数据转换下面点击“数据格式转换”，弹出如图1-2所示的对话框，在这里包含了规则数据的导入、测量成果数据的转换及各种速度普数据格式的转换，在鼠标选择一项后点击右键弹出“格式说明”菜单，点击格式说明就弹出对话框，详细的解释每种转换格式要转换的数据格式，这时用户根据自己要转的数据去选择数据格式。下面以规则数据导入和第一种速度谱格式转换为例进行说明: 规则数据导入:如T0工作站解释数据导入，数据是以空格分隔的5列规则数据。 gaosp 20517402.94 5315864.97 464.61 ts gaosp 20517419.63 5316044.47 460.06 ts gaosp 20517497.34 5315981.54 460.05 ts 工区名 X坐标 Y坐标 T0值 层位名 如图1-2所示，点击规则数据导入，右边出现输入文件对话框，点击”输入文件”，选择要导入的T0数据，如图所示， 导入起始行为1，如果文件头还有其它数据，根据行号修改导入起始行，如行号010才是需要导入的数据开始行，那么导入起始行为10，然后点击”下一步”，如图1,3，我们选择xyz散点，然后对应各数据行，如果坐标需要变换，如公里变成米，给定X,X*0.001,不需要变换默认1，然后确定，完成数据导入。 图1-2规则数据导入 图1-3规则数据导入 导入后的数据格式如下: layer M ts 20517402.94,5315864.97,464.61 20517419.63,5316044.47,460.06 20517497.34,5315981.54,460.05 以上数据格式为双狐网格化模块的标准输入数据格式。 同样如果导入的是测线或井位等其它数据，根据数据项选择对应关系即可完成。 速度谱格式转换: 输入数据格式如下: CONT BM92045 JD LINE LINE 1-320 AREA test INFO COM1 COM2 PNMO 0 0 0 1 SPNT 44301 382. 617329.2 328038.2 528 VELF 44301 100 1796 200 1796 300 1809 400 1821 500 1837 VELF 600 1867 700 1918 800 1970 900 2007 1000 2029 „„ VELF 5600 4321 5700 4349 5800 4367 5900 4381 6000 4392 SPNT 44350 431. 616662.1 329065.6 528VELF 44350 100 1797 200 1797 300 1810 400 1822 500 1837 VELF 600 1864 700 1913 800 1967 900 2010 1000 2038 首先分析要转换数据的格式，对比已经提供的数据格式要求，选择左边的格式，右边会列出 格式的输入数据样式，或鼠标左键点中左边“速度谱格式转换”，点右键，再点击“格式说明”， 会列出所有提供的数据格式。比如以上数据与格式一一致，注意:文件头可以忽略，或者在记事 本中删除。选择输入输出文件，点击“转换”完成。 图1-4 速度谱格式转换 界面中“x=x*”表示在转换过程中将x坐标乘以一个系数，如乘1000可以将公里单位变成米，其它类似，加坐标带号是在坐标的前面加上坐标带号，交换xy坐标是将x、y坐标交换一下位置，如通常在野外的x坐标在室内要当y坐标，y坐标要当x坐标，这时就要交换位置。 创建文件列表:是将所选文件的路径及文件名存成文件名列表。 图1-5 创建文件列表 其它数据格式转换:使用下面的模块，使用方法与前面格式转换模块类似。 图1-6 数据格式转换 去处断层范围内的点:是将曲线范围内的散点数据去掉。 图1-7 图1-8 断层范围内点去除前后 曲线处理:实现计算交点、闭合曲线变单线、区分断层上下盘的功能。 图1-9 测线交点计算 (二)、数据应用 按线名提取速度谱: 速度谱数据由于不同时期处理的数据差异较大，因此对速度谱数据的整理要求更高。 提取速度谱(如图1-10所示):是按测线名在速度库中提取对应测线的原始速度谱，点击“速度谱”按钮选择库文件路径，测线名可以输入也可以打开测线名文件来获得，选择需要的速度谱，然后存到指定目录。 图1-10 提取速度谱 速度谱与测量成果匹配:是完成速度谱与测量成果之间的坐标匹配(如图1-11所示)，因为测量成果中有坐标和桩号，而速度谱中只有桩号，因此只有通过匹配来获得速度谱中的桩号位置处的坐标。测量成果的格式为最终格式(见双狐速度研究软件标准数据格式)，而速度谱的格式为单条速度谱文件格式(见双狐速度研究软件标准数据格式)。应用时首先点击“测量成果”按钮，选中测量成果文件，然后点击“速度谱”按钮，选中速度谱文件，可以一次选多个文件，选中的速度谱文件将把文件名列在下面的列表中，如果有错，可以选中错误项，然后用“清除选项”来删除，如果没错就可以点击“将列表中文件名存入文件”按钮将速度谱文件的文件名存到一个文件中，如果下次应用，就从“打开文件名列表”中选择，就不用重新创建。匹配时首先看匹配关系，有哪些线能匹配上，哪些不能匹配，缺少哪些线，看时可以以线方式和点方式两种方式来看，如果不能完全匹配，可点击“查看测量成果”或“查看速度谱”来查看到底是测量成果没有，还是桩号不匹配等原因引起的。同时如果发现有错还可以对其进行修改。如果不匹配，其中原因主要有两个方面，一是测线名与速度谱名字不一样，二是测线的桩号与速度谱的桩号不对应。 图1-11 速度谱与测量成果匹配 按范围输出速度场:从已知叠加速度场中提取指定范围内的速度场数据(如图1-12所示)， 可以按坐标范围、给定点坐标和给定多边形范围提取。 图1-12 提取速度场 Vsp数据提取: 当Vsp井比较多而且又放在同一个文件夹时，可以通过查询来寻找需要的数据并输出到一个指定的目录(如图1-13所示)。 图1-13 提取vsp数据 在提取数据时，根据指定的基准面信息对输出VSP进行基准面统一校正，基准面可以是变化的浮动基准面或统一基准面。 (三)、测量成果管理 测量成果文件管理:测量成果数据的管理可以用数据库方式进行管理，也可以用文件方式进行管理(如图1-14所示)，用文件方式管理时需要同一坐标带号，将整理好的(双狐格式)测量成果存入一个库文件，应用时在这个库文件中按范围进行查询或按测线名查询(如图1-15所示)。 图1-14 追加测量成果 图1-15 测量成果文件查询 给线编CDP号:将线类数据按步长编简易CDP号。 图1-16 编CDP号 桩号转坐标:将桩号剖面等桩号数据按测量成果中的桩号和坐标关系转成平面坐标数据。 图1-17 桩号转坐标 从测量成果中提取高程:是从测量成果中提取高程数据(如图1-18所示)，测量成果文件格式为双狐测量成果格式 图1-18 提取高程 生成三维测网:根据两个顶点、原点方位角、任意三点等参数生成三维测网。 图1-19 生成三维测网 三维测网批量生成: 生成三维测网(如图1-20所示)是利用三维测网的四点坐标批量生成三维测网。输入文件为线 号与坐标对于关系文件。例如: YM21 14 1 1 558197.00 4539113.00(工区名 带号 线号 道号 X坐标 Y坐标) YM21 14 1 704 547563.00 4553106.00 YM21 14 860 704 564661.00 4566100.00 YM21 14 860 1 575295.00 4552107.00 YM17 14 1 1 573960.00 4551093.00 YM17 14 1 624 564534.00 4563492.00 YM17 14 861 624 581651.00 4576503.00 图1-20 测网生成 第二节 常用工具 软件的“常用工具” 部分界面如图2-1所示，包括“图形工具”、“数据处理工具”、“数据库查询工具”和“其他工具”四个部分。 图2-1 常用工具界面 (一)、图形工具 图形工具中包括包括图形编辑、双狐图件加载到GoogleEarth、逆断层数据分离和二维、三维网格化。图形编辑的使用在坐标工具使用手册中介绍很详细，这里不再介绍，下面介绍其余几个模块。 图件转GoogleEarth格式:首先用dfdraw打开图件，指定投影参数，然后保存文件。 图2-2 指定图形投影参数 如图2,3所示，输入双狐图形文件，然后确定完成转换。转换后将生成文件加载到GoogleEarth，启动GoogleEarth，然后在文件中打开转换后的文件，自动定位到指定位置，如图2,4所示。然后在地标框中分类列出了所有信息，并可对元素属性进行修改。 图2-3 转GoogleEarth格式 图2-4 加载到GoogleEarth 逆断层数据分离:(如图2-5): 图2-5 数据分离 分离前将解释数据和断层数据合并到一个文件中，断层类别为layer m faultage，逆断层成对出现，将上下盘数据命名为正负，如下图所示，2号断层上盘pline.2，下盘pline.-2，解释数据下盘的要与上盘断层相交但不能与下盘断层相交，同样上盘解释数据要与下盘断层相交但不能与上盘断层相交，解释数据格式如图所示。分离后数据如右边所示(如图2-6)。 图2-6 逆断层数据示意图 平面网格化: 网格化主要是由散点数据形成空间曲面模型，也就是自动勾绘等值线图，其软件模块为dfCreatef(x,y).exe模块，界面如图2-7所示，这是一个通用的工具模块，在研究过程中使用非常多，而且非常重要，用好这一模块将对速度研究起到决定性的作用。 1(点击 “输入文件”获取散点文件，同时在输入断层、边界、及输出文件。 2(网格参数:,由上一步读取文件后自动找出该散点文件的有效坐标范围，同时计算出常规情况下的‘网格数目’及‘网格间隔’，在此，用户可以根据工区的实际地质情况及勘探要求的精度等情况手工给定一个范围，同时可以在这个范围的基础上放大或缩小，人工干涉‘网格数’、‘网格间隔’ 及‘网格加密次数’数值大小。调配的原则是:在‘网格加密次数’值不变的情况下，‘网格间隔’值越大，‘网格数’越小，确定的T0曲面就越粗，相反条件下，确定的T0曲面就越精细。另外，在不改变‘网格数’、‘网格间隔’的情况下，通过调整‘网格加密次数’值的大小来提高曲面精度，当加密次数为1(表示2的1次方)时，总的网格数为81*2*2*2x68*2*2*2， 当数值增大时网格数增大，曲面精度提高，反之，曲面精度差。通常情况下，如果断层较密，那么初始的“网格数目”给大一些，如果断层少或没有，那么给小一些。 “精度”表示生成曲面的迭代过程中迭代曲面与散点之间的最大允许误差，小于这个值时迭代终止。也可以通过减小它来提高精度。 “平滑”可以根据资料品质的好坏情况及应用要求确定合适的平滑度。数值越大，表示平滑作用越大，数值越小平滑作用越小，数值为0表示没有平滑，完全过数据点。 填好参数后，点击界面工具栏“～“符号，执行网格计算。 图2-7 平面网格化 注意:在执行过程中，会弹出图2-8所示的迭代监控界面，图中的“8.200137e-002”表示收敛精度，如果越来越小，表示迭代收敛，如果越来越大，如变成“8.200137e+002”，表示迭代不收敛，应立即停止，调整网格参数重新执行，这时在显示数据的地方点一下，弹出图2-9，点击结束，这样就终止了计算回到参数控制界面。 图2-8 迭代参数控制 图2-9 3(生成等值线: 图 2-10 等值线生成 图2-10所示，首先确定输出文件的文件名，改变等值线在图上显示的疏密程度。当“步长”和“隔几根上数”数值为20和5时表示等值线间隔为20，每5根等值线一加粗上数。在确定完相应的参数值后，点击“～”按钮，进入DfDraw功能模块(如图2-11所示)，由它实现数据的可视化功能及图形编辑输出，如等值线图的编辑及各种修饰等，是一套最完备的绿色环保型工业制图软件。 图 2-11 DfDraw图形编辑界面 三维网格化: 三维空间网格化的界面如图2-12。 图 2-12 三维空间网格化 做三维空间网格化要求输入四列数据，数据格式如下:在网格化维数处选择“3D:u=f(x,y,z)”。 14778414.268100,4321348.622300,2682.680000,2862.140000 14778414.268100,4321348.622300,3792.680000,3492.350000 14778414.268100,4321348.622300,4712.680000,4319.980000 14778414.268100,4321348.622300,4922.680000,4530.570000 14778414.268100,4321348.622300,5922.680000,5029.060000 14778414.268100,4321348.622300,6922.680000,6030.940000 14778199.160500,4320371.949100,82.470000,1255.000000 在图2-13中对空间网格化的参数进行设置，begin p、end p和delt p分别是x、y、z的最大最小范围，number分别是初始网格数目，网格加密次数是在做网格化时对网格的加密次数，如加密次数等于3，那么x方向的最终网格数目,21×2×2×2(即2的3次方)，同样y和z方向也是一样。 注意:在给参数时不能给的太大，因为是三维网格化，网格数太大生成的文件就很大，如果初始网格都给30，在加密2次，那么最终形成的网格节点数就,30×4×30×4×30×4,1728000，生成的文件数据量达到21.6M，如果初始网格数给50，加密次数给3，那么形成的网格节点数为64000000，换算一下形成的文件数据量为800M，这个数据量在目前微机上处理比较慢。因此在应用时要注意参数不能给太多。 图 2-13 参数设置 在设置可信度时由于是空间上的平滑，平滑系数不能给太多，速度平滑通常给10到20左右，根据资料情况来定。 图2-14 可信度设置 要看网格化结果，可以通过DFVf(xyz)OutPut.exe模块来看，具体使用在后面有详细介绍。 图形修饰:将网格化结果对应类别和修饰，自动成图。在网格化结果文件中合并断层、井位、测 网等基础数据，然后保存文件，在等值线文件中打开保存的文件，选择修饰文件，然后在类别对应中选择等值线文件中的类别对应修饰文件中的类别，然后添加对应关系。然后根据需要选择图框、图例、责任表。确定即可。 5万(修饰符号).dfd 图2-15 图形修饰 图2-16 结果对比 (二)、数据处理工具 坐标数据运算:(如图2-17所示)是对xyz的数值进行加减乘除等数学运算。 图2-17 坐标数据运算 曲面散点运算f(x,y):(如图2-18所示)是完成散点数据与双狐曲面文件之间的各种数学运算。 图2-18 曲面散点运算 曲面散点运算dfg: 是完成散点数据与双狐曲面文件之间的各种数学运算,输入曲面文件为dfSurface网格化生成的结果。 曲面与多个散点同时运算f(x,y): 是完成多个散点数据文件与双狐曲面文件之间的加、减、乘、除运算。例如将所有T0层位从统一面减到cmp面，只需要将cmp面时间形成曲面，各T0层位为散点，然后散点减曲面即可。 图2-19 曲面散点运算 曲面文件输成散点:(如图2-20所示)是将双狐二进制曲面文件以散点(x,y,z)方式输出，输出 时可以通过修改每几个点输出一个来减少数据点，同事可以按X、Y方向线类输出。 图2-20 曲面文件输出散点 坐标投影换算:(如图2-21所示)是对xyz文件进行坐标投影换算(使用详见，《坐标工具》)。 图2-21 坐标投影换算 有时启动该模块时需要选择参数文件，如图2,22，参数文件在安装的bin目录下面。 图2-22 坐标换带 数据筛点:(如图2-23所示)对数据量过大而又没有必要的文件进行数据点筛稀，便于应用。 图2-23 数据筛点 等值线输出xyz文件:(如图2,24所示)将双狐等值线文件输出散点(x,y,z)。 图2-24 输出xyz文件 逆断层曲面插值:(如图2-25所示)断层从曲面进行插值，插值出有Z值的断层。 图2-25 逆断层插值 (三)、其他工具 数据排序:(如图2-26所示)将散点数据点按指定的列进行排序。 图2-26 数据点排序 T0层位划分小层:(如图2-27所示)由已知层位的曲面数据插值生成小层数据。 图2-27 T0层位划分小层 时间曲面计算倾角体:(如图2-28所示)由已知的各个层位，通过逐层求X、Y方向偏导数的方式计算倾角，然后形成导数体，顶底延伸是对在z,0时倾角为0，最下一层向下延伸计算倾角，去掉边界效应。 图2-28 时间曲面计算倾角体 由参数生成等时间层位:(如图2-29所示)通过给定坐标范围和网格参数及时间间隔，自动生成等间隔的等时间层位。 图2-29 等时间层位生成 图2-30 参数文件格式 从分层文件中提取数据:(如图2-31所示)从解释中加载的分层数据文件中提取某一层的钻井分层数据。数据格式看帮助。 图2-31 曲线插值计算:(如图2-32所示)输入曲线，选择插值方式和节点个数，可以得出不同插值方式的曲线。可以看生成试验数据。 指定点插值:指定x坐标，插值该点的y值; 所有线拟合一条:原始曲线有多条线时，拟合成一条; 图2-32 曲线插值计算 平面上去异常点:(如图2-33所示)根据Z值与坐标点给定半径范围内的Z值的平均值的比例来去留数据点，如果大于给的范围(如2,)就认为是奇异点，在平面上去掉。 图2-33平面上去异常点 文件曲线分离:(如图2-34所示)把文件中线名相同的曲线分离成一个单独的文件。 图2-34文件曲线分离 生成斜坐标网:(如图2-35所示)根据原点、角度、步长生成斜坐标网。 图2-35 斜坐标网生成 (四)、数据库查询工具 数据库查询工具:界面如图2-36所示。 图2-36 数据库查询 例如从Access库表中查询井坐标，形成井位图流程如下: 图2-37 数据库查询流程 第三节 速度分析 速度分析包括速度校正、速度场分析、井速度分析及其它工具四部分。(如图3-1所示) 图3-1 速度分析 (一)、速度校正 浅层建模后时间速度校正:由于浅层速度模型存在误差，尤其在山地地表变化复杂的地区，误差往往比较大，因此造成静校正存在误差，使得叠加速度和T0时间也存在误差，如果通过重新建模，能获得比较准确的浅层模型和静校正量，那么在不对资料进行重新处理的情况下，如何用新老校正量和原有速度谱进行时间速度的校正,本模块包含需要的数据格式转换和二维校正、三维校正模块，校正采用流程方式执行。 图3-2 时间速度校正 图3-3 数据格式转换 图3-4 二维校正流程生成 图3-5 二维校正执行流程 速度空白段充填:由于速度谱的第一个谱点时间不均匀，而且通常时间比较大，在100,300毫秒之间，这样在速度建场时，第一点以上的速度只能用第一个谱点的叠加速度充填，而第一个谱点速度通常误差也比较大，因此会对变速成图带来误差，如果通过速度反演，能获得比较准确的浅层速度模型，因此就用反演的速度来充填速度谱中的速度，可以减小速度场的误差。 图3-6 反演速度充填 根据T0闭合差修正速度谱校正量: (如图3-7所示)在二维剖面解释过程中，由于时间剖面存在闭合差，而时间剖面的时间和速度谱的时间是对应关系，因此速度谱中的时间也要消除时间闭合差，通常是应用单条线的时间和T0时间平面图计算每条线时间和T0的时间差，然后用这个时间差作为各自速度谱的时间校正量，对速度谱进行校正，使速度谱的时间和T0图对应。 图3-7 根据T0闭合差修正速度谱校正量 生成速度闭合差校正流程:(如图3-8所示)生成叠加速度闭合差校正及转为均方跟速度的流程，根据测量成果把原始叠加速度分为两个方向，然后沿测线计算时间倾角或从segy地震剖面自动追踪计算时间倾角，如果测线比较稀，可以在Dfdraw中复制测线，加密求倾角的控制点。 图3-8 叠加速度计算均方根速度流程生成 图3-9 叠加速度计算均方根速度执行流程 速度闭合差校正:(如图3-10所示)由于二维速度谱在主测线和联络测线之间存在速度闭合差，这个闭合差尤其在山地资料中非常突出，会对变速成图造成较大的误差，因此要对速度谱的闭合差进行校正。校正时首先根据测线名把速度谱分为两个方向，同时输入两个方向的T0时间倾角，根据时间倾角计算各自的加权系数，然后进行加权校正。具体校正方式有三种可以看帮助。 图3-10 速度谱闭合差校正 速度谱统一面与浮动面校正:是将基于统一基准面的速度场校正到以浮动基准面为准，或将基于浮动面的校正到统一面(选“反向校正”)(如图3-11所示)。 图3-11 速度谱统一面与浮动面校正 速度谱基准面校正:是根据校正量将速度谱校正到基准面，速度谱的起算面通常是cmp面，无论是统一面成图还是浮动基准面成图，速度谱的基准面与T0时间的基准面之间都存在差异，这个差异用时间来表示就是速度谱的校正量，因此成图时要用校正量对速度谱进行校正。使速度谱的基准面和T0时间的基准面一致。 -12输入文件为文件名列表，然后输入地表充填速度，就可以将速度谱校正到基准面(如图3所示)。选输入格式看数据格式。 图3-12 速度谱基准面校正 层速度替换校正:在某一层的层速度已知的情况下，可以通过层位顶底时间为界，将这一套地层的层速度充填到速度谱中去，然后重新计算速度谱，可以通过给点坐标的方式指定那些谱点要替换，其它谱点不变，或给定多边形范围，范围内的替换，范围外的不变(如图3-13所示)。 图3-13 层速度替换 vsp校正层速度:(如图3-14所示)用vsp时间速度曲线来校正某一层的层速度。 图3-14 vsp校正层速度 T0漂移问题。可以看原理简介。 (二)、速度场分析 速度谱编辑:(如图3-15所示)是检查速度谱是否有人为因素造成的错误，如果有将其改正，然后转换替换原文件。 图3-15 速度谱编辑 速度场剖面分析:(如图3-16所示)该模块的功能是实现叠加速度场的沿测线的剖面分析。输入文件为叠加速度场文件，在“输出测线索引文件”处给定输出文件名，然后确定，这时输出速度场的测线，如图3-17所示，从“选择”菜单中点击“选项”，然后选中“移动选择”，这时鼠标移动到一条测线上，点击鼠标左键，弹出如图3-18所示的单条速度剖面分析模块，这时可以对这条速度剖面进行各种方式的分析。可以以等值线方式或曲线方式进行叠加速度、层速度及平均速度等分析，当以曲面或等值线方式进行分析时，点击“形成曲面”按钮，然后对输出文件进行网格化，如果要编辑等值线，将编辑后的等值线曲线名称变为z值，然后保存净曲线，用dfcreatef(x,y).exe进行网格化，保存曲面文件，然后点击“输出速度剖面曲线”进行转换。以曲线方式进行分析时首先给定速度及时间的显示范围和网格步长，然后点击“显示or编辑”按钮，以线方式编辑，然后输出编辑结果替换原来的文件。 图3-16 速度场剖面分析 图3-17 速度场测线分析 图3-18 单剖面速度分析 在图3-16中，点击“获取任意线”，显示测网，在测网中画任意线，然后保存出任意线，这样就可以确定，然后对任意线剖面进行分析。 速度剖面平面分析:(如图3-19所示)在文件中打开dfv格式速度谱，点击任意一条线，上面窗口显示叠加速度，下面窗口显示层速度，可以交换编辑，还可以形成等值线，编辑等值线，然后修改速度谱，同样可以形成速度切片，进行编辑修改和数据统计等。 图3-19 速度剖面平面分析 叠加速度场曲线拟合:(如图3-20所示)根据叠加速度谱的变化趋势，拟合一条速度曲线，可以拟合全区的速度谱，也可以给定坐标点，按半径范围拟合，或者按给的点周围4点速度谱拟合一个曲线，输出时可以叠加速度、层速度和平均速度。 图3-20 叠加速度场曲线拟合 叠加速度谱倾角校正:(如图3-21所示)对叠加速谱做倾角校正，T0-dx..f(x,y)、T0-dy..f(x,y) 分别是各层T0曲面沿X、Y方向计算的时间偏导数，然后进行三维网格化的结果。 图3-21 叠加速度倾角校正 空间速度场转换:(如图3-22所示)用叠加速度进行各种速度转换和提取。如果提取两层之间的层速度，需要输入层位顶底界面的时间曲面，如果沿层提取速度，需要输入时间层位曲面(曲面都是DFCreatef(x,y).exe生成)，提取切片，需要输入切片时间。 图3-22 空间速度场转换 空间速度场分析:(如图3-23所示)是针对三维立体网格化结果，对数据体提取切片、按剖面输出、沿层位输出、按任意坐标点输出以及作插值等操作。 图3-23 速度场分析 空间速度场校正:(如图3-24所示)由于空间速度场与vsp的速度通常存在误差，为了保证速度场的精度，还可以使用vsp速度对速度场进行校正，校正方式有体校正和逐层校正，体校正是所有误差点都参与每个网格节点的校正。逐层校正是算每个时间点的误差，然后由这些误差在这个时间平面做平面校正。 图3-24 空间速度场校正 空间速度场转segy:(如图3-25所示)把三维网格化的结果数据转换成segy数据体。 图3-25 空间速度场转segy 由时间层位和层速度建立速度场:(如图3-26所示)根据解释的时间模型和层速度建立速度场，层速度可以是固定速度或平面层速度变化曲面，时间模型可以考虑逆断层，如果有逆断层，需要输入断面。 图3-26 由时间层位和层速度建立速度场 叠加速度场中批量提取沿层速度:(如图3-27所示)从叠加速度场中批量提取沿层速度。 图3-27 叠加速度场中批量提取沿层速度 图3-28 叠加速度场中批量提取沿层速度流程 (三)、井速度分析 井斜校正:(如图3-29所示)对斜井沿井轨迹计算各点的坐标及垂直深度。 图3-29 井斜校正 vsp测井速度分析:(如图3-30所示)是对vsp测井数据进行速度分析，以曲线方式输出各种速度。 图3-30 vsp测井速度分析 声波速度曲线坐标求取转换:(如图3-31所示)由于做速度场校正的时间速度曲线在校正速度场时必须有校正井点的坐标，本模块根据文件名和井名自动寻找井的坐标。 图3-31 vsp曲线坐标求取 vsp速度基准面校正:(如图3-32所示)将vsp速度从井口补心高校正到基准面。 图3-32 vsp基准面校正 Vsp中提取层速度:(如图3-33所示)从Vsp数据中计算任意两点的层速度。 图3-33 vsp中提取层速度 Vsp中提取层两深度点之间层速度:(如图3-34所示)从Vsp数据中计算任意两点的层速度。 -34 vsp中提取层速度 图3 提取声波时差数据:(如图3-35所示)从声波测井文件中提取用于速度分析的深度和时差数据。 图3-35 提取声波时差数据 测井曲线拼接:(如图3-36所示)将同一时期测的多段测井曲线进行拼接。 图3-36 曲线拼接 计算声波初始时间:(如图3-37所示)由于声波时差影响因素比较多，往往在初始一段误差 较大，在计算声波速度时，初始时间存在较大误差，因此可以由速度场来计算声波曲线的初始时深关系。如果应用声波测井做速度分析，通常速度场是基于基准面的，而声波如果是基于井口补心高的，需要做一个基准面校正，将声波数据从井口补心高校正到基准面，这样方便对比分析。 图3-37 计算声波初始时间 声波测井速度分析:(如图3-38所示)对声波测井数据进行速度分析，以曲线方式输出各种速度，由于输入的测井数据单位有的用米，有的用英尺，因此在应用前要确定单位。由于声波测井受各种因素的影响，误差较大，但是声波测井能反映细小的速度变化，而vsp测井比较稳定，精度较高，但测量间距大，不能反映速度的变化细微变化，因此可以用vsp对声波测井进行分段校正(如图3-39所示)。 图3-38 声波测井速度分析 图3-39 声波时差校正 时深曲线校正:(图3-40)由于声波计算的时间速度曲线因为初始时深关系误差较大，如果已知某层的T0时间，又有这层的准确分层，因此可以计算这个分层深度的时深关系，用这个时深关系对声波的时深关系做一个校正。 图3-40 时深曲线校正 (四)、其它工具 叠加速度场坐标运算:(如图3-41所示)对叠加速度场的坐标进行加减乘除运算。 图3-41 叠加速度场坐标运算 V0和β拟合:(如图3-42所示)由时间深度或时间均方跟速度拟合V0和β。 图3-42 V0和β拟合 根据测线名分离速度谱:(如图3-43所示)根据测线名将速度谱分为两个方向的文件。 图3-43 分离速度谱 根据点坐标去除速度异常点:(如图3-44所示)对于速度谱中的异常点，可以通过输入点的坐标进行去除。 图3-44根据点坐标去除速度异常点 叠加速度场横向加密:在已知叠加速度场中，给定坐标，生成这些点的叠加速度场，对原速度场进行加密(如图3-45所示)。 图3-45 叠加速度场加密 叠加速度场纵向加密:(如图3-46所示)叠加速度谱的时间采样点通常都比较少，有的几百毫秒一个点，有时层位比较薄，只有几十毫秒，因此在提取速度时存在一定误差，因此可以先对叠加速度谱点时间点进行加密，然后去计算。 图3-46 叠加速度场加密 由时间和速度重新生成速度谱:(如图3-47所示)通过已知的层速度(速度曲面或固定速度)和 已知的时间层位，通过插值计算和dix公式重新生成叠加速度场。 图3-47 由时间和速度重新生成速度谱 生成可信速度谱充填不可信速度流程:(图 3-48所示)根据剖面的质量确定那些速度谱是可信的，那些是不可信的，然后用速度谱的平面坐标把不可信速度谱的平面坐标保存出来，用“根据点坐标去除速度谱异常点模块”把速度场分成2个文件，分别是可信叠加速度场和不可信叠加速度场。 用可信度较高资料品质较好的速度充填品质较差速度不可信的区域，输入T0从CMP面起算。具体的流程生成可参考前面的流程生成模式和数据格式。 图 3-48 生成可信速度谱充填不可信速度流程 第四节 变速空校 变速空校包括时深转换、平均速度法变速、层位控制法变速、模型层析法变速空校及其它相关操作(如图4-1所示)。 图4-1 变速空校 (一)、时深转换 时间剖面转深度(层速度充填):(如图4-2所示)通过各层层速度将时间剖面直接转换成深度，还可以形成平均速度剖面，然后由平均速度剖面和时间转深度。格式看帮助。 图4-2 时间剖面转深度(层速度充填) 时间剖面转深度(平均速度剖面转):(如图4-3所示)通过平均速度剖面将二维时间模型直接转换成深度。格式说明看帮助。 图4-3 时间剖面转深度(平均速度剖面转) 时间剖面转深度(Vsp平均速度曲线):(如图4-4所示)通过Vsp平均速度曲线将时间剖面直接转换成深度，格式说明可以看帮助。 图4-4时间剖面转深度(Vsp平均速度曲线) 时间剖面转深度(三维平均速度体):(如图4-5所示)通过三维平均速度体将时间剖面转换成深度。格式说明可以看帮助。注意时间模型的文件名除了扩展名以外，要与测线中的线名完全一致，否则不能匹配。 图4-5 时间剖面转深度(三维平均速度体) 由时间segy剖面转深度segy剖面:(如图4-6所示)有单条线的测量成果和单条时间剖面(.section文件)和单条速度剖面转换成深度剖面。 图4-6 由时间segy剖面转深度segy剖面 将时间segy应用双狐解释系统，建立工区，生成inline目录下的*.Section文件。 将速度谱进行单剖面速度分析，形成平均速度剖面，注意速度谱的cdp号与segy剖面一致。 曲线关系转换:(如图4-7所示)基于曲线XY关系，已知X计算X对应的Y，比如输入时深曲线和时间散点计算深度，或输入时间速度曲线和时间由“(x,y,z*Z*0.0005)”计算深度等。 图4-7 vsp时深转换 公式转换:(如图4-8所示)当时深关系可以用公式来表示时，那么就可以通过公式来进行时深转换。 图4-8 公式时深转换 生成层速度曲线模板:(如图4-9所示)根据已有井的时间层速度生成一个时间层速度坐标模板。 图4-9 生成层速度曲线模板 (二)、变速成图流程生成 成图开始流程:要对一个工区进行成图时，可以生成该流程，集成了数据格式转换、整理等模块 和成图方法流程的生成。 图4-10 开始流程生成 量板法成图: 针对资料品质相对较差以及速度谱可信度相对较低地区运用井的层速度充填各层T0模型变速成图。 在应用时将基础数据准备好，然后就可以通过图4-11的界面生成批处理流程，应用人员按照流程执行就可，只是在中间网格化时需要控制参数。 图4-11 量板变速流程生成 量板法思路:在资料品质相对较差速度谱可信度相对较低地区，运用井的层速度充填各层T0建立模型并转换成平均速度场，得出最后各层构造图。 下面来介绍量板法数据准备和变速流程: 数据准备:各层解释T0和断层，工区内各井的层速度。 流程生成:(如图4-11所示)创建断层文件名列表(如图4-12所示)和网格参数以及输出流程的文件路径。在创建列表过程中选择对应的T0散点文件和对应的断层文件，然后按上下关系依次添加，保存列表。网格参数可以用dfCreatef(x,y).exe网格化一个范围比较大的T0文件，然后在文件\另存为下给文件名，扩展名用“ini”,保持网格化参数，再通过读参数文件来获取参数。 图4-12 创建断层文件名列表 量板法变速流程:(如图4-13所示) 图4-13量板法变速流程 第一步:网格化控制层: 网格化各层控制层，不需要计算倾角，T0划分的小层数1。由于模型是层速度和解释T0建立模型，所以不需要计算倾角和划分T0。 第二步:逆断层断面生成:由解释的T0曲面插值出断面。详细步骤流程上有说明，注意流程上说明要注意的问题。 第三步:生成控制层位量板:在做这一步之前可以看看第四步流程执行程序的帮助文件。 -14所示)由第三步的列表文件自动建立平均速度场。第四步:生成平均速度场:(如图4 图4-14 生成平均速度场 第五步:平均速度三维网格化: 第六步:应用:对任意给定的t0时间层位，要做变速，就可以用时间层位到上面形成的速度体中计算深度进行成图。 图4-15 任意层位成图 平均速度法成图: 平均速度法变速空校通常有两种方式成图，一种是由叠加剖面解释成图，一种是由偏移剖面成图。 在应用时将基础数据准备好，然后就可以通过图4-16的界面生成批处理流程，应用人员按照流程执行就可，只是在中间网格化时需要控制参数。 图4-16 平均速度法变速流程生成 如果没有断层和边界，文件名为空，如果基准面为统一基准面(如:1300m)，那么在文件名框中直接输入基准面的值1300，如果没有井分层做井校，文件名为空。 由图4-16所示的文件参数，生成的叠加剖面成图操作流程如图4-17所示:执行思路是应用速度场和叠加t0图先做时深转换，然后对深度进行偏移归位，再做基准面校正及井校。 图4-17 叠加剖面成图流程 为了应用方便，在执行之前需要将各种数据的单位统一成米(m)单位，通过看数据确认数据的单位，如果有公里单位的那么可以通过单位换算工具(dfplus.exe)或借助dfdraw的校位功能及其它方式将坐标换算成米(m)单位，然后用生成的文件替换原文件，并将文件名与源文件名一致。 生成的偏移剖面成图操作流程如图4-18所示:执行思路是应用叠加平均速度场和偏移t0图先做初始时深转换，得到一个初始深度，然后用初始深度倾角对时间反偏移，用反偏移时间再求速度，根据时间偏移量和求出的速度计算偏移速度，用偏移速度和偏移时间做时深转换，然后做基准面校正及井校。 图 4-18 偏移剖面成图流程 下面以偏移剖面成图流程为例来详细解释中间各步骤的具体操作及生成的过程文件。带“**”号的为必须执行的步骤，其余可执行可不执行。 步骤详解: 第1,8步:为数据准备和检查。首先将T0数据转换成xyz散点或带数值的等值线，速度场转换成双狐dfv格式，断层可以在dfdraw中显示，边界在dfdraw中形成闭合，井分层数据准备 好，所有坐标数据都统一成米(m)单位。数据格式如下: T0文件数据格式(正断层): Layer M 0 Layer ViewAndEdit 13510283.01,4560064.35,2760 13510293.97,4560086.82,2764 13510304.93,4560109.29,2769.6 13510315.89,4560131.75,2775.2 。。。。。。 T0文件数据格式(带逆断层，由dfSurfaceXP.exe分离得到): Layer M 0 Layer ViewAndEdit Pline.bc250\-1 (bc250为线名，\-1断层好) 13510283.01,4560064.35,2760 13510293.97,4560086.82,2764 13510304.93,4560109.29,2769.6 13510315.89,4560131.75,2775.2 。。。。。。 速度场dfv格式: STATION 9240-VEL 122.0000 13703649.0000 4364467.0000 0.00 0.00 20 285.00,1928.42 585.00,1958.78 885.00,2100.93 。。。。。。 断层文件格式: Pline.2 13703209.5,4376919.7 13707790.4,4374629.2 13712589.4,4372338.8 Pline.3 13691975.4,4371684.4 。。。。。。 边界文件格式: Layer M range Layer ViewAndEdit Pline 13677687.4,4393607.2 13757089.5,4373211.3 13751308.9,4347362 13662635.9,4370702.7 13669616.3,4395461.4 13677687.4,4393607.2 c 钻井分层文件格式: 图 4-19 井分层数据格式 h0为井口高程，如果分层已从海拔0起算，h0,0，dh为目的层的埋深。 第9步:调用下面的模块将叠加速度体转换成平均速度体，输出文件为vs-vav.xytv。 图 4-20 速度转换 第10步:调用dfDimension.exe.模块对vs-vav.xytv进行体网格华，选3d，输出vs-vav.f(x,y)文件。在填写参数时，第3行是时间的范围，通常根据作图的深度初步给个范围，大于最大时间500毫秒足够了(如最大时间4500毫秒，我们给5000)，最小通常给0，可以减少数据量，第1行是x坐标，第2行是y坐标，根据坐标范围给定deltp的值，使初始网格数目number合适，通常10