用VB6.0實(shí)現(xiàn)抗壓強(qiáng)度成果的格式化輸出

芮學(xué)來

(淮南礦業(yè)集團(tuán) 地質(zhì)勘探工程處物探隊(duì)，安徽 淮南 232052）

0 引言

利用數(shù)字測井資料計(jì)算巖石抗壓強(qiáng)度是煤田測井所具有的一個(gè)基本功能。但現(xiàn)有測井處理軟件如CLGIS、LOGANALYST等均只能提供抗壓強(qiáng)度成果曲線，不能提供地質(zhì)部門需要的分層抗壓強(qiáng)度成果表。該表提供了各分層的巖石名稱、深度、厚度及抗壓強(qiáng)度的最小值、平均值和最大值。以往這項(xiàng)工作只能完全由人工來完成，其做法是：先將測井抗壓強(qiáng)度曲線打印出來，然后用鉛筆、直尺分別量出每層抗壓強(qiáng)度的最小值，最大值，并根據(jù)曲線中心線大致估算出該層的平均值。最后統(tǒng)計(jì)成表，用WORD或EXCELL打印出來。因此這是一項(xiàng)耗時(shí)耗力的工作，不僅效率極低，而且計(jì)算結(jié)果不準(zhǔn)確，迫切需要一種高效的方法取代之。

為此，我們通過對測井抗壓強(qiáng)度數(shù)據(jù)文件和測井?dāng)?shù)據(jù)庫柱狀文件格式的解析，用VB6.0語言開發(fā)了這個(gè)抗壓強(qiáng)度成果格式化輸出程序。該程序在輸入抗壓強(qiáng)度文件名、柱狀文件名及需要輸出的深度范圍，即可自動(dòng)計(jì)算并輸出各分層的巖石名稱、深度、厚度、抗壓強(qiáng)度最小值、平均值及最大值等。該程序不僅操作方便，容易使用，而且大大提高了工作效率。

1 原始數(shù)據(jù)的獲取

我們現(xiàn)在使用的測井處理軟件有LOGSYS、CLGIS和LOGANALYST。它們對單孔數(shù)字測井資料的處理流程大致相同。即：1）建立測井?dāng)?shù)據(jù)庫，根據(jù)設(shè)置不同，庫內(nèi)可存放10條、20條、30條或更多條測井曲線；2）將現(xiàn)場采集到的測井?dāng)?shù)據(jù)讀入數(shù)據(jù)庫，不同測井參數(shù)分別以相應(yīng)測井曲線的形式存放；3）對測井?dāng)?shù)據(jù)進(jìn)行預(yù)處理；4）選用至少三種測井參數(shù)對煤巖層進(jìn)行分層定厚，并生成測井柱狀文件；5）根據(jù)設(shè)計(jì)要求計(jì)算相關(guān)物性參數(shù)如巖石抗壓強(qiáng)度等，計(jì)算結(jié)果也以相應(yīng)參數(shù)曲線的形式存放于數(shù)據(jù)庫中；6）輸出相關(guān)成果圖、表、曲線等。

為自動(dòng)計(jì)算并輸出巖石抗壓強(qiáng)度成果表，我們需要巖石抗壓強(qiáng)度原始數(shù)據(jù)和鉆孔柱狀文件。在上述測井資料處理過程中所得抗壓強(qiáng)度參數(shù)是以抗壓強(qiáng)度曲線的形式存放于測井?dāng)?shù)據(jù)庫中，該數(shù)據(jù)不能被直接引用。但可通過編程或使用LOGANALYST軟件的數(shù)據(jù)導(dǎo)出功能將其從數(shù)據(jù)庫中導(dǎo)出，形成一個(gè)文本文件。該文件由11行文件頭及隨后由深度和抗壓強(qiáng)度數(shù)據(jù)組成的文件體，抗壓強(qiáng)度數(shù)據(jù)以5cm為采樣間隔，1米20個(gè)點(diǎn)。如下所示：

井號,ZJFJ

起始深度,330.00

終止深度,965.10

采樣間隔,0.050

空數(shù)據(jù)表示,

曲線數(shù),1

數(shù)據(jù)點(diǎn)數(shù),12703

曲線1代號,KYQD

曲線1名稱,

曲線1單位,

深度,KYQD

330.00 ,171.326920

330.05 ,173.547119

330.10 ,173.522171

330.15 ,173.459915

330.20 ,173.385437

330.25 ,175.607376

……

測井柱狀文件存放測井分層定厚結(jié)果，它以數(shù)據(jù)庫名稱重復(fù)為名，以GOL為后綴，存放于數(shù)據(jù)庫目錄內(nèi)。該文件內(nèi)容可以被直接使用。文件內(nèi)容與格式如下所示：

抗壓強(qiáng)度數(shù)據(jù)文件是我們進(jìn)行分層抗壓強(qiáng)度計(jì)算的源數(shù)據(jù)，而具體怎么分層，每層的名稱、每層有多厚，則要用到測井柱狀文件中的數(shù)據(jù)。二才缺一不可。

2 程序的編制

根據(jù)程序功能要求及輸入輸出關(guān)系，設(shè)計(jì)程序框圖。然后用Visual Basic 6.0強(qiáng)大的編程功能，對一些常用的功能如對話框、提示框等，進(jìn)行程序的編寫，由于篇幅關(guān)系，程序代碼此略。當(dāng)代碼調(diào)試完畢，用VB6.0“文件”菜單中的“生成程序”功能，生成一個(gè)可獨(dú)立執(zhí)行程序，取名“分層抗壓強(qiáng)度格式化輸出程序”。以后即可隨時(shí)運(yùn)行。

3 程序的應(yīng)用

本程序使用簡單方便，具體使用步驟如下：

1）對目的層段進(jìn)行分層定厚，生成測井柱狀文件XXXXXXXX.GOL，XXXX為數(shù)據(jù)庫文件名；

2）利用測井處理軟件計(jì)算并生成抗壓強(qiáng)度源數(shù)據(jù)，然后使用LOGANALYST軟件的“文件導(dǎo)出”功能，將源數(shù)據(jù)以文本格式導(dǎo)出，取名kyqd.txt。

為方便操作，我們一般將上述兩個(gè)文件置于當(dāng)前數(shù)據(jù)庫目錄。

3）運(yùn)行程序，界面如圖1所示：

圖1

分別輸入鉆孔名稱、開始深度、結(jié)束深度、抗壓強(qiáng)度文件名及路徑、GOL文件名及路徑。程序默認(rèn)輸出對象是打印機(jī)。也可輸出到文件，或同時(shí)輸出到兩個(gè)對象。當(dāng)點(diǎn)選輸出到文件時(shí)，程序會(huì)提示輸入存放結(jié)果的文件名，路徑默認(rèn)與抗壓強(qiáng)度數(shù)據(jù)源文件相同。

4）單擊“打印成果表”，程序自動(dòng)運(yùn)行。打印輸出如下圖。對于輸出到文件內(nèi)容，則可用Excell表格軟件進(jìn)行隨意編輯和排版，更方便地質(zhì)部門使用。

表1 張集礦中央?yún)^(qū)新副檢孔分層抗壓強(qiáng)度計(jì)算成果表

4 結(jié)語

單孔抗壓強(qiáng)度數(shù)據(jù)以1000米深鉆孔為例，5厘米采樣間隔，1米20個(gè)點(diǎn)，1000米共有20000個(gè)數(shù)據(jù)。這么大的數(shù)據(jù)量用人工去分層計(jì)算抗壓強(qiáng)度平均值、最小值、最大值，幾乎是不可能。所以只有采取如圖1所示那樣大致估算的方法，即使如此，僅一個(gè)鉆孔也需要一個(gè)人2至3天的時(shí)間才能做完。而本程序可以根據(jù)測井柱狀圖中的分層結(jié)果，自動(dòng)從抗壓強(qiáng)度源文件中取數(shù)、計(jì)算并輸出成果表，只需幾分鐘。若不打印只輸出到文件，則瞬間完成，極大提高了工作效率。

[1]尹乾.中文Visual Basic 6.0實(shí)用教程[M].北京希望電子出版社.