CN106483575B - A kind of image-characterization methods of fault activities - Google Patents
A kind of image-characterization methods of fault activities Download PDFInfo
- Publication number
- CN106483575B CN106483575B CN201510543884.XA CN201510543884A CN106483575B CN 106483575 B CN106483575 B CN 106483575B CN 201510543884 A CN201510543884 A CN 201510543884A CN 106483575 B CN106483575 B CN 106483575B
- Authority
- CN
- China
- Prior art keywords
- fault activities
- fault
- activities
- intensity
- intensity map
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of image-characterization methods of fault activities, the method comprises the steps of: that obtaining fault activities to be analyzed is walking upward different statistics points in the fault activities intensity of specific one or multiple tectonic epochs；It calculates separately the acquisition different statistics points and walks upward position in the fault activities to be analyzed；The fault activities intensity map of the corresponding tectonic epochs is generated based on the fault activities intensity and the position.Compared with prior art, image-characterization methods of the invention can be with the concrete condition of more specific intuitive characterization fault activities intensity；Meanwhile present invention employs the modes of quantitatively characterizing, so that characterization result is more in line with practical geological condition.
Description
Technical field
The present invention relates to geological exploration fields, in particular relate to a kind of image-characterization methods of fault activities.
Background technique
Active procedure is broken in sedimentary basin and Gas Accumulation is closely bound up, and oil gas is usually dredged via fracture and enters circle
Among closing, also destruction can be played to oil-gas reservoir in subsequent activity.Therefore, the activity intensity of fracture is shown by stages secondaryly
It shows and, it is particularly important to the sealing process of oil gas for analysis fracture, so far, characterize the technical side of fault activities intensity
There are mainly two types of methods.
One: fault growth index (Growth index) method
Fault growth Index Definition is the decline disc thickness for being broken the same stratigraphic unit in two sides by C.E.Thorsen (1963)
With rise the ratio between disc thickness.He thinks the size for comparing different times fault growth index on growth index figure, it will be appreciated that
Activity intensity of the tomography in different times.Fault growth index can generally illustrate:
(1) fracture starts the movable time；
(2) fault activities most strong period, i.e. downthrow block stratum increase the maximum epoch；
(3) the newest epoch that the latter stage of fault activities, i.e. downthrow block stratum increase.
According to the viewpoint of C.E.Thorsen, the power of the size reflection fault activities of fault growth index, growth index etc.
In 1, it is broken inactive；Greater than 1, fault activities, growth index is bigger, and Activity of Faults is stronger.
The shortcomings that fault growth index method:
In practical geological activity, fault growth index is not only related with fault activities rate, also related with deposition rate.
Therefore the actual geological condition of reflection that fault growth index method can not be stringent.In addition, fault growth index map is only to disconnected
Some location point split cannot reflect fracture overall condition at figure.
Two: balanced cross section method
Balanced cross section method is all to be recovered to deformation structure on the section of vertical configuration trend by geometry principle
A kind of analogue technique of reasonable undeformed state.It can carry out quantitative, sxemiquantitative analysis interpretation to stratum structural evolution,
It is the important means of structural evolution quantitative analysis, is used widely in fields such as oil-gas exploration, basin modellings.
The shortcomings that balanced cross section method:
Balanced cross section is entire structure section shortening amount in the horizontal direction and shortening at the quantitative data that figure method provides
Rate, shortening amount not only have the contribution of fracture, more there is the contribution of fold.Therefore the crack conditions that balanced cross section reflects at figure method are simultaneously
Fitting actual conditions that cannot be stringent.Also, this drawing methods do not mark the horizontal shortening amount of each fracture.Together
Sample, due to the trend at figure section perpendicular to construction, this method also just for some location point in fracture at figure,
Rather than figure is fragmented into entire.
To sum up, in the prior art, the characterization result of fault activities intensity is unsatisfactory.Therefore, in the prior art
Fault activities characterizing method there are the problem of, need a kind of new fault activities characterizing method.
Summary of the invention
Fault activities characterizing method in the prior art there are aiming at the problem that, the present invention provides a kind of figures of fault activities
As characterizing method, the method is comprised the steps of:
It obtains fault activities to be analyzed and is walking upward different statistics points in the fracture of specific one or multiple tectonic epochs
Activity intensity；
It calculates separately the acquisition different statistics points and walks upward position in the fault activities to be analyzed；
The fault activities intensity map of the corresponding tectonic epochs is generated based on the fault activities intensity and the position.
In one embodiment, the fault activities intensity map is constituted by multiple column blocks are transversely arranged, in which:
The fault activities intensity map is walking upward distance as abscissa using fault activities to be analyzed, and the fracture is lived
The length of fatigue resistance map is the total length of the fault activities to be analyzed；
Each corresponding statistics point of the column block, the column block is on the fault activities intensity map
The corresponding corresponding statistics point in position is walking upward position in the fault activities to be analyzed.
In one embodiment, using fault activities intensity value as the ordinate of the fault activities intensity map, the column
The height of shape block corresponds to the statistics point in the fault activities intensity of the current tectonic epochs.
In one embodiment, it during generating the fault activities intensity map, is adopted for the different tectonic epochs
With unified ordinate ratio.
In one embodiment, the color of the column block, the different face are determined based on the specific colour code rule of correspondence
Color corresponds to the value of the different fault activities intensity, and the color of the column block corresponds to the statistics point in the current structure
Make the fault activities intensity of phase.
In one embodiment, it during generating the fault activities intensity map, is adopted for the different tectonic epochs
With the unified colour code rule of correspondence.
In one embodiment, the method also includes following steps:
For the same fault activities, longitudinal combination arranges multiple and different fault activities intensity maps to obtain
More phase fault activities maps.
In one embodiment, during combining the fault activities intensity map of the different configuration phase, using unified
Abscissa and keep the consistent of all fault activities intensity map abscissa ratios.
In one embodiment, during combining the fault activities intensity map of the different configuration phase, more phases
With the interval of specific range between the two neighboring fault activities intensity map on fault activities map.
In one embodiment, during combining the fault activities intensity map of the different configuration phase, in the longitudinal direction
The fault activities intensity map of the different configuration phase is arranged sequentially in time.
Compared with prior art, image-characterization methods of the invention can be with more specific intuitive characterization fault activities intensity
Concrete condition；Meanwhile present invention employs the modes of quantitatively characterizing, so that characterization result is more in line with practical geological condition.
Other feature or advantage of the invention will illustrate in the following description.Also, Partial Feature of the invention or
Advantage will be become apparent by specification, or be appreciated that by implementing the present invention.The purpose of the present invention and part
Advantage can be realized or be obtained by step specifically noted in the specification, claims and drawings.
Detailed description of the invention
Attached drawing is used to provide further understanding of the present invention, and constitutes part of specification, with reality of the invention
It applies example and is used together to explain the present invention, be not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is an embodiment method execution flow chart according to the present invention；
Fig. 2 is more phase fault activities maps according to an embodiment of the invention；
Fig. 3 is more phase fault activities maps according to another embodiment of the present invention.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, implementation personnel of the invention whereby
Can fully understand that how the invention applies technical means to solve technical problems, and reach technical effect realization process and according to
The present invention is embodied according to above-mentioned realization process.As long as each embodiment it should be noted that do not constitute conflict, in the present invention
And each feature in each embodiment can be combined with each other, be formed by technical solution protection scope of the present invention it
It is interior.
For more specific intuitive characterization fault activities, the invention proposes a kind of characterization image sides of fault activities
Method.Next the implementation procedure of method of the invention is specifically described based on flow chart, step shown in the flow chart of attached drawing can
To be executed in the computer system comprising such as a group of computer-executable instructions.Although showing each step in flow charts
Logical order but in some cases can be with the steps shown or described are performed in an order that is different from the one herein.
The main purpose of method of the invention is specific intuitive characterization fault activities intensity, therefore as shown in Figure 1, is executed
Method of the invention first has to execute step S110, obtains breaking strength.Geological activity can be divided into according to the time of generation
Multiple tectonic epochs, it is often movable in multiple tectonic epochs for same fracture.In the present embodiment, not for tectonic epochs
Fault activities are distinguished with (the specific activity time of fault activities is different).To be analyzed according to the determination of specific application demand
Specific one or multiple tectonic epochs, the fault activities intensity in embodiments above-mentioned specific one or multiple tectonic epochs.
When characterizing breaking strength, the statistics point characteristic manner that the present embodiment uses is moved towards for fault activities
On the breaking strengths of multiple statistics points characterized that (ideally, statistics point number is more, just to the descriptions of fault activities
It is more specific.In practice, but it is considered that statistics point selections needs according to Geological data collecting and the actual conditions of processing with
And data characterization requires to determine.).
In summary characterization require, in the step S110 of the present embodiment, obtain fault activities to be analyzed walk it is upward
Difference statistics point is in the fault activities intensity of specific one or multiple tectonic epochs.
Fault activities intensity due to walking statistics points different upwards in fault activities is different, in order to specific in detail
Fault activities are described, it is necessary to determine the detail location of statistics point.Therefore in the present embodiment, step S120 is also executed, is obtained
Statistics point position step is taken, the different statistics points of acquisition is calculated separately and walks upward position in fault activities to be analyzed.
Step S140 can be executed based on step S110 and step S120, constructs fault activities intensity map, i.e. base
In the fault activities intensity map of fault activities intensity and the corresponding specific tectonic epochs of statistics point position generation on statistics point.
In the present embodiment, fault activities intensity map is constituted by multiple column blocks are transversely arranged.Specifically, with to be analyzed
Fault activities are walking the abscissa that upward distance is fault activities intensity map, enable the length of fault activities intensity map corresponding
The total length of fault activities to be analyzed；The corresponding statistics point of each column block on fault activities intensity map, column
The corresponding corresponding statistics point in position of the block on fault activities intensity map is being walked upwards in the fault activities to be analyzed
Position.In this way, fault activities intensity map, which is composed, can intuitively show that fault activities walk the position of upward each statistics point.
Further, it is living that the fracture accordingly counted on point is described by each column block on fault activities intensity map
Fatigue resistance, fault activities intensity map, which is composed, can intuitively show that the intensity of fault activities moves towards situation.
In the present embodiment, it uses in a manner of characterization intensity.The specific colour code rule of correspondence is established, with difference
Color correspond to different fault activities intensity values.The color of column block is determined based on the colour code rule of correspondence.In this way, column block
Color can correspond to statistics point in the fault activities intensity of current tectonic epochs.
As shown in Fig. 2, the color bar 200 of bottom is exactly the colour code rule of correspondence table constructed, 6 horizontal colors of top
Item is exactly the fault activities intensity map for respectively corresponding 6 different configuration phases.
It in other embodiments of the invention, can also be strong using the fault activities of other modes embodiments statistics point
Degree.In an alternative embodiment of the invention, using fault activities intensity value as the ordinate of fault activities intensity map, column block
The fault activities intensity of the corresponding statistics point of height in the current tectonic epochs.
It is often movable in multiple tectonic epochs for same fracture.In order to comprehensively describe same fracture in multiple structures
The activity condition in the phase is made, combines the fault activities intensity map of multiple and different tectonic epochs with longitudinal arrangement in the present embodiment
Mode generate more phase fault activities maps.
Next the generating process of more phase fault activities maps is specifically described.Since more phase fault activities maps are based on more
The fault activities intensity map composite construction of a different tectonic epochs.Therefore as shown in Figure 1, the present embodiment generation more phases
During fault activities map, step S130 is executed after executing step S110 and S120, determines that more phase fractures are lived
The quantity of activity intensity map is broken in cardon spectrum.Meanwhile in step s 130 also it needs to be determined that constituting more phase fault activities figures
Spectrum the corresponding tectonic epochs of each fault activities intensity map (usually in multiple continuous tectonic epochs same fault activities it is disconnected
It splits activity intensity map and constitutes more phase fault activities maps).
Next the fault activities intensity map for constructing more phase fault activities maps can be generated.In the present embodiment,
More phase fault activities maps are constructed by the way of arranging while generating (in an alternative embodiment of the invention, can also adopt
With first generating, all fault activities intensity maps are then unified to be arranged in a manner of constructing more phase fault activities maps).Meanwhile being
It changes with time trend convenient for characterization fault activities intensity, in the mistake of the fault activities intensity map of combination different configuration phase
Cheng Zhong arranges the fault activities intensity map of different configuration phase sequentially in time in the longitudinal direction.
Specifically, step S140 is first carried out, earliest/late-tectonic fault activities intensity map is generated.Then it executes
Earliest/late-tectonic fault activities intensity map is placed on more phases by step S150, assembled arrangement fault activities intensity map
The most lower/upper side of fault activities map.Next step S160 is executed, map quantity is verified, that is, verifies more phase fault activities figures
Whether the fault activities intensity map quantity placed in spectrum is for confirmation generates and has disposed all fault activities intensity maps
Spectrum.When step S160 judgement does not generate and disposes all fault activities intensity map time spectrums, execution step S140, raw again
At the fault activities intensity map of next (early one phase an of phase/evening) tectonic epochs.Step S150 is executed again, by newly-generated fracture
Activity intensity map is placed in the above/below of upper fault activities intensity map in more phase fault activities maps.
Step S160, S140 and S150 are repeated until completing to the fault activities intensity map of possessive construction phase
It generates to dispose.Step S170 can be executed at this time, generate the more phase fault activities maps being finally completed.
During disposing fault activities intensity map, for the ease of distinguishing the fault activities intensity of different tectonic epochs
Map retains the interval of specific range on more phase fault activities maps between two neighboring fault activities intensity map.
Further, fault activities intensity of the same statistics point in the different configuration phase is compared for convenience, in the present embodiment
In, during combining the fault activities intensity map of different configuration phase, using unified abscissa and keep all fractures
The one of activity intensity map abscissa ratio show guarantee multiple column blocks of longitudinal same abscissa it is corresponding be with unifying
Enumeration.
In addition, in the present embodiment, using unification for the fault activities intensity map of different tectonic epochs in construction
The colour code rule of correspondence.Same statistics point thus can further be compared in the fault activities intensity of different configuration phase.Equally,
In the embodiment of the height characterization fault activities intensity based on column block, for the ease of comparing same statistics point in different configuration
The fault activities intensity of phase uses unified ordinate ratio for the fault activities intensity map of different tectonic epochs.
By taking an application examples as an example, more phase fault activities maps for ultimately generating as shown in Fig. 2, bottom color bar 200
For colour code rule of correspondence table, 6 lateral chromatism colour bars of top are exactly the fault activities intensity map for respectively corresponding 6 different configuration phases
Spectrum.6 lateral chromatism colour bars respectively correspond the construction 1-6 phase from morning to night from bottom to top.
Next based on a specific application example come detailed description of the present invention effect.It is fractured into Tarim Basin
Example, as shown in figure 3, constructing the fracture in end of middle Ordovician Period, more phase fault activities of latest Ordovician and the Silurian Period in last three periods
Map (Fig. 3).As seen from Figure 3:
It is 100-200m, the activity at the middle part (at 60km) of fracture in the average turn-off of end of middle Ordovician Period fracture mass activity
Property is stronger, and turn-off is up to 500m；
It is broken in latest Ordovician substantially inactive；
In Silurian end fracture mass activity enhancing, average turn-off reaches 200-300m, and strong zone of action is by being broken
Near being transferred at 90km at the 60km in center to east side, fracture is in the least active (turn-off of end of middle Ordovician Period script at 90km
For 0m), become the most active (turn-off rises to 600m) to Silurian end.
It is possible thereby to prove, compared to the prior art, image-characterization methods of the invention can be with more specific intuitive table
Levy the concrete condition of fault activities intensity；Meanwhile present invention employs the modes of quantitatively characterizing, so that characterization result is more in line with
Practical geological condition.
While it is disclosed that embodiment content as above but described only to facilitate understanding the present invention and adopting
Embodiment is not intended to limit the invention.Method of the present invention can also have other various embodiments.Without departing substantially from
In the case where essence of the present invention, those skilled in the art make various corresponding changes or change in accordance with the present invention
Shape, but these corresponding changes or deformation all should belong to scope of protection of the claims of the invention.
Claims (5)
1. a kind of image-characterization methods of fault activities, which is characterized in that the method comprises the steps of:
It obtains fault activities to be analyzed and is walking upward different statistics points in the fault activities of specific one or multiple tectonic epochs
Intensity；
It calculates separately the acquisition different statistics points and walks upward position in the fault activities to be analyzed；
The fault activities intensity map of the corresponding tectonic epochs is generated based on the fault activities intensity and the position；And
For the same fault activities, longitudinal combination arranges multiple and different fault activities intensity maps to obtain more phases
Fault activities map, and during combining the fault activities intensity map of different configuration phase, more phase fractures are lived
Cardon composes the interval with specific range between the upper two neighboring fault activities intensity map；
Wherein, the fault activities intensity map is constituted by multiple column blocks are transversely arranged, in which:
The fault activities intensity map is walking upward distance as abscissa using fault activities to be analyzed, and the fault activities are strong
The length for spending map is the total length of the fault activities to be analyzed；
Each corresponding statistics point of the column block, position of the column block on the fault activities intensity map
The corresponding corresponding statistics point is walking upward position in the fault activities to be analyzed；And the fault activities intensity map
Using fault activities intensity value as the ordinate of the fault activities intensity map, the height of the column block corresponds to the statistics
The fault activities intensity of the point in the current tectonic epochs；Or
Determine the color of the column block based on the specific colour code rule of correspondence, the different colors corresponds to different described disconnected
The value of activity intensity is split, the color of the column block corresponds to the statistics point in the fault activities of the current tectonic epochs
Intensity.
2. the method according to claim 1, wherein during generating the fault activities intensity map, needle
Unified ordinate ratio is used to the different tectonic epochs.
3. the method according to claim 1, wherein during generating the fault activities intensity map, needle
To the different tectonic epochs using the unified colour code rule of correspondence.
4. the method according to claim 1, wherein in the fault activities intensity map of combination different configuration phase
During spectrum, using unified abscissa and the consistent of all fault activities intensity map abscissa ratios is kept.
5. according to the method described in claim 4, it is characterized in that, in the fault activities intensity map for combining the different configuration phase
During spectrum, the fault activities intensity map of different configuration phase is arranged sequentially in time in the longitudinal direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510543884.XA CN106483575B (en) | 2015-08-28 | 2015-08-28 | A kind of image-characterization methods of fault activities |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510543884.XA CN106483575B (en) | 2015-08-28 | 2015-08-28 | A kind of image-characterization methods of fault activities |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106483575A CN106483575A (en) | 2017-03-08 |
| CN106483575B true CN106483575B (en) | 2019-04-30 |
Family
ID=58235317
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510543884.XA Active CN106483575B (en) | 2015-08-28 | 2015-08-28 | A kind of image-characterization methods of fault activities |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106483575B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114647004B (en) * | 2022-02-25 | 2024-03-15 | 中国地质大学（北京） | Method for confirming sliding directions of underground sliding fracture at different periods |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103245972A (en) * | 2012-02-07 | 2013-08-14 | 中国石油天然气集团公司 | Method for determining complex geologic structure in two-dimensional space |
| CN103514630A (en) * | 2013-10-16 | 2014-01-15 | 北京石油化工学院 | Fault structure three-dimensional modeling method |
-
2015
- 2015-08-28 CN CN201510543884.XA patent/CN106483575B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103245972A (en) * | 2012-02-07 | 2013-08-14 | 中国石油天然气集团公司 | Method for determining complex geologic structure in two-dimensional space |
| CN103514630A (en) * | 2013-10-16 | 2014-01-15 | 北京石油化工学院 | Fault structure three-dimensional modeling method |
Non-Patent Citations (4)
| Title |
|---|
| The process zone: A microstructural view of fault growth;Jan M.Vermilye et al.;《JOURNAL OF GEOPHYSICAL RESEARCH》;19980610;第103卷(第B6期);第12223-12237页 |
| 一种断裂活动强度计算方法及其应用;卢异等;《天然气地球科学》;20100831;第21卷(第4期);第612-616页 |
| 徐深地区徐西深大断裂特征与火山岩天然气成藏条件;梁正中等;《天然气地球科学》;20090630;第20卷(第3期);第374-375页，图5 |
| 黄骅坳陷中区新生代断裂系统及其成因分析;何书等;《地质科学》;20080731;第43卷(第3期);第536-538页，图3 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106483575A (en) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pollastro | Total petroleum system assessment of undiscovered resources in the giant Barnett Shale continuous (unconventional) gas accumulation, Fort Worth Basin, Texas | |
| Cao et al. | The Permian hybrid petroleum system in the northwest margin of the Junggar Basin, northwest China | |
| Hu et al. | Crustal stress pattern in China and its adjacent areas | |
| CN106873032B (en) | A kind of granite hill reservoir crack prediction method | |
| Ellero et al. | Transpressional structuring of the High Atlas belt, Morocco | |
| Westman et al. | Passive seismic tomography for three-dimensional time-lapse imaging of mining-induced rock mass changes | |
| CN104267431A (en) | Method for obtaining minor fault for seismic tectonic interpretation | |
| CN106202673B (en) | Determine the method and device of oil reservoir moisture percentage and oily recovery percent of reserves relationship | |
| CN106338765B (en) | Method and system for determining the composite elastic parameter oil-gas recognition factor | |
| CN106483575B (en) | A kind of image-characterization methods of fault activities | |
| CN104181595A (en) | Novel method for quantitative recognition of fault associated cracks of complex tension structure system | |
| Palombo | To what extent could functional diversity be a useful tool in inferring ecosystem responses to past climate changes? | |
| CN106483579B (en) | A method of for obtaining fault activities situation | |
| Zhai et al. | A new tool for multi-cluster & multi-well hydraulic fracture modeling | |
| CN105986819A (en) | Method and device used for automatic processing and comprehensive interpretation of logging information | |
| CN113803044B (en) | Method and system for integrally designing unconventional reservoir volume fracturing and well distribution scheme | |
| Xiang et al. | Hydrocarbon migration and accumulation along the fault intersection zone—A case study on the reef-flat systems of the No. 1 slope break zone in the Tazhong area, Tarim Basin | |
| Manchanda et al. | An integrated approach to development optimization using monitor wells and hydraulic fracture diagnostics in the Permian Basin | |
| CN106483560B (en) | A method of obtaining the normal throw of fracture | |
| CN104298856A (en) | Tunnel advance geology forecast method based on surrounding rock deformation monitoring and numerical simulation | |
| Joseph et al. | An Investigation of Candidate Mechanisms for Hydraulic Fracture Swarming through High-Fidelity Numerical Modeling | |
| Moradi et al. | Learning from Behavioral Frac Maps: A Montney Case Study in Integration of Modern Microseismic and Production Data Analyses | |
| Zhu et al. | Hydraulic fracture propagation in Denver-Julesburg basin constrained by cross-well distributed strain measurements | |
| Lang et al. | Macro-characterisation of cave damage for palaeoseismological investigations in regions of low strain: A case study from central–western North Island (Waitomo caves), New Zealand | |
| Shoemaker et al. | Integrating 3D seismic and geomechanical properties with microseismic acquisition and fracturing parameters to optimize completion practices within the Wolfcamp Shale Play of the Midland Basin |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |