测井解释4－砂泥岩解释5－测井解释软件

《测井解释4－砂泥岩解释5－测井解释软件》由会员分享，可在线阅读，更多相关《测井解释4－砂泥岩解释5－测井解释软件（38页珍藏版）》请在装配图网上搜索。

1、单击以编辑母版标题样式,,单击以编辑母版文本样式,,第二级,,第三级,,第四级,,第五级,,,,*,第五节 测井解释软件,1,,POR程序的解释原理,POR程序是从美国Atlas公司引进的单孔隙度测井泥质砂岩分析程序,简单实用，输入的测井曲线数目少，在地质情况比较简单的情况下得到比较好的解释结果， 在国内普遍使用或针对地区具体条件作了改进使用。,一、POR程序的特点,二、POR程序的解释原理,2,,对测井资料进行定量解释之前，首先要计算出V,sh.,,该程序最多采用5种测井值计算。,,5,种测井值排列顺序,GR 、SP 、RT 、CNL 、NLL（中子寿命）,上式中的i =1、2、3、4、5,

2、五条,曲线采用一个公式,1、计算V,sh,3,,SH,i,: 第i条曲线相对值,,SHFG： 测井值求V,sh,的标识符,,最终的 V,sh,= min(V,sh1,,V,sh2,,…,V,sh5,）,2、求孔隙度,采用单矿物含水泥质岩石模型计算,F，在实际计算中只进行泥质校正，,PFG,为求,F的标识符。,PFG=1，,F= F,D,-V,sh,,F,Dsh,,PFG=2，,F=（ F,S,/CP）-V,sh,,F,Dsh,,PFG=3，,F= F,N,-V,sh,,F,Nsh,当SH>SHCT(,泥岩截止值,) 为泥岩,求,泥岩孔隙度,为:,F= F,,(1-,V,sh,)

3、 用此式把泥岩与砂岩孔隙度区分开。,4,,3、求Sw,Swop：求Sw的标识符,,Swop=2,，由阿尔奇公式求Sw： a=1 n=2,,m=1.87+0.019/,F， F,>0.1 取 m=2.1 ，m>4 取m=4,Swop=1,， 用simandoux公式的简化形式求Sw ：,Swop=3,，由阿尔奇公式求Sw ： a=0.62 、m=2.15 、n=2,5,,4. 求k 采用Timur公式,5．求其他地质参数,（1）、,F,w、,F,xo,F,w=,F,,Sw,F,xo =,F,,Sxo,（2）、 经验公式求S,hr,.,,S,hr,=SRHM,,(1-S,

4、W,),,SRHM: 为残余油气饱和度和含油气饱和度地区经验系数，隐含为0.5,6,,（3） 、冲洗带残余油气相对体积及质量m,hr.,,,V,hr,=S,hr,,F,,m,hr.=,V,hr,,h,,,当,,h,准确可靠时，可判断油气界面,,油层 V,hr,=m,hr,因为体积,,质量,,气层 V,hr,>> m,hr,V,hr,=m,hr,（4）累计孔隙度（PF）和累计油气厚度（HF）,F,i_,——第i个点的,F,DH,——测井曲线采样间隔（常为0.125m或0.1m）,7,,见P244图4-50,PF、HF表示从某一深度开始累计得到的纯孔隙厚度和纯油气厚度，在解释成果图上，通

5、常在某些深度位置上用短线表示，每相邻短线之间累计孔隙厚度或累计油气厚度为1m或1ft,如果井段短线越多，说明地层孔隙越发育或油气越多，如处理井段为N，该井控制面积为S，则处理井段油气体积Vh=NS,（5）、求BULK（出砂指数）,,指导采油作业的 一个参数,,BULK =13400*,r,b,/,D,t,2,D,t——,s/ft,,BULK：出砂指数（10,6,lb/in=7.04,10,8,Kg/m,2,）,8,,当,BULK >3,,正常求产方式下采油不出砂,,,BULK <3,, 正常求产方式下出砂，应减小油咀生产,可不出砂或少出砂。,,BULK是指采油作业的一个参数,1< BULK

6、 <10,至少一种孔隙度测井、GR 、RT、CAL、SP，真正的POR程序有42个输入参数，15个输出参数。,见附录P566或原程序,三．输入资料,9,,2．地层特性,四 、成果显示及实例,测井数据处理成果图是一张随深度变化的地质参数曲线图，需要时也可采用数据表的形式，数据处理成果图由以下几个部分组成：,打印深度标记和解释结论,在砂泥岩地层中，通常用V,sh,和k曲线的变化说明。V,sh,线性刻度放在左侧第一道，k采用对数刻度。,3、油气分析,S,w,V,hr,m,hr,线性刻度 、,,V,hr,、,,m,hr,取决于参数形式提供油气密度的准确性，POR程序中显示这两条曲线。,1、深度标识区

7、,10,,5．地层体积分析（岩性分析）,,位于第四道：Vsh、Vma、,F,4、孔隙度分析（也称流体分析）,,位于第三道,11,,P227 图4-44a、b POR,程序的流程图,演示过程：用POR程序处理实际资料,P229 图POR数字处理成果图,12,,SAND2 程序,1、SAND2是一种改进的泥质砂岩分析程序，用多种方法求V,sh ，,据具体情况选用适当的数值作为采用值。,2、在考虑泥质和油气影响的情况下，用中子密度交会图求,F,和Vsh。,3、还可采用一种经验方法求粉砂指数,一、SAND2程序特点,13,,4、当有CAL曲线，程序对井径扩大的井段进行自动处理。,绝大部分的参数符

8、号与POR相同，只有少许不同（5个），输出参数有19个，15个与POR相同。,二、输入资料和输出参数,1．输入资料： 必须有,2．输入输出参数,见原程序,14,,三．解释方法原理,1．求V,sh,,（1）用GR算出的为V,sh,,,（2）选,CNL、NLL、SP、RT中的一条用GR的公式求出的为SH,。,（3）计算SHP,,当Neut=0时，用SP求SHP,当Neut=1时，用Neu求SHP,15,,（4）合理确定V,sh,,①、,V,sh,=min(V,sh,，SHP),,②、SHDF=0或100（SHDF代表泥质含量差值）,2．求,F,和Vsh,,SAND2程序在考虑泥质和油气对中子密度测

9、井影响的情况下，用迭代法计算,F,和Vsh。,（1）计算公式,,①、对密度和中子进行泥质校正求,F,的公式与POR相同,16,,②、求油气密度,中子则是 PORN=(Neu-Lcor-0.5Nsh(SH))*0.01,,测井值%,骨架值%,泥岩%,小数,SH=Vsh,SRH：冲洗带残余油气饱和度,17,,③、计算地层孔隙度和Vsh,对于含油气泥质砂岩的冲洗带：,经变换得含水泥质砂岩模型：,18,,得方程组的解为：,其中A,1,、 A,2,、B,1,、 B,2,、 C,1,、 C,2,、 D为：,19,,当,F,hr,=0 即地层不含油气，这里用的泥岩是相当于泥岩点。,,因砂岩骨架点，泥岩点和水

10、点构成一个表示砂泥岩特性的三角形，而系数A,1,,B,1,,C,1,和A,2,,B,2,,C,2,决定了三角形中一个点（x,y）的孔隙度和泥质，故这些系数在SAND2中称为交会三角形系数。并有一个确定这些系数的子程序MAT.,,A,1,=XXM(1) B,1,=XXM(2) C,1,=XXM(3),20,,A,2,=YM（1） B,2,=YM(2) C,2,=YM(3),,Ф,xo,=POR V,sh,=SH Ф,hr=,GAS,,Ф有效=Ф,xo,+Ф,hr,（2）、进行油气迭代校正,1)、V,sh,=V,sh,(GR) Ф,N,、Ф,D,只经泥质校正,3）、

