CN112685971B - 亲水油藏的参数解释方法及装置 - Google Patents

亲水油藏的参数解释方法及装置 Download PDF

Info

Publication number
CN112685971B
CN112685971B CN202011570188.5A CN202011570188A CN112685971B CN 112685971 B CN112685971 B CN 112685971B CN 202011570188 A CN202011570188 A CN 202011570188A CN 112685971 B CN112685971 B CN 112685971B
Authority
CN
China
Prior art keywords
parameter
oil phase
water
phase
polymer
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
Application number
CN202011570188.5A
Other languages
English (en)
Other versions
CN112685971A (zh
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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202011570188.5A priority Critical patent/CN112685971B/zh
Publication of CN112685971A publication Critical patent/CN112685971A/zh
Application granted granted Critical
Publication of CN112685971B publication Critical patent/CN112685971B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Instructional Devices (AREA)
  • Treatment Of Water By Ion Exchange (AREA)

Abstract

本申请公开一种亲水油藏的参数解释方法及装置。该方法包括:获得亲水油藏的第一参数集中各参数的参数值;设定亲水油藏的第二参数集中各参数的参数值;利用预设的聚合物驱模型求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值;根据前述求解结果计算油产量和水产量,再计算含水率,之后对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合,如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足预设精度要求,则重新设定第二参数集中至少一个参数的参数值,并执行后续步骤。本申请公开的方案能够提高亲水油藏的参数解释精度。

Description

亲水油藏的参数解释方法及装置
技术领域
本申请属于石油开采技术领域,尤其涉及亲水油藏的参数解释方法及装置。
背景技术
聚合物驱是指向地层中注入聚合物进行驱油的一种增产措施,已经在油田得到大规模工业化应用。聚合物驱的数值试井过程是分析和推算油藏参数的过程。
申请人发现,目前针对亲水油藏的参数解释,其解释结果往往存在较大的偏差。
发明内容
有鉴于此,本申请的目的在于提供一种亲水油藏的参数解释方法及装置,以提高亲水油藏的参数解释精度。
为实现上述目的,本申请提供如下技术方案:
本申请提供一种亲水油藏的参数解释方法,所述方法包括:
获得所述亲水油藏的第一参数集中各参数的参数值,所述第一参数集包括:油相的相对渗透率曲线、水相的相对渗透率曲线、参考油相地层体积系数、参考水相地层体积系数、油相重度、水相重度、压缩系数、厚度、深度、油相粘度、纯水粘度、以及静止时的聚合物水溶液粘度;
设定所述亲水油藏的第二参数集中各参数的参数值,所述第二参数集包括:渗透率下降系数、参考孔隙度、不可及孔隙体积、聚合物浓度初始分布、水饱和度初始分布、一价阳离子浓度初始分布、井底压力初始分布、绝对渗透率、水相粘浓曲线、吸附曲线、以及剪切变稀曲线;
利用预设的聚合物驱模型、所述亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,其中,所述聚合物驱模型包括水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程,并且,所述水相流动方程、所述油相流动方程、所述聚合物流动方程和所述一价阳离子流动方程中均包含渗透率下降系数;
根据求解出的所述油相压力、所述水相饱和度、所述聚合物浓度和所述一价阳离子浓度的参数值计算油产量和水产量;
根据所述油产量和所述水产量计算含水率;
对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合;
如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足所述预设精度要求,则重新设定所述第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果。
可选的,
所述油相流动方程为:
Figure BDA0002862606590000021
所述水相流动方程为:
Figure BDA0002862606590000022
所述聚合物流动方程为:
Figure BDA0002862606590000023
所述一价阳离子流动方程为:
Figure BDA0002862606590000024
其中,
Figure BDA0002862606590000025
为梯度符号,fro为油相流动能力的校正系数,frw为水相流动能力的校正系数,K为绝对渗透率,Kro为油相相对渗透率,Krw为水相相对渗透率,Rk为渗透率下降系数,μo为油相粘度,μw为流动时的聚合物水溶液粘度,Bo为油相地层体积系数,Bw为水相地层体积系数,po为油相压力,pw为水相压力,γo为油相重度,γw为水相重度,D为深度,φ为孔隙度,φp为不可及孔隙体积,So为油相饱和度,Sw为水相饱和度,qo为油相流量,qw为水相流量,qp为聚合物流量,qcl为一价阳离子流量,Cp为聚合物浓度,Ccl为一价阳离子浓度,Cpads为聚合物吸附浓度,t为时间。
可选的,所述对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合,包括:
对求解出的聚合物浓度的参数值与实测得到的聚合物浓度的参数值进行拟合,对求解出的油相压力的参数值与实测得到的油相压力的参数值进行拟合,对求解出的含水率的参数值与实测得到的含水率的参数值进行拟合。
本申请还提供一种亲水油藏的参数解释装置,包括:
参数值确定单元,用于获得所述亲水油藏的第一参数集中各参数的参数值,所述第一参数集包括:油相的相对渗透率曲线、水相的相对渗透率曲线、参考油相地层体积系数、参考水相地层体积系数、油相重度、水相重度、压缩系数、厚度、深度、油相粘度、纯水粘度、以及静止时的聚合物水溶液粘度;
参数值设定单元,用于设定所述亲水油藏的第二参数集中各参数的参数值,所述第二参数集包括:渗透率下降系数、参考孔隙度、不可及孔隙体积、聚合物浓度初始分布、水饱和度初始分布、一价阳离子浓度初始分布、井底压力初始分布、绝对渗透率、水相粘浓曲线、吸附曲线、以及剪切变稀曲线;
参数值求解单元,用于利用预设的聚合物驱模型、所述亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,其中,所述聚合物驱模型包括水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程,并且,所述水相流动方程、所述油相流动方程、所述聚合物流动方程和所述一价阳离子流动方程中均包含渗透率下降系数;
产量求解单元,用于根据求解出的所述油相压力、所述水相饱和度、所述聚合物浓度和所述一价阳离子浓度的参数值计算油产量和水产量;
含水率求解单元,用于根据所述油产量和所述水产量计算含水率;
拟合单元,用于对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合;
处理单元,用于如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足所述预设精度要求,则重新设定所述第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果。
可选的,所述拟合单元具体用于:
对求解出的聚合物浓度的参数值与实测得到的聚合物浓度的参数值进行拟合,对求解出的油相压力的参数值与实测得到的油相压力的参数值进行拟合,对求解出的含水率的参数值与实测得到的含水率的参数值进行拟合。
由此可见,本申请的有益效果为:
本申请公开的亲水油藏的参数解释方法和装置,首先获得亲水油藏的第一参数集中各参数的参数值,设定亲水油藏的第二参数集中各参数的参数值,之后利用预设的聚合物驱模型、亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,其中,聚合物驱模型中的水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程均包含渗透率下降系数,之后根据求解出的油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值计算油产量和水产量,再根据油产量和水产量计算含水率,之后对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合,如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足预设精度要求,则重新设定第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果。本申请公开的技术方案中,聚合物驱模型中的水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程均包含渗透率下降系数,这使得聚合物驱模型更加符合亲水油藏的流动规律,因此,能够提高亲水油藏的参数解释精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请公开的一种亲水油藏的参数解释方法的流程图;
图2为本申请公开的油相和水相的相对渗透率曲线的示意图;
图3为本申请公开的一种亲水油藏的参数解释装置的结构示意图。
具体实施方式
本申请提供一种亲水油藏的参数解释方法及装置,以提高亲水油藏的参数解释精度。
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
参见图1,图1为本申请公开的一种亲水油藏的参数解释方法的流程图。该方法包括:
S1:获得亲水油藏的第一参数集中各参数的参数值。
其中,第一参数集包括:油相的相对渗透率曲线、水相的相对渗透率曲线、参考油相地层体积系数、参考水相地层体积系数、油相重度、水相重度、压缩系数、厚度、深度、油相粘度、纯水粘度、以及静止时的聚合物水溶液粘度。
需要说明的是,第一参数集中各个参数的参数值是通过实测或者试验手段确定的,在参数解释过程中,第一参数集中各个参数的参数值是不调整的。
S2:设定亲水油藏的第二参数集中各参数的参数值。
其中,第二参数集包括:渗透率下降系数、参考孔隙度、不可及孔隙体积、聚合物浓度初始分布、水饱和度初始分布、一价阳离子浓度初始分布、井底压力初始分布、绝对渗透率、水相粘浓曲线、吸附曲线、以及剪切变稀曲线。需要说明的是,第二参数集中各个参数是待解释的参数。
S3:利用预设的聚合物驱模型、亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值。
其中,聚合物驱模型包括水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程。并且,水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程中均包含渗透率下降系数。
申请人经过大量研究发现:针对亲水油藏,注入的聚合物会使亲水油藏中的水相流体和油相流体的流动都变慢。因此,本申请所使用的聚合物驱流动模型,水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程中均包含渗透率下降系数,这使得本申请进行参数解释过程所使用的聚合物驱流动模型更加符合亲水油藏的流动规律。
S4:根据求解出油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值计算油产量和水产量。
S5:根据油产量和水产量计算含水率。
S6:对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合。
可选的,对求解出的聚合物浓度的参数值与实测得到的聚合物浓度的参数值进行拟合,对求解出的油相压力的参数值与实测得到的油相压力的参数值进行拟合,对求解出的含水率的参数值与实测得到的含水率的参数值进行拟合。
S7:如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足预设精度要求,则重新设定第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果。
本申请公开的亲水油藏的参数解释方法,首先获得亲水油藏的第一参数集中各参数的参数值,设定亲水油藏的第二参数集中各参数的参数值,之后利用预设的聚合物驱模型、亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,其中,聚合物驱模型中的水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程均包含渗透率下降系数,之后根据求解出的油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值计算油产量和水产量,再根据油产量和水产量计算含水率,之后对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合,如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足预设精度要求,则重新设定第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果。本申请公开的亲水油藏的参数解释方法,聚合物驱模型中的水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程均包含渗透率下降系数,这使得聚合物驱模型更加符合亲水油藏的流动规律,因此,能够提高亲水油藏的参数解释精度。
在一个实施例中,
油相流动方程为:
Figure BDA0002862606590000081
水相流动方程为:
Figure BDA0002862606590000082
聚合物流动方程为:
Figure BDA0002862606590000083
一价阳离子流动方程为:
Figure BDA0002862606590000084
其中:
Figure BDA0002862606590000085
为梯度符号。
fro为油相流动能力的校正系数,根据前一周期求解出的油相饱和度、聚合物浓度以及预先构建的油相的相对渗透率校正系数曲线确定。
frw为水相流动能力的校正系数,根据前一周期求解出的水相饱和度、聚合物浓度以及预先构建的水相的相对渗透率校正系数曲线确定。
也就是说,预先构建油相的相对渗透率校正系数曲线和水相的相对渗透率校正系数的曲线;在当前周期内,根据前一周期求解出的油相饱和度和水相饱和度,并结合油相的相对渗透率校正系数曲线和水相的相对渗透率校正系数的曲线,就可以确定当前周期的油相流动能力的校正系数fro和水相流动能力的校正系数frw
K为绝对渗透率,绝对渗透率K可能是一个值,也可能是一个分布。
Kro为油相相对渗透率,根据油相的相对渗透率曲线以及前一周期求解出的油相饱和度确定。
Krw为水相相对渗透率,根据水相的相对渗透率曲线以及前一周期求解出的水相饱和度确定。其中,水相饱和度和油相饱和度的和为1。
如图2所示,曲线L1为油相的相对渗透率曲线,曲线L2为水相的相对渗透率曲线,根据前一周期求解出的水相饱和度Sw以及水相的相对渗透率曲线就可以确定水相相对渗透率,根据前一周期求解出的水相饱和度Sw以及油相的相对渗透率曲线就可以确定油相的相对渗透率曲线。
Rk为渗透率下降系数。
μo为油相粘度。
μw为流动时的聚合物水溶液粘度,根据纯水粘度、静止时的聚合物水溶液粘度、当前设定的剪切变稀曲线、当前设定的水相粘浓曲线、以及前一周期求解出的一价阳离子浓度确定。
Bo为油相地层体积系数,是根据参考油相地层体积系数确定的。
Bw为水相地层体积系数,是根据参考水相地层体积系数确定的。
po为油相压力。
pw为水相压力。其中,油相压力po和水相压力pw是求解出的,在忽略毛管力的情况下,油相压力po和水相压力pw相等。
γo为油相重度。
γw为水相重度。
D为深度。
φ为孔隙度,根据压缩系数和当前设定的参考孔隙度确定。
φp为不可及孔隙体积,聚合物的分子较大,对于部分很小的孔隙,聚合物是无法进入的。
So为油相饱和度。
Sw为水相饱和度。其中,油相饱和度So和水相饱和度Sw是求解出的,且油相饱和度So和水相饱和度Sw的和为1。
qo为油相流量,根据油相粘度、参考油相地层体积系数、油相的相对渗透率曲线、厚度、深度、当前设定的聚合物浓度初始分布、当前设定的水饱和度初始分布、当前设定的井底压力初始分布、以及当前设定的绝对渗透率确定。
qw为水相流量,根据纯水粘度、静止时的聚合物水溶液粘度、参考水相地层体积系数、水相的相对渗透率曲线、厚度、深度、当前设定的聚合物浓度初始分布、当前设定的水饱和度初始分布、当前设定的一价阳离子浓度初始分布、当前设定的井底压力初始分布、当前设定的绝对渗透率、当前设定的水相粘浓曲线、当前设定的吸附曲线、以及当前设定的剪切变稀曲线确定。
qp为聚合物流量,根据纯水粘度、静止时的聚合物水溶液粘度、参考水相地层体积系数、水相的相对渗透率曲线、厚度、深度、当前设定的聚合物浓度初始分布、当前设定的水饱和度初始分布、当前设定的一价阳离子浓度初始分布、当前设定的井底压力初始分布、当前设定的绝对渗透率、当前设定的水相粘浓曲线、当前设定的吸附曲线、以及当前设定的剪切变稀曲线确定。
qcl为一价阳离子流量,根据纯水粘度、静止时的聚合物水溶液粘度、参考水相地层体积系数、水相的相对渗透率曲线、厚度、深度、当前设定的聚合物浓度初始分布、当前设定的水饱和度初始分布、当前设定的一价阳离子浓度初始分布、当前设定的井底压力初始分布、当前设定的绝对渗透率、当前设定的水相粘浓曲线、当前设定的吸附曲线、以及当前设定的剪切变稀曲线确定。
Cp为聚合物浓度,是求解出的。
Ccl为一价阳离子浓度,是求解出的。
Cpads为聚合物吸附浓度,根据当前设定的吸附曲线以及前一周期求解出的聚合物浓度确定。
t为时间。
本申请上述公开了亲水油藏的参数解释方法,相应的,本申请还公开亲水油藏的参数解释装置,说明书中关于两者的描述可以相互参考。
参见图3,图3为本申请公开的一种亲水油藏的参数解释装置的结构示意图。该装置包括:
参数值确定单元100,用于获得亲水油藏的第一参数集中各参数的参数值。
其中,第一参数集包括:油相的相对渗透率曲线、水相的相对渗透率曲线、参考油相地层体积系数、参考水相地层体积系数、油相重度、水相重度、压缩系数、厚度、深度、油相粘度、纯水粘度、以及静止时的聚合物水溶液粘度。
参数值设定单元200,用于设定亲水油藏的第二参数集中各参数的参数值。
其中,第二参数集包括:渗透率下降系数、参考孔隙度、不可及孔隙体积、聚合物浓度初始分布、水饱和度初始分布、一价阳离子浓度初始分布、井底压力初始分布、绝对渗透率、水相粘浓曲线、吸附曲线、以及剪切变稀曲线。
参数值求解单元300,用于利用预设的聚合物驱模型、亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值。
其中,聚合物驱模型包括水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程,并且,水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程中均包含渗透率下降系数。
产量求解单元400,用于根据求解出的油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值计算油产量和水产量。
含水率求解单元500,用于根据油产量和水产量计算含水率。
拟合单元600,用于对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合。
处理单元700,用于如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足预设精度要求,则重新设定第二参数集中至少一个参数的参数值,并触发参数值求解单元300求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,直至得到解释结果。
本申请公开的亲水油藏的参数解释装置,聚合物驱模型中的水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程均包含渗透率下降系数,这使得聚合物驱模型更加符合亲水油藏的流动规律,因此,能够提高亲水油藏的参数解释精度。
在一个实施例中,
油相流动方程为:
Figure BDA0002862606590000121
水相流动方程为:
Figure BDA0002862606590000122
聚合物流动方程为:
Figure BDA0002862606590000123
一价阳离子流动方程为:
Figure BDA0002862606590000124
在一个实施例中,拟合单元600具体用于:
对求解出的聚合物浓度的参数值与实测得到的聚合物浓度的参数值进行拟合,对求解出的油相压力的参数值与实测得到的油相压力的参数值进行拟合,对求解出的含水率的参数值与实测得到的含水率的参数值进行拟合。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的参数解释装置而言,由于其与实施例公开的参数解释方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (4)

1.一种亲水油藏的参数解释方法,其特征在于,所述方法包括:
获得所述亲水油藏的第一参数集中各参数的参数值,所述第一参数集包括:油相的相对渗透率曲线、水相的相对渗透率曲线、参考油相地层体积系数、参考水相地层体积系数、油相重度、水相重度、压缩系数、厚度、深度、油相粘度、纯水粘度、以及静止时的聚合物水溶液粘度;
设定所述亲水油藏的第二参数集中各参数的参数值,所述第二参数集包括:渗透率下降系数、参考孔隙度、不可及孔隙体积、聚合物浓度初始分布、水饱和度初始分布、一价阳离子浓度初始分布、井底压力初始分布、绝对渗透率、水相粘浓曲线、吸附曲线、以及剪切变稀曲线;
利用预设的聚合物驱模型、所述亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,其中,所述聚合物驱模型包括水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程,并且,所述水相流动方程、所述油相流动方程、所述聚合物流动方程和所述一价阳离子流动方程中均包含渗透率下降系数;
根据求解出的所述油相压力、所述水相饱和度、所述聚合物浓度和所述一价阳离子浓度的参数值计算油产量和水产量;
根据所述油产量和所述水产量计算含水率;
对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合;
如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足所述预设精度要求,则重新设定所述第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果;
所述油相流动方程为:
Figure FDA0003710547140000021
所述水相流动方程为:
Figure FDA0003710547140000022
所述聚合物流动方程为:
Figure FDA0003710547140000023
所述一价阳离子流动方程为:
Figure FDA0003710547140000024
其中,
Figure FDA0003710547140000025
为梯度符号,fro为油相流动能力的校正系数,frw为水相流动能力的校正系数,K为绝对渗透率,Kro为油相相对渗透率,Krw为水相相对渗透率,Rk为渗透率下降系数,μo为油相粘度,μw为流动时的聚合物水溶液粘度,Bo为油相地层体积系数,Bw为水相地层体积系数,po为油相压力,pw为水相压力,γo为油相重度,γw为水相重度,D为深度,φ为孔隙度,φp为不可及孔隙体积,So为油相饱和度,Sw为水相饱和度,qo为油相流量,qw为水相流量,qp为聚合物流量,qcl为一价阳离子流量,Cp为聚合物浓度,Ccl为一价阳离子浓度,Cpads为聚合物吸附浓度,t为时间。
2.根据权利要求1所述的方法,其特征在于,所述对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合,包括:
对求解出的聚合物浓度的参数值与实测得到的聚合物浓度的参数值进行拟合,对求解出的油相压力的参数值与实测得到的油相压力的参数值进行拟合,对求解出的含水率的参数值与实测得到的含水率的参数值进行拟合。
3.一种亲水油藏的参数解释装置,其特征在于,包括:
参数值确定单元,用于获得所述亲水油藏的第一参数集中各参数的参数值,所述第一参数集包括:油相的相对渗透率曲线、水相的相对渗透率曲线、参考油相地层体积系数、参考水相地层体积系数、油相重度、水相重度、压缩系数、厚度、深度、油相粘度、纯水粘度、以及静止时的聚合物水溶液粘度;
参数值设定单元,用于设定所述亲水油藏的第二参数集中各参数的参数值,所述第二参数集包括:渗透率下降系数、参考孔隙度、不可及孔隙体积、聚合物浓度初始分布、水饱和度初始分布、一价阳离子浓度初始分布、井底压力初始分布、绝对渗透率、水相粘浓曲线、吸附曲线、以及剪切变稀曲线;
参数值求解单元,用于利用预设的聚合物驱模型、所述亲水油藏的第一参数集和第二参数集中各参数的参数值,求解油相压力、水相饱和度、聚合物浓度和一价阳离子浓度的参数值,其中,所述聚合物驱模型包括水相流动方程、油相流动方程、聚合物流动方程和一价阳离子流动方程,并且,所述水相流动方程、所述油相流动方程、所述聚合物流动方程和所述一价阳离子流动方程中均包含渗透率下降系数;
产量求解单元,用于根据求解出的所述油相压力、所述水相饱和度、所述聚合物浓度和所述一价阳离子浓度的参数值计算油产量和水产量;
含水率求解单元,用于根据所述油产量和所述水产量计算含水率;
拟合单元,用于对求解出的聚合物浓度、油相压力和含水率中至少一个参数的参数值与对应参数的实测的参数值进行拟合;
处理单元,用于如果拟合结果满足预设精度要求,则将当前设定的第二参数集中各参数的参数值作为解释结果,如果拟合结果不满足所述预设精度要求,则重新设定所述第二参数集中至少一个参数的参数值,并执行后续步骤,直至得到解释结果;
所述油相流动方程为:
Figure FDA0003710547140000031
所述水相流动方程为:
Figure FDA0003710547140000032
所述聚合物流动方程为:
Figure FDA0003710547140000033
所述一价阳离子流动方程为:
Figure FDA0003710547140000041
其中,
Figure FDA0003710547140000042
为梯度符号,fro为油相流动能力的校正系数,frw为水相流动能力的校正系数,K为绝对渗透率,Kro为油相相对渗透率,Krw为水相相对渗透率,Rk为渗透率下降系数,μo为油相粘度,μw为流动时的聚合物水溶液粘度,Bo为油相地层体积系数,Bw为水相地层体积系数,po为油相压力,pw为水相压力,γo为油相重度,γw为水相重度,D为深度,φ为孔隙度,φp为不可及孔隙体积,So为油相饱和度,Sw为水相饱和度,qo为油相流量,qw为水相流量,qp为聚合物流量,qcl为一价阳离子流量,Cp为聚合物浓度,Ccl为一价阳离子浓度,Cpads为聚合物吸附浓度,t为时间。
4.根据权利要求3所述的装置,其特征在于,所述拟合单元具体用于:
对求解出的聚合物浓度的参数值与实测得到的聚合物浓度的参数值进行拟合,对求解出的油相压力的参数值与实测得到的油相压力的参数值进行拟合,对求解出的含水率的参数值与实测得到的含水率的参数值进行拟合。
CN202011570188.5A 2020-12-26 2020-12-26 亲水油藏的参数解释方法及装置 Active CN112685971B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011570188.5A CN112685971B (zh) 2020-12-26 2020-12-26 亲水油藏的参数解释方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011570188.5A CN112685971B (zh) 2020-12-26 2020-12-26 亲水油藏的参数解释方法及装置

Publications (2)

Publication Number Publication Date
CN112685971A CN112685971A (zh) 2021-04-20
CN112685971B true CN112685971B (zh) 2022-12-02

Family

ID=75452002

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011570188.5A Active CN112685971B (zh) 2020-12-26 2020-12-26 亲水油藏的参数解释方法及装置

Country Status (1)

Country Link
CN (1) CN112685971B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096300A (zh) * 2016-06-20 2016-11-09 合肥工业大学 低渗透油气藏的参数解释方法及解释系统
CN106761621A (zh) * 2017-02-08 2017-05-31 中海石油(中国)有限公司 一种聚合物驱三层窜流油藏试井解释参数的获取方法
CN107066679A (zh) * 2017-03-09 2017-08-18 中海石油(中国)有限公司 一种用于聚合物驱双层窜流油藏试井分析系统及方法
CN109441415A (zh) * 2018-12-19 2019-03-08 克拉玛依新科澳石油天然气技术股份有限公司 基于邻井干扰的聚合物驱油藏测试井的试井解释方法
CN110578503A (zh) * 2019-09-20 2019-12-17 合肥鼎盛锦业科技有限公司 聚合物驱数值试井中确定油藏初始水饱和度的方法及系统
US10815778B1 (en) * 2019-12-02 2020-10-27 China University Of Petroleum (East China) Method for obtaining formation parameters of gas hydrate reservoir through well testing interpretation

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8510089B2 (en) * 2010-08-31 2013-08-13 Chevron U.S.A., Inc. Computer-implemented systems and methods for forecasting performance of polymer flooding of an oil reservoir system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096300A (zh) * 2016-06-20 2016-11-09 合肥工业大学 低渗透油气藏的参数解释方法及解释系统
CN106761621A (zh) * 2017-02-08 2017-05-31 中海石油(中国)有限公司 一种聚合物驱三层窜流油藏试井解释参数的获取方法
CN107066679A (zh) * 2017-03-09 2017-08-18 中海石油(中国)有限公司 一种用于聚合物驱双层窜流油藏试井分析系统及方法
CN109441415A (zh) * 2018-12-19 2019-03-08 克拉玛依新科澳石油天然气技术股份有限公司 基于邻井干扰的聚合物驱油藏测试井的试井解释方法
CN110578503A (zh) * 2019-09-20 2019-12-17 合肥鼎盛锦业科技有限公司 聚合物驱数值试井中确定油藏初始水饱和度的方法及系统
US10815778B1 (en) * 2019-12-02 2020-10-27 China University Of Petroleum (East China) Method for obtaining formation parameters of gas hydrate reservoir through well testing interpretation

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Marching cube"s modification method in reservoir simulation model;Wang Haiyang等;《Computer Engineering and Applications》;20070228;第43卷(第4期);全文 *
基于卷积神经网络的径向复合油藏自动试井解释方法;李道伦等;《石油勘探与开发》;20200630;第47卷(第3期);全文 *
碱-聚合物复合驱油藏流线数值试井解释模型及其应用;姚军等;《石油学报》;20081115(第06期);全文 *
聚合物驱油藏数值试井解释方法研究;张德志等;《石油钻探技术》;20060625(第03期);全文 *

Also Published As

Publication number Publication date
CN112685971A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
EP3391025A1 (en) Method and device for determining gas permeability of a subsurface formation
CN105332698B (zh) 测定残余油饱和度的方法和油水相对渗透率关系的方法
CN105223116A (zh) 一种基于核磁共振谱系数法计算束缚水饱和度的方法
CN106204304B (zh) 一种砾岩油藏聚驱相对渗透率曲线的确定方法
CN107525746A (zh) 一种表征致密油藏低速非达西渗流特征的方法及装置
CN109932297A (zh) 一种致密砂岩储层渗透率的计算方法
CN111625917B (zh) 确定多孔介质相对渗透率曲线的方法及装置
CN112685971B (zh) 亲水油藏的参数解释方法及装置
Lin et al. Structural stability for a Brinkman fluid
NO321567B1 (no) Fremgangsmate og innretning for rask maling av resistivitetsindeksen til faste prover, slik som stein
CN110309611B (zh) 基于气水厚度分布的气水两相渗流规律预测方法及系统
CN110578503A (zh) 聚合物驱数值试井中确定油藏初始水饱和度的方法及系统
CN105844011B (zh) 一种基于毛管模型的渗透率计算方法
CN108489878B (zh) 一种基于数值模拟迭代消除末端效应的相渗曲线校正方法
Ren et al. Correlation of colloid collision efficiency with hydraulic conductivity of silica sands
Hagoort et al. Numerical solution of the material balance equations of compartmented gas reservoirs
CA2490264A1 (en) Liquid extrusion porosimeter and method
CN105823706B (zh) 气润湿性测试方法
CN106285658B (zh) 一种毛管力曲线校正的方法及装置
CN110857625B (zh) 碳酸盐岩有封闭水体缝洞单元储量获取方法、装置及设备
CN110674463A (zh) 一种确定油水两相相对渗透率曲线的方法和装置
Templeton et al. A study of gravity counterflow segregation
CN115788399A (zh) 高倍数水驱岩心油水相对渗透率曲线测量装置及方法
CN108717498A (zh) 一种油藏相渗曲线模型及相渗曲线计算方法
CN213337216U (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