CN116467967A - 一种油气藏配产数据确定和能源井一体化模拟方法及装置 - Google Patents
一种油气藏配产数据确定和能源井一体化模拟方法及装置 Download PDFInfo
- Publication number
- CN116467967A CN116467967A CN202310436435.XA CN202310436435A CN116467967A CN 116467967 A CN116467967 A CN 116467967A CN 202310436435 A CN202310436435 A CN 202310436435A CN 116467967 A CN116467967 A CN 116467967A
- Authority
- CN
- China
- Prior art keywords
- data
- target
- well
- processing
- model
- 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
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 147
- 238000000034 method Methods 0.000 title claims abstract description 107
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 82
- 238000012545 processing Methods 0.000 claims abstract description 122
- 230000001052 transient effect Effects 0.000 claims abstract description 100
- 238000012512 characterization method Methods 0.000 claims abstract description 66
- 230000008569 process Effects 0.000 claims description 34
- 238000004590 computer program Methods 0.000 claims description 19
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 239000004215 Carbon black (E152) Substances 0.000 claims description 16
- 229930195733 hydrocarbon Natural products 0.000 claims description 16
- 150000002430 hydrocarbons Chemical class 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 15
- 238000013461 design Methods 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 7
- 239000012530 fluid Substances 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 5
- 230000010354 integration Effects 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 18
- 238000004891 communication Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000003190 augmentative effect Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本说明书涉及油气藏开发生产技术领域,尤其涉及一种油气藏配产数据确定和能源井一体化模拟方法及装置。其中油气藏配产数据确定方法包括利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据;利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程,表征数据由饱和度数据确定;联立流动表现方程和垂直提升表现方程,求解得到交点数据;以及基于地面管网模型针对与终端分离器对应的参数和交点数据进行处理,得到目标流压数据和目标配产量数据。利用本说明书实施例,实现了在气藏地下渗流、井筒多相管流和地面管网集输一体化的基础上,确定目标流压数据和目标配产量数据。
Description
技术领域
本说明书涉及油气藏开发生产技术领域,尤其涉及一种油气藏配产数据确定和能源井一体化模拟方法及装置。
背景技术
目前在确定用于油气藏配产的目标配产量数据均是通过在气藏地下渗流的技术基础上确定的。由于没有考虑气井举升能力、地面管网集输能力,使得确定的目标配产量数据不准确,无法合理地指导采油采气工作。
如何在气藏地下渗流、井筒举升和地面管网集输三者的基础上,确定合理的目标配产量数据是现有技术中亟需解决的问题。
发明内容
为解决现有技术中的问题,本说明书实施例提供了一种油气藏配产数据确定和能源井一体化模拟方法及装置。基于气藏数值模拟模型、井筒模型和地面管网模型,实现了确定合理的目标配产量数据。
为了解决上述技术问题,本说明书的具体技术方案如下:
一方面,本说明书实施例提供了一种油气藏配产数据确定方法,包括,
利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据,所述第一瞬态井数据与第一采集时刻相对应;
利用井筒模型针对表征数据和所述流动表现方程进行处理,得到垂直提升表现方程,所述表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由所述饱和度数据确定子表征数据;
联立所述流动表现方程和所述垂直提升表现方程,求解得到交点数据;以及
基于地面管网模型针对与终端分离器对应的参数和所述交点数据进行处理,得到目标流压数据和目标配产量数据。
进一步,该利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据进一步包括:
采集与所述能源井的井距一半范围内的多个网格分别对应的第一瞬态数据;
确定所述第一瞬态数据的平均值为稳态数据;以及
利用所述气藏数值模拟模型针对所述稳态数据进行处理,得到所述流动表现方程和所述饱和度数据。
进一步,该第一瞬态井数据进一步包括气藏压力、渗透率、储层厚度、井筒半径、流体状态方程和供给面积。
进一步,该气藏数值模拟模型包括拟稳态气井产能公式,所述拟稳态气井产能公式包括:
其中,Q表征气产量,K表征平均渗透率,h表征平均储层厚度,Tsc表征标准条件下的温度,Pr表征平均气藏压力,Pwf表征井底流压,Psc表征标准条件下压力,μ表征地下气体粘度,Z表征气体偏差因子,T表征地下温度,CA表征Dietz形状因子,eγ表征Euler常数,Rw表征井筒半径,以及A表征供给面积。
进一步,该饱和度数据包括含水饱和度和含气饱和度,所述子表征数据包括含水率和气液比。
进一步,该在所述得到目标流压数据和目标配产量数据之后,还包括:
利用所述气藏数值模拟模型对所述目标流压数据、所述目标配产量数据和在第二采集时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据,所述第二采集时刻为所述第一采集时刻之后的时刻;
利用所述井筒模型针对目标表征数据和所述目标流动表现方程进行处理,得到目标垂直提升表现方程,所述目标表征数据包括由第二瞬态井数据确定的目标平均气藏压力数据和由所述目标饱和度数据确定目标子表征数据;
联立所述目标流动表现方程和所述目标垂直提升表现方程,求解得到目标交点数据;以及
基于所述地面管网模型针对与终端分离器对应的参数和所述目标交点数据进行处理,得到指导流压数据和指导配产量数据。
另一方面,本说明书实施例还提供了一种能源井一体化模拟方法,包括,
根据接收到的能源井数值模拟请求,确定模拟阶段信息;
在确定所述模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;以及
在确定所述模拟阶段信息指示运行阶段的情况下,采用如上油气藏配产数据确定方法中任一所述的方法,确定目标流压数据和目标配产量数据。
进一步,在确定所述模拟阶段信息指示运行阶段的情况下,还包括采用如上所述的油气藏配产数据确定方法,确定指导流压数据和指导配产量数据。
进一步,该利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度进一步包括,
基于所述气藏数值模型针对所述瞬态井数据进行处理,得到井底流压、井底产量和井底温度;
利用所述井筒模型针对所述井底流压、所述井底产量和所述井底温度进行处理,得到井口压力、井口产量和井口温度;以及
利用所述地面管网模型针对所述井口压力、所述井口产量和所述井口温度进行处理,得到所述终端分离器的压力、所述终端分离器的产量和所述终端分离器的温度。
另一方面,本说明书实施例还提供了一种油气藏配产数据确定装置,包括,
第一处理单元,用于利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据,所述第一瞬态井数据与第一采集时刻相对应;
第二处理单元,用于利用井筒模型针对表征数据和所述流动表现方程进行处理,得到垂直提升表现方程,所述表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由所述饱和度数据确定子表征数据;
第一求解单元,用于联立所述流动表现方程和所述垂直提升表现方程,求解得到交点数据;以及
第三处理单元,用于基于地面管网模型针对与终端分离器对应的参数和所述交点数据进行处理,得到目标流压数据和目标配产量数据。
进一步,还包括,
第四处理单元,用于利用所述气藏数值模拟模型对所述目标流压数据、所述目标配产量数据和在第二采集时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据,所述第二采集时刻为所述第一采集时刻之后的时刻;
第五处理单元,用于利用所述井筒模型针对目标表征数据和所述目标流动表现方程进行处理,得到目标垂直提升表现方程,所述目标表征数据包括由第二瞬态井数据确定的目标平均气藏压力数据和由所述目标饱和度数据确定目标子表征数据;
第二求解单元,用于联立所述目标流动表现方程和所述目标垂直提升表现方程,求解得到目标交点数据;以及
第六处理单元,用于基于所述地面管网模型针对与终端分离器对应的参数和所述目标交点数据进行处理,得到指导流压数据和指导配产量数据。
另一方面,本说明书实施例还提供了一种能源井一体化模拟装置,包括,
第一确定单元,用于根据接收到的能源井数值模拟请求,确定模拟阶段信息;
第七处理单元,用于在确定所述模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;以及
第八处理单元,用于在确定所述模拟阶段信息指示运行阶段的情况下,采用如油气藏配产数据确定装置中任一所述的装置,确定目标流压数据和目标配产量数据。
进一步,还包括,
第九处理单元,用于在确定所述模拟阶段信息指示运行阶段的情况下,还包括采用如上所述的如油气藏配产数据确定装置,确定指导流压数据和指导配产量数据。
另一方面,本说明书实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述的方法。
另一方面,本说明书实施例还提供了一种计算机可读存储介质,其上存储有计算机指令,该计算机指令被处理器执行时实现上述的方法。
另一方面,本说明书实施例还提供了一种计算机程序产品,包括计算机程序/指令,所述计算机程序/指令被处理器执行时实现的方法。
利用本说明书实施例,在从能源井采集得到的第一瞬态井数据时,利用气藏数值模拟模型对该第一瞬态井数据进行处理,得到流动表现方程和饱和度数据;利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程,该表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由饱和度数据确定子表征数据;联立流动表现方程和垂直提升表现方程,求解得到交点数据;以及基于地面管网模型针对与终端分离器对应的参数和交点数据进行处理,得到目标流压数据和目标配产量数据。进而实现了仅可处理稳态数据的井筒模型和地面管网模型与处理瞬态井数据的气藏数值模拟模型的联合,从而实现了在气藏地下渗流、井筒多相管流和地面管网集输三者的基础上,确定合理的目标流压数据和目标配产量数据,进一步提高油气井产能。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1所示为本说明书实施例一种油气藏配产数据确定方法的实施系统示意图;
图2所示为本说明书实施例一种油气藏配产数据确定方法的流程图;
图3所示为本说明书实施例一种流动表现方程和饱和度数据确定方法的流程图;
图4A所示为本说明书另一实施例的一种油气藏配产数据确定方法的流程图;
图4B所示为本说明书另一实施例的一种油气藏配产数据确定方法的流程图;
图5所示为本说明书实施例的一种能源井一体化模拟方法的流程图;
图6所示为本说明书实施例的一种端分离器的参数确定方法的流程图;
图7所示为本说明书另一实施例的一种能源井一体化模拟方法的流程图;
图8A所示为本说明书实施例一种油气藏配产数据确定装置的结构示意图;
图8B所示为本说明书另一实施例的一种油气藏配产数据确定装置的结构示意图;
图9A所示为本说明书实施例一种能源井一体化模拟装置的结构示意图;
图9B所示为本说明书另一实施例的一种能源井一体化模拟装置的结构示意图;
图10为本说明书实施例一种计算机设备的结构示意图。
【附图标记说明】
101、采集终端;
102、服务器;
103、用户终端;
810、第一处理单元;
820、第二处理单元;
830、第一求解单元;
840、第三处理单元;
850、第四处理单元;
860、第五处理单元;
870、第二求解单元;
880、第六处理单元;
910、第一确定单元;
920、第七处理单元;
930、第八处理单元;
940、第九处理单元;
1002、计算机设备;
1004、处理设备;
1006、存储资源;
1008、驱动机构;
1010、输入/输出模块;
1012、输入设备;
1014、输出设备;
1016、呈现设备;
1018、图形用户接口;
1020、网络接口;
1022、通信链路;
1024、通信总线。
具体实施方式
下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本说明书一部分实施例,而不是全部的实施例。基于本说明书中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本说明书保护的范围。
需要说明的是,本说明书的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本说明书的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、装置、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1所示为本说明书实施例一种油气藏配产数据确定方法的实施系统示意图,可以包括:采集终端101、服务器102和用户终端,采集终端101、服务器102和用户终端之间通过网络进行通信,网络可以包括局域网(Local Area Network,简称为LAN)、广域网(WideArea Network,简称为WAN)、因特网或其组合,并连接至网站、用户设备(例如计算设备)和后端系统。采集终端101可以为能源井内的传感器,用于采集对应的第一瞬态井数据。在得到第一瞬态井数据之后,将该第一瞬态井数据发送至服务器102。服务器102在接收到该第一瞬态井数据时,利用气藏数值模拟模型对该第一瞬态井数据进行处理,得到流动表现方程和饱和度数据;利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程,该表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由饱和度数据确定子表征数据;联立流动表现方程和垂直提升表现方程,求解得到交点数据;以及基于地面管网模型针对与终端分离器对应的参数和交点数据进行处理,得到目标流压数据和目标配产量数据,并将该目标流压数据和目标配产量数据发送至用户终端103,以用于指导下一时刻或下一采集周期的能源采集工作。
进一步,在服务器102还可以利用气藏数值模拟模型对目标流压数据、目标配产量数据和在当前时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据;利用井筒模型针对目标表征数据和目标流动表现方程进行处理,得到目标垂直提升表现方程,目标表征数据包括由第二瞬态井数据确定的目标平均气藏压力数据和由目标饱和度数据确定目标子表征数据;联立目标流动表现方程和目标垂直提升表现方程,求解得到目标交点数据;以及基于地面管网模型针对与终端分离器对应的参数和目标交点数据进行处理,得到指导流压数据和指导配产量数据,并将该指导流压数据和指导配产量数据发送至用户终端103,以用于指导再下一时刻或再下一采集周期的能源采集工作。
进一步,用户终端103可以将用户输入的能源井数值模拟请求发送至服务器102。服务器102在接收到该能源井数值模拟请求时,根据接收到的能源井数值模拟请求,确定模拟阶段信息;在确定模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度,并发送至用户终端103;以及在确定模拟阶段信息指示运行阶段的情况下,采用如上方法,确定目标流压数据和目标配产量数据,或,指导流压数据和指导配产量数据,并发送至用户终端103。
可选地,服务器102和可以是云计算系统的节点(图中未显示),或者每个服务器102可以是单独的云计算系统,包括由网络互连并作为分布式处理系统工作的多台计算机。
在一个可选的实施例中,用户终端103可以包括电子设备不限于智能手机、采集设备、台式计算机、平板电脑、笔记本电脑、智能音箱、数字助理、增强现实(AR,AugmentedReality)/虚拟现实(VR,Virtual Reality)设备、智能可穿戴设备等类型的电子设备。可选的,电子设备上运行的操作系统可以包括但不限于安卓系统、IOS系统、Linux、Windows等。
此外,需要说明的是,图1所示的仅仅是本说明书提供的一种应用环境,在实际应用中,还可以包括多个采集终端101、多个服务器102和多个用户终端103,本说明书不做限制。
如图2所示为本说明书实施例一种油气藏配产数据确定方法的流程图。在本图中描述了油气藏配产数据确定过程,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的系统或装置产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行。具体的如图2所示,方法可以包括:
S210,利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据;
S220,利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程;
S230,联立流动表现方程和垂直提升表现方程,求解得到交点数据;
S240,基于地面管网模型针对与终端分离器对应的参数和交点数据进行处理,得到目标流压数据和目标配产量数据。
利用本说明书实施例,在从能源井采集得到的第一瞬态井数据时,利用气藏数值模拟模型对该第一瞬态井数据进行处理,得到流动表现方程和饱和度数据;利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程,该表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由饱和度数据确定子表征数据;联立流动表现方程和垂直提升表现方程,求解得到交点数据;以及基于地面管网模型针对与终端分离器对应的参数和交点数据进行处理,得到目标流压数据和目标配产量数据。进而实现了仅可处理稳态数据的井筒模型和地面管网模型与处理瞬态井数据的气藏数值模拟模型的联合,从而实现了在气藏地下渗流、井筒多相管流和地面管网集输三者的基础上,确定合理的目标流压数据和目标配产量数据,进一步提高油气井产能。
根据本说明书的一个实施例,气藏数值模拟模型为基于气藏三维网格地质模型和数值模拟模型进行运算,确定气藏从地下到井底的流动特征的模型,例如,油气藏建模仿真模型例如可以为基于油藏数值模拟软件(CMG)构建的数值模拟模型。CMG是由加拿大计算机模拟软件集团(Computer Modelling Group ltd.)开发的油藏数值模拟软件,用来模拟各种复杂地质条件下油气水在地下的运移、物理化学过程。气藏数值模拟模型是一个以时间为变量的瞬态模型,该气藏数值模拟模型处理的数据均为瞬态的数据。
井筒模型可以为井筒多相管流模型,例如为基于多物理场仿真软件(COMSOL)构建的确定井底到井口的生产特征的模型。该井筒模型为不以时间为变量的稳态模型。地面管网模型可以为地面集输管网模型,例如为基于黑油模拟软件(IMEX)构建的确定井口到地面管线的生产特征的模型。该地面管网模型为不以时间为变量的稳态模型。由于气藏数值模拟模型、井筒模型和地面管网模型三个模型处理的状态并不相同,所以三个模型无法相接,本说明书为了使得以上三个模型可以相接,做了如下改进。
基于采集终端针对至少一个能源井采集第一瞬态井数据,该第一瞬态井数据即为,需要输入气藏数值模拟模型的瞬态的数据,第一瞬态井数据与第一采集时刻相对应,该第一采集时刻例如可以为初始采集时刻。针对每个能源井的第一瞬态井数据进行稳态处理,得到稳态数据。稳态处理具体为针对第一瞬态井数据去时间维度的处理。进一步,利用气藏数值模拟模型针对该稳态数据进行处理,得到流动表现方程(流动表现关系曲线,IPR曲线)和饱和度数据。由此,气藏数值模拟模型针对稳态的数据进行处理,输出稳态的结果(IPR曲线和饱和度数据)。需要说明的是,该处理过程可参见现有技术,例如将第一瞬态井数据输入CMG软件,得到流动表现方程和饱和度数据。需要注意的是,若存在多个能源井,每个能源井均有与之对应的流动表现方程和饱和度数据。
基于该饱和度数据,确定表征数据。该表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由饱和度数据确定子表征数据。利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程(垂直提升表现曲线,VLP曲线)。需要说明的是,若存在多个能源井,每个能源井均有多个与之对应的垂直提升表现方程。
在得到流动表现方程和垂直提升表现方程之后,则针对每个能源井,均存在一个与之对应的流动表现方程和多个与之对应的垂直提升表现方程。
进一步,由于地面管网模型的输入为能源井井口的压力产量关系数据,而无论是气藏数值模拟模型还是井筒模型,均无法直接输出能源井井口的压力产量关系数据(P-Q曲线)。由此,为了得到能源井井口的压力产量关系数据,针对每个能源井,联立流动表现方程和垂直提升表现方程,得到解,将得到的解作为交点数据,该交点数据即表征能源井井口的压力产量关系数据。针对每个能源井,联立流动表现方程和垂直提升表现方程,得到解具体可以为,绘制一条IPR曲线和多条VLP曲线,则得到多个交点,每个交点均为一个解,从而得到多个交点数据。
每个能源井,在地面均有匹配的终端分离器。每个终端分离器均有与之对应的参数,该参数例如可以包括终端分离器的入口压力和最大处理产量。实际的地面管网通常是由多口能源井共同建立的一个多分支节点的管网拓扑模型,能源井和能源井之间会在地面管线节点处互相制约,并且终端分离器也会有入口压力和最大处理能力的约束,那么地面管网的模拟就成为了一个非线性问题的模拟过程。地面管网模型具体可以仅包括求解算法,也可以同时包括求解算法和优化算法,该求解算法例如可以为非线性模拟牛顿-拉夫森的求解算法,该优化算法例如可以为序列二次规划(SQP)优化算法。
由此,可以利用地面管网模型针对该参数(入口压力和最大处理产量)和交点数据进行处理,得到目标流压数据和目标配产量数据。需要说明的是,该利用地面管网模型针对该参数(入口压力和最大处理产量)和交点数据进行处理,得到目标流压数据和目标配产量数据的处理过程可以参见现有技术。具体地,可以采用IMEX软件针对该参数(入口压力和最大处理产量)和交点数据进行处理,得到目标流压数据和目标配产量数据。该确定的目标流压数据和目标配产量数据可以用于指导下一时刻或下一采集周期的能源采集工作(采油或采气),即将目标流压数据和目标配产量数据作为对应能源井的采集目标输入能源采集系统,能源采集系统基于该目标流压数据和目标配产量数据,调整采集措施,以达到该采集目标。
根据本说明书的另一个实施例,第一瞬态井数据力图可以包括气藏压力、渗透率、储层厚度、井筒半径、流体状态方程和供给面积。
根据本说明书的另一个实施例,饱和度数据包括含水饱和度和含气饱和度,子表征数据包括含水率和气液比。基于含水饱和度和含气饱和度,确定含水率和气液比的方法如现有的任一确定含水率和气液比的方法。
需要说明的是,子表征数据例如可以为气水比,第一瞬态数据还包括含水饱和度数据。与每个网格对应的函数含水饱和度数据结合气水相渗数据,得到子表征数据。用户预先将气水相渗数据输入气藏数值模拟模型。
图3所示为本说明书实施例一种流动表现方程和饱和度数据确定方法的流程图。在本图中描述了流动表现方程和饱和度数据的确定过程,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。具体的如图3所示,方法可以包括:
S311,采集与能源井的井距一半范围内的多个网格分别对应的第一瞬态数据;
S312,确定第一瞬态数据的平均值为稳态数据;
S313,利用气藏数值模拟模型针对稳态数据进行处理,得到流动表现方程和饱和度数据。
根据本说明书的另一个实施例,从基于至少一个能源井构建的气藏三维网格地质模型中,确定与待分析的能源井的井距一半范围内的多个网格,并采集与该多个网格对应的第一瞬态数据。气藏三维网格地质模型包括多个三维网格,和与每个三维网格对应的第一瞬态数据。基于此,针对每个能源井,确定与之对应的多个第一瞬态数据。
针对多个第一瞬态数据,求平均值,并将该平均值作为稳态数据。例如,第一瞬态数据包括气藏压力和渗透率,则稳态数据包括平均气藏压力和平均渗透率,该平均气藏压力为多个第一瞬态数据中的气藏压力的平均值,该平均渗透率为多个第一瞬态数据中的渗透率的平均值。
在确定了稳态数据之后,将该稳态数据输入气藏数值模拟模型,得到流动表现方程和饱和度数据。
根据本说明书的另一个实施例,气藏数值模拟模型包括拟稳态气井产能公式,该拟稳态气井产能公式如以下公式(1)所示。
其中,Q表征气产量,K表征平均渗透率,h表征平均储层厚度,Tsc表征标准条件下的温度,Pr表征平均气藏压力,Pwf表征井底流压,Psc表征标准条件下压力,μ表征地下气体粘度,Z表征气体偏差因子,T表征地下温度,GA表征Dietz形状因子,eγ表征Euler常数,Rw表征井筒半径,以及A表征供给面积。需要说明的是,表征数据包括由第一瞬态井数据确定的平均气藏压力数据即为公式(1)中的Pr。
图4A所示为本说明书另一实施例的一种油气藏配产数据确定方法的流程图;图4B所示为本说明书另一实施例的一种油气藏配产数据确定方法的流程图。在本图中描述了另一种油气藏配产数据确定的过程,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。具体的如图4A所示,方法可以包括:
S450,利用气藏数值模拟模型对目标流压数据、目标配产量数据和在第二采集时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据;
S460,利用井筒模型针对目标表征数据和目标流动表现方程进行处理,得到目标垂直提升表现方程;
S470,联立目标流动表现方程和目标垂直提升表现方程,求解得到目标交点数据;
S480,基于地面管网模型针对与终端分离器对应的参数和目标交点数据进行处理,得到指导流压数据和指导配产量数据。
利用本说明书实施例,在上一时刻(第一采集时刻)确定了目标流压数据和目标配产量数据之后,若以该目标流压数据和目标配产量数据作为当前时刻的采集目标,则由于未考虑当前时刻的用于采集能源的装置的情况,导致采集目标配置不合理。为了解决该技术问题,本说明书在对第二采集时刻确定目标配置时,在考虑上一时刻(第一采集时刻)确定的目标流压数据和目标配产量数据之外,还考虑当前时刻的用于采集能源的装置的情况,由此得到更合理的采集目标。
根据本说明书的另一个实施例,第二瞬态井数据与图2中的第一瞬态井数据相似,存在采集时刻不同的差异。同时,需要说明的时,第二采集时刻为第一采集时刻之后的时刻,第二采集时刻与第一采集时刻相隔一个采集步长。在确定与第三采集时刻(与第二采集时刻相隔一个采集步长的后一个采集时刻)对应的指导流压数据和指导配产量数据时,需要将与第二采集时刻对应的指导流压数据、指导配产量数据、在第三采集时刻能源井采集得到的第三瞬态井数据输入气藏数值模拟模型中,以确定与第三采集时刻对应的指导流压数据和指导配产量数据。
进一步,将目标流压数据、目标配产量数据第二瞬态井数据输入气藏数值模拟模型,得到目标流动表现方程和目标饱和度数据。需要说明的是,需要说明的是,该处理过程可参见现有技术,例如将第二瞬态井数据、目标流压数据和目标配产量数据输入CMG软件,得到流动表现方程和饱和度数据。进而执行以上S460~S480,具体执行以上S460~S480与执行S220~S240相似,在此不再赘述。
在确定指导流压数据和指导配产量数据之后,以该指导流压数据和指导配产量数据作为对应能源井当前时刻的采集目标输入能源采集系统,能源采集系统基于该指导流压数据和指导配产量数据,调整采集措施,以达到该采集目标。
具体地,如图4B所示,在确定油气藏的平均压力Pr0和含水饱和度Sw0稳态数据之后,利用气藏数值模拟模型针对该稳态数据进行处理,得到井控流动表现方程(井控IPR)、饱和度数据气(水比GWR0)。进一步,基于井筒模型,针对井控流动表现方程(井控IPR0,IPR0为由井底流压Pwf0和流体温度T构成的曲线)、饱和度数据气进行处理,得到垂直提升表现方程(VLP曲线)。绘制一条IPR曲线和多条VLP曲线,则得到多个交点数据(基于节点分析,确定的由井口压力Pt0和流体温度T0构成的交点数据)。进一步,基于地面管网模型,针对该多个交点数据和终端分离器对应的参数进行处理(基于N-P迭代求解和/或SQP迭代优化进行处理),得到目标流压数据和目标配产量数据。在下一个时刻(第二采集时刻)时,气藏数值模拟模型的压力场和饱和度场都已经发生了变化,此时每口能源井井控区的平均压力、平均饱和度也发生变化,需要将新的井控平均压力和平均饱和度提取出来,修正上一时刻(第一采集时刻)确定的IPR曲线,利用新的IPR曲线完成当前时刻的各个能源井的井底流压求解或配产优化。再次赋予气藏数值模拟模型,上一时刻(第一采集时刻)的目标流压数据和目标配产量数据作为气藏数值模拟模型对应能源井的下一个时刻(第二采集时刻)的定压或定产条件,开始下一个时刻(第二采集时刻)的模拟计算。
图5所示为本说明书实施例的一种能源井一体化模拟方法的流程图。在本图中描述了能源井一体化模拟的过程,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。具体的如图5所示,方法可以包括:
S510,根据接收到的能源井数值模拟请求,确定模拟阶段信息;
S520,在确定模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;
S530,在确定模拟阶段信息指示运行阶段的情况下,确定目标流压数据和目标配产量数据。
利用本说明书的实施例,在地面管网设计阶段,需要确定终端分离器的选购标准,且此时进行数值模拟的方法为线性联合模拟。而地面管网搭建完毕,终端分离器的参数已知的情况下,仅可通过非线性联合模拟才可以实现对沿程管线和能源井的数据进行模拟,从而形成一体化模拟仿真。若对两个阶段不加以区分,则会导致联合模拟的结果不合理,从而误导实际作业。
根据本说明书的另一个实施例,能源井数值模拟请求包括模拟阶段信息。该模拟阶段信息可以指示地面管网设计阶段和运行阶段。在指示地网管网设计阶段时,该模拟阶段信息例如可以为终端分离器的购买参数信息,例如,终端分离器的压力标识、终端分离器的产量标识和终端分离器的温度标识。在指示运行阶段时,该模拟阶段信息例如可以为配置的采集参数信息,例如,目标流压数据标识和目标配产量数据标识。
在确定模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理。具体地,能源井数值模拟请求例如可以包括能源井采集得到的瞬态井数据。例如还可以,在确定能源井数值模拟请求之后,基于该能源井数值模拟请求,采集对应的瞬态井数据。在确定瞬态井数据之后,获取气藏数值模拟模型、井筒模型和地面管网模型,并利用该气藏数值模拟模型、井筒模型和地面管网模型分别对瞬态井数据进行处理,得到与终端分离器的压力标识对应的终端分离器的压力、与终端分离器的产量标识对应的终端分离器的产量和与终端分离器的温度标识对应的终端分离器的温度。例如,进一步,还可以根据得到终端分离器的压力、终端分离器的产量和终端分离器的温度,确定与之对应的终端分离器标识,以推荐用户购买。
在确定模拟阶段信息指示运行阶段的情况下,确定瞬态井数据,并将该瞬态井数据作为S210~S240中的第一瞬态井数据进行处理,得到与目标流压数据标识对应的目标流压数据和与目标配产量数据标识对应的目标配产量数据。
图6所示为本说明书实施例的一种端分离器的参数确定方法的流程图。在本图中描述了端分离器的参数确定的过程,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。具体的如图6所示,方法可以包括:
S621,基于气藏数值模型针对瞬态井数据进行处理,得到井底流压、井底产量和井底温度;
S622,利用井筒模型针对井底流压、井底产量和井底温度进行处理,得到井口压力、井口产量和井口温度;
S623,利用地面管网模型针对井口压力、井口产量和井口温度进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度。
根据本说明书的另一个实施例,基于气藏数值模型针对与每个能源井对应的瞬态井数据进行处理,得到与每个能源井对应的井底流压、井底产量和井底温度。基于井筒模型针对与每个能源井对应的井底流压、井底产量和井底温度进行处理,得到与每个该能源井对应的井口压力、井口产量和井口温度。
基于地面管网模型针对与每个该能源井对应的井口压力、井口产量和井口温度进行处理,得到与每个该能源井对应的终端分离器的压力、终端分离器的产量和终端分离器的温度,进一步,基于该终端分离器的压力、终端分离器的产量和终端分离器的温度,由此实现了在地网管网设计阶段的线性联合模拟,以指导终端分离器的选购。
图7所示为本说明书另一实施例的一种能源井一体化模拟方法的流程图。在本图中描述了能源井一体化模拟的过程,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。具体的如图7所示,方法可以包括:
S731,在确定模拟阶段信息指示运行阶段的情况下,确定目标运行阶段;
S732,目标运行阶段是否指示首次模拟;
S733,在确定目标运行阶段指示首次模拟的情况下,确定目标流压数据和目标配产量数据;
S734,在确定目标运行阶段不是指示首次模拟的情况下,确定指导流压数据和指导配产量数据。
根据本说明书的另一个实施例,运行阶段例如可以指示首次模拟和非首次模拟两类信息。
在运行阶段为首次模拟的情况下,则表示当前非线性联合模拟没有上一时刻的目标流压数据和目标配产量数据,由此,基于S210~S240实现首次模拟,以得到目标流压数据和目标配产量数据。
在运行阶段为非首次模拟的情况下,则表示当前非线性联合模拟有上一时刻的目标流压数据和目标配产量数据,由此,基于S450~S480实现非首次模拟,以得到指导流压数据和指导配产量数据。进一步,还可以以该指导流压数据和指导配产量数据作为当前时刻的目标流压数据和目标配产量数据,以进行下一时刻的非首次模拟。
图8A所示为本说明书实施例一种油气藏配产数据确定装置的结构示意图。如图8A所示,包括,
第一处理单元810,用于利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据,所述第一瞬态井数据与第一采集时刻相对应;
第二处理单元820,用于利用井筒模型针对表征数据和流动表现方程进行处理,得到垂直提升表现方程,该表征数据包括由第一瞬态井数据确定的平均气藏压力数据和由饱和度数据确定子表征数据;
第一求解单元830,用于联立流动表现方程和垂直提升表现方程,求解得到交点数据;以及
第三处理单元840,用于基于地面管网模型针对与终端分离器对应的参数和交点数据进行处理,得到目标流压数据和目标配产量数据。
由于上述装置解决问题的原理与上述方法相似,因此上述装置的实施可以参见上述方法的实施,重复之处不再赘述。
图8B所示为本说明书另一实施例的一种油气藏配产数据确定装置的结构示意图。如图8B所示,包括,
第四处理单元850,用于利用气藏数值模拟模型对目标流压数据、目标配产量数据和在第二采集时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据,所述第二采集时刻为所述第一采集时刻之后的时刻;
第五处理单元860,用于利用井筒模型针对目标表征数据和目标流动表现方程进行处理,得到目标垂直提升表现方程,目标表征数据包括由第二瞬态井数据确定的目标平均气藏压力数据和由目标饱和度数据确定目标子表征数据;
第二求解单元870,用于联立目标流动表现方程和目标垂直提升表现方程,求解得到目标交点数据;以及
第六处理单元880,用于基于地面管网模型针对与终端分离器对应的参数和目标交点数据进行处理,得到指导流压数据和指导配产量数据。
由于上述装置解决问题的原理与上述方法相似,因此上述装置的实施可以参见上述方法的实施,重复之处不再赘述。
图9A所示为本说明书实施例一种能源井一体化模拟装置的结构示意图,如图8B所示,包括,
第一确定单元910,用于根据接收到的能源井数值模拟请求,确定模拟阶段信息;
第七处理单元920,用于在确定模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;以及
第八处理单元930,用于在确定模拟阶段信息指示运行阶段的情况下,确定目标流压数据和目标配产量数据。
由于上述装置解决问题的原理与上述方法相似,因此上述装置的实施可以参见上述方法的实施,重复之处不再赘述。
图9B所示为本说明书另一实施例的一种能源井一体化模拟装置的结构示意图。如图9B所示包括,
第一确定单元910,用于根据接收到的能源井数值模拟请求,确定模拟阶段信息;
第七处理单元920,用于在确定模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;
第八处理单元930,用于在确定模拟阶段信息指示运行阶段的情况下,确定目标流压数据和目标配产量数据;以及
第九处理单元940,用于在确定模拟阶段信息指示运行阶段的情况下,确定指导流压数据和指导配产量数据。
由于上述装置解决问题的原理与上述方法相似,因此上述装置的实施可以参见上述方法的实施,重复之处不再赘述。
如图10所示为本说明书实施例一种计算机设备的结构示意图,本说明书中的装置可以为本实施例中的计算机设备,执行上述本说明书的方法。计算机设备1002可以包括一个或多个处理设备1004,诸如一个或多个中央处理单元(CPU),每个处理单元可以实现一个或多个硬件线程。计算机设备1002还可以包括任何存储资源1006,其用于存储诸如代码、设置、数据等之类的任何种类的信息。非限制性的,比如,存储资源1006可以包括以下任一项或多种组合:任何类型的RAM,任何类型的ROM,闪存设备,硬盘,光盘等。更一般地,任何存储资源都可以使用任何技术来存储信息。进一步地,任何存储资源可以提供信息的易失性或非易失性保留。进一步地,任何存储资源可以表示计算机设备1002的固定或可移除部件。在一种情况下,当处理设备1004执行被存储在任何存储资源或存储资源的组合中的相关联的指令时,计算机设备1002可以执行相关联指令的任一操作。计算机设备1002还包括用于与任何存储资源交互的一个或多个驱动机构1008,诸如硬盘驱动机构、光盘驱动机构等。
计算机设备1002还可以包括输入/输出模块1010(I/O),其用于接收各种输入(经由输入设备1012)和用于提供各种输出(经由输出设备1014)。一个具体输出机构可以包括呈现设备1016和相关联的图形用户接口(GUI)1018。在其他实施例中,还可以不包括输入/输出模块1010(I/O)、输入设备1012以及输出设备1014,仅作为网络中的一台计算机设备。计算机设备1002还可以包括一个或多个网络接口1020,其用于经由一个或多个通信链路1022与其他设备交换数据。一个或多个通信总线1024将上文所描述的部件耦合在一起。
通信链路1022可以以任何方式实现,例如,通过局域网、广域网(例如,因特网)、点对点连接等、或其任何组合。通信链路1022可以包括由任何协议或协议组合支配的硬连线链路、无线链路、路由器、网关功能、名称服务器等的任何组合。
本说明书实施例还提供一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,计算机程序被处理器执行时实现上述方法。
本说明书实施例还提供一种计算机程序产品,计算机程序产品包括计算机程序,计算机程序被处理器执行时实现上述方法。
本领域内的技术人员应明白,本说明书的实施例可提供为方法、系统、或计算机程序产品。因此,本说明书可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本说明书可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书是参照根据本说明书实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上的具体实施例,对本说明书的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本说明书的具体实施例而已,并不用于限定本说明书的保护范围,凡在本说明书的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本说明书的保护范围之内。
Claims (16)
1.一种油气藏配产数据确定方法,其特征在于,包括:
利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据,所述第一瞬态井数据与第一采集时刻相对应;
利用井筒模型针对表征数据和所述流动表现方程进行处理,得到垂直提升表现方程,所述表征数据包括由所述第一瞬态井数据确定的平均气藏压力数据和由所述饱和度数据确定子表征数据;
联立所述流动表现方程和所述垂直提升表现方程,求解得到交点数据;以及
基于地面管网模型针对与终端分离器对应的参数和所述交点数据进行处理,得到目标流压数据和目标配产量数据。
2.根据权利要求1所述的方法,其特征在于,所述利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据包括:
采集与所述能源井的井距一半范围内的多个网格分别对应的第一瞬态数据;
确定所述第一瞬态数据的平均值为稳态数据;以及
利用所述气藏数值模拟模型针对所述稳态数据进行处理,得到所述流动表现方程和所述饱和度数据。
3.根据权利要求1所述的方法,其特征在于,所述第一瞬态井数据包括气藏压力、渗透率、储层厚度、井筒半径、流体状态方程和供给面积。
4.根据权利要求2所述的方法,其特征在于,所述气藏数值模拟模型包括拟稳态气井产能公式,所述拟稳态气井产能公式包括:
其中,Q表征气产量,K表征平均渗透率,h表征平均储层厚度,Tsc表征标准条件下的温度,Pr表征平均气藏压力,Pwf表征井底流压,Psc表征标准条件下压力,μ表征地下气体粘度,Z表征气体偏差因子,T表征地下温度,CA表征Dietz形状因子,eγ表征Euler常数,Rw表征井筒半径,以及A表征供给面积。
5.根据权利要求1所述的方法,其特征在于,所述饱和度数据包括含水饱和度和含气饱和度,所述子表征数据包括含水率和气液比。
6.根据权利要求1所述的方法,其特征在于,在所述得到目标流压数据和目标配产量数据之后,还包括:
利用所述气藏数值模拟模型对所述目标流压数据、所述目标配产量数据和在第二采集时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据,所述第二采集时刻为所述第一采集时刻之后的时刻;
利用所述井筒模型针对目标表征数据和所述目标流动表现方程进行处理,得到目标垂直提升表现方程,所述目标表征数据包括由所述第二瞬态井数据确定的目标平均气藏压力数据和由所述目标饱和度数据确定目标子表征数据;
联立所述目标流动表现方程和所述目标垂直提升表现方程,求解得到目标交点数据;以及
基于所述地面管网模型针对与终端分离器对应的参数和所述目标交点数据进行处理,得到指导流压数据和指导配产量数据。
7.一种能源井一体化模拟方法,其特征在于,包括,
根据接收到的能源井数值模拟请求,确定模拟阶段信息;
在确定所述模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;以及
在确定所述模拟阶段信息指示运行阶段的情况下,采用如权利要求1-6中任一所述的方法,确定目标流压数据和目标配产量数据。
8.根据权利要求7所述的方法,其特征在于,在确定所述模拟阶段信息指示运行阶段的情况下,还包括采用如权利要求6所述的方法,确定指导流压数据和指导配产量数据。
9.根据权利要求7所述的方法,其特征在于,所述利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度包括:
基于所述气藏数值模型针对所述瞬态井数据进行处理,得到井底流压、井底产量和井底温度;
利用所述井筒模型针对所述井底流压、所述井底产量和所述井底温度进行处理,得到井口压力、井口产量和井口温度;以及
利用所述地面管网模型针对所述井口压力、所述井口产量和所述井口温度进行处理,得到所述终端分离器的压力、所述终端分离器的产量和所述终端分离器的温度。
10.一种油气藏配产数据确定装置,其特征在于,包括:
第一处理单元,用于利用气藏数值模拟模型对能源井采集得到的第一瞬态井数据进行处理,得到流动表现方程和饱和度数据,所述第一瞬态井数据与第一采集时刻相对应;
第二处理单元,用于利用井筒模型针对表征数据和所述流动表现方程进行处理,得到垂直提升表现方程,所述表征数据包括由所述第一瞬态井数据确定的平均气藏压力数据和由所述饱和度数据确定子表征数据;
第一求解单元,用于联立所述流动表现方程和所述垂直提升表现方程,求解得到交点数据;以及
第三处理单元,用于基于地面管网模型针对与终端分离器对应的参数和所述交点数据进行处理,得到目标流压数据和目标配产量数据。
11.根据权利要求10所述的装置,其特征在于,还包括:
第四处理单元,用于利用所述气藏数值模拟模型对所述目标流压数据、所述目标配产量数据和在第二采集时刻能源井采集得到的第二瞬态井数据进行处理,得到目标流动表现方程和目标饱和度数据,所述第二采集时刻为所述第一采集时刻之后的时刻;
第五处理单元,用于利用所述井筒模型针对目标表征数据和所述目标流动表现方程进行处理,得到目标垂直提升表现方程,所述目标表征数据包括由所述第二瞬态井数据确定的目标平均气藏压力数据和由所述目标饱和度数据确定目标子表征数据;
第二求解单元,用于联立所述目标流动表现方程和所述目标垂直提升表现方程,求解得到目标交点数据;以及
第六处理单元,用于基于所述地面管网模型针对与终端分离器对应的参数和所述目标交点数据进行处理,得到指导流压数据和指导配产量数据。
12.一种能源井一体化模拟装置,其特征在于,还包括:
第一确定单元,用于根据接收到的能源井数值模拟请求,确定模拟阶段信息;
第七处理单元,用于在确定所述模拟阶段信息指示地面管网设计阶段的情况下,利用气藏数值模拟模型、井筒模型和地面管网模型针对能源井采集得到的瞬态井数据进行处理,得到终端分离器的压力、终端分离器的产量和终端分离器的温度;以及
第八处理单元,用于在确定所述模拟阶段信息指示运行阶段的情况下,采用如权利要求10-11中任一所述的装置,确定目标流压数据和目标配产量数据。
13.根据权利要求12所述的装置,其特征在于,还包括:
第九处理单元,用于在确定所述模拟阶段信息指示运行阶段的情况下,还包括采用如权利要求11所述的装置,确定指导流压数据和指导配产量数据。
14.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1-9中任一项的方法。
15.一种计算机可读存储介质,其特征在于,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器运行时执行上述权利要求1-9任一项的方法。
16.一种计算机程序产品,包括计算机程序/指令,其特征在于,所述计算机程序/指令被处理器执行时实现根据权利要求1-9中任一项的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310436435.XA CN116467967A (zh) | 2023-04-21 | 2023-04-21 | 一种油气藏配产数据确定和能源井一体化模拟方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310436435.XA CN116467967A (zh) | 2023-04-21 | 2023-04-21 | 一种油气藏配产数据确定和能源井一体化模拟方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116467967A true CN116467967A (zh) | 2023-07-21 |
Family
ID=87180366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310436435.XA Pending CN116467967A (zh) | 2023-04-21 | 2023-04-21 | 一种油气藏配产数据确定和能源井一体化模拟方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116467967A (zh) |
-
2023
- 2023-04-21 CN CN202310436435.XA patent/CN116467967A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2711167C (en) | Executing a utility in a distributed computing system based on an integrated model | |
GB2521545B (en) | Analysis of enhanced oil recovery (EOR) processes for naturally-fractured reservoirs | |
CA2704230C (en) | Systems and methods for the determination of active constraints in a network using slack variables | |
CN110096263A (zh) | Web前端代码自动生成方法及装置 | |
US10152041B2 (en) | Method and apparatus for enabling model driven navigation | |
CA2871407C (en) | System and method for reservoir simulation using on-demand data | |
AU2014398210B2 (en) | Multi-stage linear solution for implicit reservoir simulation | |
US10970434B2 (en) | Variable equivalency on connection in a process simulation | |
US8311743B2 (en) | Method for generating an estimation of incremental recovery from a selected enhanced oil recovery (EOR) process | |
CN116882218B (zh) | 一种油藏数值模拟方法、装置、计算机设备及存储介质 | |
CN105653746B (zh) | 基于嵌入式离散裂缝模型的压裂井建模及模拟方法 | |
US10719299B2 (en) | Physical connection of network components in a graphical computer environment with one-way signal propagation | |
CN116467967A (zh) | 一种油气藏配产数据确定和能源井一体化模拟方法及装置 | |
CN110309597B (zh) | 基于构型界面的阻流带模型确定方法、装置及存储介质 | |
CN111852466B (zh) | 一种用于页岩气井规模配产与管网运行优化的方法 | |
WO2017099716A1 (en) | Visual flow analyzer for exploratory hypotheses in upstream oil and gas process data | |
CN116167246A (zh) | 一种能源井配产数据确定方法、装置及计算机设备 | |
Ignatyev et al. | The features of building the Integrated Model for development of two gas-condensate formations of Urengoyskoe field | |
CN117217131B (zh) | 一种油藏数值模拟方法、装置、设备和存储介质 | |
WO2017048266A1 (en) | Fluid flow network simulation methods and systems employing two-phase envelopes with interpolated values | |
Jia et al. | An Effective Numerical Model for Fracture-Stimulated Condensate Reservoir Production History Matching, Surveillance, and Prediction | |
CN117034717A (zh) | 一种用于陆相页岩油高效提产的单簇点压裂改造方法 | |
Wilson | Use of Emulator Methodology for Uncertainty-Reduction Quantification | |
CN114297847A (zh) | 一种用于油气开采的压裂参数优选方法、装置及设备 | |
CN116737582A (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 |