CN110487794A - 一种基于数字图像技术的岩石强度的均质度识别方法 - Google Patents

一种基于数字图像技术的岩石强度的均质度识别方法 Download PDF

Info

Publication number
CN110487794A
CN110487794A CN201910889764.3A CN201910889764A CN110487794A CN 110487794 A CN110487794 A CN 110487794A CN 201910889764 A CN201910889764 A CN 201910889764A CN 110487794 A CN110487794 A CN 110487794A
Authority
CN
China
Prior art keywords
homogeneity
rock
color
image
intensity
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
CN201910889764.3A
Other languages
English (en)
Other versions
CN110487794B (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.)
Angang Group Mining Co Ltd
Original Assignee
Angang Group Mining 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 Angang Group Mining Co Ltd filed Critical Angang Group Mining Co Ltd
Priority to CN201910889764.3A priority Critical patent/CN110487794B/zh
Publication of CN110487794A publication Critical patent/CN110487794A/zh
Application granted granted Critical
Publication of CN110487794B publication Critical patent/CN110487794B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/32Polishing; Etching
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及一种基于数字图像技术的岩石均质度识别方法,其特征在于,包括试样制作与图像采集,通过对岩石表面图像的拍摄,获取岩石表面颜色信息,利用颜色的分布特征求解颜色均质度,并拟合得出颜色均质度与强度均质度间的函数关系,并利用该函数关系求解岩石强度均质度。该方法建立了岩石表面颜色均质度与岩石强度均质度之间的拟合关系,基于图像识别技术可方便地计算得出同类型、同批次岩石试样的强度均质度。可避免采用室内基本力学试验测试方法所带来的操作复杂性,以及有效的降低强度均质度的计算误差。

Description

一种基于数字图像技术的岩石强度的均质度识别方法
技术领域
本发明属于图像处理技术领域,具体涉及一种基于数字图像技术的岩石强度的均质度识别方法。
背景技术
均质度是描述岩石材料内部均匀程度的重要参数,对岩石的破裂过程及力学参数等方面均有显著影响,也是考虑岩石非均质性数值模拟计算需要确定的重要参数之一。基于岩石介质微元体强度服从Weibull随机分布的假设,许多学者提出了均质度的计算方法。徐涛等基于岩石应变强度理论,推导出利用单轴压缩全应力应变曲线的峰值点确定均质度计算公式,杨圣奇等在此基础之上引入损伤比例系数,计算了大理岩的均质度。曹文贵基于现有岩石损伤软化统计本构模型研究,通过引进Mohr-Coulomb强度准则和岩石应力-应变全程曲线峰值应变与围压的关系建立起特定围压下均质度与围压的解析表达式。罗荣等将岩石细观单元随机参数的变异系数定义为岩石均质度,并推导了计算方法。这些研究成果促进了岩石损伤软化统计本构模型本身的完善,加深了人们对均质度的认识,但是这些方法都是通过间接手段测量岩石强度的均质度。到目前为止依然没有直接测量岩石的均质度的有效手段。
岩石是由成岩矿物胶结而成,在胶结过程中,可能出现一些微裂纹和缺陷。岩石胶结矿物的不同和岩石内部的微裂纹、缺陷的存在是造成岩石非均质性的最主要原因。在对岩石照片进行分析时,可以发现岩石不同的胶结矿物和岩石内部的微裂纹、缺陷对应着不同的颜色特征,并且不同的胶结矿物对应着不同的强度,因此基于数字图像技术的岩石均质度识别是可行的。基于以上推论,通过岩石表面的数字图像获取岩石表面颜色信息,利用颜色的分布特征求解颜色均质度,并拟合得出颜色均质度与强度均质度间的函数关系,利用该函数关系即可达到基于数字图像技术识别岩石强度均质度的目的。
发明内容
本发明的目的是提供一种基于数字图像技术的岩石强度的均质度识别方法,通过岩石表面的数字图像获取岩石表面颜色信息,利用颜色的分布特征求解颜色均质度,并拟合其与强度均质度间的函数关系,用该函数关系即可达到基于数字图像技术识别岩石强度均质度的目的。
本发明的目的是通过下述技术方案实现的:
本发明的一种基于数字图像技术的岩石强度的均质度识别方法,其特征在于,一种基于数字图像技术的岩石均质度识别方法,其特征在于,包括试样制作与图像采集,通过对岩石表面图像的拍摄,获取岩石表面颜色信息,利用颜色的分布特征求解颜色均质度,并拟合得出颜色均质度与强度均质度间的函数关系,并利用该函数关系求解岩石强度均质度,具体步骤如下:
步骤一:试样制作与图像采集,用于建立岩石图像颜色信息数据库;
步骤二:采用图像离散矩阵函数对岩石试件图像颜色信息进行提取;
步骤三:通过拟合获取试件图像颜色均质度m;
步骤四:拟合颜色均质度与强度均质度的函数关系,通过拟合关系函数,即可在已知岩石表面颜色均质度的情况下,求解任一同类岩石试件的强度均质度m'。
进一步地,所述的步骤一:试样制作与图像采集包括下列步骤:
1.1采用钻孔取芯的方式获取所需识别岩石的深部岩石试样,将现场钻取的岩芯表面呈螺纹状凹凸不平,岩芯也长短不一,上下表面不平整的岩芯进行二次加工,加工成长度一致的圆柱体试样,同时采用砂纸将试样表面打磨平整;
1.2二次加工后的圆柱体试样,采用多角度、同亮度平移拍摄,并将拍摄的多张图像拼接成一张包含有许多像素的图像,作为下一步颜色信息提取所用图像。
进一步地,所述的步骤二:采用图像离散矩阵函数对岩石试件图像颜色信息进行提取是;
将步骤1.2中和拍摄的多张图像拼接成一张包含有许多像素的图像作为一个像素矩阵看待,每一个像素都包含三个颜色分量,根据红色(R)、绿色(G)和蓝色(B)三基色,将图像离散为一个矩阵函数,即
式中:i为图像的像素高度,j为图像的像素宽度,
利用数字图像处理技术提取数字图像中颜色分量对应的像素个数,并绘制出岩石颜色分量(R)的分布曲线。
进一步地,所述步骤三:通过拟合获取试件图像颜色均质度m是:
对绘制的岩石颜色分量(R)的分布曲线进行以颜色均质度m为自变量的Weibull曲线拟合,选取拟合度最高的Weibull曲线的均质度m作为岩石颜色均质度m。
Weibull分布函数的形式为:
φ(α)=m/α0(α/α0)m-1exp(-(α/α0)m) (2)
式中:α为随机变量,对于颜色分量的分布曲线,其为颜色分量变量;m为统计分布函数的形态参数,也即均质度系数;α0是尺度参数,可通过随机变量α的均质E(α)及均质度m利用式α0=E(α)/Г(1+1/m)求出,其中Г()为gamma函数。
进一步地,所述步骤四:拟合颜色均质度与强度均质度m'的函数关系是:
基于岩石应变强度理论和岩石强度的随机统计分布假设,采用损伤力学理论,强度均质度m'可由下式计算:
m'=(ln(Eεcc))-1 (3)
式中:σc、εc分别为单轴加载条件下全程应力应变曲线的峰值点对应的应力及应变;E为弹性模量,
根据同一批岩石试样进行基本力学测试,获取σc、εc、E等参数,即可求出该批次岩石试样的强度均质度m'。
同现有技术相比,本发明的优点是:
该方法建立了岩石表面颜色均质度与岩石强度均质度之间的拟合关系,基于图像识别技术可方便地计算得出同类型、同批次岩石试样的强度均质度。可避免采用室内基本力学试验测试方法所带来的操作复杂性,以及有效的降低强度均质度的计算误差。
附图说明
图1a为大块岩石加工成尺寸为Φ50×100mm的岩石试件图。
图1b从不同的方向对试件表面进行拍摄,并将4张图像拼接成一个整体图。
图2某铁矿岩石试样颜色分布曲线及Weibull分布拟合曲线图。
图3为不同求解方法下均质度的拟合关系图。
具体实施方式
下面结合实施例对本发明作进一步说明。
如图1-图3所示,本发明的一种基于数字图像技术的岩石强度的均质度识别方法,其特征在于,包括试样制作与图像采集,通过对岩石表面图像的拍摄,获取岩石表面颜色信息,利用颜色的分布特征求解颜色均质度,并拟合得出颜色均质度与强度均质度间的函数关系,并利用该函数关系求解岩石强度均质度,具体步骤如下:
步骤一:试样制作;
1)采用钻孔取芯的方式获取所需识别岩石的深部岩石试样,将现场钻取的岩芯表面呈螺纹状凹凸不平,岩芯也长短不一,上下表面不平整的岩芯进行二次加工,加工成长度一致的圆柱体试样,同时采用砂纸将试样表面打磨平整;
2)二次加工后的圆柱体试样,采用多角度、同亮度平移拍摄,并将拍摄的多张图像拼接成一张包含有许多像素的图像,作为下一步颜色信息提取所用图像。
取某铁矿的大块岩石,工程现场普遍采用钻孔取芯的方式获取深部岩石试样,所以获取的岩石试样一般为圆柱体。现场钻取的岩芯表面多呈螺纹状凹凸不平,岩芯也长短不一,上下表面不平整,需要对现场取芯进行二次加工。一般将岩芯截成长度一致的圆柱体试样,同时采用砂纸将试样表面打磨平整,避免表面不平导致图像分析出现误差。针对于加工后的圆柱体,
以基于图像识别技术获取某铁矿岩石试样强度均质度为例,介绍图像采集方式:
将取某铁矿的大块岩石加工成尺寸为Φ50×100mm的岩石试件,如图1(a)。为了充分利用岩石试件表面的图像信息,从4个不同的方向对试件表面进行拍摄,并将4张图像拼接成一个整体,作为下一步颜色信息提取所用图像,如图1(b)所示。
步骤二:岩石试件图像颜色信息提取;
将步骤1中拼接的岩石表面图像,利用数字图像处理技术提取数字图像中颜色分量对应的像素个数,并绘制出颜色分量(R)的分布曲线,如图2中黑色曲线所示。
图像中包含有许多像素,可将图像本身作为一个像素矩阵看待。每一个像素都包含3个颜色分量,红色(R)、绿色(G)、蓝色(B),即三基色。因此图像可离散为一个矩阵函数,即
式中:i为图像的像素高度,j为图像的像素宽度。
利用数字图像处理技术提取数字图像中颜色分量对应的像素个数并绘制出颜色分量(R)的分布曲线。
步骤三3:通过拟合获取图像颜色均质度m;
图像颜色均质度m获取方法:对绘制的岩石颜色分量(R)的分布曲线进行以颜色均质度m为自变量的Weibull曲线拟合,选取拟合度最高的Weibull曲线的均质度m作为岩石颜色均质度m。
Weibull分布函数的形式为:
φ(α)=m/α0(α/α0)m-1exp(-(α/α0)m) (2)
式中:α为随机变量,对于颜色分量的分布曲线,其为颜色分量变量;m为统计分布函数的形态参数,也即均质度系数;α0是尺度参数,可通过随机变量α的均质E(α)及均质度m利用式α0=E(α)/Г(1+1/m)求出,其中Г()为gamma函数。
拟合结果如图2中红色曲线所示。其拟合系数R2高达0.98,说明颜色分量(R)呈Weibull分布,颜色均质度的拟合具有事实依据。按照以上步骤计算出其他岩石试件的颜色均质度m,部分结果如表1所示。
表1颜色均质度
步骤四:通过室内试验获取岩石试件的强度均质度m'
对同一批岩石试样进行基本力学测试,获取σc、εc、E等参数,根据式(5)可求出该批次岩石试样的强度均质度,
基于岩石应变强度理论和岩石强度的随机统计分布假设,采用损伤力学理论,强度均质度m'可由下式计算:
m'=(ln(Eεcc))-1 (3)
式中:σc、εc分别为单轴加载条件下全程应力应变曲线的峰值点对应的应力及应变;E为弹性模量,
步骤五:建立颜色均质度与强度均质度之间对应关系,通过该关系即可实现基于数字图像技术的岩石强度的均质度识别,而无需进行岩石试件基本力学测试。
图3为各岩石试件颜色均质度及强度均质度的对应关系图。采用指数函数拟合二者间的关系,拟合结果如下:
m'=0.069e0.6084m
通过上式即可在已知岩石表面颜色均质度m的情况下,求解任一取该铁矿的岩石试件的强度均质度m',而无需进行岩石试件基本力学测试。

