4. 实验数据预处理与质量检查
Langmuir拟合的可靠性并不取决于拟合软件本身,而首先取决于进入拟合程序的那一组数据是否干净、可比、基准一致。同一块煤样,若原始等温线数据未做干燥无灰基校正、压力单位混用、或把解吸支与吸附支数据混在一起,即使拟合相关系数R²高达0.999,得到的Langmuir体积V_L和Langmuir压力P_L也缺乏物理意义。本节围绕原始数据整理、异常值识别与剔除、单位统一与换算、吸附量基准校正四个环节,说明如何把实验室记录的原始数据整理成可直接用于Langmuir拟合的高质量数据集。
一、原始数据整理:先建立可追溯的数据表
煤岩吸附等温线实验通常采用容量法( volumetric )或重量法,在多个压力平衡点记录吸附量。原始记录往往分散在仪器导出文件、天平读数、温压传感器日志和人工记录本中。整理的第一步不是计算,而是建立一张字段完整、一行对应一个平衡点的数据表。建议至少包含以下字段:样品编号、煤阶/采样位置、实验温度、平衡压力(含单位)、吸附量原始读数(含基准说明)、样品质量、水分Mad、灰分Aad、实验日期、数据支别(吸附支/解吸支)、备注。
其中“数据支别”常被忽略。吸附等温线在升压段测得吸附支,降压段测得解吸支,两者在相同压力下吸附量往往不重合,形成滞后环。Langmuir拟合通常针对吸附支进行,若两支数据混排,拟合结果会被系统性拉偏。整理时应把两支数据分行标注,或至少用独立列区分,避免后续误用。
二、异常值识别与剔除:用物理约束和统计判据双重筛查
异常值来源主要有三类:仪器波动(压力传感器漂移、天平零点漂移)、操作失误(平衡时间不足、阀门切换错误)、记录错误(单位写错、小数点错位)。识别时不能只看数值大小,应结合物理约束判断。例如:吸附量随压力单调不减(吸附支),若出现明显回落,则该点可疑;低压段吸附量应趋近于零,若首点吸附量异常偏高,可能是死体积校正不当;同一压力下重复测量偏差超过设定阈值(如5%),应标记复核。
统计判据可作为辅助。对同一压力附近的重复点计算相对偏差,或对整条曲线做残差检查:先用粗略Langmuir模型拟合全部点,计算各点残差,残差绝对值超过2倍标准差的点列为候选异常点。但剔除必须谨慎——煤样在高压段吸附量增长趋缓是正常现象,不能因“看起来平”就剔除。剔除记录应保留原始值、剔除理由和处理方式,保证数据可追溯。一般建议单条等温线剔除点不超过总点数的10%,否则应重新检查实验过程而非简单删点。
三、单位统一与换算:压力、吸附量、温度三线对齐
Langmuir方程中P_L与压力同单位,V_L与吸附量同单位,单位不统一会直接导致参数数量级错误。常见压力单位有MPa、kPa、bar、atm,换算关系为1 MPa=1000 kPa,1 bar=0.1 MPa,1 atm≈0.101325 MPa。吸附量单位常见cm³/g、mL/g、mmol/g、m³/t。若采用mmol/g,需用摩尔体积22.414 L/mol(标准状态)换算:1 cm³/g≈0.0446 mmol/g。温度统一为K或℃,但Langmuir拟合中温度主要作为实验条件记录,不进入方程本身,需注意同一组拟合数据应来自同一温度。
换算时应坚持“先统一、后计算”的原则:在数据表中新增“压力_MPa”“吸附量_cm³·g⁻¹”等标准列,原始列保留不动,所有换算公式写在表头注释或单独说明中。这样既避免反复换算引入误差,也便于复核。
四、吸附量基准校正:从原始读数到干燥无灰基
煤样吸附量通常以“每克样品”为基准记录,但样品中含有水分和灰分,二者不参与甲烷吸附,却计入样品质量,导致吸附量被稀释。Langmuir拟合要求吸附量以干燥无灰基(daf)表示,才能在不同煤样之间比较。校正分两步:先扣除水分和灰分质量,得到干燥无灰基质量;再将原始吸附量除以该质量。
设原始样品质量为m(g),空气干燥基水分为Mad(%),灰分为Aad(%),则干燥无灰基质量m_daf=m×(1−Mad/100−Aad/100)。若原始吸附量V_raw以cm³/g(原样)记录,则V_daf=V_raw×m/m_daf=V_raw/(1−Mad/100−Aad/100)。例如某样品m=50.00 g,Mad=1.20%,Aad=8.50%,则m_daf=50.00×(1−0.012−0.085)=45.15 g,校正系数为50.00/45.15≈1.107,即干燥无灰基吸附量约为原样基准的1.107倍。这一系数看似不大,但当灰分较高(如20%以上)时,校正系数可达1.25以上,对V_L的影响不可忽略。
五、形成可复核的数据处理记录
完成上述四步后,应形成一份简短的数据处理记录,至少包含:原始数据文件名与来源、剔除点清单及理由、单位换算公式与系数、Mad/Aad取值来源、校正系数计算结果、最终拟合数据点数。这份记录不必复杂,但必须能让另一位研究者从原始文件复现出同一组拟合输入。实践中常见的问题是校正系数只算了一次却用于不同样品,或单位换算在表格中手工完成导致个别行遗漏,这些都可以通过“标准列+记录表”的方式避免。
需要强调的是,预处理不是对数据的“美化”。剔除异常点必须有物理或操作依据,校正必须基于实测的工业分析数据,单位换算必须双向可逆。经过这样处理后的数据集,才适合进入Langmuir拟合环节,得到的V_L和P_L才具备与煤阶、煤岩组成对比的物理意义。