说明：本文华算科技主要介绍 EIS 等效电路拟合得到的电阻、电容和弛豫时间参数为什么经常与真实电化学过程偏离，包括拟合在求解什么、同一条阻抗谱为什么能对应多套参数、测量数据中的非稳态与非线性如何让结果偏移、CPE 与弥散效应为什么让界面电容难以还原，以及哪些做法能让拟合参数贴近界面真实过程。
电化学阻抗谱的测量，是向电极施加一个小幅值的交流电压或电流扰动，并在从高到低的一系列频率上记录体系的响应。每个频率点都得到一个复数阻抗，其中实部代表能量耗散，虚部代表能量存储。
把所有频率点的复数阻抗排列起来，就构成 Nyquist 图与 Bode 图两种常见谱图。这两类谱图描述的是同一组阻抗数据，只是把实部、虚部和频率的关系用不同方式画了出来。
在这些谱图上，高频段容纳溶液电阻和界面电容这类响应很快的过程，低频段容纳电荷转移、吸附与扩散这类响应较慢的过程。
要从一条谱线里提取物理量，常规做法是把它拟合到一个等效电路模型上，最常见的是由溶液电阻、电荷转移电阻和双电层电容组成的Randles 电路，扩散占主导时再串入一个 Warburg元件。整个拟合的目标，就是为这些电路元件找到一组能复现实测谱线的数值。
图1. 碳钢在 1.0 M HCl 中加入与未加缓蚀剂时的 Nyquist 阻抗谱。DOI：10.1038/s41598-025-12838-2
在理想的Nyquist 图上，高频与实轴的交点给出溶液电阻，半圆的直径约等于电荷转移电阻，半圆顶点的频率与电荷转移电阻和电容的乘积相关，由此还能估出双电层电容的量级。这几个几何特征，构成了人们对等效电路参数的初步估计，也是拟合算法迭代时所需要的初值来源。
图2. 拟合碳钢/盐酸界面阻抗所用的等效电路模型。DOI：10.1038/s41598-025-12838-2
拟合本身依靠复数非线性最小二乘法，算法不断调整每个电路元件的取值，使模型阻抗与实测阻抗在所有频率点上的加权残差最小，从而求解 Rs、Rct、Q、n 等参数。
拟合得到的 Rct 和 Cdl 等参数是否具有物理意义，取决于电路模型、参数范围和数据质量，任何一处偏差都会让这些数值与真实的电荷转移和双电层状态产生错位，谱线吻合得再好也无济于事。
同一段被压扁的半圆，可以用好几种等效电路来描述。电阻并联电容、电阻并联 CPE、两个 RC 串联，甚至嵌套的梯形网络，都能生成形状几乎一致的阻抗曲线。
把这些电路分别拿去拟合同一组数据，软件给出的残差都很小，χ² 也都处在可接受的范围。仅凭谱线的吻合程度，几乎没办法从这几种电路里挑出唯一正确的那一个。
图3. 拟合碳钢/3.5% NaCl 体系阻抗数据的两种等效电路形式。DOI：10.1038/s41598-022-24793-3
选哪一种电路形式，取决于研究者对界面结构的物理理解，而拟合软件只优化数值，对电路本身是否合理不作评估。用单个 RC去套多孔电极，或给简单界面强行加上两个时间常数，软件照样会返回一组数学上成立、物理上空洞的参数。一旦电路形式选错，后面所有元件值都失去了原本的物理含义。
即便电路形式固定，参数之间仍然存在很强的相关。半圆的位置主要由溶液电阻决定，直径由电荷转移电阻决定，顶点频率由电荷转移电阻和双电层电容一起决定；当某个参数偏大时，算法可以用另一个参数来补偿，仍然得到接近的谱线。
CPE 的 Q 与指数 n 之间、Warburg 系数与低频电阻之间，都存在类似的参数耦合。
图4. NiCo2O4 修饰电极在碱性体系中的 Nyquist 阻抗谱。DOI：10.1038/s41598-023-49692-z
当电路元件的数目超过谱线能分辨的时间常数个数时，目标函数会出现多个相近的极小值，这就是常说的过参数化。
这时参数的置信区间变宽，相关系数接近 1，χ² 还在继续下降，元件值却随初值大幅漂移。拟合结果于是从一个稳定的解，变成一片彼此等价、却各自对应不同物理图像的解。
阻抗谱的可靠拟合，建立在线性、稳态、因果三个前提之上。扰动幅值通常控制在5 到 10 毫伏，一旦幅值过大，电流响应就会偏离线性。这时谱线里会混入谐波失真，使复数非线性最小二乘赖以成立的线性假设不再满足，拟合自然失去了根基。
图5. 碳钢在 3.5% NaCl 中不同缓蚀剂浓度下的 Nyquist 阻抗谱。DOI：10.1038/s41598-022-24793-3
稳态要求体系在整个测量期间保持稳定，但低频扫描耗时数十秒到数分钟，期间若发生电极腐蚀、SOC 变化或温度波动，低频数据就采自不同的体系状态。
Kramers–Kronig变换用来检验数据的自洽性，实部与虚部的 KK 关系是线性稳态因果体系的必然结果。
当残差沿频率呈现结构化趋势而非随机散布时，通常对应低频非稳态或高频干扰，此时再合理的电路也会把这些偏差吸收成虚假的元件参数。
频率窗口的截断会让某些参数欠定。低频扫描若没有覆盖扩散尾部，Warburg系数和低频电阻就缺乏足够约束；每十倍频程的采样点过少时，半圆和扩散段的特征点不足以支撑可靠拟合，算法对噪声的敏感度随之上升。测量频率范围本身就决定了哪些参数能被稳定求解，窗口之外的过程无法被这一次测量约束。
图6. 碳钢/盐酸体系的 Bode 模值谱与相位谱。DOI：10.1038/s41598-025-12838-2
高频端的问题主要来自硬件。电极引线、参比电极内阻和电缆之间的寄生电感会在高频产生负虚部，在 Nyquist 图上表现为实轴下方的尾部，与电化学界面过程无关。
如果不排除这些点，拟合会错误地把电感效应分配到溶液电阻或双电层电容上，使高频元件值整体偏移，连带影响串联在其后的中低频参数。
恒定相位元件CPE 的阻抗写作 Z=1/[Q(jω)n]，指数 n 偏离 1 的程度刻画了界面弥散的强弱。这种弥散来自表面粗糙、电流分布不均，以及双电层弛豫时间的分散，几乎所有真实电极都或多或少带有这种非理想性。
当 n 等于 1 时 CPE 退化为理想电容，n 等于 0.5 时接近 Warburg阻抗，介于两者之间时，Q 的量纲已经偏离法拉，不再是一个直接可读的电容值。
图7. 碳钢腐蚀体系相位角随频率的变化谱。DOI：10.1038/s41598-022-24793-3
同一个n 值在不同体系里来源各不相同，有的源于几何粗糙，有的源于多孔结构内的电流深度分布，还有的源于电极表面活性位点的不均匀。单看 Q 和 n 两个数字，很难还原界面的真实电容，也很难分辨弥散究竟来自哪一种物理过程。
把CPE 的 Q 换算成有效电容，需要借助Brug公式或 Hsu–Mansfeld公式，这是把恒定相位元件折算为等效电容的两种常见做法。两者都用到电阻和指数 n，配合体系的几何假设，算出一个等效的电容数值。
图8. 石墨烯电极/电解液界面的电容组成与电子态占据。DOI：10.1038/s41598-025-04357-x
这两套公式基于不同的弥散机制假设，对同一组参数会给出不同的结果，而且一旦电阻拟合出现偏差，有效电容的计算也跟着偏移。在多孔电极中，弥散来自孔道内的电位分布，要用传输线模型才能描述其阻抗特征。
在半导体或碳电极上，界面处还会叠加空间电荷电容和量子电容，三者串联后测得的总电容由最小的那一项主导，这样把拟合的 Q 等同于双电层电容，就会高估或低估真实的界面电容，据此算出的活性面积和缺陷密度也跟着一起失准。
测谱之前的稳态控制，是整个拟合可信的基础。需要在开路电位稳定之后再开始测量，对电池类体系固定 SOC 和温度，腐蚀体系则浸泡到接近平衡的状态。
扰动幅值要压在线性区内，通常取5 到 10 毫伏，低频扫描尽量深入到 10 毫赫兹以下来覆盖扩散过程，每十倍频程采集 6 到 10 个数据点，保证半圆和扩散段都有足够多的取样。
图9. NiCo2O4 修饰电极在析氧电位下的阻抗谱。DOI：10.1038/s41598-023-49692-z
测完之后用Kramers–Kronig 检验数据是否自洽，残差随机散布对应体系在测量过程中维持稳态。在相同条件下重复测量，让两条谱线基本重合，再把高频寄生电感产生的伪影从数据里剔除。经过这样一轮把关，送入拟合的就是一组干净、可重复、满足因果性的阻抗数据。
电路形式应当由体系的物理过程决定，而非凭经验套用固定模板。可以通过分布弛豫时间分析判断需要多少个 RC 或 RQ 单元，再为每个参数设定合理的初值和取值范围。
拟合之后检查置信区间，如果区间过宽或相关系数接近 1，说明参数冗余了。判断拟合质量不只看 χ² 的数量级，10-4 以下算较好，更要看残差是否随机、元件值是否处在物理上合理的取值范围。
把同一条件下测得的多条谱一起拟合，让元件值随实验条件单调变化、且与电化学动力学一致，参数才更有参考价值。举例来说，随着过电位升高，电荷转移电阻应当单调减小，对应电荷转移过程加快；这种参数随条件的变化趋势，比单个数据点的绝对值更有说服力，也更能抵御噪声和个别异常点对结果带来的干扰。
