跳至正文

EIS/KK 一致性检验

采用 Lin-KK 算法检验 EIS 实部与虚部的一致性,为后续拟合与 DRT 分析提供数据质量参考

EIS/KK 一致性检验

Kramers-Kronig 关系用于检验 EIS 实部与虚部的一致性。偏离可能与非线性、非稳态、有限频段或测量误差有关,但该检验不能单独识别具体原因。

KK 检验适用于等效电路拟合、DRT 分析或样品比较前的数据质量评估。其结果不用于选择等效电路,也不足以验证材料机理。

适用场景

以下情形建议先进行 KK 检验:

  • 准备进行等效电路拟合,希望先确认数据质量
  • Nyquist 或 Bode 图中有少量可疑频率点,但不确定是否应剔除
  • 低频端可能存在漂移、样品状态变化或测试时间过长
  • 同一批样品需要比较,想先排除明显不稳定的数据
  • 拟合参数不稳定,想确认问题来自数据质量还是模型选择

运行前应使用 EIS 绘图:Nyquist 和 Bode 检查原始 Nyquist/Bode 曲线。

输入数据

输入可为仪器导出的 EIS 数据文件夹或多个原始文件。支持常见文本、CSV、Excel、EC-Lab .mpr、Gamry .dta 和 VersaStudio .par 等格式。

如果已经有标准 EIS CSV,也可以直接输入。标准数据至少需要包含:

列名含义
freq_hz频率
z_real_ohm阻抗实部
z_imag_ohm阻抗虚部

分析流程

选择数据后,界面显示所有样品的原始 Nyquist 图。图表支持悬停查看频率与阻抗值、缩放、平移和双击复位,横纵坐标采用相同比例,以避免谱形因坐标缩放而失真。

程序根据各样品最低频端的谱形提供低频边界建议,并列出最低频率、低频末段 -Z'' 的起止值及对应建议。若所有样品的建议一致,界面预选对应模式;若样品间不一致,界面提示按表中建议分组计算。自动建议仅用于辅助设置,最终选择仍应依据原始 Nyquist/Bode 谱形。

程序仅读取原始数据,分析结果写入独立输出目录。

结果解释应同时考虑报告分级、残差的频率分布及原始阻抗谱,不应仅依据单一汇总指标。

方法与低频边界设置

这里采用 Lin-KK 算法 [1],以 ZZ'(实部)为拟合目标,并以按 Z|Z| 归一化的 ZZ'' 相对残差作为主要筛查指标。

低频边界有两种模式:

  • 常规谱,有低频截距:最低频端已达到或接近实轴截距,不加入串联电容项。
  • 无低频截距/低频电容尾:最低频端仍呈持续延伸的电容性尾部,加入串联电容项。

串联电容项会增加辅助模型自由度,因此程序不根据残差大小自动选择低频边界模型。用户选择记录于 kk_summary.csvboundary_modeadd_cap 两列。

输出结果

结果保存于输入目录下的 kk_validation_output 文件夹。

每个成功分析的样品会输出:

  • {sample}_kk_fit.csv:实验阻抗、Lin-KK 拟合阻抗和实部/虚部相对残差。
  • {sample}_kk_fit.png:Nyquist 对比图和残差图。

汇总文件包括:

  • kk_summary.csv:所有样品的点数、MmuRMSE(Z'')RMSE(Z')、最大虚部残差、低频边界模型和结论。

点击「确认生成」后,还可以额外导出:

  • kk_validation.opju:Origin 工程,包含 Nyquist 对比和残差图。

怎么看结果

报告中的结论分为三类:

结论含义
低残差虚部相对残差低于工程阈值;该分级不构成数据有效性的充分证据
需复核虚部相对残差超过工程阈值;需检查是否存在频段相关的系统偏离
无法判断当前结果不足以支持分级,可能与有效频点不足、频率范围有限或数值异常有关

程序以虚部相对残差 RMSE(Z'') 为主要筛查指标,实部残差仅作参考。当 RMSE(Z'') ≤ 1% 且最大单点虚部相对残差 ≤ 5% 时,结果标记为低残差。上述阈值属于产品工程标准,并非文献规定的合格判据。残差是否随频率在零点附近随机分布,仍是解释 KK 一致性的重要依据。

结果解释还应检查 {sample}_kk_fit.png

  • Nyquist 对比图中,Lin-KK 拟合曲线应大体跟随实验点。
  • 残差图中,残差最好围绕 0 随机波动,而不是在某个频率区间持续偏正或偏负。
  • 如果低频端残差明显偏离,常见原因包括样品漂移、测试时间过长、非线性响应或低频点不稳定。
  • 如果只有个别频率点偏离,建议结合原始实验记录判断是否存在接触扰动、仪器切档或数据导出异常。

方法说明

Lin-KK 算法使用一组按时间常数分布的 RC 元件拟合阻抗谱,并比较拟合值与实测值 [1]。KK 关系见 [2, 3],EIS 数据分析方法见 [4]。

KK 残差偏高可能与非线性、非稳态、有限频段或测量误差有关。该检验不能单独识别偏离来源,也不能验证等效电路或反应机理的正确性。

后续处理

  • 低残差:结合残差图和实验记录完成复核后,再进行 EIS/等效电路拟合EIS/DRT 分析 分析。
  • 需复核:使用 EIS 绘图:Nyquist 和 Bode 检查低频漂移、离群点和测试稳定性。
  • 无法判断:检查数据格式、频率范围和有效点数;当前结果不应直接解释为材料行为。

参考文献

  1. Schönleber, M., Klotz, D., and Ivers-Tiffée, E. (2014). A method for improving the robustness of linear Kramers-Kronig validity tests. Electrochim. Acta 131, 20-27. https://doi.org/10.1016/j.electacta.2014.01.034.
  2. Boukamp, B.A. (1995). A linear Kronig-Kramers transform test for immittance data validation. J. Electrochem. Soc. 142, 1885-1894. https://doi.org/10.1149/1.2044210.
  3. Agarwal, P., Orazem, M.E., and Garcia-Rubio, L.H. (1992). Measurement models for electrochemical impedance spectroscopy. I. Demonstration of applicability. J. Electrochem. Soc. 139, 1917-1927. https://doi.org/10.1149/1.2069522.
  4. Orazem, M.E., and Tribollet, B. (2017). Electrochemical Impedance Spectroscopy, 2nd ed. (John Wiley & Sons). https://doi.org/10.1002/9781119363682.