测井实验报告

1、目录一、实验背景2二、实验目的2三、实验主要内容2四、实验具体过程2（一）建立测井曲线道2（二）导入数据3（三）选定目标井段3（四）导入泥质含量、孔隙度、渗透率和饱和度曲线道3（五）编辑调整页面4（六）实验成果图4五、实验结果分析解释4（一）对测井曲线图的大致认识4（二）测井曲线的特征4（三）测井资料综合解释方法原理5（四）实验成果资料分析解释6六、实验总结6七、体会与建议7八、附录8附录1：实验成果图图18附录2：实验成果图图29附录3：泥质含量计算10附录4：孔隙度计算11附录5：渗透率计算12附录6：饱和度计算13一、实验背景作为一名地质学专业的大三学生，本学期我们跟随刘之的老师修读了测

2、井方法原理这门课程，在理论上对测井的各种方法及其原理有了大致了解，但这种了解还是浅层次的，需要实验的佐证，因此在学校和老师的安排下我们进行了此次的上机实验，对实际测井资料进行上机分析，以求达到实践与理论的结合。二、实验目的 1.通过对实际测井资料的分析，使我们深入了解测井数据格式、测井曲线类型，测井曲线的处理流程等，从而为测井资料处理奠定基础。2. 通过上机操作，使我们深入了解到物探专业基本软件carbon的使用方法和基本功能，从而为我们以后物探工作打下坚实的基础。三、实验主要内容1.认识测井曲线基本格式和基本内容。2.使用自然伽马、自然电位、井径及微电阻率测井曲线进行岩性定性识别。3.使用声

3、波时差及微电阻率曲线进行储层物性定性分析。4.使用电阻率测井曲线进行储层含油性定性分析。5.开展分析井段泥质含量、孔隙度、渗透率及含油饱和度定量计算（选择代表性深度点数据计算）6.依据上述识别和划分结果撰写课程上机实验报告。四、实验具体过程（一）建立测井曲线道1.安装卡奔3.0软件2.打开intellexpore，建立9个测井曲线道，即自然电位测井（sp）、自然伽马测井（gr）、井径测井（cal），深双侧向电阻率测井（rlld）、浅双侧向电阻率测井（rlls）、微球形聚焦测井（rmsfl）、声波测井（ac）、密度测井（den）、中子测井（cnl）。（二）导入数据在给定的多口井的实际测井资料中选

4、择一口井以excel方式打开，做一定调整后导入intellexpore中，我选择了yi42井。（三）选定目标井段本次实验做两段40m含砂性较好的井段进行资料分析解释，仔细观察gr、sp曲线对岩性进行定性分析，两者出现负异常时应该显示为砂岩层段，正异常则为泥岩层段，于是我选取了2960-3000m和3090-3130m这两个井段。（四）导入泥质含量、孔隙度、渗透率和饱和度曲线道1.导入泥质含量根据gr曲线计算泥质含量vsh，使用c语言程序实现，具体程序参见附录一。2.导入孔隙度根据ac曲线利用c语言编程计算孔隙度por，具体程序参见附录二。3.导入渗透率利用c语言编程计算渗透率perm，具体程序

5、参见附录三。4.导入饱和度利用c语言编程计算饱和度so，具体程序参见附录四。5.建立sand曲线道复制por曲线道并做一定调整建立sand曲线道。（五）编辑调整页面根据excel中数据的大小并结合测井实际要求确定每一个曲线道的刻度值，并对图头文字，曲线颜色等做一定美化。（六）实验成果图见图1（p8）、图二（p9）。五、实验结果分析解释（一）对测井曲线图的大致认识1.gr、sp、cal反映岩性2.rlld、rlls、rmsfl反映含油性3.ac、den、cnl反映孔渗性（二）测井曲线的特征表1 砂泥岩剖面测井曲线特征地层测井曲线储集层砂岩非储集层泥岩自然电位明显异常泥岩基线自然伽马低高井径缩径扩

6、径电阻率高阻低阻声波<300m/s>300m/s（三）测井资料综合解释方法原理1.岩性分析sp：在淡水泥浆的砂泥岩剖面井中，以大段泥岩层部分的自然电位曲线为基线，此时sp曲线出现负异常的井段都可以认为是渗透性岩层，纯砂岩井段出现最大的负异常，泥质的砂岩负异常幅度变低，而且随着泥质含量的增多而异常幅度下降。gr：在砂泥岩剖面中，砂岩显示出最低值，粘土（泥岩、页岩）显示出最高值，而粉砂岩泥质砂岩介于中间，并随着岩层中泥质含量增加曲线幅度增大。cal：泥岩层和某些松散岩层常常由于钻井时泥浆的浸泡和冲刷造成井壁坍塌，使实际井径大于钻头直径，出现井径扩大；渗透性岩层，常常由于泥浆滤液向岩层中

7、渗透，在井壁上形成泥饼，使实际井径小于钻头直径，出现井径缩小；而在致密岩层处，井径一般变化不大，实际井径接近钻头直径。2.物性分析以及气层识别储层物性反映的是储层质量的好坏，决定了油区的丰度和储量。测井计算反映储层物性的参数主要有孔隙度、渗透率、泥质含量以及粒度中值等。显然储层孔隙度高、渗透率大、泥质含量低则物性好，相反，储层孔隙度低、渗透率小、泥质含量高则物性差。ac：当岩层中含气时，声波时差将显著增大。另外，在泥浆侵入不深的高孔隙度疏松砂岩地层中，油层的声波时差也相应增大，一般比水%。den、cnl组合：天然气相对于地层水和石油而言，其密度很低，密度测井时，其密度测井反映值也较低，故计算出

8、的孔隙度比实际孔隙度偏大，而在中子测井曲线上气层表现为低孔隙度，因此二者曲线重叠即可识别气层。若重叠部分较大，则含气的可能性较大。3.含油性分析rlld、rlls组合：将深、浅侧向视电阻率曲线重叠绘制，观察两条曲线幅度的相对关系，在渗透层井段会出现幅度差。深侧向曲线幅度大于浅侧向曲线幅度，叫做正幅度差（意味着泥浆低侵），这种井段一般可认为是含油气井段；反之当深侧向曲线幅度小于浅侧向曲线幅度时，称为负幅度差（意味着高侵），这种井段可认为是含水井段。rlld、rlls、rmsfl组合：利用深、浅侧向和微球形聚焦三条曲线重叠绘制可以直观判断可动油、气情况，在地层水和泥浆的矿化度基本相同，且侵入不深的

9、情况下，若储集层中有可动油气则三条曲线幅度有下列关系：rmsfl<rlls<rlld；若无可动油气则三条曲线重合。（四）实验成果资料分析解释1.2960m3000m井段分析解释(见p8图1)岩性分析：该段整体来说砂质夹层和泥质砂岩比较多，而纯砂岩很少，仅在2972m和2995m处出现两段不足1m厚的纯砂岩。这一段的平均泥质含量在70%左右。物性分析及气层识别：该段整体来说声波曲线比较稳定，所以此段物性较好。在2965.162968.16m，2979.682981.84m，2984.322985.6m，2992.762994.32m这四段den和cnl有重叠，判断可能为气层。含油性分

10、析：该段整体来说rlldrllsrmsfl，所以认为此段整体基本都含可动油气。2.30903130m井段分析解释（见p9图2）岩性分析：该段整体泥质含量平均在80%左右，30903090.92m为砂质夹层；3091.923095.96m，31133115.8m和3121.843126.28m这三段sp和gr都比较小，泥质含量较低，为泥质砂岩层；3111.53113.08mgr和sp都显示最强负异常，为纯砂岩段；余下段全部为泥岩层，gr与sp显示极为明显的正异常。物性分析及气层识别：声波曲线幅度变化较大，大概出现了三次周波跳跃，分别在3100.53105m，3106.53110.5m和31183

11、122m处。含油性分析：在30903098.88m、3111.483115.92m和3121.963126.68m这三段rlld>rlls>rmsfl，说明它们可能是油层，不存在rlld<rlls的层段。六、实验总结此次通过上机实验，一方面我们加强了对测井软件carbon的熟练运用，另一方面对测井实际工作有了一定了解，还加深了对课本知识的理解。学会了利用测井曲线划分储集层和非储集层，并定性、定量分析解释储层岩性、物性以及含油性。七、体会与建议1.实践与理论相结合才能实现对理论的更高认识。2.测井资料解释与分析需要极好的观察力与洞察力，我还需要好好锻炼。3.在实验的过程中学会善

12、用课本以及网络上的各种资源提高了我做事的独立性。4.在条件允许的情况下，希望学校以后对于测井这样的课程多安排一些上机等实验。八、附录附录1：实验成果图图1图1 yi42井29603000m测井处理成果图附录2：实验成果图图2图2 yi42井30903130m测井处理成果图附录3：泥质含量计算void main()object v;v = newtrace(1, "cvrtrace", "vsh");v.left = 0;v.right = 1;v.linecolor = rgb(255, 0, 0);v.linewidth = 5;v.linestyle

13、 = 0;v.style |= 1 << 19;v.fillindex = 1;float d1, d2, m1, m2;d1 = input("计算开始深度：", 2960);d2 = input("计算结束深度：", 3000);m1 = tracemax("gr",d1,d2);m2 = tracemin("gr",d1,d2); float c = 3.7; float vp("vsh");for(depth =d1; depth <d2; depth +=0.125)

14、 float val = (gr-m2)/(m1-m2);vp = 0.83 * val * val + 1.7 * val;if(vp<0)vp=0; if(vp>1)vp=1;play();备注：将d1和d2分别改为3090、3130即可用于计算30903130m井段泥质含量。附录4：孔隙度计算void main()object v;v = newtrace(1, "cvrtrace", "por");play();v.left = 0;v.right = 1;v.linecolor = rgb(0, 127, 127);v.linewi

15、dth = 5;v.linestyle = 0;v.style |= 1 << 19;v.fillindex = 2;float ma, mf; ma = input("岩石骨架的声波时差:", 182);mf = input("泥浆滤液的声波时差:", 620);float vp("por");for(depth =2960; depth <3000; depth +=0.125) vp = (ac - ma)/(mf - ma)-vsh*(287 - ma)/(mf - ma);if(vp<0.01)vp=

16、0.01; if(vp>0.35)vp=0.35;play();备注：将语句for(depth =2960; depth <3000; depth +=0.125)，更改为for(depth =3090; depth <3130; depth +=0.125) ，将语句vp = (ac - ma)/(mf - ma)-vsh*(287 - ma)/(mf - ma)，更改为vp = (ac - ma)/(mf - ma)-vsh*(314 - ma)/(mf - ma)，即可用于计算30903130m井段的孔隙度。附录5：渗透率计算/ y4井岩心分析渗透率与孔隙度相关性比较好

17、/ 关系式：perm = 0.0014*e(0.3482*por)，(e = 2.71828182845904)/用此关系式计算渗透率void main()object v;v = newtrace(1, "cvrtrace", "perm");v.left = 0.01;v.right = 1000;v.linecolor = rgb(0, 127, 0);v.linewidth = 5;v.linestyle = 0;v.style |= 2048;float vp("perm");for(depth =2960; depth &

18、lt;3000; depth +=0.125) vp = 0.0014*pow(2.718, 0.3482*por*100);if(vp<0.001)vp=0.001; if(vp>1000)vp=1000;play();备注：将depth分别改为3090、3130即可用于计算30903130m井段的渗透率。附录6：饱和度计算/y4井为例,计算饱和度应首先有名为“por_den”的孔隙度曲线道（先执行“孔隙度计算b”方法即可）。/阿尔齐公式计算含油饱和度/公式：so = 1- (a*b*rw)/(rt * porm)n/ rw:地层水电阻率， rt:电阻率, por:孔隙度；/ a = 0.8, b = 1.1467； m = 1.8701, n = 2.1025void main()object v;v = newtrace(0, "cvrtrace", "so_auto");v.left = 0;v.right = 100;v.linecolor = rgb(255, 0, 0);v.linewidth = 5;v.linestyle = 0;float a,b,m,n;a = 0.8, b