CN110287545A - 一种计算油藏内部任意部位原始含水饱和度的方法及系统 - Google Patents

一种计算油藏内部任意部位原始含水饱和度的方法及系统 Download PDF

Info

Publication number
CN110287545A
CN110287545A CN201910474989.2A CN201910474989A CN110287545A CN 110287545 A CN110287545 A CN 110287545A CN 201910474989 A CN201910474989 A CN 201910474989A CN 110287545 A CN110287545 A CN 110287545A
Authority
CN
China
Prior art keywords
oil reservoir
original water
parameter
function
water saturation
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
Application number
CN201910474989.2A
Other languages
English (en)
Other versions
CN110287545B (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.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
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 Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201910474989.2A priority Critical patent/CN110287545B/zh
Publication of CN110287545A publication Critical patent/CN110287545A/zh
Application granted granted Critical
Publication of CN110287545B publication Critical patent/CN110287545B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种计算油藏内部任意部位原始含水饱和度的方法及系统，所述方法包括：获取油藏参数；基于所述油藏参数根据油藏原始含水饱和度模型得到所述油藏原始含水饱和度，本发明可得到油藏内部任意部位的原始含水饱和度。

Description

一种计算油藏内部任意部位原始含水饱和度的方法及系统
技术领域
本发明涉及石油勘探与开发技术领域，尤其涉及一种计算油藏内部任意部位原始含水饱和度的方法及系统。
背景技术
目前，通常采用J函数法建立油藏原始含水饱和度模型，从而计算油藏原始含水饱和度分布。当油藏储层非均质性较强时，可以根据流动单元指数(FZI)将油藏划分为若干个流动单元。对于各个流动单元，通过拟合各个流动单元的J函数，可建立每个流动单元的原始含水饱和度模型，从而得到油藏原始含水饱和度的分布情况。通常做法中，J(Sw)＝a(Sw)b，对于每个块特定的流动单元，a和b都是定值。但在实际作业过程中，取心和J函数曲线的数量是有限的，因此，无法获取油藏内部任意部位的J函数曲线，从而无法获取油藏内部任意部位的拟合的J函数的a和b，进而无法精确建立油藏内部任意部位的原始含水饱和度模型。因次，对于没有J函数曲线的岩心或油藏内部任意部位，无法精确计算其原始含水饱和度分布。
发明内容
本发明的一个目的在于提供一种计算油藏内部任意部位原始含水饱和度的方法，以得到油藏内部任意部位的原始含水饱和度。本发明的另一个目的在于提供一种计算油藏内部任意部位原始含水饱和度的系统。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。常规方法基于有限的压汞曲线来获得J函数曲线；拟合所述J函数曲线，获得拟合的J函数，进而获得拟合的J函数的参数a和b。只能计算有J函数曲线的流动单元的原始含水饱和度，不能精确建立油藏内部任意部位的原始含水饱和度模型，无法精确计算油藏内部任意部位的原始含水饱和度分布。本发明与常规方法的最大区别在于，通过建立拟合的J函数的参数a与流动单元指数的对应关系和拟合的J函数的参数b与流动单元指数的对应关系，解决了上述问题。
为了达到以上目的，本发明一方面公开了一种计算油藏内部任意部位原始含水饱和度的方法，包括：
获取油藏参数；
基于所述油藏参数根据油藏原始含水饱和度模型得到所述油藏内部任意部位的原始含水饱和度。
优选地，所述方法进一步包括形成所述油藏原始含水饱和度模型的步骤。
优选地，所述形成所述油藏原始含水饱和度模型的步骤具体包括：
确定不同岩心对应的拟合的J函数的常量参数；
形成所述常量参数与岩心的流动单元指数的对应关系；
根据所述对应关系得到所述油藏内部任意部位的原始含水饱和度模型。
优选地，所述确定不同岩心对应的拟合的J函数的常量参数具体包括：
获取每块岩心对应的J函数曲线；
拟合每条J函数曲线，得到拟合的J函数；
获取所述拟合的J函数的常量参数。
优选地，所述常量参数包括第一参数和第二参数，所述形成所述常量参数与流动单元指数的对应关系具体包括：
分别形成表示流动单元指数与所述第一参数的对应关系的第一对应关系和表示流动单元指数与所述第二参数的对应关系的第二对应关系。
优选地，所述根据所述对应关系得到所述油藏内部任意部位原始含水饱和度模型具体包括：
根据J函数的定义、拟合的J函数与原始含水饱和度的关系、地层条件下与实验室条件下毛管力的关系以及毛管力与油水密度差的关系得到原始含水饱和度与拟合的J函数的常量参数的对应关系；
根据流动单元指数与拟合的J函数的常量参数的对应关系以及原始含水饱和度与拟合的J函数的常量参数的对应关系得到原始含水饱和度与流动单元指数的对应关系，所述原始含水饱和度与流动单元指数的对应关系作为所述油藏原始含水饱和度模型。
本发明另一方面还公开了一种计算油藏内部任意部位原始含水饱和度的系统，包括参数获取单元和原始含水饱和度求解单元；
所述参数获取单元用于获取油藏参数；
所述含水饱和度求解单元用于基于所述油藏参数根据油藏原始含水饱和度模型得到所述油藏内部任意部位的原始含水饱和度。
优选地，所述系统进一步包括模型建立单元；
所述模型建立单元用于形成所述油藏原始含水饱和度模型。
优选地，所述模型建立单元进一步用于确定不同岩心对应的J函数的常量参数，形成所述常量参数与岩心的流动单元指数的对应关系，并根据所述对应关系得到所述油藏原始含水饱和度模型。
优选地，所述模型建立单元进一步用于获取每块岩心对应的J函数曲线，通过拟合，获得拟合的J函数，进而获得拟合的J函数的常量参数。
优选地，所述常量参数包括第一参数和第二参数，所述模型建立单元进一步用于分别形成表示流动单元指数与所述第一参数的对应关系的第一对应关系和表示流动单元指数与所述第二参数的对应关系的第二对应关系。
优选地，所述模型建立单元进一步用于根据J函数的定义、拟合的J函数与原始含水饱和度的关系、地层条件下与实验室条件下毛管力的关系以及毛管力与油水密度差的关系得到原始含水饱和度与拟合的J函数的常量参数的对应关系，根据流动单元指数与拟合的J函数的常量参数的对应关系以及原始含水饱和度与拟合的J函数的常量参数的对应关系得到原始含水饱和度与流动单元指数的对应关系，所述原始含水饱和度与流动单元指数的对应关系作为所述油藏原始含水饱和度模型。
本发明还公开了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质，其上存储有计算机程序，该程序被处理器执行时实现如上所述方法。
本发明通过直接获取地层的油藏参数，将油藏参数输入油藏原始含水饱和度模型中即可得到油藏内部任意部位的原始含水饱和度。本发明油藏原始含水饱和度模型需要的输入参数均为可直接从已有油藏地质模型任意部位获取的油藏参数，则通过获取油藏内部任意部位的油藏参数即可得到该部位的含水饱和度，从而可在有限取心的情况下全面了解油藏内部任意部位的含水饱和度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体实施例的流程图之一；
图2示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体实施例的流程图之二；
图3示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体实施例的流程图之三；
图4示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体实施例的流程图之四；
图5示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体实施例的流程图之五；
图6示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体实施例的流程图之六；
图7示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体算例中J函数的曲线图；
图8示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体算例中第一对应关系的关系曲线图；
图9示出本发明一种计算油藏内部任意部位原始含水饱和度的方法具体算例中第二对应关系的关系曲线图；
图10示出本发明一种计算油藏内部任意部位原始含水饱和度的系统具体实施例的结构图之一；
图11示出本发明一种计算油藏内部任意部位原始含水饱和度的系统具体实施例的结构图之二；
图12示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
根据本发明的一个方面，本实施例公开了一种计算油藏内部任意部位原始含水饱和度的方法。如图1所示，本实施例中，所述方法包括：
S100：获取油藏参数。
S200：基于所述油藏参数根据油藏原始含水饱和度模型得到所述油藏原始含水饱和度。在一个具体例子中，油藏参数可包括束缚水饱和度、水油密度差、油水界面以上高度、流动单元指数、油藏条件下油水界面张力、油藏条件下油水接触角、渗透率和孔隙度。在得到油藏原始含水饱和度模型后，通过获取油藏内部任意部位的油藏参数并输入到油藏原始含水饱和度模型中即可得到该部位的含水饱和度。
本发明通过直接获取地层的油藏参数，将油藏参数输入油藏原始含水饱和度模型中即可得到油藏内部任意部位的含水饱和度。本发明油藏原始含水饱和度模型需要的输入参数均为可直接从已有油藏地质模型任意部位获取的油藏参数，则通过获取油藏内部任意部位的油藏参数即可得到该部位的含水饱和度，从而可在有限取心的情况下全面了解油藏内部任意部位的含水饱和度。
在优选的实施方式中，如图2所示，所述方法进一步可包括形成所述油藏原始含水饱和度模型的步骤S000。在实际作业中，取心的数量有限，通过有限的取心可得到J函数曲线；根据J函数曲线可得到这些岩心的J函数曲线；拟合所述J函数曲线，得到拟合的J函数，进而得到拟合的J函数的常量参数；通过得到的拟合的J函数的常量参数可得到这些岩心的含水饱和度。拟合的J函数的常量参数与岩心的流动单元指数存在对应关系，该对应关系可类推至其油藏内部任意部位，从而通过已有的拟合的J函数可得到能够求解油藏内部任意部位的原始含水饱和度模型，以实现通过简单的油藏参数获取得到油藏内部任意部位的原始含水饱和度的目的。
在优选的实施方式中，如图3所示，所述S000具体可包括：
S010：确定不同岩心对应的拟合的J函数的常量参数。拟合的J函数通常为幂函数，包括两个常量参数。通过压汞实验得到岩心的压汞曲线后，处理压汞曲线得到J函数曲线，拟合该J函数曲线得到一个幂函数，该幂函数即为拟合的J函数，该幂函数的两个常量参数即为拟合的J函数的常量参数。
S020：形成所述常量参数与岩心的流动单元指数的对应关系。不同岩心会得到不同的拟合的J函数，由此，岩心与拟合的J函数一一对应，拟合的J函数的常量参数与对应的岩心的流动单元指数存在对应关系，通过拟合多个岩心的常量参数和流动单元指数之间的关系，可形成对于油藏内部任意部位均适用的常量参数与流动单元指数的对应关系。
S030：根据所述对应关系得到所述油藏原始含水饱和度模型。通过拟合的J函数的常量参数与流动单元指数的对应关系可将现有的油藏原始含水饱和度模型中的常量参数用流动单元指数表示，进而得到与拟合的J函数的常量参数无关的油藏原始含水饱和度模型，通过油藏参数可直接得到原始含水饱和度，大大简化了油藏原始含水饱和度的计算过程，同时可以精确确定油藏内部任意部位的原始含水饱和度。
在一个具体例子中，所述油藏原始含水饱和度模型可通过以下数学表达式表示：
其中，Sw为油水界面以上任意部位对应的原始含水饱和度；Swi为束缚水饱和度；Δρ为水油密度差，g/cm3；H为油水界面以上任意高度，m；FZI为流动单元指数；σres为油藏条件下油水界面张力，mN/m；θres为油藏条件下油水接触角度；K为渗透率，mD；为孔隙度。
在优选的实施方式中，如图4所示，所述S010具体可包括：
S011：获取每块岩心对应的J函数曲线。
S012：拟合J函数曲线，得到幂形式的拟合的J函数。其中，J函数曲线的拟合可通过例如Excel等的多种方法实现。当采用Excel拟合时，可通过Excel表示拟合的J函数，根据拟合的J函数与原始含水饱和度的关系，J＝aSwn b，最后得到拟合的J函数曲线的纵坐标为J，横坐标为标准化饱和度Swn，同时可确定a和b两个常量参数的值。
S013：获取所述拟合的J函数的常量参数。即当J＝aSwn b时，从拟合的J函数中提取a和b两个常量参数。
在优选的实施方式中，如图5所示，所述拟合的J函数为幂函数，所述常量参数包括第一参数和第二参数。所述S020具体可包括：
S021：分别形成表示流动单元指数与所述第一参数的对应关系的第一对应关系和表示流动单元指数与所述第二参数的对应关系的第二对应关系。
具体的，拟合的J函数为＝aSwn b时，则常量参数包括第一参数a和第二参数b，对于两个参数，需分别形成流动单元指数与每个参数的对应关系。形成表示流动单元指数与所述第一参数的对应关系的第一对应关系后，第一参数可仅通过流动单元指数对应表示，表示流动单元指数与所述第二参数的对应关系的第二对应关系后，第二参数可仅通过流动单元指数对应表示。将现有的油藏原始含水饱和度计算模型中的第一参数和第二参数由流动单元指数表示，得到本实施例的油藏原始含水饱和度模型。其中，建立第一参数a与FZI之间的函数关系式及第二参数b与FZI之间的函数关系式的方法有很多，例如，可以通过现有的Matlab、SPSS软件来实现，也可以使用Excel自带的拟合功能来实现。
在优选的实施方式中，如图6所示，所述S030具体可包括：
S031：根据J函数的定义、拟合的J函数与原始含水饱和度的关系、地层条件下与实验室条件下毛管力的关系以及毛管力与油水密度差的关系得到原始含水饱和度与J函数的常量参数的对应关系。
S032：根据流动单元指数与拟合的J函数的常量参数的对应关系以及原始含水饱和度与拟合的J函数的常量参数的对应关系得到原始含水饱和度与流动单元指数的对应关系，所述原始含水饱和度与流动单元指数的对应关系作为所述油藏原始含水饱和度模型。
具体的，在一个具体例子中，J函数的定义为：
其中，J称为J函数，无因次；Pclab为实验室条件下毛管力，psi；σlab为实验室条件下油水界面张力，mN/m；θlab为实验室条件下油水接触角度。
拟合的J函数与原始含水饱和度的关系为：
J＝aSwn b
其中，
Swn为标准化饱和度。
地层条件下与实验室条件下毛管力转化公式为：
其中，Pcres为油藏条件下毛管力，psi。
毛管力Pcres与油水密度差Δρ之间存在如下关系：
Pcres＝1.4214ΔρH
在一个具体算例中，选择一条J函数曲线，按照J＝aSwn b的形式拟合得到拟合的J函数，如图7所示。获取J函数的常量参数底数a和指数b分别为：a＝0.3771，b＝-0.937。同理，可得到其他J含函数曲线对应的底数a和指数b。进一步地，确定每对a和b对应的FZI，如表1所示。
表1
FZI a b
0.126 0.3771 -0.937
0.103 0.4464 -0.897
1.955 0.1234 -1.744
0.657 0.2081 -1.412
0.675 0.2169 -1.312
0.412 0.1815 -1.415
0.25 0.4682 -1.137
0.23 0.1965 -1.282
0.177 0.5264 -1.852
0.183 0.2378 -1.182
1.747 0.1357 -2.030
采用Excel自带的拟合功能以a为纵坐标、以FZI为横坐标将一一对应的a和FZI拟合得到第一对应关系的关系曲线，如图8所示，表达式为：a＝0.1663(FZI)-0.414；以b为纵坐标、以FZI为横坐标将一一对应的b和FZI拟合得到第二对应关系的关系曲线，如图9所示，表达式为：b＝-0.332ln(FZI)-1.6255。将得到的两个表达式代入现有的油藏原始含水饱和度模型中得到本发明的可求解油藏内部任意部位原始含水饱和度的油藏原始含水饱和度模型。对于油藏内部任意部位，输入油藏参数至本发明的油藏原始含水饱和度模型即可得到所述任意部位原始含水饱和度。
基于相同原理，本实施例还公开了一种计算油藏内部任意部位原始含水饱和度的系统。如图10所示，所述系统包括参数获取单元11和原始含水饱和度求解单元12。
其中，所述参数获取单元11用于获取油藏参数。
所述含水饱和度求解单元12用于基于所述油藏参数根据油藏原始含水饱和度模型得到所述油藏原始含水饱和度。
在一个具体例子中，油藏参数可包括束缚水饱和度、水油密度差、油水界面以上高度、流动单元指数、油藏条件下油水界面张力、油藏条件下油水接触角、渗透率和孔隙度。在得到油藏原始含水饱和度模型后，通过获取油藏内部任意部位的油藏参数并输入到油藏原始含水饱和度模型中即可得到该岩心的含水饱和度。
本发明通过直接获取地层的油藏参数，将油藏参数输入油藏原始含水饱和度模型中即可得到油藏内部任意部位的含水饱和度。本发明油藏原始含水饱和度模型需要的输入参数均为可直接从已有油藏地质模型的任意部位获取的油藏参数，则通过获取油藏内部任意部位的油藏参数即可得到该部位的含水饱和度，从而可在有限取心的情况下全面了解油藏内部任意部位的含水饱和度。
在优选的实施方式中，如图11所示，所述系统进一步包括模型建立单元10。所述模型建立单元10用于形成所述油藏原始含水饱和度模型。在实际作业中，取心的数量有限，通过有限的取心可得到压汞曲线；根据压汞曲线可得到这些岩心的J函数曲线；拟合所述J函数曲线，得到拟合的J函数，进而得到拟合的J函数的常量参数；通过得到的拟合的J函数的常量参数，可得到这些岩心的含水饱和度。拟合的J函数的常量参数与岩心的流动单元指数存在对应关系，该对应关系可类推至其油藏内部任意部位，从而通过已有的拟合的J函数可得到能够求解油藏内部任意部位的原始含水饱和度模型，以实现通过简单的油藏参数获取得到油藏内部任意部位的原始含水饱和度的目的。
在优选的实施方式中，所述模型建立单元10进一步用于确定不同岩心对应的拟合的J函数的常量参数，形成所述常量参数与流动单元指数的对应关系，并根据所述对应关系得到所述油藏原始含水饱和度模型。
拟合的J函数为幂函数，包括两个常量参数。通过压汞实验得到岩心的压汞曲线后，由压汞曲线，得到J函数曲线，拟合所述J函数曲线得到一个幂函数，该幂函数即为拟合的J函数，该幂函数的两个常量参数即为拟合的J函数的常量参数。不同岩心会得到不同的拟合的J函数，不同岩心具有不同流动单元指数，由此，流动单元指数与拟合的J函数一一对应，拟合的J函数的常量参数与对应的岩心的流动单元指数存在对应关系，通过拟合多块岩心的常量参数和流动单元指数之间的关系，可形成对于油藏内部任意部位均适用的常量参数与流动单元指数的对应关系。通过拟合的J函数的常量参数与流动单元指数的对应关系可将现有的油藏原始含水饱和度模型中的常量参数用流动单元指数表示，进而得到与拟合的J函数的常量参数无关的油藏原始含水饱和度模型，通过油藏内部任意部位的油藏参数可直接得到原始含水饱和度，大大简化了油藏原始含水饱和度的计算过程，同时可以精确确定藏内部任意部位的原始含水饱和度。
在一个具体例子中，所述油藏原始含水饱和度模型可通过以下数学表达式表示：
其中，Sw为油水界面以上任意部位对应的原始含水饱和度；Swi为束缚水饱和度；Δρ为水油密度差，g/cm3；H为油水界面以上任意高度，m；FZI为流动单元指数；σres为油藏条件下油水界面张力，mN/m；θres为油藏条件下油水接触角度；K为渗透率，mD；为孔隙度。
在优选的实施方式中，所述模型建立单元10进一步用于获取每块岩心对应的J函数曲线，拟合每条J函数曲线得到拟合的J函数，并获取所述拟合的J函数的常量参数。
其中，J函数的拟合可通过例如Excel等的多种方法实现。当采用Excel拟合时，可通过Excel表示得到的拟合的J函数，根据拟合的J函数与原始含水饱和度的关系，J＝aSwn b，最后得到拟合的J函数曲线的纵坐标为J，横坐标为标准化饱和度Swn，同时可确定a和b两个常量参数的值。即当J＝aSwn b时，从拟合得到的J函数中提取a和b两个常量参数。
在优选的实施方式中，所述常量参数包括第一参数和第二参数，所述模型建立单元10进一步用于分别形成表示流动单元指数与所述第一参数的对应关系的第一对应关系和表示流动单元指数与所述第二参数的对应关系的第二对应关系。具体的，拟合的J函数为＝aSwn b，则常量参数包括第一参数a和第二参数b，对于两个参数，需分别形成流动单元指数与每个参数的对应关系。形成表示流动单元指数与所述第一参数的对应关系的第一对应关系后，第一参数可仅通过流动单元指数对应表示，表示流动单元指数与所述第二参数的对应关系的第二对应关系后，第二参数可仅通过流动单元指数对应表示，将现有的油藏原始含水饱和度计算模型中的第一参数和第二参数由流动单元指数表示，得到本实施例的油藏原始含水饱和度模型。其中，建立第一参数a与FZI之间的函数关系式及第二参数b与FZI之间的函数关系式的方法有很多，例如，可以通过现有的Matlab、SPSS软件来实现，也可以使用Excel自带的拟合功能来实现。
在优选的实施方式中，所述模型建立单元10进一步用于根据J函数的定义、拟合的J函数与原始含水饱和度的关系、地层条件下与实验室条件下毛管力的关系以及毛管力与油水密度差的关系得到原始含水饱和度与J函数的常量参数的对应关系，根据流动单元指数与拟合的J函数的常量参数的对应关系以及原始含水饱和度与拟合的J函数的常量参数的对应关系得到原始含水饱和度与流动单元指数的对应关系，所述原始含水饱和度与流动单元指数的对应关系作为所述油藏原始含水饱和度模型。
具体的，在一个具体例子中，J函数的定义为：
其中，J称为J函数，无因次；Pclab为实验室条件下毛管力，psi；σlab为实验室条件下油水界面张力，mN/m；θlab为实验室条件下油水接触角度。
拟合的J函数与原始含水饱和度的关系为：
J＝aSwn b
其中，
Swn为标准化饱和度。
地层条件下与实验室条件下毛管力转化公式为：
其中，Pcres为油藏条件下毛管力，psi。
毛管力Pcres与油水密度差Δρ之间存在如下关系：
Pcres＝1.4214ΔρH
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备，具体的，计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述的由客户端执行的方法，或者，所述处理器执行所述程序时实现如上所述的由服务器执行的方法。
下面参考图12，其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图12所示，计算机设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口606。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或流程和/或方框图一个方框或方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或流程和/或方框图一个方框或方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或流程和/或方框图一个方框或方框中指定的功能的步骤。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (14)

1.一种计算油藏内部任意部位原始含水饱和度的方法，其特征在于，包括：
获取油藏参数；
基于所述油藏参数根据预设的油藏原始含水饱和度模型得到所述油藏原始含水饱和度。
2.根据权利要求1所述的计算油藏内部任意部位原始含水饱和度的方法，其特征在于，所述方法进一步包括形成所述油藏原始含水饱和度模型的步骤。
3.根据权利要求2所述的计算油藏内部任意部位原始含水饱和度的方法，其特征在于，所述形成所述油藏原始含水饱和度模型的步骤具体包括：
确定不同岩心对应的J函数的常量参数；
形成所述常量参数与流动单元指数的对应关系；
根据所述对应关系得到所述油藏原始含水饱和度模型。
4.根据权利要求3所述的计算油藏内部任意部位原始含水饱和度的方法，其特征在于，所述确定不同岩心对应的J函数的常量参数具体包括：
获取每块岩心对应的J函数曲线；
拟合所述J函数曲线，得到拟合的J函数；
获取所述拟合的J函数的常量参数。
5.根据权利要求3所述的计算油藏内部任意部位原始含水饱和度的方法，其特征在于，所述常量参数包括第一参数和第二参数，所述形成所述常量参数与岩心的流动单元指数的对应关系具体包括：
分别形成表示流动单元指数与所述第一参数的对应关系的第一对应关系和表示流动单元指数与所述第二参数的对应关系的第二对应关系。
6.根据权利要求3所述的计算油藏内部任意部位原始含水饱和度的方法，其特征在于，所述根据所述对应关系得到所述油藏原始含水饱和度模型具体包括：
根据J函数的定义、拟合的J函数与原始含水饱和度的关系、地层条件下与实验室条件下毛管力的关系以及毛管力与油水密度差的关系得到原始含水饱和度与拟合的J函数的常量参数的对应关系；
根据流动单元指数与拟合的J函数的常量参数的对应关系以及原始含水饱和度与拟合的J函数的常量参数的对应关系得到原始含水饱和度与流动单元指数的对应关系，所述原始含水饱和度与流动单元指数的对应关系作为所述油藏原始含水饱和度模型。
7.一种计算油藏内部任意部位原始含水饱和度的系统，其特征在于，包括参数获取单元和含水饱和度求解单元；
所述参数获取单元用于获取油藏参数；
所述含水饱和度求解单元用于基于所述油藏参数根据油藏原始含水饱和度模型得到所述油藏原始含水饱和度。
8.根据权利要求7所述的计算油藏内部任意部位原始含水饱和度的系统，其特征在于，所述系统进一步包括模型建立单元；
所述模型建立单元用于形成所述油藏原始含水饱和度模型。
9.根据权利要求8所述的计算油藏内部任意部位原始含水饱和度的系统，其特征在于，所述模型建立单元进一步用于确定不同岩心对应的J函数的常量参数，形成所述常量参数与岩心的流动单元指数的对应关系，并根据所述对应关系得到所述油藏原始含水饱和度模型。
10.根据权利要求9所述的计算油藏内部任意部位原始含水饱和度的系统，其特征在于，所述模型建立单元进一步用于获取每块岩心对应的J函数曲线，拟合每块岩心的所述J函数曲线得到拟合的J函数，并获取所述拟合的J函数的常量参数。
11.根据权利要求9所述的计算油藏内部任意部位原始含水饱和度的系统，其特征在于，所述常量参数包括第一参数和第二参数，所述模型建立单元进一步用于分别形成表示流动单元指数与所述第一参数的对应关系的第一对应关系和表示流动单元指数与所述第二参数的对应关系的第二对应关系。
12.根据权利要求9所述的计算油藏内部任意部位原始含水饱和度的系统，其特征在于，所述模型建立单元进一步用于根据J函数的定义、J函数与原始含水饱和度的关系、地层条件下与实验室条件下毛管力的关系以及毛管力与油水密度差的关系得到原始含水饱和度与J函数的常量参数的对应关系，根据流动单元指数与J函数的常量参数的对应关系以及原始含水饱和度与J函数的常量参数的对应关系得到原始含水饱和度与流动单元指数的对应关系，所述原始含水饱和度与流动单元指数的对应关系作为所述油藏原始含水饱和度模型。
13.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1-6任一项所述方法。
14.一种计算机可读介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-6任一项所述方法。
CN201910474989.2A 2019-06-03 2019-06-03 一种计算油藏内部任意部位原始含水饱和度的方法及系统 Active CN110287545B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910474989.2A CN110287545B (zh) 2019-06-03 2019-06-03 一种计算油藏内部任意部位原始含水饱和度的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910474989.2A CN110287545B (zh) 2019-06-03 2019-06-03 一种计算油藏内部任意部位原始含水饱和度的方法及系统

Publications (2)

Publication Number Publication Date
CN110287545A true CN110287545A (zh) 2019-09-27
CN110287545B CN110287545B (zh) 2023-02-10

Family

ID=68002994

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910474989.2A Active CN110287545B (zh) 2019-06-03 2019-06-03 一种计算油藏内部任意部位原始含水饱和度的方法及系统

Country Status (1)

Country Link
CN (1) CN110287545B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856583A (zh) * 2020-06-12 2020-10-30 中国石油天然气股份有限公司 基于孔隙度约束条件下的储层含水饱和度定量预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202763A (zh) * 2016-07-18 2016-12-07 西南石油大学 双重介质储层岩石含水饱和度的计算方法
CN106204302A (zh) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 一种计算原始含水饱和度的方法与应用
US20170292373A1 (en) * 2016-04-07 2017-10-12 Baker Hughes Incorporated Method to estimate the influence of pore-size distribution on phase equilibrium of multi-component hydrocarbon systems in unconventional shale gas and oil reservoirs
CN107701180A (zh) * 2017-09-19 2018-02-16 中国石油天然气集团公司 一种基于密闭取心的原始油藏含水饱和度计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170292373A1 (en) * 2016-04-07 2017-10-12 Baker Hughes Incorporated Method to estimate the influence of pore-size distribution on phase equilibrium of multi-component hydrocarbon systems in unconventional shale gas and oil reservoirs
CN106204302A (zh) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 一种计算原始含水饱和度的方法与应用
CN106202763A (zh) * 2016-07-18 2016-12-07 西南石油大学 双重介质储层岩石含水饱和度的计算方法
CN107701180A (zh) * 2017-09-19 2018-02-16 中国石油天然气集团公司 一种基于密闭取心的原始油藏含水饱和度计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
司马立强 等: "用J函数法求取碳酸盐岩储层饱和度方法探讨", 《岩性油气藏》 *
廖敬 等: "毛管压力曲线平均化及J函数处理", 《特种油气藏》 *
黄新波: "分流动单元平均毛管压力曲线归一化方法及应用", 《石油地质与工程》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856583A (zh) * 2020-06-12 2020-10-30 中国石油天然气股份有限公司 基于孔隙度约束条件下的储层含水饱和度定量预测方法

Also Published As

Publication number Publication date
CN110287545B (zh) 2023-02-10

Similar Documents

Publication Publication Date Title
Fenton A fifth-order Stokes theory for steady waves
CN105891089B (zh) 储层渗透率确定方法及装置
CN106121641B (zh) 油藏中贼层的识别方法及装置
CN106556877B (zh) 一种地磁通化方法及装置
CN110058310A (zh) 地震岩石的地震波速度频散及衰减预测方法及装置
CN109829217A (zh) 压裂性裂缝油藏产能模拟方法及装置
CN110489307A (zh) 接口异常调用监测方法及装置
CN104932027A (zh) 基于核磁共振测井的储层分类方法
CN104374679B (zh) 一种确定储层渗透率贡献值的方法
CN108090656A (zh) 一种确定砂体连通性的方法及装置
CN110782187B (zh) 基于含水率的枯竭油气藏型储气库动态评价方法及系统
CN110287545A (zh) 一种计算油藏内部任意部位原始含水饱和度的方法及系统
CN109597964A (zh) 理论油藏渗流参数值确定方法、系统、设备、可读介质
CN111208566B (zh) 基于sca模型的孔缝参数反演方法、装置及存储介质
CN103279651B (zh) 一种储层参数预测方法及装置
CN113392372B (zh) 气井全生命周期无阻流量确定方法及系统
CN112419493B (zh) 页岩储层三维属性模型建立方法及装置
CN109031437A (zh) 随钻核磁共振自旋回波校正方法与装置
CN110095399A (zh) 储层不可及孔隙体积确定方法及装置
CN107015289A (zh) 圈闭地质资源量确定方法及装置
CN111999331A (zh) 基于核磁测试的毛管压力确定方法及装置
CN109342273A (zh) 毛管力对pvt影响的临界孔隙半径确定方法及装置
CN111485874A (zh) 确定水平井储层含水饱和度的方法及装置
CN112052631B (zh) 完钻新井目标井段初期无阻流量参数确定方法和装置
CN108510778A (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