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

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],以 (实部)为拟合目标,并以按 归一化的 相对残差作为主要筛查指标。
低频边界有两种模式:
- 常规谱,有低频截距:最低频端已达到或接近实轴截距,不加入串联电容项。
- 无低频截距/低频电容尾:最低频端仍呈持续延伸的电容性尾部,加入串联电容项。
串联电容项会增加辅助模型自由度,因此程序不根据残差大小自动选择低频边界模型。用户选择记录于 kk_summary.csv 的 boundary_mode 和 add_cap 两列。
输出结果
结果保存于输入目录下的 kk_validation_output 文件夹。
每个成功分析的样品会输出:
{sample}_kk_fit.csv:实验阻抗、Lin-KK 拟合阻抗和实部/虚部相对残差。{sample}_kk_fit.png:Nyquist 对比图和残差图。
汇总文件包括:
kk_summary.csv:所有样品的点数、M、mu、RMSE(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 检查低频漂移、离群点和测试稳定性。
- 无法判断:检查数据格式、频率范围和有效点数;当前结果不应直接解释为材料行为。
参考文献
- 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.
- 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.
- 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.
- Orazem, M.E., and Tribollet, B. (2017). Electrochemical Impedance Spectroscopy, 2nd ed. (John Wiley & Sons). https://doi.org/10.1002/9781119363682.