GB T 18340.3-2001 地质样品有机地化测试 重馏分石油饱和烃族组成质谱分析方法.pdf

1、ICS 73. 010 D 04 G昌中华人民共和国国家标准GB/ T 18340.3- 2001 地质样品有机地化测试重馆分石油饱和怪族组成质谱分析方法Organic geochemical measurements of geological samples Method for hydrocarbon types analysis of gas-oil saturates fractions by mass spectrometry 2001-03-19发布2001-10-01实施国家质量技术监督局发布GB / T 18340. 3- 2001 前本标准是在地质矿产石油地质海洋地质局企业

2、标准重情分油饱和蛙族组成质谱分析方法的基础上，参考石油化工科学研究院质谱法测定重锢分油饱和煌的短类组成分析方法以及国内外有关资料，按照先进性、实用性与可靠性相结合的原则，并经过对比测试的实践检验后制定的。本标准确定了重馆分石油饱和短族组成的质谱分析方法、族组成定量计算过程和分析精密度要求。本标准的附录A和附录B都是标准的附录。本标准出国土资源部提出。本标准由全国地质矿产标准化技术委员会岩矿测试标准物质和方法技术委员会归口。本标准起草单位:中国新星石泊公司实验地质研究院、上海医药工业研究院、上海染料研究所。本标准主要起草人:周锦南、王惠敏、刘德晨。1 范围中华人民共和国国家标准地质样品有机地化测

3、试重锢分石油饱和短族组成质谱分析方法Organic geochemical measurements of geological samples 肌1ethodfor hydrocarbon types analysis of gas-oil saturates fractions by mass spectrometry GB/ T 18340. 3-2001 本标准规定了重馆分石油饱和短族组成的质谱分析方法、定量计算过程和分析精密度等。本标准造用于重锢分石油(沸程范围300C540C)饱和怪锢分中的7种饱和短及1种芳香煌的短族组成的定量分析。试样中应不含烯怪，芳怪含量也要20 5 4 %

4、GB/ T 18340. 3-2001 附录A(标准的附录重锢分石油饱和姬族组成计算的FORTRAN源程序A1 主程序SHAC 本方法用于测定饱和石油馆分中七种饱和短族及一种芳香短族的组成C PROGRAM SHA C如1AINPROGRA孔4LOGICAL EXISTS DIMENSION S (10) , YNO(日，YIO(日，XO(5),A(8) ,A2(8) ,RMC(8 , 8, 2) ,RA(2) CHARACTER铃32FIL，FILl，DEV1，DEV2，TIT养64，CDAT提11.CTIM秘8COMMON H (1000) ,HDI(1000) DATA YNO/. 34

5、7.606,1. 250 , 2. 439 , 4.000/ , YIO/. 0364 ,. 0505 , .0735. 铃.1061，.138/.XO/16. , 20. , 24. , 28. .32. / C 从磁盘读取质谱，进行平均和减谱处理，将干净的样品谱存盘CALL RMSCFIL ,TIT ,CDAT ,CTIM) C 校正重同位素的影响CALL DEISO C 计算.z:;71 , .z:; 69. .z:; 109 , .z:; 149. .z:; 189. .z:; 229和2:269J=O L = 71 K=4 CALL SU如1(J.L.K.S) L=69 K=6 CA

6、LL SUM(J,L.K ,S) L=109 K=7 CALL SUM(J ,L.K.S) L=149 K=8 CALL SUM (J .L.K,S) L = 189 K =9 CALL SUM (J .L,K ,S) L=229 K=lO CALL SU岛1(J, L ,K ,S) L=269 K=l1 5 6 GB/ T 18340. 3- 2001 CALL SUM(J , L ,K ,S) C 计算2.;91之1L=91 K = 5 CALL SU肌if(J, L , K ,S) C 计算L;91之2L=117 K=4 CALL SUM(J ,L ,K ,S) C 计算2.;91之3L

7、=129 K =4 CALL SU岛1(J,L ,K ,S) C 计算完整的2.;91 S(8) =S(8) +S(9) +S(10) SlGMA=O.O 003801=1,8 380 SIGMA=SlGMA+S(I) WRITEC祷，400)400 FORMAT(/) WRITE(养，410)口，S(1)，S(l)铃100.0/SIGMA，I=1，8)410 FORMATC3H S( , 11 , 2H ) = ,F12. 2, 8X ,F10. 4) C 寻找最强的链烧短(一环环皖)峰WRITE(祷，440)440 FORMATOHO，T9，碳数，T21，链烧峰强，T39，一环环烧峰强，)

8、 PMAX=O.O PBMAX=O.O 00 500 1=16,32 Pl =HDl(2+ 14关1)PIB=HDl(14祷1)WR1TE(祷，450刀，凹，P1B450 FORMATOX ,T10,I2,T 20 ,F 9. 2,T40,F9. 2) 1F (P孔1AX.GT. P1 )GOTO 460 PMAX=P1 11 = 1 460 IF(PBMAX. GT. PIB)GOTO 500 PB沁1AX=P1BIlB= 1 500 CONTINUE 1FCP孔1AX.LE. O. o. ANO. PB岛1AX.LE. O. O)GOTO 555 WRITE(餐，530)口，I1B530

9、FORMATCO由链烧(一环环烧)分子峰确定的平均碳数，13，C ,12 , ) C 另一种确定平均碳数方法，即摩尔数平均法SC=O.O SCB=16铃HDI(224)SI=O.O SIB=HDI(224) IAC=O IACB=O DO 5401=16,35 P1=HDI(14铃1+2)P1B=HDI(14铃1)SC=SC+I铃PlSCB=SCB-I铃PlBSI=SI+Pl 540 SIB=SIB+P1B GB/ T 18340. 3- 2001 IF(SI. NE. o. 0) IAC=SC/ SI +0. 5 IF(SIB. NE. o. 0) IACB=SCB/SIB+O. 5 WRI

10、TE(祷，550)IAC, IACB 550 FORMAT(由摩尔法确定的平均碳数，,13 , ( 12 , ) ) 555 WRITE(笋，560)IAC 560 FORMAT(输入平均碳数(16-32)，默认值为，, 13 ,: , ) READ(祷，570)IACl 570 FORMATCI2) IF(lAC1. GE. 16. AND. IAC1. LE. 32) IAC= IAC1 C 计算2;69/2;71和127峰/226峰比值IF(lAC. GT. 16) GOTO 680 WRITE(祷，580)580 FORMAT(。你要计算2:69/2:71和127峰/226峰比值?Y

11、/ N: , ) READ(祷， (A 1) , ) ANS IF(ANS. NE. Y. AND. ANS. NE. y) GOTO 680 RA (1) =INTC100赞S(2)/S(I)+0.5)/100.0RA(2)=INT(10祷HDI(127)/HDI(226)+0.5)/ 10.0 WRITE( * ,600) RA 600 FORMAT(/4X, 2: 69 / 2: 71= ,F5. 2，T41，127峰/226峰=，F4.1)FIL1 =NCI6 610 NF=LEN TRIM(FIL1) FIL1 =FIL1C1 :NF)/ /. TXT C 测试是否已有老文件INQU

12、IRE(FILE=FIL1 ,EXIST=EXISTS) IF(. NOT. EXISTS) GOTO 660 WRITE(费，620)FIL1 620 FORMAT(/已有，A，老文件。覆盖/重命名/退出O/R/E:, ) READ(铃，(A)ANS IF(ANS. EQ. 0 . OR. ANS. EQ. 0) THEN OPEN(2,FILE=FIL1 ,STATUS=OLD) 7 8 GB/T 18340. 3- 2001 CLOSE(2 ,STATUS= DELETE) ELSE IF(ANS. EQ. R. OR. ANS. EQ. r) THEN WRITE(祷，640)640

13、FORMA T (输入新文件名(不含扩展名):,) READ(铸，(A)FIL1 GOTO 610 ELSE STO END 1F 660 CALL PRTN16(F1Ll,F1L ,T1T ,CDAT ,CT1M ,RA) WR1TE(祷，670)F1L1 670 FORMATo正十六镜校正值已存入，A，文件中0) GOTO 1500 C 灵敏度因于插值计算680 XX=1AC CALL IPL(XO ,YNO , 5,XX ,YN) CALL 1PL (XO , YIO , 5, XX , Yl) C 计算SET值，确定链烧怪的正异构程度、按比例混合使用两个逆矩阵SET=O.O !当AM=

14、O时作异构烧处理AM=YN保HDI(14铃1AC)+Y1铃HDI(14铃1AC-29)IF(AM. NE. O. 0) SET= YN铃HDI(14铃IAC)/AMWRITE(祷，700)SET 700 FORMAT(O正异构因子SET=,F7. 3) C 计算各类短含量12 = IAC- 15 CALL RDA(8, 8, 2,I2 ,RMC) D0720 1=1 ,8 !取用正烧短逆矩阵计算各类;怪A(l)=O.O DO 710 J=1 ,8 710 A(1)=A(l)+S(J)祷RMC(1，J,1) 1F(A(l). LT. O. 0) A(I)=O.。720 CONT1NUE DO 8

15、20 1=1,8 !取用异烧短逆矩阵计算各类怪A2(l)=0.。DO 810 J =1,8 810 A2(l)=A2(l)+S(J)铃RMC(I，J，2)1F(A2(l). LT. O. 0) A2(l) = 0.0 820 CONT1NUE DO 900 1=1 , 8 !按SET值将两个矩阵计算结果混合900 A (I) =A(I)铃SET+A2(l)祷(1.。一SET)C 将各类煌含量归一化至100.0并报告各类短体积百分含量SA= O.O DO 11001=1 ,8 GB/ T 18340.3- 2001 1100 SA=SA+A(I) DO 1200 1=1,8 1200 A(I)

16、=A(I)祷100.0/SAC 储存分析结果WRITE( * ,1220) 1220 FORMAT(你要保存分析结果吗? Y /NJ ,) READ(铃，(A1)ANS1F(ANS. EQ. N. OR. ANS. EQ. n ) GOTO 1280 LC = LEN TR1M (FIL) LC2 = LEN TRIM(TIT) OPEN(UNIT=2 , FILE=SH5. LST , ACCESS = APPEND , STATUS = UN KNOWN) WRITE(2, 1240) F1L (1 :LC) , TIT(1 :LC2) , IAC 1240 FORMAT(1X ,A, ,

17、A , C ,12) WRITE(2, 1260) (A(I) ，1=1，的1260 FORMAT(1X , 7(F5. 2, , ),F 5. 2) CLOSE (UNIT = 2) C 打印分析结果1280 DEV1=PRN WRITE(铃，1290)1290 FORMAT(O选择输出设备或文件名铃不含扩展名。若是打印机，直接回车J:,) READ(善，(A)DEV2IF(DEV2. EQ. ) GOTO 1400 DEV1=DEV2 NF= LEN TRIM(DEV1) DEV1 = DEV1 (1 :NF)/ / . TXT C 测试是否已有老文件INQUIRE (F1LE = DEV

18、1 ,EXIST. = EX1STS) 1F(. NOT. EXISTS) GOTO 1400 WRITE(势，1300)DEV1 1300 FORMAT(/已有，A，老文件。覆盖/重命名/退出O/R/EJ:, ) READ(祷，(A)ANS 1F(ANS. EQ. 0 . OR. ANS. EQ. 0) THEN OPEN(2,F1LE= DEV1 ,STATUS= OLD) CLOSE(2，STAUS=DELETE) ELSE IF(ANS. EQ. R. OR. ANS. EQ. r) THEN GOTO 1280 ELSE STOP ENDIF 1400 WRITE仆，1420)142

19、0 FORMAT(你要打印几份?(1一的，默认值为1J:,) 9 GB/ T 18340. 3-2001 READ(畴(11) NP IF(NP. EQ. 0) NP=l DO 1480 N = l ,NP CALL PRT2(DEV1 ,FIL ,TIT ,CDAT ,CTIM , IAC ,SET , A) 1480 CONTINUE 1500 STOP赞誉赘重馆分石油饱和怪族组成分析完毕铃祷祷END A2 子程序PRT210 SUBROUTINE PRT2(DEV2 , W1 , TIT1 ,CDAT1 ,CTIM1 , 13 ,SET1 , A1) C 本子程序用于饱和怪族组成分析结果

20、的打印C DIMENSION A1 (8) , B(7) , LINE(51) , LC(5) CHARACTER赞32W1 , DEV2 , TIT1铃64，F2120保12，CDAT1赞11，CTIM1赞8DATA LINE/ 1H 1，50赞1H-/, LC/ 1H 1，4美1H-/OPEN(8,FILE=DEV2 ,STATUS= UNKNOWN) WRITE(8 , 1800) W1 ,CDAT1 ,CTIM1 , TIT1 1800 FORMAT(/ /分析文件名:，A，T33，日期:, A , 2X , A/ 畴样品说明:,A / ) RL=O.O DO 19001= 1 ,8

21、IF(A1 (1). GT. RL) RL=Al(l) 1900 CONTINUE ML=RL/10十lML1 = ML/ 10 ML2=ML一孔1Ll势10WRITEC8 ,2000) 13 ,SET1 2000 FORMAT(/平均碳数:，13，T33，正异构因子:，F5.3)WRITEC8 , 2100) (1祷20，I=0，ML/2)2100 FORMAT(/2X，23H-一一一一/祷4X，11HHydrocarbon, 4X , 6HVolume/ 铃7X, 4Htype , 7X , 7Hpercent ;6(1 3, 7X) WRITE怡，2110)2110 FORMAT(3X，

22、23H一一一一一-， ) F2120= ( / /CHAR(MLl + 48) / / CHAR(ML2+48)/ / (5A1) , ) WRITE(8 , F2120) (LC() , = 1，日，I=1.ML)WRITE(8, (1H 1/ ) DO 22001= 1,7 LA=Al(l)铃0.5+1.5 2200 WRITE(8 , 2300) 1 - 1 , Al(l) , (LINE() , = 1 , LA) 2300 FORMAT(6X,Il, 1X , 4HRing , 7X , F5. 1, 2X , 51A1) LA=A1(8)铃0.5+1.5 WRITE(8 , 240

23、0)A1 (肘，(LINE () , = 1 , LA) GB/ T 18340. 3- 2001 2400 FORMAT(3X, 12HMonoaromatic , 4X ,F5. 1, 2X , 51A1 ) C 除去单环芳短井重新归一化SA1=0.0 00 25001=1 ,7 2500 SA1=SA1+Al(l) 00 26001=1 ,7 2600 B(l)=Al(l)祷100.0/SA11F(A1(8). LT. O. 1) GOTO 2700 RL=O.。00 26101=1,7 1F(B(l). GT. RL) RL=B(l) 2610 CONTINUE ML=RL/ 10+1

24、 ML1=ML/ 10 ML2=ML-MLl祷10WR1TE(8,2650) 2650 FORMAT(/ /除去单环芳怪并重新归一化，/) WR1TE(8, 2100) (l提20，1=0，ML/2)WR1TE(8,2110) F2120=气，/CHAR(ML1 +48)/CHAR(ML2+ 48)/ (5A1 ) , ) WR1TE(8, F2120) (LC(J)，=l，日.1=l.ML)WR1TE(8. (1H 1/) DO 2680 1=1 . 7 LA=B(I)祷0.5+1.5 2680 WR1TE(8, 2300) 1- 1,B(I).(LINE(J), = 1.LA) C 计算样

25、品中H/C原子数比值2700 HC=2+(2铃B(1)-2铃B(3)-4祷B(4)-6赞B(5)-8祷B(6)-10铃B(7)/13/100WR1TE(8.2800) HC 2800 FORMAT(/本饱和怪中H/C=.F7. 3) 1F(A1(8). LE. 5. 0) GOTO 3100 WRITE(8,3000) 3000 FORMA T (/20X.请注意门，/ / 3X. 养本结果可能超差，这是因为试样中含有超过5.0%的单环芳炬。，/) 3100 WRITE(8.3200) 3200 FORMAT(/5X.铸铃铸分析报告结束餐善铃/)CLOSE (8) RETURN ENO A3

26、子程序SUMSUBROUT1NE SUM(J1 ,L1 .K1 .S1) C 子程序SUM用于计算峰强加和COMMON H (1000) .HDI(1000) 11 DIMENSION Sl(lO) 1=1+1 Sl(l ) = O.O DO 5000 l=l ,Kl GB/ T 18340. 3- 2001 Sl (1) =Sl (1) + HDI (L1) 5000 L1=L1+14 RETURN END A4 子程序IPLSUBROUTINE IPL(XO ,YO ,N ,X,Y) C 子程序IPL用于插值计算DIMENSIO N XO(N) , YO(N) Y=O.O DO 6300

27、I = l ,N P = l. O DO 6200 = l ,N IF(l- )6100 ,6200 ,6100 6100 P=P祷(X- XO (J) ) / (XO(l) - XO (J) 6200 CONTINUE 6300 Y=Y+P铃YO(l)RETURN END A5 子程序RDASUBROUTINE RDA (lI, ,KK ,LL,RMC1 ) C 本子程序用于读取存盘数组DlMENSION RMC1 (II，刀，KK)OPEN(UNIT=l ,FILE= CMI. DAT ,STATUS=OLD ,MODE=READ) DO 7000 N = l,LL 7000 READ(l

28、, 7100) ( ( (RMC1 (l , ,K) , = 1 ,J) ,1= 1 ,II) ,K= 1 ,KK) 7100 FORMAT(2(8(/8(2X ,F7.4) D7200 WRITE(7, 7100) ( (RMCICI, ,K) , = 1,) ,1= 1 ,II) ,K = 1 ,KK) CLOSE (UNIT = 1) RETURN END A6 子程序PRT1612 SUBROUTINE PRTN16(FIL2,W1,TIT1 ,CDAT1 ,CTIM1 ,RA) C 本子程序用于将正十六烧的两个比值存入文件DlMENSION RA(2) CHARACTER祷32FIL

29、2 ,W1 ，TIT1铃64，CDAT1铃11，CTIM1铃8OPEN (8 ,FILE= FIL2 ,status = new) G/ T . 18340. 3- 2001 WRITE饵，1700)W1 ,CDAT1 ,CTIM1 ,TIT1 1700 FORMAT(j /分析文件名:，A，T33，6H日期:，A，2X，A/祷样品说明:,A / ) WRITEC8.l720) RA 1720 FORMAT(j/ 4X , L: 69 / L: 71= ,F5. 2，T41，127峰/226峰=,F4. 1) WRITE(8,1740) 1740 FORMAT(j /) CLOSE (8) E

30、ND A7 子程序RMSSUBROUTINE RMS(FIL , TIT ,CDAT ,CTIM) C 子程序RMS用于读取存于磁盘上的质谱数据，气LOGICAL EXISTS C INTEGER提2IB(lOOO) INTEGER铃4IA(1000) ,IC(lOOO) CHARACTER铃32FIL ,FILl,FIL2 ,BFIL ,FIL3 ,FIL4 ,FIL5 ,FIL6 ,SFIL , 祷TIT铃64，ST1铃128，CDAT铃11，CTIM铸8COMMON H Cl OOO) ,HDICl OOO) WRITEC * ,20) 20 . FORMAT(j /lX，输入样品文件名

31、(不含文件扩展名):,) READC善，CA)FIL 30 40 LC= LEN TRIM CFIL) FIL1 = FILCl :LC)/ /飞CTLWRITE( * ,40) FORMAT(输入样品谱范围(起始，终了):, ) READ(祷，(215), ERR=30) IBG ,IED 读样品谱控制文件OPEN Cl ,FILE= FIL1 ,STATUS= OLD) READ(l,(jA) CDAT READCl,(A) ) CTIM READCl, 08,14) )INO , IR1 READCl, (214)M1，M2 READ(l, (j A) ) TIT CLOSE(l) WR

32、ITE(铃，60)60 FORMAT(你要减本底谱吗?Y /NJ尸，)READ(链，(A1)ANS IF(ANS. EQ. N. OR. ANS; EQ.-亘在下GOTO200 WRITE(铃，80)FIL(l :LC)、I .i 80 FORMAT(使用的本底谱是/f,A , : 4.8飞对吗?Y /NJ ,) READ(祷，(A1)ANSl. IF(ANSl. EQ. N. OR. ANSl. EQ. n) GOT 100 13 14 GB/ T 18340. 3-2001 FIL4=FILCl :LC)/ /. RAW IR2=IR1 M3=M1 M4=M2 IBG1 =4 IED1=8

33、 GOTO 180 100 WRITE(势.120)120 FORMAT(/输入本底谱文件名: ) READ(畴，(A)BFIL LC1=LEN TRI扣1BFIL)FIL3 = BFIL(1 :LC1)/ / . CTL FIL4=BFILCl :LC1)/ /. RAW WRITE(祷.160)160 FORMAT输入本底谱范围(开始，终了): ) READ(祷.(215) IBG1.IED1 OPEN(1 ,FILE=FIL3.STATUS=OLD) READ(1, (/ / /18.14) INO.IR2 READ(1 . (214) M3.M4 CLOSE(l ) C 本底谱平均18

34、0 CALL SAVERCFIL4.IBG1 .IED1.M3,M4 .IC) 200 FIL2= FIL(1 :LC)/ /. RAW C 样品谱平均CALL SAVERCFIL2.IBG.IED.M1.M2.IA) C 样品谱与本底谱相减DO 300 1=1.758 IA(l) =IA(l) - ICCD IF (lA C 1). L T . 0) IA C 1) = 0 300 CONTINUE C 减谱存盘320 WRITE(祷.340)340 FORMAT(/输入减谱文件名: ) READ(祷， (A) , ) SFIL LC3=LEN TRIM CSFIU FIL5 = SFILC

