WIS測(cè)井?dāng)?shù)據(jù)格式

1、-作者xxxx-日期xxxxWIS測(cè)井?dāng)?shù)據(jù)格式【精品文檔】測(cè)井曲線wis文件格式轉(zhuǎn)換為ASCII文本格式小軟件的開(kāi)發(fā)（C與Python結(jié)合開(kāi)發(fā)） 分類(lèi)： C/C+ 2011-03-25 12:45 627人閱讀 評(píng)論(2) 收藏 舉報(bào) 在油田開(kāi)發(fā)地質(zhì)研究工作中，測(cè)井曲線是必不可少的數(shù)據(jù)文件之一。Forward軟件是測(cè)井軟件中非常優(yōu)秀的產(chǎn)品，F(xiàn)orward中使用的是wis二進(jìn)制格式，可是很多其他類(lèi)型的軟件不能直接識(shí)別wis格式，所以需要我們將wis格式文本格式。Forward軟件有一個(gè)transfer小軟件，不能實(shí)現(xiàn)批量處理，每次只能轉(zhuǎn)換一個(gè)文件，但實(shí)際工作中通常一次要轉(zhuǎn)換上百個(gè)文件，這就有點(diǎn)不

2、方便了，因此，我在業(yè)余時(shí)間，閱讀了Forward軟件中自帶的wis格式說(shuō)明（如下，第一章內(nèi)容），利用C和Python編寫(xiě)了一個(gè)可以批量處理的小軟件。第1章 WIS測(cè)井?dāng)?shù)據(jù)文件格式WellBase平臺(tái)底層數(shù)據(jù)文件使用WIS格式，WIS文件分為文件頭、對(duì)象入口記錄和對(duì)象數(shù)據(jù)三部分。數(shù)據(jù)存放以塊為單位。WIS文件能存放三種類(lèi)型的對(duì)象，根據(jù)對(duì)象的類(lèi)型分為通道對(duì)象、表對(duì)象和流對(duì)象。通道對(duì)象一般用來(lái)存放采集數(shù)據(jù)和處理結(jié)果（如測(cè)井曲線），表對(duì)象用來(lái)存放二維表數(shù)據(jù)（如解釋結(jié)論），流對(duì)象用來(lái)存放二進(jìn)制數(shù)據(jù)塊（如解釋參數(shù)，用戶(hù)數(shù)據(jù)）。1.1 WIS文件結(jié)構(gòu)1.1.1 文件標(biāo)識(shí)WIS文件標(biāo)識(shí)從文件偏移零開(kāi)始，為10

3、個(gè)字節(jié)的字符。當(dāng)前版本的標(biāo)識(shí)為。1.1.2 文件頭結(jié)構(gòu)頭結(jié)構(gòu)緊接文件標(biāo)識(shí)。描述WIS文件的公共信息。結(jié)構(gòu)定義如下：typedef struct tagWIS_HEADWORDMachineType;WORDMaxObjectNumber;WORDObjectNumber;WORDBlockLen;DWORDEntryOffset;DWORDDataOffset;DWORDFileSize;DWORDTimeCreate;charReserved32;WIS_HEAD;偏移字節(jié)數(shù)描述02機(jī)器類(lèi)型=1 為PC； =2為SUN； =3為IBM； =4為HP。22允許記錄的最大對(duì)象數(shù)。缺省為512個(gè)，

4、該值可以在文件產(chǎn)生時(shí)給出。42當(dāng)前記錄的對(duì)象總數(shù)（包括刪除和拋棄的對(duì)象）。62塊長(zhǎng)。WIS文件對(duì)象占用的磁盤(pán)空間以塊為單位，該值指示每一數(shù)據(jù)塊的字節(jié)數(shù)。84對(duì)象入口記錄從文件開(kāi)始的偏移量。124對(duì)象數(shù)據(jù)記錄從文件開(kāi)始的偏移量。164WIS文件的字節(jié)數(shù)大小。204WIS文件產(chǎn)生的時(shí)間。2432保留字節(jié)。1.1.3 對(duì)象入口對(duì)象入口描述每個(gè)對(duì)象的公共信息，開(kāi)始位置由頭結(jié)構(gòu)給出。每個(gè)對(duì)象的描述信息前后相連。結(jié)構(gòu)定義如下：typedef struct tagWIS_OBJECT_ENTRYcharName16;longStatus;shortAttribute;shortSubAttribute;DW

