CN112682033B - 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 - Google Patents

孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 Download PDF

Info

Publication number
CN112682033B
CN112682033B CN202011376593.3A CN202011376593A CN112682033B CN 112682033 B CN112682033 B CN 112682033B CN 202011376593 A CN202011376593 A CN 202011376593A CN 112682033 B CN112682033 B CN 112682033B
Authority
CN
China
Prior art keywords
dimensionless
yield
hole
parameters
normalized
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.)
Active
Application number
CN202011376593.3A
Other languages
English (en)
Other versions
CN112682033A (zh
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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202011376593.3A priority Critical patent/CN112682033B/zh
Publication of CN112682033A publication Critical patent/CN112682033A/zh
Application granted granted Critical
Publication of CN112682033B publication Critical patent/CN112682033B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种孔洞裂缝型碳酸盐岩储集体参数确定方法及装置，该方法包括：筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；将储集体中流体的物性参数值，输入至双孔单渗数学模型中，解析获得无因次产量拉氏解；进行反演获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；确定边界参数；通过Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制Blasingame图版；确定孔洞裂缝型碳酸盐岩储集体的物性参数。本发明可以准确有效地确定物性参数和边界参数。

Description

孔洞裂缝型碳酸盐岩储集体参数确定方法及装置
技术领域
本发明涉及碳酸盐油藏开发技术领域，尤其涉及一种孔洞裂缝型碳酸盐岩储集体参数确定方法及装置。
背景技术
近年来，随着我国油气勘探开发程度越来越高，开发难度也逐渐增大。总体上，勘探开发由陆相向海相转移，由碎屑岩储集体向碳酸盐岩储集体转移，由东部地区向西部偏远地区转移，由国内区块向海外区块转移。塔里木盆地是国内主要的缝洞型碳酸盐岩油藏，其地质储量十分丰富，已达到16.7亿吨，且随着勘探开发技术的不断进步，原油探明储量仍在增长，具有广阔的开发前景。
孔洞-裂缝型碳酸盐岩储集体由于成因独特，所以储集体性质以及生产动态特征与砂岩油藏的特征差异较大。而且，孔洞-裂缝型碳酸盐岩储集体非均质性强，储集体连通性差，砂岩中常用的物质平衡方法、压力不稳定试井方法等适用性较差。因此在开发过程中缺少一种准确的有效进行孔洞裂缝型碳酸盐岩储集体参数的确定方法。
发明内容
本发明实施例提出一种孔洞裂缝型碳酸盐岩储集体参数确定方法，用以准确有效地确定孔洞裂缝型碳酸盐岩储集体的物性参数和边界参数，该方法包括：
从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；
将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；
对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；
基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；
对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；
基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。
本发明实施例提出一种孔洞裂缝型碳酸盐岩储集体参数确定装置，用以准确有效地确定孔洞裂缝型碳酸盐岩储集体的物性参数和边界参数，该装置包括：
筛选模块，用于从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；
解析模块，用于将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；
第一图版绘制模块，用于对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；
边界参数确定模块，用于基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；
第二图版绘制模块，用于对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；
物性参数确定模块，用于基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。
本发明实施例还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述孔洞裂缝型碳酸盐岩储集体参数确定方法。
本发明实施例还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述孔洞裂缝型碳酸盐岩储集体参数确定方法的计算机程序。
在本发明实施例中，从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。在上述过程中，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制的双孔单渗油井产量不稳定分析曲线图版的准确度高，从而可准确地确定孔洞裂缝型碳酸盐岩储集体的边界参数；之后，根据无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制的孔洞裂缝型碳酸盐岩储集体的Blasingame图版的准确度高，从而使得最后确定的孔洞裂缝型碳酸盐岩储集体的物性参数的准确度高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中孔洞裂缝型碳酸盐岩储集体参数确定方法的流程图；
图2为本发明实施例中双孔单渗数学模型对应的渗流示意图；
图3为本发明实施例中双孔单渗油井产量不稳定分析曲线图版的示意图；
图4为本发明实施例中无因次泄油半径对产量不稳定分析曲线图版的示意图；
图5为本发明实施例中窜流系数对产量不稳定分析曲线图版的示意图；
图6为本发明实施例中弹性储容比对产量不稳定分析曲线图版的示意图；
图7为本发明实施例中无因次泄油半径对Blasingame产量不稳定分析曲线图版的示意图；
图8为本发明实施例中窜流系数对Blasingame产量不稳定分析曲线图版的示意图；
图9为本发明实施例中弹性储容比对Blasingame产量不稳定分析曲线图版的示意图；
图10为本发明实施例中孔洞裂缝型碳酸盐岩储集体参数确定装置的示意图；
图11为本发明实施例中孔洞裂缝型碳酸盐岩储集体参数确定装置的另一示意图；
图12为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
图1为本发明实施例中孔洞裂缝型碳酸盐岩储集体参数确定方法的流程图，如图1所示，该方法包括：
步骤101，从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；
步骤102，将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；
步骤103，对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；
步骤104，基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；
步骤105，对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；
步骤106，基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。
在本发明实施例中，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制的双孔单渗油井产量不稳定分析曲线图版的准确度高，从而可准确地确定孔洞裂缝型碳酸盐岩储集体的边界参数；之后，根据无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制的孔洞裂缝型碳酸盐岩储集体的Blasingame图版的准确度高，从而使得最后确定的孔洞裂缝型碳酸盐岩储集体的物性参数的准确度高。
具体实施时，由于本发明实施例针对的是孔洞裂缝型碳酸盐岩储集体，在步骤101中，从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体。上述单井控制的孔洞裂缝型碳酸盐岩储集体为中溶蚀孔洞发育，且溶蚀孔洞直径和其分布密度均较小，而天然裂缝发育且间距大的储集体。在这种情况下可认为溶蚀孔洞的存在改变了基质的储渗能力，即将溶蚀孔洞的孔隙叠加到基质上，形成基质孔隙，裂缝直接向井筒供液，而基质孔隙向裂缝发生窜流，因此，其对应的数学模型为双孔单渗数学模型。或当碳酸盐岩储集体溶蚀孔洞不发育，即对于裂缝型储集体，只发育均匀小尺度分布的裂缝，若考虑基质会向裂缝发生窜流，此时也可用双孔单渗数学模型对储集体进行数学描述。
图2为本发明实施例中双孔单渗数学模型对应的渗流示意图，含有溶蚀孔洞的改进的基质系统(简称基质孔隙)是主要的储集空间，裂缝是主要的渗流通道。其基本假设条件如下：(1)单井以定产量生产，储集体初始压力为pi；(2)储集体岩石和流体微可压缩，且其压缩系数为常数；(3)流体的流动满足达西定律；(4)忽略重力和毛管力的影响。
在一实施例中，流体的物性参数包括单井的流压；
所述方法还包括：将单井的油压输入至统一多相流计算模型中，计算所述单井的流压。
在上述实施例中，相当于对生产井(单井)的油压进行折算，使得最后得到的流压非常准确。另外，如果有产水的话建议可以使用多相流动的统一流动模型进行计算。
在一实施例中，统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管横截面积；dL为管柱长度。
在一实施例中，孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型采用如下公式表示：
初始条件为：pfD(rD,0)＝pmD(rD,0)＝0内边界条件为：
三种外边界条件为：
无限大：
圆形定压：
圆形封闭：
无量纲量定义如下：
其中：kf为裂缝的渗透率，μm2；km为基质孔隙的渗透率，μm2；φf为裂缝孔隙度，无因次；φm为基质孔隙度，无因次；Ctf是裂缝总压缩系数，MPa-1；Ctm是基质孔隙的总压缩系数，MPa-1；pf是单井的流压，MPa；pm是基质孔隙系统在某一时刻某一点的地层压力，MPa；ω是弹性储容比，无因次；λ是窜流系数，无因次；α是形状因子，无因次；pmD是无量纲基质压力；pfD是无量纲裂缝压力。
在一实施例中，所述双孔单渗油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述孔洞裂缝型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
在步骤102中，将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解。在边界类型为无限大时，无因次井底压力拉氏解如下：
在边界类型为圆形定压时，无因次井底压力拉氏解如下：
在边界类型为圆形封闭时，无因次井底压力拉氏解如下：
式中，拟稳态窜流时：
非稳态窜流时：
利用Duhamel原理，根据上述无因次井底压力拉氏解，可得无因次产量拉氏解为：
在步骤103中，对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间tD、无因次产量qD、无因次产量积分qDi和无因次产量积分导数qDid。
其中，NpD为油井地面产量，m3/d。
之后，绘制双孔单渗油井产量不稳定分析曲线图版，图3为本发明实施例中双孔单渗油井产量不稳定分析曲线图版的示意图，上述曲线可大体分为以下五个阶段：
阶段Ⅰ)：裂缝早期线性流阶段；
阶段Ⅱ)：裂缝径向流阶段，在压力波未波及基质孔隙之前，流体沿裂缝发生径向流动，无因次产量随无因次物质平衡拟时间变化的曲线(以下简称无因次产量曲线)下降变缓慢，无因次产量积分随无因次物质平衡拟时间变化的曲线(以下简称无因次产量积分曲线)与无因次产量积分导数随无因次物质平衡拟时间变化的曲线(以下简称因次产量积分导数曲线)无近似平行；
阶段Ⅲ)：基质孔隙向裂缝的窜流阶段，无因次产量积分导数曲线下凹，主要受基质孔隙的窜流系数和弹性储容比的影响，同时，对于拟稳态窜流，凹子深窄，对于非稳态窜流，凹子浅长；
阶段Ⅳ)：整个储集体(裂缝和基质孔隙)的径向流阶段，压力波未到达外边界，无因次产量积分曲线与无因次产量积分导数曲线近似平行；
阶段Ⅴ)：边界控制流阶段，其中对于圆形封闭这一边界类型，无因次产量积分曲线及其产量积分导数曲线重合为斜率为“-1”的直线。
在一实施例中，所述边界参数包括边界距离、是否有无水体能量补充；
所述物性参数包括油井动态储量、水体大小、水侵速度、水侵量、窜流系数和弹性储容比中的至少一种。
基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数的具体过程如下：
(1)确定敏感性参数，例如图3中，确定敏感性参数为无因次泄油半径reD、窜流系数λ、弹性储容比ω；
(2)从双孔单渗油井产量不稳定分析曲线图版中分离出敏感性参数对产量不稳定分析曲线图版，例如，图4为本发明实施例中无因次泄油半径对产量不稳定分析曲线图版的示意图，图5为本发明实施例中窜流系数对产量不稳定分析曲线图版的示意图，图6为本发明实施例中弹性储容比对产量不稳定分析曲线图版的示意图。由图4可以看出，reD值越大，全系统流体径向流持续时间就越长，边界控制流阶段出现的时间越晚。由图5可以看出，λ值主要影响产量积分导数曲线“凹子”出现的早晚，即影响基质孔隙向裂缝中窜流发生的时间，λ值越小，凹子越靠后，则基质孔隙发生窜流的时间越晚。由图6可以看出，ω值主要影响产量积分导数曲线“凹子”的深浅和持续时间，ω值越大，基质孔隙储容比越小，凹子越浅越窄，相应裂缝储容能力越大，流体在流动过程中压力损耗越小，则裂缝的径向流持续时间越长，发生窜流的时间相应越晚，发生窜流的持续时间则越短。
(3)利用敏感性参数对产量不稳定分析曲线图版确定边界距离和是否有无水体能量补充。例如，在图4-图6中，边界距离可以确定为400m，且可以确定有水体能量补充。
在一实施例中，采用如下的Blasingame变换公式，对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；为无量纲井底流压
在上述实施例中，无因次物质平衡拟时间tD、无因次产量qD、无因次产量积分qDi和无因次产量积分导数qDid分别变换为了无因次规整化物质平衡拟时间tDd、无因次规整化产量qDd、无因次规整化产量积分qDdi和无因次规整化产量积分导数qDdid。通过上述处理，实现了边界控制流阶段的归一化处理，这样减少了参数的不确定性。
针对步骤106，根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版。
然后，可以从孔洞裂缝型碳酸盐岩储集体的Blasingame图版中分离出敏感性参数对对Blasingame产量不稳定分析曲线图版，例如，图7为本发明实施例中无因次泄油半径对Blasingame产量不稳定分析曲线图版的示意图，图8为本发明实施例中窜流系数对Blasingame产量不稳定分析曲线图版的示意图，图9为本发明实施例中弹性储容比对Blasingame产量不稳定分析曲线图版的示意图。
由图7可以看出，当其他参数一定时，改变无因次泄油半径reD的大小，边界控制阶段的斜率为“-1”的归一化曲线不会发生左右移动，无因次泄油半径reD越大，只会使曲线的不稳定流阶段曲线越往下，基质孔隙向裂缝开始发生窜流的物质平衡时间越靠前，即凹子越靠前，但凹子形状不会发生改变，全区径向流阶段持续时间越长。
由图8可以看出，窜流系数λ值越小，无因次规整化产量积分导数曲线“凹子”越靠右，即基质孔隙向裂缝中的窜流阶段发生的越晚，λ对边界控制流阶段无影响，与图5类似。
由图9可以看出，弹性储容比ω的大小主要影响基质孔隙向裂缝窜流阶段的深度和宽度，对窜流发生的开始时间无影响。弹性储容比ω值越大，基质储容比越小，凹子越浅越窄，则相应的全系统径向流开始的时间越靠前，边界控制流阶段出现的越早。
可见，当弹性储容比ω值一定时，边界控制阶段规整化曲线一定，则改变reD仅影响不稳定流阶段的曲线的位置，而λ值的大小也仅影响不稳定流阶段基质向裂缝中窜流阶段发生的早晚。
最后，基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，通过典型曲线拟合确定孔洞裂缝型碳酸盐岩储集体的动态控制储量、窜流系数、弹性储容比，其中水侵量计算公式为：
其中，Wei＝Nw×pi×Cw，
其中：Wei为水体潜在的最大水侵量，m3；We为累计水侵量，m3；pi为原始油藏压力，MPa；paq为水体压力，MPa；p为当前油藏的平均压力，MPa；J为水侵指数，m3/d/MPa；t为生产时间，d；Nw为水体大小，m3；Cw为水体的压缩系数，MPa-1；f为水体侵入系数，取值为0～1；θ为水体的覆盖范围，取值为0～2π；μw为水体粘度，mPa·s。
综上所述，在本发明实施例提出的方法中，从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。在上述过程中，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制的双孔单渗油井产量不稳定分析曲线图版的准确度高，从而可准确地确定孔洞裂缝型碳酸盐岩储集体的边界参数；之后，根据无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制的孔洞裂缝型碳酸盐岩储集体的Blasingame图版的准确度高，从而使得最后确定的孔洞裂缝型碳酸盐岩储集体的物性参数的准确度高。
本发明实施例还提出一种孔洞裂缝型碳酸盐岩储集体参数确定装置，其原理与孔洞裂缝型碳酸盐岩储集体参数确定方法类似，这里不再赘述。
图10为本发明实施例中孔洞裂缝型碳酸盐岩储集体参数确定装置的示意图，该装置包括：
筛选模块1001，用于从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；
解析模块1002，用于将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；
第一图版绘制模块1003，用于对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；
边界参数确定模块1004，用于基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；
第二图版绘制模块1005，用于对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；
物性参数确定模块1006，用于基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。
在一实施例中，流体的物性参数包括单井的流压；
所述装置还包括计算模块1007，如图11所示为本发明实施例中孔洞裂缝型碳酸盐岩储集体参数确定装置的另一示意图，计算模块1007用于：将单井的油压输入至统一多相流计算模型中，计算所述单井的流压。
在一实施例中，统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管横截面积；dL为管柱长度。
在一实施例中，孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型采用如下公式表示：
初始条件为：pfD(rD,0)＝pmD(rD,0)＝0内边界条件为：
三种外边界条件为：
无限大：
圆形定压：
圆形封闭：
无量纲量定义如下：
其中：kf为裂缝的渗透率，μm2；km为基质孔隙的渗透率，μm2；φf为裂缝孔隙度，无因次；φm为基质孔隙度，无因次；Ctf是裂缝总压缩系数，MPa-1；Ctm是基质孔隙的总压缩系数，MPa-1；pf是单井的流压，MPa；pm是基质孔隙系统在某一时刻某一点的地层压力，MPa；ω是弹性储容比，无因次；λ是窜流系数，无因次；α是形状因子，无因次；pmD是无量纲基质压力；pfD是无量纲裂缝压力。
在一实施例中，所述双孔单渗油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述孔洞裂缝型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
在一实施例中，第二图版绘制模块具体用于：
采用如下的Blasingame变换公式，对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；为无量纲井底流压。
在一实施例中，所述边界参数包括边界距离、是否有无水体能量补充；
所述物性参数包括油井动态储量、水体大小、水侵速度、水侵量、窜流系数和弹性储容比中的至少一种。
综上所述，在本发明实施例提出的装置中，从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数。在上述过程中，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制的双孔单渗油井产量不稳定分析曲线图版的准确度高，从而可准确地确定孔洞裂缝型碳酸盐岩储集体的边界参数；之后，根据无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制的孔洞裂缝型碳酸盐岩储集体的Blasingame图版的准确度高，从而使得最后确定的孔洞裂缝型碳酸盐岩储集体的物性参数的准确度高。
本申请的实施例还提供一种计算机设备，图12为本发明实施例中计算机设备的示意图，该计算机设备能够实现上述实施例中的孔洞裂缝型碳酸盐岩储集体参数确定方法中全部步骤，所述电子设备具体包括如下内容：
处理器(processor)1201、存储器(memory)1202、通信接口(CommunicationsInterface)1203和总线1204；
其中，所述处理器1201、存储器1202、通信接口1203通过所述总线1204完成相互间的通信；所述通信接口1203用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输；
所述处理器1201用于调用所述存储器1202中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的孔洞裂缝型碳酸盐岩储集体参数确定方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质，能够实现上述实施例中的孔洞裂缝型碳酸盐岩储集体参数确定方法中全部步骤，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的孔洞裂缝型碳酸盐岩储集体参数确定方法的全部步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种孔洞裂缝型碳酸盐岩储集体参数确定方法，其特征在于，包括：
从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；
将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；
对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；
基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；
对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得双孔单渗模型的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；
基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数；
流体的物性参数包括单井的流压；
所述方法还包括：将单井的油压输入至统一多相流计算模型中，计算所述单井的流压；
统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管横截面积；dL为管柱长度。
2.如权利要求1所述的孔洞裂缝型碳酸盐岩储集体参数确定方法，其特征在于，所述双孔单渗油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述孔洞裂缝型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
3.如权利要求1所述的孔洞裂缝型碳酸盐岩储集体参数确定方法，其特征在于，采用如下的Blasingame变换公式，对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；为无量纲井底流压。
4.如权利要求1所述的孔洞裂缝型碳酸盐岩储集体参数确定方法，其特征在于，所述边界参数包括边界距离、是否有无水体能量补充；
所述物性参数包括油井动态储量、水体大小、水侵速度、水侵量、窜流系数和弹性储容比中的至少一种。
5.一种孔洞裂缝型碳酸盐岩储集体参数确定装置，其特征在于，包括：
筛选模块，用于从单井连通的碳酸盐岩储集体中，筛选出单井控制的孔洞裂缝型碳酸盐岩储集体；
解析模块，用于将孔洞裂缝型碳酸盐岩储集体中流体的物性参数值，输入至孔洞裂缝型碳酸盐岩储集体的双孔单渗数学模型中，解析获得无因次产量拉氏解；
第一图版绘制模块，用于对不同边界参数下的无因次产量拉氏解进行反演，获得无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，并绘制双孔单渗油井产量不稳定分析曲线图版；
边界参数确定模块，用于基于双孔单渗油井产量不稳定分析曲线图版，确定孔洞裂缝型碳酸盐岩储集体的边界参数；
第二图版绘制模块，用于对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，并绘制孔洞裂缝型碳酸盐岩储集体的Blasingame图版；
物性参数确定模块，用于基于孔洞裂缝型碳酸盐岩储集体的Blasingame图版，确定孔洞裂缝型碳酸盐岩储集体的物性参数；
流体的物性参数包括单井的流压；
所述装置还包括计算模块，用于：将单井的油压输入至统一多相流计算模型中，计算所述单井的流压；
统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管横截面积；dL为管柱长度。
6.如权利要求5所述的孔洞裂缝型碳酸盐岩储集体参数确定装置，其特征在于，所述双孔单渗油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述孔洞裂缝型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
7.如权利要求5所述的孔洞裂缝型碳酸盐岩储集体参数确定装置，其特征在于，第二图版绘制模块具体用于：
采用如下的Blasingame变换公式，对无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；为无量纲井底流压。
8.如权利要求5所述的孔洞裂缝型碳酸盐岩储集体参数确定装置，其特征在于，所述边界参数包括边界距离、是否有无水体能量补充；
所述物性参数包括油井动态储量、水体大小、水侵速度、水侵量、窜流系数和弹性储容比中的至少一种。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一项所述方法的计算机程序。
CN202011376593.3A 2020-11-30 2020-11-30 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 Active CN112682033B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011376593.3A CN112682033B (zh) 2020-11-30 2020-11-30 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011376593.3A CN112682033B (zh) 2020-11-30 2020-11-30 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置

