环境科学  2017, Vol. 38 Issue (5): 1769-1774   PDF    
沧州市春季NMHCs空间分布特征
段菁春1 , 周雪明2,3 , 张鹤丰1 , 谭吉华2,3 , 胡京南1 , 柴发合1     
1. 中国环境科学研究院, 环境基准与风险评估国家重点实验室, 北京 100012;
2. 中国科学院大学资源与环境学院, 北京 100049;
3. 中国科学院生态环境研究中心, 北京 100085
摘要: 2015年春季在沧州市城区、郊区和潜在污染源附近选择了15个采样点进行同期采样.研究表明沧州市NMHCs总体上市区高于近郊及远郊区县;市区以高新区NMHCs浓度最高;郊县采样点除河间市略高外,其它采样点的浓度均明显低于市区浓度;机动车的道路排放是沧州市NMHCs的重要来源之一;沧州大化和沧州炼油在停产期间未对市区NMHCs产生明显影响;大港油田采油三厂采取了较完善的油气回收措施,未对市区NMHCs产生明显影响;平均来看,沧州市NMHCs中烷烃占65%,烯烃占16%,芳烃占19%;臭氧生成潜势(ozone formation potential,OFPs)主要来源于二甲苯(19%)、乙烯(14%)、甲苯(11%)、丙烯(5%)、异戊烷(5%)和异戊烯(5%)等;气溶胶生成潜势(formation potential of secondary organic aerosol,SOAFPs)主要来源于甲苯(28%)、蒎烯(28%)、二甲苯(16%)、乙苯(9%)和苯(9%)等.
关键词: 非甲烷烃 (NMHCs)      沧州      空间变化      臭氧生成潜势      挥发性有机物     
Spatial Distribution Characteristics of NMHCs in Spring in Cangzhou City
DUAN Jing-chun1 , ZHOU Xue-ming2,3 , ZHANG He-feng1 , TAN Ji-hua2,3 , HU Jing-nan1 , CHAI Fa-he1     
1. State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China;
2. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China;
3. Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
Abstract: Simultaneous collections of non-methane hydrocarbons (NMHCs) were carried out at 15 sampling sites including urban, suburb and potential pollution areas in Cangzhou City in spring 2015. The results showed that NMHCs were generally higher in urban areas than those in suburb and rural areas; the highest concentration of NMHCs was observed at Cangzhou High-tech zone (urban area); the concentrations of NMHCs were significantly lower at rural sites than in most urban sites except Hejian site; vehicular emissions were the main sources of NMHCs in Cangzhou; Cangzhou chemical fertilizer plant and Cangzhou oil refinery had no significant influence on urban NMHCs during their shutdown period; Dagang Oilfield, with better oil and gas recovery systems, did not have a significant impact on urban NMHCs. In general, alkanes, alkenes and aromatics accounted for 65%, 16% and 19% of NMHCs in Cangzhou City, respectively; xylene (19%), ethylene (14%), toluene (11%), propylene (5%), isopentane (5%) and isopentene (5%) were the most dominant contributors to ozone formation potential; aerosol formation potential was mainly derived from toluene (28%), pinene (28%), xylene (16%), ethylbenzene (9%) and benzene (9%).
Key words: non-methane hydrocarbons (NMHCs)      Changzhou      spatial variation      ozone formation potential      volatile organic compounds (VOCs)     

非甲烷烃 (non-methane hydrocarbons, NMHCs) 是环境空气中最主要的挥发性有机物 (volatile organic compounds, VOCs)[1].在阳光紫外线照射下, NMHCs不仅可以通过氢氧 (·OH) 等自由基发生大气光化学反应, 氧化后再与氮氧化物 (NOx) 反应生成臭氧, 某些NMHCs如甲苯、乙苯、二甲苯和三甲苯等还会通过光化学反应生成高浓度的二次有机气溶胶 (secondary organic aerosol.SOA)[2~4], 同时NMHCs中的部分化合物如苯和甲苯还会对人体产生直接的毒害作用[5]. NMHCs来源复杂, 既有生物来源, 如植物排放[6]; 也有人为来源[7~10], 如机动车排放、油品挥发、工业燃烧排放和涂料生产及应用排放等.

