CN110057990A - A kind of pH bearing calibration of multi-parameter water quality section plotter - Google Patents
A kind of pH bearing calibration of multi-parameter water quality section plotter Download PDFInfo
- Publication number
- CN110057990A CN110057990A CN201910177098.0A CN201910177098A CN110057990A CN 110057990 A CN110057990 A CN 110057990A CN 201910177098 A CN201910177098 A CN 201910177098A CN 110057990 A CN110057990 A CN 110057990A
- Authority
- CN
- China
- Prior art keywords
- data
- water quality
- section
- ratio
- parameter water
- 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.)
- Granted
Links
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/18—Water
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The present invention relates to a kind of pH bearing calibrations of multi-parameter water quality section plotter, including the following steps: data acquisition: the pH section S of the observation each website in sea area is obtained using multi-parameter water quality section plotter;Data quality control: large amplitude singular value is removed by rejecting the temperature-sensitive data in instrument continuous mode, inverse pressure correction data, calculating Local standard deviation, quality control is carried out to S, Quality Control results are indicated with Q;Data are average: carrying out sliding average to erect-position downlink data each in Q, processing result is indicated with A;Data correction: on the basis of the seawater pH value for the desk-top pH meter measurement calibrated by pH standard buffer solution, processing is corrected to A.
Description
Technical field
The present invention relates to Marine Sciences fields, and in particular to a kind of pH correction of multi-parameter water quality section plotter (abbreviation RBR)
Method.
Background technique
PH indicates the negative logarithm (pH=log10 [H+]) of water middle acid substance hydrogen ion (H+) concentration, is used as water environment
The important indicator of health, is of great significance to the existence and procreation of marine organisms.About the measurement of seawater pH value, conventional chemical
Fado acquires water sample using hydrophore, and uses the pH value of desk-top pH meter measurement seawater.Due to desk-top pH meter probe price compared with
It is low, be easily changed, therefore measurement result is more accurate, but this method takes time and effort and can only obtain single-point pH data, has larger
Limitation.The appearance of multi-parameter water quality section plotter is acquisition high frequency and continuous cross-sectional data provides possibility, has become at present
For one of most widely used sea investigation instrument equipment in the world.The instrument is (deep comprising temperature, conductance (salinity), pressure
Degree), dissolved oxygen, effective photosynthesis irradiation (PAR), chlorophyll, turbidity, multiple probes such as pH, can at most measure 13 simultaneously
A underwater parameter.Due to its compact, be very suitable for coastal waters water body fast investigation, large area investigation, water quality profiling observation,
Boat observation, towed body installation observation, buoy long-term observation, seabed long term monitoring etc. are walked, but instrument works long hours or positioned at thermohaline
When being measured at spring layer, it is understood that there may be data wander or response lag phenomenon.It therefore, can be by two methods during field investigation
It combines, on the basis of the seawater pH value of desk-top pH meter measurement, school is carried out to the pH section of multi-parameter water quality section plotter observation
Just, the two mutual supplement with each other's advantages can obtain higher frequency and accurate pH section.
PH about multi-parameter water quality section plotter is corrected, and forefathers mostly use " linear regression method ", and considers multi-parameter water
The substantial measurement errors of matter section plotter pH probe, high pH is higher, and low pH is lower, i.e., slope changes in continuous mode, therefore should
The relative error of method is larger, therefore the present invention carries out school using the pH data that " ratio method " observes multi-parameter water quality section plotter
Just, it can obtain that precision is higher, pH section of better quality, provide good data for research seawater physicochemical property and support.
Summary of the invention
In view of this, improving the present invention provides a kind of pH bearing calibration of multi-parameter water quality section plotter because of instrument response
It drifts about at any time the evaluated error of generation, obtains the pH section of higher precision, provide higher frequency to study the physicochemical property of seawater
And accurate data basis.Technical scheme is as follows:
A kind of pH bearing calibration of multi-parameter water quality section plotter, including the following steps:
1) the pH section S of the observation each website in sea area data acquisition: is obtained using multi-parameter water quality section plotter;
2) data quality control: by rejecting the temperature-sensitive data in instrument continuous mode, inverse pressure correction data, calculating part
Standard deviation removal large amplitude singular value and etc., quality control is carried out to S, Quality Control results are indicated with Q;
3) data are average: carrying out sliding average to erect-position downlink data each in Q, processing result is indicated with A;
4) Data correction: on the basis of the seawater pH value for the desk-top pH meter measurement calibrated by pH standard buffer solution,
Processing is corrected to A, the method is as follows:
Step 1: removal spring layer data utilize the temperature of multi-parameter water quality section plotter simultaneous observation, salinity, the velocity of sound, the depth of water
Data calculate the buoyancy frequency N of each erect-position, and according to logN2> -3 determine the position of epilimnion in seawater, remove spring layer in A
The pH data at place, processing result indicate with T, the calculation formula of buoyancy frequency are as follows:
In formula, c0For the velocity of sound, ρ is density of sea water, and z is the depth of water, and g is acceleration of gravity, g2/co 2One magnitude is very small, suddenly
Slightly disregard;
Step 2: the seawater pH value of desk-top pH meter measurement and the ratio of pH in corresponding T are calculated, and the ratio is defined as
Ratio has:
In formula, pHbFor the seawater pH value of desk-top pH meter measurement, pHrRefer to the seawater of corresponding multi-parameter water quality section plotter measurement
PH value;
Step 3: it is assumed that Ratio approximation Normal Distribution, calculates the probability density function of Ratio, and to discrete larger
Ratio carry out 1~3 times of standard deviation and reject, processing result is indicated with R;
Step 4: data are observed with pH in corresponding T to R and carry out linear fit, the school of each erect-position pH section in calculating observation sea area
Positive coefficient, the method is as follows: enable y=ax+b, parameter a and b are used to adjust the slope and intercept of calibration curve, and x indicates multi-parameter water
The pH section A of matter section plotter observation, i.e., comprising whole pH data including spring layer, the corresponding correction ratio of the different pH of y expression;
Step 5: A is multiplied with y one-to-one correspondence, obtains complete pH correction section.
The pH bearing calibration of multi-parameter water quality section plotter of the invention, including data acquisition, data quality control, sliding are flat
Processing and Data correction.The pH section of the observation each website in sea area is obtained with multi-parameter water quality section plotter first;It connects down
To carry out data quality control, including removal temperature-sensitive data, inverse pressure correction data and deglitch three steps of processing;Then it uses
PH section after " moving average method " controls quality carries out average output;Finally using the pH value of desk-top pH meter measurement seawater as base
Standard averagely exports section to the pH of multi-parameter water quality section plotter observation and is corrected.The feature of the method for the present invention has: 1) multi-parameter
The sample frequency of water quality section plotter is 6Hz, available pH section more abundant;2) by pH aligning step, can improve
Because instrument response is drifted about the evaluated error of generation at any time, effectively guarantee the quality of data,It is provided for research seawater physicochemical property
Higher frequency and the support of accurate data.
Detailed description of the invention
Fig. 1 is pH bearing calibration flow chart
Fig. 2 is part of test results, in which:
(a) the spring layer sectional view of B09 erect-position
(b) probability density function and histogram frequency distribution diagram of Ratio
(c) Ratio and pHrLinear correlation curve
(d)pHbWith pHrLinear correlation curve
Fig. 3 is that yellow Bohai Sea in summer in 2017 shares the big mapping of flight number correction front and back bottom pH, in which:
(a) before correcting after (b) correction
Specific embodiment
The invention mainly comprises: data acquisition, data quality control, data are average and Data correction.Fig. 1 is provided
The block diagram of mentioned method.Include:
1, data acquisition: the pH section of the observation each website in sea area is obtained using multi-parameter water quality section plotter, and is indicated with S.
Instrument has used NIST (the National Institute of Standards of pH=4.01,7.00 and 9.21 before laying
And Technology, USA) standard buffer solution at room temperature to pH probe calibrated;Guarantee to lay during laying
Shipboard windward, is uniformly transferred, so as to instrument stabilizer measurement.The sample frequency of multi-parameter water quality section plotter is 6Hz, pH probe measurement
Range 1~13, precision 0.1 can obtain pH section more abundant;
2, data quality control: in order to exclude in observation process since ship motion, winch wirerope are swung and instrument is visited
Head unstability bring external interference needs to carry out quality control to the original pH section S that multi-parameter water quality section plotter measures,
Quality Control results are indicated with Q.Detailed process is as follows:
Step 1: removal temperature-sensitive data.Multi-parameter water quality section plotter need to stop about 1min at surface layer (~1m) after entering water, complete
At instrument temperature-sensitive process, this process is mostly false data, should give removal;
Step 2: inverse pressure correction data.Work as Pi-Pi-1≤ 0, reject PiPH data at the depth of water, in formula, i is a certain website
Data sequence, PiFor the pressure at i, Pi-1For the pressure at i-1;
Step 3: deglitch processing.PH section is calculated in the first derivative of corresponding depth and the Local standard deviation of derivative,
The signal that the two ratio is more than a certain given threshold is removed as spike, and is substituted with the average value of two end signal of spike.It is former
It manages as follows:
In formula, i is the data sequence of a certain website, ziFor the depth of water of corresponding position, xiIt is measured for multi-parameter water quality section plotter
Corresponding position pH data, DiFor the first derivative of corresponding depth, s is the standard deviation of corresponding depth first derivative.So thatWherein lim is variable, using the two ratio be more than lim this threshold value signal as spike removal, and with spike two
The average signal at end replaces.
3, data are average: in order to accurately indicate measurement result, random error caused by accidental fluctuation are eliminated, often to dynamic
Test data makees smooth and filtering processing.This patent is handled (processing result is indicated with A) to Q using " moving average method ",
0.2m physical phenomenon below is filtered out, to research and analyse in next step.The basic principle of moving average method:
yj=fj+ejJ=1,2 ..., N
In formula, yjFor dynamic test data, fiFor smoother measurement result, eiFor random error.
4, Data correction: to improve because instrument response is drifted about the evaluated error of generation at any time, the pH of higher precision is obtained
Section, mentioned method is on the basis of the seawater pH value that desk-top pH meter measures, to the pH section A of multi-parameter water quality section plotter observation
It is corrected.Detailed process is as follows:
Step 1: removal spring layer data.Thermohaline spring layer nearby due to water Inversion phenomenon suddenly and change dramatically, multi-parameter water quality
There may be response lag phenomenons for section plotter, therefore, it is necessary first to be removed to the pH data near each erect-position spring layer.
Using data such as the temperature of multi-parameter water quality section plotter simultaneous observation, salinity, the velocity of sound, the depth of waters, each erect-position is calculated
Buoyancy frequency N, and according to logN2> -3 determine the position of epilimnion in seawater, remove the pH data in A at spring layer, processing knot
Fruit is indicated with T.The calculation formula of buoyancy frequency are as follows:
In formula, c0For the velocity of sound, ρ is density of sea water, and z is the depth of water, and g is acceleration of gravity.g2/co 2One magnitude is very small, can
It ignores.
Step 2: the seawater pH value of desk-top pH meter measurement and the ratio of pH in corresponding T are calculated, and the ratio is defined as
Ratio has:
In formula, pHbFor the seawater pH value of desk-top pH meter measurement, pHrRefer to the seawater of corresponding multi-parameter water quality section plotter measurement
PH value.
Step 3: it is assumed that Ratio approximation Normal Distribution, calculates the probability density function of Ratio, and to discrete larger
Ratio carry out 1~3 times of standard deviation and reject, processing result is indicated with R.Based on above-mentioned it is assumed that will have 68% numeric distribution
There is the range within 1 standard deviation in distance average, about 95% numeric distribution has within 2 standard deviations in distance average
Range, about 99.7% numeric distribution have the range within 3 standard deviations in distance average.The calculation formula of probability density function
Are as follows:
It is denoted as X~N (μ, σ2),
Wherein, x is ratio R atio defined in step 2, and μ is the expectation of Ratio, and σ is the standard deviation of Ratio, σ2For
The variance of Ratio.
Step 4: data are observed with pH in corresponding T to R and carry out linear fit, the school of each erect-position pH section in calculating observation sea area
Positive coefficient.Y=ax+b is enabled, parameter a and b are used to adjust the slope and intercept of calibration curve, and x indicates that multi-parameter water quality section plotter is seen
The corresponding correction ratio of the different pH of pH section A, y expression of survey.This patent is right by taking yellow Bohai Sea in summer in 2017 shares flight number as an example
Discrete biggish Ratio makees 2 times of standard deviation rejecting processing, fitting gained unitary linearity curve are as follows:
Y=-0.0991x+1.8037 R2=0.8225
The Ratio and pH in regression curve it can be seen from R=0.9069rGood relationship, therefore use linear fitting
PH plane section correction suitable for multi-parameter water quality section plotter.
Step 5: A is multiplied with y one-to-one correspondence, obtains complete pH correction section.That is:
Y=A × y
In formula, A is the pH section of multi-parameter water quality section plotter observation before correcting, and Y is multi-parameter water quality section plotter after correction
The pH section of observation.
Linear fit is carried out to two sets of data of correction front and back, to examine whether pH correction front and back changes investigation sea area pH's
Spatial distribution, matched curve are as shown in Figure 2 d;Fig. 3 shows that the Bohai Sea Huanghai Sea summer in 2017 shares the big EDS maps of flight number bottom pH,
It can thus be seen that the spatial distribution of correction front and back pH does not change, but the quality of data significantly improves, and more meets normal
The scientific law that seawater pH is about 7.9~8.1.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, all in essence of the invention
Within mind and principle, any modification, equivalent substitution, improvement and etc. done be should all be included in the protection scope of the present invention.
Claims (1)
1. a kind of pH bearing calibration of multi-parameter water quality section plotter, including the following steps:
1) the pH section S of the observation each website in sea area data acquisition: is obtained using multi-parameter water quality section plotter;
2) data quality control: by rejecting the temperature-sensitive data in instrument continuous mode, inverse pressure correction data, calculating local standard
Difference removal large amplitude singular value and etc., quality control is carried out to S, Quality Control results are indicated with Q;
3) data are average: carrying out sliding average to erect-position downlink data each in Q, processing result is indicated with A;
4) Data correction: on the basis of the seawater pH value for the desk-top pH meter measurement calibrated by pH standard buffer solution, to A into
Row correction process, the method is as follows:
Step 1: removal spring layer data, using the temperature of multi-parameter water quality section plotter simultaneous observation, salinity, the velocity of sound, bathymetric data,
The buoyancy frequency N of each erect-position is calculated, and according to logN2> -3 determine the position of epilimnion in seawater, remove in A at spring layer
PH data, processing result indicate with T, the calculation formula of buoyancy frequency are as follows:
In formula, c0For the velocity of sound, ρ is density of sea water, and z is the depth of water, and g is acceleration of gravity, g2/co 2One magnitude is very small, ignores not
Meter;
Step 2: calculating the seawater pH value of desk-top pH meter measurement and the ratio of pH in corresponding T, and the ratio be defined as Ratio,
Have:
In formula, pHbFor the seawater pH value of desk-top pH meter measurement, pHrRefer to the seawater pH value of corresponding multi-parameter water quality section plotter measurement;
Step 3: it is assumed that Ratio approximation Normal Distribution, calculates the probability density function of Ratio, and to discrete biggish
Ratio carries out 1~3 times of standard deviation and rejects, and processing result is indicated with R;
Step 4: data are observed with pH in corresponding T to R and carry out linear fit, the correction system of each erect-position pH section in calculating observation sea area
Number, the method is as follows: enable y=ax+b, parameter a and b are used to adjust the slope and intercept of calibration curve, and x indicates that multi-parameter water quality cuts open
The pH section A of face instrument observation, i.e., comprising whole pH data including spring layer, the corresponding correction ratio of the different pH of y expression;
Step 5: A is multiplied with y one-to-one correspondence, obtains complete pH correction section.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910177098.0A CN110057990B (en) | 2019-03-08 | 2019-03-08 | PH correction method of multi-parameter water quality profiler |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910177098.0A CN110057990B (en) | 2019-03-08 | 2019-03-08 | PH correction method of multi-parameter water quality profiler |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110057990A true CN110057990A (en) | 2019-07-26 |
CN110057990B CN110057990B (en) | 2021-09-14 |
Family
ID=67316743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910177098.0A Expired - Fee Related CN110057990B (en) | 2019-03-08 | 2019-03-08 | PH correction method of multi-parameter water quality profiler |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110057990B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111680054A (en) * | 2020-06-08 | 2020-09-18 | 中国科学院地理科学与资源研究所 | Time sequence data storage method and device based on multi-parameter water quality monitor |
CN117405975A (en) * | 2023-12-14 | 2024-01-16 | 深圳鹏城新能科技有限公司 | Method, system and medium for detecting insulation resistance of PV panel |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7348783B1 (en) * | 2006-08-18 | 2008-03-25 | Chung Yuan Christian University | Multi-functional pH meter and fabrication thereof |
JP2011220980A (en) * | 2010-04-14 | 2011-11-04 | Central Res Inst Of Electric Power Ind | Measuring method of concentration of fluorescence dyestuff in sample solution |
CN103542875A (en) * | 2012-07-09 | 2014-01-29 | 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 | Arrangement and Method for Calibrating at Least Two Sensors in Parallel |
CN107607605A (en) * | 2017-08-08 | 2018-01-19 | 安吉元融仪器仪表检测有限公司 | The calibration method of accurate pH meter |
CN108368471A (en) * | 2015-10-30 | 2018-08-03 | 豪夫迈·罗氏有限公司 | The discriminating of the calibration deviation of PH measuring devices |
CN109060611A (en) * | 2018-07-09 | 2018-12-21 | 天津大学 | A kind of suspended sediment data processing method based on laser particle analyzer |
CN109141380A (en) * | 2018-09-19 | 2019-01-04 | 天津大学 | A kind of hydrographic data processing method based on CTD instrument |
-
2019
- 2019-03-08 CN CN201910177098.0A patent/CN110057990B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7348783B1 (en) * | 2006-08-18 | 2008-03-25 | Chung Yuan Christian University | Multi-functional pH meter and fabrication thereof |
JP2011220980A (en) * | 2010-04-14 | 2011-11-04 | Central Res Inst Of Electric Power Ind | Measuring method of concentration of fluorescence dyestuff in sample solution |
CN103542875A (en) * | 2012-07-09 | 2014-01-29 | 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 | Arrangement and Method for Calibrating at Least Two Sensors in Parallel |
CN108368471A (en) * | 2015-10-30 | 2018-08-03 | 豪夫迈·罗氏有限公司 | The discriminating of the calibration deviation of PH measuring devices |
CN107607605A (en) * | 2017-08-08 | 2018-01-19 | 安吉元融仪器仪表检测有限公司 | The calibration method of accurate pH meter |
CN109060611A (en) * | 2018-07-09 | 2018-12-21 | 天津大学 | A kind of suspended sediment data processing method based on laser particle analyzer |
CN109141380A (en) * | 2018-09-19 | 2019-01-04 | 天津大学 | A kind of hydrographic data processing method based on CTD instrument |
Non-Patent Citations (1)
Title |
---|
张开远 等: "基于最小二乘法的pH值温度补偿系统设计", 《传感器与微系统》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111680054A (en) * | 2020-06-08 | 2020-09-18 | 中国科学院地理科学与资源研究所 | Time sequence data storage method and device based on multi-parameter water quality monitor |
CN111680054B (en) * | 2020-06-08 | 2023-08-04 | 中国科学院地理科学与资源研究所 | Time sequence data storage method and device based on multi-parameter water quality monitor |
CN117405975A (en) * | 2023-12-14 | 2024-01-16 | 深圳鹏城新能科技有限公司 | Method, system and medium for detecting insulation resistance of PV panel |
CN117405975B (en) * | 2023-12-14 | 2024-03-22 | 深圳鹏城新能科技有限公司 | Method, system and medium for detecting insulation resistance of PV panel |
Also Published As
Publication number | Publication date |
---|---|
CN110057990B (en) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110146895B (en) | Acoustic velocity profile inversion method based on inverted multi-beam echometer | |
Mackas et al. | Least squares multiple tracer analysis of water mass composition | |
RiSeR et al. | Salinity in Argo: A modern view of a changing ocean | |
CN110673108B (en) | Airborne marine laser radar signal processing method based on iteration Klett | |
CN110057990A (en) | A kind of pH bearing calibration of multi-parameter water quality section plotter | |
CN115342814B (en) | Unmanned ship positioning method based on multi-sensor data fusion | |
CN115327518A (en) | Ocean laser radar multiple scattering attenuation correction method based on analytical model | |
Macaulay et al. | Field measurements of acoustic absorption in seawater from 38 to 360 kHz | |
Shao et al. | Verification of echosounder measurements of thickness and spatial distribution of kelp forests | |
CN109060611B (en) | Suspended sand mass concentration data processing method based on laser particle analyzer | |
CN116242584A (en) | Floating ocean platform along with ship wave measuring device based on BP neural network | |
CN110020409A (en) | Surface air temperature Structural Observed data Analysis method based on self-adaptive kernel density estimation algorithm | |
CN116864026A (en) | Ocean dissolved oxygen concentration reconstruction model construction method based on Argo warm salt profile | |
笹野大輔 et al. | Testing a new quick response oxygen sensor," RINKO" | |
RU2769550C1 (en) | Method for completing the part of the profile of vertical distribution of the speed of sound measured from the sea surface to the bottom | |
Beaudoin et al. | Streamlining sound speed profile pre-processing: Case studies and field trials | |
RU2747739C1 (en) | Verification method for sea water salinity measurements by automatic measuring instruments | |
CN113723182A (en) | SAR image ship detection method under limited training sample condition | |
Ciuffardi et al. | Deep-water hydrodynamic observations around a cold-water coral habitat in a submarine canyon in the eastern Ligurian Sea (Mediterranean Sea) | |
CN112114367A (en) | Method for detecting electric interface in sea area by using active electric field | |
CN114818232A (en) | Reconstruction method and system of seawater space temperature profile | |
CN104777526B (en) | Method for correcting ASCAT inversion wind speed | |
CN110568417A (en) | Airborne laser sounding echo signal processing method based on multiple physical fields | |
CN118114030B (en) | Internal solitary wave parameter extraction method based on single-point single-depth current meter data | |
Fujiki et al. | Improvement in estimation stability and accuracy of wave directional spectrum by a Bayesian method for swell observation |
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 | ||
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: 20210914 |