CN108682020A - 岩心微米ct孔隙结构重构方法 - Google Patents

Abstract

本发明属于非常规储层微观技术领域，具体为岩心微米CT孔隙结构重构方法，利用FIB‑SEM实验，标定微米CT结果中灰度所对应的实际孔隙度，从而建立微米CT图像灰度和孔隙度之间的关系，然后依据该关系，给微米CT图像的每个像素赋予一定的孔隙度，建立孔隙结构的软分割方法。本发明提供的岩心微米CT孔隙结构重构方法，提出综合FIB‑SEM和CT数据的数字岩心孔隙结构重构的软分割方法，建立分数式数字岩心孔隙结构模型，进一步提高数字岩心孔隙结构模型与实际岩心的匹配度。

Description

岩心微米CT孔隙结构重构方法

技术领域

本发明属于非常规储层微观技术领域，具体为岩心微米CT孔隙结构重构方法。

背景技术

随着常规油气资源的逐渐枯竭，页岩气、页岩油、致密油气等非常规油气已经成为当今 油气领域的研究热点之一。但是实际非常规油气开发表明，非常规油气表现出独特的产能递 减规律：初期产量高，衰减快；后期衰减慢，产量低。为了分析页岩气油气的流动规律，需 要开展数字岩心研究，以表征非常规储层微观三维孔-缝结构。

由于非常规储层中，20纳米到10微米左右的孔-缝均有发育，但常规微米CT的分辨率大于 100纳米，因此常规CT中的一个三维像素块，会出现3种情况：①全部都是固体(图4a中a4区)、 ②全部都是空隙(图4a中a3区)、③部分固体部分空隙(图4a中a1和a2区)。常规微米CT数 据(图4b)也有3种情况：(1)黑色区域(图4b中b3像素，对应于图4a中a3区域)；(2)白色区域(图4b中b4像素，对应于图4a中a4区域)；(3)灰色区域(图4b中b1和b2像素，对应 于图4a中a1和a2区域)。

图4a为岩心孔隙-固体分布；图4b为与a相对应的CT扫描图像。图4b中每个格为一个像素， 像素边长对应CT的分辨率(约0.7微米-20微米)。对图4b中的CT结果采用统一阈值进行硬分 割，是以整个像素为单位，b1像素和b2像素要么整体被判定为孔隙(孔隙度为Fp＝1)，要么 整体被判定为固体(孔隙度Fp＝0)，而实际上b1像素和b2像素所对应的岩心a1和a2区域是孔 隙和固体混合的，单个像素块的实际孔隙度Fp是从0到1变化的分数。

目前基于CT的数字岩心孔隙结构重构，普遍采用硬分割方法：对于CT中的任一三维像素， 确定单一分割阈值Th,三维CT数据体中灰度值大于Th的为孔-缝，赋值为0；而三维CT数据体 中灰度值小于Th的为固体部分，赋值为1。

采用硬分割方法进行数字岩心孔隙结构重构，对于部分孔隙部分固体的像素存在误差： 按照目前的这种孔隙分割方法，

图4b中b1和b2像素只能全部划分为孔隙，或者全部划分为固体，这与实际情况不符。

发明内容

针对硬分割方法进行数字岩心孔隙结构重构的问题，即由于微米CT分辨率有限，无法区 分分辨率以下的孔隙结构，采用统一的一个阈值对每个CT像素进行硬分割，CT像素要么全部 划定为孔隙，其像素孔隙度Fp＝1；要么全部划分为固体，其像素孔隙度Fp＝0。而实际岩心中 对应于CT像素块大小的区域，孔隙是以一定比例分数存在的，并非不是0，就是1，而是一个 分数式的孔隙度。

为了解决以上问题，本发明提出软分割方法，不再将微米CT像素只划分为孔隙和固体两 个类别，而是根据CT图像的灰度，结合高分辨率FIB-SEM实验，通过建立CT灰度和FIB-SEM 实验中实际岩心孔隙度之间的关系，给CT像素赋予分数式的像素孔隙度Fp。

具体技术方案为：

利用FIB-SEM实验，标定微米CT结果中灰度所对应的实际孔隙度，从而建立微米CT图像 灰度和孔隙度之间的关系，然后依据该关系，给微米CT图像的每个像素赋予一定的孔隙度， 建立孔隙结构的软分割方法。

岩心微米CT孔隙结构重构方法，包括以下步骤：

(1)微米CT岩心样品制备

将采集到的原始岩石样品进行表面机械剖光，然后在需要开展微米CT分析的区域，钻取 微米CT样品CT_SAMPLE。该CT样品半径CTR毫米,轴长CTZ毫米；微米CT岩心样品中机械 剖光面的一端标记为样品A端，另外一端为B端；

(2)微米CT岩心样品表面离子剖光

对微米CT岩心样品A端，进行氩离子表面剖光，剖光区域为整个A端端面；然后对离子 剖光后的A端端面进行SEM扫描，拍摄A端端面的表面结构图片；

(3)打定位孔和CT扫描

根据A端SEM扫描图片结果，在A端表面确定一个边长为SEMQ微米的正方形区域SQU； 在该正方形区域SQU内，选取3个定位点PA、PB、PC，要求这3个点构成一个直角坐标关系， 针对打了定位孔的CT样品CT_SAMPLE进行微米CT扫描，扫描图像精度为CTN微米，扫描结 果为3维像素矩阵CT_DATA；

所述的步骤3，包括以下步骤：

(31)根据A端SEM扫描图片结果，在A端表面确定一个边长为SEMQ微米的正方形区域 SQU；

(32)然后在该正方形区域SQU内，选取3个定位点PA、PB、PC，要求这3个点构成一个直角坐标关系；

(33)利用聚焦离子束在3个定位点的位置打3种不同形状的定位孔：

PA标定点打圆柱形孔，孔的直径为3*CTN微米，轴长3*CTN微米；PB标定点打正方体孔， 正方体的边长为3*CTN微米；PC标定孔打三棱柱孔，表面为正三角形，边长为3*CTN微米， 纵向深3*CTN微米；

(34)针对打了定位孔的CT样品CT_SAMPLE进行微米CT扫描，扫描图像精度为CTN微米，扫描结果为3维像素矩阵CT_DATA。

(4)FIB-SEM扫描

针对打了定位孔的CT样品CT_SAMPLE进行FIB-SEM扫描，扫描区域为SQU区域，纵向切入深度为SEMQ微米，横向FIB烧蚀间隔为FIBN纳米；FIB-SEM扫描的精度为FIBN纳米；FIB-SEM扫描后，得到一个边长为SEMQ微米的立方体的岩心三维像素矩阵FIB-SEM_DATA；

(5)微米CT和FIB-SEM数据匹配

步骤(3)中CT扫描结果中可以观察到3个定位孔，CT数据中的3个定位孔为PA^CT、PB^CT、 PC^CT，同时步骤(4)FIB-SEM扫描结果中同样可以观察到3个定位孔，FIB-SEM数据中3个定 位孔为PA^FIB、PB^FIB、PC^FIB；PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB是一一对应的，根据 PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB，将CT扫描数据CT_DATA和FIB-SEM扫描数据 FIB-SEM_DATA进行位置匹配；

所述的步骤(5)，包括以下步骤：

(51)确定CT数据中3个定位孔的中心，然后以PA^CT孔的中心为坐标系原点，以PA^CT-PB^CT为x轴，以PA^CT-PC^CT为y轴，以垂直于剖光面方向为z轴，建立CT数据坐标系CT_SYSTEM；

(52)确定FIB-SEM数据中3个定位孔的中心，然后以PA^FIB孔的中心为坐标系原点，以 PA^FIB-PB^FIB为x轴，以PA^FIB-PC^FIB为y轴，以垂直于剖光面方向为z轴，建立FIB-SEM数据坐标 系FIB_SYSTEM；

(53)PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB是一一对应的，CT数据坐标系CT_SYSTEM 与FIB-SEM数据坐标系FIB_SYSTEM中坐标只相差一个缩放系数SF＝CTN/FIBN*10⁶；

(54)在FIB-SEM扫描区域与CT扫描重合区域，CT数据中任一像素P^CT在CT数据坐标系 CT_SYSTEM中的坐标为P^CT_X、P^CT_Y和P^CT_Z；由于FIB-SEM比CT数据的分辨率高，因此 CT数据中任一像素P^CT对应的FIB-SEM数据中的一个像素块P^FIB，P^FIB的x轴范围为P^CT_X*SF 至(P^CT_X+1)*SF,y轴范围为P^CT_Y*SF至(P^CT_Y+1)*SF，z轴范围为P^CT_Z*SF至(P^CT_Z+1) *SF。

(6)FIB-SEM孔隙度分析

通过步骤(5)，在FIB-SEM高分辨率数据中建立与CT数据一一对应的像素块P^FIB；对于P^FIB中的FIB-SEM数据，由于其分辨率较高，采用硬分割方式，确定P^FIB中的FIB-SEM的孔隙部分 孔隙度FP^FIB；

所述的步骤(6)，包括以下步骤：

(61)确定一个初始分割阈值Th*^,FIB，P^FIB中的FIB-SEM图像中小于Th*^,FIB的为孔隙，大 于Th*^,FIB为固体骨架，然后不断调整分割阈值，直到分割完毕的孔隙边缘与P^FIB中的原始 FIB-SEM图像的孔隙边缘相吻合位置，并以此时的分割阈值为最优阈值Th^FIB；

(62)确定了最优分割阈值之后，对P^FIB中FIB-SEM图像进行分割，统计划分孔隙的像素 和为Por^FIB，统计P^FIB的总像素数目为Pix^FIB,则P^FIB的孔隙度为FP^FIB＝Por^FIB/Pix^FIB。

(7)微米CT数据灰度和FIB-SEM孔隙度关系分析

通过步骤(5)，建立CT数据P^CT与FIB-SEM高分辨率数据P^FIB之间的一一对应关系，通过 步骤(6)可以得到P^FIB的孔隙度FP^FIB，从而可以利用最小二乘法建立P^CT的灰度G^CT与孔隙度 FP^FIB之间的多项式拟合关系FP^FIB＝f(G^CT)；

(8)微米CT软分割建立孔隙结构模型

应用步骤(7)拟合出的CT图像像素灰度与孔隙率之间的关系FP^FIB＝f(G^CT)，带入CT图 像灰度数值，计算CT扫描数据中其他数据的孔隙率，得出全部CT扫描数据中任一像素的孔隙 率FP；建立一个与CT扫描数据CT_DATA相同大小并且一一对应的像素矩阵FP_DATA，FP_DATA的数值即为相应CT_DATA中相应像素的孔隙率FP；FP_DATA即为最后输出的软分 割孔隙结构重构模型。

本发明提供的岩心微米CT孔隙结构重构方法，提出综合FIB-SEM和CT数据的数字岩心孔 隙结构重构的软分割方法，建立分数式数字岩心孔隙结构模型，进一步提高数字岩心孔隙结 构模型与实际岩心的匹配度。

附图说明

图1为实施例中CT扫描样品A端面定位点和定位孔示意图；

图2为实施例中CT数据坐标系CT_SYSTEM；

图3为实施例中FIB-SEM数据坐标系FIB_SYSTEM；

图4a为现有技术中岩心孔隙-固体分布；

图4b为与图4a相对应的CT扫描图像。

具体实施方式

结合附图说明本发明的具体技术方案。

(1)微米CT岩心样品制备

所选微米CT仪器的分辨率为CTN微米(一般为1微米)，CT仪器所需样品为圆柱形柱塞, 柱塞样品半径要求为CTR毫米，轴长要求为CTZ毫米。将采集到的原始岩石样品进行表面机 械剖光，然后在需要开展微米CT分析的区域，钻取微米CT样品CT_SAMPLE。该CT样品半径 CTR毫米,轴长CTZ毫米。微米CT岩心样品中机械剖光面的一端标记为样品A端，另外一端为B端。

(2)微米CT岩心样品表面离子剖光

对微米CT岩心样品A端，进行氩离子表面剖光，剖光区域为整个A端端面。然后对离子 剖光后的A端端面进行SEM扫描，拍摄A端端面的表面结构图片。

(3)打定位孔和CT扫描

(31)根据A端SEM扫描图片结果，在A端表面确定一个边长为SEMQ微米的正方形区域 SQU，如图1中FIB-SEM扫描区域；

(32)然后在该正方形区域SQU内，选取3个定位点PA、PB、PC,要求这3个点构成一个直角坐标关系，见图1：

(33)利用聚焦离子束(FIB)在3个定位点的位置打3种不同形状的定位孔见图1：

PA标定点打圆柱形孔，每个孔的直径为3*CTN微米,轴长3*CTN微米；PB标定点打正方 体孔，正方体的边长为3*CTN微米；PC标定孔打三棱柱孔，表面为正三角形，边长为3*CTN微米，纵向深3*CTN微米。由于PA、PB、PC孔的尺寸均大于CT分辨率CTN微米，从而保证 在微米CT中能够显示出这3个标定孔。

(34)针对打了定位孔的CT样品CT_SAMPLE进行微米CT扫描，扫描图像精度为CTN微米，扫描结果为3维像素矩阵CT_DATA。

(4)FIB-SEM扫描

针对打了定位孔的CT样品CT_SAMPLE进行FIB-SEM扫描，扫描区域为SQU区域，纵向切入深度为SEMQ微米，横向FIB烧蚀间隔为FIBN纳米(FIBN一般为10纳米)。FIB-SEM扫 描的精度为FIBN纳米。FIB-SEM扫描后，得到一个边长为SEMQ微米的立方体的岩心三维像 素矩阵FIB-SEM_DATA。

(5)微米CT和FIB-SEM数据匹配

由于步骤(3)中CT扫描结果中可以观察到3个定位孔，CT数据中的3个定位孔为PA^CT、 PB^CT、PC^CT，同时FIB-SEM扫描结果中同样可以观察到3个定位孔，FIB-SEM数据中3个定位 孔为PA^FIB、PB^FIB、PC^FIB。由于PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB都是实际岩心中3 个定位孔PA、PB、PC的扫描结果，因此PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB是一一对 应的，可以根据PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB，将CT扫描数据CT_DATA和FIB-SEM 扫描数据FIB-SEM_DATA进行位置匹配，具体匹配步骤如下：

(51)确定CT数据中3个定位孔的中心，然后以PA^CT孔的中心为坐标系原点，以PA^CT-PB^CT为x轴，以PA^CT-PC^CT为y轴，以垂直于剖光面方向为z轴，建立CT数据坐标系CT_SYSTEM，见图2。

(52)确定FIB-SEM数据中3个定位孔的中心，然后以PA^FIB孔的中心为坐标系原点，以 PA^FIB-PB^FIB为x轴，以PA^FIB-PC^FIB为y轴，以垂直于剖光面方向为z轴，建立FIB-SEM数据坐标 系FIB_SYSTEM，见图3。

(53)由于PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB是一一对应的，因此CT数据坐标系 CT_SYSTEM与FIB-SEM数据坐标系FIB_SYSTEM中坐标只相差一个缩放系数 SF＝CTN/FIBN*10⁶。

(54)在FIB-SEM扫描区域与CT扫描重合区域，如图1中灰色区域，CT数据中任一像素 P^CT在CT数据坐标系CT_SYSTEM中的坐标为P^CT_X、P^CT_Y和P^CT_Z。由于FIB-SEM比CT数据 的分辨率高，因此CT数据中任一像素P^CT对应的FIB-SEM数据中的一个像素块P^FIB，P^FIB的x轴 范围为P^CT_X*SF至(P^CT_X+1)*SF，y轴范围为PCT_Y*SF至(PCT_Y+1)*SF，z轴范围为 PCT_Z*SF至(PCT_Z+1)*SF。

(6)FIB-SEM孔隙度分析

通过步骤(5)，可以在FIB-SEM高分辨率数据中建立与CT数据一一对应的像素块P^FIB。 对于P^FIB中的FIB-SEM数据，由于其分辨率较高，采用硬分割方式，确定P^FIB中的FIB-SEM的 孔隙部分。首先确定一个初始分割阈值Th*^,FIB，P^FIB中的FIB-SEM图像中小于Th*^,FIB的为孔隙， 大于Th*^,FIB为固体骨架，然后不断调整分割阈值，直到分割完毕的孔隙边缘与P^FIB中的原始 FIB-SEM图像的孔隙边缘相吻合位置，并以此时的分割阈值为最优阈值Th^FIB。确定了最优分 割阈值之后，对P^FIB中FIB-SEM图像进行分割，统计划分孔隙的像素和为Por^FIB，统计P^FIB的总 像素数目为Pix^FIB，则P^FIB的孔隙度为FP^FIB＝Por^FIB/Pix^FIB。

(7)微米CT数据灰度和FIB-SEM孔隙度关系分析

通过步骤(5)，可以建立CT数据P^CT与FIB-SEM高分辨率数据P^FIB之间的一一对应关系。 通过步骤(6)可以得到P^FIB的孔隙度FP^FIB，从而可以利用最小二乘法建立P^CT的灰度G^CT与孔 隙度FP^FIB之间的多项式拟合关系FP^FIB＝f(G^CT)。

(8)微米CT软分割建立孔隙结构模型

应用步骤(7)拟合出的CT图像像素灰度与孔隙率之间的关系FP^FIB＝f(G^CT)，带入CT图 像灰度数值，计算CT扫描数据中其他数据(即不在FIB-SEM扫描区域与CT扫描重合区域)的 孔隙率，得出全部CT扫描数据中任一像素的孔隙率FP。建立一个与CT扫描数据CT_DATA相 同大小并且一一对应的像素矩阵FP_DATA，FP_DATA的数值即为相应CT_DATA中相应像素的孔隙率FP。FP_DATA即为最后输出的软分割孔隙结构重构模型。

例如基于蔡司微米CT仪得到的岩心微米CT图像分析中，一般以灰度值为50-60作为孔隙的 划分标准。CT图像中灰度小于50-60的为孔隙，大于50-60的固体骨架。而实际上灰度值大于 50-60的部分，并不全都是固体骨架，而是存在一部分孔隙。通过本技术方案，可以给出灰度 大于50-60部分的像素所对应的孔隙率，不再将灰度大于50-60的像素全部划分为固体骨架。 并且本技术方案是综合了FIB-SEM高精度扫描数据建立的定量分析方法。

Claims (4)

1.岩心微米CT孔隙结构重构方法，其特征在于，包括以下步骤：

(1)微米CT岩心样品制备

将采集到的原始岩石样品进行表面机械剖光，然后在需要开展微米CT分析的区域，钻取微米CT样品CT_SAMPLE。该CT样品半径CTR毫米，轴长CTZ毫米；微米CT岩心样品中机械剖光面的一端标记为样品A端，另外一端为B端；

(2)微米CT岩心样品表面离子剖光

对微米CT岩心样品A端，进行氩离子表面剖光，剖光区域为整个A端端面；然后对离子剖光后的A端端面进行SEM扫描，拍摄A端端面的表面结构图片；

(3)打定位孔和CT扫描

根据A端SEM扫描图片结果，在A端表面确定一个边长为SEMQ微米的正方形区域SQU；在该正方形区域SQU内，选取3个定位点PA、PB、PC，要求这3个点构成一个直角坐标关系，针对打了定位孔的CT样品CT_SAMPLE进行微米CT扫描，扫描图像精度为CTN微米，扫描结果为3维像素矩阵CT_DATA；

(4)FIB-SEM扫描

针对打了定位孔的CT样品CT_SAMPLE进行FIB-SEM扫描，扫描区域为SQU区域，纵向切入深度为SEMQ微米，横向FIB烧蚀间隔为FIBN纳米；FIB-SEM扫描的精度为FIBN纳米；FIB-SEM扫描后，得到一个边长为SEMQ微米的立方体的岩心三维像素矩阵FIB-SEM_DATA；

(5)微米CT和FIB-SEM数据匹配

步骤(3)中CT扫描结果中可以观察到3个定位孔，CT数据中的3个定位孔为PA^CT、PB^CT、PC^CT，同时步骤(4)FIB-SEM扫描结果中同样可以观察到3个定位孔，FIB-SEM数据中3个定位孔为PA^FIB、PB^FIB、PC^FIB；PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB是一一对应的，根据PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB，将CT扫描数据CT_DATA和FIB-SEM扫描数据FIB-SEM_DATA进行位置匹配；

(6)FIB-SEM孔隙度分析

通过步骤(5)，在FIB-SEM高分辨率数据中建立与CT数据一一对应的像素块P^FIB；对于P^FIB中的FIB-SEM数据，由于其分辨率较高，采用硬分割方式，确定P^FIB中的FIB-SEM的孔隙部分孔隙度FP^FIB；

(7)微米CT数据灰度和FIB-SEM孔隙度关系分析

通过步骤(5)，建立CT数据P^CT与FIB-SEM高分辨率数据P^FIB之间的一一对应关系，通过步骤(6)可以得到P^FIB的孔隙度FP^FIB，从而可以利用最小二乘法建立P^CT的灰度G^CT与孔隙度FP^FIB之间的多项式拟合关系FP^FIB＝f(G^CT)；

(8)微米CT软分割建立孔隙结构模型

应用步骤(7)拟合出的CT图像像素灰度与孔隙率之间的关系FP^FIB＝f(G^CT)，带入CT图像灰度数值，计算CT扫描数据中其他数据的孔隙率，得出全部CT扫描数据中任一像素的孔隙率FP；建立一个与CT扫描数据CT_DATA相同大小并且一一对应的像素矩阵FP_DATA，FP_DATA的数值即为相应CT_DATA中相应像素的孔隙率FP；FP_DATA即为最后输出的软分割孔隙结构重构模型。

2.根据权利要求1所述的岩心微米CT孔隙结构重构方法，其特征在于，所述的步骤3，包括以下步骤：

(31)根据A端SEM扫描图片结果，在A端表面确定一个边长为SEMQ微米的正方形区域SQU；

(32)然后在该正方形区域SQU内，选取3个定位点PA、PB、PC，要求这3个点构成一个直角坐标关系；

(33)利用聚焦离子束在3个定位点的位置打3种不同形状的定位孔：

PA标定点打圆柱形孔，孔的直径为3*CTN微米，轴长3*CTN微米；PB标定点打正方体孔，正方体的边长为3*CTN微米；PC标定孔打三棱柱孔，表面为正三角形，边长为3*CTN微米，纵向深3*CTN微米；

(34)针对打了定位孔的CT样品CT_SAMPLE进行微米CT扫描，扫描图像精度为CTN微米，扫描结果为3维像素矩阵CT_DATA。

3.根据权利要求1所述的岩心微米CT孔隙结构重构方法，其特征在于，所述的步骤(5)，包括以下步骤：

(51)确定CT数据中3个定位孔的中心，然后以PA^CT孔的中心为坐标系原点，以PA^CT-PB^CT为x轴，以PA^CT-PC^CT为y轴，以垂直于剖光面方向为z轴，建立CT数据坐标系CT_SYSTEM；

(52)确定FIB-SEM数据中3个定位孔的中心，然后以PA^FIB孔的中心为坐标系原点，以PA^FIB-PB^FIB为x轴，以PA^FIB-PC^FIB为y轴，以垂直于剖光面方向为z轴，建立FIB-SEM数据坐标系FIB_SYSTEM；

(53)PA^CT、PB^CT、PC^CT和PA^FIB、PB^FIB、PC^FIB是一一对应的，CT数据坐标系CT_SYSTEM与FIB-SEM数据坐标系FIB_SYSTEM中坐标只相差一个缩放系数SF＝CTN/FIBN*10⁶；

(54)在FIB-SEM扫描区域与CT扫描重合区域，CT数据中任一像素P^CT在CT数据坐标系CT_SYSTEM中的坐标为P^CT_X、P^CT_Y和P^CT_Z；由于FIB-SEM比CT数据的分辨率高，因此CT数据中任一像素P^CT对应的FIB-SEM数据中的一个像素块P^FIB，P^FIB的x轴范围为P^CT_X*SF至(P^CT_X+1)*SF，y轴范围为P^CT_Y*SF至(P^CT_Y+1)*SF，z轴范围为P^CT_Z*SF至(P^CT_Z+1)*SF。

4.根据权利要求1所述的岩心微米CT孔隙结构重构方法，其特征在于，所述的步骤(6)，包括以下步骤：

(61)确定一个初始分割阈值Th*^,FIB，P^FIB中的FIB-SEM图像中小于Th*^,FIB的为孔隙，大于Th*^，FIB为固体骨架，然后不断调整分割阈值，直到分割完毕的孔隙边缘与P^FIB中的原始FIB-SEM图像的孔隙边缘相吻合位置，并以此时的分割阈值为最优阈值Th^FIB；

(62)确定了最优分割阈值之后，对P^FIB中FIB-SEM图像进行分割，统计划分孔隙的像素和为Por^FIB，统计P^FIB的总像素数目为Pix^FIB，则P^FIB的孔隙度为FP^FIB＝Por^FIB/Pix^FIB。