CN113009592A - 一种砾岩地层岩石研磨性参数的评价方法和校正方法 - Google Patents

Abstract

一种砾岩地层岩石研磨性参数的评价方法和校正方法。所述岩石研磨性参数的评价方法，包括：根据岩样的测井数据，计算所述岩样的研磨性参数；所述测井数据包括密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据和声波测井数据。本发明利用密度测井资料、电阻率测井资料、中子孔隙度测井资料、声波测井资料和自然伽马测井资料直接评价砾岩地层的岩石研磨性参数，计算简便、成本低，对于砾岩地层钻头选型、高效破岩钻进、钻井工程降本增效具有重要意义。

Description

一种砾岩地层岩石研磨性参数的评价方法和校正方法

技术领域

本发明涉及地质勘探领域，具体的涉及了一种砾岩地层岩石研磨性参数的评价方法和校正方法。

背景技术

在钻井过程中，钻头受到岩石的磨损而变钝，从而在使用一段时间之后需要更换。岩石的研磨性参数是定量反映岩石磨损钻头切削齿能力的重要参数。目前随着砾岩油藏的开发，尤其是深部砾岩地层强烈的非均质性将导致钻头受到剧烈、不规则载荷的冲击，不同砾石体积含量、粒径、组分的砾岩研磨性差异较大，常用的CAI评价方法及等效石英含量法很难在砾岩地层适用，但岩石研磨性是影响钻头磨损的一个重要因素，对于钻头选型、钻井参数选择、降低钻井成本具有重要意义。

发明内容

本发明实施例的目的是一种砾岩地层岩石研磨性参数的评价方法和校正方法。岩石研磨性是影响钻头磨损的一个重要因素。本发明利用密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据和声波测井数据准确地评价砾岩地层的岩石研磨性，计算简便、成本低、便于优选钻头，提高钻井效率。

为了实现上述目的，本发明实施例提供了一种砾岩地层岩石研磨性参数的评价方法和校正方法，该方法包括：包括：根据岩样的测井数据，计算所述岩样的研磨性参数；所述测井数据包括密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据和声波测井数据。

可选的，所述密度测井数据包括密度测井值、砾岩地层密度测井值最大值、砾岩地层密度测井值最小值；所述声波测井数据包括声波时差测井值、声波时差测井值最大值、声波时差井值最小值；所述电阻率测井数据包括电阻率测井值；所述中子孔隙度测井数据包括中子孔隙度测井值；所述自然伽马测井数据包括自然伽马测井值。

可选的，所述根据岩样的测井数据，计算所述岩样的研磨性参数，包括：根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量；根据所述岩样的中子孔隙度测井值和电阻率测井值，计算砾石粒径中值；根据所述岩样的砾石体积含量、砾石粒径中值、声波时差测井值和自然伽马测井值，计算所述岩样的研磨性参数。

可选的，所述根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量，包括：根据所述岩样的密度测井数据，计算归一化密度测井值；根据所述岩样的声波测井数据，计算归一化声波时差测井值；根据所述归一化密度测井值和所述归一化声波时差测井值，计算砾石体积含量。

可选的，所述根据所述岩样的密度测井数据，计算归一化密度测井值，包括

其中，ρ*为归一化密度测井值，ρ为密度测井值，ρ_max为砾岩地层密度测井最大值，ρ_min为砾岩地层密度测井最小值。

可选的，所述根据所述岩样的声波测井数据，计算归一化声波时差测井值，还包括

其中，Δt*为归一化声波时差测井值，

Δt为声波时差测井值，

Δt_max为声波时差测井最大值，

Δt_min为声波时差测井最小值。

可选的，所述根据所述归一化密度测井值和所述归一化声波时差测井值，计算砾石体积含量，包括

V＝a·ρ*-b·Δt*+c

其中，V为砾石体积含量，

a，b，c为回归系数，

p*为归一化密度测井值，

Δt*为归一化声波时差测井值。

可选的，所述根据所述岩样的中子孔隙度测井值和深电阻率测井值，计算砾石粒径中值，包括

其中，D为砾石粒径中值，

d，e为回归系数，

RT为深电阻率测井值，

Φ为中子孔隙度测井值。

可选的，所述根据所述岩样的砾石体积含量、砾石粒径中值、声波时差测井值和自然伽马测井值，计算所述岩样的研磨性参数，包括

其中，ω为岩石的研磨性参数，V为砾石体积含量，D为砾石粒径中值，f，g为回归系数，Δt为声波时差测井值，GR为自然伽马测井值。

相应的，本发明实施例还提供一种砾岩地层岩石研磨性参数的校正方法，包括：根据砾岩地层岩石研磨性参数的评价方法，获取岩样的研磨性参数测量值；利用多个岩样相对应的多个岩样的研磨性参数测量值，拟合形成岩样的研磨性参数校正值与岩样的研磨性参数测量值之间的线性关系式；将岩样的研磨性参数测量值代入所述线性关系式，得到岩样的研磨性参数校正值。

通过上述技术方案，本发明通过根据岩样的测井数据，计算所述岩样的研磨性参数；所述测井数据包括密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据和声波测井数据。该方法直接利用密度测井资料、电阻率测井资料、中子孔隙度测井资料、自然伽马测井资料、声波测井资料评价砾岩地层岩石研磨性参数，计算简便，学习成本极低，推广性强，计算模型得出的研磨性参数与实际砾岩研磨性参数高度相关，评价结果准确可靠，对于砾岩钻头选型、高效破岩钻进、钻井工程降本增效具有重要意义。

本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：

图1和图2是本发明的一种砾岩地层岩石研磨性参数的评价方法的流程示意图；

图3是本发明的砾岩地层岩石研磨性参数的示意图。

具体实施方式

以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。

图1是本发明的一种砾岩地层岩石研磨性参数的评价方法的流程示意图。根据岩样的测井数据，计算所述岩样的研磨性参数；所述测井数据包括密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据和声波测井数据。如图1所示，本发明实施例的砾岩地层岩石研磨性参数的评价方法可以包括以下步骤：

步骤S101，获取测井数据中的密度测井数据。密度测井为根据伽马射线与地层的康普顿效应测定地层密度的测井方法，密度测井是划分煤层、划分致密岩层中的裂隙带，以及研究渗透性岩层的孔隙度的有效方法。根据密度测井资料获取密度测井值，用于对岩样的研磨性极值进行分析计算。所述密度测井数据包括密度测井值、砾岩地层密度测井值最大值、砾岩地层密度测井值最小值。

步骤S102，获取测井数据中的电阻率测井数据。电阻率测井是地球物理测井中常见的测井方法，它根据岩石导电性的差别，测量底层的电阻率，岩石的电阻率与岩性、储集物性、含油性有密切的关系，可以通过研究岩石的电阻率的差异区分岩性、划分储集层并评价含油性、进行底层对比，为井内研究钻井地质剖面特性的重要参数。电阻率测井是以岩石、矿石电性为基础的一组测井方法，包括普通电阻率测井、微电极、聚焦测井(深、浅三侧向、深、浅七侧向、双侧向、微侧向等)等测井方法。所述电阻率测井数据包括深电阻率测井值。根据电阻率测井资料获取电阻率测井值，用于对岩样研磨性参数进行分析计算。

步骤S103，获取测井数据中的中子孔隙度测井数据。中子孔隙度是用在中子刻度井中刻度过的中子测井仪器测出的地层孔隙度，实质上是等效含氢指数。用在中子刻度井中刻度过的中子测井仪器测出的地层孔隙度，如实际孔隙度为零的石膏的中子孔隙度为49％。所述中子孔隙度测井数据包括中子孔隙度测井值。

步骤S104，获取测井数据中的声波测井数据。声波在不同介质中传播时，速度、幅度及频率的变化等声学特性也不相同。声波测井就是利用岩石的这些声学性质来研究钻井的地质剖面，判断固井质量的一种测井方法。其中可以通过声波速度测井，得到孔隙度和力学参数；通过声幅测井，研究固井质量；通过声波电视测井，观察井壁情况、裂缝；通过噪声测井，了解井下流体的流动。声波在岩石中的传播方式包括纵波和横波，在井下，纵波和横波都能在底层传播，而泥浆中只能传播纵波。所述声波测井数据包括声波时差测井值、声波时差测井值最大值、声波时差井值最小值。

步骤S105，获取测井数据中的自然伽马测井数据。自然伽马测井是沿井身测量岩层的天然伽马射线强度的方法。岩石一般都含有不同数量的放射性元素，并且不断地放出射线。测井结果就有可能划分出钻孔的地质剖面、确定砂泥岩剖面中砂岩泥质含量和定性地判断岩层的渗透性。伽马射线是原子核衰变裂解时释放的射线之一，穿透能力极强，从液体到金属的大部分物质都能穿过，正是由于具有此特性，使其在石油工业等方面得到广泛的应用。岩石中主要含有铀(U)、钍(Th)、钾(K)等放射性元素，在沉积岩中这些放射性元素主要反映泥质含量的变化，在火山岩、花岗性风化层及某些盐类沉积，自然伽马测量值显著增高，常做为识别这类岩石类型的重要曲线标志。

步骤S106，根据岩样的测井数据，计算所述岩样的研磨性参数。所述测井数据包括密度测井数据、电阻率测井数据、中子孔隙度测井数据、声波测井数据和自然伽马测井数据。在钻进过程中，钻头在轴压和转速的作用下破碎岩石，与此同时，钻头本身也受岩石的研磨面变钝。岩石磨损钻头的能力称为岩石的研磨性。钻头被磨损，不仅增加了钻头的消耗，同时也降低了碎岩的效率。岩石研磨性与钻头寿命、生产效率、钻探成本直接相关，所以是选择钻头、设计钻头、确定规程参数、制定生产定额的重要参数。

图2是步骤S105的一种具体实施方式。按照该具体实施方式，计算所述岩样的研磨性参数，包括：

根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量；根据所述岩样的中子孔隙度测井值和电阻率测井值，计算砾石粒径中值；根据所述岩样的砾石体积含量、砾石粒径中值、声波时差测井值和自然伽马测井值，计算所述岩样的研磨性参数。

步骤S201，获取密度测井数据。所述密度测井数据包括密度测井值、砾岩地层密度测井值最大值、砾岩地层密度测井值最小值。所述根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量，包括：根据所述岩样的密度测井数据，计算归一化密度测井值，包括

其中，ρ*为归一化密度测井值，ρ为密度测井值，ρ_max为砾岩地层密度测井最大值，ρ_min为砾岩地层密度测井最小值。

步骤S202，获取声波测井数据。所述声波测井数据包括声波时差测井值、声波时差测井值最大值、声波时差井值最小值。所述根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量，包括：根据所述岩样的声波测井数据，计算归一化声波时差测井值，还包括

其中，Δt*为归一化声波时差测井值，Δt为声波时差测井值，Δt_max为声波时差测井最大值，Δt_min为声波时差测井最小值。

步骤S203，根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量。包括：根据所述归一化密度测井值和所述归一化声波时差测井值，计算砾石体积含量，包括

V＝a·ρ*-b·Δt*+c

其中，V为砾石体积含量，a，b，c为回归系数，p*为归一化密度测井值，Δt*为归一化声波时差测井值。

步骤S204，根据所述岩样的中子孔隙度测井值和深电阻率测井值，计算砾石粒径中值，包括

其中，D为砾石粒径中值，d，e为回归系数，RT为深电阻率测井值，Φ为中子孔隙度测井值。

步骤S205，根据所述岩样的砾石体积含量、砾石粒径中值、声波时差测井值和自然伽马测井值，计算所述岩样的研磨性参数，包括

其中，ω为岩石的研磨性参数，V为砾石体积含量，D为砾石粒径中值，f，g为回归系数，Δt为声波时差测井值，GR为自然伽马测井值。图3是本发明的砾岩地层岩石研磨性参数的示意图。如图所示，直接利用密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据、声波测井数据评价砾岩地层岩石研磨性参数，直观的反应出砾岩地层岩石的研磨性分布，对于砾岩地层钻头选型、高效破岩钻进、钻井工程降本增效具有重要意义。

本申请人主要通过观察砾岩岩心样品及统计方法获取砾岩样品的砾石体积含量及粒径中值，并取得工程资料中密度测井资料、电阻率测井资料、中子孔隙度测井资料、声波测井资料、自然伽马测井资料，利用多元回归的方法分析建立数学模型，实现对砾岩地层岩石研磨性参数的测定。具体步骤为：

步骤1，获取直径为100mm的干燥砾岩样品。所述砾岩样品为圆柱体砾岩岩心。干燥处理可选将所述岩样放入110℃的烤箱烘烤24h。

步骤2，将所述岩样放入岩石研磨性参数测量装置，每端面测量3组岩心的岩石研磨性参数级值，研磨环损失体积及岩石损失体积之比即为研磨性参数，并记录被研磨体积的砾石体积含量及粒径中值。获得得岩石研磨性参数、砾石体积含量、砾石粒径中值结果如下表1所示：

表1

步骤3，根据获取砾岩样品的砾石体积含量数据，建立砾石体积含量和密度测井值及声波时差测井值之间的计算模型；具体方法如下：

利用多元回归的方法分析砾石体积含量与密度测井值及声波时差测井值，建立如下数学模型：

其中，V为砾石体积含量，ρ为密度测井值，ρmax为砾岩地层密度测井最大值，ρmin为砾岩地层密度测井最小值，Δt为声波时差测井值，Δtmax为声波时差测井最大值，Δtmin为声波时差测井最小值，回归系数a，b为1.15，c为-0.014。

步骤4，根据获取砾岩样品的砾石粒径中值，建立砾岩样品砾石粒径中值与对应中子孔隙度测井值、电阻率测井值之间的计算模型。

具体方法如下：利用多元非线性回归方法分析砾石粒径中值与中子孔隙度测井值及电阻率测井值，建立如下数学模型：

其中，D为砾石粒径中值，RT为深电阻率测井值，Φ为中子测井值，回归系数d为38，e为5。

步骤5，建立砾石体积含量、砾石粒径中值及声波时差测井值与砾岩岩石研磨性参数级值之间的关系模型，利用多元非线性回归方法分析得到砾岩岩石研磨性参数计算模型：

其中，ω为岩石研磨性参数，GR为自然伽马测井值，Δt为声波时差测井值，回归系数f为0.0573，g为0.59。

上述各公式各系数取决于区块砾岩特性及多元回归数据，并非唯一的取值。

根据所确定的砾岩地层的砾石体积含量、砾石粒径中值及测井资料中的声波时差测井值、自然伽马测井值，带入上式即可得到考虑砾岩地层的岩石研磨性参数。相对于现有的针对均质性岩石的岩石研磨性评价方法，本发明计算出了砾岩的砾石体积含量及砾石粒径中值，考虑了砾石对岩石研磨性的影响，准确性强。本发明能够直接利用测井资料评价砾岩地层岩石研磨性参数，计算简便，学习成本极低，推广性强，计算模型得出的研磨性参数与实际砾岩研磨性参数高度相关，评价结果准确可靠，对于砾岩钻头选型、高效破岩钻进、钻井工程降本增效具有重要意义。

本发明还提供了一种砾岩地层岩石研磨性参数的校正方法，包括：根据上述砾岩地层岩石研磨性参数的评价方法，获取岩样的研磨性参数测量值；利用多个岩样相对应的多个岩样的研磨性参数测量值，拟合形成岩样的研磨性参数校正值与岩样的研磨性参数测量值之间的线性关系式；将岩样的研磨性参数测量值代入所述线性关系式，得到岩样的研磨性参数校正值。实现对岩样的研磨性参数进行校正。

以上结合附图详细描述了本发明实施例的可选实施方式，但是，本发明实施例并不限于上述实施方式中的具体细节，在本发明实施例的技术构思范围内，可以对本发明实施例的技术方案进行多种简单变型，这些简单变型均属于本发明实施例的保护范围。

另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施例对各种可能的组合方式不再另行说明。

此外，本发明实施例的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施例的思想，其同样应当视为本发明实施例所公开的内容。

Claims (10)

1.一种岩石研磨性参数的评价方法，用于评价砾岩地层，其特征在于，包括：

根据岩样的测井数据，计算所述岩样的研磨性参数；

所述测井数据包括密度测井数据、电阻率测井数据、中子孔隙度测井数据、自然伽马测井数据和声波测井数据。

2.根据权利要求1所述的方法，其特征在于，

所述密度测井数据包括密度测井值、砾岩地层密度测井值最大值、砾岩地层密度测井值最小值；

所述声波测井数据包括声波时差测井值、声波时差测井值最大值、声波时差井值最小值；

所述电阻率测井数据包括电阻率测井值；

所述中子孔隙度测井数据包括中子孔隙度测井值；

所述自然伽马测井数据包括自然伽马测井值。

3.根据权利要求2所述的方法，其特征在于，所述根据岩样的测井数据，计算所述岩样的研磨性参数，包括：

根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量；

根据所述岩样的中子孔隙度测井值和电阻率测井值，计算砾石粒径中值；

根据所述岩样的砾石体积含量、砾石粒径中值、声波时差测井值和自然伽马测井值，计算所述岩样的研磨性参数。

4.根据权利要求3所述的方法，其特征在于，所述根据所述岩样的密度测井数据和声波测井数据，计算砾石体积含量，包括：

根据所述岩样的密度测井数据，计算归一化密度测井值；

根据所述岩样的声波测井数据，计算归一化声波时差测井值；

根据所述归一化密度测井值和所述归一化声波时差测井值，计算砾石体积含量。

5.根据权利要求4所述的方法，其特征在于，所述根据所述岩样的密度测井数据，计算归一化密度测井值，包括

其中，ρ*为归一化密度测井值，

ρ为密度测井值，

ρ_max为砾岩地层密度测井最大值，

ρ_min为砾岩地层密度测井最小值。

6.根据权利要求4所述的方法，其特征在于，所述根据所述岩样的声波测井数据，计算归一化声波时差测井值，还包括

其中，Δt*为归一化声波时差测井值，

Δt为声波时差测井值，

Δt_max为声波时差测井最大值，

Δt_min为声波时差测井最小值。

7.根据权利要求4所述的方法，其特征在于，所述根据所述归一化密度测井值和所述归一化声波时差测井值，计算砾石体积含量，包括

V＝a·ρ*-b·Δt*+c

其中，V为砾石体积含量，

a，b，c为回归系数，

ρ*为归一化密度测井值，

Δt*为归一化声波时差测井值。

8.根据权利要求3所述的方法，其特征在于，所述根据所述岩样的中子孔隙度测井值和电阻率测井值，计算砾石粒径中值，包括

其中，D为砾石粒径中值，

d，e为回归系数，

RT为电阻率测井值，

Φ为中子孔隙度测井值。

9.根据权利要求3所述的方法，其特征在于，所述根据所述岩样的砾石体积含量、砾石粒径中值、声波时差测井值和自然伽马测井值，计算所述岩样的研磨性参数，包括

其中，ω为岩石的研磨性参数，

V为砾石体积含量，

D为砾石粒径中值，

f，g为回归系数，

Δt为声波时差测井值，

GR为自然伽马测井值。

10.一种砾岩地层岩石研磨性参数的校正方法，其特征在于，包括：根据权利要求1至9中任一权利要求所述的砾岩地层岩石研磨性参数的评价方法，获取岩样的研磨性参数测量值；

利用多个岩样相对应的多个岩样的研磨性参数测量值，拟合形成岩样的研磨性参数校正值与岩样的研磨性参数测量值之间的线性关系式；

将岩样的研磨性参数测量值代入所述线性关系式，得到岩样的研磨性参数校正值。

Images