11、用交会三角形系数求Ф,x0,,V,sh,,,,Ф,有效,=Ф,x0,+Ф,hr,2）、赋初值S,hr,=0.5 GAS=0 ρ,h0,=1,,然后计算DHYO、DHYG,,而ρ,h,=[MAX(DHYO,DHYG)+ ρ,h0,]/2,21,,4）求Sw,,Simandoux公式、阿尔齐中A、M不同，两种选择,5）、,Sw>0.7,,地层无油气或含油气很少,,不作油气校正,,,GASO=GAS ，令,r,h,=1,6）、,Sw<0.7,,再计算S,hr,=SRHM(1-Sw)，,F,hr,=S,hr,F,,,GAS=,F,hr,,7）、当,|GAS-GASO（=上次的,F,hr,）|<0

12、.005,(,根据地区资料决定)，则油气校正结束。,当,|GAS-GASO|,,0.005,，则应继续作油气校正,此时:GAS=(GAS+GASO）/2，重新计算X、Y，进行油气校正。,22,,重复上述步骤直至｜GAS-GASO|<0.005为止,,当不需要作油气校正（ＳＷ>0.7或,,｜GAS-GASO|<0.005）调整不合理的结果,再计算,DENA=（1-,F,-Vsh）,r,ma,+,r,f,f,+V,sh,r,sh,,F,NA,=（,1-,F,-Vsh）,F,N,ma,,+,F,N,f,f,+,F,N,sh,V,sh,23,,３．计算泥质的粉砂指数（用经验公式计算）,,P(x,y)

13、为测井资料点,24,,Z随资料点P 的坐标而变，同时又与交会图几个特征点的坐标有关，故可设想它是一个与粉砂指数有关的综合参数，采用统计方法确定这一关系，假设ＳＩ与Ｚ的关系为：ＳＩ＝a,1,+a,2,Z+a,3,Z,2,25,,砂岩骨架点:A=1，Z=0，SI=1得a,1,=1,,粘土点:A=(D-D,clay,)/(D-,Ф,nclay,)，Z=1，SI=0得,,a,2,=1+a,3,,泥岩点令Z=0.5（位于泥岩线中点，SI=0.354则a,3,=0.5833,,计算经验公式 ∴SI=1-1.5833Z+0.5833Z,2,4、求Vclay和重新计算Sw,Simandoux公式,Vclay

14、=Vsh（1-SI）,R,w1,、R,mf1,为地层温度下R,w,、R,mf,，求S,x0,的公式与S,w,形式同Rt=Rxo、.R,w1,=R,mf1,26,,实例演示,5、求K,K=1.861,10,7,,F,5.795,/ Vsh,F：为小数,Vsh,：为小数,SAND2程序同样计算了,F,w,、,F,f,、,F,xo,、 V,hr,、m,hr,、 PF、 HF 。,四．输出成果,,与POR程序相同，分五道,,与POR程序不同，对油气和泥质都进行迭代校正，用,F,N——,r,b,计算,F。,27,,CRA处理方法,一、CRA的功能,1、用 六种方法求,V,sh,2、用交会图技术求Ф及

15、两种岩石成分,28,,3、求,二、输入资料、输出参数,,至少，两种孔隙度资料、Ｒt、一种泥质指示测井资料，若有其他资料效果更好。,,输出参数见书,三、方法原理,,１、求Vsh用六种方法求,,Vsh＝min（Vsh,１,，Vsh,２,，Vsh,３,，……Vsh,５,）略,29,,２．求孔隙度和岩性成分,,C1,C2,C3,C4为四种矿物成分,C,1,C,2,C,3,C,4,F,N,r,b,a,b,c,点落在哪个三角形就由那两种矿物构成，（请插入图）,,最终据地质情况确定矿物的成分和属性，但最多不超过四种矿物,,以第一个三角形为例，讲解求Ф和岩性的原理,30,,设孔隙度为Ф ，矿物为Ｃ1和Ｃ2的相

16、对体积Ｖc1, Ｖc2,相应的骨架,r,c1,、,r,c2,可以求出其解：,31,,CRA中通过调用LITH子程序求V,1,V,2,V,3,，而LITH又通过调用子程序COEF来计算三角形系数。,具体步骤：（1）、求三角形系数,,（2）、对Ф进行泥质校正,,ρ,bc,=(ρ,b,-0.8V,swh,) ρ,sh,)/(1-0.8V,sh,）,,（形式类同）Ф,NC,=,,△t,C,=,（3）、对每个资料点都用三个交会三角形求解,,VM,1,(I)=TA,1,(I)*x+TB,1,(I)*y+TC1(I),,VM,2,(I)=TA,2,(I)*x+TB,2,(I)*y+TC2(I),,VM,3

17、,(I)=100- VM,1,(I)- VM,2,(I),,用循环方法对每个资料点求解,32,,3、选择解结果,,作为解释结果只能是其中某一个三角形的解，程序从第一个三角形起自动挑选合理的解释参数。,,例如第一个三角形，据以下解释结果或转到别的三角形判断,,①、VM,1,(1), VM,2,(1), VM,3,(1)>0资料点在第一个三角形内，以此为它的解释结果,,②、VM,1,(1)<0,(即Ф<0)资料点落在第一个三角形的左下方a区,,此时调整解释结果VM,1,(1)=0,,,VM,2,(1)=,33,,VM,2,(1)=,VM,1,(1)仍为原来计算的数值,,③、VM,1,(1)>0 而

18、VM,3,(1)<0资料点落在第一个三角形的上方仍由第一个三角形得，调整解释结果VM,2,(1)=100 -VM,2,(1),,④、VM,2,(1)<0，资料点落在另一个三角形的下方，即可能在第二，第三个三角形内，该点在第一个三角形内无解。,,第二，第三个三角形的检查与上相似，只是在第三个三角形时VM,3,(3)<0时，仍由第三个三角形的解，此时VM,3,(3)=100- VM,1,(3)，VM,2,(3）=0,34,,⑥求,Ф,和C,1,、C,2,、C,3,、C,4,,,以上按纯岩石模型计算的，需恢复成泥质双矿物：,,Ф=Ф(1-SH)=C,0,(1-SH), Ф<0.1%,,则,Ф=0.1

19、%,Ф>1-SH,则,Ф=1-SH,,C,1,=C,1,(1-SH)，C,2,、C,3,、C,4,形式类同,,如SH>SHCT 则,Ф=Ф(1-SH),⑤、求出资料点所在地层的矿物成分的相对体积,,当用△,1,求体积时 C,0,=VM,1,(1)， C,1,=VM,2,(1)，C,2,=VM,3,(1)，C3=0,当用△,3,求体积时 C,0,=VM,3,(1)， C,1,=0 ，C,2,=0，C3=VM,3,(3),当用△,2,求体积时 C,0,=VM,2,(1)，C,1,=0 ，C,2,=VM,2,(2)，C3= C,0,=VM,3,(2),35,,3、求其他地质参数,,⑴求S,w,和S,

20、xo,与POR方法相同,,⑵,Ф,次（缝）,,Ф,2,=Ф,次,=,Ф-Ф,s,/cp,,Ф,次,<0，,Ф,2,=0,Ф,次,>0.8,Ф则Ф,次=,0.8Ф,,（3）、求K用Timur公式,,（5）求PF和HF,,（6）求ρ,GA,（视颗粒密度）,,△t,GA,（视颗粒时差）,四、成果显示,,与其他分析程序基本相同,,五 NCRA求Ф和两种矿物成分,,利用交会图(ρ,b,-Ф,N,)或(Ф,N,-,△t)迭代求孔隙度，然后用点至直线间的距离求两种矿物的岩石体积。,36,,实例演示,37,,总结,掌握测井解释软件的基本结构、每个程序的特点、功能，能够根据测井资料，正确的建立解释模型，选择程序所需的解释参数，对测井资料进行处理，直到得到满意的结果为止。,作业,1、POR程序的输出参数有哪些？必须有哪些测井曲线才能POR程序？,2、SAND2程序与POR程序在解释原理方面有何不同？,3、CRA、SAND2、POR程序各适用于何种岩性剖面？,38,,