CN106772614A - 一种重力流水道砂体中优质储层的预测方法 - Google Patents
一种重力流水道砂体中优质储层的预测方法 Download PDFInfo
- Publication number
- CN106772614A CN106772614A CN201710111926.1A CN201710111926A CN106772614A CN 106772614 A CN106772614 A CN 106772614A CN 201710111926 A CN201710111926 A CN 201710111926A CN 106772614 A CN106772614 A CN 106772614A
- Authority
- CN
- China
- Prior art keywords
- flow channel
- gravity
- sandbody
- reservoir
- porosity
- 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 26
- 239000004576 sand Substances 0.000 title claims abstract description 22
- 230000005484 gravity Effects 0.000 title claims abstract description 20
- 239000011435 rock Substances 0.000 claims abstract description 21
- 230000002452 interceptive effect Effects 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 6
- 239000004744 fabric Substances 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 230000003667 anti-reflective effect Effects 0.000 claims 1
- 239000000463 material Substances 0.000 claims 1
- 239000011148 porous material Substances 0.000 claims 1
- 238000005553 drilling Methods 0.000 abstract description 8
- 201000004569 Blindness Diseases 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract description 3
- 230000003993 interaction Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种重力流水道砂体中优质储层的预测方法,涉及石油地质工程领域。该方法包括获取多个单井的重力流水道砂体的测井资料,得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度;并进一步得到每个所述单井的重力流水道砂体对应的波阻抗;在此基础上建立孔隙度‑波阻抗交互图版,并在该图版上确定重力流水道砂体中的优质储层;对目标重力流水道砂体储层进行波阻抗反演,结合所述孔隙度‑波阻抗交互图版所确定的所述重力流水道砂体的优质储层信息,得到所述目标重力流水道砂体储层的定量分布。本发明实施例提供的预测方法可定量预测重力流水道砂体优质储层的分布,从而克服了钻井的盲目性,提高钻井成功率。
Description
技术领域
本发明属于石油地质工程领域,特别涉及一种重力流水道砂体中优质储层的预测方法。
背景技术
在中国东部陆相盆地中,大面积分布着重力流水道砂体。由于重力流是沉积物-流体混合物在重力作用下整体顺斜坡向下运动、以悬移方式搬运为主的高密度流体流动,其形成属于事件性沉积作用,在洪水、地震、海啸、火山喷发等阵发性因素直接或间接诱发下导致沉积物的滑动形成重力流沉积。重力流水道砂体中碎屑颗粒大小混杂,分选差,储层物性差异大,如何应用地震预测优质储层难度较大。
目前国内外石油地质专家一般采用地震资料中的地震属性预测重力流水道砂体储层的分布,对重力流水道砂体储层进行定性描述和分析。
在实现本发明的过程中,本发明人发现现有技术中至少存在以下问题:
现有的技术只能区分出储层,不能定量的区别优质储层和差储层,同时,目前还没有一种用来预测重力流水道砂体中优质储层的有效方法,钻井作业比较盲目,钻井成功率低。
发明内容
鉴于此,本发明提供一种能够定量预测重力流水道砂体中优质储层的预测方法。
具体而言,包括以下的技术方案:
一种重力流水道砂体中优质储层的预测方法,所述方法包括:
获取多个单井的重力流水道砂体的测井资料,得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度;
将每个所述单井的所述岩石密度和所述纵波速度相乘,得到每个所述单井的重力流水道砂体对应的波阻抗;
根据多个所述单井的重力流水道砂体对应的所述孔隙度和所述波阻抗,建立孔隙度-波阻抗交互图版,并在该图版上确定重力流水道砂体中的优质储层;
对目标重力流水道砂体储层进行波阻抗反演,并结合所述孔隙度-波阻抗交互图版所确定的所述重力流水道砂体的优质储层信息,得到所述目标重力流水道砂体储层的定量分布。
进一步地,所述根据多个所述单井的重力流水道砂体对应的所述孔隙度和所述波阻抗,建立孔隙度-波阻抗交互图版,并在该图版上确定重力流水道砂体中的优质储层包括:根据所述孔隙度-波阻抗交互图版,分别选用主频对所述单井进行正演,来确定所述重力流水道砂体中的优质储层。
优选地,所选用的所述主频为20Hz、30Hz、40Hz。
进一步地,所述获取单井的重力流水道砂体的测井资料,得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度中,所述测井资料包括声波时差。
优选地,所述孔隙度的计算方式为:Φ=(Δt-Δtma)/[(Δtf-Δtma)*Cp],其中:Φ为孔隙度;Δt为声波时差;Δtma为岩石骨架声波时差值;Δtf为岩石孔隙流体声波时差值;Cp为压实系数。
优选地,所述孔隙度的计算方式为:Φ=(ρma-ρb)/(ρma-ρf),其中:Φ为孔隙度;ρma是岩石密度骨架值;ρb是岩石密度测量值;ρf是岩石孔隙流体密度值。
进一步地,所述纵波速度的计算方式为:V=1×106/△t,其中:V为纵波速度。
进一步地,所述孔隙度-波阻抗交互图版包括:当孔隙度大于12%时,波阻抗值小于9900g.cm-3.m.s-1;当孔隙度小于8%时,波阻抗值大于11600g.cm-3.m.s-1。
优选地,所述对目标重力流水道砂体储层进行波阻抗反演,并结合所述孔隙度-波阻抗交互图版所确定的重力流水道砂体的优质储层信息,得到所述目标重力流水道砂体储层的定量分布包括:
在所述目标重力流水道砂体储层的地震数据体中解释所述目标重力流水道砂体储层的顶底界面;
提取所述顶底界面内的地震属性,定性预测所述目标重力流水道砂体储层中优质储层的分布情况;
根据所述分布情况,进行波阻抗约束反演,并结合所述孔隙度-波阻抗交互图版,定量预测所述目标重力流水道砂体储层中优质储层的分布。
本发明实施例提供的技术方案的有益技术效果如下:
本发明实施例提供了一种能够定量预测重力流水道砂体中优质储层的预测方法。具体地,首先通过运用单井的重力流水道砂体储层的测井资料,得到单井的重力流水道砂体对应的波阻抗,建立孔隙度-波阻抗交互图版,并确定重力流水道砂体中的优质储层,再进行井震结合,通过目标重力流水道砂体储层的地震属性和测井约束反演进行优质储层的定性与定量预测。本发明实施例提供的预测方法可定量预测重力流水道砂体优质储层的分布,克服了钻井的盲目性,提高钻井成功率。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种重力流水道砂体中优质储层的预测方法的方法流程图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。除非另有定义,本发明实施例所用的所有技术术语均具有与本领域技术人员通常理解的相同的含义。
本发明提供一种重力流水道砂体优质储层的预测方法,参见图1,该方法包括如下步骤:
步骤101:获取多个单井的重力流水道砂体的测井资料,得到每个单井的重力流水道砂体对应的孔隙度Φ、岩石密度和纵波速度V;
根据已钻单井的测井资料(包括声波时差△t等),求取各单井中重力流水道砂体对应的孔隙度Φ和纵波速度V,优选地,可以用两种方式计算孔隙度Φ:
一种是采用声波时差孔隙度的计算方式:
Φ=(Δt-Δtma)/[(Δtf-Δtma)*Cp],其中:Δt为声波时差;Δtma为岩石骨架声波时差值;Δtf为岩石孔隙流体声波时差值;Cp为压实系数。
另一种是采用密度测井孔隙度计算方式:
Φ=(ρma-ρb)/(ρma-ρf),其中:ρma是岩石密度骨架值;ρb是岩石密度测量值;ρf是岩石孔隙流体密度值。
纵波速度V=1×106/Δt,声波时差Δt单位为μs/m。
由此,根据多个单井的重力流水道砂体储层的孔隙度Φ和纵波速度V的值,可建立反映孔隙度与纵波速度的对应关系的地质模型,如下:
孔隙度Φ大于12%的重力流水道砂体对应低速层,纵波速度V=4000-4350m.s-1;孔隙度Φ小于8%的重力流水道砂体对应高速层,纵波速度V=5400-5500m.s-1;即优质储层对应低速层,差储层对应高速层。
步骤102:将每个单井的岩石密度和纵波速度相乘,得到每个单井的重力流水道砂体对应的波阻抗;
用岩石密度的测量值ρb与纵波速度V相乘,求取各单井重力流水道砂体对应的波阻抗(impedanee)。
步骤103:根据多个单井的重力流水道砂体对应的孔隙度Φ和波阻抗,建立孔隙度-波阻抗交互图版,并在该图版上确定重力流水道砂体中的优质储层;
根据多个单井的重力流水道砂体储层的孔隙度Φ和波阻抗的值,可建立反映孔隙度与波阻抗的线性关系交互图版,由关系图版可得出,孔隙度Φ大于12%的中高孔隙度区,波阻抗值小于9900g.cm-3.m.s-1;孔隙度Φ小于8%的致密区,波阻抗值大于11600g.cm- 3.m.s-1,从而可知,中低波阻抗对应的砂体储层物性较好,为优质储层,高阻抗对应的砂体储层物性相对较差,为致密砂岩发育区。
进一步地,在孔隙度-波阻抗交互图版的基础上,分别选用主频对各单井进行正演,以进一步确定重力流水道砂体中的优质储层。结果表明,在厚度相近的情况下,中波阻抗对应的砂体储层物性较好,高阻抗对应的砂体储层物性相对较差。
且由于在地震采集时,主要目的层的主频一般在20-30Hz之间,而提高主频可以更精确的推测砂体积尖灭点的位置,但是太高了又容易造成假象,因此优选地,用20Hz、30Hz和40Hz进行正演。
步骤104:对目标重力流水道砂体储层进行波阻抗反演,并结合孔隙度-波阻抗交互图版所确定的重力流水道砂体的优质储层信息,得到目标重力流水道砂体储层的定量分布。
具体地,在目标重力流水道砂体储层的地震数据体中解释目标重力流水道砂体储层的顶底界面;
提取顶底界面内的地震属性,定性预测目标重力流水道砂体储层中优质储层的分布情况;
根据分布情况,进行波阻抗约束反演,并结合孔隙度-波阻抗交互图版,定量预测所述目标重力流水道砂体储层中优质储层的分布。
本发明实施例提供的预测方法可定量预测重力流水道砂体优质储层的分布,经5口井钻探实施后,优质储层钻遇率由原来的50%提高到100%,单井日产油由原来的0.01-5.2吨提高到10.5--40.8吨,克服了钻井的盲目性,提高钻井成功率。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案,并不用以限制本发明。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种重力流水道砂体中优质储层的预测方法,其特征在于,所述方法包括:
获取多个单井的重力流水道砂体的测井资料,得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度;
将每个所述单井的所述岩石密度和所述纵波速度相乘,得到每个所述单井的重力流水道砂体对应的波阻抗;
根据多个所述单井的重力流水道砂体对应的所述孔隙度和所述波阻抗,建立孔隙度-波阻抗交互图版,并在该图版上确定重力流水道砂体中的优质储层;
对目标重力流水道砂体储层进行波阻抗反演,并结合所述孔隙度-波阻抗交互图版所确定的所述重力流水道砂体的优质储层信息,得到所述目标重力流水道砂体储层的定量分布。
2.根据权利要求1所述的方法,其特征在于,所述根据多个所述单井的重力流水道砂体对应的所述孔隙度和所述波阻抗,建立孔隙度-波阻抗交互图版,并在该图版上确定重力流水道砂体中的优质储层包括:根据所述孔隙度-波阻抗交互图版,分别选用主频对所述单井进行正演,来确定所述重力流水道砂体中的优质储层。
3.根据权利要求2所述的方法,其特征在于,所选用的所述主频为20Hz、30Hz、40Hz。
4.根据权利要求1所述的方法,其特征在于,所述获取单井的重力流水道砂体的测井资料,得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度中,所述测井资料包括声波时差。
5.根据权利要求4所述的方法,其特征在于,所述孔隙度的计算方式为:Φ=(Δt-Δtma)/[(Δtf-Δtma)*Cp],其中:Φ为孔隙度;Δt为声波时差;Δtma为岩石骨架声波时差值;Δtf为岩石孔隙流体声波时差值;Cp为压实系数。
6.根据权利要求4所述的方法,其特征在于,所述孔隙度的计算方式为:Φ=(ρma-ρb)/(ρma-ρf),其中:Φ为孔隙度;ρma是岩石密度骨架值;ρb是岩石密度测量值;ρf是岩石孔隙流体密度值。
7.根据权利要求4所述的方法,其特征在于,所述纵波速度的计算方式为:V=1×106/△t,其中:V为纵波速度。
8.根据权利要求1所述的方法,其特征在于,所述孔隙度-波阻抗交互图版包括:当孔隙度大于12%时,波阻抗值小于9900g.cm-3.m.s-1;当孔隙度小于8%时,波阻抗值大于11600g.cm-3.m.s-1。
9.根据权利要求1所述的方法,其特征在于,所述对目标重力流水道砂体储层进行波阻抗反演,并结合所述孔隙度-波阻抗交互图版所确定的重力流水道砂体的优质储层信息,得到所述目标重力流水道砂体储层的定量分布包括:
在所述目标重力流水道砂体储层的地震数据体中解释所述目标重力流水道砂体储层的顶底界面;
提取所述顶底界面内的地震属性,定性预测所述目标重力流水道砂体储层中优质储层的分布情况;
根据所述分布情况,进行波阻抗约束反演,并结合所述孔隙度-波阻抗交互图版,定量预测所述目标重力流水道砂体储层中优质储层的分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111926.1A CN106772614A (zh) | 2017-02-28 | 2017-02-28 | 一种重力流水道砂体中优质储层的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111926.1A CN106772614A (zh) | 2017-02-28 | 2017-02-28 | 一种重力流水道砂体中优质储层的预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106772614A true CN106772614A (zh) | 2017-05-31 |
Family
ID=58959446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710111926.1A Pending CN106772614A (zh) | 2017-02-28 | 2017-02-28 | 一种重力流水道砂体中优质储层的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772614A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107165627A (zh) * | 2017-07-18 | 2017-09-15 | 中国石油大学(华东) | 一种坝砂有利储层分布范围预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454685A (zh) * | 2013-08-09 | 2013-12-18 | 中国石油天然气股份有限公司 | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 |
CN105089632A (zh) * | 2015-08-04 | 2015-11-25 | 中国海洋石油总公司 | 一种高温高压储层co2流体纵波时差骨架参数的获取方法 |
CN106368691A (zh) * | 2015-07-24 | 2017-02-01 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
-
2017
- 2017-02-28 CN CN201710111926.1A patent/CN106772614A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454685A (zh) * | 2013-08-09 | 2013-12-18 | 中国石油天然气股份有限公司 | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 |
CN106368691A (zh) * | 2015-07-24 | 2017-02-01 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
CN105089632A (zh) * | 2015-08-04 | 2015-11-25 | 中国海洋石油总公司 | 一种高温高压储层co2流体纵波时差骨架参数的获取方法 |
Non-Patent Citations (5)
Title |
---|
PU XIUGANG ET AL.: "Gravity flow sedimentation and tight oil exploration in lower first member of Shahejie Formation in slope area of Qikou Sag, Bohai Bay Basin", 《PETROLEUM EXPLORATION AND DEVELOPMENT》 * |
孔德博等: "歧南斜坡沙一段中部重力流水道优质砂体分布与演化", 《成都理工大学学报(自然科学版)》 * |
宋春刚: "G油田古近系重力流水道相储层预测研究", 《中外能源》 * |
赵忠军等: "基于井校地震属性层序切片砂体预测技术研究", 《第十一届宁夏青年科学家论坛论文集》 * |
陈楠: "塔中地区志留系柯坪塔格组沉积相与砂体展布规律", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107165627A (zh) * | 2017-07-18 | 2017-09-15 | 中国石油大学(华东) | 一种坝砂有利储层分布范围预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Satter et al. | Reservoir engineering: the fundamentals, simulation, and management of conventional and unconventional recoveries | |
CN105469159B (zh) | 定量预测油气有利聚集区的方法 | |
CN104912550B (zh) | 一种核磁共振测井资料定量计算储层产液剖面方法 | |
CN109441422A (zh) | 一种页岩气井间距优化开采方法 | |
CN108868748A (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
CN105467464A (zh) | 近岸水下扇砂砾岩有效连通体划分和对比方法 | |
CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
CN105467466B (zh) | 一种基于多尺度信息的致密储层成岩相预测方法 | |
Zeng et al. | The effect of multi-scale faults and fractures on oil enrichment and production in tight sandstone reservoirs: A case study in the southwestern Ordos Basin, China | |
Wang et al. | Technologies for the benefit development of low-permeability tight sandstone gas reservoirs in the Yan'an Gas Field, Ordos Basin | |
Delorme et al. | An efficient model to simulate reservoir stimulation and induced microseismic events on 3D discrete fracture network for unconventional reservoirs | |
Lothe et al. | Storage resources for future European CCS deployment; A roadmap for a Horda CO2 storage hub, offshore Norway | |
Yang et al. | Numerical investigation on the spewing mechanism of earth pressure balance shield in a high‐pressure water‐rich sand stratum | |
Giamminonni et al. | Condensate banking phenomenon evaluation in heterogeneous low permeability reservoirs | |
CN106772614A (zh) | 一种重力流水道砂体中优质储层的预测方法 | |
Suarez et al. | Fracturing-to-Production Simulation Approach for Completion Optimization in the Vaca Muerta Shale | |
Burns et al. | The challenges of shale gas exploration and appraisal in Europe and North Africa | |
Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
Chiarella et al. | Stratigraphic and sedimentological characterization of the Late Cretaceous post-rift intra Lange Sandstones of the Gimsan Basin and Grinda Graben (Halten Terrace, Norwegian Sea) | |
CN106646594A (zh) | 基于“断‑砂”耦合关系的浅层岩性油气藏定量预测方法 | |
CN115407402A (zh) | 一种多因子联合砂砾岩优质储层预测方法 | |
Tan et al. | A new methodology for modelling of sand wormholes in field scale thermal simulation | |
Abbaszadeh et al. | Development of an integrated reservoir model for a naturally fractured volcanic reservoir in China | |
Carr et al. | Characterization of a Vertically Compartmentalized Reservoir in a Supergiant Field, Burgan Formation, Greater Burgan Field, Kuwait; Part 2: Evolving Concepts of Reservoir Performance | |
Rodriguez et al. | Application of modern reservoir characterization in mature fields to unravel hidden reserves |
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: 20170531 |