CN108487904A - 一种基于图版的消除末端效应的相渗曲线校正方法 - Google Patents
一种基于图版的消除末端效应的相渗曲线校正方法 Download PDFInfo
- Publication number
- CN108487904A CN108487904A CN201810115163.2A CN201810115163A CN108487904A CN 108487904 A CN108487904 A CN 108487904A CN 201810115163 A CN201810115163 A CN 201810115163A CN 108487904 A CN108487904 A CN 108487904A
- Authority
- CN
- China
- Prior art keywords
- phase
- mutually
- end effect
- rock core
- percolation 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.)
- Granted
Links
- 230000000694 effects Effects 0.000 title claims abstract description 114
- 238000005325 percolation Methods 0.000 title claims abstract description 77
- 230000008030 elimination Effects 0.000 title claims abstract description 7
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 7
- 239000011435 rock Substances 0.000 claims abstract description 118
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 77
- 238000002474 experimental method Methods 0.000 claims abstract description 47
- 238000006073 displacement reaction Methods 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 38
- 230000035699 permeability Effects 0.000 claims abstract description 38
- 238000012512 characterization method Methods 0.000 claims abstract description 28
- 238000012937 correction Methods 0.000 claims abstract description 19
- 230000007246 mechanism Effects 0.000 claims abstract description 11
- 238000004519 manufacturing process Methods 0.000 claims abstract description 9
- 239000012530 fluid Substances 0.000 claims description 25
- 239000004519 grease Substances 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 20
- 238000004088 simulation Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 6
- 230000000149 penetrating effect Effects 0.000 claims description 6
- 238000000205 computational method Methods 0.000 claims description 4
- 238000009533 lab test Methods 0.000 claims description 4
- 238000009736 wetting Methods 0.000 claims description 4
- 230000035515 penetration Effects 0.000 abstract description 3
- 230000005514 two-phase flow Effects 0.000 abstract description 3
- 238000011161 development Methods 0.000 description 8
- 230000018109 developmental process Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000001035 drying Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012332 laboratory investigation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- 238000005406 washing Methods 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
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于图版的消除末端效应的相渗曲线校正方法,主要包括非稳态法测定岩石相对渗透率、考虑末端效应的数值模拟模型建立、考虑末端效应的相渗表征参数校正图版建立、相渗校正四部分。本发明基于末端效应的产生机理,考虑毛管力对岩心中油水两相渗流规律的影响,利用虚拟网格模拟末端效应的影响区域,并利用线性相渗对末端效应影响区域进行表征,首次建立了考虑末端效应的一维岩心两相驱替数值模拟模型,对末端效应的影响效果进行了模拟,为渗透率较低的低渗透、超低渗透岩心、致密岩心、页岩岩心的非稳态法相渗测定实验提供了较为切实可行的误差校正方法,使得相渗曲线更加贴近油藏实际,可以反映油藏中油水两相真实的渗流规律。
Description
技术领域
本发明涉及油气田开发油气渗流技术领域,特别涉及一种适用末端效应影响较大的致密岩心非稳态法测定相对渗透率曲线的测定及校正方法。
背景技术
在油气田开发领域的室内实验研究中,非稳态法是测定油水相对渗透率比较常用的方法。该方法通常是用溶剂洗净岩样,烘干后抽空饱和水,求得孔隙度。然后用油驱替水,得到束缚水饱和度,接着进行水驱油实验。在实验中记录不同驱替时间的岩心两端压差、驱出的油量和水量。用实验前测定的绝对渗透率、油水粘度、岩样的长度和截面积就可以确定相对渗透率。但是进行该实验时必须采用较高的驱动速度,使相似准数Luμw≥1(L-岩样长度,cm;u-注入速度,cm/min;μw-水的粘度,厘泊),以克服末端效应。但对于致密岩心相渗测定过程中,由于致密岩心渗透率较低,即使增大驱动压差,驱动速度也难以达到克服末端效应的要求,反而对实验设备的承压能力提出了更高的要求。为了保证致密岩心非稳态法所测定的相渗曲线可以反映岩心的实际油水渗透能力,避免相渗曲线的误差在数值模拟工作中被放大,消除末端效应对相渗曲线造成的影响,需要对常规非稳态法所测量的相渗曲线进行校正。
常规的不稳定法测量油水相对渗透率是以水驱油基本理论为基础,利用传统的JBN相渗计算方法进行相渗计算,该方法假设在水驱油过程中,油水饱和度在岩心中的分布是时间和距离的函数,没有考虑末端效应对岩心饱和度分布的影响,如图1所示。通过在水驱油过程中准确测量出恒定驱替压差下油、水流量,就可由贝克莱-列维尔特非活塞式驱油理论计算出岩心出口断面上任意时刻的含水饱和度及有效渗透率,具有测定速度快、设备简单、操作方便等优点。图2是非稳态法测量岩石相对渗透率的仪器。
非稳态法油水相渗测定的步骤如下:
1、将岩样洗净烘干、计量尺寸及干重。
2、将岩样抽空饱和水,称重法测量孔隙度,放入恒温箱内的岩心夹持器中,恒温至实验温度,一般放置时间大于4小时,水测岩石渗透率kw。
3、用油驱至束缚水饱和度,一般10PV以上,老化一定时间,一般超过10天,在实验温度下测定束缚水下油相相对渗透率。
4、以恒定的驱替压差进行水驱油,测量不同时间岩心两端的压力差、累积产油、累积产水、水驱30PV后,测定残余油下水的相对渗透率。
5、根据非稳态法油水相渗测定的基本数据和贝克莱-列维尔特驱油机理可以推导得到油水相渗的计算公式,这种方法计算的相渗曲线传统称为JBN相渗曲线,具体计算公式如下:
式中:
kro——出口端饱和度下的油相相对渗透率;
krw——出口端饱和度下的水相相对渗透率;
Swe——出口端含水饱和度;
——无因次累积注水量;
Vt——累积注水量,cm3;
——无因次累积采油量;
Vp——岩样孔隙体积,cm3;
fo(Swe)——出口端含油率;
fw(Swe)——出口端含水率;
Swi——共存水或束缚水饱和度;
μo——油相粘度,mPa·s;
μw——水相粘度,mPa·s;
I——任意时刻与初始时刻的流动能力比;
k——岩石绝对渗透率,μm2;
A——岩样渗流截面积,cm2;
L——岩样长度,cm;
Q(t)——t时刻出口端产液量,cm2/s;
Δp(t)——t时刻岩样两端出口压差,10-1MPa;
针对渗透率较低的致密岩心,由于传统的非稳态法实验中驱替速度达不到克服末端效应的要求,同时致密岩心中毛管力效应更为明显。水驱油过程中一旦见水,其驱替过程中的饱和度剖面如图3所示,与常规的稳态法相渗测量方法的假设相差较大,使得末端效应对相渗测定的影响误差较大。目前还没有较好的方法针对非稳态法致密岩心相渗测定过程中末端效应的影响进行校正。
与室内驱替实验相比,数值模拟技术是一种更为直观的描述基本渗流机理和现象的基本方法。通过数值模拟计算描述油水渗流机理,在给定的限制条件下计算油水的运动规律和状态,同时也是可进行一维油藏岩心油水驱替模拟的一种有效方法。而目前常规的一维油藏油水驱替数值模拟模型中无法考虑末端效应的影响,并不能针对有末端效应影响的一维岩心水驱油过程中的流动进行精确模拟,不能反映出有末端效应影响下一维岩心油水流动规律。
发明内容
基于上述技术问题,本发明提供一种基于图版的消除末端效应的相渗曲线校正方法。
本发明所采用的技术解决方案是:
一种基于图版的消除末端效应的相渗曲线校正方法,包括以下步骤:
(1)非稳态法测定岩石相对渗透率
利用非稳态法相渗测定基本仪器,按照规范流程对岩心相对渗透率进行测定;
(2)考虑末端效应的数值模拟模型建立
考虑毛管力的影响,根据达西定律和状态方程,可以得到油水两相的连续性方程:
辅助方程有:
So+Sw=1
pc=po-pw
初始条件有:
p(x,0)=pi
sw(x,0)=swc
岩心两端定压驱替,其边界条件有:
p|x=0=pin
p|x=L=pout
渗流微分方程
结合毛管力计算公式,利用有限差分形式对上式进行差分离散,可得:
式中:
λ——总流动系数,λ=λo+λw,其中λo=kkro/μw,λw=kkrw/μw
pw——水相压力,0.1MPa;
po——油相压力,0.1MPa;
pc——毛管力,0.1MPa;
qwv——水相源汇项;
qov——油相源汇项;
φ——岩心孔隙度;
sw——含水饱和度;
so——含油饱和度;
pi——岩心初始压力,0.1MPa;
swc——岩心束缚水饱和度;
pin——岩心入口端压力,0.1MPa;
pout——岩心出口端压力,0.1MPa;
qv——表示油水的总流量,qv=qo+qw
在网格划分过程中,在岩心末端加上一段距离的虚拟网格,来模拟网格末端效应的作用区域,并且在此区域内,考虑末端效应的产生机理,没有毛管力的影响,在末端虚拟网格中其油水两相的相对渗透能力符合线性相渗的规律;
根据以上的考虑毛管力的渗流微分方程和末端效应影响区域的虚拟网格和相渗曲线,建立了考虑末端效应影响的一维岩心油水两相驱替数值模拟模型,并利用IMPES方法进行压力和饱和度的求解,得到不同驱替时间岩心末端的出水量和出油量,利用传统的JBN相渗曲线计算方法,得到有末端效应影响的相渗曲线,并与输入的没有末端效应影响的相渗曲线进行比较,前者称为实验相渗,后者称为油藏相渗;
(3)考虑末端效应的相渗表征参数校正图版建立
根据步骤(2)所建立的数值模拟模型,输入该数值模拟模型的典型相渗曲线,该相渗曲线为理想情况下没有末端效应影响下的相渗曲线,基于步骤(2)所建立的数值模拟模型得到的产油、产水数据,并利用JBN方法计算得到的相渗曲线为有末端效应影响的相渗曲线;考虑不同的岩心长度和驱替压差,可以得到不同实验条件下的油藏相渗和实验相渗,相渗曲线根据下列公式进行拟合,每条相渗曲线利用3个参数进行拟合,定量表征相渗曲线形态;表1为相渗曲线拟合参数及拟合公式统计表;
表1
其中为归一化含水饱和度,其计算公式为:
Ewn——油藏相渗表达式中润湿相流体(水)的表征参数;
Awn——油藏相渗表达式中润湿相流体(水)的常数;
Krwend——润湿相(水)相渗曲线端点处相对渗透率;
Enwn——油藏相渗表达式中非润湿相流体(油)的表征参数;
Anwn——油藏相渗表达式中非润湿相流体(油)的常数;
Krnwend——非润湿相(油)相渗曲线端点处相对渗透率;
Ewe——实验相渗表达式中润湿相流体(水)的表征参数;
Awe——实验相渗表达式中润湿相流体(水)的常数;
Enwe——实验相渗表达式中非润湿相流体(油)的表征参数;
Anwe——实验相渗表达式中非润湿相流体(油)的常数;
(4)相渗校正
利用步骤(3)建立的参数校正图版对室内所测相渗曲线进行校正,消除末端效应导致的较大误差;首先利用测定基本仪器进行非稳态法岩心相对渗透率测量,并对室内实验测得的相渗曲线进行参数拟合,得到油水实验相渗曲线的拟合参数Enwe,Anwe,Ewe,Awe;根据室内非稳态法岩心相渗测定实验中岩心的长度和定压测量时所用的测量压差进行读图,得到对应参数的比值,进而确定油藏相渗的油水相渗曲线表征参数Enwn,Anwn,Ewn,Awn,从而确定油水相渗曲线的形态。
本发明的有益技术效果是:
第一,基于末端效应的产生机理,本发明考虑毛管力对岩心中油水两相渗流规律的影响,利用虚拟网格模拟末端效应的影响区域,并利用线性相渗对末端效应影响区域进行表征,首次建立了考虑末端效应的一维岩心两相驱替数值模拟模型,对末端效应的影响效果进行了模拟。
第二,结合室内实验的具体实施条件,根据不同岩心长度、不同的驱替压差进行了校正,揭示了不同岩心长度和驱替压差下末端效应的影响程度和影响规律,并建立了致密岩心考虑末端效应影响的相渗参数校正图版,可在传统非稳态相渗实验测定的基础上,根据校正图版对致密岩心的油水相渗曲线进行校正,减小测量过程中由于末端效应导致的误差。
第三,与以往为减小末端效应影响而改进的实验装置相比,本发明不需要耐压程度较高的实验仪器,不需要高精度的实验仪器设备,保证了实验条件的切实可行,有力地降低了实验成本,保障了实验操作过程中的安全。
第四,本方法可以进一步推广,为渗透率较低的低渗透、超低渗透岩心、致密岩心、页岩岩心的非稳态法相渗测定实验提供了较为切实可行的误差校正方法,使得相渗曲线更加贴近油藏实际,可以反映油藏中油水两相真实的渗流规律,有力地保证了后期开发方案设计及优化过程中数值模拟工作的精确度和可信度。
附图说明
下面结合附图与具体实施方式对本发明作进一步说明:
图1所示为传统非稳态法相渗测定中假设岩心中含水饱和度的分布示意图;其中Sor-残余油饱和度,Swf-水驱前缘含水饱和度,Swc-束缚水饱和度,t1、t2、t3、t4、t5分别代表不同时刻的含水饱和度剖面;
图2所示为非稳态法测量岩石相对渗透率的仪器;1-高压平流泵,2-手动计量泵,3-六通筏,4-带活塞的水容器,5-带活塞的油容器,6-压力表,7-岩心夹持器,8-油水分离器,9-压力传感器组,10-压力显示仪,11-压力记录仪;
图3所示为末端效应影响下的致密岩心中含水饱和度的实际分布示意图;其中Sor-残余油饱和度,Swf-水驱前缘含水饱和度,Swc-束缚水饱和度,t1、t2、t3、t4、t5分别代表不同时刻的含水饱和度剖面;
图4所示为本发明所提供的考虑末端效应的致密岩心非稳态法相渗测定及校正方法流程图;
图5所示为数值模拟模型中末端效应影响区域及末端虚拟网格示意图;
图6所示为数值模拟模型中末端效应影响区域线性相渗示意图;
图7所示为油相相渗表征参数Enw校正图版;
图8所示为油相相渗表征参数Anw校正图版;
图9所示为水相相渗表征参数Ew校正图版;
图10所示为水相相渗表征参数Aw校正图版;
图11所示为实验实例中1#岩心非稳态法测得的实验相渗;
图12所示为利用本发明校正后的消除末端效应影响的油藏相渗;
图13所示为校正前后岩心实验相渗和油藏相渗的对比图。
具体实施方式
本发明考虑末端效应的影响机理及效果,首次建立了考虑末端效应的一维岩心两相驱替数值模拟模型,并通过模拟计算首次建立了致密岩心考虑末端效应影响的相渗参数校正图版,可在传统非稳态相渗实验测定的基础上,对致密岩心的油水相渗曲线进行校正,减小测量过程中由于末端效应导致的误差,以在后期的开发方案设计和优选过程中开发动态预测更贴近油藏实际。
本发明的目的在于提供一种考虑末端效应影响的相渗测量误差的校正方法,减小室内非稳态岩心相渗测定实验中,尤其是致密岩心中由于末端效应导致的误差。主要解决的技术问题如下:一、根据传统的室内岩心相渗非稳态测定实验数据,传统的JBN相渗计算方法其饱和度剖面的基本假设如图1所示,而实际岩心驱替实验中由于末端效应的影响,其饱和度剖面分布与假设相差较大,导致所测量的相渗曲线有较大误差。二、目前在实验中针对减小末端效应的措施主要为两种方法,一是提高流速,但在具体实验操作中可以发现,对于中高渗岩心来说,传统的增加流速驱动可以降低末端效应的影响,而对于致密岩心,由于其渗透率较低,难以达到消除末端效应的流速。二是三段岩心法,通过在岩心末端加一段人造岩心或天然岩心,来消除末端效应的影响;但对于致密岩心相渗测定实验来说,由于流速较低,流速计量需要十分精确,三段岩心法无法保证相渗测定实验中的精确度。需要在考虑末端效应的影响下对相对渗透率进行校正,而目前在实验操作过程中没有切实可用的校正方法。三、目前利用数值模拟技术对于末端效应影响下的一维岩心油水两相驱替模拟模型尚未成熟,末端效应对相渗曲线的影响程度并不明确。需要考虑末端效应的产水机理及影响效果,建立考虑末端效应影响的一维油水两相驱替数值模拟模型,明确末端效应对相渗曲线的影响。四、传统的非稳态法相渗测定实验流程难以达到消除末端效应的要求,本发明基于数值模拟技术,对末端效应对相渗影响效果进行表征,通过对传统方法测得的相渗进行改进,从而消除末端效应对相渗曲线测定导致的误差。通过本发明技术,对传统的非稳态法岩心相渗实验所得的相渗曲线进行校正,减小末端效应对相渗曲线导致的误差,使得相渗曲线可以代表实际油藏中的油水相对渗透能力,以保证在大量的数值模拟工作准确性,为后期开发历史拟合、开发方案预测、开发方案优选奠定基础。
本发明所提供的考虑末端效应的非稳态法相渗测定及校正方法主要流程如图4所示,主要包括非稳态法测定岩石相对渗透率、考虑末端效应的数值模拟模型建立、考虑末端效应的相渗表征参数校正图版建立、相渗校正四部分。
1、利用非稳态法相渗测定基本仪器,按照规范流程对岩心相渗进行测定。
2、考虑末端效应影响的一维油水两相驱替数值模拟模型。根据末端效应的产生机理,由于毛管力的不连续性,导致岩心末端润湿相积聚,含水饱和度上升,如图3所示。考虑毛管力的影响,根据达西定律和状态方程,可以得到油水两相的连续性方程:
辅助方程有:
So+Sw=1
pc=po-pw
初始条件有:
p(x,0)=pi
sw(x,0)=swc
岩心两端定压驱替,其边界条件有:
p|x=0=pin
p|x=L=pout
渗流微分方程:
结合毛管力计算公式,利用有限差分形式对上式进行差分离散,可得:
式中:
pw——水相压力,0.1MPa;
po——油相压力,0.1MPa;
pc——毛管力,0.1MPa;
qwv——水相源汇项,
qov——油相源汇项,
φ——岩心孔隙度;
sw——含水饱和度;
so——含油饱和度;
pi——岩心初始压力,0.1MPa;
swc——岩心束缚水饱和度;
pin——岩心入口端压力,0.1MPa;
pout——岩心出口端压力,0.1MPa;
qv——表示油水的总流量,qv=qo+qw
考虑末端效应的产生机理,因为末端效应影响范围主要是岩心驱替的末端,在网格划分过程中,在岩心末端加上一段距离的虚拟网格,来模拟网格末端效应的作用区域,并且在此区域内,考虑末端效应的产生机理,没有毛管力的影响,在末端虚拟网格中其油水两相的相对渗透能力符合线性相渗的规律。
根据以上的考虑毛管力的渗流微分方程和末端效应影响区域的虚拟网格和相渗曲线,建立了考虑末端效应影响的一维岩心油水两相驱替数值模拟模型,并利用IMPES方法进行压力和饱和度的求解,可以得到不同驱替时间岩心末端的出水量和出油量,利用传统的JBN相渗曲线计算方法,可以得到有末端效应影响的相渗曲线,并可以与输入的没有末端效应影响的相渗曲线进行比较。为便于区分,前者称为实验相渗,后者称为油藏相渗。
3、油藏相渗和实验相渗的参数校正图版。根据步骤2所建立的数值模拟模型,输入该数值模拟模型的典型相渗曲线,具体为真实反映不同饱和度下渗透能力的相渗曲线,该相渗曲线为理想情况下没有末端效应影响下的相渗曲线(油藏相渗)。而通过步骤2建立的数值模拟模型和JBN方法计算得到的相渗曲线为有末端效应影响的相渗曲线(实验相渗)。考虑不同的岩心长度和驱替压差,可以得到不同实验条件下的油藏相渗和实验相渗,相渗曲线根据下列公式进行拟合,每条相渗曲线利用3个参数进行拟合,定量表征相渗曲线形态。将对应的参数比值根据模拟条件(岩心长度和压差)进行画图,如图7-10所示。
表1为相渗曲线拟合参数及拟合公式统计表。
表1
其中为归一化含水饱和度,其计算公式为:
Ewn——油藏相渗表达式中润湿相流体(水)的表征参数;
Awn——油藏相渗表达式中润湿相流体(水)的常数;
Krwend——润湿相(水)相渗曲线端点处相对渗透率;
Enwn——油藏相渗表达式中非润湿相流体(油)的表征参数;
Anwn——油藏相渗表达式中非润湿相流体(油)的常数;
Krnwend——非润湿相(油)相渗曲线端点处相对渗透率;
Ewe——实验相渗表达式中润湿相流体(水)的表征参数;
Awe——实验相渗表达式中润湿相流体(水)的常数;
Enwe——实验相渗表达式中非润湿相流体(油)的表征参数;
Anwe——实验相渗表达式中非润湿相流体(油)的常数;
4、利用参数校正图版对室内所测相渗曲线进行校正,消除末端效应导致的较大误差。首先利用室内仪器装置进行非稳态法岩心相对渗透率测量,并对室内实验测得的相渗曲线进行参数拟合,得到油水实验相渗曲线的拟合参数Enwe,Anwe,Ewe,Awe。根据室内非稳态法岩心相渗测定实验中岩心的长度和定压测量时所用的测量压差进行读图,得到对应参数的比值,进而确定油藏相渗的油水相渗曲线表征参数Enwn,Anwn,Ewn,Awn,从而确定油水相渗曲线的形态。
下面结合具体应用实例对本发明进行说明:
参考测量岩石相对渗透率规范流程对某区块1#致密岩心进行洗油、烘干、测定孔隙度等预处理工作,该岩心的基本参数如表2(岩心基本参数统计表)所示。
表2
利用图2所示的非稳态法测量岩石相对渗透率的仪器及规范流程对岩心进行相渗测定,实验中采取的定压驱替,两端压差为0.1MPa,并根据JBN方法对数据进行处理,绘制岩心的相渗曲线,如图11所示。并按照表1中相应的拟合公式利用最小二乘法对油水相渗曲线进行参数拟合。确定的拟合参数如下表3(实验相渗拟合参数统计表)所示。
表3
Enwe | Anwe | Ewe | Awe | Krnwend | Krwend |
140.5 | 1.167 | 0.06061 | 0.9223 | 1 | 0.1595 |
建立考虑末端效应的一维岩心非稳态油水两相驱替数值模拟模型,并建立油藏相渗和实验相渗的参数校正图版,如图7-图10所示。
根据岩心长度4.128cm确定相渗参数校正图版的横坐标、选择驱替压差0.1MPa的曲线,读图得到油藏相渗和实验相渗对应表征参数的比值,即图版的纵坐标。末端效应校正图版读取值统计表如下表4。
表4
根据图版求取油藏相渗的表征参数Ewn,Awn,Enwn,Anwn,例如:
校正后的油藏相渗表征参数统计表,见表5:
表5
Enwn | Anwn | Ewn | Awn |
238.136 | 1.525 | 0.041 | 0.518 |
根据油藏相渗的表征参数绘制油藏相渗,如图12所示,完成了从实验相渗到油藏相渗的校正,达到了消除末端效应影响的目的。
本发明的有益效果在于:第一,基于末端效应的产生机理,本发明考虑毛管力对岩心中油水两相渗流规律的影响,利用虚拟网格模拟末端效应的影响区域,并利用线性相渗对末端效应影响区域进行表征,首次建立了考虑末端效应的一维岩心两相驱替数值模拟模型,对末端效应的影响效果进行了模拟。第二,结合室内实验的具体实施条件,根据不同岩心长度、不同的驱替压差进行了技术,揭示了不同岩心长度和驱替压差下末端效应的影响程度和影响规律,并建立了致密岩心考虑末端效应影响的相渗参数校正图版,可在传统非稳态相渗实验测定的基础上,根据校正图版对致密岩心的油水相渗曲线进行校正,减小测量过程中由于末端效应导致的误差。第三,与以往为减小末端效应影响而改进的实验装置相比,本发明不需要耐压程度较高的实验仪器,不需要高精度的实验仪器设备,保证了实验条件的切实可行,有力地降低了实验成本,保障了实验操作过程中的安全。第四,本方法可以进一步推广,为渗透率较低的低渗透、超低渗透岩心、致密岩心、页岩岩心的非稳态法相渗测定实验提供了较为切实可行的误差校正方法,使得相渗曲线更加贴近油藏实际,可以反映油藏中油水两相真实的渗流规律,有力地保证了后期开发方案设计及优化过程中数值模拟工作的精确度和可信度。
上述方式中未述及的部分采取或借鉴已有技术即可实现。
以上对本发明的原理及实施方式进行了阐述,并提供了使用案例,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (1)
1.一种基于图版的消除末端效应的相渗曲线校正方法,其特征在于包括以下步骤:
(1)非稳态法测定岩石相对渗透率
利用非稳态法相渗测定基本仪器,按照规范流程对岩心相对渗透率进行测定;
(2)考虑末端效应的数值模拟模型建立
考虑毛管力的影响,根据达西定律和状态方程,可以得到油水两相的连续性方程:
辅助方程有:
So+Sw=1
pc=po-pw
初始条件有:
p(x,0)=pi
sw(x,0)=swc
岩心两端定压驱替,其边界条件有:
p|x=0=pin
p|x=L=pout
渗流微分方程
结合毛管力计算公式,利用有限差分形式对上式进行差分离散,可得:
式中:
λ——总流动系数,λ=λo+λw,其中λo=kkro/μw,λw=kkrw/μw
pw——水相压力,0.1MPa;
po——油相压力,0.1MPa;
pc——毛管力,0.1MPa;
qwv——水相源汇项;
qov——油相源汇项;
φ——岩心孔隙度;
sw——含水饱和度;
so——含油饱和度;
pi——岩心初始压力,0.1MPa;
swc——岩心束缚水饱和度;
pin——岩心入口端压力,0.1MPa;
pout——岩心出口端压力,0.1MPa;
qv——表示油水的总流量,qv=qo+qw
在网格划分过程中,在岩心末端加上一段距离的虚拟网格,来模拟网格末端效应的作用区域,并且在此区域内,考虑末端效应的产生机理,没有毛管力的影响,在末端虚拟网格中其油水两相的相对渗透能力符合线性相渗的规律;
根据以上的考虑毛管力的渗流微分方程和末端效应影响区域的虚拟网格和相渗曲线,建立了考虑末端效应影响的一维岩心油水两相驱替数值模拟模型,并利用IMPES方法进行压力和饱和度的求解,得到不同驱替时间岩心末端的出水量和出油量,利用传统的JBN相渗曲线计算方法,得到有末端效应影响的相渗曲线,并与输入的没有末端效应影响的相渗曲线进行比较,前者称为实验相渗,后者称为油藏相渗;
(3)考虑末端效应的相渗表征参数校正图版建立
根据步骤(2)所建立的数值模拟模型,输入该数值模拟模型的典型相渗曲线,该相渗曲线为理想情况下没有末端效应影响下的相渗曲线,通过步骤(2)建立的数值模拟模型和JBN方法计算得到的相渗曲线为有末端效应影响的相渗曲线;考虑不同的岩心长度和驱替压差,可以得到不同实验条件下的油藏相渗和实验相渗,相渗曲线根据下列公式进行拟合,每条相渗曲线利用3个参数进行拟合,定量表征相渗曲线形态;表1为相渗曲线拟合参数及拟合公式统计表;
表1
其中为归一化含水饱和度,其计算公式为:
Ewn——油藏相渗表达式中润湿相流体的表征参数;
Awn——油藏相渗表达式中润湿相流体的常数;
Krwend——润湿相相渗曲线端点处相对渗透率;
Enwn——油藏相渗表达式中非润湿相流体的表征参数;
Anwn——油藏相渗表达式中非润湿相流体的常数;
Krnwend——非润湿相相渗曲线端点处相对渗透率;
Ewe——实验相渗表达式中润湿相流体的表征参数;
Awe——实验相渗表达式中润湿相流体的常数;
Enwe——实验相渗表达式中非润湿相流体的表征参数;
Anwe——实验相渗表达式中非润湿相流体的常数;
(4)相渗校正
利用步骤(3)建立的参数校正图版对室内所测相渗曲线进行校正,消除末端效应导致的较大误差;首先利用测定基本仪器进行非稳态法岩心相对渗透率测量,并对室内实验测得的相渗曲线进行参数拟合,得到油水实验相渗曲线的拟合参数Enwe,Anwe,Ewe,Awe;根据室内非稳态法岩心相渗测定实验中岩心的长度和定压测量时所用的测量压差进行读图,得到对应参数的比值,进而确定油藏相渗的油水相渗曲线表征参数Enwn,Anwn,Ewn,Awn,从而确定油水相渗曲线的形态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810115163.2A CN108487904B (zh) | 2018-02-06 | 2018-02-06 | 一种基于图版的消除末端效应的相渗曲线校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810115163.2A CN108487904B (zh) | 2018-02-06 | 2018-02-06 | 一种基于图版的消除末端效应的相渗曲线校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108487904A true CN108487904A (zh) | 2018-09-04 |
CN108487904B CN108487904B (zh) | 2020-05-19 |
Family
ID=63344336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810115163.2A Active CN108487904B (zh) | 2018-02-06 | 2018-02-06 | 一种基于图版的消除末端效应的相渗曲线校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108487904B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110424944A (zh) * | 2019-07-22 | 2019-11-08 | 中国海洋石油集团有限公司 | 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法 |
CN110598167A (zh) * | 2019-10-11 | 2019-12-20 | 中国石油化工股份有限公司 | 低渗透油藏油水相对渗透率实验数据的处理方法 |
CN111220509A (zh) * | 2020-01-20 | 2020-06-02 | 中国石油天然气股份有限公司 | 一种考虑渗透率时变性的油水相对渗透率曲线修正方法 |
CN112179826A (zh) * | 2020-09-14 | 2021-01-05 | 中国石油大学(华东) | 一种基于时域反射技术的测量高温高压岩心动态毛管力的装置及实验方法 |
CN113295594A (zh) * | 2021-05-20 | 2021-08-24 | 中国石油大学(北京) | 一种基于分流模型与ct扫描的相对渗透率测定系统及方法 |
CN113295592A (zh) * | 2021-05-20 | 2021-08-24 | 中国石油大学(北京) | 一种低渗透岩石的相对渗透率测定系统及方法 |
US11408811B2 (en) | 2020-02-04 | 2022-08-09 | Saudi Arabian Oil Company | Methods and systems for determining residual fluid saturation of a subsurface formation |
CN115452681A (zh) * | 2022-11-11 | 2022-12-09 | 中国石油大学(华东) | 一种储层油气相对渗透率确定方法 |
CN117664804A (zh) * | 2024-02-01 | 2024-03-08 | 中国地质大学(北京) | 一种考虑端面效应的相对渗透率曲线校正方法及其应用 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4622643A (en) * | 1983-10-21 | 1986-11-11 | Mobil Oil Corporation | Method for determining consistent water relative permeability values from dynamic displacement data |
CN102915406A (zh) * | 2012-10-26 | 2013-02-06 | 中国石油大学(华东) | 径向流条件下油水相对渗透率曲线的计算方法 |
CN103573234A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种确定完整油水相对渗透率曲线的方法 |
CN104102802A (zh) * | 2013-04-03 | 2014-10-15 | 中国石油化工股份有限公司 | 不同驱替压力梯度条件下油水相对渗透率曲线的表征方法 |
CN104179499A (zh) * | 2013-05-27 | 2014-12-03 | 中国石油化工股份有限公司 | 考虑油藏参数时变的数值模拟方法 |
CN105239976A (zh) * | 2015-11-04 | 2016-01-13 | 中国海洋石油总公司 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
CN107346518A (zh) * | 2016-05-06 | 2017-11-14 | 中国石油化工股份有限公司 | 致密低渗透油藏油水两相流最大渗流阻力梯度的获取方法 |
-
2018
- 2018-02-06 CN CN201810115163.2A patent/CN108487904B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4622643A (en) * | 1983-10-21 | 1986-11-11 | Mobil Oil Corporation | Method for determining consistent water relative permeability values from dynamic displacement data |
CN103573234A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种确定完整油水相对渗透率曲线的方法 |
CN102915406A (zh) * | 2012-10-26 | 2013-02-06 | 中国石油大学(华东) | 径向流条件下油水相对渗透率曲线的计算方法 |
CN104102802A (zh) * | 2013-04-03 | 2014-10-15 | 中国石油化工股份有限公司 | 不同驱替压力梯度条件下油水相对渗透率曲线的表征方法 |
CN104179499A (zh) * | 2013-05-27 | 2014-12-03 | 中国石油化工股份有限公司 | 考虑油藏参数时变的数值模拟方法 |
CN105239976A (zh) * | 2015-11-04 | 2016-01-13 | 中国海洋石油总公司 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
CN107346518A (zh) * | 2016-05-06 | 2017-11-14 | 中国石油化工股份有限公司 | 致密低渗透油藏油水两相流最大渗流阻力梯度的获取方法 |
Non-Patent Citations (1)
Title |
---|
雷霄等: "一种改进的油水相渗实验末端效应校正方法_", 《中国海上油气》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110424944A (zh) * | 2019-07-22 | 2019-11-08 | 中国海洋石油集团有限公司 | 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法 |
CN110424944B (zh) * | 2019-07-22 | 2022-08-23 | 中国海洋石油集团有限公司 | 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法 |
CN110598167A (zh) * | 2019-10-11 | 2019-12-20 | 中国石油化工股份有限公司 | 低渗透油藏油水相对渗透率实验数据的处理方法 |
CN110598167B (zh) * | 2019-10-11 | 2022-10-28 | 中国石油化工股份有限公司 | 低渗透油藏油水相对渗透率实验数据的处理方法 |
CN111220509A (zh) * | 2020-01-20 | 2020-06-02 | 中国石油天然气股份有限公司 | 一种考虑渗透率时变性的油水相对渗透率曲线修正方法 |
US11408811B2 (en) | 2020-02-04 | 2022-08-09 | Saudi Arabian Oil Company | Methods and systems for determining residual fluid saturation of a subsurface formation |
CN112179826A (zh) * | 2020-09-14 | 2021-01-05 | 中国石油大学(华东) | 一种基于时域反射技术的测量高温高压岩心动态毛管力的装置及实验方法 |
CN113295594B (zh) * | 2021-05-20 | 2022-07-12 | 中国石油大学(北京) | 一种基于分流模型与ct扫描的相对渗透率测定系统及方法 |
CN113295592B (zh) * | 2021-05-20 | 2022-06-07 | 中国石油大学(北京) | 一种低渗透岩石的相对渗透率测定系统及方法 |
CN113295592A (zh) * | 2021-05-20 | 2021-08-24 | 中国石油大学(北京) | 一种低渗透岩石的相对渗透率测定系统及方法 |
CN113295594A (zh) * | 2021-05-20 | 2021-08-24 | 中国石油大学(北京) | 一种基于分流模型与ct扫描的相对渗透率测定系统及方法 |
CN115452681A (zh) * | 2022-11-11 | 2022-12-09 | 中国石油大学(华东) | 一种储层油气相对渗透率确定方法 |
CN117664804A (zh) * | 2024-02-01 | 2024-03-08 | 中国地质大学(北京) | 一种考虑端面效应的相对渗透率曲线校正方法及其应用 |
CN117664804B (zh) * | 2024-02-01 | 2024-04-26 | 中国地质大学(北京) | 一种考虑端面效应的相对渗透率曲线校正方法及其应用 |
Also Published As
Publication number | Publication date |
---|---|
CN108487904B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108487904A (zh) | 一种基于图版的消除末端效应的相渗曲线校正方法 | |
CN108896599B (zh) | 一种测试气水相对渗透率曲线的系统及方法 | |
US7406857B2 (en) | Electronic humidity chamber for vapor desorption to determine high capillary pressures | |
Walls et al. | Effects of pressure and partial water saturation on gas permeability in tight sands: experimental results | |
CN208171813U (zh) | 一种多功能的渗透率测试装置 | |
CN106814018A (zh) | 一种致密岩石气相相对渗透率测量装置及方法 | |
CN210264648U (zh) | 一种多功能的岩心驱替装置 | |
CN104359819A (zh) | 低渗致密岩心气水相对渗透率测定装置及测定方法 | |
CN109184644B (zh) | 一种考虑聚合物非牛顿性和渗流附加阻力的早期注聚效果评价方法 | |
CN105181739B (zh) | 一种检测固体推进剂体膨胀系数的方法 | |
CN107121370A (zh) | 致密油层水相圈闭损害实验评价方法 | |
CN109682850A (zh) | 一种在线渗吸实验用核磁共振测试装置及实验方法 | |
CN112986097B (zh) | 一种确定致密储层稳态法相对渗透率曲线的实验测量方法 | |
CN104297126A (zh) | 低渗透储层气体渗流启动压力梯度测量装置及测量方法 | |
CN108489878A (zh) | 一种基于数值模拟迭代消除末端效应的相渗曲线校正方法 | |
CN207586089U (zh) | 一种稠油启动压力梯度以及渗流规律测量装置 | |
CN110309611B (zh) | 基于气水厚度分布的气水两相渗流规律预测方法及系统 | |
Farahani et al. | Sensitivity of unsteady-state gas-water relative permeability to experimental artefacts and interpretation techniques; case study from a gas reservoir in south Iran | |
CN109085104A (zh) | 一种致密气藏的气体相对渗透率的确定方法 | |
Christoffersen et al. | Gas/oil capillary pressure of chalk at elevated pressures | |
CN111963149B (zh) | 一种考虑滞地液量增压的压裂后地层压力求取方法 | |
CN106526079B (zh) | 一种研究致密砂岩孔喉结构动态变化的方法 | |
CN108153944A (zh) | 一种应力敏感储层产水气井天然气单井控制储量确定方法 | |
CN110095584B (zh) | 一种储层油水饱和度校正方法 | |
Liang et al. | Acquisition and analysis of transient data through unsteady-state core flooding experiments |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230601 Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Patentee after: CHINA PETROLEUM & CHEMICAL Corp. Patentee after: CHINA PETROLEUM & CHEMICAL CORPORATION EXPLORATION & PRODUCTION Research Institute Address before: 266580 No. 66 Changjiang West Road, Huangdao economic and Technological Development Zone, Qingdao, Shandong Patentee before: CHINA University OF PETROLEUM (EAST CHINA) |