5、ORDPosition;DWORDBlockNum;DWORDTimeCreate;DWORDTimeWrite;charReserved32;WIS_OBJECT_ENTRY;偏移字節(jié)數(shù) 描述016對(duì)象的名稱(chēng)，以零結(jié)尾的字符串。164對(duì)象的狀態(tài)：=0為正常； =1為拋棄； =2為刪除。202 對(duì)象的主屬性：=1為通道對(duì)象； =2為表對(duì)象； =3為流對(duì)象。222對(duì)象的子屬性，描述對(duì)應(yīng)主屬性的子屬性。244對(duì)象數(shù)據(jù)體從文件開(kāi)始處的偏移量。284對(duì)象數(shù)據(jù)體占用磁盤(pán)的塊數(shù)。324對(duì)象產(chǎn)生的時(shí)間。364對(duì)象最近修改的時(shí)間。4032保留字節(jié)。1.1.4 對(duì)象數(shù)據(jù)體對(duì)象數(shù)據(jù)體記錄各個(gè)對(duì)象的具體特性及數(shù)據(jù)。

6、根據(jù)不同的主屬性分三種類(lèi)型。對(duì)象數(shù)據(jù)體在WIS文件中的位置由對(duì)象入口指定。1.1.5 通道對(duì)象通道對(duì)象用來(lái)存放采集和計(jì)算結(jié)果數(shù)據(jù)(如測(cè)井曲線)。分為通道信息和通道數(shù)據(jù)兩部分。WIS文件將在一定時(shí)空內(nèi)對(duì)某一采集或計(jì)算的物理信息數(shù)據(jù)集統(tǒng)稱(chēng)為通道數(shù)據(jù)。通道信息描述通道數(shù)據(jù)的存放形式，分為基本信息和維信息，基本信息描述信息的基本物理含義，維信息描述信息的時(shí)空特性，可以等間隔(連續(xù))或非等間隔(離散)。最大允許有四維信息，通道信息共占用一個(gè)塊空間，結(jié)構(gòu)定義如下： typedef struct tagWIS_CHANNLEcharUnit8;charAliasName16;charAliasUnit16;

7、WORDRepCode;WORDCodeLen;floatMinVal;floatMaxVal;WORDReserved;WORDNumOfDimension;WIS_CHANNEL_DIMENSION DimInfo4;WIS_CHANNEL;偏移字節(jié)數(shù)描述08對(duì)象的單位，以零結(jié)尾的字符串。816對(duì)象的別名，以零結(jié)尾的字符串。2416單位的別名，以零結(jié)尾的字符串。402對(duì)象數(shù)據(jù)類(lèi)型，參見(jiàn)。422數(shù)據(jù)類(lèi)型的長(zhǎng)度。444對(duì)象的最小值(測(cè)井曲線缺省左刻度值)。484對(duì)象的最大值(測(cè)井曲線缺省右刻度值)。522保留字節(jié)。542對(duì)象維信息數(shù)。564*56對(duì)象維信息。通道維信息結(jié)構(gòu)定義如下：typede

8、f struct tagWIS_CHANNLE_DIMENSIONcharName8;charUnit8;charAliasName16;floatStartVal;floatDelta;DWORDSamples;DWORDMaxSamples;DWORDSize;WORDRepCode;WORDReserved;WIS_CHANNEL_DIMENSION;偏移字節(jié)數(shù)描述08維的名稱(chēng)，以零結(jié)尾的字符串。88維的單位，以零結(jié)尾的字符串。1616維的別名，以零結(jié)尾的字符串。324維的開(kāi)始值。364維的采集或計(jì)算增量。對(duì)于離散數(shù)據(jù)，該值為0，數(shù)據(jù)中記錄該維的值。404維的數(shù)據(jù)采樣點(diǎn)數(shù)。如果該值為0，

9、采樣點(diǎn)數(shù)為可變值，數(shù)據(jù)中記錄該值。對(duì)于第一維數(shù)據(jù)，該值不能為0。444維的數(shù)據(jù)采樣最大點(diǎn)數(shù)。該值僅當(dāng)采樣點(diǎn)數(shù)信息為0（可變采樣點(diǎn)）時(shí)有效，該維信息在數(shù)據(jù)中所占用的字節(jié)數(shù)通過(guò)該值計(jì)算。48 4該維上每一采樣點(diǎn)所占用的字節(jié)數(shù)。522維的數(shù)據(jù)類(lèi)型，參見(jiàn)。542保留字節(jié)。通道數(shù)據(jù)從通道描述信息的下一塊開(kāi)始。下面為一個(gè)包含深度和時(shí)間維的物理信息數(shù)據(jù)體的存放順序。第一維為深度，第二維為時(shí)間。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代表深度值，當(dāng)深度維信息結(jié)構(gòu)中的Delta為零時(shí)，

10、記錄此值。N2代表當(dāng)前深度點(diǎn)上的時(shí)間采樣點(diǎn)數(shù)，當(dāng)時(shí)間維信息結(jié)構(gòu)中的采樣點(diǎn)數(shù)為零時(shí)，記錄此值。B1，B2， ，BN代表時(shí)間值，當(dāng)時(shí)間維信息結(jié)構(gòu)中的Delta為零時(shí)，記錄此值。X1，X2， ，N代表物理信息的值。1.1.6 表對(duì)象表對(duì)象用來(lái)存放二維表數(shù)據(jù)，分為表信息和表數(shù)據(jù)體兩個(gè)部分。表信息由不同的表項(xiàng)組成，每一表項(xiàng)稱(chēng)為字段。表信息結(jié)構(gòu)定義如下：typedef struct tagWIS_TABLEDWORDRecordCount;DWORDFieldCount;WIS_TABLE_FIELD *pField;WIS_TABLE; 偏移字節(jié)數(shù)描述04表的記錄數(shù)。44表的字段數(shù)。84指向字段信息結(jié)構(gòu)

11、的指針。字段信息結(jié)構(gòu)定義如下：typedef struct tagWIS_TABLE_FIELDcharName32;WORDRepCode;WORDLength;DWORDReserved;WIS_TABLE_FIELD;偏移字節(jié)數(shù)描述032 字段的名稱(chēng)，以零結(jié)尾的字符串。322 字段值的浮點(diǎn)類(lèi)型，參見(jiàn)。342 字段值的長(zhǎng)度。表數(shù)據(jù)體(記錄)從表信息記錄的下一塊開(kāi)始。1.1.7 流對(duì)象流對(duì)象用來(lái)存放二進(jìn)制數(shù)據(jù)塊。開(kāi)始為4個(gè)字節(jié)的無(wú)符號(hào)長(zhǎng)整形數(shù)，代表數(shù)據(jù)流的長(zhǎng)度。接著為該流的二進(jìn)制值。第二章 WIS測(cè)井?dāng)?shù)據(jù)格式轉(zhuǎn)換文本格式C程序代碼代碼如下：view plaincopy to clipboar

12、dprint?1. /*2. *Copyright(c)2011,SunYunqiang3. *Allrightsreserved.4. *Software:wis2txt_win5. *Version:1.06. *Author:SunYunqiang7. *EnglishName:AlexSun8. *ReleaseDate:18/03/20119. */10. 11. /*/12. /*ThisProgramisforconvertingwisformatfiletotxtformatfile*/13. /*Command:args1args2*/14. /*Description:*/

13、15. /*wis2txt_win.exe:theexecutablefileascommand*/16. /*args1:thefirstargumentthatisthewisfilename*/17. /*args2:thesecondargumentthatisthetxtfilename*/18. /*/19. 20. #include 21. #include 22. #include 23. #include 24. #include 25. #include 26. #include 27. 28. #defineROW65536 29. #defineCOL512 30. 3

14、1. 32. /*/33. /*Belowblocksarethestructsofwisfile*/34. /*/35. typedefstructtagWIS_HEAD36. 37. WORDMachineType;/0-PC1-SUN2-IBM3-HP 38. WORDMaxObjectNumber;39. WORDObjectNumber;40. WORDBlockLen;41. DWORDEntryOffset;42. DWORDDataOffset;43. DWORDFileSize;44. time_tTimeCreate;45. charReserved32;46. WIS_H

15、EAD;47. 48. typedefstructtagWIS_OBJECT_ENTRY49. 50. charName16;51. longStatus;52. shortAttribute;/1-通道對(duì)象2-表對(duì)象3-流對(duì)象 53. shortSubAttribute;/1-曲線對(duì)象2-波形對(duì)象3-地層測(cè)試對(duì)象4-時(shí)深時(shí)對(duì)象 54. DWORDPosition;/對(duì)象數(shù)據(jù)體從文件開(kāi)始處的偏移量 55. DWORDBlockNum;56. time_tTimeCreate;57. time_tTimeWrite;58. charReserved32;59. WIS_OBJECT_ENTRY;

16、60. 61. typedefstructtagWIS_CHANNEL_DIMENSION62. 63. charName8;64. charUnit8;65. charAliasName16;66. floatStartVal;67. floatDelta;68. DWORDSamples;69. DWORDMaxSamples;70. DWORDSize;71. WORDRepCode;72. WORDReserved;73. WIS_CHANNEL_DIMENSION;74. 75. typedefstructtagWIS_CHANNEL76. 77. charUnit8;78. cha

17、rAliasName16;79. charAliasUnit16;80. WORDRepCode;81. WORDCodeLen;82. floatMinVal;83. floatMaxVal;84. WORDReserved;85. WORDNumOfDimension;86. WIS_CHANNEL_DIMENSIONDimInfo4;87. WIS_CHANNEL;88. 89. typedefstructtagWIS_STREAM90. 91. DWORDLength;92. DWORDOffset;93. WIS_STREAM;94. 95. typedefstructtagWIS_

18、TABLE_FIELD96. 97. charName32;98. WORDRepCode;99. WORDLength;100. DWORDReserved;101. WIS_TABLE_FIELD;102. 103. typedefstructtagWIS_TABLE104. 105. DWORDRecordCount;106. DWORDFieldCount;107. WIS_TABLE_FIELD*pField;108. WIS_TABLE;109. 110. 111. typedefstructtagWIS_TABLE_DEFAULT_FILED112. 113. charName1

19、6;114. charAlias16;115. charUnit8;116. charType8;117. WORDLength;118. WORDCount;119. charDefVal6412;120. WIS_TABLE_DEFAULT_FIELD;121. 122. 123. typedefstructtagWIS_DEFAULT_TABLE124. 125. charName16;126. charAlias16;127. charAttrb8;128. DWORDFieldCount;129. WIS_TABLE_DEFAULT_FIELD*pField;130. WIS_DEF

20、AULT_TABLE;131. 132. 133. 134. /*曲線的信息結(jié)構(gòu)體*/135. typedefstructtagCURVE_HEAD136. 137. charName16;/曲線的名字 138. charUnit8;/曲線的單位 139. floatstartMD;/曲線的起始深度 140. floatstopMD;/曲線的結(jié)束深度 141. intpointNum;/數(shù)據(jù)樣點(diǎn)數(shù) 142. intposition;/數(shù)據(jù)位置 143. structtagCURVE_HEAD*next;144. CURVE_HEAD;145. 146. /*曲線的數(shù)據(jù)結(jié)構(gòu)體*/147. typ

21、edefstructtagCURVE_DATA148. 149. charName16;/曲線的名字 150. charUnit8;/曲線的單位 151. floatDepthROW;/曲線的深度值 152. floatValueROW;/曲線的值 153. structtagCURVE_DATA*next;154. CURVE_DATA;155. 156. 157. 158. 159. /*/160. /*BelowblockistheMainprogramofprocessingwisformatfile*/161. /*/162. 163. intmain(intargc,char*ar

22、gv)164. 165. FILE*wisfile;/輸入文件 166. FILE*txtfile;/輸出文件 167. char*infilename=argv1;168. char*outfilename=argv2;169. /char*infilename=demoa.wis; 170. /char*outfilename=demoa.txt; 171. 172. if(argc!=3)173. 174. printf(ERROR:Commandparametersisnotdefinedcorrectly!/n);175. exit(1);176. 177. 178. 179. wi

23、sfile=fopen(infilename,rb);180. if(wisfile=NULL)181. 182. printf(ERROR:Readfileterminated!/n);183. exit(1);184. 185. txtfile=fopen(outfilename,w);186. if(txtfile=NULL)187. 188. printf(ERROR:Writefileterminated!/n);189. exit(1);190. 191. 192. /*Processwisfile,Readitandtransferittotext*/193. 194. /*wi

24、s文件標(biāo)識(shí)符從文件偏移0開(kāi)始，為10個(gè)字節(jié)的字符*/195. charwisid10;196. fread(wisid,10,1,wisfile);197. /printf(文件標(biāo)識(shí)符:%s/n,wisid); 198. fprintf(txtfile,ASCIIWIS%c%c%c/n,wisid4,wisid5,wisid6);199. 200. /*頭文件緊接文件標(biāo)識(shí)*/201. WIS_HEAD*wishead;202. wishead=(WIS_HEAD*)malloc(sizeof(WIS_HEAD);203. fseek(wisfile,10,SEEK_SET);204. frea

25、d(wishead,sizeof(WIS_HEAD),1,wisfile);205. /printf(對(duì)象入口記錄從文件開(kāi)始的偏移量EntryOffset:%d/n,wishead-EntryOffset); 206. /printf(對(duì)象數(shù)據(jù)記錄從文件開(kāi)始的偏移量DataOffset:%d/n,wishead-DataOffset); 207. /printf(當(dāng)前記錄的對(duì)象總數(shù)ObjectNumber:%d/n,wishead-ObjectNumber); 208. 209. 210. /*對(duì)象入口，位置由頭結(jié)構(gòu)中EntryOffset參數(shù)指定最多512條曲線*/211. /*curveP

26、osition數(shù)組保存曲線數(shù)據(jù)開(kāi)始的位置WIS_OBJECT_ENTRY-Position的值*/212. /*curveName數(shù)組保存每條曲線的名稱(chēng)*/213. /*curveNum變量文件中保存曲線的數(shù)目*/214. intcurvePosition512;215. intcurveName51216;216. intcurveNum=0;/curveNumistheidentiferofthecurvesnumber 217. fseek(wisfile,wishead-EntryOffset,SEEK_SET);218. /fseek(wisfile,10+sizeof(WIS_HE

27、AD),SEEK_SET); 219. while(1)220. 221. WIS_OBJECT_ENTRY*objectEntry;222. objectEntry=(WIS_OBJECT_ENTRY*)malloc(sizeof(WIS_OBJECT_ENTRY);223. fread(objectEntry,sizeof(WIS_OBJECT_ENTRY),1,wisfile);224. 225. if(objectEntry-Attribute=0)226. 227. break;228. 229. elseif(objectEntry-Attribute=1)230. 231. /p

28、rintf(通道類(lèi)型:%d/n,objectEntry-Attribute); 232. /printf(No.%d/n,curveNum); 233. /printf(曲線對(duì)象的名稱(chēng)Name:%s/n,objectEntry-Name); 234. /fprintf(outfile,%s,objectEntry-Name);/打印曲線名字 235. /printf(對(duì)象數(shù)據(jù)體絕對(duì)偏移量Position:%d/n/n,objectEntry-Position); 236. curvePositioncurveNum=objectEntry-Position;237. inti=0;238. f

29、or(i=0;iNamei;241. 242. 243. curveNum+;244. 245. else246. 247. continue;248. 249. 250. 251. printf(/n%s文件的曲線數(shù)目:%d/n/n,infilename,curveNum);252. 253. 254. 255. /*/256. /*Belowblockreadeverycurveheadinformationandputitinalink*/257. /*/258. 259. /*創(chuàng)建鏈表data_node開(kāi)始讀取數(shù)據(jù)*/260. /*num變量每條曲線的行數(shù)*/261. CURVE_HE

30、AD*first_head=NULL;262. floatdelta=0.0;263. 264. CURVE_HEAD*data_head=NULL;265. intnum=0;266. for(num=0;numposition=curvePositionnum;271. 272. fseek(wisfile,curvePositionnum,SEEK_SET);273. WIS_CHANNEL*channel;274. channel=(WIS_CHANNEL*)malloc(sizeof(WIS_CHANNEL);275. fread(channel,sizeof(WIS_CHANNEL

31、),1,wisfile);276. /printf(對(duì)象的單位Unit:%s/n,channel-Unit); 277. /printf(對(duì)象的別名AliasName:%s/n,channel-AliasName); 278. /printf(單位的別稱(chēng)AliasUnit:%s/n,channel-AliasUnit); 279. /printf(對(duì)象的最小值MinVal:%f/n,channel-MinVal); 280. /printf(對(duì)象的最大值MaxVal:%f/n,channel-MaxVal); 281. /printf(對(duì)象維信息數(shù)NumOfDimension:%d/n,cha

32、nnel-NumOfDimension); 282. /printf(維的開(kāi)始值StartVal:%f/n,channel-DimInfo0.StartVal); 283. /printf(維的增量值Delta:%f/n,channel-DimInfo0.Delta); 284. 285. data_head-pointNum=channel-DimInfo0.MaxSamples;/深度采樣點(diǎn)的個(gè)數(shù) 286. data_head-startMD=channel-DimInfo0.StartVal;/起始深度 287. delta=channel-DimInfo0.Delta;/間隔點(diǎn)值 28

33、8. data_head-stopMD=channel-DimInfo0.StartVal+delta*channel-DimInfo0.MaxSamples;/結(jié)束深度 289. 290. 291. inti=0;292. for(i=0;iNamei=curveNamenumi;295. 296. for(i=0;iUniti=channel-Uniti;299. 300. 301. printf(No.%d/t,num+1);/曲線編號(hào) 302. printf(%s/t,data_head-Name);/曲線名 303. printf(%s/t,data_head-Unit);/曲線單位

34、 304. printf(%d/t,data_head-pointNum);/采樣點(diǎn)數(shù) 305. printf(%9.3f/t,delta);/采樣點(diǎn)間隔 306. printf(%9.3f/t,data_head-startMD);/起始深度 307. printf(%9.3f/n,data_head-stopMD);/結(jié)束深度 308. 309. 310. 311. data_head-next=first_head;312. first_head=data_head;313. 314. 315. 316. 317. /*/318. /*Belowblockisfordecidingthe

35、minandmaxDepth,andthemaxlinenum*/319. /*/320. 321. /從這里開(kāi)始判斷所有曲線中最小深度和最大深度值 322. floatminDepth=0.0;323. floatmaxDepth=0.0;324. minDepth=first_head-startMD;325. maxDepth=first_head-stopMD;326. for(data_head=first_head;data_head!=NULL;data_head=data_head-next)327. 328. if(minDepthdata_head-startMD)329.

36、 330. minDepth=data_head-startMD;331. 332. if(maxDepthstopMD)333. 334. maxDepth=data_head-stopMD;335. 336. 337. /printf(最小深度值：%f/n,minDepth); 338. /printf(最大深度值：%f/n,maxDepth); 339. 340. /利用最小最大深度值和間隔點(diǎn)值，確定輸出的采樣點(diǎn)數(shù) 341. floatdepthROW;342. intlineNum=0;343. lineNum=(maxDepth-minDepth)/delta;344. /print

37、f(采樣點(diǎn)個(gè)數(shù)：%d/n,lineNum); 345. 346. introw=0;347. for(row=0;rowlineNum;row+)348. 349. depthrow=minDepth+delta*row;350. 351. 352. 353. 354. /*/355. /*Belowblockputeverypointvalueintoalinkasacurve*/356. /*/357. 358. /從這里開(kāi)始將數(shù)據(jù)值放入曲線各自對(duì)應(yīng)的鏈表里 359. 360. CURVE_DATA*first_node=NULL;361. CURVE_DATA*data_node=NULL;362. 363. data_head=first_head;364. 365. for(num=0;numcurveNum;num+)366. 367. data_node=(CURVE_DATA*)malloc(sizeof(CURVE_DATA);368. inti=0;369. for(i=0;iNamei=data_head-Namei;372. 373. for(i=0;iUniti=data_head-Uniti;376. 377. /printf(