CN105569641A - 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 - Google Patents
一种改善近井压降的水平井瞬态压力快速计算模型建立方法 Download PDFInfo
- Publication number
- CN105569641A CN105569641A CN201510956689.XA CN201510956689A CN105569641A CN 105569641 A CN105569641 A CN 105569641A CN 201510956689 A CN201510956689 A CN 201510956689A CN 105569641 A CN105569641 A CN 105569641A
- Authority
- CN
- China
- Prior art keywords
- horizontal well
- well
- pressure
- pressure drop
- alpha
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
Abstract
本发明公开了一种改善近井压降的水平井瞬态压力快速计算模型建立方法,包括以下步骤:S1.计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解S2.将水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解分解为地层的平面流动部分和顶底封闭边界效应的镜像迭加部分S3.计算水平井近井流动压降;S4.用水平井近井流动压降替换水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解中的顶底封闭边界效应的镜像迭加部分得到水平井的井筒压力计算模型本发明通过分解水平井的压降响应组成,采用近井径向流机制合理描述水平井的近井流动压降,真实反映垂向渗透率、井筒位置、正负表皮及井筒半径的影响。
Description
技术领域
本发明涉及水平井瞬态压力技术领域,特别是涉及一种改善近井压降的水平井瞬态压力快速计算模型建立方法。
背景技术
水平井技术已经成为油气藏高效开发的主要手段,瞬态压力分析是评价水平井产能及油藏的基本途径。由于试井解释要求实时计算和快速响应,应用渗流数学模型解析解仍然是试井评价的主要工具。Ozkan&Raghavan(1991)根据叠加原理通过点源响应沿水平井井筒轴线积分,获得了线源水平井的压力响应计算模型,该模型计算简单、应用广泛,已经成为复杂结构水平井压力响应计算的基本方法。
尽管Ozkan点源迭加水平井压降解整体上反映出水平井的压降动态,但是对于近井地带的压降(对应典型曲线早期)描述却表现出失真:模型对井筒半径的变化不敏感,不能描述有效井径扩大的措施效果;模型对垂向渗透率变化不敏感,时常导致试井解释估计的垂向水平渗透率比过低,与实际的岩心资料严重不符。另外,该模型在顶底封边边界效应的镜像迭加部分的计算量大,历史拟合耗时过高。
发明内容
本发明的目的在于克服现有技术的不足,提供一种改善近井压降的水平井瞬态压力快速计算模型建立方法,通过分解水平井的压降响应组成,采用近井径向流机制合理描述水平井的近井流动压降,真实反映垂向渗透率、井筒位置、正负表皮及井筒半径的影响。
本发明的目的是通过以下技术方案来实现的:一种改善近井压降的水平井瞬态压力快速计算模型建立方法,包括以下步骤:
S1.计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解包括以下子步骤:
S11.定义无因次压力PD、无因次时间tD和无因次坐标(xD,yD,zD),其中:
其中,k是平均渗透率,q是油井产量;φ是孔隙度;Ct是流体压缩系数;μ是流体粘度;B是地层体积系数;Pi是油藏初始压力;P(x,y,z,t)是在空间位置(x,y,z)处t时间的油藏压力;h是油层有效厚度;Lh是水平井半长度;L是无因次长度参考量,取L=Lh;kx,ky和kz分别地层的方向渗透率;x,y,z为空间坐标。
S12.计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解
其中,s为拉普拉斯空间变量;u=s×f(s),f(s)是双重孔隙介质窜流函数,均质是双重孔隙介质的特例,均质油藏取f(s)=1;井筒源点位置为(xwD,ywD,zwD),观测点位置在(xD,yD,zD),当观测点与源点位置相同时即为井筒压力;K0为虚宗变量第二类修正贝塞尔函数;α为积分变量;
S2.将水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解分解为地层的平面流动部分和顶底封闭边界效应的镜像迭加部分其中:
S3.计算水平井近井流动压降,包括以下子步骤:
S31.将根据无因次压力PD的定义,储层厚度为h、井筒半径为rwD、全打开的直井在半径rD处的井筒压力为:
S32.将S31中的压力响应转换为视厚度2Lh的水平井在半径rD处的压力响应
S33.水平井近井流动压降等于从径向流边界rh到井筒rw的流动压降
其中,
S4.用水平井近井流动压降替换水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解中的顶底封闭边界效应的镜像迭加部分得到水平井的井筒压力计算模型
定义xD=xwD=0,yD=ywD=0,则水平井的井筒压力计算简化模型为:
当水平井井筒位于储层顶部或底部时,水平井近井带半径向流范围内的压降,等于从径向流边界re到井筒rw的流动压降
其中,
处于偏心位置的水平井的近井流动压降采用中心位置压降与顶底部位置压降的线性插值表示:
其中,
在考虑水平井的井储和表皮效应时,采用Everdingen和Hurst方法,得到水平井的井筒压力
式中,是不考虑井储和表皮效应的拉普拉斯空间无因次压力;Sh是基于产层厚度h定义的水平井拟表皮系数;CD是无因次井储系数。
水平井表皮系数是以水平井长度为基础的井筒机械表皮Sm,它与以储层厚度为基础的等效拟表皮Sh之间存在等效的压降关系:
对正负井筒机械表皮处理分别为:
当Sm≥0时,采用拟表皮系数形式,取
当Sm<0时,采用有效井筒半径形式,取
本发明的有益效果是:
(1)本发明通过分解水平井的压降响应组成,采用近井径向流机制合理描述水平井的近井流动压降,真实反映垂向渗透率、井筒位置、正负表皮及井筒半径的影响;
(2)本发明提出的水平井瞬时压力快速计算模型的压降动态已完整的反映出水平井的多种流动机制,计算简化,计算速度提高约10倍;
(3)本发明综合垂向渗透率、井筒位置、正负表皮及井筒存储系数的影响,突出近井流动压降的参数敏感性,为准确预测水平井的产能及评价地层特性提供有效工具。
附图说明
图1为本发明一种改善近井压降的水平井瞬态压力快速计算模型建立方法的流程图;
图2为水平井模型示意图;
图3为水平井井筒附近径向流区域示意图;
图4为水平井井筒附近径向流简化模型示意图;
图5为水平井位于储层顶部的半径向流简化模型示意图;
图6为水平井井筒附近径向流简化模型示意图;
图7为无限大油藏均匀流率水平井压降典型曲线;
图8为早期近井径向流示意图;
图9为地层线性流示意图;
图10为系统拟径向流示意图;
图11为垂直-水平渗透率比的影响曲线;
图12为负表皮系数的影响曲线;
图13为正表皮系数的影响曲线;
图14为水平井长度的影响曲线;
图15为井筒位置的影响曲线;
图16为井筒存储系数的影响曲线。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,一种改善近井压降的水平井瞬态压力快速计算模型建立方法,包括以下步骤:
S1.考虑水平井实在各向异性油藏中,在在x方向和y方向为无限大,其中kx、ky和kz分别为地层的方向渗透率。假设油藏中为单相流体渗流,流体为微可压缩,其水平井的长度为2Lh,油藏厚度为h,油藏的上下为封闭边界(z=0和z=h),水平井平行于上下边界,忽略重力作用影响,油藏的初始压力均匀,其油藏和水平井的物理模型示意图如图2所示。计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解包括以下子步骤:
S11.S11.定义无因次压力PD、无因次时间tD和无因次坐标(xD,yD,zD),其中:
其中,k是平均渗透率,单位为μm2,q是油井产量,单位为m3/d;φ是孔隙度,小数;Ct是流体压缩系数,单位为MPa-1;μ是流体粘度,单位为mPa.s;B是地层体积系数,无因次;Pi是油藏初始压力,单位为MPa;P(x,y,z,t)是在空间位置(x,y,z)处t时间的油藏压力单位为,MPa;h是油层有效厚度,单位为m;Lh是水平井半长度,单位为m;L是无因次长度参考量,取L=Lh;kx,ky和kz分别地层的方向渗透率,x,y,z为空间坐标,定义如图1所示。
S12.水平井压力响应关系由Ozkan和Raghavan导出,计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解
其中,s为拉普拉斯空间变量;u=s×f(s),f(s)是双重孔隙介质窜流函数,均质是双重孔隙介质的特例,均质油藏取f(s)=1;井筒源点位置为(xwD,ywD,zwD),观测点位置在(xD,yD,zD),当观测点与源点位置相同时即为井筒压力;K0为虚宗变量第二类修正贝塞尔函数;α为积分变量;
垂向渗透率影响转换为均质化的等效厚度h’的无因次水平井储层厚度hD:
S2.将水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解分解为地层的平面流动部分和顶底封闭边界效应的镜像迭加部分地层的平面流动部分类似于均匀流率的压裂井模型响应,该部分积分可转换为解析关系,计算准确、快速;顶底封闭边界效应的镜像迭加部分反映早期的垂向流动和向井筒收缩的径向流;该部分需要迭加的次数过多(>500),是压力响应计算中最耗时的部分,其中:
水平井近井径向流区域示意图见图3,将近井径向流简化为图4所示模型:假设水平井位于油藏中部,流动初期是近井径向流(相当于储层厚度等于水平井长度的直井),当压降传播到上下边界时,近井区的不稳定过程结束,表现为流线收缩的拟表皮压降(井筒到径向半径的压降)。
S3.计算水平井近井流动压降,包括以下子步骤:
S31.根据无因次压力PD的定义,储层厚度为h、井筒半径为rwD、全打开的直井在半径rD处的压力响应为:
S32.将S31中的压力响应转换为视厚度2Lh的水平井在半径rD处的压力响应
S33.水平井近井流动压降等于从径向流边界rh到井筒rw的流动压降
其中,
S4.用水平井近井流动压降替换水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解中的顶底封闭边界效应的镜像迭加部分得到水平井的井筒压力计算模型
定义xD=xwD=0,yD=ywD=0,则水平井的井筒压力计算模型为:
取无因次参考长度量L=Lh,中心位置的水平井的压力响应:
当水平井井筒偏离储层中心位置时,近井带的流线收缩区域扩大,存在比中心位置更高的附加压降。极端情况是井筒位于储层顶部或底部,流动机制如图5所示。产量为q的井在半径向流条件下的压降相当于2倍产量的全径向流井的压降,根据图5的参数定义,当水平井井筒位于储层顶部或底部时,水平井近井带半径向流范围内的压降,等于从径向流边界re到井筒rw的流动压降
其中,
如图6所示,处于偏心位置的水平井的近井流动压降采用中心位置压降与顶底部位置压降的线性插值表示:
其中,
水平井的压力响应可简化为:
在考虑水平井的井储和表皮效应时,采用Everdingen和Hurst方法,得到水平井的井筒压力
式中,是不考虑井储和表皮效应的拉普拉斯空间无因次压力;Sh是基于产层厚度h定义的水平井拟表皮系数;CD是无因次井储系数。
真实的水平井表皮系数是以水平井长度为基础的井筒机械表皮Sm,它与以储层厚度为基础的等效拟表皮Sh之间存在等效的压降关系:
对正负井筒机械表皮处理分别为:
当Sm≥0时,采用拟表皮系数形式,取
当Sm<0时,采用有效井筒半径形式,取
对水平井的井筒压力(14)式进行数值反演获得实空间压力响应,图7给出了无限大均质油藏均匀流率水平井的双对数压力和压力导数曲线,从图中可以清楚的看到水平井有几个流动阶段:早期纯井储控制期(斜率为1)、近井径向流动的早期径向流期(导数水平段)、外围向水平井流动的地层线性流期(导数斜率1/2)、压降漏斗近似圆形的系统拟径向流期(导数水平段)。本瞬态压力模型已完整的反映出水平井如图8、图9和图10所示各流动期的流动机制,新开发的模型没有镜像迭加、计算耗时少,与原始模型(4)式的计算对比发现,计算速度提高约10倍。
垂直水平渗透率比kz/k影响压力波的垂向传播速度、流线收缩的垂向流动压降,图11反映出kz/k由大变小时的压降典型曲线,kz/k下降导致压降增加,早期导数水平段上升表现出近井渗透率下降,同时垂向压力传播速度低,储层等效厚度放大,早期径向流时间延长。
图12和图13反映出表皮系数对压降动态的影响。表皮为负增长时,井筒有效半径扩大、压降下降,同时近井径向流区域变小、早期径向流存在时间变短(如图10所示)。表皮为正增长时,压降增大、压力导数驼峰上升,井储系数或表皮系数较大时可完全掩盖早期径向流特征。
图14反映出水平井长度的影响,厚度-半长比h/Lh的增加对应水平井长度的下降,导致压降增加。图15反映出水平井井筒位置的影响,当井筒偏离中心位置时,压降小幅度增加。图16反映出井筒存储系数影响,无因次井储系数大时出现压力导数驼峰,并可掩盖早期径向流水平段特征。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (6)
1.一种改善近井压降的水平井瞬态压力快速计算模型建立方法,其特征在于:包括以下步骤:
S1.计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解包括以下子步骤:
S11.定义无因次压力PD、无因次时间tD和无因次坐标(xD,yD,zD),其中:
其中,k是平均渗透率,q是油井产量;φ是孔隙度;Ct是流体压缩系数;μ是流体粘度;B是地层体积系数;Pi是油藏初始压力;P(x,y,z,t)是在空间位置(x,y,z)处t时间的油藏压力;h是油层有效厚度;Lh是水平井半长度;L是无因次长度参考量,取L=Lh;kx,ky和kz分别地层的方向渗透率;x,y,z为空间坐标。
S12.计算水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解
其中,s为拉普拉斯空间变量;u=s×f(s),f(s)是双重孔隙介质窜流函数,均质是双重孔隙介质的特例,均质油藏取f(s)=1;井筒源点位置为(xwD,ywD,zwD),观测点位置在(xD,yD,zD),当观测点与源点位置相同时即为井筒压力;K0为虚宗变量第二类修正贝塞尔函数;α为积分变量;
S2.将水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解分解为地层的平面流动部分和顶底封闭边界效应的镜像迭加部分其中:
S3.计算水平井近井流动压降,包括以下子步骤:
S31.将根据无因次压力PD的定义,储层厚度为h、井筒半径为rwD、全打开的直井在半径rD处的井筒压力为:
S32.将S31中的压力响应转换为视厚度2Lh的水平井在半径rD处的压力响应
S33.水平井近井流动压降等于从径向流边界rh到井筒rw的流动压降
其中,rwD=rw/L,rhD=rh/L=hD/2;
S4.用水平井近井流动压降P~rhD替换水平井在无限大油藏和均匀流率分布条件下的拉普拉斯空间解中的顶底封闭边界效应的镜像迭加部分得到水平井的井筒压力计算模型
2.根据权利要求1所述的一种改善近井压降的水平井瞬态压力快速计算模型建立方法,其特征在于:定义xD=xwD=0,yD=ywD=0,则水平井的井筒压力计算简化模型为:
3.根据权利要求2所述的一种改善近井压降的水平井瞬态压力快速计算模型建立方法,其特征在于:当水平井井筒位于储层顶部或底部时,水平井近井带半径向流范围内的压降,等于从径向流边界re到井筒rw的流动压降
其中,rwD=rw/L,reD=re/L=hD。
4.根据权利要求1所述的一种改善近井压降的水平井瞬态压力快速计算模型建立方法,其特征在于:处于偏心位置的水平井的近井流动压降采用中心位置压降与顶底部位置压降的线性插值表示:
其中,
5.根据权利要求1所述的一种改善近井压降的水平井瞬态压力快速计算模型建立方法,其特征在于:在考虑水平井的井储和表皮效应时,采用Everdingen和Hurst方法,得到水平井的井筒压力
式中,是不考虑井储和表皮效应的拉普拉斯空间无因次压力;Sh是基于产层厚度h定义的水平井拟表皮系数;CD是无因次井储系数。
6.根据权利要求5所述的一种改善近井压降的水平井瞬态压力快速计算模型建立方法,其特征在于:水平井表皮系数是以水平井长度为基础的井筒机械表皮Sm,它与以储层厚度为基础的等效拟表皮Sh之间存在等效的压降关系:
对正负井筒机械表皮处理分别为:
当Sm≥0时,采用拟表皮系数形式,取rwD=rw/L;
当Sm<0时,采用有效井筒半径形式,取Sh=0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510956689.XA CN105569641B (zh) | 2015-12-17 | 2015-12-17 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510956689.XA CN105569641B (zh) | 2015-12-17 | 2015-12-17 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105569641A true CN105569641A (zh) | 2016-05-11 |
CN105569641B CN105569641B (zh) | 2021-09-28 |
Family
ID=55880197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510956689.XA Active CN105569641B (zh) | 2015-12-17 | 2015-12-17 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105569641B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106246171A (zh) * | 2016-09-09 | 2016-12-21 | 西南石油大学 | 部分连通断层边界的无限大双重介质油藏数学建模方法 |
CN106354957A (zh) * | 2016-10-31 | 2017-01-25 | 延安大学 | 一种基于贝塞尔函数的分支水平井开采方法 |
CN107563027A (zh) * | 2017-08-21 | 2018-01-09 | 西南石油大学 | 用于分析多缝洞单元串联结构的试井解释模型与方法 |
CN108442923A (zh) * | 2018-03-29 | 2018-08-24 | 西南石油大学 | 致密油气藏水平井椭圆形泄流区域长、短半轴的确定方法 |
CN109611088A (zh) * | 2018-12-14 | 2019-04-12 | 西南石油大学 | 一种诊断边水早期水侵的试井分析方法 |
CN110826250A (zh) * | 2019-11-19 | 2020-02-21 | 成都北方石油勘探开发技术有限公司 | 用于水平井的钻井液伤害计算方法 |
CN111950112A (zh) * | 2019-05-16 | 2020-11-17 | 中国石油化工股份有限公司 | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 |
CN112196510A (zh) * | 2020-10-22 | 2021-01-08 | 中国石油天然气集团有限公司 | 一种分段压裂大斜度气井试井图版生成方法及系统 |
CN113586042A (zh) * | 2021-07-09 | 2021-11-02 | 中国石油天然气股份有限公司 | 孔隙型碳酸盐岩动态渗透率确定方法及装置 |
CN113719271A (zh) * | 2021-11-03 | 2021-11-30 | 中法渤海地质服务有限公司 | 一种试井设计参数矫正方法 |
CN114109379A (zh) * | 2021-11-30 | 2022-03-01 | 中海油田服务股份有限公司 | 稠油油藏地层的复合增油措施的效果预测方法及装置 |
CN116629154A (zh) * | 2023-05-24 | 2023-08-22 | 西南石油大学 | 分形复合气藏压裂井瞬态压力计算方法、系统和设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1263831A1 (ru) * | 1985-03-06 | 1986-10-15 | Всесоюзный научно-исследовательский институт гидрогеологии и инженерной геологии | Способ определени гидравлических характеристик прискважинной зоны |
CN101446196A (zh) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | 三重介质油藏分支水平井的试井分析方法及装置 |
CN102220865A (zh) * | 2011-05-13 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种检测灰岩地层孔隙压力的方法 |
CN103590824A (zh) * | 2013-10-21 | 2014-02-19 | 中国石油天然气股份有限公司 | 经过多段压裂改造后的致密气藏水平井的产能计算方法 |
CN104500016A (zh) * | 2014-11-17 | 2015-04-08 | 成都理工大学 | 利用页岩气藏压裂施工压降段数据分析储层物性的新方法 |
-
2015
- 2015-12-17 CN CN201510956689.XA patent/CN105569641B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1263831A1 (ru) * | 1985-03-06 | 1986-10-15 | Всесоюзный научно-исследовательский институт гидрогеологии и инженерной геологии | Способ определени гидравлических характеристик прискважинной зоны |
CN101446196A (zh) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | 三重介质油藏分支水平井的试井分析方法及装置 |
CN102220865A (zh) * | 2011-05-13 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种检测灰岩地层孔隙压力的方法 |
CN103590824A (zh) * | 2013-10-21 | 2014-02-19 | 中国石油天然气股份有限公司 | 经过多段压裂改造后的致密气藏水平井的产能计算方法 |
CN104500016A (zh) * | 2014-11-17 | 2015-04-08 | 成都理工大学 | 利用页岩气藏压裂施工压降段数据分析储层物性的新方法 |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106246171A (zh) * | 2016-09-09 | 2016-12-21 | 西南石油大学 | 部分连通断层边界的无限大双重介质油藏数学建模方法 |
CN106246171B (zh) * | 2016-09-09 | 2019-09-20 | 西南石油大学 | 部分连通断层边界的无限大双重介质油藏数学建模方法 |
CN106354957A (zh) * | 2016-10-31 | 2017-01-25 | 延安大学 | 一种基于贝塞尔函数的分支水平井开采方法 |
CN106354957B (zh) * | 2016-10-31 | 2019-06-14 | 延安大学 | 一种基于贝塞尔函数的分支水平井开采方法 |
CN107563027A (zh) * | 2017-08-21 | 2018-01-09 | 西南石油大学 | 用于分析多缝洞单元串联结构的试井解释模型与方法 |
CN108442923A (zh) * | 2018-03-29 | 2018-08-24 | 西南石油大学 | 致密油气藏水平井椭圆形泄流区域长、短半轴的确定方法 |
CN109611088A (zh) * | 2018-12-14 | 2019-04-12 | 西南石油大学 | 一种诊断边水早期水侵的试井分析方法 |
CN111950112A (zh) * | 2019-05-16 | 2020-11-17 | 中国石油化工股份有限公司 | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 |
CN110826250A (zh) * | 2019-11-19 | 2020-02-21 | 成都北方石油勘探开发技术有限公司 | 用于水平井的钻井液伤害计算方法 |
CN112196510A (zh) * | 2020-10-22 | 2021-01-08 | 中国石油天然气集团有限公司 | 一种分段压裂大斜度气井试井图版生成方法及系统 |
CN112196510B (zh) * | 2020-10-22 | 2022-11-04 | 中国石油天然气集团有限公司 | 一种分段压裂大斜度气井试井图版生成方法及系统 |
CN113586042A (zh) * | 2021-07-09 | 2021-11-02 | 中国石油天然气股份有限公司 | 孔隙型碳酸盐岩动态渗透率确定方法及装置 |
CN113586042B (zh) * | 2021-07-09 | 2023-09-26 | 中国石油天然气股份有限公司 | 孔隙型碳酸盐岩动态渗透率确定方法及装置 |
CN113719271A (zh) * | 2021-11-03 | 2021-11-30 | 中法渤海地质服务有限公司 | 一种试井设计参数矫正方法 |
CN114109379A (zh) * | 2021-11-30 | 2022-03-01 | 中海油田服务股份有限公司 | 稠油油藏地层的复合增油措施的效果预测方法及装置 |
CN114109379B (zh) * | 2021-11-30 | 2023-08-18 | 中海油田服务股份有限公司 | 稠油油藏地层的复合增油措施的效果预测方法及装置 |
CN116629154A (zh) * | 2023-05-24 | 2023-08-22 | 西南石油大学 | 分形复合气藏压裂井瞬态压力计算方法、系统和设备 |
CN116629154B (zh) * | 2023-05-24 | 2024-01-09 | 西南石油大学 | 分形复合气藏压裂井瞬态压力计算方法、系统和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN105569641B (zh) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105569641A (zh) | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 | |
Lu et al. | Nonlinear consolidation of vertical drains with coupled radial–vertical flow considering well resistance | |
CN105787174A (zh) | 一种基于响应面法的高堆石坝瞬变流变参数反演方法 | |
CN103995947A (zh) | 改进的煤层底板突水脆弱性评价方法 | |
CN111861107A (zh) | 一种基于强度折减法的滑坡动态稳定性评价方法 | |
CN103835708A (zh) | 一种油气田类比开发指标定量预测方法 | |
CN106204304A (zh) | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 | |
CN105654482B (zh) | 基于数字图像技术的土石混合体表征单元体尺度的确定方法 | |
CN109829217A (zh) | 压裂性裂缝油藏产能模拟方法及装置 | |
Jeng et al. | Three-dimensional consolidation of a porous unsaturated seabed under rubble mound breakwater | |
CN113238007A (zh) | 一种粗糙裂缝内暂堵模拟的实验方法及装置 | |
CN108843296B (zh) | 一种基于多因素影响下的单井重复压裂效果预测方法 | |
Li et al. | Shear behaviour of aeolian sand with different density and confining pressure | |
Luo et al. | Hydro-mechanical coupling experiments on suffusion in sandy gravel foundations containing a partially penetrating cut-off wall | |
Dastpak et al. | Soil erosion due to defective pipes: a hidden hazard beneath our feet | |
Vu et al. | Numerical modelling of soft ground improvement by vacuum preloading considering the varying coefficient of permeability | |
Deshpande et al. | Review and assessment of the theories of stable alluvial channel design | |
Duan et al. | Numerical investigation of the wave/current–induced responses of transient soil around a square mono-pile foundation | |
CN106600075A (zh) | 一种考虑地层各向异性的水平井产能预测新方法 | |
Pan et al. | Analytical Modeling with Laboratory Data and Observations of the Mechanisms of Backward Erosion Piping Development | |
Echeverribar et al. | Numerical simulation of 2D real large scale floods on GPU: the Ebro River | |
CN106874645A (zh) | 一种考虑地层各向异性的水平井产能预测新方法 | |
Ky et al. | Mechanical stability of river banks submitted to fluctuations of the water level | |
Wang et al. | A hydraulic analysis of an online pressure compensating emitter using CFD-CSD technology | |
Chen et al. | Solution methods of hydraulic conductivity in low permeable saturated soil during liquid releasing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |