CN110598167B - 低渗透油藏油水相对渗透率实验数据的处理方法 - Google Patents
低渗透油藏油水相对渗透率实验数据的处理方法 Download PDFInfo
- Publication number
- CN110598167B CN110598167B CN201910970956.7A CN201910970956A CN110598167B CN 110598167 B CN110598167 B CN 110598167B CN 201910970956 A CN201910970956 A CN 201910970956A CN 110598167 B CN110598167 B CN 110598167B
- Authority
- CN
- China
- Prior art keywords
- oil
- water
- relative permeability
- permeability
- core
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 159
- 230000035699 permeability Effects 0.000 title claims abstract description 149
- 238000003672 processing method Methods 0.000 title abstract description 4
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000002347 injection Methods 0.000 claims abstract description 21
- 239000007924 injection Substances 0.000 claims abstract description 21
- 239000012530 fluid Substances 0.000 claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000002474 experimental method Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000000605 extraction Methods 0.000 claims abstract description 9
- 239000011435 rock Substances 0.000 claims description 49
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 239000007788 liquid Substances 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 3
- 239000008346 aqueous phase Substances 0.000 abstract 1
- 239000012071 phase Substances 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 229910052582 BN Inorganic materials 0.000 description 1
- PZNSFCLAULLKQX-UHFFFAOYSA-N Boron nitride Chemical compound N#B PZNSFCLAULLKQX-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Data Mining & Analysis (AREA)
- Immunology (AREA)
- Animal Husbandry (AREA)
- Computational Mathematics (AREA)
- Pathology (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Biochemistry (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Agronomy & Crop Science (AREA)
- Algebra (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Dispersion Chemistry (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供一种低渗透油藏油水相对渗透率实验数据的处理方法,包括:开展低渗透岩心最小启动压力梯度实验,获取地层油通过该低渗透岩心的最小启动压力值;开展非稳态条件油水相对渗透率实验;计算相邻两个时刻注入采出能力比,该值记作为Z值;根据JBN方法,计算含水饱和度及油、水相对渗透率曲线;计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线;计算不同含水饱和度下的水相相对渗透率曲线。该低渗透油藏油水相对渗透率实验数据的处理方法综合考虑了低渗透油藏流体流动的非线性渗流特征和启动压力梯度存在的问题,为明确低渗透油藏渗流规律及完善低渗透油藏非线性渗流数值模拟技术奠定基础。
Description
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种低渗透油藏油水相对渗透率实验数据的处理方法。
背景技术
研究结果表明,由于低渗透多孔介质孔隙几何形状的复杂性、流体性质的特殊性和同时受到流动条件的制约,使得流体在低渗透多孔介质中渗流特征不再符合传统的达西渗流规律,表现为明显的非线性渗流特征,同时具有一定的启动压力梯度。
油水两相相对渗透率曲线被广泛应用于油田开发动态分析、开发方案优化及油藏数值模拟技术中,其综合反映了储层中油水两相的渗流规律和基本特征。目前,油水相对渗透率的主要获取方法:(1)理论模型;(2)岩心实验测试。而岩心实验测试法较真实的反映流体在油藏中的流动特征而被研究人员广泛的认可。目前,实验测试油水相对渗透率的主要方法主要是稳态法和非稳态法。使用稳态法测试油水相对渗透率可以直接获取可靠的含水饱和度。然而,稳态法不适合油田开发的实际应用,测试油水相对渗透率需要更多的时间和更高的成本。测试油水相对渗透率的非稳态法是基于水驱油原理,描述了驱替过程中含水饱和度分布与时间和距离的关系。非稳态法可以分为等速法和恒压法。利用等速法测量低渗透岩心的油水相对渗透率时,由于岩心的孔隙度和渗透率低,导致驱替压力传播缓慢。较高的流速可以导致在注入端有较高的压力而出口端没有流体输出,而较低的流速可以导致较低的驱替压力差和较长的测试时间。
随着低渗透率油藏开发的发展,非稳态法的应用变得尤为重要,被研究人员广泛的采用。目前,使用非稳态法测试油水相对渗透率时,应用JBN方法处理实验数据。而JBN方法的主要理论基础是达西渗流规律和BL水驱油方程。当流体的渗流规律不再遵循达西渗流时,应用JBN方法计算低渗透油藏的油水相对渗透率必然产生较大的误差。
在申请号:201410221806.3的中国专利申请中公开了一种利用煤层气井生产数据测量气水相渗曲线的方法,通过物质平衡方法得到平均储层压力与平均含水饱和度,进而得到绝对渗透率与相对渗透率变化;利用束缚水饱和度和孔隙度的关系,预测束缚水饱和度的变化;将绝对渗透率和端点标定后的相对渗透率在同一状态下耦合,得到区域平均有效渗透率动态预测模型,进而得到有效渗透率曲线。该专利是针对煤层气藏中气水渗流规律建立的气水相对渗透率的处理方法,不适用于油藏。
在申请号:201810257741.6的中国专利申请中提供一种基于Harris模型的油相相对渗透率曲线表征方法,包括:选取能反映区块或油藏物性特征的典型油水相对渗透率曲线数据;对油相相对渗透率和含水饱和度进行标准化处理,判断标准化后的油相相对渗透率与含水饱和度之间的关系;基于Harris模型构建油相相对渗透率与含水饱和度关系曲线的函数表达式;应用最小二乘法求解构建油相相对渗透率与含水饱和度关系曲线的函数表达式中的系数。该专利是在已有油相相对渗透率曲线数据基础上,基于Harris模型构建油相相对渗透率与含水饱和度理论模型。该专利中仅建立了油相相对渗透率理论模型,未建立水相相对渗透率模型。
为此我们发明了一种新的低渗透油藏油水相对渗透率实验数据的处理方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种综合考虑了低渗透油藏流体流动的非线性渗流特征和启动压力梯度存在的问题的低渗透油藏油水相对渗透率实验数据的处理方法。
本发明的目的可通过如下技术措施来实现:低渗透油藏油水相对渗透率实验数据的处理方法,该低渗透油藏油水相对渗透率实验数据的处理方法包括:步骤1,开展低渗透岩心最小启动压力梯度实验,获取地层油通过该低渗透岩心的最小启动压力值;步骤2,开展非稳态条件油水相对渗透率实验,记录不同时刻、不同驱替压力下的液量、油量、水量及累积液量、油量和水量这些实验数据信息;步骤3,计算相邻两个时刻注入采出能力比,该值记作为Z值;步骤4,根据JBN方法,计算含水饱和度及油、水相对渗透率曲线;步骤5,计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线;步骤6,计算不同含水饱和度下的水相相对渗透率曲线。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,测量岩心长度、横截面积孔隙度和渗透率这些常规岩心的基本参数及流体粘度和密度这些流体的常规参数,开展低渗透岩心最小启动压力梯度实验,获取地层油通过该低渗透岩心的最小启动压力值。
在步骤3中,计算相邻两个时刻注入采出能力比Z值的计算公式为:
其中:式中,Q(t)为t时刻水的累积注入量,cm3;φ为低渗岩心的孔隙度,无因次;A为岩心的横截面积,cm2;L岩心的长度,cm;Wi为岩心的累积注入倍数,无因次;v为记录的岩心出口端的液体流出速度,cm/s;μo为地层油(或模拟油)的粘度,mPa·s;k为岩心的绝对渗透率,μm2;加载在岩心两端的驱替压力梯度,0.1MPa/cm;Z为注入采出能力比,无因次。
在步骤4中,根据JBN方法,计算含水饱和度swe,油相对渗透率曲线krod(swe)和水相对渗透率曲线krwd(swe),其计算公式为:
式中,Z为注入采出能力比,无因次;swe为岩心出口端的含水饱和度;fo(swe)为岩心出口端含水饱和度为swe时的含油率,无因次;krwd(swe)为应用JBN方法计算的岩心出口含水饱和度为swe时的水相相对渗透率,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率,无因次。
在步骤5中,根据步骤1中获得的最小启动压力梯度Go、步骤3中的Z值及步骤4中得到的油相相对渗透率值krod(swe),计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线kro(swe)。
在步骤5中,计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线kro(swe)的计算公式为
式中,Go为地层油通过低渗岩心的最小启动压力梯度,0.1MPa/cm;Z为注入采出能力比,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率,无因次;kro(swe)为考虑非线性渗流特征后且岩心出口端含水饱和度为swe时的油相相对渗透率,无因次;krw(swe)岩心出口端含水饱和度为swe时的水相相对渗透率,无因次,加载在岩心两端的驱替压力梯度,0.1MPa/cm;I为单位量纲。
在步骤6中,根据步骤4中得到的水相、油相相对渗透率值及步骤5中得到考虑非线性渗流特征的油相相对渗透率值,计算不同含水饱和度下的水相相对渗透率曲线。
在步骤6中,计算不同含水饱和度下的水相相对渗透率曲线的计算公式为:
式中,Go为地层油通过低渗岩心的最小启动压力梯度,0.1MPa/cm;加载在岩心两端的驱替压力梯度,0.1MPa/cm;krwd(swe)为应用JBN方法计算的岩心出口含水饱和度为swe时的水相相对渗透率,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率,无因次;kro(swe)为考虑非线性渗流特征后且岩心出口端含水饱和度为swe时的油相相对渗透率,无因次;krw(swe)岩心出口端含水饱和度为swe时的水相相对渗透率,无因次。
流体在低渗透多孔介质中渗流特征不再符合传统的达西渗流规律,表现为明显的非线性渗流特征,同时具有一定的启动压力梯度。目前,实验室内广泛使用的是非稳态的JBN方法处理油水相对渗透率实验数据。而JBN方法的主要理论基础是达西渗流规律和BL水驱油方程。当低渗透储层中流体的规律不再遵循传统的达西渗流规律情况下,应用JBN方法计算低渗透油藏的油水相对渗透率必然产生较大的误差。
本发明中的低渗透油藏油水相对渗透率实验数据的处理方法,综合考虑了低渗透油藏流体流动的非线性渗流特征和存在启动压力梯度的问题,而目前实验室内广泛使用的JBN方法均没有考虑启动压力梯度及非线性渗流特征。本发明中提出的低渗透油藏油水相对渗透率处理方法,解决了现有JBN方法在处理低渗透油藏水油相对渗透率实验数据处理方面存在的技术难题,为明确低渗透油藏渗流规律及完善低渗透油藏非线性渗流数值模拟技术奠定基础。
附图说明
图1为本发明的低渗透油藏油水相对渗透率实验数据的处理方法的一具体实施例的流程图;
图2为本发明的一具体实施例中基于JBN方法和考虑非线性渗流特征的油水相对渗透率曲线对比图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
本发明的低渗透油藏油水相对渗透率实验数据的处理方法包括了以下步骤:
在步骤1中,测量岩心长度、横截面积孔隙度和渗透率等常规岩心的基本参数及流体粘度和密度等流体的常规参数,开展低渗透岩心最小启动压力梯度实验,获取地层油(或模拟油)通过该低渗透岩心的最小启动压力值Go。
在步骤2中,开展非稳态条件油水相对渗透率实验,记录不同时刻、不同驱替压力下的液量、油量、水量及累积液量、油量和水量等实验数据信息。
在步骤3中,计算相邻两个时刻注入采出能力比,该值记作为Z值,其计算公式为
在步骤4中,根据JBN方法,计算含水饱和度swe,油相对渗透率曲线krod(swe)和水相对渗透率曲线krwd(swe),其计算公式为
在步骤5中,根据步骤1中获得的最小启动压力梯度Go、步骤3中的Z值及步骤4中得到的油相相对渗透率值krod(swe),计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线kro(swe),其计算公式为
步骤6,根据步骤4中得到的水相、油相相对渗透率值及步骤5中得到考虑非线性渗流特征的油相相对渗透率值,计算不同含水饱和度下的水相相对渗透率曲线,其计算公式为
式中,Go为地层油通过低渗岩心的最小启动压力梯度,0.1MPa/cm;Q(t)为t时刻水的累积注入量,cm3;φ为低渗岩心的孔隙度,无因次;A为岩心的横截面积,cm2;L岩心的长度,cm;Wi为岩心的累积注入倍数,无因次;v为记录的岩心出口端的液体流出速度,cm/s;μo为地层油(或模拟油)的粘度,mPa·s;k为岩心的绝对渗透率,μm2;加载在岩心两端的驱替压力梯度,0.1MPa/cm;Z为注入采出能力比,无因次;swe为岩心出口端的含水饱和度;fo(swe)为岩心出口端含水饱和度为swe时的含油率,无因次;krwd(swe)为应用JBN方法计算的岩心出口含水饱和度为swe时的水相相对渗透率,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率,无因次;kro(swe)为考虑非线性渗流特征后且岩心出口端含水饱和度为swe时的油相相对渗透率,无因次;krw(swe)岩心出口端含水饱和度为swe时的水相相对渗透率,无因次。
以下为应用本发明的一具体实施例,如图1所示,图1为本发明的低渗透油藏油水相对渗透率实验数据的处理方法的流程图。
在步骤101中,测量岩心长度、横截面积孔隙度和渗透率等常规岩心的基本参数及流体粘度和密度等流体的常规参数;开展地层油通过低渗透岩心的最小启动压力梯度实验,获取地层油(或模拟油)通过该低渗透岩心的最小启动压力值,实验用岩样基本参数见表1所示,实验用模拟油及注入水参数见表2所示。
表1岩样基本参数表
表2模拟油及注入水基本参数表
在步骤102中,开展非稳态条件油水相对渗透率实验,记录不同时刻、不同驱替压力下的瞬时指标(液量、油量、含水)及累积指标(液量、油量)等实验数据信息,计算结果见表3。
表3非稳态实验数据表
表4相邻两个时刻注入采出能力比Z
其中:
式中:角标t1和角标t2为两个相邻的时刻。
在步骤104中,根据JBN方法,计算岩心出口端不同时刻的含水饱和度swe,油相对渗透率曲线krod(swe)和水相对渗透率曲线krwd(swe)。计算结果见边5。
岩心出口端含水饱和度的计算方法为
swe=swi+Vo-fo(swe)·Wi (2)
岩心出口端油相对渗透率曲线krod(swe)的计算公式为
岩心出口端水相相对渗透率曲线krwd(swe)的计算公式为
表5 JBN方法计算油水相对渗透率表
swe | krod(JBN) | krwd(JBN) |
0.505 | 1 | 0 |
0.6053 | 0.1697 | 0.0182 |
0.6488 | 0.0194 | 0.1008 |
0.6655 | 0.0030 | 0.0947 |
0.6759 | 0.0039 | 0.1346 |
0.6992 | 0.0014 | 0.1489 |
0.7087 | 0.0007 | 0.1019 |
0.7246 | 0.0003 | 0.1661 |
在步骤105中,根据步骤101中获得的样品油通过低渗岩心的最小启动压力梯度Go、步骤3中的Z值及步骤104中得到的油相相对渗透率值krod(swe),计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线kro(swe),其计算公式为
步骤106,根据步骤101中获得的样品油通过低渗岩心的最小启动压力梯度Go、步骤104中基于JBN方法得到的水相相对渗透率值krwd(swe)、油相相对渗透率值krod(swe)及步骤105中得到考虑非线性渗流特征的油相相对渗透率值kro(swe),计算不同含水饱和度下的水相相对渗透率值krw(swe),其计算公式为
根据步骤105和步骤106计算的考虑非线性渗流特征的油水相对渗透率见表6所示。基于JBN方法和考虑非线性渗流特征的油水相对渗透率曲线对比如图2所示。
表6考虑非线性渗流特征的油水相对渗透率表
sw | kro | krw |
0.505 | 1 | 0 |
0.6053 | 0.1605 | 0.0171 |
0.6488 | 0.0183 | 0.0943 |
0.6655 | 0.0028 | 0.0888 |
0.6759 | 0.0037 | 0.1269 |
0.6992 | 0.0014 | 0.1403 |
0.7087 | 0.0007 | 0.0962 |
0.7246 | 0.0003 | 0.1565 |
式中,Go为地层油通过低渗岩心的最小启动压力梯度,0.1MPa/cm;Q(t)为t时刻水的累积注入量,cm3;φ为低渗岩心的孔隙度,无因次;A为岩心的横截面积,cm2;L岩心的长度,cm;Wi为岩心的无因次累积注入倍数,无因次;v为记录的岩心出口端的液体流出速度,cm/s;μo为地层油(或模拟油)的粘度,mPa·s;μw为模拟水的粘度,mPa·s;k为岩心的绝对渗透率,μm2;加载在岩心两端的驱替压力梯度,0.1MPa/cm;Z为注入采出能力比,无因次;swe为岩心出口端的含水饱和度;fo(swe)为岩心出口端含水饱和度为swe时的含油率,无因次;krwd(swe)为应用JBN方法计算的岩心出口含水饱和度为swe时的水相相对渗透率,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率,无因次;kro(swe)为考虑非线性渗流特征的且岩心出口端含水饱和度为swe时的油相相对渗透率,无因次;krw(swe)岩心出口端含水饱和度为swe时的水相相对渗透率,无因次;swi为岩心束缚水饱和度;Vo为无因次累积采油量,无因次;Qo(t)为t时刻的岩心出口端累积油量,cm3。
Claims (2)
1.低渗透油藏油水相对渗透率实验数据的处理方法,其特征在于,该低渗透油藏油水相对渗透率实验数据的处理方法包括:
步骤1,开展低渗透岩心最小启动压力梯度实验,获取地层油通过该低渗透岩心的最小启动压力值;
步骤2,开展非稳态条件油水相对渗透率实验,记录不同时刻、不同驱替压力下的液量、油量、水量及累积液量、油量和水量这些实验数据信息;
步骤3,计算相邻两个时刻注入采出能力比,记作为Z值;
步骤4,根据JBN方法,计算含水饱和度及油、水相对渗透率曲线;
步骤5,计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线;
步骤6,计算不同含水饱和度下的水相相对渗透率曲线;
在步骤3中,计算相邻两个时刻注入采出能力比Z值的计算公式为:
式中,Q(t)为t时刻水的累积注入量,cm3;φ为低渗岩心的孔隙度,无因次;A为岩心的横截面积,cm2;L岩心的长度,cm;Wi为岩心的累积注入倍数,无因次;v为记录的岩心出口端的液体流出速度,cm/s;μo为地层油或模拟油的粘度,mPa·s;k为岩心的绝对渗透率,μm2;加载在岩心两端的驱替压力梯度,0.1MPa/cm;Z为注入采出能力比,无因次;
式中:角标t1和角标t2为两个相邻的时刻;
在步骤4中,根据JBN方法,计算含水饱和度swe,油相对渗透率曲线krod(swe)和水相对渗透率曲线krwd(swe),油相对渗透率曲线krod(swe)的计算公式为
式中,Z为注入采出能力比,无因次;swe为岩心出口端的含水饱和度;fo(swe)为岩心出口端含水饱和度为swe时的含油率,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率曲线,无因次;
水相对渗透率曲线krwd(swe)的计算公式为
swe=swi+Vo-fo(swe)·Wi (5)
swi为岩心束缚水饱和度;krwd(swe)为应用JBN方法计算的岩心出口含水饱和度为swe时的水相相对渗透率,无因次;
在步骤5中,根据步骤1中获得的最小启动压力梯度Go、步骤3中的Z值及步骤4中得到的油相相对渗透率曲线krod(swe),计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线kro(swe);
计算不同含水饱和度下的考虑非线性渗流特征的油相相对渗透率曲线kro(swe)的计算公式为
式中,Go为地层油通过低渗岩心的最小启动压力梯度,0.1MPa/cm;Z为注入采出能力比,无因次;krod(swe)为应用JBN方法计算的含水饱和度为swe时的油相相对渗透率,无因次;加载在岩心两端的驱替压力梯度,0.1MPa/cm;I为单位量纲;
在步骤6中,根据步骤4中得到的水相相对渗透率曲线、油相相对渗透率曲线及步骤5中得到考虑非线性渗流特征的油相相对渗透率曲线,计算不同含水饱和度下的水相相对渗透率曲线;
计算不同含水饱和度下的水相相对渗透率曲线的计算公式为:
2.根据权利要求1所述的低渗透油藏油水相对渗透率实验数据的处理方法,其特征在于,在步骤1中,测量岩心长度、横截面积孔隙度和渗透率这些常规岩心的基本参数及流体粘度和密度这些流体的常规参数,开展低渗透岩心最小启动压力梯度实验,获取地层油通过该低渗透岩心的最小启动压力值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910970956.7A CN110598167B (zh) | 2019-10-11 | 2019-10-11 | 低渗透油藏油水相对渗透率实验数据的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910970956.7A CN110598167B (zh) | 2019-10-11 | 2019-10-11 | 低渗透油藏油水相对渗透率实验数据的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110598167A CN110598167A (zh) | 2019-12-20 |
CN110598167B true CN110598167B (zh) | 2022-10-28 |
Family
ID=68866841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910970956.7A Active CN110598167B (zh) | 2019-10-11 | 2019-10-11 | 低渗透油藏油水相对渗透率实验数据的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110598167B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111208048A (zh) * | 2020-01-17 | 2020-05-29 | 中国石油天然气股份有限公司 | 一种基于相渗测试的贾敏效应动态变化定量表征方法 |
CN111220509A (zh) * | 2020-01-20 | 2020-06-02 | 中国石油天然气股份有限公司 | 一种考虑渗透率时变性的油水相对渗透率曲线修正方法 |
CN111535787B (zh) * | 2020-04-09 | 2022-02-22 | 中国石油化工股份有限公司 | 高含水油藏动态渗流界面的判识模型及判识界限构建方法 |
CN112163379B (zh) * | 2020-09-21 | 2022-02-15 | 西南石油大学 | 一种基于孔隙网络模型的非稳态两相相对渗透率计算方法 |
CN112081588B (zh) * | 2020-10-26 | 2022-02-15 | 西南石油大学 | 一种判断相对渗透率曲线中化学体系黏度变化的方法 |
CN112081589B (zh) * | 2020-10-26 | 2022-03-01 | 西南石油大学 | 标准化等渗点变量获取化学体系相对渗透率曲线的方法 |
CN112685970B (zh) * | 2020-12-25 | 2022-08-30 | 中国石油大学(华东) | 一种水驱油藏流动单元渗流界面定量表征方法及系统 |
CN113295594B (zh) * | 2021-05-20 | 2022-07-12 | 中国石油大学(北京) | 一种基于分流模型与ct扫描的相对渗透率测定系统及方法 |
CN113295592B (zh) * | 2021-05-20 | 2022-06-07 | 中国石油大学(北京) | 一种低渗透岩石的相对渗透率测定系统及方法 |
CN114330155B (zh) * | 2021-12-06 | 2024-05-24 | 西安石油大学 | 考虑离子水浓度下岩心自发渗吸油水相对渗透率的预测方法 |
CN115234216B (zh) * | 2022-05-25 | 2023-09-19 | 东营市永昇能源科技有限责任公司 | 一种测井识别版图的建立方法及其应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4638447A (en) * | 1983-10-21 | 1987-01-20 | Mobil Oil Corporation | Method for determining consistent oil relative permeability values from dynamic displacement data |
CN107346518A (zh) * | 2016-05-06 | 2017-11-14 | 中国石油化工股份有限公司 | 致密低渗透油藏油水两相流最大渗流阻力梯度的获取方法 |
CN108487904A (zh) * | 2018-02-06 | 2018-09-04 | 中国石油大学(华东) | 一种基于图版的消除末端效应的相渗曲线校正方法 |
-
2019
- 2019-10-11 CN CN201910970956.7A patent/CN110598167B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4638447A (en) * | 1983-10-21 | 1987-01-20 | Mobil Oil Corporation | Method for determining consistent oil relative permeability values from dynamic displacement data |
CN107346518A (zh) * | 2016-05-06 | 2017-11-14 | 中国石油化工股份有限公司 | 致密低渗透油藏油水两相流最大渗流阻力梯度的获取方法 |
CN108487904A (zh) * | 2018-02-06 | 2018-09-04 | 中国石油大学(华东) | 一种基于图版的消除末端效应的相渗曲线校正方法 |
Non-Patent Citations (1)
Title |
---|
不同油藏条件下相渗曲线分析;阳晓燕;《科技技术与工程》;20121231(第14期);第3340-3343页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110598167A (zh) | 2019-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110598167B (zh) | 低渗透油藏油水相对渗透率实验数据的处理方法 | |
CN108896599B (zh) | 一种测试气水相对渗透率曲线的系统及方法 | |
CN104101564B (zh) | 一种非稳态高温高压测试低渗透岩心启动压力梯度的方法 | |
CN105547961B (zh) | 衰竭式开发砂岩凝析气藏储层中反凝析油饱和度确定方法 | |
Shubao et al. | Dynamic effect of capillary pressure in low permeability reservoirs | |
CN107462936B (zh) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 | |
CN109339774B (zh) | 基于边界层的低渗透变形介质非线性渗流数值模拟方法 | |
CN109488282B (zh) | 可动储量物性下限的确定方法 | |
Gao et al. | Apparent permeability and gas flow behavior in carboniferous shale from the Qaidam Basin, China: an experimental study | |
CN204330547U (zh) | 一种煤岩渗透率测定参数检测装置 | |
You et al. | Pore structure and limit pressure of gas slippage effect in tight sandstone | |
Zeng et al. | Experimental and theoretical characterization of the natural gas migration and accumulation mechanism in low-permeability (tight) sandstone cores | |
CN113825996A (zh) | 用于在脉冲衰减实验中确定岩心渗透率的方法和系统 | |
CN112031719A (zh) | 一种基于流动系数下启动压力的油藏开发方式优选方法 | |
CN115559715A (zh) | 一种超高压低渗透气藏储层产水评价方法 | |
Wang et al. | Gas-phase relative permeability characterization on tight-gas samples | |
CN110095584B (zh) | 一种储层油水饱和度校正方法 | |
CN116383573B (zh) | 一种基于多区相变传质渗流耦合的凝析气产能评价方法 | |
Cho et al. | Water and Oil Relative Permeability of Middle Bakken Formation: Experiments and Numerical Modeling | |
CN111208048A (zh) | 一种基于相渗测试的贾敏效应动态变化定量表征方法 | |
Tale et al. | Estimating transition time from deep filtration of particles to external cake during produced water re-injection and disposal | |
Fang et al. | Experimental study of gas flow characteristics in micro-/nano-pores in tight and shale reservoirs using microtubes under high pressure and low pressure gradients | |
Ge et al. | Predicting the relative permeability of water phase based on theory of coupled electricity-seepage and capillary bundle model | |
Meng et al. | Experimental Study on Gas–Water Seepage in High-Rank Coal and Its Research Significance | |
CN109556996A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |