在MAPGIS下快速实现柱状剖面图的绘制.pdf

第2 6 卷 第2 期2 0 0 4 年 5 月 物探 化探 计 算技 术 V o 1 ． 2 6 N o ． 2 M a y 2 0 0 4 C OMP UTI NG TE C HNI QUE S F OR GE OP HYS I C AL AND GE O C HE MI C AL E XP L ORATI ON 文章编 号 1 0 0 1 1 7 4 9 2 0 0 4 0 2 --0 1 7 3 0 4 在 MAP GI S下快速 实现柱 状剖面 图的绘制 李沙园，陈昕华，戈永怡 中国石化勘探 开发研 究院舍肥化探 所 ，安徽合肥2 3 0 0 2 2 摘 要 将 油 气化 探 数 据 用 F ORTRAN 短 程 序 生成 MAP GI S点、 线 、 面 明码 文 件 后 ， 再利 用 MAP GI S图形编 辑软件 对 图形进 行编辑加 工 ， 达到快速制 作油 气化探 柱状 剖面 图的 目的 。使 用 该方法 能充分 利用 MAP GI S系统 资源， 降低 编制绘 图程序的难度 。 关键词 MAP GI S；制作柱状 剖面 图；油气化探 中图分类号 P 6 3 2 文献标 识码 A A F AS T DRAW I NG M ETHOD FoR CoLUM NAR S ECTI ON US IN G M APGI S LI Sh a y u a n ，CHEN Xi n h u a ，GE Yo n g y i P e t r o l e u m Ge o c h e mi c a l P r o s p e c t i n g I n s t i t u t e ，S I NO PE C， 崩An h u i 2 3 0 0 2 2，C h i n a Ab s t r a c t Th i s a r t i c l e i n t r o d u c e s a me t h od t O c o n v e r t d a t a o f g e o c h e mi c a l e x p l o r a t i o n f o r o i l a n d g a s i n t o M AP GI S p o i n t ，l i n e a n d s u r f a c e c od e f i l e s u s i n g a s h o r t F ORTRAN p r o g r a m a n d t h e n t h e c o l u mn a r s e c t i o n o f g e o c h e mi c a l e x p l o r a t i o n c a n b e d r a wn b y M APGI S e d i t s o f t wa r e ． Th e me t h od wi t h M APGI S r e s o u r c e s c a n r e d u c e t h e d i f f i e u l t i e s e n c o u n t e r e d i n c o n v e n t i o n a l p r o g r a m d e s i g n f o r c o l u mn a r s e c t i o n d r a wi n g ． Ke y wo r d s M AP GI S；c o l u mn a r ma p d r a wi n g；g e o c h e mi c a l e x p l o r a t i o n f o r o i l a n d g a s 0 引言 获 2 0 0 2年度 国家科 技进 步二等奖 的 MAP GI S软件 ， 作 为 国产 工具 型 GI S软件 ， 其 强大 的功能 ， 良好 的技术 支持 ， 已经 获得普遍 的认可 ， 使用 该系统软件 编辑 的图件具有 方便 、 速 度快 、 精度 高 的特 点 ， 在地学 领域 具有 广大 的用户群 。 油气化探柱状 剖 面图是将 岩性 剖 面和多个化探 指标浓度 曲线在 同一个 图上按 其取 样深 度表示 出来 ， 并在 图上标注地层 岩性 说 明、 化探采样点 及样号 。笔者 利用 自编程序 ， 把 已有 的地层岩性 数据及化探数 据 文件转换 成 MAP GI S点 、 线 、 面明码文件 ， g N, l MAP GI S图形 编辑软件 对图形 进行编 辑加工 ， 最后形 成 一 幅理 想的柱状剖 面图 。该其方 法简便 、 易行且速度快 ， 在实 际应用 中取得 了较好 的效 果 。 1 作图的基本思想 为了表示不 同的岩性 ， 需 要编程形 成很多 图库 ， 以便绘 图时调用 。 其工作量 较大且有 一定难度 ， 又不直 观 。而 MAP GI S本 身带有近千 种子图库 符号库 、 线型库 、 图案库 ， 要 编辑添加 子图也很 方便 、 直观 。因此 收稿 日期 2 O O 3 一O 9 一O 2 维普资讯 1 7 4 物探 化探 计算技 术 2 6卷 利用 MAP GI S本身的子 图库 ， 再根据 需要 编辑添 加一些子 图 ， 就可作 为岩 性 图库 ， 而使 绘 图程序 简单化 ， 稍有 编程能力即可实 现 。 2实现化探柱状 剖面图的绘制 由上所述 ， 在做图前先要编辑添加一些子图， 作为岩性图库时， 可利用 MA P G I S本身的子图库就能方 便地制作添加 。图 1 是我们制作的部分图案。 2 ． 1 化探柱状 剖面 图的绘制过程 用此 方法绘 制油气化探柱状 剖面 图的过 程见图 2所示 。 序 号 图 案 岩 性 l J 固J I 白 云 质 灰 岩 与 l 9 9 l 1 1 l 藻团块灰岩互层 I l l l 2 0 o I I l I 灰岩 l 1 1 f 2 0 l 白云质灰岩 I 1 l l I 够 l 国 1 0 l 2 0 2 簟 粒灰 岩 l I 母 I J f I I I I 2 0 3 l I 、 ， I I 含 生 物 屑 灰 岩 陶辑绘图敌据文件 l 程 序 睦序转换明码文件 I 【形成点强 、 面文件 l 编辑点 线 、 面文件 I 形成工程文件 I 输 出 图 1 部 分 岩性 图例 圈 2 简单绘 图漉程 Fi g． 1 Th e p a r t s o f l i t h o l o g y l e g e n d s Fi g ． 2 Th e g e n e r a l d r a wi ng flo wc h a r t 从上述 流程 图可 以看 出 ， 所 绘制 的剖 面 图是将化 探数 据用 程序 生成 MAP GI S明码文 件后 再转 换 为 点 、 线、 面文件 ， 以图形 叠加 的方式制 作 的。 2 ． 2编程实现 将化探数据转换成 MA P G I S明码文件的 F O R T R AN程序主要部分如下 。 参数说 明 n s c a l e 比例尺 y i z l i o u t f i l e 2 na mel na me 2 h wd 每 层厚度 累积厚度 ．wa l 文件名 ．wa p文件名 ．wa t 文件名 图的长度 图的宽度 岩性分层 数 化探采样 数 第 i 列宽度 第 i 层 岩性代号 ， 和 区参 数 中的图案对应 略 程序 头 ， 定 义数组 、 变量 、 读数据 i x 0 0 C C i x l 一1 x3 0．3 x00．0 i y0 0 d d n s e a l e ／ 1 0 0 0 ． 0 y 0 h ／ d d x d wd 打开 MAP GI S点 、 线、 面文件并写 文件 头 o p e n 5， f i l e o u t f i l e 2 o p e n 6， f i l e n a me 1 o p e n 7． f i l e n a me 2 wr i t e 7， ’ 1 x， a ’ ’ W M AP9 0 2 2 ’ wr i t e 6， ’ 1 x， a W M AP9 0 2 3 wr i t e 5， ’ 1 x， a ’ W M AP9 0 2 1’ 读数据 xl x O a 1 xl 1 xl一 3 O． 1 i o ． n n a C 维普资讯 2期 李沙 园等 在 MAP GI S下快速 实现柱状剖 面图的绘制 1 7 5 x2 xl一 一 20 yl y0 nn 0 y d一 0 ． nn n 0 d o 6 0 0 i 一 1 ．n1 y d y d y i C 判 断是否是 相 同的岩性 ， 如果 两层 相 同则 归 并为一层 ， 否则 进一步判断后 画图 i f i c i 一 i c i 1 ． n e ． 0 t h e n C 判断 厚度 是否 ≥ 2 ram，2 ram 以上可 以作 岩 性 图 i f y d ． g t ． 2 ．*d dt h e n n n n n n n 1 i a n n n 一 i n n n i l 1 n t n t 1 h dz l i ／ d d y 2 y 0 -- h d C 在面文件 中画线 ． wa p wr i t e 6 ，’ 1 x， a ’ 1 ，0 ，1 ，1 ，1 0 ．0 0 0 0 0 0，1 0 ． 0 0 00 00， 1， 0， 0’ wr i’ t e 6 ， 1 x， i l， ’ ， ’ i 1 i x 0， i x 0 wr i t e 6 ， ’ ’ 1 x， i 6 ， ’ ， ’ ， i 2 n n， i x0 wr i t e 6 ， 4 O O 5 wnt e 6 ， 2 0 0 xl ， y l wr i t e 6 ， 2 0 0 x2， y l wr t t e 6 ， 2 0 0 x 2， y 2 wr i t e 6 ， 2 0 0 xl ， y 2 wr t t e 6 ， 2 0 0 xl ， yl wr t t e 6 ， 1 x， i 6， ’ ， ’ ， f 3 ． 1 n n， x 3 e n d i f C 画线 ． wa d wr i t e 5，’ 1 x， a 1 ，0，1 ，1，1 0 ．0 0 0 0 0 0，1 0 ． 00 000 0， 1， 0， 0 wr i t e 5， 4 0 O 5 wr i t e 5， 2 0 0 xl ， yl wr i t e 5 ， 2 0 0 x 2， y l wr i t e 5 ， 2 0 0 x 2， y 2 wnt e 5 ， 2 0 0 xl ， y 2 wnt e 5 ， 2 0 0 x l， y l wr i t e 5， 1 x， i 6， ’ ， ’ ， f 3 ． 1 n t ， x 3 y d一 0 ． yl y2 e n d i f 6 00 c ont i nue C 造 区用 I C n n 来 表示不 同的岩性 wr i t e 6， 4 0 0 i x1 wr i t e 6 ， 4 0 0 n n n d o 7 0 0 i 一 1， n n n n n i a i wr i t e 6， 1 x， a， i 4， a， i 6 ， a 9， ， i c n n ， ， 3． 00 00 00， 3．0 0 000 0， 1， 0， 0， 0， ， i ， ， 0， 0 wr i t e 6， ’ 1 x a ’ ’ 2 wr i t e 6， 4 0 0 i wr i t e 6 ， ’ 1 x， i l ’ i x 0 7 0 0 c o n t i nu e z 0 0． c 对 应 的化 探 指标 曲线 图 ， 为 了 图形 的美 观 ， 用 对数数据作 图 ， 标 注为实际值 d o 8 0 0 i k k．n nt 一 1 wr i t e 5 ，’ 1 x， a’ ’ 1， 0，1，1 ，1 0 ．0 0 0 0 0 0，1 0 ． 00 000 0， 1， 0， 0’ wr i t e 5， 4 0 0 n 2 xl x i yl y 0 t i t n t 1 i f z 2 i k k 1 ． I t ． 0 ． t h e n z 0一 一 1 ． e n d i f C 画刻度 d o l 1 0 0 k 1， 4 x 2 一 x l k*1 0 ． wr i t e 7， 1 x， f 1 1 ． 6， ’ ， ， f 1 1 ． 6， a x 2， y 0， ， 0， 1 ， 5 6， 1 ． 5 0 0 0 0 0， 1 ． 5 0 0 0 0 0， 0 ． 0 0 0 0 0 0， 1， 1， 1 ， 1 ， 0 ’ wr i t e 7， ’ ’ 1 x， f l 1 ． 6， ， ’ ， f 1 1 ． 6 ， a ， f 8 ．1 ， a x 2 2 ．， y 0 2 ．， ’ ， 0， 0， ’ ， z 0 a b s z 2 i k k 1 *k， ’ ， 1 ． 5 0 0 0 0 0， 1 ．5 00 000， 0， 0． 000 00 0， 1， 1， 0， 0， 1， 2， 0 1 1 00 c o nt i nu e C 画化探 曲线 d o 1 0 0 0 j 1 ， n 2 x 2 x l s j ， i k k 1 --z 0 ／ a b s z 2 i k k 1 * 1 0 ． h d z j ／ d d yl y 0-- h d wr i t e 5， 2 0 0 x 2， y l 维普资讯 1 7 6 物探 化探计 算技 术 2 6卷 1 00 0 8 00 ⋯⋯ 标 注化探样点 、 画地 质界 、 系、 统 、 组 、 段分界线及标注 岩层厚度 可 以看 出程序 比较简单 。将形 成的 MAP GI S明码 文件 ， 转 换成 MAP GI S认 可 的点 、 线 、 面文 件 ， 并 利 用 编辑 子系统 ， 对 图形进行修 饰 、 加工 ， 就可转换成 工程文件输 出。 2 ． 3 运行环境 本程 序 在 F o r t r a n P o we r S a t i o n 4 ． 0下 编号 ， 在 Wi n d o ws 9 8以上 操 作 系统 运 行 ， 形 成 的 图形 文 件在 MAP GI S 5 ． 0以上版本均 可运行 。 3 实例显示 下 面给 出用 上 述 方 法绘 制 的 某 地 区 地质 剖面化 探指标纵 向曲线 图的实例 。由 于原 图较 大 ， 只给 出部分 图件 内容 如 图 3 所示 。 4 作用及 意义 该 方法的应用 ， 利 用 了 MAP GI S平 台 的 自身功能， 结合 F O R T R A N程序 ， 直接 地层单位 岩性 酸解烃甲烷 奠解烃重烃 厚度 化探点及编号 界 系 统 组 段 代号 剖面 { m 由 下 T 1口 马 ● ● ● ● ● 执 - 松 由 多 组 下 T l m 第 蓍 耋 I 生 五 ● 叠 段 T ￡ 血 第 霾 ● 四 统 段 T lz ‘ 第 蘑 i 界 三 系 段 T lg3 转换 油气 化探数据 。制作容 易 ， 方法简便 ， 圈 3 地质剖面与化探指标纵向曲线图示例 缩减 了时间， 降低 了成本及编程难 度， 适 F i g ． 3 T h e d e m o n s t r a t i o n o f g e o lo g ic a l s e c t io n a n d v e r t ic a l 合于时间紧任务急的应用项 目， 尤其适用 c u r v e o f g e o c h e m ic a l e x p lo r a t io n in d ic a t o r 于 没 有 制 作 岩 性 剖 面 图 相 应 的 程 序 ， 而 有 MAP GI S软 件 的 人 或 单 位 。 由 于利 用 F OR TRAN 编 程 和 MAP G I S平台作图， 图形精度高 ， 为在 MAP G I S平台下制作油气化探柱状剖面图提供 了一种可选择的途 径 。 5 建议 该 方法是针对应用 项 目而提 出的 ， 偏 重实用 ， 具成 图快 的特 点 。 用别 的语 言代替 F ORTR AN语 言编程 可 ， 进一步加强对 字符 的操作 能力 ， 而且对一些 岩性符号 的表示 ， 用线子 图代替 面图案将更 准确 。 参考文献 [ 1 ] 中国地质大学信息工程学院 ． MAP GI S平台地理信息用户教程[ M] ． 武汉 武汉大学 出版社 ， 1 9 9 8 ． [ 2 3 桂 良进 ， 王军 ， 董波 ． F o rtr a n P o w e r S t a t i o n 4 ． 0使 用 与编 程[ M ] ． 北 京 北 京航 空 航天 大学 出版 社 2 0 0 0 ． 1- 3 3 敖 文波 ， 周 蒂 ， 胡光 道 ． MAP GI S下 重磁 数 据直 接成 图的 实现 及 意K E J ] ． 物探 化 探计 算 技术 ， 2 0 0 1 ， 2 3 3 2 8 [ 4 3 曹静 平 ， 王文联 ．如 何应 用 1 2 0 0 0 0 0化 探 数据 库构 造 GI S图层 [ J 3 ． 物 探 化探 计算 技 术 ． 1 9 9 8 ， 2 0 2 1 3 0 [ 5 3 张滨生 ．利用 MAP GI S快速制作点位图的一种方法[ 3 3 ． 物探化探计算技术， 2 0 0 1 ， 2 3 1 9 0 ． 作者简介 李沙园 1 9 5 6 一 ， 女， 高级工程师， 1 9 8 0年毕业于北京大学数 学系应用数学专业， 现从事油气 化探数 据处理及信 息管理 工作 。 Lr； 一 维普资讯