浅地层剖面仪SEGY型数据文件的读取

浅地层剖面仪SEGY型数据文件的读取 () 文章编号 :10012909 X 20090420063209 浅地层剖面与测深数据综合成图法 在舟山大陆连岛工程中的应用 丁维凤 ,刘 奎 ,韩富江 ,陆晓鸣 ,钱万民 ( )国家海洋局 第二海洋研究所 ,浙江 杭州 310012 () 摘 要 :利用 2006 年舟山市大陆连岛工程 三期岱山跨海大桥工程物探调查 , 、讨论了海底 水深剖面图的生成 、浅地层剖面资料解释 、时深转换及其地质剖面图生成过程中存在的问题 ,提出了 实际调查航迹的坐标投影 、浅地层剖面反射层位人机交互解释拾取 、时深转换层速度的输入等方法 , 并应用这些方法快速有效地生成了工程调查数据地质剖面图 。 关键词 :水深剖面 ;浅地层剖面 ;时深转换 ;地质剖面图 中图分类号 : P714 . 8文献标识码 : A 定位同步工作 ,以获得准确的调查位置数据 。但由于 引言0 实际调查航迹与设计测线之间存在固有偏差 ,同时受 海底浅地层剖面探测与水深测量在海洋地质调 浅地层剖面资料解释方法和时深转换不同层位声波 查与研究工作中是不可缺少的重要调查内容 。水深 旅行速度差异等因素影响 ,要生成精确的地质解释剖 测量可提供实时 、连续的海底水深地形变化 ,而浅地 面图 ,需要进行专门研究 。层剖面探测可提供直观 、连续的海底地质结构特征与 实际工作中 ,调查船测量航迹与设计测线间常存 构造状况 ,为大陆架第四纪地质 、浅海灾害地质等研 在很大的偏差 ,因此在室内解释 、绘制海底地层剖面 [ 1 - 2 ] 究提供重要参考。根据海底探测目标以及声波 图前 ,必须先将实际测量航迹向设计测线上投影 , 以 频率作用范围的不同 ,水深测量与浅地层剖面设备的 解决剖面图位置的偏差问题 。室内绘制海底地质剖 换能器发射声波频率范围与强度存在很大差异 ,一般 面图的传统方法是 :人工在浅地层剖面记录纸上读取 高频声波用于测量中 、浅海水深或海底面扫测 ,低频 反射层位厚度 ,并手工绘制到经潮位改正后的水深剖 [ 3 ] 声波则用 于 测 量 深 海 水 深 或 探 测 浅 层 地 层 信 息。 面图上 ,然后经扫描 ,再经 CAD 数字化等系列工作 , () 所以单波束 或多波束测量设备采用高分辨率 、高灵 才能完成一幅地质剖面图 。可见 ,这种传统的绘制方 () 敏度的高频声波 一般几十至几百千赫,而浅地层剖 法不适合对长距离测线的剖面解释 ; 同时 ,手工绘制 () 面仪采用穿透力强的低频声波 一般 1,20 k Hz,两 的所有层位厚度都是由同一个声波在海水中的传播 者结合使用 ,就可获得海底水深地形与海底地质剖面 () 速度 一般取 1 500 m/ s计算所得 ,而由于岩性等的 及地质属性等信息 。 不同 ,声波的传播速度必然有一定的差异 ,因次 ,利用 ( 为了获得可靠的地质解释剖面图 ,单波束 或多 同一速度值所计算的层位厚度必然给深度成图带来 很大误差 。目前 , 国 际 上 大 多 采 用 专 业 的 商 业 软 件 ) 波束水深测量与浅地层剖面探测必须与 GPS 导航 收稿日期 :2008203218 ( ) :中国海洋大学海底科学与探测技术教育部重点实验室资助项目 2008205;国家海洋局第二海洋研究所基本科研业务费专项资助项 目 项 ( )目 SZ0812 、J T0901 ( ) 作者简介 :丁维凤 1978 - ,男 ,江西德安县人 ,高级工程师 ,主要从事海洋地球物理勘查与研究工作 。 (如 Trito n 公司的 SB2Int e rp ret e r 解释软件和 Che sa2 野外调查布置了 3 条测线 : 北测线 、中心轴线和南测 )() p ea ke 公司的 So na r WizMap + SB P 采集处理软件等线 表 1 、图 1。区域内海水深度为 6,20 m ,局部可 ? 进行计算和成图 ,但这些专业软件只能简单 、机械地 达 55 m 。 图 2 为长白支线现场调查航迹图 ,包括北测线 、 拾取反射层位 ,不能灵活地进行人机交互式解释 , 无 (法将解释结果与其它资料 如水深资料和侧扫声呐资 中心轴线 、南测线和与之垂直的 2 条联络测线 。野外 ) 料等进行有效结合 ;同时 ,这些软件也都是采用同一 调查对 3 条测线都同时进行了单波束水深测量与浅 () 个速度 声波在海水中的传播速度进行时深转换的 , 地层剖面探测 ,对中心轴线还作了 Boo me r 震源单道 未提供不同层位层速度的输入界面 ,因此其解释结果 地震探测 。 也一定存在较大误差 。 表 1 长白支线 3 条设计测线的端点位置 ( ) 本文利用 2006 年舟山市大陆连岛工程 三期岱 Tab. 1 The coo r dinate s of t hree Cha ngbai 山跨海大桥工程物探调查资料 ,对实际调查航迹向设 pla nned surveying line s 计测线投影方法 、浅地层剖面人机交互处理解释技术 测线 端点 X/ m Y/ m 以及时深转换方法进行了研究 ,并将解释结果与水深 C12QD1 北测线 3 341 190 . 134 41 405 394 . 793 剖面准确结合生成了该区的地质剖面图 。 C12QD2 3 342 077 . 389 41 402 131 . 070 1 资料与方法 中心轴线 C22QD1 3 340 900 . 641 41 405 316 . 094 C22QD2 3 341 787 . 896 41 402 052 . 370 舟山市大陆连岛工程是连接舟山本岛和岱山岛 C32QD1 3 340 611 . 148 41 405 237 . 394 南测线 的一项重要工程 ,三期工程阶段的长白支线东起长白 3 341 498 . 402 41 401 973 . 670 C32QD2 岛 ,西接舟山本岛至岱山岛的西线 ,全长约 3. 3 km , 图 1 长白支线测线的位置 Fig. 1 The locatio n of t he Cha ngbai surveyi ng line s ( ) ?杭州国海海洋工程勘测设计研究院. 舟山大陆连岛工程 三期岱山跨海大桥工程可行性研究专题浅地层剖面探测报告 ,2006 . 丁维凤 等 :浅地层剖面与测深数据综合成图法在舟山大陆连岛工程中的应用 ?65 ? 图 2 长白支线的调查航迹图 Fig. 2 The navigatio n t rack of Changbai surveying line s 表 2 提取的轴线剖面部分数据和计算结果 2 结果与讨论 Tab. 2 The p art of hydro grap hic data and co mp uti ng re sult s of ce nt er sur veyi ng li ne 2 . 1 长白支线轴线的水深剖面 Event 水面 水深 水底 将 H yp ack 导航软件 输出 的导 航 定位 文件 按 一 时间 X/ m Y/ m 号 高程/ m / m 高程/ m 定时间间隔提取 GPS 位置与水深数据 ,并根据 1985 13 : 12 :21 105 3 341 418 . 03 403 423 . 34 1 . 81 8 . 23 - 6 . 42 ( 国家高程 基准 面 对提 取的 水 深数 据作 潮 位改 正 表 13 : 12 :35 106 3 341 409 . 14 403 445 . 38 1 . 77 8 . 17 - 6 . 40 ) 2;将改正后的水深数据调入 CAD 等成图软件 ,根据 13 : 13 :23 107 3 341 379 . 97 403 526 . 19 1 . 75 7 . 98 - 6 . 23 的成图比例尺 ,用三角网格法插值计算出 D TM 13 : 14 :08 108 3 341 360 . 27 403 606 . 58 1 . 76 7 . 75 - 5 . 99 () 网格数据 ,并沿设计测线切出水深剖面图 图 3。 13 : 14 : 20 109 3 341 356 . 52 403 628 . 63 1 . 79 7 . 62 - 5 . 83 2 . 2 浅地层剖面人机交互解释与层位时深转换 2 . 2 . 1 浅地层剖面人机交互解释 浅海测深必需使用 采集的浅剖资料中水深数据的准确率在风浪影响下 高频率声波测量 ,由于 Chi rp2 也不高 。因此 ,在绘制地质剖面图时 ,不能以浅地层 III 浅地层剖面仪发射的 2 ,7 k Hz 的低频声波对海 剖面资料上的反射海底面作为地质剖面图中的海底 底软弱表层反应不敏感 ,再加上 Chi rp III 拖体的入水 面 ,而应以经波浪与潮位改正后的测深数据作为海底 深度难以人为控制等诸多因素的影响 ,浅地层剖面资 面 ,将浅地层剖面反射层的厚度数据与该测深数据结 料上显示的海底深度数据的精度并不高 。另外 ,现场 合后再成图 。作业时 ,浅地层剖面仪大多未接入波浪改正仪器 ,所 图 3 长白支线中心轴线的水深剖面图 Fig. 3 The bat hymet ric p rofile of Cha ngbai center surveying line 数字记录的浅地层剖面资料中 ,除包含有发射 、 据不会带来很大误差 ,但对不同岩性反射地层的浅地 接收的实时位置外 ,还存储有由导航软件按固定距离 层剖面数据却会带来较大误差 。 据已有的地质和钻输出的定位 Eve nt 号 ,通过该 Eve nt 号 , 既可随时确 探资料显示 ,调查区上部的第 ( 定剖面的航迹位置 ,还可与其它资料结合 如水深剖四系为亚粘土 、粘土 、细砂与砾石等 ,其下为侏罗系的 ) 面和侧扫声呐资料等进行准确定位 。在进行资料解 凝灰岩 等 , 凝 灰 岩 地 层 的 声 波 旅 行 速 度 远 远 超 过 ? 释时 ,利用浅地层剖面资料上的定位 Eve nt 号 ,可在 1 500 m/ s 。表 3 为调查区不同岩性地层声波的传 层位形态变化帄缓处只拾取 Eve nt 号位置 ,在层位形 播速度 。由表 3 可见 ,调查区不同岩性地层的声波传 () 态变化剧烈处加密拾取 不考虑反射波偏移等因素。 播速度是不同的 ,且多数岩性地层中声波的传播速度 另外 ,对于长距离测线 ,在计算机屏幕上人工用鼠标 都大于 1 500 m/ s 。以图 6 中的两个反射层为例 ,若 点击拾取 ,必然会影响拾取的精度与速度 ,若能开发 拾取的两地层声波传播速度都取 1 500 m/ s 计算的 计算机自动拾取功能软件 ,则可以加快拾取进度 , 提 话 , Eve nt 标号为 2 278 处拾取的第二层特征点深度 高解释的精度 。最后 ,为了能将浅地层剖面解释数据 值为 66 . 10 m ; 若考虑第二层以上的粘土 、细砂与砾 与测深数据结合成图 ,需计算出浅地层剖面各反射层 石层的岩性特征 ,取声波帄均传播速度为 1 850 m/ s , (相对于海底面反射层的厚度 该厚度需要经过时深转 第一层以上声波传播速度仍沿用 1 500 m/ s 进行计 ) 换获得,计算方法有两种 : 一种是由软件提供上 、下 算 ,则得到的特征点深度值为 80 . 39 m , 两值相差达 层位同位置拾取的功能 ,即使各反射层拾取的特征点 14 . 29 m ,且越向地层的深处 ,该差值越大 ,相对误差 与海底面反射层拾取位置保持一致 ,然后将拾取特征 可达 18 %左右 。所以取同一个声波传播速度值进行 点按上 、下相 同位 置 的深 度值 相 减即 可 得 到 厚 度 值 时深转换 ,其解释结果的误差是不可忽视的 。 石油地() 图 4a;另一种是由软件自动插值计算拾取特征点上 震勘探剖面的时深转换借助于确定的地 海底面反射层相同位置处的深度值 ,再将插值点与拾 震波速模型 ,而波速模型建立的关键是追踪经钻孔验 [ 4 ] 取点的深度值相减就可得到拾取特征点的厚度值证了的层位 ,利用钻孔标定的标准层深度和反射 () 图 4 b。 地震剖面上该层对应的双程走时确定时深转换模型 , 按照上述方法 ,本文研究开发了浅地层剖面定位 并利用声测井 、中子密度测井和 V SP 资料合成地震 [ 6 ] Eve nt 号拾取 、计算机自动拾取 、等 pi ng 间隔拾取 、手 图作检验。可见 ,若采用地震波速模型方法对声学 工拉直线多点拾取与层位厚度的自动计算等系列功 地层剖面记录作时深转换 ,将很难得到用于合成声学 [ 4 - 5 ] 能的解释处理软件。图 5 和图 6 分别为利用该 记录剖面的测井数据 ,同时 ,合成的模型还需要不断 套开发软件对浅地层时间剖 面和 Boo me r 震源 单 道 的验证与改进 ,运算起来非常繁琐 。 地震时间剖面资料进行人机交互解释的结果 。为了避开上述问题 ,本文只对解释拾取的反射层 位在输出时作时深转换 ,未对浅地层剖面记录本身作 转换 。利用这种方法计算 ,需要用户明确各解释层的 层速度以及目的层以上已有的解释层位 。图 7 为反 射层解释拾取点处的厚度与深度计算示意图 。各层 相对于第一层的相对厚度计算公式分别为 H= 0 ( )1 1 () ( ) 4 上 、下层同位置拾取 a与上层插值计算 b的特征点 图 ( ) ( ) 2H= V t- t/ 2 +H 2 2 2 11 1Fig. 4 The sketch map of calculating ( )) t he picking point sπ t hickness ( H3 - H t/ 2 + H= V t2+ 122 3 3 3 各层的深度计算公式分别为 2 . 2 . 2 解释层位的时深转换 ( ) t / 2 4 D= V 1 1 1 目前很多商业软件都没有对时深转换功能作专 DDH = + ( )2125 门的研究与设计 ,都是简单地采用了声波在海水中的 ( ) 6 3 DHD = + 1 3 帄均传播速度作时深转换的 。该方法虽对测深的数 式中 : H、H和 H分别为 3 个层位相对于第一层的 1 2 3 ?浙江省地球物理勘查院. 岱山跨海大桥工程可行性研究专题工程地质物探勘察报告 ,2007 . 丁维凤 等 :浅地层剖面与测深数据综合成图法在舟山大陆连岛工程中的应用 ?67 ? 相对厚度 ; D、D和 D分别为 3 个层位的深度 ; V 、 8 为浅地层剖面交互解释拾取各解释层层速度 图 1 2 3 1 V 和 V 分别为 3 个层位的声波传播速度 ; t、t、t、与层位序号定义的界面 ,计算所得的解释层拾取特征 2 3 1 11 2 t和 t分别为声波到达 3 个层位的双程旅行时 。点的深度和相对于第一层的相对厚度值列于表4 中 。 22 3 2 . 2 . 3 声学地层剖面的偏移处理 偏移处理可使倾斜地层面上的反射波 、断层面上 的断面波 、弯曲界面上的回转波以及断点和尖灭点上 的绕射波归 位和 收敛 , 得 到地 下 反射 界 面 的 真 实 位 置 、构造形态和清晰可辨的断点与尖灭 ,这对提高剖 [ 7 ] 面资料的分辨率具有重要作用。 调查区浅地层时间剖面资料显示 ,该区浅部反射 地层倾斜角度很小 ,界面反射波组连续 ,不存在活动 断层与地层尖灭等地质现象 ,可不对其进行偏移处理 () 图 5。而 Boo me r 震源单道地震时间剖面上第二解 释层的反射界面起伏较大 ,反射波组不连续 ,应对其 () 进行偏移处理 图 6。由于时间方面的问题 ,本文在调 查时未对 Boo mer 单道地震资料作专门的偏移处理 。 ?表 3 不同岩性地层声波的传播速度 Ta b. 3 The so nic p ropagating velocit y in ?t he diff erent st rat ums of surveying a rea - 1 项目 岩性 ) ( 纵波速度/ m ?s 海水 — 1 450,1 500 ( )1 300,1 550 第四系 Q 淤泥质亚粘土 亚粘土 、粘土 1 500,1 900 粉细砂 1 600,1 900 砾石 熔结凝1 800,2 300 灰岩 凝灰质4 000,4 500 ()第四系上侏罗系 J 3 砾岩 4 000,5 000 霏细斑岩 4 000,4 800 表 4 浅地层时间剖面交互解释拾取的部分输出结果 Ta b. 4 The pa rt of picking point s abo ve sub2bo t to m interactive picking 相对于第一层的层位特征点厚度/ m 层位特征点深度/ m (())号 Event 纬度 N? 经度 E? 第二层 第三层 第四层 第五层 第二层 第三层 第四层 第五层 105 30 . 187 821 121 . 998 144 4 . 54 10 . 69 20 . 02 17 . 39 9 . 76 15 . 90 25 . 24 22 . 60 106 30 . 187 699 121 . 998 665 4 . 54 10 . 88 19 . 88 18 . 12 9 . 58 15 . 91 24 . 91 23 . 16 107 30 . 187 533 121 . 999 137 4 . 64 10 . 69 19 . 80 17 . 60 9 . 49 15 . 54 24 . 65 22 . 45 108 30 . 187 447 121 . 999 642 4 . 73 10 . 40 19 . 52 17 . 10 9 . 58 15 . 25 24 . 36 21 . 95 109 30 . 187 333 122 . 000 163 4 . 91 10 . 49 19 . 50 16 . 86 9 . 58 15 . 16 24 . 16 21 . 53 并计算出投影点在设计测线上的帄面坐标以及该投 3 海底地质剖面图的生成 影点上的改正水深值 ,然后将表 4 中的厚度值与该水 为解决浅地层剖面 、Boo mer 震源单道地震数据 ( 深值相加 , 得 到 解 释 层 位 各 特 征 点 的 高 程 图 9 、表 ) 与设计测线之间的位置偏差问题 ,在生成海底地质剖 5,最后将这些高程与帄面坐标数据调入 A uto CAD 面图前 ,需将拾取的浅地层时间剖面与 Boo me r 震源 中 ,通过专门编制的 L i sp 语言模块 ,绘制出最终的地 () 单道地震位置坐标向设计测线投影 。投影方法为 : 将 质剖面图 图 10。 丁维凤 等 :浅地层剖面与测深数据综合成图法在舟山大陆连岛工程中的应用 ?69 ? 图 9 特征点坐标向设计测线上的投影 Fig. 9 The sketch map of picked t ype point s p rojectio n to pla ned line 表 5 浅地层剖面和 Boomer 震源单道地震的解释层位特征点高程 Ta b. 5 The layer sπ dep t h of integrating interp reted layer s wit h hydro grap h data 层位特征点高程/ m 投影后 X/ m 投影后 Y/ m 第一层 第二层 第三层 第四层 第五层 第六层 3 341 392 . 630 403 506 . 338 - 6 . 42 - 10 . 96 - 17 . 11 - 26 . 44 - 23 . 81 - 106 . 03 3 341 378 . 999 403 556 . 479 - 6 . 40 - 10 . 94 - 17 . 28 - 26 . 28 - 24 . 52 - 105 . 15 3 341 366 . 234 403 603 . 433 - 6 . 23 - 10 . 87 - 16 . 92 - 26 . 03 - 23 . 83 - 104 . 87 3 341 353 . 274 403 651 . 106 - 5 . 99 - 10 . 72 - 16 . 39 - 25 . 51 - 23 . 09 - 103 . 23 3 341 339 . 704 403 701 . 022 - 5 . 83 - 10 . 74 - 16 . 32 - 25 . 33 - 22 . 69 - 102 . 14 图 10 地质剖面图 Fig. 10 The re sult of geolo gical p rofile 成 D TM 网格数据的方法 , 快速 、有效地解决了水深 4 结语 测量实际调查航迹与设计测线间的偏差问题 ,并生成 地质剖面图是海洋地质调查中的重要成果图件 。 了沿设计测线的水深剖面 。利用浅地层时间剖面与 为了准确地生成地质剖面解释图件 ,通常需要逐步解 Boo me r 震源单道地震数据的定位 Eve nt 号人机交互决调查航迹的位置投影 、水深剖面的生成 、浅地层剖 拾取 、计算机自动拾取 、等 pi ng 间隔拾取 、手工拉直 面与 Boo mer 震源单道地震反射层位的解释拾取 、层 线多点拾取 与层 位厚 度的 自 动计 算 等 诸 多 方 法 , 快 位的时深转换 、水深剖面与反射层位的结合等一系列 速 、准确地解释了浅地层剖面与 Boo me r 震源单道地 震数据中的反射层位 。通过研究 、设计的浅地层时间 问题 。 本文采用目前成熟的专业软件 ,通过三角网格生剖面与 Boo mer 震源单道地震解释 层 位层 速度 输 入 ( ) 术研究[J ] . 声学学报 ,2001 ,26 4:365 - 371 . 等方法 ,解决了时深转换时误差大 、运算繁琐等棘手 [ 2 ] 夏美永 ,米晓利. 浅海剖面仪在海洋工程中的应用 [ J ] . 物探装备 , 的时深转换问题 。另外将拾取的浅地层时间剖面与( ) 2002 ,12 1: 52 - 54 . Boo me r 震源单道地震位置坐标向设计测线投影 , 并 [ 3 ] 吴自银 ,郑玉龙 ,初凤友 , 等. 海底浅表层信息声探测技术研究现 插值计算出设计测线上对应的帄面坐标 、提取 D TM ( ) 状及发展[J ] . 地球科学进展 ,2005 ,20 11:1 210 - 1 217 . 网格数据中的水深高程 ,利用 A uto CAD 中专门编制 [ 4 ] 丁维凤 ,罗进华 ,来向华 , 等. 浅地层剖面交互拾取解释技术研究 的 L i sp 语言模块 ,实现了将水深剖面 、浅地层剖面与 ( ) [J ] . 海洋科学 ,2008 ,32 9:1 - 6 . Boo me r 震源单道地震反射层位数据结合 ,生成了较 [ 5 ] 罗进华 ,丁维凤 ,潘国富. 能量比值法在浅地层剖面海底面自动追 ( ) 踪中的应用[J ] . 物探与化探技术 ,2009 ,32 1: 20 - 24 . 准确的地质剖面图 。 [ 6 ] 杨文采 ,杨午阳 ,程振炎. 中国大陆科学钻探孔区的地震波速模型 ( ) [J ] . 地球物理学报 ,2006 ,49 2:477 - 489 . [ 7 ] 黄德济 ,贺振华 ,包吉山. 地震勘探资料数字处理 [ M ] . 北京 : 地质 参考文献 : 出版社 ,1990 . [ 1 ] 季文赟 ,林亦俊 ,张叔英. 地质声呐记录的图像处理与地层识别技 丁维凤 等 :浅地层剖面与测深数据综合成图法在舟山大陆连岛工程中的应用 ?71 ? Geol ogic prof ile chart ing by integra t ing sub2bottom prof ile an d ba thymetr ic da ta D IN G Wei2f e ng , L IU Kui , HA N Fu2jia ng , L U Xiao2mi ng , Q IA N Wa n2mi n ( )S econ d I ns t i t ute o f Oce a no g r a p h y , S O A , H a n g z hou 310012 , C hi n a Abstract : Itπs ver y i mpo r t a nt fo r hydro grap h sur veyi ng sync hro no u sl y wit h sub2bo t to m e xp lo ri ng fo r ocea n geolo gical i nve sti gatio n . Whe n c ha r ti ng t he fi nal geolo gical p rofile t he re a re so me t ro u ble so me p ro ble ms to be sol ve d. Fi r st , t he act ual navi gatio nπs t rac k i s dep a rt ure f ro m de si gned li ne s w he n spo t sur veyi ng , a nd t he sub2 bo t to m p rofileπs diff e re nt i nt e rp ret atio n met ho d ma ke t he diff ere nt i nt e rp ret e d re sultπs p reci sio n ; Seco nd , w he n ti me2dep t h co nver sio n we mu st co n si de t he diff e re nt st rat u msπ velocit y beca u se of t he so nic diff e re nt p rop a gati ng velocit y i n diff ere nt st rat ums , w hile e xi sti ng co mmercial sof t s t reat t hi s p ro ble m wit h o ne velocit y ; Thi r d , w he n cha r ti ng t he geolo gical p rofile , we mu st i nt e grat e t he hydro grap hπs dat a wit h t he st rat u msπ t hic k ne sse s co mi ng f ro m su b2bo t to m p rofile . Ba sed o n t he act ual i nve sti gat ed dat a , t he a ut ho r a nal yse s a bo ve p ro ble mπs p ro duci ng rea so n , a nd t a ke vali d met ho ds to sol ve t he se p ro ble m s. To sol ve t he fi r st p ro ble m , a ut ho r u se t he i deal co mme rcial sof t to ma ke f a st l y t he hydro grap h p rofile a nd bat hymet ric dat a alo ng t he p la nne d sur veyi ng li ne , w hile develop e d my self a i ndep e nde nt sof t to i nt e ractive i nt e rp reti ng t he sub2bo t to m p rofile , t hi s sof t off e r s eve nt ma r k locatio n picki ng a nd picki ng ra ndo ml y each pi ng a nd li ne picki ng a nd a uto picki ng et c . To sol ve t he seco nd p ro ble m , a ut ho r off e r s a sof t i nt erf ace to defi ne diff ere nt st rat u mπs velocit y a nd t he sof t co mp ut e s t he i nt e rp ret ed st rat umπs t hick ne ss to sea be d ref lecti ng layer a nd dep t h u si ng defi ni ng velocit y , at la st , t he sof t o utp ut s each st rat umπs t hick ne ss a nd dep t h dat a acco r di ng to u se rπs defi ni ng file na me . The t hi r d p ro ble m ca n be sol ve d wit h i nt egrati ng t he o utp ut st rat umπs t hic k ne ss wit h bat hymet ric dat a , b ut befo re i nt e grati ng t he t wo dat a , t he picked sub2bo t to mπs t rack po sitio n mu st be p roject ed to sur vey p la nned li ne , o t her wi se , t here i s po sitio n er ro r bet wee n pic ke d su b2bo t to mπs t rack po sitio n a nd hydro grap h p rofileπs po sitio n . Wit h all of t he a bo ve met ho ds , we i nt erp ret e s f a st l y t he su b2bo t to m p rofile a nd p ro duce t he hydro grap hy p rofile , at la st we buil d t he geolo gical p rofile wit h t he t wo dat a . All t he i nt e rp ret ed dat a a nd geolo gical p rofile a re sub mit t e d to t he clie nt s , a nd we a re app ro ved by t he clie nt s fo r all of t he se efficie nt met ho ds. Key words : bat hymet ric p rofile ; sub2bo t to m p rofile ; ti me2dep t h co nver sio n ; geolo gical p rofile