CN105740574B - A kind of ancient landform section restoration methods - Google Patents

A kind of ancient landform section restoration methods Download PDF

Info

Publication number
CN105740574B
CN105740574B CN201610121264.1A CN201610121264A CN105740574B CN 105740574 B CN105740574 B CN 105740574B CN 201610121264 A CN201610121264 A CN 201610121264A CN 105740574 B CN105740574 B CN 105740574B
Authority
CN
China
Prior art keywords
well
thickness
ancient landform
karst
obtaining
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
CN201610121264.1A
Other languages
Chinese (zh)
Other versions
CN105740574A (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 Beijing
Original Assignee
China University of Petroleum Beijing
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 Beijing filed Critical China University of Petroleum Beijing
Priority to CN201610121264.1A priority Critical patent/CN105740574B/en
Publication of CN105740574A publication Critical patent/CN105740574A/en
Application granted granted Critical
Publication of CN105740574B publication Critical patent/CN105740574B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention relates to a kind of ancient landform section restoration methods, it is characterised in that: 1) establishes sequence stratigraphy framework: carrying out stratigraphic unit division to the target zone of target ancient karst oil gas field, sequence stratigraphy framework whens establishing etc.;2) it obtains ancient landform section and restores data: 1. obtaining superstratum thickness data;2. obtaining sedimentation thickness data;3. obtaining residual thickness data;3) ancient landform sectional view restores: 1. obtaining hatching and well location point;2. obtaining karst palaeotectonic setting line;3. obtaining sedimentation relative contour;4. obtaining stratum residual thickness line;5. completing ancient landform section according to above-mentioned steps to restore, ancient landform sectional view is generated, and indicates every graph element.The present invention is restored by ancient landform sectional view, is improved three-dimensional, restrictive, the quantitative, intuitive of ancient landform research, is enhanced geomorphic unit and portray precision.

Description

