杨凯

（中国人民武装警察部队黄金第二支队，内蒙古呼和浩特 010010）

一种基于FORTRAN的多测深剖面水平切面提取方法

杨凯

（中国人民武装警察部队黄金第二支队，内蒙古呼和浩特 010010）

测深是勘察地下电阻率或其他物理特征分布的重要地球物理方法，主要分为电阻率测深、激电测深、电磁测深等，一般只需得到测深剖面的一二维断面图就可达到一定的勘探目的，但有些时候为了从不同角度、不同侧重点来研究异常，还要有不同深度的水平电阻率等值线图。基于此，通过FORTRAN语言编程，利用二维插值原理来统一提取多条测深剖面相应深度的电阻率数据，进而绘制水平电阻率等值线图。

Fortran 90；测深；二维插值；水平切面

FORTRAN语言是世界上最早出现的高级程序设计语言，其以严谨、规范、高效的特点普遍应用于数值分析、工程计算和科学研究等领域［1］。利用FORTRAN可以对物探数据进行批量、高效的处理。本文就是笔者在实践工作中探索出的根据二维差值原理利用FORTRAN编程实现从多条测深剖面数据中提取相应深度水平切面的过程。

1 基本原理

根据给定矩形域上n×m个结点上的函数值，利用二元插值公式计算出指定插值点处的函数值。

已知矩形域上n×m个结点的2个方向上的坐标分别为：

其相应的函数值为：

计算插值点（u，v）为中心，在X方向上，前后各取4个坐标：

在Y方向上，前后也取4个坐标：

然后利用二元插值公式：

计算插值点（u，v）处的函数值。

根据以上原理即可通过二元插值在电阻率等值断面图GRD网格化文件中提取出所需深度的电阻率值。……