一种基于CT扫描的煤岩逆向化建模力学参数确定方法技术

本发明专利技术公开了一种基于CT扫描的煤岩逆向化建模力学参数确定方法，运用游标卡尺及电子称测出煤岩试件的基本尺寸a及煤岩试件质量m并求出煤岩试件体积V，利用公式V＝a3、计算出煤岩试件密度，采用超声波探伤仪探测出穿过煤岩试件横波和纵波用时ts、tp，并利用公式求出煤岩试件弹性模量E，CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件对煤岩试件内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi，得出多个煤岩试件基本参数后，运用公式Mi＝mi1+mi2+……+mj＝ρi1Vi1+ρi2Vi2+……+ρijV得出不同组份材料密度ρij，最终运用公式得出不同组份材料的弹性模量Ei，达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。

【技术实现步骤摘要】

本专利技术属于数学计算方法领域，涉及煤岩逆向化建模力学参数确定方法，尤其涉 及一种基于CT扫描的煤岩逆向化建模力学参数确定方法。
技术介绍
煤岩在细观尺度上是非均匀材料，其物理力学性质在空间上是非连续性。长期以 来的理论研宄是建立在均质性的假设之上研宄岩石的力学性质及变化过程，而非均质性对 于煤岩的强度、应力分布及破坏形式起着重要的作用。在煤岩非均匀性理论分析上，常采用 统计分布规律来描述。在对煤岩破坏形式的数值计算中，大部分学者是在煤岩非均质性理 论分析的基础之上来探讨非均质性对煤岩破坏形式的影响，例如Weibull分布以及以此为 基础探讨煤岩的非均质性对其破坏形式影响的数值计算。介于煤岩试件中所含的矿物成份 比较复杂，为了真实的展现煤岩试件材料的非均质性，一部分学者通过工业CT及逆向化三 维重建软件来实现煤岩试件非均质性的真实展现。但是，目前基于逆向化来探讨非均质性 对煤岩破裂形式影响的数值计算，并不能够准确获取各组份材料的参数，数值计算结果的 准确性有待考证。基于上述情况，迫切需要一种简单可靠，基于CT扫描的煤岩逆向化建模 力学参数确定方法，有效可提高煤岩内部裂表征精度的基于CT扫描的受载煤岩损伤本构 方程构建方法，达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。
技术实现思路
本专利技术运用游标卡尺及电子称测出煤岩试件的基本尺寸a及煤岩试件质量m并求 出煤岩试件体积V，利用公式V=a3、P= $计算出煤岩试件密度，采用超声波探伤仪探测出 穿过煤岩试件横波和纵波用时ts、tp，并利用公式求出煤岩试件弹性模量E，CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软 件对煤岩试件内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi，得出多个 煤岩试件基本参数后，运用公式Mi=mn+mi2+......+mj=p nVn+p i2Vi2+......+P i」V得出不 同组份材料密度Pw最终运用公5导出不同组份材料的弹 性模量Ep达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。为了 实现上述目的，本专利技术采用如下技术方案： 一种基于CT扫描的煤岩逆向化建模力学参数确定方法，其特征在于，包括如下步 骤： (1)运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤 岩试件体积V，结合基本的计算公V=a3、P= ^7，计算得出煤岩试件的密度P。 (2)采用超声波探伤仪，对煤岩试件进行超声波实验，分别以横波、纵波穿过煤岩 试件。记录下横波及纵波穿过煤岩试件的时间，时间分别为：ts、tp，并计算出横波及纵波波 速戈W用公式得出煤岩试件弹性模量E。 (3)利用CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件利用二 维灰度图像构建三维重建模型，对三维重建模型内部不同矿物成份的材料进行分组以及统 计不同组份材料体积Vi。 (4)结合煤岩试件质量m、不同组份材料体积Vi及煤岩试件弹性模量E，首先运用 公式Mi=mn+mi2+......+mj=p nVn+p i2Vi2+......+P 计算出不同组份材料的密度P "， 得出不同组份材料密度Pw结合公3:出不同组份材料 的弹性模量Ep 与现有相关方法相比，本专利技术具有如下优点： 本专利技术提供了一种新的基于CT扫描的煤岩逆向化建模力学参数确定方法，测量 出煤岩试件基本尺寸、煤岩试件质量并求出煤岩试件体积，然后对煤岩试件进行超声波实 验，利用相关公式推导出煤岩试件不同组份材料的参数，为基于CT扫描的逆向化三维重构 煤岩试件数值计算提供精准参数。【附图说明】 图1为本专利技术的煤岩试件体积求解过程图； 图2为本专利技术的煤岩试件质量求解过程图； 图3为本专利技术的煤岩试件弹性模量求解过程图； 图4为本专利技术的不同组份材料体积求解过程图。 图中：1_煤岩试件；2-CT扫描系统；3-二维灰度图像；4-三维重建软件；5-三维 重建模型；6-不同组份材料体积Vi;7_游标卡尺；8-煤岩试件体积V;9_电子称；10-煤岩 试件质量m;11_超声波探伤仪；12-煤岩试件弹性模量E。【具体实施方式】 结合图1、图2、图3、图4所示，一种基于CT扫描的煤岩逆向化建模力学参数确 定方法，包括煤岩试件1、CT扫描系统2、二维灰度图像、三维重建软件4、三维重建模型 5、不同组份材料体积Vi6、游标卡尺7、煤岩试件体积、电子称9、煤岩试件质量mlO、超声 波探伤仪11、煤岩试件弹性模量E12,运用游标卡尺7及电子称9测出煤岩试件1的基本 尺寸及煤岩试件质量mlO并求出煤岩试件体积V8,利用公式V=a3、P= ^算出煤岩试 件密度，采用超声波探伤仪11探测出穿过煤岩试件1横波和纵波用时ts、tp，并利用公式十算出煤岩试件弹性模量E12,CT扫描系统2扫描 煤岩试件1得出二维灰度图像3,运用三维重建软件4对煤岩试件1内部不同矿物成份的材 料进行分组以及统计不同组份材料体积Vi6,得出多个煤岩试件1基本参数后，运用公式Mi =mil+mi2+......+mj=pnVn+pi2Vi2+......+PUV得出不同组份材料密度p最终运用公式辱出不同组份材料的弹性模量Ep 以某矿立方体非标准煤岩试件为例，阐述本专利技术具体步骤如下： (1)运用游标卡尺7及电子称9测出煤岩试件1的基本尺寸a、煤岩试件质量mlO 并求出煤岩试件体积V8,结合基本的计算公式V=a3、P= ^，计算得出煤岩试件1的密度 P〇 (2)采用超声波探伤仪11，对煤岩试件1进行超声波实验，分别以横波、纵波穿过 煤岩试件1。记录下横波及纵波穿过煤岩试件1的时间，时间分别为：ts、tp，并计算出横波 及纵波波速为I岩试件弹性模量E12。 (3)利用CT扫描系统2扫描煤岩试件1得出二维灰度图像3,运用三维重建软件 4利用二维灰度图像3构建三维重建模型5,对三维重建模型5内部不同矿物成份的材料进 行分组以及统计不同组份材料体积Vi6。 (4)结合煤岩试件质量mlO、不同组份材料体积Vi6及煤岩试件弹性模量E12,首 先运用公式Mi=mn+mi2+......+mj=pnVn+pi2Vi2+......+PUV计算出不同组份材料的密 度Pw得出不同组份材料密度Pw结合公3:出不同组 份材料的弹性模量Ep 当然，以上说明仅仅为本专利技术的较佳实施例，本专利技术并不限于列举上述实施例，应 当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替换、明 显变形形式，均落在本说明书的实质范围之内，理应受到本专利技术的保护。【主权项】1. 一种基于CT扫描的煤岩逆向化建模力学参数确定方法，其特征在于，包括如下步 骤： (1) 运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤岩试 件体积V，结合基本的计算公V=a3、，计算得出煤岩试件的密度P。 (2) 采用超声波探伤仪，对煤岩试件进行超声波实验，分别以横波、纵波穿过煤岩试件。 记录下横波及纵波穿过煤岩试件的时间，时间分别为：ts、tp，并计算出横波及纵波波速为1得出煤岩试件弹性模量E。 (3) 利用CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件利用二维灰 度图像构建三维重建模型，对三维重建模型内本文档来自技高网...

【技术保护点】
一种基于CT扫描的煤岩逆向化建模力学参数确定方法，其特征在于，包括如下步骤：(1)运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤岩试件体积V，结合基本的计算公V＝a3、计算得出煤岩试件的密度ρ。(2)采用超声波探伤仪，对煤岩试件进行超声波实验，分别以横波、纵波穿过煤岩试件。记录下横波及纵波穿过煤岩试件的时间，时间分别为：ts、tp，并计算出横波及纵波波速为vs=ats,vp=atp,]]>利用公式E=ρvs2(3vp2-4vs2)vp2-vs2]]>得出煤岩试件弹性模量E。(3)利用CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件利用二维灰度图像构建三维重建模型，对三维重建模型内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi。(4)结合煤岩试件质量m、不同组份材料体积Vi及煤岩试件弹性模量E，首先运用公式Mi＝mi1+mi2+……+mj＝ρi1Vi1+ρi2Vi2+……+ρijV计算出不同组份材料的密度ρij，得出不同组份材料密度ρij，结合公式求出不同组份材料的弹性模量Ei。

【技术特征摘要】

【专利技术属性】
技术研发人员：赵毅鑫，张通，许多，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京;11