CN106547034B - A kind of method of determining compact reservoir rock brittleness index - Google Patents
A kind of method of determining compact reservoir rock brittleness index Download PDFInfo
- Publication number
- CN106547034B CN106547034B CN201610982279.7A CN201610982279A CN106547034B CN 106547034 B CN106547034 B CN 106547034B CN 201610982279 A CN201610982279 A CN 201610982279A CN 106547034 B CN106547034 B CN 106547034B
- Authority
- CN
- China
- Prior art keywords
- rock
- young
- modulus
- angle
- brittleness index
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Abstract
The present invention provides a kind of method quickly calculating rock brittleness based on well-log information, to realize the quick calculating of compact reservoir rock brittleness under formation condition, promotes fine and close oil-gas reservoir Efficient Development.This method is more complicated according to rock rupture form mainly from rock rupture form and fracture energy angle, and fracture can smaller, the bigger principle of rock brittleness index.Reflect the ability of rock resistance to deformation, the rate of deformation respectively using Young's modulus, dilative angle, can be very good the Brittleness of description rock, while Young's modulus and dilative angle can also reflect the feature of different phase on stress-strain curve of rock.The weight coefficient that Young's modulus and dilative angle are calculated using grey correlation theory, to set up the new method for calculating rock brittleness index.This method is the computational methods based on well-log information fast prediction brittleness index, improves the reasonability of rock brittleness evaluation, promotes the Efficient Development of fine and close oil-gas reservoir.
Description
Technical field：
The present invention relates to unconventionaloil pool exploration and development fields, and in particular to one kind calculating fine and close storage based on well-log information
The method of layer rock brittleness index.
Background technology：
The fine and close oil and gas reserves in China is abundant, has good potentiality to be exploited, but development difficulty is larger.By shale gas reservoir
The inspiration of effective exploitation, fine and close oil-gas reservoir have carried out the practice of seam net pressure break.In stitching net pressure break, the brittleness of rock is to judge storage
The important parameter to form seam net ability is split in lamination.At present petroleum industry oil and gas development field using most mineral brittleness indexs and
The method that dynamic rock mechanics parameters calculate rock brittleness is all the summary of U.S.'s shale gas pressure break experience.And rock mechanics field
Rock brittleness be mainly based upon again rock indoor experiment stress-strain diagram and mechanical characteristics calculated, without and oil gas
Development field combines, it is difficult to be used for field application.Therefore, fine and close oil gas pressure break is badly in need of a kind of based on rock brittleness mechanics table
Now with the convenient rock brittleness computational methods of field application.
Invention content：
Deficiency of the present invention for the rock brittleness computational methods of fine and close oil gas field application at present, it is proposed that one kind is based on
The method that rock fracture energy, ruptured type and well-log information quickly calculate rock brittleness index causes to realize under formation condition
The close brittle quick calculating of reservoir rock, foundation is provided for compact reservoir design of hydraulic fracturing scheme.
The microcosmic fracture of rock can be the specific manifestation of rock brittleness with the ruptured type of macroscopic view.It is proposed by the present invention
A method of compact reservoir rock brittleness index is calculated, reflects rock resistance to deformation respectively using Young's modulus, dilative angle
The rate of ability and deformation has the principle of the Brittleness of description rock well.By assigning the corresponding weights of two parameter
To establish new rock brittleness computational methods, and the parameter value for calculating brittleness index is obtained using well-log information.Master of the present invention
Include the following steps：
1. Well logging Data prepares.Well-log information includes full wave train acustic logging or dipole acoustic and natural gamma
The data such as well logging include at least shear wave slowness, compressional wave time difference, density, porosity, the Natural Gamma-ray Logging Curves of fractured interval；
2. calculating separately out the Young's modulus of rock using shear wave slowness, compressional wave time difference, density curve, gamma ray curve
E, Poisson's ratio ν, internal friction angleShale content VshSection uses formula (1)-(5) to calculate respectively：
In formula, EdFor kinetic Youngs modulus, νdFor dynamic Poisson's ratio, ρ is density, VtFor velocity of longitudinal wave, VsFor shear wave speed
Degree,For internal friction angle, φ is porosity, and SH is natural gamma relative value, GR, GRmax、GRminRespectively target zone natural gamma
Value, the natural gamma value on pure shale stratum, the natural gamma value of pure lithologic character stratum, VshFor shale content, GCUR is to have with the age
The empirical coefficient of pass, old stratum take 3.7, and new stratum takes 2.
3. the rock-mechanics property parameter and shale content section that are obtained using step 2 be calculated compressive strength of rock and
Confining pressure section uses formula (6) and (7) to calculate respectively：
σc=(0.0045+0.0035Vsh)Ed (6)
In formula, σcFor the compression strength of rock, PcFor confining pressure, H is the depth of reservoir, and α is effective stress coefficient.
As shale content VshWhen > 0.8, effective stress factor alpha=0.6；As shale content VshWhen < 0.2, effective stress
Factor alpha=0.9；As 0.2 < V of shale contentshWhen < 0.8, effective stress coefficient
4. the dilative angle section of the compressive strength of rock and calculating of surrounding rock pressure rock that are obtained using step 3 is counted with formula (8)
It calculates；
In formula, ψ is the dilative angle of rock,For internal friction angle, σcFor compressive strength of rock, PcFor confining pressure.
5. establishing new rock brittleness index number evaluation method using Young's modulus and dilative angle the two parameters, pass through grey
Relevance theory calculates the weight coefficient of Young's modulus and dilative angle；
6. the weight coefficient that the dilative angle and step 5 that are obtained with Young's modulus, step 4 that step 2 obtains obtain, calculates
Rock brittleness index profile is obtained, rock brittleness index B is calculated with formula (9)I。
BI=W1En+W2ψn (9)
In formula, BIFor brittleness index, EnFor normalized Young's modulus, ψnFor normalized dilative angle, W1For Young's modulus
Shared weight, W2For weight shared by dilative angle.
Wherein, EnAnd ψnCalculation formula it is as follows：
In formula, E, Emax, EminThe respectively Young's modulus of stratum arbitrary point, fractured interval stratum maximum Young's modulus, pressure
Parting section stratum minimum Young's modulus；ψ, ψmax, ψminRespectively the dilative angle of stratum arbitrary point, fractured interval stratum maximum are cut swollen
Angle, fractured interval stratum minimum dilative angle.
Compared with the existing method for calculating rock brittleness index, the method have the advantages that：The present invention is based on rocks
Fracture energy and ruptured type, it is proposed that one kind can utilize well-log information fast prediction brittleness index computational methods, by rock
The brittleness of stone combines with the ruptured type of laboratory experiment rock, provides more intuitive judgement for subsurface rock brittleness, is
Fine and close oil-gas reservoir Fracturing Project optimization provides foundation, is conducive to improve fine and close oil gas fracturing effect.
Description of the drawings：
Fig. 1 is the step block diagram that the present invention calculates rock brittleness.
Fig. 2 is Dagang Oilfield compact oil reservoir 4050m~4065m well section well-log information schematic diagrames of the present invention
Fig. 3 is to calculate compact oil reservoir 4050m~4065m well section brittleness index diagrammatic cross-sections according to the method for the present invention
Specific implementation mode：
Implementation method of the present invention is described in further detail below, it is specific as follows：
By taking a certain compact oil reservoir reservoir of Dagang Oilfield aperture leakage as an example, which is 4050~4065m, the well section
Layer average pressure is 45MPa.For at buried depth of strata 4050m, brittleness index is calculated using the method for the present invention.
1. according to step 1, the shear wave slowness, compressional wave time difference, density, porosity and the natural gamma that prepare fractured interval are surveyed
Well curve, Well logging Data are as shown in Fig. 2；
2. calculating separately the Young's modulus of rock, Poisson's ratio, internal friction angle, shale using related logging data to contain
Amount, it is E=36580.5MPa, Poisson that the Young's modulus at buried depth 4050m, which is calculated, by formula (1)-(5) in step 2
Than being for v=0.258, internal friction angleShale content is Vsh=0.209；
3. it is σ to calculate separately to obtain the compressive strength of rock at buried depth 4050m according to formula (6)-(7) in step 3c
=191.44MPa, confining pressure Pc=59.85MPa, according to the formula (8) in step 4 be calculated rock dilative angle be ψ=
5.69°；
4. the weight coefficient according to the two parameters of calculating Young's modulus and dilative angle in step 5 is respectively 0.426 He
0.574, and it is that rock brittleness index is B at 4050m that buried depth, which is calculated, according to formula (9)-(11) in step 6I=
0.824。
Similarly, method through the invention can be calculated the brittleness that buried depth is 4050~4065m compact reservoir rocks and refer to
Number section, as a result as shown in Fig. 3.
Claims (5)
1. a kind of method calculating compact reservoir rock brittleness index, includes the following steps：
(1) Well logging Data prepares：Well-log information includes that full wave train acustic logging or dipole acoustic and natural gamma are surveyed
Well data include at least shear wave slowness, compressional wave time difference, density, porosity, Natural Gamma-ray Logging Curves；
(2) the Young mould of rock is calculated separately out using shear wave slowness, compressional wave time difference, density, porosity, gamma ray curve
Amount, Poisson's ratio, internal friction angle, shale content；
(3) utilize the obtained rock-mechanics property parameter of step (2) and shale content section be calculated compressive strength of rock and
Confining pressure section；
(4) the dilative angle section of the compressive strength of rock and calculating of surrounding rock pressure rock obtained using step (3)；
(5) it utilizes Young's modulus and dilative angle the two parameters to establish new rock brittleness index number evaluation method, is closed by grey
Join the weight coefficient of theoretical calculation Young's modulus and dilative angle；
(6) weight coefficient that the dilative angle and step (5) that Young's modulus, the step (4) for using step (2) to obtain obtain obtain,
Rock brittleness index profile is calculated.
2. according to a kind of method of the above-mentioned calculating compact reservoir rock brittleness index of claim 1, it is characterised in that：The step
Suddenly the Young's modulus in (2), Poisson's ratio, internal friction angle, the calculation formula of shale content are：
In formula, EdFor kinetic Youngs modulus, νdFor dynamic Poisson's ratio, ρ is density, VtFor velocity of longitudinal wave, VsFor shear wave velocity,For
Internal friction angle, φ are porosity, and SH is natural gamma relative value, GR, GRmax、GRminRespectively target zone natural gamma value, pure
The natural gamma value of the natural gamma value of mudstone stratum, pure lithologic character stratum, VshFor shale content, GCUR is warp related with the age
Coefficient is tested, old stratum takes 3.7, and new stratum takes 2.
3. according to a kind of method of the above-mentioned calculating compact reservoir rock brittleness index of claim 1, it is characterised in that：The step
Suddenly the compressive strength of rock in (3) and calculating of surrounding rock pressure formula are：
σc=(0.0045+0.0035Vsh)Ed (6)
In formula, σcFor the compression strength of rock, PcFor confining pressure, H is the depth of reservoir, and α is effective stress coefficient；
As shale content VshWhen > 0.8, effective stress factor alpha=0.6；As shale content VshWhen < 0.2, effective stress factor alpha
=0.9；As 0.2 < V of shale contentshWhen < 0.8, effective stress coefficient
4. according to a kind of method of the above-mentioned calculating compact reservoir rock brittleness index of claim 1, it is characterised in that：The step
Suddenly the dilative angle calculation formula in (4) is：
In formula, ψ is the dilative angle of rock,For internal friction angle, σcFor compressive strength of rock, PcFor confining pressure.
5. according to a kind of method of the above-mentioned calculating compact reservoir rock brittleness index of claim 1, it is characterised in that：The step
Suddenly the rock brittleness index B in (6)ICalculation formula be：
BI=W1En+W2ψn (9)
In formula, BIFor brittleness index, EnFor normalized Young's modulus, ψnFor normalized dilative angle, W1Shared by Young's modulus
Weight, W2For weight shared by dilative angle；
Wherein, EnAnd ψnCalculation formula it is as follows：
In formula, E, Emax, EminThe respectively Young's modulus of stratum arbitrary point, stratum maximum Young's modulus, stratum minimum Young mould
Amount；ψ, ψmax, ψminThe respectively dilative angle of stratum arbitrary point, stratum maximum dilative angle, stratum minimum dilative angle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610982279.7A CN106547034B (en) | 2016-11-09 | 2016-11-09 | A kind of method of determining compact reservoir rock brittleness index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610982279.7A CN106547034B (en) | 2016-11-09 | 2016-11-09 | A kind of method of determining compact reservoir rock brittleness index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106547034A CN106547034A (en) | 2017-03-29 |
| CN106547034B true CN106547034B (en) | 2018-11-02 |
Family
ID=58395393
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610982279.7A Active CN106547034B (en) | 2016-11-09 | 2016-11-09 | A kind of method of determining compact reservoir rock brittleness index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106547034B (en) |
Families Citing this family (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107045580B (en) * | 2017-04-27 | 2020-04-21 | 中国石油大学(华东) | Shale mechanical parameter rapid calculation method based on digital rock core |
| CN107747916B (en) * | 2017-10-27 | 2019-07-30 | 辽宁工程技术大学 | A kind of measuring method of strain localization band dilative angle |
| CN108519281B (en) * | 2018-02-27 | 2020-08-11 | 中国石油天然气股份有限公司 | Method, device and system for determining brittleness index of rock |
| CN108875115B (en) * | 2018-03-12 | 2022-05-10 | 中国石油天然气集团有限公司 | Method for determining rock strength |
| CN108805158B (en) * | 2018-04-16 | 2021-06-01 | 北京师范大学 | Compact oil reservoir lithogenous phase dividing method |
| CN109164499B (en) * | 2018-08-30 | 2020-04-17 | 中国石油化工股份有限公司 | Rock brittleness index calculation and explanation method considering sedimentary diagenesis process |
| CN109917102A (en) * | 2019-01-28 | 2019-06-21 | 西南石油大学 | A kind of brittleness evaluation method considering rock rupture complexity |
| CN109778880B (en) * | 2019-02-14 | 2020-12-15 | 长安大学 | Construction method of geocell retaining wall |
| CN109828031B (en) * | 2019-02-15 | 2021-07-23 | 西南石油大学 | Rock brittleness evaluation method and device |
| CN110058323A (en) * | 2019-04-03 | 2019-07-26 | 天津科技大学 | A kind of tight sand formation brittleness index calculation method |
| CN110288233B (en) * | 2019-06-26 | 2023-08-08 | 重庆科技学院 | Deep shale gas compressibility evaluation method based on fuzzy gray correlation method |
| CN110566171A (en) * | 2019-07-15 | 2019-12-13 | 西南石油大学 | Ultrahigh pressure tight fractured sandstone gas reservoir sand production prediction method |
| CN110426751A (en) * | 2019-08-09 | 2019-11-08 | 中国石油天然气股份有限公司 | A method of shear wave slowness is predicted using well-log information |
| CN112649858A (en) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | Shale brittleness prediction method and system based on core test |
| CN110618456A (en) * | 2019-10-24 | 2019-12-27 | 西南石油大学 | Method for determining years of multi-period active fracture of sedimentary basin |
| CN111396057A (en) * | 2020-03-04 | 2020-07-10 | 中国石油大学(华东) | Deep shale reservoir brittleness evaluation method and computer readable storage medium |
| CN111366452B (en) * | 2020-03-26 | 2021-01-29 | 北京科技大学 | Method for measuring energy storage level of self-energy-storage rock mass |
| CN113756793B (en) * | 2020-06-05 | 2023-09-26 | 中国石油天然气股份有限公司 | Shale oil exploitation mode determining method, device and equipment and readable storage medium |
| CN113818862A (en) * | 2020-06-05 | 2021-12-21 | 中国石油天然气股份有限公司 | Method, device and equipment for dividing shale oil exploitation area and readable storage medium |
| CN113033021A (en) * | 2021-04-21 | 2021-06-25 | 北京艾迪博科油气技术有限公司 | Method and device for reforming and evaluating tight reservoir |
| CN113189647B (en) * | 2021-04-30 | 2022-03-11 | 西南石油大学 | Method for predicting formation brittleness index of transverse isotropic shale |
| CN113189648B (en) * | 2021-04-30 | 2022-03-11 | 西南石油大学 | Method for predicting brittleness index of orthotropic shale |
| CN113625358B (en) * | 2021-08-11 | 2024-01-19 | 中国石油化工股份有限公司 | Method for judging influence degree of rock components on physical properties of tight sandstone reservoir |
| CN114235579B (en) * | 2021-12-17 | 2023-12-19 | 广东石油化工学院 | Triaxial testing device and method for compression resistance and shear strength of natural gas hydrate core |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103344705A (en) * | 2013-06-25 | 2013-10-09 | 中国石油大学(北京) | Method of measuring rock brittleness index through applying acoustic emission energy values |
| CN104677734A (en) * | 2015-02-12 | 2015-06-03 | 中国石油大学(华东) | Multi-scale information integration tight formation brittleness index measuring method |
| CN105221141A (en) * | 2014-06-23 | 2016-01-06 | 中国石油化工股份有限公司 | A kind of mud shale brittleness index Forecasting Methodology |
| CN105221140A (en) * | 2014-06-20 | 2016-01-06 | 中国石油化工股份有限公司 | A kind ofly determine that shale formation can the method for pressure break sex index |
| CN105675635A (en) * | 2015-12-31 | 2016-06-15 | 中国石油天然气股份有限公司 | Method for determining relative content of components of compact rocks and brittleness index of compact rocks, and apparatus thereof |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US216609A (en) * | 1879-06-17 | Improvement in hotel and restaurant waiters | ||
| NO863745L (en) * | 1985-09-19 | 1987-03-20 | Seislith Dev Inc | SEISMIC INVESTIGATION OF LAW. |
| GB0017227D0 (en) * | 2000-07-14 | 2000-08-30 | Schlumberger Ind Ltd | Fully coupled geomechanics in a commerical reservoir simulator |
-
2016
- 2016-11-09 CN CN201610982279.7A patent/CN106547034B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103344705A (en) * | 2013-06-25 | 2013-10-09 | 中国石油大学(北京) | Method of measuring rock brittleness index through applying acoustic emission energy values |
| CN105221140A (en) * | 2014-06-20 | 2016-01-06 | 中国石油化工股份有限公司 | A kind ofly determine that shale formation can the method for pressure break sex index |
| CN105221141A (en) * | 2014-06-23 | 2016-01-06 | 中国石油化工股份有限公司 | A kind of mud shale brittleness index Forecasting Methodology |
| CN104677734A (en) * | 2015-02-12 | 2015-06-03 | 中国石油大学(华东) | Multi-scale information integration tight formation brittleness index measuring method |
| CN105675635A (en) * | 2015-12-31 | 2016-06-15 | 中国石油天然气股份有限公司 | Method for determining relative content of components of compact rocks and brittleness index of compact rocks, and apparatus thereof |
Non-Patent Citations (5)
| Title |
|---|
| "Micro-fracture characteristics of tight sandstone reservoirs and its evaluation by capillary pressure curves: A case study of Permian sandstones in Ordos Basin, China";Peng Zhu 等;《Journal of Natural Gas Science and Engineering》;20150822(第27期);第90-97页 * |
| "不同围压下煤岩变形与剪胀扩容模型";王云飞 等;《水文地质工程地质》;20150131;第42卷(第1期);第106-111页 * |
| "基于岩石物理实验的页岩脆性测井评价方法";徐赣川 等;《天然气工业》;20141231;第38-45页 * |
| "复杂储层岩石脆性分析及应用研究";许孝凯 等;《测井技术》;20150831;第39卷(第4期);第486-490页 * |
| "岩石剪胀角模型与验证";赵星光 等;《岩石力学与工程学报》;20100521;第29卷(第5期);第970-981页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106547034A (en) | 2017-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106547034B (en) | A kind of method of determining compact reservoir rock brittleness index | |
| CN105156103B (en) | A kind of multiple dimensioned shale reservoir three-dimensional compressibility evaluation method of landwaste-rock core-wellbore-reservoir | |
| Taleghani et al. | Numerical simulation of hydraulic fracture propagation in naturally fractured formations using the cohesive zone model | |
| CN105134156B (en) | A kind of modeling method for tight sandstone reservoir three-dimensional compressibility model | |
| Guo et al. | Evaluation of fracability and screening of perforation interval for tight sandstone gas reservoir in western Sichuan Basin | |
| CN108009705A (en) | A kind of shale reservoir compressibility evaluation method based on support vector machines technology | |
| CN105484710B (en) | A method of selection Deviated Well Fracturing perforation orientation | |
| Lu et al. | Engineering geological characteristics and the hydraulic fracture propagation mechanism of the sand-shale interbedded formation in the Xu5 reservoir | |
| Suppachoknirun et al. | Hydraulic fracturing and production optimization in eagle ford shale using coupled geomechanics and fluid flow model | |
| Gao et al. | Fracability evaluation of lacustrine shale in the Yanchang Formation of southeastern Ordos Basin | |
| Igor et al. | Specifics of mechanical and strength rock properties estimation for wells drilling and exploitation | |
| Ligtenberg et al. | Fractures in the Dutch Rotliegend—An Overview | |
| Giroldi* et al. | Applications of broadband seismic inversion in the assessment of drilling and completion strategies: A case study from eastern Saudi Arabia | |
| Ge et al. | Analytical modeling on 3D stress redistribution and fault reactivation during hydraulic fracturing stimulation | |
| Wan et al. | Numerical simulation and experimental study of near-wellbore fracture initiation mechanism on sandstone coal interbedding | |
| Al Mulhim | Optimization of hydraulic fracture treatments and landing zone intervals within the Eagle Ford | |
| Li et al. | Numerical Investigation of multi-well, multi-stage hydraulic fracture height growth in laminated shale reservoirs using finite-discrete element method | |
| Liu et al. | Fracturing Optimization Design of Fractured Volcanic Rock in Songliao Basin Based on Numerical Research and Orthogonal Test | |
| Liu et al. | Fracability evaluation method and influencing factors of the tight sandstone reservoir | |
| Hernandez et al. | Drill Cutting Reinjection Feasibility Study in a Critical Environment, Apaika Field, Ecuador | |
| Haggerty et al. | Investigating the Impact of Dynamic Effective Stress While Perforating Low-Strength Formations | |
| Bauer et al. | Numerical model of extended leak-off test (XLOT) | |
| Lv et al. | Hydraulic fracture initiation while staged fracturing for horizontal wells | |
| Hudyma et al. | Variability of weathered limestone properties in light of drilled shaft design in Florida | |
| Hashemian et al. | Study of the relationship between fracture intensity distribution and geomechanical parameters in Sarvak reservoir, in an oilfield located in Dezful Embayment |
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 |