35、l :LC3)/ /. CTL FIL6=SFIL(1 :LC3)/ / . RA W C 测试是否有同名文件INQUIRE(FILE=FIL5.EXIST=EXISTS) IF(. NOT . EXISTS) GOTO 400 WRITE仆，360)SFILCl :LC) 360 FORMA T ( ，文件，A，巳存在!覆盖/重命名/跳过?O/R/S:. ) READ(祷.(A)ANS GB/ T 18340. 3- 2001 IFCANS. EQ. 0. OR. ANS. EQ. 0) THEN OPEN(2 ,FILE=FIL5,STATUS=OLD ) CLOSEC2 ,STATUS=

36、DELETE ) OPENC2 ,FILE=FIL6,STATUS=OLD ,ERR =400) CLOSEC2 ,STATUS= DELETE) ELSE IFCANS. EQ. R. OR. ANS. EQ.气)THEN GOTO 320 ELSE GOTO 500 ENDIF C 产生减谱控制文件400 OPEN (1 ,FILE=FILl ,ST ATUS= OLD ) OPEN (2 ,FILE=FIL5 ,ST ATUS= NEW) D04201=1,5 READC1 (A) ) ST1 LC2= LEN TRIMCSTl) 420 WRITEC2.,CA) STI C1 :LC2

37、) READ(l, (14) ) N1 WRITE(2. (U) 1 READC1, (A) ST1 DO 4401=128，1，一1IF (ST1 (1 : 1). NE.门.AND. ST1CI:1). NE. CHARCO) GOTO 460 440 CONTINUE 460 LC2=I TIT=STl (1 :LC2) WRITEC2 ,480) ST1 (1 :LC2) ,F IL C1 :LC) ,IBG ,IED.BFILCl :LC口，IBG1.IED1480 FORMATCA, ,A, ( .14, .14.)一，,A , C ,14, . .14, ) ) WRITEC2,

38、 (1口)1 CLOSE(l) CLOSE(2) C 产生减谱数据文件IB=IA OPENCUNIT=2.FILE=FIL6,STATUS=NEW .FORM=, BINARY ) WRITE(2) (lB(J) ,=Ml,M2) CLOSE CUNIT = 2) 500 H=IA RETURN END A8 子程序SAVERSUBROUTINE SAVERCFIL2 ,IBG ,IED,M1.M2 ,IL) C 读取多幅质谱并将它们平均INTEGER赞2IA(1000)15 GB/ T 18340. 3- 2001 INTEGER铃4IL (1000) ,14 CHARACTER铃32FIL

39、2 IL=O IR1=M2- M1+ 1 OPENCUNIT=2,FILE=FIL2,STATUS=OLD ,FORM=BINARY , 铃RECL=IR1祷2，ACCESS=DIRECT) DO 220 I = IBG ,IED READ(2 ,REC = l) CIA (J) ,J = M1 ,M2) DO 220 J=M1 ,M2 14 = IA(J) IF(I4. LT. 0) 14=14十65536220 IL(J)=IL(J)+I4 260 CLOSE (UNIT= 2) IED=I-1 IL = IL/ (IED- IBG+ 1) RETURN END A9 子程序DEISO16

40、 SUBROUTINE DEISO C C 本子程序计算单同位素峰强，假设所有离子的Z值在+2至一11之间C 分子通式为C(N)H(2N+Z): C COMMON H (lOOO) ,HDI(1000) DIMENSION NCARB(1000) ,NHYD(lOOO) DO 10 1 = 12 , 758 NCARB(I)= 0 10 NHYD(l)=O DO 20 K=12, 758 NCARB(K) = (K十11)/14NHYD(K) = K -12铃NCARB(K)IF(NHYD(K). LT. 0) NHYD(K) = 0 20 CONTINUE DO 30 K =14,758 H

41、DI (K) = H (K) - HDI (K -1)铃(.010811铃FLOATCNCARB(K-1)+.00015赞赞FLOAT(NHYD(K-1)HDI(K)二HDI(K)+HDI(K-2)势已00005844祷FLOAT(NCARB(K一2)祷势(l-NCARB(K -2)十.1125E-7祷FLOATCNHYD (K -2)赞(1-NHYDCK一2) 铃-.162165E-5传FLOAT(NCARB(K-2)势NHYD(K一2)IF(HDI(K). LT. O. 0) HDI(K) = O. 0 30 CONTINUE RETURN END B1 GB/T 18340. 3- 20

42、01 附录B(标准的附录)饱和坦族组成计算用道矩阵系组(CMI. DAT) C16C32饱和炬E异构造矩阵系数O. 5344 一0.0292一O.0066 0.0215 0. 0299 O. 0000 一O.0610 0.3404 一O.2146 一O.1162 一0.0362O. 0000 一O.0039 0.0170 O. 8491 一0.6968一O.3420 0.0000 O. 0000 一O.0004 0.0115 1. 7220 一1.3545 O. 0000 O. 0001 O. 0004 0.0039 一O.0138 3.2594 O. 0000 O. 0000 O. 0000

43、 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 -0. 0007 一O.0029 一0.0237一O.1566 一O.3494 O. 0000 0.6543 一0.0358-0.0081 O. 0264 0.0366 O. 0000 一O.0866 0.3416 一0.2143一O.1171 一0.0377O. 0000 一O.0053 0.0172 0.8492 -0.6968 一0.3420O. 0000 O. 0001 一O.0004 0.0115 1. 7220 -1. 3

44、545 O. 0000 O. 0000 O. 0004 0.0039 一0.01383. 2594 0.0000 O. 0000 0.0000 O. 0000 O. 0000 0.0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 0.0001 一0.0029一O.0237 一O.1565 -0.3493 O. 0000 O. 5243 一0.0311一O.0075 O. 02 27 0. 0322 O. 0000 一0.06600.3403 一0.2130一0.1164-0.0385 O. 0000 一0.00380. 0

45、154 O. 8375 一0.6826一0.3318O. 0000 0.0000 一O.0004 O. 0095 1. 6824 一1.3111 O. 0000 0.0001 O. 0004 O. 0039 一0.01473.1247 O. 0000 0.0000 O. 0000 O. 0000 O. 0000 0. 0000 O. 0000 0.0000 O. 0000 O. 0000 O. 0000 O. 0000 O. 0000 -0. 0007 一0.0027一O.0220 一O.1514 一0.3331O. 0000 0.6435 -0. 0382 一O.0092 O. 0279 0

46、. 0395 O. 0000 一0.0942O. 3418 一0.2125一0.1176一0.0403O. 0000 一0.00540. 0155 O. 8375 一0.6826一0.3319O. 0000 0.0000 一0.0002O. 0090 1. 6825 一1.3111 O. 0000 0. 0000 0. 0004 O. 0040 一0.01473. 1247 O. 0000 0. 0000 0. 0000 O. 0000 O. 0000 O. 0000 O. 0000 0.0000 O. 0000 O. 0000 O. 0000 0. 0000 O. 0000 0. 0000

47、一O.0027 -0. 0220 -0. 1514 一0.3331O. 0000 0. 5175 一0.0338一O.0085 O. 0234 O. 0344 O. 0000 O. 0000 一0.01510.0000 一0.01120.0000 一0.00480.0000 0.0152 0.0000 一0.0485O. 0000 0.0000 O. 0000 0.0000 0.0000 O. 3521 O. 0000 一0.01850.0000 一0.01010.0000 -0. 0046 0. 0000 0.0152 O. 0000 -0.0485 O. 0000 0.0000 O. 00

48、00 0.0000 O. 0000 O. 3521 O. 0000 一0.0163O. 0000 一0.0121O. 0000 一0.0052O. 0000 0.0166 O. 0000 一0.05270.0000 0.0000 O. 0000 0.0000 O. 0000 O. 3612 O. 0000 一0.0200O. 0000 一0.0112O. 0000 一O.0052 O. 0000 0.0166 O. 0000 一O.0527 O. 0000 0.0000 O. 0000 O. 0000 0.0000 O. 3612 0.0000 一0.017817 GB/ T 18340.3- 2001 一O.0720 - 0. 0039 0.3403一0.2091-0.1183一0.04040.0138 0.8183一O.6626 -