A kind of ancient landform section restoration methods
Technical field
The present invention relates to a kind of ancient landform section restoration methods, belong to exploration of oil and gas field field.
Background technique
Ancient karst oil gas field is a kind of important oil gas field type, widely distributed in China, such as: Bohai gulf basin Renqiu oil Field, Tahe Oilfield of The Tarim Basin take turns southern oil gas field and Ordos Basin is pacified the border region gas field etc..The formation of such oil gas field and ground Shell lifting is closely related with long-term weathering and erosion, rain drop erosion, chemolysis and eluviation, is distributed strictly by karst The control of ancient landform development characteristics.Meticulous depiction ancient landform unit and its regularity of distribution, to grasp karst reservoir space symbiosis And Hydrocarbon Enrichment Regularity has great practice significance.
Related palaeogeomorphic research reports more that end result is mainly landforms plan view at present, and geomorphological profile is discussed State less and rarely seen geomorphic profile.This not only causes geomorphologic map to lack constraint and verifying, and geomorphic unit portrays shortage directly The property seen, there are randomnesss for the criteria for classifying.
Summary of the invention
In view of the above-mentioned problems, the object of the present invention is to provide it is a kind of it is easy intuitive, the ancient landform section of strong operability is extensive Compound method.
To achieve the above object, the present invention takes following technical scheme: a kind of ancient landform section restoration methods comprising with Lower step:
1) sequence stratigraphy framework is established
To the target zone of target ancient karst oil gas field, stratigraphic unit division is carried out, sequence stratigraphy framework whens establishing etc.;
2) it obtains ancient landform section and restores data
1. obtaining superstratum thickness data: it is thick to obtain each well superstratum for the layering data implemented according to Strata Comparison Degree;
2. obtaining sedimentation thickness data: stratum is retained complete well by the layering data implemented according to Strata Comparison It filters out, determines sedimentation thickness before each well weathering karst;
3. obtaining residual thickness data: the layering data implemented according to Strata Comparison determines each well stratum residual thickness;
3) ancient landform sectional view restores
1. obtaining hatching and well location point: choosing geomorphological profile position, sequence reads section well underground coordinate, calculates between well Actual range obtains hatching;
2. obtaining karst palaeotectonic setting line: using actual range is abscissa between the section well well, karst palaeotectonic setting Line number obtains karst palaeotectonic setting line according to being ordinate;
3. obtaining sedimentation relative contour: using actual range is abscissa between the section well well, sedimentation thickness line number According to for ordinate, sedimentation relative contour is obtained;
4. obtaining stratum residual thickness line: using actual range is abscissa between the section well well, stratum residual thickness line number According to for ordinate, stratum residual thickness line is obtained;
5. completing ancient landform section according to above-mentioned steps to restore, ancient landform sectional view is generated, and indicates every graph element.
In the step 3), karst palaeotectonic setting line data calculation formula:
Y '=- Y+Ymax+Hmax
In formula, Y ' is karst palaeotectonic setting line number evidence;Y is superstratum thickness;YmaxFor superstratum thickness maximum Value;HmaxFor sedimentation maximum thickness.
In the step 3), sedimentation relative contour data calculation formula:
H '=Y '-H
In formula, Y ' is karst palaeotectonic setting line number evidence;H is sedimentation thickness;H ' is sedimentation relative contour data.
In the step 3), stratum residual thickness line data calculation formula:
H '=H '+h
In formula, H ' is sedimentation relative contour data;H is stratum residual thickness;H ' is stratum residual thickness line number evidence.
In the step 2), sedimentation thickness is obtained by trend analysis method before each well weathering karst.
Ymax+HmaxThe desirable integer more than or equal to sum of the two.
The invention adopts the above technical scheme, which has the following advantages: 1, the present invention is extensive by ancient landform sectional view It is multiple, three-dimensional, restrictive, the quantitative, intuitive of ancient landform research is improved, geomorphic unit is enhanced and portrays precision.2, originally Inventive method is directly used in the recovery of ancient landform section, and then is applied to palaeokarst-related reservoir genesis mechanism, rich accumulation of oil and gas and advises at hiding Numerous aspects such as rule, Favorable Areas screening and well site deployment.
Detailed description of the invention
Fig. 1 is flow chart of the invention;
Fig. 2 is the ancient landform spread figure of the embodiment of the present invention;
Fig. 3 is the weathering crust superstratum thickness chart of the embodiment of the present invention;
Fig. 4 is the weathering crust sedimentation thickness chart of the embodiment of the present invention;
Fig. 5 is the weathering crust stratum residual thickness figure of the embodiment of the present invention;
Fig. 6 is the recovery figure of the hatching of the embodiment of the present invention;
Fig. 7 is the recovery figure of the karst palaeotectonic setting line of the embodiment of the present invention;
Fig. 8 is the recovery figure of the sedimentation relative contour of the embodiment of the present invention;
Fig. 9 is the recovery figure of the stratum residual thickness line of the embodiment of the present invention;
Figure 10 is the ancient landform section recovery figure of the embodiment of the present invention.
Specific embodiment
The present invention is described in detail below with reference to the accompanying drawings and embodiments.
As shown in Figure 1, ancient landform section restoration methods of the invention the following steps are included:
1 establishes sequence stratigraphy framework
First to the target zone of target ancient karst oil gas field (i.e. weathering crust stratum), using conventional Stratigraphic Division And Correlation method Carry out stratigraphic unit division, establish etc. whens sequence stratigraphy framework;
2, which obtain ancient landform section, restores data
2.1 obtain superstratum thickness data: it is thick to obtain each well superstratum for the layering data implemented according to Strata Comparison Degree is according to Y;
2.2 acquisition sedimentation thickness datas: ancient karst will lead to missing and stratum exposure difference, i.e. stratum at the top of stratum Reserving degree is different.According to the layering data that Strata Comparison is implemented, stratum is retained into complete well screen and is elected, passes through trend point Analysis method determines sedimentation thickness H before each well weathering karst;
2.3 obtain residual thickness data: the layering data implemented according to Strata Comparison determines each well stratum residual thickness h;
Wherein, the target zone of ancient karst oil gas field is commonly referred to as weathering crust stratum in geology.
3 ancient landform sectional views restore
3.1 obtain hatching and well location point: choosing geomorphological profile position, sequence reads section well underground coordinate, calculates well Between actual range, obtain hatching;
3.2 obtain karst palaeotectonic setting line: due to the superstratum and mesh stratum stratum of ancient karst oil gas field target zone Between in filling up relationship, therefore impression method is generallyd use, according to weathering crust superstratum thickness to karst palaeotectonic setting Restored.General superstratum thickness is small, then is in karst paleostructure high-order bit;Conversely, being then in karst paleostructure lower curtate Position.
Using actual range between the section well well as abscissa, karst palaeotectonic setting line number is ordinate according to Y ', obtains rock Molten palaeotectonic setting line, karst palaeotectonic setting line number is according to Y ' calculation formula:
Y '=- Y+Ymax+Hmax
In formula, Y ' is karst palaeotectonic setting line number evidence;Y is superstratum thickness;YmaxFor superstratum thickness maximum Value;HmaxFor sedimentation maximum thickness;
Wherein, it is contemplated that mapping is beautiful rationally, Ymax+HmaxThe desirable integer more than or equal to sum of the two.
3.3 obtain sedimentation relative contour: using actual range is abscissa between the section well well, sedimentation thickness line number It is ordinate according to H ', obtains sedimentation relative contour, sedimentation relative contour data H ' calculation formula:
H '=Y '-H
In formula, Y ' is karst palaeotectonic setting line number evidence;H is sedimentation thickness;H ' is sedimentation relative contour data;
3.4 obtain stratum residual thickness line: using actual range is abscissa between the section well well, stratum residual thickness line number It is ordinate according to h ', obtains stratum residual thickness line, stratum residual thickness line number is according to h ' calculation formula:
H '=H '+h
In formula, H ' is sedimentation relative contour data;H is stratum residual thickness;H ' is stratum residual thickness line number evidence;
3.5, which complete ancient landform section according to above-mentioned steps, restores, and generates ancient landform sectional view, and indicates every graph and want Element.
Below by a specific embodiment, technical effect of the invention is further illustrated.
Embodiment
As shown in Fig. 2, gas field of pacifying the border region is a typical ancient rock using 330 wellblock of gas field east Shan of pacifying the border region as research object Molten gas field, in caledonian movement advanced stage, the gas field is integrally lifted, and makes to expose under earth's surface after the regression of Eopaleozoic Middle Ordovician Paleozoic group carbonate rock is subject to weathering and erosion, rain drop erosion and the chemolysis of 130Ma, eluviation, forms Ma Jiagou Five sections of top crust reservoirs of weathering of group.
1) sequence stratigraphy framework is established
Using conventional Stratigraphic Division And Correlation method carry out the division of weathering crust stratigraphic unit, establish etc. whens sequence stratigraphy framework;
2) ancient landform section restores data preparation
It obtains superstratum thickness data: data being layered according to the 330 wellblock drilling well of east Shan, it is thick to obtain each well superstratum Degree generates superstratum thickness chart (as shown in Figure 3) according to Y.Wherein, local area superstratum is the Carboniferous System;
It obtains sedimentation thickness data: data being layered according to the 330 wellblock drilling well of east Shan, stratum is retained into complete well It screens, by trend analysis, determines sedimentation thickness H before each well weathering karst, generate sedimentation thickness chart (as shown in Figure 4);
It obtains residual thickness data: data being layered according to the 330 wellblock drilling well of east Shan, determines each well stratum residual thickness h, It generates stratum residual thickness figure (as shown in Figure 5);
3) ancient landform sectional view restores
It obtains hatching and well location point: choosing in Fig. 2, B-B1 profile position, sequence reads section well underground coordinate, meter Actual range (as shown in table 1) between calculation well, generates hatching (as shown in Figure 6);
Actual range between 1 NS2 section well of table
Obtain karst palaeotectonic setting line: using actual range between well as abscissa, karst palaeotectonic setting line number is according to Y ' Ordinate obtains karst palaeotectonic setting line (as shown in Figure 7);
Due to Ymax+Hmax=70.88+60.2=131.08 (m), it is contemplated that mapping is beautiful rationally, take C=200 (> 131.08);
Obtain sedimentation relative contour: using actual range between well as abscissa, sedimentation relative contour data H ' is vertical sits Mark obtains sedimentation relative contour (as shown in Figure 8);
Obtain stratum residual thickness line: using actual range between well as abscissa, stratum residual thickness line number is vertical sit according to h ' Mark obtains stratum residual thickness line (as shown in Figure 9);
The ancient landform section for completing B-B1 section restores, and indicates every graph element (as shown in Figure 10), specific data As shown in table 2.
2 ancient landform sectional view of table restores data
Above description is only embodiments herein and the explanation to institute's application technology principle.Those skilled in the art answer Work as understanding, invention scope involved in the application, however it is not limited to technical solution made of the specific combination of above-mentioned technical characteristic, It should also cover in the case where not departing from the inventive concept simultaneously, carry out any group by above-mentioned technical characteristic or its equivalent feature The other technical solutions closed and formed.Such as features described above and (but being not limited to) disclosed herein have the skill of similar functions Art feature is replaced mutually and the technical solution that is formed.

Claims (6)

1. a kind of ancient landform section restoration methods comprising following steps:
1) sequence stratigraphy framework is established
To the weathering crust stratum of target ancient karst oil gas field, stratigraphic unit division is carried out, sequence stratigraphy framework whens establishing etc.;
2) it obtains ancient landform section and restores data
1. obtaining superstratum thickness data: the layering data implemented according to Strata Comparison obtains each well superstratum thickness;
2. obtaining sedimentation thickness data: stratum is retained complete well screen and selected by the layering data implemented according to Strata Comparison Out, sedimentation thickness before each well weathering karst is determined;
3. obtaining residual thickness data: the layering data implemented according to Strata Comparison determines each well stratum residual thickness;
3) ancient landform sectional view restores
1. acquisition hatching and well location point: choosing geomorphological profile position, sequence reads section well underground coordinate, calculates reality between well Distance obtains hatching;
2. obtaining karst palaeotectonic setting line: using actual range is abscissa between the section well well, karst palaeotectonic setting line number According to for ordinate, karst palaeotectonic setting line is obtained;
3. obtaining sedimentation relative contour: using actual range between the section well well as abscissa, sedimentation relative contour data are Ordinate obtains sedimentation relative contour;
4. obtaining stratum residual thickness line: using actual range between the section well well as abscissa, residual thickness line number evidence in stratum is Ordinate obtains stratum residual thickness line;
5. completing ancient landform section according to above-mentioned steps to restore, ancient landform sectional view is generated, and indicates every graph element.
2. a kind of ancient landform section restoration methods as described in claim 1, it is characterised in that: in the step 3), karst is ancient Tectonic setting line data calculation formula:
Y '=- Y+Ymax+Hmax
In formula, Y ' is karst palaeotectonic setting line number evidence;Y is superstratum thickness;YmaxFor superstratum maximum thickness;Hmax For sedimentation maximum thickness.
3. a kind of ancient landform section restoration methods as described in claim 1, it is characterised in that: in the step 3), stratum is heavy Product relative contour data calculation formula:
H '=Y '-H
In formula, Y ' is karst palaeotectonic setting line number evidence;H is sedimentation thickness;H ' is sedimentation relative contour data.
4. a kind of ancient landform section restoration methods as described in claim 1, it is characterised in that: in the step 3), stratum is residual Stay relative contour data calculation formula:
H '=H '+h
In formula, H ' is sedimentation relative contour data;H is stratum residual thickness;H ' is stratum residual thickness line number evidence.
5. a kind of ancient landform section restoration methods as described in claim 1, it is characterised in that: in the step 2), each well wind Sedimentation thickness is obtained by trend analysis method before changing karst.
6. a kind of ancient landform section restoration methods as claimed in claim 2, it is characterised in that: Ymax+HmaxIt is desirable to be greater than or wait In the integer of sum of the two.
CN201610121264.1A 2016-03-03 2016-03-03 A kind of ancient landform section restoration methods Expired - Fee Related CN105740574B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610121264.1A CN105740574B (en) 2016-03-03 2016-03-03 A kind of ancient landform section restoration methods

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610121264.1A CN105740574B (en) 2016-03-03 2016-03-03 A kind of ancient landform section restoration methods

