近红外光谱法测定毛竹综纤维素的含量研究

第２７卷第１期　２００７年２月　林　产　化　学　与　工　业　Ｃｈｅｍｉｓｔｒｙ　ａｎｄ　Ｉｎｄｕｓｔｙ　ｏｆ　Ｆｏｒｅｓｔｒ　Ｐｒｏｄｕｃｔｓ　ＶｏＩ＿２７　Ｎｏ．１　Ｆｅｂ．２００７　近红外光谱法测定毛竹综纤维素的含量研究　　，李改云ｒ，王戈　，黄安民　（１．中国林业科学研究院木材工业研究所，北京１０００９１；２．国际竹藤网络中心，北京１００１０２）　摘要：研究了用近红外光谱（ＮＩＲ）结合多变量统计分析技术对毛竹综纤维素含量的快速测定。　用常规实验室方－：－ｋ￣４定了５４个竹材样品的综纤维素含量，用近红外光谱仪采集相应样品的光谱，　ＪＩＡＮＧ　Ｚｅ．ｈｕｉ　对原始光谱进行二阶导数和２５点平滑预处理后，从５４个竹材样品中挑选４１个代表性的样品建模，　选择１０１ｌ一１６７５　ａｍ和１９３０—２４８８　ａｍ波段区间，用偏最小二乘法（ＰＩ＿Ｓ１）和完全交互验证方式建立　毛竹综纤维素含量的预测模型。结果表明，毛竹综纤维素含量和近红外光谱之间存在非常好的相关性，预测模型的相关　系数（Ｒ　）为０．９５，预测模型的标准偏差（ＳＥＰ）为０．７６％。　关键词：毛竹；综纤维素；近红外光谱；偏最小二乘法　中图分类号：ＴＱ９１　文献标识码：Ａ　文章编号：０２５３—２４１７（２００７）０１—００１５一ｏ４　Ｓｔｕｄｙ　ｏｎ　Ｄｅｔｅｒｍｉｎａｔｉｏｎ　ｏｆ　Ｈｏｌｏｃｅｌｌｕｌｏｓｅ　Ｃｏｎｔｅｎｔ　ｏｆ　Ｂａｍｂｏｏ　ｂｙ　Ｎｅａｒ—ｉｎｆｒａｒｅｄ　Ｒｅｆｌｅｃｔａｎｃｅ　Ｓｐｅｃｔｒｏｓｃｏｐｙ　ＪＩＡＮＧ　Ｚｅ．ｈｕｉ　，ＬＩ　Ｇａｉ．ｙｕｎ　，ＷＡＮＧ　Ｇｅ　，ＨＵＡＮＧ　Ａｎ．ｍｉｎ　（１．Ｒｅｓｅａｒｃｈ　Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｗｏｏｄ　Ｉｎｄｕｓｔｒｙ，ＣＡＦ，Ｂｅｉｊｉｎｇ　１０００９１，Ｃｈｉｎａ；　２．Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｃｅｎｔｅｒ　ｆｏｒ　Ｂａｍｂｏｏ　ａｎｄ　Ｒａｔｔｎ，Ｂｅｉａｊｉｎｇ　１００１０２，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ａ　ｒａｐｉｄ　ｄｅｔｅｒｍｉｎａｔｉｏｎ　ｏｆ　ｈｏｌｏｃｅｌｌｕｌｏｓｅ　ｃｏｎｔｅｎｔ　ｏｆ　ｂａｍｂｏｏ（Ｐｈｙｌｌｅｓｔａｃｈｙｓ　ｈｅｔｅｒｏｃｙｃｌａ（Ｃａｒｔ．）Ｍｉｆｏｒｄ　ＣＶ．ｐｕ，ｂｅｓｃｅｎｓ）ｈａｓ　ｂｅｅｎ　ｄｅｖｅｌｏｐｅｄ　ｂｙ　ｎｅａｒ—ｉｎｆｒｒｅｄ　ａｒｅｆｌｅｃｔａｎｃｅ　ｓｐｅｃｔｒｏｓｃｏｐｙ（ＮＩＲ）ｃｏｕｐｌｅｄ　ｂｙ　ｍｕｌｔｉｖａｒｉａｔｅ　ｓｔａｔｉｓｔｉｃａｌ　ｎａａｌｙｓｉｓ　ｔｅｃｈｎｉｑｕｅ．Ｔｈｅ　ｈｏｌｏ－　ｃｅｌｌｕｌｏｓｅ　ｃｏｎｔｅｎｔｓ　ｏｆ　５４　ｓａｍｐｌｅｓ　ｗｅｒｅ　ａｎａｌｙｚｅｄ　ｗｉｔｈ　ｔｒａｄｉｔｉｏｎａｌ　ｗｅｔ　ｃｈｅｍｉｃａｌ　ｔｅｃｈｎｉｑｕｅ．ｔｈｅ　ｓｐｅｃｔｒａ　ｏｆ　ｗｈｉｃｈ　ｗｅｒｅ　ｃｏｌｌｅｃｔｅｄ　ｂｙ　ＮＩＲ．Ａｆｔｅｒ　ｔｈｅ　ｓｐｅｃｔａ　ｗｅｒｅ　ｃｏｎｖｅ￣ｅｄ　ｔｒｏ　ｔｈｅ　ｓｅｃｏｎｄ　ｄｅｒｉｖａｔｉｖｅ　ｓｐｅｃｔｒａ　ｗｉｔｈ　ａ　ｇａｐ　ｏｆ　２５　ｎｉｎ，ｔｈｅ　ｐｒｅｄｉｃｔｉｏｎ　ｍｏｄｅｌ　ｆｏｒ　ｈｅ　ｈｏｌｏｃｅｌ—ｔ　ｌｕｌｏｓｅ　ｃｏｎｔｅｎｔ　ｏｆ　ｂａｍｂｏｏ　ｓａｍｐｌｅｓ　ＷａＳ　ｅｓｔａｂｌｉｓｈｅｄ　ｕｓｉｎｇ　ｐａｒｔｉｌａ　ｌｅａｓｔ－ｓｕａｑｒｅ　ｒｅｇｒｅｓｓｉｏｎ（ＰＩ＿Ｓ１）ａｎｄ　ｆｕｌｌ　ｃｒｏｓｓ　ｖａｌｉｄａｔｉｏｎ　ｉｎ　ｔｈｅ　ｒａｎｇｅ　ｏｆ　１０１１—１６７５　ｎｍ　ａｎｄ　１９３０—２４８８　ａｍ．Ｈｉ．ｓｈ　ｃｏｒｒｅｌａｔｉｏｎｓ　ｗｅｒｅ　ｏｂｔａｉｎｅｄ　ｂｅｔｗｅｅｎ　ｔｈｅ　ｐｒｅｄｉｃｔｅｄ　ＮＩＲ　ｒｅｓｕｌｔｓ　ａｎｄ　ｔｈｏｓｅ　ｏｂｔｉｎｅｄ　ｆａｒｏｍ　ｔｒａｄｉｔｉｏｎａｌ　ｃｈｅｍｉｃａｌ　ｍｅｔｈｏｄ．Ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｆｉｃｉｅｎｔ　ｏｆ　ｐｒｅｄｉｃｔｉｏｎ　ｍｏｄｅｌ　ｉｓ　０．９５　ａｎｄ　ｓｔａｎｄａｒｄ　ｅｒｒｏｒ　ｏｆ　ｐｒｅｄｉｃｔｉｏｎ　（ＳＥＰ）ｉｓ０．７６％．　Ｋｅｙ　ｗｏｒｄｓ：Ｐｈｙｌｌｅｓｔａｃｈｙｓ　ｈｅｔｅｒｏｃｙｃｌａ（Ｃａｒｔ．）Ｍｉｔｆｏｒｄ　ＣＶ．ｐｕｂｅｓｃｅｎｓ；ｈｏｌｏｅｅｌｌｌｏｓｅ；ｕｎｅａｒ－ｉｆｎｒａｒｅｄ　ｅｆｒｌｃｔｅａｎｃｅ　ｓｐｅｃｔｒｏｓｃｏｐｙ；ＰＬＳ１　近红外光谱区的吸收主要是分子或原子振动基频在２０００　ｃｍ　以上的倍频、合频吸收，记录的主要　是有机物中ｃ—Ｈ、Ｏ—Ｈ、Ｎ—Ｈ、Ｓ—Ｈ、Ｐ—Ｈ等含氢基团的倍频与合频吸收带。这些含氢基团的吸收　频率特征性强，受分子内外环境影响较小，光谱特性比较稳定。但倍频与合频吸收带吸收强度弱，比基　频吸收带宽得多，使得有机物近红外光谱谱带重叠严重，从而使近红外光谱的图谱解析异常困难，复杂　样品近红外光谱分析的信息提取主要依靠化学计量学方法与计算机数据处理…。综纤维素是纤维原　料除去抽出物和木质素后所留下的部分，主要包括纤维素和半纤维素，纤维素主要由葡萄糖基以／３－１，４　苷键连接起来的链状高分子化合物，而半纤维素是由不同糖基构成的复合聚糖的混合物　，综纤维素　中大量的ｃ—Ｈ和Ｏ—Ｈ含氢基团使其在近红外区域具有丰富的吸收。目前测量植物纤维原料中综纤　维素含量的方法主要是依据国家标准，首先用苯一醇溶液抽提样品，然后用亚氯酸钠将木质素氯化和氧　收稿日期：２００６—０６—２６　基金项目：国家林业局“９４８”引进项目“人工林木材的ＮＩＲ材性预测及增值利用”（２００３—０４－２７）　作者简介：（１９３８一），女，江苏扬州人，教授，博士生导师，主要从事木材和竹材材性、加工与利用研究　通讯作者：李改云。助理研究员。博士生，主要从事木材化学与胶黏剂研究；Ｅ・ｎｌａｉｌ．．１ｉｇｙ＠ｆｏｒｅｓｔｒｙ．ａｃ．ｃａ。　维普资讯 http://www.cqvip.com

ｌ６　林产化学与工业　第２７卷　化后除去，测量一个样品需２—３　ｄ，操作手续繁琐，成本高，无法实现大批量样品的快速测定。近红外光　谱法是一种快速、无损的检测方法，国外已有用近红外技术测定人工林蓝桉（Ｅｕｃａｌｙｐｔｕｓ　６　糖、甘露糖和半乳糖含量及３种落叶松（　Ｌ．）和　亮果桉（Ｅ．ｎｉｔｅｎｓ）的纤维素含量和木质素含量、火炬松（Ｐｉｎｕｓ　ｔａｅｄａ　Ｌ．）的木质素、抽提物、葡萄糖、木　ｓｐ．）心材抽提物的研究报道　』。中国林科院近年来已开　始从事近红外技术在木材科学上的应用研究工作，并取得一定的研究成果　。但国内外目前尚无用近　红外光谱法预测竹材中综纤维素含量的相关报道。本研究依据我国制浆造纸工业国标测定竹材中的综　纤维素含量，利用近红外光谱技术和化学计量学软件建立了近红外光谱数据和竹子中综纤维素含量的　相关模型，对实现综纤维素含量的快速、准确测定和推进近红外光谱技术在我国木材科学与造纸领域中　的应用具有重要意义。　１　实验部分　１．１试剂与仪器　原料：毛竹（Ｐｈｙｌｌｅｓｔａｃｈｙｓ　ｈｅｔｅｒｏｃｙｃｌａ（Ｃａｒｒ．）Ｍｉｔｆｏｒｄ　ＣＶ．ｐｕｂｅｓｃｅｎｓ）试材采自浙江省杭州市萧山区。　取２、４和６年生毛竹各３棵，取０．５—２．５　ｍ和４．５—６．５　ｍ高的两个竹段，风干，每个竹段分成竹青、竹　肉和竹黄３部分，切成小薄片后粉碎、过筛，取通过０．４５　ｍｍ筛孔而截留在０．３０　ｍｍ筛孔的竹粉，放置　于样品瓶中，供分析和采谱用。试剂：亚氯酸钠（７４．６％），工业纯；苯、乙醇、冰醋酸、丙酮、硫酸均为分　析纯。仪器：ＬａｂＳｐｅｃ＠Ｐｒｏ近红外光谱仪，检测器：３５０—２５００　ｎｍ，低噪声５１２阵元ＰＤＡ；光谱采样间　隔：１．４　ｎｍ（３５０—１０５０　ｎｍ），２　ｎｍ（１０００—２５００　ｎｍ）；光谱分辨率：３　ｎｍ（７００　ｎｍ）和１０　ｎｍ（１４００　ｎｍ，　２ｌ００　ｎｍ），美国ＡＳＤ公司生产。　１．２综纤维素含量的测定　按照国标ＧＢ／Ｔ　２６７７．１０—１９９５测定竹粉中综纤维素的含量。同时做两份平行测定，取其算术平均　值作为测定结果。　１．３　ＮＩＲ光谱采集　近红外光谱仪和光谱采集均在装有空调的恒温（２０　ｏＣ±２　ｏＣ）室内进行，样品在恒温室内敞口放置　２周后采谱。利用美国ＡＳＤ公司开发的Ｉｎｄｉｃｏ数据采集和预处理软件采集竹粉的近红外光谱。先用杯　光源监测器对商用聚四氟乙烯白板进行空白校准后，再对样品的近红外光谱在全光谱范围内（３５０—　２５００　ｎｍ）进行采集，每次扫描３０次并自动平均为一个光谱，每个样品装样３次。　１．４模型建立与评价　完成对竹子粉末样品的光谱数据进行微分和平滑预处理后，使用瑞典ＣＡＭＯ公司的多变量统计分　析软件Ｕｎｓｃｒａｍｂｌｅｒ对每个样品的３个光谱平均后用偏最小二乘法和完全交互验证方式在样品综纤维　素含量和光谱数据之间建立模型，选择不同谱区、主成分维数及选择建模样品等以确立最佳建模条件。　用内部交叉验证方式评价模型的优劣，相关系数越高，标准偏差越小，模型质量越好。　２结果与讨论　２．１毛竹综纤维素含量的变异性　样品的选择对模型的建立具有非常大的影响，为了使选择的样品具有代表性，考虑到毛竹生长过程　和构造特点，按照实验部分所列方案制备了样本，表ｌ对５４个竹子样本分别按照竹龄、横向和纵向进行　了统计。如表中所述，毛竹作为一种天然生长的有机材料，随着竹龄的增长，竹杆材质生长经历增进期、　稳定期和下降期３个阶段，构成细胞壁主要成分之一的综纤维素含量也发生变化，４年生毛竹的综纤维　素含量达到较大值，６年生毛竹的综纤维素含量略有下降。竹壁的内、中、外３层因构造的不同，综纤维　素含量差别较大，从内到外，综纤维素含量逐渐增加。综纤维素在竹杆纵向部位分布比较均匀，为了使　样品的综纤维素含量不同，本试验选取０．５—２．５　ｍ和４．５—６．５　ｍ两个竹段。结果表明，随着竹杆部位　从０．５—２．５　ｍ升高到４．５—６．５　ｍ，综纤维素含量增加。　维普资讯 http://www.cqvip.com

第１期　，等：近红外光谱法测定毛竹综纤维素的含量研究　１７　建立近红外预测模型的化学值来自参比分析方法，预测模型的精确度与准确度依赖于参比分析的　精确度与准确度　。本试验５４个毛竹样品的综纤维素含量由１个经验丰富的分析人员依照ＧＢ／Ｔ　２６７７．１０—１９９５测定，每个样品作平行试验，５４个样品平行试验的平均标准偏差为０．１７，说明综纤维素　含量的实验室数据具有高度重复性。　表１毛竹综纤维素分析结果　Ｔａｂｌｅ　１　ＡｎＭ￣ｉｃＭ　ｒｅｓｕｌｔｓ　ｏｆ　ｈｏｌｏｃｅｌｌｕｌｏｓｅ　ｃｏｎｔｅｎｔｓ　ｆｏｒ　ｔｈｅ　５４　ｓａｍｐｌｅｓ　辍蹄鑫¨　Ｏ　２　Ｏ　２．２光谱预处理　图１是５４个竹粉样品的原始近红外光谱谱图。因为近红外光谱仪采集的原始数据中除包含与样　本物质化学结构相关的信息外，还包含其他干扰信息，影响建模的准确性。在建模之前，对近红外光谱　的原始数据进行处理可以获得高信噪比、低背景干扰的分析信号。对光谱进行微分处理可以消除基线　漂移或平缓背景干扰的影响，提供比原始光谱更高的分辨率和更清晰的光谱轮廓变化。光谱数据的平　滑处理可有效平滑高频噪声，提高信噪比…。作者对原始近红外光谱数据进行二阶微分处理的同时，　采用２５点平滑处理，扣除仪器背景或漂移对信号的影响，较好的净化了谱图信息，同时去掉高频噪音对　信号的干扰，最终提高了光谱的分辨率。　Ｉ　Ｌ。　一Ｉｌ　．酞　ｊｌ　入　ｒＩ　●　．　Ｖ’　．　竹　Ｉ．　２０００　２５Ｏ０　０　５００　１０００　１５Ｏ０　波长／ｎｍ　波长／ｎｍ　图１　５４个竹子粉末样品的ＮＩＲ原始谱图　Ｆｉｇ．１　Ｒａｗ　ＮＩＲ　ｓｐｅｃｔｒａ　ｆｏｒ　５４　ｐｏｗｄｅｒｅｄ　ｂａｍｂｏｏ　ｓａｍｐｌｅｓ　图２　５４个竹子粉末样品的二阶导数近红外光谱　Ｆｉｇ．２　Ｔｈｅ　ｓｅｃｏｎｄ　ｄｅｒｉｖａｔｉｖｅ　ＮＩＲ　ｓｐｅｃｔｒａ　ｄａｔａ　ｏｆ　５４　ｏｗｄｅｒｐｅｄ　ｂａｍｂｏｏ　ｓａｍｐｌｅｓ　２．３波长和样品的选择　竹材在３５０—１０１０　ａｍ的可见光和短波区域内，光谱信息量较少（图１）。在预处理后的谱图和相关　系数图中（图２和图３），对相关性存在严重的干扰，建立模型时应消除此波段的影响。在１６７６—　１８８０　ａｍ波长区域，样品的吸光度变化较小，信噪比差，分辨率较低，而在相关系数图中的贡献也非常　小。１９２３　ａｍ附近为水的合频吸收峰，对模型的精度有负面影响，删除１６７６—１９２９　ａｍ波段可屏蔽水分　子对模型的影响，提高模型的准确度。在近红外的长波区域（２４８９—２５００　ａｍ），主要是一些干扰信息，　经过二阶微分处理后干扰信息被消除。消除３５０—１０１０　ａｍ，１６７６—１９２９　ａｍ和２４８９—２５００　ａｍ这三段光　维普资讯 http://www.cqvip.com

１８　林产化学与工业　第２７卷　谱后，即选择１０１１～１６７５　ｎｍ和１９３０～２４８８　ｎｍ波段区间，提高了模型的精度，获得了很好的预测效果。　代表性样品的选择是近红外分析中很重要的一个环节。足量样品是建立模型的前提，但近红外数学　模型包含的样品数也不是越多越好，样品数越多引入　误差的机会越多…。本试验首先依据毛竹生长过程　２０ｏ０　４０００　赧　０　ｍ　ｈ　。　和构造特点，选取了５４个毛竹样品，并对其综纤维素　含量进行了分析，然后按照综纤维素的化学值挑选出　４１个样品用于建模。利用偏最ｄｘ－－乘法和全交互验　证方式建立预测模型，选择主成分数为６，模型的相关　系数（　）达到０．９５，预测模型的标准偏差（ＳＥＰ）为　０．７６％　回　一２０ｏ０　３结论　对竹材的漫反射近红外原始光谱进行二阶微分　和２５点平滑处理后，从５４个竹材样品中挑选４１个代　表陛的样品，在合适的光谱区间内，用偏最小二乘法　建立了毛竹综纤维素含量的预测模型，取得了很好的预测效果，证明用近红外光谱技术可以实现对竹材　样品中综纤维素含量的快速测定。　参考文献：　［１］严衍禄．近红外光谱分析基础与应用［Ｍ］．北京：中国轻工业出版社，２００５：９，９８—１０３，１２２．　［２］杨淑蕙．植物纤维化学［Ｍ］．北京：中国轻工业出版社，２００１：１６３，２１４．　［３］ＫＥＬＬＥＹ　Ｓ　Ｓ，ＲＩＡＬＳ　Ｔ　Ｇ，ＳＮＥＬＬ　Ｒ，ｅｔ　ａ１．Ｕｓｅ　ｏｆ　ｎｅａｒ　ｉｎｆｒａｒｅｄ　ｓｐｅｃｔｒｏｓｃｏｐｙ　ｔｏ　ｍｅａｓｕｒｅ　ｔｈｅ　ｃｈｅｍｉｃａｌ　ａｎｄ　ｍｅｃｈａｎｉｃａｌ　ｐｒｏｐｏｍｅｓ　ｏｆ　ｓｏｈｄ　ｗｏｏｄ［Ｊ］．Ｗｏｏｄ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ，２００４，３８（４）：２５７．　［４］ＳＣＨＩＭＬＥＣＫ　Ｌ　Ｒ，ＷＲＩＧＨＴ　Ｐ　Ｊ，ＭＩＣＨＥＬＬ　Ａ　Ｊ，ｅｔ　ａ１．Ｎｅａｒ－ｉｎｆｒａｒｅｄ　ｓ￣ｃｔｒａ　ａｎｄ　ｃｈｅｍｉｃａｌ　ｃｏｍｐｏｓｉｔｉｏｎｓ　ｏｆＥ　ｇｌｏｂｕｌｅｓ　ａｎｄ　Ｅ　ｎ／ｔｅｎｓ　ｐｌａｎｔａ－　ｔｉｏｎ　ｗｏｏｄｓ［Ｊ］．Ａｐｐｉｔａ　Ｊｏｕｒｎａｌ，１９９７，５０（１）：４０－４６．　［５］ＹＡＭＡＤＡ　Ｔ，ＹＥＨ　Ｔ　Ｆ，ＣＨＡＮＧ　Ｈ　Ｍ，ｅｔ　ａ１．Ｒａｐｉｄ　ａｎａｌｙｓｉｓ　ｏｆ　ｔｒａｎｓｇｅｎｉｃ　ｔｒｅｅｓ　ｕｓｉｎｇ　ｔｒａｎｓｍｉｔｔａｎｃｅ　ｎｅａｒ—ｉｎｆｒａｒｅｄ　ｓｐｅｃｔｒｏｓｃｏｐｙ（ＮＩＲ）［Ｊ］．　Ｈｏｌｚｆｏｒｓｅｈｕｎｇ，２００６，６０（１）：２４—２８．　［６］ＧＩＥＲＨＮＧＥＲ　Ｎ，ＳＣＨＷＡＮＮＩＮＧＥＲ　Ｍ，ＨＩＮＴＥＲＳＴＯＩＳＳＥＲ　Ｂ，ｅｔ　ａ１．Ｒａｐｉｄ　ｄｅｔｅｒｍｉｎａｔｉｏｎ　ｏｆ　ｈｅａｒｔｗｏｏｄ　ｅｘｔｒａｅｔｉｖｅｓ　ｉｎ　Ｆｏｕｒｉｅｒ　ｔｒａｎｓｆｏｒｍ　ｎｅａｒ　ｉｆｒｎａｒｅｄ　ｓｐｅｃｔｒｏｓｃｏｐｙ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｎｅａｒ　Ｉｎｆｒａｒｅｄ　Ｓｐｅｃｔｒｏｓｃｏｐｙ，２００２，１０（３）：２０３—２１４．　ｓｐ．ｂｙ　ｍｅａｎｓ　ｏｆ　［７］，黄安民，费本华，等．利用近红外光谱和ｘ射线衍射技术分析木材微纤丝角［Ｊ］．光谱学与光谱分析，２００６，２６（７）：１２３０一　】２３３．　关于《林产化学与工业》改双月刊的紧急通知　《林产化学与工业》自１９８１年创刊以来，得到了广大读者和作者的关心与大力支持，在此我们表示　衷心的感谢。由于大家的厚爱，本刊的稿源非常充实，而刊物容量有限，因而刊出周期较长。为了缩短　稿件的发表周期，更进一步提高刊物质量，更快、更好、更全面地为广大读者和作者服务，经国家林业局、　江苏省新闻出版局批准，《林产化学与工业》定于２００７年起由季刊改为双月刊，逢双月月末出版，定价　１５．００元，全年９０．００元。欢迎广大读者及时到邮局订阅，同时欢迎广大作者、读者继续踊跃投稿，我们　将竭诚为您服务！