一份色譜實驗記錄(二十）

上一篇 / 下一篇 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
      1 BH                    字符型                            8                                      否
      2 NAME                  字符型                           10                                      否
      3 CS                    數值型                            8                                      否
      4 AREA1                 數值型                            8                                      否
      5 AREA2                 數值型                            8                                      否
      6 AREA3                 數值型                            8                                      否
      7 AREA4                 數值型                            8                                      否
      8 AREA5                 數值型                            8                                      否
      9 AREA6                 數值型                            8                                      否
     10 AREA7                 數值型                            8                                      否
     11 AREA8                 數值型                            8                                      否
     12 AREA9                 數值型                            8                                      否
     13 AREA10                數值型                            8                                      否
     14 AREA11                數值型                            8                                      否
     15 AREA12                數值型                            8                                      否
     16 AREA13                數值型                            8                                      否
     17 AREA14                數值型                            8                                      否
     18 AREA15                數值型                            8                                      否
     19 AREA16                數值型                            8                                      否
     20 AREA17                數值型                            8                                      否
     21 AREA18                數值型                            8                                      否
     22 AREA19                數值型                            8                                      否
     23 AREA20                數值型                            8                                      否
     24 AREA21                數值型                            8                                      否
     25 AREA22                數值型                            8                                      否
     26 AREA23                數值型                            8                                      否
     27 AREA24                數值型                            8                                      否
     28 AREA25                數值型                            8                                      否
     29 AREA26                數值型                            8                                      否
     30 TF                    數值型                            5        2                             否
** 總計 **                                                     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"

* 順式不飽和脂肪酸-----------------------------------------------------------

* 反式不飽和脂肪酸

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)

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)甲酯 '