Claims (5)

1.一种基于数字图像技术的岩石均质度识别方法,其特征在于,包括试样制作与图像采集,通过对岩石表面图像的拍摄,获取岩石表面颜色信息,利用颜色的分布特征求解颜色均质度,并拟合得出颜色均质度与强度均质度间的函数关系,并利用该函数关系求解岩石强度均质度,具体步骤如下:
步骤一:试样制作与图像采集,用于建立岩石图像颜色信息数据库;
步骤二:采用图像离散矩阵函数对岩石试件图像颜色信息进行提取;
步骤三:通过拟合获取试件图像颜色均质度m;
步骤四:拟合颜色均质度与强度均质度的函数关系,通过拟合关系函数,即可在已知岩石表面颜色均质度的情况下,求解任一同类岩石试件的强度均质度m'。
2.根据权利要求1所述的基于数字图像技术的岩石均质度识别方法,其特征在于,所述的步骤一:试样制作与图像采集包括下列步骤:
1.1 采用钻孔取芯的方式获取所需识别岩石的深部岩石试样,将现场钻取的岩芯表面呈螺纹状凹凸不平,岩芯也长短不一,上下表面不平整的岩芯进行二次加工,加工成长度一致的圆柱体试样,同时采用砂纸将试样表面打磨平整;
1.2 二次加工后的圆柱体试样,采用多角度、同亮度平移拍摄,并将拍摄的多张图像拼接成一张包含有许多像素的图像,作为下一步颜色信息提取所用图像。
3.根据权利要求2所述的基于数字图像技术的岩石均质度识别方法,其特征在于,所述的步骤二:采用图像离散矩阵函数对岩石试件图像颜色信息进行提取是;
将步骤1.2中和拍摄的多张图像拼接成一张包含有许多像素的图像作为一个像素矩阵看待,每一个像素都包含三个颜色分量,根据红色(R)、绿色(G)和蓝色(B)三基色,将图像离散为一个矩阵函数,即
式中:i为图像的像素高度,j为图像的像素宽度,
利用数字图像处理技术提取数字图像中颜色分量对应的像素个数,并绘制出岩石颜色分量(R)的分布曲线。
4.根据权利要求3所述的基于数字图像技术的岩石均质度识别方法,其特征在于,所述步骤三:通过拟合获取试件图像颜色均质度m是:
对绘制的岩石颜色分量(R)的分布曲线进行以颜色均质度m为自变量的Weibull曲线拟合,选取拟合度最高的Weibull曲线的均质度m作为岩石颜色均质度m。
Weibull分布函数的形式为:
φ(α)=m/α0(α/α0)m-1exp(-(α/α0)m) (2)
式中:α为随机变量,对于颜色分量的分布曲线,其为颜色分量变量;m为统计分布函数的形态参数,也即均质度系数;α0是尺度参数,可通过随机变量α的均质E(α)及均质度m利用式α0=E(α)/Г(1+1/m)求出,其中Г()为gamma函数。
5.根据权利要求1所述的基于数字图像技术的岩石均质度识别方法,其特征在于,所述步骤四:拟合颜色均质度与强度均质度m'的函数关系是:
基于岩石应变强度理论和岩石强度的随机统计分布假设,采用损伤力学理论,强度均质度m'可由下式计算:
m'=(ln(Eεcc))-1 (3)
式中:σc、εc分别为单轴加载条件下全程应力应变曲线的峰值点对应的应力及应变;E为弹性模量,
根据同一批岩石试样进行基本力学测试,获取σc、εc、E等参数,即可求出该批次岩石试样的强度均质度m'。
CN201910889764.3A 2019-09-20 2019-09-20 一种基于数字图像技术的岩石强度的均质度识别方法 Active CN110487794B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910889764.3A CN110487794B (zh) 2019-09-20 2019-09-20 一种基于数字图像技术的岩石强度的均质度识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910889764.3A CN110487794B (zh) 2019-09-20 2019-09-20 一种基于数字图像技术的岩石强度的均质度识别方法

Publications (2)

Publication Number Publication Date
CN110487794A true CN110487794A (zh) 2019-11-22
CN110487794B CN110487794B (zh) 2021-10-22

Family

ID=68558794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910889764.3A Active CN110487794B (zh) 2019-09-20 2019-09-20 一种基于数字图像技术的岩石强度的均质度识别方法

Country Status (1)

Country Link
CN (1) CN110487794B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112798592A (zh) * 2020-12-28 2021-05-14 山东大学 基于岩相学特征分析的岩石强度预测系统及方法
CN114399486A (zh) * 2022-01-11 2022-04-26 上海飞机制造有限公司 一种用于获取喷丸工艺中的弹丸散落分布函数的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130308831A1 (en) * 2012-05-18 2013-11-21 Ingrain, Inc. Method And System For Estimating Rock Properties From Rock Samples Using Digital Rock Physics Imaging
CN104849276A (zh) * 2015-05-20 2015-08-19 辽宁工程技术大学 一种基于像素统计的花岗岩三维细观结构重构方法
CN106501206A (zh) * 2016-01-21 2017-03-15 中国石油化工股份有限公司胜利油田分公司投资发展处 一种基于显微红外光谱技术的储层岩石非均质性表征方法
CN107782640A (zh) * 2017-10-20 2018-03-09 东北大学 基于数字图像进行岩石试件含水均匀性检测和扩散系数计算的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130308831A1 (en) * 2012-05-18 2013-11-21 Ingrain, Inc. Method And System For Estimating Rock Properties From Rock Samples Using Digital Rock Physics Imaging
CN104849276A (zh) * 2015-05-20 2015-08-19 辽宁工程技术大学 一种基于像素统计的花岗岩三维细观结构重构方法
CN106501206A (zh) * 2016-01-21 2017-03-15 中国石油化工股份有限公司胜利油田分公司投资发展处 一种基于显微红外光谱技术的储层岩石非均质性表征方法
CN107782640A (zh) * 2017-10-20 2018-03-09 东北大学 基于数字图像进行岩石试件含水均匀性检测和扩散系数计算的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
朱泽奇等: "基于数字图像处理的非均质岩石材料破坏过程模拟", 《岩土力学》 *
李昂等: "软硬互层状岩体局部区域的细观非均匀特性研究", 《岩石力学与工程学报》 *
罗荣等: "岩石非均质度对其力学性能的影响研究", 《岩土力学》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112798592A (zh) * 2020-12-28 2021-05-14 山东大学 基于岩相学特征分析的岩石强度预测系统及方法
CN114399486A (zh) * 2022-01-11 2022-04-26 上海飞机制造有限公司 一种用于获取喷丸工艺中的弹丸散落分布函数的方法

Also Published As

Publication number Publication date
CN110487794B (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN103969168B (zh) 一种疏松矿物横截面孔隙率的定量测定方法
CN102012356B (zh) 一种沥青混凝土路面集料级配的快速检测方法
CN109359416B (zh) 一种反映真实工程土石混合体分布的颗粒流数值模拟方法
Polat et al. The correlation between aggregate shape and compressive strength of concrete: Digital image processing approach
CN110487794A (zh) 一种基于数字图像技术的岩石强度的均质度识别方法
CN108956424A (zh) 一种页岩中孔隙定量表征的方法
CN102749046B (zh) 岩体结构面直剪试验中剪切面积的测量方法
CN104700097B (zh) 一种基于矿物标准色数据库的矿物组分识别方法
Reedy et al. Image analysis in quantitative particle studies of archaeological ceramic thin sections
Song et al. Advances in measuring air-void parameters in hardened concrete using a flatbed scanner
CN100596330C (zh) 冶金矿相智能识别方法
CN113158490B (zh) 致密砂岩储层渗透率计算模型的建立方法
CN110487696A (zh) 一种致密砂岩全视域的孔隙、物性检测方法和装置
Guo et al. Texture analysis of the microstructure of internal curing concrete based on image recognition technology
WO2023246459A1 (zh) 一种确定软岩内摩擦角的方法
Wei et al. The colorimetry method in assessing fire-damaged concrete
CN106570269A (zh) 一种基于数字图像的页岩破裂模式分析方法
CN111307070B (zh) 基于数字图像处理对混凝土粗集料的棱角度测量方法
CN111766304B (zh) 一种基于压缩试验判断脆性岩石宏细观行为联系的方法
CN116502557B (zh) 基于网格分析法的储层压裂后可动流体饱和度计算方法
Santos et al. Study of size effect using digital image correlation
CN116291330A (zh) 一种基于岩石破碎块质量分形的可压裂性评价方法
Kukela et al. Simplified method of assessment of weathering on historical stone monuments: An example of El-Merdani mosque, Cairo, Egypt
CN110059648A (zh) 一种基于增补分割的复杂组成烧结矿的非结构信息库的建立方法
Roy et al. Shape parameters of recycled aggregate and determination its area of mortar and aggregate using imaging techniques

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