CN102313713B - 基于中红外光谱的植物15n示踪同位素丰度快速检测方法 - Google Patents
基于中红外光谱的植物15n示踪同位素丰度快速检测方法 Download PDFInfo
- Publication number
- CN102313713B CN102313713B CN 201110196570 CN201110196570A CN102313713B CN 102313713 B CN102313713 B CN 102313713B CN 201110196570 CN201110196570 CN 201110196570 CN 201110196570 A CN201110196570 A CN 201110196570A CN 102313713 B CN102313713 B CN 102313713B
- Authority
- CN
- China
- Prior art keywords
- sample
- wave number
- mean value
- representative sample
- number place
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000700 radioactive tracer Substances 0.000 title claims abstract description 21
- 238000001228 spectrum Methods 0.000 title abstract description 16
- 238000001514 detection method Methods 0.000 title abstract description 5
- 238000002835 absorbance Methods 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 31
- 239000000843 powder Substances 0.000 claims abstract description 11
- IOLCXVTUBQKXJR-UHFFFAOYSA-M potassium bromide Chemical compound [K+].[Br-] IOLCXVTUBQKXJR-UHFFFAOYSA-M 0.000 claims abstract description 8
- 238000002360 preparation method Methods 0.000 claims abstract description 6
- 230000003595 spectral effect Effects 0.000 claims description 13
- 239000002245 particle Substances 0.000 claims description 7
- 238000004433 infrared transmission spectrum Methods 0.000 claims description 6
- 238000002156 mixing Methods 0.000 claims description 6
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 239000004570 mortar (masonry) Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 9
- 241000196324 Embryophyta Species 0.000 description 29
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 23
- 230000000155 isotopic effect Effects 0.000 description 22
- 229910052757 nitrogen Inorganic materials 0.000 description 11
- 238000005259 measurement Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 238000002329 infrared spectrum Methods 0.000 description 9
- 239000000618 nitrogen fertilizer Substances 0.000 description 7
- 230000035945 sensitivity Effects 0.000 description 7
- 238000000862 absorption spectrum Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 6
- 238000004949 mass spectrometry Methods 0.000 description 6
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000004993 emission spectroscopy Methods 0.000 description 4
- 239000003337 fertilizer Substances 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 235000009024 Ceanothus sanguineus Nutrition 0.000 description 3
- 240000003553 Leptospermum scoparium Species 0.000 description 3
- 235000015459 Lycium barbarum Nutrition 0.000 description 3
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 description 3
- 238000012271 agricultural production Methods 0.000 description 3
- 229940037003 alum Drugs 0.000 description 3
- 125000004429 atom Chemical group 0.000 description 3
- 238000003947 neutron activation analysis Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012113 quantitative test Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 244000269722 Thea sinensis Species 0.000 description 2
- 239000004202 carbamide Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000005445 isotope effect Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 125000004433 nitrogen atom Chemical group N* 0.000 description 2
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 241001096696 Alseodaphnopsis hainanensis Species 0.000 description 1
- 244000144730 Amygdalus persica Species 0.000 description 1
- 241000345998 Calamus manan Species 0.000 description 1
- 241001164374 Calyx Species 0.000 description 1
- 244000025254 Cannabis sativa Species 0.000 description 1
- 241001081178 Cryptocarya Species 0.000 description 1
- 240000007687 Dacrycarpus imbricatus Species 0.000 description 1
- 240000003537 Ficus benghalensis Species 0.000 description 1
- 241000238631 Hexapoda Species 0.000 description 1
- 244000048199 Hibiscus mutabilis Species 0.000 description 1
- 235000003973 Hibiscus mutabilis Nutrition 0.000 description 1
- 238000004566 IR spectroscopy Methods 0.000 description 1
- 241000193986 Kalmia latifolia Species 0.000 description 1
- 244000007392 Lambertia formosa Species 0.000 description 1
- 235000004461 Lambertia formosa Nutrition 0.000 description 1
- 241000633855 Litsea pungens Species 0.000 description 1
- 240000007163 Livistona chinensis Species 0.000 description 1
- 235000007688 Lycopersicon esculentum Nutrition 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 241000712582 Neolitsea aurata Species 0.000 description 1
- 240000007817 Olea europaea Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 1
- 235000006040 Prunus persica var persica Nutrition 0.000 description 1
- 240000003768 Solanum lycopersicum Species 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 244000228451 Stevia rebaudiana Species 0.000 description 1
- 235000006092 Stevia rebaudiana Nutrition 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 244000098338 Triticum aestivum Species 0.000 description 1
- 229920001807 Urea-formaldehyde Polymers 0.000 description 1
- 241000607479 Yersinia pestis Species 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000012136 culture method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 235000020774 essential nutrients Nutrition 0.000 description 1
- 230000004720 fertilization Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 229910052500 inorganic mineral Chemical class 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N iron(III) oxide Inorganic materials O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 150000002605 large molecules Chemical group 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 239000011707 mineral Chemical class 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 235000015277 pork Nutrition 0.000 description 1
- 239000011591 potassium Substances 0.000 description 1
- 229910052700 potassium Inorganic materials 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 235000012950 rattan cane Nutrition 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000035943 smell Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000001845 vibrational spectrum Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Description
技术领域
本发明属于信息检测技术领域,涉及一种基于中红外光谱的植物15N示踪同位素丰度快速检测方法。
背景技术
氮是植物必需的营养元素,是作物施肥三要素之首,氮肥占我国化肥消费比重达60%左右。目前我国占世界7%的耕地上消耗着全球35%的氮肥,氮肥的过量施用给我国的农业生产和生态环境带来了巨大危害。中国每年因不合理施肥造成近1000万吨的氮素流失,直接经济损失约300亿元。而且由于氮肥利用率低下,有60%左右进入环境。因而过量施氮已成为威胁中国长期粮食安全和环境安全的重大问题。减少氮肥施用量的关键是提高氮肥利用率,而15N同位素示踪(标记)法是目前唯一能直接反映作物对氮肥实际利用情况的方法。15N同位素示踪技术的关键是要对15N同位素丰度的动态变化进行跟踪和检测。但是目前15N同位素分析测量仪器昂贵、成本高、操作繁琐、耗时长、需要大量化学试剂,难以满足农业生产经济、高效的要求。
目前,15N同位素丰度测量的方法有质谱法、15N发射光谱法、核磁共振法和中子活化分析等,其中应用较多的是质谱法。
质谱法是通过测量各种物质的原子或原子团的荷质比来确定被测原子或原子团的性质和质量,具有高灵敏度及精确度的特征,广泛用于同位素的实验室测量。但是这种方法所采用的仪器价格昂贵、系统庞大、操作复杂,样品制备繁琐,操作过程引入污染。且最终的结果只能给出元素的量,不能测定元素的化学形态,这限制了质谱法在农业生产中的应用。
15N发射光谱法需要将样品转化成氮气(N2),然后用适当方法来激发氮分子发光,得到氮的光谱。随着组成氮分子的氮原子质量数不同而使氮的光谱波长发生位移,利用这一现象来测定15N同位素丰度。与质谱法相比,15N发射光谱法不需要昂贵的仪器设备,操作相对简便,并能测定含氮量为微克级的样品。但是,有些元素的原子很难激发,而且有机物一般是大分子,组成复杂,激发后谱线杂乱难以分析,限制了15N发射光谱法的应用。
核磁共振法是利用1-102兆赫的电磁波照射置于强磁场中的样品,样品中的某些原子核在特定的磁场强度下,与照射电磁波发生共振现象,产生强弱不同的吸收信号。以此来判明某些核素在分子中所处的位置及某些功能团上核素的相对数目。但是由于其造价非常昂贵,而且测定灵敏度较低,故极少用于15N同位素的肥料利用率研究。
中子活化分析是通过鉴别和测试样本因辐照感生的放射性核素的特征辐射,进行元素和核素分析的放射分析化学方法。该法的优点是灵敏度极高,准确度和精密度也很高,可测定元素范围广;但同样存在仪器价格昂贵,分析周期较长,操作技术复杂的缺陷,在15N同位素丰度的测量上鲜见报道。而且一般情况下,中子活化分析只能给出元素的量,不能测定元素的化学形态。
中红外吸收光谱法就是利用物质对中红外光区电磁辐射的选择性吸收来捕捉分子中有关基团结构定性、定量信息的分析方法。中红外吸收光谱法以其高度的特征性,固、液、气态样品均可用,且用量少、不破坏样品、分析速度快,已成为分析化学中结构定性、定量分析最成熟的手段之一。中红外吸收光谱特别适合于捕捉光谱同位素效应,因为绝大多数的有机化合物和很多无机化合物的化学键转动-振动的基频都出现在中红外区。由于同位素的核质量不同而使分子的转动和振动初始能态发生变化,由之引起分子光谱的谱线位移,所以每一种同位素分子都有其专属的转动-振动光谱。
物质对中红外辐射的吸收强度与物质含量的关系符合朗伯—比尔定律,所以中红外光谱可用于定量分析。但是中红外光子能量较低,传统光谱仪的检测器灵敏度也很低,限制了中红外光谱在定量分析中的广泛应用。随着计算机技术的深入发展,傅里叶变换红外光谱-(FTIR)技术逐渐成为红外光谱分析的重要工具。与传统色散型光谱仪相比,FTIR不需要分光,信号强,灵敏度很高;由于采用干涉仪,与扫描单色仪相比扫描速度提高了数百倍,而且干涉仪中没有狭缝的限制,其辐射通量比色散型仪器大得多。基于此类优点,FTIR特别适用于测量弱信号光谱,而且具有很高的灵敏度。凭借其高灵敏度、高分辨率、高信噪比、扫描迅速、光通量大、联机操作性强及高度计算机化的特点,FTIR推动经典中红外光谱技术的定性、定量分析能力取得了突破性进展。
由于朗伯-比尔定律具有可加性,对于有几个组分的混合物来说,如果每个组分都符合朗伯-比尔定律,那么混合物的中红外光谱是各个纯成分中红外光谱的叠加。由于中红外光谱的谱带较多,选择余地较大,有利于排除干扰找到对应每个组分的且不受其它组分吸收峰干扰的“独立峰”,基于“独立峰”的筛选和指认,定量分析可以通过直接测定混合物的中红外光谱来实现,从而简化了繁琐的制样过程。
目前基于中红外光谱技术的稳定同位素分析还处于初探阶段,大部分研究局限于光谱同位素效应的定性解析,少量的定量分析研究对象仅以简单气体化合物为主。而对于15N示踪植物体内15N同位素丰度的中红外光谱测量国内外还未见报道。然而,植物体内15N同位素丰度快速检测方法的研究,对于研究植物的营养、植物的生长生理、以及植物与环境的关系都具有极其重要的推动作用。尤其对于研究植物对氮素的吸收、运转、同化和代谢等过程以及与生态环境相互作用机理上将起到巨大的作用。相对传统15N同位素丰度检测法的破坏性实验而言,中红外光谱测量具有无损、快速、多组分同时检测的优点,可以准确而全面地对作物整个生长周期的氮吸收、利用的动态情况进行跟踪,并且可以减少因植物体间个体差异明显,采样点、采样频率和采样时机等随机性和偶然性因素所带来的影响。同时,中红外光谱技术能有效克服传统15N同位素分析法测量费用昂贵、分析测试耗时、采样数目偏少和检测周期长等缺陷。
发明内容
本发明提出了植物15N示踪同位素丰度的中红外光谱快速检测方法,为快速有效地监测植物体内15N同位素丰度的动态变化规律提供理论依据,有望实现植物15N示踪同位素丰度的快速、低成本测量。
发明提供一种快速测量植物15N示踪同位素丰度的方法,它包括:
①首先把待测植物样本烘干,然后在红外灯照射下把样本放入玛瑙研钵研磨成粒径小于20μm的粉末,按照1:20~49(质量比)的比例与干燥的溴化钾粉末混匀,使用压片机把混匀后的粉末制成约2mm的薄片,完成样本制备。
②把制备好的样本固定在支架上,放入中红外光谱仪的样品室,进行波数4000cm-1~400cm-1范围,分辨率为1cm-1的中红外透射光谱扫描,获取样本的中红外透射光谱。
③把样本的中红外透射光谱转换成吸光度,基于特征波数1235cm-1~1071cm-1和3475cm-1~3465cm-1处的吸光度,采用11点平均法提取得到16个特征变量,把这16个特征变量代入下式,即可计算得到植物样本的15N示踪同位素丰度。
Y=3.676e-03+3.543e-08x1+1.080e-07x2-4.259e-08x3-1.643e-07x4-2.913e-08x5
+3.106e-07x6-5.966e-08x7-3.084e-08x8-3.130e-07x9+3.647e-07x10-3.066e-07x11
+2.091e-08x12+9.052e-08x13+5.094e-08x14-9.086e-08x15+5.143e-08x16
其中:Y代表模型对样本的15N示踪同位素丰度的预测值;
x1代表样本在1235cm-1~1225cm-1的11个波数处吸光度的平均值;
x2代表样本在1224cm-1~1214cm-1的11个波数处吸光度的平均值;
x3代表样本在1213cm-1~1203cm-1的11个波数处吸光度的平均值;
x4代表样本在1202cm-1~1192cm-1的11个波数处吸光度的平均值;
x5代表样本在1191cm-1~1181cm-1的11个波数处吸光度的平均值;
x6代表样本在1180cm-1~1170cm-1的11个波数处吸光度的平均值;
x7代表样本在1169cm-1~1159cm-1的11个波数处吸光度的平均值;
x8代表样本在1158cm-1~1148cm-1的11个波数处吸光度的平均值;
x9代表样本在1147cm-1~1137cm-1的11个波数处吸光度的平均值;
x10代表样本在1136cm-1~1126cm-1的11个波数处吸光度的平均值;
x11代表样本在1125cm-1~1115cm-1的11个波数处吸光度的平均值;
x12代表样本在1114cm-1~1104cm-1的11个波数处吸光度的平均值;
x13代表样本在1103cm-1~1093cm-1的11个波数处吸光度的平均值;
x14代表样本在1092cm-1~1082cm-1的11个波数处吸光度的平均值;
x15代表样本在1081cm-1~1071cm-1的11个波数处吸光度的平均值;
x16代表样本在3475cm-1~3465cm-1的11个波数处吸光度的平均值。
本发明所测量的同位素丰度是相对丰度,即N15原子数目占所有氮原子(包括N15原子和N14原子)数目的百分含量。
本发明可测定植物包括水稻、油菜、番茄、小麦,茶叶、粗叶木、山矾、光叶山矾、华南省藤、山月桂、高山蒲葵、裂叶棕梅、锈毛新木姜子、九节、海南罗伞、平滑琼楠、中华厚壳桂、柏拉木、露兜草、蜜花山矾、斜茎盘算子、海南木莲、雄鸡树、海南杨桐、披针山矾、降真香、鸡毛松、残枝蒲桃、毛宁、白肉榕、大萼木姜子、燕尾葵、毛荔枝、油丹、桃榄中的15N示踪同位素丰度。
本发明与背景技术相比,具有以下特点:
(1)快速,中红外光谱技术具有特征性强、灵敏度高的优点,避免了传统15N同位素丰度测量的分离、提取等繁琐、耗时的样品制备过程。而其多通道优点,扫描速度快,可在1s内完成整个中红外波段范围(4000cm-1~400cm-1)的扫描,更大大加快了测量时间。
(2)低成本,中红外光谱仪与传统分析法所用的质谱仪、核磁共振仪等仪器设备相比,成本大大降低。
(3)具有良好的经济效益,传统的测量手段在取样、制样、测定等方面需要耗费大量的人力、财力、物力,本测量方法简单、使用方便,可以快速、准确的测量植物15N示踪同位素丰度,故具有良好的经济效益。
附图说明
图1是本发明方法流程图;
具体实施方式
以下结合附图对本发明作进一步说明。
如图1所示,本发明方法具体包括以下步骤:
一、建立中红外光谱与植物15N示踪同位素丰度间的定量关系模型
①.15N同位素示踪下茶树的栽培。以茶树品种龙井43为研究对象,采用盆栽法,每盆栽4株1年生无性系茶苗,栽后修剪,进行正常的水分管理和病虫防治及防冻措施。用15N标记尿素为氮源,采用丰度分别为8%,5%,4%和2%尿素形成4个水平的氮示踪梯度和一个空白试验对照(即施用非15N标记的普通尿素),每个梯度设计处理4盆样品,每盆搭配施入0.6g磷和0.6g钾,肥料于春季分两次施下,间隔1月。
②.收集15N示踪茶树样本。施肥后30天、60天、90天、120天分别取样,每次取各盆中一株,4次重复。取样时,按盆连同根系取出整株植株,并拣净落在土壤中的根系,然后分株按根、茎、叶等器官取样,共得到67个样本。于80℃烘箱内烘至恒重,在红外灯照射下用玛瑙研钵研磨成粒径小于20μm的粉末,然后按1:49的比例把粉末与溴化钾磨细混匀,用钢铲把磨细混匀的粉末约0.05g移入压片模具,把压片模具水平放置在压片机座上,施加20t/cm2的压力,得到厚度约2mm的样本薄片。
③.获取样本的中红外吸收光谱。把制备好的样本薄片固定在支架上,放入中红外光谱仪的样本室,中红外光谱仪的光谱采集范围是4000cm-1~400cm-1,光谱分辨率为1cm-1,光谱采集时室温为25℃左右,湿度小于65%,每个样本的光谱是64次扫描的平均值,采集样本的透射光谱。把样本的透射光谱,按照公式“吸光度=2-log(透过率%)”的数学变换得到样本的中红外吸收光谱如图2所示。随后用质谱法测量样本的15N同位素丰度。
④.基于粒子群优化法捕捉植物N15同位素丰度的中红外光谱特征波数。基于粒子群优化法对示踪植物的N15同位素丰度的特征波数进行优化和筛选,由于4000cm-1~400cm-1范围共3601个波数变量,故设定对应的优化矩阵为1行3601列的二进制代码矩阵,优化矩阵中向量取值为0或1,0代表所对应的波数为非特征波数,1代表所对应的波数为特征波数。优化矩阵中所有为1的向量对应的波数成为特征波数组合,以这些特征波数的吸光度作为自变量,N15同位素丰度作为因变量,建立多元线性回归模型。粒子群优化的目标函数是多元线性回归模型的均方根误差最小。当均方根误差最小时对应粒子群优化法寻找得到的最优特征波数。研究发现,对于植物N15同位素丰度测量的特征波数为1235cm-1~1071cm-1和3475cm-1~3465cm-1。
⑤.建立样本中红外光谱与15N同位素丰度之间的定量关系模型。基于样本的中红外吸收光谱,获得样本在特征波数处(1235cm-1~1071cm-1和3475cm-1~3465cm-1)的吸光度,采用11点平均法提取得到16个特征变量,依据朗伯-比尔定律,应用多元线性回归算法建立这16个特征变量与15N同位素丰度的定量关系模型。模型的数学表达式如下所示:
Y=3.676e-03+3.543e-08x1+1.080e-07x2-4.259e-08x3-1.643e-07x4-2.913e-08x5
+3.106e-07x6-5.966e-08x7-3.084e-08x8-3.130e-07x9+3.647e-07x10-3.066e-07x11
+2.091e-08x12+9.052e-08x13+5.094e-08x14-9.086e-08x15+5.143e-08x16
其中:Y代表模型对样本的15N同位素丰度的预测值;
x1代表样本在1235cm-1~1225cm-1的11个波数处吸光度的平均值;
x2代表样本在1224cm-1~1214cm-1的11个波数处吸光度的平均值;
x3代表样本在1213cm-1~1203cm-1的11个波数处吸光度的平均值;
x4代表样本在1202cm-1~1192cm-1的11个波数处吸光度的平均值;
x5代表样本在1191cm-1~1181cm-1的11个波数处吸光度的平均值;
x6代表样本在1180cm-1~1170cm-1的11个波数处吸光度的平均值;
x7代表样本在1169cm-1~1159cm-1的11个波数处吸光度的平均值;
x8代表样本在1158cm-1~1148cm-1的11个波数处吸光度的平均值;
x9代表样本在1147cm-1~1137cm-1的11个波数处吸光度的平均值;
x10代表样本在1136cm-1~1126cm-1的11个波数处吸光度的平均值;
x11代表样本在1125cm-1~1115cm-1的11个波数处吸光度的平均值;
x12代表样本在1114cm-1~1104cm-1的11个波数处吸光度的平均值;
x13代表样本在1103cm-1~1093cm-1的11个波数处吸光度的平均值;
x14代表样本在1092cm-1~1082cm-1的11个波数处吸光度的平均值;
x15代表样本在1081cm-1~1071cm-1的11个波数处吸光度的平均值;
x16代表样本在3475cm-1~3465cm-1的11个波数处吸光度的平均值。二、基于植物15N示踪同位素丰度的中红外光谱测量模型对未知样本的15N示踪同位素丰度进行测定
基于所建立的植物15N示踪同位素丰度的中红外光谱测量模型对67个样本的预测相关系数为0.719,均方根误差为2.752e-09,标准均方根误差为7.488e-05%。具体预测结果如表1所示。说明该测量模型能够实现植物15N示踪同位素丰度的快速检测。
表1、模型对于67个样本的预测结果
样本号 | 真实值 | 预测值 |
(1) | 0.003675989106297490 | 0.003675988578237590 |
(2) | 0.003675990036688740 | 0.003675994336605070 |
(3) | 0.003675984906964000 | 0.003675988879986110 |
(4) | 0.003675997027195990 | 0.003675994713790710 |
(5) | 0.003675994940102100 | 0.003675992827862500 |
(6) | 0.003675987949594850 | 0.003675990740768610 |
(7) | 0.003675984906964000 | 0.003675986465997990 |
(8) | 0.003675987019203600 | 0.003675986818037930 |
(9) | 0.003675990967079990 | 0.003675993833690880 |
(10) | 0.003675994940102100 | 0.003675990539602930 |
(11) | 0.003675986088812350 | 0.003675990112125870 |
(12) | 0.003675993054173880 | 0.003675993355922400 |
(13) | 0.003675994940102100 | 0.003675986893475050 |
(14) | 0.003675989106297490 | 0.003675989408046010 |
(15) | 0.003675994940102100 | 0.003675994437187910 |
(16) | 0.003675993054173880 | 0.003675996071659030 |
(17) | 0.003675994940102100 | 0.003675995241850610 |
(18) | 0.003675994009710850 | 0.003675993984565140 |
(19) | 0.003675996096804740 | 0.003675994261167940 |
(20) | 0.003675994009710850 | 0.003675992098636920 |
(21) | 0.003675996096804740 | 0.003675995795056220 |
(22) | 0.003675997957587240 | 0.003675993909128010 |
(23) | 0.003676000044681130 | 0.003675995744764800 |
(24) | 0.003675993054173880 | 0.003675989282317460 |
(25) | 0.003675994940102100 | 0.003675992098636920 |
(26) | 0.003675996096804740 | 0.003675992450676860 |
(27) | 0.003675994940102100 | 0.003675994562916450 |
(28) | 0.003675991897471250 | 0.003675993531942360 |
(29) | 0.003675994009710850 | 0.003675994889810680 |
(30) | 0.003675997027195990 | 0.003675998284481460 |
(31) | 0.003675997027195990 | 0.003675997932441530 |
(32) | 0.003675997027195990 | 0.003675998359918590 |
(33) | 0.003675998887978490 | 0.003675995116122060 |
(34) | 0.003675993054173880 | 0.003675991847179830 |
(35) | 0.003675993054173880 | 0.003675993582233780 |
(36) | 0.003675991897471250 | 0.003675993481650950 |
(37) | 0.003675994940102100 | 0.003675993909128010 |
(38) | 0.003675990036688740 | 0.003675994814373550 |
(39) | 0.003676000975072380 | 0.003676001754589370 |
(40) | 0.003675997027195990 | 0.003675993682816620 |
(41) | 0.003676000044681130 | 0.003675999089144170 |
(42) | 0.003675996096804740 | 0.003675996247678990 |
(43) | 0.003675991897471250 | 0.003675992199219760 |
(44) | 0.003675998887978490 | 0.003675996297970410 |
(45) | 0.003676000044681130 | 0.003676001050509510 |
(46) | 0.003675994009710850 | 0.003675996826030310 |
(47) | 0.003675997027195990 | 0.003675996323116120 |
(48) | 0.003675997957587240 | 0.003675994738936420 |
(49) | 0.003675991897471250 | 0.003675992048345500 |
(50) | 0.003675993054173880 | 0.003675996976904570 |
(51) | 0.003675993054173880 | 0.003675992199219760 |
(52) | 0.003675990967079990 | 0.003675995719619090 |
(53) | 0.003675986088812350 | 0.003675991998054080 |
(54) | 0.003675987019203600 | 0.003675989508628840 |
(55) | 0.003675990967079990 | 0.003675991822034120 |
(56) | 0.003675990967079990 | 0.003675992727279660 |
(57) | 0.003675990036688740 | 0.003675990891642870 |
(58) | 0.003675994940102100 | 0.003675994286313650 |
(59) | 0.003675994009710850 | 0.003675995820201930 |
(60) | 0.003675991897471250 | 0.003675998862832780 |
(61) | 0.003675993054173880 | 0.003675994060002260 |
(62) | 0.003675991897471250 | 0.003675991092808540 |
(63) | 0.003675986088812350 | 0.003675989156588910 |
(64) | 0.003675994009710850 | 0.003675993607379490 |
(65) | 0.003675987949594850 | 0.003675986365415150 |
(66) | 0.003675994940102100 | 0.003675995870493350 |
(67) | 0.003675990967079990 | 0.003675987974740560 |
Claims (1)
1.基于中红外光谱的植物15N示踪同位素丰度快速检测方法,其特征在于该方法包括以下步骤:
步骤1.把待测植物样本烘干,然后在红外灯照射下把样本放入玛瑙研钵研磨成粒径小于20μm的粉末,按照质量比1:20~49的比例与干燥的溴化钾粉末混匀,使用压片机把混匀后的粉末制成约2mm的薄片,完成样本制备;
步骤2.把制备好的样本固定在支架上,放入中红外光谱仪的样品室,进行波数4000cm-1~400cm-1范围,分辨率为1cm-1的中红外透射光谱扫描,获取样本的中红外透射光谱;
步骤3.把样本的中红外透射光谱转换成吸光度,基于特征波数1235cm-1~1071cm-1和3475cm-1~3465cm-1处的吸光度,采用11点平均法提取得到16个特征变量,把这16个特征变量代入下式,即可计算得到植物样本的15N示踪同位素丰度Y,
Y=3.676e-03+3.543e-08x1+1.080e-07x2-4.259e-08x3-1.643e-07x4-2.913e-08x5
+3.106e-07x6-5.966e-08x7-3.084e-08x8-3.130e-07x9+3.647e-07x10-3.066e-07x11
+2.091e-08x12+9.052e-08x13+5.094e-08x14-9.086e-08x15+5.143e-08x16
其中:Y代表模型对样本的15N示踪同位素丰度的预测值;
x1代表样本在1235cm-1~1225cm-1的11个波数处吸光度的平均值;
x2代表样本在1224cm-1~1214cm-1的11个波数处吸光度的平均值;
x3代表样本在1213cm-1~1203cm-1的11个波数处吸光度的平均值;
x4代表样本在1202cm-1~1192cm-1的11个波数处吸光度的平均值;
x5代表样本在1191cm-1~1181cm-1的11个波数处吸光度的平均值;
x6代表样本在1180cm-1~1170cm-1的11个波数处吸光度的平均值;
x7代表样本在1169cm-1~1159cm-1的11个波数处吸光度的平均值;
x8代表样本在1158cm-1~1148cm-1的11个波数处吸光度的平均值;
x9代表样本在1147cm-1~1137cm-1的11个波数处吸光度的平均值;
x10代表样本在1136cm-1~1126cm-1的11个波数处吸光度的平均值;
x11代表样本在1125cm-1~1115cm-1的11个波数处吸光度的平均值;
x12代表样本在1114cm-1~1104cm-1的11个波数处吸光度的平均值;
x13代表样本在1103cm-1~1093cm-1的11个波数处吸光度的平均值;
x14代表样本在1092cm-1~1082cm-1的11个波数处吸光度的平均值;
x15代表样本在1081cm-1~1071cm-1的11个波数处吸光度的平均值;
x16代表样本在3475cm-1~3465cm-1的11个波数处吸光度的平均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110196570 CN102313713B (zh) | 2011-07-14 | 2011-07-14 | 基于中红外光谱的植物15n示踪同位素丰度快速检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110196570 CN102313713B (zh) | 2011-07-14 | 2011-07-14 | 基于中红外光谱的植物15n示踪同位素丰度快速检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102313713A CN102313713A (zh) | 2012-01-11 |
CN102313713B true CN102313713B (zh) | 2013-04-03 |
Family
ID=45427070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110196570 Active CN102313713B (zh) | 2011-07-14 | 2011-07-14 | 基于中红外光谱的植物15n示踪同位素丰度快速检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102313713B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102928396B (zh) * | 2012-10-29 | 2014-10-22 | 浙江大学 | 基于拉曼光谱的尿素同位素丰度的快速检测方法 |
CN103267740A (zh) * | 2012-12-20 | 2013-08-28 | 江苏大学 | 秸秆发酵过程特征波数软仪表装置及其构造方法 |
CN105424640A (zh) * | 2015-11-27 | 2016-03-23 | 浙江大学 | 一种茶叶中美术绿添加含量的检测方法 |
CN106417194A (zh) * | 2016-12-09 | 2017-02-22 | 广州微因生物科技有限公司 | 一种用于昆虫蛋白质定量和示踪的n15稳定同位素标记方法 |
CN109580413B (zh) * | 2017-09-28 | 2021-04-23 | 宁海德宝立新材料有限公司 | 一种二元混合物的红外光谱分析方法及其应用 |
CN107741417B (zh) * | 2017-09-28 | 2019-11-08 | 上海氘峰医疗器械有限公司 | 一种原位快速检测细胞生物过程的方法 |
CN112305007B (zh) * | 2020-10-30 | 2022-12-09 | 上海化工研究院有限公司 | 一种利用核磁氢谱或氘谱测定氘标记化合物的氘同位素丰度的方法 |
CN113324870B (zh) * | 2021-06-18 | 2022-02-01 | 河北省农林科学院粮油作物研究所 | 一种定量作物非根器官氮素吸收及其氮营养贡献的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2312743A (en) * | 1996-05-01 | 1997-11-05 | Japan Radio Co Ltd | Analysis of carbon isotopes |
CN101303279A (zh) * | 2008-06-26 | 2008-11-12 | 华东理工大学 | 用于近红外检测低丰度组分的在线富集装置及其使用方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5696378A (en) * | 1996-03-20 | 1997-12-09 | Baylor University | High accuracy determination of chlorine content by isotope dilution flame infrared emission spectrometry (ID-FIRE) |
-
2011
- 2011-07-14 CN CN 201110196570 patent/CN102313713B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2312743A (en) * | 1996-05-01 | 1997-11-05 | Japan Radio Co Ltd | Analysis of carbon isotopes |
CN101303279A (zh) * | 2008-06-26 | 2008-11-12 | 华东理工大学 | 用于近红外检测低丰度组分的在线富集装置及其使用方法 |
Non-Patent Citations (6)
Title |
---|
Standardizing Infra-red Measures of Bone Mineral Crystallinity: an Experimental Approach;Todd A. Surovell 等;《Journal of Archaeological Science》;20011231;第28卷;第633-642页 * |
Todd A. Surovell 等.Standardizing Infra-red Measures of Bone Mineral Crystallinity: an Experimental Approach.《Journal of Archaeological Science》.2001,第28卷第633-642页. |
唐强 等.基于高光谱反射率的冬小麦生长后期氮素丰度监测研究.《光谱学与光谱分析》.2010,第30卷(第11期),第3061-3066页. |
基于高光谱反射率的冬小麦生长后期氮素丰度监测研究;唐强 等;《光谱学与光谱分析》;20101130;第30卷(第11期);第3061-3066页 * |
我国各地氮同位素丰度的质谱测量;王建耀 等;《化学通报》;19881231(第8期);第34-35页 * |
王建耀 等.我国各地氮同位素丰度的质谱测量.《化学通报》.1988,(第8期),第34-35页. |
Also Published As
Publication number | Publication date |
---|---|
CN102313713A (zh) | 2012-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102313713B (zh) | 基于中红外光谱的植物15n示踪同位素丰度快速检测方法 | |
CN102426153B (zh) | 一种基于冠层高光谱指数的小麦植株水分监测方法 | |
CN102928396B (zh) | 基于拉曼光谱的尿素同位素丰度的快速检测方法 | |
CN104502283A (zh) | 估测大豆产量和地上部干重的一套两波段高光谱指数和预测模型 | |
CN110567892B (zh) | 一种基于临界氮浓度的夏玉米氮素高光谱预测方法 | |
CN103293111A (zh) | 一种土壤背景干扰下小麦叶层氮含量光谱监测模型及建模方法 | |
CN105445229A (zh) | 基于氮平衡光谱法的小麦春季氮肥施用方法及其追氮量模型的构建方法 | |
CN103954567A (zh) | 基于连续统去除法的土壤盐分测定方法 | |
Li et al. | Use of fluorescence-based sensors to determine the nitrogen status of paddy rice | |
CN101762569A (zh) | 一种畜禽粪便工厂化堆肥发酵过程的无损监测方法 | |
CN105784628A (zh) | 一种利用中红外光谱检测土壤有机质化学组成的方法 | |
CN101836561A (zh) | 一种海滨盐土蓖麻产量预测模型及其构建方法和应用 | |
Wang et al. | Solar-induced chlorophyll fluorescence as an indicator for determining the end date of the vegetation growing season | |
CN107748144A (zh) | 快速测定土壤碳氮及其稳定同位素的中红外光谱检测系统 | |
CN113670913B (zh) | 水稻氮素含量反演高光谱植被指数构建方法 | |
CN105067666A (zh) | 一种利用热分析检测植物生物质toc的方法 | |
Tang et al. | Soil respiration and carbon balance in a Moso bamboo (Phyllostachys heterocycla (Carr.) Mitford cv. Pubescens) forest in subtropical China | |
CN103674929B (zh) | 光谱分析用于植物幼苗的体检方法 | |
CN103575680A (zh) | 一种评估有机肥质量指标的光谱学方法 | |
CN112362812A (zh) | 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 | |
AU2018101606A4 (en) | A method for identifying meconopsis quintuplinervia regel from different geographical origins | |
CN111595806A (zh) | 一种利用中红外漫反射光谱监测土壤碳组分的方法 | |
CN102230894B (zh) | 15n示踪尿素同位素丰度的红外光谱快速检测方法 | |
Gao et al. | Phylogeography of Ulmus elongata based on Fourier transform-infrared spectroscopy (FTIR), thermal gravimetric and differential thermal analyses | |
Liu et al. | Adaptive mechanism exploration of camellia limonia in karst high calcium environment by integrative analysis of metabolomics and metagenomics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130403 |