专利名称：色谱柱短化后在程序升温下气相色谱的保留时间的预测方法
技术领域：
本发明涉及一种色谱柱短化后在程序升温下气相色谱保留时间的预测方法，属于色谱分析技术领域。
背景技术：
色谱保留时间是色谱分析定性的基础，也是鉴别化合物最佳分离条件的一个重要的依据。目前保留时间的预测方法，都是通过恒温数据对相同柱长的同一根色谱柱在一定的温度、压力或流量程序条件下的保留时间进行预测。但随着色谱柱的频繁使用致使柱效损坏，及在实验操作过程中不同型号色谱柱间的更换不当导致柱子的折断，使色谱柱的长度明显的短化。我们测的三种样品在短化27%的HP-5毛细管柱上保留时间，与在原长HP-5 毛细管柱测得的三种样品在恒温下的保留时间有着明显的变化，这种差异在保留时间越长的样品上就越发明显。如萘在80°C恒温条件下柱长改变前后的保留时间差值大于lOmin， 乙酸丁酯在30°C恒温条件下柱长改变前后的保留时间差值大于7min。这种差异对于保留时间短的样品也是不能忽略的，如乙醇在30°C恒温条件下柱长改变前后的保留时间差值大于lmin。在高温下如250°C下三种样品柱长改变前后的保留时间差值也都近似lmin。即所有样品的保留时间较柱子折断前至少减少了 40%，这说明两个问题第一，两次测量结果差异极大；第二，保留时间的减少不只受柱长的影响，还受到柱内流量、压力变化的影响。因此当毛细管柱的长度发生短化后，需要重新测量各样品在短柱中恒温条件下的保留时间，利用新的恒温数据对样品程序升温保留时间进行预测。因此每当色谱柱长度明显变短后，都需要舍弃以前长柱的恒温数据，而重新测量短柱的恒温数据进行程序升温保留时间预测。这样不仅增加了最佳色谱条件选择的试验次数，并且延长了分析时间，给色谱分析工作带来很大不便。在这种情况下能否继续使用原长色谱柱的恒温数据，在程序升温条件下预测组分在短化后的色谱柱上的保留时间，具有十分重要的意义。其技术关键包括两个方面，一是方法的建立，二是短化后的色谱柱塔板数的校正。有关文献报道，根据保留时间锁定的方法可以解决色谱柱变短带来的保留时间明显变化的问题。即通过测定不同压力下的保留时间，找到压力与保留时间之间的关系并进行复杂地校正，达到在相应的条件下锁定某一化合物的保留时间的目的。但是这种方法成本太高，不能广泛普及应用。
发明内容
本发明的目的在于提供一种色谱柱短化后在程序升温下气相色谱保留时间的预测方法，该方法具有过程简单、实用、准确和可信度高的特点。本发明是通过以下技术方案加以实现的，一种色谱柱短化后在程序升温下气相色谱的保留时间的预测方法，该方法以非线性塔板理论为基础，针对HP-6890气相色谱仪配置的HP-5毛细管色谱柱，运用色谱柱短化前的长柱原有的恒温保留时间，预测程序升温下组分在短化后色谱柱中的气相色谱保留时间，其特征在于包括以下过程1)确定原长色谱柱在任意温度点对应的死时间和容量因子(1)确定原长色谱柱任意温度点对应的死时间t^i的值将在温度30°c、50°c、10(rc、15(rc、20(rc和250°C六个温度点下测得甲烷的保留时间作为死时间，将六个温度点的死时间与相应的温度值分别代入式1得到包含六个方程的方程组，通过解方程组，得到系数a'、b'、c'和d'的值tMJ_! = a' T3+b' T2+c' T+d' 式 1将确定后系数a'、b'、c'和d'的值再代入式1中，得到在原长色谱柱上任意温度点与死时间的关系；(2)确定原长色谱柱任意温度点对应的容量因子k'的值①采用步骤1)中(1)确定的系数a'、b'、c'和d'的值，通过式1计算出原长柱在30°c、50°c、10(rc、15(rc、20(rc和250°C六个温度点下的死时间；②将组分在温度30°c、50°c、10(rc、15(rc、20(rc和250°C六个温度点下的保留时间，及步骤1) (2)中①确定的相应温度点下的死时间，通过式2得到组分六个温度点下相应的各个容量因子k'k' = tE/tm"l式 2式2中tK为组分在对应温度点的保留时间，tm为对应温度点的死时间，k'为组分在对应温度点的容量因子；③将六个温度点下的各自容量因子与相应的温度分别代入式3得到包含六个方程的方程组，通过解方程组，得到系数a、b、c和d的值Ink' = aT3+bT2+cT+d式 3将确定后系数a、b、c和d的值再代入式3中，得到该组分在原长色谱柱上任意温度点与容量因子的关系；2)在恒温条件下，以苯为标准物校正短化后色谱柱的塔板数η (1)计算短化后色谱柱长1与原长色谱柱长L的比值r ；(2)利用步骤(1)中的比值r，确定短化后色谱柱的塔板数nx nx = N*{l-r-0. 01*(x_l)}式 4式4中x为校正短柱的塔板数需要计算的次数，其中1 < χ且χ为整数，nx为第χ次计算时对应的短化后色谱柱的塔板数，N为原长色谱柱的固有的塔板数，为已知；(3)计算恒温温度T下，短化后色谱柱的塔板数为nx时，流动相每次跳跃经历的时间At:At = tm/nx式 5式5中tm为指定温度T下原长色谱柱的死时间，已由式1计算得到，其中T值在60°C _200°C任意设定；(4)计算恒温温度T下，流动相在短化后色谱柱内跳跃j次经历的时间、tj = j* Δ t式 6其中nx彡j且j为整数；
(5)确定恒温温度T下，标准物苯在原长色谱柱的容量因子过程与步骤1)中O)的方法相同；(6)根据步骤( 中的温度T下标准物苯在原长色谱柱的容量因子，计算流动相在短化后色谱柱内跳跃j次后，苯在第0块，第1块，---，第nx_l块塔板中流动相和固定相的浓度k' = CSiJ/CMiJ式 7Csij+CMi」=Csicj-D+Cm(H)(计式 8式7中β为色谱柱的相比，对于已选定的色谱柱而言，它为已知值，Csij和Cmu分别为在、时刻，组分在第i块塔板中固定相和流动相里的浓度，式8中=Csi(H)为ty时刻，组分在第i块塔板中固定相的浓度，Cm(H)(H)为tH时刻，组分在第块塔板中流动相的浓度；其中，0彡 i 彡 nx-l，Cmoo = 1 μ g/ml ；(7)通过步骤(6)中确定的组分在各块塔板中流动相的浓度，判断组分在短化色谱柱中流动相里的最大浓度将、时刻的组分在第nx_l块塔板中流动相的平衡浓度CM(I1)j值，与时刻的组分在第nx_l块塔板中流动相的平衡浓度(H)值，及与tm时刻的组分在第nx_l块塔板中流动相的平衡浓度Cn^i -ι) 0+i)的值进行比较当CM( ,-i)(j-i)< Cm(Bi-I) j> Cm( -i)(j+i)时，停止计算；此时对应的流动相完成j次跳跃的时间tj即为温度T下标准物苯的校正保留时间 tpE ；当不满足不等式<^(>1,-1)()-1)<(叫 ,-1)〗>(叫 「1)(]+1)，则继续重复步骤2)中(6)-(7) 的计算，直至达到满足CMK-I)(κ)< CM(nx-i)j> Cm(^1)(H)不等式为止；(8)测定恒温温度T下，标准物苯在短化后色谱柱的保留时间teK (9)计算恒温温度T下，标准物苯的校正保留时间tpK与实测保留时间U的相对
偏差E%= (teE-tpE)/teE*100式 9当标准物苯的校正保留时间tpK与实测保留时间teK的相对偏差的绝对值小于时，对应的塔板数为短化后色谱柱的校正塔板数η ；当标准物苯的校正保留时间tpK与实测保留时间teK的相对偏差绝对值大于时，重复步骤2)中(2)-(9)的计算，直至标准物苯的校正保留时间tpK与实测保留时间teK 的相对偏差E%绝对值小于；3)利用步骤1)所确定的原长色谱柱的任意温度点对应的死时间和容量因子，及步骤2)确定的短化后色谱柱的校正塔板数n，在程序升温下，预测组分在短化色谱柱中的保留时间(1)确定程序升温需要的总时间tt。tal ①当采用单阶程序升温时用式10确定程序升温需要的总时间tt。tal ttotal = (Tf-T0)/ν式 10式10中Tf为设定的终止温度，该温度在30°C _250°C任意设定，T0为起始温度，该温度值在20°C -30°C任意设定，
ν为升温速率，该值在0°C /min-30°C /min任意设定；②当采用多阶程序升温时用式11确定程序升温需要的总时间tt。tal ttotal = th+t!+t2式 11式11中th为起始温度保持时间，该值在0min-20min任意设定；tl为第一阶段程序升温需要的时间，由式12确定tl = (Tm-T0Vv1式 12式12中Tm为第一阶段程序升温的终止温度，该温度值在30°C -70°C任意设定，V1为第一阶段程序升温的升温速率，其值在0°C /min-20°C /min任意设定；t2为第二阶段程序升温需要的时间，由式13确定t2 = (Tf-Tm) /v2式 13式13中Tf为第二阶段程序升温的终止温度，该温度值在70°C -250°C任意设定，V2为第二阶段程序升温的升温速率，该值在0°C /min-30°C /min任意设定；(2)确定程序升温条件下，流动相跳跃j次经历的时间、由式14确定Q=Xpi η式 14式14中m为流动相在色谱柱内的跳跃次数，其值在η < m < 10000间任意设定； τ j为流动相第j次跳跃经历的时间，通过式15计算τ j = t^^/n (j = 1，……m)式 15式15中‘为流动相第j_l次跳跃后所在温度下的死时间，它通过程序升温条件下，流动相跳跃j-ι次经历的时间y对应的短化色谱柱的温度IV1加以确定，确定过程如下①确定ty时刻短化色谱柱的温度IV1 (一 )对单阶程序升温时th时刻短化色谱柱的温度IV1确定当ty < ttotal时，tjM时刻的短化色谱柱温度Tj+通过式16计算Tj^1 = v*t+T0式 16当V1彡ttotal时，V1时刻短化色谱柱温度IV1 = Tf ；( 二)对多阶程序升温时th时刻短化色谱柱温度IV1确定当t」_i彡th时，V1时刻的短化色谱柱温度IV1 = T0 ；当th < tjM ( ti+th时，tjM时刻的短化色谱柱温度IV1由式17计算Tj^1 = V1* (t-th) +T0 式 17当ti+th < tjM < ttotal时，V1时刻短化色谱柱的温度IV1由式18计算I^1 = v2* (t-tfth)+Tm 式 18当V1彡tttotal时，V1时刻短化色谱柱的温度IV1 = Tf ；②根据色谱柱的温度IV1，确定ty时刻短化色谱柱的死时间t^ 将步骤；3)⑵中①计算得到的IV1值，代入式1计算得到t^i的值；③将t^i的值代入式14和式15得到、的值(3)确定、时刻短化色谱柱的温度Tj 过程与步骤3)⑵中①的方法相同；(4)确定温度Tj对应的长柱的容量因子k'的值将&的值，代入式3中，得到短化色谱柱的温度为L时，相应的长柱的容量因子值；(5)运用温度L对应的长柱的容量因子值，借助式7和式8，分别计算组分在短化色谱柱中第0块，第1块，一，第n-1块塔板中流动相和固定相的浓度；(6)通过步骤(5)中确定的组分在各块塔板中流动相的浓度，判断组分在短化色谱柱中流动相里的最大浓度将、时刻的组分在第Π-1块塔板中流动相的平衡浓度CM(n_⑴值，与ty时刻的组分在第n-1块塔板中流动相的平衡浓度Cm^)(η)值，及与tm时刻的组分在第n-1块塔板中流动相的平衡浓度Cm^) +1)的值进行比较当CM^kh^CMkDpCMkD&d 时，停止计算；此时对应的流动相完成j次跳跃的时间、即为相应程序升温条件下，组分在短化后色谱柱中的保留时间tK;当不满足不等式(^^……〈(临-们冗临-…州)，贝丨」继续重复步骤3)中的计算，直至达
到 两足 Cm(^1) (j—:0〈CM(n—山冗恤―d (j+1) 不等式为止；4)根据步骤；3)中(1)确定的程序升温条件，测定组分在短化后色谱柱上的保留时间；5)预测值的验证采用经步骤3)得到保留时间的预测值减去步骤4)相应组分保留时间的实测值得到偏差值，偏差值除以保留时间的实测值即为组分保留时间预测的相对偏差，当相对偏差的绝对值小于2%时，认为预测精度达到要求。有益效果本发明具有以下优点本发明运用组分在长柱中恒温下的保留时间，预测了其在短柱中程序升温下的保留时间。不仅节约了实验成本，而且降低了分析工作的强度，缩短了分析的时间，为短柱的气相色谱最佳分离条件的选择提供了可行途径。除此之外，本发明再现了组分在色谱柱两相内真实的、动态的分配过程，因此不需要加入许多额外的参数进行预测，预测过程清晰、明了、适合大范围地普及应用。
图1为本发明预测在程序升温下组分在短化后色谱柱中保留时间的计算模拟流程图。
具体实施例方式实施例1仪器HP6890气相色谱仪，氢火焰离子化检测器，6890气相色谱工作站；色谱柱非极性的HP-5 (5%苯基甲基聚硅氧烷)柱；原长，HP-5 (Crosslinked 5% PH ME SiIoxane，30mX0. 32mmX0. 25 μ m)；短化后，HP-5 (Crosslinked 5% PH ME Siloxane, 22mX0. 32mmX0. 25 μ m)；条件检测器的温度为250°C，进样口温度为250°C ；载气使用高纯氮气(纯度不低于99. 999%)，恒流操作模式，即载气在柱出口处， 质量流量保持恒定，为lml/min ；
进样方式分流进样，分流比为50 1，每次的进样量为0.2 μ 1，初始浓度Cmcici = 1 μ g/ml ；程序升温条件程序升温起始温度Ttl = 30°C、升温速率ν = 5°C /min、终止温度乙 =250 0C ；样品萘、乙酸丁酯、甲醇、甲烷、苯；(一)调用甲烷在原长 30m 的 HP-5 柱中 30°C、50°C、100°C、150°C、200°C和 250°C 六个温度点下的保留时间，它们分别为2. 78min、2. 63min、2. 38min、2. 17min、2. OOmin和
1.85min(甲烷的保留时间看作死时间)；(二)将步骤(一)中六个温度点下的死时间与相应的温度代入式1，得到式1 的系数 a'、b'、c'和 d'的值，它们分别为-0. 00000003188、0. 00002185、-0. 08045 和
2.994。将系数af、b'、c'和d'的值代入式1，得到在原长30m的HP-5柱中死时间与温度的关系tMJ_! = -0. 00000003188 X T3+0. 00002185 X T2-O. 08045 X Τ+2. 994 ；(三)分别调用乙酸丁酯、甲醇、标准物苯和萘在原长30m的HP_5柱中30°C、50°C、 100°C、150°C、200°C和250°C六个温度点下的保留时间，它们分别对应为20. 26min、9. 10min、3. 44min、2. 53min、2. 21min 和 2. 02min ；3. 35min、3. 01min、2. 60min、2. 34min、2. 14min 和 1. 99min ；6. 74min、4. 53min、2. 89min、2. 41min、2. 17min 和 2. OOmin ；由于萘在低温下保留时间大于1个小时，分析时间太长，因此只调用了萘在原长30m的HP-5柱中100°C、150°C、200 V和250°C四个温度点下的保留时间，它们分别为 12. 60min、4. 09min、2. 63min 和 2. 18min ；(四)将步骤(一)中六个温度点下的死时间与乙酸丁酯在相应的温度下的保留时间分别代入式2，得到乙酸丁酯在六个温度点下的容量因子，它们分别为6. 286,2. 462、 0.4441,0. 1668,0. 1053 和 0. 09162 ；(五)按照步骤(四)的计算过程，分别得到甲醇和标准物苯在六个温度点下的容量因子它们分别为 0. 2036,0. 1452,0. 0908,0. 07719,0. 07025 和 0. 07432 ；1. 425,0. 7230,0. 2122,0. 1127,0. 08275 和 0. 07946 ；同样按照步骤(四)的计算过程，得到萘在四个温度点下的容量因子，它们分别为 4.292,0. 8825,0. 3130 和 0. 1768 ；(六)将步骤(四)中乙酸丁酯在六个温度点下的容量因子及相应的温度代入式 3，得到乙酸丁酯的容量因子与温度关系的系数a、b、c和d的值，它们分别为a(乙酸丁酯) =-0. 0000003, b(乙酸丁酯)=0.0005，c(乙酸丁酯)=-0.2507，和 d(乙酸丁酯)= 43. 696 ；(七)按照步骤(六)的计算过程，分别得到甲醇、标准物苯和萘的容量因子与温度关系的系数a、b、c和d的值，它们分别为a(甲醇)=-0.0000002，b(甲醇)=0.0002， c(甲醇)=-0. 1136，和 d(甲醇)=15. 989 ；a(苯)=-0. 0000002，b(苯)=0. 0004， c (苯)=-0. 193，和 d(苯)=32. 127 及 a(萘)=-0. 0000001，b (萘)=0. 0003，c (萘) =-0. 2006，和 d (萘)=42. 197 ；(八)用标准物苯在短柱100°C下的实测保留时间，校正短化后色谱柱的塔板数η 将以下参数的值原长色谱柱长L = 30m，短化后色谱柱长1 = 22m, T = 100°C，tm =2. 38min (100°C ), teE = 1. Hmin (苯在短柱100°C下的实测保留时间)、N = 10000、组分的浓度 Cmoo = 1 μ g/ml 及 a (苯)=-0. 00000003，b (苯)=0. 00005，c (苯)=-0. 0288， 和d (苯)=7. 7105输入，程序自动按照步骤幻计算，得到短化后色谱柱的校正塔板数η =1000 ；(九)将以下参数的值程序起始温度Ttl= 30°C、升温速率ν = 5 0C /min、终止温度Tf = 250°C、组分的浓度Cmcici= 1 μ g/ml、短化后色谱柱的校正塔板数η = 1000、 a ‘ =-0. 00000003188，b ‘ = 0. 00002185, c ‘ = -0. 08045 和 d ‘ = 2. 994 及 a(萘) =-0. 0000001，b(萘)=0.0003，c(萘)=-0. 2006，和 d(萘)=42. 197 ;a(乙酸丁酯) =-0. 0000003, b(乙酸丁酯)=0. 0005, c(乙酸丁酯)=-0. 2507 和 d(乙酸丁酯)= 43.696 ；a(甲醇)=-0.0000002，b(甲醇)=0.0002，c(甲醇)=-0. 1136 和 d(甲醇) =15. 989输入；程序自动按照步骤幻进行计算；为了验证预测值，需在同样单阶程序升温条件下，测量三种组分在变短后2 !的HP-5色谱柱上的真实保留时间。三种物质在上述程序升温条件下预测值与实测值的误差分析如表1所示表1升温速率ν = 5°C /min条件下三种化合物保留时间的实测值与预测值得比较
化合物保留时间(min)偏差相对偏差(％)实测值预测值萘16.3916.610.221.34乙酸丁酯6.376.400.030.47甲醇2.012.010.000.00实施例2本实施例的过程和条件与实施例1相同，不同的是程序升温条件程序升温起始温度Ttl = 30°C、升温速率ν = 250C /min、终止温度 Tf = 250 0C ；(一 )只需将实施例1中步骤(九)的输入的程序升温速率的参数数值改为升温速率ν = 250C /min ；其余参数值和预测过程与实施例1完全相同；(二)为了验证预测值，需在程序升温起始温度Ttl = 30°C、升温速率ν = 25 0C / min、终止温度Tf = 2500C的条件下，测量三种组分在变短后2 !的HP-5色谱柱上的真实保留时间。三种物质在上述程序升温条件下预测值与实测值的误差分析如表2所示表2升温速率ν = 250C /min条件下三种化合物保留时间的实测值与预测值得比较
1权利要求
1. 一种色谱柱短化后在程序升温下气相色谱的保留时间的预测方法，该方法以非线性塔板理论为基础，针对HP-6890气相色谱仪配置的HP-5毛细管色谱柱，运用色谱柱短化前的长柱原有的恒温保留时间，预测程序升温下组分在短化后色谱柱中的气相色谱保留时间，其特征在于包括以下过程1)确定原长色谱柱在任意温度点对应的死时间和容量因子(1)确定原长色谱柱任意温度点对应的死时间^^的值将在温度30°C、50°C、100°C、150°C、200°C和250°C六个温度点下测得甲烷的保留时间作为死时间，将六个温度点的死时间与相应的温度值分别代入式1得到包含六个方程的方程组，通过解方程组，得到系数a'、b'、c'和d'的值 V1 = a' T3+b' T2+c' T+d' 式 1将确定后的系数a'、b'、(和(1'的值再代入式1中，得到在原长色谱柱上任意温度点与死时间的关系；(2)确定原长色谱柱任意温度点对应的容量因子k'的值①采用步骤1)中(1)确定的系数a'、b'、c'和d'的值，通过式1计算出原长柱在 30°c、50°c、100°c、15(rc、20(rc和250°C六个温度点下的死时间；②将组分在温度30V、50°C、100°C、150°C、200 V和250 V六个温度点下的保留时间， 及步骤1) (2)中①确定的相应温度点下的死时间，通过式2得到组分六个温度点下相应的各个容量因子k'k' = tE/tm-l式 2式2中tK为组分在对应温度点的保留时间， tm为对应温度点的死时间， k'为组分在对应温度点的容量因子；③将六个温度点下的各自容量因子与相应的温度分别代入式3得到包含六个方程的方程组，通过解方程组，得到系数a、b、c和d的值In k' = aT3+bT2+cT+d式 3将确定后的系数a、b、c和d的值再代入式3中，得到该组分在原长色谱柱上任意温度点与容量因子的关系；2)在恒温条件下，以苯为标准物校正短化后色谱柱的塔板数η(1)计算短化后色谱柱长1与原长色谱柱长L的比值r；(2)利用步骤(1)中的比值r，确定短化后色谱柱的塔板数nx: nx = N*{l-r-0. 01*(x-l)}式 4式4中X为校正短柱的塔板数需要计算的次数，其中1 < χ且χ为整数， ηχ为第χ次计算时对应的短化后色谱柱的塔板数， N为原长色谱柱的固有的塔板数，为已知；(3)计算恒温温度T下，短化后色谱柱的塔板数为ηχ时，流动相每次跳跃经历的时间At At = tm/nx式 5式5中tm为指定温度T下原长色谱柱的死时间，已由式1计算得到， 其中T值在60°C _200°C任意设定；(4)计算恒温温度T下，流动相在短化后色谱柱内跳跃j次经历的时间、 tj = j* Δ t式 6其中j且j为整数；(5)确定恒温温度T下，标准物苯在原长色谱柱的容量因子 过程与步骤1)中(2)的方法相同；(6)根据步骤(5)中的温度T下标准物苯在原长色谱柱的容量因子，计算流动相在短化后色谱柱内跳跃j次后，苯在第0块，第1块，一-，第nx_l块塔板中流动相和固定相的浓度k' = CSiJ/CMiJ式 7Csij+。Mij — CSi(j-l)+CM(i-l) (j-1) 式 8式7中β为色谱柱的相比，对于已选定的色谱柱而言，它为已知值， Csij和Cmu分别为在、时刻，组分在第i块塔板中固定相和流动相里的浓度， 式8中=Csi(H)为y时刻，组分在第i块塔板中固定相的浓度， CMa-Do-D为ty时刻，组分在第i-Ι块塔板中流动相的浓度； 其中，0 ( i ( nx-l，Cmoo = 1 μ g/ml ；(7)通过步骤(6)中确定的组分在各块塔板中流动相的浓度，判断组分在短化色谱柱中流动相里的最大浓度将、时刻的组分在第nx_l块塔板中流动相的平衡浓度Cm(Y1) j值，与V1时刻的组分在第nx_l块塔板中流动相的平衡浓度CM(Y1)(HL)值，及与tm时刻的组分在第nx_l块塔板中流动相的平衡浓度-1) (j+i)的值进行比较当C1m( -ι)(Η)< CM(Bi-I)j> Cm( -i)(j+i)时，停止计算；此时对应的流动相完成j次跳跃的时间、即为温度T下标准物苯的校正保留时间tpK ； 当不满足不等式CMk-DG-i^CMdDjXTMh-Dij+D，则继续重复步骤2)中(6)-(7)的计算，直至达到满足cM(n,-i)G-i)<j> <^Μ( ,-ι)( +ι)不等式为止；(8)测定恒温温度T下，标准物苯在短化后色谱柱的保留时间teK(9)计算恒温温度T下，标准物苯的校正保留时间tpK与实测保留时间U的相对偏差E% E%= (teE-tpE)/teE*100式 9当标准物苯的校正保留时间tpK与实测保留时间teK的相对偏差的绝对值小于时，对应的塔板数为短化后色谱柱的校正塔板数η ；当标准物苯的校正保留时间tpK与实测保留时间teK的相对偏差绝对值大于时，重复步骤2)中(2)-(9)的计算，直至标准物苯的校正保留时间tpK与实测保留时间U的相对偏差E %绝对值小于；3)利用步骤1)所确定的原长色谱柱的任意温度点对应的死时间和容量因子，及步骤 2)确定的短化后色谱柱的校正塔板数n，在程序升温下，预测组分在短化色谱柱中的保留时间(1)确定程序升温需要的总时间tt。tal ①当采用单阶程序升温时用式10确定程序升温需要的总时间tt。tal ttotal = (Tf-T0)A式 10式10中Tf为设定的终止温度，该温度在30°C _250°C任意设定， T0为起始温度，该温度值在20°C -30°C任意设定， ν为升温速率，该值在0°C /min-30°C /min任意设定； ②当采用多阶程序升温时用式11确定程序升温需要的总时间tt。tal ttotal = Vt^t2式 11式11中th为起始温度保持时间，该值在0min-20min任意设定； 、为第一阶段程序升温需要的时间，由式12确定 ti = (Tffl-T0)A1式 12式12中Tm为第一阶段程序升温的终止温度，该温度值在30°C _70°C任意设定， V1为第一阶段程序升温的升温速率，其值在0°C /min-20°C /min任意设定； t2为第二阶段程序升温需要的时间，由式13确定 t2 = (Tf-Tffl)/v2式 13式13中Tf为第二阶段程序升温的终止温度，该温度值在70°C _250°C任意设定，V2为第二阶段程序升温的升温速率，该值在0°C /min-30°C /min任意设定；(2)确定程序升温条件下，流动相跳跃j次经历的时间、由式14确定 t) = 式 η式14中m为流动相在色谱柱内的跳跃次数，其值在η < m < 10000间任意设定；τ j 为流动相第j次跳跃经历的时间，通过式15计算 τ j = tnjXj = 1, ......m)式 15式15中‘为流动相第j-Ι次跳跃后所在温度下的死时间，它通过程序升温条件下， 流动相跳跃j_l次经历的时间ty对应的短化色谱柱的温度IV1加以确定，确定过程如下①确定ty时刻短化色谱柱的温度IV1(一)对单阶程序升温时ty时刻短化色谱柱的温度IV1确定当ty < ttotal时，tj^时刻的短化色谱柱温度IV1,通过式16计算 Tj^1 = v*t+TQ式 16当y彡ttotal时，tj^时刻短化色谱柱温度IV1 = Tf ；(二)对多阶程序升温时ty时刻短化色谱柱温度IV1确定 当ty ( th时，tjM时刻的短化色谱柱温度IV1 = T0 ；当th < tj^ ( ti+th时，tj^时刻的短化色谱柱温度IV1由式17计算 IV1 = Vl*(t-th)+TQ式 17当ti+th < tjM < ttotal时，tj^时刻短化色谱柱的温度IV1由式18计算 Tw = v2*(t-t「th)+Tm式 18当y彡ttotal时，tj^时刻短化色谱柱的温度IV1 = Tf ；②根据色谱柱的温度IV1，确定ty时刻短化色谱柱的死时间tMH 将步骤3)⑵中①计算得到的IV1值，代入式1计算得到t^的值；③将t^的值代入式14和式15得到tj的值(3)确定、时刻短化色谱柱的温度Tj 过程与步骤3) (2)中①的方法相同；(4)确定温度L对应的长柱的容量因子k'的值将L的值，代入式3中，得到短化色谱柱的温度为L时，相应的长柱的容量因子值(5)运用温度L对应的长柱的容量因子值，借助式7和式8，分别计算组分在短化色谱柱中第0块，第1块，一，第n-1块塔板中流动相和固定相的浓度；(6)通过步骤(5)中确定的组分在各块塔板中流动相的浓度，判断组分在短化色谱柱中流动相里的最大浓度将、时刻的组分在第n-1块塔板中流动相的平衡浓度CM(n_⑴值，与ty时刻的组分在第n-1块塔板中流动相的平衡浓度Cm(m)(h)值，及与tm时刻的组分在第n-1块塔板中流动相的平衡浓度 ^M (n-1) (j+1) 的值进行比较当。Μ(η—l)(j—1)〈。Μ(η—1”〉。Μ(η—1)G+1)时，停止计算此时对应的流动相完成j次跳跃的时间、即为相应程序升温条件下，组分在短化后色谱柱中的保留时间tK;当不满足不等式 ^M (n-1) (j-1) 〈C.dAC.d&d，则继续重复步骤⑴中的计算，直至达到满足 Cm(n-1) (j-1)〈。Μ(η-1)」·〉。Μ(η-1) (j+1) 不等式为止；4)根据步骤幻中(1)确定的程序升温条件，测定组分在短化后色谱柱上的保留时间；5)预测值的验证采用经步骤3)得到保留时间的预测值减去步骤4)相应组分保留时间的实测值得到偏差值，偏差值除以保留时间的实测值即为组分保留时间预测的相对偏差，当相对偏差的绝对值小于2%时，认为预测精度达到要求。
全文摘要
本发明公开了一种色谱柱短化后在程序升温下气相色谱保留时间的预测方法。该方法过程包括首先确定原长色谱柱在任意温度点对应的死时间和容量因子、并以苯为标准物校正短化后色谱柱的塔板数；再次在程序升温下，预测组分在短化后色谱柱中的保留时间；然后根据预测步骤中确定的程序升温，测定组分在短化后的色谱柱上的保留时间；最后计算预测组分的保留时间的相对偏差，当相对偏差的绝对值小于2％时，认为预测精度达到要求。本发明优点节约了实验成本，降低了分析工作的强度，为短柱的气相色谱最佳分离条件的选择提供了可行途径，预测过程清晰、明了、适合大范围地普及应用。
文档编号G01N30/86GK102507813SQ201110286568
公开日2012年6月20日 申请日期2011年9月26日 优先权日2011年9月26日
发明者张兰兰, 敖敏, 皇甫旭丹, 范国樑 申请人:天津大学