一种基于ct扫描的煤岩逆向化建模力学参数确定方法
【技术领域】
[0001] 本发明属于数学计算方法领域，涉及煤岩逆向化建模力学参数确定方法，尤其涉 及一种基于CT扫描的煤岩逆向化建模力学参数确定方法。
【背景技术】
[0002] 煤岩在细观尺度上是非均匀材料，其物理力学性质在空间上是非连续性。长期以 来的理论研宄是建立在均质性的假设之上研宄岩石的力学性质及变化过程，而非均质性对 于煤岩的强度、应力分布及破坏形式起着重要的作用。在煤岩非均匀性理论分析上，常采用 统计分布规律来描述。在对煤岩破坏形式的数值计算中，大部分学者是在煤岩非均质性理 论分析的基础之上来探讨非均质性对煤岩破坏形式的影响，例如Weibull分布以及以此为 基础探讨煤岩的非均质性对其破坏形式影响的数值计算。介于煤岩试件中所含的矿物成份 比较复杂，为了真实的展现煤岩试件材料的非均质性，一部分学者通过工业CT及逆向化三 维重建软件来实现煤岩试件非均质性的真实展现。但是，目前基于逆向化来探讨非均质性 对煤岩破裂形式影响的数值计算，并不能够准确获取各组份材料的参数，数值计算结果的 准确性有待考证。基于上述情况，迫切需要一种简单可靠，基于CT扫描的煤岩逆向化建模 力学参数确定方法，有效可提高煤岩内部裂表征精度的基于CT扫描的受载煤岩损伤本构 方程构建方法，达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。
【发明内容】
[0003] 本发明运用游标卡尺及电子称测出煤岩试件的基本尺寸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扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。为了 实现上述目的，本发明采用如下技术方案：
[0004] 一种基于CT扫描的煤岩逆向化建模力学参数确定方法，其特征在于，包括如下步 骤：
[0005] (1)运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤 岩试件体积V，结合基本的计算公V=a3、P= ^7，计算得出煤岩试件的密度P。
[0006] (2)采用超声波探伤仪，对煤岩试件进行超声波实验，分别以横波、纵波穿过煤岩 试件。记录下横波及纵波穿过煤岩试件的时间，时间分别为：ts、tp，并计算出横波及纵波波 速戈
W用公式
得出煤岩试件弹性模量E。
[0007] (3)利用CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件利用二 维灰度图像构建三维重建模型，对三维重建模型内部不同矿物成份的材料进行分组以及统 计不同组份材料体积Vi。
[0008] (4)结合煤岩试件质量m、不同组份材料体积Vi及煤岩试件弹性模量E，首先运用 公式Mi=mn+mi2+......+mj=p nVn+p i2Vi2+......+P 计算出不同组份材料的密度P "， 得出不同组份材料密度Pw结合公3
:出不同组份材料 的弹性模量Ep
[0009] 与现有相关方法相比，本发明具有如下优点：
[0010] 本发明提供了一种新的基于CT扫描的煤岩逆向化建模力学参数确定方法，测量 出煤岩试件基本尺寸、煤岩试件质量并求出煤岩试件体积，然后对煤岩试件进行超声波实 验，利用相关公式推导出煤岩试件不同组份材料的参数，为基于CT扫描的逆向化三维重构 煤岩试件数值计算提供精准参数。
【附图说明】
[0011] 图1为本发明的煤岩试件体积求解过程图；
[0012] 图2为本发明的煤岩试件质量求解过程图；
[0013] 图3为本发明的煤岩试件弹性模量求解过程图；
[0014] 图4为本发明的不同组份材料体积求解过程图。
[0015] 图中：1_煤岩试件；2-CT扫描系统；3-二维灰度图像；4-三维重建软件；5-三维 重建模型；6-不同组份材料体积Vi;7_游标卡尺；8-煤岩试件体积V;9_电子称；10-煤岩 试件质量m;11_超声波探伤仪；12-煤岩试件弹性模量E。
【具体实施方式】
[0016] 结合图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
[0017] 以某矿立方体非标准煤岩试件为例，阐述本发明具体步骤如下：
[0018] (1)运用游标卡尺7及电子称9测出煤岩试件1的基本尺寸a、煤岩试件质量mlO 并求出煤岩试件体积V8,结合基本的计算公式V=a3、P= ^，计算得出煤岩试件1的密度 P〇
[0019] (2)采用超声波探伤仪11，对煤岩试件1进行超声波实验，分别以横波、纵波穿过 煤岩试件1。记录下横波及纵波穿过煤岩试件1的时间，时间分别为：ts、tp，并计算出横波 及纵波波速为
I岩试件弹性模量E12。
[0020] (3)利用CT扫描系统2扫描煤岩试件1得出二维灰度图像3,运用三维重建软件 4利用二维灰度图像3构建三维重建模型5,对三维重建模型5内部不同矿物成份的材料进 行分组以及统计不同组份材料体积Vi6。
[0021] (4)结合煤岩试件质量mlO、不同组份材料体积Vi6及煤岩试件弹性模量E12,首 先运用公式Mi=mn+mi2+......+mj=pnVn+pi2Vi2+......+PUV计算出不同组份材料的密 度Pw得出不同组份材料密度Pw结合公3:出不同组 份材料的弹性模量Ep
[0022] 当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应 当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替换、明 显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。
【主权项】
1. 一种基于CT扫描的煤岩逆向化建模力学参数确定方法，其特征在于，包括如下步 骤： (1) 运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤岩试 件体积V，结合基本的计算公V=a3、，计算得出煤岩试件的密度P。 (2) 采用超声波探伤仪，对煤岩试件进行超声波实验，分别以横波、纵波穿过煤岩试件。 记录下横波及纵波穿过煤岩试件的时间，时间分别为：ts、tp，并计算出横波及纵波波速为1得出煤岩试件弹性模量E。 (3) 利用CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件利用二维灰 度图像构建三维重建模型，对三维重建模型内部不同矿物成份的材料进行分组以及统计不 同组份材料体积Vi。 (4) 结合煤岩试件质量m、不同组份材料体积Vi及煤岩试件弹性模量E，首先运用公式 Mi=mil+mi2+......+mj= p nVn+p i2Vi2+......+ P"V计算出不同组份材料的密度P得出 不同组份材料密度Pm结合公式求出不同组份材料的弹 性模量Ep
【专利摘要】本发明公开了一种基于CT扫描的煤岩逆向化建模力学参数确定方法，运用游标卡尺及电子称测出煤岩试件的基本尺寸a及煤岩试件质量m并求出煤岩试件体积V，利用公式V＝a3、计算出煤岩试件密度，采用超声波探伤仪探测出穿过煤岩试件横波和纵波用时ts、tp，并利用公式求出煤岩试件弹性模量E，CT扫描系统扫描煤岩试件得出二维灰度图像，运用三维重建软件对煤岩试件内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi，得出多个煤岩试件基本参数后，运用公式Mi＝mi1+mi2+……+mj＝ρi1Vi1+ρi2Vi2+……+ρijV得出不同组份材料密度ρij，最终运用公式得出不同组份材料的弹性模量Ei，达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。
【IPC分类】G01N3/00, G01N29/07, G01N23/04
【公开号】CN104990788
【申请号】CN201510374180
【发明人】赵毅鑫, 张通, 许多
【申请人】中国矿业大学(北京)
【公开日】2015年10月21日
【申请日】2015年7月1日