测井曲线wis文件格式转换为ASCII文本格式小软件的开发（C与Python结合开发）

2024-04-04 11:28:18

在油田开发地质研究工作中，测井曲线是必不可少的数据文件之一。Forward软件是测井软件中非常优秀的产品，Forward中使用的是wis二进制格式，可是很多其他类型的软件不能直接识别wis格式，所以需要我们将wis格式文本格式。Forward软件有一个transfer小软件，不能实现批量处理，每次只能转换一个文件，但实际工作中通常一次要转换上百个文件，这就有点不方便了，因此，我在业余时间，阅读了Forward软件中自带的wis格式说明（如下，第一章内容），利用C和Python编写了一个可以批量处理的小软件。

第1章 WIS测井数据文件格式

WellBase平台底层数据文件使用WIS格式，WIS文件分为文件头、对象入口记录和对象数据三部分。数据存放以块为单位。

WIS文件能存放三种类型的对象，根据对象的类型分为通道对象、表对象和流对象。通道对象一般用来存放采集数据和处理结果（如测井曲线），表对象用来存放二维表数据（如解释结论），流对象用来存放二进制数据块（如解释参数，用户数据）。

1.1 WIS文件结构

1.1.1 文件标识

WIS文件标识从文件偏移零开始，为10个字节的字符。当前版本的标识为WIS 1.0。

1.1.2 文件头结构

头结构紧接文件标识。描述WIS文件的公共信息。结构定义如下：

typedef struct tagWIS_HEAD

WORD MachineType;

WORD MaxObjectNumber;

WORD ObjectNumber;

DWORD EntryOffset;

DWORD DataOffset;

DWORD FileSize;

DWORD TimeCreate;

char Reserved[32];

偏移 字节数 描述

0 2 机器类型=1 为PC； =2为SUN； =3为IBM；

2 2 允许记录的最大对象数。缺省为512个，该值可以在文件产生时给出。

4 2 当前记录的对象总数（包括删除和抛弃的对象）。

6 2 块长。WIS文件对象占用的磁盘空间以块为单位，该值指示每一数据块的字节数。

8 4 对象入口记录从文件开始的偏移量。

12 4 对象数据记录从文件开始的偏移量。

16 4 WIS文件的字节数大小。

20 4 WIS文件产生的时间。

24 32 保留字节。

1.1.3 对象入口

对象入口描述每个对象的公共信息，开始位置由头结构给出。每个对象的描述信息前后相连。结构定义如下：

typedef struct tagWIS_OBJECT_ENTRY

short Attribute;

short SubAttribute;

DWORD Position;

DWORD BlockNum;

DWORD TimeCreate;

DWORD TimeWrite;

char Reserved[32];

}WIS_OBJECT_ENTRY;

偏移 字节数 描述

0 16 对象的名称，以零结尾的字符串。

16 4 对象的状态：=0为正常； =1为抛弃； =2为删除。

20 2 对象的主属性：=1为通道对象； =2为表对象； =3为流对象。

22 2 对象的子属性，描述对应主属性的子属性。

24 4 对象数据体从文件开始处的偏移量。

28 4 对象数据体占用磁盘的块数。

32 4 对象产生的时间。

36 4 对象最近修改的时间。

40 32 保留字节。

1.1.4 对象数据体

对象数据体记录各个对象的具体特性及数据。根据不同的主属性分三种类型。对象数据体在WIS文件中的位置由对象入口指定。

1.1.5 通道对象

通道对象用来存放采集和计算结果数据(如测井曲线)。分为通道信息和通道数据两部分。

WIS文件将在一定时空内对某一采集或计算的物理信息数据集统称为通道数据。通道信息描述通道数据的存放形式，分为基本信息和维信息，基本信息描述信息的基本物理含义，维信息描述信息的时空特性，可以等间隔(连续)或非等间隔(离散)。最大允许有四维信息，通道信息共占用一个块空间，结构定义如下：

typedef struct tagWIS_CHANNLE

char AliasName[16];

char AliasUnit[16];

WORD NumOfDimension;

WIS_CHANNEL_DIMENSION DimInfo[4];

}WIS_CHANNEL;

偏移字节数描述

0 8 对象的单位，以零结尾的字符串。

8 16 对象的别名，以零结尾的字符串。

24 16 单位的别名，以零结尾的字符串。

40 2 对象数据类型，参见3.2.1。

42 2 数据类型的长度。

44 4 对象的最小值(测井曲线缺省左刻度值)。

48 4 对象的最大值(测井曲线缺省右刻度值)。

52 2 保留字节。

54 2 对象维信息数。

56 4*56 对象维信息。

通道维信息结构定义如下：

typedef struct tagWIS_CHANNLE_DIMENSION

char AliasName[16];

float StartVal;

DWORD MaxSamples;

}WIS_CHANNEL_DIMENSION;

偏移字节数描述

0 8 维的名称，以零结尾的字符串。

8 8 维的单位，以零结尾的字符串。

16 16 维的别名，以零结尾的字符串。

32 4 维的开始值。

36 4 维的采集或计算增量。对于离散数据，该值为0，数据中记录该维的值。

40 4 维的数据采样点数。如果该值为0，采样点数为可变值，数据中记录该值。对于第一维数据，该值不能为0。

44 4 维的数据采样最大点数。该值仅当采样点数信息为0（可变采样点）时有效，该维信息在数据中所占用的字节数通过该值计算。

48 4 该维上每一采样点所占用的字节数。

52 2 维的数据类型，参见3.2.1。

54 2 保留字节。

通道数据从通道描述信息的下一块开始。

下面为一个包含深度和时间维的物理信息数据体的存放顺序。第一维为深度，第二维为时间。

[A1]+[N2]+[B1]+X1+[B2]+X2+···+[BN]+XN+

[A2]+[N2]+[B1]+X1+[B2]+X2+···+[BN]+XN+

[AN]+[N2]+[B1]+X1+[B2]+X2+···+[BN]+XN

A1，A2，··· ，AN代表深度值，当深度维信息结构中的Delta为零时，记录此值。

N2代表当前深度点上的时间采样点数，当时间维信息结构中的采样点数为零时，记录此值。

B1，B2，··· ，BN代表时间值，当时间维信息结构中的Delta为零时，记录此值。

X1，X2，··· ，N代表物理信息的值。

1.1.6 表对象

表对象用来存放二维表数据，分为表信息和表数据体两个部分。表信息由不同的表项组成，每一表项称为字段。表信息结构定义如下：

typedef struct tagWIS_TABLE

DWORD RecordCount;

DWORD FieldCount;

WIS_TABLE_FIELD *pField;

}WIS_TABLE;

偏移字节数描述

0 4 表的记录数。

4 4 表的字段数。

8 4 指向字段信息结构的指针。

字段信息结构定义如下：

typedef struct tagWIS_TABLE_FIELD

DWORD Reserved;

}WIS_TABLE_FIELD;

偏移字节数描述

0 32 字段的名称，以零结尾的字符串。

32 2 字段值的浮点类型，参见3.2.1。

34 2 字段值的长度。

表数据体(记录)从表信息记录的下一块开始。

1.1.7 流对象

流对象用来存放二进制数据块。开始为4个字节的无符号长整形数，代表数据流的长度。接着为该流的二进制值。

第二章 WIS测井数据格式转换文本格式C程序代码

代码如下：

/* * Copyright (c) 2011, Sun Yunqiang * All rights reserved. * Software: wis2txt_win * Version: 1.0 * Author: Sun Yunqiang * English Name: Alex Sun * Release Date: 18/03/2011 */ /** ***************************************************************** */ /** This Program is for converting wis format file to txt format file */ /** Command: wis2txt_win.exe args1 args2 */ /** Description: */ /** wis2txt_win.exe: the executable file as command */ /** args1: the first argument that is the wis file name */ /** args2: the second argument that is the txt file name */ /** ***************************************************************** */ #include <windows.h> #include <stdio.h> #include <stdlib.h> #include <string.h> #include <ctype.h> #include <time.h> #include <math.h> #define ROW 65536 #define COL 512 /* **************************************** */ /* Below blocks are the structs of wis file */ /* **************************************** */ typedef struct tagWIS_HEAD { WORD MachineType; // 0-PC 1-SUN 2-IBM 3-HP WORD MaxObjectNumber; WORD ObjectNumber; WORD BlockLen; DWORD EntryOffset; DWORD DataOffset; DWORD FileSize; time_t TimeCreate; char Reserved[32]; }WIS_HEAD; typedef struct tagWIS_OBJECT_ENTRY { char Name[16]; long Status; short Attribute; // 1-通道对象 2-表对象 3-流对象 short SubAttribute; // 1-曲线对象 2-波形对象 3-地层测试对象 4-时深时对象 DWORD Position; // 对象数据体从文件开始处的偏移量 DWORD BlockNum; time_t TimeCreate; time_t TimeWrite; char Reserved[32]; }WIS_OBJECT_ENTRY; typedef struct tagWIS_CHANNEL_DIMENSION { char Name[8]; char Unit[8]; char AliasName[16]; float StartVal; float Delta; DWORD Samples; DWORD MaxSamples; DWORD Size; WORD RepCode; WORD Reserved; }WIS_CHANNEL_DIMENSION; typedef struct tagWIS_CHANNEL { char Unit[8]; char AliasName[16]; char AliasUnit[16]; WORD RepCode; WORD CodeLen; float MinVal; float MaxVal; WORD Reserved; WORD NumOfDimension; WIS_CHANNEL_DIMENSION DimInfo[4]; }WIS_CHANNEL; typedef struct tagWIS_STREAM { DWORD Length; DWORD Offset; }WIS_STREAM; typedef struct tagWIS_TABLE_FIELD { char Name[32]; WORD RepCode; WORD Length; DWORD Reserved; } WIS_TABLE_FIELD; typedef struct tagWIS_TABLE { DWORD RecordCount; DWORD FieldCount; WIS_TABLE_FIELD *pField; }WIS_TABLE; typedef struct tagWIS_TABLE_DEFAULT_FILED { char Name[16]; char Alias[16]; char Unit[8]; char Type[8]; WORD Length; WORD Count; char DefVal[64][12]; } WIS_TABLE_DEFAULT_FIELD; typedef struct tagWIS_DEFAULT_TABLE { char Name[16]; char Alias[16]; char Attrb[8]; DWORD FieldCount; WIS_TABLE_DEFAULT_FIELD *pField; } WIS_DEFAULT_TABLE; /* 曲线的信息结构体*/ typedef struct tagCURVE_HEAD { char Name[16]; // 曲线的名字 char Unit[8]; // 曲线的单位 float startMD; // 曲线的起始深度 float stopMD; // 曲线的结束深度 int pointNum; // 数据样点数 int position; // 数据位置 struct tagCURVE_HEAD *next; } CURVE_HEAD; /* 曲线的数据结构体 */ typedef struct tagCURVE_DATA { char Name[16]; // 曲线的名字 char Unit[8]; // 曲线的单位 float Depth[ROW]; // 曲线的深度值 float Value[ROW]; // 曲线的值 struct tagCURVE_DATA *next; } CURVE_DATA; /* ************************************************************* */ /* Below block is the Main program of processing wis format file */ /* ************************************************************* */ int main(int argc, char *argv[]) { FILE *wisfile; // 输入文件 FILE *txtfile; // 输出文件 char *infilename = argv[1]; char *outfilename = argv[2]; // char *infilename = "demoa.wis"; // char *outfilename = "demoa.txt"; if(argc != 3 ) { printf("ERROR: Command parameters is not defined correctly!/n"); exit(1); } wisfile = fopen(infilename, "rb"); if(wisfile == NULL) { printf("ERROR: Read file terminated!/n"); exit(1); } txtfile = fopen(outfilename, "w"); if(txtfile == NULL) { printf("ERROR: Write file terminated!/n"); exit(1); } /* Process wis file, Read it and transfer it to text */ /** wis文件标识符从文件偏移0开始，为10个字节的字符 */ char wisid[10]; fread(wisid, 10, 1, wisfile); // printf("文件标识符: %s/n", wisid); fprintf(txtfile, "ASCII WIS %c%c%c/n", wisid[4], wisid[5], wisid[6]); /** 头文件紧接文件标识 */ WIS_HEAD *wishead; wishead = (WIS_HEAD *)malloc(sizeof(WIS_HEAD)); fseek(wisfile, 10, SEEK_SET); fread(wishead, sizeof(WIS_HEAD), 1, wisfile); // printf("对象入口记录从文件开始的偏移量EntryOffset: %d/n", wishead->EntryOffset); // printf("对象数据记录从文件开始的偏移量DataOffset: %d/n", wishead->DataOffset); // printf("当前记录的对象总数ObjectNumber: %d/n", wishead->ObjectNumber); /** 对象入口，位置由头结构中EntryOffset参数指定最多512条曲线 */ /** curvePosition数组保存曲线数据开始的位置 WIS_OBJECT_ENTRY->Position的值 */ /** curveName数组保存每条曲线的名称 */ /** curveNum变量文件中保存曲线的数目 */ int curvePosition[512]; int curveName[512][16]; int curveNum = 0; // curveNum is the identifer of the curves number fseek(wisfile, wishead->EntryOffset, SEEK_SET); // fseek(wisfile, 10+sizeof(WIS_HEAD), SEEK_SET); while(1) { WIS_OBJECT_ENTRY *objectEntry; objectEntry = (WIS_OBJECT_ENTRY *)malloc(sizeof(WIS_OBJECT_ENTRY)); fread(objectEntry, sizeof(WIS_OBJECT_ENTRY), 1, wisfile); if(objectEntry->Attribute == 0) { break; } else if(objectEntry->Attribute == 1) { // printf("通道类型: %d/n",objectEntry->Attribute); // printf("No. %d/n", curveNum); // printf("曲线对象的名称Name: %s/n", objectEntry->Name); // fprintf(outfile, "%s ", objectEntry->Name); // 打印曲线名字 // printf("对象数据体绝对偏移量Position: %d/n/n", objectEntry->Position); curvePosition[curveNum] = objectEntry->Position; int i=0; for(i=0; i<16; i++) { curveName[curveNum][i] = objectEntry->Name[i]; } curveNum++; } else { continue; } } printf("/n%s文件的曲线数目: %d/n/n", infilename, curveNum); /* ****************************************************************** */ /* Below block read every curve head information and put it in a link */ /* ****************************************************************** */ /** 创建链表 data_node 开始读取数据 */ /** num变量每条曲线的行数 */ CURVE_HEAD *first_head = NULL; float delta = 0.0; CURVE_HEAD *data_head = NULL; int num = 0; for(num = 0; num < curveNum ; num++) // num相当于第几条曲线 { data_head = (CURVE_HEAD *)malloc(sizeof(CURVE_HEAD)); data_head->position = curvePosition[num]; fseek(wisfile, curvePosition[num], SEEK_SET); WIS_CHANNEL *channel; channel = (WIS_CHANNEL *)malloc(sizeof(WIS_CHANNEL)); fread(channel, sizeof(WIS_CHANNEL), 1, wisfile); // printf("对象的单位Unit: %s/n", channel->Unit); // printf("对象的别名AliasName: %s/n", channel->AliasName); // printf("单位的别称AliasUnit: %s/n", channel->AliasUnit); // printf("对象的最小值MinVal: %f/n", channel->MinVal); // printf("对象的最大值MaxVal: %f/n", channel->MaxVal); // printf("对象维信息数NumOfDimension: %d/n", channel->NumOfDimension); // printf("维的开始值StartVal: %f/n", channel->DimInfo[0].StartVal); // printf("维的增量值Delta: %f/n", channel->DimInfo[0].Delta); data_head->pointNum = channel->DimInfo[0].MaxSamples; // 深度采样点的个数 data_head->startMD = channel->DimInfo[0].StartVal; // 起始深度 delta = channel->DimInfo[0].Delta; // 间隔点值 data_head->stopMD = channel->DimInfo[0].StartVal + delta * channel->DimInfo[0].MaxSamples; // 结束深度 int i = 0; for(i=0; i<16; i++) { data_head->Name[i] = curveName[num][i]; } for(i=0; i<8; i++) { data_head->Unit[i] = channel->Unit[i]; } printf("No.%d/t", num+1); // 曲线编号 printf("%s/t", data_head->Name); // 曲线名 printf("%s/t", data_head->Unit); // 曲线单位 printf("%d/t", data_head->pointNum); // 采样点数 printf("%9.3f/t", delta); // 采样点间隔 printf("%9.3f/t", data_head->startMD); // 起始深度 printf("%9.3f/n", data_head->stopMD); // 结束深度 data_head->next = first_head; first_head = data_head; } /* *********************************************************************** */ /* Below block is for deciding the min and max Depth, and the max line num */ /* *********************************************************************** */ // 从这里开始判断所有曲线中最小深度和最大深度值 float minDepth = 0.0; float maxDepth = 0.0; minDepth = first_head->startMD; maxDepth = first_head->stopMD; for(data_head = first_head; data_head!= NULL; data_head=data_head->next) { if(minDepth>data_head->startMD) { minDepth = data_head->startMD; } if(maxDepth < data_head->stopMD) { maxDepth = data_head->stopMD; } } // printf("最小深度值：%f/n", minDepth); // printf("最大深度值：%f/n", maxDepth); // 利用最小最大深度值和间隔点值，确定输出的采样点数 float depth[ROW]; int lineNum = 0; lineNum = (maxDepth - minDepth) / delta; // printf("采样点个数：%d/n", lineNum); int row=0; for(row=0; row<lineNum; row++) { depth[row] = minDepth + delta * row; } /* ******************************************************** */ /* Below block put every point value into a link as a curve */ /* ******************************************************** */ // 从这里开始将数据值放入曲线各自对应的链表里 CURVE_DATA *first_node = NULL; CURVE_DATA *data_node = NULL; data_head = first_head; for(num = 0; num < curveNum ; num++) { data_node = (CURVE_DATA *)malloc(sizeof(CURVE_DATA)); int i = 0; for(i=0; i<16; i++) { data_node->Name[i] = data_head->Name[i]; } for(i=0; i<8; i++) { data_node->Unit[i] = data_head->Unit[i]; } // printf("曲线名称：%s/n", data_node->Name); // printf("曲线单位：%s/n", data_node->Unit); int startLine = 0; for(row=0; row<lineNum; row++) { if( abs(depth[row]-data_head->startMD)<= delta && depth[row] >= data_head->startMD ) { startLine = row; break; } } // printf("起始点：%d/n", startLine); // fseek(wisfile, data_head->position, SEEK_SET); fseek(wisfile, data_head->position+2*sizeof(WIS_CHANNEL)+464, SEEK_SET); // printf("测试样点数%d/n", data_head->pointNum); for(row=0; row<lineNum; row++) { float data[1]; data[0] = -9999.000; if(row>=startLine && (row-startLine) < data_head->pointNum) { fread(data, sizeof(float), 1, wisfile); data_node->Depth[row] = data_head->startMD + delta * (row-startLine); data_node->Value[row] = data[0]; } else { data_node->Depth[row] = minDepth + delta * row; data_node->Value[row] = -9999.000; } // printf("读取值：%f/n", data_node->Value[row]); } data_node->next = first_node; first_node = data_node; if(data_head->next == NULL) { break; } data_head = data_head->next; // printf("下一条曲线名：%s/n", data_head->Name); } /* **************************************************** */ /* Below block print all point value for every curve */ /* The important thing is how to decide the right depth */ /* **************************************************** */ // 从这里开始写入文本 // Print Curve Name fprintf(txtfile, "DEPTH "); for(data_node=first_node; data_node!=NULL; data_node=data_node->next) { fprintf(txtfile, "%9s ", data_node->Name); } fprintf(txtfile, "/n"); // printf("/n"); // Print Curve Unit fprintf(txtfile, "M "); for(data_node=first_node; data_node!=NULL; data_node=data_node->next) { fprintf(txtfile, "%9s ", data_node->Unit); } fprintf(txtfile, "/n"); // printf("/n"); // Print curve value for(row=0; row<lineNum; row++) { if(depth[row] > maxDepth) { break; } fprintf(txtfile, "%-9.4f ", depth[row]); for(data_node=first_node; data_node!=NULL; data_node=data_node->next) { // fprintf(txtfile, "%-9.4f ", data_node->Depth[row]); fprintf(txtfile, "%9.3f ", data_node->Value[row]); } fprintf(txtfile, "/n"); } // */ // Process Over fclose(wisfile); fclose(txtfile); printf("/n%s Complete!/n", infilename); // system("PAUSE"); // Program Over return 0; }

最后，利用Python写好界面，并调用以上编译好的exe程序，这样就可以同时处理多个wis文件，界面如下：

右侧消息列表会显示每个Wis文件处理的信息！

当然转换后的结果文件会存到所选取的路径中。转换为文本后的曲线在进行处理或者运用就方便多了！

测井曲线wis文件格式转换为ASCII文本格式小软件的开发（C与Python结合开发）相关推荐

福昕PDF编辑器把PDF转换为富文本格式Word
福昕PDF编辑器把PDF转换为富文本格式Word 今天主要是一个记录我喜欢用latex写paper,但是有时候需要提交Word版本,所以我昨天刚想到我写大论文还是用latex写,暂时提交就把PDF转 ...
正则表达式实现将html文本转换为纯文本格式(将html字符串转换为纯文本方法)
Regex regex = new Regex("<.+?>", RegexOptions.IgnoreCase); string strOutput = regex. ...
android 文件流转换为zip文件_PDF文件如何转换为TXT文本？教您一招高效转换的方法...
我们见过很多的文件格式,如果说最常见的文件格式想必要属于TXT文本格式了. 大多数人都用过TXT文本格式临时记录几件事或者临时记录几句话,毕竟TXT文件是微软在操作系统中附带的一种文本格式. 如果说T ...
html转换txt文本格式方法,hthtml转txtml转换txt文本格式方法
html 转换txt文本格式,让记事本打开. 首先网页html文档直接可以使用记事本打开,如果需要转换转化为txt文本格式,大家只需要将扩展名改为.txt文本格式扩展名即可. html转化txt截图- ...
Word自定义快捷键：将粘贴内容设为纯文本格式
这两天做了一些文档的工作,文档中间复制粘贴格式不同很麻烦,仅仅靠Ctrl+V,无法实现纯文本粘贴,因为word里自带的粘贴快捷方式"Ctrl+V"默认是带原来文本格式的粘贴,如需转 ...
matlab数据转换为tecplot格式[ASCII码格式下的plt文件]
matlab数据转换为tecplot格式[ASCII码格式下的plt文件] 1 tecplot能读取的的常见ASCII格式 1.1 文件格式与文件头 1.2 zone相关的格式 2 一维数据 3 二维 ...
las格式测井曲线_邹榕，等：顺北和托甫台区块奥陶系断裂结构单元测井响应特征初探...
引用格式:邹榕,徐中祥,张晓明,等.顺北和托甫台区块奥陶系断裂结构单测井响应特征初探[J].油气藏评价与开发,2020,10(2):18-23.ZOUR, XU Z X, ZHANG X M, et ...
python 文件格式转换_python实现txt文件格式转换为arff格式
本文实例为大家分享了python实现txt文件格式转换为arff格式的具体代码,供大家参考,具体内容如下将文件读取出来的时候默认都是字符型的,所以有转换出来有点问题,但是还是可以用的. 文件要求第一 ...
日期格式化为yyyymmdd_Excel小技巧——如何将多行日期快速转换为数字文本
Hello~大家好,本来计划这一期给大家介绍一下如何根据身份证号快速获得性别信息的小技巧,但是在上一节内容如何将数字快速转换成日期格式发布后,就有朋友私信问我,可不可以再介绍一下将日期快速转换成数字文 ...
python写文件格式转换程序_python实现txt文件格式转换为arff格式
本文实例为大家分享了python实现txt文件格式转换为arff格式的具体代码,供大家参考,具体内容如下将文件读取出来的时候默认都是字符型的,所以有转换出来有点问题,但是还是可以用的. 文件要求第一 ...

最新文章

热门文章