07-第7章-测井数据处理-total.ppt

第四节 CRA程序 CRA是Atlas公司的复杂岩性分析程序。主要功能包括： 6种方法计算泥质含量（类似POR程序）； 交会图技术确定孔隙度及两种岩石成分； 计算次生孔隙度、含水饱和度、渗透率、视颗粒密度等（详见参考书11.6节）。 一、概述 第四节 CRA程序 CRA计算泥质含量的方法与POR程序基本相同。 二、计算泥质含量 第四节 CRA程序 CRA采程序设有C1、C2、C3和C4四种矿物成分（即标准四矿物：石英、方解石、白云石和石膏） ，按其在交会图上的位置，可与水点构成三个三角形，由上之下顺序称为第一、第二、第三个三角形。 资料点落入某三角形内，就认为它是那两种矿物组成的岩石。 交会图可用任意两种孔隙度测井组合（一般采用中子—密度）。 可以根据地区地质情况指定指定两种矿物，也可以通过程序自动确定。 二、计算孔隙度及矿物成分的方法原理 第四节 CRA程序 下面以中子—密度交会图和第一个三角形（石英、方解石）为例，说明计算孔隙度和矿物成分的方法。 二、计算孔隙度及矿物成分的方法原理 二、计算孔隙度及矿物成分的方法原理 含水双矿物纯岩石密度、中子响应方程为： 把上述方程组写成数学上的规格化形式，令： 则上述方程组变成： （4-4-1） （4-4-2） 解之得： 其中， （4-4-3） （4-4-4） 系数A1、B1、C1和A2、B2、C2称为三角形系数； 上述三角形系数公式形式为第一个三角形的； 可以推广到第二、第三个三角形，但需作变化？？ 第四节 CRA程序 1、计算三角形系数： 按照设定的四种矿物，与水点构成三个三角形。三次调用三角形系数子程序得到这三个三角形的系数。 三、CRA计算孔隙度和矿物成分的流程 问题：若岩石含泥质上述计算是否存在问题？ （i=1,2,3；分别表示第1、第2和第3个三角形的系数。） 第四节 CRA程序 2、对孔隙度测井资料做泥质校正： 由于方程组（4-4-1） 是按照纯岩石模型得到的，如用（4-4-4）式 计算泥质岩石孔隙度和矿物成分，就必须把孔隙度测井资料结果校正到纯岩石的响应。 三、CRA计算孔隙度和矿物成分的流程 第四节 CRA程序 2、对孔隙度测井资料做泥质校正： 如右图，根据泥质岩石体积模型 可写出密度测井响应方程为： 三、CRA计算孔隙度和矿物成分的流程 若使上式满足（4-4-1）（即纯岩石模型），必须有： 第四节 CRA程序 2、对孔隙度测井资料做泥质校正： CRA程序中实际采用经验公式进行泥质校正，即： 三、CRA计算孔隙度和矿物成分的流程 3、对每个资料点都用三个交会三角形求解： 其中， i =1,2,3；分别表示第1、第2和第3个三角形的结果。 4、选择解释结果： （1）若VM1(1)、VM2(1)、VM3(1)>0，说明资料点落入第一个三角形内，输出其结果； （2）若VM1(1)<0，说明资料点落入第一个三角形的左下方（a区）。 令：VM1(1)=0， VM2(1)/(1- VM2(1))，VM3(1)/(1- VM3(1))； （3）若VM1(1) 、VM2(1) >0，而VM3(1)<0，说明资料点落在第一个三角形上方（b区）。 令VM3(1)=0， VM2(1) = 1-VM1(1) （4）若VM2(1) <0，说明资料点落在第一个三角形右下方（c区），也可能落入第二或第三个三角形内，或落在第三个三角形右下方。 进一步判断是否落入第二个三角形中（？）； 否则，再判断是否落入第三个三角形中（？） 若VM2(3) <0，则令VM2(3) =0， VM3(3)=1- VM1(1) b a c 判断标准： 以第一个三角形为例： V1>0,V2>0,V3>0,资料点落在三角形内 V1<0，资料点落在a区 V1>0,V3<0,资料点落在b区 V2<0,资料点落在第一个三角形下方的c区，表明该点在第一个三角形无解，再去判断第二个三角形是否有解 b a c A点：X=1,Y=2.52 V1>0,V2>0,V3<0，落在b区 B点：X=4,Y=2.6 V1>0,V2>0,V3>0，落在第一个三角形内 C点： X=7,Y=2.72 V2<0,落在c区，经计算落在第二个三角形内 D点： X=12,Y=2.84 V2<0,落在c区，经计算落在第二个三角形下方 例：矿物1取砂岩，矿物2取灰岩，矿物3取白云岩 A B C D 1 2 第七章 作业 7-1、思考如何把POR及SAND2程序结合起来， 编程处理实际