Publications (2)

Publication Number Publication Date
CN112682033A CN112682033A (zh) 2021-04-20
CN112682033B true CN112682033B (zh) 2023-08-22

Family

ID=75446995

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011376593.3A Active CN112682033B (zh) 2020-11-30 2020-11-30 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置

Country Status (1)

Country Link
CN (1) CN112682033B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016192077A1 (zh) * 2015-06-04 2016-12-08 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 一种致密气压裂水平井数值试井模型建立求解方法
CN110984973A (zh) * 2019-11-29 2020-04-10 中国石油集团川庆钻探工程有限公司 缝洞型碳酸盐岩气藏单井控制储量的确定方法
AU2020100337A4 (en) * 2020-01-09 2020-04-16 China University Of Geosciences, Beijing Method and system for generating unsteady-state well test chart of highly-deviated well in gas reservoir
CN111581786A (zh) * 2020-04-19 2020-08-25 东北石油大学 用于分析缝洞串联模式双孔复合储层的试井解释模型
CN111794740A (zh) * 2020-08-07 2020-10-20 西南石油大学 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170191348A1 (en) * 2016-01-04 2017-07-06 Saudi Arabian Oil Company Modeling to characterize fractures network in homogeneous petroleum reservoirs
CN110334868B (zh) * 2019-07-08 2020-12-08 西南石油大学 一种耦合流体流动与地质应力预测最优焖井时间的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016192077A1 (zh) * 2015-06-04 2016-12-08 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 一种致密气压裂水平井数值试井模型建立求解方法
CN110984973A (zh) * 2019-11-29 2020-04-10 中国石油集团川庆钻探工程有限公司 缝洞型碳酸盐岩气藏单井控制储量的确定方法
AU2020100337A4 (en) * 2020-01-09 2020-04-16 China University Of Geosciences, Beijing Method and system for generating unsteady-state well test chart of highly-deviated well in gas reservoir
CN111581786A (zh) * 2020-04-19 2020-08-25 东北石油大学 用于分析缝洞串联模式双孔复合储层的试井解释模型
CN111794740A (zh) * 2020-08-07 2020-10-20 西南石油大学 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法

Also Published As

Publication number Publication date
CN112682033A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN105840187B (zh) 致密性油藏水平井分段压裂产能计算方法
CN104834807B (zh) 一种基于分形理论的应力敏感储层相对渗透率计算方法
CN110096718B (zh) 一种获取碳酸盐岩油藏中溶洞的体积的方法
CN110107277B (zh) 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法
CN116306385B (zh) 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质
CN112377182B (zh) 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置
CN109388866A (zh) 一种沥青质沉淀油藏流固耦合等效数值模拟方法
Cui et al. Comprehensive proppant settling model in hydraulic fractures of unconventional gas reservoir considering multifactorial influence
Sun et al. A new numerical well test method of multi-scale discrete fractured tight sandstone gas reservoirs and its application in the Kelasu Gas Field of the Tarim Basin
Suzuki et al. Characterization of 3D printed fracture networks
CN112682033B (zh) 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置
Meng et al. Study of water Huff-n-Puff in low-permeability oil reservoirs with horizontal fractures: a case study of Chang 6 reservoir in Yanchang, China
CN112377184B (zh) 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置
Zhu et al. Multicluster fractures propagation during temporary plugging fracturing in naturally fractured reservoirs integrated with dynamic perforation erosion
Zhang et al. Analytical method for performance evaluation of fractured horizontal wells in tight reservoirs
CN112377181B (zh) 定容型碳酸盐岩储集体的参数确定方法及装置
Xulin et al. Numerical simulation for coupling flow and geomechanics in embedded discrete fracture model based on XFEM-MBEM
CN113530534A (zh) 储层改造方法及装置
CN112377185B (zh) 单洞与多重介质复合型储层参数敏感性分析方法及装置
Mei et al. Analytical Time-Dependent Shape Factor for Counter-Current Imbibition in Fractal Fractured Reservoirs
CN114086933B (zh) 一种试井解释方法、装置、存储介质以及计算机设备
CN113586023B (zh) 页岩油藏压后闷井时间的确定方法及设备
CN110671102A (zh) 一种气井临界出砂压差的确定方法及系统
CN114086934B (zh) 一种试井解释方法、装置、存储介质以及计算机设备
CN112668136B (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