CN107045050B - The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat - Google Patents

The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat Download PDF

Info

Publication number
CN107045050B
CN107045050B CN201710248685.5A CN201710248685A CN107045050B CN 107045050 B CN107045050 B CN 107045050B CN 201710248685 A CN201710248685 A CN 201710248685A CN 107045050 B CN107045050 B CN 107045050B
Authority
CN
China
Prior art keywords
water
value
pollutant
flow rate
model
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
Application number
CN201710248685.5A
Other languages
Chinese (zh)
Other versions
CN107045050A (en
Inventor
丁晓雯
康博
翟爱丰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN201710248685.5A priority Critical patent/CN107045050B/en
Publication of CN107045050A publication Critical patent/CN107045050A/en
Application granted granted Critical
Publication of CN107045050B publication Critical patent/CN107045050B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/18Water
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02Alarms for ensuring the safety of persons
    • G08B21/12Alarms for ensuring the safety of persons responsive to undesired emission of substances, e.g. pollution alarms
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/18Status alarms
    • G08B21/182Level alarms, e.g. alarms responsive to variables exceeding a threshold
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/20Controlling water pollution; Waste water treatment

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Toxicology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Water Treatment By Sorption (AREA)

Abstract

The invention discloses a kind of grading forewarning system methods that flow rate threshold is utilized under pop-up threat, which comprises the following steps: 1) divides to pollutant warning grade, using hierarchy model (A)2) flow rate threshold is calculated using the pollutant concentration monitor value c detected, determines dispersal pattern model (1) of the pollutant in river,(1) 3) when pop-up threat occurs, pollutant concentration value detected by the detector positioned at reservoir is transferred to early warning hierarchy model and flow rate threshold computation model, divides to pollutant degree of danger, determines warning level.

Description

The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat
Technical field
The invention belongs to environmental projects and ecological engineering technology field, and in particular to stream is utilized under a kind of pop-up threat The grading forewarning system method of fast threshold value simultaneously copes with pop-up threat using this method.
Background technique
In recent years, Three Gorges Reservoir upstream establishes more medium-sized chemical plant, and considerably increase river water source ground contamination can It can property.In addition, data are shown, oneself has been more than 3500 tons to the annual petroleum transportation in Three Gorges, therefore how research ensures water head site The water quality safety of drinking water is highly important.According to reservoir area of Three Gorges potable water source district Water Environment Status quo, commented in conjunction with water quality safety Sentence system, constructs reservoir area of Three Gorges Water Quality Safety of Drinking Water Source early-warning and predicting model framework.Based on water quality model to pollutant It is diffused simulation, establishes water quality safety forecast submodel;Based on water monitoring data, warning grade division is carried out, is established Water quality safety early warning submodel;Using MATLAB powerful calculating and graphic capability, reservoir area of Three Gorges potable water source district is constructed Water quality safety early-warning and predicting model.
Summary of the invention
To achieve the goals above, the technical solution of the present invention is as follows: utilizing flow rate threshold under a kind of pop-up threat Grading forewarning system method, comprising the following steps:
1) pollutant warning grade is divided, using hierarchy model (A).
Wherein, c- monitor value;c1Background value;c2Standard value;c3Pollutant average value
The monitor value is automatic control station surveyed Mixing Coefficient in Rectangular Channels water concentration numerical value when being occurred by accident.The back Scape value is the average value in water source area pollutant concentration many years monitoring data.The standard value is in water environment quality standard (GB3838-2002) meet the requirement of third level in.
When the functional value EW (C) of the hierarchy model (A) obtained according to pollutant concentration monitor value c value is less than 0, drink Water head site starts water quality safety primary warning;When EW (C) is between 0 and 1, it is pre- that potable water source district starts water quality safety middle rank It is alert;When Early-warning Model functional value EW (C) is greater than 1, potable water source district starts the advanced early warning of water quality safety, stops water intaking.
2) flow rate threshold is calculated using the pollutant water concentration monitor value c detected
Determine dispersal pattern model of the pollutant in river, with formula (l } indicate
Assuming that pollution sources are located between two boundaries there are two boundaries, consider the reflex on boundary, can lead at this time The reflex that the virtual source assumed carrys out simulating boundary is crossed, model (1) is changed to formula (2)
c≤c1
ux>0
Constraint condition: b < B
The calculation formula of each parameter is as follows in model:
Wherein: mono- monitor value of c, mg/l;
c1One national water standard concentration value, mg/l;
c0One water body background concentration value, mg/l;
Mono- pollutant quality of M, g;
The mono- river depth of water of h, m;
Mono- limiting time of t, s;
The diffusion coefficient in the direction Dx-X;
The diffusion coefficient in the direction Dy-Y;
Mono- attenuation coefficient of K, d-1
uX-X direction flow velocity, m/s;
Mono- pollution sources of b are at a distance from bank, m;
Mono- river width of B, m;
Mono- hydraulic gradient of i.
Using the model determination in the case of pop-up threat, by determining that the depth of water and limiting time etc. constrain item Part is meeting national water standard (under the requirement of water environment quality standard (GB 3838-2002 >), when calculating restriction Interior to reach flow rate of water flow required for national water standard, when flow velocity is greater than flow velocity door screen value, river water can pass through self-cleaning Ability reaches three classes water standard in 5 days.If flow velocity is less than flow velocity threshold value, can intuitively it judge within the crash time River cannot solve the problems, such as water pollution by itself self-purification capacity.Then need to provide a response item for scheduling aspect Part.
3) when pop-up threat occurs, pollutant concentration value detected by the detector positioned at reservoir is transferred to early warning Hierarchy model and flow rate threshold computation model, divide pollutant degree of danger, the stream for determining warning level, while obtaining Fast threshold value is compared with flow velocity actual in reservoir, if VIt is practical> VThreshold, reservoir passes through self-cleaning can reach water quality mark in 5 days Quasi- value, if VIt is practical< VThreshold, river cannot solve the problems, such as water pollution by itself self-purification capacity, should take other measures such as object Reason, chemistry and biological control measure are handled.
Effect of the invention, the present invention utilize existing numerical model and Visualization Platform, and inverting typical case's potable water source district exists Suddenly accident scene under different water environment situations, and grade classification, water are carried out to the pollution level of unexpected environmental accident When matter contamination accident happens suddenly, concentration value measured by monitoring station is carried out by the monitoring data of water quality early-warning model, pollutant concentration Analysis, can determine warning level, subsequent urgency signal can be published to relevant departments, take the measure of corresponding warning level. And consider that concentration issues different alarms from flow velocity, for different grades of sudden water pollution event, in conjunction with related physical, Chemistry and biological control measure, are administered.
Detailed description of the invention
Fig. 1 is grading forewarning system schematic diagram.
The relational graph of pollutant concentration and flow velocity when Fig. 2 is flood season.
The relational graph of pollutant concentration and flow velocity when Fig. 3 is the water storage phase
Fig. 4 is the relational graph of pollutant concentration and flow velocity when falling the phase of disappearing.
Fig. 5 is the combination figure for being classified alarm and flow rate threshold.
Fig. 6 is that reservoir area of Three Gorges water quality safety forecasts submodel figure.
Specific embodiment
With reference to the accompanying drawing, the present invention is described in more detail.
The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat, comprising the following steps:
1) pollutant warning grade is divided, the present Research and its hair of the current Water Quality Evaluation method of comprehensive analysis Exhibition trend, selection operation is simple, and operation is convenient and meets the progress warning level division of the method for early warning of reservoir area of Three Gorges actual conditions. Inventor creates hierarchy model (A).
Wherein, c- monitor value;c1Background value;c2Standard value;c3Pollutant average value.The monitor value is to pass through accident Automatic control station surveyed Mixing Coefficient in Rectangular Channels water concentration numerical value when generation.The background value is more in water source area pollutant concentration The average value of year monitoring data.The standard value meets third level in water environment quality standard (GB3838-2002) Requirement.
Fig. 1 is grading forewarning system schematic diagram, as the Early-warning Model functional value EW obtained according to pollutant concentration monitor value c value (C) less than 0 when, potable water source district starts water quality safety primary warning.When EW (C) is between 0 and 1, potable water source district is opened Dynamic water quality safety middle rank early warning, shows water head site there are risk of exceeding criterion, relevant departments carry out Emergency management preparation, consider to stop It only fetches water, carries out drinking water deposit work, it is ensured that the just common water of resident.When Early-warning Model functional value EW (C) be greater than 1, drinking water Start to source the advanced early warning of water quality safety, shows that water head site is exceeded, the safety of the water quality safety of potable water source district to the mankind It is formed dangerous, it is necessary to stop water intaking, relevant departments take urgent measure improvement at once, strive for making water head site drinking water quality as early as possible Reach safe concentration.
2) flow rate threshold is calculated using the pollutant water concentration monitor value c detected.
First determine dispersal pattern model of the pollutant in river, can with formula (1 } indicate
Assuming that pollution sources are located between two boundaries there are two boundaries, consider the reflex on boundary, can lead at this time The reflex that the virtual source assumed carrys out simulating boundary is crossed, model (1) is changed to formula (2)
c≤c1
ux>0
ux>0
Constraint condition: b < B
The calculation formula of each parameter is as follows in model:
Wherein: mono- monitor value of c, mg/l;
c1One national water standard concentration value, mg/l;
c.One water body background concentration value, mg/l;
Mono- pollutant quality of M, g;
The mono- river depth of water of h, m;
Mono- limiting time of t, s;
The diffusion coefficient in the direction Dx-X;
The diffusion coefficient in the direction Dy-Y;
Mono- attenuation coefficient of K, d-1
uX-X direction flow velocity, m/s;
Mono- pollution sources of b are at a distance from bank, m;
Mono- river width of B, m;
Mono- hydraulic gradient of i.
The present invention is for simulating the strain of the Changjiang river Yichang and return water head site Oil spills accident, wherein data used for 2015 times Hydrographic data is calculated, and flood season is May to August part, and the water storage phase is September And October, and the phase that falls that disappears is April November to next year. Corresponding river width B is 2000m, and the decay coefficient K of petroleum is 0.0000015, and contamination accident chooses river center, i.e. b= 1000m, national three classes water standard cl=O.OSmg/l.By calculating, flood season depth of water h=133.04m, Dx=61.11, Dy= l.ss;Water storage phase depth of water h=149.33m, Dx=72.7, Dy=1.84;Disappear and falls phase depth of water h=153.74m, Dx=75.91, Dy =1.92.
According to Three Gorges hydrological variation, being divided into three periods is calculated, with the feasibility of confirmatory experiment.Fig. 1, figure 2 and Fig. 3 respectively represents flood season, water storage phase and the flow velocity threshold value for falling in period phase and guaranteeing drinking water safety that disappears.In required scheduling It in limit, is calculated by simulation, can obtain the river flow thresholds of the guarantee water quality reaching standard in each period, such as this simulation case In, the emergency scheduling time be 5 days, therefore in the calculating simulation calculating time of flow velocity threshold value be limited with 5 days, in simulation process Pollutant selects petroleum, quality 600t.By calculate, the flow velocity threshold value in flood season, water storage phase and the phase that falls that disappears be distributed as 0.017m/s, 0.014m/s and 0.013mls.When the flow velocity in the period is greater than flow velocity threshold value, river water can be by from net energy Power reaches three classes water standard in 5 days.If flow velocity is less than flow velocity threshold value, can intuitive judgment go out river within the crash time Stream cannot solve the problems, such as water pollution by itself self-purification capacity.Need to provide a response condition for scheduling aspect.
3) Fig. 5 is the combination figure for being classified alarm and flow rate threshold.When pop-up threat occurs, positioned at the detector of reservoir Detected pollutant concentration value is transferred to early warning hierarchy model and flow rate threshold computation model, to pollutant degree of danger into Row divides, and the flow rate threshold for determining warning level, while obtaining is compared with flow velocity actual in reservoir, if VIt is practical> VThreshold, water Library can reach water standard value by self-cleaning in 5 days, if VIt is practical< VThreshold, river cannot solve water by itself self-purification capacity The problem of matter pollutes, should take other measures such as physics, chemistry and biological control measure to be handled.
Physical method: absorption method goes the organic contaminations such as benzene homologues, phenols, pesticide in water removal using adsorbent materials such as active carbons Object.Activated carbon application is extensive at present, can cope with 60 Some Organic Pollutants.Active carbon is divided into Powdered Activated Carbon (PAC) and granular Active carbon (GAC).In water pollution in Songhuajiang River accident urban water supply emergency processing in 2005, PAC absorption source water nitre is formd Base benzene and processing water factory's sand filter increase the emergency handling process of GAC filtering layer dual safety barrier newly, and nitro phenenyl concentration is full after processing Sufficient water standard.
Chemical method: the chemical emergency disposal measure for burst water pollution accident is mainly oxidative decomposition process.Oxygenolysis Method uses the oxidants such as potassium permanganate, ozone by organic pollutants oxidation removal.
Bioanalysis: biological treatment be by the phagocytosis of microbe matter, by the effects of metabolism by it It decomposes and converts, to achieve the purpose that remove organic matter.
Technical solution of the present invention is described in detail in above-described embodiment.It is apparent that the present invention is not limited being retouched The embodiment stated.Based on the embodiments of the present invention, those skilled in the art can also make a variety of variations accordingly, but appoint What is equal with the present invention or similar variation shall fall within the protection scope of the present invention.

Claims (4)

1. utilizing the grading forewarning system method of flow rate threshold under a kind of pop-up threat, which comprises the following steps:
1) pollutant warning grade is divided, using hierarchy model (A),
Wherein, c- monitor value;c1Background value;c2Standard value;c3Pollutant average value;
The monitor value is the surveyed Mixing Coefficient in Rectangular Channels water concentration numerical value of automatic control station, the background value when being occurred by accident It is in the average value of water source area pollutant concentration many years monitoring data, the standard value is in water environment quality standard (GB3838-2002) meet the requirement of third level in;
2) flow rate threshold is calculated using the pollutant water concentration monitor value c detected, determines diffusion of the pollutant in river Pattern model is indicated with formula (l)
Model (1) is changed to formula (2)
c≤c1
ux>0
Constraint condition: b < B
The calculation formula of each parameter is as follows in model:
Wherein: c-monitor value, mg/l;
c1- country water standard concentration value, mg/l;
C-water body background concentration value, mg/l;
M-pollutant quality, g;
H-river the depth of water, m;
T-limiting time, s;
DxThe diffusion coefficient of -X direction;
DyThe diffusion coefficient of -Y direction;
K-attenuation coefficient, d-1
uX-X direction flow velocity, m/s;
B-pollution sources are at a distance from bank, m;
B-river width, m;
I-hydraulic gradient;
3) when pop-up threat occurs, pollutant concentration value detected by the detector positioned at reservoir is transferred to early warning classification Model and flow rate threshold computation model, divide pollutant degree of danger, determine warning level;
In step 1), when the functional value EW (c) of the hierarchy model (A) obtained according to pollutant concentration monitor value c value is less than 0 When, potable water source district starts water quality safety primary warning;When EW (c) is between 0 and 1, potable water source district starts water quality peace Full middle rank early warning;When Early-warning Model functional value EW (c) is greater than 1, potable water source district starts the advanced early warning of water quality safety, stops taking Water.
2. method according to claim 1, which is characterized in that in step 2), when flow velocity is greater than flow rate threshold, river water Body can reach three classes water standard by self-purification capacity in 5 days;If flow velocity is less than flow rate threshold, can intuitively it judge The river within the crash time cannot solve the problems, such as water pollution by itself self-purification capacity out, then need to mention for scheduling aspect For response condition.
3. method according to claim 1, which is characterized in that actual stream in the flow rate threshold and reservoir obtained in step 3) Speed is compared, if VIt is practical> VThreshold, reservoir can reach water standard value by self-cleaning in 5 days, if VIt is practical< VThreshold, river cannot Water pollution is solved the problems, such as by itself self-purification capacity, should be taken at other measures such as physics, chemistry and biological control measure Reason.
4. method according to claim 3, which is characterized in that the physical measure is absorption method, is gone using absorbent charcoal material Benzene homologues, phenols, insecticide pollution in water removal;The chemistry measure is using the oxidative decomposition process of potassium permanganate, using ozone Oxidant method;The biological control measure is to be decomposed it by metabolism by the phagocytosis of microbe matter And conversion, achieve the purpose that remove organic matter.
CN201710248685.5A 2017-04-17 2017-04-17 The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat Expired - Fee Related CN107045050B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710248685.5A CN107045050B (en) 2017-04-17 2017-04-17 The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710248685.5A CN107045050B (en) 2017-04-17 2017-04-17 The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat

Publications (2)

Publication Number Publication Date
CN107045050A CN107045050A (en) 2017-08-15
CN107045050B true CN107045050B (en) 2019-03-12

Family

ID=59544939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710248685.5A Expired - Fee Related CN107045050B (en) 2017-04-17 2017-04-17 The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat

Country Status (1)

Country Link
CN (1) CN107045050B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107403034B (en) * 2017-06-28 2021-08-17 西交利物浦大学 Water quality pollution diffusion early warning visualization method
CN111524035B (en) * 2020-05-13 2021-05-14 水利部交通运输部国家能源局南京水利科学研究院 Hydrodynamic regulation and control threshold determination method for improvement of river network water environment in plain city
CN115684516B (en) * 2022-08-29 2024-03-26 邯郸市亿润工程咨询有限公司 Hydraulic engineering sewage treatment detection method and device, electronic equipment and medium

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105976086A (en) * 2016-04-26 2016-09-28 天津大学 Early warning and plan generating method for sudden water pollution event of water conveyance project

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105976086A (en) * 2016-04-26 2016-09-28 天津大学 Early warning and plan generating method for sudden water pollution event of water conveyance project

Also Published As

Publication number Publication date
CN107045050A (en) 2017-08-15

Similar Documents

Publication Publication Date Title
Varol Spatio-temporal changes in surface water quality and sediment phosphorus content of a large reservoir in Turkey
CN107045050B (en) The grading forewarning system method of flow rate threshold is utilized under a kind of pop-up threat
Zhai et al. Integrated approach of hydrological and water quality dynamic simulation for anthropogenic disturbance assessment in the Huai River Basin, China
Muñoz-Barbosa et al. Relationship between metal enrichments and a biological adverse effects index in sediments from Todos Santos Bay, northwest coast of Baja California, México
Zhu et al. Spatial distribution and contamination assessment of heavy metals in surface sediments of the Caofeidian adjacent sea after the land reclamation, Bohai Bay
Issakhov et al. Numerical modeling of water pollution by products of chemical reactions from the activities of industrial facilities at variable and constant temperatures of the environment
López et al. Estuaries as filters for riverine microplastics: Simulations in a large, coastal-plain estuary
Xu et al. Impact of cascade reservoirs on continuity of river water temperature: A temperature trend hypothesis in river
Paz-Alberto et al. Climate change vulnerability and disaster risk assessment using remote sensing technology and adaptation strategies for resiliency and disaster risk management in selected coastal municipalities of Zambales, Philippines
Park et al. Impact of climate change on the persistent turbidity issue of a large dam reservoir in the temperate monsoon region
Long et al. Saltwater intrusion induced by a complete neap tide and its effect on nutrients variation in the estuary of Pearl River, China
Li et al. Status and trends of sediment metal pollution in Bohai Sea, China
Chen et al. Research on marine disaster prevention and mitigation information platform system based on big data
Zhou et al. Study on health risk assessment of potentially toxic elements in the soil around landfill site in Shannan City, Tibet
Hanslík et al. Observed half-lives of 3H, 90Sr and 137Cs in hydrosphere in the Vltava River basin (Bohemia)
Zhao et al. Water quality assessment, possible origins and health risks of toxic metal (loid) s in five cascade reservoirs in the upper Mekong
Liu et al. Using network to enhance the insights on correlation and pollution assessment of co-occurring metals in marine sediments, the East China Sea
Jacob et al. Environmental isotopes to test hypotheses for fluid mud (mud bank) generation mechanisms along the southwest coast of India
Qiao et al. Fast and optimal decision for emergency control of sudden water pollution accidents in long distance water diversion projects
de Castro Barcellos et al. The role of the Arroio Pavuna river in the transport of particulate heavy metals to Jacarepagua lagoon, Brazil
Devkota et al. Response characteristics of the Perdido and Wolf Bay system to inflows and sea level rise.
Tian et al. Water quality's responses to water energy variability of the Yangtze River
Tri et al. Using numerical modelling in the simulation of mass fish death phenomenon along the Central Coast of Vietnam
Nikanorov et al. Tendencies of long-term changes in water quality of water bodies in the South of Russia
Ziying et al. Research of risk assessment system on tailings pond water pollution

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

Granted publication date: 20190312

CF01 Termination of patent right due to non-payment of annual fee