[发明专利]预测多阶程序升温条件下伸舌色谱峰形的方法

权利要求书

1.一种预测多阶程序升温条件下伸舌色谱峰形的方法，该方法是以非线性塔板理论为基础，针对HP6890气相色谱仪及非极性的HP-5色谱柱，预测过程中采用任意恒定的虚拟死时间，其特征在于包括以下过程：

1)虚拟死时间τ的设定：色谱柱设定温度变化范围为30-250℃，测定待测化合物在T₁＝30℃、T₂＝50℃、T₃＝100℃、T₄＝150℃、T₅＝200℃和T₆＝250℃六个恒温下的保留时间t_R1、t_R2、t_R3、t_R4、t_R5和t_R6，确定其中最小保留时间值，以小于该最小保留时间值的任意一个时间值均可作为虚拟的死时间；

2)恒温条件下保留因子k的计算：根据步骤1)测定的六个恒温下的保留时间t_R1、t_R2、t_R3、t_R4、t_R5、t_R6和已经确定的虚拟死时间τ，采用式1，计算六个恒温下对应的保留因子：k₁、k₂、k₃、k₄、k₅和k₆，

k＝(t_R-τ)/τ                         式1

式1中：k为保留因子，

       τ为虚拟的死时间，

       t_R为对应各个温度点的保留时间；

3)恒温条件下峰形参数(peak shape parameter，psp)的计算：

峰形参数psp定义为基于非线性塔板理论建立的待测化合物在同一块塔板上固定相与流动相中的浓度之间的函数关系中的系数，函数关系如下式：

CS=psp·CM2+(k+1)·CM]]>式2

C_S+C_M＝C_n                            式3

式2中：psp为峰形参数值；

       C_S、C_M为待测化合物在同一块塔板上固定相和流动相中的浓度；

式3中：C_n为待测化合物在同一块塔板上固定相与流动相中浓度之和；

以下述计算机程序运行，计算不同温度条件下的峰形参数值：

首先在程序中输入已知参数：恒温条件下色谱柱温度T、待测化合物的初始浓度C_M00、待测化合物在该温度条件下的保留因子k、理论塔板数N以及跳跃次数l；再利用式2和式3，输入三个不同的峰形参数值，分别对比计算机程序运行获得的峰形与实验峰形，得出峰形参数与峰形大小之间变化趋势的关系；不断调整输入的峰形参数值的大小，直到获得的该恒温条件下的峰形与实验峰形完全吻合时，该数值便为该待测化合物在该温度条件下的实际峰形参数值；

利用第i次跳跃对应的保留因子k_i，以及设定的峰形参数值psp，根据式4和式5，分别计算待测化合物在第1、2、3...n...N块塔板中固定相和流动相的浓度：

CSni=psp·CMni2+(k+1)·CMni]]>式4

C_Sni+C_Mni＝C_ni                        式5

式4中：C_Sni和C_Mni分别为跳跃次数为i时，待测化合物在第n块塔板中固定相和流动相中的浓度；

式5中：C_ni为跳跃次数为i时，待测化合物在第n块塔板中的总浓度，它由式6确定：

C_ni＝C_S(n-1)i+C_M(n-1)(i-1)           式6

C_M00＝1μg/mL

式6中：C_M00为待测化合物的起始浓度，

       C_S(n-1)i为跳跃次数为i时，待测化合物在第n-1块塔板的固定相中的浓度，

       C_M(n-1)(i-1)为跳跃次数为i-1时，待测化合物在第n-1块塔板的流动相中的浓度；

当待测化合物跳跃至最后一块塔板时，待测化合物仍按照式4、式5关系进行分配，当进行下一次跳跃时，最后一块塔板上流动相中的待测化合物流出，记录该时刻流出的浓度值，固定相中的待测化合物与上一塔板流动相带入的待测化合物重新进行分配，如此，直至达到设定的跳跃次数，保证待测化合物全部流出色谱柱，由此得到待测化合物在每一次跳跃时对应流出浓度值的数据点，并由这些数据点构成的图形即为预测的伸舌色谱峰形；

在上述计算机程序中分别设定T₁、T₂、T₃、T₄、T₅和T₆六种温度数值，并由此得到相应的峰形参数值psp₁、psp₂、psp₃、psp₄、psp₅和psp₆；

4)确定待测化合物在程序升温条件下色谱柱内任意一次跳跃时，对应的温度Ti：

(1)以式6确定每次跳跃需要的时间Δτ，

Δτ＝τ/(N-1)                        式6

式6中：τ为虚拟的死时间，由步骤1)已确定，

        N为色谱柱固有的理论塔板数；

(2)待测化合物在色谱柱内跳跃i次，共需时间t_i，由式7计算，

t_i＝i×Δτ                           式7

式7中，i是跳跃次数；

(3)在多阶程序升温中，计算多阶程序升温的总时间t：

t＝t_h1+t₁+t_h2+t₂                        式8

式8中：t_h1为起始温度的保持时间，经验数值为：1-5min，

       t_h2为第一阶段程序升温的终止温度的保持时间，经验数值为：1-5min，

       t₁为第一阶段程序升温的需要时间，

       t₂为第二阶段程序升温的需要时间，

       t₁和t₂分别由式6和式7计算得到：

t₁＝(T_m-T₀)/r₁                         式9

t₂＝(T_f-T_m)/r₂                         式10

式9中：T_m为第一阶段程序升温的终止温度，

       T₀为起始温度，

       r₁为第一阶段程序升温的升温速率；

式10中：T_f为第二阶段程序升温的终止温度，

        r₂为第二阶段程序升温的升温速率，

其中，r₁和r₂的经验取值范围为5-30℃/min；

(4)确定待测化合物第i次跳跃时，此时色谱柱对应的温度T_i：

当t_i＜t_h1，则柱温T_i＝T₀，

当t_h1＜t_i＜(t_h1+t₁)，则柱温T_i＝r₁×(t_i-t_h1)+T₀，

当(t_h1+t₁)≤t_i≤(t_h1+t₁+t_h2)，则柱温T_i＝T_m，

当(t_h1+t₁+t_h2)＜t_i＜t，则柱温T_i＝r₂×(t_i-t₁-t_h1-t_h2)+T_m，

当t_i＞t，则柱温T_i＝T_f；

5)多阶程序升温条件下保留因子k及峰形参数psp与温度关系的确定：

(1)利用六个恒温温度值T₁、T₂、T₃、T₄、T₅和T₆以及步骤2)中求得的相应的保留因子k₁、k₂、k₃、k₄、k₅和k₆，通过回归分析，得到保留因子k与温度T之间的函数关系式为式11：

ln k＝aT³+bT²+cT+d           式11

其中参数a、b、c和d均为定值；

由此计算出程序升温中任意温度点T_i对应的保留因子k_i；

(2)利用六个恒温温度值T₁、T₂、T₃、T₄、T₅和T₆以及步骤3)中获得的相应的峰形参数值psp₁、psp₂、psp₃、psp₄、psp₅和psp₆，通过回归分析，得到峰形参数psp与温度T之间的函数关系式为式12：

ln psp＝a′T³+b′T²+c′T+d′            式12

式12中，参数a′、b′、c′和d′均为定值；

由此计算出程序升温中任意温度点T₁对应的峰形参数psp；

6)多阶程序升温条件下待测化合物预测峰形的获得：

在计算机程序中输入以下数据：柱相比β、待测化合物起始浓度C_M00、理论塔板数N、跳跃次数n、六个恒温温度值、虚拟的死时间τ、保留因子和峰形参数与温度关系式中的各参数、多阶程序升温的初始温度、第一阶段终止的温度、第二阶段终止的温度、第一阶段的升温速率、第二阶段的升温速率、初始温度下保持的时间和第一阶段程序升温终止温度下保持的时间；通过步骤3)中所述的计算机运行程序，得到待测化合物在该程序升温条件下的预测的伸舌色谱峰形。