CN113610292A - 一种致密气藏压裂水平井动态储量的确定方法和装置 - Google Patents

一种致密气藏压裂水平井动态储量的确定方法和装置 Download PDF

Info

Publication number
CN113610292A
CN113610292A CN202110885050.2A CN202110885050A CN113610292A CN 113610292 A CN113610292 A CN 113610292A CN 202110885050 A CN202110885050 A CN 202110885050A CN 113610292 A CN113610292 A CN 113610292A
Authority
CN
China
Prior art keywords
pressure
dimensionless
horizontal well
denotes
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.)
Pending
Application number
CN202110885050.2A
Other languages
English (en)
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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202110885050.2A priority Critical patent/CN113610292A/zh
Publication of CN113610292A publication Critical patent/CN113610292A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Animal Husbandry (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Agronomy & Crop Science (AREA)
  • Development Economics (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请提供了一种致密气藏压裂水平井动态储量的确定方法和装置,其中,该方法包括:根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;将无因次压力敏感系数和无量纲压力代入渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;获取目标生产井的产量,根据目标生产井的产量计算目标生产井的井底压力;根据井底压力和基于应力敏感的压裂水平井动态储量预测模型,计算目标生产井基于压力敏感的动态储量。通过上述方案解决了现有的致密气藏压裂水平井动态储量预测准确度较低的问题,达到了提升致密气藏压裂水平井动态储量预测准确度的技术效果。

Description

一种致密气藏压裂水平井动态储量的确定方法和装置
技术领域
本申请属于地质勘探技术领域,尤其涉及一种致密气藏压裂水平井动态储量的确定方法和装置。
背景技术
随着常规油气勘探资源越来越少,且勘探开发理论与技术的不断创新,非常规油气资源逐渐成为油气勘探开发增储上产的重要途径。其中,致密气的勘探开发是非常规油气资源勘探开发的主要方向。
然而,致密气储层因为低孔(孔隙度小于10%)、低渗(渗透率低于0.1mD)的特点,导致对致密气储层的天然能量开发不具备工业化规模开发的基础,一般所采用的开发方式是“水平井+水力压裂”。
动态储量是油气田开发方案设计参考的重要指标之一,动态储量对水平井部署、生产井合理配产、井网调整等开发方案具有重要的指导作用。然而,如何对致密气藏压裂水平井动态储量进行准确有效的预测,目前尚未提出有效的解决方案。
发明内容
本申请目的在于提供一种致密气藏压裂水平井动态储量的确定方法和装置,可以实现对致密气藏压力水平动态储量的准确预测。
本申请提供一种致密气藏压裂水平井动态储量的确定方法和装置是这样实现的:
一种致密气藏压裂水平井动态储量的确定方法,所述方法包括:
根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
在一个实施方式中,根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程,包括:
建立如下的渗流方程:
Figure BDA0003193713060000021
其中,γ表示压力敏感系数,r表示无量纲距离,pi表示气藏原始压力,p表示压力,φ表示孔隙度,μi表示初始时刻气体粘度,Cgi表示压缩系数,i表示初始状态,ki表示初始渗透率。
在一个实施方式中,对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力,包括:
按照如下公式对压力敏感系数进行无因次化处理得到无因次压力敏感系数:
Figure BDA0003193713060000022
其中,γD表示无因次压力敏感系数,q表示产量,psc表示外边界压力,T表示温度,h表示地层厚度,Tsc表示内边界处温度;
按照如下公式对所述压力参数进行无量纲处理得到无量纲压力:
Figure BDA0003193713060000023
其中,pD表示无量纲压力。
在一个实施方式中,所述基于应力敏感的压裂水平井动态储量预测模型表示为:
Figure BDA0003193713060000024
其中,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数,tD表示无因次时间。
在一个实施方式中,在得到基于应力敏感的压裂水平井动态储量预测模型之后还包括:
为所述基于应力敏感的压裂水平井动态储量预测模型建立初始条件、内边界条件和外边界条件。
在一个实施方式中,所述初始条件表示为:
pD(rD,0)=0
所述内边界条件表示为:
Figure BDA0003193713060000031
所述外边界条件表示为:
Figure BDA0003193713060000032
其中,reD表示外边界的无量纲距离,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数。
一种致密气藏压裂水平井动态储量的确定装置,包括:
建立模块,用于根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
处理模块,用于对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
代入模块,用于将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
获取模块,用于获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
计算模块,用于根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
在一个实施方式中,所述建立模块具体用于建立如下的渗流方程:
Figure BDA0003193713060000033
其中,γ表示压力敏感系数,r表示无量纲距离,pi表示气藏原始压力,p表示压力,φ表示孔隙度,μi表示初始时刻气体粘度,Cgi表示压缩系数,i表示初始状态,ki表示初始渗透率。
一种电子设备,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现上述方法的步骤。
一种计算机可读存储介质,其上存储有计算机指令,所述指令被执行时实现上述方法的步骤。
本申请提供的致密气藏压裂水平井动态储量的确定方法和装置,通过引入压力敏感系数来建立渗流方程,并对压力敏感系数进行无因次处理,得到基于应力敏感的压裂水平井动态储量预测模型,从而得到了考虑压力敏感系数的动态储量预测模型,使得计算得到的动态储量更为准确,从而解决了现有的致密气藏压裂水平井动态储量预测准确度较低的问题,达到了提升致密气藏压裂水平井动态储量预测准确度的技术效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的致密气藏压裂水平井动态储量的确定方法一种实施例的方法流程图;
图2是本申请提供的压裂水平井的示意图;
图3是本申请提供的苏9-11H井的生产状况图;
图4是本申请提供的一种致密气藏压裂水平井动态储量的确定方法的电子设备的硬件结构框图;
图5是本申请提供的致密气藏压裂水平井动态储量的确定装置一种实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请提供的致密气藏压裂水平井动态储量的确定方法一种实施例的方法流程图。虽然本申请提供了如下述实施例或附图所示的方法操作步骤或装置结构,但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中,这些步骤的执行顺序或装置的模块结构不限于本申请实施例描述及附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置或终端产品应用时,可以按照实施例或者附图所示的方法或模块结构连接进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至分布式处理环境)。
具体的,如图1所示,上述的致密气藏压裂水平井动态储量的确定方法可以包括如下步骤:
步骤101:根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
步骤102:对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
步骤103:将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
步骤104:获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
步骤105:根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
在上例中,通过引入压力敏感系数来建立渗流方程,并对压力敏感系数进行无因次处理,得到基于应力敏感的压裂水平井动态储量预测模型,从而得到了考虑压力敏感系数的动态储量预测模型,使得计算得到的动态储量更为准确,从而解决了现有的致密气藏压裂水平井动态储量预测准确度较低的问题,达到了提升了致密气藏压裂水平井动态储量预测准确度的技术效果。
在本例中,根据致密储层渗流机理,可以建立如下的基于应力敏感的致密气藏的渗流方程:
Figure BDA0003193713060000051
其中,γ表示压力敏感系数,r表示无量纲距离,p(r,t=0)=pi,其中,t表示时间,pi表示气藏原始压力,p表示压力,φ表示孔隙度,μi表示初始时刻气体粘度,Cgi表示压缩系数,i表示初始状态,ki表示初始渗透率。
针对上述的渗流方程,可以对其中的压力敏感系数进行无因次化处理得到无因次压力敏感系数,对压力参数进行无量纲处理得到无量纲压力。
具体的,可以按照如下公式对压力敏感系数进行无因次化处理得到无因次压力敏感系数:
Figure BDA0003193713060000052
其中,γD表示无因次压力敏感系数,q表示产量,psc表示外边界压力,T表示温度,h表示地层厚度,Tsc表示内边界处温度;
可以按照如下公式对所述压力参数进行无量纲处理得到无量纲压力:
Figure BDA0003193713060000061
其中,pD表示无量纲压力;
将上述得到的无因次压力敏感系数和无量纲压力,代入上述的渗流方程,可以得到如下的基于应力敏感的压裂水平井动态储量预测模型:
Figure BDA0003193713060000062
其中,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数,tD表示无因次时间。
针对上述的基于应力敏感的压裂水平井动态储量预测模型,可以建立初始条件、内边界条件和外边界条件。
其中,初始条件可以表示为:
pD(rD,0)=0
内边界条件可以表示为:
Figure BDA0003193713060000063
外边界条件可以表示为:
Figure BDA0003193713060000064
其中,reD表示外边界的无量纲距离,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数。
在根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程的时候,可以根据实际区块的地质条件,假定在均质气藏中心内有口压裂水平井,模型的边界条件为封闭圆形,裂缝方向为垂直横向,水平段的长度为L,水平段方向处在Y轴方向,而裂缝面与水平段方向垂直。模型假设水平井压裂裂缝为有限导流,每条横切水平井的裂缝具有相同属性,即裂缝的长度、宽度、高度和导流能力都保持一致,流体通过裂缝流入水平井筒,且总流量保持不变。
基于此,可以假设地层满足如下条件:气藏原始压力为pi,地层厚度为h;气藏各向异性,即,在x、y方向基质渗透率分别为kx、ky,基质孔隙度为Φm,裂缝渗透率kf,水平井段长度为L,气藏顶底和圆形边界均为不渗透边界;气体的渗流过程符合达西定律,不考虑毛管力和重力的影响;假定裂缝半长为Lf,其方向与井筒方向垂直,完全穿透储层;气井以一定产量q生产。可以在这种假定前提下,建立基于应力敏感的致密气藏的渗流方程的时候。
下面结合一个具体实施例对上述方法进行说明,然而值得注意的是,该具体实施例仅是为了更好地说明本申请,并不构成对本申请的不当限定。
针对无法对致密气藏压裂水平井动态储量进行准确预测的问题,在本例中提供了一种基于应力敏感效应的致密气藏压裂水平井动态储量预测方法,可以包括如下步骤:
步骤1:基于致密储层渗流机理,建立考虑应力敏感压裂水平井动态储量预测模型;
具体的:
1)根据实际区块的地质条件,假定在均质气藏中心内有口压裂水平井,模型的边界条件为封闭圆形,裂缝方向为垂直横向,水平段的长度为L,如图2所示,水平段方向处在Y轴方向,而裂缝面与水平段方向垂直。模型假设水平井压裂裂缝为有限导流,每条横切水平井的裂缝具有相同属性,即裂缝的长度、宽度、高度和导流能力都保持一致,流体通过裂缝流入水平井筒,且总流量保持不变。
假设条件为:
1.1气藏原始压力为pi,地层厚度为h;
1.2气藏各向异性,即,在x、y方向基质渗透率分别为kx、ky,基质孔隙度为Φm,裂缝渗透率kf,水平井段长度为L,气藏顶底和圆形边界均为不渗透边界;
1.3气体的渗流过程符合达西定律,不考虑毛管力和重力的影响;
1.4假定裂缝半长为Lf,其方向与井筒方向垂直,完全穿透储层;
1.5气井以一定产量q生产。
无量纲拟压力:
Figure BDA0003193713060000071
Figure BDA0003193713060000072
其中,pD表示无因次压力,μ表示气体粘度,单位为mPa·s,T表示温度,单位为℃,k表示渗透率,单位为10-3μm3,z表示气体偏差因子。
无量纲物质平衡拟时间:
Figure BDA0003193713060000081
其中,tD表示无因次时间;Ct表示综合性压缩系数,单位为MPa-1;φ表示孔隙度;km表示基质渗透率,单位为10-3μm3;tca表示为生产过程中的某一时刻,单位为s;μgi表示初始时刻气体粘度,单位为mPa·s。
无量纲距离:
Figure BDA0003193713060000082
无因次井筒半径:
Figure BDA0003193713060000083
其中,rw表示井筒,单位为mm。
无量纲水平井半长:
Figure BDA0003193713060000084
其中,xD表示x方向的无因次长度,yD表示y方向的无因次长度。
2)建立连续性方程、状态方程、达西定律方程、辅助方程和边界、初始条件。
①连续性方程:
Figure BDA0003193713060000085
②状态方程:
pM=ZρRT (7)
其中,M表示气体质量,kg;z表示气体偏差因子;ρ表示气体密度,g/L;R表示气体常数,T表示温度,单位为℃。
③达西定律:
Figure BDA0003193713060000086
④孔隙的压缩性:
Figure BDA0003193713060000091
⑤气体等温压缩系数:
Figure BDA0003193713060000092
初始条件
p(r,t=0)=pi (11)外边界条件:
Figure BDA0003193713060000093
内边界条件:
Figure BDA0003193713060000094
由于致密气藏的渗透率比较小,认为在基质中流动为拟稳态流动,渗流方程为:
Figure BDA0003193713060000095
其中,γ表示应力敏感系数。
将应力敏感系数无因次化处理:
Figure BDA0003193713060000096
无量纲压力为:
Figure BDA0003193713060000097
将上面无量纲压力和应力敏感系数引入渗流方程,将渗流方程(14)转化为以下形式:
Figure BDA0003193713060000098
初始条件:
pD(rD,0)=0 (16)
内边界条件
Figure BDA0003193713060000101
外边界条件
Figure BDA0003193713060000102
步骤2:选取典型生产井,利用产量计算出井底压力,根据井底压力,对流压和产量数据进行变换,进行图版拟合,计算出考虑应力敏感的动态储量和不考虑应力敏感的动态储量。
如图3所示,为苏9-11H井的生产状况图,根据产量计算出井底压力,对流压和产量数据进行变换,进行图版拟合可以得到拟合图,不考虑应力敏感动态储量为41.06×106m3,考虑压力敏感拟合得到的动态储量35.42×106m3,应力敏感造成储量损失为13.73%。
本申请上述实施例所提供的方法实施例可以在移动终端、计算机终端或者类似的运算装置中执行。以运行在电子设备上为例,图4是本申请提供的一种致密气藏压裂水平井动态储量的确定方法的电子设备的硬件结构框图。如图4所示,电子设备10可以包括一个或多个(图中仅示出一个)处理器02(处理器02可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)、用于存储数据的存储器04、以及用于通信功能的传输模块06。本领域普通技术人员可以理解,图4所示的结构仅为示意,其并不对上述电子装置的结构造成限定。例如,电子设备10还可包括比图4中所示更多或者更少的组件,或者具有与图4所示不同的配置。
存储器04可用于存储应用软件的软件程序以及模块,如本申请实施例中的致密气藏压裂水平井动态储量的确定方法对应的程序指令/模块,处理器02通过运行存储在存储器04内的软件程序以及模块,从而执行各种功能应用以及数据处理,即实现上述的应用程序的致密气藏压裂水平井动态储量的确定方法。存储器04可包括高速随机存储器,还可包括非易失性存储器,如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中,存储器04可进一步包括相对于处理器02远程设置的存储器,这些远程存储器可以通过网络连接至电子设备10。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
传输模块06用于经由一个网络接收或者发送数据。上述的网络具体实例可包括电子设备10的通信供应商提供的无线网络。在一个实例中,传输模块06包括一个网络适配器(Network Interface Controller,NIC),其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中,传输模块06可以为射频(Radio Frequency,RF)模块,其用于通过无线方式与互联网进行通讯。
在软件层面,上述装置可以如图5所示,包括:
建立模块501,用于根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
处理模块502,用于对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
代入模块503,用于将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
获取模块504,用于获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
计算模块505,用于根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
在一个实施方式中,上述建立模块501具体可以用于建立如下的渗流方程:
Figure BDA0003193713060000111
其中,γ表示压力敏感系数,r表示无量纲距离,pi表示气藏原始压力,p表示压力,φ表示孔隙度,μi表示初始时刻气体粘度,Cgi表示压缩系数,i表示初始状态,ki表示初始渗透率。
在一个实施方式中,对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力,可以包括:
按照如下公式对压力敏感系数进行无因次化处理得到无因次压力敏感系数:
Figure BDA0003193713060000112
其中,γD表示无因次压力敏感系数,q表示产量,psc表示外边界压力,T表示温度,h表示地层厚度,Tsc表示内边界处温度;
按照如下公式对所述压力参数进行无量纲处理得到无量纲压力:
Figure BDA0003193713060000113
其中,pD表示无量纲压力。
在一个实施方式中,所述基于应力敏感的压裂水平井动态储量预测模型可以表示为:
Figure BDA0003193713060000121
其中,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数,tD表示无因次时间。
在一个实施方式中,在得到基于应力敏感的压裂水平井动态储量预测模型之后还可以包括:
为所述基于应力敏感的压裂水平井动态储量预测模型建立初始条件、内边界条件和外边界条件。
在一个实施方式中,所述初始条件可以表示为:
pD(rD,0)=0
所述内边界条件可以表示为:
Figure BDA0003193713060000122
所述外边界条件可以表示为:
Figure BDA0003193713060000123
其中,reD表示外边界的无量纲距离,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数。
本申请的实施例还提供能够实现上述实施例中的致密气藏压裂水平井动态储量的确定方法中全部步骤的一种电子设备的具体实施方式,所述电子设备具体包括如下内容:处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线;其中,所述处理器、存储器、通信接口通过所述总线完成相互间的通信;所述处理器用于调用所述存储器中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的致密气藏压裂水平井动态储量的确定方法中的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
步骤1:根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
步骤2:对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
步骤3:将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
步骤4:获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
步骤5:根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
从上述描述可知,本申请实施例通过引入压力敏感系数来建立渗流方程,并对压力敏感系数进行无因次处理,得到基于应力敏感的压裂水平井动态储量预测模型,从而得到了考虑压力敏感系数的动态储量预测模型,使得计算得到的动态储量更为准确,从而解决了现有的致密气藏压裂水平井动态储量预测准确度较低的问题,达到了提升了致密气藏压裂水平井动态储量预测准确度的技术效果。
本申请的实施例还提供能够实现上述实施例中的致密气藏压裂水平井动态储量的确定方法中全部步骤的一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的致密气藏压裂水平井动态储量的确定方法的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
步骤1:根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
步骤2:对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
步骤3:将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
步骤4:获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
步骤5:根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
从上述描述可知,本申请实施例通过引入压力敏感系数来建立渗流方程,并对压力敏感系数进行无因次处理,得到基于应力敏感的压裂水平井动态储量预测模型,从而得到了考虑压力敏感系数的动态储量预测模型,使得计算得到的动态储量更为准确,从而解决了现有的致密气藏压裂水平井动态储量预测准确度较低的问题,达到了提升了致密气藏压裂水平井动态储量预测准确度的技术效果。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已,并不用于限制本说明书实施例。对于本领域技术人员来说,本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本说明书实施例的权利要求范围之内。

Claims (10)

1.一种致密气藏压裂水平井动态储量的确定方法,其特征在于,所述方法包括:
根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
2.根据权利要求1所述的方法,其特征在于,根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程,包括:
建立如下的渗流方程:
Figure FDA0003193713050000011
其中,γ表示压力敏感系数,r表示无量纲距离,pi表示气藏原始压力,p表示压力,φ表示孔隙度,μi表示初始时刻气体粘度,Cgi表示压缩系数,i表示初始状态,ki表示初始渗透率。
3.根据权利要求2所述的方法,其特征在于,对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力,包括:
按照如下公式对压力敏感系数进行无因次化处理得到无因次压力敏感系数:
Figure FDA0003193713050000012
其中,γD表示无因次压力敏感系数,q表示产量,psc表示外边界压力,T表示温度,h表示地层厚度,Tsc表示内边界处温度;
按照如下公式对所述压力参数进行无量纲处理得到无量纲压力:
Figure FDA0003193713050000013
其中,pD表示无量纲压力。
4.根据权利要求3所述的方法,其特征在于,所述基于应力敏感的压裂水平井动态储量预测模型表示为:
Figure FDA0003193713050000021
其中,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数,tD表示无因次时间。
5.根据权利要求4所述的方法,其特征在于,在得到基于应力敏感的压裂水平井动态储量预测模型之后还包括:
为所述基于应力敏感的压裂水平井动态储量预测模型建立初始条件、内边界条件和外边界条件。
6.根据权利要求5所述的方法,其特征在于:
所述初始条件表示为:
pD(rD,0)=0
所述内边界条件表示为:
Figure FDA0003193713050000022
所述外边界条件表示为:
Figure FDA0003193713050000023
其中,reD表示外边界的无量纲距离,pD表示无量纲压力,rD表示无量纲距离,γD表示无因次压力敏感系数。
7.一种致密气藏压裂水平井动态储量的确定装置,其特征在于,包括:
建立模块,用于根据致密储层渗流机理,建立基于应力敏感的致密气藏的渗流方程;
处理模块,用于对压力敏感系数进行无因次化处理得到无因次压力敏感系数,对所述压力参数进行无量纲处理得到无量纲压力;
代入模块,用于将所述无因次压力敏感系数和所述无量纲压力代入所述渗流方程,得到基于应力敏感的压裂水平井动态储量预测模型;
获取模块,用于获取目标生产井的产量,根据所述目标生产井的产量计算所述目标生产井的井底压力;
计算模块,用于根据所述井底压力和所述基于应力敏感的压裂水平井动态储量预测模型,计算所述目标生产井基于压力敏感的动态储量。
8.根据权利要求7所述的装置,其特征在于,所述建立模块具体用于建立如下的渗流方程:
Figure FDA0003193713050000031
其中,γ表示压力敏感系数,r表示无量纲距离,pi表示气藏原始压力,p表示压力,φ表示孔隙度,μi表示初始时刻气体粘度,Cgi表示压缩系数,i表示初始状态,ki表示初始渗透率。
9.一种电子设备,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现权利要求1至6中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机指令,所述指令被执行时实现权利要求1至6中任一项所述方法的步骤。
CN202110885050.2A 2021-08-03 2021-08-03 一种致密气藏压裂水平井动态储量的确定方法和装置 Pending CN113610292A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110885050.2A CN113610292A (zh) 2021-08-03 2021-08-03 一种致密气藏压裂水平井动态储量的确定方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110885050.2A CN113610292A (zh) 2021-08-03 2021-08-03 一种致密气藏压裂水平井动态储量的确定方法和装置

Publications (1)

Publication Number Publication Date
CN113610292A true CN113610292A (zh) 2021-11-05

Family

ID=78339192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110885050.2A Pending CN113610292A (zh) 2021-08-03 2021-08-03 一种致密气藏压裂水平井动态储量的确定方法和装置

Country Status (1)

Country Link
CN (1) CN113610292A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114169158A (zh) * 2021-11-30 2022-03-11 中国石油大学(北京) 一种致密砂岩气藏多层合采产能预测方法、装置及设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106484933A (zh) * 2015-08-31 2017-03-08 中国石油化工股份有限公司 一种用于确定页岩气井井控动态储量的方法及系统
CN107609986A (zh) * 2017-08-25 2018-01-19 中国石油天然气股份有限公司 气井动态储量的确定方法和装置
CN108266184A (zh) * 2018-01-14 2018-07-10 常州大学 一种基于生产数据的气井单井控制储量估算方法
CN110318741A (zh) * 2018-03-28 2019-10-11 中国石油化工股份有限公司 一种异常高压、高含水致密低渗气藏剩余气分布描述方法
CN111582532A (zh) * 2019-02-19 2020-08-25 中国石油天然气股份有限公司 应力敏感性油藏水平井产液能力预测方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106484933A (zh) * 2015-08-31 2017-03-08 中国石油化工股份有限公司 一种用于确定页岩气井井控动态储量的方法及系统
CN107609986A (zh) * 2017-08-25 2018-01-19 中国石油天然气股份有限公司 气井动态储量的确定方法和装置
CN108266184A (zh) * 2018-01-14 2018-07-10 常州大学 一种基于生产数据的气井单井控制储量估算方法
CN110318741A (zh) * 2018-03-28 2019-10-11 中国石油化工股份有限公司 一种异常高压、高含水致密低渗气藏剩余气分布描述方法
CN111582532A (zh) * 2019-02-19 2020-08-25 中国石油天然气股份有限公司 应力敏感性油藏水平井产液能力预测方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王少军: "《累积产量图版法预测致密气藏动态储量》", 石油天然气学报, vol. 35, no. 5, pages 83 - 87 *
田冷: "《低渗透气田砂岩储层应力敏感试井模型研究"》", 石油钻探技术, vol. 35, no. 6, pages 89 - 92 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114169158A (zh) * 2021-11-30 2022-03-11 中国石油大学(北京) 一种致密砂岩气藏多层合采产能预测方法、装置及设备
CN114169158B (zh) * 2021-11-30 2022-11-22 中国石油大学(北京) 一种致密砂岩气藏多层合采产能预测方法、装置及设备

Similar Documents

Publication Publication Date Title
CN111101929B (zh) 一种计算油气藏平均地层压力的方法、装置及系统
CN107366534B (zh) 粗化渗透率的确定方法和装置
CN110414723B (zh) 基于微震事件的裂缝油气藏历史拟合的方法、装置及系统
CN105760564A (zh) 分析油层套管失效的方法和装置
CN111507537B (zh) 一种水驱气藏的储量和水侵量的预测方法及装置
CN107806916B (zh) 储层改造体积的确定方法和装置
CN113610292A (zh) 一种致密气藏压裂水平井动态储量的确定方法和装置
CN108227032B (zh) 一种地应力的确定方法及装置
CN109681184B (zh) 页岩气井生产动态预测方法及设备
CN108131127B (zh) 一种泡沫油型超重油油田的生产气油比的获取方法及装置
Li et al. Coupled CO2 sequestration simulation using Abaqus and Eclipse
CN110083852B (zh) 井底压力的确定方法、装置及存储介质
CN108681793B (zh) 深层油藏采油指数预测方法及装置
CN108319743B (zh) 古老海相碳酸盐岩油气资源丰度预测方法及装置
CN108363114B (zh) 致密油甜点区评价方法及装置
CN108952676A (zh) 一种页岩气藏非均质性评价方法及其装置
CN111323814B (zh) 基于岩石物理模板定量确定砂岩储层含水饱和度的方法
CN112613171A (zh) 一种煤储层渗透率确定方法及装置
CN108629463B (zh) 地应力变化预测方法和装置
CN114580143B (zh) 矿井涌水量动态预测方法及相关装置
CN113392372B (zh) 气井全生命周期无阻流量确定方法及系统
CN112081565B (zh) 水驱油藏体积波及系数的确定方法及装置
CN108868724A (zh) 一种确定凝析气井气举增产油气量的方法及其装置
CN113445989A (zh) 致密油藏压裂水平井的产能预测方法、介质、终端和装置
CN108108512B (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