CN109799531A - A kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties - Google Patents

A kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties Download PDF

Info

Publication number
CN109799531A
CN109799531A CN201910211333.1A CN201910211333A CN109799531A CN 109799531 A CN109799531 A CN 109799531A CN 201910211333 A CN201910211333 A CN 201910211333A CN 109799531 A CN109799531 A CN 109799531A
Authority
CN
China
Prior art keywords
frequency
grid
coherence properties
vertex
frequency dividing
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
CN201910211333.1A
Other languages
Chinese (zh)
Other versions
CN109799531B (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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201910211333.1A priority Critical patent/CN109799531B/en
Publication of CN109799531A publication Critical patent/CN109799531A/en
Application granted granted Critical
Publication of CN109799531B publication Critical patent/CN109799531B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

A kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties, is under the jurisdiction of oil gas geophysics reservoir prediction field, it is desirable to provide a kind of fracture reservoir analytical technology for accurately identifying different scale crack and prediction distribution rule.The following steps are included: 1. extracting the frequency dividing coherence properties of purpose of horizon different frequency, and carry out enhancing processing;2. dividing coherence properties data to gridding treated purpose of horizon, the average value and variance on 4 vertex of each grid are calculated;3. the range of statistical variance value takes the intermediate value of the range, the average value of the corresponding all grids of the value is counted, the threshold values for dividing seam unit and non-seam unit is set according to statistical result, quickly and effectively obtains the distribution situation in crack;4. to being the grid for stitching unit in two different frequency purpose of horizon frequency dividing coherence properties data, vertex average value away from nearest non-seam unit is assigned to 4 vertex of the grid in upper frequency, the part being overlapped in the results of fracture prediction of upper frequency with the results of fracture prediction of lower frequency can be effectively removed;5. utilizing the frequency dividing coherence properties data through above-mentioned steps treated purpose of horizon different frequency, by plane at figure, fracture reservoir prediction is realized.The present invention is based on purpose of horizon to divide coherence properties data, divides seam unit and non-seam unit, compares the distribution situation for stitching unit in different frequency data, remove intersection, complete to different scale crack accurately identify and the Accurate Prediction of the regularity of distribution.

Description

A kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties
Technical field
It is a kind of based on poststack 3-D seismics the present invention relates to the fracture reservoir forecasting problem in oil gas geophysics field The fracture reservoir prediction technique of attributive analysis, for provide it is a kind of accurately identify different scale crack and prediction distribution rule split Stitch reservoir prediction analytical technology.
Background technique
Poststack 3-D seismics attributive analysis is the important method of crack identification and prediction, detects crack master using seismic properties There are coherent body, curvature, inclination angle etc. according to the discontinuity of seismic waveform.Since crack has multiple dimensioned property, how to difference Crack carries out accurate characterization, the prediction to the identification in different scale crack, the regularity of distribution, is the emphasis and difficulty of research all the time Point.Based on many scholars at home and abroad propose some technology and methods, such as Zhang Guangzhi (2011) is converted by Curvelet, Different reconstruction coefficients are given in bent wave zone, the seismic data cube of prominent different frequency bands and different direction are obtained, in conjunction with phase Stem body edge detecting technology predicts fracture and Fractured Zone and its trend;Guo Rui etc. (2014) proposes that one kind is based on The crack analysis method of high frequency recovery and multi-resolution decomposition, widens the frequency band of three dimensional seismic data, and seismic channel is decomposed into correspondence The data volume of different scale (corresponding with reconstruct wavelet dominant frequency), analyzes crack using relevant and curvature attributes, portrays different rulers Spend lower formation fracture development behavior and distribution characteristics;Wu Peng etc. (2011) is carried out high-precision using the coherence properties technology that small echo divides Spend the explanation in crack, geology body characteristics.Mainly rear 3D seismic data carries out multi-resolution decomposition to above method through the stack, then In conjunction with poststack 3-D seismics attribute analysis technology relevant to crack, the crack of different scale and prediction distribution rule are identified.But Above method all has a problem that: the poststack 3D seismic data based on different decomposition scale carries out FRACTURE PREDICTION, obtains As a result there is intersection, cause true to the identification in different scale crack and the forecasting inaccuracy of the regularity of distribution.Base of the present invention as a result, In poststack 3-D seismics attribute analysis technology, it is contemplated that FRACTURE PREDICTION knot of the crack in the seismic data based on different decomposition scale Distribution characteristics having the same in fruit, the present invention will based on this, accurately identify crack and the prediction distribution rule of different scale Rule.
Summary of the invention
The present invention is intended to provide a kind of fracture reservoir forecast analysis technology based on earthquake frequency dividing coherence properties, effective to solve Certainly different scale crack is accurately identified and the problem of prediction distribution rule.Using spectra decomposition to seismic data into Row processing is analyzed and extracts frequency dividing attribute;The frequency dividing coherence properties of different frequency are obtained in conjunction with seismic coherence attribute analysis technology Body;Using the seismic horizon explained, the frequency dividing coherence properties of purpose of horizon different frequency are extracted;To the frequency dividing phase of purpose of horizon Dry attribute carries out gridding processing, and statisticallys analyze the vertex average value and variance of all grids, is drawn according to statistic analysis result Divide grid property;The distribution situation that unit grid is stitched in omparison purpose layer position different frequency data, finds out intersection, i.e. certain net Lattice are seam unit in the data of different frequency, are handled this kind of grid in higher frequency data;After processing Purpose of horizon different frequency frequency dividing coherence properties, carry out fracture reservoir prediction.
Specific steps of the invention include:
(1) time frequency analysis based on continuous wavelet transform is carried out to poststack 3D seismic data first, analysis is found out effectively Seismic band range generates a series of single-frequency data volume of different frequencies using frequency filtering technology in the range.
(2) on obtained single-frequency data volume, coherence properties extraction is carried out using third generation coherent algorithm, obtains different frequencies The frequency dividing coherence properties body of rate;
(3) enhancing processing is carried out to the frequency dividing coherence properties body of different frequency;
(4) interpreted good seismic horizon is utilized, for specific on the frequency dividing coherence properties body handled through step (3) Purpose of horizon, extract the frequency dividing coherence properties of corresponding layer position;
(5) gridding processing is carried out to the frequency dividing coherence properties data of purpose of horizon, calculates the flat of 4 vertex of each grid Mean value and variance;The variance on the vertex of all grids is counted, its distribution is obtained;It takes in variance distribution Value finds out the grid that all variances are the value, counts the vertex average value of these grids, sets one further according to the result of statistics Threshold value to stitch unit, otherwise is non-seam unit when the vertex average value of grid is greater than threshold value;
(6) distribution for comparing the seam unit and non-seam unit of two different frequency purpose of horizon frequency dividing coherence properties data, when Grid is seam unit in the result of two frequencies, then carries out assignment processing to 4 vertex of the grid in upper frequency, tax Value is the non-seam unit vertex average value nearest apart from the grid.
(7) FRACTURE PREDICTION is carried out using the frequency dividing coherence properties through treated purpose of horizon different frequency.
The present invention is a kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties, has following features, mainly It shows themselves in that
(1) coherence properties data are divided to gridding treated purpose of horizon earthquake, calculates 4 tops of each grid The average value and variance, the range of statistical variance value of point, take the intermediate value of the range, count being averaged for the corresponding all grids of the value Value sets the threshold values for dividing seam unit and non-seam unit according to statistical result, can quickly and effectively obtain the distribution feelings in crack Condition;
(2) to being the grid for stitching unit in two different frequency purpose of horizon frequency dividing coherence properties data, to higher-frequency 4 vertex of the grid carry out assignment in rate, and assigned value is the vertex average value of the non-seam unit nearest apart from the grid, can It is divided with effectively removing in the results of fracture prediction based on upper frequency earthquake frequency dividing coherence properties with based on lower frequency earthquake The forecasting inaccuracy of the part that the results of fracture prediction of coherence properties is overlapped, the identification and the regularity of distribution that solve different scale crack is true The problem of.
Specific embodiment
A kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties, specific implementation step are as follows:
(1) it extracts purpose of horizon and divides coherence properties
Step 1: small echo divides attributes extraction:
1. using small echo to carry out time frequency analysis to seismic data for the continuous wavelet transform of Morlet small echo, by seismic data Body transforms to frequency domain from time-domain;
2. calculating the frequency spectrum at each seismic channel all time points, the earthquake for generating frequency altogether is then rearranged according to frequency Data;
Step 2: the single-frequency seismic data for the different frequency that step 1 obtains is calculated using third generation coherent algorithm The coherence properties body of different frequency single-frequency seismic data cube;
Step 3: enhancing processing is carried out to the coherence properties body for the different frequency that step 2 is calculated;
Step 4: interpreted good seismic horizon is utilized, for spy on the frequency dividing coherence properties body handled through step 4 Fixed purpose of horizon extracts the frequency dividing coherence properties of corresponding layer position;
(2) FRACTURE PREDICTION is carried out based on purpose of horizon frequency dividing coherence properties
Step 1: the distribution characteristics in crack in analysis purpose of horizon different frequency frequency dividing coherence properties data:
1. the frequency dividing coherence properties data to purpose of horizon different frequency carry out gridding processing, and calculate each grid 4 The average value and variance on a vertex;
2. counting to the variance on the vertex of all grids, its distribution is obtained;The intermediate value of variance distribution is taken, The grid that all variances are the value is found out, the vertex average value of these grids is counted, sets a threshold further according to the result of statistics Value to stitch unit, otherwise is non-seam unit, obtains purpose of horizon different frequency when the vertex average value of grid is greater than threshold value Divide the distribution characteristics in crack in coherence properties data;
Step 2: point of the seam unit of the frequency dividing coherence properties data of omparison purpose two different frequency of layer position and non-seam unit Cloth then carries out at assignment 4 vertex of the grid in upper frequency when grid is seam unit in the result of two frequencies Reason, assigned value are the non-seam unit vertex average value nearest apart from the unit grid.Grid 1 and grid are calculated using following formula The distance between 2:
In formula: a, b are respectively the vertical survey line number and cross-track number of the bottom right vertex of grid 1, and i, j are respectively grid 2 The vertical survey line number and cross-track number of bottom right vertex, Δ x are vertical survey line spacing, and Δ y is cross-track spacing.
Step 3: utilize the frequency dividing coherence properties data through treated purpose of horizon different frequency, by plane at Figure, obtains the graphical effect of FRACTURE PREDICTION, to realize the prediction of fracture reservoir.

Claims (3)

1. a kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties, it is characterised in that using a kind of completely new side Method, such as following steps:
(1) using poststack 3D seismic data as research object, using small echo is the continuous wavelet transform of Morlet small echo to earthquake Data carry out time frequency analysis, and seismic data is transformed to frequency domain from time-domain, are generated using frequency filtering technology a series of The single-frequency data volume of different frequency obtains the earthquake frequency dividing coherence properties body of different frequency using third generation coherent algorithm;
(2) enhancing processing is carried out to the earthquake frequency dividing coherence properties body of different frequency, the interpreted good seismic horizon of utilization mentions Take the frequency dividing coherence properties of corresponding layer position;
(3) gridding processing is carried out to the frequency dividing coherence properties data of purpose of horizon, calculates the average value on each 4 vertex of grid With variance, the variance on the vertex of all grids is counted, obtains its distribution, the intermediate value of the range is taken, finds out variance For the grid of the value, the vertex average value of these grids is counted, a threshold value is set further according to the result of statistics, when the top of grid It to stitch unit, otherwise is non-seam unit when point average value is greater than the threshold value;
(4) distribution for comparing the seam unit and non-seam unit of two different frequency purpose of horizon frequency dividing coherence properties data, works as grid It is seam unit in the result of two frequencies, then assignment processing is carried out to 4 vertex of the grid in upper frequency, assigned value is The non-seam unit vertex average value nearest apart from the grid, using formula (1) calculate between grid 1 and 2 central point of grid away from From:
In formula: a, b are respectively the vertical survey line number and cross-track number of the bottom right vertex of grid 1, and i, j are respectively the bottom right of grid 2 The vertical survey line number and cross-track number on vertex, Δ x are vertical survey line spacing, and Δ y is cross-track spacing;
(5) FRACTURE PREDICTION is carried out using the frequency dividing coherence properties through treated purpose of horizon different frequency.
2. a kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties according to claim 1, feature exist In: coherence properties data are divided to gridding treated purpose of horizon earthquake, calculate being averaged for 4 vertex of each grid Value and variance, the range of statistical variance value take the intermediate value of the range, count the average value of the corresponding all grids of the value, according to Statistical result setting divides the threshold values of seam unit and non-seam unit, quickly and effectively obtains the distribution situation in crack.
3. a kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties according to claim 1, feature exist In: to being the grid for stitching unit in two different frequency purpose of horizon frequency dividing coherence properties data, to the net in upper frequency 4 vertex of lattice carry out assignment, and assigned value is the vertex average value of the non-seam unit nearest apart from the grid, can effectively go Except in the results of fracture prediction based on upper frequency earthquake frequency dividing coherence properties and based on lower frequency earthquake frequency dividing coherence properties Results of fracture prediction be overlapped part, calculate grid element center point between distance use formula 1.
CN201910211333.1A 2019-03-20 2019-03-20 Fracture reservoir prediction method based on seismic frequency division coherence attribute Expired - Fee Related CN109799531B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910211333.1A CN109799531B (en) 2019-03-20 2019-03-20 Fracture reservoir prediction method based on seismic frequency division coherence attribute

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910211333.1A CN109799531B (en) 2019-03-20 2019-03-20 Fracture reservoir prediction method based on seismic frequency division coherence attribute

Publications (2)

Publication Number Publication Date
CN109799531A true CN109799531A (en) 2019-05-24
CN109799531B CN109799531B (en) 2022-04-01

Family

ID=66563690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910211333.1A Expired - Fee Related CN109799531B (en) 2019-03-20 2019-03-20 Fracture reservoir prediction method based on seismic frequency division coherence attribute

Country Status (1)

Country Link
CN (1) CN109799531B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111239817A (en) * 2020-02-05 2020-06-05 成都理工大学 Method for improving fault likelihood attribute analysis precision
CN111399056A (en) * 2020-04-29 2020-07-10 西南石油大学 Method for predicting crack strength based on divided azimuth filtering
CN112415586A (en) * 2020-11-16 2021-02-26 北京孚梅森石油科技有限公司 Slurry leakage risk assessment method and device
CN112649863A (en) * 2019-10-12 2021-04-13 中国石油化工股份有限公司 Frequency division seismic attribute data optimization method and system
CN116520419A (en) * 2023-05-05 2023-08-01 中国矿业大学(北京) Hot fluid crack channel identification method

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0655861A1 (en) * 1993-11-26 1995-05-31 Koninklijke Philips Electronics N.V. Image composition method and imaging apparatus for performing said method
US6407404B1 (en) * 1999-03-15 2002-06-18 Denso Corporation Apparatus for the examining defect of monolithic substrate and method for examining the same
US7254091B1 (en) * 2006-06-08 2007-08-07 Bhp Billiton Innovation Pty Ltd. Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs
US20100259541A1 (en) * 2009-04-08 2010-10-14 Canon Kabushiki Kaisha Image processing apparatus and image processing method
US20100260412A1 (en) * 2009-04-13 2010-10-14 Canon Kabushiki Kaisha Image processing apparatus and control method thereof
CN102306376A (en) * 2009-11-03 2012-01-04 蒋慧琴 Method for adaptive medical image enhancement
CN103777244A (en) * 2012-10-25 2014-05-07 中国石油化工股份有限公司 Method for quantitatively analyzing earthquake crack attribute volume
CN104749617A (en) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 Multi-scale fractured reservoir forward model establishing method
US20170358143A1 (en) * 2016-06-09 2017-12-14 Adobe Systems Incorporated Controlling Patch Usage in Image Synthesis
US20170358128A1 (en) * 2016-06-09 2017-12-14 Jakub Fiser Illumination-Guided Example-Based Stylization of 3D Renderings
CN107909083A (en) * 2017-11-10 2018-04-13 中国科学院地理科学与资源研究所 A kind of hough transform extracting method based on outline optimization
CN108931814A (en) * 2017-05-26 2018-12-04 中国石油化工股份有限公司 A method of the basement rock FRACTURE PREDICTION based on the fusion of more attributes
CN109117851A (en) * 2018-07-06 2019-01-01 航天星图科技(北京)有限公司 A kind of video image matching process based on lattice statistical constraint
CN109389172A (en) * 2018-10-11 2019-02-26 中南大学 A kind of radio-signal data clustering method based on printenv grid

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0655861A1 (en) * 1993-11-26 1995-05-31 Koninklijke Philips Electronics N.V. Image composition method and imaging apparatus for performing said method
US6407404B1 (en) * 1999-03-15 2002-06-18 Denso Corporation Apparatus for the examining defect of monolithic substrate and method for examining the same
US7254091B1 (en) * 2006-06-08 2007-08-07 Bhp Billiton Innovation Pty Ltd. Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs
US20100259541A1 (en) * 2009-04-08 2010-10-14 Canon Kabushiki Kaisha Image processing apparatus and image processing method
US20100260412A1 (en) * 2009-04-13 2010-10-14 Canon Kabushiki Kaisha Image processing apparatus and control method thereof
CN102306376A (en) * 2009-11-03 2012-01-04 蒋慧琴 Method for adaptive medical image enhancement
CN103777244A (en) * 2012-10-25 2014-05-07 中国石油化工股份有限公司 Method for quantitatively analyzing earthquake crack attribute volume
CN104749617A (en) * 2013-12-26 2015-07-01 中国石油化工股份有限公司 Multi-scale fractured reservoir forward model establishing method
US20170358143A1 (en) * 2016-06-09 2017-12-14 Adobe Systems Incorporated Controlling Patch Usage in Image Synthesis
US20170358128A1 (en) * 2016-06-09 2017-12-14 Jakub Fiser Illumination-Guided Example-Based Stylization of 3D Renderings
CN108931814A (en) * 2017-05-26 2018-12-04 中国石油化工股份有限公司 A method of the basement rock FRACTURE PREDICTION based on the fusion of more attributes
CN107909083A (en) * 2017-11-10 2018-04-13 中国科学院地理科学与资源研究所 A kind of hough transform extracting method based on outline optimization
CN109117851A (en) * 2018-07-06 2019-01-01 航天星图科技(北京)有限公司 A kind of video image matching process based on lattice statistical constraint
CN109389172A (en) * 2018-10-11 2019-02-26 中南大学 A kind of radio-signal data clustering method based on printenv grid

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
于敏捷: "分频地震属性在裂缝储层预测中的应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 *
伍 鹏, 陈学华等: "基于小波分频的地震资料相干属性对比分析", 《长江大学学报(自然科学版)》 *
陈波,孙德胜,朱筱敏等: "利用地震数据分频相干技术检测火山岩裂缝_陈波", 《石油地球物理勘探》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112649863A (en) * 2019-10-12 2021-04-13 中国石油化工股份有限公司 Frequency division seismic attribute data optimization method and system
CN111239817A (en) * 2020-02-05 2020-06-05 成都理工大学 Method for improving fault likelihood attribute analysis precision
CN111399056A (en) * 2020-04-29 2020-07-10 西南石油大学 Method for predicting crack strength based on divided azimuth filtering
CN111399056B (en) * 2020-04-29 2020-12-08 西南石油大学 Method for predicting crack strength based on divided azimuth filtering
CN112415586A (en) * 2020-11-16 2021-02-26 北京孚梅森石油科技有限公司 Slurry leakage risk assessment method and device
CN112415586B (en) * 2020-11-16 2023-10-03 北京孚梅森石油科技有限公司 Slurry leakage risk assessment method and device
CN116520419A (en) * 2023-05-05 2023-08-01 中国矿业大学(北京) Hot fluid crack channel identification method
CN116520419B (en) * 2023-05-05 2024-01-30 中国矿业大学(北京) Hot fluid crack channel identification method

Also Published As

Publication number Publication date
CN109799531B (en) 2022-04-01

Similar Documents

Publication Publication Date Title
CN109799531A (en) A kind of fracture reservoir prediction technique based on earthquake frequency dividing coherence properties
CN107976713B (en) A kind of method and device of the lower removal sedimentation setting of higher-dimension seismic data input
CN105510964A (en) Seismic recognition method of low-order strike-slip faults in complex structural areas
CN104267431B (en) The craven fault acquisition methods that a kind of seismotectonics is explained
CN103728659B (en) A kind of method improving detection of karst cave precision
CN105334534A (en) Low order fault interpretation method based on construction mode guidance
WO2010096759A1 (en) Dynamic extraction of extrema-based geometric primitives in 3d seismic voxel volume
CN103576191B (en) A kind of method adopting seismic properties identification tomography
CN104020492A (en) Edge-preserving filtering method of three-dimensional earthquake data
CN105866832B (en) The recognition methods for the low-grade fault that Deep Tight Sandstone is hidden and device
CN102831588B (en) De-noising processing method for three-dimensional seismic images
CN103364835A (en) Stratum structure self-adaption median filtering method
CN105445801B (en) A kind of processing method for eliminating 2-d seismic data random noise
CN105259579A (en) A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes
CN104181591A (en) Multi-dimension seismic attribute fusion method and device
CN105445787B (en) A kind of relevant crack prediction method of top optimization direction daughter
CN107942405A (en) The method for predicting thin sand-mud interbed sand body cumulative thickness
CN108226997A (en) A kind of seismic facies analysis method based on earthquake data before superposition
CN103592684A (en) Massive seismic data compression method and device for preserving spatial attribute information
CN1325935C (en) Superiority frequency band coherent handling method used for fine tomography explanation
Cano-Domingo et al. Segmentation and characteristic extraction for Schumann Resonance transient events
CN110954958A (en) Crack and fault prediction method and system
CN108680950A (en) A kind of desert seismic signal method for detecting position based on Self-adaptive Block Matching
Bravo et al. Volume curvature attributes to identify subtle faults and fractures in carbonate reservoirs: Cimarrona Formation, Middle Magdalena Valley Basin, Colombia.
CN109061726B (en) Method and device for identifying multiples

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