GB T 18340.4-2001 地质样品有机地化测试 重馏份石油芳香烃族组成质谱谱分析方法.PDF

《GB T 18340.4-2001 地质样品有机地化测试 重馏份石油芳香烃族组成质谱谱分析方法.PDF》由会员分享，可在线阅读，更多相关《GB T 18340.4-2001 地质样品有机地化测试 重馏份石油芳香烃族组成质谱谱分析方法.PDF（26页珍藏版）》请在麦多课文档分享上搜索。

1、GB/T 18340.4-2001前台二习本标准是在地质矿产石油地质海洋地质局企业标准重馏分油芳香烃族组成质谱分析方法的基础上，参考石油化工科学研究院质谱法测定重馏分油芳香烃的烃类组成分析方法以及国内外有关资料，按照先进性、实用性与可靠性相结合的原则，并经过对比测试的实践检验后制定的。本标准确定了重馏分石油芳香烃族组成的质谱分析方法、定量计算过程和分析精密度要求。本标准的附录A是标准的附录。本标准由国土资源部提出。本标准由全国地质矿产标准化技术委员会岩矿测试标准物质和方法技术委员会归口。本标准起草单位：中国新星石油公司实验地质研究院、上海医药工业研究院、上海染料研究所。本标准主要起草人：周锦南

2、、王惠敏、刘德晨。中华人民共和国国家标准地质样品有机地化测试重馏分石油芳香烃族组成质谱分析方法Organic geochemical measurements of geological samples-Method for aromatic types analysis of gas-oils aromaticfractions by mass spectrometrycB/T 18340.4-2001范围本标准规定了重馏分石油芳香烃族组成的质谱分析方法、定量计算过程和分析精密度等。本标准适用于重馏分石油（沸程范围300C -540 C）芳香烃馏分中的18种芳香烃及3种芳香唾吩的烃族组成的定

3、量分析。试样中必须不含烯烃，总硫量必须低于总重量的100，非芳烃含量应50oa注：在总计21种烃族中，只有15种已定名，其余6种为相应于各自同系物的未鉴定芳烃。2方法提要被测油样不经分离直接作质谱测定，得到该样品的混合质谱图。根据质谱数据算出I VI类芳烃的特征峰组强度加和，每一特征峰组包括多同位素分子离子峰系列和单同位素分子离子峰减1系列。用根据纯芳烃化合物导出的7X7逆阵计算出7类芳烃的相对含量。又根据基线法把每类烃分为标称型(0型）、第一重叠型(1型）和第二重叠型(2型）。最后求出7类3型共计21种芳香烃族的相对含量，烃族含量以体积百分数计。3仪器及设备3.1质谱计：扇形磁场质谱计或磁场

4、一静电场双聚集质谱计。3.2进样系统：任何能升温至125C-350 C、无样品损失、无污染的样品汽化导人装置。可以是储样罐，也可以是直接蒸发进样器，用后者并配合高速数据系统可做微量分析。3.3数据处理系统：能对谱图作数据实时采集、相加、相减、列表、作图等处理的质谱数据系统，带FORTRAN 77, MS-FORTRAN V5.0,MS-FORTRAN POWERSTATION V4.0或其他更高版本的FORTRAN程序编译器。3.4微量注射器：1 )UL, 10 aI，和100 uI.o4试剂及材料全氟煤油（PFK）或其他已知准确质量的参考物：质量数校准用。正十六烷：仪器调节和校准用，色谱纯。

5、氯仿：溶剂，样品稀释和实验器具清洗用，分析纯。J卫0亡Q；月目月月分析条件5.1质量分辨率：800a国家质It技术监督局2001一03一19批准2001一10一01实施GB/T 18340.4-20015.25.3电子轰击能量：70 eV o电离室温度：250Co进样器加热温度：125 C -350C o质量扫描方式：指数扫描。质量扫描范围：60u-760u,电子倍增器增益：）100离子通道真空度：蕊1 X 10-Pa.月卜叹U，1八0亡Jl阮J工J亡J一J6校准和其他准备工作6.1各种仪器校准和分析工作均应在仪器运转达到稳定时进行。6.2用PFK（或其他已知准确质量的参考物）校准数据系统的质

6、量标尺，使在m/z60u-760u范围内的质量偏差不超过士0. 2u，否则应建立新的质量校正表。6.3用正十六烷校准质谱裂解模型，使Z69/Z71比值在0.18士0. 03之间，m/zl27与m/z226的峰高比在1.4士0.1之间。若比值超过此范围，应重新调整离子源内离子推斥极电压等相关的工作参数。艺69与习71的定义为：名69=m/z69-I-m/z83-m/z97+m/zll1+m/zl25+m/zl39的总峰强；习71=m/z71+m/z85+m/z99+m/zll3的总峰强。上述各质谱峰强是指已作了重同位素校正的“单同位素”峰强。6.4计算芳香烃族组成用的FORTRAN程序载于本标准

7、的附录A（标准的附录），它由一个主程序（附录A1和四个子程序（附录A2-A5）组成。由于质谱数据储存方式和计算机系统兼容性的缘故，允许使用者对程序中的某些语句作必要的修改。在进样分析前，要用相应的编译程序对其进行编译和链接，使之成为一个完整的可执行程序。注：在本标准附录A（标准的附录）所提供的程序里，质谱数据存取格式是按中国科学院科学仪器中心的软件产品KYKY-GC/MS-DSI质谱数据系统的数据存盘格式编写的。如果用户所使用的数据系统与此不同，则应对RMS. FOR（附录A4)和SAVER. FOR（附录A5)这两个子程序中的质谱数据输人输出语句进行相应的修改。分析步骤7.1使用精馏切割法获

8、取沸点为300 C-540 C的重馏分石油，并用柱层析法分离得到芳香烃馏分。7.2检查质谱计各部分的工作状态，确定它们都在第6节规定的分析条件下运行且是稳定的。7.3样品质谱分析对于储样罐与直接蒸发两种不同的进样方式，取用不同的分析步骤。7.3.1储样罐进样方式7.31.1启动质谱扫描，在m/z60u-760u范围内，记录下约10幅本底离子流质谱。73.1.2关闭储样罐抽气阀，用微量注射器吸取3 pL-5 t.1.试样注人罐内，被汽化的试样通过漏孔或漏阀导人离子源。注意控制进人离子源的样品流量，保持离子流放大器工作在线性区域。7.3.1.3在不间断地采集质谱的同时，观察总离子流的变化。当总离子

9、流达极大后，持续采集10幅以上的样品质谱，然后结束样品谱的采集。7.3.1.4打开储样罐抽气阀，用泵抽去被测试样。待残余本底质谱总离子流低于上次样品总离子流的1/50后。可准备作下一次进样。7.3.2直接蒸发进样方式7.3.2.1用氯仿作溶剂，将样品稀释至10 tg/jA1.-100 Ag/l.l.左右的浓度。用1 til，的微量注射器吸取0. 1 FAL-1 tAL溶液注人直接进样用专用柑祸，待溶剂自然挥发干后，样品柑涡内将留有约5 )g-20 jig量的芳香烃试样。GB/T 18340.4-20017.3.2.2将柑祸装到直接进样杆上，再把进样杆装到离子源的直接进样口上。打开前级真空抽气阀

10、，将进样杆预抽至10 Pa量级的前级真空。7.3-2.3启动质谱扫描，在m/z60u-760u范围内，记录下约10幅本底离子流质谱。7.3-2.4关闭前级真空抽气阀，打开离子源高真空隔离阀，把进样杆推人离子源，位于进样杆顶端的盛样柑祸则进人电离室。加热柑祸使样品均衡地蒸发，直至样品全部蒸发完。在此期间质谱采集不间断地进行，记录下全部样品质谱。采集的样品谱总数以30-100幅为最佳。注：若质谱计带有直接蒸发进样总离子流控制装置，则使用该装置完成样品燕发过程。总离子流控制装置可保证样品均衡地蒸发，并保持离子流放大器工作在线性区域。7.3-2.5将进样杆拉出离子源，关闭高真空隔离阀，让高真空泵不断地

11、抽去被测试样品在离子源内的残余物。待本底质谱总离子流低于上次样品总离子流的1乃0后，可准备下一次进样。8计算8.1选择样品质谱范围烃族组成计算用的样品质谱范围是从总离子流图上来选择的。对于储样罐与直接蒸发两种不同的进样方式，使用不同的选择方法。8.1.，储样罐进样方式选择总离子流强度较高而且较稳定的8-20幅样品质谱，但不要选用质谱采集时间过长、甚至总离子流已明显下跌的那些谱。注：不选用质谱采集时间过长的那些谱，是因为在进样系统中为控制样品流量使用了分子漏孔，这将导致产生质量歧视效应和同位素分馏效应，造成谱图的变形。这二种效应对计算结果的不利影响将随时间的增加而逐渐加剧。8.1.2直接蒸发进样

12、方式选择样品受热蒸发总离子流开始上升的谱作为样品起始谱，选择总离子流下降到接近本底的谱作为样品终了谱。8.2选择进样前的仪器本底离子流质谱，总数5-8幅即可。8.3用计算机计算芳香烃烃族组成调用由附录A（标准的附录）所产生的计算机可执行程序计算芳香烃烃族组成，该程序对储样罐与直接蒸发两种进样方式均适用。按程序的提示依次输人相关信息并加载样品质谱数据，计算机即开始进行下列各项计算。8.3.1对选定范围的样品谱与本底谱分别进行平均计算。8.3.2将样品平均谱与本底平均谱相减，求得扣除f本底的真实的样品谱。8.3.3开始计算芳香烃样品中21种烃族的体积百分含量。这21种烃族的类别、类型、名称以及分子

13、通式见表to有关芳香烃烃族组成计算步骤的详细说明，参见附录A（标准的附录）中的多条程序注释。9精密度同一样品两次平行分析，其测定结果的允许双差D见表注：双差D的定义为：DA一B式中：A第一次（基本分析）测量值；B一一第二次（检查分析）测量值。GB/T 18340.4-2001表1重馏分石油芳香烃的类别、类型、名称以及分子通式类别类型名称分子式I0烷基苯CH2。1苯并唾吩CH2.,oS2环烷菲C, H2,-20皿0环烷苯C,H2” s1花类C,H2,z22未鉴定皿0二环烷苯CHz”一101类C, H2.-Z42未鉴定N0禁类CH2 i21二苯并唾吩CH2 -,6S2未鉴定V0危及二苯并吠喃C,H

14、2a；及CH2，。（）1类CM 2,282未鉴定研0药类CH2-is1二苯并葱CH2,rgo2未鉴定IN0菲类CH2， -is1蔡苯并唆吩CH2， -22S2未鉴定表2重馏分石油芳香烃族组成质谱平行测定的允许双差D质量分数范围允许双差525.20320G5GB/T 18340.4-2001附录A（标准的附录）计算重馏分石油芳香烃族组成用的FORTRAN源程序A1主程序ATAPROGRAM ATA本程序是计算重馏分石油芳香烃族组成用的主程序程序中的H(M）代表质量M的多同位素峰高程序中的HDI (M）代表质量M的单同位素峰高CCCCCDIMENSION MS(8),HT(8)CHARACTER，

15、32 FIL,TIT二64,CDAT * 11,CTIM * 8COMMON H(1000),HDI(1000)质谱的读入、平均、相减和储存CALL RMS(FIL,TIT,CDAT,CTIM)校正重同位素影响CALL DEISO计算和报告芳香烃族组成CALL AROMTC(FIL,TIT,CDAT,CTIM)STOPENDA2子程序DEISOSUBROUTINE DEISO本子程序计算单同位素峰强，假设所有离子的Z值在2至一11之间分子通式为C(N)H(2N十Z).COMMON H(1000),HDI(1000)DIMENSION NCARB(1000),NHYD(1000)DO 10 I二

16、12,758NCARB(D=010 NHYD(I)=0DO 20 K=12,758NCARB (K )(K11)/14NHYD(K）二K一12，NCARB(K)IF(NHYD(K). LT. 0) NHYD(K）020 CONTINUEDO 30K14，758HDI (K）=H(K）一HDI(K一1）二（010811，FI.OAT(NCARB(K一1），00015二，FI.OAT (NHYD(K一I）CCCCGB/T 18340.4-200130HDI ( K）二HDI ( K）十HDI(K一2）二（00005844二FI-OAT(NCARB(K一2)(I一NCARB(K一2)）十1125E一

17、7FLOAT(NHYD(K一2)（1一NHYD(K一2)）一162165E一5FI.OAT(NCARB(K一2)，NHYD(K一2)）IF (HDI(K). LT. 0. 0) HDI (K）二0.0CONTINUERETURNENDA3子程序AROMTCSUBROUTINE AROMTC(W,TITLE,CDATI,CTIM1)本子程序给出了石油重芳烃馏分计算过程详细的计算过程描述见Anal Chem 41,1548-54(1969).CCCCCOMMON H(1000),HDI(1000)DIMENSION AIN(7，7),BA(7),BB(7),SR(1000)CHARACTER AN

18、S1，W*32，DEV，32，F25018，TITLECTIMI二8DATA AIN/-I-1. 8094，一1952，十0124，一0027，一0015，一0一1601，十2. 0479，一2806，一0401，十0082,-1-.一0943，一2287，十2. 3024，一4936，一0601，一1一0292，十0033，一0580,+1.9404，一1337，一一0022，一0033，一0026，一0195，十1. 9773，一一0420，十0026，一0018，一0151，一，0584，十2.164，CDAT111，0028，0000，0089，0043，0123，4193，9904/

19、NC 1=LENNC2=I_EN一2346，一1069，一，0267，一0019，一0057，一TRIM(W)TRIM(TITLE）CC“开方”数组初始化C2132DO 2132 I12，750SR (I）0.0加和芳烃馏分相关谱峰以用于7X7矩阵处理A6,A7,A8.的值先是代表Z值为6,7,8.的峰强加和然后在A6,A7,A8.中又加人了相应于分子离子峰一1系列的奇数Z值峰强A6二A6+A7，A8=A8十A9.CCCCCC2106A6二0. 0DO 2106 M78，750，14A6=A6+H(M）A70. 0GB/T 18340.4-200121072108210921102111211

20、221132114211521162117DO 2107 M91，750,14A7二A7十HDI (M )M=M十A7IF(A6.I_T. 1.0) A6=1.0A8二0.0DO 2108 M104，750，14A8=A8+H(M）A9=0. 0DO 2109M=117，750，14A9=A9十HDI(M)A8=A8+A9IF(A8. LT. 1. 0) A8二1.0A10一0.0DO 2110 M二130,750,14A10=A10+H(M）All0.0DO 2111 M129，750，14All二All+HDI(M）A10二A10十AllIF(Al0. LT. l. 0) A10二1.0A

21、12=0. 0DO 2112 M二128,750,14A12二A12+H(M)A13=0. 0DO 2113 M=141，750，14A13=Al3+HDI(M)A12一A12+A13IF (A12. LT. 1. 0) A12二1.0A14=0. 0DO 2114M=154，750，14A14=A14+H(M)A15=0. 0DO 2115 M167，750,14A15二A15十HDI (M)A14=Al4+Al5IF(A14. LT. 1. 0) A14=1. 0A16=0. 0DO 2116 M=166，750，14A16=Alb+H(M)A17二0.0DO 2117 M二179,750

22、,14A17=A17+HDI(M）A16=A16十A17Gs/T 18340.4-200121182119IF(A16.I_T. 1.0) A16=1.0A18=0. 0DO 2118 M178，750,14A18二A18+H(M)A19=0. 0DO 2119 M191，750，14A19=A19十HDI (M）A18=A18+A19IF(A18. LT. 1. 0) A18=1. 0对质量数175,176,189,190,200,213处的不恰当的离子峰强作校正校正HDI(175)CCCCCCD1175=HDI (161）一（HDI (161）一HDI(203)/3.0IF(HDI（175

23、). GE. CDI175）GO TO 1046CD1175二HDI（175)CC校正H(176)C1046 CH176=H(162）一（H(162）一H(204)/3.0IF(H(176). GE. CH176）GO TO 1048CH176=H(176)CC校正HDI(189)C1048 CD1189=CD1175一（CDI175一HDI (203) /2. 0IF (HDI (189). GE. CDI189) GO TO 1049CDI189二HDI(189)CC校正H(190)C1049 CH190=CH196一（CH196一H(204)/2.0IF (H (190). GE. CH

24、 190) GO TO 2101CH190二H(190)CC校正H(200)C2101 CH200=(H(186)H(214)/2.0IF (H (200). GE. CH200) GO TO 2102CH200=H(200)CC校正HDI(213)GB/r 18340.4-20012102CD1213二（HDI(199)+HDI(227)/2.0IF(HDI(213). GE. CDI213) GO TO 2103CDI213二HDI(213)下述语句校正A6,A10加和2103 A6=A6一（HDI(175)）+HDI(189)+H(176)+H(190)）二十CDI175+CDI189+

25、CH176+CH190A10=A10一（H(200)+HDI(213)）+CH200十CDI213CC把A6,A8.，加和值重定义为数组元素变量，然后以芳烃逆矩阵AIN(I,J)乘之CCC21242125212621272128BA (1）=A6BA (2）二A8BA(3)=A10BA(4)=A12BA(5)=A14BA(6)=A16BA(7)=A18DO 2125）1，7BB(J)=0. 0DO 2124 I=1，7BB (J）=BB(J）+BA(1)、AIN(J,I)CONTINUEDO 2127）1，7IF(BB(J)）2126，2127，2127BB(J)=0.0CONTINUEAA6

26、=BB(l）AA8=BB(2)AA10=BB(3)AA12=BB(4)AA14二BB (5)AA16=BB(6)AA18=BB(7)SUMAA=0. 0DO 2128）1，7SUMAA二SUMAA+BB(J）计算出的AA6,AA8.即相应于Z=-6,-8.等类芳烃的分量大小SUMAA是AA值的加和，它代表所计算芳烃的总量下面用基线法求解Z=-6重叠型CCCCCGB/T 18340.4-2001213021312133A7二A7一HDI(175）一HDI (189) +CDI175+CDII89HDI (175) =CDI175HDI(189)=CD1189DO 2130M105，750，14I

27、F (HDI (M）2130，2131，2130CONTINUEMM=M一14SLOPE= ( ( (0. 72，HDI(105)），二0.5）一（HDI (MM )）二，0.5)/(90.71一（1000. 0/FLOAT(MM)），二2)B= (0. 72，HDI (105）0.5一90.71SLOPEDO 2133 M二147,MM，14REALM=MSR (M) =SLOPE二（1000. 0/REALM)，2+B以上是Z二一6重叠型，烷基苯峰强的开方储存于数组SR (I）中以下是解Z=-8重叠型CCCC213421352136DO 2134 M二215，750，14IF (HDI (

28、M）2134，2135，2134CONTINUEMN二M一14SLOPE= (0. 66HDI(173)X0.5）一（HDI(MN)）二0-5)/(34.12一（1000. 0/FLOAT (MN )），二2)B二（0. 66、HDI(173) * * 0. 5一34. 12SLOPEDO 2136 M215 ,MN，14REALM=MSR (M）=SLOPE，(1000. 0/REALM）2十B以下是解Z一10重叠型CCC213721382139All=All一HDI(213）十CDI213HDI (213) =CDI213DO 2137 M241，750，14IF(HDI(M)）2137，

29、21382137CONTINUEMO=M一14SLOPE二（(HDI(185)）二0.5一（HDI(MO)）二0.5)/(29. 22一（1000. 0/FLOAT (MO )），二2)B=HDI (185 )二0. 5一29. 22SLOPEDO 2139 M241, MO，14REALM=MSR (M）=SLOPE（1000. 0/REALM）2+BcB/T 18340.4-2001以下是解z一12重叠型CC214021412142DO 2140 M197,750，14IF(HDI(M)）2140，2141，2140CONTINUEMP=M一14SLOPE= ( ( (0. 25HDI(1

30、83)）二0.5）一（HDI (MP) )、0.5)/(29.86一（1000. 0/FLOAT(MP)）二，2)B二（0.25二HDI(183）二二0.5一29.86二SLOPEDO 2142 M197,MP，14REALM=MSR (M)=SLOPE(1000. 0/REALM)，2+B以下是解z二一14重叠型CCC214321442145DO 2143 M=265，750，14IF (HDI (M)）2143,2144，2143CONTINUEMQ=M一14SLOPE二（(0.64二HDI(251)，0.5）一（HDI(MQ)）二二0.5)/（15.87一（1000. 0/FLOAT (

31、MQ )）2)B二（0.64HDI (251)0.5一15.87SLOPEDO 2145 M二265 , MQ，14REAI,M=MSR (M)=SLOPE二（1000. 0/REALM )，2+B以下是解Z一16重叠型CCC214621472148DO 2146 M=291，750，14IF (HDI (M)）2146，2147，2146CONTINUEMR二M一14SLOPE =（(0.7HDI(277)），0. 5）一（HDI (MR)），0-5)/(13.03一（1000. o/FI.OAT(MR)）二二2)B= (0. 7，HDI(277) * * 0. 5一13. 03 * SLO

32、PEDO 2148 M二291,MR，14REALM二MSR(M)=SLOPE(1000. 0/REALM），2+B以下是解Z一18重叠型CCCDO 2149 M247，750,14IF (HDI(M)）2149,2150，2149GB/T 18340.4-20012149 CONTINUE2150 MS二M一14SLOPE= (0. 58，HDI(233)）二二0. 5）一（HDI (MS) )，0.5)/，(18.42一（1000. 0/FLOAT(MS)）二，2)B=(0. 58 * HDI (233)）0.5一18.42SLOPEDO 2151 M二247,MS,14REALM二M21

33、51 SR(M）二SLOPE(1000. 0/REALM )二2+B至此已计算了开方数组对某些谱在零峰高区仍得到斜率和截距。如发生这一情况，将会把错误带人SR数组以下语句使HDI=O. 0时SR也为零CCCCC21522153DO 2153 I=12，750IF (HDI (I)）2152，2152，2153SR (I）二0.0CONTINUESR数组平方，得到各个Z值的未校准标称型峰高CCCDO 2154 I二12，7502154 SR(I)=SR(I) * * 2CC对某些非线性开方关系的SR (I）值进行校准CSR(147)=SR(147)，1.44SR(197）二SR(197)3. 1

34、0SR (211）二SR(211），2.52SR(225）二SR(225) * 2. 07SR (239) =SR (239) * 1.83SR(253)=SR(253)1.59SR(267）二SR(267)，1.39SR (281）=SR(281）1.28SR(295)=SR(295) * 1. 26SR(309)=SR(309)1.14SR(323）二SR(323)1. 06SR(265)=SR(265）二1.42SR (279）二SR(279) * 1. 24SR(293)=SR(293)，1.12SR (307) =SR(307)1.06SR(291）二SR(291）二1. 24SR(

35、305)=SR(305) * 1. 15SR(319)=SR(319)1.07GB/T 18340.4-2001SR(333）二SR(333) * 1.06SR(347)=SR(347）二1.05SR(361）=SR(361）* l. 03SR(247）二SR(247)1.61SR(261）一SR(261），1.50SR(275)=SR(275)，1.44SR(289)=SR(289)1.37SR (303) =SR (303）二1.28SR(317)=SR(317) * 1. 28SR(331）=SR(331），l. 21SR(345) =SR(345）二1.10SR(359)=SR(359

36、)，1.09SR(373)=SR(373)1.07SR(387)=SR(387)，1.05CC如发现有SR(M）值超过对应的HDI ( M）则进行校正C21552156DO 2156 M二12,750IF (SR (M）一HDI(M) 2156,2156,2155SR (M)=HDI(M)CONTINUE计算A7的各个部分即A6A,AIOS,A20A各个类型CCC21572158215921602161A6A=0. 0DO 2157 M=91，133，14A6A=A6A十HDI (M)DO 2158 M=147,MM，14A6A=A6A+SR(M)AIOS=0. 0DO 2159 M二147，

37、189，14AIDS=AlOS十HDI ( M）一SR(M)AIOS=AIOS/0. 75A20A=A7一A6A一A 1 OSIF(A20A) 2160,2161,2161A 20A二0.0A l OS二A7一A6ACONTINUE计算分量A6A,AIOS和A20ACCCTRASH二（A6一AA6 * 0. 5579)IF(TRASH. LT. 0. 0) TRASH=O.(A7/A6)0GB/T 18340.4-2001A7二A7一TRASHA6A二A6A一TRASHIF(A6A. LT. 0. 0) A6A=0. 0IF(A6A. EQ. 0. 0) A7二A l OS十A20AIF (A

38、7. EQ. 0. 0) THENA6A=0. 0AIOS=0. 0A20A=0. 0ELSEA6A二（A6A/A7）二AA6Al0S= (AlOS/A7）二AA6A20A二（A20A/A7)AA6END IFcc计算A9的各个部分即A8A,A22A和A36AcA8A=0. 0DO 2162 M二117，201，142162 A8A=A8A十HDI(M)DO 2163 M215,MN，142163 A8A=A8A十SR(M）A22A二0.0DO 2164 M215，257，142164 A22A二A22A+HDI (M）一SR(M)A22A二A22A/0. 75A36A=A9一A8A一A22A

39、IF(A36A) 2165，2166，21662165 A36A=0.0A22A=A9一A8A2166 CONTINUECC计算分量A8A,A22A和A36ACTRASH=(A8一AA80.4997）二（A9/A8)IF(TRASH. LT. 0. 0) TRASH二0.0A9二A9一TRASHA8A=A8A一TRASHIF(A8A. LT. 0. 0) A8A0.0IF(A8A. EQ. 0. 0) A9二A22A+A36AIF (A9. EQ. 0. 0) THENA8A0.0A22A二0.0A36A=0. 0Gs/T 18340.4-2001ELSEA8A=(A8A/A9)AA8A22A

40、=(A22A/A9)，AA8A36A(A36A/A9)，AA8END IFCC计算All的各个部分即A10A,A24A和A38ACA10A=0. 0DO 2167 M129,227，142167 A10A二AIOA+HDI(M)DO 2168 M二241,MO,142168 Al0A二Al0A十SR(M)A24A=0. 0DO 2169 M二241，283，142169 A24A=A24A十HDI(M）一SR(M)A24A二A24A/0. 75A38A二All一A 10A一A24AIF(A38A) 2170,2171，21712170 A38A二0.0A24A=All一Al0A2171 CON

41、TINUECC计算分量AlOA,A24A和A38ACTRASH二（A10一AA10，0.4435）二（All/A10)IF (TRASH. LT. 0. 0) TRASH二0.0All=All一TRASHAIOA二A 10A一TRASHIF(AlOA. LT. 0. 0) A10A=0. 0IF(AlOA. EQ. 0.0) All=A24A十A38AIF (A11. EQ. 0. 0) THENAIOA二0.0A24A=0. 0A38A二0.0ELSEA 10A二（AlOA/All）二AA10A24A=(A24A/All），AA10A38A=(A38A/All）AA10END IFCC计算A

42、13的各个部分即A12A,A16S和A26ACGs/T 18340.4-200121722173217421752176A12A=0. 0DO 2172 M141，183，14A12A=Al2A+HDI (M）DO 2173 M=197,MR，14A12A=Al2A十SR(M）A16S=0. 0DO 2174 M197，225，14A 16S二A16S十HDI ( M）一SR (M）A16S=Al6S/0. 625A26A=Al3一A12A一A16SIF(A26A) 2175,2176，2176A26A二0.0A16S=AI3一A12ACONTINUE计算分量A12A,A16S和A26ACCC

43、TRASH=(A12一AA120. 5192)(A13/A12)IF (TRASH. LT. 0. 0) TRASH= 0. 0A13=A13一TRASHA12A=Al2A一TRASHIF(A12A. LT. 0. 0) A12A=0. 0IF(A12A. EQ. 0. 0) A13一A16S+A26AIF (A13. EQ. 0. 0) THENA12A=0. 0A16S=0. 0A26A=0. 0ELSEA12A= (A12A/A13）二AA12A16S= (A16S/A13)，AA12A26A=(A26A/Al3)，AA12END IFCC计算A15的各个部分即A14A,A28A和A42

44、ACA14A=0. 0DO 2177 M=167，251，142177 A14A=Al4A+HDI(M)DO 2178 M=265,MQ，142178 A14A=Al4A+SR(M)A28A=0. 0DO 2179 M二265，307，142179 A28A二A28A+HDI(M）一SR (M)GB/T 18340.4-200121802181A28A=A28A/0. 75A42A=A15一A14A一A28AIF(A42A) 2180,2181，2181A42A二0.0A28A=Al5一A14ACONTINUE计算分量A14A,A28A和A42ACCCTRASH= (A14一AA14二0.50

45、75)(A15/A14)IF (TRASH. LT. 0. 0) TRASH= 0. 0A15=A15一TRASHA14A=Al4A一TRASHIF(A14A. LT. 0. 0) A14A二0,0IF(A14A. EQ. 0. 0) A15二A28A+A42AIF (A15. EQ. 0. 0) THENA14A=o. 0A28A二0.0A42A=0. 0ELSEA14A二（A14A/A15)AA14A28A二（A28A/A15）二AA14A42A=(A42A/A15)，AA14END IFCC计算A17的各个部分即A16A,A30A和A44ACA16A一0.0DO 2182 M二179，2

46、77，142182 A16A=Al6A+HDI(M）DO 2183 M291,MR，142183 A16A=Al6A一SR(M）A30A=0. 0DO 2184 M二291，333，142184 A30A二A30A十HDI ( M）一SR(M)A30A=A30A10. 75A44A=Al7一A 16A一A30AIF(A44A) 2185，2186，21862185 A44A=0.0A30A=Al7一A16A2186 CONTINUECC计算分量A16A,A30A和A44AGs/T 18340.4一2001TRASH= (A16一AA160.4910)(Al7/A16)IF(TRASH. LT.

47、 0. 0) TRASH=0. 0A17=A17一TRASHA16A二A16A一TRASHIF(Al6A.I.T.0.0) A16A=0.0IF(Al6A. EQ. 0. 0) A17=A30A一A44AIF(Al7. EQ. 0. 0) THENA16A=0. 0A30A=0. 0A44A=0. 0ELSEA16A二（A16A/Al7)，AA16A30A二（A30A/Al7)AA16A44A=(A4A/Al7)，AA16END IFCC计算A19的各个部分即A18A,A22S,A32ACA18A=0. 0DO 2187M=191，233，142187 A18A=Al8A+HDI(M)DO 2

48、188 M247,MS，142188 A18A=Al8A+SR(M)A22S二0.0DO 2189 M二247,289，142189 A22S二A22S+HDI(M）一SR(M)A22S二A22S/0. 75A32A=Al9一A18A一A22SIF(A32A) 2190,2191，21912190 A32A=0.0A22S二A19一A18A2191 CONTINUECC计算分量A18A,A22S和A32ACTRASH= (A18一AA18二0. 5073)(A19/Al8)IF(TRASH.I.T. 0.0) TRASH二0.0A19二A19一TRASHA 18A二A18A一TRASHIF (A18A. LT. 0. 0) A 18A二0.0IF(A18A. EQ. 0. 0) A19=A22S+A32AIF (A19. EQ. 0. 0) THENGB/T 18340.4-2001A18A=0. 0A22S=0. 0A32A=0. 0ELSEA18A二（A18A/A19)A22S= (A22S/A19）二A32A= (A32A/A19)END IF* AA18AA18二AA18