CN113155873A - 一种基于x线ct设备技术测算木材等效扩散系数的方法 - Google Patents
一种基于x线ct设备技术测算木材等效扩散系数的方法 Download PDFInfo
- Publication number
- CN113155873A CN113155873A CN202110391745.5A CN202110391745A CN113155873A CN 113155873 A CN113155873 A CN 113155873A CN 202110391745 A CN202110391745 A CN 202110391745A CN 113155873 A CN113155873 A CN 113155873A
- Authority
- CN
- China
- Prior art keywords
- wood
- dimensional image
- ray
- measuring
- calculating
- 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
- 239000002023 wood Substances 0.000 title claims abstract description 42
- 238000009792 diffusion process Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000005516 engineering process Methods 0.000 title claims abstract description 18
- 239000011148 porous material Substances 0.000 claims abstract description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 23
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000010586 diagram Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 238000009877 rendering Methods 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 238000012800 visualization Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 5
- 238000004445 quantitative analysis Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 230000008034 disappearance Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 238000004458 analytical method Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 abstract description 3
- 238000002591 computed tomography Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (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)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于X线CT设备技术测算木材等效扩散系数的方法，通过采用X‑CT技术测量木材真实的孔隙结构进行测算水分流，从而进行扩散等效模拟，不损坏试件，可呈现3维结构体，并且在这基础上，可以将描述水分扩散现象的控制方程引入几何体中，直接构建水分扩散模型，通过研究孔隙结果对等效扩散系数的影响规律，极大程度的保证模拟计算精度。
Description
技术领域
本发明涉及木材科学领域，本发明提供一种基于X线CT设备技术测算木材等效扩散系数的方法。
背景技术
扩散现象，是指在浓度梯度作用下，液体或气体在材料内部迁移的过程，在木材科学领域，扩散性质是极其重要的表征水分在木材内部迁移能力的物性参数之一；其深入研究可完善木材与水分关系的基础理论，对木材干燥、木材改性等工的艺优化具有理论指导意义。
关于木材水分扩散性质的研究，主要集中在木材扩散系数求解及影响因素分析上，即在有限的检测条件下，以实验数据为基础，借助扩散理论构建水分迁移模型，进而逆求木材水分扩散系数；在此基础上，改变木材温度与含水率状态，建立扩散系数与外部因素的联系，该方法被称为逆分析法，目前多数研究者采用此手段计算木材水分扩散系数。
但是，上述方法的缺点是，不能通过研究孔隙结果对等效扩散系数的影响规律；同时，实验过程繁琐，检测精度极大地受限于实验条件设备精度。
发明内容
本发明的目的在于提供一种基于X线CT设备技术测算木材等效扩散系数的方法，已解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：.一种基于X线CT设备技术测算木材等效扩散系数的方法，包括有以下步骤：
S1扫描：将木材样件放置进入X线CT设备进行扫描得出二维图像扫描数据体；
S2三维图像重建：将二维图像扫描数据体经Avizo软件中进行可视化处理得出三维图像，其中，可视化处理由以下工序组成：
a)立体渲染工序，将二维图像扫描进行三维体渲染得出三维图像；
b)图像裁切设定工序，通过在三维图像上选取部分区域进行详细分析，从而达到减少图像处理计算量；
c)噪声过滤，对选取后的图像进行优化，过滤由于环境与设备带来的影响；
d)阈值分割，选取合适括值对木材孔隙与基质进行分割；
e)标记，通过标记确定每个3维孔隙结构；
f)量化分析，定量计算诸如孔隙率、孔体积及孔径分布等参数；
S3等效扩散计算：将S2得出的三维图像结构体部分模拟图进行设定孔隙与木材基层导热系数，在三维图像结构体部分模拟图上选择计算区域与选择计算方向。
作为本发明的进一步方案，所述木材样件呈圆柱形，且所述三维图像结构体模拟图呈圆柱形。
作为本发明的进一步方案，所述S2的工序b中所述图像裁切选取部分区域裁切尺寸为200*200*200pixel。
作为本发明的进一步方案，所述S3包括以下步骤：
S3.1通过选中conducting materials的“0”代表孔与“1”代表木材基质；根据需要选中水分流方向，其中含有X,Y,Z，三个方向；
S3.2边界条件选择温度的流入与流出：且所述流入温度设定298K，所述流出温度设定273K；
S3.3设置基质与孔隙内空气的导热系数：基质导热系数值设定为0.4396，孔隙的导热系数值0.0236W m-1 K-1，其它选项为默认值；
S3.4在浓度分布模拟结果中，进行可视化选择。首先在display中，分别选择层(slice)，等高面(isosurface)与等高线(isocontour)三个选项；
S3.5计算水分流流动方向，即在浓度流数据上，选择compute——gradient计算，直接单击应用，生成浓度的梯度分布；在此数据体上，单击display-illuminatedstreamlines，在此选项中选择，流线数num line为300，线长length为120，透明度opacity为0.5，消失因子fade factor为0.9，步尺寸step size为500000，分布distribute选择proportional，其它选项为默认值。
与现有技术相比，本发明公开了一种基于X线CT设备技术测算木材等效扩散系数的方法：通过采用X线CT设备技术测量木材真实的孔隙结构进行测算水分流，从而进行扩散等效模拟，不损坏试件，可呈现3维结构体，并且在这基础上，可以将描述水分扩散现象的控制方程引入几何体中，直接构建水分扩散模型，通过研究孔隙结果对等效扩散系数的影响规律，极大程度的保证模拟计算精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明的圆柱形三维图像示意图。
图2是本发明的部分区域裁剪示意图。
图3是本发明的水分流模拟示意图。
具体实施方式
下面结合附图对本发明作作为本发明的进一步方案说明：
参考图1至图2，一种基于X线CT设备技术测算木材等效扩散系数的方法，包括有以下步骤：
S1扫描：将木材样件放置进入X线CT设备进行扫描得出二维图像扫描数据体；
S2三维图像重建：将二维图像扫描数据体经Avizo软件中进行可视化处理得出三维图像，其中，可视化处理由以下工序组成：
a)立体渲染工序，将二维图像扫描进行三维体渲染得出三维图像；
b)图像裁切设定工序，通过在三维图像上选取部分区域进行详细分析，从而达到减少图像处理计算量；
c)噪声过滤，对选取后的图像进行优化，过滤由于环境与设备带来的影响；
d)阈值分割，选取合适括值对木材孔隙与基质进行分割；
e)标记，通过标记确定每个3维孔隙结构；
f)量化分析，定量计算诸如孔隙率、孔体积及孔径分布等参数；
S3等效扩散计算：将S2得出的三维图像结构体部分模拟图进行设定孔隙与木材基层导热系数，在三维图像结构体部分模拟图上选择计算区域与选择计算方向。
作为本发明的进一步方案，所述木材样件呈圆柱形，且所述三维图像结构体模拟图呈圆柱形。
作为本发明的进一步方案，所述S2的工序b中所述图像裁切选取部分区域裁切尺寸为200*200*200pixel。
作为本发明的进一步方案，所述S3包括以下步骤：
S3.1通过选中conducting materials的“0”代表孔与“1”代表木材基质；根据需要选中水分流方向，其中含有X,Y,Z，三个方向；
S3.2边界条件选择温度的流入与流出：且所述流入温度设定298K，所述流出温度设定273K；
S3.3设置基质与孔隙内空气的导热系数：基质导热系数值设定为0.4396，孔隙的导热系数值0.0236W m-1 K-1，其它选项为默认值；
S3.4在浓度分布模拟结果中，进行可视化选择。首先在display中，分别选择层(slice)，等高面(isosurface)与等高线(isocontour)三个选项；
S3.5计算水分流流动方向，即在浓度流数据上，选择compute——gradient计算，直接单击应用，生成浓度的梯度分布；在此数据体上，单击display-illuminatedstreamlines，在此选项中选择，流线数num line为300，线长length为120，透明度opacity为0.5，消失因子fade factor为0.9，步尺寸step size为500000，分布distribute选择proportional，其它选项均为默认值。
实施例一
其中，在slice中，选择合适的方向角度，即orientastion，xy-zx-yz；在colormap中将显示颜色设定为graylevel，根据在窗口中图像的视觉效果，进行具体数值设定，范围在298-273之间；其它选项为默认值。
其中，在isocontour中，选择合适的方向角度，即orientastion，xy-zx-yz；将阈值范围设定为298-273K之间，显示等高线条数为30；图像质量设定为高，分辨率设定300；线宽设定2；其它选项为默认值，单击apply。
在isosurface中，选择合适的方向角度，即orientastion，xy-zx-yz；点击auto-refresh自动更新；阈值设定为min-max覆盖整个区域；在colormap中将显示颜色设定为graylevel；其它选项为默认值，单击apply。
等效导热系数的各向异性表征thermal conductivity tensor simulation，这里我们在此模块中重复上述步骤1，但此时，必须将3个方向的水分流都选中，单击运行，这里直接计算出如下矩阵数据的形式，代表材料内部水分扩散性质具有矢量性。
再重复上述2与3步骤，将再次出现模拟结果，但此模拟结果是3个方向的矢量结合，不具有单一方向性。
最终，所获数据的形式由3*3的特征值与特征向量表达。
本发明公开了一种基于X线CT设备技术测算木材等效扩散系数的方法：通过采用X线CT设备技术测量木材真实的孔隙结构进行测算水分流，从而进行扩散等效模拟，不损坏试件，可呈现3维结构体，并且在这基础上，可以将描述水分扩散现象的控制方程引入几何体中，直接构建水分扩散模型，通过研究孔隙结果对等效扩散系数的影响规律，极大程度的保证模拟计算精度。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变形，本发明的范围由所附权利要求及其等同物限定。
Claims (5)
1.一种基于X线CT设备技术测算木材等效扩散系数的方法，其特征在于：包括有以下步骤：
S1扫描：将木材样件放置进入X线CT设备进行扫描得出二维图像扫描数据体；S2三维图像重建：将二维图像扫描数据体经Avizo软件中进行可视化处理得出三维图像，其中，可视化处理由以下工序组成：
a)立体渲染工序，将二维图像扫描进行三维体渲染得出三维图像；
b)图像裁切设定工序，通过在三维图像上选取部分区域进行详细分析，从而达到减少图像处理计算量；
c)噪声过滤，对选取后的图像进行优化，过滤由于环境与设备带来的影响；
d)阈值分割，选取合适括值对木材孔隙与基质进行分割；
e)标记，通过标记确定每个3维孔隙结构；
f)量化分析，定量计算诸如孔隙率、孔体积及孔径分布等参数；
S3等效扩散计算：将S2得出的三维图像结构体部分模拟图进行设定孔隙与木材基层导热系数，在三维图像结构体部分模拟图上选择计算区域与选择计算方向。
2.根据权利要求1所述一种基于X线CT设备技术测算木材等效扩散系数的方法，其特征在于：所述木材样件呈圆柱形，且所述三维图像结构体模拟图呈圆柱形。
3.根据权利要求1所述一种基于X线CT设备技术测算木材等效扩散系数的方法，其特征在于：所述S2的工序f中所述量化分析选取孔径分析。
4.根据权利要求1所述一种基于X线CT设备技术测算木材等效扩散系数的方法，其特征在于：所述S2的工序b中所述图像裁切选取部分区域裁切尺寸为200*200*200pixel。
5.根据权利要求1所述一种基于X线CT设备技术测算木材等效扩散系数的方法，其特征在于，所述S3包括以下步骤：
S3.1通过选中conducting materials的“0”代表孔与“1”代表木材基质；根据需要选中水分流方向，其中含有X,Y,Z，三个方向；
S3.2边界条件选择温度的流入与流出：且所述流入温度设定298K，所述流出温度设定273K；
S3.3设置基质与孔隙内空气的导热系数：基质导热系数值设定为0.4396，孔隙的导热系数值0.0236W m-1K-1，其它选项为默认值；
S3.4在浓度分布模拟结果中，进行可视化选择，首先在display中，分别选择层(slice)，等高面(isosurface)与等高线(isocontour)三个选项；
S3.5计算水分流流动方向，即在浓度流数据上，选择compute——gradient计算，直接单击应用，生成浓度的梯度分布；在此数据体上，单击display-illuminated streamlines，在此选项中选择，流线数numline为300，线长length为120，透明度opacity为0.5，消失因子fade factor为0.9，步尺寸step size为500000，分布distribute选择proportional，其它选项均为默认值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110391745.5A CN113155873A (zh) | 2021-04-13 | 2021-04-13 | 一种基于x线ct设备技术测算木材等效扩散系数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110391745.5A CN113155873A (zh) | 2021-04-13 | 2021-04-13 | 一种基于x线ct设备技术测算木材等效扩散系数的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113155873A true CN113155873A (zh) | 2021-07-23 |
Family
ID=76890034
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110391745.5A Pending CN113155873A (zh) | 2021-04-13 | 2021-04-13 | 一种基于x线ct设备技术测算木材等效扩散系数的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113155873A (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140003573A1 (en) * | 2012-06-28 | 2014-01-02 | Hitachi, Ltd. | X-Ray CT System for Measuring Three Dimensional Shapes and Measuring Method of Three Dimensional Shapes by X-Ray CT System |
| CN103884628A (zh) * | 2014-03-14 | 2014-06-25 | 大连理工大学 | 一种应用ct测量多孔介质中二氧化碳扩散系数的装置与方法 |
| CN108535162A (zh) * | 2018-03-12 | 2018-09-14 | 广东省测试分析研究所（中国广州分析测试中心） | 一种基于x射线显微ct表征木材孔隙结构的方法 |
| CN109883928A (zh) * | 2019-03-29 | 2019-06-14 | 上海市城市建设设计研究总院（集团）有限公司 | 基于ct扫描和三维重构的水泥混凝土渗透能力模拟方法 |
-
2021
- 2021-04-13 CN CN202110391745.5A patent/CN113155873A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140003573A1 (en) * | 2012-06-28 | 2014-01-02 | Hitachi, Ltd. | X-Ray CT System for Measuring Three Dimensional Shapes and Measuring Method of Three Dimensional Shapes by X-Ray CT System |
| CN103884628A (zh) * | 2014-03-14 | 2014-06-25 | 大连理工大学 | 一种应用ct测量多孔介质中二氧化碳扩散系数的装置与方法 |
| CN108535162A (zh) * | 2018-03-12 | 2018-09-14 | 广东省测试分析研究所（中国广州分析测试中心） | 一种基于x射线显微ct表征木材孔隙结构的方法 |
| CN109883928A (zh) * | 2019-03-29 | 2019-06-14 | 上海市城市建设设计研究总院（集团）有限公司 | 基于ct扫描和三维重构的水泥混凝土渗透能力模拟方法 |
Non-Patent Citations (2)
| Title |
|---|
| LIANGYAN GUO 等: "Pore Structure Characterization of Oak via X-ray Computed Tomography", 《BIORESOURCES》 * |
| 何新婷等: "气体在微孔介质中的等效扩散系数预测", 《工程热物理学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AlRatrout et al. | Automatic measurement of contact angle in pore-space images | |
| Ni et al. | Quantitative 3D spatial characterization and flow simulation of coal macropores based on μCT technology | |
| Roose et al. | Challenges in imaging and predictive modeling of rhizosphere processes | |
| CN108763711B (zh) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 | |
| Georgiadis et al. | Pore-scale micro-computed-tomography imaging: Nonwetting-phase cluster-size distribution during drainage and imbibition | |
| Soulaine et al. | The impact of sub-resolution porosity of X-ray microtomography images on the permeability | |
| Ramandi et al. | Digital rock analysis for accurate prediction of fractured media permeability | |
| Montminy et al. | The 3D structure of real polymer foams | |
| Gerke et al. | Improving watershed-based pore-network extraction method using maximum inscribed ball pore-body positioning | |
| Bock et al. | Geometric classification of open-cell metal foams using X-ray micro-computed tomography | |
| Mascini et al. | Event-based contact angle measurements inside porous media using time-resolved micro-computed tomography | |
| CN104751473B (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
| Cı́slerová et al. | CT derived porosity distribution and flow domains | |
| CN106093083B (zh) | 一种采用数字岩心模拟建立岩电关系的方法 | |
| Yu et al. | Response of pore structures to long‐term fertilization by a combination of synchrotron radiation X‐ray microcomputed tomography and a pore network model | |
| Pot et al. | Effect of resolution, reconstruction settings, and segmentation methods on the numerical calculation of saturated soil hydraulic conductivity from 3D computed tomography images | |
| Jing et al. | Reconstruction and seepage simulation of a coal pore-fracture network based on CT technology | |
| Callow et al. | Optimal X-ray micro-CT image based methods for porosity and permeability quantification in heterogeneous sandstones | |
| Wang et al. | Three-dimensional reconstruction and morphologic characteristics of porous metal fiber sintered sheet | |
| De Boever et al. | Comparison between traditional laboratory tests, permeability measurements and CT-based fluid flow modelling for cultural heritage applications | |
| Zhao et al. | Soil pore identification with the adaptive fuzzy C-means method based on computed tomography images | |
| Ferreira et al. | An analysis of three XCT-based methods to determine the intrinsic permeability of soil aggregates | |
| CN114283254A (zh) | 基于核磁共振数据的岩心数字化孔隙网络模型构建方法 | |
| CN113155873A (zh) | 一种基于x线ct设备技术测算木材等效扩散系数的方法 | |
| Kang et al. | Construction of complex digital rock physics based on full convolution network |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210723 |
|
| RJ01 | Rejection of invention patent application after publication |