CN104990788A - 一种基于ct扫描的煤岩逆向化建模力学参数确定方法 - Google Patents
一种基于ct扫描的煤岩逆向化建模力学参数确定方法 Download PDFInfo
- Publication number
- CN104990788A CN104990788A CN201510374180.4A CN201510374180A CN104990788A CN 104990788 A CN104990788 A CN 104990788A CN 201510374180 A CN201510374180 A CN 201510374180A CN 104990788 A CN104990788 A CN 104990788A
- Authority
- CN
- China
- Prior art keywords
- coal rock
- rock specimen
- coal
- component materials
- different component
- 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.)
- Pending
Links
- 239000003245 coal Substances 0.000 title claims abstract description 104
- 239000011435 rock Substances 0.000 title claims abstract description 87
- 238000002591 computed tomography Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 title claims abstract description 16
- 239000000463 material Substances 0.000 claims abstract description 37
- 239000000203 mixture Substances 0.000 claims abstract description 8
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 7
- 239000011707 mineral Substances 0.000 claims abstract description 7
- 230000007547 defect Effects 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009172 bursting Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于CT扫描的煤岩逆向化建模力学参数确定方法,运用游标卡尺及电子称测出煤岩试件的基本尺寸a及煤岩试件质量m并求出煤岩试件体积V,利用公式V=a3、计算出煤岩试件密度,采用超声波探伤仪探测出穿过煤岩试件横波和纵波用时ts、tp,并利用公式求出煤岩试件弹性模量E,CT扫描系统扫描煤岩试件得出二维灰度图像,运用三维重建软件对煤岩试件内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi,得出多个煤岩试件基本参数后,运用公式Mi=mi1+mi2+……+mj=ρi1Vi1+ρi2Vi2+……+ρijV得出不同组份材料密度ρij,最终运用公式得出不同组份材料的弹性模量Ei,达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。
Description
技术领域
本发明属于数学计算方法领域,涉及煤岩逆向化建模力学参数确定方法,尤其涉及一种基于CT扫描的煤岩逆向化建模力学参数确定方法。
背景技术
煤岩在细观尺度上是非均匀材料,其物理力学性质在空间上是非连续性。长期以来的理论研究是建立在均质性的假设之上研究岩石的力学性质及变化过程,而非均质性对于煤岩的强度、应力分布及破坏形式起着重要的作用。在煤岩非均匀性理论分析上,常采用统计分布规律来描述。在对煤岩破坏形式的数值计算中,大部分学者是在煤岩非均质性理论分析的基础之上来探讨非均质性对煤岩破坏形式的影响,例如Weibull分布以及以此为基础探讨煤岩的非均质性对其破坏形式影响的数值计算。介于煤岩试件中所含的矿物成份比较复杂,为了真实的展现煤岩试件材料的非均质性,一部分学者通过工业CT及逆向化三维重建软件来实现煤岩试件非均质性的真实展现。但是,目前基于逆向化来探讨非均质性对煤岩破裂形式影响的数值计算,并不能够准确获取各组份材料的参数,数值计算结果的准确性有待考证。基于上述情况,迫切需要一种简单可靠,基于CT扫描的煤岩逆向化建模力学参数确定方法,有效可提高煤岩内部裂表征精度的基于CT扫描的受载煤岩损伤本构方程构建方法,达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。
发明内容
本发明运用游标卡尺及电子称测出煤岩试件的基本尺寸a及煤岩试件质量m并求出煤岩试件体积V,利用公式V=a3、计算出煤岩试件密度,采用超声波探伤仪探测出穿过煤岩试件横波和纵波用时ts、tp,并利用公式求出煤岩试件弹性模量E,CT扫描系统扫描煤岩试件得出二维灰度图像,运用三维重建软件对煤岩试件内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi,得出多个煤岩试件基本参数后,运用公式Mi=mi1+mi2+……+mj=ρi1Vi1+ρi2Vi2+……+ρijV得出不同组份材料密度ρij,最终运用公式得出不同组份材料的弹性模量Ei,达到为基于CT扫描的逆向化三维重构煤岩试件提供精确计算参数的目的。为了实现上述目的,本发明采用如下技术方案:
一种基于CT扫描的煤岩逆向化建模力学参数确定方法,其特征在于,包括如下步骤:
(1)运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤岩试件体积V,结合基本的计算公V=a3、计算得出煤岩试件的密度ρ。
(2)采用超声波探伤仪,对煤岩试件进行超声波实验,分别以横波、纵波穿过煤岩试件。记录下横波及纵波穿过煤岩试件的时间,时间分别为:ts、tp,并计算出横波及纵波波速为利用公式得出煤岩试件弹性模量E。
(3)利用CT扫描系统扫描煤岩试件得出二维灰度图像,运用三维重建软件利用二维灰度图像构建三维重建模型,对三维重建模型内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi。
(4)结合煤岩试件质量m、不同组份材料体积Vi及煤岩试件弹性模量E,首先运用公式Mi=mi1+mi2+……+mj=ρi1Vi1+ρi2Vi2+……+ρijV计算出不同组份材料的密度ρij,得出不同组份材料密度ρij,结合公式求出不同组份材料的弹性模量Ei。
与现有相关方法相比,本发明具有如下优点:
本发明提供了一种新的基于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、煤岩试件质量m10、超声波探伤仪11、煤岩试件弹性模量E12,运用游标卡尺7及电子称9测出煤岩试件1的基本尺寸及煤岩试件质量m10并求出煤岩试件体积V8,利用公式V=a3、算出煤岩试件密度,采用超声波探伤仪11探测出穿过煤岩试件1横波和纵波用时ts、tp,并利用公式计算出煤岩试件弹性模量E12,CT扫描系统2扫描煤岩试件1得出二维灰度图像3,运用三维重建软件4对煤岩试件1内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi6,得出多个煤岩试件1基本参数后,运用公式Mi=mi1+mi2+……+mj=ρi1Vi1+ρi2Vi2+……+ρijV得出不同组份材料密度ρij,最终运用公式得出不同组份材料的弹性模量Ei。
以某矿立方体非标准煤岩试件为例,阐述本发明具体步骤如下:
(1)运用游标卡尺7及电子称9测出煤岩试件1的基本尺寸a、煤岩试件质量m10并求出煤岩试件体积V8,结合基本的计算公式V=a3、计算得出煤岩试件1的密度ρ。
(2)采用超声波探伤仪11,对煤岩试件1进行超声波实验,分别以横波、纵波穿过煤岩试件1。记录下横波及纵波穿过煤岩试件1的时间,时间分别为:ts、tp,并计算出横波及纵波波速为利用公式煤岩试件弹性模量E12。
(3)利用CT扫描系统2扫描煤岩试件1得出二维灰度图像3,运用三维重建软件4利用二维灰度图像3构建三维重建模型5,对三维重建模型5内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi6。
(4)结合煤岩试件质量m10、不同组份材料体积Vi6及煤岩试件弹性模量E12,首先运用公式Mi=mi1+mi2+……+mj=ρi1Vi1+ρi2Vi2+……+ρijV计算出不同组份材料的密度ρij,得出不同组份材料密度ρij,结合公式求出不同组份材料的弹性模量Ei。
当然,以上说明仅仅为本发明的较佳实施例,本发明并不限于列举上述实施例,应当说明的是,任何熟悉本领域的技术人员在本说明书的教导下,所做出的所有等同替换、明显变形形式,均落在本说明书的实质范围之内,理应受到本发明的保护。
Claims (1)
1.一种基于CT扫描的煤岩逆向化建模力学参数确定方法,其特征在于,包括如下步骤:
(1)运用游标卡尺及电子称测出煤岩试件的基本尺寸a、煤岩试件质量m并求出煤岩试件体积V,结合基本的计算公V=a3、计算得出煤岩试件的密度ρ。
(2)采用超声波探伤仪,对煤岩试件进行超声波实验,分别以横波、纵波穿过煤岩试件。记录下横波及纵波穿过煤岩试件的时间,时间分别为:ts、tp,并计算出横波及纵波波速为 利用公式 得出煤岩试件弹性模量E。
(3)利用CT扫描系统扫描煤岩试件得出二维灰度图像,运用三维重建软件利用二维灰度图像构建三维重建模型,对三维重建模型内部不同矿物成份的材料进行分组以及统计不同组份材料体积Vi。
(4)结合煤岩试件质量m、不同组份材料体积Vi及煤岩试件弹性模量E,首先运用公式Mi=mi1+mi2+……+mj=ρi1Vi1+ρi2Vi2+……+ρijV计算出不同组份材料的密度ρij,得出不同组份材料密度ρij,结合公式求出不同组份材料的弹性模量Ei。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510374180.4A CN104990788A (zh) | 2015-07-01 | 2015-07-01 | 一种基于ct扫描的煤岩逆向化建模力学参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510374180.4A CN104990788A (zh) | 2015-07-01 | 2015-07-01 | 一种基于ct扫描的煤岩逆向化建模力学参数确定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104990788A true CN104990788A (zh) | 2015-10-21 |
Family
ID=54302629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510374180.4A Pending CN104990788A (zh) | 2015-07-01 | 2015-07-01 | 一种基于ct扫描的煤岩逆向化建模力学参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104990788A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108663287A (zh) * | 2018-04-25 | 2018-10-16 | 中国地质大学(北京) | 一种利用ct图像精确计算煤岩密度的方法 |
CN110308166A (zh) * | 2019-07-05 | 2019-10-08 | 中国矿业大学 | 利用x射线测定煤岩中水分含量与分布的实验装置及方法 |
CN105973706B (zh) * | 2016-06-07 | 2020-07-03 | 中国矿业大学(北京) | 一种基于工业ct的煤岩体多尺度力学特性分析方法 |
CN112557424A (zh) * | 2020-11-25 | 2021-03-26 | 东南大学 | 一种预测水泥浆体有效弹性模量的方法 |
CN113390906A (zh) * | 2021-06-11 | 2021-09-14 | 内蒙古科技大学 | 一种细观尺度下低渗透性煤岩体致裂增透效果的评价方法 |
CN113588419A (zh) * | 2021-07-27 | 2021-11-02 | 中国科学院光电技术研究所 | 一种准确、快捷测量反应烧结SiC弹性模量的方法 |
CN116183875A (zh) * | 2023-04-21 | 2023-05-30 | 煤炭科学研究总院有限公司 | 岩体波速的确定方法、装置和电子设备 |
-
2015
- 2015-07-01 CN CN201510374180.4A patent/CN104990788A/zh active Pending
Non-Patent Citations (7)
Title |
---|
北京工商大学物理教研室主编: "《大学物理实验》", 31 January 2010, 北京:机械工业出版社 * |
吴建宝等: "《大学物理实验教程》", 30 November 2013, 北京:清华大学出版社 * |
王文婕: ""煤层冲击倾向性对冲击地压的影响机制研究"", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
胡彦智等: ""基于微焦点CT的煤岩材料变形破坏特征数值模拟研究"", 《北京力学会第21届学术年会暨北京振动工程学会第22届学术年会论文集》 * |
赵毅鑫等: ""基于微焦点CT扫描的煤岩材料应变局部化数值模拟研究"", 《塑性力学新进展-2011年全国塑性力学会议论文集》 * |
郑剑: "《高能固体推进剂性能及配方设计专家系统》", 30 June 2014, 北京:国防工业出版社 * |
陈世涛等: "《大学物理实验教程 第2版》", 31 January 2014, 成都:西南交通大学出版社 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973706B (zh) * | 2016-06-07 | 2020-07-03 | 中国矿业大学(北京) | 一种基于工业ct的煤岩体多尺度力学特性分析方法 |
CN108663287A (zh) * | 2018-04-25 | 2018-10-16 | 中国地质大学(北京) | 一种利用ct图像精确计算煤岩密度的方法 |
CN108663287B (zh) * | 2018-04-25 | 2021-02-26 | 中国地质大学(北京) | 一种利用ct图像精确计算煤岩密度的方法 |
CN110308166A (zh) * | 2019-07-05 | 2019-10-08 | 中国矿业大学 | 利用x射线测定煤岩中水分含量与分布的实验装置及方法 |
CN110308166B (zh) * | 2019-07-05 | 2020-06-05 | 中国矿业大学 | 利用x射线测定煤岩中水分含量与分布的实验装置及方法 |
WO2021003996A1 (zh) * | 2019-07-05 | 2021-01-14 | 中国矿业大学 | 利用x射线测定煤岩中水分含量与分布的实验装置及方法 |
CN112557424A (zh) * | 2020-11-25 | 2021-03-26 | 东南大学 | 一种预测水泥浆体有效弹性模量的方法 |
CN113390906A (zh) * | 2021-06-11 | 2021-09-14 | 内蒙古科技大学 | 一种细观尺度下低渗透性煤岩体致裂增透效果的评价方法 |
CN113588419A (zh) * | 2021-07-27 | 2021-11-02 | 中国科学院光电技术研究所 | 一种准确、快捷测量反应烧结SiC弹性模量的方法 |
CN116183875A (zh) * | 2023-04-21 | 2023-05-30 | 煤炭科学研究总院有限公司 | 岩体波速的确定方法、装置和电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104990788A (zh) | 一种基于ct扫描的煤岩逆向化建模力学参数确定方法 | |
US20220136945A1 (en) | Method for evaluating anchor bolt embedment in concrete | |
Li et al. | Comparison of two methods in acoustic emission source location using four sensors without measuring sonic speed | |
CN203981165U (zh) | 一种新型河工动床模型地形测量装置 | |
CN104833599A (zh) | 测定混凝土损伤演变规律的试验装置和定量描述方法 | |
CN105158342A (zh) | 一种超声水浸无损评价残余应力的方法 | |
CN105158796B (zh) | 确定toc含量的方法和装置 | |
CN104976970A (zh) | 一种超声相控阵楔块自检测系统及其方法 | |
CN103091342B (zh) | 一种对岩芯样品进行ct扫描分析处理的方法 | |
CN104634876A (zh) | 一种超声扫描显微镜检测金属材料内部夹杂物的方法 | |
CN110346453A (zh) | 一种混凝土结构内缺陷小排列反射回波快速检测方法 | |
CN110795793A (zh) | 一种隧道围岩快速分级设备系统及其操作方法 | |
Resende et al. | Rock joint topography: three-dimensional scanning and numerical analysis | |
CN206235600U (zh) | 一种用于校准超声波检测仪系统延迟时间的消声水槽装置 | |
CN104165795A (zh) | 一种古建筑木梁的剩余抗弯承载力测定方法 | |
CN105044215A (zh) | 一种非破坏性的材料声速现场测量方法 | |
Shah et al. | Photogrammetry and Monte Carlo Simulation based statistical characterization of rock mass discontinuity parameters | |
Draudvilienė et al. | Validation of dispersion curve reconstruction techniques for the A0 and S0 modes of Lamb waves | |
CN104573357A (zh) | 一种预应力通道注浆密实度检测方法 | |
Borosnyói | NDT assessment of existing concrete structures: spatial analysis of rebound hammer results recorded in-situ | |
CN102735374B (zh) | 一种确定具有轴向静载荷的弹性杆中入射波大小的方法 | |
Bragov et al. | A combined approach to dynamic testing of structural materials | |
CN110231217B (zh) | 一种坦克扭力轴力学性能退化的共线混频超声检测方法 | |
CN105891339A (zh) | 岩土工程介质的信号处理方法、装置及系统 | |
Tan et al. | Stochastic three-dimensional joint geometry model and the properties of REV for a jointed rock mass |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151021 |
|
RJ01 | Rejection of invention patent application after publication |