Publications (2)

Publication Number Publication Date
CN105740574A CN105740574A (en) 2016-07-06
CN105740574B true CN105740574B (en) 2019-02-01

Family

ID=56249080

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610121264.1A Expired - Fee Related CN105740574B (en) 2016-03-03 2016-03-03 A kind of ancient landform section restoration methods

Country Status (1)

Country Link
CN (1) CN105740574B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107015290B (en) * 2017-03-13 2018-05-08 西北大学 A kind of method that reworked garden basin primary deposit looks are recovered
CN108734779B (en) * 2017-04-17 2021-07-30 中国石油化工股份有限公司 Ancient landform restoration method
CN107167853A (en) * 2017-06-01 2017-09-15 陕西延长石油(集团)有限责任公司研究院 A kind of fine method for recovering Palaeokarst Landform
CN107589470B (en) * 2017-08-24 2019-05-03 西南石油大学 Palaeokarst Landform restoration methods in small area
CN110794477B (en) * 2018-08-03 2022-03-01 中国石油天然气股份有限公司 Method for identifying overtaking stratum pinch-out line by combining well and seismic
CN110361781B (en) * 2019-08-14 2020-10-09 西南石油大学 Sedimentary-phase ancient landform restoration method based on seismic reflection amplitude analysis
CN111983678B (en) * 2020-07-22 2024-02-23 中海石油(中国)有限公司深圳分公司 Method for rapidly evaluating development potential of deep water sand body
CN112130210B (en) * 2020-09-03 2024-07-26 中国石油天然气集团有限公司 Well position optimization method and device based on ancient and modern geological elements
CN113009591B (en) * 2020-11-03 2023-04-07 中国石油天然气股份有限公司 Ancient landform restoration method
CN113312695B (en) * 2021-05-26 2022-03-11 成都理工大学 Ancient landform restoration method based on trend surface analysis
CN113640870B (en) * 2021-07-27 2024-02-23 中国石油天然气股份有限公司西南油气田分公司勘探开发研究院 Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface
CN114415233B (en) * 2022-01-20 2024-06-25 北京师范大学 Method for recovering morphology of stratum of mined fossil energy

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103424773A (en) * 2012-05-25 2013-12-04 中国石油化工股份有限公司 S Stratification position flattening method-based ancient landform restoration method
CN103853882A (en) * 2014-02-13 2014-06-11 中铁工程设计咨询集团有限公司 Drawing method and system for engineering geological profile
CN103901471A (en) * 2014-03-31 2014-07-02 成都理工大学 Method for recovering depositional sequence prototype section

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103424773A (en) * 2012-05-25 2013-12-04 中国石油化工股份有限公司 S Stratification position flattening method-based ancient landform restoration method
CN103853882A (en) * 2014-02-13 2014-06-11 中铁工程设计咨询集团有限公司 Drawing method and system for engineering geological profile
CN103901471A (en) * 2014-03-31 2014-07-02 成都理工大学 Method for recovering depositional sequence prototype section

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"鄂尔多斯盆地中部气田奥陶系古地貌研究";代金友等;《石油学报》;20050531;第20卷(第3期);第37-39,43页
"鄂尔多斯盆地延长探区奥陶纪末古地貌恢复与储层预测";韩敏强等;《石油与天然气地质》;20111031;第32卷(第54期);第760-767页

