CN112081565A - 水驱油藏体积波及系数的确定方法及装置 - Google Patents
水驱油藏体积波及系数的确定方法及装置 Download PDFInfo
- Publication number
- CN112081565A CN112081565A CN201910504252.0A CN201910504252A CN112081565A CN 112081565 A CN112081565 A CN 112081565A CN 201910504252 A CN201910504252 A CN 201910504252A CN 112081565 A CN112081565 A CN 112081565A
- Authority
- CN
- China
- Prior art keywords
- water
- determining
- water content
- volume sweep
- sweep coefficient
- 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
- 238000000034 method Methods 0.000 title claims abstract description 65
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 236
- 238000004519 manufacturing process Methods 0.000 claims abstract description 44
- 239000000463 material Substances 0.000 claims abstract description 29
- 230000009466 transformation Effects 0.000 claims abstract description 23
- 230000006870 function Effects 0.000 claims description 31
- 230000035699 permeability Effects 0.000 claims description 24
- 238000012545 processing Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 6
- 238000011161 development Methods 0.000 abstract description 25
- 230000000694 effects Effects 0.000 abstract description 10
- 239000003921 oil Substances 0.000 description 104
- 238000004891 communication Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
-
- 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
-
- 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
- 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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Control Of Non-Electrical Variables (AREA)
Abstract
本申请提供一种水驱油藏体积波及系数的确定方法及装置，方法包括：根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。本申请能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发的效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
Description
技术领域
本发明涉及油层开发技术领域，具体涉及一种水驱油藏体积波及系数的确定方法及装置。
背景技术
水驱油开发过程中，水驱体积波及系数是确定水驱开发效果及采收率的重要依据，对油藏制定后期调整方向和调整措施具有重要意义。
目前水驱体积波及系数的确定方法主要有室内实验、数值模拟及水驱特征曲线法。室内实验方法受岩心平面及纵向层段影响，实验结果并不能完全反映油藏的水驱渗流规律，且实验过程需要耗费大量的人力及物力。数值模拟是通过建立油藏模型，然后模拟水驱过程，因此数值模拟需要耗费大量的时间成本及人工成本，且油藏建模及数值模拟的参数不确定性因素较多，对模拟结果的精度影响较大。水驱特征曲线法计算过程简单且容易掌握，主要是利用甲、乙、丙、丁水驱特征曲线获得曲线斜率及截距，然后通过曲线斜率及截距直接计算体积波及系数与含水率或累积产油量的关系，由于井控地质储量对体积波及系数的计算影响非常大，水驱特征曲线法通常采取容积法计算井控地质储量或者以最大可采储量代替井控地质储量，容积法计算误差较大，而最大可采储量通常远小于井控地质储量，从而导致体积波及系数偏大或者偏小。
因此，亟需一种高精度的体积波及系数来指导油田高效开发，提高油田水驱开发得效率。
发明内容
针对现有技术中的问题，本发明提供一种水驱油藏体积波及系数的确定方法及装置，能够更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。
为解决上述技术问题，本发明提供以下技术方案：
第一方面，本发明提供一种水驱油藏体积波及系数的确定方法，包括：
根据含水率和含水饱和度的变换曲线建立含水率模型；
基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；
通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
进一步的，所述根据含水率和含水饱和度的变换曲线建立含水率模型，包括：
根据含水率和含水饱和度的变换曲线，确定基于Logistic模型的含水率微分方程；
对所述含水率微分方程进行预处理得到含水率模型。
进一步的，所述对所述含水率微分方程进行预处理得到含水率模型，包括：
对所述含水率微分方程进行变量分离以及积分处理得到含水率模型。
进一步的，所述基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量，包括：
基于所述含水率模型和所述物质平衡方程，确定含水率与累计产油量之间关系的第一表达式；
根据所述含水率模型和所述莱文莱特函数，确定含水率与相对渗透率之间关系的第二表达式；
通过生产动态数据确定所述第一表达式的曲线斜率和截距以及通过油水相对渗透率曲线确定所述第二表达式的曲线斜率和截距；
基于所述第一表达式的曲线斜率和截距以及所述第二表达式的曲线斜率和截距确定井控储量。
第二方面，本发明提供一种水驱油藏体积波及系数的确定装置，包括：
建模单元，用于根据含水率和含水饱和度的变换曲线建立含水率模型；
井控储量单元，用于基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；
体积波及系数单元，用于通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
进一步的，所述建模单元，包括：
微分子单元，用于根据含水率和含水饱和度的变换曲线，确定基于Logistic模型的含水率微分方程；
建模子单元，用于对所述含水率微分方程进行预处理得到含水率模型。
进一步的，所述建模子单元，包括：
处理模块，用于对所述含水率微分方程进行变量分离以及积分处理得到含水率模型。
进一步的，所述井控储量单元，包括：
第一子单元，用于基于所述含水率模型和所述物质平衡方程，确定含水率与累计产油量之间关系的第一表达式；
第二子单元，用于根据所述含水率模型和所述莱文莱特函数，确定含水率与相对渗透率之间关系的第二表达式；
第三子单元，用于通过生产动态数据确定所述第一表达式的曲线斜率和截距以及通过油水相对渗透率曲线确定所述第二表达式的曲线斜率和截距；
第四子单元，用于基于所述第一表达式的曲线斜率和截距以及所述第二表达式的曲线斜率和截距确定井控储量。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的水驱油藏体积波及系数的确定方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的水驱油藏体积波及系数的确定方法的步骤。
由上述技术方案可知，本发明提供一种水驱油藏体积波及系数的确定方法及装置，通过根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数，能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的水驱油藏体积波及系数的确定装置的一种通信结构示意图。
图2为本发明的水驱油藏体积波及系数的确定装置的另一种通信结构示意图。
图3为本发明实施例中的水驱油藏体积波及系数的确定方法的流程示意图。
图4为本发明实施例提供的水驱油藏体积波及系数的确定方法中含水率、含水率导数和含水饱和度之间的关系曲线示意图。
图5为本发明实施例提供的水驱油藏体积波及系数的确定方法中油水相对渗透率的曲线示意图。
图6为本发明实施例提供的水驱油藏体积波及系数的确定方法中W与T的关系曲线示意图。
图7为本发明实施例提供的水驱油藏体积波及系数的确定方法中Y与X的关系曲线示意图。
图8为本发明实施例中的水驱油藏体积波及系数的确定装置的结构示意图。
图9为本发明实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
考虑到现有的对水驱油藏体积波及系数的确定方法中存在的准确低和误差大的问题。本发明提供一种水驱油藏体积波及系数的确定方法、水驱油藏体积波及系数的确定装置、电子设备及计算机可读存储介质，通过根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数，能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
基于上述内容，本发明还提供水驱油藏体积波及系数的确定装置，该装置可以为一种服务器A1，参见图1，该服务器A1可以与客户端设备B1通信连接，用户可以将水驱油藏资料及其他相关数据输入所述客户端设备B1，所述客户端设备B1可以在线将水驱油藏资料及其他相关数据发送至服务器A1，所述服务器A1可以在线接收所述客户端设备B1发送的水驱油藏资料及其他相关数据，而后离线或在线根据水驱油藏资料获取对应的水驱油藏体积波及系数，根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。而后，所述服务器A1在线将体积波及系数发送至所述客户端设备B1，使得用户经由所述客户端设备B1获知体积波及系数数据。
进一步来说，上述服务器A1还可以与一水驱油藏资料采集设备C1通信连接，参见图2，该水驱油藏资料采集设备C1可以直接自目标区域获取水驱油藏资料及其他相关数据，也可以与一数据库D1通信连接，自该数据库D1中获取对应的水驱油藏资料及其他相关数据。而后，所述水驱油藏资料采集设备C1将水驱油藏资料及其他相关数据发送至所述服务器A1。
可以理解的是，所述客户端设备B1可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中，进行水驱油藏体积波及系数的确定的部分可以在如上述内容所述的服务器A1侧执行，即，如图1所示的架构，也可以所有的操作都在所述客户端设备B1中完成。具体可以根据所述客户端设备B1的处理能力，以及用户使用场景的限制等进行选择。本发明对此不作限定。若所有的操作都在所述客户端设备B1中完成，所述客户端设备B1还可以包括处理器，用于进行水驱油藏体积波及系数的确定的具体处理。
上述的客户端设备可以具有通信模块(即通信单元)，可以与远程的服务器进行通信连接，实现与所述服务器的数据传输。例如，通信单元可以将上述用户输入的水驱油藏资料及其他相关数据发送至服务器，以便服务器根据这些水驱油藏资料及其他相关数据进行水驱油藏体积波及系数的确定。通信单元还可以接收服务器返回的确定结果。所述服务器可以包括任务调度中心一侧的服务器，其他的实施场景中也可以包括中间平台的服务器，例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备，也可以包括多个服务器组成的服务器集群，或者分布式装置的服务器结构。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信，包括在本发明提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然，所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol，远程过程调用协议)、REST协议(Representational State Transfer，表述性状态转移协议)等。
为了能够有效提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率，本发明提供一种水驱油藏体积波及系数的确定方法的实施例，参见图3，具体包含有如下内容：
S101：根据含水率和含水饱和度的变换曲线建立含水率模型；
可以理解的是，含水率可以是描述储层产液性质最直接的参数，一般可以根据生产动态资料定量计算含水率。含水率是一个有限的值，随含水饱和度的增加而逐渐增加，数值大小介于0到0.98之间。
对于任意一个确定的油田，水驱开发过程中含水率fw(Sw)与含水饱和度Sw普遍具有图4所示的规律特征。随着含水饱和度的增加，含水率导数fw’(Sw)首先逐渐增加，达到高峰后又逐渐下降。结合含水率fw(Sw)曲线及导数fw’(Sw)曲线的变化规律，含水率fw(Sw)随含水饱和度Sw的变化关系普遍符合S型曲线(又称为Logistic曲线)增长规律，则根据含水率和含水饱和度的变换曲线建立含水率模型，具体包括：
S1011：根据含水率和含水饱和度的变换曲线，确定基于Logistic模型的含水率微分方程；
其中，基于Logistic模型的含水率微分方程为：
式中：fw(Sw)为含水率，％；Sw为含水饱和度，％；fwmax是含水率的最大值，一般为98％，％；β是模型常数。
S1012对所述含水率微分方程进行预处理得到含水率模型。
在本步骤中，对所述含水率微分方程进行变量分离以及积分处理得到含水率模型，具体包括：
分离变量且积分：
式中：Swi是初始含水饱和度，％；fwi(Swi)是初始含水率，％。
对公式(3)进行简化得到含水率模型：
S102：基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；
在本步骤中，在含水率模型基础上，结合物质平衡方程和莱文莱特函数确定井控储量，能够提高井控储量的精度，具体包括：
S1021：基于所述含水率模型和所述物质平衡方程，确定含水率与累计产油量之间关系的第一表达式；
在本步骤中，确定油井生产动态数据和物质平衡方程，其中，物质平衡方程为：
式中，Np是含水饱和度从Swi变化到Sw时的累积产油量，104t；N是井控储量，104t。
结合含水率模型和物质平衡方程，即将公式(5)代入公式(4)中，则：
令：
则公式(6)为：
W＝a+bT (8)。
需要说明的是，公式(6)即为第一表达式。
S1022：根据所述含水率模型和所述莱文莱特函数，确定含水率与相对渗透率之间关系的第二表达式；
在本步骤中，将含水率模型(公式4)变换为：
确定莱文莱特函数：
式中：Qw是油藏的年产水量，Qo是油藏的年产油量，Kro是油相相对渗透率，无因次；Krw是水相相对渗透率，无因次；μw是水粘度，mpa.s；μo是油粘度，mpa.s。
结合含水率模型和莱文莱特函数，即，将公式(10)代入公式(9)中：
令：
则公式(11)为：
Y＝a′+β′X (13)；
需要说明的是，公式(13)即为第二表达式。
S1023：通过生产动态数据确定所述第一表达式的曲线斜率和截距以及通过油水相对渗透率曲线确定所述第二表达式的曲线斜率和截距；
在本步骤中，将生产动态数据中的多项含水率和每项含水率对应的累计产油量代入公式(7)中，确定多个参数W和多个参数T，将多个参数W和多个参数T代入公式(8)中确定曲线斜率b和截距a的具体值。
根据相对渗透率曲线，采用不同含水饱和度Sw(Sw≥Swi)下的油相相对渗透率值、水相相对渗透率值代入公式(12)中求取多个参数X和多个参数Y的值。将不同含水饱和度的参数X和参数Y的值代入公式(13)作曲线，确定曲线斜率β′和截距a′。
曲线斜率b和截距a是通过公式(8)及油藏动态数据拟合获得，曲线斜率β′和截距a′是通过公式(13)及相对渗透率曲线静态数据拟合获得。公式(8)及公式(13)都是由公式(4)变换而来，理论上截距a的值等于截距a′的值。实际上，相对渗透率曲线作为静态数据往往与油藏实际生产动态存在一定的误差，因此需要引入常数C对a′、β′值进行修正，则公式(13)为：
Y＝Ca′+Cβ′X＝a+β”X (14)；
因此，对β′值进行修正后的β”值为：
S104：基于所述第一表达式的曲线斜率和截距以及所述第二表达式的曲线斜率和截距确定井控储量。
在本步骤中，将公式(15)中的β”作为更新后的β代入公式(7)中得到井控储量N为：
S103：通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
在本步骤中，累积产油量Np可表示为：
Np＝NRo (17)；
采出程度Ro可表示为：
Ro＝EvEd (18)；
式中：Ev为体积波及系数，小数；Ed为驱油效率，小数。
驱油效率Ed为：
Ed＝(1-Swi-Sor)/(1-Swi) (19)；
式中：Swi为初始含水饱和度，Sor为残余油饱和度。
将公式(17)和公式(19)代入公式(18)，确定体积波及系数Ev为：
本实施例提供的水驱油藏体积波及系数的确定方法能够提高水驱油藏体积波及系数的精度，根据该精度更高的体积波及系数应用在驱油效率的计算中，可以更准确的确定油藏水驱的规律，提高驱油的效率，进而有效提高油田水驱开发的准确性，并缩短开发期，提高开发效率。
将本实施例中的体积波及系数应用在采出程度的计算中，可以提高油藏的采出程度的计算精度，进而可以根据采出程度确定采收率，更好的指导油田高效的开发，而且对油藏后期的方向调整和措施调整具有重要意义。
从上述描述可知，本实施例提供的水驱油藏体积波及系数的确定方法，通过根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数，能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
进一步地，在计算得到与所述待测含水率相对应的水驱体积波及系数之后，可以根据所述待测含水率以及与所述待测含水率相对应的水驱体积波及系数，绘制相应的图版，确定所述待测油藏的水驱规律以及能达到的开采程度，这样可以更好的分析油田水驱开发效果，指导油田高效开发。同时，通过确定待测油藏的最终波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
为进一步地说明本方案，本发明提供了一种水驱油藏体积波及系数的确定方法的实例，具体包含有如下内容：
C井的油藏为块状边底水油藏，地层平均有效厚度22.63m，平均含油面积8.8525km2，平均孔隙度20.02％，平均有效渗透率284.4×10-3μm2，原始平均含油饱和度76％，含水饱和度Swi为0.24，原始地层压力13.9MPa，原始地层温度28℃，地下原油粘度2.46mPa.s，地下原油密度0.885g/cm3，油的体积系数1.032，地下水粘度0.51mPa.s，残余油饱和度Sor约0.21。C井生产动态数据如表1，油藏相对渗透率曲线如图5所示。
表1 C井生产动态数据
根据生产动态数据，结合公式(7)和公式(8)作如图6所示的W与T的关系曲线，从而求得参数a为-3.743，参数b为0.9563。
根据相对渗透率曲线，结合公式(12)和公式(13)作如图7所示的Y与X的关系曲线，从而求得参数a′为-3.572，参数β′为15.873。
根据公式(15)及(16)计算出井控储量N为13.21864×104t。最后，根据公式(20)可求得不同累积产量对应的体积波及系数，如表2所示。
表2体积波及系数的预测结果
从上述描述可知，本申请根据水驱开发过程中含水率、含水率导数随含水饱和度的变化规律特征，建立了含水率的S型曲线增长模型。在含水率模型基础上，结合油井生产动态数据及油藏参数、相对渗透率曲线静态数据，在合理确定井控储量基础上，能够简单、快速、有效的确定水驱油藏体积波及系数。
本发明实施例提供一种能够实现所述水驱油藏体积波及系数的确定方法中全部内容的水驱油藏体积波及系数的确定装置的具体实施方式，参见图8，水驱油藏体积波及系数的确定装置具体包括如下内容：
建模单元10，用于根据含水率和含水饱和度的变换曲线建立含水率模型；
井控储量单元20，用于基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；
体积波及系数单元30，用于通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
进一步的，所述建模单元10，包括：
微分子单元，用于根据含水率和含水饱和度的变换曲线，确定基于Logistic模型的含水率微分方程；
建模子单元，用于对所述含水率微分方程进行预处理得到含水率模型。
进一步的于，所述建模子单元，包括：
处理模块，用于对所述含水率微分方程进行变量分离以及积分处理得到含水率模型。
进一步的，所述井控储量单元20，包括：
第一子单元，用于基于所述含水率模型和所述物质平衡方程，确定含水率与累计产油量之间关系的第一表达式；
第二子单元，用于根据所述含水率模型和所述莱文莱特函数，确定含水率与相对渗透率之间关系的第二表达式；
第三子单元，用于通过生产动态数据确定所述第一表达式的曲线斜率和截距以及通过油水相对渗透率曲线确定所述第二表达式的曲线斜率和截距；
第四子单元，用于基于所述第一表达式的曲线斜率和截距以及所述第二表达式的曲线斜率和截距确定井控储量。
本发明提供的水驱油藏体积波及系数的确定装置的实施例具体可以用于执行上述实施例中的水驱油藏体积波及系数的确定方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
从上述描述可知，本发明实施例提供的一种水驱油藏体积波及系数的确定装置，通过根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数，能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
本发明的实施例还提供能够实现上述实施例中的水驱油藏体积波及系数的确定方法中全部步骤的一种电子设备的具体实施方式，参见图9，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的水驱油藏体积波及系数的确定方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
从上述描述可知，本发明实施例提供的电子设备，通过根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数，能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
本发明的实施例还提供能够实现上述实施例中的水驱油藏体积波及系数的确定方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的水驱油藏体积波及系数的确定方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
从上述描述可知，本发明实施例提供的计算机可读存储介质，通过根据含水率和含水饱和度的变换曲线建立含水率模型；基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数，能够提高水驱体积波及系数的精度，可以更好的指导油田高效开发，确定油藏水驱规律，提高油田水驱开发得效率。进一步地，通过确定待测油藏的体积波及系数，可以准确的获取待测油藏所能够达到的最佳开采效果。
虽然本发明提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。术语“上”、“下”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
本发明的说明书中，说明了大量具体细节。然而能够理解的是，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。类似地，应当理解，为了精简本发明公开并帮助理解各个发明方面中的一个或多个，在上面对本发明的示例性实施例的描述中，本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该公开的方法解释呈反映如下意图：即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如权利要求书所反映的那样，发明方面在于少于前面公开的单个实施例的所有特征。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。本发明并不局限于任何单一的方面，也不局限于任何单一的实施例，也不局限于这些方面和/或实施例的任意组合和/或置换。而且，可以单独使用本发明的每个方面和/或实施例或者与一个或更多其他方面和/或其实施例结合使用。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (10)
1.一种水驱油藏体积波及系数的确定方法，其特征在于，包括：
根据含水率和含水饱和度的变换曲线建立含水率模型；
基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；
通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
2.根据权利要求1所述的水驱油藏体积波及系数的确定方法，其特征在于，所述根据含水率和含水饱和度的变换曲线建立含水率模型，包括：
根据含水率和含水饱和度的变换曲线，确定基于Logistic模型的含水率微分方程；
对所述含水率微分方程进行预处理得到含水率模型。
3.根据权利要求2所述的水驱油藏体积波及系数的确定方法，其特征在于，所述对所述含水率微分方程进行预处理得到含水率模型，包括：
对所述含水率微分方程进行变量分离以及积分处理得到含水率模型。
4.根据权利要求1所述的水驱油藏体积波及系数的确定方法，其特征在于，所述基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量，包括：
基于所述含水率模型和所述物质平衡方程，确定含水率与累计产油量之间关系的第一表达式；
根据所述含水率模型和所述莱文莱特函数，确定含水率与相对渗透率之间关系的第二表达式；
通过生产动态数据确定所述第一表达式的曲线斜率和截距以及通过油水相对渗透率曲线确定所述第二表达式的曲线斜率和截距；
基于所述第一表达式的曲线斜率和截距以及所述第二表达式的曲线斜率和截距确定井控储量。
5.一种水驱油藏体积波及系数的确定装置，其特征在于，包括：
建模单元，用于根据含水率和含水饱和度的变换曲线建立含水率模型；
井控储量单元，用于基于所述含水率模型，结合物质平衡方程和莱文莱特函数确定井控储量；
体积波及系数单元，用于通过所述井控储量、累计产油量、初始含水饱和度和残余油饱和度确定体积波及系数。
6.根据权利要求5所述的水驱油藏体积波及系数的确定装置，其特征在于，所述建模单元，包括：
微分子单元，用于根据含水率和含水饱和度的变换曲线，确定基于Logistic模型的含水率微分方程；
建模子单元，用于对所述含水率微分方程进行预处理得到含水率模型。
7.根据权利要求6所述的水驱油藏体积波及系数的确定装置，其特征在于，所述建模子单元，包括：
处理模块，用于对所述含水率微分方程进行变量分离以及积分处理得到含水率模型。
8.根据权利要求5所述的水驱油藏体积波及系数的确定装置，其特征在于，所述井控储量单元，包括：
第一子单元，用于基于所述含水率模型和所述物质平衡方程，确定含水率与累计产油量之间关系的第一表达式；
第二子单元，用于根据所述含水率模型和所述莱文莱特函数，确定含水率与相对渗透率之间关系的第二表达式；
第三子单元，用于通过生产动态数据确定所述第一表达式的曲线斜率和截距以及通过油水相对渗透率曲线确定所述第二表达式的曲线斜率和截距；
第四子单元，用于基于所述第一表达式的曲线斜率和截距以及所述第二表达式的曲线斜率和截距确定井控储量。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述的水驱油藏体积波及系数的确定方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至4任一项所述的水驱油藏体积波及系数的确定方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910504252.0A CN112081565B (zh) | 2019-06-12 | 2019-06-12 | 水驱油藏体积波及系数的确定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910504252.0A CN112081565B (zh) | 2019-06-12 | 2019-06-12 | 水驱油藏体积波及系数的确定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112081565A true CN112081565A (zh) | 2020-12-15 |
| CN112081565B CN112081565B (zh) | 2023-01-10 |
Family
ID=73733215
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910504252.0A Active CN112081565B (zh) | 2019-06-12 | 2019-06-12 | 水驱油藏体积波及系数的确定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112081565B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116861714A (zh) * | 2023-09-05 | 2023-10-10 | 西南石油大学 | 一种确定缝洞型油藏水驱波及程度的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2480584C1 (ru) * | 2011-10-26 | 2013-04-27 | федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Пермский национальный исследовательский политехнический университет" | Способ оперативного прогнозирования основных показателей разработки нефтяных залежей |
| CN106150477A (zh) * | 2015-04-23 | 2016-11-23 | 中国石油化工股份有限公司 | 一种确定单井控制储量的方法 |
| CN107194180A (zh) * | 2017-05-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 油藏水驱体积波及系数的确定方法、装置及系统 |
| CN109424362A (zh) * | 2017-08-29 | 2019-03-05 | 中国石油化工股份有限公司 | 计算底水油藏单井控制原油储量的方法及系统 |
-
2019
- 2019-06-12 CN CN201910504252.0A patent/CN112081565B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2480584C1 (ru) * | 2011-10-26 | 2013-04-27 | федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Пермский национальный исследовательский политехнический университет" | Способ оперативного прогнозирования основных показателей разработки нефтяных залежей |
| CN106150477A (zh) * | 2015-04-23 | 2016-11-23 | 中国石油化工股份有限公司 | 一种确定单井控制储量的方法 |
| CN107194180A (zh) * | 2017-05-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 油藏水驱体积波及系数的确定方法、装置及系统 |
| CN109424362A (zh) * | 2017-08-29 | 2019-03-05 | 中国石油化工股份有限公司 | 计算底水油藏单井控制原油储量的方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| 孙恩慧等: "计算底水油藏单井水驱波及体积系数的新方法", 《石油化工应用》 * |
| 孙来喜等: "基于Logistic模型及水驱曲线的一种综合预测模型", 《矿物岩石》 * |
| 赵艳武等: "水驱油田特高含水期含水率预测模型", 《特种油气藏》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116861714A (zh) * | 2023-09-05 | 2023-10-10 | 西南石油大学 | 一种确定缝洞型油藏水驱波及程度的方法 |
| CN116861714B (zh) * | 2023-09-05 | 2023-11-24 | 西南石油大学 | 一种确定缝洞型油藏水驱波及程度的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112081565B (zh) | 2023-01-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108898218A (zh) | 一种神经网络模型的训练方法、装置、及计算机设备 | |
| CN109709607B (zh) | 一种预测薄层砂岩储层厚度方法及装置 | |
| CN111101929A (zh) | 一种计算油气藏平均地层压力的方法、装置及系统 | |
| US10235478B2 (en) | Pseudo-phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation | |
| CN104216884A (zh) | 一种基于时间衰减的协同过滤系统及方法 | |
| US20130332120A1 (en) | System and method for aggregating reservoir connectivities | |
| CN112081565B (zh) | 水驱油藏体积波及系数的确定方法及装置 | |
| US10060228B2 (en) | Pseudo phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation in order to rank multiple petro-physical realizations | |
| CN108681793B (zh) | 深层油藏采油指数预测方法及装置 | |
| CN111830559A (zh) | 一种地震资料拓频处理的方法及装置 | |
| CN111983679A (zh) | 基于沉积参数的相控型白云岩储层地震预测方法和装置 | |
| CN111434889A (zh) | 带凝析气顶油藏中注入水配比的获取方法及系统 | |
| CN112329232B (zh) | 缝洞型油藏生产动态表征方法、装置、设备及存储介质 | |
| CN113610292A (zh) | 一种致密气藏压裂水平井动态储量的确定方法和装置 | |
| CN114201668A (zh) | 推荐系统的数据处理方法、装置、电子设备和存储介质 | |
| CN113033857B (zh) | 预测直线递减生产井初始年产量的方法和装置 | |
| CN113033858B (zh) | 预测原油技术可采储量的方法和装置 | |
| US10178034B1 (en) | Iterative autocorrelation function calculation for streamed data using components | |
| CN113033855B (zh) | 预测油井技术极限递减生产日历年数的方法和装置 | |
| CN112668137B (zh) | 重力流砂岩有效厚度的校正方法及校正装置 | |
| CN113027414B (zh) | 预测原油经济可采储量的方法和装置 | |
| CN113033856A (zh) | 预测直线递减生产井初始递减率的方法和装置 | |
| CN111292205A (zh) | 一种司法数据解析方法、装置、设备及存储介质 | |
| US10235414B1 (en) | Iterative kurtosis calculation for streamed data using components | |
| CN112052631B (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 |