CN112696194A - 一种确定特高含水油藏可动油饱和度的方法 - Google Patents

一种确定特高含水油藏可动油饱和度的方法 Download PDF

Info

Publication number
CN112696194A
CN112696194A CN201911007527.6A CN201911007527A CN112696194A CN 112696194 A CN112696194 A CN 112696194A CN 201911007527 A CN201911007527 A CN 201911007527A CN 112696194 A CN112696194 A CN 112696194A
Authority
CN
China
Prior art keywords
water
oil
saturation
relative permeability
determining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201911007527.6A
Other languages
English (en)
Inventor
王明
许寻
熊运斌
莫冰
张蕴曼
杨瑞敏
谢爱华
席伟安
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201911007527.6A priority Critical patent/CN112696194A/zh
Publication of CN112696194A publication Critical patent/CN112696194A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes 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)
  • Geophysics (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明属于油气田开发技术领域,具体涉及一种确定特高含水油藏可动油饱和度的方法。本发明的确定特高含水油藏可动油饱和度的方法主要包括以下步骤:根据高倍数水驱后的目标储层岩心样品的油相相对渗透率曲线以及水相相对渗透率曲线确定高倍数水驱后的残余油饱和度,然后根据高倍数水驱前的残余油饱和度以及高倍数水驱后的残余油饱和度确定可动油饱和度。本发明的方法精度高,可以反映特高含水期的渗流特征,并且可以应用于油田现场,为注水开发油藏方案的调整提供依据。

Description

一种确定特高含水油藏可动油饱和度的方法
技术领域
本发明属于油气田开发技术领域,具体涉及一种确定特高含水油藏可动油饱和度的方法。
背景技术
油田经过长期的注水开发,有的油田现进入到特高含水开发期,寻找地下剩余油的分布规律并制定相应的挖潜措施已成为当务之急。剩余油包括残余油及可动油,其中残余油指水驱后仍然不能采出而残留于油层孔隙中的原油,可动油是指在一定的采油工艺条件下在储层中可以流动且部分地从油层中产出的油。虽然水驱后的可动油不能完全从储层中生产出来,但水驱后剩余油挖潜的主要目的是降低水驱后的可动油饱和度,因此可动油饱和度是油藏下一步调整挖潜的对象,也是油藏评价中的一个重要参数。
目前水驱后可动油饱和度研究主要采用室内岩心,结合核磁共振技术、常规压汞技术、油驱水和水驱油实验,只能得到岩心剩余可动油百分数。张树宝等在《关于剩余油饱和度分布的研究方法》(断块油气田,2000年第03期,19~22页)一文中公开了一种确定可动油饱和度的方法,包括以下步骤:在半对数坐标上作出油水相对渗透率比值与含水饱和度的关系曲线,进行线性回归,得到Kro/Krw=ae-bSw;根据油水分流理论求含水率fw;根据含水率求含水饱和度Sw;根据含水饱和度求水驱剩余油饱和度;根据剩余油饱和度求水驱可动油饱和度。在该方法的计算过程中存在相互求取迭代的缺陷,另外精度低。并且上述方法不能直接应用于油田现场,不便于油藏工程方案编制实施。
发明内容
本发明的目的在于提供一种确定特高含水油藏可动油饱和度的方法,该方法能够提高得到的可动油饱和度的精度。
为实现上述目的,本发明的确定特高含水油藏可动油饱和度的方法采用的技术方案为:
一种确定特高含水油藏可动油饱和度的方法,包括以下步骤:
(1)将目标储层的岩心样品洗油、干燥、然后饱和水,再饱和油制造束缚水,得到一定含水饱和度的岩心样品;将制造束缚水后的岩心样品进行高倍数水驱,测定岩心样品高倍数水驱后的油相相对渗透率、水相相对渗透率以及含水饱和度;所述高倍数水驱的驱替倍数至少为1000pv;
(2)重复步骤(1)测定不同含水饱和度下岩心样品水驱后的油相相对渗透率、水相相对渗透率数据以及含水饱和度,然后对式(Ⅰ)进行拟合,确定系数a、b、c,得到岩心样品的相对渗透率曲线:
ln(Kro/Krw)=aSw 2+bSw+c (Ⅰ)
式(Ⅰ)中:Kro为油相相对渗透率,Krw为水相相对渗透率,Sw为含水饱和度,a、b和c为无单位系数;
(3)根据拟合后的公式(Ⅰ)对不同含水饱和度下的油相相对渗透率以及水相相对渗透率进行校正,确定特高含水油藏水驱后的残余油饱和度Sorh=1-Sw
(4)利用高倍数水驱后的残余油饱和度Sorh和特高含水油藏高倍数水驱前的残余油饱和度Sorq确定特高含水油藏的可动油饱和度Sd:Sd=Sorq-Sorh
特高含水期表现出的渗流规律、特征不同于中高含水阶段,特高含水期的相对渗透率比值曲线明显呈现非线性,因此特高含水阶段水驱特征曲线都出现“上翘”现象,当水驱特征曲线出现“上翘”现象时,传统油藏工程的方法不再适用。本发明利用目标储层在高倍数水驱前后的相对渗透率数据,定量表征了高倍数水驱前后的残余油饱和度,从而确定特高含水油藏可动油饱和度。本发明的确定可动油饱和度的方法精度高,适用于特高含水油藏,能够直接应用于油田现场,便于油藏工程方案编制实施。
步骤(4)中所述特高含水油藏高倍数水驱前的残余油饱和度Sorq由以下方法确定:
1)获取目标储层的物性参数和流体特征参数;所述物性参数包括渗透率和孔隙度;所述流体特征参数包括地层原油粘度和地层水粘度;
2)利用所述物性参数和流体特征参数,测定目标储层岩心样品在不同含水饱和度下的油相相对渗透率数据和水相相对渗透率数据;
3)根据步骤2)所得油相相对渗透率和水相相对渗透率确定相对渗透率曲线,然后利用相对渗透率曲线确定特高含水油藏高倍数水驱前的残余油饱和度Sorq
优选的,步骤(1)中所述高倍数水驱的驱替倍数为1000pv。
优选的,步骤2)中测定不同含水饱和度下的油相相对渗透率和水相相对渗透率时的水驱倍数为20~50pv。
优选的,所用岩心样品的孔隙度为0.2~0.3。
附图说明
图1为本发明的实施例中高倍数水驱前归一化后的油水相相对渗透率图;
图2为本发明的实施例中高倍数水驱后归一化后的油水相相对渗透率图;
图1和图2中,曲线ⅰ为油相相对渗透率曲线,曲线ⅱ为水相相对渗透率曲线。
具体实施方式
下面结合具体实施例对本发明的技术方案作进一步说明。
实施例1
本实施例采用濮城油田西区沙二下储层为例,采用以下方法确定储层的可动油饱和度:
一、确定高倍数水驱前的残余油饱和度
(1)利用濮城油田沙二下储层的岩心样品(渗透率k为115mD),测定其孔隙度为0.214,长度为5cm,半径为2.5cm;
(2)利用濮城油田沙二下储层的油水样品,测定地层原油粘度μo为6.9mPa·s,地层水粘度μw为1.34mPa·s;
(3)将经过吸油、烘干后的岩心样品放入岩心夹持器内饱和水,然后采用与地层原油粘度相同的模拟油驱水制造束缚水,测定束缚水状态下的油相相对渗透率。然后水驱(水驱倍数为20~50pv),在水驱过程中记录见水时间、见水时的累积产油量、累积产液量和驱替压差。
(4)利用以下公式可得到岩心样品在不同饱和度下的油相相对渗透率和水相相对渗透率:
Figure BDA0002243207610000031
Figure BDA0002243207610000032
Figure BDA0002243207610000033
Figure BDA0002243207610000034
Figure BDA0002243207610000035
其中:fo(Sw)—含油率的数值,用小数表示;
Figure BDA0002243207610000036
—无因次累积采油量的数值,以孔隙体积(V,
Figure BDA0002243207610000037
r为岩心样品的半径,L为岩心样品的长度,
Figure BDA0002243207610000038
为岩心样品的孔隙度。)的倍数表示;
Figure BDA0002243207610000039
—无因次累积采液量的数值,以孔隙体积的倍数表示;
Kro—油相相对渗透率的数值,用小数表示;
Krw—水相相对渗透率的数值,用小数表示;
I—相对注入能力的数值,又称流动能力比;
Q(t)—t时刻岩样出口端面产液流量的数值,单位为立方厘米每秒(cm3/s);
Qo—初始时刻岩样出口端面产液流量的数值,单位为立方厘米每秒(cm3/s);
Δpo—初始驱动压差的数值,单位为兆帕(MPa);
Δp(t)—t时刻驱动压差的数值,单位为兆帕(MPa);
Swe—岩样出口端面含水饱和度的数值,用小数表示;
Sws—束缚水饱和度的数值,用小数表示。
通过对储层中岩心样品的相渗实验,测试岩心样品在不同含水饱和度下的油相相对渗透率和水相相对渗透率,测试结果如表1所示。
表1高倍数水驱前相渗实验结果
Figure BDA0002243207610000041
Figure BDA0002243207610000051
(5)根据表1中的相对渗透率数据得相对渗透率曲线(如图1所示)。根据表1、图1可知,高倍数水驱前的含水饱和度Sw1为0.7288,则高倍数水驱前的残余油饱和度Sorq=1-Sw1=0.2712。
二、确定高倍数水驱后的残余油饱和度
(1)将上述岩心样品洗油后烘干,然后饱和水,之后饱和油制造束缚水。
(2)对制造束缚水后的岩心样品进行高倍数水驱,驱替倍数为1000pv,在水驱过程中记录见水时间、见水时的累积产油量、累积产液量和驱替压差,然后根据确定高倍数水驱前的残余油饱和度的步骤(4)中的各公式得到不同含水饱和度下的油相相对渗透率和水相相对渗透率如表2所示。
表2高倍数水驱后相渗实验结果
驱替倍数/pv 含水饱和度/% 油相相对渗透率(K<sub>ro</sub>) 水相相对渗透率(K<sub>rw</sub>)
0.34 0.385 1.0000 0.0000
1.25 0.436 0.4980 0.001
2.84 0.569 0.1613 0.012
5.66 0.664 0.0100 0.022
8.04 0.677 0.0072 0.025
11.22 0.698 0.0038 0.028
24.92 0.719 0.0021 0.032
30.68 0.732 0.0014 0.034
107.56 0.74 0.0011 0.035
353.79 0.755 0.0007 0.038
588.10 0.762 0.0005 0.039
822.41 0.767 0.0005 0.04
1001.92 0.778 0.0001 0.042
(3)将表(2)中的数据,依据式(Ⅰ)通过数据拟合可得反映特高含水期渗流特征的相渗公式(Ⅱ)。
ln(Kro/Krw)=aSw 2+bSw+c (Ⅰ)
ln(Kro/Krw)=-131.692Sw 2+164.615Sw-53.747 (Ⅱ)
(4)依据步骤(3)中的公式(Ⅱ)得到校正后的特高含水期的相对渗透率曲线数据如表3所示,校正后的特高含水期的相渗曲线如图2所示。
表3校正后的相对渗透率数据
含水饱和度/% 油相相对渗透率(K<sub>ro</sub>) 水相相对渗透率(K<sub>rw</sub>)
0.385 1 0
0.436 0.348 0.001
0.569 0.122 0.012
0.664 0.01 0.022
0.677 0.007 0.025
0.698 0.004 0.028
0.719 0.0025 0.032
0.732 0.002 0.034
0.74 0.0015 0.035
0.755 0.0009 0.038
0.762 0.0008 0.039
0.767 0.0006 0.04
0.778 0.0002 0.042
0.823 0.0001 0.042
根据图2可知,高倍数水驱后的含水饱和度Sw2为0.823,高倍数水驱后的残余油饱和度Sorh=1-Sw2=0.177。
三、确定可动油饱和度
根据高倍水驱前后的残余油饱和度,可知可动油饱和度Sd=Sorq-Sorh=0.2712-0.177=0.0942,即采用水驱方式可以流动且从油层中产出的油的饱和度为0.0942。

Claims (10)

1.一种确定特高含水油藏可动油饱和度的方法,其特征在于,包括以下步骤:
(1)将目标储层的岩心样品洗油、干燥、然后饱和水,再饱和油制造束缚水,得到一定含水饱和度的岩心样品;将制造束缚水后的岩心样品进行高倍数水驱,测定岩心样品高倍数水驱后的油相相对渗透率、水相相对渗透率以及含水饱和度;所述高倍数水驱的驱替倍数至少为1000pv;
(2)重复步骤(1)测定不同含水饱和度下岩心样品水驱后的油相相对渗透率、水相相对渗透率数据以及含水饱和度,然后对式(Ⅰ)进行拟合,确定系数a、b、c,得到岩心样品的相对渗透率曲线:
ln(Kro/Krw)=aSw 2+bSw+c (Ⅰ)
式(Ⅰ)中:Kro为油相相对渗透率,Krw为水相相对渗透率,Sw为含水饱和度,a、b和c为无单位系数;
(3)根据拟合后的公式(Ⅰ)对不同含水饱和度下的油相相对渗透率以及水相相对渗透率进行校正,确定特高含水油藏水驱后的残余油饱和度Sorh=1-Sw
(4)利用高倍数水驱后的残余油饱和度Sorh和特高含水油藏高倍数水驱前的残余油饱和度Sorq确定特高含水油藏的可动油饱和度Sd:Sd=Sorq-Sorh
2.根据权利要求1所述的确定特高含水油藏可动油饱和度的方法,其特征在于,步骤(4)中所述特高含水油藏高倍数水驱前的残余油饱和度Sorq由以下方法确定:
1)获取目标储层的物性参数和流体特征参数;所述物性参数包括渗透率和孔隙度;所述流体特征参数包括地层原油粘度和地层水粘度;
2)利用所述物性参数和流体特征参数,测定目标储层岩心样品在不同含水饱和度下的油相相对渗透率和水相相对渗透率;
3)根据步骤2)所得油相相对渗透率和水相相对渗透率确定相对渗透率曲线,然后利用相对渗透率曲线确定特高含水油藏高倍数水驱前的残余油饱和度Sorq
3.根据权利要求1或2所述的确定特高含水油藏可动油饱和度的方法,其特征在于,步骤(1)中所述高倍数水驱的驱替倍数为1000pv。
4.根据权利要求2所述的确定特高含水油藏可动油饱和度的方法,其特征在于,步骤2)中测定不同含水饱和度下的油相相对渗透率和水相相对渗透率时的水驱倍数为20~50pv。
5.根据权利要求1、2或4所述的确定特高含水油藏可动油饱和度的方法,其特征在于,所用岩心样品的孔隙度为0.2~0.3。
6.根据权利要求1、2或4所述的确定特高含水油藏可动油饱和度的方法,其特征在于,不同含水饱和度下的油相相对渗透率Kro根据以下公式计算得到:
Figure FDA0002243207600000021
其中fo(Sw)为含油率,
Figure FDA0002243207600000022
为无因此累积采油量的数值,I为相对注入能力的数值。
7.根据权利要求6所述的确定特高含水油藏可动油饱和度的方法,其特征在于,所述含油率fo(Sw)根据以下公式计算得到:
Figure FDA0002243207600000023
其中
Figure FDA0002243207600000024
为无因此累积采油量的数值,
Figure FDA0002243207600000025
为无因次累积采油量的数值。
8.根据权利要求6所述的确定特高含水油藏可动油饱和度的方法,其特征在于,所述相对注入能力的数值I根据以下公式计算得到:
Figure FDA0002243207600000026
其中Q(t)为t时刻岩样出口端面产液流量的数值,Qo为初始时刻岩心样品出口端面产液流量的数值,Δpo为初始驱动压差的数值,Δp(t)为t时刻驱动压差的数值。
9.根据权利要求1或2任一项所述的确定特高含水油藏可动油饱和度的方法,其特征在于,所述不同含水饱和度下的水相相对渗透率Krw根据以下公式计算得到:
Figure FDA0002243207600000027
fo(Sw)为含油率,μw为地层水粘度,μo为地层原油粘度,Kro为油相相对渗透率。
10.根据权利要求1所述的确定特高含水油藏可动油饱和度的方法,其特征在于,式(Ⅰ)拟合后为:ln(Kro/Krw)=-131.692Sw 2+164.615Sw-53.747,其中Kro为油相相对渗透率,Krw为水相相对渗透率,Sw为含水饱和度。
CN201911007527.6A 2019-10-22 2019-10-22 一种确定特高含水油藏可动油饱和度的方法 Pending CN112696194A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911007527.6A CN112696194A (zh) 2019-10-22 2019-10-22 一种确定特高含水油藏可动油饱和度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911007527.6A CN112696194A (zh) 2019-10-22 2019-10-22 一种确定特高含水油藏可动油饱和度的方法

Publications (1)

Publication Number Publication Date
CN112696194A true CN112696194A (zh) 2021-04-23

Family

ID=75504716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911007527.6A Pending CN112696194A (zh) 2019-10-22 2019-10-22 一种确定特高含水油藏可动油饱和度的方法

Country Status (1)

Country Link
CN (1) CN112696194A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916740A (zh) * 2021-09-27 2022-01-11 中国海洋石油集团有限公司 一种测量中高渗岩心水驱轻油相渗的实验方法及装置
CN117648523A (zh) * 2024-01-29 2024-03-05 成都英沃信科技有限公司 一种有水气藏动态储量及水侵常数计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4349737A (en) * 1980-07-14 1982-09-14 Standard Oil Company (Indiana) Determination of movable oil saturations
US20080157584A1 (en) * 2006-12-29 2008-07-03 Kieschnick John A System and method for identifying productive gas shale formations
CN106600693A (zh) * 2016-12-14 2017-04-26 西安石油大学 含油饱和度的四维建模方法
CN110095584A (zh) * 2018-01-31 2019-08-06 中国石油化工股份有限公司 一种储层油水饱和度校正方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4349737A (en) * 1980-07-14 1982-09-14 Standard Oil Company (Indiana) Determination of movable oil saturations
US20080157584A1 (en) * 2006-12-29 2008-07-03 Kieschnick John A System and method for identifying productive gas shale formations
CN106600693A (zh) * 2016-12-14 2017-04-26 西安石油大学 含油饱和度的四维建模方法
CN110095584A (zh) * 2018-01-31 2019-08-06 中国石油化工股份有限公司 一种储层油水饱和度校正方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916740A (zh) * 2021-09-27 2022-01-11 中国海洋石油集团有限公司 一种测量中高渗岩心水驱轻油相渗的实验方法及装置
CN117648523A (zh) * 2024-01-29 2024-03-05 成都英沃信科技有限公司 一种有水气藏动态储量及水侵常数计算方法
CN117648523B (zh) * 2024-01-29 2024-04-05 成都英沃信科技有限公司 一种有水气藏动态储量及水侵常数计算方法

Similar Documents

Publication Publication Date Title
CN106651610B (zh) 一种浅层超低渗砂岩油藏注水开发动态分析方法
CN104018829B (zh) 一种利用煤层气井生产数据测量气水相渗曲线的方法
CN104834807B (zh) 一种基于分形理论的应力敏感储层相对渗透率计算方法
CN110306960B (zh) 双孔双渗介质储层的压裂液渗吸增产模拟方法及装置
CN107038268B (zh) 一种确定非均质储层五点井网水驱波及系数的方法
CN111353205B (zh) 用于致密气藏产水气井地层压力和动态产能的计算方法
CN103912269A (zh) 页岩气储层地层破裂压力梯度录井确定方法
CN112696194A (zh) 一种确定特高含水油藏可动油饱和度的方法
CN112031719A (zh) 一种基于流动系数下启动压力的油藏开发方式优选方法
CN111577264B (zh) 裂缝孔隙型油藏水平井产能预测方法及装置
CN110095584B (zh) 一种储层油水饱和度校正方法
CN111535803B (zh) 一种油田化学调剖调驱药剂合理注入压力预测方法
CN116201538B (zh) 一种基于采出程度的全生命周期储层损害评价方法
CN105300865A (zh) 一种钻井液滤饼渗透率测试方法
CN113136185A (zh) 一种低温高矿化度油藏用有机堵水冻胶
CN106501300A (zh) 一种高含水致密凝析气藏非平衡相变的测试方法
CN112989721B (zh) 一种致密油藏体积压裂水平井改造体积的快速计算方法
CN113312798B (zh) 补充水驱油实验驱替倍数与驱油效率早期缺失数据的方法
CN109932498A (zh) 一种分析酸蚀裂缝壁面稳定性的实验方法
RU2648135C1 (ru) Способ разработки нефтяного месторождения
CN116337705A (zh) 一种确定砾岩储层气驱下限实验方法
CN111157073B (zh) 聚合物溶液在多孔介质中的滞留信息的测定方法及***
CN110965992B (zh) 一种确定地层含气原油粘度的方法
Sensoy Use of nanoparticles for maintaining shale stability
CN106677764B (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