Also Published As

Publication number Publication date
CN105740574A (en) 2016-07-06

Similar Documents

Publication Publication Date Title
CN105740574B (en) A kind of ancient landform section restoration methods
Sharp et al. Stratigraphic architecture and fracture-controlled dolomitization of the Cretaceous Khami and Bangestan groups: an outcrop case study, Zagros Mountains, Iran
CN104632200B (en) Fluvial facies multilayer and thick-layer oil reservoir fine modeling method
CN106875471B (en) Three-dimensional visual modeling method for coal-series water-containing or water-resisting layer
CN103899285B (en) Multilamellar sandstone oil reservoir near limit Water-cut Period is rotated water drive method
CN108828666B (en) Method for recovering buried history of stratum near normal fault
CN105425316A (en) Ancient landform restoration method based on impression method
Lindenmaier et al. Structure and genesis of the Cubango Megafan in northern Namibia: implications for its hydrogeology
Frost III et al. Syndepositional deformation in a high-relief carbonate platform and its effect on early fluid flow as revealed by dolomite patterns
CN109783944A (en) Three-dimensional Flow model of element method for building up, device and equipment
Sonntag et al. Sedimentological controls on the fracture distribution and network development in Mesaverde Group sandstone lithofacies, Uinta Basin, Utah, USA
Tinker et al. Reservoir characterization of a Permian giant: Yates field, west Texas
CN106338761A (en) Three-dimensional embedded entity modeling method
CN105719339A (en) Shale gas reservoir lamellation fracture three-dimensional modeling method
CN104574513A (en) Representation method for accurately depicting three-dimensional distribution of interlayers based on geometrical morphology
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
Wilson et al. Regional geologic characterization of the Grayburg-San Andres reservoir for salt water disposal management, Midland Basin, Texas
Brantferger et al. Impact and lessons of using high permeability streaks in history matching a giant offshore Middle East carbonate reservoir
Ciftci et al. Outcrop-based three-dimensional modeling of the Tensleep sandstone at Alkali Creek, Bighorn Basin, Wyoming
CN109188525A (en) A kind of acquisition methods and system of marine facies rammell buried depth data
Faulds et al. Cenozoic evolution of the abrupt Colorado Plateau–Basin and Range boundary, northwest Arizona: A tale of three basins, immense lacustrine-evaporite deposits, and the nascent Colorado River
Raynaud et al. Geological and topographical study of the original hills at the base of Fourth Dynasty Egyptian monuments of the Memphite plateau
La Fontaine et al. Integrated outcrop and subsurface geomodeling of the Turonian Wall Creek Member of the Frontier Formation, Powder River Basin, Wyoming, USA
Le et al. Integrated Stratigraphic Controls for Flow Simulation of the Wall Creek Member in the Frontier Formation: Western Powder River Basin, Wyoming
French et al. Accommodation-controlled Systems-tractspecific Facies Partitioning and Resulting Geometric Development of Reservoir Grainstone Ramp-crest Shoal Bodies

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: 20190201

Termination date: 20210303

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