CN112377181A - 定容型碳酸盐岩储集体的参数确定方法及装置 - Google Patents
定容型碳酸盐岩储集体的参数确定方法及装置 Download PDFInfo
- Publication number
- CN112377181A CN112377181A CN202011376514.9A CN202011376514A CN112377181A CN 112377181 A CN112377181 A CN 112377181A CN 202011376514 A CN202011376514 A CN 202011376514A CN 112377181 A CN112377181 A CN 112377181A
- Authority
- CN
- China
- Prior art keywords
- dimensionless
- yield
- normalized
- well
- integral
- 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
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 137
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011435 rock Substances 0.000 title claims abstract description 32
- 239000000463 material Substances 0.000 claims abstract description 91
- 238000004088 simulation Methods 0.000 claims abstract description 70
- 239000003129 oil well Substances 0.000 claims abstract description 38
- 238000004458 analytical method Methods 0.000 claims abstract description 37
- 239000012530 fluid Substances 0.000 claims abstract description 24
- 230000000704 physical effect Effects 0.000 claims abstract description 22
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000012216 screening Methods 0.000 claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims description 41
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 36
- 238000004364 calculation method Methods 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 16
- 239000000126 substance Substances 0.000 claims description 13
- 230000009545 invasion Effects 0.000 claims description 12
- 239000007788 liquid Substances 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000009469 supplementation Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 19
- 239000003921 oil Substances 0.000 description 10
- 238000011161 development Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 239000010779 crude oil Substances 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000013589 supplement Substances 0.000 description 4
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
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
- 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
-
- 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
-
- 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)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (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为本发明实施例中无因次泄油半径对Blasingame产量不稳定分析曲线图版；
图6为本发明实施例中定容型碳酸盐岩储集体的参数确定装置的示意图；
图7为本发明实施例中定容型碳酸盐岩储集体的参数确定装置的另一示意图；
图8为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
图1为本发明实施例中定容型碳酸盐岩储集体的参数确定方法的流程图，如图1所示，该方法包括：
步骤101，从单井连通的碳酸盐岩储集体中，筛选出单井控制的定容型碳酸盐岩储集体；
步骤102，根据单井控制的定容型碳酸盐岩储集体中流体的物性参数值，计算单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数；
步骤103，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制单井的单重介质油井产量不稳定分析曲线图版；
步骤104，基于单井的单重介质油井产量不稳定分析曲线图版，确定定容型碳酸盐岩储集体的边界参数；
步骤105，对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数；
步骤106，根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制定容型碳酸盐岩储集体的Blasingame图版；
步骤107，基于定容型碳酸盐岩储集体的Blasingame图版，确定定容型碳酸盐岩储集体的物性参数。
在本发明实施例中，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制的不同边界条件下的单井的单重介质油井产量不稳定分析曲线图版的准确度高，从而可准确地确定定容型碳酸盐岩储集体的边界参数；之后，根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制的定容型碳酸盐岩储集体的Blasingame图版的准确度高，从而使得最后确定的定容型碳酸盐岩储集体的物性参数的准确度高。
具体实施时，由于本发明实施例针对的是定容型碳酸盐岩储集体，在步骤101中，从单井连通的碳酸盐岩储集体中，筛选出单井控制的定容型碳酸盐岩储集体。上述单井控制的定容型碳酸盐岩储集体为仅有小尺度裂缝发育，基质的储渗能力可以忽略的储集体。定容型碳酸盐岩储集体可用单重介质模型来描述。图2为本发明实施例中单重介质模型的渗流示意图，其基本假设条件如下：(1)单井以定产量生产，储集体初始压力为pi；(2)储集体岩石和流体微可压缩，且其压缩系数为常数；(3)流体的流动满足达西定律；(4)忽略重力和毛管力的影响。
在一实施例中，流体的物性参数包括单井的流压；
所述方法还包括：将油井的油压输入至统一多相流计算模型中，计算所述单井的流压。
在上述实施例中，相当于对生产井(单井)的油压进行折算，使得最后得到的流压非常准确。另外，如果有产水的话建议可以使用多相流动的统一流动模型进行计算。
在一实施例中，统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管的横截面积；dL为管柱长度。
在步骤102中，根据单井控制的定容型碳酸盐岩储集体中流体的物性参数值，计算单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，计算公式如下：
其中：qD为无因次产量；qDi为无因次产量积分；qDid为无因次产量积分导数；tD为无因次物质平衡拟时间；r为任一点到井的距离，m；rw为油井半径，m；re为泄油半径，m；k为储层渗透率，μm2；φ为储层孔隙度，无因次；Ct是储层总压缩系数，MPa-1；h为储层厚度，m；μ为原油粘度，mPa·s；B为原油体积系数，无因次；NpD为油井地面产量，m3/d；t为生产时间，h；pi为原始地层压力，MPa；p为单井的流压，也称为储层在某一时刻某一点的地层压力，MPa；reD为无因次泄油半径。
在一实施例中，所述单重介质油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述定容型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
在步骤103中，根据单井的无因次物质平衡拟时间tD、无因次产量qD、无因次产量积分qDi和无因次产量积分导数qDid，绘制上述单井的单重介质油井产量不稳定分析曲线图版，图3为本发明实施例中单重介质油井产量不稳定分析曲线图版的示意图，通过图3可以看出无因次产量积分随无因次物质平衡拟时间变化的曲线(以下简称无因次产量积分曲线)和无因次产量积分导数随无因次物质平衡拟时间变化的曲线(以下简称无因次产量积分导数曲线)具有明显特征来展示原油的流动特征，而无因次产量随无因次物质平衡拟时间变化的曲线(以下简称无因次产量曲线)不具有明显的特征，则确定边界参数时，不需要考虑无因次产量曲线。
在一实施例中，所述边界参数包括边界距离、边界类型和是否有无水体能量补充中的至少一种；
所述物性参数包括油井动态储量、水体大小、水侵速度和水侵量中的至少一种。
图3中，将定容型碳酸盐岩储集体的渗流分为三个阶段：
阶段I)：早期线性流阶段，表现为平面线形流，该阶段一般出现的较早，若考虑井储和表皮的影响，该阶段容易被覆盖，无因次产量积分曲线与无因次产量积分导数曲线斜率为-1/2；
阶段II)：中晚期径向流阶段，为压力波传到外边界之前的不稳定阶段，无因次产量积分曲线与无因次产量积分导数曲线下降速度变缓，且近似平行，以上2个阶段属于不稳定流阶段；
阶段III)：边界控制流阶段，曲线主要受边界类型的影响，若边界类型为圆形封闭，无因次产量积分曲线与无因次产量积分导数曲线重合为斜率为“-1”的直线；若边界类型为圆形定压，无因次产量曲线和无因次产量积分曲线变水平，无因次产量积分导数曲线下降变为0；若边界类型为无限大，无因次产量积分曲线与无因次产量积分导数曲线持续中期径向流。
对于步骤104，基于单井的单重介质油井产量不稳定分析曲线图版，确定定容型碳酸盐岩储集体的边界参数的具体过程如下：
(1)从单重介质油井产量不稳定分析曲线图版中确定边界类型，例如，前述中国，从图3中可以确定3种边界类型，分别为圆形封闭、圆形定压和无限大。
(2)确定敏感性参数，例如图3中，确定敏感性参数为无因次泄油半径reD；
(3)从单重介质油井产量不稳定分析曲线图版中分离出敏感性参数对产量不稳定分析曲线图版，例如，图4为本发明实施例中无因次泄油半径对产量不稳定分析曲线图版的示意图，由图4可以看出，无因次泄油半径reD主要影响中后期径向流阶段的持续长短，无因次泄油半径越大，压力波传递到边界的时间越晚，无因次产量积分曲线与无因次产量积分导数曲线的边界反应阶段越靠右，即边界控制流出现的越晚，利用该部分曲线即可准确的确定边界距离和是否有无水体能量补充。例如，在图4中，边界距离可以确定为400m，且可以确定有水体能量补充。
在一实施例中，采用如下的Blasingame变换公式，对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；
为无量纲井底流压。
在上述实施例中，无因次物质平衡拟时间tD、无因次产量qD、无因次产量积分qDi和无因次产量积分导数qDid分别变换为了无因次规整化物质平衡拟时间tDd、无因次规整化产量qDd、无因次规整化产量积分qDdi和无因次规整化产量积分导数qDdid。通过上述处理，实现了边界控制流阶段的归一化处理，这样减少了参数的不确定性。
针对步骤106，根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制定容型碳酸盐岩储集体的Blasingame图版。
然后，可以从定容型碳酸盐岩储集体的Blasingame图版中分离出敏感性参数对对Blasingame产量不稳定分析曲线图版，例如，图5为本发明实施例中无因次泄油半径对Blasingame产量不稳定分析曲线图版，由图5可以看出，无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线(以下简称为无因次规整化产量曲线)在边界控制流阶段汇聚成一条斜率为-1的直线，其值不会随着无因次泄油半径reD的改变而左右移动，可见，经过Blasingame变换的处理，reD不再影响边界控制流阶段的影响，仅影响不稳定流阶段(阶段I和阶段II)，这对于储集体物性参数的确定减少了较多的不确定性。尤其当含有水体能量补充时，若实际油井的生产数据规整化处理后，最终没有拟合上归一化的边界控制流阶段，则表示生产数据已经受到了水体的影响，水体已开始侵入储集体，此时可通过结合水侵模型处理实际数据来获得定容型碳酸盐岩储集体的Blasingame图版。
针对步骤107，基于定容型碳酸盐岩储集体的Blasingame图版，通过图版拟合，确定定容型碳酸盐岩储集体的物性参数，如reD、渗透率、以及水体大小和水侵速度。
其中，水体大小的公式为
其中，
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图版的准确度高，从而使得最后确定的定容型碳酸盐岩储集体的物性参数的准确度高。
本发明实施例还提出一种定容型碳酸盐岩储集体的参数确定装置，其原理与定容型碳酸盐岩储集体的参数确定方法类似，这里不再赘述。
图6为本发明实施例中定容型碳酸盐岩储集体的参数确定装置的示意图，该装置包括：
筛选模块601，用于从单井连通的碳酸盐岩储集体中，筛选出单井控制的定容型碳酸盐岩储集体；
第一计算模块602，用于根据单井控制的定容型碳酸盐岩储集体中流体的物性参数值，计算单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数；
第一图版绘制模块603，用于根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制单井的单重介质油井产量不稳定分析曲线图版；
边界参数确定模块604，用于基于单井的单重介质油井产量不稳定分析曲线图版，确定定容型碳酸盐岩储集体的边界参数；
第二计算模块605，用于对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数；
第二图版绘制模块606，用于根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制定容型碳酸盐岩储集体的Blasingame图版；
物性参数确定模块607，用于基于定容型碳酸盐岩储集体的Blasingame图版，确定定容型碳酸盐岩储集体的物性参数。
在一实施例中，流体的物性参数包括单井的流压；
所述装置还包括第三计算模块608，如图7所示为本发明实施例中定容型碳酸盐岩储集体的参数确定装置的另一示意图，第三计算模块608用于：将油井的油压输入至统一多相流计算模型中，计算所述单井的流压。
在一实施例中，统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管横截面积；dL为管柱长度。
在一实施例中，所述单重介质油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述定容型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
在一实施例中，第二计算模块具体用于：
采用如下的Blasingame变换公式，对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；
为无量纲井底流压。
在一实施例中，所述边界参数包括边界距离、边界类型和是否有无水体能量补充中的至少一种；
所述物性参数包括油井动态储量、水体大小、水侵速度和水侵量中的至少一种。
综上所述，在本发明实施例提出的装置中，从单井连通的碳酸盐岩储集体中，筛选出单井控制的定容型碳酸盐岩储集体；根据单井控制的定容型碳酸盐岩储集体中流体的物性参数值，计算单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数；根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制单井的单重介质油井产量不稳定分析曲线图版；基于单井的单重介质油井产量不稳定分析曲线图版，确定定容型碳酸盐岩储集体的边界参数；对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数；根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制定容型碳酸盐岩储集体的Blasingame图版；基于定容型碳酸盐岩储集体的Blasingame图版，确定定容型碳酸盐岩储集体的物性参数。在上述过程中，根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制的不同边界条件下的单井的单重介质油井产量不稳定分析曲线图版的准确度高，从而可准确地确定定容型碳酸盐岩储集体的边界参数；之后，根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制的定容型碳酸盐岩储集体的Blasingame图版的准确度高，从而使得最后确定的定容型碳酸盐岩储集体的物性参数的准确度高。
本申请的实施例还提供一种计算机设备，图8为本发明实施例中计算机设备的示意图，该计算机设备能够实现上述实施例中的定容型碳酸盐岩储集体的参数确定方法中全部步骤，所述电子设备具体包括如下内容：
处理器(processor)801、存储器(memory)802、通信接口(CommunicationsInterface)803和总线804；
其中，所述处理器801、存储器802、通信接口803通过所述总线804完成相互间的通信；所述通信接口803用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输；
所述处理器801用于调用所述存储器802中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的定容型碳酸盐岩储集体的参数确定方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质，能够实现上述实施例中的定容型碳酸盐岩储集体的参数确定方法中全部步骤，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的定容型碳酸盐岩储集体的参数确定方法的全部步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (14)
1.一种定容型碳酸盐岩储集体的参数确定方法，其特征在于，包括：
从单井连通的碳酸盐岩储集体中，筛选出单井控制的定容型碳酸盐岩储集体；
根据单井控制的定容型碳酸盐岩储集体中流体的物性参数值，计算单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数；
根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制单井的单重介质油井产量不稳定分析曲线图版；
基于单井的单重介质油井产量不稳定分析曲线图版，确定定容型碳酸盐岩储集体的边界参数；
对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数；
根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制定容型碳酸盐岩储集体的Blasingame图版；
基于定容型碳酸盐岩储集体的Blasingame图版，确定定容型碳酸盐岩储集体的物性参数。
2.如权利要求1所述的定容型碳酸盐岩储集体的参数确定方法，其特征在于，流体的物性参数包括单井的流压；
所述方法还包括：将油井的油压输入至统一多相流计算模型中，计算所述单井的流压。
3.如权利要求2所述的定容型碳酸盐岩储集体的参数确定方法，其特征在于，统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管的横截面积；dL为管柱长度。
4.如权利要求1所述的定容型碳酸盐岩储集体的参数确定方法，其特征在于，所述单重介质油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述定容型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
5.如权利要求1所述的定容型碳酸盐岩储集体的参数确定方法，其特征在于，采用如下的Blasingame变换公式，对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；
为无量纲井底流压。
6.如权利要求1所述的定容型碳酸盐岩储集体的参数确定方法，其特征在于，所述边界参数包括边界距离、边界类型和是否有无水体能量补充中的至少一种；
所述物性参数包括油井动态储量、水体大小、水侵速度和水侵量中的至少一种。
7.一种定容型碳酸盐岩储集体的参数确定装置，其特征在于，包括：
筛选模块，用于从单井连通的碳酸盐岩储集体中，筛选出单井控制的定容型碳酸盐岩储集体；
第一计算模块，用于根据单井控制的定容型碳酸盐岩储集体中流体的物性参数值，计算单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数；
第一图版绘制模块，用于根据单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数，绘制单井的单重介质油井产量不稳定分析曲线图版；
边界参数确定模块，用于基于单井的单重介质油井产量不稳定分析曲线图版，确定定容型碳酸盐岩储集体的边界参数；
第二计算模块，用于对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数；
第二图版绘制模块，用于根据单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化产量积分和无因次规整化产量积分导数，绘制定容型碳酸盐岩储集体的Blasingame图版；
物性参数确定模块，用于基于定容型碳酸盐岩储集体的Blasingame图版，确定定容型碳酸盐岩储集体的物性参数。
8.如权利要求7所述的定容型碳酸盐岩储集体的参数确定装置，其特征在于，流体的物性参数包括单井的流压；
所述装置还包括第三计算模块，用于：将油井的油压输入至统一多相流计算模型中，计算所述单井的流压。
9.如权利要求8所述的定容型碳酸盐岩储集体的参数确定装置，其特征在于，统一多相流计算模型采用如下公式表示：
其中，τS为段塞体剪切应力，τWF为气芯中液膜的剪切应力，τWG为气芯中气泡的剪切应力，LS为段塞体长度，LU为段塞单元长度，LF为段塞单元中液膜的长度，SF为液膜的周长，SG为段塞单元中气芯的周长；dp为流压与油压的差值；ρ为流体密度；U为流体体积；g为重力加速度；θ为管柱倾角；d为管柱横截面积；A为油管的横截面积；dL为管柱长度。
10.如权利要求7所述的定容型碳酸盐岩储集体的参数确定装置，其特征在于，所述单重介质油井产量不稳定分析曲线图版包括无因次产量随无因次物质平衡拟时间变化的曲线、无因次产量积分随无因次物质平衡拟时间变化的曲线和无因次产量积分导数随无因次物质平衡拟时间变化的曲线；
所述定容型碳酸盐岩储集体的Blasingame图版包括无因次规整化产量随无因次规整化物质平衡拟时间变化的曲线、无因次规整化产量积分随无因次规整化物质平衡拟时间变化的曲线和无因次规整化产量积分导数随无因次规整化物质平衡拟时间变化的曲线。
11.如权利要求7所述的定容型碳酸盐岩储集体的参数确定装置，其特征在于，第二计算模块具体用于：
采用如下的Blasingame变换公式，对单井的无因次物质平衡拟时间、无因次产量、无因次产量积分和无因次产量积分导数进行Blasingame变换，获得单井的无因次规整化物质平衡拟时间、无因次规整化产量、无因次规整化累积产量、无因次规整化产量积分和无因次规整化产量积分导数：
其中，tDd为无因次规整化物质平衡拟时间；qDd为无因次规整化产量；NpDd为无因次规整化累积产量；qDdi为无因次规整化产量积分；qDdid为无因次规整化产量积分导数；tD为无因次物质平衡拟时间；reD为无因次泄油半径；L-1表示拉氏变换；
为无量纲井底流压。
12.如权利要求7所述的定容型碳酸盐岩储集体的参数确定装置，其特征在于，所述边界参数包括边界距离、边界类型和是否有无水体能量补充中的至少一种；
所述物性参数包括油井动态储量、水体大小、水侵速度和水侵量中的至少一种。
13.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6任一项所述方法。
14.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至6任一项所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011376514.9A CN112377181B (zh) | 2020-11-30 | 2020-11-30 | 定容型碳酸盐岩储集体的参数确定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011376514.9A CN112377181B (zh) | 2020-11-30 | 2020-11-30 | 定容型碳酸盐岩储集体的参数确定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112377181A true CN112377181A (zh) | 2021-02-19 |
| CN112377181B CN112377181B (zh) | 2023-12-22 |
Family
ID=74589229
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011376514.9A Active CN112377181B (zh) | 2020-11-30 | 2020-11-30 | 定容型碳酸盐岩储集体的参数确定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112377181B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106703796A (zh) * | 2016-12-30 | 2017-05-24 | 中国石油天然气股份有限公司 | 一种获取油藏的动态储量及水体大小的方法及装置 |
| CN111236931A (zh) * | 2020-01-09 | 2020-06-05 | 中国地质大学（北京） | 一种气藏大斜度井非稳态试井图版生成方法和系统 |
| CN111794740A (zh) * | 2020-08-07 | 2020-10-20 | 西南石油大学 | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 |
-
2020
- 2020-11-30 CN CN202011376514.9A patent/CN112377181B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106703796A (zh) * | 2016-12-30 | 2017-05-24 | 中国石油天然气股份有限公司 | 一种获取油藏的动态储量及水体大小的方法及装置 |
| CN111236931A (zh) * | 2020-01-09 | 2020-06-05 | 中国地质大学（北京） | 一种气藏大斜度井非稳态试井图版生成方法和系统 |
| CN111794740A (zh) * | 2020-08-07 | 2020-10-20 | 西南石油大学 | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 于清艳等: "钻遇大型溶洞的油井参数及储层参数动态评价方法", vol. 41, no. 3, pages 647 - 653 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112377181B (zh) | 2023-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106599449B (zh) | 一种溶洞体积计算的试井解释方法 | |
| CN105243182B (zh) | 致密油压裂水平井的动态储量计算方法 | |
| Wang et al. | Streamline approach for history matching production data | |
| CN110096718B (zh) | 一种获取碳酸盐岩油藏中溶洞的体积的方法 | |
| CN104881522B (zh) | 一种基于流体势场的剩余油表征方法及装置 | |
| CN116306385B (zh) | 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质 | |
| CN109884269A (zh) | 岩心自发渗吸采收率预测方法及系统 | |
| CN110245454A (zh) | 一种低渗气藏产能计算方法 | |
| CN112377182B (zh) | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
| CN110096669B (zh) | 一种获取碳酸盐岩油藏中裂缝体积的方法 | |
| Nakashima et al. | Near-well upscaling for three-phase flows | |
| AU2014352632B2 (en) | Slug flow initiation in fluid flow models | |
| AU2013405168B2 (en) | Optimizing flow control device properties on a producer well in coupled injector-producer liquid flooding systems | |
| WO2019015749A1 (en) | FLUID FLOW PRIMING IN FLUID FLOW MODELS | |
| CN112836438B (zh) | 一种非均质边水气藏气井见水时间预测方法及装置 | |
| CN112377181A (zh) | 定容型碳酸盐岩储集体的参数确定方法及装置 | |
| Soto et al. | Pressure jumps during drainage in macroporous soils | |
| CN112682033B (zh) | 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
| Salimi et al. | The influence of heterogeneity, wetting, and viscosity ratio on oil recovery from vertically fractured reservoirs | |
| CN113887110A (zh) | 油井产能预测方法、装置、电子设备及计算机存储介质 | |
| CN112377184A (zh) | 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置 | |
| CN112270430A (zh) | 潜山油藏底水锥进范围预测方法和装置 | |
| CN115217467A (zh) | 多重介质底水气藏水平井底水突破时间的确定方法、装置 | |
| CN106952169B (zh) | 一种缝洞型油藏流固耦合模型的构建方法 | |
| CN113530534A (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 |