CN106873032B - A kind of granite hill reservoir crack prediction method - Google Patents
A kind of granite hill reservoir crack prediction method Download PDFInfo
- Publication number
- CN106873032B CN106873032B CN201710083193.5A CN201710083193A CN106873032B CN 106873032 B CN106873032 B CN 106873032B CN 201710083193 A CN201710083193 A CN 201710083193A CN 106873032 B CN106873032 B CN 106873032B
- Authority
- CN
- China
- Prior art keywords
- ditch
- ridge
- line
- buried hill
- top surface
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000010438 granite Substances 0.000 title claims abstract description 12
- 206010017076 Fracture Diseases 0.000 claims abstract description 19
- 208000010392 Bone Fractures Diseases 0.000 claims abstract description 18
- 238000011161 development Methods 0.000 claims abstract description 9
- 230000007704 transition Effects 0.000 claims description 6
- 230000007306 turnover Effects 0.000 claims description 3
- 239000013049 sediment Substances 0.000 claims description 2
- 239000000284 extract Substances 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 3
- 238000011160 research Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
The present invention relates to a kind of granite hill reservoir crack prediction methods, which is characterized in that includes the following steps：1) paleogeographic reconstruction is carried out to target area buried hill top surface；2) to the buried hill top surface after evening up for the last time, i.e. buried hill ancient landform, interpolation generates ancient landform contour map；3) buried hill ancient landform " ditch " and " ridge " are identified；4) " ditch " and " ridge " distributed areas are portrayed and is extracted；5) fracture development area is predicted.The present invention is directed to Granite Reservoir FRACTURE PREDICTION, it identifies and extracts using " ditch " and " ridge " that buried hill ancient landform is relatively high, relatively low is represented, come Prediction of fracture area, and Vector triangle is used in the process and identifies ancient landform " ditch " and " ridge " one by one, the prediction regulation and standardization for making " ditch " and " ridge " makes the more traditional Geological working method of FRACTURE PREDICTION workload be greatly reduced.
Description
Technical field
The present invention relates to a kind of crack prediction methods, especially a kind of to be directed to granite hill reservoir crack prediction method.
Background technology
The country has the history of decades about the research of fracture-type reservoir, is basically limited to before the 1990s
Traditional geological method, only from the recognition methods in the source investigation crack of wellbore itself；And in the world, main oil company encloses
Or carrying out a large amount of tackling key problem research around crack.Recent years, petrochina exploration and development research institute constitute a packet
The multidisciplinary crack research project group for including geology, well logging, earthquake and reservoir engineering expert, by summarizing domestic and international more than ten cracks
Property oil reservoir exploration and development practice, explore the crack identification and forecast of distribution geologic assessment technology of a set of synthesis.Generally speaking,
Crack characterization Research Thinking is gradually intended to the synthesis of multi-specialized, multiple dimensioned, more technological means, that is, outcrop is applied to study,
Crack observation, the description of rock core crack, identifying fracture in well logging, Seismic Fracture identification, Rock Mechanics Test, stress field research etc. are a variety of
Method, in conjunction with a variety of professions such as earthquake, well logging, geology, oil reservoir, use is qualitative, Semi-qualitative, the multiple means such as quantitative, integrates macro
A variety of scales such as sight, microcosmic take the number of ways such as numerical simulation, physical analogy, comprehensive to determine that description and prediction Fractured are stored up
Layer spread.But since fractue spacing is extremely complex, research difficulty it is big, description and prediction be petroleum industrial circle difficult point, heat
Point.Wherein, the objective description of intrinsic fracture development degree is the basis of crack description and prediction.
Invention content
In view of the above-mentioned problems, the object of the present invention is to provide a kind of granite hill reservoir crack prediction methods.
To achieve the above object, the present invention takes following technical scheme：A kind of granite hill reservoir crack prediction method,
It includes the following steps：
1) buried hill top surface paleogeographic reconstruction
Paleogeographic reconstruction is carried out to target area buried hill top surface using tradition stratum flattening technique：
Target area buried hill top surface is successively evened up from top to bottom, gradually reflects buried hill top surface landforms by newly to old change
Change；The interfaces at different levels of flattening are each stratum in buried hill overlying sediments stratum " group ", " section " top surface, and lowest part interface (represents most
Buried hill top surface after evening up in early days) reflects landform shape of the buried hill before receiving to deposit, and this change procedure reflects respectively
Period buried hill top surface the changing of physiognomy situation, that is, reflect the historical process of Landform Evolution.
2) ancient landform contour in buried hill is made picture
To the buried hill top surface after evening up for the last time, i.e. buried hill ancient landform, interpolation generates ancient landform contour map.
3) buried hill ancient landform " ditch " and " ridge " identification
Relatively a height of " ridge " on the ancient landform of buried hill, relatively low is " ditch ".Relatively high, relatively low judgement uses triangle
Shape rule：(1) adjacent three points of connection landform turnover, form triangle；(2) with minimum with horizontal line angle in triangle
On the basis of side, when the 3rd vertex is located on the side, the position of third point represents palaeogeomorphic " ridge "；(3) when the 3rd
When a vertex is located under the side, the position of third point represents palaeogeomorphic " ditch ".It, can be by using above-mentioned Vector triangle
One identifies buried hill palaeogeomorphic " ditch " and " ridge ".
4) " ditch " and " ridge " distributed areas portraying and extracting
Buried hill ancient landform is a rugged curved surface, the step of portraying and extract on curved surface " ditch " and " ridge " distributed areas
It is as follows：
A. main " ditch " and " ridge " side of trend in working region (region in crack to be predicted, hereinafter referred to as working region) is determined
To；
B. since working region is on one side, perpendicular to main " ditch " and " ridge " trend, the 1st arbitrary line is drawn；
C. each terrain transition point on the arbitrary line is determined；
D. judge that each terrain transition point is " ditch " or " ridge " on the arbitrary line using Vector triangle；
E. it is parallel to the 1st arbitrary line, the 2nd arbitrary line is drawn with 50m distances, step a~d is repeated and determines that the 2nd is appointed
" ditch " and " ridge " anticipated on line.The 3rd is drawn successively, the 4th arbitrary line ... determine respectively until working region another side
" ditch " and " ridge " on each arbitrary line；
F. " ditch " and " ridge " on above-mentioned arbitrary line is added on the paleogeomorphological map of buried hill top surface；
G. using buried hill top surface paleogeomorphological map as auxiliary, " ditch " and " ridge " on the arbitrary line of each item is sequentially connected, is formed
" ditch " and " ridge " line distribution map；
H. " ditch " and " ridge " distributed areas are determined：Centered on " ditch " line, both sides are with the center line of " ditch " and " ridge " line
Boundary limits ranging from buried hill top surface ancient landform " ditch " distributed areas；Centered on " ridge " line, both sides are with " ridge " and " ditch " line
Center line is boundary, limits ranging from buried hill top surface ancient landform " ridge " distributed areas.
5) fracture developing zone is predicted
Buried hill top surface landforms " ditch " distributed areas development hundreds of microns more than large fracture, quantity be generally 2~8/
m；Develop hundreds of microns small crack below, generally tens of~hundreds of/m of quantity, " ridge " distributed areas in " ridge " distributed areas
As advantageous fracture reservoir development area.
The invention adopts the above technical scheme, which has the following advantages：1, the present invention is directed to Granite Reservoir crack
Prediction identifies and extracts using " ditch " and " ridge " that buried hill ancient landform is relatively high, relatively low is represented, comes Prediction of fracture area,
And Vector triangle is used in the process and identifies ancient landform " ditch " and " ridge " one by one, makes the prediction specification of " ditch " and " ridge "
Change, standardization.2, the present invention by way of arbitrary line, realize in working region " ditch " and " ridge " distribution portray with
Extraction, realizes the prediction of fracture developing zone, has stronger operability, to user almost without skill requirement.3, of the invention
Due to using procedure, standardized working method, the more traditional Geological working method of FRACTURE PREDICTION workload being made to be greatly reduced.
Description of the drawings
Fig. 1 is the Reservoir Fractures result figure described based on paleostructure ditch ridge in the embodiment of the present invention.
In figure："+" indicates that paleostructure high point, "/" indicate that arbitrary line, " " indicate that triangle fixed point (ditch), " " indicate triangle
Shape vertex (ridge),Indicate paleostructure contour,Indicate paleostructure ditch line,Indicate paleostructure crestal line,Indicate paleostructure ditch region,Indicate paleostructure ridge region.
Specific implementation mode
The present invention is described in detail below with reference to the accompanying drawings and embodiments.It should be appreciated, however, that the offer of attached drawing is only
For a better understanding of the present invention, they should not be interpreted as limitation of the present invention.
A kind of granite hill reservoir crack prediction method provided by the invention comprising following steps are (by taking the oil fields P as an example
It illustrates)：
1) buried hill top surface paleogeographic reconstruction
Paleogeographic reconstruction is carried out using tradition stratum flattening technique.
P oilfield buried hills top surface is successively evened up from top to bottom using one section of Guantao group, bright hypomere as interface is evened up, bright hypomere
Landforms after interface is evened up represent buried hill top surface, the presedimentary landform shape of group in Guantao, i.e. buried hill described in the present embodiment
Ancient landform.
2) ancient landform contour in buried hill is made picture
To buried hill ancient landform, interpolation generates ancient landform contour map.
3) buried hill ancient landform " ditch " and " ridge " identification
Relatively a height of " ridge " on the ancient landform of buried hill, relatively low is " ditch ".Relatively high, relatively low judgement uses triangle
Shape rule：(1) adjacent three points of connection landform turnover, form triangle；(2) with minimum with horizontal line angle in triangle
On the basis of side, when the 3rd vertex is located on the side, the position of third point represents palaeogeomorphic ridge；(3) when the 3rd
When vertex is located under the side, the position of third point represents palaeogeomorphic ditch.Using the Vector triangle, can identify one by one
Go out P oilfield buried hills ancient landform " ditch " and " ridge ".
4) " ditch " and " ridge " range portraying and extracting
Steps are as follows for the pickup of P oilfield buried hills ancient landform " ditch " " ridge " distributed areas：
A. determine that " ditch " and " ridge " main in working region moves towards direction；
B. since working region is on one side, perpendicular to main " ditch " and " ridge " trend, the 1st arbitrary line R1 is drawn；
C. each terrain transition point (totally three points) on the arbitrary line is determined；
D. judge that 3 terrain transition points are " ditch " or " ridge ", triangular apex on the arbitrary line using Vector triangle
(ditch) is indicated with open circles, and triangular apex (ridge) is indicated with filled circles；
E. it is parallel to the 1st arbitrary line, the 2nd arbitrary line is drawn with 50m distances, the 2nd is determined arbitrarily by a~d methods
" ditch " and " ridge " on line.The 3rd is drawn successively, the 4th arbitrary line ... determine each respectively until working region another side
" ditch " and " ridge " on arbitrary line；
F. " ditch " and " ridge " on above-mentioned arbitrary line is added on paleogeomorphological map；
G. using paleogeomorphological map as auxiliary, " ditch " and " ridge " on the arbitrary line of each item is sequentially connected (and by filled circles with
Open circles distribution connects), (ditch line is open circles line, and crestal line is that filled circles connect for formation " ditch " and " ridge " line distribution map
Line)；
H. ditch ridge distribution is determined：Centered on " ditch " line, both sides are limited using the center line of " ditch " and " ridge " line as boundary
Fixed ranging from buried hill ancient landform " ditch " region (i.e. * fill areas)；Centered on " ridge " line, both sides are with " ridge " and " ditch " line
Center line is boundary, is limited ranging from buried hill ancient landform " ridge " region (oblique dotted line fill area).
5) FRACTURE PREDICTION
Ancient landform " ditch " region is this oil field large fracture development area, and quantity is generally 2~8/m；" ridge " region is this oil
The owner of farmland wants reservoir fissure development area, and to develop hundreds of microns small crack below as main feature, fracture spacing is big, and quantity is general
For tens of~hundreds of/m.
The various embodiments described above are merely to illustrate the present invention, wherein the structure of each component, connection type and manufacture craft etc. are all
It can be varied from, every equivalents carried out based on the technical solution of the present invention and improvement should not exclude
Except protection scope of the present invention.
Claims (3)
1. a kind of granite hill reservoir crack prediction method, which is characterized in that include the following steps：
1) paleogeographic reconstruction is carried out to target area buried hill top surface；
2) to the buried hill top surface after evening up for the last time, i.e. buried hill ancient landform, interpolation generates ancient landform contour map；
3) buried hill ancient landform " ditch " and " ridge " are identified；
4) " ditch " and " ridge " distributed areas are portrayed and is extracted, included the following steps：
A. determine that main " ditch " and " ridge " moves towards direction in the region in crack to be predicted；
B. since the region in crack to be predicted is on one side, perpendicular to main " ditch " and " ridge " trend, the 1st arbitrary line is drawn；
C. each terrain transition point on the arbitrary line is determined；
D. judge that each terrain transition point is " ditch " or " ridge " on the arbitrary line using Vector triangle；
E. it is parallel to the 1st arbitrary line, the 2nd arbitrary line is drawn with 50m distances, step a~d is repeated and determines the 2nd arbitrary line
On " ditch " and " ridge "；Draw successively the 3rd, the 4th arbitrary line ... until crack to be predicted region another side, and respectively
Determine " ditch " and " ridge " on each arbitrary line；
F. " ditch " and " ridge " on above-mentioned arbitrary line is added on the paleogeomorphological map of buried hill top surface；
G. using buried hill top surface paleogeomorphological map as auxiliary, " ditch " and " ridge " on the arbitrary line of each item is sequentially connected, is formed " ditch "
" ridge " line distribution map；
H. " ditch " and " ridge " distributed areas are determined：Centered on " ditch " line, both sides are using the center line of " ditch " and " ridge " line as boundary, institute
Limit ranging from buried hill top surface ancient landform " ditch " distributed areas；Centered on " ridge " line, both sides are with the center of " ridge " and " ditch " line
Line is boundary, limits ranging from buried hill top surface ancient landform " ridge " distributed areas；
5) fracture development area is predicted.
2. a kind of granite hill reservoir crack prediction method as described in claim 1, which is characterized in that above-mentioned steps 1) in
Paleogeographic reconstruction is carried out using tradition stratum flattening technique, i.e., target area buried hill top surface is successively evened up from top to bottom, gradually instead
Buried hill top surface landforms are reflected by newly to old variation；The interfaces at different levels of flattening be each stratum in buried hill overlying sediments stratum " group ",
" section " top surface, the buried hill top surface after lowest part interface is evened up reflect landform shape of the buried hill before receiving to deposit.
3. a kind of granite hill reservoir crack prediction method as described in claim 1, which is characterized in that above-mentioned steps 3) in
Identify buried hill palaeogeomorphic " ditch " and " ridge " one by one using Vector triangle：(1) adjacent three points of connection landform turnover, shape
At triangle；(2) by triangle on the basis of the side of horizontal line angle minimum, when the 3rd vertex is located on the side,
The position of three points represents palaeogeomorphic " ridge "；(3) when the 3rd vertex is located under the side, the position generation of third point
Table palaeogeomorphic " ditch ".
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710083193.5A CN106873032B (en) | 2017-02-16 | 2017-02-16 | A kind of granite hill reservoir crack prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710083193.5A CN106873032B (en) | 2017-02-16 | 2017-02-16 | A kind of granite hill reservoir crack prediction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106873032A CN106873032A (en) | 2017-06-20 |
| CN106873032B true CN106873032B (en) | 2018-09-28 |
Family
ID=59167025
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710083193.5A Active CN106873032B (en) | 2017-02-16 | 2017-02-16 | A kind of granite hill reservoir crack prediction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106873032B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107664772B (en) * | 2017-07-28 | 2019-08-13 | 核工业北京地质研究院 | A kind of prediction technique for sandstone-type uranium mineralization with respect favorab1e sandbadies and Favorable Zones |
| CN108267783B (en) * | 2018-01-02 | 2019-08-06 | 中国石油天然气集团公司 | A kind of method, apparatus and system of determining buried-hill trap |
| CN109283578B (en) * | 2018-01-22 | 2020-07-10 | 中国石油天然气股份有限公司 | Method and device for determining the position of a hydrocarbon supply window of a submerged hill |
| CN114594529B (en) * | 2020-12-07 | 2023-12-01 | 中国石油化工股份有限公司 | Method for predicting reservoir of buried hill based on azimuth seismic attribute change rate |
| CN113433589A (en) * | 2021-08-03 | 2021-09-24 | 成都理工大学 | Weathered crust reservoir bottom interface identification method based on mathematical statistics |
| CN114019565B (en) * | 2021-11-02 | 2023-10-13 | 中海石油(中国)有限公司 | Rapid prestack crack prediction method, rapid prestack crack prediction device, computer storage medium and computer storage device |
| CN114265116B (en) * | 2022-03-02 | 2022-05-17 | 中海油研究总院有限责任公司 | Quantitative statistical method and device for channel and ridge amplitude of granite buried hill weathering crust |
Citations (2)
| 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 |
| CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
-
2017
- 2017-02-16 CN CN201710083193.5A patent/CN106873032B/en active Active
Patent Citations (2)
| 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 |
| CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
Non-Patent Citations (5)
| Title |
|---|
| Fine buried hill reservoir characterization technology application research;Jiannan Zhou et al.;《SPG/SEG Beijing 2016 International Geophysical Conference》;20161231;第29-32页 * |
| P 潜山油田基于组合地震属性的裂缝预测研究;聂妍等;《SPG/SEG 北京2016国际地球物理会议》;20160930;第374-377页、第363-365页 * |
| 花岗岩潜山储层裂缝建模表征方法—以渤海花岗岩潜山A油田为例;张雨晴等;《吉林大学学报(地球科学版)》;20160930;第46卷(第5期);第1313-1315页，图3-4 * |
| 鄂尔多斯盆地中部气田剥蚀脊与沟槽;代金友等;《石油勘探与开发》;20051231;第32卷(第6期);第29-31页 * |
| 靖边古潜台剥蚀特征精细刻画及对气藏发育的影响;罗国文等;《石油化工应用》;20120731;第31卷(第7期);第33-37页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106873032A (en) | 2017-06-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106873032B (en) | A kind of granite hill reservoir crack prediction method | |
| CN104992468B (en) | Fracture-cavity type carbonate oil-gas reservoir three-dimensional geological modeling method | |
| CN104564041B (en) | Hyposmosis clastic reservoir rock efficiency evaluation method based on exploitation permeability limits | |
| US9026420B2 (en) | Method and system of updating a geocellular model | |
| CN102748016A (en) | Geologic history period sandstone reservoir porosity evolution recovery method | |
| CN109709301B (en) | Method, device and system for determining distribution of fracture pore type tight sandstone reservoir | |
| CN104166163A (en) | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube | |
| CN105626023A (en) | Well test determination method for vertical fracturing fracture azimuth of low-permeability oil reservoir | |
| CN105651962B (en) | Diagenesis phase identification method | |
| CN110410027A (en) | It is a kind of to carry out the continuous evaluation method of rock-mass quality and system using drilling core | |
| Beddingfield et al. | Polygonal impact craters on Dione: Evidence for tectonic structures outside the wispy terrain | |
| CN104597487B (en) | A kind of fault description method | |
| Humphris et al. | A new radiocarbon chronology for ancient iron production in the Meroe region of Sudan | |
| CN105888659A (en) | Method and device for determining lithologic oil-gas reservoir forming probability | |
| CN115185015B (en) | Deep lithium beryllium ore investigation method | |
| CN109583016A (en) | A kind of space geometry form of fracture hole body is determining and quantitative calculation method | |
| CN102562051A (en) | Method for processing original data of gas-measuring plate | |
| CN107218969B (en) | A kind of phase of determining oil-gas migration time, the time, power and direction method | |
| CN109446744A (en) | A kind of stand growth simulating method and system for considering space structure and being interacted with growth | |
| CHEN et al. | 3D geological modeling of dual porosity carbonate reservoirs: A case from the Kenkiyak pre-salt oilfield, Kazakhstan | |
| CN104749618B (en) | Mud shale low-angle dip crack poststack probability quantitatively characterizing method | |
| CN106504319B (en) | Reservoir Three Dimensional Contrast map generalization method and device between well | |
| Howland | Long-distance trade and social complexity in Iron Age Faynan, Jordan | |
| CN106483560B (en) | A method of obtaining the normal throw of fracture | |
| CN115437013A (en) | Oil reservoir small slip distance slip fracture analysis method, medium and equipment |
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 | ||
| CB02 | Change of applicant information |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant after: China Offshore Oil Group Co., Ltd. Applicant after: CNOOC research institute limited liability company Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant before: China National Offshore Oil Corporation Applicant before: CNOOC Research Institute |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |