亚星游戏官网-www.yaxin868.com

专利名称：基于中红外光谱的植物¹⁵N示踪同位素丰度快速检测方法
技术领域：
本发明属于信息检测技术领域，涉及一种基于中红外光谱的植物1N示踪同位素丰度快速检测方法。
背景技术：
氮是植物必需的营养元素，是作物施肥三要素之首，氮肥占我国化肥消费比重达 60%左右。目前我国占世界7%的耕地上消耗着全球35%的氮肥，氮肥的过量施用给我国的农业生产和生态环境带来了巨大�：�。中国每年因不合理施肥造成近1000万吨的氮素流失，直接经济损失约300亿元。而且由于氮肥利用率低下，有60%左右进入环境。因而过量施氮已成为威胁中国长期粮食安全和环境安全的重大问题。减少氮肥施用量的关键是提高氮肥利用率，而15N同位素示踪(标记)法是目前唯一能直接反映作物对氮肥实际利用情况的方法。1N同位素示踪技术的关键是要对1N同位素丰度的动态变化进行跟踪和检测。但是目前15N同位素分析测量仪器昂贵、成本高、操作繁琐、耗时长、需要大量化学试剂，难以满足农业生产经济、高效的要求。目前，15N同位素丰度测量的方法有质谱法、1M发射光谱法、核磁共振法和中子活化分析等，其中应用较多的是质谱法。质谱法是通过测量各种物质的原子或原子团的荷质比来确定被测原子或原子团的性质和质量，具有高灵敏度及精确度的特征，广泛用于同位素的实验室测量。但是这种方法所采用的仪器价格昂贵、系统庞大、操作复杂，样品制备繁琐，操作过程引入污染。且最终的结果只能给出元素的量，不能测定元素的化学形态，这限制了质谱法在农业生产中的应用。1N发射光谱法需要将样品转化成氮气( )，然后用适当方法来激发氮分子发光， 得到氮的光谱。随着组成氮分子的氮原子质量数不同而使氮的光谱波长发生位移，利用这一现象来测定1N同位素丰度。与质谱法相比,1 发射光谱法不需要昂贵的仪器设备，操作相对简便，并能测定含氮量为微克级的样品。但是，有些元素的原子很难激发，而且有机物一般是大分子，组成复杂，激发后谱线杂乱难以分析，限制了 1N发射光谱法的应用。核磁共振法是利用I-IO2兆赫的电磁波照射置于强磁场中的样品，样品中的某些原子核在特定的磁场强度下，与照射电磁波发生共振现象，产生强弱不同的吸收信号。以此来判明某些核素在分子中所处的位置及某些功能团上核素的相对数目。但是由于其造价非常昂贵，而且测定灵敏度较低，故极少用于15N同位素的肥料利用率研究。中子活化分析是通过鉴别和测试样本因辐照感生的放射性核素的特征辐射，进行元素和核素分析的放射分析化学方法。该法的优点是灵敏度极高，准确度和精密度也很高， 可测定元素范围广；但同样存在仪器价格昂贵，分析周期较长，操作技术复杂的缺陷，在1N 同位素丰度的测量上鲜见报道。而且一般情况下，中子活化分析只能给出元素的量，不能测定元素的化学形态。中红外吸收光谱法就是利用物质对中红外光区电磁辐射的选择性吸收来捕捉分子中有关基团结构定性、定量信息的分析方法。中红外吸收光谱法以其高度的特征性，固、 液、气态样品均可用，且用量少、不破坏样品、分析速度快，已成为分析化学中结构定性、定量分析最成熟的手段之一。中红外吸收光谱特别适合于捕捉光谱同位素效应，因为绝大多数的有机化合物和很多无机化合物的化学键转动-振动的基频都出现在中红外区。由于同位素的核质量不同而使分子的转动和振动初始能态发生变化，由之引起分子光谱的谱线位移，所以每一种同位素分子都有其专属的转动-振动光谱。物质对中红外辐射的吸收强度与物质含量的关系符合朗伯一比尔定律，所以中红外光谱可用于定量分析。但是中红外光子能量较低，传统光谱仪的检测器灵敏度也很低，限制了中红外光谱在定量分析中的广泛应用。随着计算机技术的深入发展，傅里叶变换红外光谱-(FTIR)技术逐渐成为红外光谱分析的重要工具。与传统色散型光谱仪相比，FIlR不需要分光，信号强，灵敏度很高；由于采用干涉仪，与扫描单色仪相比扫描速度提高了数百倍，而且干涉仪中没有狭缝的限制，其辐射通量比色散型仪器大得多。基于此类优点，FTIR 特别适用于测量弱信号光谱，而且具有很高的灵敏度。凭借其高灵敏度、高分辨率、高信噪比、扫描迅速、光通量大、联机操作性强及高度计算机化的特点，FIlR推动经典中红外光谱技术的定性、定量分析能力取得了突破性进展。由于朗伯-比尔定律具有可加性，对于有几个组分的混合物来说，如果每个组分都符合朗伯-比尔定律，那么混合物的中红外光谱是各个纯成分中红外光谱的叠加。由于中红外光谱的谱带较多，选择余地较大，有利于排除干扰找到对应每个组分的且不受其它组分吸收峰干扰的“独立峰”，基于“独立峰”的筛选和指认，定量分析可以通过直接测定混合物的中红外光谱来实现，从而简化了繁琐的制样过程。目前基于中红外光谱技术的稳定同位素分析还处于初探阶段，大部分研究局限于光谱同位素效应的定性解析，少量的定量分析研究对象仅以简单气体化合物为主。而对于 1N示踪植物体内15N同位素丰度的中红外光谱测量国内外还未见报道。然而，植物体内1N 同位素丰度快速检测方法的研究，对于研究植物的营养、植物的生长生理、以及植物与环境的关系都具有极其重要的推动作用。尤其对于研究植物对氮素的吸收、运转、同化和代谢等过程以及与生态环境相互作用机理上将起到巨大的作用。相对传统1N同位素丰度检测法的破坏性实验而言，中红外光谱测量具有无损、快速、多组分同时检测的优点，可以准确而全面地对作物整个生长周期的氮吸收、利用的动态情况进行跟踪，并且可以减少因植物体间个体差异明显，采样点、采样频率和采样时机等随机性和偶然性因素所带来的影响。同时，中红外光谱技术能有效克服传统1~同位素分析法测量费用昂贵、分析测试耗时、采样数目偏少和检测周期长等缺陷。

发明内容
本发明提出了植物1N示踪同位素丰度的中红外光谱快速检测方法，为快速有效地监测植物体内1N同位素丰度的动态变化规律提供理论依据，有望实现植物1N示踪同位素丰度的快速、低成本测量。发明提供一种快速测量植物1N示踪同位素丰度的方法，它包括
①首先把待测植物样本烘干，然后在红外灯照射下把样本放入玛瑙研钵研磨成粒径小于20 μ m的粉末，按照1:20 49 (质量比)的比例与干燥的溴化钾粉末混勻，使用压片机把混勻后的粉末制成约2mm的薄片，完成样本制备。②把制备好的样本固定在支架上，放入中红外光谱仪的样品室，进行波数4000 cnT1 400 CnT1范围，分辨率为1 CnT1的中红外透射光谱扫描，获取样本的中红外透射光
■i並曰ο③把样本的中红外透射光谱转换成吸光度，基于特征波数1235 cnT1 1071 cnT1
和3475 cm—1 3465 cm—1处的吸光度，采用11点平均法提取得到16个特征变量，把这16 个特征变量代入下式，即可计算得到植物样本的1N示踪同位素丰度。
] 其中Γ代表模型对样本的1EN示踪同位素丰度的预测值；A代表样本在1235cm 1广-1225cm 1的11个波· 处吸光度的平均值；5代表样本在12 cm 1214cm的11个波_女处吸光度的平均值；巧代表样本在1213cm 广一 1203cm的11个波娄 处吸光度的平均值；X4代表样本在1202cm 1192cm的11个波_女处吸光度的平均值；巧代表样本在1191cm 广一 1181cm的11个波娄 处吸光度的平均值；々代表样本在1180cm 广一 1170cm的11个波娄 处吸光度的平均值；X7代表样本在1169cm 1159cm的11个波_女处吸光度的平均值；而代表样本在1158cm 广一 1148cm的11个波娄 处吸光度的平均值；巧代表样本在1147cm 广一 1137cm的11个波娄 处吸光度的平均值；%代表样本在1136cm 1126cm1的11个波I发处吸光度的平均值%代表样本在1125cm 1115cm1的11个波I发处吸光度的平均值%代表样本在1114cm 1104cm1的11个波I发处吸光度的平均值%代表样本在1103cm 1093cm1的11个波I发处吸光度的平均值%代表样本在1092cm 1082cm1的11个波I发处吸光度的平均值%代表样本在1081cm 1071cm1的11个波I发处吸光度的平均值%代表样本在；3475cm 3465cm1的11个波I发处吸光度的平均值。
本发明所测量的同位素丰度是相对丰度，即N15原子数目占所有氮原子(包括N15原子和N14原子)数目的百分含量。
本发明可测定植物包括水稻、油菜、番茄、小麦，茶叶、粗叶木、山矾、光叶山矾、华
南省藤、山月桂、高山蒲葵、裂叶棕梅、锈毛新木姜子、九节、海南罗伞、平滑琼楠、中华厚壳桂、柏拉木、露兜草、蜜花山矾、斜茎盘算子、海南木莲、雄鸡树、海南杨桐、披针山矾、降真香、鸡毛松、残枝蒲桃、毛宁、白肉榕、大萼木姜子、燕尾葵、毛荔枝、油丹、桃榄中的1N示踪同位素丰度。本发明与背景技术相比，具有以下特点
(1)快速，中红外光谱技术具有特征性强、灵敏度高的优点，避免了传统15N同位素丰度测量的分离、提取等繁琐、耗时的样品制备过程。而其多通道优点，扫描速度快，可在Is内完成整个中红外波段范围(4000 cnT1 400 cnT1)的扫描，更大大加快了测量时间。(2)低成本，中红外光谱仪与传统分析法所用的质谱仪、核磁共振仪等仪器设备相比，成本大大降低。(3)具有良好的经济效益，传统的测量手段在取样、制样、测定等方面需要耗费大量的人力、财力、物力，本测量方法简单、使用方便，可以快速、准确的测量植物1N示踪同位素丰度，故具有良好的经济效益。


图1是本发明方法流程图2是1N示踪茶树样本的中红外吸收光谱图。
具体实施例方式以下结合附图对本发明作进一步说明。如图1所示，本发明方法具体包括以下步骤
一、建立中红外光谱与植物1N示踪同位素丰度间的定量关系模型 ①.1N同位素示踪下茶树的栽培。以茶树品种龙井43为研究对象，采用盆栽法，每盆栽4株1年生无性系茶苗，栽后修剪，进行正常的水分管理和病虫防治及防冻措施。用1M 标记尿素为氮源，采用丰度分别为8%，5%，4%和洲尿素形成4个水平的氮示踪梯度和一个空白试验对照(即施用非1N标记的普通尿素)，每个梯度设计处理4盆样品，每盆搭配施入 0. 6g磷和0. 6g钾，肥料于春季分两次施下，间隔1月。②.收集1N示踪茶树样本。施肥后30天、60天、90天、120天分别取样，每次取各盆中一株，4次重复。取样时，按盆连同根系取出整株植株，并拣净落在土壤中的根系， 然后分株按根、茎、叶等器官取样，共得到67个样本。于80°C烘箱内烘至恒重，在红外灯照射下用玛瑙研钵研磨成粒径小于20 μ m的粉末，然后按1:49的比例把粉末与溴化钾磨细混勻，用钢铲把磨细混勻的粉末约0. 05g移入压片模具，把压片模具水平放置在压片机座上， 施加20t/cm2的压力，得到厚度约2mm的样本薄片。③.获取样本的中红外吸收光谱。把制备好的样本薄片固定在支架上，放入中红外光谱仪的样本室，中红外光谱仪的光谱采集范围是4000 cnT1 400 cm"1,光谱分辨率为 1 cnT1，光谱采集时室温为25°C左右，湿度小于65%，每个样本的光谱是64次扫描的平均值， 采集样本的透射光谱。把样本的透射光谱，按照公式“吸光度=2_log (透过率％) ”的数学变换得到样本的中红外吸收光谱如图2所示。随后用质谱法测量样本的15N同位素丰度。
④.基于粒子群优化法捕捉植物N15同位素丰度的中红外光谱特征波数。基于粒子群优化法对示踪植物的N15同位素丰度的特征波数进行优化和筛�。�由于4000 cm—1 400 cm-1范围共3601个波数变量，故设定对应的优化矩阵为1行3601列的二进制代码矩阵，优化矩阵中向量取值为0或1，0代表所对应的波数为非特征波数，1代表所对应的波数为特征波数。优化矩阵中所有为1的向量对应的波数成为特征波数组合，以这些特征波数的吸光度作为自变量，N15同位素丰度作为因变量，建立多元线性回归模型。粒子群优化的目标函数是多元线性回归模型的均方根误差最小。当均方根误差最小时对应粒子群优化法寻找得到的最优特征波数。研究发现，对于植物N15同位素丰度测量的特征波数为1235 cm—1 1071 CnT1 和 3475 cnT1 3465 cnT1。⑤.建立样本中红外光谱与15N同位素丰度之间的定量关系模型。基于样本的中红外吸收光谱，获得样本在特征波数处(1235 cnT1 1071 cnT1和3475 cnT1 3465 cnT1) 的吸光度，采用11点平均法提取得到16个特征变量，依据朗伯-比尔定律，应用多元线性回归算法建立这16个特征变量与1N同位素丰度的定量关系模型。模型的数学表达式如下所示
Y = 3.6<76eUB + 3.5433^8 + !,OSOe^T2 - 4.25 "^ - IMle^7Xi - 2.913e^% + 3.1060-07 — 5.966e^8x7 — 3.03 . — 3.130^"07 +3.647e^x10 — 3.066^^7 X11 + 2.09 ^8 +9.0520^8 + 5.094^-08% — 9.086^^38 +5.14 -°8
其中1代表模型对样本的15N同位素丰度的预测值；
巧代表样本在1235cm 1 1225cm 1的11个波· 处吸光度的平均值h代表样本在12 cm 1214cm的11个波_女处吸光度的平均值巧代表样本在1213cm 1203cm的11个波娄 处吸光度的平均值X4代表样本在1202cm 1192cm的11个波_女处吸光度的平均值A代表样本在1191cm 1181cm的11个波娄 处吸光度的平均值A代表样本在1180cm 1170cm的11个波娄 处吸光度的平均值巧代表样本在1169cm 1159cm的11个波_女处吸光度的平均值巧代表样本在1158cm 1148cm的11个波娄 处吸光度的平均值、代表样本在1147cm 1137cm的11个波娄 处吸光度的平均值%代表样本在1136cm 1126cm1的11个波I发处吸光度的平均值&代表样本在1125cm 1115cm1的11个波I发处吸光度的平均值%代表样本在1114cm 1104cm1的11个波I发处吸光度的平均值%代表样本在1103cm 1093cm1的11个波I发处吸光度的平均值￥代表样本在1092cm 1082cm1的11个波I发处吸光度的平均值b代表样本在1081cm 1071cm1的11个波I发处吸光度的平均值~代表样本在3475 cm—1 3465 cm—1的11个波数处吸光度的平均值。二、基于植物1N示踪同位素丰度的中红外光谱测量模型对未知样本的1N示踪同位素丰度进行测定
基于所建立的植物1~示踪同位素丰度的中红外光谱测量模型对67个样本的预测相关系数为0.719，均方根误差为2.75 ,，标准均方根误差为7. 488e_°5%。具体预测结果如表 ι所示。说明该测量模型能够实现植物1N示踪同位素丰度的快速检测。表1、模型对于67个样本的预测结果
权利要求
1.基于中红外光谱的植物1N示踪同位素丰度快速检测方法，其特征在于该方法包括以下步骤步骤1.把待测植物样本烘干，然后在红外灯照射下把样本放入玛瑙研钵研磨成粒径小于20 μ m的粉末，按照质量比1 20 49的比例与干燥的溴化钾粉末混勻，使用压片机把混勻后的粉末制成约2mm的薄片，完成样本制备；步骤2.把制备好的样本固定在支架上，放入中红外光谱仪的样品室，进行波数4000 cnT1 400 CnT1范围，分辨率为1 CnT1的中红外透射光谱扫描，获取样本的中红外透射光谱；步骤3.把样本的中红外透射光谱转换成吸光度，基于特征波数1235 cnT1 1071 cm—1和；3475 cm—1 ；3465 cm—1处的吸光度，采用11点平均法提取得到16个特征变量，把这16个特征变量代入下式，即可计算得到植物样本的15N示踪同位素丰度 F，-β8 1 m -O r., r\ t r, -OS+3. , +1.080/ — x-s - 1.643e X4 -+ 3.106 -07 — 5.966e^8j — 3.08 ^8 — 3.130^7 +3.647g^7x + 2.09Ie^z12 +9.052i ^8%+ 5.09 ^8 - 9.086e^8%+5.14 ^8, 其中.〗代表模型对样本的1M示踪同位素丰度的预测值；的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值； 1的11个波数处吸光度的平均值。X511巧代表样本在1235cm 1 “-1225cm&代表样本在12 cm 1214cm而代表样本在1213cm 广一 1203cmX4代表样本在1202cm 1192cm 代表样本在1191cm 广一 1181cm4代表样本在1180cm 广一 1170cm巧代表样本在1169cm 1159cm巧代表样本在1158cm 广一 1148cmA代表样本在1147cm 广一 1137cm%代表样本在1136cm 1126cmh代表样本在1125cm 1115cm%代表样本在1114cm 1104cm￥代表样本在1103cm 1093cm%代表样本在1092cm 1082cm%代表样本在1081cm 1071cm 代表样本在3475cm 3465cm
全文摘要
本发明涉及一种基于中红外光谱的植物15N示踪同位素丰度快速检测方法。传统的检测方法速度慢且操作复杂。本发明首先把待测植物样本研磨成粉末样本，将粉末样本与溴化钾后压片，完成样本制备。然后将制备好的样本进行的中红外透射光谱扫描，得到样本的中红外透射光谱，再把样本中红外透射光谱转换成吸光度。最后将样本在特征波数处的吸光度，代入公式，即可得到样本的同位素丰度。本发明使用方便，可以快速、有效的测量植物15N示踪同位素丰度，故具有良好的经济效益。
文档编号G01N21/35GK102313713SQ201110196570
公开日2012年1月11日 申请日期2011年7月14日 优先权日2011年7月14日
发明者何勇, 李晓丽, 聂鹏程 申请人:浙江大学