xca演算法
1. qtnyjl.sys 是個什麼文件
MODFastTracker、StarTrekker、NoiseTracker(等等)音樂模塊文件；Microsoft多計劃電子表格；Amiga/PC磁軌文件MOVQuickTimeforWindows電影MP2第二層MPEG音頻文件MP3第三層MPEG音頻文件MPAMPEG相關文件，MIME「mpeg類型」MPE，MPEG，MPGMPEG動畫文件MPPMicrosoft工程文件；CAD繪圖文件格式MPRMicrosoftFoxPro菜單（已編譯）MRIMRI掃描文件MSA魔術陰影檔案MSDLManchester的場景描述語言MSGMicrosoft郵件消息MSIWindows安裝器包MSNMicrosoft網路文檔；DescentMission文件MSPMicrosoftPaint（畫圖）點陣圖文件；WindowsInstaller路徑文件MSTWindows安裝器傳輸文件MTMMulti跟蹤器音樂模塊（MOD）文件MULUltima在線MUS音樂MUS10Mus10聲音MVBMicrosoft多媒體查看器文件MWPLotusWordPro97SmartMaster文件NANNanoscope文件(RawGrayscale)NAPNAP元文件NCBMicrosoftDeveloperStudio文件NCDNorton改變目錄NCFNetWare命令文件；LotusNotes內部剪切板NDO3D低多邊形建模器，NendonetCDF網路公用數據表單NFF中性文件格式NFTNetObjectFusion模板文件NILNorton游標庫文件（EasyIcons-兼容）NISTNISTSphere聲音NLBOracle7數據NLMNetWare可裝載模塊NLS用於本地化的國家語言支持文件（例如，Uniscape）NLUNortonLiveUpdatee-mail觸發器文件NODNetObjectFusion文件NSFLotusNotes資料庫NSONetObjectFusion文檔文件NSTNoiseTracker音樂模塊（MOD）文件NS2LotusNotes資料庫（第二版）NTFLotusNotes資料庫模板NTXCA-Clipper索引文件NWCNoteworthyComposer歌曲文件NWSMicrosoftOutlookExpress新聞消息（MIMERFC822）O01台風聲音文件OBDMicrosoftOffice活頁夾OBJ對象文件OBZMicrosoftOffice活頁夾向導OCXMicrosoft對象鏈接與嵌入定製控制項ODSMicrosoftOutlookExpress郵箱文件OFF3D網狀物對象文件格式OFNMicrosoftOfficeFileNew文件OFTMicrosoftOutlook模板OKTOktalyzer音樂模塊（MOD）文件OLBOLE對象庫OLEOLE對象OOGL面向對象圖形庫OPL組織者編程語言源文件——Psion/SymbianOPOOPL輸出可執行文件OPTMicrosoftDeveloperStudio文件OPXOPL擴展DLL（動態鏈接庫）ORAOracle7配置文件ORCOracle7腳本文件ORGLotusOrganizer文件OR2LotusOrganizer2文件OR3LotusOrganizer97文件OSSMicrosoftOffice查找文件OSTMicrosoftExchange/Outlook離線文件OTLSuperNoteTab模板文件OUTC語言輸出文件（工程設計器）文件P10TektronixPlot10繪圖文件P65PageMaker6.5文件P7CDigitalID文件（MIME）PABMicrosoft個人地址簿PACSBStudioⅡ包PAKQuakeWAD文件PAL壓縮文件PARTGo！Zilla部分下載文件PASPascal源代碼PATDataCADHatch模式文件；CorelDRAW模式；高級GravisUltrasound/Forte技術；碎片文件PBDPowerBuilder動態庫，作為本地DLL的一個替代物PBFTurtleBeach的Pinnacle銀行文件PBKMicrosoftPhoneBook（電話簿）PBL用於訶owerBuilder開發環境中的PowerBuilder動態庫PBM可導出點陣圖PBRPowerBuilder資源PCDKodakPhoto-CD映像；P-Code編譯器測試腳本，由Microsoft測試與MicrosoftVisual測試PCEMapsEudora郵箱名字的DOS文件名PCLHewlett-Packard列印機控制語言文件（列印機備用點陣圖）PCM聲音文件格式；OKIMSM6376合成晶元PCM格式PCPSymantecLiveUpdatePro文件PCSPICS動畫文件PCTMacintoshPICT籜?募?PCXZsoftPC畫筆點陣圖PDB3ComPalmPilot資料庫文件PDD可以用PaintShopPro或其他圖像處理軟體打開的圖形圖像PDFAdobeAcrobat可導出文檔格式文件（可用Web瀏覽器顯示）；Microsoft系統管理伺服器包定義文件；NetWare列印機定義文件PDPBroderbund的PrintShopDeluxe文件PDQPatton&PattonFlowerchartingPDQLite文件PDS攝影圖像文件（該文件格式的來源不清楚）PFAladdin系統對私人文件進行加密的文件PFA類型1字體（ASCⅡ）PFB類型1字體（二進制）PFCPF組件PFM列印機字體尺度PGD良好隱私（PrettyGoodPrivacy，PGP）虛擬磁碟文件PGLHP繪圖儀繪圖文件PGM可輸出灰度圖（點陣圖）PGP用良好隱私（PGP）演算法加密文件PH由Microsoft幫助文件編譯器產生的臨時文件PHP，PHP3包含有PHP腳本的HTML網頁PHTML包含有PHP腳本的HTML網頁；由Perl分析解釋的HTMLPICPC畫圖點陣圖；Lotus圖片；MacintoshPICT繪圖PICTMacintoshPICT圖形文件PIF程序信息文件；IBMPIF繪圖文件PIGLucasArts的DarkForcesWAD文件PINEpicPinball數據文件PIX內置系統點陣圖PJMKS源完整性文件PJX，PJTMicrosoftVisualFoxPro工程文件PKGMicrosoftDeveloperStudio應用程序擴展（與DLL文件類似）PKRPGP的公用鑰匙環PLPerl程序PLG由REND386/AVRIL使用的文件格式PLIOracle7數據描述PLMDiscorderTracker2模塊PLSDisorderTracker2抽樣文件；MPEGPlayList文件（由WinAmp使用）PLTHPGL繪圖儀繪圖文件；AutoCADplot繪圖文件；Gerber標志製作軟體PM5Pagemaker5.0文件PM6Pagemaker6.0文件PNG可移植的網路圖形點陣圖；PaintShopPro瀏覽器目錄PNT，PNTGMacPaint圖形文件POGDescent2PIG文件擴展POLWindowsNT策略文件POPVisualdBASE上托文件POTMicrosoftPowerpoint模塊POV視頻射線跟蹤器暫留PP4PicturePublisher4點陣圖PPAMicrosoftPowerpoint內插器PPFTurtleBeach的Pinnacle程序文件PPM可移植的象素映射點陣圖PPPParsonPowerPublisher；SerifPagePlus桌面出版預設輸出PPSMicrosoftPowerpoint幻燈片放映PPTMicrosoftPowerpoint演示文稿PQIPowerQuest驅動器圖像文件PRC3COMPalmPiltt資源（文本或程序）文件PRELotusFreelance演示文稿PRFWindows系統文件，Macromedia導演設置文件PRGdBASEClipper和FoxPro程序源文件；WAVmaker程序PRJ3DStudio（DOS）工程文件PRN列印表格（用空格分隔的文本）；DataCADWindows列印機文件PRPOberson的Prospero數據轉換產品保存的工程文件PRSHarvardGraphicsforWindows演示文件PRT列印格式化文件；Pro/ENGINEER元件文件PRVPsiMailInternet提供者模板文件PRZLotusFreelanceGraphics97文件PSPostscript格式化文件（PostScript列印機可讀文件）點陣圖文件PSIPSIONa-Law聲音文件PSMProtrackerStudio模型格式；Epic游戲的源數據文件PSPPaintShopPro圖像文件PSTMicrosoftOutlook個人文件夾文件PTDPro/ENGINEER表格文件PTMPolytracker音樂模塊（MOD）文件PUBVenturaPublisher出版物；MicrosoftPublisher文檔PWDMicrosoftPocketWord文檔PWLWindows95口令列表文件PWPPhotoworks圖像文件（能被Photoworks瀏覽的一系列文件）PWZMicrosoftPowerpoint向導PXLMicrosoftPocketExcel電子表格PY來自Yahoo的電子消息；Python腳本文件PYCPython腳本文件QADPFQuickArt文檔QBWQuickBooksforWindows文件QDT來自QuickenUK的QuickBooks數據文件，帳目/稅/貨單程序QD3DApple的QuickDraw3D元文件格式QFLFAMILYLAWYER文檔QICMicrosoft備份文件QIFQuickTime相關圖像（MIME）；Quicken導入文件QLBQuick庫QMQualityMotion文件QRYMicrosoft查詢文件QSTQuakeSpyTab文件QT，QTMQuickTime電影QTI，QTIFQuickTime相關圖像QTPQuickTime優先文件QTSMacPICT圖像文件；QuickTime相關圖像QTXQuickTime相關圖像QWSymantecQ&AWrite程序文件QXDQuarkXPress文件RPegasus郵件資源文件RARealAudio聲音文件RAMRealAudio元文件RARRAR壓縮檔案（EugeneRoshall格式）RASSun光柵圖像點陣圖RAWRAW文件格式（點陣圖）；Raw標識的PCM數據RBH由RoboHELP維持的RBH文件，它加入到一個幫助工程文件的信息中RDF資源描述框架文件（涉及XML和元數據）RDLDescent注冊水平文件REC錄音機宏；RapidComm聲音文件REG注冊表文件REPVisualdBASE報表文件RESMicrosoftVisualC++資源文件RFT可修訂的表單文本（IBM的DCA一部分或文檔內容框架結構一部分）RGB，SGISilicon圖形RGB文件RLERun-Length編碼的點陣圖RL2Descent2注冊水平文件RMRealAudio視頻文件RMDMicrosoftRegMaid文檔RMFRichMap格式（3D游戲編輯器使用它來保存圖）RMIM1D1音樂ROM基於盒式磁帶的家庭游戲模擬器文件（來自Atari2600、Colecovision、Sega、Nintendo等盒式磁帶里的ROM完全拷貝，在兩個模擬器之間不可互修改）ROVRescueRover數據文件RPMRedHat包管理器包（用於Linux）報表文件RRSAcegameRoadRash保存的文件RSLBorland的Paradox7報表RSMWinWayResumeWriter恢復文件RTFRichText格式文檔RTKRoboHELP使用的用來模擬Windows幫助的搜索功能RTMRealTracker音樂模塊（MOD）文件RTSRealAudio的RTSL文檔；RoboHELP對復雜操作進行加速RULInstallShield使用的擴展名RVPMicrosoftScan配置文件（MIME）Rxx多卷檔案上的RAR壓縮文件（xx=1～99間的一個數字）S匯編源代碼文件S3IScreamTrackerv3設備S3MScreamTrackerv3的聲音模塊文件SAMAmi專業文檔；8位抽樣數據SAV游戲保存文件SB原始帶符號位元組（8位）數據SBKCreativeLabs的Soundfont1.0Bank文件；(Soundblaster)/EMUSonndFontv1.xBank文件SBLShockwaveFlash對象文件SC2MicrosoftSchele+7文件格式；SAS目錄（Windows95/NT、OS/2、Mac）SC3SimCity3000保存的游戲文件SCCMicrosoftSourceSafe文件SCDMatrix/ImaproSCODL幻燈片圖像；MicrosoftSchele+7SCFWindowsExplorer命令文件SCHMicrosoftSchele+1SCIScanVecInspire本地文件格式SCNTrueSpace2場景文件SCP撥號網路腳本文件SCRWindows屏幕保護；傳真圖像；腳本文件SCTSAS目錄（DOS）；ScitexCT點陣圖；MicrosoftFoxPro表單SCT01SAS目錄（UNIX）SCVScanVecCASmate本地文件格式SCXMicrosoftFoxPro表單文件SDSoundDesigner1聲音文件SD2SoundDesigner2展平文件/數據分X指令；SAS資料庫（Windows95/NT、OS/2、Mac）SDF系統數據文件格式—LegacyUnisys（Sperry）格式SDKRolandS—系列軟盤映像SDLSmartDraw庫文件SDRSmartDraw繪圖文件SDS原始Midi抽樣轉儲標准文件SDTSmartDraw模板SDV分號分隔的值文件SDWLotusWordPro圖形文件；原始帶符號的DWORD（32位）數據SDX由SDX壓縮的Midi抽樣轉儲標准文件SEA自解壓檔案（StuffltforMacintosh或其他軟體使用的文件）SEP標簽圖像文件格式（TIFF）點陣圖SESCoolEditSession文件（普通數據聲音編輯器文件）SFIRCAM聲音文件格式SF2EmuSoundfontv2.0文件；CreativeLabs的Soundfont2.0Bank文件（SoundBlaster）SFDSoundStage聲音文件數據SFISoundStage聲音文件信息SFRSonicFoundrySample資源SFWSeattle電影工程（損壞的JPEG）SFXRAR自解壓檔案SGML標准通用標簽語言SHBCorelShow演示文稿；文檔快捷文件SHG熱點點陣圖SHP3DStudio（DOS）形狀文件；被一些應用程序用於多部分交互三角形模型的3D建模SHSShellscrap文件；據載用於發送「口令盜竊者」SHTML含有伺服器端包括（SSI）的HTML文件SHWCorelShow演示文稿SIG符號文件SITMac的StuffIt檔案文件SIZOracle7配置文件SKAPGP秘鑰SKLMacromedia導演者資源文件SLPACT的保存布局擴展名SLBAutodeskSlide庫文件格式SLDAutodeskSlide文件格式SLKSymbolicLink（SYLK）電子表格SM3DataCAD標志文件SMPSamplevision格式；AdLibGold抽樣文件SNDNeXT聲音；Mac聲音資源；原始的未符號化的PCM數據；AKAIMPC系列抽樣文件SNDRSounder聲音文件SNDTSndtool聲音文件SOUSBStudioⅡ聲音SPDSpeech數據文件SPLShockwaveFlash對象；DigiTrakker抽樣SPPACKSPPack聲音抽樣SPRITEAcorn的點陣圖格式SQC結構化查詢語言（SQR）普通代碼文件SQLInFORMixSQL查詢；通常被資料庫產品用於SQL查詢（腳本、文本、二進制）的文件擴展名SQR結構化查詢語言（SQR）程序文件SSDO1SAS數據集合（UNIX）SSDSAS資料庫（DOS）SSF可用的電子表格文件STAtariST磁碟映像STLSterolithography文件STM.shtml的短後綴形式，含有一個服務端包括（SSI）的HTML文件；ScreamTrackerV2音樂模塊（MOD）文件STR屏幕保護文件STYVenturaPublisher風格表SVXAmiga8SVX聲音；互交換文件格式，8SVX/16SVSW原始帶符號字（16位）數據SWA在Macromedia導演文件（MP3文件）中的Shockwave聲音文件SWFShockwaveFlash對象SWPDataCAD交換文件SYS系統文件SYWYamahaSY系列波形文件T64Commodore64模擬器磁帶映像文件TABGuitar表文件TAR磁帶檔案TAZUNIXgzip/tape檔案TBKAsymetrixToolbook交互多媒體文件TCL用TCL/TK語言編寫的腳本TDBThumbsPlus資料庫TDDDImagine和TurboSilver射線跟蹤器使用的文件格式TEX正文文件TGATarga點陣圖TGZUNIXgzip/tap檔案文件THEMEWindows95桌面主題文件THNGraphicsWorkShopforWindows速寫TIF，TIFF標簽圖像文件格式（TIFF）點陣圖TIG虎形文件，美國政府用於分發地圖TLBOLE類型庫TLE兩線元素集合（NASA）TMPWindows臨時文件TOCEudora郵箱內容表TOLKodak照片增強器TOSAtari16/32和32/32計算機操作系統文件
2. 用C語言編程 在屏幕上顯示圍棋棋盤
－－
以怎麼不用 MS VC＋＋6 ？？
編C語言一般都用啊
還有\xCD是C語言
Ascii代碼 表示一個圖形
用循環語句啊！
等我一下 晚上寫出來 貼上哦
晚上10：00 上完課回來
呵呵 分給不給 隨便 你的問題讓我熟練了許多
分不重要了
我終於寫出來了
調試了30分鍾啊！
純粹的循環 ＋ 判斷
－－ 我就這點水平了 演算法可能不好
如果有更好的告訴我 上QQ 或發網路消息給我
網路的顯示有問題 最好
用源代碼格式化工具格式化一下 網路有搜
#include <stdio.h>
void main()
{
int n,m,x,y,temp,j=0,i=0 ;
printf("輸入棋盤的橫格數，豎格數：\n");
scanf("%d,%d",&x,&y);
printf("%d,%d\n",x,y);
n=2*x+1 ;
m=2*y+1 ;
printf("%d,%d\n",n,m);
while(j<m)
{
j++;
temp=j%2 ;
i=0 ;
if(j==1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("\xC9");
continue ;
}
if(n==i)
{
printf("\xBB\n");
continue ;
}
if(i%2!=0)
{
printf("\xCB");
continue ;
}
if(i%2==0)
{
printf("\xCD\xCD\xCD");
continue ;
}
}
continue ;
}
if(m==j)
{
while(i<n)
{
i++;
if(i==1)
{
printf("\xC8");
continue ;
}
if(n==i)
{
printf("\xBC\n");
continue ;
}
if(i%2!=0)
{
printf("\xCA");
continue ;
}
if(i%2==0)
{
printf("\xCD\xCD\xCD");
continue ;
}
}
continue ;
}
if(temp==0)
{
while(i<n)
{
i++;
if(i%2!=0)
{
printf("\xBA");
if(n==i)printf("\n");
continue ;
}
if(i%2==0)
{
printf(" ");
continue ;
}
}
continue ;
}
if(temp!=0&&j!=1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("\xCC");
continue ;
}
if(n==i)
{
printf("\xB9\n");
continue ;
}
if(i%2!=0)
{
printf("\xCE");
continue ;
}
if(i%2==0)
{
printf("\xCD\xCD\xCD");
continue ;
}
}
continue ;
}
}
}
100%編譯通過
可以畫出任意大小的棋盤
如 5*3 5*5
只要你的屏幕能夠顯示
以上是用圖像庫畫的棋盤
下面是用 __ 和 |
畫的棋盤 其實很簡單 替換一下就可以了
#include <stdio.h>
void main()
{
int n,m,x,y,temp,j=0,i=0 ;
printf("輸入棋盤的橫格數，豎格數：\n");
scanf("%d,%d",&x,&y);
printf("%d,%d\n",x,y);
n=2*x+1 ;
m=2*y+1 ;
printf("%d,%d\n",n,m);
while(j<m)
{
j++;
temp=j%2 ;
i=0 ;
if(j==1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("__");
continue ;
}
if(n==i)
{
printf("__\n");
continue ;
}
if(i%2!=0)
{
printf("__");
continue ;
}
if(i%2==0)
{
printf("__");
continue ;
}
}
continue ;
}
if(m==j)
{
while(i<n)
{
i++;
if(i==1)
{
printf("__");
continue ;
}
if(n==i)
{
printf("__\n");
continue ;
}
if(i%2!=0)
{
printf("__");
continue ;
}
if(i%2==0)
{
printf("__");
continue ;
}
}
continue ;
}
if(temp==0)
{
while(i<n)
{
i++;
if(i%2!=0)
{
printf("|");
if(n==i)printf("\n");
continue ;
}
if(i%2==0)
{
printf(" ");
continue ;
}
}
continue ;
}
if(temp!=0&&j!=1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("__");
continue ;
}
if(n==i)
{
printf("__\n");
continue ;
}
if(i%2!=0)
{
printf("__");
continue ;
}
if(i%2==0)
{
printf("__");
continue ;
}
}
continue ;
}
}
}
3. 急求高手賜教，Matlab循環問題，在線等
%通過定義全局變數來解決這個問題..
function F=myfun(x)
global P; %將P設置為全局變數, 這樣matlab就會在已有變數中搜尋之.
F=[x(1)+x(2)-1-P;
x(1)-x(2)-P];
end
---------------------------
此外, 在主函數中也修改
global P;
for P=1:10
x0 = [-6;-5];
x = fsolve(@myfun,x0) %使用fsolve 函數求解方程，options默認，等價形式還可以寫成x-fosolve('myfun',x0)
end
-----------------------------------------
測試通過:
Optimization terminated: first-order optimality is less than options.TolFun.
x =
1.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
2.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
3.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
4.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
5.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
6.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
7.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
8.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
9.5000
0.5000
Optimization terminated: first-order optimality is less than options.TolFun.
x =
10.5000
0.5000
4. 如何用C語言繪出一個圍棋盤
－－
以怎麼不用 MS VC＋＋6 ？？
編C語言一般都用啊
還有\xCD是C語言
Ascii代碼 表示一個圖形
用循環語句啊！
等我一下 晚上寫出來 貼上哦
晚上10：00 上完課回來
呵呵 分給不給 隨便 你的問題讓我熟練了許多
分不重要了
我終於寫出來了
調試了30分鍾啊！
純粹的循環 ＋ 判斷
－－ 我就這點水平了 演算法可能不好
如果有更好的告訴我 上QQ 或發網路消息給我
網路的顯示有問題 最好
用源代碼格式化工具格式化一下 網路有搜
#include <stdio.h>
void main()
{
int n,m,x,y,temp,j=0,i=0 ;
printf("輸入棋盤的橫格數，豎格數：\n");
scanf("%d,%d",&x,&y);
printf("%d,%d\n",x,y);
n=2*x+1 ;
m=2*y+1 ;
printf("%d,%d\n",n,m);
while(j<m)
{
j++;
temp=j%2 ;
i=0 ;
if(j==1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("\xC9");
continue ;
}
if(n==i)
{
printf("\xBB\n");
continue ;
}
if(i%2!=0)
{
printf("\xCB");
continue ;
}
if(i%2==0)
{
printf("\xCD\xCD\xCD");
continue ;
}
}
continue ;
}
if(m==j)
{
while(i<n)
{
i++;
if(i==1)
{
printf("\xC8");
continue ;
}
if(n==i)
{
printf("\xBC\n");
continue ;
}
if(i%2!=0)
{
printf("\xCA");
continue ;
}
if(i%2==0)
{
printf("\xCD\xCD\xCD");
continue ;
}
}
continue ;
}
if(temp==0)
{
while(i<n)
{
i++;
if(i%2!=0)
{
printf("\xBA");
if(n==i)printf("\n");
continue ;
}
if(i%2==0)
{
printf(" ");
continue ;
}
}
continue ;
}
if(temp!=0&&j!=1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("\xCC");
continue ;
}
if(n==i)
{
printf("\xB9\n");
continue ;
}
if(i%2!=0)
{
printf("\xCE");
continue ;
}
if(i%2==0)
{
printf("\xCD\xCD\xCD");
continue ;
}
}
continue ;
}
}
}
100%編譯通過
可以畫出任意大小的棋盤
如 5*3 5*5
只要你的屏幕能夠顯示
以上是用圖像庫畫的棋盤
下面是用 __ 和 |
畫的棋盤 其實很簡單 替換一下就可以了
#include <stdio.h>
void main()
{
int n,m,x,y,temp,j=0,i=0 ;
printf("輸入棋盤的橫格數，豎格數：\n");
scanf("%d,%d",&x,&y);
printf("%d,%d\n",x,y);
n=2*x+1 ;
m=2*y+1 ;
printf("%d,%d\n",n,m);
while(j<m)
{
j++;
temp=j%2 ;
i=0 ;
if(j==1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("__");
continue ;
}
if(n==i)
{
printf("__\n");
continue ;
}
if(i%2!=0)
{
printf("__");
continue ;
}
if(i%2==0)
{
printf("__");
continue ;
}
}
continue ;
}
if(m==j)
{
while(i<n)
{
i++;
if(i==1)
{
printf("__");
continue ;
}
if(n==i)
{
printf("__\n");
continue ;
}
if(i%2!=0)
{
printf("__");
continue ;
}
if(i%2==0)
{
printf("__");
continue ;
}
}
continue ;
}
if(temp==0)
{
while(i<n)
{
i++;
if(i%2!=0)
{
printf("|");
if(n==i)printf("\n");
continue ;
}
if(i%2==0)
{
printf(" ");
continue ;
}
}
continue ;
}
if(temp!=0&&j!=1)
{
while(i<n)
{
i++;
if(i==1)
{
printf("__");
continue ;
}
if(n==i)
{
printf("__\n");
continue ;
}
if(i%2!=0)
{
printf("__");
continue ;
}
if(i%2==0)
{
printf("__");
continue ;
}
}
continue ;
}
}
}
5. 岩石密度及地震波速計算
1.程序功能
程序SEIVEL可根據用戶提供的全岩化學成分和密度，或礦物成分和含量，計算結晶岩的地震波速VP和VS。在後一種情況下，可同時計算岩石密度。
2.方法原理
由全岩化學成分和岩石密度計算地震波速的方法據Simmons（1964）、Anderson和Samis（1970）。在該法中，根據實驗資料，將地震波速擬合為全岩平均原子量、CaO的質量分數和岩石密度的函數。在缺少岩石密度資料、但已知各礦物的成分及含量的情況下，可採用Duffy等（1989）和Jackson等（1990）關於常見造岩礦物端員組分的密度資料，按照礦物的相對含量加權，計算出全岩的密度和地震波速。對於含有輝石類礦物的岩石，可以採用馬鴻文等（1998，未發表資料）有關計算輝石族晶胞體積和密度的方法。通過對比計算發現，與採用全岩化學成分和岩石密度的方法相比，採用Duffy等（1989）和Jackson等（1990）的礦物端員組分密度計演算法，計算的VP和VS波速的誤差分別達0.183～0.216和0.106；而採用馬鴻文等（1998，未發表資料）計算輝石族礦物的密度，則VP和VS波速的計算誤差分別減小為0.009～0.029和0.016，精度顯著提高。
3.程序結構
結晶岩熱力學軟體
4.使用說明
（1）輸入格式
程序運行過程中，按照屏幕提示，依次提供以下參數：
Option 選擇計算方法
IFN/OFN 輸入/輸出文件名
樣品名稱和氧化物含量由輸入文件讀入。
選擇由全岩化學成分和岩石密度計算地震波速時（Option=1），原始數據輸入格式為：A6,F6.3,13F6.2。
各變數排列順序依次為：Sample（樣品號）、Densty（密度）、SiO2、TiO2、Al2O3、Cr2O3、Fe2O3、FeO、MnO、NiO、MgO、CaO、Na2O、K2O、P2O5。
選擇由礦物化學成分、含量和端員組分密度計算岩石密度和地震波速時（Option=2,3），原始數據輸入格式為：A6,I3（樣品號，礦物相數）；A3,F6.3,13F6.2。
各變數排列順序依次為：minerl（礦物代號）、Xmi（礦物含量，摩爾分數）、SiO2O、TiO2、Al2O3、Cr2O3、Fe2O3、FeO、MnO、NiO、MgO、CaO、Na2O、K2O、P2O5。
各樣品最末一行按相同格式依次排列樣品號、岩石密度（若已知，否則可輸入0.000）及全岩各氧化物含量。
每個樣品的礦物相數和每次計算的樣品個數不限。
程序中使用的礦物代號如下：
qtz石英
plg斜長石
opx斜方輝石
cpx單斜輝石
olv橄欖石
spn尖晶石
grt石榴石
礦物代號一律採用小寫字母。
（2）輸出格式
全部計算結果輸出到文件OFN中。內容包括：樣品號，密度，室溫下的VP和VS數值，計算的p=1.0～2.6GPa、t=400～1400℃下的VP和Vs值以列表形式輸出。
多個樣品的計算結果按相同格式輸出。
5.程序文本
結晶岩熱力學軟體
結晶岩熱力學軟體
30 format（A6,F6.3,13F6.2）
else if（Job.eq.2.or.Job.eq.3）then
read（3,36,ERR=120,END=100）Sample,m
36 format（A6,I3）
Densty=0.0
do i=1,m
read（3,38）minerl,Xmi,（OXD（j）,j=1,n）
38 format（A3,F6.3,13F6.2）
sum=0.0
do j=1,n
OXD（j）=Ncat（j）*OXD（j）/GFW（j）
sum=sum+OXD（j）
end do
do j=1,n
OXD（j）=OXD（j）/sum
end do
if（OXD（6）.gt.0.or.OXD（9）.gt.0）XFe=OXD（6）/（OXD（6）+OXD（9）
if（minerl.eq.』qtz』）then
Densty=Densty+2.648*Xmi
else if（Minerl.eq.』plg』）then
XAn=OXD（10）/（OXD（10）+OXD（11）
Densty=Densty+（2.617+0.145*XAn）*Xmi
else if（Minerl.eq.』opx』）then
if（Job.eq.2）then
Densty=Densty+（3.204+0.799*XFe）*Xmi
else if（Job.eq.3）then
JFe3＝』Y'
call PYNORM（n,CMP,OXD,5,JFe3,Pgfw）
Dpyx=Pgfw*4/（Avog*Volpyx（n,OXD）
DenSty=DenSty+Dpyx*Xmi
end if
else if（Minerl.eq.』cpx』）then
if（Job.eq.2）then
XJd=OXD（11）/（OXD（10）+OXD（11）
Densty=Densty+（3.280+0.377*XFe）
$*（1.0-XJd）+3.320*XJd）*Xmi
else if（Job.eq.3）then
JFe3＝』Y'
call PYNORM（n,CMP,OXD,6,JFe3,Pgfw）
Dpyx=Pgfw*4/（Avog*Volpyx（n,OXD）
Densty=Densty+Dpyx*Xmi
end if
else if（Minerl.eq.'olv'）then
Densty=Densty+（3.222+1.182*XFe）*Xmi
else if（Minerl.eq.'spn'）then
Densty=Densty+（3.582+0.676*XFe）*XMi
else if（Minerl.eq.'grt'）then
XCa=OXD（10）/（OXD（6）+OXD（9）+OXD（10）
Densty=Densty+（3.562+0.758*XFe）*（1.0-XCa）+3.595*XCa）*Xmi
end if
end do
read（3,38）Rock,DST,（OXD（j）,j=1,n）
end if
****CalCulate CaO weight fraction and mean atomic weight****************
sum=0.0
do j=1,n
sum=sum+OXD（j）
end do
CaOwtf=OXD（10）/sum
do j=1,n
OXD（j）=OXD（j）/GFW（j）
end do
Oxy=0.0
Cat=0.0
do j=1,n
Oxy=Oxy+Noxy（j）*OXD（j）
OXD（j）=Ncat（j）*OXD（j）
Cat=Cat+OXD（j）
end do
sum=Oxy+Cat
do j=1,n
OXD（j）=OXD（j）/sum
end do
Oxy=Oxy/sum
MAW=15.9996*Oxy
do j=1,n
MAW=MAW+ATW（j）*OXD（j）
end do
****Calculate Vn ＆ Vsat room temperature*******************************
Vp=-0.98+0.70*（21.0-MAW）+2.76*Densty+4.600*CaOwtf
Vs=-0.63+0.21*（21.0-MAW）+1.56*Densty+0.016*CaOwtf
write（*,70）Sample,Densty,Vp,Vs
write（4,70）Sample,Densty,Vp,Vs
70 format（/5X，』Sample：』，A6，』Density＝』，F6.3,
$』 Vp＝』，F6.3，』Vs＝』，F6.3，』at room temperature』/）
75 format（1X，』P（GPa）/T（C） 400 500 600 700 800 900』，
$』 1000 1100 1200 1300 1400』）
****Calculate Vp ＆ Vs at 1.0-2.6（GPa）and 400-1400（C）*******************
do i=1,17
do j=1,11
Vptp（i,j）=Vp+（j*100+300-20）*dVpdtp+（i-1）*dVpdpt
Vstp（i,j）=Vs+（j*100+300-20）*dVsdtp+（i-1）*dVsdpt
end do
end do
write（4,*）』Seismic velocity Vp（km/s）』
write（4,75）
do i=1,17
write（4,96）（i+9.0）/10.0,（Vptp（i,j）,j=1,11）
end do
write（4,*）』Seismic velocity Vs（km/s）』
write（4,75）
do i=1,17
write（4,96）（i+9.0）/10.0,（Vstp（i,j）,j=1,11）
end do
96 format（4X,F4.1,3X,11F6.3）
goto 25
100write（*,110）OFN
110format（/5X，』Edit』，A10，』to look over the results！』）
goto 130
120 write（*,*）』File read error,data skipped！』
130end
*************************************************************************
subrout ine PYNORM（n,CMP,OXD,IM,JFe3,Pgfw）
*For deriving pyroxene formulas from chemical analyses with anion-
*based,hydrogen-equivalent method（Jackson et al,1967;US Geol SURVEY
*PROF PAPER 575C,P.C23-31）.
*Written by MA Hongwen,Nov.1988
*China Universi ty of Geosc i ences,Be i j ing 100083
*Link: MINORM----READEL
*----Fe3cal
*CMP（13）:SiO2O Ti02O A1203 Cr2O3 Fe2O3 FeO MnO NiO
*MgO CaO Na2O K2O Li2O
*n--component numbers EQW--equivalent weights NVL--cation valances
*OXD（n）--oxide（wt%）（input）,Coef（output）
*V.04/06/97
parameter（m=13,Nox=12）
character CMP（n）*6,ELM（m）*6,JFe3
dimension EQW（m）,NVL（m）,NCT（m）,OXD（n）
call READEL（n,CMP,ELM,EQW,NVL,NCT）
do j=1,m
EQW（j）=EQW（j）/（NVL（j）*NCT（j）
end do
do j=1,m
OXD（j）=OXD（j）/EQW（j）
end do
Csum=0
do j=1,m
Csum=Csum+OXD（j）
end do
Fact=Nox/Csum
do j=1,m
OXD（j）=Fact*OXD（j）/NVL（j）
end do
if（JFe3.eq.』y』.or.JFe3.eq.』Y』）then
call Fe3cal（m,OXD,NVL,IM,Nox）
end if
Pgfw=O
do j=1,m
Pgfw=Pgfw+NVL（j）*EQW（j）*OXD（j）
end do
end
*************************************************************************
6.計算實例
吉林省汪清地區尖晶石二輝橄欖岩包體（WQ023）（馬鴻文等，1994，未發表資料）的地震波速VP和VS值計算。輸出文件中同時列出採用程序中所提供的3種計算方法所得到的計算結果。對比3種計算結果可以發現，採用馬鴻文等（1998，未發表資料）計算輝石族礦物的密度，則VP和VS速的計算誤差僅相當於採用礦物端員組分密度計算地震波速（Duffyet al.，1989；Jackson et al.，1990）誤差的1/7。因此，在岩石密度未知的情況下，建議採用方法3 （Option=3）計算岩石密度和地震波速。
輸入文件：exam81.dat
結晶岩熱力學軟體
輸出文件：exam82.dat
結晶岩熱力學軟體
結晶岩熱力學軟體
Sample:1YQ023 Density=3.328 Vp=8.194 Vs=4.523 at room temperature
（Option No.2）
Seismic velocity Vp（km/s）
結晶岩熱力學軟體
結晶岩熱力學軟體
結晶岩熱力學軟體
6. 脈沖中子測井原理
（一）碳氧比伽馬能譜測井的譜分析和數據處理
圖3-57是碳氧比能譜數據採集時序示意圖，對應每一脈沖中子重復周期設置三個數據採集時間門：非彈性門與中子發射持續時間對應，主要測量由快中子非彈性散射產生的伽馬射線；同時還記錄上一周期剩餘的俘獲輻射和活化產生的伽馬射線；為得到「凈譜」，需從總譜中扣除本底，本底門設置在非彈性門之後，測量本底的近似值；非彈性譜門測量俘獲伽馬能譜。
1.伽馬能譜的解析
測井得到的中子非彈性散射伽馬譜和俘獲伽馬譜都是由多種核素生成的混合譜，解析就是從混合譜中將每種核素的貢獻分離出來，方法和自然伽馬能譜處理類似。以中子非彈性散射伽馬譜為例，設第j種核素快中子非彈性散射截面為σj，單位體積岩石中該種元素（穩定核素的豐度為常數）的原子數為nj，它對i道計數率yi的貢獻yij應與乘積xj=σjnj成正比，即
地球物理測井
式中：αij為第j種核素對i值計數率的響應系數；yi中包含所有核素的貢獻和統計及測量誤差εi，有
地球物理測井
若道數為m，有貢獻的核素為s種，m＞s，並考慮到各道計數率精度差別很大，用加權最小二乘法得到一矩陣方程：
地球物理測井
式中：A為響應系數組成的m×s階矩陣；AT為A的轉置矩陣；W為權因子ωi組成的對角矩陣，i=1，2，…，m；X為由s種核素的xj組成的列矩陣；Y為m道計數率組成的列矩陣。
對X求解，得：
地球物理測井
令
地球物理測井
有
地球物理測井
式中，E為s×m階矩陣。
從演算法上看，E是一數字濾波器，通過它能從實測混合譜中將每種元素的貢獻提取出來。對12C、16O、28Si、40Ca分別得到xC、xO、xSi、xCa，並定義碳氧化和硅鈣比分別為
地球物理測井
xj與地層中j種元素的非彈伽馬發射率成正比，也叫產額系數。用同樣的方法對俘獲伽馬譜進行解析，可獲得xH、xCl、xSi、xCa、xFe、xS、xn等參數，它們都是相應元素的俘獲輻射產額系數。
用上述方法對全譜進行解析，充分利用已獲取的信息會遇到兩點困難：①道計數率低，統計精度不高；②元素的標准譜難以獲得，如在原油中測得碳譜，原油對伽馬射線的散射和吸收與地層不同，同樣鈣和硅元素的標准譜也不易測定。若從全譜中選定幾個特徵道區（能窗），其積分計數率將會有較好的統計精度，再用礦物的標准譜代替元素的標准譜對儀器進行刻度，會更接近地層的實際情況。
2.碳氧比能譜測井的探測深度和環境影響
（1）碳氧比能譜測井的探測深度
實驗數據如圖3-58。當源距為55.88 cm時，探測深度為21.336 cm；源距增大到68.58 cm時，探測深度增大到 28.448 cm。由此可見：①碳氧比能譜測井的探測深度較小，井的影響不能忽視；②探測深度與源距有關，用較小源距研究井眼的影響和尋求井眼環境的校正方法。
圖3-58 碳氧比能譜測井的探測深度
（2）碳氧比能譜測井的環境影響
碳氧比能譜測井儀器的源距不同，井眼和地層條件不同，探測深度也不盡相同。其探測深度一般不超過30 cm，國內外模型實驗都證明了這一點。若水侵入油層深度超過20 cm，用碳氧比很難區別油層和水層。在裸眼井中，侵入帶一般都超過這一范圍。在已射孔的套管井，除侵入影響外，還有套管和水泥環的影響，情況更為復雜。
對未射孔的套管井，為使侵入帶消失，需要等適當的時間。此時井眼中流體的類型直接影響測得的碳氧比值，而水泥環對碳氧比及硅鈣比都有影響。
3.響應方程
碳氧比能譜測井的主要用途，是在孔隙水的礦化度低、不穩定或未知的條件下，在套管井中測定地層的含油飽和度，特別是測定注水開發油層的剩餘油飽和度。在其他條件相同情況下，當含油飽和度高時，單位體積地層中碳原子數較多而氧原子數較少，或者說碳氧原子數比值較高。作為用碳氧比求含油飽和度的基礎，先計算單位體積地層中的碳和氧原子數。
（1）單位體積地層中的碳和氧原子數及其比值
設a為每立方厘米原油中碳原子的數目，b為每立方厘米岩石骨架中碳原子的數目，c為每立方厘米淡水中氧原子的數目，d為每立方厘米岩石骨架中氧原子的數目。若原油密度為0.87 g/cm3，分子式為CnH2n，可以算得：
地球物理測井
每立方厘米淡水中氧原子的數目為
地球物理測井
對純砂岩地層，岩石骨架中不含碳，b=0；而每立方厘米岩石骨架中氧原子的數目為
地球物理測井
對石灰岩地層，每立方厘米岩石骨架中碳原子的數目為
地球物理測井
氧原子的數目為
地球物理測井
純砂岩地層中，孔隙度為φ，含油飽和度為So，則每立方厘米岩石的碳原子數為
地球物理測井
每立方厘米岩石的氧原子數為
地球物理測井
碳氧原子數比的響應方程為
地球物理測井
由式（3-122）可見，給定孔隙度φ，碳氧原子數比與含油飽和度So有單值關系，由此式可繪制出關系曲線如圖3-59所示。
從式3-122和圖3-59可以看出：當孔隙度大時，曲線的斜率大，測定含油飽和度的靈敏度高；對孔隙度相同的地層，含油飽和度高時靈敏度高；孔隙度高和含油飽和度也高的地層對碳氧比測井有利，可達到較高的精度。低孔隙度和高含水地層對測井不利，得不到理想的效果。
純石灰岩地層中，碳氧原子數比為
地球物理測井
相應圖形見圖3-60。
與圖3-59相比，圖3-60的不同之處有：①當含油飽和度為零時，碳氧原子數比為0.333，比孔隙度為35%和含油飽和度高達90%的純砂岩地層還要高。②當含油飽和度達到20%時，孔隙度不同的各條曲線交於一點，將曲線簇分成兩部分。③當含油飽和度小於20%時，對應於同一含油飽和度，孔隙度大的地層碳氧原子數比值低。④當含油飽和度大於20%時，對應於同一含油飽和度，孔隙度大的地層碳氧原子數比值高。
由以上分析可知，識別岩性對碳氧比能譜測井定量解釋非常重要。
圖3-59 純砂岩碳氧原子數比（COR）與含油飽和度的關系
圖3-60 純石灰岩碳氧原子數比（COR）與含油飽和度的關系
（2）產額（系數）比和含油飽和度模型
地層的碳氧產額比為
地球物理測井
式中：A=σC/σO，即截面比。令
地球物理測井
式中：nC1、nC2為單位體積地層中岩石骨架和孔隙流體中的碳原子數；nO1、nO為單位體積地層中岩石骨架和孔隙流體中的氧原子數；BC、BO為井內流體對碳和氧測量結果的影響。
顯然，這幾個量分別與地層中骨架、油和水的相對體積以及井液中的持油率或持水率成正比，有
地球物理測井
式中：Yw為井液持水率；KC1、KC2為碳的非彈性散射伽馬產額對岩石骨架和油的相對體積的靈敏度；KO1、KO2為氧的非彈性散射伽馬產額對岩石骨架和水的相對體積的靈敏度；KC3、KO3為碳或氧的非彈性散射伽馬產額對井眼中持油或持水率的靈敏度。
在岩性和孔隙度已知的情況下，對單探測器儀器求含油飽和度或含水飽和度，除需通過刻度井確定式（3-125）的六個系數外，還需測定持水率，或用實驗方法測定井液影響校正曲線。但是，對雙探測器儀器，可利用長、短源距探測器探測范圍的差別（圖3-61）來補償井液的影響。
雙探測器儀器的解釋模型是一組聯立方程。長、短源距探測器的產額比分別為
地球物理測井
解方程組可得Sw和Yw。
（3）其他產額比和岩性、孔隙度、泥質和礦化度響應
地層碳氧比主要反映含油飽和度，也可稱碳氧化能譜測井的含油飽和度響應，簡稱飽和度響應或含油飽和度指數。類似的比值有四個。
圖3-61 長、短源距探測器探測范圍示意圖
岩性指數：
地球物理測井
純碳酸鹽岩岩性指數近於零，純砂岩岩性指數近於一。因受套管外水泥環的影響，即使是純砂岩，測出的岩性指數也小於一。岩性指數幾乎不受孔隙度、含油飽和度和地層水礦化度的影響。用硅鈣非彈性散射伽馬產額或俘獲伽馬產額比，都能指示岩性，並可用以校正碳酸鹽岩的碳氧化。
孔隙度指數：
地球物理測井
式中各個元素的產額由俘獲伽馬譜求出，孔隙度指數可定性指示孔隙度的大小。
泥質指數：
地球物理測井
式中各個元素的產額由俘獲伽馬譜求出。
在裸眼井中，泥質指數從零到大於1；而對套管井，該指數可達1.5～2.5。
礦化度指數：
地球物理測井
在有利條件下，這一俘獲伽馬產額比可定性指示地層水礦化度。
（4）各元素產額曲線的比較和監測剩餘油飽和度
圖3-62是石灰岩地層，從深度A～B為油層，B以下為水層，裸眼完井，井徑為16 cm；井中充滿原油，井眼內油水界面在深度D處。從xC、xO、xCl和xH四條曲線上，都能將這兩個界面分出來；鈣的非彈性散射和俘獲伽馬產額xCa均與零線接近，而非彈性散射鈣產額高，正確地指示出岩性為石灰岩；地層和井筒中氯的影響，，使C～D和D以下井段xCl升高，xH和鈣俘獲產額xCa降低；鐵的非彈性散射和俘獲產額都近於零。對元素產額曲線的變化特點有深刻理解之後，各種比值曲線和由比值導出的飽和度曲線的變化規律和解釋方法也就不難理解了。
圖3-62 元素產額曲線圖
若用SOl和SO2分別表示在裸眼井用電測井確定的原始含油飽和度和用碳氧比測井測出的剩餘油飽和度，比較這兩條曲線就可觀察到原油采出的程度和油水界面的變化。
（二）脈沖中子孔隙度測井
脈沖中子孔隙度測井是用同位素中子源的中子孔隙度測井的替代方法，比較典型的儀器是APS測井儀。
1.APS陣列脈沖中子測井儀
圖3-63為APS陣列脈沖中子測井儀示意圖。脈沖中子發生器發射14 MeV的中子，由五個3H計數管記錄超熱中子和熱中子。儀器有貼井壁裝置，中子計數管的背後用碳化硼屏蔽起來，以消減井液的影響。最大的一對超熱中子探測器，由短源距和長源距超熱中子計數管組成，用與補償中子測井類似的計數率比值法求地層的中子孔隙度。在上述兩個計數器之間有三個計數器組成陣列，離源較近的一對是超熱中子探測器，離源較遠的一個是熱中子探測器。用成對超熱中子探測器測量：①與時間無關的超熱中子計數率，高解析度薄層超熱中子測井曲線；②與時間無關的超熱中子計數率，即中子脈沖間隔中的超熱中子計數率時間分布，其衰減常數是快中子慢化時間的量度，與地層含氫指數相關。矩陣中的熱中子探測器，測量熱中子計數率的時間分布，求地層的熱中子宏觀截面Σ或熱中子壽命τ。
2.中子慢化時間與孔隙度的關系
圖3-64為不同孔隙度石灰岩超熱中子計數率衰減曲線。可以看出，孔隙度大的地層計數率衰減快，孔隙度小的地層衰減慢。圖3-65 給出慢化時間的倒數與孔隙度的關系，石灰岩、白雲岩和砂岩三種不同岩性數據點偏離不明顯，即對岩性不敏感。實驗表明，用比值法求出的超熱中子孔隙度受岩性影響較大。這是因為，能量較高的快中子最初的一二次碰撞所佔時間非常短，對慢化時間貢獻很小，慢化時間主要是由能量已降低的中子與氫核的彈性碰撞決定的；用計數率比值法求孔隙度，通過測量中子慢化長度來求地層孔隙度，最初幾次碰撞對中子慢化長度影響很大，因而受岩性影響較大。
（三）熱中子壽命測井
熱中子壽命測井，也稱熱中子衰減時間測井。用脈沖中子源向地層發射能量為14 MeV的中子，測量熱中子或俘獲伽馬計數率隨時間的衰減，算出地層的熱中子宏觀俘獲截面或壽命。在地層水礦化度高時，可求出地層含水飽和度。
圖3-63 APS結構示意圖
1.岩石的熱中子壽命和宏觀俘獲截面
熱中子壽命τ是指熱中子從產生的瞬時起到被俘獲的時刻止經過的平均時間。由計算可知，它等於原有的熱中子已有63.2%被俘獲而剩下的還有36.8%所經歷的時間。在常遇地層中，熱中子壽命τ主要與含氯量有關。熱中子壽命τ與宏觀俘獲截面Σ的關系為
地球物理測井
式中：v為熱中子速度，cm/s。
圖3-64 石灰岩超熱中子計數率衰減曲線
圖3-65 慢化時間與孔隙度的關系
熱中子速度與環境溫度有關，即
地球物理測井
式中：T為絕對溫度。
若熱中子壽命τ以μs為單位，並將25℃時的熱中子速度2.2×105 cm/s代入式（3-131），有
地球物理測井
測井時，通常選用10-3cm-1作為宏觀俘獲截面的單位，記作cu，於是有
地球物理測井
單一化合物的宏觀俘獲截面可用式（3-135）計算，即
地球物理測井
式中：ρ為密度，g/cm3；ni為化合物分子中第i種原子的個數；σi為第i種原子核的微觀俘獲截面；m為相對分子質量。
純岩石的熱中子宏觀俘獲截面為
地球物理測井
式中：Σma、Σw、Σh分別為岩石骨架、地層水和烴的熱中子宏觀俘獲截面。
當地層含泥質時，式（3-136）變成：
地球物理測井
式中：Vsh、Σsh分別為泥質的相對體積和熱中子宏觀俘獲截面。
表3-15給出幾種礦物的宏觀俘獲截面和壽命值。
從式（3-136）、式（3-137）和表3-14中數據可見：①高礦化度地層水熱中子宏觀俘獲截面比石英、白雲石和方解石等孔隙性岩石骨架礦物大一個數量級，是淡水或原油俘獲截面的2～3倍，因而一般儲層的宏觀俘獲截面主要決定於高礦化度地層水的相對體積。②高礦化度地層水的熱中子宏觀俘獲截面和壽命與原油有明顯區別，因而用中子壽命測井可測定含水飽和度。③地層中熱中子俘獲截面非常大的某些元素，如硼對中子壽命測井有嚴重的影響。④地層骨架礦物俘獲截面與孔隙流體有明顯區別，中子壽命測井對孔隙度敏感。⑤粘土礦物的俘獲截面大，泥質含量對中子壽命測井有較大影響。
表3-14 不同礦物的熱中子宏觀俘獲截面和熱中子壽命
2.熱中子壽命和宏觀俘獲截面的測定
本章第三節曾給出用擴散方程表示的中子數守衡定律，即
地球物理測井
式中：v為中子速度；Φ為中子通量。
測定熱中子壽命時，開始計數的時間比中子發射時間要滯後一些，熱中子產生項S已為零。通過選擇源距，擴散項的影響也可減小，必要時可做適當校正。這樣，式（3-138）簡化為
地球物理測井
積分此式，得：
地球物理測井
式中：τ為中子壽命；Φ0為初值。
因中子或伽馬計數率N與中子通量成正比，用N代替Φ後（3-140）仍正確。測井時，測量的計數率Nt包括三部分射線源的貢獻：①井內介質對熱中子計數率或俘獲伽馬計數率的貢獻N1；②地層對熱中子或俘獲伽馬計數率的貢獻N2；③井內介質和地層生成的穩定的背景值N3。這部分按式（3-140）隨時間衰減。用公式表示：
地球物理測井
或
地球物理測井
式中：N01、N02和N3是常數。
圖3-66是按式3-142計算得到的關系曲線，用以說明中子壽命測井計數率衰減曲線的組成和各分量的特點。計算時取∑1=52.5 cu和∑2=12 cu。與實測曲線相比，圖中未顯示統計漲落，而其他特點是相同的。在設計數據採集時序時，應考慮圖中顯示的這些特點。中子壽命測井的主要用途是求地層的含水飽和度。由式（3-137）可得到含水飽和度：
圖3-66 中子壽命測井原理圖
圖3-67 石灰岩熱中子計數率衰減曲線
地球物理測井
式中：Σ為測井值；Σma、Σh和Σsh分別為骨架、烴和泥質常數；Σw為地層水的宏觀俘獲截面，對原狀地層Σw是常數，而對注水開發油田它是變數；Vsh為泥質體積含量；φ為孔隙度。
孔隙度不同時，衰減曲線的斜率不同。圖3-67是不同孔隙度的石灰岩地層的熱中子計數率衰減曲線，這說明求准孔隙度對用中子壽命測井確定含水飽和度是很重要的。
定量解釋可信系數c應大於0.5，計算公式為
地球物理測井
中子壽命測井還可以監測油水界面、測定可動流體相對體積和剩餘油飽和度變化。