CN107918150A - 基于海上油田大井距下的单河道识别方法 - Google Patents

基于海上油田大井距下的单河道识别方法 Download PDF

Info

Publication number
CN107918150A
CN107918150A CN201610884783.3A CN201610884783A CN107918150A CN 107918150 A CN107918150 A CN 107918150A CN 201610884783 A CN201610884783 A CN 201610884783A CN 107918150 A CN107918150 A CN 107918150A
Authority
CN
China
Prior art keywords
river
single river
recognition methods
wide spacing
offshore oilfield
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
CN201610884783.3A
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.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
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 China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201610884783.3A priority Critical patent/CN107918150A/zh
Publication of CN107918150A publication Critical patent/CN107918150A/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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • 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/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种基于海上油田大井距下的单河道识别方法，该方法采用“井震结合、模式指导与动态约束”相结合，包括：步骤1，分频段对地震属性进行优选，重构地震数据体，从重构的数据体提取与砂体厚度相关性最好的属性，预测复合河道分布；步骤2，进行分频带反演，得到反演数据体；步骤3，根据经验公式确定单河道带的规模，初步划分单河道；以及步骤4，结合运用反演剖面高程差异、正演剖面以及动态资料分析进一步划分单河道。该基于海上油田大井距下的单河道识别方法技术思路清楚、应用简单，突破了以往方法的局限性，与以往相比具有更好的操作性，具有创新性、实用性，利于推广。

Description

基于海上油田大井距下的单河道识别方法
技术领域
本发明涉及石油开发中的地质研究技术领域，特别是涉及到基于海上油田大井距下的单河道识别方法。
背景技术
曲流河作为陆相湖盆沉积体系的主要构成单元，是我国已发现油田的重要储层类型。据统计，河流相储层石油储量占我国已开发油田动用储量的40％以上。自上个世纪60-70年代以来，我国大部分油田经历了十几年甚至几十年的注水开发，目前已处于高含水采油阶段，开发矛盾比较突出。但是，在钻新井时仍然能发现有高产井，这是储层内部非均质性的突出表现，即储层内部构型单元界面形成的渗流屏障及构型单元渗流差异导致的剩余油富集。在很多密井网小井距的陆上油田已经开展了构型侧积体级别的研究，并取得了一定的效果。但是海上油田一般井距大，资料少，很多方法不能适应，因此发明一种适合海上油田大井距的单河道识别方法对砂体连通模式的认识至关重要。
发明内容
本发明的目的是提供一种曲流河单河道划分的方法，实现了基于海上油田大井距下的单河道划分。
本发明的目的可通过如下技术措施来实现：基于海上油田大井距下的单河道识别方法，该基于海上油田大井距下的单河道识别方法包括：步骤1，分频段对地震属性进行优选，重构地震数据体，从重构的数据体提取与砂体厚度相关性最好的属性，预测复合河道分布；步骤2，进行分频带反演，得到反演数据体；步骤3，根据公式确定单河道带的规模，初步划分单河道；以及步骤4，结合运用反演剖面高程差异、正演剖面以及动态资料分析进一步划分单河道。
本发明的目的还可通过如下技术措施来实现：
在步骤1中，通过地震资料分析，优选与河道砂体相关性最好的属性，然后按照不同的中心频率进行小波分解，得到具有一定主频和频宽的分频地震数据体。
在步骤1中，在对不同频带地震数据体分析的基础上，选取与河道砂体相关性较好的频带地震数据体进行融合重构，按不同方式组合可以得到多个重构地震数据体。
在步骤1中，分别从原始地震数据和重构地震数据提取振幅、频率、相位类多种地震属性，并分析各地震属性与单井解释的砂体厚度之间的相关关系，找到与砂体厚度相关性最好的属性，根据属性规律刻画复合河道带的砂体边界。
在步骤2中，对研究区各测井曲线进行分析，挑选出能较好反演岩性的测井曲线，并进行归一化处理；选取中心频段对原始地震数据进行小波分频，进而得到分频反演数据体和各连井反演剖面。
在步骤3中，根据经验公式计算出满岸宽度和单河道宽度，
单河道规模通过单井测算完整正旋回厚度，恢复压实作用后利用Leeder公式求取满岸宽度：log(w)＝1.54×log(h)+0.83
其中，w为河流满岸宽度，m；h为河流满岸深度，m；
求取满岸宽度后可再通过Lorenz公式计算单河道宽度：
其中，Wm为单河道宽度，m；W2为活动河道宽度，m。
在步骤4中，在单河道宽度的指导下，根据分频反演剖面协助划分相邻的两口井应该分别归属哪一条单河道内。
在步骤4中，在单河道宽度的指导下以及概念模型正演响应特征分析的基础上，根据正演剖面特征分析，协助划分单河道带。
在步骤4中，在单河道边界的井，分析注水井与生产井的产量变化，得到邻井的连通状况，进而根据连通状况判断单河道的划分位置。
本发明中的基于海上油田大井距下的单河道识别方法采用“井震结合、模式指导与动态约束”相结合，通过多个步骤，多个方法的结合，满足了大井距条件下单河道砂体识别的需要。该方法技术思路清楚、应用简单，突破了单独一种方法识别单河道的局限，克服了大井距下完钻井资料少划分单河道的困难，与以往相比具有更好的效果，具有创新性、实用性，利于推广。为海上大井距单河道识别提供了切实可行的方法。该方法在胜利埕岛油田埕北11井区获得应用，成功划分馆上段4砂组、5砂组的各小层单河道分布，为砂体连通关系研究提供了地质基础。
附图说明
图1为本发明的基于海上油田大井距下的单河道识别方法的一具体实施例的流程图；
图2为本发明的一具体实施例中规模估算得到的单河道分布平面图；
图3为本发明的一具体实施例中对子井①的对比剖面与反演剖面图；
图4为本发明的一具体实施例中对子井②的对比剖面图；
图5为本发明的一具体实施例中对子井②的正演剖面图；
图6为本发明的一具体实施例中对子井④的对比剖面与反演剖面图；
图7为本发明的一具体实施例中连井③的对比剖面图；
图8为本发明的一具体实施例中连井③的开发曲线图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的基于海上油田大井距下的单河道识别方法的流程图。该基于海上油田大井距下单河道识别方法采用“井震结合、模式指导与动态约束”相结合，包括：
步骤101，分频段对地震属性进行优选，分析各频段数据体与砂体发育的关系，优选出较好反映砂体发育的频段，选取特定主频的分频数据体进行融合，从而重构出适合预测特定砂体厚度的“重构地震数据体”。从重构的数据体提取与砂体厚度相关性最好的属性，预测复合河道分布。
通过地震资料分析，优选与河道砂体相关性最好的属性，然后按照不同的中心频率进行小波分解，得到具有一定主频和频宽的分频地震数据体。在对不同频带地震数据体分析的基础上，选取与河道砂体相关性较好的频带地震数据体进行融合重构，按不同方式组合可以得到多个重构地震数据体。分别从原始地震数据和重构地震数据提取振幅、频率、相位类多种地震属性，并分析各地震属性与单井解释的砂体厚度之间的相关关系，找到与砂体厚度相关性最好的属性，根据属性规律刻画复合河道带的砂体边界。
步骤102，进行分频带反演，得到反演数据体。对研究区各测井曲线进行分析，挑选出能较好反演岩性的测井曲线，并进行归一化处理。选取中心频段对原始地震数据进行小波分频，进而得到分频反演数据体和各连井反演剖面。
步骤103，根据经验公式确定单河道带的规模，初步划分单河道。根据前人的经验公式计算出满岸宽度和单河道宽度。
单河道规模通过单井测算完整正旋回厚度，恢复压实作用后利用Leeder公式求取满岸宽度：
其中，w为河流满岸宽度，m；h为河流满岸深度，m；
求取满岸宽度后可再通过Lorenz公式计算单河道宽度：
其中，Wm为单河道宽度，m；W2为活动河道宽度，m。
步骤104，结合运用反演剖面高程差异、正演剖面以及动态资料分析进一步划分单河道。在单河道宽度的指导下，根据分频反演剖面协助划分相邻的两口井应该分别归属哪一条单河道内。
在单河道宽度的指导下以及概念模型正演响应特征分析的基础上，根据正演剖面特征分析，协助划分单河道带。
在单河道边界的井，可以分析注水井与生产井的产量变化，得到邻井的连通状况，进而根据连通状况判断单河道的划分位置。
在应用本发明的一具体实施例中，包括了以下步骤：
在步骤1，通过地震资料分析，埕北11井区为曲流河沉积，振幅类属性与河道砂体的分布具有较好的正相关性，对振幅类属性按照10～70Hz(间距为10Hz)7个中心频率进行小波分解，得到具有一定主频和频宽的分频地震数据体。
对不同频带地震数据体分析，中心频率为10Hz时，地震属性高值区凌乱分布，并无明显规律；中心频率提高到20Hz时，属性已经呈现出一条较为明显的蛇曲状高值带，且在属性高值带以外零星分布的高值区很少；到70Hz时，属性规律又变差。因此，分别选取20、30、40、50、60Hz分频数据体融合得到重构“地震数据体A”，20、30、40、50Hz分频数据体融合得到重构“地震数据体B”，25、35、45、55Hz分频数据体融合得到重构“地震数据体C”，选取30、40、50、60Hz分频数据体融合得到重构“地震数据体D”。
分别从原始地震数据和重构地震数据提取振幅、频率、相位类多种地震属性，并分析各地震属性与单井解释的砂体厚度之间的相关关系，找到与砂体厚度相关性最好的属性，根据属性规律刻画复合河道带的砂体边界。地震属性与砂体厚度相关性分析表如表1所示。
表1 Ng541单层地震属性与砂体厚度相关性分析表
在步骤2，对研究区各测井曲线进行分析，挑选出能较好反演岩性的测井曲线，并进行归一化处理。选取中心频段对原始地震数据进行小波分频，进而得到分频反演数据体和各连井反演剖面。
在步骤3，根据前人的经验公式计算出满岸宽度和单河道宽度。
单河道规模通过单井测算完整正旋回厚度，恢复压实作用后利用Leeder公式求取满岸宽度：log(w)＝1.54×log(h)+0.83
其中，w为河流满岸宽度，m；h为河流满岸深度，m；
求取满岸宽度后可再通过Lorenz公式计算单河道宽度：
其中，Wm为单河道宽度，m；W2为活动河道宽度，m。
图2为Ng54 1单层根据公式估算得到的3条单河道，单河道1估算宽度约610m，单河道2估算宽度约800m，单河道3估算宽度约640m。
在步骤4，在单河道宽度的指导下，根据分频反演剖面协助划分相邻的两口井应该分别归属哪一条单河道内。图3为具体实例Ng54 1单层对子井11NA-3与11F-1的对比剖面和反演剖面。根据这个剖面可以确定单河道在此处的边界位置。
图4为具体实例Ng54 1单层对子井11NB-6与CB11E-6的对比剖面，图5为具体实例Ng54 1单层对子井11NB-6与CB11E-6的正演剖面。根据这两个剖面可以确定单河道在此处的边界位置。
图6为具体实例Ng54 1单层对子井11NB-7与11H-4的对比剖面和反演剖面。根据这个剖面可以确定单河道在此处的边界位置。
图7为具体实例Ng54 1连井11NC-1、11NB-3与11NB-1的对比剖面，图8为具体实例Ng54 1连井11NC-1、11NB-3与11NB-1的开发曲线。在Ng54 1层11NB-3给11NC-1井注水，开发曲线显示注水受效明显，说明此处两井连通性好，应该为同一单河道内，而11NB-3给11NB-1井注水，开发曲线显示注水不受效，说明此处两井连通性差，应该为不同的单河道内，根据单井动态分析可以判断此处单河道边界应在11NB-3与11NC-1之间。
以此类推可以进一步确定单河道2与单河道3的边界位置，这样就可以把在Ng54 1层的复合河道划分为3条单河道了。

Claims (9)

1.基于海上油田大井距下的单河道识别方法，其特征在于，该基于海上油田大井距下的单河道识别方法包括：
步骤1，分频段对地震属性进行优选，重构地震数据体，从重构的数据体提取与砂体厚度相关性最好的属性，预测复合河道分布；
步骤2，进行分频带反演，得到反演数据体；
步骤3，根据公式确定单河道带的规模，初步划分单河道；以及
步骤4，结合运用反演剖面高程差异、正演剖面以及动态资料分析进一步划分单河道。
2.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤1中，通过地震资料分析，优选与河道砂体相关性最好的属性，然后按照不同的中心频率进行小波分解，得到具有一定主频和频宽的分频地震数据体。
3.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤1中，在对不同频带地震数据体分析的基础上，选取与河道砂体相关性较好的频带地震数据体进行融合重构，按不同方式组合可以得到多个重构地震数据体。
4.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤1中，分别从原始地震数据和重构地震数据提取振幅、频率、相位类多种地震属性，并分析各地震属性与单井解释的砂体厚度之间的相关关系，找到与砂体厚度相关性最好的属性，根据属性规律刻画复合河道带的砂体边界。
5.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤2中，对研究区各测井曲线进行分析，挑选出能较好反演岩性的测井曲线，并进行归一化处理；选取中心频段对原始地震数据进行小波分频，进而得到分频反演数据体和各连井反演剖面。
6.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤3中，根据经验公式计算出满岸宽度和单河道宽度，
单河道规模通过单井测算完整正旋回厚度，恢复压实作用后利用Leeder公式求取满岸宽度：log(w＝1.54×log h+0.83
其中，w为河流满岸宽度，m；h为河流满岸深度，m；
求取满岸宽度后可再通过Lorenz公式计算单河道宽度：
其中，Wm为单河道宽度，m；W2为活动河道宽度，m。
7.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤4中，在单河道宽度的指导下，根据分频反演剖面协助划分相邻的两口井应该分别归属哪一条单河道内。
8.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤4中，在单河道宽度的指导下以及概念模型正演响应特征分析的基础上，根据正演剖面特征分析，协助划分单河道带。
9.根据权利要求1所述的基于海上油田大井距下的单河道识别方法，其特征在于，在步骤4中，在单河道边界的井，分析注水井与生产井的产量变化，得到邻井的连通状况，进而根据连通状况判断单河道的划分位置。
CN201610884783.3A 2016-10-10 2016-10-10 基于海上油田大井距下的单河道识别方法 Pending CN107918150A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610884783.3A CN107918150A (zh) 2016-10-10 2016-10-10 基于海上油田大井距下的单河道识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610884783.3A CN107918150A (zh) 2016-10-10 2016-10-10 基于海上油田大井距下的单河道识别方法

Publications (1)

Publication Number Publication Date
CN107918150A true CN107918150A (zh) 2018-04-17

Family

ID=61892454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610884783.3A Pending CN107918150A (zh) 2016-10-10 2016-10-10 基于海上油田大井距下的单河道识别方法

Country Status (1)

Country Link
CN (1) CN107918150A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109083640A (zh) * 2018-09-26 2018-12-25 中国石油化工股份有限公司 海上河流相油藏边部储量有效动用方法
CN110412649A (zh) * 2019-06-17 2019-11-05 中国石油化工股份有限公司 一种单期水下分流河道的识别方法
CN111997604A (zh) * 2020-10-16 2020-11-27 陕西延长石油（集团）有限责任公司 一种单期河道边界的识别方法
CN112651547A (zh) * 2020-12-11 2021-04-13 中国石油天然气股份有限公司 一种曲流河砂岩储层规模预测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502969A (zh) * 2014-12-30 2015-04-08 中国石油化工股份有限公司 河道砂岩性油藏的识别方法
CN104698496A (zh) * 2013-12-05 2015-06-10 中国石油化工股份有限公司 一种小砂体边界判识方法及小砂体空间定量描述方法
CN105116447A (zh) * 2015-08-14 2015-12-02 中国海洋石油总公司 一种基于曲率异常条带的地质河道方向判别方法
CN105182444A (zh) * 2015-09-09 2015-12-23 中国海洋石油总公司 一种高分辨率层序地层格架约束的地质统计学反演方法
CN105372703A (zh) * 2015-08-11 2016-03-02 中国石油化工股份有限公司 一种河道砂体精细识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104698496A (zh) * 2013-12-05 2015-06-10 中国石油化工股份有限公司 一种小砂体边界判识方法及小砂体空间定量描述方法
CN104502969A (zh) * 2014-12-30 2015-04-08 中国石油化工股份有限公司 河道砂岩性油藏的识别方法
CN105372703A (zh) * 2015-08-11 2016-03-02 中国石油化工股份有限公司 一种河道砂体精细识别方法
CN105116447A (zh) * 2015-08-14 2015-12-02 中国海洋石油总公司 一种基于曲率异常条带的地质河道方向判别方法
CN105182444A (zh) * 2015-09-09 2015-12-23 中国海洋石油总公司 一种高分辨率层序地层格架约束的地质统计学反演方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐洪斌 等: "《地层、岩性油气藏地震勘探方法与技术》", 31 July 2012, 石油工业出版社 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109083640A (zh) * 2018-09-26 2018-12-25 中国石油化工股份有限公司 海上河流相油藏边部储量有效动用方法
CN110412649A (zh) * 2019-06-17 2019-11-05 中国石油化工股份有限公司 一种单期水下分流河道的识别方法
CN110412649B (zh) * 2019-06-17 2021-06-18 中国石油化工股份有限公司 一种单期水下分流河道的识别方法
CN111997604A (zh) * 2020-10-16 2020-11-27 陕西延长石油（集团）有限责任公司 一种单期河道边界的识别方法
CN111997604B (zh) * 2020-10-16 2023-12-08 陕西延长石油（集团）有限责任公司 一种单期河道边界的识别方法
CN112651547A (zh) * 2020-12-11 2021-04-13 中国石油天然气股份有限公司 一种曲流河砂岩储层规模预测方法
CN112651547B (zh) * 2020-12-11 2024-05-28 中国石油天然气股份有限公司 一种曲流河砂岩储层规模预测方法

Similar Documents

Publication Publication Date Title
Sullivan et al. Architectural analysis of deep-water outcrops: Implications for exploration and development of the Diana sub-basin, western Gulf of Mexico
Fonnesu et al. Short length-scale variability of hybrid event beds and its applied significance
Sylvester et al. Seismic stratigraphy of a shelf-edge delta and linked submarine channels in the northeastern Gulf of Mexico
CN103993862B (zh) 复杂断块稀油油藏分层开发方法
CN103510947B (zh) 滩坝砂岩微相识别模式的建立方法
CN107918150A (zh) 基于海上油田大井距下的单河道识别方法
Francheteau et al. Oceanic fracture zones do not provide deep sections in the crust
CN104453834A (zh) 一种井组注采关系优化调整方法
CN101936165A (zh) 一种岩溶型碳酸盐岩储层布井方法
CN106437674A (zh) 仿水平井注水开发井网适配方法
CN104698496A (zh) 一种小砂体边界判识方法及小砂体空间定量描述方法
Zhuo et al. Seismic geomorphology, architecture and genesis of Miocene shelf sand ridges in the Pearl River Mouth Basin, northern South China Sea
CN107966731B (zh) 一种基于地震波形结构属性的河流相砂体叠置区识别方法
CN109783944A (zh) 三维流动单元模型建立方法、装置和设备
Meckel Shelf margin deltas: the key to big reserves
Taixian et al. Techniques for high-efficient development of offshore fluvial oilfields
Navarro et al. Depositional architecture and evolution of a deep-marine channel-levee complex: Isaac Formation (Windermere Supergroup), Southern Canadian Cordillera
Cummings et al. Sedimentology and stratigraphy of a thick, areally extensive fluvial–marine transition, Missisauga Formation, offshore Nova Scotia, and its correlation with shelf margin and slope strata
Zhang et al. Sedimentary characteristics and controlling factors of shelf sand ridges in the Pearl River Mouth Basin, northeast of South China Sea
Cliff et al. Exploration of the Rakhine Basin, pushing out the barriers with new 3D
Anderson et al. Late Quaternary shelf-margin delta and slope-fan complexes of the east Texas–western Louisiana margin: Variable response to eustasy and sediment supply
Brunt et al. Heterogeneity of fill within an incised channel: the Oligocene Gres du Champsaur, SE France
Pinghua et al. Comprehensive Dissection of Meandering-River Sedimentary System in Qinhuangdao 32-6 Oilfield, Bohai Bay Basin
Zhang et al. Case studies on removing strong reflections in carbonate and sand reservoir prediction using volume-based seismic waveform decomposition
Potapova et al. Integrated approach to the correlation polyfacies objects deposits

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20180417