CN113719280A - 一种低阻气层测井识别方法 - Google Patents
一种低阻气层测井识别方法 Download PDFInfo
- Publication number
- CN113719280A CN113719280A CN202111052085.4A CN202111052085A CN113719280A CN 113719280 A CN113719280 A CN 113719280A CN 202111052085 A CN202111052085 A CN 202111052085A CN 113719280 A CN113719280 A CN 113719280A
- Authority
- CN
- China
- Prior art keywords
- low
- logging
- gas
- phi
- gas layer
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 24
- 238000011161 development Methods 0.000 claims abstract description 16
- 239000011435 rock Substances 0.000 claims abstract description 16
- 238000005553 drilling Methods 0.000 claims abstract description 8
- 230000005611 electricity Effects 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 239000012530 fluid Substances 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 claims description 8
- 230000000704 physical effect Effects 0.000 claims description 8
- 238000011160 research Methods 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 8
- 230000009897 systematic effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000008021 deposition Effects 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 239000007789 gas Substances 0.000 abstract description 90
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 12
- 239000003345 natural gas Substances 0.000 abstract description 6
- 230000004888 barrier function Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于天然气勘探开发、地球物理测井技术领域，尤其是一种低阻气层测井识别方法，针对现有技术大都是运用单一、特定的方法对低阻气层进行识别，并不能完全的将低阻气层有效识别出来，对低阻气层的解释精度明显较低的问题，现提出如下方案，其包括以下步骤：S1：充分收集钻井取心、试气以及测井资料，如有必要，对测井曲线进行标准化；S2：岩电归位；S3：“四性”关系研究；S4：拾取低阻气层发育段；S5：计算三孔隙度值ФS、ФD、ФN、等效弹性模量差比值DR、计算声阻抗比值Z。本发明能对低阻气层进行识别的方法，操作方式灵活，能实现全井位、全油田广泛覆盖。
Description
技术领域
本发明涉及天然气勘探开发、地球物理测井技术领域，尤其涉及一种低阻气层测井识别方法。
背景技术
在天然气勘探过程当中，一般天然气层位都表现出电阻率为高值的特征，但是随着勘探程度的深入，遇到低阻气层(相对低电阻率气层)也具有较高产量，低阻气层指含有一定商业价值的烃类流体，但测井电阻率值较低的产气层，一般认为如果在同一气水系统内气层与纯水层的电阻率比值小于2，即电阻增大率小于2，这种气层就是低阻气层。低阻气层对于气田储量的意义越来越大。准确的识别和评价低阻气层对于油田储量升级以及开发方案的制定都具有重要意义。但是，由于低阻地层成因复杂，每个地区低阻气层的主控因素都不完全一致，导致目前提出的方法都具有针对性和局限性，都只适用于特定区块。
现有技术大都是运用单一、特定的方法对低阻气层进行识别，并不能完全的将低阻气层有效识别出来，对低阻气层的解释精度明显较低，有较大一批气层会被解释为水层，不利于气层的勘探开发，没办法实现全井位、全油田广泛覆盖。
发明内容
本发明的目的是为了解决现有技术大都是运用单一、特定的方法对低阻气层进行识别，并不能完全的将低阻气层有效识别出来，对低阻气层的解释精度明显较低，有较大一批气层会被解释为水层，不利于气层的勘探开发，没办法实现全井位、全油田广泛覆盖的缺点，而提出的一种低阻气层测井识别方法。
为了实现上述目的，本发明采用了如下技术方案：
一种低阻气层测井识别方法，包括以下步骤：
S1：充分收集钻井取心、试气以及测井资料，如有必要，对测井曲线进行标准化；
S2：岩电归位；
S3：“四性”关系研究；即研究岩性、物性、含气性与电性之间的关系，分析低阻气层测井响应特征
S4：拾取低阻气层发育段；
S5：计算三孔隙度值ФS、ФD、ФN、等效弹性模量差比值DR、计算声阻抗比值Z，将等效弹性模量差比值DR和声阻抗Z在一个曲线道内对称显示；
S6：综合识别气层，根据试气、试产资料，计算符合率。
优选的，所述S1中充分收集钻井取心、试气以及测井资料，如有必要，对测井曲线进行标准化；测井曲线可以被看作是有效信号、随机噪音以及系统误差的综合，测井曲线标准化就是消除测井曲线系统误差的过程。
优选的，所述S2中岩电归位，即将岩心深度与电测深度对应，将其恢复到真实深度的工作，整理和加载岩心数据时，必须对照电测资料，加上地质人员的判断，同时要校正钻具长度和测井深度上的系统误差，要将岩心的不同岩性与电测曲线解释的岩性仔细对应，恢复岩心所在真实深度。
优选的，所述S3中“四性”关系研究。即研究岩性、物性、含气性与电性之间的关系，分析低阻气层测井响应特征。测井曲线是储层岩性、物性、流体性质等因素的共同反映，同时受控于构造、沉积环境、地层水的变化等，并受储层非均质性的影响。示例低阻气层在测井曲线上表现为：低伽马10-50API、低岩性-密度指数PE<2.5巴/电子、高声波时差△t>230μs/m、低电阻率Rt<40Ω·m、低密度ρ<2.45g/cm3等特征。
优选的，所述S4中拾取低阻气层发育段，根据S3中总结出的该研究地区低阻气层测井响应特征，对应试气资料，拾取研究区内低阻气层发育段。
优选的，所述S5中计算三孔隙度值ФS、ФD、ФN、计算等效弹性模量差比值DR、计算声阻抗比值Z，将三孔隙度值转换成同一刻度后在一个曲线道显示，等效弹性模量差比值DR和声阻抗Z在一个曲线道内对称显示。
三孔隙度值计算公式：
式中，ФS、ФD、ФN分别为声波时差孔隙度、密度孔隙度和中子孔隙度。
等效弹性模量差比值计算公式为：
式中：DR为等效弹性模量差比值；M0为目的层饱含水时岩石的等效弹性模量；M为目的层岩石的等效弹性模量。
式中：
为纵波声波时差，μs/m；ρb为岩石体积密度，g/cm3。
声阻抗比值计算公式为：
式中：T*为流体时差比；P*为流体密度比。
将三孔隙度值ФS、ФD、ФN转换成同一刻度后在一个曲线道显示，等效弹性模量差比值DR和声阻抗比值Z在一个曲线道内对称显示。可以将实例中的低阻气层识别出来。
优选的，所述S6中综合识别气层，根据试气、试产资料，计算符合率，自然电位负异常幅度很大，自然伽马值低，曲线呈箱状，PE指数值较小，高声波时差，低密度，电阻率值较低<40Ω·m，声波时差和自然电位重叠图以及视自然电位与自然电位重叠图均有幅度差，三孔隙度重叠明显有ФS>ФD>ФN，△RM>0，故将其解释为低阻气层，如此一来，就可以将原本油田解释为气水同层层5、层6的低阻气层正确解释出来，经试气结果证实确为气层，两层日产气量合计为6534m3，不产水。
与现有技术相比，本发明的优点在于：
(1)本方案综合利用三孔隙度差值法、纵波等效弹性模量差比值法以及声阻抗比值法，对低阻气层进行识别。
(2)本发明可以根据油田实际情况做各种更改和变化，对于低阻气层的识别这一技术难点，通过优选非电阻率参数，综合运用三孔隙度差值法、纵波等效弹性模量差比值法以及声阻抗比值法，对低阻气层进行有效识别，提高气层识别的准确定，减少有效气层漏失，为计算天然气储量和后期勘探开发提供技术支撑；
(3)操作方式灵活，能实现全井位、全油田广泛覆盖。
附图说明
图1为本发明提出的低阻气层识别流程图；
图2为本发明提出的某井低阻气层解释成果图。
具体实施方式
下面将结合本实施例中的附图，对本实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本实施例一部分实施例，而不是全部的实施例。
参照图1-2，一种低阻气层测井识别方法，包括以下步骤：
S1：充分收集钻井取心、试气以及测井资料，如有必要，对测井曲线进行标准化；
S2：岩电归位；
S3：“四性”关系研究；即研究岩性、物性、含气性与电性之间的关系，分析低阻气层测井响应特征
S4：拾取低阻气层发育段；
S5：计算三孔隙度值ФS、ФD、ФN、等效弹性模量差比值DR、计算声阻抗比值Z。将等效弹性模量差比值DR和声阻抗Z在一个曲线道内对称显示；
S6：综合识别气层，根据试气、试产资料，计算符合率。
本实施例，S1中充分收集钻井取心、试气以及测井资料，如有必要，对测井曲线进行标准化；测井曲线可以被看作是有效信号、随机噪音以及系统误差的综合，测井曲线标准化就是消除测井曲线系统误差的过程。
本实施例，S2中岩电归位，即将岩心深度与电测深度对应，将其恢复到真实深度的工作，整理和加载岩心数据时，必须对照电测资料，加上地质人员的判断，同时要校正钻具长度和测井深度上的系统误差，要将岩心的不同岩性与电测曲线解释的岩性仔细对应，恢复岩心所在真实深度。
本实施例，S3中“四性”关系研究。即研究岩性、物性、含气性与电性之间的关系，分析低阻气层测井响应特征。测井曲线是储层岩性、物性、流体性质等因素的共同反映，同时受控于构造、沉积环境、地层水的变化等，并受储层非均质性的影响。如图2所示，示例低阻气层在测井曲线上表现为：低伽马10-50API、低岩性-密度指数PE<2.5巴/电子、高声波时差△t>230μs/m、低电阻率Rt<40Ω·m、低密度ρ<2.45g/cm3等特征。
本实施例，S4中拾取低阻气层发育段，根据S3中总结出的该研究地区低阻气层测井响应特征，对应试气资料，拾取研究区内低阻气层发育段。
本实施例，S5中计算三孔隙度值ФS、ФD、ФN、计算等效弹性模量差比值DR、计算声阻抗比值Z，将三孔隙度值转换成同一刻度后在一个曲线道显示，等效弹性模量差比值DR和声阻抗Z在一个曲线道内对称显示。
三孔隙度值计算公式：
式中，ФS、ФD、ФN分别为声波时差孔隙度、密度孔隙度和中子孔隙度。
等效弹性模量差比值计算公式为：
式中：DR为等效弹性模量差比值；M0为目的层饱含水时岩石的等效弹性模量；M为目的层岩石的等效弹性模量。
式中：
为纵波声波时差，μs/m；ρb为岩石体积密度，g/cm3。
声阻抗比值计算公式为：
式中：T*为流体时差比；P*为流体密度比。
如图2所示，将三孔隙度值ФS、ФD、ФN转换成同一刻度后在一个曲线道显示，等效弹性模量差比值DR和声阻抗比值Z在一个曲线道内对称显示，可以将实例中的低阻气层识别出来。
本实施例，S6中综合识别气层，根据试气、试产资料，计算符合率，图2为某井实际资料解释成果图，图2中，层3、层5、层6对应的曲线特征为：自然电位负异常幅度很大，自然伽马值低，曲线呈箱状，PE指数值较小，高声波时差，低密度，电阻率值较低<40Ω·m，声波时差和自然电位重叠图以及视自然电位与自然电位重叠图均有幅度差，三孔隙度重叠明显有ФS>ФD>ФN，△RM>0，故将其解释为低阻气层，如此一来，就可以将原本油田解释为气水同层层5、层6的低阻气层正确解释出来，经试气结果证实确为气层，两层日产气量合计为6534m3，不产水。
本实施例，考虑到电阻率参数已经不能对低阻气层有很好的反映，在原有的反映渗透性的声波时差曲线和自然电位曲线重叠图以及反映可动烃的视自然电位曲线与自然电位曲线重叠图的基础上，抛开电阻率信号，应用非电阻率测井数据，通过对数据的转换、整合和放大含气信息，综合利用三孔隙度差值法、纵波等效弹性模量差比值法以及声阻抗比值法，对低阻气层进行识别。
其原理是：当地层中含天然气时，与地层中含水相比或者，声波时差测井值增大，所以声波孔隙度增大ФS；补偿中子测井值减小，即补偿中子孔隙度减小ФD；补偿密度测井值减小，即密度孔隙度增大ФN。因而气层等效弹性模量值减小、声阻抗减小。
以上所述，仅为本实施例较佳的具体实施方式，但本实施例的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本实施例揭露的技术范围内，根据本实施例的技术方案及其发明构思加以等同替换或改变，都应涵盖在本实施例的保护范围之内。
Claims (7)
1.一种低阻气层测井识别方法，其特征在于，包括以下步骤：
S1：充分收集钻井取心、试气以及测井资料，对测井曲线进行标准化；
S2：岩电归位；
S3：“四性”关系研究；
S4：拾取低阻气层发育段；
S5：计算三孔隙度值ФS、ФD、ФN、等效弹性模量差比值DR、计算声阻抗比值Z；
S6：综合识别气层，根据试气、试产资料，计算符合率。
2.根据权利要求1所述的一种低阻气层测井识别方法，其特征在于，所述S1中充分收集钻井取心、试气以及测井资料，对测井曲线进行标准化；测井曲线可以被看作是有效信号、随机噪音以及系统误差的综合，测井曲线标准化就是消除测井曲线系统误差的过程。
3.根据权利要求1所述的一种低阻气层测井识别方法，其特征在于，所述S2中岩电归位，将岩心深度与电测深度对应，将其恢复到真实深度的工作。
4.根据权利要求1所述的一种低阻气层测井识别方法，其特征在于，所述S3中“四性”关系研究，即研究岩性、物性、含气性与电性之间的关系，分析低阻气层测井响应特征，测井曲线是储层岩性、物性、流体性质的共同反映，同时受控于构造、沉积环境、地层水的变化，并受储层非均质性的影响。
5.根据权利要求1所述的一种低阻气层测井识别方法，其特征在于，所述S4中拾取低阻气层发育段，根据S3中总结出的该研究地区低阻气层测井响应特征，对应试气资料，拾取研究区内低阻气层发育段。
6.根据权利要求1所述的一种低阻气层测井识别方法，其特征在于，所述S5中计算三孔隙度值ФS、ФD、ФN、计算等效弹性模量差比值DR、计算声阻抗比值Z，将三孔隙度值转换成同一刻度后在一个曲线道显示，等效弹性模量差比值DR和声阻抗Z在一个曲线道内对称显示。
三孔隙度值计算公式：
式中，ФS、ФD、ФN分别为声波时差孔隙度、密度孔隙度和中子孔隙度。
等效弹性模量差比值计算公式为：
式中：DR为等效弹性模量差比值；M0为目的层饱含水时岩石的等效弹性模量；M为目的层岩石的等效弹性模量。
式中：
为纵波声波时差，μs/m；ρb为岩石体积密度，g/cm3。
声阻抗比值计算公式为：
式中：T*为流体时差比；P*为流体密度比。
将三孔隙度值ФS、ФD、ФN转换成同一刻度后在一个曲线道显示，等效弹性模量差比值DR和声阻抗比值Z在一个曲线道内对称显示，将低阻气层识别出来。
7.根据权利要求1所述的一种低阻气层测井识别方法，其特征在于，所述S6中综合识别气层，根据试气、试产资料，计算符合率，自然电位负异常幅度很大，自然伽马值低，曲线呈箱状，PE指数值较小，高声波时差，低密度，电阻率值较低<40Ω·m，声波时差和自然电位重叠图以及视自然电位与自然电位重叠图均有幅度差，三孔隙度重叠明显有ФS>ФD>ФN，△RM>0，故将其解释为低阻气层。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111052085.4A CN113719280A (zh) | 2021-09-08 | 2021-09-08 | 一种低阻气层测井识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111052085.4A CN113719280A (zh) | 2021-09-08 | 2021-09-08 | 一种低阻气层测井识别方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113719280A true CN113719280A (zh) | 2021-11-30 |
Family
ID=78682680
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111052085.4A Pending CN113719280A (zh) | 2021-09-08 | 2021-09-08 | 一种低阻气层测井识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113719280A (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120109603A1 (en) * | 2009-06-22 | 2012-05-03 | Ning Li | Quantitative calculation method for oil (gas) saturation of fractured reservoir |
| US20140379265A1 (en) * | 2011-10-24 | 2014-12-25 | Repsol, S.A. | Real-time method for determining the porosity and water saturation of an underground formation using gas and mud logging data |
| CN104453851A (zh) * | 2014-11-26 | 2015-03-25 | 中国石油天然气股份有限公司 | 油田驱替物、油田驱替方法与油田驱替物的试验制备装置 |
| CN111271057A (zh) * | 2020-03-25 | 2020-06-12 | 西北大学 | 一种测井识别气层方法 |
-
2021
- 2021-09-08 CN CN202111052085.4A patent/CN113719280A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120109603A1 (en) * | 2009-06-22 | 2012-05-03 | Ning Li | Quantitative calculation method for oil (gas) saturation of fractured reservoir |
| US20140379265A1 (en) * | 2011-10-24 | 2014-12-25 | Repsol, S.A. | Real-time method for determining the porosity and water saturation of an underground formation using gas and mud logging data |
| CN104453851A (zh) * | 2014-11-26 | 2015-03-25 | 中国石油天然气股份有限公司 | 油田驱替物、油田驱替方法与油田驱替物的试验制备装置 |
| CN111271057A (zh) * | 2020-03-25 | 2020-06-12 | 西北大学 | 一种测井识别气层方法 |
Non-Patent Citations (3)
| Title |
|---|
| 张齐等: "三孔隙度重叠法和三孔隙度差值及比值法在保山盆地永铸街气田气层识别中的应用", 《石油天然气学报》, pages 90 - 93 * |
| 许风光等: "东海平湖油气田流体识别方法与应用研究", 《海洋石油》, pages 94 - 99 * |
| 阎媛子等: "鄂尔多斯盆地苏东南地区盒8段低阻气层识别", 《地球物理学进展》, pages 1326 - 1331 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8498177B2 (en) | Determining a position of a geological layer relative to a wavelet response in seismic data | |
| AU2009217610B2 (en) | Monitoring of reservoir fluid moving along flow pathways in a producing oil field using passive seismic emissions | |
| WO2011144623A1 (en) | Passive monitoring method for seismic events | |
| CN111679343B (zh) | 地震电磁复合数据采集系统和地下储层油气储量预测方法 | |
| CN105277982A (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
| CN104005760A (zh) | 基于方位各向异性弹性阻抗的裂缝检测方法 | |
| Williams et al. | In-situ stress orientations in the UK Southern North Sea: Regional trends, deviations and detachment of the post-Zechstein stress field | |
| CN104153768A (zh) | 一种评价花岗岩储层储集性能的方法 | |
| US11767756B2 (en) | Method and system for the determination of hydrocarbon accumulations | |
| US20170307771A1 (en) | Monitoring microseismic events | |
| Behm | Feasibility of borehole ambient noise interferometry for permanent reservoir monitoring | |
| US20150109887A1 (en) | Sonic adaptor for converting sonic or ultrasonic waveform data for use with a seismic-based computer program | |
| US11474272B2 (en) | Methods and systems for identifying and plugging subterranean conduits | |
| CN113719280A (zh) | 一种低阻气层测井识别方法 | |
| Li et al. | Hydraulic fracturing-induced microseismicity controlled by rock brittleness and natural fractures in Tongren, Guizhou, China | |
| CN106990433B (zh) | 一种隆起区微小侵蚀沟槽的识别方法 | |
| CN212364624U (zh) | 地震电磁复合数据采集系统 | |
| CN104516013A (zh) | 一种利用垂直地震数据进行时频分析与油气检测的方法 | |
| Opara et al. | Low velocity layer characterization in the Niger Delta: Implications for seismic reflection data quality | |
| Yankun et al. | Optically Quantifying Spatiotemporal Responses of Water Injection-Induced Strain via Downhole Distributed Fiber Optics Sensing | |
| CN113504569B (zh) | 阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质 | |
| Bo et al. | Study on Artificial Fracture Monitoring Technology of Coalbed Methane Reservoir Based on Microseismic Events Attribute | |
| Hardin | Fracture characterization from attenuation and generation of tube waves | |
| Adekunle et al. | Use of seismic attributes and well logs data in quantitative reservoir characterization of K-field, Onshore Niger delta area, Nigeria | |
| Westerhausl et al. | No clear evidence for temporal variations of tidal tilt prior to the 1999 Izmit and Düzce earthquakes in NW-Anatolia |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20211130 |
|
| WD01 | Invention patent application deemed withdrawn after publication |