WIS测井数据格式.docx

1、WIS测井数据格式测井曲线wis文件格式转换为ASCII文本格式小软件的开发（C与Python结合开发） 分类： C/C+ 2011-03-25 12:45 627人阅读 评论(2) 收藏 举报 在油田开发地质研究工作中，测井曲线是必不可少的数据文件之一。Forward软件是测井软件中非常优秀的产品，Forward中使用的是wis二进制格式，可是很多其他类型的软件不能直接识别wis格式，所以需要我们将wis格式文本格式。Forward软件有一个transfer小软件，不能实现批量处理，每次只能转换一个文件，但实际工作中通常一次要转换上百个文件，这就有点不方便了，因此，我在业余时间，阅读了For

2、ward软件中自带的wis格式说明（如下，第一章内容），利用C和Python编写了一个可以批量处理的小软件。第1章 WIS测井数据文件格式WellBase平台底层数据文件使用WIS格式，WIS文件分为文件头、对象入口记录和对象数据三部分。数据存放以块为单位。WIS文件能存放三种类型的对象，根据对象的类型分为通道对象、表对象和流对象。通道对象一般用来存放采集数据和处理结果（如测井曲线），表对象用来存放二维表数据（如解释结论），流对象用来存放二进制数据块（如解释参数，用户数据）。1.1 WIS文件结构1.1.1 文件标识WIS文件标识从文件偏移零开始，为10个字节的字符。当前版本的标识为WIS 1

3、.0。1.1.2 文件头结构头结构紧接文件标识。描述WIS文件的公共信息。结构定义如下：typedef struct tagWIS_HEADWORD MachineType;WORD MaxObjectNumber;WORD ObjectNumber;WORD BlockLen;DWORD EntryOffset;DWORD DataOffset;DWORD FileSize;DWORD TimeCreate;char Reserved32;WIS_HEAD;偏移 字节数 描述0 2 机器类型=1 为PC； =2为SUN； =3为IBM； =4为HP。2 2 允许记录的最大对象数。缺省为512

4、个，该值可以在文件产生时给出。4 2 当前记录的对象总数（包括删除和抛弃的对象）。6 2 块长。WIS文件对象占用的磁盘空间以块为单位，该值指示每一数据块的字节数。8 4 对象入口记录从文件开始的偏移量。12 4 对象数据记录从文件开始的偏移量。16 4 WIS文件的字节数大小。20 4 WIS文件产生的时间。24 32 保留字节。1.1.3 对象入口对象入口描述每个对象的公共信息，开始位置由头结构给出。每个对象的描述信息前后相连。结构定义如下：typedef struct tagWIS_OBJECT_ENTRYchar Name16;long Status;short Attribute;s

5、hort SubAttribute;DWORD Position;DWORD BlockNum;DWORD TimeCreate;DWORD TimeWrite;char Reserved32;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 对象最近修改的时间

6、。40 32 保留字节。1.1.4 对象数据体对象数据体记录各个对象的具体特性及数据。根据不同的主属性分三种类型。对象数据体在WIS文件中的位置由对象入口指定。1.1.5 通道对象通道对象用来存放采集和计算结果数据(如测井曲线)。分为通道信息和通道数据两部分。WIS文件将在一定时空内对某一采集或计算的物理信息数据集统称为通道数据。通道信息描述通道数据的存放形式，分为基本信息和维信息，基本信息描述信息的基本物理含义，维信息描述信息的时空特性，可以等间隔(连续)或非等间隔(离散)。最大允许有四维信息，通道信息共占用一个块空间，结构定义如下： typedef struct tagWIS_CHANNL

7、Echar Unit8;char AliasName16;char AliasUnit16;WORD RepCode;WORD CodeLen;float MinVal;float MaxVal;WORD Reserved;WORD NumOfDimension;WIS_CHANNEL_DIMENSION DimInfo4;WIS_CHANNEL;偏移 字节数 描述0 8 对象的单位，以零结尾的字符串。8 16 对象的别名，以零结尾的字符串。24 16 单位的别名，以零结尾的字符串。40 2 对象数据类型，参见42 2 数据类型的长度。44 4 对象的最小值(测井曲线缺省左刻度值)。48 4

8、对象的最大值(测井曲线缺省右刻度值)。52 2 保留字节。54 2 对象维信息数。56 4*56 对象维信息。通道维信息结构定义如下：typedef struct tagWIS_CHANNLE_DIMENSIONchar Name8;char Unit8;char AliasName16;float StartVal;float Delta;DWORD Samples;DWORD MaxSamples;DWORD Size;WORD RepCode;WORD Reserved;WIS_CHANNEL_DIMENSION;偏移 字节数 描述0 8 维的名称，以零结尾的字符串。8 8 维的单位，以

9、零结尾的字符串。16 16 维的别名，以零结尾的字符串。32 4 维的开始值。36 4 维的采集或计算增量。对于离散数据，该值为0，数据中记录该维的值。40 4 维的数据采样点数。如果该值为0，采样点数为可变值，数据中记录该值。对于第一维数据，该值不能为0。44 4 维的数据采样最大点数。该值仅当采样点数信息为0（可变采样点）时有效，该维信息在数据中所占用的字节数通过该值计算。48 4 该维上每一采样点所占用的字节数。52 2 维的数据类型，参见54 2 保留字节。通道数据从通道描述信息的下一块开始。下面为一个包含深度和时间维的物理信息数据体的存放顺序。第一维为深度，第二维为时间。A1+N2+

10、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_TAB

11、LEDWORD RecordCount;DWORD FieldCount;WIS_TABLE_FIELD *pField;WIS_TABLE;偏移 字节数 描述0 4 表的记录数。4 4 表的字段数。8 4 指向字段信息结构的指针。字段信息结构定义如下：typedef struct tagWIS_TABLE_FIELDchar Name32;WORD RepCode;WORD Length;DWORD Reserved;WIS_TABLE_FIELD;偏移 字节数 描述0 32 字段的名称，以零结尾的字符串。32 2 字段值的浮点类型，参见34 2 字段值的长度。表数据体(记录)从表信息记录的

12、下一块开始。1.1.7 流对象流对象用来存放二进制数据块。开始为4个字节的无符号长整形数，代表数据流的长度。接着为该流的二进制值。第二章 WIS测井数据格式转换文本格式C程序代码代码如下：view plaincopy to clipboardprint?1./*2.*Copyright(c)2011,SunYunqiang3.*Allrightsreserved.4.*Software:wis2txt_win5.*Version:1.06.*Author:SunYunqiang7.*EnglishName:AlexSun8.*ReleaseDate:18/03/20119.*/10./*/11

13、./*ThisProgramisforconvertingwisformatfiletotxtformatfile*/12./*Command:wis2txt_win.exeargs1args2*/13./*Description:*/14./*wis2txt_win.exe:theexecutablefileascommand*/15./*args1:thefirstargumentthatisthewisfilename*/16./*args2:thesecondargumentthatisthetxtfilename*/17./*/18.#include 19.#include 20.#

14、include 21.#include 22.#include 23.#include 24.#include 25.#defineROW65536 26.#defineCOL512 27./*/28./*Belowblocksarethestructsofwisfile*/29./*/30.typedefstructtagWIS_HEAD31.32.WORDMachineType;/0-PC1-SUN2-IBM3-HP 33.WORDMaxObjectNumber;34.WORDObjectNumber;35.WORDBlockLen;36.DWORDEntryOffset;37.DWORD

15、DataOffset;38.DWORDFileSize;39.time_tTimeCreate;40.charReserved32;41.WIS_HEAD;42.typedefstructtagWIS_OBJECT_ENTRY43.44.charName16;45.longStatus;46.shortAttribute;/1-通道对象2-表对象3-流对象 47.shortSubAttribute;/1-曲线对象2-波形对象3-地层测试对象4-时深时对象 48.DWORDPosition;/对象数据体从文件开始处的偏移量 49.DWORDBlockNum;50.time_tTimeCreate

16、;51.time_tTimeWrite;52.charReserved32;53.WIS_OBJECT_ENTRY;54.typedefstructtagWIS_CHANNEL_DIMENSION55.56.charName8;57.charUnit8;58.charAliasName16;59.floatStartVal;60.floatDelta;61.DWORDSamples;62.DWORDMaxSamples;63.DWORDSize;64.WORDRepCode;65.WORDReserved;66.WIS_CHANNEL_DIMENSION;67.typedefstructtag

17、WIS_CHANNEL68.69.charUnit8;70.charAliasName16;71.charAliasUnit16;72.WORDRepCode;73.WORDCodeLen;74.floatMinVal;75.floatMaxVal;76.WORDReserved;77.WORDNumOfDimension;78.WIS_CHANNEL_DIMENSIONDimInfo4;79.WIS_CHANNEL;80.typedefstructtagWIS_STREAM81.82.DWORDLength;83.DWORDOffset;84.WIS_STREAM;85.typedefstr

18、ucttagWIS_TABLE_FIELD86.87.charName32;88.WORDRepCode;89.WORDLength;90.DWORDReserved;91.WIS_TABLE_FIELD;92.typedefstructtagWIS_TABLE93.94.DWORDRecordCount;95.DWORDFieldCount;96.WIS_TABLE_FIELD*pField;97.WIS_TABLE;98.typedefstructtagWIS_TABLE_DEFAULT_FILED99.100.charName16;101.charAlias16;102.charUnit

19、8;103.charType8;104.WORDLength;105.WORDCount;106.charDefVal6412;107.WIS_TABLE_DEFAULT_FIELD;108.typedefstructtagWIS_DEFAULT_TABLE109.110.charName16;111.charAlias16;112.charAttrb8;113.DWORDFieldCount;114.WIS_TABLE_DEFAULT_FIELD*pField;115.WIS_DEFAULT_TABLE;116./*曲线的信息结构体*/117.typedefstructtagCURVE_HE

20、AD118.119.charName16;/曲线的名字 120.charUnit8;/曲线的单位 121.floatstartMD;/曲线的起始深度 122.floatstopMD;/曲线的结束深度 123.intpointNum;/数据样点数 124.intposition;/数据位置 125.structtagCURVE_HEAD*next;126.CURVE_HEAD;127./*曲线的数据结构体*/128.typedefstructtagCURVE_DATA129.130.charName16;/曲线的名字 131.charUnit8;/曲线的单位 132.floatDepthROW;

21、/曲线的深度值 133.floatValueROW;/曲线的值 134.structtagCURVE_DATA*next;135.CURVE_DATA;136./*/137./*BelowblockistheMainprogramofprocessingwisformatfile*/138./*/139.intmain(intargc,char*argv)140.141.FILE*wisfile;/输入文件 142.FILE*txtfile;/输出文件 143.char*infilename=argv1;144.char*outfilename=argv2;145./char*infilena

22、me=demoa.wis; 146./char*outfilename=demoa.txt; 147.if(argc!=3)148.149.printf(ERROR:Commandparametersisnotdefinedcorrectly!/n);150.exit(1);151.152.wisfile=fopen(infilename,rb);153.if(wisfile=NULL)154.155.printf(ERROR:Readfileterminated!/n);156.exit(1);157.158.txtfile=fopen(outfilename,w);159.if(txtfi

23、le=NULL)160.161.printf(ERROR:Writefileterminated!/n);162.exit(1);163.164./*Processwisfile,Readitandtransferittotext*/165./*wis文件标识符从文件偏移0开始，为10个字节的字符*/166.charwisid10;167.fread(wisid,10,1,wisfile);168./printf(文件标识符:%s/n,wisid); 169.fprintf(txtfile,ASCIIWIS%c%c%c/n,wisid4,wisid5,wisid6);170./*头文件紧接文件

24、标识*/171.WIS_HEAD*wishead;172.wishead=(WIS_HEAD*)malloc(sizeof(WIS_HEAD);173.fseek(wisfile,10,SEEK_SET);174.fread(wishead,sizeof(WIS_HEAD),1,wisfile);175./printf(对象入口记录从文件开始的偏移量EntryOffset:%d/n,wishead-EntryOffset); 176./printf(对象数据记录从文件开始的偏移量DataOffset:%d/n,wishead-DataOffset); 177./printf(当前记录的对象总数O