CN105301059A - Device and method for measuring gas-liquid cyclone liquid-phase standing time distribution - Google Patents

Device and method for measuring gas-liquid cyclone liquid-phase standing time distribution Download PDF

Info

Publication number
CN105301059A
CN105301059A CN201510712331.2A CN201510712331A CN105301059A CN 105301059 A CN105301059 A CN 105301059A CN 201510712331 A CN201510712331 A CN 201510712331A CN 105301059 A CN105301059 A CN 105301059A
Authority
CN
China
Prior art keywords
liquid
liquid phase
gas
phase
cyclone separator
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
Application number
CN201510712331.2A
Other languages
Chinese (zh)
Other versions
CN105301059B (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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201510712331.2A priority Critical patent/CN105301059B/en
Publication of CN105301059A publication Critical patent/CN105301059A/en
Application granted granted Critical
Publication of CN105301059B publication Critical patent/CN105301059B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Volume Flow (AREA)

Abstract

The invention belongs to the field of chemical engineering, and particularly relates to a device for measuring gas-liquid cyclone liquid-phase standing time distribution. A liquid-phase storage tank is connected with the inlet of a cyclone separator through a liquid-phase feed pipeline. A liquid-phase feed pump, a liquid-phase flow meter and a liquid-phase feed valve are arranged on the liquid-phase feed pipeline. The inlet of the cyclone separator is connected with a tracer agent injecting opening. The bottom liquid discharge opening of the cyclone separator is connected with a draught fan through a liquid discharge pipeline. A rectifying element, a conductivity probe, a liquid discharge separating tank, a liquid discharge branch gas-phase flow meter and a liquid discharge branch valve are arranged on the liquid discharge pipeline. The top gas discharge opening of the cyclone separator is connected with the draught fan through a gas discharge pipeline. A gas discharge branch flow meter and a gas discharge branch valve are arranged on the gas discharge pipeline. The conductivity probe is connected with a conductivity instrument which is connected with a computer and used for recording data to obtain liquid-phase standing time distribution. The device can measure liquid-phase standing time distribution under the condition of gas-liquid cyclone, and the defect that a conventional measuring method can not be used for dispersing the liquid phase is overcome.

Description

Measure the apparatus and method of gas-liquid eddy flow liquid phase residence time destribution
Technical field
The invention belongs to field of chemical engineering, particularly, relate to a kind of apparatus and method measuring gas-liquid eddy flow liquid phase residence time destribution, for measuring liquid phase residence time destribution under gas-liquid Swirl Condition, overcoming the limitation of conventional residence time assay method.
Background technology
Residence time destribution is the important parameter of research reactor or separation equipment performance, and residence time destribution can reflect the admixture of fluid in reaction or separation equipment, is such Optimum Design and the requisite theoretical foundation of operation.The mensuration of residence time destribution is not only widely used in Chemical Reaction Engineering and chemical separation process, and is applied to other field of design flow process.
Liquid phase residence time assay method conventional at present has rays method, fluorescence method, photoelectric method and electrical conductivity method etc.Wherein rays method is by detecting its concentration of radioactivity determination of tracer agent, and this method exists radiation pollution and requires higher to technical equipment, now seldom adopts; Fluorescence method is the fluorescence intensity utilizing fainter light detector to detect fluorescent tracer, and the preparation of this method fluorescent tracer is comparatively complicated, is mainly used in the measurement of polymkeric substance; Photoelectric method is by detecting the intensity of specific wavelength light through staining solution, measure its stain concentration, but this method light intensity easily by the impact of liquid film layer thickness, Gas-Liquid Dispersion state and liquid phase rotational flow, is not suitable for the measurement of gas liquid two-phase flow state.
Electrical conductivity method is the method for mensuration liquid phase residence time destribution the most frequently used at present.Its principle adds tracer agent (being generally electrolyte solution) in the liquid phase, by measuring conductivity determination tracer concentration in liquid phase, to obtain residence time destribution.But conductance measurement must ensure for Continuous Liquid Phase between two electrodes, be the state difficulty or ease Measurement accuracy of rotating flow, liquid throttling and spray flow for gas-liquid two-phase flow pattern.
At present, the research about gas-liquid cyclonic unit liquid phase residence time destribution rarely has report, and existing measuring liquid phase stay time distribution method is only applicable to pure liquid phase and measures, and rarely has report for the residence time destribution mensuration of liquid phase mixed flow state.During due to gas-liquid mixed flow, when especially liquid phase is disperse phase, difficulty in various degree is all caused for above-mentioned measuring method.
Summary of the invention
In order to overcome the defect that prior art exists, the invention provides a kind of apparatus and method measuring gas-liquid eddy flow liquid phase residence time destribution, for measuring liquid phase residence time destribution under gas-liquid Swirl Condition.
For achieving the above object, the present invention adopts following proposal:
Measure the device of gas-liquid eddy flow liquid phase residence time destribution, comprising: liquid phase storage tank, cyclone separator, conductivity probe, blower fan; Wherein: cyclone separator arranges liquid phase mixing entrance, top vent, bottom leakage fluid dram; Liquid phase storage tank is connected by liquid phase feeding pipeline with the entrance of cyclone separator, and liquid phase feeding pipeline is provided with liquid phase feeding pump, liquid phase stream gauge, liquid phase feeding valve successively by liquid phase storage tank to the direction of cyclone separator; The entrance of cyclone separator connects tracer agent inlet simultaneously, injects tracer agent by tracer agent inlet; The bottom leakage fluid dram of cyclone separator is connected with blower fan by tapping line, and tapping line is provided with rectifier cell, conductivity probe, discharge opeing separating tank, discharge opeing branch vapor phase stream gauge, discharge opeing branch valve successively by cyclone separator to blower fan direction; Cyclone separator top vent is connected with blower fan by exhaust line, and gas outlet to blower fan direction is provided with Exhaust branch flowmeter, Exhaust branch valve successively; Conductivity probe is connected with conductivity meter by signal wire, and the simulating signal that conductivity probe exports is converted to digital signal by conductivity meter, and is transferred to computing machine, carries out data record, to obtain liquid phase residence time destribution.
Preferably, tracer agent is generally NaCl or KCl solution, and tracer agent, by changing liquid phase conductivity, detects outlet conductivity variations by conductivity probe, obtains residence time destribution through data processing.
Relative to prior art, beneficial effect of the present invention is as follows: this device can measure liquid phase residence time destribution under gas-liquid Swirl Condition, overcomes the defect that routine measurement mode cannot be used for dispersion phase; Meanwhile, adopt branch flow to control, the liquid phase residence time destribution under different gas phase split ratio condition can be measured, to the measuring liquid phase stay time distribution of gas-liquid mixed flow system, there is applicability widely.
Accompanying drawing explanation
Fig. 1 is the device schematic diagram measuring gas-liquid eddy flow liquid phase residence time destribution;
Fig. 2 is conductivity measurement point partial enlarged drawing;
Fig. 3 is residence time measurement result example;
In figure: 1, liquid phase storage tank; 2, cyclone separator; 3, conductivity probe; 4, blower fan; 5, liquid phase feeding pump; 6, liquid phase stream gauge; 7, liquid phase feeding valve; 8, tracer agent inlet; 9, rectifier cell; 10, discharge opeing separating tank; 11, discharge opeing branch vapor phase stream gauge; 12, discharge opeing branch valve; 13, separating tank tapping valve; 14, Exhaust branch flowmeter; 15, Exhaust branch valve; 16, conductivity meter; 17, computing machine; 18, drop; 19, liquid film.
Embodiment
As shown in Figure 1 and Figure 2, measure the device of gas-liquid eddy flow liquid phase residence time destribution, comprising: liquid phase storage tank 1, cyclone separator 2, conductivity probe 3, blower fan 4; Wherein:
Cyclone separator 2 arranges liquid phase mixing entrance, top vent, bottom leakage fluid dram; The liquid phase blending agent entered by entrance is separated by cyclone separator, and partial gas phase flows out from top vent pipe, and residue gas phase and whole liquid phase flow out from bottom drain;
Liquid phase storage tank 1 is connected by liquid phase feeding pipeline with the entrance of cyclone separator 2, and this pipe liquid phase feeding pipeline is provided with liquid phase feeding pump 5, liquid phase stream gauge 6, liquid phase feeding valve 7 successively by liquid phase storage tank to the direction of cyclone separator; Liquid phase pump in liquid phase storage tank 1 is delivered to cyclone separator 2 entrance by liquid phase feeding pump 5, provides system liquid phase feeding, and liquid phase is measured by liquid phase stream gauge 6, controlled by liquid phase feeding valve 7; Liquid phase is separated by cyclone separator 2; The entrance of cyclone separator 2 connects tracer agent inlet 8 simultaneously, inject tracer agent by tracer agent inlet 8, tracer agent is generally NaCl or KCl solution, and tracer agent is by changing liquid phase conductivity, detect outlet conductivity variations by conductivity probe, obtain residence time destribution through data processing;
The bottom leakage fluid dram of cyclone separator 2 is connected with blower fan 4 by tapping line, and tapping line is provided with rectifier cell 9, conductivity probe 3, discharge opeing separating tank 10, discharge opeing branch vapor phase stream gauge 11, discharge opeing branch valve 12 successively by cyclone separator 2 to blower fan 4 direction; Dispersant liquid drop 18 rectification in discharging tube is the continuous liquid film 19 of adherent flowing by rectifier cell 9, is convenient to conductivity probe 3 and effectively measures liquid phase conductivity; Be provided with separating tank tapping valve 13 bottom discharge opeing separating tank 10, the liquid phase flowed out bottom cyclone 2 is separated with unsegregated gas phase by discharge opeing separating tank 10 again, and liquid phase is discharged by the tapping valve 13 of separating tank lower curtate, and gas phase is emptying after being inhaled into blower fan 4; Gas phase in discharging tube is measured by discharge opeing affluent-dividing gauge 11, is controlled by discharge opeing branch valve 12.
Cyclone separator 2 top vent is connected with blower fan 4 by exhaust line, and gas outlet is provided with Exhaust branch flowmeter 14, Exhaust branch valve 15 successively to blower fan 4 direction; In gas outlet, gas phase is measured by Exhaust branch flowmeter 14, is controlled by Exhaust branch valve 15.
Conductivity probe 3 is connected with conductivity meter 16 by signal wire, and the simulating signal that conductivity probe 3 exports is converted to digital signal by conductivity meter 16, and is transferred to computing machine 17, carries out data record, to obtain liquid phase residence time destribution.
Measure the method for gas-liquid eddy flow liquid phase residence time destribution, adopt the device of said determination gas-liquid eddy flow liquid phase residence time destribution, comprise the steps:
(1), start liquid phase feeder pump 5, liquid phase in liquid phase storage tank 1 is pumped into cyclone separator 2 entrance, adjustment liquid phase feeding valve 7 aperture, makes liquid phase stream gauge 6 displayed value reach experiment setting value, record liquid phase flow rate;
(2), blower fan 4 is started, gas phase is sucked by cyclone separator 2 entrance, adjustment Exhaust branch valve 15 and discharge opeing branch valve 12, to control gas phase split ratio, discharge opeing branch vapor phase stream gauge 11 and Exhaust branch flowmeter 14 is made to reach setting value, record exhaust and discharge opeing branch gas phase flow rate respectively, can calculate gas phase split ratio:
Gas phase split ratio=cyclone gas exhausting tube phase flow rate/cyclone inlet gas phase flow rate
(3), after device is stable, tracer agent is injected by tracer agent inlet 8, timing by computing machine simultaneously; Tracer agent is generally NaCl or KCl solution, injection mode can be that impulses injection or step are injected, grooved faces or plug-in type conductivity probe 3 is adopted to measure liquid phase conductivity, measuring-signal is sent to computing machine 17 through conductivity meter 16, carry out signal-obtaining and record, calculate residence time destribution by following formula
E ( t ) = C i ΣC i Δt i
Wherein, C ifor each moment conductivity of computer recording, Δ t ifor data step-length writing time, E (t) is residence time destribution probability.
Figure 3 shows that cyclone separator enters implication speed 16m/s, the liquid phase residence time destribution result measured during liquid phase volume content 1%.

Claims (5)

1. measure a device for gas-liquid eddy flow liquid phase residence time destribution, comprising: liquid phase storage tank, cyclone separator, conductivity probe, blower fan; It is characterized in that: cyclone separator arranges liquid phase mixing entrance, top vent, bottom leakage fluid dram; Liquid phase storage tank is connected by liquid phase feeding pipeline with the entrance of cyclone separator, and liquid phase feeding pipeline is provided with liquid phase feeding pump, liquid phase stream gauge, liquid phase feeding valve successively by liquid phase storage tank to the direction of cyclone separator; The entrance of cyclone separator connects tracer agent inlet simultaneously, injects tracer agent by tracer agent inlet; The bottom leakage fluid dram of cyclone separator is connected with blower fan by tapping line, and tapping line is provided with rectifier cell, conductivity probe, discharge opeing separating tank, discharge opeing branch vapor phase stream gauge, discharge opeing branch valve successively by cyclone separator to blower fan direction; Cyclone separator top vent is connected with blower fan by exhaust line, and gas outlet to blower fan direction is provided with Exhaust branch flowmeter, Exhaust branch valve successively; Conductivity probe is connected with conductivity meter by signal wire, and the simulating signal that conductivity probe exports is converted to digital signal by conductivity meter, and is transferred to computing machine, carries out data record, to obtain liquid phase residence time destribution.
2. the device of mensuration gas-liquid eddy flow liquid phase residence time destribution according to claim 1, is characterized in that: tracer agent is generally NaCl or KCl solution.
3. measure a method for gas-liquid eddy flow liquid phase residence time destribution, adopt the device of the mensuration gas-liquid eddy flow liquid phase residence time destribution described in claim 1-2, it is characterized in that, comprise the steps:
(1), start liquid phase feeder pump, liquid phase in liquid phase storage tank is pumped into cyclone separator entrance, adjustment liquid phase feeding valve opening, makes liquid phase stream gauge displayed value reach experiment setting value, record liquid phase flow rate;
(2), blower fan is started, gas phase is sucked by cyclone separator entrance, adjustment Exhaust branch valve and discharge opeing branch valve, to control gas phase split ratio, discharge opeing branch vapor phase stream gauge and Exhaust branch flowmeter is made to reach setting value, record exhaust and discharge opeing branch gas phase flow rate respectively, can calculate gas phase split ratio:
Gas phase split ratio=cyclone gas exhausting tube phase flow rate/cyclone inlet gas phase flow rate
(3), after device is stable, tracer agent is injected by tracer agent inlet, timing by computing machine simultaneously; Adopt grooved faces or plug-in type conductivity probe to measure liquid phase conductivity, measuring-signal is sent to computing machine through conductivity meter, carries out signal-obtaining and record, calculates residence time destribution by following formula
E ( t ) = C i ΣC i Δt i
Wherein, C ifor each moment conductivity of computer recording, Δ t ifor data step-length writing time, E (t) is residence time destribution probability.
4. the method for mensuration gas-liquid eddy flow liquid phase residence time destribution according to claim 3, is characterized in that, the dispersant liquid drop rectification in discharging tube is the continuous liquid film of adherent flowing by employing rectifier cell.
5. the method for the mensuration gas-liquid eddy flow liquid phase residence time destribution according to claim 3-4, it is characterized in that, conductivity probe used is grooved faces or plug-in type conductivity probe.
CN201510712331.2A 2015-10-28 2015-10-28 The apparatus and method for determining gas-liquid eddy flow liquid phase residence time destribution Expired - Fee Related CN105301059B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510712331.2A CN105301059B (en) 2015-10-28 2015-10-28 The apparatus and method for determining gas-liquid eddy flow liquid phase residence time destribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510712331.2A CN105301059B (en) 2015-10-28 2015-10-28 The apparatus and method for determining gas-liquid eddy flow liquid phase residence time destribution

Publications (2)

Publication Number Publication Date
CN105301059A true CN105301059A (en) 2016-02-03
CN105301059B CN105301059B (en) 2018-01-12

Family

ID=55198567

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510712331.2A Expired - Fee Related CN105301059B (en) 2015-10-28 2015-10-28 The apparatus and method for determining gas-liquid eddy flow liquid phase residence time destribution

Country Status (1)

Country Link
CN (1) CN105301059B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932445A (en) * 2017-03-03 2017-07-07 四川龙蟒钛业股份有限公司 A kind of method of mean residence time and distribution of the test titanium white slurry in sand mill
CN107421999A (en) * 2017-06-19 2017-12-01 天津大学 Biphase gas and liquid flow plays rotation phase-splitter
CN108256197A (en) * 2018-01-10 2018-07-06 四川大学 The determining method of water body residence time in energy dissipation tank
CN110220556A (en) * 2019-06-25 2019-09-10 西安交通大学 A kind of biphase gas and liquid flow double parameter transient measurement system
CN110672542A (en) * 2019-09-26 2020-01-10 清华大学 Device and method for on-line measuring residence time distribution in liquid-liquid and gas-liquid continuous reactor

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5390547A (en) * 1993-11-16 1995-02-21 Liu; Ke-Tien Multiphase flow separation and measurement system
US20070084340A1 (en) * 2003-05-16 2007-04-19 Jianwen Dou Adjustable gas-liquid centrifugal separator and separating method
CN101793854A (en) * 2010-02-10 2010-08-04 中北大学 Device and method for online measurement of residence time distribution of rotating packed bed
CN101839395A (en) * 2009-10-27 2010-09-22 梁法春 Gas-liquid two-phase flow uniform distribution device
CN202142228U (en) * 2011-07-06 2012-02-08 浙江大学宁波理工学院 Experimental apparatus for measuring liquid phase stay time distribution in gas-liquid reactor
CN103239894A (en) * 2012-05-21 2013-08-14 中国石油天然气股份有限公司 Oil-gas mixed transportation separating flow divider and oil-gas separating flow-dividing method
CN103585934A (en) * 2013-09-29 2014-02-19 山东百川同创能源有限公司 Reactor suitable for cracking liquid raw materials of micro fluidized bed reaction analyzer
CN103979668A (en) * 2014-06-05 2014-08-13 陕西瑞科特种设备技术有限公司 Device for reacting discharge plasma positioned on gas-liquid interface of swirl gas column
GB2511262A (en) * 2011-12-22 2014-08-27 Statoil Petroleum As Method and system for fluid separation with an integrated control system

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5390547A (en) * 1993-11-16 1995-02-21 Liu; Ke-Tien Multiphase flow separation and measurement system
US20070084340A1 (en) * 2003-05-16 2007-04-19 Jianwen Dou Adjustable gas-liquid centrifugal separator and separating method
CN101839395A (en) * 2009-10-27 2010-09-22 梁法春 Gas-liquid two-phase flow uniform distribution device
CN101793854A (en) * 2010-02-10 2010-08-04 中北大学 Device and method for online measurement of residence time distribution of rotating packed bed
CN202142228U (en) * 2011-07-06 2012-02-08 浙江大学宁波理工学院 Experimental apparatus for measuring liquid phase stay time distribution in gas-liquid reactor
GB2511262A (en) * 2011-12-22 2014-08-27 Statoil Petroleum As Method and system for fluid separation with an integrated control system
CN103239894A (en) * 2012-05-21 2013-08-14 中国石油天然气股份有限公司 Oil-gas mixed transportation separating flow divider and oil-gas separating flow-dividing method
CN103585934A (en) * 2013-09-29 2014-02-19 山东百川同创能源有限公司 Reactor suitable for cracking liquid raw materials of micro fluidized bed reaction analyzer
CN103979668A (en) * 2014-06-05 2014-08-13 陕西瑞科特种设备技术有限公司 Device for reacting discharge plasma positioned on gas-liquid interface of swirl gas column

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
杨军卫 等: "气液旋流器内液相平均停留时间研究[J/OL]", 《化工学报》 *
梁正熙 等: "电导率仪与PC-1500联机测定液相停留时间分布", 《化学工程》 *
陈鸥: "气液两相流中液相局部速度电导探针法测量实验研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础学科辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932445A (en) * 2017-03-03 2017-07-07 四川龙蟒钛业股份有限公司 A kind of method of mean residence time and distribution of the test titanium white slurry in sand mill
CN107421999A (en) * 2017-06-19 2017-12-01 天津大学 Biphase gas and liquid flow plays rotation phase-splitter
CN108256197A (en) * 2018-01-10 2018-07-06 四川大学 The determining method of water body residence time in energy dissipation tank
CN108256197B (en) * 2018-01-10 2021-10-22 四川大学 Method for determining residence time of water body in energy dissipation pool
CN110220556A (en) * 2019-06-25 2019-09-10 西安交通大学 A kind of biphase gas and liquid flow double parameter transient measurement system
CN110672542A (en) * 2019-09-26 2020-01-10 清华大学 Device and method for on-line measuring residence time distribution in liquid-liquid and gas-liquid continuous reactor

Also Published As

Publication number Publication date
CN105301059B (en) 2018-01-12

Similar Documents

Publication Publication Date Title
CN105301059A (en) Device and method for measuring gas-liquid cyclone liquid-phase standing time distribution
CN105572307B (en) The calibration tooling and scaling method of a kind of gas sensor
US20090139345A1 (en) Isokinetic sampling method and system for multiphase flow from subterranean wells
CN208296945U (en) Electronic scale Static mass method water flow standard apparatus used in parallel
CN105004761A (en) Denitration system escaping ammonia on-line continuous monitoring device and method
CN104989371A (en) Online water content analyzing device for oil well mouth
CN101413817A (en) Double-differential pressure throttle humid gas measuring method
CN108333299A (en) A kind of marine main engine discharges pollutants quick precise testing device and method
CN203657895U (en) Liquid flow verification system for flow meter
CN211235687U (en) Chromatographic analysis system for separating methanol, acetaldehyde, ethylene oxide and chloroethylene
CN204877436U (en) Online moisture analytical equipment of oil well wellhead
CN110297077B (en) Lubricating oil moisture content measuring system and method based on Laval nozzle
CN104730216B (en) A kind of analysis metering device and liquid analysis system
CN204925005U (en) Flow cell
CN104360090B (en) Soil automatic analyzer and the method with this analysis-e/or determining Trace Elements in Soil
CN206627464U (en) The detection means of gas dew point
CN101603889B (en) Pointwise scanning type micro-fluid metering device and method
CN109752507A (en) A kind of rock core displacement test oil-gas-water metering device
CN102621039B (en) Method for measuring gaseous diffusion coefficient with double-flow gaseous quasi-static process
CN108732226A (en) Thioether class gas-detecting device and method
CN204575616U (en) A kind of analysis measuring apparatus and liquid analysis system
CN107561209A (en) Sulfur hexafluoride gas acidity detection method and device
CN209085690U (en) Based on the Z-type natural gas moisture real-time measurement apparatus being mutually separated in pipe
CN208223567U (en) A kind of integral type gas-liquid biphase flowmeter
CN206177375U (en) Automatic measure device of agricultural shower nozzle droplet flow

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20180112

Termination date: 20211028

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