1、ICS 73-010 D 04 道:;中华人民共和国国家标准GB/ T 18340.6- 2001 地质样品有机地化测试汽油族组成质谱分析方法Organic geochemical measurements. of geological samples Method for hydrocarbon types analysis of gasoline by mass spectrometry 2001-03-19发布2001-10-01实施国家质量技术监督局发布GB/ T 18340.6- 2001 前去一口本标准是在地质矿产石油地质海洋地质局企业标准汽油馆分族组成质谱分析方法的基础上，参考石
2、油化工科学研究院质谱法测定汽油的短类组成分析方法以及国内外有关资料，按照先进性、实用性与可靠性相结合的原则，并经过对比测试的实践检验后制定的。本标准确定了汽油族组成的质谱分析方法、定量计算过程和分析精密度要求等。本标准的附录A和附录B都是标准的附录。本标准由国土资源部提出。本标准由全国地质矿产标准化技术委员会岩矿测试标准物质和方法技术委员会归口。本标准起草单位z中国新星石油公司实验地质研究院、上海医药工业研究院、上海染料研究所。本标准主要起草人z周锦南、王惠敏、刘德展。1 范围中华人民共和国国家标准地质样晶有机地化测试汽油族组成质谱分析方法Organic geochemical measure
3、ments of geological samples Method for hydrocarbon types analysis of gasoline by mass spectrometry GB/ T 18340.6- 2001 本标准规定了汽油族组成的质谱分析方法、定量计算过程和分析精密度等。本标准适用于测定轻锢分石油(95%馆出物温度低于196C)中的3种饱和炬及3种芳香炬。试样中烯怪的含量要20 ;:;5 % 3 GB/ T 18340. 6- 2001 附录A(标准的附录汽油族组成计算的FORTRAN源程序A1 主程序GSL4 C 本程序用于用质谱法分析汽油的族组成C 方法的基
4、础是ASTM02789-71 C PROGRA岛GSL C 主程序DIMENSION S(6) ,AS(6) COMMON H (1000) ,HDI(lOOO) CHARACTER铃320EV1，W，TIT铃64，COAT铃11，CTIM铃8CALL RMS(W ,TIT ,COAT ,CTIM) WRITE(祷，60)60 FORMAT(选择输出设备或文件名隐含为打印机J:,) REAO(头(A) )OEV1 IF (OEV1. EQ.川)OEV1=PRN CALL GS5 (DEV1 , W ,TIT ,AC1 ,AC2 ,S ,AS) C 储存和打印分析结果WRITE(饼，500)50
5、0 FORMAT(你要储存结果?Y /NJ , ) READ(势，(Al)ASKIF(ASK. EQ. N七OR.ASK. EQ. n) GOTO 700 NC1=LEN TRIM(W) NC2=LEN TRIM(TIT) OPEN(UNIT=2,FILE = GSL. LST , ACCESS = , APPEND ,STATUS=UN KNOWN) WRITE(2,540) W (1 :NC1) , TITO :NC2) 540 FORMATOX,A , ,A) WRITE(Z ,560) CAS(I) ,1= 1 ,6) 560 FORMATOX,5(F5.2,),F5.2) CLOSE
6、 (UNIT = 2) 700 WRITE( * ,760) 760 FORMATC你要打印几份?(1 - 9) : , ) REAO(祷，780，ERR=700)NP 780 FORMATCI1) 00800 N =l,NP CALL PRTl (DEV1 , W , TIT ,CDAT ,CTIM) CALL PRT3(DEVl ,ACl ,AC2 ,S ,AS) 800 CONTINUE GB/T 18340.6- 2001 1000 STOP 铃铃铃分析结束铃铃铃END A2 子程序GS5SUBROUTINE GS5(DEV1 ,W ,TIT ,AC1 ,AC2 ,S ,AS) C C
7、 本子程序用于计算汽油族组成C 相关描述见ASTMD2789-71 C DIMENSION S(的，CF1(7) ,CF2(7) ,CF3(7) ,CF4(7) ,AS(6) ,AA(6) CHARACTER铃32DEV1，W，TIT铃64COMMON H (1000) ,HDI(1000) DATA CF1/0. 0668 , 0. 0779 , 0. 089 , 0.1001 , 0.1113 , 0.1224 , 0.1335/ , 祷CF2/1.0, 0. 92 ,1. 4 ,1. 8,1. 9, 2. 0, 2.1/, 铃CF3/0.0, 0. 0769 , 0.088, 0. 09
8、91 ,0. 1102 , 0. 1212 ,0.1323/ , 铃CF4/1.0,1. 7, 2.2, 2.4, 2.7, 2.8, 2.9/ C 计算特征峰强组加和43，41，67，77，103，128S (1) = H(43)+H(57) +H(71) +H(85) +H(99) S(2) =H(41)+H(55)+H(69) +H(83) +H(97) S(3) =H(67) + H(68) + H(81) +H(82) + H(95) + H(96) S(4) =H (77) + H(78) +H(79) + H(91)寸H(92)+H(105)+H(106)+H(119)祷十H(1
9、20)+ H (133) + H (134) + H (147) + H (148) +H(l61) + H (162) S(5) =H(103)+H(104) +H(117) + H (118) +H(131)+H(132) 铃+H(145)+H(146) +H(159) +H(160) S(6) =H(128)+H(141) +H(142)十H(155)+H(156) C 屏幕输出各类短特征峰强组加和值300 , SIGMA = o. 0 DO 3201=1,6 320 SIGMA=SIGMA+S(l) WRITE(祷，(/)WRITE(祷，360)(l ,S(I) ，S(l)铃100.O
10、/SIGMA,I=l , 6) 360 FORMAT(3H S(，口，2H)=,F10. 0, 12X ,F7. 2) C 计算链皖及皖基苯平均碳数WRITE(铃，400)400 FORMAT(1HO，20X，6H链烧煌，l1X，6H烧基苯WRITE(赞，420)420 FORMAT(1H，7X，8H平均碳数，6X，4H摩尔，13X，4H摩尔/)F1=0.0 F2= 0.0 F3= 0.0 F4=0.0 UM=O.O DO 4501= 1,7 5 6 11=72+1铃1412=64+1费1413=1+5 GB/ T 18340.6- 2001 U1=H(I1)-CFl(l)铃H(I1-1)IF
11、(I. EQ. 1) U1 = U1 +0.0026铃H(84)一0.014祷(H(92)一0.0769铃H(91)怜一0.008祷(H(106)一0.088铃H(l05)一0.008祷(H(120)一0.0991铃H(119)IF(U1. LT. O. 0) U1 =0. 0 U2=H(l2)一CF3(I)提H(l2-1)IF(U2. LT. o. 0) U2=0. 0 U1=CF2(l)铃U1U2=CF4(l)提U2IF(Ul. GT. UM) UM=U1 WRlTE(铃，440)口，U1，U2440 FORMAT(9X,12 ,2(8X,F9.2) F1=F1十13祷U1F2=F2+U1
12、 F3=F3+13祷U2F4= F4+U2 450 CONTINUE IF(F2. EQ. O. 0) GOTO 510 AC1=F1/F2 510 IF(F4. EQ. O. 0) GOTO 520 AC2=F3/ F4 520 WRITE(祷，530)AC1 ,AC2 530 FORMATo链皖短及烧基苯平均碳数铃=,F5. 2, ,F5. 2) 1F(ACl. LT. 6. O. OR. ACl. GE. 11. 0) GOTO 550 1F(AC2. LT. 6. O. OR. AC2. GE. 11. 0) GOTO 550 GOTO 600 550 WRlTE(祷，560)560
13、FORMAT( $输入链烧炬及烧基苯平均碳数6.-10.J:)READ(祷(F5. 2 ,F5. 2) AC1 ,AC2 C 计算各类炬相对含量600 CALL RMO(AC1 ,S ,AS) CALL RMO(AC2 ,S ,AA) C 含量归一化至100.0，报告体积百分数DO 700 1=4,6 700 AS(I) = AA(I) SA=O.O DO 720 1= 1,6 720 SA=SA+AS(l) DO 740 1= 1,6 740 AS(I) =AS(l)祷100.0/SAGB/ T 18340.6- 2001 1000 RETURN ENO A3 子程序RMOSUBROUTIN
14、E RMO(AC,S,A1) C 本子程序用于计算各类怪相对含量C DIMENSION S(的，A1(的，A2(6),RMC(6.6) IAC=AC一5CALL ROAOAC,6,RMC) DO 6601=1,6 Al(l)=O.O 006201=1,6 620 A10)=Al(l)+RMCO，J)铃SO)IF(Al(l).LT.O.O) A1 (l) = 0.0 660 CONTINUE IF(lAC. EQ. 5) RETURN IAC=IAC+1 CALL ROAOAC,6,RMC) 007001=1,6 A2(l)=0.0 006801=1,6 680 A2(l) = A2(1) +
15、RMC(l ，J)铃SO)IF(A2(l). LT.0.0) A20)=0.0 700 CONTINUE C 对非整型碳数进行插值计算FAC=AC-INT(AC) 008001=1,6 800 Al(1) =A1 (1) + F AC铃(A2(1) - A1 (1) RETURN ENO A4 子程序RDASUBROUTINE ROA(LL ,JO,RMC) C 按平均碳数从磁盘读人相应的逆矩阵DIMENSION RMC(6，的OPEN (UNIT = 1, FILE = GSLRMX. OAT , STATUS = OLO , MOOE = REAO) 005000 L=l,LL 5000
16、REAO(l,5100 , ENO= 5900) (RMC(l.J) ,J=l , JO) .1 = l.JO) 5100 FORMAT (/6(2X,F7.4) o WRITE(8 , 5100) (RMC(!,J) .J= l , JO) ,1= 1 , JO) !炬阵打印5900 CLOSE(UNIT= l) 7 RETURN END GB /T 18340. 6- 2001 A5 子程序RMS8 SUBROUTINE RMS(FIL , TIT ,CDAT ,CTIM) C 子程序RMS用于读取存于磁盘上的质谱数据C LOGICAL EXISTS INTEGER * 2 IB (100
17、0) INTEGER择4IA(1000) ,IC(lOOO) CHARACTER赞32FIL ,FILl ,FIL2 ,BFIL,FIL3 ,FIL4 ,FIL5 ,FIL6 ,SFIL, 势TIT善64，ST1铃128，CDAT祷11，CTIM怜8COMMON H (1000) ,HDI (1000) WRITE(祷，20)20 FORMAT(/lX，输入样品文件名(不含文件扩展名):,) READ(祷，(A)FIL 30 40 LC = LEN TRIM (FIL) FIL1 = FIL (1 : LC) / /飞CTLWRITE(铃，40)FORMAT(输入样品谱范围(起始，终了):,)
18、 READ(铃，(215),ERR= 30) IBG ,IED 读样品谱控制文件OPEN(l ,FILE=FIL1 ,STATUS=OLD) READ(l, (/A) CDAT READ (1, (A) CTIM READ(1, 08,14) lNO , IRl READ(l , (214) M1 , M2 READ (1 , (/A) TIT CLOSE(l) WRITE(祷，60)60 FORMAT(你要减本底谱吗?Y /N : ,) READ(祷，(A1)ANS IF(ANS. EQ. N飞OR.ANS. EQ. n) GOTO 200 WRITE(祷，80)FIL(l :LC) 80
19、FORMAT(使用的本底谱是,A , :4-8，对吗?Y /NJ , ) READ(赞，(A1)ANSl IF(ANSl. EQ. N. OR. ANSl. EQ. n) GOTO 100 FIL4 = FIL(1 :LC)/ / . RA W IR2=IR1 M3 = Ml 如14=如12IBG1 = 4 IED1 = 8 GOTO 180 GB/ T 18340. 6- 2001 100 WRITE(祷，120)120 FORMAT(/输入本底谱文件名:,) READ(铃， (A) BFIL LC1=LEN TRIM(BFIL) FIL3= BFIL(1 :LC1)/ /仁CTLFIL4=
20、 BFIL(1 :LCl)/ /. RAW WRITE(笋，160)160 FORMAT(输入本底谱范围(开始，终了): ., ) READ(祷，(215)IBG1 ,IED1 OPEN(1 ,FILE=FIL3,STATUS=OLD) READ(1, (/ / / 18 ,14) INO ,IR2 READ(1, (汩的M3，M4CLOSE(1) C 、本底谱平均180 CALL SAVER(FIL4, IBGl,IED1 ,M3 , M4;IC) 200 FIL2= FIL(1 :LC)/ / . RAW C 样品谱平均CALL SA VER (FIL2 , IBG ,IED ，岛11，岛
21、12，IA)C 样品谱与本底谱相减DO 300 1=1,758 IA (1) = IA (1) - IC (1) IF(IA(I). LT. 0) IA(I) =0 300 CONTINUE C 减谱存盘320 WRITE(铃，340)340 FORMAT(/输入减谱文件名:,) READ(祷， (A) ) SFIL LC3= LEN TRIM(SFIL) F IL 5=SFIL(1 :LC3)/ /仁CTLFIL6= SFIL(1 :LC3)/ / . RAW C 测试是否有同名文件INQUIRE (FILE = FIL5 ,EXIST = EXISTS) IF (. NOT. EXISTS
22、) GOTO 400 WRITE(势，360)SFIL(1 :LC) 360 FORMAT(文件，A，已存在!覆盖/重命名/跳过?O/ R/ S: , ) READ(祷，(A)ANS IF (ANS. EQ. 0 . OR. ANS. EQ. 0) THEN OPEN (2 ,FILE= FIL5 ,STATUS= OLD) CLOSE(2,STATUS=DELETE ) OPEN (2 ,FILE=FIL6,STATUS=OLD ,ERR=400) CLOSE(2 ,STATUS= DELETE) ELSE IF(ANS. EQ. R. OR. ANS. EQ. r) THEN GOTO 3
23、20 9 ELSE GOTO 500 END lF C 产生减谱控制文件GB/ T 18340. 6- 2001 400 OPEN Cl .FILE= FILl.STATUS=OLD) OPEN(2.FILE=FIL5.STATUS= NEW) D04201=1.5 READCl. (A) ST1 LC2=LEN TRlt-.在(ST1)420 WRITE(2 ,(A) ST1 Cl :LC2) READ(1,(l4) N1 WRITE(2，气11)1 READ(l,(A) ST1 DO 4401=128，1.一1lF CST1 (l;I). NE. AND. ST1(I:l). NE. CH
24、AR(O) GOTO 460 440 CONTINUE 460 LC2 = I TIT=ST1 cl :LC2) WRITE(2 .480) ST1 cl :LC2) .FILCl :LC) .IBG.IED.BFILCl :LC1) .IBG1.IEDl 480 FORMA T (A，勺，. A. ( .14. .14. )一A.( .14. .14.) WRITE(2. (11) 1 CLOSECl) CLOSE(2) C 产生减谱数据文件IB=IA OPEN(UNIT=2.FILE=FIL6,STATUS= NEW .FORM=BINARY) WRITE(2) (IB(J) .J= M1
25、.M2) CLOSE (UNIT = 2) 500 H=IA RETURN END A6 子程序SAVER10 SUBROUTINE SAVER(FIL2 ,IBG ,IED.M1 ,M2.IL) C 读取多幅质谱并将它们平均INTEGER美21AClOOO)lNTEGER祷4lLCl OOO) .14 CHARACTER祷32FIL2 IL=O IR1=M2Ml+1 OPENCUNIT=2,FILE=FIL2,STATUS= OLD .FORM=BINARY 铃RECL=IR1铃2.ACCESS = , DIRECT ) DO 220 I=IBG,IED G/ T 18340.6- 2001
26、 READ(2,REC=D (lA (J) , = M1 ,M2) DO 220 J =M1 ,M2 I4=IA() IF(l4. LT. 0) 14=14+65536 220 IL()=IL()+I4 260 CLOSE(UNIT=2) IED=I- 1 IL= IL/ (IED- IBG + l) RETURN END A7 子程序PRT1SUBROUTINE PRT1 (DEV2, W1 ,TIT1 ,CDA T1 ,CTIM1) C 本子程序用于打印分析报告头部CHARACTER铃32W1 ,DEV2, TIT1铃64，CDAT1於11，CTIM1铃8OPEN(8,FILE=DEV2,
27、STATUS= NEW ) WRITE(8,1800) W1 ,CDAT1 ,CTIM1,TITl 1800 FORMAT(/ /13H分析文件名:，A，T33，6H日期:,A ,2X ,A/ 铃llH样品说明:,A/) CLOSE(8) RETURN END A8 子程序PRT3SUBROUTINE PRT3(DEV2,AC1 ,AC2 ,S ,Al) C 本子程序用于打印汽油族组成分析结果C DIMENSI0N A1 (的，S(6)CHARACTER DEV2怜32，LINE铃51，ST祷108，LC祷5DATA LINE/叫一一一一-一/DATA STO:18)/Paraffins /
28、DATA ST09:36)/ Monocycloparaffins / DATA ST(37:54)/ Dicycloparaffins / DATA STC55:72)/ Alkyl benzenes / DATA STC73:90)/ lndans/ Tetralines / DATA ST(91:108) / Naphthalenes /，LC/一一斗/RL=O.O DO 19001=1,6 IF(Al(l). GT. RL) RL= Al(l) 1900 CONTINUE ML = RL/10+1 OPEN(6,FILE= DEV2 ,STATUS=NEW) WRITE(6j2100)
29、 (I祷20，I=0，ML/2)11 12 GB /T 18340. 6-2001 2100 FORMAT(/lHO，31H一一一一一一一一一/铃4X，11HHydrocarbon, 5X , 3Hlon, 3X , 6HVolume/ 铃7X,4Htype , 7X , 6Hlnten. ,1X , 7Hpercent ,6(12 ,8X)/ ) WRITE(6,2120) (LC ,I= l ,ML) 2120 FOR孔1AT(/lX，32H一-一一一一-1，10A)D02200I= 1 ,6 LA=A1(1)头0.5+1.5 2200 WRITE(6,2300) ST(18铃。一1)+1
30、:18铃。一1)+1肘，S(l),Al(I) ,LlNE (1 :LA) 2300 FORMAT Cl HO,A ,F8. 0,F 5. 1 ,A ) WRITE(6, 2400) AC1 ,AC2 2400 FORMAT(/lX，31H-一一一一一-/祷Note:saturate average carbon number: ,F 5. 2/ 祷，aromatic average carbon number: ,F 5. 2/) RETURN END GB/T 18340.6-2001 附景B (标准的附录汽油族组成计算用逆矩阵系鼓(GSLRMX. DAT) B1 C6-C10汽油锢分逆矩阵
31、系鼓0.9016 -0.1368 0.0257 一0.00030.0000 0. 0000 一0.44711. 0285 一0.2391一0.00020.0000 0. 0000 0.0100 一0.02580.4325 0. 0000 0.0000 0. 0000 0.0017 一0.0048-0.0149 0.5117 0.0000 0. 0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 O. 72注1一0.06550.0105 一0.0100一0.01000.
32、0000 一0.25420.7283 一0.1695一0.0051一0.00350.0000 0.0167 一0.05230.4387 0.0001 0.0003 0. 0000 0.0010 一0.0044一0.01340.4576 一0.08970. 0000 0.0000 0.0000 一0.00020.0000 0.5424 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 O. 6449 一0.05840.0090 一0.0011一0.0105一0.0082一0.19020.6132 一O.1428 一0.0063一0.0029O. 0
33、006 0.0128 一0.04690.4375 0.0001 0.0003 一0.00040.0007 一0.0049一0.01250.4375 一0.0857一0.02710.0000 0. 0002 0.0004 一0.02070.5465 一0.00260.0000 0.0000 0.0000 0.0000 0.0000 0.5757 0.6043 一0.06730.0071 一0.0018一0.0095一0.0075一0.19330.6183 一0.1929一0.0130一0.00170. 0011 0.0212 一0.08220.6809 0.0003 0.0004 一0.0006
34、0.0007 一0.0040-0.0261 0.4015 一0.0787一0.02480.0001 0.0002 0.0020 一0.03610.5496 一0.0016一0.00900.0008 0.0000 0.0000 0.0001 0.5759 0. 5766 一0.15620.0606 0.0001 一0.0025一0.0070一O.1897 0.7443 一0.3315一0.0270一0.00040.0015 0.0666 -0.2792 O. 7592 0.0087 一0.0032一0.0009一0.00060.0021 一0.0201O. 3903 一0.1240一0.0238
35、0.0002 一0.00010.0029 一O.0709 0.7315 一0.0007一0.01200.0033 一0.0012一0.0006一0.0174O. 5761 13 FCON- .O叮伺俨H阁。国华人民共和国家标准地质样晶有机地化测试汽油族组成质谱分析方法GB/T 18340. 6-2001 中 中国标准出版社出版北京复兴门外三里河北街16号邮政编码:100045电话:6852394668517548 中国标准出版社秦皇岛印刷厂印刷新华书店北京发行所发行各地新华书店经售4陪开本880X 1230 1/16 印张174字数26千字2001年8月第一版2001年8月第一次印刷印数1-1000定价12.00元* 书号:155066.1-17763 网址版权专有侵权必究举报电话:(010)68533533* 578-544 科目GB/T 18340.6-2001
