[发明专利]双柱串联气相色谱保留时间的预测方法

摘要

本发明涉及一种双柱串联气相色谱保留时间的预测方法，属气相色谱技术领域。该方法过程包括：首先确定组分在独立柱1、柱2任意温度点容量因子；其次利用甲醇标准物校正柱1和柱2串联柱的任意温度点死时间；然后确定组分在串联柱中程序升温段任意时刻容量因子；最后利用上述数据对组分在串联柱上程序升温条件下的保留时间进行预测，并与实测值进行比较。本发明优点：减少实验次数，降低工作量，节约分析时间，预测精度高。此外，本发明为双柱串联气相色谱最佳分离条件的选择提供可行途径。

主权项

1.一种双柱串联气相色谱保留时间的预测方法，所述的双柱串联，是在HP 6890气相色谱仪上，以一个非极性色谱柱和一个强极性色谱柱通过死体积可忽略不计的接头直接连接构成，其特征在于包括以下步骤：第一步确定待测组分在独立柱1、柱2任意温度点的容量因子：1)确定待测组分在独立柱1任意温度点的容量因子：(1)柱1在温度30-250℃范围内设定6个温度点：T₁＝30℃、T₂＝50℃、T₃＝100℃、T₄＝150℃、T₅＝200℃和T₆＝250℃，分别测定待测组分在这6个温度点下的保留时间；(2)测定甲烷在上述6个温度点下的保留时间作为柱1各个温度点的死时间；(3)将步骤(1)确定的保留时间和步骤(2)确定的死时间分别代入式1中，得到待测组分在6个温度点下相应的各个容量因子：k′＝t_R/t_M-1             式1式1中：t_R为待测组分在对应温度点的保留时间，t_M为对应温度点的死时间，k′为待测组分的在对应温度点的容量因子；(4)将由式1计算得到的6个温度点下的各自容量因子与相应的温度值分别代入式2中，得到包含6个方程的方程组，通过解方程组，确定柱1未知系数a、b、c和d的值：ln k′＝aT³+bT²+cT+d     式2(5)将确定后的a、b、c和d再代入式2中，得到待测组分在柱1上任意温度点的容量因子k′；2)确定待测组分在独立柱2任意温度点的容量因子：(1)柱2在温度30-220℃范围内设定6个温度点：T₁’＝30℃、T₂’＝50℃、T₃’＝100℃、T₄’＝150℃、T₅’＝200℃和T₆’＝220℃；(2)上述6个温度点对应的保留时间、死时间、容量因子及待测组分在柱2上任意温度点的容量因子确定过程与步骤1)相同；第二步利用甲醇标准物校正柱1和柱2串联后的柱1、柱2任意温度点的死时间：1)以甲醇标准物在T₁’＝30℃、T₂’＝50℃、T₃’＝100℃、T₄’＝150℃、T₅’＝200℃和T₆’＝220℃6个温度点校正串联后柱1、柱2的死时间：(1)在串联柱上测定甲烷在上述6个温度点下的保留时间，并作为串联柱的死时间；(2)上述串联柱的死时间，是串联后柱1和柱2在相同温度点下的死时间之和：t_M＝t_M1+t_M2              式3式3中：t_M为串联柱的死时间，经步骤(1)已测得，t_M1为串联后柱1的死时间，为未知，t_M2为串联后柱2的死时间，为未知；(3)测定串联柱上甲醇标准物在上述6个温度点下的保留时间t_R；(4)计算串联柱上甲醇标准物在上述6个温度点下的校正保留时间t_P：以t_M1＝0，t_M2＝t_M为起始条件，计算上述6个温度点甲醇校正保留时间，计算步骤如下：①计算甲醇在串联后柱1上6个温度点的校正保留时间t_P1：t_P1＝t_M1(1+k₁’)       式4式4中：k₁’是根据第一步中1)计算出甲醇在独立柱1上每个温度点对应的容量因子；②计算甲醇在串联后柱2上6个温度点的校正保留时间t_P2：t_P2＝t_M2(1+k₂’)       式5式5中：k₂’是根据第一步中2)计算出甲醇在独立柱2上每个温度点对应的容量因子；③甲醇在6个温度点下柱1和柱2串联柱的校正保留时间t_P为t_P1+t_P2；(5)甲醇在上述6个温度点下校正保留时间tP与实测保留时间t_R的相对误差RE％，由式6计算得到：RE%=tP-tRtR×100%]]>式6(6)改变起始计算条件，将t_M1＝0，t_M2＝t_M变为t_M1+Δt，t_M2-Δt，并按等步长变化，直至使起始计算条件为：t_M1＝t_M，t_M2＝0为止，按每一次改变起始条件值并按步骤(4)-(5)计算对应6个温度点下的相对偏差，其中Δt为0.001min-0.005min任意设定；(7)根据不断变化的t_M1、t_M2，选取甲醇在这6个温度点下相对偏差值最小的死时间t_M1为串联后柱1上6个温度点校正后的死时间，及相对偏差值最小的死时间t_M2为串联后柱2上6个温度点校正后的死时间；2)确定柱1和柱2串联柱在任意温度点下的死时间：(1)将第二步中1)确定的串联后柱1在6个温度点下的死时间及对应的温度值分别代入式7中，得到包含6个方程的方程组，通过解方程组，得到串联后柱1死时间t_M1的未知系数a′、b′、c′和d′：t_M1＝a′T³+b′T²+c′T+d′ 式7将已确定的系数a′、b′、c′和d′再代回式7中，得到串联后柱1的死时间与串联柱温度的对应关系；(2)串联后柱2的死时间与串联柱温度的对应关系的确定过程与步骤(1)中确定串联后柱1死时间与串联柱温度的对应关系相同；第三步确定柱1和柱2串联柱程序升温段中任意时刻的温度值：1)确定程序升温需要的总时间t_total：(1)当采用单阶程序升温时用式8确定程序升温需要的总时间t_total：t_total＝(T_f-T₀)/v        式8式8中：T₀为起始温度，该温度值为20℃-30℃任意设定，T_f为终止温度，该温度值为30℃-220℃任意设定，v为升温速率，该值为0℃/min-30℃/min任意设定；(2)当采用多阶程序升温时用式9确定程序升温需要的总时间t_total：t_total＝t_h1+t₁+t_h2+t₂         式9式9中：t_h1为起始温度保持时间，该值为0min-20min任意设定；t₁为第一阶段程序升温需要的时间，由式10确定：t₁＝(T_m-T₀)/v₁                式10式10中：T_m为第一阶段程序的终止温度，该温度值为30℃-70℃任意设定，v₁为第一阶段程序升温的升温速率，其值为0℃/min-20℃/min任意设定；t_h2为第一阶段程序的终止温度保持时间，该值为0min-20min任意设定；t₂为第二阶段程序升温需要的时间，由式11确定：t₂＝(T_f-T_m)/v₂                式11式11中：T_f为第二阶段程序的终止温度，该温度值为70℃-220℃任意设定，v₂为第二阶段程序升温的升温速率，该值为0℃/min-30℃/min任意设定；2)确定流动相由第0块跳到第i块塔板经历的时间t_j：tj=Σj=1mτj]]>式12式12中：m为流动相在色谱柱内的跳跃次数，其值在n＜m＜10000间任意设定，其中n为柱1和柱2串联柱固有的塔板数，为已知值，且柱1固有的塔板数为q，柱2固有的塔板数为n-q；τ_j为流动相第j次跳跃经历的时间，通过式13计算：τj=tMj-1n(j=1,······m)]]>式13式13中：为流动相第j-1次跳跃后所在温度下的死时间，它通过流动相跳跃j-1次经历的时间t_j-1加以确定，确定过程如下：(1)确定t_j-1时刻柱1和柱2串联柱的温度T_j-1：①对单阶程序升温时t_j-1时刻柱1和柱2串联柱的温度T_j-1确定：当t_j-1＜t_total时：              T_j-1＝v×t_j-1+T₀                  式14当t_j-1≥t_total时：              T_j-1＝T_f②对多阶程序升温时t_j-1时刻柱1和柱2串联柱的温度T_j-1确定：当t_j-1≤t_h1时：                 T_j-1＝T₀当t_h1＜t_j-1＜t_h1+t₁时：         T_j-1＝v₁×(t_j-1-t_h1)+T₀           式15当t_h1+t₁≤t_j-1≤t_h1+t₁+t_h2时：  T_j-1＝T_m当t_h1+t₁+t_h2＜t_j-1＜t_total时：  T_j-1＝v₂×(t_j-1-t_h1-t₁-t_h2)+T_m    式16当t_j-1≥t_total时：              T_j-1＝T_f(2)根据串联柱的温度T_j-1，确定t_j-1时刻柱1和柱2串联柱的死时间；将第三步2)中(1)计算得到的T_j-1值，代入式7中，分别计算在此温度值下串联柱中柱1、柱2各自的死时间，然后相加，得到串联柱的值；(3)将的值代入式13和式12得到t_j的值；3)确定t_j时刻柱1和柱2串联柱的温度T_j：过程与第三步2)中(1)的方法相同；第四步待测组分保留时间的预测：1)将t_j时刻柱1和柱2串联柱的温度T_j带入独立柱2容量因子与温度的关系式中，计算待测组分在温度Tj时独立柱2的容量因子k_b′；2)运用温度T_j对应的独立柱2的容量因子值，通过式17和式18，分别计算待测组分在柱1和柱2串联柱中第q+1块，第q+2块，...，第n块塔板中流动相和固定相的浓度：k_b′×β_b＝C_Sij/C_Mij                式17C_Sij+C_Mij＝C_Si(j-1)+C_M(i-1)(j-1)    式18式17中：β_b为独立柱2的相比，为已知值；C_Sij和C_Mij分别为在t_j时刻，待测组分在第i块塔板中固定相和流动相里的浓度，其中q+1≤i≤n；式18中：C_Si(j-1)为t_j-1时刻，待测组分在第i块塔板中固定相的浓度；C_M(i-1)(j-1)为t_j-1时刻，待测组分在第i-1块塔板中流动相的浓度，其中，C_M00＝1μg/ml；3)判断待测组分在柱1和柱2串联柱中流动相里的最大浓度：将t_j时刻的待测组分在第n块塔板中流动相里的平衡浓度C_Mnj值，与t_j-1时刻的待测组分在第n块塔板中流动相里的平衡浓度C_Mn(j-1)值，及与t_j+1时刻的待测组分在第n块塔板中流动相里的平衡浓度C_Mn(j+1)的值进行比较：当C_Mn(j-1)<C_Mnj>C_Mn(j+1)时，停止计算；此时对应的流动相完成j次跳跃的时间tj即为保留时间t_p’；当不满足不等式C_Mn(j-1)<C_Mnj>C_Mn(j+1)，则继续重复上述第三步、第四步中的计算，直至达到满足C_Mn(j-1)<C_Mnj>C_Mn(j+1)不等式为止；4)根据第三步中1)确定的程序升温条件，测定待测组分在柱1和柱2串联柱上的保留时间t_R’；5)预测值的验证：RE%=tp′-tR′tR′×100%]]>式19根据式19，经第四步中3)得到的待测组分保留时间预测值t_p’减去待测组分实测的保留时间值t_R’，得到偏差值，用偏差值除以实测的保留时间值t_R’即为待测组分的保留时间相对偏差，以此来对预测结果进行验证。