測(cè)井巖相軟件使用說明手冊(cè)

巖石物理相分析軟件使用說明書二零零九年七月一、 軟件介紹及運(yùn)行環(huán)境測(cè)井相解釋軟件系統(tǒng)是在微機(jī)上利用VC和FORTRAN語言環(huán)境編寫的軟件系統(tǒng)。該系統(tǒng)軟件運(yùn)行環(huán)境：1） 硬件環(huán)境：PII或以上處理器的微機(jī)，64M或更大的內(nèi)存空間，800X600或更高分辨率的彩色顯示器。2） 軟件環(huán)境：Windows9X、2000、XP等微軟的操作系統(tǒng)，具備VisualC++編程環(huán)境。二、 軟件的功能2.1數(shù)據(jù)轉(zhuǎn)換1）實(shí)現(xiàn)LA（80）＞Forward（文本）和自定義文本三種文件格式之間的互相轉(zhuǎn)換:首先打開文件，通過數(shù)據(jù)文件的不同后綴確定讀入的是哪類測(cè)井?dāng)?shù)據(jù)格式：（LA（80）格式*.la、Forward（文本）格式*"乂七、自定義文本格式*.dat）；然后保存為目標(biāo)格式的數(shù)據(jù)文件即可，文件類型同樣根據(jù)文件名后綴確定。LA（80）為80條標(biāo)準(zhǔn)測(cè)井LA格式；Forward（文本）為Forward軟件標(biāo)準(zhǔn)文本格式；自定義文本數(shù)據(jù)格式如下：(geng18.dat，/**/內(nèi)的是注釋）geng18/*井名*/3/*曲線條數(shù)*/2040.0000,2325.0000/*起止深度*/0.1250/*采樣間隔*/END/*間隔字符*/#DEPTHGRCALAC2040.00083.4690027.41300250.867002040.12588.4710027.60300247.289002040.25095.1761927.60300239.789792040.375101.4576027.53900238.608602040.500107.4931927.43500241.652392040.625112.6471927.36700249.787802040.750116.4800027.31000260.840612040.875118.5320027.35800270.81122

2041.000118.7002027.47900276.277402041.000118.7002027.47900276.277402041.125 117.40440 27.51400 275.683012041.250 115.06059 27.49000 269.69260… … … … …2）將地質(zhì)相等數(shù)據(jù)插入井?dāng)?shù)據(jù)中：首先點(diǎn)擊文件一＞打開，打開原始井?dāng)?shù)據(jù)文件，可選擇將巖性類型插為單條和多條兩種情況，點(diǎn)擊插入新曲線一＞單條曲線，選擇巖心曲線或解釋巖性曲線，可分別將巖心分析巖性資料輸入為一條曲線；點(diǎn)擊插入新曲線一＞多條曲線可將巖心分析巖性資料輸入為多條曲線，每一類巖性輸入為一條曲線，這樣可便于成果圖輸出，最后點(diǎn)擊文件一＞保存即可。圖1插入測(cè)井?dāng)?shù)據(jù)及數(shù)據(jù)格式轉(zhuǎn)換界面2.2測(cè)井相解釋本部分是此軟件系統(tǒng)的核心部分，主要是通過對(duì)測(cè)井?dāng)?shù)據(jù)進(jìn)行測(cè)井相分析，包括地層分層、主成分及聚類分析、譜系圖分析、選取標(biāo)準(zhǔn)樣本層和貝葉斯逐步判別分析等功能模塊，建立地區(qū)測(cè)井相-巖性數(shù)值解釋模型。程序運(yùn)行主界面如下（圖2）。圖2測(cè)井相解釋對(duì)話框2.2.1測(cè)井曲線分層及特征參數(shù)提取功能說明：測(cè)井相分析是按層段進(jìn)行解釋的，因此，利用測(cè)井資料進(jìn)行測(cè)井相劃分前首先要根據(jù)測(cè)井曲線對(duì)井剖面進(jìn)行層段的劃分，此處采用對(duì)稱差斜率法實(shí)現(xiàn)地層自動(dòng)分層。輸入數(shù)據(jù)：解釋參數(shù)文件、樣本深度文件；輸出數(shù)據(jù)：SSBP.DATo操作流程：1）首先編輯解釋參數(shù)文件和樣本深度文件；a解釋參數(shù)文件解釋參數(shù)文件格式如上圖所示，第一行為處理井井名；第二行為處理起止深度；第三行以后為各類解釋參數(shù)，中間以逗號(hào)分隔；結(jié)尾標(biāo)志符為/*。各解釋參數(shù)含義及取值說明如下：

圖3解釋參數(shù)文件SHFG：計(jì)算泥質(zhì)含量VSH標(biāo)志符（隱含值為2）,SHFG=0不計(jì)算泥質(zhì)含量；SHFG=1用GR曲線計(jì)算泥質(zhì)含量；SHFG=2用SP曲線計(jì)算泥質(zhì)含量；SHFG=3用CNL曲線計(jì)算泥質(zhì)含量；SHFG=4用RT曲線計(jì)算泥質(zhì)含量；SHFG=5用NLL曲線計(jì)算泥質(zhì)含量；s:選用解釋曲線標(biāo)志符（隱含值0），s=1選用該條曲線進(jìn)行解釋；s=0不選i i i用該條曲線進(jìn)行解釋；（S/與曲線對(duì)應(yīng)關(guān)系見下表）s1s2s3s4s5s6S7s8s9s10s11s12s13曲線名RXO1DENCNLGRACSPUKTHCGRRXORTPED最小值MFMNDEMNCNMNGRMNACMNSPMNUMNKMNTHMNCGMNRXMNRTMNPEMN最大值MFMXDEMXCNMXGRMXACMXSPMXUMXKMXTHMXCGMXRXMXRTMXPEMXs14S15S16S17S18S19S20S21S22S23S24S25曲線名TPTSTSTDTREDALPALSALTAPNASNATNVSH最小值TPMNTSMNTTMNDRMNEDMNPLMNSLMNTLMNPNMNSNMNTNMNSHMN最大值TPMXTSMXTTMXDRMXEDMXPLMXSLMXTLMXPNMXSNMXTNMXSHMX各參數(shù)代碼說明如下（括號(hào)內(nèi)為程序默認(rèn)取值）：MFMX，MFMN：淺探測(cè)電阻率測(cè)井曲線RXO1的最大、最小值（100，0.1）；RXMX，RXMN：中探測(cè)電阻率測(cè)井曲線RXO的最大、最小值（100，0.1）；RTMX，RTMN:深探測(cè)電阻率測(cè)井曲線RT的最大、最小值（100，0.1）；DEMX，DEMN:密度測(cè)井曲線DEN的最大、最小值（100，0）；CNMX，CNMN:中子測(cè)井曲線CNL的最大、最小值（100，0）；GRMX，GRMN：自然伽瑪測(cè)井曲線GR的最大、最小值（100，0）；ACMX，ACMN:聲波測(cè)井曲線AC的最大、最小值（100，0）；SPMX,SPMN：自然電位測(cè)井曲線SP的最大、最小值(100,0)；UMX,UMN:鈾測(cè)井曲線U的最大、最小值(100,0)；KMX,KMN:鉀測(cè)井曲線K的最大、最小值(100,0)；THMX,THMN：釷測(cè)井曲線TH的最大、最小值(100,0)；CGMX,CGMN:無鈾自然伽瑪能譜測(cè)井曲線CGR的最大、最小值(100,0)；PEMX,PEMN:光電吸收截面指數(shù)測(cè)井曲線PED的最大、最小值(100,0)；SHMX,SHMN:泥質(zhì)含量曲線VSH的最大、最小值(100,0)；b樣本層深度文件圖4樣本深度文件樣本深度文件格式如上圖所示，其第一列為樣本段起始深度，第二列為樣本段終止深度，第三列為巖相類型編號(hào)，中間以逗號(hào)分隔。點(diǎn)擊分層處理彈出讀取參數(shù)卡文件對(duì)話框(如圖5)和樣本層深度文件對(duì)話框(如圖6),讀入解釋參數(shù)文件和樣本深度文件，完成計(jì)算后提示分層處理完成。圖5讀取參數(shù)卡文件對(duì)話框

圖6樣本層深度文件3）輸出文件SSBP.DAT自動(dòng)保存在工作目錄下。2.2.2主成分及聚類分析功能說明：對(duì)上面模塊輸出的參數(shù)進(jìn)行主成分分析，將各樣本參數(shù)信息綜合到少量幾個(gè)主成分中，根據(jù)主成分對(duì)樣本巖相進(jìn)行聚類分析，將相似樣本聚集到一起。輸入數(shù)據(jù)：SSBP.DAT；輸出數(shù)據(jù)：SORR.DAT（聚類分析結(jié)果文件）。操作流程：點(diǎn)擊聚類分析彈出讀取數(shù)據(jù)文件對(duì)話框（如圖7），選擇上一步運(yùn)算結(jié)果文件，計(jì)算完成后提示聚類分析完成。圖7讀取數(shù)據(jù)文件對(duì)話框2.2.3譜系圖分析功能說明：對(duì)上面模塊輸出的聚類分析結(jié)果，以譜系圖的形式直觀地顯示出來，通過剔除少量樣本將相同巖相類型的樣本層選擇為同一類測(cè)井相，為選取有效建模樣本和快速建立測(cè)井相與巖相對(duì)應(yīng)關(guān)系提供了極大的方便。輸入數(shù)據(jù)：SORR.DAT（聚類分析結(jié)果文件）；輸出數(shù)據(jù)：SOR.DAT（譜系圖分析結(jié)果文件）、EFTL.DAT（測(cè)井相類型與巖相類型對(duì)應(yīng)關(guān)系文件）。操作流程：1）點(diǎn)擊繪制譜系圖彈出選擇數(shù)據(jù)文件對(duì)話框（如圖8），將上一步生成聚類分析結(jié)果文件讀入，顯示聚類分析譜系圖（如圖9）。圖8選擇數(shù)據(jù)文件對(duì)話框圖9譜系圖顯示結(jié)果2）點(diǎn)擊譜系圖一＞調(diào)整聚類結(jié)果（如圖10），圖中標(biāo)題部分提示了總樣本個(gè)數(shù)，左邊窗口中兩列數(shù)字分別是樣本號(hào)和樣本對(duì)應(yīng)巖相類型編號(hào)，右邊窗口中的縱線表示聚類

線，可通過鼠標(biāo)方便的進(jìn)行選擇及編輯;圖10譜系圖調(diào)整聚類分析結(jié)果模塊對(duì)聚類結(jié)果的編輯：右鍵單擊右邊的分類線，在彈出菜單中單擊選擇項(xiàng)，表示標(biāo)記了該分類結(jié)果，分類線變成紅色，同時(shí)左側(cè)相應(yīng)樣本也被綠色線框在一起（如圖11）；如果認(rèn)為選擇的聚類線不合理，可以單擊取消項(xiàng)，取消對(duì)該聚類線的標(biāo)記。對(duì)樣本的編輯：如果同一類的眾多樣本中，只有少數(shù)樣本與其它樣本巖相類型不同，同樣可以通過編輯對(duì)不合理的樣本進(jìn)行剔除操作，該樣本數(shù)字變成灰色（如圖11）；如果不小心剔除了準(zhǔn)確的樣本，可以通過同樣的步驟重新選取樣本。圖11譜系圖調(diào)整示例圖3）對(duì)樣本聚類結(jié)果調(diào)整好后，可以通過單擊譜系圖一＞統(tǒng)計(jì)聚類結(jié)果，統(tǒng)計(jì)樣本總數(shù)、聚類選擇的有效樣本數(shù)和聚類的有效分類數(shù)；

4） 完成聚類分析的調(diào)整編輯后，即可得到相應(yīng)合理的聚類結(jié)果，并可得到用于下一步建模的有效樣本層編號(hào)及測(cè)井相與巖相類型的對(duì)應(yīng)關(guān)系，單擊譜系圖一＞輸出聚類結(jié)果，將上訴信息保存為SOR.DAT。5） 根據(jù)SOR.DAT，建立測(cè)井相類型與巖相類型對(duì)應(yīng)關(guān)系，將其保存為EFTL.DAT，其格式如圖12所示，第一行為測(cè)井相分類數(shù)，以后幾行均為測(cè)井相與巖相類型對(duì)應(yīng)關(guān)系，前面數(shù)字為測(cè)井相類型編號(hào)，后面數(shù)字為巖相類型編號(hào)。圖12測(cè)井相類型與巖相類型對(duì)應(yīng)關(guān)系示例2.2.4選取建模有效樣本層功能說明：根據(jù)上面模塊輸出的調(diào)整好聚類分析結(jié)果及建立的測(cè)井相類型與巖相類型對(duì)應(yīng)關(guān)系，將用于建立模型的有效樣本層提取出來保存為NSBP.DAT。輸入數(shù)據(jù)：SOR.DAT、SSBP.DAT；輸出數(shù)據(jù)：NSBP.DAT（建模有效樣本層文件）。操作流程：點(diǎn)擊建模一＞選取標(biāo)準(zhǔn)層彈出讀取數(shù)據(jù)文件對(duì)話框（如圖13），讀入SSBP.DAT文件，根據(jù)上一步輸出的SOR.DAT和EFTL.DAT選取可用于建模的有效樣本層。圖13讀取數(shù)據(jù)文件對(duì)話框2.2.5建立地區(qū)巖相解釋數(shù)值模型功能說明：根據(jù)上面模塊輸出的有效建模樣本數(shù)據(jù)，應(yīng)用貝葉斯逐步判別分析的方法建立地區(qū)沉積微相解釋數(shù)值模型，建立模型的同時(shí)輸出模型對(duì)建模數(shù)據(jù)的判別分析結(jié)果，可以依此評(píng)價(jià)模型的正確性。輸入數(shù)據(jù)：NSBP.DAT；輸出數(shù)據(jù)：DATA.DAT（各類測(cè)井相判別數(shù)值模型文件）、BSSAD.DAT（建模過程明細(xì)表）。操作流程：點(diǎn)擊建模一＞建立模型彈出讀取建模數(shù)據(jù)文件對(duì)話框（如圖14），讀入建模數(shù)據(jù)文件，運(yùn)算后提示模型建立完成。圖14讀取建模數(shù)據(jù)對(duì)話框