岩石力学大作业.docx

1、岩石力学大作业岩石力学结课大作业张宗仁1 前言：在石油勘探开发过程中，由于地层岩石其特有的性质，在钻井过程中，由于井壁围岩失稳所造成漏、喷、塌、卡等安全事故，以及开采过程中，出沙、套管损坏等问题已经造成巨大损失。利用岩石力学知识，利用测井数据分析井眼周围岩石力学性质，是解决这些问题的有效方法。本文利用给出的测井数据，对地层力学参数、地应力、地层坍塌压力与破裂压力进行分析计算，作出地层力学参数、地层主应力、地层坍塌、破裂压力剖面，分析井壁坍塌原因；研究储层段的出砂可能性。2 利用自然伽马测井数据简单分析地层岩性21 估算地层泥质含量自然伽玛测井是在井内测量岩层中自然存在的放射性核素核衰变过程中放

2、射出来的射线的强度，它可用于划分岩性，估算地层泥质含量。 由于泥质颗粒细小，具有较大的比面，使它对放射性物质有较大的吸附能力，并且沉积时间长，有充分的时间与溶液中的放射性物质一起沉积下来，所以泥质有很高的放射性。在不含放射性矿物的情况下，泥质含量的多少就决定了沉积岩石的放射性的强弱。所以有可能利用自然伽玛测井资料来估算泥质的含量，具体方法有两种1 P125：相对值法 （2-1）式中：为泥质的体积含量。GCUR为希尔奇指数，与地质时代有关，可根据取芯分析资料与自然伽玛测井值进行统计确定，对于第三系地层取，老地层取2（本作业中取2）。为泥质含量指数。 （2-2）分别表示目的层的、纯泥岩层的和纯砂岩

3、层的自然伽玛值。本测井数据中=；=。计算出的泥质含量与井深的关系见第8部分图1。22 估算地层biot系数当泥质含量时，biot系数=，当泥质含量时，biot系数=，当时，biot系数。计算出的biot系数与井深的关系见第8部分图2。3 利用测井数据计算分析地层的弹性模量、泊松比。31纵、横声波速度2声波速度测井是测量地层声波速度的测井方法，声波速度测井可测量滑行波通过地层传播的时差T，纵波时差ts和横波时差tp可从由测井公司提供的测井曲线或磁盘数据中得到，经过换算即可得到纵横波速度为： （m/s） （3-1）式中：为纵波时差，；为横波时差，。大部分的油田测井作业中，并不作全波列测井，即缺少横

4、波测井资料，因此针对某一地层，就需要借助经验公式来估算横波速度。对于大多数地层，其柏松比一般在之间，因此有： （3-2）计算出纵、横声波速度与井深的关系见第8部分图3。利用B井测井数据计算分析地层强度参数32岩石的动态弹性参数确定2岩石为各向同性无限弹性体，则根据纵波速度和横波速度计算动态泊松比和动态杨氏模量的关系为： (3-3) (3-4)其中：为纵波速度；为横波速度；为岩石密度；为动态泊松比；为动态弹性模量。应用上式以及从密度测井和声波测井得到的一系列数据，我们能确定整个井深剖面内的岩石弹性参数。33岩石的动态弹性参数与静态弹性参数的转化3根据文献3岩石的动态参数与静态参数的关系式可以表示

5、为： （3-5） （3-6）动态弹性模量与静态弹性模量井深的关系见第8部分图4；动态泊松比与静态泊松比与井深的关系见第8部分图5。4 利用测井数据计算粘聚力、内摩擦角与地层抗拉强度的连续剖面长期以来，一直有不少研究人员尝试着从测井资料获取岩石特性参数。地层声波测井反映声波在岩石中的传播速度，它与岩石的密度、孔隙度、结构强度等密切相关，它作为衡量岩石强度参数的一个重要指标，长期以来为众多学者所重视。41岩石单轴抗压强度的确定2Miller和Deere在实验基础之上建立了岩石单轴抗压和岩石弹性模量、粘土含量之间的关系； （4-1）式中：为砂岩的泥质含量；为砂岩的动态杨氏模量，MPa。单轴抗压强度与

6、井深的关系见第8部分图6Cotaes等人继Miller和Deere之后，提出了岩石固有剪切强度（）与单轴抗压强度（）之间的关系 （4-2）岩石实验表明，岩石的单轴抗压强度一般是其抗拉强度的815倍，因此，可以用下式近似计算岩石的抗拉强度： （4-3）式中：为岩石的抗拉强度，MPa。Kcl=815,本次计算取15。抗拉强度与井深的关系见第8部分图742岩石粘聚力和内摩擦角的确定Coast提出沉积岩的年粘聚力和单轴抗压强度（）的经验关系式： （4-4）式中：为岩石的动态体积压缩模量。 粘聚力与井深的关系见第8部分图7对岩石强度内摩擦角的计算，在斯伦贝谢公司推出的力学稳定性测井软件中假定所有岩石的内

7、摩擦角为30，这与实际情况不符，岩石的类型、颗粒大小等均对有很大的影响。一般岩石的值与值存在着一定的对应关系。根据玛湖油田13组岩芯的实测强度参数值，通过回归分析得到泥页岩地层内摩擦角与粘聚力间的相关关系式为： （4-5）内摩擦角与井深的关系见第8部分图85 计算地层地应力地层间或者层间的不同岩性岩石的物理特性、力学特性以及地层孔隙压力等方面的差别造成了层间或者层内应力分布的非均匀性。测井资料可以有效地计算这些特征量。根据上述测井资料解释出来的岩层物理、力学参数可以有效地解释地应力。由于0-700米的井段缺少测井数据，我们需要从下部的底层的测井数据推倒0-700岩石参数，由于在自然沉积过程中，

8、由于岩石的压实效应，岩石的密度往往随地层深度的增加而增加并且近似的成正比关系，我们不妨用700-3800米的地层密度测井数据拟和岩石密度与井深的变化关系。 (5-1)式中：H为井深，A、B为拟合数据。本次计算中拟合结果为：当H=700米处，该处的上腹岩层压力约为 (5-2)黄荣樽等假设地下岩层的地应力主要由上覆岩层压力与水平方向构造应力产生，且水平方向的构造应力与上覆岩层压力成正比4P92。 (5-3)式中：、为地应力系数；、为地应力系数；是区域地层的特性值，那么我们可根据水力压裂试验的方法反求它们的值。 (5-4)式中： 为有效应力系数； 为孔隙压力（MPa）。当井深2450m时，=59Mp

9、a，=49Mpa，由测井数据及以上内容计算：Vcl= 带入上式可算得：把上述系数带入式（5-3）可以求得地层地应力随井深的变化规律。计算结果见第8部分图9。从计算结果可以看出，在该底层段，最小水平主应力 上覆岩层压力2104MPa时，在正常生产时油层不会出砂；（2）当2104MPaB104MPa时，油层微量出砂，但油层见水后就会严重出砂，需生产中适时防砂；（3）当B2104MPa，在正常生产时油层不会出砂；（2）当2104MPaEc104MPa时，油层微量出砂；（3）当Ec104MPa时，正常生产时会严重出砂。组合模量与井深的关系见第8部分图13，从结果可以看出，该储层岩石强度较高，在储层段E

10、c都大于20000MPa，在油井正常生产时不出砂。8 计算结果及其分析81利用自然伽马测井数据简单分析地层泥质含量 图1地层泥质含量与井深关系曲线 图2 Boit系数与井深关系曲线82地层的动静态弹性模量、泊松比图3 纵横波波速 图4 地层地层动静态弹性模量 图5 地层地层动静态泊松比83地层岩石单轴抗拉强度、粘聚力、内摩擦角与地层抗拉强度的连续剖面 图6 地层岩石单轴抗压强度 图7 地层岩石内磨擦角图8 地层岩石粘聚力与抗拉强度84地层地应力图9 地层地应力当量密度85地层坍塌压力和破裂压力 图10地层坍塌压力和破裂压力 图11 井径从上图可以看出在该井在2600米3300米处出现了较为严重

11、的扩径现象，在坍塌破裂压力曲线图上我们不难发现在该井段泥浆密度低于坍塌压力，在下部井段提高泥浆密度过后，扩径明显减小，井眼质量得到改善。从计算结果可以看出不难发现发现，该井在部分井段上有井壁坍塌得可能性，从井径与井深关系曲线（第8部分图10），我们不难看出在该井700-2500米存在不同程度的坍塌或者缩径，在的（2600-2700米）井段以及（2700-3300米）井段出现了严重的扩径现象，（3300-3800米）井段扩径现象减小，一方面由于泥浆密度提高，一方面由于岩石强度增大以及地应力水平相对降低。86用B指数法预测地层出砂可能性图12 B指数图13 组合模量从结果可以看出，该储层岩石强度较

12、高，在储层段B指数基本上都大于20000MPa， Ec都大于20000MPa，在油井正常生产时不出砂。9 结论与总结本文根据声波、密度、伽玛测井资料及泥浆密度等参数，利用本课程中所学到的岩石强度参数计算公式计算出各地层的强度参数，再根据密度对井深的积分计算出上覆岩层应力，利用压裂试验数据，通过六五模式，强度参数计算了最大、最小主应力。然后用以上所得参数计算出地层的坍塌及破裂压力。分析了在井段出现扩径的原因。最后利用B指数和组合模量法对出砂情况进行预测。通过计算，该储层正常生产时，出砂的可能性较小。通过学习本课程，我对石油工程岩石力学有了更深更系统的了解，为以后继续学习和工作奠定了一定的基础，蔚宝华老师讲课认真细致，对我理解本课程的知识点帮助很大，感谢蔚老师！参考文献1 丁次乾.矿场地球物理M.中国石油大学出版社，2 楼一珊,金业权.岩石力学与石油工程 M.石油工业出版社，3 林英松 葛洪魁.岩石动静力学参数的试验研究J.岩石力学与工程学报，1998,17(2)2162224 陈勉,金衍,张广清.石油工程岩石力学M.北京：科学出版社，,287-295.