CN115718081B - 基于光谱指纹的琥珀产地溯源模型的构建方法及应用 - Google Patents
基于光谱指纹的琥珀产地溯源模型的构建方法及应用 Download PDFInfo
- Publication number
- CN115718081B CN115718081B CN202211429390.5A CN202211429390A CN115718081B CN 115718081 B CN115718081 B CN 115718081B CN 202211429390 A CN202211429390 A CN 202211429390A CN 115718081 B CN115718081 B CN 115718081B
- Authority
- CN
- China
- Prior art keywords
- amber
- vibration
- sample
- functional group
- origin
- 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
- 230000003595 spectral effect Effects 0.000 title claims abstract description 22
- 238000010276 construction Methods 0.000 title claims abstract description 20
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 38
- 125000000524 functional group Chemical group 0.000 claims abstract description 33
- 238000004458 analytical method Methods 0.000 claims abstract description 28
- 238000004519 manufacturing process Methods 0.000 claims abstract description 19
- 238000000513 principal component analysis Methods 0.000 claims abstract description 15
- 238000010606 normalization Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 150000001735 carboxylic acids Chemical class 0.000 claims description 24
- 238000010521 absorption reaction Methods 0.000 claims description 23
- 238000012360 testing method Methods 0.000 claims description 21
- 238000000103 photoluminescence spectrum Methods 0.000 claims description 19
- 238000000862 absorption spectrum Methods 0.000 claims description 15
- 150000003900 succinic acid esters Chemical class 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 11
- 229910052760 oxygen Inorganic materials 0.000 claims description 11
- 230000005284 excitation Effects 0.000 claims description 10
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 8
- 239000001301 oxygen Substances 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000005452 bending Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 238000005498 polishing Methods 0.000 claims description 5
- 238000003672 processing method Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 125000002496 methyl group Chemical group [H]C([H])([H])* 0.000 claims description 4
- 125000000325 methylidene group Chemical group [H]C([H])=* 0.000 claims description 4
- KDYFGRWQOYBRFD-UHFFFAOYSA-L succinate(2-) Chemical compound [O-]C(=O)CCC([O-])=O KDYFGRWQOYBRFD-UHFFFAOYSA-L 0.000 claims description 4
- 238000010998 test method Methods 0.000 claims description 4
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 3
- 229920001971 elastomer Polymers 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 2
- 238000012417 linear regression Methods 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 abstract description 11
- 238000010183 spectrum analysis Methods 0.000 abstract description 2
- 238000009825 accumulation Methods 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 40
- 230000009102 absorption Effects 0.000 description 21
- 238000005516 engineering process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000004566 IR spectroscopy Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 241000282376 Panthera tigris Species 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000002290 gas chromatography-mass spectrometry Methods 0.000 description 2
- 238000004949 mass spectrometry Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 239000000025 natural resin Substances 0.000 description 2
- 230000003647 oxidation Effects 0.000 description 2
- 238000007254 oxidation reaction Methods 0.000 description 2
- IOLCXVTUBQKXJR-UHFFFAOYSA-M potassium bromide Chemical compound [K+].[Br-] IOLCXVTUBQKXJR-UHFFFAOYSA-M 0.000 description 2
- 239000002994 raw material Substances 0.000 description 2
- 239000013074 reference sample Substances 0.000 description 2
- 238000001028 reflection method Methods 0.000 description 2
- 239000002893 slag Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 235000001674 Agaricus brunnescens Nutrition 0.000 description 1
- 101100313164 Caenorhabditis elegans sea-1 gene Proteins 0.000 description 1
- 235000008694 Humulus lupulus Nutrition 0.000 description 1
- 244000025221 Humulus lupulus Species 0.000 description 1
- 244000137852 Petrea volubilis Species 0.000 description 1
- 238000001237 Raman spectrum Methods 0.000 description 1
- 244000028419 Styrax benzoin Species 0.000 description 1
- 235000000126 Styrax benzoin Nutrition 0.000 description 1
- 235000008411 Sumatra benzointree Nutrition 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 239000002775 capsule Substances 0.000 description 1
- 125000002915 carbonyl group Chemical group [*:2]C([*:1])=O 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 229930004069 diterpene Natural products 0.000 description 1
- 125000000567 diterpene group Chemical group 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 150000002148 esters Chemical class 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000004817 gas chromatography Methods 0.000 description 1
- 239000010437 gem Substances 0.000 description 1
- 229910001751 gemstone Inorganic materials 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 235000019382 gum benzoic Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- -1 ivory Substances 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000004811 liquid chromatography Methods 0.000 description 1
- 238000004895 liquid chromatography mass spectrometry Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000001819 mass spectrum Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000000643 oven drying Methods 0.000 description 1
- 230000001590 oxidative effect Effects 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000000045 pyrolysis gas chromatography Methods 0.000 description 1
- 238000000009 pyrolysis mass spectrometry Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000010301 surface-oxidation reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了基于光谱指纹的琥珀产地溯源模型的构建方法,取不同产地的琥珀样品,将样品进行预处理;将处理后的样品进行红外光谱采集;对得到的红外光谱采用OPUS软件进行谱线归一化处理,然后采用Origin软件进行曲线拟合,获得峰位;确定峰位代表的官能团,运用主成分分析和判别分析对所得官能团数据进行分析,建立产地溯源模型。本发明首先收集不同国家的琥珀样品,经预处理和多种光谱分析,对样品光谱进行曲线拟合和官能团峰位识别;结合主成分分析和判别分析,主成分析得到八个主成分因子特征值,总贡献率的得分值累积达到93.03%,以此建立了一种出土琥珀产地鉴别的方法,该方法快速、准确,准确率高达100%。
Description
技术领域
本发明属于出土琥珀科技考古以及珠宝检测技术领域,具体涉及基于光谱指纹的琥珀产地溯源模型的构建方法及应用。
背景技术
琥珀属于有机宝石的一种,是天然树脂因地壳运动埋藏于地下,经地质作用转化为的树脂化石。琥珀的形成历经几千万年至上亿年,是记录地史的“时间胶囊”。与此同时,琥珀作为古今中外人类喜爱的珍宝,也反映了不同时期的东西方文化。作为人类最早的艺术媒介之一,不同民族的人们赋予了琥珀许多神秘的传说和美好的寄托。西方的古希腊神话中记载琥珀是由泪珠凝结而成的。汉代时人们认为琥珀是“虎目光沦入地所为也”,称之为虎魄。琥珀曾同金、象牙、玉等一样为皇室和权贵的奢侈品,琥珀制成的装饰品、工艺品及实用器具深受追捧和喜爱,其加工工艺也不断精进,在东西方都形成了独特的琥珀艺术。古代琥珀重要产地主要集中在欧洲波罗的海沿岸国家和亚洲的缅甸。据历代古籍记载,中国古代琥珀原料的产地包括:永昌郡(今云南、缅甸交界地区)、罽宾(今克什米尔)、安息、波斯、大秦、波罗的海等。琥珀原料从邻国甚至欧洲国家抵达古代中国的过程中,不仅涉及了贸易和交通往来,还留下了东西方琥珀工艺及文化交流的线索。
目前国内主要应用形态分析、放大观察和红外光谱技术对少量出土琥珀进行材质鉴定和溯源。覃春雷等应用形态分析、放大观察和无损的红外光谱测试对4件古代琥珀珠饰开展年代、材质和产地鉴定,初步判定4件样品年代为西汉到魏晋南北朝,且产地为波罗的海。罗武干等采用红外光谱指出河南省南阳市体育中心游泳馆汉画像石墓M18中出土的3件琥珀的产源与缅甸琥珀较为相似,推测琥珀样品属于缅甸汉代文化互动与交流的产物。国外关于出土琥珀产地溯源研究的手段相对较多,无损技术常用红外光谱、拉曼光谱,有损技术多采用了气相色谱/质谱。Colombini等利用裂解气相色谱/质谱及红外光谱对俄罗斯、波兰等博物馆所藏的20块琥珀藏品进行了产地来源和降解机制的研究,研究表明这些琥珀藏品来自波罗的海,外层多孔粗糙的外观可能与较高的浸出率和游离二萜从表面蒸发有关。
通过以上分析,现有技术存在的问题:
(1)采用单一谱学方法,没有进行多谱学方法的联用,测试得到的数据有限,产地判别困难,无法得出准确产地信息;
(2)气相色谱/质谱和液相色谱/质谱方法属于有损,不适合不可损坏的文物样品;
(3)在琥珀产地溯源的研究中,基本是通过人眼观察谱图进行判断,需要研究人员具备较为深厚的经验。而当琥珀经风化和氧化,其谱峰发生一定程度的飘移时,琥珀产地的区分就会显得较为困难,此时采用统计方法对琥珀产地进行分类研究更为科学和准确。
发明内容
本发明的目的在于,针对现有技术的上述不足,提供了基于光谱指纹的琥珀产地溯源模型的构建方法及应用。
为实现上述目的,本发明采用如下的技术方案:
本发明的第一目的提供一种基于光谱指纹的琥珀产地溯源模型的构建方法,包括以下具体步骤:
步骤S1,取不同产地的琥珀样品,将样品进行预处理,所述产地包括缅甸、波罗的海和中国辽宁抚顺;
步骤S2,将处理后的样品进行红外光谱、光致发光光谱和紫外可见吸收光谱数据采集分析;
步骤S3,采用OPUS软件对得到的红外光谱图进行谱线归一化处理,然后采用Origin软件进行曲线拟合,获得峰位;
步骤S4,确定峰位代表的官能团,运用主成分分析和判别分析对所得官能团数据进行分析,建立产地溯源模型;
所述产地溯源模型包括:根据每个样品的红外光谱,得到官能团数据,运用主成分分析初步判别影响较大的官能团数据;利用线性回归分析,判别官能团的多重共线性,VIF小于10视为不存在多重共线性;根据筛选后的不具有多重共线性的官能团数据进行判别分析;利用判别分析和交叉检验准确判别的不同产地琥珀的分布规律,
红外光谱特征信息指示的琥珀产地溯源的判别函数分别为:
Y1(x)=38.977x1-41.439x2+15.038x3+8.046x4-7.024x5+16.104x6-34.596x7+25.375x8;
Y2(x)=8.952x1+19.421x2-8.037x3+8.765x4-19.969x5-7.038x6+18.129x7+34.421x8;
其中,x1、x2、x3、x4、x5、x6、x7和x8代表八个主成分因子特征值,具体为总贡献率的得分值,x1与1737cm-1处琥珀中琥珀酸酯中C=O伸缩振动及1698cm-1处羧酸中C=O振动相关;x2与琥珀酸酯中1733cm-1处C=O伸缩振动及1693cm-1处羧酸中C=O振动的相关;x3与1710cm-1处羧酸中C=O振动及1378cm-1处甲基、亚甲基中的C-H弯曲振动相关;x4与1666cm-1处不饱和C=C伸缩振动及1162cm-1处含氧官能团C-O单键伸缩振动相关;x5与1737cm-1处琥珀酸酯中C=O伸缩振动及1704cm-1处羧酸中C=O振动的相关;x6与1716cm-1处羧酸中C=O振动与1658cm-1处不饱和C=C伸缩振动相关;x7与1627cm-1处不饱和C=C伸缩振动及1253cm-1处含氧官能团C-O单键伸缩振动;x8与1741cm-1处琥珀酸酯中C=O伸缩振动及1718cm-1处羧酸中C=O振动的相关。
进一步的,步骤S2中,采集的红外光谱数据范围为1800-400cm-1。根据不同产地测得的红外光谱,得到琥珀中对产地判别影响较大的官能团,选取对产地判别影响较大的红外光谱数据范围为1800-400cm-1。
进一步的,步骤S1中,所述样品预处理步骤包括:
步骤S11,对样品表面进行清理;
步骤S12,进行样品进行切割,制成厚度为5mm的薄片;
步骤S13,对薄片样品进行双面平行抛光。
进一步的,所述红外光谱采集条件为:背景扫描32次,样品扫描64次,分辨率为2cm-1,测试范围为400-4000cm-1。
进一步的,步骤S3中,所述红外光谱谱线处理方法为:采用红外处理软件OPUS中气氛补偿进行H2O和CO2补偿,采用凹式rubbenband校正方法进行基线校正,迭代次数为25,基线点数为100,校正范围选择为4000-400cm-1波数区域进行基线校正,使用OPUS软件的向量归一化方法进行谱线归一化处理。
本发明的第二目的是提供一种出土琥珀产地的鉴别方法,包括以下具体步骤:
步骤S1,获得待测出土琥珀样品在波长1800~400cm-1的红外光谱图;
步骤S2,按照权利要求1-5中任一项所述的构建方法,将步骤S1获得的红外光谱图带入所述构建方法中运算进行主成分分析,得到八个主成分x1、x2、x3、x4、x5、x6、x7和x8的得分值;
步骤S3,将x1~x8的得分值代入权利要求1~5中任一项所述的构建方法得到的产地溯源模型中,得到判别界值;
步骤S4,对比判别界值到组质心距离最小即可判别产地;
若当同时判别为波罗的海和缅甸,需采用光致发光光谱鉴别,其中在405nm激发下,波罗的海琥珀荧光特征峰位于500-550nm范围内,缅甸琥珀荧光特征峰位于450-500nm范围内;
若当同时判别为波罗的海和中国辽宁抚顺,需采用光致发光光谱鉴别,其中在405nm激发下,波罗的海琥珀荧光特征峰位于500-550nm范围内,中国辽宁抚顺琥珀荧光特征峰位于450-500nm范围内,
若当同时判别为缅甸和中国辽宁抚顺,需采用紫外可见吸收光谱鉴别,比较875nm-925nm范围内的是否吸收峰,有吸收峰的为中国辽宁抚顺琥珀,没有吸收峰的为缅甸琥珀。
进一步的,所述光致发光光谱采集条件为:激发光源波长405nm,测试范围210-1000nm,扫描时间5s,扫描次数3。
进一步的,所述紫外可见吸收光谱采集条件为:测试方法为透射法,测试范围为200-800nm,数据间隔1nm,探测器狭缝宽度2nm,扫描速度为266.75nm/min。
本发明的第三目的是提供一种电子设备,包括处理器和存储有计算机程序的存储器,所述计算机程序被所述处理器执行时,使得所述处理器执行上述的构建方法。
本发明的第四目的是提供一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,使得所述处理器执行上述的构建方法。
与现有技术比较,本发明提供的技术方案带来的有益效果是:
(1)本发明提供了基于光谱指纹的琥珀产地溯源模型的构建方法,取不同产地的琥珀样品,将样品进行预处理;将处理后的样品进行红外光谱采集;采用OPUS软件对得到的红外光谱图进行谱线归一化处理,然后采用Origin软件进行曲线拟合,获得峰位;确定峰位代表的官能团,运用主成分分析和判别分析对所得官能团数据进行分析,建立产地溯源模型。本发明首先收集不同国家的琥珀样品,经预处理和多种光谱分析,对样品光谱进行曲线拟合和官能团峰位识别;结合主成分分析、回归分析和判别分析,对波罗的海、缅甸和中国辽宁抚顺的琥珀进行原产地溯源,构建产地溯源模型,建立了一种出土琥珀产地的鉴别方法,该方法快速、准确,准确率高达100%;
(2)相较于气相色谱-质谱及液相色谱质谱等质谱方法,采用红外光谱、光致发光光谱和紫外可见吸收光谱方法联用可对完整样品进行无损测试,而对于已经碎裂的样品仅需要微量样品即可得到产地信息,该方法适用于出土琥珀可用于质检和文物检测中;
(3)本方法提出一种出土琥珀产地的鉴别方法,对样品的要求都较为简单,无需复杂的制样过程。结合本发明建立的琥珀产地溯源模型,和多元统计方法,极大的减弱了因氧化或老化导致的谱峰漂移带来的产地溯源困难。选取光谱学指纹参数和数学模型结合,可实现出土琥珀的快速、准确产地识别,为考证古代经济、文化和商贸交流情况打下基础。
附图说明
图1为本发明提供的基于光谱指纹的琥珀产地溯源模型的构建方法流程图;
图2为本方法实施例1提供的不同产地琥珀的红外光谱图;
图3为本方法实施例1提供的不同产地琥珀的光致发光光谱图;
图4为本方法实施例1提供的缅甸和抚顺琥珀的紫外可见吸收光谱图;
图5为本方法实施例1提供的判别分析二维散点图;
图6为本发明实施例2提供的汉代出土琥珀样品判别分析二维散点图;
图7为本发明实施例3提供的清代出土琥珀样品判别分析二维散点图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合具体实施例和附图,对本发明的具体实施方式作进一步详细描述。实施例中未注明具体技术或条件的,按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。所用试剂或仪器未注明生产厂商者,均为可以通过市购获得的常规产品。
本发明的出土琥珀是指墓葬中出土的琥珀文物。这些琥珀文物由于埋藏环境的影响,易受氧化和风化作用产生龟裂、碎裂,甚至渣化等氧化蚀变。不同程度的氧化蚀变均会对检测仪器产生信号噪声,致使峰形、峰位发生改变,给产地识别带来挑战。
实施例1
琥珀产地溯源模型的构建及出土琥珀产地的鉴别方法
步骤S1,采集波罗的海样品5个、缅甸6个和中国辽宁抚顺琥珀4个,样品采集信息表见表1所示;对样品表面进行清理,使用金刚粉线锯对所有样品进行切割,制成厚度为2~5mm左右的薄片,并借助砂纸和抛光板或抛光布对样品进行双面平行抛光。
表1.出土琥珀样品采集信息表
产地 | 产地代号 | 样品数量/个 |
波罗的海 | 1 | 10 |
缅甸 | 2 | 10 |
中国辽宁抚顺 | 3 | 10 |
步骤S2,将处理后的样品分别进行红外光谱、光致发光光谱和紫外可见吸收光谱数据采集
(1)红外光谱测试
采用傅里叶红外光谱仪(型号:德国BruckerTensor-27)对样品进行红外光谱测试。打开OPUS软件,选择反射模式,将样品抛光面放置于ATR样品台检测窗上,确保光束照射在样品表面。
设置红外光谱采集条件:背景扫描32次,样品扫描64次,分辨率为2cm-1,测试范围为400-4000cm-1。
分别测试得到三个不同产地各10个样品的谱图。
(2)光致发光光谱测试
采用显微光致发光光谱仪(型号:韩国SC11012)对样品进行测试。测试条件:激发光源波长405nm;测试范围200-1000nm;扫描时间5s;扫描次数3。测试方法为透射法,在暗室环境或在对探头和样品做遮光处理后测试,调节探头到样品距离,保证激发光聚焦在样品表面的光斑最小,以获得最强光谱信号。将样品放置于光束下,使激发光汇聚到样品表面,积分球接收到样品背面对应位置的透射光后即可得到该样品的光致发光光谱。
(3)紫外可见吸收光谱测试
使用PerkinElmer Lambda 650S型的紫外可见吸收光谱仪,测试方法选择反射法,测试范围200~800nm,数据间隔1nm,探测器狭缝宽度2nm,扫描速度为266.75nm/min。选择所制样品中透明度较高的部分进行测试。
步骤S3,采用OPUS和origin进行谱线处理,确定峰位及官能团归属
红外光谱谱线处理方法为:采用红外处理软件OPUS中气氛补偿进行H2O和CO2补偿,采用凹式rubbenband校正方法进行基线校正,迭代次数为25,基线点数为100,校正范围选择为4000-400cm-1波数区域进行基线校正。使用OPUS软件的向量归一化方法进行谱线归一化处理。利用Origin对所得谱图进行拟合,获得谱图峰位。
光致发光光谱处理方法为:对所得发光谱图利用origin绘制,得到特征发光谱图。
紫外可见吸收光谱处理方法为:对所得反射谱图利用origin进行拟合,获得吸收峰及反射峰。
步骤S4,运用多元统计分析,建立不同产地琥珀产地溯源模型
琥珀的红外光谱主要包括两个部分,4000-1500cm-1称为官能团区,1500-400cm-1称为指纹区。其中3000-2800cm-1为由C-H饱和键伸缩振动引起的谱带,-CH2、-CH3中的C-H弯曲振动引起的红外吸收峰位于1464、1380cm-1,C=O伸缩振动引起的红外吸收峰出现在1740-1695cm-1。1172-1272cm-1、1032cm-1范围内的红外吸收谱带归属为C-O伸缩振动所致,970cm-1附近的吸收峰归属于C-H面外弯曲振动。结合现有琥珀的红外光谱研究,并且利用origin谱峰拟合的谱峰,选取其中红外光谱中1800-400cm-1的红外光谱数据进行主成分分析。琥珀的产地鉴别特征红外吸收峰主要存在于1800-400cm-1范围内,红外光谱中1750-1690cm-1之间为C=O官能团伸缩振动引起的吸收峰,波罗的海琥珀在此范围内存在1732cm-1处的特征峰;由于缅甸琥珀的成熟度较高,在1724cm-1处存在一强峰,抚顺琥珀在1724cm-1、1697cm-1均有吸收峰。1300-1000cm-1的吸收峰由C-O单键振动引起,波罗的海琥珀在1250-1175cm-1存在特征的肩峰,而抚顺琥珀和缅甸琥珀在此范围内由一组强吸收峰,呈“山”字形,另外抚顺琥珀1136cm-1附近的峰较弱,1259-1033cm-1之间有多个弱吸收。波罗的海琥珀成熟度较缅甸、抚顺琥珀更低,其红外光谱中存在与双键有关的1647cm-1、887cm-1吸收峰。
在经过多个变量的线性变化组合后,由表2可知前8个主成分能够解释原始波长变量的94.04%,说明前8个主成分可代表不同产地琥珀红外光谱的主要信息。其中第一主成分特征值为0.01475,对结果的总贡献率为47.26%,第二主成分特征值为0.00601,对结果的总贡献率为19.26%,第三主成分特征值为0.00415,对结果的总贡献率为13.29%,第四主成分特征值为0.00138,对结果的总贡献率为4.42%,第五主成分特征值为0.00601,对结果的总贡献率为19.26%,第六主成分特征值为9.17433E-4,对结果的总贡献率为2.16%,第七主成分特征值为6.45063E-4,对结果的总贡献率为2.07%,第八主成分特征值为5.12397E-4,对结果的总贡献率为1.64%,第一主成分到第八主成分贡献率累计已达93.03%,表明主成分分析效果较好,可以有效地提取出红外光谱原始信息。如表2可知:
表2.
其中,x1与1737cm-1处琥珀中琥珀酸酯中C=O伸缩振动及1698cm-1处羧酸中C=O振动相关;x2与琥珀酸酯中1733cm-1处C=O伸缩振动及1693cm-1处羧酸中C=O振动的相关;x3与1710cm-1处羧酸中C=O振动及1378cm-1处甲基、亚甲基中的C-H弯曲振动相关;x4与1666cm-1处不饱和C=C伸缩振动及1162cm-1处含氧官能团C-O单键伸缩振动相关;x5与1737cm-1处琥珀酸酯中C=O伸缩振动及1704cm-1处羧酸中C=O振动的相关;x6与1716cm-1处羧酸中C=O振动与1658cm-1处不饱和C=C伸缩振动相关;x7与1627cm-1处不饱和C=C伸缩振动及1253cm-1处含氧官能团C-O单键伸缩振动;x8与1741cm-1处琥珀酸酯中C=O伸缩振动及1718cm-1处羧酸中C=O振动的相关。
以1、2、3代号分别代表波罗的海、缅甸和抚顺琥珀,以产地代表符号为应变量,共30个样品的八个主成分得分(x1-x8)为自变量进行多重共线性分析,以排除多重共线性数据。如表3结果显示所选择的八个主成分得分数据均不具有共线性,可以进行判别分析。
表3.八个主成分得分和产地之间的共线性分析结果
取30个样品的八个主成分得分进行判别分析及交叉检验,采用一起输入变量法进行分析,得到出土琥珀产地判别的函数,并绘制三个产地数据的典则函数散点图。
如图2所示,为不同产地琥珀的红外光谱谱图,纵坐标代表相对吸光度,横坐标表示不同波数(cm-1)。
从不同产地琥珀的红外光谱图可以看出,三个产地之间的红外谱峰存在区别,其中部分谱峰具有产地特征,如波罗的海琥珀的888cm-1吸收峰在其他两个产地琥珀的红外光谱中没有出现。琥珀在成熟过程中环外双键逐渐消失,羰基逐渐转化为酯,与成熟度相关的888cm-1处吸收峰强度降低直至消失。
如图3所示,为缅甸和中国辽宁抚顺琥珀的光致发光光谱,纵坐标代表归一化后的发光强度,横坐标表示不同波长(nm)。
不同产地琥珀的树源不同,埋藏环境不同,物质成分不同,其发光行为之间也存在这差异。光致发光光谱可以捕捉这种差异,并以谱图形式显示,图3中还可以看出三个产地在同一激发光源下的发光位置有所区别,其中波罗的海琥珀和其他两个产地发光峰位置差别较大,缅甸琥珀和中国辽宁抚顺琥珀发光位置存在相似性,因此可以通过光致发光光谱将波罗的海琥珀与其他两个产地的琥珀进行进一步区分。
如图4所示,为缅甸和中国辽宁抚顺琥珀的紫外可见吸收光谱,纵坐标代表反射率,横坐标表示不同波长(nm)。
由于波罗的海琥珀通过红外光谱和光致发光光谱可以较为简便的同缅甸琥珀和中国辽宁抚顺琥珀进行区分,而后两种琥珀由于地质时期相近,树源相似,内部官能团组成及发光行为较为类似,红外光谱和光致发光光谱较难对两个产地琥珀进行区分,而其紫外可见吸收光谱中两者有所差别,可作为缅甸和抚顺琥珀产地区分的依据。
图4中还可以看出,缅甸琥珀和抚顺琥珀的紫外可见吸收光谱谱线走势基本相似,具有相近位置的反射峰和吸收峰,但其中吸收率较高位置的几个吸收峰位置有细微差别,可作为产地判别依据。
表4为判别分析所得出的典型判别式函数系数,根据系数建立琥珀产地溯源的判别分析函数。
表4.
所建立的红外光谱特征信息产地溯源判别函数表述为:
Y1(x)=38.977x1-41.439x2+15.038x3+8.046x4-7.024x5+16.104x6-34.596x7+25.375x8;
Y2(x)=8.952x1+19.421x2-8.037x3+8.765x4-19.969x5-7.038x6+18.129x7+34.421x8
其中x1-x8为对三个产地琥珀样品红外光谱数据进行主成分分析后贡献率最高的8个主成分得分。x1与1737cm-1处琥珀中琥珀酸酯中C=O伸缩振动及1698cm-1处羧酸中C=O振动相关;x2与琥珀酸酯中1733cm-1处C=O伸缩振动及1693cm-1处羧酸中C=O振动的相关;x3与1710cm-1处羧酸中C=O振动及1378cm-1处甲基、亚甲基中的C-H弯曲振动相关;x4与1666cm-1处不饱和C=C伸缩振动及1162cm-1处含氧官能团C-O单键伸缩振动相关;x5与1737cm-1处琥珀酸酯中C=O伸缩振动及1704cm-1处羧酸中C=O振动的相关;x6与1716cm-1处羧酸中C=O振动与1658cm-1处不饱和C=C伸缩振动相关;x7与1627cm-1处不饱和C=C伸缩振动及1253cm-1处含氧官能团C-O单键伸缩振动;x8与1741cm-1处琥珀酸酯中C=O伸缩振动及1718cm-1处羧酸中C=O振动的相关。
如图5所示,为典则判别函数的二维散点图,可以看出该函数可以较好的将三个产地的琥珀进行区分,并且没有未被分类的样品。根据函数1、2为坐标轴,组1(波罗的海)组质心为(-7.982,-0.529),组2(缅甸)组质心为(2.671,2.663),组3(抚顺)组质心为(5.311,-2.134),表5为组质心处的函数值表,
表5.
将30个样本代入所述的判别方程,按各母体的后验概率重新分类,表6为经过交叉检验后的分类结果,所有样品都得到了正确的产地判别,判别正确率达100%,说明此判别函数可以很好的对出土琥珀进行产地鉴别。
表6.
为了更好的阐述本发明的出土琥珀产地溯源模型的优势,本发明人还进行了具体实施验证:
实施例2
出土汉代琥珀产地溯源
对出土的一件汉代金珀饰物进行了检测。该饰物保存较好,对其获取官能团的特征峰,其在405nm下的发光图谱与缅甸琥珀和抚顺琥珀相近。
为了进一步识别其产地,将经实施例1中谱线处理方法对归一化后红外光谱数据按照所述产地溯源模型进行计算,将出土样品与参考样品进行主成分分析,其前八个主成分分得分分别为x1=0.17913,x2=0.06973,x3=-0.14923,x4=-0.02389,x5=-0.00165,x6=0.01279,x7=-0.00157,x8=0.02016,代入判别方程可得Y1=2.43950434,Y2=4.55616141,与模型建立所使用的标准样品共同绘制散点图图6,分别计算与三个组质心的距离,得到其属于组2(缅甸)。
图6中可以看出该出土金珀(实心点为此样品)属于缅甸组群。该出土琥珀谱学参数和主成分分析结果均显示其与缅甸琥珀一致。
实施例3
出土清代琥珀产地溯源
对出土的一件清代琥珀饰物进行了检测。由于该件样品氧化严重,且表面产生了较厚的氧化层,并蚀化掉渣,反射法测试红外光谱较为困难。取样品表面氧化程度较低的碎渣,按照1:100比例与溴化钾粉末混匀、研磨、烘干、装模,压制为透明薄片进行透射法红外光谱测试。对所得到的光谱进行K-K转换为吸光度光谱,经与三个产地标准样品的红外谱图对比,谱图与波罗的海琥珀谱峰较为相似。将经实施例1中谱线处理方法对归一化后红外光谱数据按照所述产地溯源模型进行计算,将出土样品与参考样品进行主成分分析,其前八个主成分分得分分别为x1=-0.06535,x2=0.39779,x3=0.62141,x4=0.60344,x5=0.21039,x6=0.06617,x7=0.04802,x8=0.02016,代入判别方程可得Y1=-7.58008504,Y2=2.722631,与模型建立所使用的标准样品共同绘制散点图图7(实心点为此样品),分别计算与三个组质心的距离,得到其属于组(波罗的海)。确定该样品来自波罗的海。
在不冲突的情况下,本文中上述实施例及实施例中的特征可以相互结合。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于光谱指纹的琥珀产地溯源模型的构建方法,其特征在于,包括以下具体步骤:
S1、取不同产地的琥珀样品,将样品进行预处理,所述产地包括缅甸、波罗的海和中国辽宁抚顺;
S2、将处理后的样品进行红外光谱采集,获得其官能团特征峰;
S3、采用OPUS软件对得到的红外光谱图进行谱线归一化处理,然后采用Origin软件进行曲线拟合,获得峰位;
S4、确定峰位代表的官能团,运用主成分分析和判别分析对所得官能团数据进行分析,建立产地溯源模型;
所述产地溯源模型包括:根据每个样品的红外光谱,得到官能团数据,运用主成分分析初步判别影响较大的官能团数据;利用线性回归分析,判别官能团的多重共线性,其中VIF小于10视为不存在多重共线性;根据筛选后的不具有多重共线性的官能团数据进行判别分析;利用判别分析和交叉检验准确判别的不同产地琥珀的分布规律,
红外光谱特征信息指示的琥珀产地溯源的判别函数分别为:
Y1(x)=38.977x1-41.439x2+15.038x3+8.046x4-7.024x5+16.104x6-34.596x7+25.375x8;
Y2(x)=8.952x1+19.421x2-8.037x3+8.765x4-19.969x5-7.038x6+18.129x7+34.421x8;
其中,x1、x2、x3、x4、x5、x6、x7和x8代表八个主成分因子特征值,具体为总贡献率的得分值,x1与1737cm-1处琥珀中琥珀酸酯中C=O伸缩振动及1698cm-1处羧酸中C=O振动相关;x2与琥珀酸酯中1733cm-1处C=O伸缩振动及1693cm-1处羧酸中C=O振动的相关;x3与1710cm-1处羧酸中C=O振动及1378cm-1处甲基、亚甲基中的C-H弯曲振动相关;x4与1666cm-1处不饱和C=C伸缩振动及1162cm-1处含氧官能团C-O单键伸缩振动相关;x5与1737cm-1处琥珀酸酯中C=O伸缩振动及1704cm-1处羧酸中C=O振动的相关;x6与1716cm-1处羧酸中C=O振动与1658cm-1处不饱和C=C伸缩振动相关;x7与1627cm-1处不饱和C=C伸缩振动及1253cm-1处含氧官能团C-O单键伸缩振动;x8与1741cm-1处琥珀酸酯中C=O伸缩振动及1718cm-1处羧酸中C=O振动的相关。
2.如权利要求1所述的构建方法,其特征在于,步骤S2中,采集的红外光谱数据范围为1800-400cm-1。
3.如权利要求1所述的构建方法,其特征在于,步骤S1中,所述样品预处理步骤包括:
S11、对样品表面进行清理;
S12、进行样品进行切割,制成厚度为2~5mm的薄片;
S13、对薄片样品进行双面平行抛光。
4.如权利要求1所述的构建方法,其特征在于,所述红外光谱采集条件为:背景扫描32次,样品扫描64次,分辨率为2cm-1,测试范围为400-4000cm-1。
5.如权利要求4所述的构建方法,其特征在于,步骤S3中,所述红外光谱谱线处理方法为:采用红外处理软件OPUS中气氛补偿进行H2O和CO2补偿,采用凹式rubbenband校正方法进行基线校正,迭代次数为25,基线点数为100,校正范围选择为4000-400cm-1波数区域进行基线校正,使用OPUS软件的向量归一化方法进行谱线归一化处理。
6.一种出土琥珀产地的鉴别方法,其特征在于,包括以下具体步骤:
S1、获得待测出土琥珀样品在波长1800~400cm-1的红外光谱图;
S2、按照权利要求1-5中任一项所述的构建方法,将步骤S1获得的红外光谱图带入所述构建方法中运算进行主成分分析,得到八个主成分x1、x2、x3、x4、x5、x6、x7和x8的得分值;
S3、将x1~x8的得分值代入权利要求1~5中任一项所述的构建方法得到的产地溯源模型中,得到判别界值;
S4、对比判别界值到组质心距离最小即可判别产地;
若当同时判别为波罗的海和缅甸,需采用光致发光光谱鉴别,其中在405nm激发下,波罗的海琥珀荧光特征峰位于500-550nm范围内,缅甸琥珀荧光特征峰位于450-500nm范围内;
若当同时判别为波罗的海和中国辽宁抚顺,需采用光致发光光谱鉴别,其中在405nm激发下,波罗的海琥珀荧光特征峰位于500-550nm范围内,中国辽宁抚顺琥珀荧光特征峰位于450-500nm范围内;
若当同时判别为缅甸和中国辽宁抚顺,需采用紫外可见吸收光谱鉴别,比较875nm-925nm范围内的是否有吸收峰,有吸收峰的为中国辽宁抚顺琥珀,没有吸收峰的为缅甸琥珀。
7.如权利要求6所述的鉴别方法,其特征在于,所述光致发光光谱采集条件为:激发光源波长405nm,测试范围210-1000nm,扫描时间5s,扫描次数3。
8.如权利要求6所述的鉴别方法,其特征在于,所述紫外可见吸收光谱采集条件为:测试方法为透射法,测试范围为200-800nm,数据间隔1nm,探测器狭缝宽度2nm,扫描速度为266.75nm/min。
9.一种电子设备,包括处理器和存储有计算机程序的存储器,其特征在于,所述计算机程序被所述处理器执行时,使得所述处理器执行如权利要求1~5中任一项所述的构建方法。
10.一种计算机可读存储介质,存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,使得所述处理器执行如权利要求1~5中任一项所述的构建方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211429390.5A CN115718081B (zh) | 2022-11-15 | 2022-11-15 | 基于光谱指纹的琥珀产地溯源模型的构建方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211429390.5A CN115718081B (zh) | 2022-11-15 | 2022-11-15 | 基于光谱指纹的琥珀产地溯源模型的构建方法及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115718081A CN115718081A (zh) | 2023-02-28 |
CN115718081B true CN115718081B (zh) | 2024-04-12 |
Family
ID=85255161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211429390.5A Active CN115718081B (zh) | 2022-11-15 | 2022-11-15 | 基于光谱指纹的琥珀产地溯源模型的构建方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115718081B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114965439B (zh) * | 2022-05-07 | 2024-09-06 | 中国地质大学(武汉) | 基于人工神经网络的青白色软玉产地溯源研究方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852422A (zh) * | 2014-02-14 | 2014-06-11 | 温州市质量技术监督检测院 | 一种无损鉴定琥珀类样品的方法 |
JP2019032231A (ja) * | 2017-08-08 | 2019-02-28 | 国立大学法人秋田大学 | 琥珀砕片の種別および透明度の推定方法、推定装置、ならびに推定プログラム |
JP6679053B1 (ja) * | 2019-01-23 | 2020-04-15 | 九鬼産業株式会社 | ゴマの産地判別装置及びゴマの産地判別方法 |
CN112014346A (zh) * | 2020-09-03 | 2020-12-01 | 中国地质大学(武汉) | 一种基于红外光谱的煤炭产地溯源方法 |
CN113552074A (zh) * | 2021-07-14 | 2021-10-26 | 浙江方圆检测集团股份有限公司 | 基于UV-Vis漫反射光谱与PL光谱特征判定蓝珀产地溯源的方法 |
-
2022
- 2022-11-15 CN CN202211429390.5A patent/CN115718081B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852422A (zh) * | 2014-02-14 | 2014-06-11 | 温州市质量技术监督检测院 | 一种无损鉴定琥珀类样品的方法 |
JP2019032231A (ja) * | 2017-08-08 | 2019-02-28 | 国立大学法人秋田大学 | 琥珀砕片の種別および透明度の推定方法、推定装置、ならびに推定プログラム |
JP6679053B1 (ja) * | 2019-01-23 | 2020-04-15 | 九鬼産業株式会社 | ゴマの産地判別装置及びゴマの産地判別方法 |
CN112014346A (zh) * | 2020-09-03 | 2020-12-01 | 中国地质大学(武汉) | 一种基于红外光谱的煤炭产地溯源方法 |
CN113552074A (zh) * | 2021-07-14 | 2021-10-26 | 浙江方圆检测集团股份有限公司 | 基于UV-Vis漫反射光谱与PL光谱特征判定蓝珀产地溯源的方法 |
Non-Patent Citations (4)
Title |
---|
Effects of weathering on FTIR spectra and origin traceability of archaeological amber: The case of the Han Tomb of Haihun Marquis, China;Tong Zhao 等;Journal of Archaeological Science;20230317;全文 * |
出土琥珀产地指纹特征筛选与溯源方法研究;李 妍 等;宝石和宝石学杂志 ( 中英文 );20230731;第25卷(第4期);全文 * |
琥珀的激光拉曼光谱特征研究;吴文杰;王雅玫;;宝石和宝石学杂志;20140215(01);全文 * |
辽代的资源获取与文化互动-以陈国公主及驸马墓葬为例;赵彤 等;宝石和宝石学杂志(中英文);20221130;第24卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115718081A (zh) | 2023-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Discrimination of Ganoderma lucidum according to geographical origin with near infrared diffuse reflectance spectroscopy and pattern recognition techniques | |
Fortier et al. | Identification of cotton and cotton trash components by Fourier transform near-infrared spectroscopy | |
CN109540831B (zh) | 基于高光谱成像技术的枸杞子品种识别方法 | |
CN104251839B (zh) | 南疆红枣建模用南疆红枣样品成分的光谱分离检测方法 | |
CN115718081B (zh) | 基于光谱指纹的琥珀产地溯源模型的构建方法及应用 | |
CN103278473B (zh) | 白胡椒中胡椒碱及水分含量的测定和品质评价方法 | |
CN110376153B (zh) | 一种atr-ftir结合rbf神经网络对市售西红花产地溯源的方法 | |
CN106501208A (zh) | 一种基于近红外光谱特征的烟叶风格相似性分类方法 | |
WO2023077859A1 (zh) | 一种基于空间偏移拉曼光谱的带包装食品信号分离方法 | |
Li et al. | Multi-information based on ATR-FTIR and FT-NIR for identification and evaluation for different parts and harvest time of Dendrobium officinale with chemometrics | |
CN102937575B (zh) | 一种基于二次光谱重组的西瓜糖度快速建模方法 | |
CN113008817A (zh) | 一种基于高光谱成像技术快速鉴别苦杏仁真伪优劣的方法 | |
Sun et al. | Non-destructive detection for foreign bodies of tea stalks in finished tea products using terahertz spectroscopy and imaging | |
Cui et al. | Identification of maize seed varieties based on near infrared reflectance spectroscopy and chemometrics | |
Karolina et al. | Identification of treated Baltic amber by FTIR and FT-Raman–a feasibility study | |
Chen et al. | Rapid identification of three varieties of Chrysanthemum with near infrared spectroscopy | |
Zhuang et al. | Rapid determination of green tea origins by near-infrared spectroscopy and multi-wavelength statistical discriminant analysis | |
CN107402192A (zh) | 一种快速分析香精香料质量稳定性的方法 | |
CN104596975A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶木质素的方法 | |
CN104596979A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶纤维素的方法 | |
CN105223140A (zh) | 同源物质的快速识别方法 | |
CN109358022A (zh) | 一种快速判别烟用爆珠类型的方法 | |
CN113030007B (zh) | 基于相似度学习算法快速检验烟用香精质量稳定性的方法 | |
CN111398208B (zh) | 一种利用近红外技术对传统工艺翡翠进行快速鉴定的方法 | |
CN115901694A (zh) | 一种鉴别琥珀产地的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |