CN108875115B - 一种确定岩石强度的方法 - Google Patents

一种确定岩石强度的方法 Download PDF

Info

Publication number
CN108875115B
CN108875115B CN201810201388.XA CN201810201388A CN108875115B CN 108875115 B CN108875115 B CN 108875115B CN 201810201388 A CN201810201388 A CN 201810201388A CN 108875115 B CN108875115 B CN 108875115B
Authority
CN
China
Prior art keywords
rock
strength
content
weight coefficient
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810201388.XA
Other languages
English (en)
Other versions
CN108875115A (zh
Inventor
崔式涛
孙佩
刘春雷
李娜
王谦
段朝伟
牟明洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN201810201388.XA priority Critical patent/CN108875115B/zh
Publication of CN108875115A publication Critical patent/CN108875115A/zh
Application granted granted Critical
Publication of CN108875115B publication Critical patent/CN108875115B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种确定岩石强度的方法,包括:首先利用岩心刻度常规测井资料或ECS资料,建立参数模型,计算储层矿物含量,分析岩石强度参数与矿物含量之间的对应关系,通过引入多矿物含量构建新的函数方程,求取单轴抗压强度,固有剪切强度,最后利用得到的岩石强度参数计算方程,代入岩石矿物含量以及岩石力学参数,即可得到岩石强度参数,并结合测井资料连续评价地层特征。相比于原有的计算方法,能够更准确的计算储层强度参数,满足储层可压裂性分析的需求。经本发明所求得的岩石强度参数,与岩石力学实验测的岩石强度参数符合性好,准确度高,可以及时服务或满足测井资料解释评价的需要,为目前我国储层压裂评价提供可靠的依据。

Description

一种确定岩石强度的方法
技术领域
本发明涉及储层压裂评价领域,特别涉及一种确定岩石强度的方法。
背景技术
随着我国石油勘探领域中对复杂储层评价的不断加深,地层的可压裂性评价成为了部分地区必不可少的评价内容。计算岩石的强度是地层可压裂性评价中的基础,国内外对于岩石强度计算方法已经有很多研究,现在最常用的是Coates和Denco(1981)在Deere和Miller(1966)实验基础上建立了计算岩石抗强的关系式:
单轴抗压强度:
SC=Ed[0.008Vsh+0.0045(1-Vsh)]
固有剪切强度:
Figure BDA0001594590650000011
单轴抗拉强度:
ST=SC/12
但在实际应用中,不同地区不同地层效果不尽相同,不能满足准确计算岩石强度参数的需求。本发明通过矿物含量对岩石强度参数进行拟合,建立了一种更为精确的岩石强度计算方法,满足复杂储层评价的需求。
发明内容
为解决现有技术中存在的上述缺陷,本发明的目的在于提供一种确定岩石强度的方法,该方法通过矿物含量与岩石强度参数之间建立关系,能够较为准确的计算储层的岩石强度参数,解决了由于原有计算模型由于缺少地区参数而无法较为准确的计算岩石气囊度参数的问题,能够及时服务或满足实际生产的需要,实现了油气储层可压裂性评价的需求。
本发明是通过下述技术方案来实现的。
所述方法按照以下步骤进行操作:
步骤(1),获取岩心刻度常规测井资料或ECS资料,建立参数模型,计算储层矿物含量,分析岩石强度参数与矿物含量之间的对应关系。
步骤(2),利用计算方程求取岩石单轴抗压强度和固有剪切强度参数;
步骤(3),利用步骤(2)得到的岩石强度参数计算方程,代入岩石矿物含量以及岩石力学参数,即可得到岩石强度参数;
步骤(4),根据求得的岩石强度参数结合测井资料连续评价地层特征。
上述技术方案中,还包括进一步限定的方案:
所述步骤(1),计算储层矿物含量按照下式计算:
Figure BDA0001594590650000021
Figure BDA0001594590650000022
Figure BDA0001594590650000023
砂质含量:
VSAND=ea×P+b×N-c×M-d
泥质含量:
VSH=j×GR-f×M+g×P+h×N-i
式中:Δtf,Δtma,Δt分别为地层水、岩石骨架、岩石的声波时差,μs/ft;ρfmab分别为地层水的密度、岩石骨架密度、体积密度,g/cm3
Figure BDA0001594590650000032
Figure BDA0001594590650000033
分别为地层水、岩石骨架、岩石的中子孔隙度,%;GR为自然伽马曲线值,API,a为砂质含量的声波中子比值权重系数,b为砂质含量的中子密度比值权重系数,c为砂质含量的声波密度比值权重系数,d为砂质含量的权重系数;j为泥质含量的自然伽马值权重系数,f为泥质含量的声波密度比值权重系数,g为泥质含量的声波中子比值权重系数,h为泥质含量的中子密度比值权重系数,i为泥质含量的权重系数。
所述步骤(2),通过下述公式计算岩石强度参数:
单轴抗压强度:
Sc=E*(x1*TOC+x2*VSAND+x3*VSH+x4)
固有剪切强度:
C0=X*(y1*TOC+y2*VSH+y3*VSAND+y4*VLIME+y5)
其中:
Figure BDA0001594590650000031
式中:Sc为单轴抗压强度,MPa;E为动态杨氏模量,GPa;TOC为有机碳含量,%;VSAND为砂质含量,%;VSH为泥质含量,%;C0为固有剪切强度,MPa;VLIME为砂质含量,%;x1,x2,x3,x4,y1,y2,y3,y4,y5为待定拟合参数;VP为纵波速度;ρb为流体密度;υ为泊松比。
步骤(3)利用步骤(2)得到的岩石强度参数计算方程,代入岩石矿物含量以及岩石力学参数,即可得到岩石强度参数。
所述步骤(4),求得的岩石强度参数结合测井资料连续评价地层特征,通过下述方法进行:
通过ECS资料结合常规测井资料确定地层矿物含量,根据不同矿物含量确定岩石强度,得到全井段连续分布的抗压强度、抗张强度及固有剪切强度,进而评价连续地层的可压裂性。
相比于原有的计算方法,本发明方法能够更准确的计算储层强度参数,满足储层可压裂性分析的需求。经本发明所求得的岩石强度参数,与岩石力学实验测的岩石强度参数符合性好,准确度高,岩石强度参数中,抗张强度相对误差不大于10.81%,抗压强度相对误差不大于12.91%,固有剪切强度相对误差不大于16.25%。可以及时服务或满足测井资料解释评价的需要,为目前我国储层压裂评价提供可靠的依据。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明计算的岩石强度参数计算效果图。
具体实施方式
下面结合附图和实施例对发明作进一步的详细说明,但并不作为对发明做任何限制的依据。
本发明确定岩石强度的方法,包括以下步骤:
步骤(1),获取岩心刻度常规测井资料或ECS资料,计算储层矿物含量,分析岩石强度参数与矿物含量之间的对应关系;
计算储层矿物含量按照下式计算:
Figure BDA0001594590650000051
Figure BDA0001594590650000052
Figure BDA0001594590650000053
砂质含量:
VSAND=ea×P+b×N-c×M-d
泥质含量:
VSH=j×GR-f×M+g×P+h×N-i
式中:Δtf,Δtma,Δt分别为地层水、岩石骨架、岩石的声波时差,μs/ft;ρfmab分别为地层水的密度、岩石骨架密度、体积密度,g/cm3
Figure BDA0001594590650000054
Figure BDA0001594590650000055
分别为地层水、岩石骨架、岩石的中子孔隙度,%;GR为自然伽马曲线值,API,a为砂质含量的声波中子比值权重系数,b为砂质含量的中子密度比值权重系数,c为砂质含量的声波密度比值权重系数,d为砂质含量的权重系数;j为泥质含量的自然伽马值权重系数,f为泥质含量的声波密度比值权重系数,g为泥质含量的声波中子比值权重系数,h为泥质含量的中子密度比值权重系数,i为泥质含量的权重系数。
步骤(2),利用计算方程求取岩石单轴抗压强度和固有剪切强度参数;
通过下述公式计算岩石强度参数:
单轴抗压强度:
Sc=E*(x1*TOC+x2*VSAND+x3*VSH+x4)
固有剪切强度:
C0=X*(y1*TOC+y2*VSH+y3*VSAND+y4*VLIME+y5)
其中:
Figure BDA0001594590650000061
式中:Sc为单轴抗压强度,MPa;E为动态杨氏模量,GPa;TOC为有机碳含量,%;VSAND为砂质含量,%;VSH为泥质含量,%;C0为固有剪切强度,MPa;VLIME为砂质含量,%;x1,x2,x3,x4,y1,y2,y3,y4,y5为待定拟合参数;VP为纵波速度;ρb为流体密度;υ为泊松比。
结合岩石力学实验资料与岩石矿物组分含量资料,建立该区块的岩石强度参数计算模型。
步骤(3),利用步骤(2)得到的岩石强度参数计算方程,代入岩石矿物含量以及岩石力学参数,即可得到岩石强度参数;
求得的岩石强度参数结合测井资料连续评价地层特征,通过下述方法进行:
通过ECS资料结合常规测井资料确定地层矿物含量,根据不同矿物含量确定岩石强度,得到全井段连续分布的抗压强度、抗张强度及固有剪切强度,进而评价连续地层的可压裂性。
图1是本发明计算的岩石强度参数计算效果图。从图中可以看出,求取的抗张强度值、抗压强度值及固有剪切强度值,与岩石力学试验数据对比可见,所计算的岩石强度参数相对误差率低,准确度高,应用效果好。
Figure BDA0001594590650000062
Figure BDA0001594590650000071
从上述表1可以看出,岩石强度参数中,抗张强度相对误差不大于10.81%,抗压强度相对误差不大于12.91%,固有剪切强度相对误差不大于16.25%。
本发明通过矿物含量与岩石强度参数之间建立关系,能够较为准确的计算储层的岩石强度参数,解决了由于原有计算模型由于缺少地区参数而无法较为准确的计算岩石气囊度参数的问题,相比于原有的计算方法,能够更准确的计算储层强度参数,与岩石力学实验测的岩石强度参数符合性更好,可以及时服务或满足测井资料解释评价的需要,应用效果好。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内一种。

Claims (4)

1.一种确定岩石强度的方法,其特征在于,包括以下步骤:
步骤(1),获取岩心刻度常规测井资料或ECS资料,计算储层矿物含量,分析岩石强度参数与矿质含量之间的对应关系;
步骤(2),利用计算方程求取岩石单轴抗压强度和固有剪切强度参数;
通过下述公式计算岩石强度参数:
单轴抗压强度:
Sc=E*(x1*TOC+x2*VSAND+x3*VSH+x4)
固有剪切强度:
C0=X*(y1*TOC+y2*VSH+y3*VSAND+y4*VLIME+y5)
其中:
Figure FDA0003482801260000011
式中:Sc为单轴抗压强度,MPa;E为动态杨氏模量,GPa;TOC为有机碳含量,%;VSAND为砂质含量,%;VSH为泥质含量,%;C0为固有剪切强度,MPa;VLIME为砂质含量,%;x1,x2,x3,x4,y1,y2,y3,y4,y5分别为待定拟合参数;VP为纵波速度;ρb为流体密度;υ为泊松比;
步骤(3),利用步骤(2)得到的岩石单轴抗压强度和固有剪切强度参数计算方程,代入岩石矿物含量以及岩石力学参数,即可得到岩石强度参数;
步骤(4)根据求得的岩石强度参数结合测井资料连续评价地层特征。
2.如权利要求1所述的确定岩石强度的方法,其特征在于,所述步骤(1),计算储层矿物含量按照下式计算:
Figure FDA0003482801260000012
Figure FDA0003482801260000013
Figure FDA0003482801260000021
砂质含量:
VSAND=ea×P+b×N-c×M-d
泥质含量:
VSH=j×GR-f×M+g×P+h×N-i
式中:Δtf,Δtma,Δt分别为地层水、岩石骨架、岩石的声波时差,μs/ft;ρfmab分别为地层水的密度、岩石骨架密度、体积密度,g/cm3
Figure FDA0003482801260000022
Figure FDA0003482801260000023
分别为地层水、岩石骨架、岩石的中子孔隙度,%;GR为自然伽马曲线值,API,a为砂质含量的声波中子比值权重系数,b为砂质含量的中子密度比值权重系数,c为砂质含量的声波密度比值权重系数,d为砂质含量的权重系数;j为泥质含量的自然伽马值权重系数,f为泥质含量的声波密度比值权重系数,g为泥质含量的声波中子比值权重系数,h为泥质含量的中子密度比值权重系数,i为泥质含量的权重系数。
3.如权利要求1所述的确定岩石强度的方法,其特征在于,所述步骤(4),求得的岩石强度参数结合测井资料连续评价地层特征,通过下述方法进行:
通过ECS资料结合常规测井资料确定地层矿物含量,根据不同矿物含量确定岩石强度,得到全井段连续分布的抗压强度、抗张强度及固有剪切强度,进而评价连续地层的可压裂性。
4.如权利要求1所述的确定岩石强度的方法,其特征在于,岩石强度参数中,抗张强度相对误差不大于10.81%,抗压强度相对误差不大于12.91%,固有剪切强度相对误差不大于16.25%。
CN201810201388.XA 2018-03-12 2018-03-12 一种确定岩石强度的方法 Active CN108875115B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810201388.XA CN108875115B (zh) 2018-03-12 2018-03-12 一种确定岩石强度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810201388.XA CN108875115B (zh) 2018-03-12 2018-03-12 一种确定岩石强度的方法

Publications (2)

Publication Number Publication Date
CN108875115A CN108875115A (zh) 2018-11-23
CN108875115B true CN108875115B (zh) 2022-05-10

Family

ID=64326054

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810201388.XA Active CN108875115B (zh) 2018-03-12 2018-03-12 一种确定岩石强度的方法

Country Status (1)

Country Link
CN (1) CN108875115B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257134A (zh) * 2018-11-30 2020-06-09 中国石油天然气股份有限公司 剪切应力的连续深度处理方法
CN111596029A (zh) * 2019-09-01 2020-08-28 中国铁道科学研究院集团有限公司铁道建筑研究所 一种适用于机制骨料卵石母岩性能快速评价方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308433A (zh) * 2013-05-03 2013-09-18 中国石油天然气集团公司 一种基于孔隙演化分析评价致密砂岩储层成岩相的方法
CN104345339A (zh) * 2013-07-25 2015-02-11 中国石油天然气集团公司 一种用阵列声波测井资料确定岩石脆性系数的方法
CN106097125A (zh) * 2016-06-17 2016-11-09 延长油田股份有限公司 一种致密砂岩储层可压裂性评价方法
CN106547034A (zh) * 2016-11-09 2017-03-29 西南石油大学 一种计算致密储层岩石脆性指数的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308433A (zh) * 2013-05-03 2013-09-18 中国石油天然气集团公司 一种基于孔隙演化分析评价致密砂岩储层成岩相的方法
CN104345339A (zh) * 2013-07-25 2015-02-11 中国石油天然气集团公司 一种用阵列声波测井资料确定岩石脆性系数的方法
CN106097125A (zh) * 2016-06-17 2016-11-09 延长油田股份有限公司 一种致密砂岩储层可压裂性评价方法
CN106547034A (zh) * 2016-11-09 2017-03-29 西南石油大学 一种计算致密储层岩石脆性指数的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Effects of Grain Scale Heterogeneity on Rock Strength and the Chipping Process;Marlene C.Villeneuve 等;《International Journal of Geomechanics》;20121231;第1-60页 *
利用测井资料评价泥页岩油气"五性"指标;张晋言 等;《测井技术》;20120420;第36卷(第02期);第146-153页 *
基于多矿物模型分析的孔隙度计算方法在M地区中的应用;李权 等;《贵州师范大学学报(自然科学版)》;20161015;第34卷(第05期);第77-82页 *
柳沟庄泥云岩储层岩性与矿物组分特性研究;王俊骏 等;《石化技术》;20160228(第2期);第178-179页 *

Also Published As

Publication number Publication date
CN108875115A (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
CN109441422B (zh) 一种页岩气井间距优化开采方法
CN106761677B (zh) 页岩气水平井单井产能的测井预测方法
CN111460601B (zh) 基于岩石物理建模的正交各向异性地层地应力预测方法
CN104314563B (zh) 一种煤层气储层可压裂性的测井定量评价方法
CN111460602B (zh) 基于岩石物理建模的横观各向同性地层地应力预测方法
CN113901681B (zh) 一种全寿命周期页岩气储层双甜点三维可压性评估方法
CN107038313B (zh) 基于数值岩芯的分层地应力精细描述方法
CN105822298B (zh) 基于产气指数的页岩气层绝对无阻流量的获取方法
CN104007485A (zh) 一种建立复杂孔隙泥页岩岩石物理模型方法
CN111027818B (zh) 一种页岩油分类评价方法
CN105527652A (zh) 一种岩石脆性的测井方法和装置
CN104632201A (zh) 一种页岩油气“甜点区”有机碳含量下限值的确定方法
CN106054279A (zh) 一种煤岩脆性指数的确定方法
CN108875115B (zh) 一种确定岩石强度的方法
CN109211745A (zh) 一种富含有机质泥页岩孔隙度的演化过程的恢复方法
CN112560246B (zh) 一种目标井散点地层压力系数的预测方法
CN111764895B (zh) 一种适用于页岩气储层地质模型的测井评价方法
CN110133724B (zh) 一种应用阵列声波测井数据进行气层识别的方法
CN107808055B (zh) 一种基于双重影响的页岩气饱和度测井计算方法
Evans et al. A geological approach to permeability prediction in clastic reservoirs
Huang et al. Geophysical evaluation technology for shale gas reservoir: A case study in Silurian of Changning Area in Sichuan Basin
Zhang et al. Staged fracturing of horizontal wells in continental tight sandstone oil reservoirs: a case study of Yanchang Formation in Western Ordos Basin, China
CN112814655A (zh) 一种校正页岩气储层束缚水饱和度的方法
CN112133376A (zh) 一种饱含水地层体积模量的流体识别方法
Zheng et al. Geological model evaluation through well test simulation: a case study from the Wytch farm oilfield, southern England

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant