CN103452547B - 试井资料中续流数据的分析处理方法和系统 - Google Patents
试井资料中续流数据的分析处理方法和系统 Download PDFInfo
- Publication number
- CN103452547B CN103452547B CN201210281796.3A CN201210281796A CN103452547B CN 103452547 B CN103452547 B CN 103452547B CN 201210281796 A CN201210281796 A CN 201210281796A CN 103452547 B CN103452547 B CN 103452547B
- Authority
- CN
- China
- Prior art keywords
- speed
- time
- reservoir
- during
- curve
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 189
- 238000004458 analytical method Methods 0.000 title claims abstract description 101
- 238000003672 processing method Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 116
- 230000008569 process Effects 0.000 claims abstract description 61
- 238000005259 measurement Methods 0.000 claims abstract description 52
- 239000012530 fluid Substances 0.000 claims abstract description 50
- 238000004519 manufacturing process Methods 0.000 claims abstract description 34
- 230000035699 permeability Effects 0.000 claims description 50
- 230000015572 biosynthetic process Effects 0.000 claims description 20
- 230000006399 behavior Effects 0.000 claims description 12
- 238000007405 data analysis Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 abstract description 14
- 238000011156 evaluation Methods 0.000 description 28
- 230000018109 developmental process Effects 0.000 description 25
- 238000011161 development Methods 0.000 description 24
- 238000005516 engineering process Methods 0.000 description 21
- 238000011160 research Methods 0.000 description 21
- 230000008859 change Effects 0.000 description 20
- 238000010586 diagram Methods 0.000 description 10
- 239000003921 oil Substances 0.000 description 9
- 238000003860 storage Methods 0.000 description 7
- 239000011435 rock Substances 0.000 description 5
- 238000011109 contamination Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 241001269238 Data Species 0.000 description 2
- 238000003556 assay Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000008034 disappearance Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009969 flowable effect Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 102000003712 Complement factor B Human genes 0.000 description 1
- 108090000056 Complement factor B Proteins 0.000 description 1
- 230000005483 Hooke's law Effects 0.000 description 1
- 206010027336 Menstruation delayed Diseases 0.000 description 1
- 239000003637 basic solution Substances 0.000 description 1
- 239000010430 carbonatite Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011157 data evaluation Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- -1 lava Substances 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000000246 remedial effect Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000000153 supplemental effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Chemical & Material Sciences (AREA)
- Geophysics (AREA)
- Dispersion Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明实施例公开了一种试井资料中续流数据的分析处理方法和系统,其中,方法包括:预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型;在对储层试井过程中,在不同采集时刻分别采集井底压力,并获取等时速率VP与等时时间间隔t之间的实测等时速率曲线;选取试井解释模型参数,根据等时速率数学模型和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取VP’与模型等时时间间隔t’之间的模型等时速率曲线;将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数。本发明实施例可以实现对试井资料中续流数据的分析处理。
Description
技术领域
本发明涉及油田勘探开发中试井资料的评价分析处理技术,尤其是一种试井资料中续流数据的分析处理方法和系统。
背景技术
试井(也称为:测试)是油田勘探与开发过程中,及时获取油藏特征信息和对油气层进行定量评价分析的重要手段,其成果是进行油田的储量计算、油藏模型、制定勘探方案、开发方案和措施优选所必需的储层特征参数资料。
试井解释是对地层测试与试井资料评价分析处理的关键环节,试井工艺取得的资料都要应用试井解释理论与方法,通过对试井资料的综合评价分析处理求取储层特性参数,即:储层的各项特征参数,是一种以线性流渗流理论为基础的解释评价技术。试井解释理论与方法在漫长的发展历程中,大致经历了如下三个阶段:以压降解释理论为基础的霍纳解释方法和格林卡登图版解释方法为主体的手工解释时代,以霍纳解释理论与压降解释理论为基础的计算机模型解释阶段,和以叠加解释理论和现代试井解释理论为基础的多模型解释阶段。
从计算机技术应用于试井解释的整个发展过程中,都是以建立最新试井解释理论研究成果与计算机技术相结合的专业技术应用系统为发展主体,形成不同时期的试井解释理论与方法,实现对压降资料或压力恢复资料的高精度分析处理,但在整体的发展过程中,试井解释理论与方法都是建立于达西渗流理论与径向流理论基础上,并在此理论基础上,针对不同的储层特征,通过各类物理模型与数学模型的建立,形成了多种试井解释图版。
随着试井工艺技术、测量仪器和解释理论的不断完善和计算机技术的高速发展,试井的工艺技术方法越来越多,*获取的测试资料精度越来越高,试井解释评价分析的模型范围越来越多,建立起了一套以稳定试井和不稳定试井为核心内容的现代试井理论与方法的体系。但不论是何种试井工艺,其应用的理论基础仍然是达西渗流理论与径向流试井基础理论体系,在对测试资料分析处理和获取储层特性参数的工作中,无论是半对解释方法、改进的混合分布假说(MixtureOfDistributionHypothesis,以下简称:MDH)法或者是现代试井解释方法,其基础主要是研究径向流条件下压力的变化特征,并通过对压力变化特征的研究建立不同储层的基础模型,通过对基础模型与实测数据进行对比分析,求解测试储层的压力、渗透率、外界干扰和损害程度等特性参数。由于受到基础理论与模型建立受到假设条件的制约,上述方法的应用基础是获取的压力数据以达到消除井筒干扰和储层续流干扰为基础时的形态----径向流状态下实现求解。但在实际的应用过程中,随着油气勘探开发向低渗、特低渗储层的发展,一方面受非达西流渗流影响,使建立于径向流理论基础上的模型应用受到限制,试井解释误差增大;另一方面由于测试时间不合理、储层物性低渗等原因,导致储层导流能力差、生产过程产出少、压力传导慢等问题突出,使得测试过程中无法测取达到径向流条件下的试井数据,要获取符合径向流试井解释理论要求的资料难度极大。特别是在特低渗储层测试中,虽然通过大量的工艺技术改进、优化测试时间等措施,测试资料品质只能得到一定的改善,但取得符合径向流试井解释理论的测试资料仍然很低,大量的测试资料不能有效的进行储层评价,造成储层评价分析第一手资料缺失,使勘探开发和储层改造的针对性降低,导致测试成本和勘探开发大幅度增加,低产出和高投入的问题突出。
在试井过程中,实际测试的压力数据按照阶段从前至后依次分为井筒储集数据、续流数据(也称为过渡段数据)、径向流数据和晚期数据四部分。由于储层物性较差或实际测试时间分配不合理的原因,在低渗透和特低渗透储层的试井过程中,总有大量的测试数据达不到径向流试井解释理论条件的要求,所取得的资料处于续流阶段,通称为续流数据,这将导致应用传统的径向流试井解释理论与方法,无法对测试的续流数据进行评价分析,无法通过试井数据和资料的评价分析获取储层特性参数,导致试井资料失去应用价值。由于不能取得试井解释评价成果,无法对储层进行定量与定性评价分析,使得油田勘探与开发综合分析资料缺失,给油田的勘探开发方案制定带来难度。特别是对致密岩性储层勘探开发工作展开后,在对油田储层的评价分析工作中,试井工艺无论采取什么测试方法和补救措施,在实际的地层测试过程中,储层渗流能力低,导压性能差,要取得径向流段数据的难度极大,不仅需要大量的测试费用,而且需要很长的测试时间,大部分储层仍然无法达到获取径向流试井数据的要求。
由此,在实现本发明的过程中,发明人发现现有技术至少存在以下问题:
以径向流渗流理论为基础的试井解释理论与方法无法对续流数据进行评价分析,也没有形成相应的理论基础与评价分析方法,在长期的应用研究中,研究人员应用生产史叠加分析方法、力学校正法、早期续流时间校正方法和单井数学模型等方法,通过对细分流量、关井时间校正、生产时间补偿等,都离不开将续流数据通过一定的经验校正或力学校正,使不符合径向流理论的数据,通过校正后应用径向流理论进行分析处理,由于其理论基础的应用就存在一定的问题,而由此而建立的解释方法复杂且具极大的局限性,仅能作为经验性对比研究和个性分析评价中参考性的分析,无法满足试井资料的评价分析和矿场应用要求。
因此,应用径向流试井理论与方法无法系统全面的进行续流数据的处理和评价分析,以通过试井为手段而获取储层的特性参数资料缺失成为制约油田勘探开发的主导因素。同时,随着致密岩性、裂缝性火山岩等非达西流储层的不断增多,由于其渗流状态完全不符合径向流理论,采用以径向流试井解释理论与方法处理评价的精度降低,储层特征分析难度增大,解决续流资料和非达西流储层的评价分析成为国内外近年来重点的技术攻关研究的内容。
发明内容
本发明实施例所要解决的一个技术问题是:提供一种试井资料中续流数据的分析处理方法和系统,针对在试井过程中取得的实测资料达不到径向流条件时,实现对试井资料中续流数据的分析处理,获取储层的特征参数,以实现对非达西流储层的评价分析。
本发明实施例提供的一种试井资料中续流数据的分析处理方法,包括:
预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,所述等时速率与各种影响因素之间的关系包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系;
根据等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型;
在对储层试井过程中,在不同采集时刻分别采集井底压力,并获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线;其中,P2与P1分别为等时时间间隔t的结束时刻t2与开始时刻t1采集到的压力值;
选取试井解释模型参数,根据建立的等时速率数学模型和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取VP’与模型等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线;其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率;
将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,所述储层特征参数包括渗透率K与窜流时间tc。
本发明实施例提供的一种试井资料中续流数据的分析处理系统,包括:
采集单元,用于在对储层试井过程中,在不同采集时刻分别采集井底压力;
参数选取单元,用于选取试井解释模型参数;
存储单元,用于存储预先建立的等时速率数学模型;所述等时速率数学模型基于预先测定的各种岩性条件下的等时速率与各种影响因素之间的关系获取得到的等时速率与储层岩性、储层结构、流体特性之间的关系数据建立,所述等时速率与各种影响因素之间的关系包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系;
计算单元,用于根据参数选取单元选取的试井解释模型参数和存储单元中的等时速率数学模型,计算在不同时刻的模型等时速率VP’;其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率;
获取单元,用于根据采集单元在不同采集时刻分别采集井底压力,获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线;其中,P2与P1分别为等时时间间隔t的结束时刻t2与开始时刻t1采集到的压力值;以及获取计算单元计算出的模型等时速率VP’与模型等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线;
拟合单元,用于将获取单元获取的模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,所述储层特征参数包括渗透率K与窜流时间tc。
基于本发明上述实施例提供的试井资料中续流数据的分析处理方法和系统,以研究储层在相等时间内压力的变化速率为基础,预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系;根据等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型。基于等时速率解释理论与方法,在试井过程中在井筒储层阶段完成后,所取得的压力恢复或压力降落资料未达到稳定的径向流状态条件时,在井筒储层阶段完成后,在不同采集时刻分别采集井底压力,并获取表示等时速率VP与等时时间间隔t之间关系的实测等时速率曲线,并根据预先建立的等时速率数学模型和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取表示VP’与模型等时时间间隔t之间的对应关系的模型等时速率曲线,然后将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,包括渗透率K与窜流时间tc,从而实现了对续流数据的分析处理,可应用于全面分析测试半径范围内储层的渗流特征,实现对续流资料的解释处理,取得储层的特征参数,达到同径向流资料相同的分析处理之目的,以实现对非达西流储层的评价分析。本发明实施例是对续流数据进行评价分析的全新的续流解释理论与方法,是通过研究压力传导过程中的等时速率为分析对象的新型解释理论与方法。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明试井资料中续流数据的分析处理方法一个实施例的流程图。
图2为本发明实施例中等时速率与岩性关系的一个示意图。
图3为本发明实施例中等时速率与渗透率的关系的一个示意图。
图4为本发明实施例中等时速率与污染特征的关系的一个示意图。
图5为本发明实施例中等时速率与有限导流裂缝干扰的关系的一个示意图。
图6为本发明实施例中等时速率与无限导流裂缝干扰的关系的一个示意图。
图7为本发明方法试井资料中续流数据的分析处理另一个实施例的流程图。
图8为本发明实施例中VP~t曲线的一个具体示例。
图9为本发明实施例中第一曲线图的一个具体示例。
图10为本发明实施例中第二曲线图的一个具体示例。
图11为本发明实施例中实测P~t曲线的一个具体示例。
图12为本发明试井资料中续流数据的分析处理系统一个实施例的结构示意图。
图13为本发明试井资料中续流数据的分析处理系统另一个实施例的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
国外以线性流渗流理论为基础理论的研究相对较早,早期就形成了以砂岩储层、复合储层、裂缝性储层为主体的多种渗流解释理论,并形成了美国信息及通信技术(ICT)公司的试井解释软件、江斯顿的试井系统和美国SSI公司的workbench试井解释软件等。在后期的发展中,美国SSI公司不断完善试井解释理论与方法,并结合生产应用需求和计算机网络技术,研发成为集试井设计、试井解释和油藏数字模型为主体的试井网络应用系统和解释工作站,试井解释技术从基础理论研究到软件开发得到了快速的发展。90年代开始,英国的EPS公司和法国Kappa公司的SaPhir试井解释软件逐步以强劲的试井基础理论和解释方法研究成果占据了世界试井解释的较大空间。并通过不断的解释理论创新和解释方法个性化研究,紧跟计算机网络技术与矿场应用需要,将试井解释、试井设计和油藏数字模型技术可视化、人性化,不断的完善试井解释模型的数量,模型也由单一的均质储层向裂缝性储层、复合储层发展,井类由单一的直井,向水平井、斜井、多支井发展,并建立了气井、水井、油井、稠油井等多种复杂的解释理论模型和多类型边界影响模型,使线性流基础理论为解释基础的试井解释技术发展到较完善的水平,系统地指导石油勘探与开发的地质评价。
近年以来,随着非线性流测试资料和续流测试资料的不断增多,研究工作注重于以麦金利续流理论的深度研究与应用,建立了数字试井解释方法,并进一步完善其研究成果,但受到多重影响因素的影响,技术研究取得的突破较小,非线性流资料与续流资料解释系统的发展与解释精度仍然不能达到高精度分析的需要。
早期国内的试井解释工作主要依赖于引进国外的解释图版或系统软件,80年代由华北测试公司通过引进消化美国ICT公司的试井解释软件,帅先研制开发了WTC试井解释软件,并随着计算机技术的快速发展,研制开发了super-WTIS试井解释软件,成为国内试井解释主导应用软件。随着国内勘探开发储层类型的复杂和工艺技术的进步,国产软件模型单一和功能低下的问题表现出较大的不适用性,很快受到国外先进试井解释软件的冲击而退出应用领域。针对国内油田以低渗、特低渗为主的现状和勘探开发向低渗透火山岩储层的发展,非达西流储层和低渗透续流资料的不断增多,国内外勘探开发对象的完全不同,引进的试井解释软件系统和国外的技术已经不能满足国内的生产需要,使国内在非线性流基础理论研究工作启动较早,解释基础理论与方法的研究领先于国外的发展。但由于试井解释软件系统的研究工作长期以径向流渗流理论为基础,导致试井解释技术无法取得较大的突破,应用条件受限、精度低和分析评价应用难度大的问题突出,导致对续流数据的分析评价问题没有得到根本的解决。
目前以径向流理论为基础对续流数据的评价分析的方法,是基于径向流压力变化规律研究而形成的理论,以此理论为基础的续流资料分析方法,主要是对关井时间或产量进行叠加分析而使续流数据通过补偿转化而形成径向流段数据后,应用径向流理论进行分析处理。由于其补偿转化的方法和多个未知参数影响,建立系统的解释模型与方法难度高、误差大。而本发明实施例采用等时速率解释理论与方法,以储层在相等时间内压力的变化速率为基础,从根本上建立了自己独立的理论和方法,并将这一理论与传统的达西渗流理论相结合,以实际试验结果进行参数回归转化。
本发明的发明人在长期的研究工作中,经过对压力变化规律性特征进行的大量室内试验、矿场试验与对试井资料的数学分析研究发现,储层流体的压力的传导速率是储层特征的真实反映,压力传导速率在井筒储集阶段完成后,规律性变化特征已经形成,可以表述为:
VPL=(P2-P1)/(t2-t1)(1)
上述公式(1)中,P1表示储层流体在时间【t1,t2】的开始时刻t1的压力,P2表示储层流体在时间【t1,t2】的结束时刻t2的压力;VPL表示压力的传导速率(以下简称:压力速率)。
根据上述试验结果,本发明实施例针对不同岩性条件的岩矿进行了大量的渗流试验和导压特性试验,证明岩石的压力速率随储层特征的变化而变化,压力速率的变化特征在达西流和非达西流状态下都具有规律性变化的特征,流动相态的变化反映到压力速率特征上总是表现为压力速率的增大与减小。而将储层压力速率按同一时间间隔进行处理后建立的等时速率,其规律性变化出现的时间,远早于径向流态,在井筒储层阶段完成后,等时速率就形成了规律的变化特征,可以表述为:
VP=(P2-P1)/t(2)
上述公式(2)中,VP表示储层流体在等时时间间隔t内的压力速率(以下简称:等时速率),可以真实反映储层特征;t可以表示为t2-t1,P1表示储层流体在等时时间间隔t的开始时刻t1的压力,P2表示储层流体在等时时间间隔t的结束时刻t2的压力。上述公式(2)所表示的理论在本发明中可以称为等时速率解释理论,是试井解释基础理论与方法的研究成果,可以反映储层渗流特征。
等时速率早期规律性变化特征和随流动相态的规律性变化特征,使应用等时速率曲线的早期线性回归求解得到实现,也使非达西流的针对性分析得到理论支持。本发明将等时速率变化规律的研究与储层特征参数研究相结合,应用等时速率曲线规律性特征进行评价分析,并以此为模型分析计算储层特性参数,从而取得储层特性参数,达到了以续流的压力速率变化进行续流数据的评价分析,将规律性的压力速率变化量转化为储层渗透率,并求解其它的储层特性参数。等时速率解释理论成为对续流数据和非达西流储层针对性评价分析处理的理论基础。
图1为本发明试井资料中续流数据的分析处理方法一个实施例的流程图。如图1所示,该实施例的试井资料中续流数据的分析处理方法包括:
步骤101,预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系。
步骤102,根据等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型。
步骤103,在对储层试井过程中,例如,在井筒储层阶段完成后,在不同采集时刻分别采集井底压力,并获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线。
其中,P2与P1分别为等时时间间隔t的结束时刻t2与开始时刻t1采集到的压力值。
步骤104,选取与实际测试等时速率曲线相符合的试井解释模型参数,例如,与实测等时速率曲线较为接近的平均生产产量Q、地层流体体积系数B、地层流体粘度u、层流体的渗透率K、窜流时间tc、生产时间tP、窜流系数λ,根据建立的等时速率数学模型和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取VP’与等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线。
其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率。
步骤105,将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,该储层特征参数包括渗透率K与窜流时间tc。
本发明上述实施例提供的试井资料中续流数据的分析处理方法,以研究储层在相等时间内压力的变化速率为基础,预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系;根据等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型。基于等时速率解释理论与方法,在试井过程中在井筒储层阶段完成后,所取得的压力恢复或压力降落资料未达到稳定的径向流状态条件时,在井筒储层阶段完成后,在不同采集时刻分别采集井底压力,并获取表示等时速率VP与等时时间间隔t之间关系的实测等时速率曲线,并根据预先建立的等时速率数学模型和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取表示VP’与模型等时时间间隔t之间的对应关系的模型等时速率曲线,然后将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,包括渗透率K与窜流时间tc,从而实现了对续流数据的分析处理,可应用于全面分析测试半径范围内储层的渗流特征,实现对续流资料的解释处理,取得储层的特征参数。
本发明实施例将等时速率解释理论用于续流数据进行评价分析时,针对等时速率变化特征与储层特征的对应关系,依据储层的渗流特征建立典型的等时速率曲线,并建立等时速率与各种影响因素之间的关系,最终导出等时速率数学模型。通过对实测等时速率曲线与典型等时速率曲线的拟合对比分析,获取实测储层的各项特征参数数据。本发明人通过大量的试验室岩石渗流试验,全面测定了各种岩性条件下的等时速率与各种影响因素之间的关系,从而获得等时速率与储层的岩性、储层结构、流体特性的关系数据λQBu/K。其中,λ表示地层流体的窜流系数,为预设大于零的常数;Q表示平均生产产量;B表示地层流体体积系数;u表示粘度;K表示地层流体的渗透率,均为无因次物理量。如下图2~图6所示,依次为各种岩性条件下,等时速率与储层的岩性特征的关系示意图,等时速率与渗透率的关系示意图,等时速率与污染特性的关系示意图,等时速率与有限导流裂缝干扰的关系示意图,等时速率与无限导流裂缝干扰的关系示意图,其中,图5与图6所示的有限导流裂缝干扰与无限导流裂缝干扰属于裂缝窜流特征。图2~图6中的横坐标表示时间,纵坐标表示等时速率。
本发明实施例中,通过储层岩性特征、渗透率、污染特征、裂缝窜流特征对等时速率的影响的分析,建立了如下等时速率数学模型:
Vp=λQBu【lg(((t2(tp+t1))/((t1(tp+t2)))】/(t2-t1)/kh(4)
上述公式(3)和(4)中,VP表示模型等时速率;t1表示等时时间间隔t’的开始时刻t1’;t2表示等时时间间隔t’的结束时刻t2’;tP表示生产时间;Q表示平均生产产量,为预知取值;B表示地层流体体积系数,为预知取值;u表示地层流体粘度,为预知取值;K表示地层流体的渗透率;h表示储层厚度,可以预先测试获得;λ表示地层流体的窜流系数,为预设大于零的校正常数,与储层的岩性特征与导流特征具有相关性。上述物理量均为无因次物理量。其中,导流特征是岩性、流体特征条件下的一种导流方式的总称。
根据本发明试井资料中续流数据的分析处理方法的一个具体示例而非限制,窜流系数λ的取值根据储层的介质类型以及窜流时间tC预设为1.81~2.59。
本发明实施例中,应用等时速率解释理论对续流数据和非达西流数据进行评价分析处理时,将储层分为单介质渗流(均质)储层和多介质渗流(非均质)储层两大类。
对于单介质渗流储层,具有井筒储集特征、均质、等厚、无限大储层特征,除层的岩石弹性压缩时符合虎克定律,储层各点的压力均衡,生产产量恒定。本发明实施例中提出的单介质渗流储层等时速率数学模型为:
定解条件为:VP(r,t-ts)=VPo,VP(ω,t-ts)=0(7)
上述公式(5)~(7)中,rd表示对储层流体的无因次测试半径;td表示无因次渗流时间,为预知取值;表示地层流体粘度,为预知取值;CD表示地层总压缩系数,为预知取值;VPo表示无因次等时速率,为预知取值。
根据上述单介质渗流储层等时速率数学模型与定解条件可以导出单介质渗流储层的等时速率方程为公式(4),即为:
VP=λQBu【lg(((t2(tP+t1))/((t1(tP+t2)))/(t2-t1)】/Kh。
大量的试验证实,储层的等时速率与渗流特征的关系表现为:在以等时速率为分析前提条件下,储层的渗流过程,直接受制于储层的流体特性、生产压差和储层渗透率,是储层特性的直接反映。在测试开井过程中,对于单介质渗流储层而言,开井流动过程分为介质内渗流、介质向井筒渗流和井筒流动三个过程;关井后流动态为:关井初期介质向井筒渗流。测试开井过程中人为建立的生产压差,将形成不同渗流特征。当开井时间较长时,井筒流动压力与介质压力呈规律性降低,等时速率表现为规律性变化。
当开井时间过短,人工建立的生产压差没有得到平衡,井筒和近井带介质内压力与远井带介质内压力形成差异,关井初期的等时速率表现为非线性流规律性变化,但时间较短,进入续流后进入线性流规律性变化。当井筒内压力和近井带介质内压力与远井带介质内压力成规律性差异时,流动进入介质内渗流和介质向井筒内渗流的同步形态,代表着井筒储集过程的结束,储层进入相等的速率变化过程中,等时速率表现为线性流的规律性变化,其等时速率具有储层自身特征表现。
因此,对于单介质渗流储层而言,当开井过程中流体流动时间较长,流体的流动达到了线性流的状态,渗流特征符合等时速率方程,采用此等时速率方程进行储层分析研究是正确的。开井时间过短时,流动未达到线性流状态,渗流特征不符合等时速率方程,但井储之后便符合等时速率方程。
对于多介质渗流储层而言,开井流动过程分为孔隙渗流、裂隙渗流、裂缝窜流、裂缝向井筒窜流、井筒流动五个阶段,关井后流动态为:关井初期为孔隙(或称为:裂隙)向井筒和裂缝渗流,等时速率表现为非线性流规律性变化,当井筒与裂缝的压力达到孔隙(裂隙)或微裂缝压力时,进入裂隙渗流阶段,等时速率表现为裂隙线性流规律性变化,当裂隙的压力达到裂隙周围孔隙压力时,进入孔隙渗流状态,等时速率表现孔隙线性流规律性变化,而对于基质无储渗条件的裂缝性储层而言,则没有后期的孔隙性线性流阶段。无论是井筒储集阶段、续流阶段或径向流动阶段,等时速率取决于流动驱动力(也即:流动压差)的大小,等时速率总是与流动驱动力呈规律性的变化,储层渗流特性是储层真实流动能力的直接反映,压力速率是储层渗流能力的直接反映。本发明人在实现本发明的过程中发现,裂缝储层进入裂隙流动阶段后,渗流特征与孔隙的渗流特征虽然不近相同,但速率的变化与孔隙储层是一致的,呈现线性流动的特征,符合压力恢复等时速率方程,只是压力恢复受到裂缝非线性窜流影响,导致了在相同的时间内,实测的压力值远大于了线性流压力值,影响了压力恢复的时间,本发明实施例中将这一时间称为窜流时间,表示为tc。因此,裂缝性储层实测取得的压力值,将是孔渗模型压力时间t与窜流时间tC之和,由此推导裂缝性储层的等时速率方程则为:
VP=λQBu【lg((((t2+tc)(tP+t1+tc))/(((t1+tc)(tP+t2+tc)))】/(t2-t1)/Kh
(8)
P2-P1=λQBu【lg((((t2+tc)(tP+t1+tc))/(((t1+tc)(tP+t2+tc)))】/Kh(9)
根据储层损害分析试验获知,不论是何种污染,反映在储层渗流特征上则表现为有效渗透度的降低,它将使储层的导流和导压能力降低,导致我们在测压力恢复过程中,单位时间内实测的压力值则小于未受污染时的压力值,也就是说,受污染后要测取到储层未受污染时的压力值所需要的时间比原来的时间要长,这个时间差我们定义为tS,则测取的压力值是理论时间减去tS(tS的具体取值可以为小于或等于零的数值),由此得到,在储层受到污染后的等时速率方程为:
VP=λQBu【lg((((t2+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))】/(t2-t1)/Kh
(10)
P2-P1=λQBu【lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))】/Kh(11)
上述公式(8)与(10)即为根据等时速率数学模型设定多介质渗流储层的等时速率方程。其中,tc表示裂缝性储层导致的窜流时间,tS表示储层受污染后测取未受污染时的压力所需时间与未受污染时测取未受污染时的压力所需时间的差值。特别地,在公式(10)中,t1、t2表示关井点时刻,t2-t1表示关井总时间。其中,上述公式(10)中,tc和tS的取值为零时得到单介质渗流储层的等时速率方程。对于多介质渗流储层,通过本发明实施例获取到的窜流时间为体现污染特征的等效窜流时间,即:通过本发明实施例获取到的等效窜流时间tc=tc-tS。
在实际应用中,可以通过等时速率曲线的斜率自动计算污染时间Ts。在本发明实施例中,通过改变上述等时速率方程中的λ和tc,便可建立不同岩性、流体、渗流特征条件下的单介质渗流储层和多介质渗流储层的典型等时速率曲线模型,从而满足对续流数据和非达西流储层的评价分析需求。
本发明上述各公式中的VP即为模型等时速率VP’,t2与t1即为模型等时时间间隔t’的结束时刻t2’与开始时刻t1’。
图7为本发明方法试井资料中续流数据的分析处理方法另一个实施例的流程图。如图7所示,该实施例的试井资料中续流数据的分析处理方法包括:
步骤201,预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,该等时速率与各种影响因素之间的关系包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系。
步骤202,根据等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型。
步骤203,在对储层试井过程中,在井筒储层阶段完成后,在不同采集时刻分别采集井底压力,并获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线。
其中,P2与P1分别为采集时刻t2、t1采集到的压力值。如图8所示,为本发明实施例中VP~t曲线的一个具体示例。
步骤204,选取与实际测试等时速率曲线相符合的试井解释模型参数,根据预先建立的等时速率方程和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取VP’与模型等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线。
其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率。等时速率方程中的t1与t2分别表示等时时间间隔始点与未点,t1’与t2’分别模型表示等时时间间隔始点与未点,采用等时速率方程计算模型等时速率VP’时,t1’与t2’分别为等时速率方程中的t1与t2,计算出的VP即为VP’,此处区分t1’与t1、t2’与t2、VP’与VP的意义仅在于表示二者的具体取值未必相同,可以相同,也可以不同。
步骤205,绘制实测等时速率曲线中等时速率VP与对数lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))之间的第一曲线图,以及绘制采集的井底压力P与等时速率VP之间的第二曲线图。
如图9所示,为本发明实施例中第一曲线图的一个具体示例,其中的纵坐标表示等时速率VP,横坐标表示对数lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))。如图10所示,为本发明实施例中第二曲线图的一个具体示例,其中的纵坐标表示等时速率VP,横坐标表示井底压力P。图9与图10中,第一线段表示实测压力数据对应的关系线,第二线段表示曲线未端直线段的直线。
步骤206,分别利用曲线末端线性求解方法,对第一曲线与第二曲线进行求解,通过求取第一曲线未端的直线斜率获取实测等时速率曲线的近似渗透率K’,通过求取第二曲线未端的直线斜率获取近似地层压力Pi’,并根据实测等时速率曲线的曲率获取窜流系数λ,例如,可以依据VP=0.001858λQBu/Kh,在其它参数值已知的情况下获取窜流系数λ的取值。
步骤207,对近似渗透率K’、近似地层压力Pi’、与等效窜流时间tC-tS进行精细调整,直至模型等时速率曲线与实测等时速率曲线重合,以重合时模型等时速率曲线对应的渗透率K、地层压力Pi与tc-tS作为储层特征参数。
通过本发明上述实施例续流数据的分析处理方法获得储层特征参数后,可以进一步绘制表示不同采集时刻分别采集井底压力P与采集时刻t之间关系的实测P~t曲线,如图11所示,为本发明实施例中实测P~t曲线的一个具体示例。并根据获取的储层特征参数由公式(11)绘制模型P~t曲线,对比实测P~t曲线与模型P~t曲线是否完全拟合(即:重合),若不完全重合,可以进一步调整储层特征参数使实测P~t曲线与模型P~t曲线完全拟合,从而更精确的获得储层特征参数。
本发明实施例中,可以根据实测P~t曲线与模型P~t曲线是否完全拟合,来获取地层压力Pi;并且,可以用来分析判断模型等时速率曲线与实测等时速率曲线拟合时的模型是否一致,只有模型等时速率曲线与地层实际的类型一致,实测P~t曲线与模型P~t曲线的形态才完全一致,模型等时速率曲线与实测等时速率曲线拟合的结果才是正确的,否则便需进一步调整模型等时速率曲线的相关试井解释模型参数。
通过本发明上述实施例试井资料中续流数据的分析处理方法,获得储层特征参数后,根据本发明续流数据的分析处理方法的进一步实施例,可以基于这些储层特征参数进一步获取其它的储层特征参数。例如,利用关井总时间和开井生产时间tP,计算获取储层流体的测试半径rd;基于窜流时间tc计算异常点距离、裂缝半长;基于渗透率K和预知的基质渗透率K1计算获取窜流系数λ;利用渗透率K计算地层系数KH和流动系数KH/u。
图12为本发明试井资料中续流数据的分析处理系统一个实施例的结构示意图。该实施例的分析处理系统可用于实现本发明上述各分析处理方法实施例的流程。如图12所示,其包括采集单元301、参数选取单元302、存储单元303、计算单元304、获取单元305与拟合单元306。
其中,采集单元301,用于在对储层试井过程中,例如,在井筒储层阶段完成后,在不同采集时刻分别采集井底压力。
参数选取单元302,用于选取试井解释模型参数。
存储单元303,用于存储预先建立的等时速率数学模型。该等时速率数学模型基于预先测定的各种岩性条件下的等时速率与各种影响因素之间的关系获取得到的等时速率与储层岩性、储层结构、流体特性之间的关系数据建立。其中的等时速率与各种影响因素之间的关系包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系。
示例性地,该等时速率数学模型基于预先测定的各种岩性条件下的等时速率与各种影响因素之间的关系获取得到的等时速率与储层岩性、储层结构、流体特性之间的关系数据建立。其中,等时速率与各种影响因素之间的关系具体可以包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系。其中,等时速率与裂缝窜流特征的关系包括等时速率与有限导流裂缝干扰的关系、等时速率与无限导流裂缝干扰的关系。
计算单元304,用于根据参数选取单元302选取的试井解释模型参数和存储单元303中的等时速率数学模型,计算在不同模型时刻的模型等时速率VP’。其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率。
获取单元305,用于根据采集单元301在不同采集时刻分别采集井底压力,获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线。其中,P2与P1分别等时时间间隔t的结束时刻t2与开始时刻t1采集到的压力值。以及获取计算单元304计算出的模型等时速率VP’与模型等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线。
拟合单元306,用于将获取单元305获取的模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,该储层特征参数包括渗透率K与窜流时间tc。
本发明上述实施例提供的试井资料中续流数据的分析处理系统,以研究储层在相等时间内压力的变化速率为基础,基于等时速率解释理论与方法,在试井过程中在井筒储层阶段完成后,所取得的压力恢复或压力降落资料未达到稳定的径向流状态条件时,在井筒储层阶段完成后,在不同采集时刻分别采集井底压力,并获取表示等时速率VP与等时时间间隔t之间关系的实测等时速率曲线,并根据预先建立的等时速率数学模型和选取的试井解释模型参数,计算在不同模型时刻的模型等时速率VP’,并获取表示VP’与模型等时时间间隔t之间的对应关系的模型等时速率曲线,然后将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,包括渗透率K与窜流时间tc,从而实现了对续流数据的分析处理,可应用于全面分析测试半径范围内储层的渗流特征,实现对续流资料的解释处理,取得储层的特征参数。
根据本发明分析处理系统实施例的一个具体示例而非限制,存储单元303中存储的等时速率数学模型可以是:
Vp=λQBu【lg(((t2(tp+t1))/((t1(tp+t2)))】/(t2-t1)/kh(4)
其中,VP表示模型等时速率VP’;t1表示等时时间间隔t’的开始时刻t1’;t2表示等时时间间隔t’的结束时刻t2’;tP表示生产时间;Q表示平均生产产量;B表示地层流体体积系数;u表示粘度;K表示地层流体的渗透率;h表示储层厚度;λ表示地层流体的窜流系数,为预设大于零的常数。
其中,λ的取值可以根据储层的介质类型以及窜流时间tc预设为1.81~2.59。储层根据介质类型可以包括单介质渗流储层与多介质渗流储层。
根据本发明分析处理系统实施例的另一个具体示例而非限制,存储单元303中还可以包括根据上述公式(4)的等时速率数学模型设定的单介质渗流储层的等时速率方程,具体与公式(4)相同。
另外,存储单元303中还包括根据等时速率数学模型设定的多介质渗流储层的等时速率方程:
VP=λQBu【lg((((t2+tC-tS)(tP+t1+tC-tS))/(((t1+tC-tS)(tP+t2+tC-tS)))】/(t2-t1)/Kh(10)
其中,tc表示裂缝性储层导致的窜流时间,tS表示储层受污染后测取未受污染时的压力所需时间与未受污染时测取未受污染时的压力所需时间的差值。
根据本发明分析处理系统实施例的又一个具体示例而非限制,与本发明上述分析处理方法实施例相应地,计算单元304具体可以根据预先建立的等时速率方程和选取的试井解释模型参数,计算在不同模型时刻的模型等时速率。
图13为本发明试井资料中续流数据的分析处理系统另一个实施例的结构示意图。该实施例的分析处理系统可用于实现本发明图7所示实施例的分析处理方法流程。如图13所示,与图12所示的实施例相比,该实施例中,拟合单元306具体包括绘制子单元401、获取子单元402与调整子单元403。
其中,绘制子单元401,用于绘制获取单元305获取的实测等时速率曲线中等时速率VP与存储单元303存储的等时速率方程中的对数lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))之间的第一曲线图,以及绘制采集单元采集的井底压力P与等时速率VP之间的第二曲线图。
获取子单元402,用于分别利用曲线末端线性求解方法,对绘制子单元401绘制的第一曲线与第二曲线进行求解,获取实测等时速率曲线的近似渗透率K’和近似地层压力Pi’,并根据实测等时速率曲线的曲率获取窜流系数λ。
调整子单元403,用于对获取子单元402获取的近似渗透率K’、近似地层压力Pi’与窜流时间tc进行精细调整,并指示绘制子单元402重新绘制第一曲线图与第二曲线图,直至模型等时速率曲线与实测等时速率曲线重合,以重合时模型等时速率曲线对应的渗透率K、地层压力Pi与窜流时间tc作为储层特征参数。
根据本发明分析处理系统实施例的再一个具体示例而非限制,试井资料中续流数据的分析处理系统中除采集单元301以外的其它单元可操作性地耦合在用户设备中。其中的用户设备可以是计算机或者其它终端设备。
本说明书中各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似的部分相互参见即可。对于系统实施例而言,由于其与方法实施例基本对应,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
可能以许多方式来实现本发明的方法和系统。例如,可通过软件、硬件、固件或者软件、硬件、固件的任何组合来实现本发明的方法和系统。用于所述方法的步骤的上述顺序仅是为了进行说明,本发明的方法的步骤不限于以上具体描述的顺序,除非以其它方式特别说明。此外,在一些实施例中,还可将本发明实施为记录在记录介质中的程序,这些程序包括用于实现根据本发明的方法的机器可读指令。因而,本发明还覆盖存储用于执行根据本发明的方法的程序的记录介质。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
本发明实施例基于等时速率解释理论与方法,实现了对续流资料的解释处理,取得储层的特征参数,达到同径向流资料相同的分析处理之目的,以实现对非达西流储层的评价分析。
本发明实施例提出的等时速率解释理论与方法,是针对不同的储层类型进行分析处理的重要的方法,充分发挥其研究对象与传统解释方法的不同,应用等时速率与时间变化特征理论,以其解释理论基础为模型框架,以实际测试储层的速率特征为模型建立依据,形成了随实测储层特征而自动建立模型的强大功能,可以实现对各类储层的解释评价能力,应用于对砂岩、碳酸岩、火山岩、孔隙储层、裂缝储层,双孔储层等的解释评价。
本发明实施例提出的等时速率解释理论与方法,建立在渗流速率为储层特征的基础上,在测试过程中的各种渗流状态下,渗流速度的变化特征均代表储层的自身的渗流能力,从而实现了对续流资料的评价分析和径向流资料的评价分析。
本发明实施例提出的等时速率解释理论与方法,形成的针对续流资料的处理功能,用户可以在进行低渗储层测试时,无需考虑压力降对恢复时间的影响,最大限度的延长开井时间,扩大波及半径和测试范围,更加完善的了解低渗储层的非均质特征。同时,针对开发试井中关井时间过长影响产量的问题,用户可以依据测试目的,针对性的缩短关井时间,实现以测取续流资料为方法的短试井作业工艺,减少测试对油田产量的影响。
本发明的描述是为了示例和描述起见而给出的,而并不是无遗漏的或者将本发明限于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显然的。选择和描述实施例是为了更好说明本发明的原理和实际应用,并且使本领域的普通技术人员能够理解本发明从而设计适于特定用途的带有各种修改的各种实施例。
Claims (19)
1.一种试井资料中续流数据的分析处理方法,其特征在于,包括:
预先测定的各种岩性条件下等时速率与各种影响因素之间的关系,所述等时速率与各种影响因素之间的关系包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系;
根据等时速率与各种影响因素之间的关系,获取等时速率与储层岩性、储层结构、流体特性之间的关系数据,并建立等时速率数学模型;
在对储层试井过程中,在不同采集时刻分别采集井底压力,并获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线;其中,P2与P1分别为等时时间间隔t的结束时刻t2与开始时刻t1采集到的压力值;
选取试井解释模型参数,根据建立的等时速率数学模型和选取的试井解释模型参数,计算在不同时刻的模型等时速率VP’,并获取VP’与模型等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线;其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率;
将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,所述储层特征参数包括渗透率K与窜流时间tc。
2.根据权利要求1所述的方法,其特征在于,所述等时速率数学模型为:
VP=λQBu【lg(((t2(tP+t1))/((t1(tP+t2)))/(t2-t1)】/Kh;
其中,VP表示模型等时速率VP’;t1表示等时时间间隔t’的开始时刻t1’;t2表示等时时间间隔t’的结束时刻t2’;tP表示生产时间;Q表示平均生产产量;B表示地层流体体积系数;u表示粘度;K表示地层流体的渗透率;h表示储层厚度;λ表示地层流体的窜流系数,为预设大于零的常数。
3.根据权利要求2所述的方法,其特征在于,λ的取值根据储层的介质类型以及窜流时间tc预设为1.81~2.59。
4.根据权利要求2所述的方法,其特征在于,所述储层根据介质类型包括单介质渗流储层与多介质渗流储层。
5.根据权利要求4所述的方法,其特征在于,还包括:
根据所述等时速率数学模型设定单介质渗流储层的等时速率方程为:
VP=λQBu【lg(((t2(tP+t1))/((t1(tP+t2)))/(t2-t1)】/Kh。
6.根据权利要求4所述的方法,其特征在于,还包括:
根据所述等时速率数学模型设定多介质渗流储层的等时速率方程为:
VP=λQBu【lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))】/(t2-t1)/Kh;
其中,tc表示裂缝性储层导致的窜流时间,tS表示储层受污染后测取未受污染时的压力所需时间与未受污染时测取未受污染时的压力所需时间的差值。
7.根据权利要求5或6所述的方法,其特征在于,根据建立的等时速率数学模型和选取的试井解释模型参数,计算在不同模型时刻的模型等时速率包括:
根据预先建立的等时速率方程和选取的试井解释模型参数,计算在不同模型时刻的模型等时速率。
8.根据权利要求7所述的方法,其特征在于,将模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数包括:
绘制实测等时速率曲线中等时速率VP与对数lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))之间的第一曲线图,以及绘制采集的井底压力P与等时速率VP之间的第二曲线图;
分别利用曲线末端线性求解方法对第一曲线与第二曲线进行求解,获取实测等时速率曲线的近似渗透率K’和近似地层压力Pi’,并根据实测等时速率曲线的曲率获取窜流系数λ;
对近似渗透率K’、近似地层压力Pi’与窜流时间tc进行精细调整,直至模型等时速率曲线与实测等时速率曲线重合,以重合时模型等时速率曲线对应的渗透率K、地层压力Pi与窜流时间tc作为储层特征参数。
9.一种试井资料中续流数据的分析处理系统,其特征在于,包括:
采集单元,用于在对储层试井过程中,在不同采集时刻分别采集井底压力;
参数选取单元,用于选取试井解释模型参数;
存储单元,用于存储预先建立的等时速率数学模型;所述等时速率数学模型基于预先测定的各种岩性条件下的等时速率与各种影响因素之间的关系获取得到的等时速率与储层岩性、储层结构、流体特性之间的关系数据建立,所述等时速率与各种影响因素之间的关系包括等时速率与渗流状态的关系、等时速率与渗透率的关系、等时速率与污染的关系、等时速率与裂缝窜流特征的关系;
计算单元,用于根据参数选取单元选取的试井解释模型参数和存储单元中的等时速率数学模型,计算在不同时刻的模型等时速率VP’;其中,VP’为在等时速率数学模型下,根据选取的试井解释模型参数获取到的模型等时时间间隔t’对应的等时速率;
获取单元,用于根据采集单元在不同采集时刻分别采集井底压力,获取等时速率VP=(P2-P1)/(t2-t1)与等时时间间隔t=t2-t1之间的对应关系VP~t曲线作为实测等时速率曲线;其中,P2与P1分别为等时时间间隔t的结束时刻t2与开始时刻t1采集到的压力值;以及获取计算单元计算出的模型等时速率VP’与模型等时时间间隔t’=t2’-t1’之间的对应关系VP’~t’曲线作为模型等时速率曲线;
拟合单元,用于将获取单元获取的模型等时速率曲线与实测等时速率曲线进行拟合,获取储层特征参数,所述储层特征参数包括渗透率K与窜流时间tc。
10.根据权利要求9所述的系统,其特征在于,所述等时速率数学模型为:
VP=λQBu【lg(((t2(tP+t1))/((t1(tP+t2)))/(t2-t1)】/Kh;
其中,VP表示模型等时速率VP’;t1表示等时时间间隔t’的开始时刻t1’;t2表示等时时间间隔t’的结束时刻t2’;tP表示生产时间;Q表示平均生产产量;B表示地层流体体积系数;u表示粘度;K表示地层流体的渗透率;h表示储层厚度;λ表示地层流体的窜流系数,为预设大于零的常数。
11.根据权利要求10所述的系统,其特征在于,λ的取值根据储层的介质类型以及窜流时间tc预设为1.81~2.59。
12.根据权利要求10所述的系统,其特征在于,所述储层根据介质类型包括单介质渗流储层与多介质渗流储层。
13.根据权利要求12所述的系统,其特征在于,所述存储单元中还包括根据所述等时速率数学模型设定的单介质渗流储层的等时速率方程:
VP=λQBu【lg(((t2(tP+t1))/((t1(tP+t2)))/(t2-t1)】/Kh。
14.根据权利要求12所述的系统,其特征在于,所述存储单元中还包括根据所述等时速率数学模型设定的多介质渗流储层的等时速率方程:
VP=λQBu【lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))】/(t2-t1)/Kh;
其中,tc表示裂缝性储层导致的窜流时间,tS表示储层受污染后测取未受污染时的压力所需时间与未受污染时测取未受污染时的压力所需时间的差值。
15.根据权利要求13或14所述的系统,其特征在于,所述计算单元具体根据建立的等时速率方程和选取的试井解释模型参数,计算在不同模型时刻的模型等时速率。
16.根据权利要求15所述的系统,其特征在于,所述拟合单元包括:
绘制子单元,用于绘制获取单元获取的实测等时速率曲线中等时速率VP与存储单元存储的等时速率方程中的对数lg((((t2+tc-tS)(tP+t1+tc-tS))/(((t1+tc-tS)(tP+t2+tc-tS)))之间的第一曲线图,以及绘制采集单元采集的井底压力P与等时速率VP之间的第二曲线图;
获取子单元,用于分别利用曲线末端线性求解方法对第一曲线与第二曲线进行求解,获取实测等时速率曲线的近似渗透率K’和近似地层压力Pi’,并根据实测等时速率曲线的曲率获取窜流系数λ;
调整子单元,用于对获取子单元获取的近似渗透率K’、近似地层压力Pi’与窜流时间tc进行精细调整,并指示绘制子单元重新绘制第一曲线图与第二曲线图,直至模型等时速率曲线与实测等时速率曲线重合,以重合时模型等时速率曲线对应的渗透率K、地层压力Pi与窜流时间tc作为储层特征参数。
17.根据权利要求9至14任意一项所述的系统,其特征在于,所述试井资料中续流数据的分析处理系统中除采集单元以外的其它单元可操作性地耦合在用户设备中。
18.根据权利要求15所述的系统,其特征在于,所述试井资料中续流数据的分析处理系统中除采集单元以外的其它单元可操作性地耦合在用户设备中。
19.根据权利要求16所述的系统,其特征在于,所述试井资料中续流数据的分析处理系统中除采集单元以外的其它单元可操作性地耦合在用户设备中。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210281796.3A CN103452547B (zh) | 2012-08-09 | 2012-08-09 | 试井资料中续流数据的分析处理方法和系统 |
PCT/CN2013/080791 WO2014023201A1 (zh) | 2012-08-09 | 2013-08-05 | 试井资料中续流数据的分析处理方法和系统 |
US14/536,528 US20150066372A1 (en) | 2012-08-09 | 2014-11-07 | Method and system for analyzing and processing continued flow data in well testing data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210281796.3A CN103452547B (zh) | 2012-08-09 | 2012-08-09 | 试井资料中续流数据的分析处理方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103452547A CN103452547A (zh) | 2013-12-18 |
CN103452547B true CN103452547B (zh) | 2016-03-23 |
Family
ID=49735330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210281796.3A Active CN103452547B (zh) | 2012-08-09 | 2012-08-09 | 试井资料中续流数据的分析处理方法和系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20150066372A1 (zh) |
CN (1) | CN103452547B (zh) |
WO (1) | WO2014023201A1 (zh) |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105756660B (zh) * | 2014-12-19 | 2018-11-16 | 中石化胜利石油工程有限公司钻井工艺研究院 | 一种气井压回法压井时机的确定方法 |
CN105930932B (zh) * | 2016-04-25 | 2019-11-12 | 中石化石油工程技术服务有限公司 | 基于含气指数的页岩气层标准化无阻流量的获取方法 |
CN106021793A (zh) * | 2016-06-01 | 2016-10-12 | 中国地质大学(武汉) | 基于存储系数和渗流系数的低渗储层甜点评价方法 |
CN106194163B (zh) * | 2016-07-26 | 2019-04-12 | 中国石油集团渤海钻探工程有限公司 | 一种油水井试井资料解释自动选择方法 |
CN108490487B (zh) * | 2018-02-24 | 2019-10-18 | 中国石油化工股份有限公司 | 利用试井资料确定近井裂缝参数的新方法 |
CN110593860B (zh) * | 2018-05-25 | 2021-08-13 | 中国科学技术大学 | 一种参数计算方法、装置及设备 |
CN108729906A (zh) * | 2018-06-07 | 2018-11-02 | 西南石油大学 | 一种低渗致密气藏改进逆序修正等时试井测试方法 |
CN111255434B (zh) * | 2018-11-30 | 2023-05-26 | 中国石油天然气股份有限公司 | 气井的试井方法、装置和计算机存储介质 |
CN111594113B (zh) * | 2019-02-20 | 2022-06-17 | 中国石油化工股份有限公司 | 一种致密储层井间裂缝开度动态反演方法 |
CN111691869B (zh) * | 2019-03-15 | 2023-07-04 | 中国石油化工股份有限公司 | 一种高含水致密气藏气井产量变化趋势的预测方法 |
CN112343566B (zh) * | 2019-08-08 | 2022-12-02 | 中国石油天然气股份有限公司 | 注入聚合物井的井底流压确定方法及装置 |
CN110374572B (zh) * | 2019-08-12 | 2021-01-05 | 西南石油大学 | 非常规油气藏多段压裂水平井各改造段特性参数确定方法 |
CN110749532A (zh) * | 2019-10-30 | 2020-02-04 | 高伟 | 水利工程渗流智能监测系统及方法 |
RU2722900C1 (ru) * | 2019-12-23 | 2020-06-04 | Общество с ограниченной ответственностью "Газпром добыча Уренгой" | Способ прогнозирования длительности регистрации кривой восстановления давления скважины |
CN111537663B (zh) * | 2020-04-20 | 2022-10-04 | 中国石油天然气集团有限公司 | 一种岩性识别剂携带装置及基于其的岩性识别系统和方法 |
CN113931621B (zh) * | 2020-07-14 | 2023-08-22 | 中国石油天然气股份有限公司 | 气井积液信息的确定方法、装置及存储介质 |
CN113947005B (zh) * | 2020-07-15 | 2022-09-23 | 中国石油化工股份有限公司 | 一种基于机器学习的试井解释方法及系统 |
WO2022051626A1 (en) * | 2020-09-04 | 2022-03-10 | Saudi Arabian Oil Company | Graph framework (database methods) to analyze trillion cell reservoir and basin simulation results |
CN112211626B (zh) * | 2020-10-30 | 2022-03-11 | 西南石油大学 | 一种非均质气藏气井产能试井测试类型的优选方法 |
CN112392475B (zh) * | 2020-11-19 | 2023-07-18 | 中海石油(中国)有限公司 | 一种用于确定近临界油藏测试米采油指数的方法 |
CN113050190A (zh) * | 2021-03-03 | 2021-06-29 | 河北益坤岩土工程新技术有限公司 | 直线边界非稳定流抽水试验水文地质参数智能计算方法 |
CN113065246B (zh) * | 2021-03-26 | 2024-03-29 | 中国科学技术大学 | 一种基于射孔测试联作系统的地层参数反演方法 |
CN113027391B (zh) * | 2021-04-22 | 2022-10-28 | 中海油田服务股份有限公司 | 一种增注液波及半径的计算方法 |
CN113552037B (zh) * | 2021-06-03 | 2022-08-02 | 浙江大学 | 一种测试垃圾双孔隙度渗流参数的装置及方法 |
CN113586046B (zh) * | 2021-06-23 | 2023-07-28 | 中海油能源发展股份有限公司 | 一种跟踪评价水平井调驱效果的方法 |
CN115822562B (zh) * | 2022-12-28 | 2023-07-11 | 中海石油(中国)有限公司海南分公司 | 一种考虑层内窜流的纵向非均质气藏产能评价方法 |
CN115828787B (zh) * | 2022-12-13 | 2023-09-08 | 西南石油大学 | 一种改进的mdh试井分析方法 |
CN116006167B (zh) * | 2022-12-19 | 2024-05-14 | 新奥(天津)能源技术有限公司 | 压裂缝网检测方法、装置、介质及设备 |
CN116255116A (zh) * | 2023-02-09 | 2023-06-13 | 新疆敦华绿碳技术股份有限公司 | Co2驱底水砂岩油藏的油井的防腐措施确定方法 |
CN116341393B (zh) * | 2023-05-26 | 2023-08-04 | 中国石油大学(华东) | 一种自动非稳态试井解释方法、装置、设备及介质 |
CN116399784B (zh) * | 2023-06-07 | 2023-08-18 | 太原理工大学 | 一种模拟裂缝型热储层地热开采渗流试验的装置及方法 |
CN116455946B (zh) * | 2023-06-19 | 2023-09-19 | 安徽井上天华科技有限公司 | 一种基于云端高频井口压力生产数据分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4858130A (en) * | 1987-08-10 | 1989-08-15 | The Board Of Trustees Of The Leland Stanford Junior University | Estimation of hydraulic fracture geometry from pumping pressure measurements |
CN1438404A (zh) * | 2003-02-17 | 2003-08-27 | 大庆油田有限责任公司 | 不依赖径向流直线段的试井分析方法 |
CN101560879A (zh) * | 2008-04-15 | 2009-10-21 | 中国石油大学(北京) | 用于低渗透气藏的试井分析控制系统及方法 |
-
2012
- 2012-08-09 CN CN201210281796.3A patent/CN103452547B/zh active Active
-
2013
- 2013-08-05 WO PCT/CN2013/080791 patent/WO2014023201A1/zh active Application Filing
-
2014
- 2014-11-07 US US14/536,528 patent/US20150066372A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4858130A (en) * | 1987-08-10 | 1989-08-15 | The Board Of Trustees Of The Leland Stanford Junior University | Estimation of hydraulic fracture geometry from pumping pressure measurements |
CN1438404A (zh) * | 2003-02-17 | 2003-08-27 | 大庆油田有限责任公司 | 不依赖径向流直线段的试井分析方法 |
CN101560879A (zh) * | 2008-04-15 | 2009-10-21 | 中国石油大学(北京) | 用于低渗透气藏的试井分析控制系统及方法 |
Non-Patent Citations (1)
Title |
---|
低渗井的早期续流段试井资料解释方法;李笑萍等;《大庆石油地质与开发》;20051031;第24卷(第5期);第55-58页 * |
Also Published As
Publication number | Publication date |
---|---|
US20150066372A1 (en) | 2015-03-05 |
WO2014023201A1 (zh) | 2014-02-13 |
CN103452547A (zh) | 2013-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103452547B (zh) | 试井资料中续流数据的分析处理方法和系统 | |
CN102953726B (zh) | 一种水驱油田优势通道识别方法及装置 | |
CN105510203B (zh) | 一种确定不同温度梯度下砂岩油藏油水相对渗透率的方法 | |
CN104912551B (zh) | 一种油水相渗曲线和驱油效率的标定方法 | |
CN108661631B (zh) | 一种产量预测方法 | |
CN104750896A (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN106932324B (zh) | 一种确定高含水砂岩油藏储层渗透率变化规律的方法 | |
CN104564041A (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
CN110096718B (zh) | 一种获取碳酸盐岩油藏中溶洞的体积的方法 | |
CN105740563A (zh) | 一种成熟油田二次开发之优势通道识别技术 | |
CN104314547A (zh) | Hafd水平井压裂酸化分段优化方法及系统 | |
CN104634804A (zh) | 一种利用核磁共振t2谱确定储层相对渗透率的方法 | |
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
CN110107277B (zh) | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 | |
CN106204304A (zh) | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 | |
CN104100263A (zh) | 一种非均匀水淹层剩余油饱和度的确定方法 | |
CN104695952B (zh) | 非常规气藏地层参数的解释方法及解释系统 | |
CN111502652A (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
CN105204081B (zh) | 一种预测页岩气储层粘土矿物及其组分含量的方法 | |
CN108982320A (zh) | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 | |
CN114201932A (zh) | 一种复杂情况下的致密油藏压裂井试井模拟方法 | |
CN106501146B (zh) | 一种致密油藏物性上限确定方法及装置 | |
CN108843296B (zh) | 一种基于多因素影响下的单井重复压裂效果预测方法 | |
CN105844011A (zh) | 一种基于毛管模型的渗透率计算方法 | |
Kolin et al. | Pressure build-up test analysis of the reservoir system with the multiphase flow |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |