本发明涉及产地鉴别，具体为一种基于稳定同位素比值鉴别白术产地的方法。
背景技术：
：白术为菊科植物白术(atractylodesmacrocephalakoidz.)的干燥根茎，是著名中药“浙八味”之一，具有健脾益气、燥湿利水、止汗、安胎的功效，为常用补益类中药，俗有“北参南术”之说。白术主产于浙江、安徽、湖南、湖北等省，江西、福建、四川等省亦有栽培，全国范围内栽培白术的县市超过了200个，但品质良莠不齐。自宋代以来，浙江於潜、新昌、磐安等地由于出产的白术品质优良，逐渐发展为白术道地产区。由于道地药材自古以来就被视为优质药材，其价格通常高于非道地药材，受经济利益驱使，目前市场上充斥着假冒的道地药材，大大影响了我国中医药行业的健康发展，损害了消费者的利益。因此，对白术进行原产地鉴别意义重大。白术产地鉴别一直以来是一大难题，传统的人工鉴别法需要依赖专家的个人经验，无法大规模推广，目前常用的分子生物学方法和化学成分分析法(如：dna条形码技术、高效液相色谱指纹图谱、气相色谱指纹图谱、红外光谱指纹图谱)各有局限。其中分子生物学方法能准确地区分某种特定中药材内部不同种质间的微小差异，但若种源相同则无法区分，而非道地产区的白术其种源很可能是来自道地产区，在这种情况下分子标记技术将无能为力。化学成分分析法的准确性不够高，现有的报道显示不同地区所产白术在一些化学成分上有时会存在差异，但却无法直接根据这些数据推测白术原产地，且化学成分容易受产地以外的因素(如：种源)的影响。因此，我们亟需一种直接、可靠的方法对白术原产地进行鉴定。技术实现要素：本发明提供了一种基于稳定同位素比值鉴别白术产地的方法，该方法具有判别准确率高、检测指标少、样品前处理简单以及样品需要量少的优点。为实现上述目的，本发明提供如下技术方案：一种基于稳定同位素比值鉴别白术产地的方法，该方法包括如下步骤：s01、选取多组不同真实产地的白术样本；s02、测定每个所述白术样本的稳定同位素比值a，构建真实产地白术样本的稳定同位素比值数据库；其中，所述的稳定同位素比值a包括：δ13c、δ15n、δ2h、δ18o以及δ34s中的至少一种；s03、根据真实产地白术样本的稳定同位素比值数据库，利用统计软件对其进行正交偏最小二乘法判别分析，建立产地溯源模型；s04、选取未知产地的待测白术样本，测定得到待测白术样本的稳定同位素比值b；其中，所述的稳定同位素比值b包括：δ13c’、δ15n’、δ2h’、δ18o’以及δ34s’中的至少一种；s05、将上述得到的稳定同位素比值b作为输入值，输入至产地溯源模型，输出得到相应的产地判别结果。其中，构成生物体的主要元素有碳、氢、氧、氮、硫等，来自不同地区的农产品，其稳定同位素比值(如：δ13c,δ15n,δ2h,δ18o和δ34s)受气候条件和地理环境影响会有所不同，这种差异源自于当地的光照、温度、大气、降水、土壤等自然条件，是农产品自带的天然指纹标签，不受农产品后期加工的影响，能真实地反映它们的原产地。作为优选，在步骤s01和s02之间，还包括步骤s01.1，所述步骤s01.1为：(1)将白术样本在50℃～70℃的温度条件下放置5～7h，冷却至室温后，碾碎，过50～80目筛，得到白术样本粉末；(2)将所述的白术样本粉末用锡杯包好，置于元素分析-稳定同位素质谱仪的样品盘上即可用于步骤s02对稳定同位素比值a的测定。作为优选，步骤(1)具体为：将白术样本在60℃的温度条件下放置6h，冷却至室温后，碾碎，过60目筛，得到白术样本粉末。作为优选，用于测定稳定同位素比值a的仪器为元素分析-稳定同位素质谱仪。作为优选，在稳定同位素比值a中，δ13c、δ15n以及δ34s的测试条件为：元素分析-稳定同位素质谱仪的氦气吹扫流量为：230ml/min；氧化炉和还原炉的温度分别为1150℃和850℃；进入元素分析-稳定同位素质谱仪氦载气流量为100ml/min；在稳定同位素比值b中,δ13c’、δ15n’以及δ34s’的测试条件同上。作为优选，在稳定同位素比值a中，δ2h以及δ18o的测试条件为：元素分析-稳定同位素质谱仪的氦气吹扫流量为200ml/min；反应管温度为1400℃，进入元素分析-稳定同位素质谱仪氦载气流量为110ml/min；在稳定同位素比值b中，δ2h’以及δ18o’的测试条件同上。作为优选，在步骤s05中，统计软件选用simca-p。本发明的有益效果为：本发明通过检测碳、氢、氧、氮、硫五种元素的稳定同位素比值，结合正交偏最小二乘法判别分析，建立了一种白术产地鉴别方法，具有判别准确率高，检测指标少，样品前处理简单，样品需要量少的优点。且相较于传统的人工鉴定法，本方法无需依赖专家的个人经验，可操作性强，有利于推广。具体实施方式本实施例提供一种技术方案：一种基于稳定同位素比值鉴别白术产地的方法，该方法包括如下步骤：s01、选取多组不同真实产地的白术样本，样本数量分别为：浙江11个、陕西8个以及河北9个；s02、测定已知产地的白术样本的稳定同位素比值a，构建真实产地白术样本的稳定同位素比值数据库，见表1。其中，所述的稳定同位素比值a包括：δ13c、δ15n、δ2h、δ18o以及δ34s。s03、根据真实产地白术样本的稳定同位素比值数据库，利用统计软件对其进行正交偏最小二乘法判别分析，建立产地溯源模型。经七折交叉验证后，该模型的判别正确率见表2。s04、选取未知产地的待测白术样本，测定得到待测白术样本的稳定同位素比值b，见表3。其中，稳定同位素比值b包括：δ13c’、δ15n’、δ2h’、δ18o’以及δ34s’。s05、将s04得到的稳定同位素比值b作为输入值，输入至s03的产地溯源模型，输出得到相应的产地判别结果。表1稳定同位素比值a编号产地δ13cδ15nδ2hδ18oδ34s1陕西-25.901.92-50.2824.6510.152陕西-25.182.75-53.5625.026.563陕西-25.14-1.29-50.4825.1612.864陕西-25.601.87-50.8924.818.885陕西-24.824.24-59.1623.837.976陕西-25.302.78-63.7021.987.577陕西-25.02-1.26-66.2222.655.388陕西-25.203.92-60.3323.063.479河北-26.17-0.90-56.7826.896.1510河北-25.67-0.85-55.1026.453.9511河北-25.680.65-56.1624.976.8612河北-25.49-0.64-52.4627.027.6213河北-26.030.31-52.5825.548.5814河北-24.371.36-44.8024.63-1.9515河北-24.891.83-41.7425.481.7416河北-25.450.28-42.7625.59-1.0217河北-24.381.20-43.6825.930.9218浙江-24.374.69-45.3322.78-4.4119浙江-24.157.87-50.2524.27-2.8120浙江-23.723.48-45.8426.01-0.7621浙江-23.412.95-43.6725.16-0.3322浙江-24.443.02-42.8023.76-0.2923浙江-23.584.58-49.5424.370.6024浙江-24.144.59-49.4424.12-4.2725浙江-23.904.70-46.0424.16-1.4826浙江-24.181.27-56.9723.43-0.7427浙江-24.702.35-51.1823.57-2.1028浙江-24.593.15-52.1223.863.27表2白术产地溯源模型判别结果在本实施例中，通过对不同真实产地白术样本的稳定同位素比值数据库进行正交偏最小二乘法判别分析，建立白术产地溯源模型，经七折交叉验证后，该模型对浙江省、陕西省和河北省三个产地白术的判别正确率分别达到了100％、100％和88.89％，整体判别正确率达到96.43％，具体结果见表2。需要说明的是，在实际检验过程中，样本数量非常多，本实施例仅仅截取部分。本例中虽有1个河北样本被错判为浙江，但该模型的整体判别正确率相比其他现有手段仍是相当高的，几乎接近100％，因此可准确用于白术产地的溯源。在上述得到的白术产地溯源模型可用于准确追溯白术产地的基础上，实验人员做下一步实例：选取未知真实产地的白术盲样2个，测得其稳定同位素比值(即稳定同位素比值b)如表3所示，将其输入至s03中建立的产地溯源模型，对两组数据进行产地判别分析，结果显示，第一个样本“未知-1”的产地为浙江，第二个样本“未知-2”的产地为陕西。因此，实验人员采用上述方法判别白术的原产地时，具有判别准确率高，检测指标少以及样品前处理简单的优点。可适用于对批量样本的判别，提高实验人员的效率。表3稳定同位素比值b编号δ13cδ15nδ2hδ18oδ34s未知-1-24.433.94-50.4424.48-0.71未知-2-25.340.88-55.8622.6710.57当前第1页12