CN111060968B - 一种用于球面扩散补偿精确井控速度校正方法 - Google Patents
一种用于球面扩散补偿精确井控速度校正方法 Download PDFInfo
- Publication number
- CN111060968B CN111060968B CN201911350090.6A CN201911350090A CN111060968B CN 111060968 B CN111060968 B CN 111060968B CN 201911350090 A CN201911350090 A CN 201911350090A CN 111060968 B CN111060968 B CN 111060968B
- Authority
- CN
- China
- Prior art keywords
- speed
- well
- velocity
- rms
- points
- 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
- 238000009792 diffusion process Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000002159 abnormal effect Effects 0.000 claims abstract description 46
- 230000008859 change Effects 0.000 claims abstract description 24
- 238000012937 correction Methods 0.000 claims abstract description 9
- 238000012216 screening Methods 0.000 claims abstract description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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)
Abstract
本发明提出了一种用于球面扩散补偿精确井控速度校正方法，步骤S1，对全区的测井资料进行筛选，并根据VSP速度或者根据测井曲线生成RMS速度，找到RMS速度的异常点；步骤S2，将异常点代入到全区的井速度模型，转为含有异常速度的地层速度；步骤S3，利用井速度模型，约束地震速度的拾取，对异常点进行校正，形成井控速度场；步骤S4，将井控速度场，应用到球面扩散补偿的程序中。本发明采用补偿因子和井控速度场作为全区振幅补偿的参数，这样能够体现空间变化；能够在复杂地区的使用，其应用范围广。
Description
技术领域
本发明涉及井控补偿技术领域，特别涉及一种用于球面扩散补偿精确井控速度校正方法。
背景技术
地震信号在向地下传播的过程中，能量有一定的衰减，横向炮与炮之间能量也存在差别。影响地震振幅的因素大致有两种：一是随着传播距离的增大，地震反射信号能量迅速降低，使得深、浅层反射能量差异悬殊，不利于地震成像和地质分析，显然这种影响是非地质因素造成的，需要进行恢复或补偿，这类因素有球面扩散、吸收衰减等，它们主要造成纵向上的能量差异，其影响是时间的函数；二是由于激发、接收条件的差异对反射所造成的振幅影响，它与空间位置有密切的关系，这种影响不利于成像和储层预测，需要予以排除，它们主要造成横向上的能量差异，是空间位置的函数。因此，为了改善成像和储层预测的效果，需要对记录的地震信号应用时变或空变的增益函数g(t)合理补偿不同因素引起的振幅衰减和差异。
不同的振幅补偿方法，对应不同的增益函数。常规的振幅补偿方法，往往是通过试验来选择参数，如球面扩散补偿系数扫描，通过分析补偿前后的单炮记录及叠加剖面的能量变化，确定最佳的补偿参数。由于地震信号在均匀介质中传播时，能量衰减正比于，(r为球面半径)；地震波的振幅正比于能量的平方根，它随1/r衰减。实际上地层速度一般随深度增加，这使地震波更进一步发散，振幅也随距离衰减变快。常用补偿函数为T补偿函数，power是T补偿参数，以及基于速度的球面扩散补偿，一般是应用工区的区域速度做振幅补偿，当工区的构造起伏大，地下构造的横向变化比较大,不能使用统一的速度函数时，常规的振幅补偿方法对振幅保真是不利的。
早期井控振幅补偿方法根据VSP资料求取的确定性TAR因子进行振幅恢复，该项技术能够实现纵向能量的合理补偿。对于构造起伏较大的三维地震数据来说，早期井控技术不能体现地层的吸收、弱响应和地层粘性特征；不能体现地表变化；不能体现空间变化；简单使用统一的速度进行振幅补偿，对振幅的保真性是不利的。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种用于球面扩散补偿精确井控速度校正方法，采用补偿因子和井控速度场作为全区振幅补偿的参数，这样能够体现空间变化；能够在复杂地区的使用，其应用范围广，还能够体现地层的吸收、弱响应和地层粘性特征。
为了实现上述目的，本发明提供一种用于球面扩散补偿精确井控速度校正方法，包括以下步骤：
步骤S1，对全区的测井资料进行筛选，并根据VSP速度或者根据测井曲线生成RMS速度，找到RMS速度的异常点；
在工区构造的高点和低点，选择所有的井资料，获得其RMS速度，并查找RMS速度的变化规律，对于变化点做标记，并利用相邻井资料进行核实；这些核实后的变化点就是RMS速度的异常点；
步骤S2，将异常点代入到全区的井速度模型，转为含有异常速度的地层速度；
步骤S3，利用井速度模型，约束地震速度的拾取，对异常点进行校正，形成井控速度场；
将异常点对应的含有异常速度的地层速度，与该地层速度的相同地层相同时间同一位置所对应的正常速度进行对比，并计算其速度差，当速度差大于一定值时，则对该异常点所对应的RMS速度进行调整，直到调整后的异常点的速度差不大于一定值时为止，并将调整后的异常点与正常点形成井控速度场；
步骤S4，将井控速度场，应用到球面扩散补偿的程序中。
在上述任一方案中优选的是，在步骤S1中，采用VSP技术获取全区的测井资料时，首先将零偏和非零偏校正到统一基准面，然后在初至拾取前确定资料的极性，接着按不同震源类型拾取正确的初至位置，在拾取初至时，根据VSP速度生成RMS速度，然后根据RMS速度变化规律，找到RMS速度的异常点。
在上述任一方案中优选的是，零偏包括直井标定，直井标定是直接利用走廊叠加和大井位置的地面地震剖面进行波组特征匹配，再利用动校道集确定时深关系。
在上述任一方案中优选的是，零偏还包括斜井标定，在进行斜井标定时，如井斜很小，则进行走廊叠加，然后按照井斜对走廊进行轨迹校正，再进行标定；如井斜很大，则进行零偏成像处理，利用零偏成像进行标定。
在上述任一方案中优选的是，在步骤S1中，根据测井曲线计算测井曲线上每个RMS速度点斜率，找到与其他相邻RMS速度点斜率差别大的RMS速度点作为变化点，并对其进行标记，并利用相邻井资料进行核实；这些核实后的变化点就是RMS速度的异常点。
在上述任一方案中优选的是，步骤S2中，通过克里金法插值，获得全区的井速度模型，将异常点到全区的井速度模型中，转为含有异常速度的地层速度，并被记录下来。
在上述任一方案中优选的是，全区的井速度模型公式为：
其中，其中s为距离，t为传播时间，V为该点的速度。
在上述任一方案中优选的是，步骤S3中，速度差的计算公式为：
其中，ΔV为速度差；Vavg为正常速度；Vk为第k个点的速度。
在上述任一方案中优选的是，进行球面扩散补偿时，首先，根据VSP资料求取的确定性补偿因子，然后将补偿因子和井控速度场代入到球面扩散补偿的振幅计算公式中，进行振幅恢复。
本发明的用于球面扩散补偿精确井控速度校正方法具有以下有益效果：
1、本发明采用改进型的井控速度球面扩散补偿技术，通过利用工区内多种井资料(VSP、测井曲线等)，多个井点位置速度形成井控速度场，约束地震速度的拾取，从而形成井控约束下的全区速度场，将此速度用于全区振幅补偿，将会取得较好的效果。这种基于井速度约束的修正速度更准确，时、空变的速度更灵活适应复杂构造，更高阶的补偿函数实现更高精度的补偿。
2、本发明将异常点对应的含有异常速度的地层速度，与该地层速度的相同地层相同时间同一位置所对应的正常速度进行对比，并计算其速度差，当速度差大于一定值时，则对该异常点所对应的RMS速度进行调整，从而约束地震速度拾取。
3、本发明采用补偿因子和井控速度场作为全区振幅补偿的参数，这样能够体现空间变化；能够在复杂地区的使用，其应用范围广，还能够体现地层的吸收、弱响应和地层粘性特征。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是本发明的用于球面扩散补偿精确井控速度校正方法的流程图；
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
VSP是Vertical Seismic Profiling的缩写，中文为垂直地震剖面法，地面地震通常是将震源和检波器都置于地面进行采集；而VSP技术是将震源和检波器其中的一种置于井下进行的地震采集的技术。
常规球面扩散振幅补偿方法往往是使用一个最佳的补偿参数，基于速度的球面扩散补偿，一般是应用工区的区域速度做振幅补偿，对于构造起伏大，地下构造的横向变化比较大的三维地震数据来说，常规补偿方法对振幅保真是不利的，不能满足反射层及层间细层的振幅恢复需求。
基于此，本发明提供一种用于球面扩散补偿精确井控速度校正方法，如图1所示，包括以下步骤：
步骤S1，对全区的测井资料进行筛选，并根据VSP速度或者根据测井曲线生成RMS速度，找到RMS速度的异常点；
在工区构造的高点和低点，选择所有的井资料，获得其RMS速度，并查找RMS速度的变化规律，对于变化点做标记，并利用相邻井资料进行核实；这些核实后的变化点就是RMS速度的异常点；
采用VSP技术获取全区的测井资料时，首先将零偏和非零偏校正到统一基准面，然后在初至拾取前确定资料的极性，接着按不同震源类型拾取正确的初至位置，在拾取初至时，根据VSP速度生成RMS速度，然后根据RMS速度变化规律，找到RMS速度的异常点。
零偏包括直井标定，直井标定是直接利用走廊叠加和大井位置的地面地震剖面进行波组特征匹配，再利用动校道集确定时深关系。
零偏还包括斜井标定，在进行斜井标定时，如井斜很小，则进行走廊叠加，然后按照井斜对走廊进行轨迹校正，再进行标定；如井斜很大，则进行零偏成像处理，利用零偏成像进行标定。
采用测井曲线获取全区的测井资料时，根据测井曲线计算测井曲线上每个RMS速度点斜率，找到与其他相邻RMS速度点斜率差别大的RMS速度点作为变化点，并对其进行标记，并利用相邻井资料进行核实；这些核实后的变化点就是RMS速度的异常点。
步骤S2，将异常点代入到全区的井速度模型，转为含有异常速度的地层速度；
通过克里金法插值，获得全区的井速度模型，将异常点到全区的井速度模型中，转为含有异常速度的地层速度，并被记录下来。
全区的井速度模型公式为：
其中，其中s为距离，t为传播时间，V为该点的速度。
步骤S3，利用井速度模型，约束地震速度的拾取，对异常点进行校正，形成井控速度场；
将异常点对应的含有异常速度的地层速度，与该地层速度的相同地层相同时间同一位置所对应的正常速度进行对比，并计算其速度差，当速度差大于一定值时，则对该异常点所对应的RMS速度进行调整，直到调整后的异常点的速度差不大于一定值时为止，并将调整后的异常点与正常点形成井控速度场；
速度差的计算公式为：
其中，ΔV为速度差；Vavg为正常速度；Vk为第k个点的速度。
本发明的另一个实施例中，一定值优选1％，即ΔV的绝对值大于1％时，则对该异常点所对应的RMS速度进行调整。
步骤S4，将井控速度场，应用到球面扩散补偿的程序中。
进行球面扩散补偿时，首先，根据VSP资料求取的确定性补偿因子，然后将补偿因子和井控速度场代入到球面扩散补偿的振幅计算公式中，进行振幅恢复。
对于构造起伏较大的三维地震数据来说，简单使用统一的补偿因子值进行振幅补偿，对振幅的保真性是不利的。本发明采用改进型的井控速度球面扩散补偿技术，通过利用工区内多种井资料(VSP、测井曲线等)，多个井点位置速度形成井控速度场，约束地震速度的拾取，从而形成井控约束下的全区速度场，将此速度用于全区振幅补偿，将会取得较好的效果。这种基于井速度约束的修正速度更准确，时、空变的速度更灵活适应复杂构造，更高阶的补偿函数实现更高精度的补偿。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (7)
1.一种用于球面扩散补偿精确井控速度校正方法，其特征在于，包括以下步骤：
步骤S1，对全区的测井资料进行筛选，并根据VSP速度或者根据测井曲线生成RMS速度，找到RMS速度的异常点；
在工区构造的高点和低点，选择所有的井资料，获得其RMS速度，并查找RMS速度的变化规律，对于变化点做标记，并利用相邻井资料进行核实；这些核实后的变化点就是RMS速度的异常点；
(1)当根据VSP速度生成RMS速度找到RMS速度的异常点时：
采用VSP技术获取全区的测井资料时，首先将零偏和非零偏校正到统一基准面，然后在初至拾取前确定资料的极性，接着按不同震源类型拾取正确的初至位置，在拾取初至时，根据VSP速度生成RMS速度，然后根据RMS速度变化规律，找到RMS速度的异常点；
(2)当根据测井曲线生成RMS速度找到RMS速度的异常点时：
根据测井曲线计算测井曲线上每个RMS速度点斜率，找到与其他相邻RMS速度点斜率差别大的RMS速度点作为变化点，并对其进行标记，并利用相邻井资料进行核实；这些核实后的变化点就是RMS速度的异常点；
步骤S2，将异常点代入到全区的井速度模型，转为含有异常速度的地层速度；
步骤S3，利用井速度模型，约束地震速度的拾取，对异常点进行校正，形成井控速度场；
将异常点对应的含有异常速度的地层速度，与该地层速度的相同地层相同时间同一位置所对应的正常速度进行对比，并计算其速度差，当速度差大于一定值时，则对该异常点所对应的RMS速度进行调整，直到调整后的异常点的速度差不大于一定值时为止，并将调整后的异常点与正常点形成井控速度场；
步骤S4，将井控速度场，应用到球面扩散补偿的程序中。
2.如权利要求1所述的用于球面扩散补偿精确井控速度校正方法，其特征在于，零偏包括直井标定，直井标定是直接利用走廊叠加和大井位置的地面地震剖面进行波组特征匹配，再利用动校道集确定时深关系。
3.如权利要求1或2所述的用于球面扩散补偿精确井控速度校正方法，其特征在于，零偏还包括斜井标定，在进行斜井标定时，如井斜小，则进行走廊叠加，然后按照井斜对走廊进行轨迹校正，再进行标定；如井斜大，则进行零偏成像处理，利用零偏成像进行标定。
4.如权利要求1所述的用于球面扩散补偿精确井控速度校正方法，其特征在于，步骤S2中，通过克里金法插值，获得全区的井速度模型，将异常点到全区的井速度模型中，转为含有异常速度的地层速度，并被记录下来。
5.如权利要求1或4所述的用于球面扩散补偿精确井控速度校正方法，其特征在于，全区的井速度模型公式为：
其中，其中s为距离，t为传播时间，V为该点的速度。
6.如权利要求1所述的用于球面扩散补偿精确井控速度校正方法，其特征在于，步骤S3中，速度差的计算公式为：
其中，ΔV为速度差；Vavg为正常速度；Vk为第k个点的速度。
7.如权利要求1所述的用于球面扩散补偿精确井控速度校正方法，其特征在于，进行球面扩散补偿时，首先，根据VSP资料求取的确定性补偿因子，然后将补偿因子和井控速度场代入到球面扩散补偿的振幅计算公式中，进行振幅恢复。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911350090.6A CN111060968B (zh) | 2019-12-24 | 2019-12-24 | 一种用于球面扩散补偿精确井控速度校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911350090.6A CN111060968B (zh) | 2019-12-24 | 2019-12-24 | 一种用于球面扩散补偿精确井控速度校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111060968A CN111060968A (zh) | 2020-04-24 |
| CN111060968B true CN111060968B (zh) | 2022-10-21 |
Family
ID=70303242
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911350090.6A Active CN111060968B (zh) | 2019-12-24 | 2019-12-24 | 一种用于球面扩散补偿精确井控速度校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111060968B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114964042B (zh) * | 2022-05-20 | 2023-10-20 | 西安交通大学 | 一种曲线轮廓在线测量中数据异常点判别与辨识方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102998700A (zh) * | 2011-09-08 | 2013-03-27 | 中国石油天然气集团公司 | 一种垂直地震剖面资料的地表一致性能量补偿方法 |
| CN103089306A (zh) * | 2011-11-03 | 2013-05-08 | 大港油田集团有限责任公司 | 利用枯竭式油气藏建造地下储气库的方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101738638B (zh) * | 2008-11-24 | 2011-10-05 | 中国石油集团东方地球物理勘探有限责任公司 | 石油地震构造图等值线的填充方法 |
| US8942064B2 (en) * | 2009-06-10 | 2015-01-27 | Baker Hughes Incorporated | Sending a seismic trace to surface after a vertical seismic profiling while drilling measurement |
| CN102943669B (zh) * | 2012-10-09 | 2015-03-11 | 中国石油化工股份有限公司 | 适用于水平井轨迹设计的高精度深度域随机模拟反演方法 |
| CN104914463B (zh) * | 2014-03-14 | 2018-04-27 | 恒泰艾普集团股份有限公司 | 小尺度裂缝预测的成像方法 |
-
2019
- 2019-12-24 CN CN201911350090.6A patent/CN111060968B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102998700A (zh) * | 2011-09-08 | 2013-03-27 | 中国石油天然气集团公司 | 一种垂直地震剖面资料的地表一致性能量补偿方法 |
| CN103089306A (zh) * | 2011-11-03 | 2013-05-08 | 大港油田集团有限责任公司 | 利用枯竭式油气藏建造地下储气库的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111060968A (zh) | 2020-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
| US6917564B2 (en) | Method of processing seismic data | |
| US5999489A (en) | High vertical resolution crosswell seismic imaging | |
| US6201765B1 (en) | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data | |
| CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN108196305B (zh) | 一种山地静校正方法 | |
| AU2001275708A1 (en) | A method of processing seismic data | |
| WO2006054181A1 (en) | Method for processing at least two sets of seismic data | |
| CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
| CN104483705A (zh) | 一种三维剩余静校正方法 | |
| US20030187583A1 (en) | Method and apparatus for resolving shear wave seismic data | |
| US6311131B1 (en) | Method of and apparatus for determining the quality of seismic data | |
| CN111060968B (zh) | 一种用于球面扩散补偿精确井控速度校正方法 | |
| EP0288154B1 (en) | Method for estimating shear wave reflection data from acquired compressional wave reflection data | |
| CN113568046B (zh) | 裂缝方位约束的多维信息裂缝表征方法 | |
| CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
| CN116500679A (zh) | 一种井震标定方法、装置和相关设备 | |
| CN109839659B (zh) | 一种对叠前深度偏移剖面进行迭代优化的方法 | |
| CN110967751B (zh) | 基于地面浅井监测的微地震事件的定位方法及存储介质 | |
| CN112305604A (zh) | 一种可控震源组合激发条件下的近地表物性区分方法 | |
| WO2023123971A1 (zh) | 基于vsp的深度域地震剖面层位标定方法及装置 | |
| CN111337975B (zh) | 一种基于方差分形维数的微地震事件自动识别方法 | |
| CN104316967B (zh) | 大偏移距动校正速度时移误差补偿方法 | |
| CN117991352A (zh) | 基于纵波avo截距减梯度属性的多波数据匹配方法 |
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 |