优化玻璃pH电极在pH范围?4.5至15.5的非线性校准：使用改进的酸误差和碱误差统计处理方法

《Talanta》：Optimising Nonlinear Calibration of Glass pH Electrodes over the pH Range ?4.5 to 15.5 Using Improved Statistical Treatment of Acid and Alkaline Errors

　　

准确校准广泛使用的玻璃pH电极在高（pH < 1）和低（pH > 11）质子活度（a_H+）条件下尤其具有挑战性，此时电池电势与log(a_H+)之间的线性Nernst关系失效。通过校正赤池信息准则（AICc）和均方根误差（RMSE）从七个非线性函数候选中选出的非对称双逻辑函数，最佳地描述了从pH ?4.5至15.5、温度25°C时，涵盖酸误差和碱误差区域的电动势（EMF）–pH关系（R² = 1.00，RMSE = 2.3 mV）。参考pH值通过质子活度（Pitzer模型，MacInnes约定）赋予五个NIST/DIN缓冲液、十个H₂SO₄标准液（0.10–7.54 mol·L^-1 H₂SO₄；pH ?4.5至1）和九个NaOH标准液（0.10–6.91 mol·L^-1 NaOH；pH 12.9至15.5）。与传统的三点校准相比，该函数在pH ?4.5处将酸误差偏差降低最多3.0个pH单位，在pH 15.5处将碱误差偏差降低最多1.5个pH单位。全子集和留一法蒙特卡洛分析进一步表明，五个酸性标准液和三个碱性标准液，连同NIST/DIN可追溯缓冲液，在极端条件下将全数据集精度保持在±0.25 pH单位以内。该非线性校准模型可在近二十个pH单位内实现连续校准；当仅关注一个区域时，该模型简化为四参数逻辑函数。缩减的标准集在保持精度的同时，将校准时间至少减少50%。
新颖性声明：这是首次在单个连续非线性函数中描述玻璃pH电极的酸误差和碱误差，将可靠测量范围从pH ?4.5扩展至15.5。统计分析进一步确定了所需极端pH标准液的最小数量，将先前单区域方法的校准工作量大约减半。

**研究背景与问题**
玻璃pH电极是测量pH的标准传感器，其精度依赖于Nernst响应。然而，在强酸性（pH < 1）和强碱性（pH > 11）条件下，电极响应偏离线性，产生酸误差和碱误差。传统三点校准（如pH 4.01、7.01、10.01）仅适用于中间范围，在极端pH下会系统性地高估或低估pH。许多自然和工程系统（如酸性矿山排水、火山湖泊、水泥孔隙液、湿法冶金流程）的pH范围远超常规缓冲液范围。现有非线性校准方法（如Nordstrom等和Licht等）仅覆盖单一极端区域，且使用大量校准溶液，耗时且可能损害电极。因此，研究人员旨在开发一个连续非线性函数，同时描述酸误差和碱误差，并优化校准标准数量，以扩展可靠测量范围并提高效率。该论文发表在《Talanta》。

**研究内容与结论**
研究人员通过制备10个H₂SO₄标准液（pH ?4.5至1）和9个NaOH标准液（pH 12.9至15.5），结合5个NIST/DIN缓冲液（pH 4.01至13.00），在25°C下测量电动势（EMF），并利用Pitzer模型和MacInnes约定计算质子活度以分配参考pH。通过比较7个候选非线性模型（包括多项式、逻辑函数、Richards模型、Gompertz模型、双逻辑函数），使用校正赤池信息准则（AICc）和均方根误差（RMSE）进行选择。双逻辑函数（R2 = 1.00，RMSE = 2.3 mV）以最低AICc（60.7）和最小RMSE胜出，能同时描述酸误差和碱误差，将传统三点校准的偏差（pH ?4.5处+3.0 pH单位，pH 15.5处?1.5 pH单位）降至±0.25 pH单位以内。通过全子集分析、指纹组合分析和留一法蒙特卡洛分析，确定最小校准标准数量：酸性区域5个H₂SO₄标准液（含最极端pH ?4.5），碱性区域3个NaOH标准液（含最极端pH 15.5），在保持精度（ΔpH ≤ 0.25）的同时将校准时间减少50%以上。该模型可应用于跨极端pH的系统（如火山湖泊、酸性矿山排水），或简化至单区域四参数逻辑函数。

**关键技术方法**
（1）**质子活度计算**：采用Pitzer模型及MacInnes约定，通过PHREEQC软件计算H₂SO₄和NaOH标准液（浓度由电位滴定验证）的质子活度，分配参考pH值。
（2）**电动势测量**：在25°C下使用Metrohm Unitrode电极测量各标准液的EMF，稳定性阈值<0.5 mV·min^-1。
（3）**非线性模型选择**：通过AICc和RMSE比较7个候选函数（包括3次多项式、4参数逻辑函数、5参数逻辑函数、Richards模型、Gompertz模型、双逻辑函数），选择最优模型。
（4）**校准设计优化**：通过全子集分析、指纹组合分析和留一法蒙特卡洛分析，确定最小标准集数量及组成，利用局部斜率|dEMF/dpH|将pH误差容忍度（ΔpH ≤ 0.25或0.50）转换为RMSE阈值。

**研究结果**
**3.1 非理想溶液中的质子活度（pH < 0 和 pH > 14）**
通过Pitzer模型计算，H₂SO₄溶液中质子活度系数γ_H+从0.77增至3.15×103，pH降至?4.5；NaOH溶液中羟基活度系数γ_OH-从0.78增至2.56，pH升至15.5。EMF测量显示，酸误差校正因子（CF）在pH ?4.5处达175.5 mV，碱误差CF在pH 15.5处达?82.3 mV，表明线性Nernst模型在极端pH下显著偏离。

**3.2 pH ?4.5至15.5的非线性校准模型**
比较7个候选模型，双逻辑函数获得最低AICc（60.7）和最低平均RMSE（2.3 mV，酸性区域1.7 mV，碱性区域2.8 mV），而其他模型RMSE为4.8–7.4 mV。双逻辑函数能对称处理两个误差区域，而单S形模型（如Gompertz）在酸性或碱性区域表现不对称。

**3.3 双逻辑函数在非对称pH校准中的性能**
双逻辑函数拟合24个数据点，得到公式：EMF = 542.9×[1+exp(0.4×(pH?3.1))]^?1 ? 411×[1+exp(0.5×(pH?10.7))]^?1。中间pH 1–11区域Nernst斜率为?58.1 mV·pH^-1（理论的98.2%）。酸误差与碱误差分别对应两个逻辑项，其拐点pH 3.1和10.7反映了质子饱和与Na?干扰的物理化学机制。该模型可简化为单一4参数逻辑函数用于单区域应用。

**3.4 玻璃pH电极的酸误差和碱误差评估**
与文献数据比较，本研究的酸误差（pH ?4.5处约+3 pH单位）大于早期研究（如Corning 015电极），归因于更宽的负pH范围。碱误差（pH 15.5处约?1.5 pH单位）与Licht等报告的总碱误差一致，但Licht仅分离了K?干扰贡献，而本模型包括总偏差。

**3.5 扩展非线性pH范围的最小校准标准数量与选择**
通过全子集分析，酸性区域中5个标准的中位数RMSE为2.00 mV（接近不可约下限1.81 mV），严格容忍度（ΔpH ≤ 0.25）通过率59.9%；碱性区域中3个标准通过率91.3%。指纹组合分析表明，酸性区域最优组合包含最极端pH（?4.5）和拐点区域（pH ?2至?1）；碱性区域必然包含最极端pH 15.5。留一法分析显示，最极端标准（pH ?4.5和15.5）对参数约束最大，不可替代。推荐最小设计：酸性5个H₂SO₄标准（含pH ?4.5、?2至?1、+1）加缓冲液；碱性3个NaOH标准（含pH 15.5）加缓冲液，分别减少校准工作量50%和67%。

**总结与讨论**
该研究首次用单一连续双逻辑函数描述玻璃pH电极在pH ?4.5至15.5的响应，将酸误差和碱误差视为Nernst行为的连续偏离。这一连续性使跨极端pH系统（如火山湖泊、酸性矿山排水、湿法冶金流）可在同一标尺下测量，无需切换程序。统计优化将该精度保持仅需5个酸性标准和3个碱性标准，缩短校准时间并减少电极暴露于极端条件。讨论部分指出，该模型在pH ?4.5处的RMSE（1.7 mV）对应于±0.22 pH单位的不确定性，在pH 15.5处为±0.18 pH单位，低于严格容忍度。与传统方法相比，该方法避免了分别处理酸误差和碱误差的不一致性，并显著降低了校准成本。

**研究结论**
（翻译原文结论部分）一个单一的双逻辑函数描述了玻璃电极从pH ?4.5至15.5的完整响应，将其解构为两个S形转变，分别对应两种误差的不同物理化学起源：强酸性溶液中膜质子饱和与液接电势升高，以及强碱性溶液中膜表面的Na?干扰。将酸误差和碱误差视为Nernst行为的连续偏离，而非两个单独校正的区域，使得单一校准在近二十个pH单位内保持有效，同时在中性区域保持近Nernst行为。这种连续性是一项实际进步。跨越两个极端条件的系统（如火山湖泊、正在修复的酸性矿山排水、湿法冶金流）可以在同一标尺上测量，无需切换程序；仅关注一个极端时，该函数简化为一个逻辑项。对校准设计的统计评估进一步表明，这种精度仅需五个酸性标准液和三个碱性标准液，约为本研究中测量数量的一半。制备和测量更少的极端pH溶液可缩短整体校准时间，同样重要的是，限制了电极暴露于逐渐降解玻璃膜的浓酸和浓碱。