专利名称：一种多孔材料孔结构分形维数的计算方法
技术领域：
本发明涉及一种多孔材料孔隙结构的分析表征技术领域，尤其是涉及 一种多孔材料孔结构分形维数的计算方法。
背景技术：
多孔材料是一类重要的结构和功能一体化材料，广泛存在于自然界和 工程应用中。由于其内部大量孔隙的存在使得其具有诸多不同于致密材料 的优异性能，但是由于孔隙大小不一、随机分布，孔隙空间结构的杂乱无 序，使得多孔材料的孔结构描述和性能研究变得复杂化。传统孔结构性能 检测和研究方法是将孔隙等效为圆柱孔模型，这种圆柱孔模型与实际复杂 的孔结构相比，过于简化，难以反映多孔介质中孔的真实空间结构。
借助于分形理论来描述多孔材料的孔结构是新近发展起来的一种有 效方法。多孔材料孔结构中的孔隙、孔壁表面、孔通道等均可以视为分形
体。孔结构分形表征中釆用的分形维数主要包括孔隙体积分形维数、孔 壁表面分形维数、孔通道曲线分形维数、孔截面周长分形维数以及孔截面 分形维数等。这些分形维数分别从不同角度来描述孔结构，任何一个维数 都难以完整表述复杂孔结构的特征。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一 种多孔材料孔结构分形维数的计算方法，其能有效解决现有分形维数不能 完整表述被检测多孔材料复杂孔结构特征的实际问题。
为解决上述技术问题，本发明釆用的技术方案是 一种多孔材料孔结构 分形维数的计算方法，其特征在于该方法包括以下步骤步骤一、采用图像釆集设备摄取被检测多孔材料外表面或一个剖面的 二维孔隙结构图像，并将所采集到的二维孔隙结构图像传输至计算机进行 分析处理，获得被检测多孔材料的二维孔隙边界图像，所述二维孔隙边界
图像中各孔隙的二维孔隙边界为由多个像素点为顶点组成的平面多边形； 步骤二、用计算机确定所述平面多边形各顶点的坐标值，其坐标原点
为所述二维孔隙边界图像中的任一位置；
步骤三、用计算机处理计算得出各孔隙的形貌因子Z,其计算机处理
过程包括以下步骤
(a) 根据公式Q=L2/S计算各孔隙的不规则因子Q，其中，L和S分别
为单个孔隙的周长和面积；
(b) 计算各孔隙的锐度因子A:对单个孔隙而言，首先确定该孔隙边 界对应平面多边形的所有顶点的内角角度e，并将所述所有顶点按照内角 角度大小划分到多个角度区间中，对各角度区间内的顶点点位个数进行统 计后，计算出该孔隙的锐度因子A;
(c) 计算各孔隙的形貌因子Z: Z-;rQ+尸A，式中％、 /为系数，对于 单个孔隙而言，首先通过试验分别分析孔隙弯曲情况和孔隙曲率半径对孔 隙性能的影响程度，当孔隙弯曲情况为主要影响因素时，,取值大于/值； 当孔隙曲率半径为主要影响因素时，％取值小于/值；
步骤四、计算各孔隙形貌的分形维数以单个孔隙的形貌因子Z为标 度对象，釆用砂盒法对其进行分形分析得到综合描述被检测多孔材料孔隙 结构特征的孔形貌分形维数，计算分析过程包括以下步骤
(A)以步骤一中所述二维孔隙结构图像的中心点为中心，以所述二 维孔隙结构图像的最小边长度为边长A的正方形一，将s,正方形中包含的 所有孔隙对应的Z值累加，得到一个对应的再以所述中心点为中心 将所述正方形一的边长缩小一倍得到边长为& = *6的正方形二，将&正方
形中包含的所有孔隙对应的Z值累加，得到一个对应的W(e2);以此类推， 直到边长"对应的正方形N，"大于所述二维孔隙结构图像中最小孔隙的直径，孔隙结构分析精度越高，"取值越小；最终得到一组数据A^,);
(B) 在以ln/V,化)为纵坐标，lne,为横坐标的平面直角坐标系中，对一 组数据A^,)对应的一组数据点进行线性回归处理，得到一条直线，直线斜 率的绝对值就是被检测多孔材料孔形貌的分形维数A 。
步骤(b)中所述的计算各孔隙的锐度因子A时，根据加权法，通过 公式A-ga,n,，计算出该孔隙的锐度因子A，其中n,为各角度区间内的顶点
i二0
点位个数，a,为对应角度区间的加权系数。
步骤一所述的二维孔隙结构图像传输至计算机进行分析处理时，所述 计算机的分析处理过程包括以下步骤
(C) 将所述二维孔隙结构图像转化为多灰度级灰度图像；
(D) 对所述多灰度级灰度图像进行去噪和滤波处理；
(E) 将经去噪和滤波处理的所述多灰度级灰度图像转换为二灰度级 图像即黑白图像；
所述黑白图像由像素值为1的白色图像区域和像素值为0的黑色图像 区域组成，其中，所述被检测多孔材料的实体骨架材料对应的图像区域为 白色图像区域或黑色图像区域，所述孔隙空间对应的图像区域相应为黑色 图像区域或白色图像区域；
(F) 用计算机提取所述黑白图像中黑色图像区域和白色图像区域的 交界线即孔隙空间边界并用彩色线条进行标示，所标示出来的孔隙空间边 界为由有限个像素点为顶点组成的平面多边形，确定所述平面多边形所有 顶点的二维坐标值并对其各顶点进行连续编号，即获得被检测多孔材料的 一个或多个相互平行断面的二维孔隙边界。
步骤 一 中所述的图像釆集设备为金相显微镜或扫描电子显微镜。 步骤一中所述的二维孔隙结构图像为分辨率不低于500 x 500的清晰
彩色图像或多灰度级灰度图像。
步骤(D)中所述的滤波处理为最小值滤波和/或中值滤波和/或最大
值滤波和/或均值滤波和/或低通滤波和/或高通滤波处理。步骤(E)中所述的实体骨架材料对应的图像区域为白色图像区域，
所述孔隙空间对应的图像区域相应为黑色图像区域。
步骤三中还需对所述多边形进行简化处理，所述简化处理的方法为减 少所述多边形的边数并相应增加各边边长。
步骤(F)中确定所述平面多边形所有顶点的二维坐标值时，选取所 述黑白图像左下方的顶点为坐标原点。
步骤(F)中对所述平面多边形各顶点进行连续编号时，以二维坐标 绝对值最小的顶点为1 #顶点，然后按照逆时针顺序依次对剩余顶点进行连 续编号。
本发明与现有技术相比具有以下优点，1、所得到的孔隙边界即孔隙
空间边界反映了材料的真实孔隙结构，同时具有操作简单、成本低等优点； 2、通过数学分析方法和计算机数据分析处理技术，再结合孔隙边界数值 化后的大量结果数据，对多孔材料结构进行系统分析，利用分形理论对多 孔材料进行数学分析，进而建立孔隙边界特征与材料应用性能间的关系，
以便于指导多孔材料的制造和应用；3、适用范围广，不仅可以对已知的
材质，例如金属、陶瓷、高分子、岩土等制成多孔材料的表面或剖面孔隙 结构进行分析，也可以对骨骼、木材等具有多孔结构的材料的表面或剖面 孔隙结构进行分析。总之，本发明提出的孔形貌分形维数充分考虑了孔隙 的形状和曲率分布情况，不但能够更加科学全面地描述孔隙的真实结构， 而且物理意义更加明确，能够将分形维数与力学性能、渗透性能、过滤性 能等与多孔材料孔隙形貌有关的物理性能参数更加紧密地联系起来，对多 孔材料孔隙结构表征和性能研究具有重要推动作用。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
图1为本发明所获取的二维孔隙边界多边形结构示意图。附图标记说明
4一简化多边形。
具体实施例方式
本发明所述的一种多孔材料孔结构分形维数的计算方法包括以下步
骤
步骤一、釆用图像釆集设备摄取被检测多孔材料外表面或一个剖面的 二维孔隙结构图像，并将所釆集到的二维孔隙结构图像传输至计算机进行 分析处理，获得被检测多孔材料的二维孔隙边界图像，所述二维孔隙边界 图像中各孔隙的二维孔隙边界为由多个像素点为顶点组成的平面多边形。 其中，所釆集的二维孔隙结构图像中，所述多孔材料的实体骨架材料和孔 隙空间分别对应的图像区域之间具有明显的色彩对比。所述图像采集设备 为金相显微镜或扫描电子显微镜。实践中，实践中，还可以釆用其它光学 或电子的材料结构图像釆集设备获取上述微观孔隙结构图像。所检测的多 孔材料可以是已知的材质，例如金属、陶瓷、高分子、岩土等制成多孔材 料的表面或剖面微观结构图像，也可以是骨骼、木材等具有多孔结构的材 料的表面或剖面微观结构图像。
本步骤中，所述二维孔隙结构图像传输至计算机进行分析处理时，所 述计算机的分析处理过程包括以下步骤
(C) 将所述二维孔隙结构图像转化为多灰度级灰度图像即将RGB、 CMYK等色彩模式转换为灰度的色彩模式。
(D) 对所述多灰度级灰度图像进行去噪和滤波处理所述滤波处理 为最小值滤波和/或中值滤波和/或最大值滤波和/或均值滤波和/或低通 滤波和/或高通滤波处理。实践中，可以根据具体实际情况，相应选择最 小值滤波、中值滤波、最大值滤波、均值滤波、低通滤波和高通滤波中的 一种或几种滤波处理方式的组合方式进行滤波处理。
(E) 将经去噪和滤波处理的所述多灰度级灰度图像转换为二灰度级 图像即黑白图像；所述黑白图像由像素值为1的白色图像区域和像素值为o的黑色图像
区域组成，其中，所述被检测多孔材料的实体骨架材料对应的图像区域为 白色图像区域或黑色图像区域，所述孔隙空间对应的图像区域相应为黑色 图像区域或白色图像区域。
也就是说，本步骤中，将所述孔隙结构图像中的所有像素值分别置为
0 (黑色图像区域)或1 (白色图像区域)，本实施例中，所述实体骨架材
料对应的图像区域为白色图像区域，所述孔隙空间对应的图像区域相应为 黑色图像区域。实际操作过程中，根据不同的多孔材料以及不同的孔隙结 构图像获取方式，黑色图像区域和白黑色图像区域分别与实体骨架材料和 孔隙空间的对应关系不同。
(F)用计算机提取所述黑白图像中黑色图像区域和白色图像区域的 交界线即孔隙空间边界并用彩色线条进行标示，所标示出来的孔隙空间边 界为由有限个像素点为顶点组成的平面多边形，确定所述平面多边形所有 顶点的二维坐标值并对其各顶点进行连续编号，即获得被检测多孔材料的 一个或多个相互平行断面的二维孔隙边界。另外，还需对所述平面多边形 进行简化处理，所述简化处理的方法为减少所述平面多边形的边数并相应 增加各边边长。在确定所述平面多边形所有顶点的二维坐标值时，选取所 述黑白图像左下方的顶点为坐标原点，也就是说，以黑白图像左下方的顶 点为坐标原点，确定所述平面多边形谷顶点的二维坐标值。对所述平面多 边形各顶点进行连续编号时，以二维坐标绝对值最小的顶点为1#顶点，然 后按照逆时针顺序依次对剩余顶点进行连续编号。
结合图1,本步骤中，从转换后的黑白图像中提取孔隙空间边界即孔 隙边界，为便于识别孔隙边界，利用计算机将黑色图像区域和白色图像区 域的交界线识别出来后，并用红色或蓝色等彩色线条在所述黑白图像中标 示出来，所标示出来的孔隙空间边界即孔隙边界是由有限个像素点为顶点 组成的多边形。为降低后续计算处理的工作量，可对所标示出来的多边形 进行简化处理，其简化处理的方法是减少多边形的边数并增加多边形的边长，优选的减少多边形的边数并增加多边形边长的实施方法是每隔若干个 像素点将首尾两个像素点连接起来构成简化的多边形，本实施例中，每隔
5个像素点将首尾两个像素点连接起来构成简化多边形4。
步骤二、用计算机确定所述平面多边形各顶点的坐标值，其坐标原点 为所述二维孔隙边界图像中的任一位置。本步骤中，为方便对孔隙边界的 多边形进行数学分析处理，先对孔隙边界多边形进行数值化处理，再利用 计算机确定所述多边形所有顶点的二维坐标值，并对所述多边形的所有顶 点进行编号。本实施例中，确定所述平面多边形所有顶点的二维坐标值时， 选取所述二维孔隙结构图像即所述黑白图像左下方的顶点为坐标原点，而 对所述平面多边形各顶点进行连续编号时，以二维坐标绝对值最小的顶点 为1#顶点，然后按照逆时针顺序依次对剩余的12个顶点进行连续编号。
步骤三、用计算机处理计算得出各孔隙的形貌因子Z，其计算机处理 过程包括以下步骤
(a) 根据公式Q-L7S计算各孔隙的不规则因子Q，其中，L和S分别 为单个孔隙的周长和面积。
(b) 计算各孔隙的锐度因子A:对单个孔隙而言，首先确定该孔隙边
界对应平面多边形的所有顶点的内角角度e,并将所述所有顶点按照内角
角度大小划分到多个角度区间中，对各角度区间内的顶点点位个数进行统 计后，计算出该孔隙的锐度因子A。本步骤中所述的计算各孔隙的锐度因 子A时，根据加权法，通过公式A^ga,n,，计算出该孔隙的锐度因子A,
i=0
其中A为各角度区间内的顶点点位个数，a,为对应角度区间的加权系数。
本实施例中，角度区间的默认值为l(T (实际操作中根据计算精度的 需要确定角度区间，角度区间越�。仍礁撸扑懔吭酱�)，也就是说， 各个角度区间的起至角度根据需要进行调整，典型的区间划分方法为以0 。为起点每隔10°划分一个角度区间。加权系数a,根据多孔材料性能的关注 重点进行确定，例如，对力学和腐蚀性能，孔隙边界各处的角度大小影响 较显著，角度小的加权系数取值较大。系数^的取值应根据已有孔隙的曲率对性能影响关系的认识，例如，研究材料的力学性能时，a,取对应角度 区间初始值(弧度)的平方(例如，对30'-40"的角度区间，取30"对应的 弧度值的平方(tt/6)2)。
(c)计算各孔隙的形貌因子Z: Z=^Q+rA,式中％、 /为系数，对于 单个孔隙而言，首先通过试验分别分析孔隙弯曲情况和孔隙曲率半径对孔 隙性能的影响程度，当孔隙弯曲情况为主要影响因素时，x取值大于/值； 当孔隙曲率半径为主要影响因素时，％取值小于/值。优选地，釆用实验 的方法确定,、/值。以研究多孔材料的压缩强度为例，针对一组实验样品进 行试验时，根据上述计算方法获取每种样品的Q、 A和对应压缩强度o数据， 通过数据拟合得到o=尸.Q + /"A，该拟合关系式中的系数尸和历分别对应Z 值计算式中的^和/。
步骤四、计算各孔隙形貌的分形维数以单个孔隙的形貌因子Z为标 度对象，釆用砂盒法对其进行分形分析得到综合描述被检测多孔材料孔隙 结构特征的孔形貌分形维数，计算分析过程包括以下步骤
(A) 以步骤一中所述二维孔隙结构图像的中心点为中心，以所述二 维孔隙结构图像的最小边长度为边长A的正方形一，将s,正方形中包含的 所有孔隙对应的Z值累加，得到一个对应的A^》；再以所述中心点为中心 将所述正方形一的边长缩小一倍得到边长为^ = |^的正方形二，将^正方
形中包含的所有孔隙对应的Z值累加，得到一个对应的A^(&);以此类推， 直到边长"对应的正方形N, n大于所述二维孔隙结构图像中最小孔隙的 直径，孔隙结构分析精度越高，"取值越�。蛔钪盏玫揭蛔槭軼(《)。
(B) 在以lnW,(s,)为纵坐标，ln《为横坐标的平面直角坐标系中，对一
组数据A^,)对应的一组数据点进行线性回归处理，得到一条直线，直线斜 率的绝对值就是被检测多孔材料孔形貌的分形维数A。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡
结构变化，均仍属于本发明技术方案的保护范围内。
权利要求
1. 一种多孔材料孔结构分形维数的计算方法，其特征在于该方法包括以下步骤步骤一、采用图像采集设备摄取被检测多孔材料外表面或一个剖面的二维孔隙结构图像，并将所采集到的二维孔隙结构图像传输至计算机进行分析处理，获得被检测多孔材料的二维孔隙边界图像，所述二维孔隙边界图像中各孔隙的二维孔隙边界为由多个像素点为顶点组成的平面多边形；步骤二、用计算机确定所述平面多边形各顶点的坐标值，其坐标原点为所述二维孔隙边界图像中的任一位置；步骤三、用计算机处理计算得出各孔隙的形貌因子Z，其计算机处理过程包括以下步骤(a)根据公式Q＝L2/S计算各孔隙的不规则因子Q，其中，L和S分别为单个孔隙的周长和面积；(b)计算各孔隙的锐度因子A对单个孔隙而言，首先确定该孔隙边界对应平面多边形的所有顶点的内角角度θ，并将所述所有顶点按照内角角度大小划分到多个角度区间中，对各角度区间内的顶点点位个数进行统计后，计算出该孔隙的锐度因子A；(c)计算各孔隙的形貌因子Z∶Z＝X·Q+y·A，式中x、y为系数，对于单个孔隙而言，首先通过试验分别分析孔隙弯曲情况和孔隙曲率半径对孔隙性能的影响程度，当孔隙弯曲情况为主要影响因素时，x取值大于y值；当孔隙曲率半径为主要影响因素时，x取值小于y值；步骤四、计算各孔隙形貌的分形维数以单个孔隙的形貌因子Z为标度对象，采用砂盒法对其进行分形分析得到综合描述被检测多孔材料孔隙结构特征的孔形貌分形维数，计算分析过程包括以下步骤(A)以步骤一中所述二维孔隙结构图像的中心点为中心，以所述二维孔隙结构图像的最小边长度为边长ε1的正方形一，将ε1正方形中包含的所有孔隙对应的Z值累加，得到一个对应的N(ε1)；再以所述中心点为中心将所述正方形一的边长缩小一倍得到边长为的正方形二，将ε2正方形中包含的所有孔隙对应的Z值累加，得到一个对应的N(ε2)；以此类推，直到边长εn对应的正方形N，εn大于所述二维孔隙结构图像中最小孔隙的直径，孔隙结构分析精度越高，εn取值越�。蛔钪盏玫揭蛔槭軳(εi)；(B)在以lnNi(εi)为纵坐标，lnεi为横坐标的平面直角坐标系中，对一组数据N(εi)对应的一组数据点进行线性回归处理，得到一条直线，直线斜率的绝对值就是被检测多孔材料孔形貌的分形维数DS。
2. 按照权利要求1所述的一种多孔材料孔结构分形维数的计算方法， 其特征在于步骤(b)中所述的计算各孔隙的锐度因子A时，根据加权 法，通过公式A^l^n,，计算出该孔隙的锐度因子A,其中n,为各角度区间内的顶点点位个数，a|为对应角度区间的加权系数。
3. 按照权利要求1或2所述的 一 种多孔材料孔结构分形维数的计算方 法，其特征在于步骤一所述的二维孔隙结构图像传输至计算机进行分析 处理时，所述计算机的分析处理过程包括以下步骤(C) 将所述二维孔隙结构图像转化为多灰度级灰度图像；(D) 对所述多灰度级灰度图像进行去噪和滤波处理；(E) 将经去噪和滤波处理的所述多灰度级灰度图像转换为二灰度级 图像即黑白图像；所述黑白图像由像素值为1的白色图像区域和像素值为0的黑色图像 区域组成，其中，所述被检测多孔材料的实体骨架材料对应的图像区域为 白色图像区域或黑色图像区域，所述孔隙空间对应的图像区域相应为黑色 图像区域或白色图像区域；(F) 用计算机提取所述黑白图像中黑色图像区域和白色图像区域的 交界线即孔隙空间边界并用彩色线条进行标示，所标示出来的孔隙空间边 界为由有限个像素点为顶点组成的平面多边形，确定所述平面多边形所有 顶点的二维坐标值并对其各顶点进行连续编号，即获得被检测多孔材料的一个或多个相互平行断面的二维孔隙边界。
4. 按照权利要求1或2所述的 一 种多孔材料孔结构分形维数的计算方 法，其特征在于步骤一中所述的图像釆集设备为金相显微镜或扫描电子 显微镜。
5. 按照权利要求1或2所述的一种多孔材料孔结构分形维数的计算方 法，其特征在于步骤一中所述的二维孔隙结构图像为分辨率不低于500 x 500的清晰彩色图像或多灰度级灰度图像。
6. 按照权利要求3所述的 一 种多孔材料孔结构分形维数的计算方法， 其特征在于步骤(D)中所述的滤波处理为最小值滤波和/或中值滤波和 /或最大值滤波和/或均值滤波和/或低通滤波和/或高通滤波处理。
7. 按照权利要求3所述的 一种多孔材料孔结构分形维数的计算方法， 其特征在于步骤(E)中所述的实体骨架材料对应的图像区域为白色图 像区域，所述孔隙空间对应的图像区域相应为黑色图像区域。
8. 按照权利要求3所述的 一种多孔材料孔结构分形维数的计算方法， 其特征在于步骤三中还需对所述多边形进行简化处理，所述简化处理的方法为减少所述多边形的边数并相应增加各边边长。
9. 按照权利要求3所述的 一 种多孔材料孔结构分形维数的计算方法， 其特征在于步骤(F)中确定所述平面多边形所有顶点的二维坐标值时，选取所述黑白图像左下方的顶点为坐标原点。
10. 按照权利要求3所述的一种多孔材料孔结构分形维数的计算方法， 其特征在于步骤(F)中对所述平面多边形各顶点进行连续编号时，以 二维坐标绝对值最小的顶点为1#顶点，然后按照逆时针顺序依次对剩余顶 点进行连续编号。
全文摘要
本发明公开了一种多孔材料孔结构分形维数的计算方法，该方法包括以下步骤获得被检测多孔材料的二维孔隙边界图像，所述二维孔隙边界图像中各孔隙的二维孔隙边界为由多个像素点为顶点组成的平面多边形；用计算机确定所述平面多边形各顶点的坐标值；用计算机处理计算得出各孔隙的形貌因子Z，其计算机处理过程包括以下步骤计算各孔隙的不规则因子Q、计算各孔隙的锐度因子A；计算各孔隙的形貌因子Z；计算各孔隙形貌的分形维数D<sub>S</sub>。本发明能有效解决现有分形维数不能完整表述被检测多孔材料复杂孔结构特征的实际问题。
文档编号G01N15/08GK101413871SQ20081023261
公开日2009年4月22日 申请日期2008年12月5日 优先权日2008年12月5日
发明者奚正平, 朱纪磊, 程 李, 汤慧萍, 王建永 申请人:西北有色金属研究院
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
