CN108490491A - 一种基于波形指示反演的滩体预测方法 - Google Patents

一种基于波形指示反演的滩体预测方法 Download PDF

Info

Publication number
CN108490491A
CN108490491A CN201810184219.XA CN201810184219A CN108490491A CN 108490491 A CN108490491 A CN 108490491A CN 201810184219 A CN201810184219 A CN 201810184219A CN 108490491 A CN108490491 A CN 108490491A
Authority
CN
China
Prior art keywords
beach
data
waveform
inverting
beach body
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
Application number
CN201810184219.XA
Other languages
English (en)
Inventor
谭荣彪
梁虹
黎枫佶
游来胜
王潇
李晓娟
陈华
董伟
吴勇
龙思萍
段春晖
苟思
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BGP Inc
Original Assignee
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BGP Inc filed Critical BGP Inc
Priority to CN201810184219.XA priority Critical patent/CN108490491A/zh
Publication of CN108490491A publication Critical patent/CN108490491A/zh
Pending legal-status Critical Current

Links

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
    • G01V1/50—Analysing data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6169—Data from specific type of measurement using well-logging

Abstract

本发明提供了一种基于波形指示反演的滩体预测方法。所述预测方法包括：获得地震偏移数据、地震层位数据、钻井数据和测井数据；通过井‑震标定方法获得准确的时深关系，将深度域测井曲线转换为时间域测井曲线；分析得到敏感时间域测井曲线，通过所述敏感时间域测井曲线获得初始框架模型；分析得到适合滩体预测的反演参数，通过所述地震偏移数据和初始框架模型获得反演成果数据体；通过所述敏感时间域测井曲线获得滩体对应的响应值域，根据所述反演成果数据体和滩体对应的响应值域，提取出所述反演成果数据体中代表滩体的数据，获得滩体预测成果。根据本发明方法能够更加准确地预测滩体，能够对于储层和油气预测起到重要的指导意义。

Description

一种基于波形指示反演的滩体预测方法
技术领域
本发明属于地震资料解释技术领域，更具体地讲，涉及一种基于波形指示反演的滩体预测方法。
背景技术
滩体是有利的油气储集体，通过地震数据可以对滩体的规模、空间展布特征进行描述，其准确划分对于油气勘探意义重大。
传统的滩体预测大多采用定性的人工识别方法，随着科学技术水平的发展和地震资料采集技术的不断提高，地震数据中包含的信息更加丰富，地震属性分析方法被不断地应用于滩体预测中。滩体在地震剖面上多表现为地层增厚、滩体底波峰反射等特征，现阶段对滩体的识别主要为定性识别，主要通过对地震剖面波形特征、地震波峰属性等地震信息进行分析，达到滩体定性识别的效果。但是，目前的预测方法只能通过波形特征和波峰振幅属性等地震属性切片，定性地划分出滩体横向上的边界，对于多个滩体纵向叠置的情况，无法有效地识别滩体边界，因此解释结果存在精度不高、多解性强等问题。
发明内容
针对现有技术中存在的问题，本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的在于解决传统的定性划分滩体横向边界时，对于多个滩体纵向叠置的情况，无法有效地识别滩体边界，导致解释结果精度不高、多解性强的问题。
为了实现上述目的，本发明提供了一种基于波形指示反演的滩体预测方法，所述预测方法包括：获取地震偏移数据、地震层位数据、钻井数据和测井数据；获得准确的时深关系，将深度域测井曲线转换为时间域测井曲线；获取敏感时间域测井曲线，通过所述敏感时间域测井曲线获得初始框架模型；确认适合滩体预测的反演参数，通过所述地震偏移数据和初始框架模型获得反演成果数据体；通过所述敏感时间域测井曲线获得滩体对应的响应值域，根据所述反演成果数据体和滩体对应的响应值域，提取出所述反演成果数据体中代表滩体的数据，获得滩体预测成果，其中，所述获得初始框架模型包括：对所述敏感时间域测井曲线进行标准化处理；以目的层段顶、底地震层位为界面，根据地震反射特征和/或沉积规律确定地层接触关系。
在本发明的一个示例性实施例中，所述时深关系可以通过井震标定获得。
在本发明的一个示例性实施例中，所述获得初始框架模型步骤中还可以包括在标准化处理前对敏感时间域测井曲线中的异常值进行剔除。
在本发明的一个示例性实施例中，所述标准化处理可以采用直方图分析法。
在本发明的一个示例性实施例中，所述敏感时间域测井曲线的确定可以通过直方图分析法和/或交汇图分析法。
在本发明的一个示例性实施例中，所述反演参数可以包括波阻抗曲线、有效样本数、平滑系数和截止频率。
在本发明的一个示例性实施例中，所述有效样本数可以为3～5。
在本发明的一个示例性实施例中，所述平滑系数可以为1～3。
在本发明的一个示例性实施例中，所述截止频率中低频部分的高通频率可以为5～8Hz，高截频率可以为8～10Hz，高频部分的低截频率可以为40～60Hz，低通频率可以为60～80Hz，高通频率可以为150～300Hz，高截频率可以为250～400Hz。
与现有技术相比，本发明的有益效果包括：基于波形指示反演的滩体预测方法能够识别多个滩体纵向叠置情况下的滩体边界；能够提高解释结果的精度；能够降低多解性；能够更加准确的预测滩体；能够对于储层和油气预测起到重要的指导意义。
附图说明
图1示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的一个流程图。
图2示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的另一个流程图。
图3示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的预测剖面图。
图4示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的预测平面图。
具体实施方式
在下文中，将结合示例性实施例和附图来详细说明本发明的一种基于波形指示反演的滩体预测方法。
图1示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的流程图。如图1所示，在一个示例性实施例中，本发明的基于波形指示反演的滩体预测方法可通过如下步骤实现：
步骤(1)，获取地震偏移数据、地震层位数据、测井数据和钻井数据。例如，获取处理好的四川盆地某三维工区地震偏移数据，作为波形指示反演的基础数据；获取目的层段(例如，龙王庙组)的顶、底地震层位，作为波形指示反演的约束时窗；获取该三维工区内的钻井数据和测井数据，作为井-震标定的基础数据。该三维工区的地震偏移数据可以通过野外数据采集和叠前时间偏移处理得到。
步骤(2)，通过井-震标定方法获得准确的时深关系，将深度域测井曲线转换为时间域测井曲线。步骤(1)中获取的测井数据为深度域，需要将深度域测井数据转换为时间域测井数据用于波形指示反演。
步骤(3)，分析得到敏感时间域测井曲线，通过敏感时间域测井曲线获得初始框架模型。对该三维工区中的时间域测井曲线进行分析，确定出存在异常值的时间域测井曲线，剔除时间域测井曲线中的异常值，然后对正常的时间域测井曲线进行标准化处理，例如，可以采用直方图分析法进行标准化处理；例如，可以根据时间域测井曲线中的声波速度、密度、伽马、孔隙度、波阻抗等，选取目的层段样点，根据实际钻井情况，将样点分为滩体和非滩体两类，可以通过直方图分析方法和/或交汇图分析方法，对目的层段时间域测井曲线进行分析，确定出能够区分出滩体和非滩体的时间域测井曲线，得到敏感时间域测井曲线，例如，针对该三维工区的敏感时间域测井曲线可以为波阻抗曲线；可以目的层段顶、底地震层位为界面，可以根据地震反射特征和/或沉积规律确定地层接触关系，可以井点敏感时间域测井曲线为基础，进行全区线性内插，计算得到初始框架模型。
步骤(4)，分析得到适合滩体预测的反演参数，通过地震偏移数据和初始框架模型获得反演成果数据体。通过试验得出适合该三维工区滩体预测的反演参数，例如，反演参数可以包括输入曲线、有效样本数、平滑系数和截止频率，其中，输入曲线可以为波阻抗曲线，有效样本数可以为3～5，优选地，有效样本数为5；平滑系数可以为1～3，优选地，平滑系数为1；截止频率中低频部分的高通频率可以为5～8Hz，高截频率可以为8～10Hz，高频部分的低截频率可以为40～60Hz，低通频率可以为60～80Hz，高通频率可以为150～300Hz，高截频率可以为250～400Hz，优选地，最佳截止频率中低频部分的高通频率可以为5Hz，高截频率可以为10Hz，高频部分的低截频率可以为60Hz，低通频率可以为80Hz，高通频率可以为200Hz，高截频率可以为300Hz。通过波形指示反演获得反演成果数据体。
步骤(5)，通过敏感时间域测井曲线获得滩体对应的响应值域，根据步骤(4)获得的反演成果数据体和滩体对应的响应值域，提取出反演成果数据体中代表滩体的数据，获得滩体预测成果，用于储层含油气性的预测，如图3、图4所示。如图3所示为龙王庙组滩体预测剖面图，图中的A、B、C区域分别表示预测出的三个滩体，虚线表示滩体的边界。如图4所示为采用本方法对四川盆地某三维工区龙王庙组开展的滩体预测效果图，图中磨溪以及宝龙加数字代表井名，钻井钻遇滩体的井均位于滩体发育区，滩体预测符合率100％，滩体预测成果对于储层与含油气性解释有很好的定性指导作用。
如图2所示，在本发明的另一个示例性实施例中，基于波形指示反演的滩体预测方法可由以下步骤实现：
步骤(1)，获取地震偏移数据、地震层位数据、钻井数据和测井数据。首先，以本领域公知的方法进行野外勘探采集，得到原始采集的数据，经过本领域公知的叠前时间偏移处理方法进行地震资料处理，得到地震偏移数据，作为波形指示反演的基础数据；然后，获取目的层段的顶、底地震层位，作为波形指示反演的约束时窗；最后，获得工区内钻、测井数据，作为井-震标定的基础数据。例如，应用本方法可对四川盆地某三维工区龙王庙组开展滩体预测。
步骤(2)，井-震标定。通过本领域公知的井-震标定方法，对目的层段滩体进行标定，从而确定滩体的地震响应特征，并为下一步波形指示反演提供准确的时深关系。
步骤(3)，波形指示反演。地震有效频带一般为10～60HZ，无法直接获取高频成分，但地震波形的横向变化反映了沉积环境的变化，而相似的沉积环境具有可类比的波形特征，因此可以充分利用地震波形的横向变化开展高频成分估计。波形指示反演和传统统计学反演最大的区别就在于统计样本的筛选，传统方法是基于空间域变差函数的，只能粗略表达空间变异程度，无法体现相变特征，而且对样本分布均匀程度要求比较苛刻，限制该方法的应用领域。而波形指示反演利用沉积学基本原理，充分利用地震波形的横向变化来反应储层空间的相变特征，进而分析储层垂向岩性组合高频结构特征，更好地体现了相控的思想，是一种真正的井震结合高频模拟方法，使反演结果从完全随机到逐步确定，同时对井位分布的均匀性没有严格要求。反演步骤如下：
①数据预处理
a、对工区内井数据进行分析，找到存在异常值的测井曲线，剔除曲线中的异常值，并存储为新曲线名备用。
b、可以采用本领域公知的直方图分析法对测井曲线进行标准化处理。
②敏感参数分析
可以采用测井曲线中的声波速度、密度、伽马、孔隙度、波阻抗等曲线，选取目的层段样点，根据钻井实际情况，将样点分为滩体和非滩体两类，可以通过直方图、交汇图分析方法，对目的层段测井曲线进行分析，找出能区分出滩体和非滩体的测井曲线，例如，在本工区中波阻抗曲线能够区分出滩体和非滩体。
③框架模型建立
以目的层段顶、底地震层位为界面，可以依据地震反射特征和沉积规律设置地层接触关系，可以井点波阻抗曲线值为基础，进行全区线性内插，计算出初始框架模型。
④地震反演
通过试验，可以确定适合本工区滩体预测的反演参数，完成波形指示反演，获得反演成果数据体。反演参数可以包括波阻抗曲线、有效样本数、平滑系数和截止频率，其中，输入曲线可以为波阻抗曲线，有效样本数可以为3～5，优选地，有效样本数为5；平滑系数可以为1～3，优选地，平滑系数为1；截止频率中低频部分的高通频率可以为5～8Hz，高截频率可以为8～10Hz，高频部分的低截频率可以为40～60Hz，低通频率可以为60～80Hz，高通频率可以为150～300Hz，高截频率可以为250～400Hz，优选地，最佳截止频率中低频部分的高通频率可以为5Hz，高截频率可以为10Hz，高频部分的低截频率可以为60Hz，低通频率可以为80Hz，高通频率可以为200Hz，高截频率可以为300Hz。通过波形指示反演获得反演成果数据体。
步骤(4)，滩体预测。根据测井曲线敏感参数分析结果，以滩体对应响应值域为提取范围，提取出反演成果数据体中代表滩体的数据，在反演剖面上将滩体顶、底界对比出来，如图3所示为龙王庙组滩体预测剖面图，图中A、B、C区域分别表示预测出的三个滩体，虚线表示滩体的边界。
如图4所示为采用本方法对四川盆地某三维工区龙王庙组开展的滩体预测效果图，图中磨溪以及宝龙加数字代表井名，钻井钻遇滩体的井均位于滩体发育区，滩体预测符合率100％，滩体预测成果对于储层与含油气性解释有很好的定性指导作用。
综上所述，本发明提出了一种基于波形指示反演的滩体预测方法，以地震偏移数据为输入，采用波形指示反演技术，获得目的层段波阻抗反演数据体。根据实际钻井资料和井-震标定，获得钻井已钻遇滩体的波阻抗值域范围，从根据波形指示反演得到的波阻抗反演数据体中提取出滩体值域范围值的数据，可以有效地实现对滩体的预测，实践表明该方法对于该工区的储层和油气预测具有重要的指导意义。
尽管上面已经结合附图和示例性实施例描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。

Claims (9)

1.一种基于波形指示反演的滩体预测方法，其特征在于，所述预测方法包括以下步骤：
获取地震偏移数据、地震层位数据、钻井数据和测井数据；
获得准确的时深关系，将深度域测井曲线转换为时间域测井曲线；
获取敏感时间域测井曲线，通过所述敏感时间域测井曲线获得初始框架模型；
确认适合滩体预测的反演参数，通过所述地震偏移数据和初始框架模型获得反演成果数据体；
通过所述敏感时间域测井曲线获得滩体对应的响应值域，根据所述反演成果数据体和滩体对应的响应值域，提取出所述反演成果数据体中代表滩体的数据，获得滩体预测成果，其中，
所述获得初始框架模型包括：对所述敏感时间域测井曲线进行标准化处理；以目的层段顶、底地震层位为界面，根据地震反射特征和/或沉积规律确定地层接触关系。
2.根据权利要求1所述的基于波形指示反演的滩体预测方法，其特征在于，所述时深关系通过井震标定获得。
3.根据权利要求1所述的基于波形指示反演的滩体预测方法，其特征在于，所述获得初始框架模型步骤中还包括在标准化处理前对敏感时间域测井曲线中的异常值进行剔除。
4.根据权利要求1所述的基于波形指示反演的滩体预测方法，其特征在于，所述标准化处理采用直方图分析法。
5.根据权利要求1所述的基于波形指示反演的滩体预测方法，其特征在于，所述敏感时间域测井曲线的确定通过直方图分析法和/或交汇图分析法。
6.根据权利要求1所述的基于波形指示反演的滩体预测方法，其特征在于，所述反演参数包括波阻抗曲线、有效样本数、平滑系数和截止频率。
7.根据权利要求6所述的基于波形指示反演的滩体预测方法，其特征在于，所述有效样本数为3～5。
8.根据权利要求6所述的基于波形指示反演的滩体预测方法，其特征在于，所述平滑系数为1～3。
9.根据权利要求6所述的基于波形指示反演的滩体预测方法，其特征在于，所述截止频率中低频部分的高通频率为5～8Hz，高截频率为8～10Hz，高频部分的低截频率为40～60Hz，低通频率为60～80Hz，高通频率为150～300Hz，高截频率为250～400Hz。
CN201810184219.XA 2018-03-06 2018-03-06 一种基于波形指示反演的滩体预测方法 Pending CN108490491A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810184219.XA CN108490491A (zh) 2018-03-06 2018-03-06 一种基于波形指示反演的滩体预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810184219.XA CN108490491A (zh) 2018-03-06 2018-03-06 一种基于波形指示反演的滩体预测方法

Publications (1)

Publication Number Publication Date
CN108490491A true CN108490491A (zh) 2018-09-04

Family

ID=63341588

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810184219.XA Pending CN108490491A (zh) 2018-03-06 2018-03-06 一种基于波形指示反演的滩体预测方法

Country Status (1)

Country Link
CN (1) CN108490491A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111077578A (zh) * 2018-10-22 2020-04-28 中国石油天然气股份有限公司 岩层分布预测方法和装置
CN111856566A (zh) * 2019-04-28 2020-10-30 中国石油天然气股份有限公司 湖相滩坝砂体中薄储层预测方法及装置
CN112180440A (zh) * 2019-07-03 2021-01-05 中国石油天然气集团有限公司 一种基于avo特征分析的叠前随机反演方法及系统
CN112444857A (zh) * 2019-09-02 2021-03-05 中国石油天然气股份有限公司 一种基于叠后地震资料的膏盐地球物理预测方法
CN113138407A (zh) * 2020-01-20 2021-07-20 中国石油天然气集团有限公司 一种深层页岩气的多尺度断裂地震预测方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150316673A1 (en) * 2012-12-05 2015-11-05 Jacky M. Wiener Systems and Methods for 3D Seismic Data Depth Conversion Utilizing Artificial Neural Networks
CN105353407A (zh) * 2015-10-28 2016-02-24 中国石油化工股份有限公司 一种叠后地震波阻抗反演方法
CN107121699A (zh) * 2017-02-28 2017-09-01 中国石油化工股份有限公司 一种地震相控制下的沉积微相识别方法
CN107329171A (zh) * 2017-06-07 2017-11-07 中国石油天然气股份有限公司 深度域储层地震反演方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150316673A1 (en) * 2012-12-05 2015-11-05 Jacky M. Wiener Systems and Methods for 3D Seismic Data Depth Conversion Utilizing Artificial Neural Networks
CN105353407A (zh) * 2015-10-28 2016-02-24 中国石油化工股份有限公司 一种叠后地震波阻抗反演方法
CN107121699A (zh) * 2017-02-28 2017-09-01 中国石油化工股份有限公司 一种地震相控制下的沉积微相识别方法
CN107329171A (zh) * 2017-06-07 2017-11-07 中国石油天然气股份有限公司 深度域储层地震反演方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周影: "《地震储层预测在油气勘探中的应用》", 《中国优秀硕士论文全文数据库基础科学辑》 *
顾雯等: "《强屏蔽下薄储层高精度预测研究》", 《石油物探》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111077578A (zh) * 2018-10-22 2020-04-28 中国石油天然气股份有限公司 岩层分布预测方法和装置
CN111856566A (zh) * 2019-04-28 2020-10-30 中国石油天然气股份有限公司 湖相滩坝砂体中薄储层预测方法及装置
CN112180440A (zh) * 2019-07-03 2021-01-05 中国石油天然气集团有限公司 一种基于avo特征分析的叠前随机反演方法及系统
CN112180440B (zh) * 2019-07-03 2023-05-26 中国石油天然气集团有限公司 一种基于avo特征分析的叠前随机反演方法及系统
CN112444857A (zh) * 2019-09-02 2021-03-05 中国石油天然气股份有限公司 一种基于叠后地震资料的膏盐地球物理预测方法
CN112444857B (zh) * 2019-09-02 2023-10-27 中国石油天然气股份有限公司 一种基于叠后地震资料的膏盐地球物理预测方法
CN113138407A (zh) * 2020-01-20 2021-07-20 中国石油天然气集团有限公司 一种深层页岩气的多尺度断裂地震预测方法及系统

Similar Documents

Publication Publication Date Title
CN113759425B (zh) 井震联合评价深层古岩溶储层充填特征的方法与系统
CN108490491A (zh) 一种基于波形指示反演的滩体预测方法
CN109061765B (zh) 非均质薄砂岩互层油藏的圈闭评价方法
JP3390445B2 (ja) 地質学的特性を推定するために統計学的較正技術を使用する地震トレース解析方法
CN110318745B (zh) 一种沉积微相约束下粒径岩性测井评价方法
CN105093306B (zh) 一种地球物理勘探中储层自动解释与厚度求取方法
CN105986813A (zh) 准致密储层快速评价方法及准致密储层多指标评价方法
CN114114459B (zh) 一种相控约束下的深层-超深层碳酸盐岩薄储层预测方法
CN105240006A (zh) 一种适用于火山岩储层的油水层识别方法
CN108957540B (zh) 一种高效提取复杂储层中衰减品质因子的方法
CN112963145B (zh) 一种碳酸盐岩储层气井产能的预测方法
CN112505754B (zh) 基于高精度层序格架模型的井震协同划分沉积微相的方法
CN104459790A (zh) 含油气性盆地有效储层的分析方法和装置
CN116168224A (zh) 基于成像砾石含量的机器学习岩相自动识别方法
CN106033127A (zh) 基于横波速度变化率的地应力方位地震预测方法
CN115390155A (zh) 针对非均质储层的测井解释方法、装置、电子设备及介质
CN112528106A (zh) 火山岩岩性识别方法
CN109707378A (zh) 一种基于泥浆侵入特性及纵向对比的低阻油层识别方法
CN112946782B (zh) 一种致密油气储渗体地震精细刻画方法
US11719851B2 (en) Method and system for predicting formation top depths
CN106353813A (zh) 基于阵列声波测井的流体性质识别方法
CN113514884A (zh) 一种致密砂岩储层预测方法
CN109581489B (zh) 嵌套式地震相的提取方法及系统
CN112731556B (zh) 一种裂缝发育区的预测方法及用于预测裂缝发育区的计算机存储介质
CN115932968B (zh) 一种基于地震振幅比值属性的碳酸盐岩薄储层预测方法

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180904