共计 8702 个字符，预计需要花费 22 分钟才能阅读完成。

分析了好久的BET数据，每次在科学指南针或者一些高校科研院所，或者企业里拿到测试结果，都会第一眼看看有没有原始记录数据。慢慢发现，市面上的 BET 分析软件、仪器厂商软件和报告格式非常不统一。不同软件的默认拟合区间、厚度方程和结果展示方式各不相同，让不同样品之间的横向对比对我来说变成一件苦差事。有人说直接看输出的PDF和excel里面的结果不就可以了么？然而实际并非如此，软件内部的参数可调性极高，每种软件中都有所规定的固定的一套判断处理方式。为此我开发了一一套BET综合分析软件。

BET综合分析软件下载（点击图片直接下载）

目前支持的厂家和设备（仍在github中实时更新）：

厂家 / 软件体系	目前已支持或验证的设备 / 版本	支持导入格式	当前支持状态
Micromeritics	TriStar II 3020	SMP	支持等温线、BET、Langmuir、t-Plot、BJH；保留 TriStar II 3020 的默认点数和常数差异
Micromeritics	MicroActive for TriStar II Plus	SMP、XLS、XLSX、XLSM	支持原始 SMP 和官方 Excel；BET 可读取官方选点区间，BJH 按 MicroActive 口径继续修正
Micromeritics	3Flex 3500 / Flex 6.03	SMP、XLS	支持 3Flex SMP 手动等温线点表和 Flex 6.03 官方 XLS 报表
Micromeritics	ASAP 2460	SMP、XLS、XLSX、XLSM	支持 ASAP 等温线和官方 Excel 导出；BET 默认区间含 ASAP 2460 特定修正
Micromeritics	ASAP 2020 Plus	SMP、XLS、XLSX、XLSM	支持原始或导出数据读取，并进入统一分析流程
MicrotracBEL	BELSORP / BELMaster	DAT、XLS、XLSX、XLSM	支持 BELMaster 等温线导入，并进入 BET / Langmuir / t-Plot / BJH 分析
Quantachrome	Autosorb iQ / QuadraSorb / NovaWin 11.02 文本型 Excel 报告	QPS、XLSX	支持 QPS 原始数据和 NovaWin 文本型 Excel 报告读取
贝士德	BSD-660MC，软件版本 V.9.1.15.0 Date 26.04.28	XLS、XLSX、XLSM	支持官方 Excel；BET、Langmuir、t-Plot 已按 BSD 报表口径复现，BJH 读取官方逐点表
精微高博	JWGB 官方 Excel，多 sheet 报表	XLSX	支持 Info、Isotherm、BET、Langmuir、t-Plot、BJH 等 sheet 读取与统一分析

软件最大优势：软件支持多样品等温线叠加显示，可以同时观察不同样品的吸附/脱附曲线差异；BET、Langmuir、t-Plot、BJH、DH、HK、DFT （可实时调整正则化系数同时应用全部样品，实时预览）等分析图也支持多样品对比。这样不仅能看到某一个样品的结果，还能快速判断一组样品在比表面积、孔容、孔径分布和微孔结构上的变化趋势。可后台自动检查软件更新！

算法分析（极为干货）

BET 分析软件真正难的地方，并不是把一条等温线画出来，而是把不同仪器、不同软件版本、不同默认选点规则、不同厚度方程和不同孔径分布修正方式，放到同一个可追溯的算法框架里。我的 BET 综合分析软件正是围绕这个目标开发的：让原始数据、拟合区间、计算公式、参数来源、厂家默认规则和最终结果之间的关系变得透明。

一、BET 算法：软件到底如何计算比表面积？

BET 模型来自 Brunauer、Emmett 和 Teller 的多分子层吸附理论。它比 Langmuir 模型更进一步，因为 Langmuir 只考虑单分子层吸附，而 BET 认为第一层吸附和后续多层吸附的能量不同，因此可以描述相对压力较低区域内的多层吸附行为。

BET 原始方程可写为：

V = Vm * C * x / [(1 - x) * (1 + (C - 1) * x)]

其中：x = P/P0，V 是相对压力 x 下的吸附量，Vm 是单层吸附量，C 是 BET 常数，反映第一层吸附热与液化热之间的差异。

为了让软件能稳定拟合，BET 方程通常被线性化：

x / [V * (1 - x)] = 1 / (Vm * C) + [(C - 1) / (Vm * C)] * x

令：

Y = x / [V * (1 - x)]

则软件对 Y 和 x 做线性回归：

Y = slope * x + intercept

拟合完成后，BET 参数由以下公式得到：

Vm = 1 / (slope + intercept)

C = slope / intercept + 1

比表面积计算为：

S_BET = Vm * NA * sigma * 1e-18 / 22414

其中 S_BET 单位为 m²/g，Vm 单位为 cm³/g STP，NA 为阿伏伽德罗常数，sigma 是吸附质分子截面积。对于 77 K 氮气，常用 sigma = 0.162 nm²。软件会优先读取文件中的吸附质参数，如果文件没有提供，则使用默认氮气参数。

BET 选点为什么比公式更关键？

BET 的计算结果高度依赖拟合区间。相同公式，如果选点不同，比表面积和 C 常数可能相差明显。因此我的软件不是简单固定 0.05–0.30，而是把选区作为核心算法对象。

当前软件在自动选区时主要考虑：

1. 有效点约束：吸附量必须大于 0，相对压力必须满足 0 < P/P0 < 1，并且至少有 3 个有效点。

2. C 常数约束：BET 常数 C 必须为正。若 C <= 0，软件会标记为警告，因为这通常表示所选区间不符合 BET 物理意义。

3. Rouquerol 型约束：软件检查 V * (1 - P/P0) 在候选区间内是否随压力增加而单调不下降。这个约束对应 Rouquerol 判据中的重要思想：BET 选区应位于单层容量形成前后的合理区域，而不是随便拿一段线性看起来很好的数据。

4. 线性质量：候选区间需要有较高 R²，但软件不会只看 R²。因为吸附等温线中很多局部区间都可能线性很好，但不一定具有正确物理意义。

5. 压力区间偏好：软件默认优先搜索低压段，一般搜索到 0.35，必要时放宽到 0.50；评分时会轻微偏好接近常见 BET 区间的范围，例如低端接近 0.15、高端接近 0.30，并惩罚过窄区间。

不同软件版本的 BET 差异

来源/版本	当前软件中的处理方式	为什么要区别处理
TriStar II 3020	默认取前 7 个吸附点作为 BET 初始区间，并沿用历史口径的阿伏伽德罗常数 6.023e23。	老版本软件在默认点数和物理常数上与现代软件略有差异，若全部统一成新常数，会造成与原报告不一致。
MicroActive for TriStar II Plus	优先读取官方文件或 Excel 中保存的 BET 区间；若无官方区间，则在 0.08–0.30 内按 Rouquerol 型规则截取。	MicroActive 通常保存了软件选择或人工选择的点号，复现原报告时应优先尊重这些信息。
ASAP 2460 / ASAP 2020 Plus	优先使用官方拟合区间；ASAP 2460 若官方下限为 0 或异常低值，软件会按当前实现修正为 0.08 起算。	ASAP 系列报告中低压起点的处理方式与普通固定区间不同，需要避免把过低点错误纳入 BET 区间。
3Flex / Flex 6.03	优先读取官方区间；Langmuir 可能使用完整吸附分支；t-Plot 默认厚度范围使用 0.35–0.50 nm。	Flex 报告体系中拟合区间、自由空间修正和 BJH 输出表都有自己的报表口径。
BSD / JWGB	优先读取官方 Excel 中的区间或点号；t-Plot 的纵轴口径按厂家报表使用吸附量 STP，再换算表面积。	这些国产软件的报表口径和 Micromeritics 并不完全相同，不能直接套用同一个 t-Plot 纵轴定义。

二、Langmuir：它和 BET 到底有什么区别？

Langmuir 模型假设吸附只发生在均一表面上的单分子层，所有吸附位点等价，且吸附分子之间没有相互作用。其基本形式为：

V = Vm * b * x / (1 + b * x)