沧州市位于中国空气污染最严重京津冀地区的东南部, 冬季霾频发, 以PM2.5为代表大气污染物浓度居高不下, 夏季近地面臭氧污染严重, 这均与沧州市NMHCs污染具有一定联系[1].根据文献[11], 沧州市2014年PM2.5平均浓度为88 μg·m-3, 超过现行环境空气质量标准1.52倍.沧州市2015年机动车保有量已超过150万辆.虽然沧州市在推广新能源汽车、提升燃油品质、开展加油站油气回收利用等方面采取了相当多的措施, 机动车VOCs排放仍不容忽视.除了机动车排放外, 沧州市还存在大量石油开采、炼油和化工等排放VOCs的工业企业.近年来, 国内研究者虽然对北京、广州、上海等大城市大气VOCs, 特别是NMHCs的污染特征, 时空变化规律开展了较多的研究[12~15], 但对沧州市这样的地级市NMHCs的时空分布特征报道较少.本研究通过对沧州市15个采样点的同期采样, 开展沧州市NMHCs空间分布和臭氧生成潜势的初步研究, 该研究对于了解沧州市NMHCs污染特征, 以期为环境保护部门制定VOCs控制对策具有重要参考意义.

1 材料与方法 1.1 采样点选取

为获得沧州市NMHCs空间分布特征, 同时初步了解主要NMHCs污染源的排放特征以及对环境大气的影响, 本研究共选取15个采样点.其中市环境保护局、沧县城建局、电视转播站、省质控站等4个采样点位于主城区内, 代表城市污染特征; 经济开发区和高新区2个采样点, 代表工业开发区污染特征; 在东南西北这4个方向选取黄骅市、泊头市、河间市和青县4个采样点, 代表沧州市远郊区县污染特征; 利用采样车以巡回采样方式在沧州市潜在污染源沧州大化和沧州炼油的下风向, 大港油田采油三厂作业区内, 城市重要交通十字路口, 以及沧州近郊农田内设置大化门口、沧炼北、采油三厂和新华路口这4个潜在污染源采样点, 代表潜在污染源可能对沧州市污染的影响; 选取沧州西站西1个采样点, 位于沧州市近郊农业区, 远离特点污染源, 代表沧州市区背景污染水平, 具体采样点信息详见图 1表 1.

图 1 沧州市NMHCs采样点位置分布示意 Fig. 1 Map of sampling sites of NMHCs in Cangzhou City

表 1 沧州市NMHCs罐采样采样点位及采样时间 Table 1 Information of location and data during the NMHCs sampling in Cangzhou City

1.2 样品采集

采样日期为2015年4月27日, 在沧州市10个环境空气质量监测站点同期进行采样, 采样时间为上午09:00;另有5个非固定采样点采取巡回采样的方式进行采样, 采样时间为上午08:13~11:15.非甲烷烃 (NMHCs) 用容积为2L的美国NUPRO公司97-300系列采样罐采集, 采集前用美国ENTECH公司3100型打样罐清洗仪用高纯氮气 (>99.999%) 清洗5次后抽成真空.采样时缓慢打开采样阀门, 尽量使空气均匀进入罐中, 当采样罐中的气压与空气大气压平衡时关闭阀门, 每个样品采样时间大约3min.样品采集后送中国科学院广州地球化学所有机地球化学国家重点实验室进行定量分析, 分析工作在一周内完成.

1.3 样品分析

NMHCs样品的分析在预浓缩系统/色谱/质谱上根据美国EPA的TO-14法测定完成.分析条件:色谱柱为HP-VOC毛细管柱 (60 mL×0.32 mm i. d. ×1.8 μm), 柱温40℃保留2 min, 然后以6℃·min-1由40℃升至230℃, 再保留10 min.氦气流速为1.2 mL·min-1, 分流模式为不分流进样, 进样量为250 mL. NMHCs通过对照标样与样品的保留时间以及质谱图进行定性并通过色谱图峰面积计算定量.混合标样 (Supelco TO-14 Calibration Mix) 稀释成7个不同浓度做标准曲线, 标准曲线的回归系数R2大于0.995.

2 结果与讨论 2.1 沧州市15个采样点NMHCs浓度分布及组成特征

图 2可以看出, 沧州市NMHCs以新华路口机动车道路边最高, 表明机动车的排放是沧州市NMHCs的重要来源之一; 采样点大化门口和沧炼北距离沧州市区重要NMHCs污染源沧州大化和沧州炼油较近, 虽然采样点NMHCs略有升高, 但与城区其它采样点相比变化不大, 这主要是因为采样期间两企业均停产; 采样点采油三厂是潜在的油气挥发污染源之一, 但实际NMHCs浓度较低, 这表明大港油田采油三厂采取了较完善的油气回收措施; 郊县采样点青县、河间市、泊头市和黄骅市, 除河间市略高外, 其它郊县采样点的浓度均明显低于沧州市区的采样点浓度, 表明沧州市的NMHCs可能主要来自于市区; 在沧州市区的几个采样点中高新区NMHCs浓度最高, 表明这里可能存在涉VOCs的污染源排放.

图 2 沧州市NMHCs浓度分布特征 Fig. 2 Spatial characteristics of NMHCs in Cangzhou City

从总体上看 (图 2), 沧州市NMHCs中烷烃、烯烃、芳烃和炔烃的构成比例较为接近, 其中烷烃约占53.1%, 烯烃12.6%, 芳烃15.6%和炔烃18.7%.

图 3可以看出, 沧州市NMHCs中, 烷烃以乙烷、丙烷、异丁烷、正丁烷、异戊烷和正戊烷等为主; 烯烃以乙烯、丙烯、异丁烯和异戊烯等为主; 芳烃以苯、甲苯、二甲苯和乙苯等为主.

图 3 沧州市NMHCs物种构成特征 Fig. 3 Species composition characteristics of NMHCs in Cangzhou City

2.2 沧州市NMHCs的臭氧生成潜势

为了评估NMHCs臭氧生成潜势 (ozone formation potential, OFPs), 最常用的方法是利用NMHCs组分的最大臭氧益增活性 (MIR) 计算臭氧生成潜势的方法[16~18].臭氧生成潜势可以用VOCs物种的浓度与MIR系数的乘积来衡量, 计算公式为:

(1)

式中, [VOC]i为观测到大气中某种VOC浓度; MIRi为该VOC化合物在臭氧最大增量反应中的臭氧生成系数.

图 4表 2可以看出沧州市OFPs以新华路口机动车道路边最高, 表明机动车的排放对沧州市臭氧的形成有重要影响; 工业点源大化门口和沧炼北与城区其它采样点OFPs相比变化不大, 这主要是因为采样期间两企业均停产; 潜在污染点源采油三厂采取了较完善的油气回收措施, OFPs也较低; 郊县采样点青县、河间市、泊头市和黄骅市, 除河间市略高外, 其它郊县采样点的OFPs浓度均明显低于沧州市区; 沧州市区除电视转播站和火车西站二个点较低外, 其他几个点OFPs都比较高, 由于市区受机动车影响较大, 说明市区机动车对臭氧的生成有重要影响.从沧州市臭氧生成潜势分布来看, OFPs主要来源于芳烃贡献, 其次是烯烃和烷烃, 炔烃对OFPs也有少量贡献.从图 5可以看出沧州市的臭氧生成潜势 (OFPs) 最主要的物种来源由大到小依次为:二甲苯 (19%)、乙烯 (14%)、甲苯 (11%)、丙烯 (5%)、异戊烷 (5%) 和异戊烯 (5%) 等.

表 2 沧州市芳烃、烷烃、烯烃、炔烃、臭氧生成潜势以及气溶胶生成潜势/μg·m-3 Table 2 Concentrations of aromatics, alkanes, alkenes, alkyne, OFPs and SOAFPs in Cangzhou City/μg·m-3

图 4 沧州市臭氧生成潜势 (OFPs) 的空间分布 Fig. 4 Spatial characteristics of OFPs in Cangzhou City

图 5 沧州市OFPs物种构成特征 Fig. 5 Species composition characteristics of OFPs in Cangzhou City

2.3 沧州市NMHCs的气溶胶生成潜势

大气中的VOCs可以通过·OH、NO3和O3等氧化剂的氧化作用生成SOA. SOA不仅会恶化空气质量危害人体健康, 而且会降低大气能见度并影响全球气候变化.国内外研究者提出多种方法评估SOA对大气颗粒的贡献率[19~24], 但由于不同的方法各有其不足之处, 因此现有的方法尚不能准确评估SOA的浓度水平及贡献率.国内研究表明, 利用元素碳示踪法评估, 广州SOA对PM2.5的贡献约为4.2%~6.8%[25].本研究通过气溶胶生成系数 (fractional aerosol coefficient, FAC) 对SOA的生成潜势 (formation potential of secondary organic aerosol, SOAFPs) 进行估算[26~30], 气溶胶生成系数FAC是SOA生成潜势的一种表达方式.对于可以形成SOA的VOC组分i, 其FAC可以定义为:

(2)

这样, 前体VOC的环境初始浓度乘以相应的FAC值就可以得到生成的SOA浓度.目前使用比较多的一套FAC值是Grosjean等[26, 27]在综合大量烟雾箱实验数据和大气化学动力学数据的基础上提出[28~30].这种方法的优点在于它能从VOCs的排放清单或环境浓度直接估算出SOA的浓度, 并能反映各SOA前体物的相对贡献[4].

图 6可以看出沧州市市区气溶胶生成潜势除电视转播站和火车西站二个点较低外, 其他几个点SOAFPs都比较高.其中以新华路口机动车道路边最高, 表明机动车的排放对沧州市二次有机气溶胶的形成有重要影响; 郊县采样点青县、河间市、泊头市和黄骅市, 除河间市略高外, 其它郊县采样点的SOAFPs浓度均明显低于沧州市区.结果还表明沧州市气溶胶生成潜势主要来源于芳烃和烯烃, 烷烃和炔烃贡献微小.从图 7可以看出沧州市的气溶胶生成潜势 (SOAFPs) 最主要的物种来源由大到小依次为:甲苯 (28%)、蒎烯 (28%)、二甲苯 (16%)、乙苯 (9%) 和苯 (9%) 等.

图 6 沧州市气溶胶生成潜势 (SOAFPs) 的空间分布 Fig. 6 Spatial characteristics of SOAFPs in Cangzhou City

图 7 沧州市SOAFPs物种构成特征 Fig. 7 Species composition characteristics of SOAFPs in Cangzhou City

3 结论

(1) 沧州市区以高新区NMHCs浓度最高, 表明这里可能存在涉VOCs的污染源排放; 远郊区县除河间市略高外, 其它采样点的浓度均明显低于沧州市区浓度; 机动车的排放是沧州市NMHCs的重要来源之一; 停产期间未发现沧州大化和沧州炼油对市区NMHCs的明显影响; 大港油田采油三厂未对市区NMHCs产生明显影响.

(2) 沧州市NMHCs中, 烷烃以乙烷、丙烷、异丁烷、正丁烷、异戊烷和正戊烷等为主; 烯烃以乙烯、丙烯、异丁烯和异戊烯等为主; 芳烃以苯、甲苯、二甲苯和乙苯等为主.

(3) 沧州市的臭氧生成潜势 (OFPs) 从物种上主要源于二甲苯 (1%)、乙烯 (14%)、甲苯 (11%)、丙烯 (5%)、异戊烷 (5%) 和异戊烯 (5%) 等.

(4) 沧州市的气溶胶生成潜势 (SOAFPs) 从物种来源上看主要源于甲苯 (28%)、蒎烯 (28%)、二甲苯 (16%)、乙苯 (9%) 和苯 (9%) 等.

(5) 建议沧州市加强对机动车排放水平的管控, 同时加强高新区涉VOCs企业的排放管理.

致谢: 感谢中国科学院广州地球化学研究所王新明老师为采样提供的采样罐和对NMHCs的分析定量; 感谢沧州市环境保护局和环境监测中心为本研究提供的采样帮助.
参考文献
[1] Duan J C, Tan J H, Yang L, et al. Concentration, sources and ozone formation potential of volatile organic compounds (VOCs) during ozone episode in Beijing[J]. Atmospheric Research, 2008, 88(1): 25–35. DOI: 10.1016/j.atmosres.2007.09.004
[2] Tan J H, Duan J C, Chen D H, et al. Chemical characteristics of haze during summer and winter in Guangzhou[J]. Atmospheric Research, 2009, 94(2): 238–245. DOI: 10.1016/j.atmosres.2009.05.016
[3] 谭吉华, 赵金平, 段菁春, 等. 广州秋季灰霾污染过程大气颗粒物有机酸的污染特征[J]. 环境科学, 2013, 34(5): 1982–1987. Tan J H, Zhao J P, Duan J C, et al. Pollution characteristics of organic acids in atmospheric particles during haze periods in autumn in Guangzhou[J]. Environmental Science, 2013, 34(5): 1982–1987.
[4] 吕子峰, 郝吉明, 段菁春, 等. 北京市夏季二次有机气溶胶生成潜势的估算[J]. 环境科学, 2009, 30(4): 969–975. Lv Z F, Hao J M, Duan J C, et al. Estimate of the formation potential of secondary organic aerosol in Beijing Summertime[J]. Environmental Science, 2009, 30(4): 969–975.
[5] Tang J H, Chan L Y, Chan C Y, et al. Characteristics and diurnal variations of NMHCs at urban, suburban, and rural sites in the Pearl River Delta and a remote site in South China[J]. Atmospheric Environment, 2007, 41(38): 8620–8632. DOI: 10.1016/j.atmosenv.2007.07.029
[6] 张钢锋, 谢绍东. 基于树种蓄积量的中国森林VOC排放估算[J]. 环境科学, 2009, 30(10): 2816–2822. Zhang G F, Xie S D. Estimation of VOC Emission from forests in China based on the volume of tree species[J]. Environmental Science, 2009, 30(10): 2816–2822. DOI: 10.3321/j.issn:0250-3301.2009.10.002
[7] 邵敏, 赵美萍, 白郁华, 等. 燕山石化地区NMHC的特征研究[J]. 环境化学, 1994, 13(1): 40–45. Shao M, Zhao M P, Bai Y H, et al. A study of NMHC character in atmosphere of yanshan petrochemical region[J]. Environmental Chemistry, 1994, 13(1): 40–45.
[8] 沈旻嘉, 郝吉明, 王丽涛. 中国加油站VOC排放污染现状及控制[J]. 环境科学, 2006, 27(8): 1473–1478. Shen M J, Hao J M, Wang L T. VOC emission situation and control measures of gas station in China[J]. Environmental Science, 2006, 27(8): 1473–1478.
[9] 魏巍, 王书肖, 郝吉明. 中国涂料应用过程挥发性有机物的排放计算及未来发展趋势预测[J]. 环境科学, 2009, 30(10): 2809–2815. Wei W, Wang S X, Hao J M. Estimation and forecast of volatile organic compounds emitted from paint uses in China[J]. Environmental Science, 2009, 30(10): 2809–2815. DOI: 10.3321/j.issn:0250-3301.2009.10.001
[10] 马永亮, 谭吉华, 贺克斌, 等. 佛山灰霾期挥发性有机物的污染特征[J]. 环境科学, 2011, 32(12): 3549–3554. Ma Y L, Tan J H, He K B, et al. Characteristics of volatile organic compounds during haze episode in Foshan City[J]. Environmental Science, 2011, 32(12): 3549–3554.
[11] 沧州市环境保护局. 2014年沧州市环境状况公报[R]. 沧州: 沧州市环境保护局, 2015.
[12] 张新民, 柴发合, 岳婷婷, 等. 天津武清大气挥发性有机物光化学污染特征及来源[J]. 环境科学研究, 2012, 25(10): 1085–1091. Zhang X M, Chai F H, Yue T T, et al. Photochemical characteristics and sources of volatile organic compounds in Wuqing, Tianjin[J]. Research of Environmental Sciences, 2012, 25(10): 1085–1091.
[13] 吴方堃, 王跃思, 安俊琳, 等. 北京奥运时段VOCs浓度变化, 臭氧产生潜势及来源分析研究[J]. 环境科学, 2010, 31(1): 10–16. Wu F K, Wang Y S, An J L, et al. Study on concentration, ozone production potential and sources of VOCs in the atmosphere of Beijing during olympics period[J]. Environmental Science, 2010, 31(1): 10–16.
[14] Pang X B, Mu Y J, Zhang Y J, et al. Contribution of isoprene to formaldehyde and ozone formation based on its oxidation products measurement in Beijing, China[J]. Atmospheric Environment, 2009, 43(13): 2142–2147. DOI: 10.1016/j.atmosenv.2009.01.022
[15] Xu J, Ma J Z, Zhang X L, et al. Measurements of ozone and its precursors in Beijing during summertime:impact of urban plumes on ozone pollution in downwind rural areas[J]. Atmospheric Chemistry and Physics, 2011, 11(23): 12241–12252. DOI: 10.5194/acp-11-12241-2011
[16] 周雪明, 项萍, 段菁春, 等. 佛山市冬夏季非甲烷烃污染特征研究[J]. 环境科学, 2016, 37(11): 4124–4132. Zhou X M, Xiang P, Duan J C, et al. Pollution characteristics of non-methane hydrocarbons during winter and summer in Foshan City[J]. Environmental Science, 2016, 37(11): 4124–4132.
[17] Carter W P L. Development of ozone reactivity scales for volatile organic compounds[J]. Air & Waste, 1994, 44(7): 881–899.
[18] So K L, Wang T. C3-C12 non-methane hydrocarbons in subtropical Hong Kong:spatial-temporal variations, source-receptor relationships and photochemical reactivity[J]. Science of the Total Environment, 2004, 328(1-3): 161–174. DOI: 10.1016/j.scitotenv.2004.01.029
[19] Schauer J J, Rogge W F, Hildemann L M, et al. Source apportionment of airborne particulate matter using organic compounds as tracers[J]. Atmospheric Environment, 1996, 30(22): 3837–3855. DOI: 10.1016/1352-2310(96)00085-4
[20] Zheng M, Cass G R, Schauer J J, et al. Source apportionment of PM2.5 in the southeastern United States using solvent-extractable organic compounds as tracers[J]. Environmental Science & Technology, 2002, 36(11): 2361–2371.
[21] Pandis S N, Harley R H, Cass G R, et al. Secondary organic aerosol formation and transport[J]. Atmospheric Environment. Part A. General Topics, 1992, 26(13): 2269–2282. DOI: 10.1016/0960-1686(92)90358-R
[22] Strader R, Lurmann F, Pandis S N. Evaluation of secondary organic aerosol formation in winter[J]. Atmospheric Environment, 1999, 33(29): 4849–4863. DOI: 10.1016/S1352-2310(99)00310-6
[23] Hildemann L M, Rogge W F, Cass G R, et al. Contribution of primary aerosol emissions from vegetation-derived sources to fine particle concentrations in Los Angeles[J]. Journal of Geophysical Research, 1996, 101(D14): 19541–19549. DOI: 10.1029/95JD02136
[24] Turpin B J, Huntzicker J J. Identification of secondary organic aerosol episodes and quantitation of primary and secondary organic aerosol concentrations during SCAQS[J]. Atmospheric Environment, 1995, 29(23): 3527–3544. DOI: 10.1016/1352-2310(94)00276-Q
[25] Duan J C, Tan J H, Cheng D X, et al. Sources and characteristics of carbonaceous aerosol in two largest cities in Pearl River Delta Region, China[J]. Atmospheric Environment, 2007, 41(14): 2895–2903. DOI: 10.1016/j.atmosenv.2006.12.017
[26] Grosjean D. In situ organic aerosol formation during a smog episode:estimated production and chemical functionality[J]. Atmospheric Environment. Part A. General Topics, 1992, 26(6): 953–963. DOI: 10.1016/0960-1686(92)90027-I
[27] Grosjean D, Seinfeld J H. Parameterization of the formation potential of secondary organic aerosols[J]. Atmospheric Environment (1967), 1989, 23(8): 1733–1747. DOI: 10.1016/0004-6981(89)90058-9
[28] Barthelmie R J, Pryor S C. Secondary organic aerosols:formation potential and ambient data[J]. Science of the Total Environment, 1997, 205(2-3): 167–178. DOI: 10.1016/S0048-9697(97)00200-3
[29] Dechapanya W, Russell M, Allen D T. Estimates of anthropogenic secondary organic aerosol formation in Houston, Texas[J]. Aerosol Science and Technology, 2004, 38(S1): 156–166.
[30] Kourtidis K, Ziomas I. Estimation of secondary organic aerosol (SOA) production from traffic emissions in the city of Athens[J]. Global NEST Journal, 1999, 1(1): 33–38.