CN110807237A - 获取剩余可动水体比例的方法和装置 - Google Patents
获取剩余可动水体比例的方法和装置 Download PDFInfo
- Publication number
- CN110807237A CN110807237A CN201810879636.6A CN201810879636A CN110807237A CN 110807237 A CN110807237 A CN 110807237A CN 201810879636 A CN201810879636 A CN 201810879636A CN 110807237 A CN110807237 A CN 110807237A
- Authority
- CN
- China
- Prior art keywords
- water
- gas field
- field well
- formation
- proportion
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 300
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000009545 invasion Effects 0.000 claims abstract description 52
- 239000008398 formation water Substances 0.000 claims description 90
- 230000015572 biosynthetic process Effects 0.000 claims description 79
- 238000004422 calculation algorithm Methods 0.000 claims description 30
- 239000011435 rock Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 7
- 230000007423 decrease Effects 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 10
- 239000012071 phase Substances 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 2
- 229910052753 mercury Inorganic materials 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 239000008346 aqueous phase Substances 0.000 description 1
- 230000005465 channeling Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
- Pipeline Systems (AREA)
Abstract
本发明提供一种获取剩余可动水体比例的方法和装置,通过满足水流能量方程的气田井所在的地层水的参数和气田井所在的储层的参数,计算得到气田井当前剩余可动水体比例,进而可以基于气田井的剩余可动水体比例指导气田井的排水采气。由于上述方法无需再获取气田井的可动水体量,因此,可以避免现有技术中因强水侵均质模型无法准确刻画气田井的储层,导致通过强水侵均质模型无法准确获取气田井的可动水体量,进而无法通过气田井的可动水体量,指导气田井的排水采气的情况。
Description
技术领域
本发明涉及气田开发技术领域,尤其涉及一种获取剩余可动水体比例的方法和装置。
背景技术
在气田开发过程中,气田井出水后,储层变为气液两相渗流,水很快占据渗流通道而阻塞、封闭气相流动,导致气田井所在气田产能下降、采收率降低。目前,可以通过气田井的可动水体量,指导气田井的排水采气,引导水相流动,由此减弱、解除储层气相流动阻力,来提高出水后气田井所在气田的产能。
现有技术中,主要通过强水侵均质模型获取气田井的可动水体量。该强水侵均质模型在假设气田井所在的储层均质、气田井的水侵入量的大小满足物质平衡原理(即气田井的采气亏空量与“水侵入量和地层空隙变化量之和”相等)的场景下,计算气田井的可动水体量。
然而,由于强水侵均质模型无法准确的刻画气田井所在的储层,导致通过强水侵均质模型无法准确获取气田井的可动水体量,进而无法通过气田井的可动水体量,指导气田井的排水采气。
发明内容
本发明提供了一种获取剩余可动水体比例的方法和装置,能够获取气田井当前剩余可动水体比例,进而可以基于气田井的剩余可动水体比例指导气田井的排水采气。
本发明第一方面提供一种获取剩余可动水体比例的方法,该方法包括:
获取气田井所在的地层水的参数和所述气田井所在的储层的参数,所述地层水的参数包括:所述地层水在初始地层压力下的体积系数、所述地层水在当前地层压力下的体积系数、所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数、所述地层水中入侵水的当前地层压力,以及,所述入侵水的初始地层压力,所述入侵水包括:边水和/或底水;所述储层的参数包括:所述储层中岩石的有效压缩系数、所述储层的束缚水饱和度;所述地层水的参数和所述储层的参数满足水流能量方程;
根据所述地层水的参数和所述储层的参数,计算能够表征所述气田井当前的水体可排性的所述气田井当前剩余可动水体比例。
可选的,所述根据所述地层水的参数和所述储层的参数,计算所述气田井当前剩余可动水体比例,包括:
其中,所述R为所述气田井当前剩余可动水体比例,所述Bwi为所述地层水在原始地层压力下的体积系数,所述BW为所述地层水在当前地层压力下的体积系数,所述CW为所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数,所述Cf为所述储层中岩石的有效压缩系数,所述Swc为所述储层的束缚水饱和度,所述ΔP为所述入侵水的压力变化量,所述ΔP为所述入侵水的初始压力与所述入侵水的当前压力之差。
可选的,所述方法还包括:
获取所述气田井在多个时刻的剩余可动水体比例;
根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线。
可选的,所述根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线,包括:
根据所述气田井在多个时刻的剩余可动水体比例,采用预设算法,绘制所述气田井的剩余可动水体比例趋势曲线,所述预设算法包括下述至少一项:回归算法、差值算法、预测算法。
可选的,所述绘制所述气田井的剩余可动水体比例趋势曲线之后,所述方法还包括:
根据所述气田井的剩余可动水体比例趋势曲线,确定所述气田井的水体可排性。
本发明第二方面提供一种获取剩余可动水体比例的装置,该装置包括:
第一获取模块,用于获取气田井所在的地层水的参数和所述气田井所在的储层的参数,所述地层水的参数包括:所述地层水在初始地层压力下的体积系数、所述地层水在当前地层压力下的体积系数、所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数、所述地层水中入侵水的当前地层压力,以及,所述入侵水的初始地层压力,所述入侵水包括:边水和/或底水;所述储层的参数包括:所述储层中岩石的有效压缩系数、所述储层的束缚水饱和度;所述地层水的参数和所述储层的参数满足水流能量方程;
计算模块,用于根据所述地层水的参数和所述储层的参数,计算所述气田井当前剩余可动水体比例。
其中,所述R为所述气田井当前剩余可动水体比例,所述Bwi为所述地层水在原始地层压力下的体积系数,所述BW为所述地层水在当前地层压力下的体积系数,所述CW为所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数,所述Cf为所述储层中岩石的有效压缩系数,所述Swc为所述储层的束缚水饱和度,所述ΔP为所述入侵水的压力变化量,所述ΔP为所述入侵水的初始压力与所述入侵水的当前压力之差。
可选的,所述装置,还包括:
第二获取模块,用于获取所述气田井在多个时刻的剩余可动水体比例;
绘制模块,用于根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线。
可选的,所述绘制模块,具体用于根据所述气田井在多个时刻的剩余可动水体比例,采用预设算法,绘制所述气田井的剩余可动水体比例趋势曲线,所述预设算法包括下述至少一项:回归算法、差值算法、预测算法。
可选的,所述装置还包括:
确定模块,用于在所述绘制模块绘制所述气田井的剩余可动水体比例趋势曲线之后,根据所述气田井的剩余可动水体比例趋势曲线,确定所述气田井的水体可排性。
本发明第三方面提供一种获取剩余可动水体比例的装置,该装置包括:至少一个处理器和存储器;
所述存储器存储计算机执行指令;所述至少一个处理器执行所述存储器存储的计算机执行指令,以执行如第一方面和第一方面任一可能的实施方式所述的方法。
本发明第四方面提供一种计算机可读存储介质,该计算机可读存储介质中存储有程序指令,所述程序指令被处理器执行时实现第一方面和第一方面任一可能的实施方式所述的方法。
本发明第五方面提供一种计算机程序产品,该计算机程序产品包括计算机程序或指令,当所述计算机程序或指令被运行时,实现如第一方面和第一方面任一可能的实施方式所述的方法。
本发明提供的获取剩余可动水体比例的方法和装置,通过满足水流能量方程的气田井所在的地层水的参数和气田井所在的储层的参数,计算得到气田井当前剩余可动水体比例,进而可以基于气田井的剩余可动水体比例指导气田井的排水采气。由于上述方法无需再获取气田井的可动水体量,因此,可以避免现有技术中因强水侵均质模型无法准确刻画气田井的储层,导致通过强水侵均质模型无法准确获取气田井的可动水体量,进而无法通过气田井的可动水体量,指导气田井的排水采气的情况。
附图说明
图1为本发明实施例提供的气田井的示意图;
图2为本发明实施例提供的一种剩余可动水体比例的方法的流程示意图;
图3为本发明实施例提供的气田井的可动水体图;
图4为本发明提供的一种获取剩余可动水体比例的装置的结构示意图;
图5为本发明提供的另一种获取剩余可动水体比例的装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的气田井的示意图。如图1所示,在气田开发过程中,气田井1出水后,储层变为气液两相渗流,水4通过水侵通道3入侵气藏2,占据渗流通道而阻塞、封闭气相流动,导致气田井1所在气田产能下降、采收率降低。
现有技术中,主要通过强水侵均质模型获取气田井的可动水体量,以通过气田井的可动水体量,指导气田井的排水采气,引导水相流动,由此减弱、解除储层气相流动阻力,来提高出水后气田井所在气田的产能。该强水侵均质模型在假设气田井所在的储层均质、气田井的水侵入量的大小满足物质平衡原理(即气田井的采气亏空量与“水侵入量和地层空隙变化量之和”相等)的场景下,计算气田井的可动水体量。
实际中大多气田井所在的储层的非均质性很强,强水均质模型无法准确的刻画出储层的裂缝和断层。尤其是针对泄流特征复杂的气田井,以及,针对碳酸盐岩、致密砂岩等裂缝气藏的气田井,强水均质模型更无法准确的刻画出储层的裂缝和断层。但是,由于入侵水往往是通过裂缝、断层等入侵,导致通过上述强水侵均质模型无法准确获取气田井的可动水体量,进而无法通过气田井的可动水体量指导气田井的排水采气。这里所说的入侵水可以是边水或底水。
考虑到上述准确获取气田井的可动水体量比较难的问题,本发明提供了一种获取剩余可动水体比例的方法,通过使用满足水流能量方程的气田井的地层水的参数和储层的参数,可以准确的获取到气田井当前剩余可动水体比例。由于气田井当前剩余可动水体比例可以间接的反映出气田井当前的可动水体量,且水流能量方程不限定气田井所在储层的泄流特征,因此可以避免因气田井所在的储层的非均质性很导致无法准确的获取气田井的可动水体量的问题。
下面以具体地实施例对本申请的技术方案进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。
图2为本发明实施例提供的一种剩余可动水体比例的方法的流程示意图。本实施例涉及的是通过满足水流能量方程的气田井的地层水的参数和储层的参数,计算得到气田井的当前剩余可动水体比例的过程。本实施例的执行主体为任一具有处理能力的电子设备,例如:计算机、服务器、终端设备等。如图2所示,该方法包括:
S101、获取气田井所在的地层水的参数和气田井所在的储层的参数,地层水的参数包括:地层水在初始地层压力下的体积系数、地层水在当前地层压力下的体积系数、地层水从原始地层压力下降到当前地层压力时的平均压缩系数、地层水中入侵水的当前地层压力,以及,地层水中入侵水的初始地层压力;入侵水包括:边水和/或底水;储层的参数包括:储层中岩石的有效压缩系数、储层的束缚水饱和度;地层水的参数和储层的参数满足水流能量方程。
S102、根据地层水的参数和储层的参数,计算气田井当前剩余可动水体比例。
具体的,上述地层水的参数中的地层水在初始地层压力下的体积系数Bwi、地层水在当前地层压力下的体积系数Bw、地层水从原始地层压力下降到当前地层压力时的平均压缩系数CW、储层参数中的储层中岩石的有效压缩系数Cf,可以由测试人员通过查阅油气藏专业的工具书获取后,输入至电子设备。由于属于同一气藏的多个气田井的参数相同,因此,上述测试人员可以将与上述气田井属于同一气藏的其他气田井的这些参数作为该气田井的参数,输入至电子设备。
上述储层参数中的储层的束缚水饱和度Swc可以由电子设备通过该气田井的压汞数据分析获得,或,使用属于同一气藏的其他气田井的压汞数据分析获得,具体可以参见现有技术,对此不再赘述。
上述地层水的参数中的地层水中入侵水的当前地层压力,以及,地层水中入侵水的初始地层压力,可以采用如下两种方式获取:
第一种方式:测试人员通过现有的气田井压力测试,测试气田井地层水中入侵水的初始地层压力和地层水中入侵水的当前地层压力,并将气田井地层水中入侵水的初始地层压力和地层水中入侵水的当前地层压力输入至电子设备,以使电子设备获取到这两个地层压力。其中,这里所说的初始地层压力是指在气田井投产之前的地层压力,此时,入侵水还未向气田井入侵。
第二种方式:电子设备可以通过与气田井同属于同一气藏的其他气田井(简称邻井)的地层压力,以及下述公式(1)计算得到气田井的地层水中入侵水的当前地层压力或入侵水的初始地层压力。其中,公式(1)例如可以如下:
其中,z为与渗流数学模型无因次时间相对应的拉普拉斯空间自变量,表示气田井中的实际气体受到压缩后与理想气体受到同样的压力压缩后在体积上的偏差;Pf为拉普拉斯空间中气田井所在的储层的裂缝的无因次压力;ω为气田井中裂缝孔隙双重介质的储层的储容比,无因次;λ为气田井中裂缝孔隙双重介质的储层的窜流系数,无因次;x为研究单向线性渗流问题的一维坐标系中的无因次坐标值,表示邻井与该气田井的相对位置。
上述公式(1)中,设定气田井的气水交界处x等于0,水区外边界处(即气田井所在气藏的含气外边界)x等于1,x0为邻井处的无因次坐标值,P0为邻井处的无因次压力,则上述公式(1)中的参数满足如下条件:
在本实施例中,通过上述方式获取到的气田井所在的地层水的参数和气田井所在的储层的参数,满足水流能量方程。这里所说的水流能量方程不限定气田井所在储层的泄流特征,例如可以空隙型、裂缝型、孔洞型及复合型等复杂渗流气藏的气田井所在的储层。该水流能量方程可以如下述公式(5)所示:
其中,We为侵入气田井的水量,单位为104×m3;Wp为气田井的累积排水量,单位为104×m3;W为气田井的水体储量,单位为104×m3;CW的单位为MPa-1;Cf的单位为MPa-1;ΔP为地层水中入侵水的压力变化量,即地层水中入侵水的初始压力与地层水中入侵水的当前压力之差,单位为MPa。
We为侵入气田井的水量,Bw为地层水在当前地层压力下的体积系数,因此,通过将We与Bw相除,得到当前地层压力下侵入气田井的水的体积。通过将当前地层压力下侵入气田井的水的体积与气田井的累积排水量Wp相加,可以得到气田井的已入侵水体量。通过将气田井的已入侵水体量与气田井的水体储量相除,即可得到气田井的水体储量中当前已入侵水体量比例。进而通过将1与气田井当前已入侵水体量比例相减,即可得到气田井当前剩余可动水体比例R。即,随气田井的开采自由流动且能够侵入到气田井的气区的水体比例。
基于上述原理,电子设备可以采用下述公式(6)计算气田井当前剩余可动水体比例R,具体地:
由于气田井的泄流特征复杂,参数We、Wp、W可能无法准确获取。因此,电子设备通过获取满足水流能量方程的气田井的地层水的参数和气田井的储层的参数,可以利用上述公式(5)所示的储层的参数、地层水的参数与参数We、Wp、W的关系,将上述公式(6)中的参数We、Wp、W进行转换,得到下述公式(7),具体如下:
通过上述公式(7),电子设备可以根据满足水流能量方程的气田井的地层水的参数和气田井的储层的参数,准确的得到气田井当前剩余可动水体比例R。由于气田井当前剩余可动水体比例R可以间接的反映出气田井当前的可动水体量,因此,后续可以基于气田井的剩余可动水体比例R对气田井的排水采气进行指导,无需再获取气田井的可动水体量。
本发明提供的获取剩余可动水体比例的方法,通过满足水流能量方程的气田井所在的地层水的参数和气田井所在的储层的参数,计算得到气田井当前剩余可动水体比例,进而可以基于气田井的剩余可动水体比例指导气田井的排水采气。由于上述方法无需再获取气田井的可动水体量,因此,可以避免现有技术中因强水侵均质模型无法准确刻画气田井的储层,导致通过强水侵均质模型无法准确获取气田井的可动水体量,进而无法通过气田井的可动水体量,指导气田井的排水采气的情况。
随着气田井的生产、地层压力逐渐降低、地层水的入侵水不断入侵,气田井所在的地层水的参数和所述气田井所在的储层的参数随时间变化,使得气田井的剩余可动水体比例也不断变化。因此,在上述实施例的基础上,电子设备在采用上述方式获取到气田井当前剩余可动水体比例R之后,可以继续采用上述方式获取该气田井在多个时刻的剩余可动水体比例R。然后,电子设备可以根据该气田井在多个时刻的剩余可动水体比例R,绘制气田井的剩余可动水体比例的趋势曲线。该趋势曲线为时间函数曲线,通过该趋势曲线,可以宏观的反映出气田井随时间变化的剩余可动水体比例。其中,上述多个时刻中任意两个相邻时刻的时间间隔可以相同,也可以不同。这里所说的时间间隔的大小具体可以根据用户的需求设定。
具体实现时,电子设备可以根据气田井在多个时刻的剩余可动水体比例,采用预设算法,绘制气田井的剩余可动水体比例趋势曲线。这里所说的预设算法例如可以包括下述至少一项:回归算法、差值算法、预测算法。通过这种方式,可以使用较少的时刻的剩余可动水体比例,得到绘制气田井的剩余可动水体比例趋势曲线,提高了绘制气田井的剩余可动水体比例趋势曲线的效率。
在得到上述气田井的剩余可动水体比例趋势曲线之后,电子设备可以根据气田井的剩余可动水体比例趋势曲线,确定气田井的水体可排性。
具体实现时,电子设备可以结合气田井的剩余可动水体比例趋势曲线,以及,气田井排水线,绘制气田井的可动水体图。进而,电子设备可以根据气田井的可动水体图,确定气田井的水体可排性。
图3为本发明实施例提供的气田井的可动水体图。如图3所示,气田井的可动水体图包括:气田井的剩余可动水体比例趋势曲线、无水采气期、气田井排水线等。可选的,气田井的可动水体图还可以根据实际需求在当前基础上新增或删减信息,本申请对此不限定。
在气田井的可动水体图中,气田井排水线用于确定气田井的水体可排性。其中,气田井排水线的数量具体可以根据可排性等级划分的数量确定。排水线的位置可以根据气田井排水的难易,以及,其他气藏经验类比等方法确定。
以可排性等级包括可排性强、可排性中等、可排性差为例,则上述排水线包括排水线1、排水线2、排水线3。排水线1、排水线2、排水线3的位置如图3所示。当气田井的剩余可动水体比例趋势曲线位于排水线1和排水线2之间时,气田井的水体可排性为可排性强。当气田井的剩余可动水体比例趋势曲线位于排水线2和排水线3之间时,气田井的水体可排性为可排性中等。当气田井的剩余可动水体比例趋势曲线位于排水线3之外时,气田井的水体可排性为可排性差。假定上述方式所绘制的气田井的可动水体比例趋势曲线如图3所示,则该气田井的的水体可排性为可排性强,可作为排水井排水保护气藏。
本发明提供的获取剩余可动水体比例的方法,可以根据气田井的剩余可动水体比例趋势曲线,确定气田井的水体可排性,以基于气田井的水体可排性对气田井的排水采气进行指导,提高了指导气田井排水采气的效率。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
图4为本发明提供的一种获取剩余可动水体比例的装置的结构示意图。该获取剩余可动水体比例的装置可以通过软件、硬件或者两者的结合实现前述电子设备的部分或者全部。如图4所示,该获取剩余可动水体比例的装置可以包括:第一获取模块11和计算模块12。其中,
第一获取模块11,用于获取气田井所在的地层水的参数和所述气田井所在的储层的参数,所述地层水的参数包括:所述地层水在初始地层压力下的体积系数、所述地层水在当前地层压力下的体积系数、所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数、所述地层水中入侵水的当前地层压力,以及,所述入侵水的初始地层压力,所述入侵水包括:边水和/或底水;所述储层的参数包括:所述储层中岩石的有效压缩系数、所述储层的束缚水饱和度;所述地层水的参数和所述储层的参数满足水流能量方程;
计算模块12,用于根据所述地层水的参数和所述储层的参数,计算所述气田井当前剩余可动水体比例。
可选的,所述计算模块12,具体用于根据计算所述气田井的剩余可动水体比例;其中,所述R为所述气田井当前剩余可动水体比例,所述Bwi为所述地层水在原始地层压力下的体积系数,所述BW为所述地层水在当前地层压力下的体积系数,所述CW为所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数,所述Cf为所述储层中岩石的有效压缩系数,所述Swc为所述储层的束缚水饱和度,所述ΔP为所述入侵水的压力变化量,所述ΔP为所述入侵水的初始压力与所述入侵水的当前压力之差。
继续参照图4,可选的,在一些实施例中,上述装置还可以包括:第二获取模块13和绘制模块14。其中,
第二获取模块13,用于获取所述气田井在多个时刻的剩余可动水体比例;
绘制模块14,用于根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线。可选的,所述绘制模块14,具体用于根据所述气田井在多个时刻的剩余可动水体比例,采用预设算法,绘制所述气田井的剩余可动水体比例趋势曲线,所述预设算法包括下述至少一项:回归算法、差值算法、预测算法。
继续参照图4,在该实现方式下,上述装置还可以包括:
确定模块15,用于在所述绘制模块14绘制所述气田井的剩余可动水体比例趋势曲线之后,根据所述气田井的剩余可动水体比例趋势曲线,确定所述气田井的水体可排性。
本发明提供的获取剩余可动水体比例的装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
图5为本发明提供的另一种获取剩余可动水体比例的装置的结构示意图。如图5所示,该装置可以包括:至少一个处理器21和存储器22。图5示出的是以一个处理器为例的装置,其中,
存储器22,用于存放程序。具体地,程序可以包括程序代码,所述程序代码包括计算机操作指令。存储器22可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatile memory),例如至少一个磁盘存储器。
处理器21用于执行所述存储器22存储的计算机执行指令,以实现上述实施例中的获取剩余可动水体比例的方法,其实现原理和技术效果类似,在此不再赘述。
其中,处理器21可能是一个中央处理器(Central Processing Unit,简称为CPU),或者是特定集成电路(Application Specific Integrated Circuit,简称为ASIC),或者是被配置成实施本发明实施例的一个或多个集成电路。
可选的,在具体实现上,如果通信接口、存储器22和处理器21独立实现,则通信接口、存储器22和处理器21可以通过总线相互连接并完成相互间的通信。所述总线可以是工业标准体系结构(Industry Standard Architecture,简称为ISA)总线、外部设备互连(Peripheral Component,简称为PCI)总线或扩展工业标准体系结构(Extended IndustryStandard Architecture,简称为EISA)总线等。所述总线可以分为地址总线、数据总线、控制总线等,但并不表示仅有一根总线或一种类型的总线。
可选的,在具体实现上,如果通信接口、存储器22和处理器21集成在一块芯片上实现,则通信接口、存储器22和处理器21可以通过内部接口完成相同间的通信。
本发明还提供了一种计算机可读存储介质,该计算机可读存储介质可以包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random AccessMemory)、磁盘或者光盘等各种可以存储程序代码的介质。具体的,该计算机可读存储介质中存储有程序指令,程序指令用于上述实施例中的方法。
本发明还提供了一种计算机程序产品,该计算机程序产品包括计算机程序或指令,当所述计算机程序或指令被运行时,可以实现上述实施例中的方法。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种获取剩余可动水体比例的方法,其特征在于,包括:
获取气田井所在的地层水的参数和所述气田井所在的储层的参数,所述地层水的参数包括:所述地层水在初始地层压力下的体积系数、所述地层水在当前地层压力下的体积系数、所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数、所述地层水中入侵水的当前地层压力,以及,所述入侵水的初始地层压力,所述入侵水包括:边水和/或底水;所述储层的参数包括:所述储层中岩石的有效压缩系数、所述储层的束缚水饱和度;所述地层水的参数和所述储层的参数满足水流能量方程;
根据所述地层水的参数和所述储层的参数,计算所述气田井当前剩余可动水体比例。
3.根据权利要求1所述的方法,其特征在于,所述方法还包括:
获取所述气田井在多个时刻的剩余可动水体比例;
根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线。
4.根据权利要求3所述的方法,其特征在于,所述根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线,包括:
根据所述气田井在多个时刻的剩余可动水体比例,采用预设算法,绘制所述气田井的剩余可动水体比例趋势曲线,所述预设算法包括下述至少一项:回归算法、差值算法、预测算法。
5.根据权利要求3所述的方法,其特征在于,所述绘制所述气田井的剩余可动水体比例趋势曲线之后,所述方法还包括:
根据所述气田井的剩余可动水体比例趋势曲线,确定所述气田井的水体可排性。
6.一种获取剩余可动水体比例的装置,其特征在于,包括:
第一获取模块,用于获取气田井所在的地层水的参数和所述气田井所在的储层的参数,所述地层水的参数包括:所述地层水在初始地层压力下的体积系数、所述地层水在当前地层压力下的体积系数、所述地层水从原始地层压力下降到当前地层压力时的平均压缩系数、所述地层水中入侵水的当前地层压力,以及,所述入侵水的初始地层压力,所述入侵水包括:边水和/或底水;所述储层的参数包括:所述储层中岩石的有效压缩系数、所述储层的束缚水饱和度;所述地层水的参数和所述储层的参数满足水流能量方程;
计算模块,用于根据所述地层水的参数和所述储层的参数,计算所述气田井当前剩余可动水体比例。
8.根据权利要求6所述的装置,其特征在于,所述装置,还包括:
第二获取模块,用于获取所述气田井在多个时刻的剩余可动水体比例;
绘制模块,用于根据所述气田井在多个时刻的剩余可动水体比例,绘制所述气田井的剩余可动水体比例趋势曲线。
9.根据权利要求8所述的装置,其特征在于,所述绘制模块,具体用于根据所述气田井在多个时刻的剩余可动水体比例,采用预设算法,绘制所述气田井的剩余可动水体比例趋势曲线,所述预设算法包括下述至少一项:回归算法、差值算法、预测算法。
10.根据权利要求8所述的装置,其特征在于,所述装置还包括:
确定模块,用于在所述绘制模块绘制所述气田井的剩余可动水体比例趋势曲线之后,根据所述气田井的剩余可动水体比例趋势曲线,确定所述气田井的水体可排性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810879636.6A CN110807237B (zh) | 2018-08-03 | 2018-08-03 | 获取剩余可动水体比例的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810879636.6A CN110807237B (zh) | 2018-08-03 | 2018-08-03 | 获取剩余可动水体比例的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110807237A true CN110807237A (zh) | 2020-02-18 |
CN110807237B CN110807237B (zh) | 2022-11-01 |
Family
ID=69486847
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810879636.6A Active CN110807237B (zh) | 2018-08-03 | 2018-08-03 | 获取剩余可动水体比例的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110807237B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104018829A (zh) * | 2014-05-23 | 2014-09-03 | 中国地质大学(北京) | 一种利用煤层气井生产数据测量气水相渗曲线的方法 |
CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
CN106484933A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于确定页岩气井井控动态储量的方法及系统 |
CN107575207A (zh) * | 2017-10-19 | 2018-01-12 | 陕西延长石油(集团)有限责任公司研究院 | 一种预测油田水驱波及半径的方法 |
-
2018
- 2018-08-03 CN CN201810879636.6A patent/CN110807237B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104018829A (zh) * | 2014-05-23 | 2014-09-03 | 中国地质大学(北京) | 一种利用煤层气井生产数据测量气水相渗曲线的方法 |
CN106484933A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于确定页岩气井井控动态储量的方法及系统 |
CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
CN107575207A (zh) * | 2017-10-19 | 2018-01-12 | 陕西延长石油(集团)有限责任公司研究院 | 一种预测油田水驱波及半径的方法 |
Non-Patent Citations (2)
Title |
---|
张翠兰等: "XL气田茅口组有水气藏水侵量计算及水侵分析", 《内江科技》 * |
彭远进等: "一种裂缝型有水气藏物质平衡新模型", 《天然气工业》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110807237B (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20170218729A1 (en) | Method and system for hydraulic fracturing based on skin factor analysis | |
CN109858701B (zh) | 裂缝性边水气藏水侵量的定量识别方法及系统 | |
US8725478B2 (en) | Reservoir upscaling method with preserved transmissibility | |
CN107451311B (zh) | 一种油水相对渗透率曲线计算方法及装置 | |
US10282385B2 (en) | Estimation of oil recovery in foam flooded hydrocarbon reservoirs | |
Abbaszadeh et al. | Experimentally-based empirical foam modeling | |
Daripa et al. | Modeling and simulation of surfactant–polymer flooding using a new hybrid method | |
Osei-Bonsu et al. | Quantitative analysis of phase topology evolution during three-phase displacements in porous media | |
CN110807237B (zh) | 获取剩余可动水体比例的方法和装置 | |
Garcia et al. | Flow instabilities during injection of CO2 into saline aquifers | |
CN108536982B (zh) | 一种多裂缝油水混合交叉驱替的评价方法 | |
Amrollahinasab et al. | Simultaneous interpretation of SCAL data with different degrees of freedom and uncertainty analysis | |
CN111577264A (zh) | 裂缝孔隙型油藏水平井产能预测方法及装置 | |
Saafan et al. | Inferring capillary pressure curve from 2D rock images based on fractal theory in low-permeability sandstone: a new integrated approach | |
CN111626530B (zh) | 确定压裂页岩气藏井控储量的方法及装置 | |
CN116291406A (zh) | 一种海上疏松砂岩油藏大孔道识别方法和系统 | |
CN111027892B (zh) | 气井的储层渗透率的确定方法、装置和服务器 | |
CN114925547A (zh) | 一种堵水剂用量确定方法、装置、电子设备和存储介质 | |
Fager et al. | Evaluation of directly simulated WAG hysteresis at pore scale and its effect on injectivity index | |
US20230400401A1 (en) | Method and system for automated in-situ capillary pressure and wettabilty characterization of porous media | |
CN110863825B (zh) | 区分特征区域的方法和装置 | |
Thomas et al. | Effect of scaleup and aggregation on the use of well tests to identify geological properties | |
CN107133880B (zh) | 一种估算开发井泄油面积的方法 | |
US11808148B2 (en) | Systems and methods for back-allocation of oil produced by waterflooding | |
CN112983373B (zh) | 储层的改造系统及储层的改造方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |