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，计算公式为
地球物理测井
中子寿命测井还可以监测油水界面、测定可动流体相对体积和剩余油饱和度变化。