CN110487987B - 流域农业面源污染的frn-cssi联合溯源示踪解析方法 - Google Patents
流域农业面源污染的frn-cssi联合溯源示踪解析方法 Download PDFInfo
- Publication number
- CN110487987B CN110487987B CN201910728202.0A CN201910728202A CN110487987B CN 110487987 B CN110487987 B CN 110487987B CN 201910728202 A CN201910728202 A CN 201910728202A CN 110487987 B CN110487987 B CN 110487987B
- Authority
- CN
- China
- Prior art keywords
- frn
- silt
- soil
- types
- source
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Processing Of Solid Wastes (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明提供一种流域农业面源污染的FRN‑CSSI联合溯源示踪解析方法,包括以下步骤:(1)FRN示踪和CSSI溯源技术的样品采集;(2)FRN示踪和CSSI溯源技术的样品分析;(3)流域不同土地利用类型/植被类型的输沙速率;(4)流域不同植被类型泥沙来源贡献比例;(5)泥沙及其污染物来源解析。本发明采用FRN和CSSI联合溯源技术,可以定量流域泥沙污染物的土地利用来源及其输入水系的负荷通量法,根据不同土地利用类型的土壤侵蚀产沙速率和土地利用面积,可以定量辨识和划分流域内土壤侵蚀及其污染的关键控制区域。
Description
技术领域
本发明属于污染监控技术领域,具体涉及一种流域农业面源污染的FRN-CSSI联合溯源示踪解析方法。
背景技术
降雨径流引起农业源物质(泥沙、氮磷、农药等)向地表水系输送,已经成为全球河流、湖泊、水库中面源污染的最重要贡献,也正威胁着流域水体质量。为了有效防控输入河湖库的泥沙污染负荷,必须首先辨识泥沙污染的来源,弄清污染物进入水体的路径。传统面源污染来源解析方法,如地球化学指纹技术、生物指纹、SWAT等模型解析及同位素示踪等均不能辨识不同土地利用/植被类型的污染物来源。联合FRN示踪和CSSI溯源技术,可以定量解析流域泥沙污染的土地利用来源及其输入水系的负荷通量,相比传统的农业面源污染来源解析方法具有显著优势,可有效防控输入河湖库的泥沙污染负荷。
散落的环境放射性核素(FRN,fallout radionuclides)示踪技术是利用环境放射性核素定量评价不同土地利用类型土壤侵蚀产沙速率的方法(Ritchie and Mc Henry,1990;Walling et al.,1993;Li et al.,2015)。全球普遍应用的环境放射性核素包括137Cs、210Pb和7Be等,分别来源于上世纪50-70年代全球大气核爆试验、土壤/岩石中的释放和宇宙射线。这些核素均随降雨散落在地表,很快被地表土壤细颗粒吸附固定,只有在地表径流、风力等外力作用下发生运移。FRN在全球广泛分布,一般在北半球的含量高于南半球。由于137Cs、210Pb和7Be半衰期不同,分别为30年、22年和53天,因此采集一次FRN土壤样品,可以确定不同时间(50年、100年和1个月/次降雨事件)和空间尺度土壤流失和泥沙堆积速率(Zapata,2003)。目前FRN示踪技术应用主要在两个方面:一是确定不同土地利用类型或土地利用管理变化引起的土壤流失再分布速率,二是用于辨析河库泥沙污染物沉积的时间和泥沙的侵蚀类型来源。
特定单体化合物稳定性同位素(CSSI,compound-specific stable isotope)分析技术是根据不同土地利用类型中的植物特有的土壤有机标志物(如脂肪酸和正构烷烃),辨识不同土地利用类型对河库泥沙贡献比例的生物指纹技术(Gibbs,2008;Mabit,et al.,2013;Reiffarth,et al.,2014;Tolosa et al.,2013)。与传统示踪方法比较,CSSI技术在辨识湖泊、水库的泥沙来源方面具有专一、快速、动态等优点。Gibbs(2008)以新西兰Mahurangi河流域为对象,调查分析了流域主要土地利用类型的特定单体化合物稳定性碳同位素。研究结果表明,河口三角洲沉积的泥沙主要来源于流域内的松树林土壤,认为砍伐松树林导致土壤侵蚀产沙是河口泥沙沉积加剧的主要动因。Blake等(2012)利用CSSI技术研究了英国西南部River Otter支流泥沙来源,发现占流域面积的65%的草地对流域泥沙的总负荷贡献最大,但草地单位面积的侵蚀产沙量(0.13±0.02t ha-1)却显著小于冬小麦地(0.44±0.15t ha-1)。Hancock等(2013)应用CSSI方法评价了澳大利亚昆士兰州Logan和Albert Rivers流域的泥沙污染来源,指出应用脂肪酸中的δ13C值可以辨别流域侵蚀产沙的来源。Alewell等(2016)应用CSSI技术研究了瑞士中部低地的河流中泥沙的植被来源类型,发现林地和农地的侵蚀产沙对河流水体的贡献有明显差异。
然而,FRN示踪技术可以确定不同土地利用类型土壤侵蚀速率的时空变异,但不能按照土地利用类型区分流域水系中泥沙污染的来源。相反,CSSI技术可以用于辨识地表水体中泥沙污染物土地利用来源的贡献比例,但无法了解不同土地利用类型土壤的侵蚀产沙速率。联合FRN示踪和CSSI溯源技术,可以定量了解流域泥沙污染的土地利用来源及其输入水系的负荷通量。目前国内外学者通常分开应用FRN示踪和CSSI溯源,难于定量解析河湖水库流域泥沙污染物的来源类型,不能解释流域侵蚀产生的泥沙及其携带的污染物输送过程与水体污染负荷的源汇关系。
发明内容
本发明提供一种流域农业面源污染的FRN-CSSI联合溯源示踪解析方法,可以了解上游泥沙污染来源区与下游水体污染汇集区的定量关系,为流域水系面源污染的防控提供科技支撑。FRN-CSSI联合溯源示踪技术能解决如下问题:
(1)从流域尺度了解上游不同土地利用或作物类型土壤向下游水体输送的污染通量:
(2)能够用于评价流域环境治理措施的有效性。
流域农业面源污染的FRN-CSSI联合溯源示踪解析方法,包括以下步骤:
(1)FRN示踪和CSSI溯源技术的样品采集
在流域景观地,选择未扰动地点作为参考点,采集FRNs(如137Cs、210Pb和7Be)土壤剖面样品,采样深度超过核素分布深度,以确保采集到土壤剖面所有核素;在各子流域不同土地利用类型采集土壤样品,采集深度超过核素分布深度;
流域内CSSI样品采集分为两类;一类是采集代表泥沙来源的子流域不同土地利用类型土壤;一类是采集流域出口低洼处的沉积泥沙样品;泥沙来源和混合样的采样方法参考Gibbs(2014)的操作步骤,所有样品采集深度为表层0-2cm,每个采样点采集10个样品混合为一个样品;
(2)FRN示踪和CSSI溯源技术的样品分析
采集的FRN样品经风干过2mm筛后进行核素活度测定;参考点和不同土地利用类型土壤核素测定采用无源效率刻度伽玛能谱仪(BE5030宽能型HPGe探测器);同时需要测定采集的土壤样品容重;
采集的不同土地利用类型土壤和流域泥沙样品一部分经二氯甲烷提取、皂化后,调至酸性过无水硫酸钠分离柱,然后采用三氟化硼甲醇溶液甲酯化后溶解在正己烷溶剂,利用气相色谱-同位素质谱联用(GC-IRMS)测定脂肪酸中的δ13C丰度;此外,一部分土壤和流域泥沙样品经酸化后,直接用元素分析仪测定土壤bulkδ13C丰度;
(3)流域不同土地利用类型的输沙速率
根据参考点和各研究点的FRN活度、容重和采样深度,计算其面积活度(公式2);
式中:A-面积活度,Bq/m2;Cn-第n层质量活度,Bq/kg;ρn-第n层容重,g/cm3;hn-第n层深度,cm;
不同土地利用类型需要选择不同的侵蚀模型(Zapata,2002):农地可以采用比例模型或质量平衡模型,灌木地、林地和草地选用轮廓形状模型或扩散迁移模型;根据流域内各土地利用类型选择相应的侵蚀模型,输入参考点和各土地利用类型的FRN面积活度,估算不同土地利用类型泥沙污染物输出速率(tha-1yr-1);
(4)流域不同土地利用类型泥沙来源贡献比例
通过方差分析选择合适的脂肪酸来确定泥沙中不同土地利用类型来源的碳同位素比例;基于不同土地利用类型筛选的脂肪酸δ13C和bulkδ13C的结果,应用Iso-Source混合模型(Gibbs,2014)确定泥沙混合物中不同土地利用类型来源的碳同位素比例;将泥沙中的有机碳含量和模型输出碳同位素比例转换成不同土地利用类型泥沙来源贡献(公式2);
式中:Sn%-泥沙第n个源土壤贡献,In-使用同位素混合模型估算的混合物中第n个源土壤的同位素比例(%),Cn%-第n个源土壤有机碳含量;
(5)泥沙及其污染物来源解析
基于FRN估算的不同土地利用类型土壤泥沙污染物输出速率和CSSI解析的入河泥沙污染的土地利用类型贡献比例,定量计算流域泥沙污染的土地利用来源及其输入水系的泥沙通量(公式3);
Fn=Rn×Sn%×Sn(3)
式中:Fn-不同土地利用类型泥沙通量(t yr-1);Rn-泥沙第n个源土地利用类型泥沙污染物输出速率(t ha-1yr-1);Sn%-泥沙第n个源土壤贡献;Sn-泥沙第n个源土地利用面积。
本发明提供的流域农业面源污染的FRN-CSSI联合溯源示踪解析方法,可以定量了解流域泥沙污染的土地利用来源及其输入水系的泥沙通量;根据不同土地利用类型的土壤侵蚀产沙速率和土地利用面积,可以定量辨识和划分流域内土壤侵蚀及其污染的关键控制区域。相比传统的农业面源污染来源解析方法,本发明具有显著优势,可有效防控输入河湖库的泥沙污染负荷。
附图说明
图1为实施例的流域泥沙及其污染物来源解析结果。
具体实施方式
结合实施例说明本发明的技术方案。
以某流域为例,流域出口泥沙来自3个子流域,每个流域土地利用类型有差异,在每个子流域不同土地利用类型采集FRN和CSSI样品并分析,基于FRN估算的不同土地利用类型土壤泥沙污染物输出速率、CSSI解析的入河泥沙污染的土地利用类型贡献比例,可以定量计算流域泥沙污染的土地利用来源及其输入水系的泥沙通量,如图1。
流域农业面源污染的FRN-CSSI联合溯源示踪解析方法的应用会带来如下有益效果:
(1)有助于从流域尺度定量上游不同土地利用类型或作物类型土壤向下游水体输送的农业源污染负荷,有助于提升农业面源污染防控技水平。联合FRN示踪和CSSI溯源技术,可以定量了解流域泥沙污染的土地利用来源及其输入水系的泥沙通量。根据不同土地利用类型的土壤侵蚀产沙速率和土地利用面积,可以定量辨识和划分流域内土壤侵蚀及其污染的关键控制区域。以上述流域为例,联合FRN示踪和CSSI溯源技术,可以定量计算流域泥沙污染的土地利用来源及其输入水系的泥沙通量(表1)。结果显示农地是该流域下游水体污染物的主要来源。通过在农地采取保护性耕作、降低农药化肥施用量等方法减少流域面源污染。
表1FRN-CSSI联合示踪的不同流域输沙通量
(2)可定量评价环境治理措施的有效性。在农业面源污染严重的流域,经过环境治理后,需要定量评价治理效果。联合FRN示踪和CSSI溯源技术,确定不同土地利用类型土壤侵蚀速率,辨识地表水体中泥沙污染物土地利用来源的贡献比例,但定量计算流域泥沙污染的土地利用来源及其输入水系的泥沙通量。通过治理前后流域农业面源污染的来源贡献调查,比较环境治理措施实施前后氮磷等污染物来源的泥沙通量变化,评价不同治理措施的有效性。
Claims (2)
1.流域农业面源污染的FRN-CSSI联合溯源示踪解析方法,其特征在于,包括以下步骤:
(1)环境放射性核素示踪和单体化合物稳定性同位素溯源技术的样品采集
在流域景观地,选择未扰动地点作为参考点,采集FRN土壤剖面样品,采样深度超过核素分布深度,以确保采集到土壤剖面所有核素;在各子流域不同土地利用类型采集土壤样品,采集深度超过核素分布深度;
流域内CSSI样品采集分为两类:一类是采集代表泥沙来源的子流域不同土地利用类型土壤;一类是采集流域出口低洼处的沉积泥沙样品;所有样品采集深度为表层0-2cm,每个采样点采集10个样品混合为一个样品;
(2)FRN示踪和CSSI溯源技术的样品分析
采集的FRN样品经风干过2mm筛后进行核素活度测定;参考点和不同土地利用类型土壤核素测定采用无源效率刻度伽玛能谱仪;同时需要测定采集的土壤样品容重;
采集的不同土地利用类型土壤和流域泥沙样品一部分经二氯甲烷提取、皂化后,调至酸性过无水硫酸钠分离柱,然后采用三氟化硼甲醇溶液甲酯化后溶解在正己烷中,利用气相色谱-同位素质谱联用测定脂肪酸中的δ13C丰度;另一部分样品经酸化后,用元素分析仪测定土壤bulkδ13C丰度;
(3)流域不同土地利用类型的输沙速率
根据参考点和各研究点的FRN活度、容重和采样深度,使用公式1计算其面积活度:
式中:A-面积活度,Bq/m2;Cn-第n层质量活度,Bq/kg;ρn-第n层容重,g/cm3;hn-第n层深度,cm;
不同土地利用类型需要选择不同的侵蚀模型;根据流域内各土地利用类型选择相应的侵蚀模型,输入参考点和各土地利用类型的FRN面积活度,估算不同土地利用类型泥沙污染物输出速率,t ha-1yr-1;
(4)流域不同土地利用类型泥沙来源贡献比例
通过方差分析选择合适的脂肪酸来确定泥沙中不同土地利用类型来源的碳同位素比例;基于不同土地利用类型筛选的脂肪酸δ13C和bulkδ13C的结果,应用Iso-Source混合模型确定泥沙混合物中不同土地利用类型来源的碳同位素比例;根据公式2将泥沙中的有机碳含量和模型输出碳同位素比例转换成不同土地利用类型泥沙来源贡献:
式中:Sn%-泥沙第n个源土壤贡献;In-使用同位素混合模型估算的混合物中第n个源土壤的同位素比例,%;Cn%-第n个源土壤有机碳含量;
(5)泥沙及其污染物来源解析
基于FRN估算的不同土地利用类型土壤泥沙污染物输出速率和CSSI解析的入河泥沙污染的土地利用类型贡献比例,根据公式3定量计算流域泥沙污染的土地利用来源及其输入水系的泥沙通量:
Fn=Rn×Sn%×Sn (3)
式中:Fn-不同土地利用类型泥沙通量,t yr-1;Rn-泥沙第n个源土地利用类型泥沙污染物输出速率,t ha-1yr-1;Sn%-泥沙第n个源土壤贡献;Sn-泥沙第n个源土地利用面积。
2.根据权利要求1所述的流域农业面源污染的FRN-CSSI联合溯源示踪解析方法,其特征在于,步骤(3)中,所述的不同土地利用类型需要选择不同的侵蚀模型:农地采用比例模型或质量平衡模型;灌木地、林地和草地选用轮廓形状模型或扩散迁移模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910728202.0A CN110487987B (zh) | 2019-08-08 | 2019-08-08 | 流域农业面源污染的frn-cssi联合溯源示踪解析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910728202.0A CN110487987B (zh) | 2019-08-08 | 2019-08-08 | 流域农业面源污染的frn-cssi联合溯源示踪解析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110487987A CN110487987A (zh) | 2019-11-22 |
CN110487987B true CN110487987B (zh) | 2020-08-28 |
Family
ID=68550215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910728202.0A Expired - Fee Related CN110487987B (zh) | 2019-08-08 | 2019-08-08 | 流域农业面源污染的frn-cssi联合溯源示踪解析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110487987B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113536598B (zh) * | 2021-08-13 | 2022-04-26 | 谢军 | 一种基于CAESAR-Lisflood模型的流域泥沙溯源模拟方法 |
CN114354830B (zh) * | 2021-12-27 | 2023-10-10 | 中国环境科学研究院 | 有机氯污染物溯源方法及装置、计算机设备和存储介质 |
CN114460270B (zh) * | 2022-01-19 | 2024-04-19 | 广西大学 | Cssi联合实时监测定量解析侵蚀泥沙重金属入河负荷的方法 |
CN114384224B (zh) * | 2022-01-19 | 2022-08-05 | 中国农业科学院农业环境与可持续发展研究所 | 基于多同位素联合示踪的流域氮磷污染物解析方法与*** |
CN114460271A (zh) * | 2022-01-19 | 2022-05-10 | 广西大学 | 基于cssi溯源与实时监测的流域入河泥沙来源负荷辨识方法 |
CN114580167A (zh) * | 2022-03-02 | 2022-06-03 | 广西大学 | 一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2159198A1 (en) * | 2008-09-01 | 2010-03-03 | Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH) | Method for the degradation of pollutants in water and/ or soil |
CN104933300B (zh) * | 2015-06-03 | 2018-02-06 | 中国农业科学院农业资源与农业区划研究所 | 流域农业面源污染物河道削减系数计算方法 |
CN107655961B (zh) * | 2017-09-26 | 2020-07-21 | 临沂大学 | 基于沉积物同位素分析的计算农业面源重金属流失负荷值的方法 |
KR101991634B1 (ko) * | 2018-10-26 | 2019-06-20 | 건국대학교 산학협력단 | 성분별 안정동위원소 분석을 이용한 유기농 우유 및 일반 우유의 판별 방법 |
-
2019
- 2019-08-08 CN CN201910728202.0A patent/CN110487987B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110487987A (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110487987B (zh) | 流域农业面源污染的frn-cssi联合溯源示踪解析方法 | |
Qiao et al. | Quantitative analysis of the factors influencing spatial distribution of soil heavy metals based on geographical detector | |
Guzmán et al. | Sediment tracers in water erosion studies: current approaches and challenges | |
Schuller et al. | Using 137Cs and 210Pbex and other sediment source fingerprints to document suspended sediment sources in small forested catchments in south-central Chile | |
Fox et al. | The use of carbon and nitrogen isotopes to study watershed erosion processes 1 | |
Mukundan et al. | Sediment fingerprinting to determine the source of suspended sediment in a southern Piedmont stream | |
Iwagami et al. | Temporal changes in dissolved 137Cs concentrations in groundwater and stream water in Fukushima after the Fukushima Dai-ichi Nuclear Power Plant accident | |
Tsuji et al. | Distribution of dissolved and particulate radiocesium concentrations along rivers and the relations between radiocesium concentration and deposition after the nuclear power plant accident in Fukushima | |
Wallbrink et al. | Quantifying the contributions of sediment, sediment-P and fertiliser-P from forested, cultivated and pasture areas at the landuse and catchment scale using fallout radionuclides and geochemistry | |
Mizugaki et al. | Estimation of suspended sediment sources using 137Cs and 210Pbex in unmanaged Japanese cypress plantation watersheds in southern Japan | |
Wallbrink | Quantifying the erosion processes and land-uses which dominate fine sediment supply to Moreton Bay, Southeast Queensland, Australia | |
Palazón et al. | Combining catchment modelling and sediment fingerprinting to assess sediment dynamics in a Spanish Pyrenean river system | |
Wakiyama et al. | Land use types control solid wash-off rate and entrainment coefficient of Fukushima-derived 137Cs, and their time dependence | |
Huon et al. | Suspended sediment source and propagation during monsoon events across nested sub-catchments with contrasted land uses in Laos | |
Xiong et al. | Soil erosion and chemical weathering in a region with typical karst topography | |
Loughran et al. | A survey of soil erosion in Australia using caesium‐137 | |
Palazón et al. | Evaluating the importance of surface soil contributions to reservoir sediment in alpine environments: a combined modelling and fingerprinting approach in the Posets-Maladeta Natural Park | |
Chen et al. | Quantifying sediment source contributions in an agricultural catchment with ephemeral and classic gullies using 137Cs technique | |
Schimmack et al. | Can 239+ 240Pu replace 137Cs as an erosion tracer in agricultural landscapes contaminated with Chernobyl fallout? | |
Huang et al. | Using reservoir deposits to quantify the source contributions to the sediment yield in the Black Soil Region, Northeast China, based on the fingerprinting technique | |
Tian et al. | Assessing sediment yield and sources using fingerprinting method in a representative catchment of the Loess Plateau, China | |
CN114384224B (zh) | 基于多同位素联合示踪的流域氮磷污染物解析方法与*** | |
Wilken et al. | Understanding the role of water and tillage erosion from 239+ 240 Pu tracer measurements using inverse modelling | |
Liu et al. | Evaluation of soil erosion rates in the hilly-gully region of the Loess Plateau in China in the past 60 years using global fallout plutonium | |
Smith et al. | Modelling particle residence times in agricultural river basins using a sediment budget model and fallout radionuclide tracers |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Yu Hanqing Inventor after: Liu Wenxiang Inventor before: Yu Hanqing |
|
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200828 Termination date: 20210808 |