一份色谱实验记录(二十)

上一篇 / 下一篇  2013-11-04 20:31:54

    在我接触过的色谱项目中,反式脂肪酸这个项目是比较特殊的,它的特殊在于两点。第一点就在于,它不是我们熟悉的含量检测,而是指标检测,需要在检测一系列组分后,再进行一个看上去比较复杂的运算,最终得到一个类似于比例关系的数值。这与我们通常接触到的色谱项目的很大区别,没有做过这个项目的人多不理解,有些人谈到反式脂肪酸时张嘴就来,但是并不一定能搞明白这个项目的真正含义,还有些人一谈到地沟油检测时,一张嘴就是要上反式脂肪酸,这些人多是完全没有自己做过反式脂肪酸的检测,仅仅是听过这五个汉字罢了。

    第二点,就是这个项目在定义存在着一定的语言模糊性,说到项目的模糊性,类似地沟油的检测也是这个问题,只是很多人没有意识到这个对检测工作的影响有多么严重,地沟油的检测以后再谈。 反式脂肪酸定义的模糊性表现在,我们谈到反式脂肪酸时,觉得是面对一个很清晰的名词,其实它并不是,我们并没有确定这个项目中到底需要是检测多少种组份,无论是检测5 种脂肪组分还是检测30多种脂肪组分的方法,都说自己在检测反式脂肪,这个项目好像自己就有生命力似的。在2000年左右看美国油脂协会的检测方法时,看到的项目是8-10种组分,而在两年后,我在一篇论文上就看到了检测15种组分方法,我自己参考的反式脂肪酸方法,是一份国标申请稿,检测的是26种组分,在我建成项目并完成程序后一段时间,居然又看到了30多种组分的论文,这个项目的内容居然会象细菌一样自己繁殖,太神奇了。之后再过几年,我看到一个检测婴儿奶粉中式脂肪酸的国标,居然又只有5 种组分,呃,还会收缩。检测5 种组分和几十种组分,这难度可差太远了,连100米的柱子都不用。

     做色谱,尤其是做气相,会接触不同规格的色谱柱。我用过的最长的柱子就是传说中的脂肪酸检测专用柱,100米长,这种柱子不多见,我只找到两个厂家的产品,一个是瓦里安的,一个是HP-88。当时一样各买了一根,对比实验也做了,分享效果差不太多。但是后来发现瓦里安的那根柱子居然是两根拼接的,在柱架子上面有一个联接头,把两根连在一起。做色谱的人最看不得柱子上有些不可预测的节点问题,更别说断柱重接了,于是看那个联接头怎么看怎么不爽,以至后来就不太愿意用它。

 

(2)

    反式脂肪的检测难点也很特殊,首先是个概念理解问题,就是你要明白这个项目在做什么,你又需要做什么,它的检测技术的难度反而不高,真得不高,一根比较贵的专用柱,有钱谁都能买不要啥技术,一台正常的气相色谱仪,检测器FID就够了,在带入日常工作后,真正难点就是计算问题,组分项目太多,过程很罗嗦,算起来很麻烦,抄结果要累死。所以,在看明白这个方法后,我就萌生了一个想法,编一个计算反式脂肪结果的程序。这个程序与测量不确定度计算的程序的风格类似,也是不仅是为了单纯计算出结果,而是在完成计算的同时,把实验记录也同时完成了。

    (很幸运,最近在清理几个旧电脑时,居然找到了一些以前编的老程序,其中就有这个反式脂肪的程序,由于是给自己用,就没有编什么界面,一切从简。我编程序有两个特点,第一是不喜欢设密码,第二是不喜欢写注释,但是这个程序我却写了很多注释,没办法,参数太多了,不写明白自己肯定很快就蒙了,能不能编完程序都难说)。

   下面附一个数据庫结构和一个子程序,有点程序基础的人,直接COPY略微修改一下就能用了。

 

  

表结构:                       D:\备份资料\程序\FAT\FAT.DBF
数据记录数:                   11      
最近更新的时间:               01/08/08
代码页:                       936    
   字段  字段名                类型                           宽度   小数位      索引   排序         Nulls
      BH                    字符型                                                               
      NAME                  字符型                           10                                     
      CS                    数值型                                                               
      AREA1                 数值型                                                               
      AREA2                 数值型                                                               
      AREA3                 数值型                                                               
      AREA4                 数值型                                                               
      AREA5                 数值型                                                               
      AREA6                 数值型                                                               
     10  AREA7                 数值型                                                               
     11  AREA8                 数值型                                                               
     12  AREA9                 数值型                                                               
     13  AREA10                数值型                                                               
     14  AREA11                数值型                                                               
     15  AREA12                数值型                                                               
     16  AREA13                数值型                                                               
     17  AREA14                数值型                                                               
     18  AREA15                数值型                                                               
     19  AREA16                数值型                                                               
     20  AREA17                数值型                                                               
     21  AREA18                数值型                                                               
     22  AREA19                数值型                                                               
     23  AREA20                数值型                                                               
     24  AREA21                数值型                                                               
     25  AREA22                数值型                                                               
     26  AREA23                数值型                                                               
     27  AREA24                数值型                                                               
     28  AREA25                数值型                                                               
     29  AREA26                数值型                                                               
     30  TF                    数值型                                                            
** 总计 **                                                     240

 

 

 

 

 

 

(3)

      反式脂肪酯计算程序    
                 
注:共26种,由t1 - t26命名
峰面积    由 a1-a26 输入数组a(26)  
换算系数  由 q1-q26 固定数组q(26)
校准因子  由 f1-f26 固定数组f(26)
样品名称   由 n1-n26 固定数组n(26)
**

   饱和脂肪酸 -----------------------------------------------------
*1     C4:0      丁酸(C4:0)甲酯                  Butyric/丁酸 
*2     C6:0      己酸(C6:0)甲酯                  Caproic/己酸 
*3     C8:0      辛酸(C8:0)甲酯                  Octanoic/辛酸
*4     C10:0  奎酸(C10:0)甲酯                 Capric/癸酸 
*5     C12:0  十二烷酸(C12:0)甲酯             Lauric/月桂酸 
*6     C14:0 十四烷酸(C14:0)甲酯              Tetradecanoic/豆寇酸 
*7     C16:0 十六烷酸甲酯,C16:0 棕榈酸甲酯    Palmitic/棕榈酸 
*8     C18:0  十八烷酸甲酯,C18:0 硬酯酸甲酯   Stearic/硬脂酸
*9     C20:0 二十烷酸甲酯,C20:0 花生酸甲酯    Arachidic/花生酸 
*10    C22:0 二十二烷酸甲酯,C22:0山嵛酸甲酯   Behenic/山嵛酸 
*11    C24:0 二十四烷酸甲酯,C24:0 木蜡酸甲酯  Lignoceric/掬焦油酸 

顺式不饱和脂肪酸-----------------------------------------------------------

*12     C16:1 9cis            9-顺-十六碳烯酸甲酯(C16:1 9cis),棕榈油酸甲酯      Palmitoleic/棕榈油酸 
*13     C18:1 6cis            6-顺-十八碳烯酸甲酯(C18:1 6cis),油酸甲酯          Petroselinic/岩芹酸 
*14     C18:1 9cis           9-顺-十八烯酸甲酯,(C18:1 9cis)油酸甲酯           Oleic/油酸 
*15     C18:1 11cis           11-顺-十八烯酸甲酯,(C18:1 11cis)异油酸甲酯         Vaccenic/异油酸 
*16     C18:2 9cis12cis        顺,顺-9,12-十八烯酸甲酯(C18:2 9cis12cis)(亚油酸甲酯)   Linoleic/亚油酸
*17     C18:3 9cis12cis15cis 顺,顺,顺-9,12,15-十八烯酸甲酯(亚麻酸甲酯)(C18:3 9cis12cis15cis)   Linolenic/亚麻酸 
*18     C20:1 11cis          11-顺-二十碳烯酸甲酯,(C20:1 11cis)花生一烯酸甲酯                Eicosenoic/花生烯酸 
*19     C22:1 13cis           13-顺-二十二烯酸甲酯,(C22:1 13cis)顺芥子酸甲酯                  Erucic/芥酸 

反式不饱和脂肪酸

*20     C16:1 9trans        反棕榈烯酸甲酯(C16:1 9trans)              Palmitelaidic/棕榈油酸 
*21     C18:1 6 trans         6-反-十八烯酸甲酯(C18:1 6 trans)       Petroselaidic/岩芹烯酸 
*22     C18:1 9 trans         反油酸甲酯(C18:1 9 trans)                Olaidic/油烯酸 
*23     C18:1 11 trans         异油酸(C18:1 11 trans)                 Transvaccenic/异油酸 
*24     C18:2 9 trans 12 trans  反亚油酸甲酯(C18:2 9trans12trans)       Linolaidic/反-十八碳二烯酸 
*25     C20:1 11 trans        11-反-二十烯酸甲酯(C20:1 11trans)        Trans Eicosenoic/反-花生烯酸
*26     C22:1 13 trans        13-反-二十二烯酸甲酯(C22:1 13trans)     Brassidic/芸苔酸 

 

dimension a(26),q(26),f(26),n(26),x(26)
stor 0 to a
stor 0 to x
f(7)=0.003576
f(1)=1.5106*f(7)
f(2)=1.2836*f(7)
f(3)=1.1702*f(7)
f(4)=1.1021*f(7)
f(5)=1.0567*f(7)
f(6)=1.0243*f(7)
f(7)=1.0000*f(7)
f(8)=0.9811*f(7)
f(9)=0.9660*f(7)
f(10)=0.9536*f(7)
f(11)=0.9433*f(7)
f(12)=0.9925*f(7)
f(13)=0.9745*f(7)
f(14)=0.9745*f(7)
f(15)=0.9745*f(7)
f(16)=0.9678*f(7)
f(17)=0.9612*f(7)
f(18)=0.9600*f(7)
f(19)=0.9482*f(7)
f(20)=0.9925*f(7)
f(21)=0.9745*f(7)
f(22)=0.9745*f(7)
f(23)=0.9745*f(7)
f(24)=0.9678*f(7)
f(25)=0.9600*f(7)
f(26)=0.9482*f(7)

q(1)=0.9868
q(2)=0.9897
q(3)=0.9915
q(4)=0.9928
q(5)=0.9937
q(6)=0.9945
q(7)=0.9950
q(8)=0.9955
q(9)=0.9959
q(10)=0.9962
q(11)=0.9965
q(12)=0.9950
q(13)=0.9955
q(14)=0.9955
q(15)=0.9955
q(16)=0.9954
q(17)=0.9954
q(18)=0.9959
q(19)=0.9962
q(20)=0.9950
q(21)=0.9955
q(22)=0.9955
q(23)=0.9955
q(24)=0.9954
q(25)=0.9959
q(26)=0.9962


n(1)='  C4:0( Butyric/丁酸)      丁酸(C4:0)甲酯 '
n(2)='  C6:0( Caproic/己酸)      己酸(C6:0)甲酯 '
n(3)='  C8:0( Octanoic/辛酸 )    辛酸(C8:0)甲酯 '               
n(4)='  C10:0( Capric/癸酸)      奎酸(C10:0)甲酯 '                
n(5)='  C12:0( Lauric/月桂酸)    十二烷酸(C12:0)甲酯'            
n(6)='  C14:0(Tetradecanoic/豆寇酸)    十四烷酸(C14:0)甲酯'             
n(7)='  C16:0( Palmitic/棕榈酸      十六烷酸甲酯,C16:0 棕榈酸甲酯'
n(8)='  C18:0( Stearic/硬脂酸)         十八烷酸甲酯,C18:0 硬酯酸甲酯'
n(9)='  C20:0( Arachidic/花生酸)       二十烷酸甲酯,C20:0 花生酸甲酯'   
n(10)='  C22:0( Behenic/山嵛酸      二十二烷酸甲酯,C22:0山嵛酸甲酯'
n(11)='  C24:0( Lignoceric/掬焦油酸)   二十四烷酸甲酯,C24:0 木蜡酸甲酯'
n(12)='  C16:1 9cis  ( Palmitoleic/棕榈油酸    9-顺-十六碳烯酸甲酯(C16:1 9cis),棕榈油酸甲酯'   
n(13)='  C18:1 6cis  Petroselinic/岩芹酸    6-顺-十八碳烯酸甲酯(C18:1 6cis),油酸甲酯 '  
n(14)='  C18:1 9cis  ( Oleic/油酸)               9-顺-十八烯酸甲酯,(C18:1 9cis)油酸甲酯'  
n(15)='  C18:1 11cis ( Vaccenic/异油酸)          11-顺-十八烯酸甲酯,(C18:1 11cis)异油酸甲酯'       
n(16)='  C18:2 9cis12cis( Linoleic/亚油酸 )      顺,顺-9,12-十八烯酸甲酯(C18:2 9cis12cis)(亚油酸甲酯)'  
n(17)='  C18:3 9cis12cis15cis( Linolenic/亚麻酸)顺,顺,顺-9,12,15-十八烯酸甲酯(亚麻酸甲酯)(C18:3 9cis12cis15cis)'  
n(18)='  C20:1 11cis  Eicosenoic/花生烯酸)      11-顺-二十碳烯酸甲酯,(C20:1 11cis)花生一烯酸甲酯 '           
n(19)='  C22:1 13cis( Erucic/芥酸)             13-顺-二十二烯酸甲酯,(C22:1 13cis)顺芥子酸甲酯 '           
n(20)='  C16:1 9trans (Palmitelaidic/棕榈油酸   反棕榈烯酸甲酯(C16:1 9trans)'     
n(21)='  C18:1 6 trans( Petroselaidic/岩芹烯酸)  6-反-十八烯酸甲酯(C18:1 6 trans)'    
n(22)='  C18:1 9 trans(  Olaidic/油烯酸 )      反油酸甲酯(C18:1 9 trans) '   
n(23)='  C18:1 11 trans( Transvaccenic/异油酸)  异油酸(C18:1 11 trans)'      
n(24)='  C18:2 9 trans 12 trans(  Linolaidic/反-十八碳二烯酸 )   反亚油酸甲酯(C18:2 9trans12trans)'  
n(25)='  C20:1 11 trans( Trans Eicosenoic/反-花生烯酸)         11-反-二十烯酸甲酯(C20:1 11trans)'    
n(26)='  C22:1 13 trans(  Brassidic/芸苔酸)       13-反-二十二烯酸甲酯(C22:1 13trans)'

set color to n/bg,*r/w,GR*


************************************************
bh1=space(8)

@ 15,16 say replicate('-',55)
@ 13,15 say ' 2--修改旧样品的数据,如名称,峰面积等 '      
@ 16,15 say '样品编号: ' get BH1

 
read
   if bh1=space(8)
   retu
   endif
use fat
locate  for bh=bh1

   if eof()
   retu
   endif

 

@ 20,5 say '修改样品名称: ' get name
 
read

clear
i=1
@ 2+i,40 say '修改以下项目的样品峰面积(未出峰的峰面积为 0 )'

 

   do while i<=26

     store str(INT(i)) to iii
    store 'Area'+Ltrim(iii) to areax
        @ 3+i,5 say '第'+alltrim(str(i))+'种 '+n(i) get &areax
    i=i+1
    enddo
    read
   
   
    ******************   
   i=1
   do while i<=26
 
    store str(INT(i)) to iii
    store 'Area'+Ltrim(iii) to areax
     stor &areax to  a(i)
   
    i=i+1
    enddo
 
 
