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

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

Info

Publication number
CN108875115A
CN108875115A CN201810201388.XA CN201810201388A CN108875115A CN 108875115 A CN108875115 A CN 108875115A CN 201810201388 A CN201810201388 A CN 201810201388A CN 108875115 A CN108875115 A CN 108875115A
Authority
CN
China
Prior art keywords
rock
strength
content
parameter
weight coefficient
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.)
Granted
Application number
CN201810201388.XA
Other languages
English (en)
Other versions
CN108875115B (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 Petroleum and Natural Gas Co Ltd
China Petroleum Logging Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
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 Petroleum and Natural Gas Co Ltd, China Petroleum Logging Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
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

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)]
固有剪切强度:
单轴抗拉强度:
ST=SC/12
但在实际应用中,不同地区不同地层效果不尽相同,不能满足准确计算岩石强度参数的需求。本发明通过矿物含量对岩石强度参数进行拟合,建立了一种更为精确的岩石强度计算方法,满足复杂储层评价的需求。
发明内容
为解决现有技术中存在的上述缺陷,本发明的目的在于提供一种确定岩石强度的方法,该方法通过矿物含量与岩石强度参数之间建立关系,能够较为准确的计算储层的岩石强度参数,解决了由于原有计算模型由于缺少地区参数而无法较为准确的计算岩石气囊度参数的问题,能够及时服务或满足实际生产的需要,实现了油气储层可压裂性评价的需求。
本发明是通过下述技术方案来实现的。
所述方法按照以下步骤进行操作:
步骤(1),获取岩心刻度常规测井资料或ECS资料,建立参数模型,计算储层矿物含量,分析岩石强度参数与矿物含量之间的对应关系。
步骤(2),利用计算方程求取岩石单轴抗压强度和固有剪切强度参数;
步骤(3),利用步骤(2)得到的岩石强度参数计算方程,代入岩石矿物含量以及岩石力学参数,即可得到岩石强度参数;
步骤(4),根据求得的岩石强度参数结合测井资料连续评价地层特征。
上述技术方案中,还包括进一步限定的方案:
所述步骤(1),计算储层矿物含量按照下式计算:
砂质含量:
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 分别为地层水、岩石骨架、岩石的中子孔隙度,%;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)
其中:
式中: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资料,计算储层矿物含量,分析岩石强度参数与矿物含量之间的对应关系;
计算储层矿物含量按照下式计算:
砂质含量:
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 分别为地层水、岩石骨架、岩石的中子孔隙度,%;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)
其中:
式中: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是本发明计算的岩石强度参数计算效果图。从图中可以看出,求取的抗张强度值、抗压强度值及固有剪切强度值,与岩石力学试验数据对比可见,所计算的岩石强度参数相对误差率低,准确度高,应用效果好。
从上述表1可以看出,岩石强度参数中,抗张强度相对误差不大于10.81%,抗压强度相对误差不大于12.91%,固有剪切强度相对误差不大于16.25%。
本发明通过矿物含量与岩石强度参数之间建立关系,能够较为准确的计算储层的岩石强度参数,解决了由于原有计算模型由于缺少地区参数而无法较为准确的计算岩石气囊度参数的问题,相比于原有的计算方法,能够更准确的计算储层强度参数,与岩石力学实验测的岩石强度参数符合性更好,可以及时服务或满足测井资料解释评价的需要,应用效果好。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内一种。

Claims (5)

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

Cited By (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
MARLENE C.VILLENEUVE 等: "Effects of Grain Scale Heterogeneity on Rock Strength and the Chipping Process", 《INTERNATIONAL JOURNAL OF GEOMECHANICS》 *
张晋言 等: "利用测井资料评价泥页岩油气"五性"指标", 《测井技术》 *
李权 等: "基于多矿物模型分析的孔隙度计算方法在M地区中的应用", 《贵州师范大学学报(自然科学版)》 *
王俊骏 等: "柳沟庄泥云岩储层岩性与矿物组分特性研究", 《石化技术》 *

Cited By (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 中国铁道科学研究院集团有限公司铁道建筑研究所 一种适用于机制骨料卵石母岩性能快速评价方法

Also Published As

Publication number Publication date
CN108875115B (zh) 2022-05-10

Similar Documents

Publication Publication Date Title
CN105469159B (zh) 定量预测油气有利聚集区的方法
Datta-Gupta Streamline Simulation: A Technology Update (includes associated papers 71204 and 71764)
Bell et al. In situ stress magnitude and orientation estimates for Cretaceous coal-bearing strata beneath the plains area of central and southern Alberta
CN104698492B (zh) 一种计算异常地层压力的方法
CN106223941B (zh) 基于测井资料的页岩气储层有机孔隙度确定方法
CN106054279B (zh) 一种煤岩脆性指数的确定方法
CN104007485A (zh) 一种建立复杂孔隙泥页岩岩石物理模型方法
Turner Rock-socketed shafts for highway structure foundations
Hoseinpour et al. Determination of the mud weight window, optimum drilling trajectory, and wellbore stability using geomechanical parameters in one of the Iranian hydrocarbon reservoirs
Yang et al. Natural fractures and mechanical properties in a horn river shale core from well logs and hardness measurements
LeCain Air-injection testing in vertical boreholes in welded and nonwelded tuff, Yucca Mountain, Nevada
Zhang et al. Postfracturing permeability prediction for CBM well with the analysis of fracturing pressure decline
CN109458176A (zh) 碳酸盐岩储层压力的预测方法及其应用
CN108875115A (zh) 一种确定岩石强度的方法
Ocak Estimating the modulus of elasticity of the rock material from compressive strength and unit weight
CN101315428A (zh) 基于横波速度的高孔隙碎屑岩沉积地层的高分辨率勘探方法
Samani et al. Characterizing a heterogeneous aquifer by derivative analysis of pumping and recovery test data
CN108595834A (zh) 一种基于多地质因素的煤层顶底板动力分区评价方法
Carr et al. Characterization of a Vertically Compartmentalized Reservoir in a Supergiant Field, Burgan Formation, Greater Burgan Field, Kuwait; Part 2: Evolving Concepts of Reservoir Performance
Haws et al. Applications of pressure-interference data in reservoir characterization studies, Big Horn basin, Wyoming
Geldon Results and interpretation of preliminary aquifer tests in boreholes UE-25c# 1, UE-25c# 2, and UE-25c# 3, Yucca Mountain, Nye County, Nevada
Johnson Constructing a Niobrara reservoir model using outcrop and downhole data
Geldon Hydraulic tests of Miocene volcanic rocks at Yucca Mountain and Pahute Mesa and implications for groundwater flow in the Southwest Nevada Volcanic Field, Nevada and California
Hossaini Some aspects of the strength characteristics of intact and jointed rocks
CN113763194B (zh) 页岩单井含气丰度预测方法及系统

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