测井课程设计

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

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

3、和基本内容。2. 使用自然伽马、自然电位、井径及微电阻率测井曲线进行岩性定性识别。3. 使用声波时差及微电阻率曲线进行储层物性定性分析。4. 使用电阻率测井曲线进行储层含油性定性分析。5. 开展分析井段泥质含量、 孔隙度、渗透率及含油饱和度定量计算 （选择代表性深度 点数据计算）6. 依据上述识别和划分结果撰写课程上机实验报告。四、实验具体过程（一）建立测井曲线道1. 安装卡奔 3.0 软件2. 打开IntellExpore ，建立9个测井曲线道，即自然电位测井（SP）、自然伽马测井（GR、井径测井（CAL,深双侧向电阻率测井（RLLD、浅双侧向电阻率测井（RLLS、 微球形聚焦测井（RMSF

4、、声波测井（AC、密度测井（DE、中子测井（CNL。（二、导入数据在给定的多口井的实际测井资料中选择一口井以 EXCEL方式打开，做一定调整后导入IntellExpore 中,我选择了 YI42 井。（三、选定目标井段本次实验做两段40m含砂性较好的井段进行资料分析解释，仔细观察 GR SP曲线对 岩性进行定性分析， 两者出现负异常时应该显示为砂岩层段， 正异常则为泥岩层段， 于是 我选取了 2960-3000m和3090-3130m这两个井段。（四、导入泥质含量、孔隙度、渗透率和饱和度曲线道1. 导入泥质含量根据GR曲线计算泥质含量Vsh,使用C语言程序实现，具体程序参见附录一。2. 导入孔

5、隙度根据AC曲线利用C语言编程计算孔隙度POR具体程序参见附录二。3. 导入渗透率利用C语言编程计算渗透率PERM具体程序参见附录三。4. 导入饱和度利用C语言编程计算饱和度So,具体程序参见附录四。5. 建立Sand曲线道复制POR曲线道并做一定调整建立Sand曲线道（五）编辑调整页面根据EXCEL中数据的大小并结合测井实际要求确定每一个曲线道的刻度值，并对图头 文字，曲线颜色等做一定美化。（六）实验成果图见图 1（ P8）、图二（P9）。五、实验结果分析解释（一）对测井曲线图的大致认识1. GR、SP、CALS映岩性2. RLLD RLLS RMSF反映含油性3. AC、DEN CNL反映

6、孔渗性（二）测井曲线的特征表1砂泥岩剖面测井曲线特征测井曲线层储集层一砂岩非储集层一泥岩自然电位明显异常泥岩基线自然伽马低高井径缩径扩径电阻率高阻低阻声波300 卩 m/s（三）测井资料综合解释方法原理1. 岩性分析SP:在淡水泥浆的砂泥岩剖面井中，以大段泥岩层部分的自然电位曲线为基线，此时SP 曲线出现负异常的井段都可以认为是渗透性岩层，纯砂岩井段出现最大的负异常，泥 质的砂岩负异常幅度变低，而且随着泥质含量的增多而异常幅度下降。GR在砂泥岩剖面中，砂岩显示出最低值，粘土（泥岩、页岩）显示出最高值，而粉 砂岩泥质砂岩介于中间，并随着岩层中泥质含量增加曲线幅度增大。CAL泥岩层和某些松散岩层常

7、常由于钻井时泥浆的浸泡和冲刷造成井壁坍塌，使实 际井径大于钻头直径，出现井径扩大；渗透性岩层， 常常由于泥浆滤液向岩层中渗透， 在 井壁上形成泥饼，使实际井径小于钻头直径，出现井径缩小； 而在致密岩层处，井径一般 变化不大，实际井径接近钻头直径。2. 物性分析以及气层识别储层物性反映的是储层质量的好坏， 决定了油区的丰度和储量。 测井计算反映储层物 性的参数主要有孔隙度、渗透率、泥质含量以及粒度中值等。显然储层孔隙度高、渗透率 大、泥质含量低则物性好，相反，储层孔隙度低、渗透率小、泥质含量高则物性差。AC 当岩层中含气时， 声波时差将显著增大。 另外，在泥浆侵入不深的高孔隙度疏松 砂岩地层中，

8、油层的声波时差也相应增大，一般比水 %。DEN CNL组合：天然气相对于地层水和石油而言，其密度很低，密度测井时，其密 度测井反映值也较低， 故计算出的孔隙度比实际孔隙度偏大， 而在中子测井曲线上气层表 现为低孔隙度， 因此二者曲线重叠即可识别气层。 若重叠部分较大， 则含气的可能性较大。3. 含油性分析RLLD、 RLLS 组合：将深、浅侧向视电阻率曲线重叠绘制，观察两条曲线幅度的相对 关系，在渗透层井段会出现幅度差。 深侧向曲线幅度大于浅侧向曲线幅度， 叫做正幅度差 （意味着泥浆低侵） ，这种井段一般可认为是含油气井段；反之当深侧向曲线幅度小于浅 侧向曲线幅度时，称为负幅度差（意味着高侵）

9、 ，这种井段可认为是含水井段。RLLD RLLS RMSF组合：禾用深、浅侧向和微球形聚焦三条曲线重叠绘制可以直观判断可动油、气情况， 在地层水和泥浆的矿化度基本相同，且侵入不深的情况下，若储集层中有可动油气则三条曲线幅度有下列关系：RMSFLvRLLSvRL；若无可动油气则三条曲线重合。（四）实验成果资料分析解释1.2960m 3000m井段分析解释（见P8图1） 岩性分析：该段整体来说砂质夹层和泥质砂岩比较多，而纯砂岩很少，仅在2972m和2995m处出现两段不足1m厚的纯砂岩。这一段的平均泥质含量在70流右。 物性分析及气层识别：该段整体来说声波曲线比较稳定，所以此段物性较好。在2965

10、.16 2968.16m, 2979.68 2981.84m, 2984.32 2985.6m, 2992.76 2994.32m 这四 段DEN和CNL有重叠，判断可能为气层。 含油性分析：该段整体来说RLL RLLS RMSFL所以认为此段整体基本都含可动 油气。2.3090 3130m井段分析解释（见 P9图2）岩性分析：该段整体泥质含量平均在 80%左右，3090 3090.92m为砂质夹层； 3091.92 3095.96m, 3113 3115.8m 和 3121.84 3126.28m这三段 SP和 GR都比较小， 泥质含量较低，为泥质砂岩层；3111.5 3113.08mGF和

11、SP都显示最强负异常，为纯砂岩 段；余下段全部为泥岩层，GR与 SP显示极为明显的正异常。物性分析及气层识别：声波曲线幅度变化较大， 大概出现了三次周波跳跃， 分别在3100.5 3105m 3106.5 3110.5m 和 3118 3122m处。含油性分析：在 3090 3098.88m 3111.48 3115.92m和 3121.96 3126.68m这三 段RLLDRLLSRMSF说明它们可能是油层，不存在 RLLDvRLLS勺层段。六、实验总结此次通过上机实验，一方面我们加强了对测井软件Carbon的熟练运用，另一方面对测井实际工作有了一定了解， 还加深了对课本知识的理解。 学会了

12、利用测井曲线划分储集 层和非储集层，并定性、定量分析解释储层岩性、物性以及含油性。七、体会与建议1. 实践与理论相结合才能实现对理论的更高认识。2. 测井资料解释与分析需要极好的观察力与洞察力，我还需要好好锻炼。3. 在实验的过程中学会善用课本以及网络上的各种资源提高了我做事的独立性。4. 在条件允许的情况下，希望学校以后对于测井这样的课程多安排一些上机等实验。八、附录附录 1：实验成果图图 10GR 15000騎厂 - - /*/ /F：L; -ui-T境亠 毎0 諮胡 2988 蜃_RLLD1 100RLLSwoRMSFL1 - 100AC450 L Nc -5 I 1I I P j a

13、亠.1501dcills00Ta85图1 YI42 井2960 3000m测井处理成果图附录2:实验成果图图2图2 YI42 井3090-3130m测井处理成果图附录 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 = 0;v.style |= 1 19;v.FillIndex = 1;float d1, d2, m1, m2;d1 = Input( 计算开始深度：

14、, 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) float val = (GR-m2)/(m1-m2); vp = 0.83 * val * val + 1.7 * val;if(vp1)vp=1;Play();备注：将di和d2分别改为3090、3130即可用于计算3090-3130m井段泥质含量附录 4：孔隙度计算void main(

15、)object v;v = NewTrace(1, CvrTrace, Por);Play();v.Left = 0;v.Right = 1;v.LineColor = RGB(0, 127, 127);v.LineWidth = 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

16、= (AC - ma)/(mf - ma)-Vsh*(287 - ma)/(mf - ma);if(vp0.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) 即可用于计算3090-3130m井段的孔隙

17、度。附录 5：渗透率计算/ y4 井岩心分析渗透率与孔隙度相关性比较好/关系式：PERM = 0.0014*eA(0.3482*Por) , (e = 2.704)/ 用此关系式计算渗透率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 3000; dept

18、h +=0.125) vp = 0.0014*Pow(2.718, 0.3482*POR*100);if(vp1000)vp=1000;Play();备注：将depth分别改为3090、3130即可用于计算3090 3130m井段的渗透率附录 6：饱和度计算y4井为例,计算饱和度应首先有名为“ POR_DE”的孔隙度曲线道(先执行“孔隙 度计算B”方法即可)。/ 阿尔齐公式计算含油饱和度 / 公式：So = 1- (a*b*Rw)/(Rt * PorAm)An/ 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,