其中 b 是 Langmuir 吸附平衡常数。线性化后：

x / V = 1 / (Vm * b) + x / Vm

软件对 x/V 和 x 做线性拟合：

x / V = slope * x + intercept

因此：

Vm = 1 / slope

b = slope / intercept

S_Langmuir = Vm * NA * sigma * 1e-18 / 22414

Langmuir 和 BET 最大的区别在于：Langmuir 是单层模型，BET 是多层模型。对于非微孔材料的氮气物理吸附，比表面积通常以 BET 为标准；Langmuir 面积常作为辅助对比。很多软件中的 Langmuir 计算确实没有 BET 那么复杂，原因不是软件偷懒，而是 Langmuir 模型本身没有 Rouquerol 选区、C 常数和多层吸附判据这些细节。

在我的软件中，Langmuir 默认区间一般为 0.05–0.30，但如果数据文件保存了官方区间，则优先使用官方区间。对于某些 Flex 数据，Langmuir 可能按完整吸附分支拟合，这与 BET 的严格选区逻辑不同。

三、t-Plot：如何分离外比表面积和微孔体积？

t-Plot 的思想是：如果样品表面是非微孔开放表面，吸附量与吸附膜厚度 t 应近似线性；如果存在微孔填充，低压段会出现额外吸附量，线性外推截距可用于估算微孔体积。

软件先把相对压力转换为统计膜厚：

t = f(P/P0)

再将吸附量换算为液体体积：

V_liq = V_STP * density_conversion_factor

氮气默认换算因子为：

density_conversion_factor = 0.0015468

然后拟合：

V_liq = slope * t + intercept

当 V_liq 单位为 cm³/g，t 单位为 nm 时：

S_external = slope * 1000

V_micropore = max(0, intercept)

这里乘以 1000 的原因是单位换算：V = S * t，m²/g 与 nm 转成 cm³/g 时有 1 m² = 10⁴ cm²、1 nm = 10⁻⁷ cm，因此 V(cm³/g) = S(m²/g) * t(nm) * 10⁻³。

软件支持的厚度方程

Harkins-Jura：

t(nm) = 0.1 * [13.99 / (0.034 - log10(P/P0))]^0.5

Halsey：

t(nm) = 0.1 * 3.54 * [-5 / ln(P/P0)]^0.333

Kruk-Jaroniec-Sayari：

t(nm) = 0.1 * [60.65 / (0.03071 - log10(P/P0))]^0.3968

Broekhoff-de Boer：软件通过数值求解厚度 t，满足：

log10(P/P0) = -16.11 / t² + 0.1682 * exp(-0.1137 * t)

这里 t 先以 Å 计算，再换算成 nm。

Carbon Black STSA：

t(nm) = 0.1 * [2.98 + 6.45x + 0.88x²]

此外，软件还支持参考厚度曲线，使用参比点插值。这一点很重要，因为不同厂家软件使用的厚度曲线并不完全一样，t-Plot 的外比表面积和微孔体积对厚度方程非常敏感。

四、BJH：介孔孔径分布如何从等温线中计算出来？

BJH 用于介孔分布分析，核心思想是 Kelvin 方程加吸附膜厚度修正。软件对吸附或脱附分支逐点处理，先计算 Kelvin 半径，再加上吸附膜厚得到孔半径。

Kelvin 方程可写为：

rK = -A / ln(P/P0)

其中 A 与液氮表面张力、液体摩尔体积、温度和气体常数有关。软件中氮气 77.35 K 的默认等效因子为 0.953 nm，并按温度做比例修正：

A_T = A_77.35 * 77.35 / T

孔半径：

rp = rK + t

孔径：

D = 2 * rp

吸附量先转换成液体体积：

V_liq = V_STP * density_conversion_factor

每个压力区间的孔体积增量不是简单相减，因为在压力下降或上升时，已经存在的孔壁吸附膜厚度也在变化。因此标准 BJH 需要壁面修正：

volume_step = V_liq(high) - V_liq(low)

wall_volume = 已打开孔中吸附膜厚度变化贡献

core_volume = volume_step - wall_volume

delta_Vpore = core_volume * (rp / rK)^2

孔面积由圆柱孔近似得到：

delta_A = 4000 * delta_V / D

其中 delta_V 单位为 cm³/g，D 单位为 nm，得到 delta_A 单位为 m²/g。

常见显示方式包括：

dV/dlogD = delta_V / |log10(D_high) - log10(D_low)|

dV/dD = delta_V / |D_high - D_low|

软件还提供平滑处理：在 logD 网格上用 Akima 插值，并使用 9 点平滑导数计算微分孔容。这样做的目的是减少逐点实验噪声对微分孔径分布曲线的放大。

BJH 中不同软件的差异

Micromeritics、Quantachrome、BSD、JWGB 等软件在 BJH 上的差别可能比 BET 更明显。差异包括 Kelvin 因子、液体体积换算、最小孔径截断、吸附/脱附分支、厚度方程、标准修正或 Faas 修正、是否读取官方逐点表、是否平滑微分曲线等。

五、DH：它和 BJH 的区别是什么？

DH，即 Dollimore-Heal 方法，同样用于孔径分布分析。它和 BJH 一样基于毛细凝聚、Kelvin 半径和吸附膜厚度，但在区间平均孔径、孔长或孔壁处理上与 BJH 有不同近似。

在当前软件实现中，DH 与 BJH 共用同一套基础计算框架：

rK = -A / ln(P/P0)

D = 2 * (rK + t)

V_liq = V_STP * density_conversion_factor

差别主要在区间孔径的定义上。BJH 标准修正时会使用更接近 Micromeritics 报表口径的平均孔径；DH 模式下当前软件使用：

D_interval = 0.5 * (D_high + D_low)

并采用更低的默认最小孔径下限，例如当前 DH 默认最小孔径为 0.75 nm。如果文件中含有官方 DH 表，软件会优先读取官方 DH 表作为结果来源。

六、HK：微孔孔径为什么不能只靠 Kelvin 方程？

Horvath-Kawazoe 方法主要用于微孔，尤其是小于 2 nm 的孔。微孔中吸附不是简单的毛细凝聚，孔壁势场会显著影响吸附分子，因此 Kelvin 方程不再充分。

HK 方法的核心是把相对压力与吸附势能联系起来：

ln(P/P0) = -DeltaG / (R * T)

软件中对 slit、cylinder、sphere 三种孔型分别计算吸附势，并通过数值方法反求孔宽。以 slit pore 为例，软件使用的形式可概括为：

ln(P/P0) = [NA / (R*T)] * [A / (sigma^4 * (L - 2d0))] * F(L, d0, sigma)

其中 L 是核到核孔宽，d0 与吸附质和吸附剂分子直径有关，sigma 与零势能直径有关，A 是吸附剂-吸附质相互作用参数，F 包含吸引项和排斥项。

软件允许选择吸附剂和吸附质预设，例如 Zeolite、Aluminophosphate、Carbon-Graphite、N2、Ar、CO2 等。每组预设包含分子直径、零势能直径、极化率、磁化率和表面密度。相互作用参数可以手动输入，也可以由这些物性参数计算：

A = density_adsorptive * dispersion_adsorptive + density_adsorbent * dispersion_sample

HK 的输出不是介孔 BJH 那种 Kelvin 孔径，而是由吸附势反推出的微孔宽度。软件会逐点求解孔宽，并由液体体积变化得到累计孔容和微分孔容。

七、DFT / NLDFT：为什么这是更高级的孔径分布方法？

DFT 和 NLDFT 不再把孔看成简单圆柱或简单 slit 的 Kelvin 几何问题，而是用统计热力学模型描述不同孔宽下的吸附等温线。软件中的数学形式可以写成离散反问题：

q_i = sum_j K(P_i, w_j) * V_j

其中 q_i 是第 i 个压力点的吸附量，K(P_i, w_j) 是孔宽 w_j 在压力 P_i 下的理论吸附核，V_j 是对应孔宽区间的孔体积增量。

这是一个典型的病态反问题，不能直接普通最小二乘，否则结果会强烈振荡。因此软件采用非负正则化反演：

min ||K * V - q||² + lambda * ||D2 * V||²

约束条件：

V_j >= 0

其中 D2 是二阶差分矩阵，用于惩罚孔径分布曲线过度振荡；lambda 是正则化强度。lambda 越小，结果越贴近实验点但可能噪声更大；lambda 越大，曲线越平滑但可能掩盖真实窄峰。

当前软件默认正则化参数为 0.316，并提供从 0 到 10 的候选值用于诊断。软件会输出模型拟合吸附量、残差 RMS 和分布粗糙度，帮助判断 DFT 结果是否可信。

如果存在官方 DFT/NLDFT 模型核文件，软件优先使用模型文件；如果模型文件缺失，则使用平滑解析核作为备用，以保证界面可以继续运行。严格发表时，应优先使用与材料类型、孔几何和吸附质相匹配的官方或文献模型核。

八、为什么统一软件必须保留“厂家差异”？

很多人会认为统一软件就是把所有数据都套同一个公式。实际上这会带来新的问题：不同厂家软件的默认结果之所以不同，不只是界面不同，而是背后存在具体算法差异。

这些差异包括：

1. BET 默认选点不同。 有的软件保存官方点号，有的软件按压力范围，有的软件根据 Rouquerol 判据截断。

2. 物理常数不同。 例如氮气分子截面积、阿伏伽德罗常数、液体体积换算因子、Kelvin 因子都可能存在历史版本差异。

3. t-Plot 纵轴定义不同。 有的软件用液体体积，有的软件报表中以 STP 吸附量为纵轴，再换算外比表面积。

4. BJH 输出来源不同。 有的软件可以导出逐点官方 BJH 表，有的软件只能从等温线重算；有的软件对低孔径区间、平滑和修正方式有专门处理。

5. DFT 模型核不同。 DFT/NLDFT 结果强依赖模型核，slit、cylinder、carbon、oxide、77 K N2、87 K Ar 等都不能混用。

因此，这套软件的设计不是“强行消除差异”，而是分成两层：第一层尽量复现来源软件的默认规则，第二层允许用户用统一规则重新计算。这样既能核对原报告，也能进行横向比较。

九、总结：专业 BET 软件不应该只给一个数字

一个专业的 BET 分析软件，不应该只输出“比表面积 = 多少 m²/g”。它应该告诉用户：用了哪些原始点，采用了哪个压力区间，C 常数是否合理，是否满足 Rouquerol 型判据，厚度方程是什么，孔径分布来自官方表还是重新计算，是否经过平滑，DFT 使用了哪个模型核和正则化参数。

这也是我开发这套 BET 综合分析软件的核心目标：让每一个结果都能回到数据、公式、参数和软件版本。只有这样，不同样品之间、不同仪器之间、不同厂商软件之间的比较才真正有意义。

参考文献

[1] Brunauer, S.; Emmett, P. H.; Teller, E. Adsorption of Gases in Multimolecular Layers. https://doi.org/10.1021/ja01269a023

[2] Langmuir, I. The Adsorption of Gases on Plane Surfaces of Glass, Mica and Platinum. https://doi.org/10.1021/ja02242a004

[3] Thommes, M. et al. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution. https://doi.org/10.1515/pac-2014-1117

[4] Barrett, E. P.; Joyner, L. G.; Halenda, P. P. The Determination of Pore Volume and Area Distributions in Porous Substances. https://doi.org/10.1021/ja01145a126

[5] Lippens, B. C.; de Boer, J. H. Studies on pore systems in catalysts. https://doi.org/10.1016/0021-9517(65)90307-6

[6] Dollimore, D.; Heal, G. R. An Improved Method for the Calculation of Pore Size Distribution from Adsorption Data. https://doi.org/10.1002/jctb.5010140302

[7] Horvath, G.; Kawazoe, K. Method for the Calculation of Effective Pore Size Distribution in Molecular Sieve Carbon. https://doi.org/10.1252/jcej.16.470

[8] Lastoskie, C.; Gubbins, K. E.; Quirke, N. Pore size distribution analysis of microporous carbons: a density functional theory approach. https://doi.org/10.1021/j100120a035