********************
i=1
af3=0.0000
do while i<=26
af3=af3+a(i)*f(i)
i=i+1
enddo

********************

clear
xq3=0
i=1
@ 1+i,40 say '显示各个项目的脂肪酸甲酯含量'
do while i<=26
x(i)=a(i)*f(i)/af3
@ 3+i,1 say '第'+alltrim(str(i))+'种 '+n(i)+'的脂肪酸甲酯含量: '+str(x(i)*100,6,2)+'%'
xq3=xq3+x(i)*q(i)
i=i+1
enddo
wait

clear

tg=0.0000
i=1
@ 1+i,40 say '显示各个项目的甘油三酯含量'
do while i<=26
@ 3+i,3 say '第'+alltrim(str(i))+'种 '+n(i)+'的甘油三酯的含量  '+str(100*x(i)*q(i)/xq3,6,2)+'%'
  if i>=20
  tg=tg+q(i)*x(i)/xq3
  endif
i=i+1
enddo
 
@ 31,5 say '样品编号: ' + BH
@ 32,5 say '样品名称: '+ name
* @ 33,5 say  '进样次数: '+ str(cs)+'次'
if af3=0
repl tf with 0
@ 34,40 say '总甘油三酯的含量: 0.00%'
else
repl tf with tg*100

@ 34,40 say '总甘油三酯的含量'+str(tg*100,7,3)+'%'

endif

wait ''
return

 

 

(4)

set colo off
set path to d:\程序
zt1="font '宋体',18"
zt2="font '宋体',10"
zt3="font '宋体',13"

 

      反式脂肪酯计算程序    
                 
注:共26种,由t1 - t26命名
峰面积    由 a1-a26 输入数组a(26)  
换算系数  由 q1-q26 固定数组q(26)
校准因子  由 f1-f26 固定数组f(26)
样品名称   由 n1-n26 固定数组n(26)
**

   饱和脂肪酸 -----------------------------------------------------
*1     C4:0      丁酸(C4:0)甲酯                  Butyric/丁酸 
*2     C6:0      己酸(C6:0)甲酯                  Caproic/己酸 
*3     C8:0      辛酸(C8:0)甲酯                  Octanoic/辛酸
*4     C10:0  奎酸(C10:0)甲酯                 Capric/癸酸 
*5     C12:0  十二烷酸(C12:0)甲酯             Lauric/月桂酸 
*6     C14:0 十四烷酸(C14:0)甲酯              Tetradecanoic/豆寇酸 
*7     C16:0 十六烷酸甲酯,C16:0 棕榈酸甲酯    Palmitic/棕榈酸 
*8     C18:0  十八烷酸甲酯,C18:0 硬酯酸甲酯   Stearic/硬脂酸
*9     C20:0 二十烷酸甲酯,C20:0 花生酸甲酯    Arachidic/花生酸 
*10    C22:0 二十二烷酸甲酯,C22:0山嵛酸甲酯   Behenic/山嵛酸 
*11    C24:0 二十四烷酸甲酯,C24:0 木蜡酸甲酯  Lignoceric/掬焦油酸 

顺式不饱和脂肪酸-----------------------------------------------------------

*12     C16:1 9cis            9-顺-十六碳烯酸甲酯(C16:1 9cis),棕榈油酸甲酯      Palmitoleic/棕榈油酸 
*13     C18:1 6cis            6-顺-十八碳烯酸甲酯(C18:1 6cis),油酸甲酯          Petroselinic/岩芹酸 
*14     C18:1 9cis           9-顺-十八烯酸甲酯,(C18:1 9cis)油酸甲酯           Oleic/油酸 
*15     C18:1 11cis           11-顺-十八烯酸甲酯,(C18:1 11cis)异油酸甲酯         Vaccenic/异油酸 
*16     C18:2 9cis12cis        顺,顺-9,12-十八烯酸甲酯(C18:2 9cis12cis)(亚油酸甲酯)   Linoleic/亚油酸
*17     C18:3 9cis12cis15cis 顺,顺,顺-9,12,15-十八烯酸甲酯(亚麻酸甲酯)(C18:3 9cis12cis15cis)   Linolenic/亚麻酸 
*18     C20:1 11cis          11-顺-二十碳烯酸甲酯,(C20:1 11cis)花生一烯酸甲酯                Eicosenoic/花生烯酸 
*19     C22:1 13cis           13-顺-二十二烯酸甲酯,(C22:1 13cis)顺芥子酸甲酯                  Erucic/芥酸 

反式不饱和脂肪酸

*20     C16:1 9trans        反棕榈烯酸甲酯(C16:1 9trans)              Palmitelaidic/棕榈油酸 
*21     C18:1 6 trans         6-反-十八烯酸甲酯(C18:1 6 trans)       Petroselaidic/岩芹烯酸 
*22     C18:1 9 trans         反油酸甲酯(C18:1 9 trans)                Olaidic/油烯酸 
*23     C18:1 11 trans         异油酸(C18:1 11 trans)                 Transvaccenic/异油酸 
*24     C18:2 9 trans 12 trans  反亚油酸甲酯(C18:2 9trans12trans)       Linolaidic/反-十八碳二烯酸 
*25     C20:1 11 trans        11-反-二十烯酸甲酯(C20:1 11trans)        Trans Eicosenoic/反-花生烯酸
*26     C22:1 13 trans        13-反-二十二烯酸甲酯(C22:1 13trans)     Brassidic/芸苔酸 

 

dimension a(26),q(26),f(26),n(26),x(26)
stor 0 to a
stor 0 to x
f(7)=0.003576
* f(7)即f(16:0)的实值
f(1)=1.5106*f(7)
f(2)=1.2836*f(7)
f(3)=1.1702*f(7)
f(4)=1.1021*f(7)
f(5)=1.0567*f(7)
f(6)=1.0243*f(7)
f(7)=1.0000*f(7)
f(8)=0.9811*f(7)
f(9)=0.9660*f(7)
f(10)=0.9536*f(7)
f(11)=0.9433*f(7)
f(12)=0.9925*f(7)
f(13)=0.9745*f(7)
f(14)=0.9745*f(7)
f(15)=0.9745*f(7)
f(16)=0.9678*f(7)
f(17)=0.9612*f(7)
f(18)=0.9600*f(7)
f(19)=0.9482*f(7)
f(20)=0.9925*f(7)
f(21)=0.9745*f(7)
f(22)=0.9745*f(7)
f(23)=0.9745*f(7)
f(24)=0.9678*f(7)
f(25)=0.9600*f(7)
f(26)=0.9482*f(7)

q(1)=0.9868
q(2)=0.9897
q(3)=0.9915
q(4)=0.9928
q(5)=0.9937
q(6)=0.9945
q(7)=0.9950
q(8)=0.9955
q(9)=0.9959
q(10)=0.9962
q(11)=0.9965
q(12)=0.9950
q(13)=0.9955
q(14)=0.9955
q(15)=0.9955
q(16)=0.9954
q(17)=0.9954
q(18)=0.9959
q(19)=0.9962
q(20)=0.9950
q(21)=0.9955
q(22)=0.9955
q(23)=0.9955
q(24)=0.9954
q(25)=0.9959
q(26)=0.9962


n(1)='  C4:0( Butyric/丁酸)      丁酸(C4:0)甲酯 '
n(2)='  C6:0( Caproic/己酸)      己酸(C6:0)甲酯 '
n(3)='  C8:0( Octanoic/辛酸 )    辛酸(C8:0)甲酯 '               
n(4)='  C10:0( Capric/癸酸)      奎酸(C10:0)甲酯 '         

TAG:

学海无涯乐作舟 引用 删除 zzl   /   2013-12-04 09:41:06
5
 

评分:0

我来说两句

显示全部

:loveliness::handshake:victory::funk::time::kiss::call::hug::lol:'(:Q:L;P:$:P:o:@:D:(:)

Open Toolbar