CN112377185B - 单洞与多重介质复合型储层参数敏感性分析方法及装置 - Google Patents
单洞与多重介质复合型储层参数敏感性分析方法及装置 Download PDFInfo
- Publication number
- CN112377185B CN112377185B CN202011380695.2A CN202011380695A CN112377185B CN 112377185 B CN112377185 B CN 112377185B CN 202011380695 A CN202011380695 A CN 202011380695A CN 112377185 B CN112377185 B CN 112377185B
- Authority
- CN
- China
- Prior art keywords
- yield
- dimensionless
- oil well
- physical parameters
- hole
- 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
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 25
- 239000003129 oil well Substances 0.000 claims abstract description 123
- 238000004458 analytical method Methods 0.000 claims abstract description 88
- 238000013178 mathematical model Methods 0.000 claims abstract description 19
- 230000035945 sensitivity Effects 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 17
- 230000000704 physical effect Effects 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 25
- 238000004519 manufacturing process Methods 0.000 description 22
- 239000011159 matrix material Substances 0.000 description 12
- 230000003628 erosive effect Effects 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000005260 corrosion Methods 0.000 description 7
- 230000007797 corrosion Effects 0.000 description 7
- 239000011148 porous material Substances 0.000 description 7
- 239000012530 fluid Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 239000010779 crude oil Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 239000000758 substrate Substances 0.000 description 5
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000003921 oil Substances 0.000 description 4
- 230000002093 peripheral effect Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000005465 channeling Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种单洞与多重介质复合型储层参数敏感性分析方法及装置，该方法包括：确定目标单洞与多重介质复合型储层的初始物性参数；将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版。本发明可以分析单洞与多重介质复合型储层的参数敏感性，准确度高。
Description
技术领域
本发明涉及缝洞型碳酸盐岩油藏开发技术领域，尤其涉及一种单洞与多重介质复合型储层参数敏感性分析方法及装置。
背景技术
我国缝洞型碳酸盐岩油藏资源丰富，具有广阔的勘探开发前景。缝洞型介质通常具有复杂的内部结构，不仅包含基岩、裂缝而且还有不同尺度空间的溶洞，这些溶洞的空间大小从毫米级跨越到密集。同时，大量的测井、岩心和露头资料显示：受后期构造运动影响，裂缝和溶蚀孔洞充填比较严重加剧了储层的非均质性。
缝洞型介质结构的复杂性和多尺度性使得如何认识和正确描述缝洞型介质中的流体流动规律成为了一项具有挑战性的工作，其主要难点在于流体在介质中流动不仅有渗流而且还存在大空间自由流动。国内外学者针对裂缝性油藏提出三重介质模型，考虑了裂缝系统、基岩系统和溶洞系统间的物质交换，但该模型仍是传统连续性介质模型并没有体现出缝洞型介质中多尺度耦合流动的特点。国内学者采用离散缝洞网络模型在离散裂缝模型的基础上增加了溶洞系统使其能适用于缝洞型介质的研究，对于洞穴和多重介质复合型储层的研究相对较少，尤其缺乏一种单洞与多重介质复合型储层敏感性分析的方法。
发明内容
本发明实施例提出一种单洞与多重介质复合型储层参数敏感性分析方法，用以分析单洞与多重介质复合型储层的参数敏感性，准确度高，该方法包括：
确定目标单洞与多重介质复合型储层的初始物性参数；
将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；
对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；
基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；
从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。
本发明实施例提出一种单洞与多重介质复合型储层参数敏感性分析装置，用以分析单洞与多重介质复合型储层的参数敏感性，准确度高，该装置包括：
初始物性参数确定模块，用于确定目标单洞与多重介质复合型储层的初始物性参数；
解析模块，用于将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；
反演模块，用于对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；
图版绘制模块，用于基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；
敏感性分析模块，用于从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。
本发明实施例还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述单洞与多重介质复合型储层参数敏感性分析方法。
本发明实施例还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述单洞与多重介质复合型储层参数敏感性分析方法的计算机程序。
在本发明实施例中，确定目标单洞与多重介质复合型储层的初始物性参数；将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。在上述过程中，本发明实施例采用了哑铃状缝洞型储层的数学模型，该模型准确度高，使得解析获得的油井井底的无因次产量拉氏解的准确度高，通过反演后获得的结果可绘制油井产量不稳定分析曲线图版的准确度高，即最后分解出的不同物性参数对油井产量不稳定影响分析曲线图版的准确度高，从而可以更准确地进行不同物性参数对油井产量的敏感性影响的分析。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中单洞与多重介质复合型储层参数敏感性分析方法的流程图；
图2为本发明实施例中单洞与多重介质复合型储层数学模型对应的物理模型的示意图；
图3为本发明实施例中单洞与多重介质复合型储层数学模型对应的物理模型的渗流示意图；
图4为本发明实施例中单洞与三重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版的示意图；
图5为本发明实施例中无量纲大型溶洞半径对油井产量不稳定影响分析曲线图版的示意图；
图6为本发明实施例中大型溶洞的弹性储容比对油井产量不稳定影响分析曲线图版的示意图；
图7为本发明实施例中无量纲泄油半径对油井产量不稳定影响分析曲线图版的示意图；
图8为本发明实施例中单洞与双重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版的示意图；
图9为本发明实施例中单洞与单重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版的示意图；
图10为本发明实施例中单洞与多重介质复合型储层参数敏感性分析装置的示意图；
图11为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
图1为本发明实施例中单洞与多重介质复合型储层参数敏感性分析方法的流程图，如图1所示，该方法包括：
步骤101，确定目标单洞与多重介质复合型储层的初始物性参数；
步骤102，将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；
步骤103，对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；
步骤104，基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；
步骤105，从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。
在本发明实施例提出的方法中，本发明实施例采用了哑铃状缝洞型储层的数学模型，该模型准确度高，使得解析获得的油井井底的无因次产量拉氏解的准确度高，通过反演后获得的结果可绘制油井产量不稳定分析曲线图版的准确度高，即最后分解出的不同物性参数对油井产量不稳定影响分析曲线图版的准确度高，从而可以更准确地进行不同物性参数对油井产量的敏感性影响的分析。
具体实施时，多重介质是指裂缝、基质和溶蚀孔洞，图2为本发明实施例中单洞与多重介质复合型储层数学模型对应的物理模型的示意图，基于地震、测井、生产动态等资料，可以获得图2对应的物理模型，其中储层考虑为三维圆柱复合模型，大型溶洞为内部圆柱体，圆柱高度为h，半径为R，外部为多重介质，与圆柱体同心；大型溶洞考虑为自由流，外围多重介质考虑为渗流。图3为本发明实施例中单洞与多重介质复合型储层数学模型对应的物理模型的渗流示意图，在渗流时，溶蚀孔洞和基质通过裂缝向大型溶洞供液，油井从大型溶洞中采出原油。
在一实施例中，单洞与多重介质复合型储层数学模型采用如下公式表示：
其中，kf、km、kv分别为裂缝、基质和溶蚀孔洞的渗透率，μm2；φf、φmφv分别为裂缝、基质和溶蚀孔洞的孔隙度，无因次；φV为大型溶洞的孔隙度，无因次；Ctf、Ctm、Ctv分别为裂缝、基质和溶蚀孔洞的总压缩系数，MPa-1；CtV为大型溶洞的总压缩系数，MPa-1；pf、pm、pv分别为裂缝、基质和溶蚀孔洞在某一时刻某一点的地层压力，MPa；pV为大型溶洞的地层压力，MPa；pi为原始地层压力，MPa；ωf、ωm、ωv分别为裂缝、基质和溶蚀孔洞的弹性储容比，无因次；ωV为大型溶洞的弹性储容比，无因次；λm、λv分别为基质和溶蚀孔洞的窜流系数，无因次；αm、αv分别为基质、溶蚀孔洞的形状因子，无因次；R为大型溶洞的半径，m；r为任一点到井的距离，m；rw为井半径，m；re为泄油半径，m；h为储层厚度，m；μ为原油粘度，mPa·s；B为原油体积系数，无因次；q为油井产量，m3/d；t为生产时间，d；pfD,mD,vD,VD分别为无量纲基质压力、无量纲裂缝压力、无量纲溶蚀孔洞压力和无量纲大型溶洞压力；tD为无量纲生成时间；rD为任一点到井的无量纲距离；RD为无量纲大型溶洞半径；reD为无量纲泄油半径；pwD为无量纲井底流压。
在步骤102中，将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解，具体过程如下：
联立单洞与多重介质复合型储层数学模型中的公式，且通过拉普拉斯变换，可得无因次Laplace空间解表达式：
其中：
利用Duhamel原理，井底无因次产量拉氏解为：
在步骤103中，对所述无因次产量拉氏解进行数值反演，获得无因次产量qD、无因次产量积分qDi、无因次产量积分导数qDid、无因次时间tD。
之后，即可基于无因次产量qD、无因次产量积分qDi、无因次产量积分导数qDid、无因次时间tD，绘制单洞与多重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版。
由于在本发明实施例中，多重介质是指裂缝、基质和溶蚀孔洞，可分别考虑有三重介质、双重介质和一重介质的情况。
在一实施例中，所述油井产量不稳定分析曲线图版包括无因次产量随无因次时间变化曲线、无因次产量积分随无因次时间变化曲线、无因次产量积分导数随无因次时间变化曲线。
图4为本发明实施例中单洞与三重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版的示意图，该图版以无因次产量qD的对数、无因次产量积分qDi的对数和无因次产量积分导数qDid的对数作纵坐标，以无因次时间tD的对数为横坐标绘制，在图4中，ωv＝0.1，ωf＝0.01，ωV＝0.2，RD＝20，reD＝1000000，λm＝10-9，λv＝10-7。
从图4可以看出，单洞与三重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版可指导目标单洞与多重介质复合型储层的流体流动规律分析，其流体流动规律可以分为以下七个阶段：
阶段Ⅰ：为大型溶洞反应阶段，相当于一个扩大的“井筒”，表现出井筒储集效应，无因次产量积分与产量积分导数曲线重合为斜率为“-1”的直线；
阶段Ⅱ：为外围三重介质裂缝线性流与径向流反应阶段，主要与多重介质中裂缝的性质有关；
阶段Ⅲ：为外围三重介质中溶蚀孔洞向裂缝窜流阶段，无因次产量积分导数曲线出现第一个“凹子”；
阶段Ⅳ：为三重介质裂缝与溶蚀孔洞整体的径向流阶段，无因次产量积分与产量积分导数曲线平行；
阶段Ⅴ：为三重介质中基质向裂缝窜流阶段，无因次产量积分导数曲线出现第二个“凹子”；
阶段Ⅵ：为外围三重介质整体径向流阶段；
阶段Ⅶ：为边界控制流阶段，对于封闭边界，产量积分曲线与产量积分导数曲线重合为斜率为“-1”的直线。
在步骤105中，要进行敏感性分析，其方法有多种，下面给出其中一个实施例。
在一实施例中，从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，包括：
从物性参数中确定待分析的关键物性参数；
从油井产量不稳定分析曲线图版中，分解出关键物性参数对油井产量不稳定影响分析曲线图版；
分析每个关键物性参数对油井产量不稳定影响分析曲线图版，在确定该关键物性参数为敏感性物性参数时，分析该关键物性参数对油井产量的影响。
以图4为例，基于该图进行敏感性分析时，首先确定待分析的关键物性参数为无量纲大型溶洞半径RD、大型溶洞的弹性储容比ωV和无量纲泄油半径reD，然后从油井产量不稳定分析曲线图版中，分解出这3个关键物性参数对油井产量不稳定影响分析曲线图版，其中，图5为本发明实施例中无量纲大型溶洞半径对油井产量不稳定影响分析曲线图版的示意图，图6为本发明实施例中大型溶洞的弹性储容比对油井产量不稳定影响分析曲线图版的示意图，图7为本发明实施例中无量纲泄油半径对油井产量不稳定影响分析曲线图版的示意图，从图5中可以看出，无量纲大型溶洞半径RD越大，无因次产量积分与产量积分导数曲线相交的直线越靠上，反映的井筒储集效应阶段持续时间越长。从图6中可以看出，大型溶洞的弹性储容比ωV越大，无因次产量积分与产量积分导数曲线相交的直线越靠上，反映的井筒储集效应阶段持续时间越长，但是ωV的影响比RD小得多，因此，确定无量纲大型溶洞半径RD为敏感性物性参数。从图7中可以看出，无量纲泄油半径reD越大，对流体流动规律的阶段Ⅵ至阶段Ⅶ影响较大，也属于敏感性物性参数。之后，可根据图5和图7的图版，分析目标单洞与多重介质复合型储层的油井产量。
与图4同理，可以得到图8表示的本发明实施例中单洞与双重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版的示意图，以及图9表示的本发明实施例中单洞与单重介质复合型储层数学模型对应的油井产量不稳定分析曲线图版的示意图，图8和图9以无因次产量qD的对数、无因次产量积分qDi的对数和无因次产量积分导数qDid的对数作纵坐标，以无因次时间tD的对数为横坐标绘制，在图8中，ωv＝0，ωf＝0.01，ωV＝0.2，RD＝20，reD＝1000000，λm＝10-9，λv＝10-7。在图9中，ωv＝0，ωf＝0.8，ωV＝0.2，RD＝20，reD＝1000000，λm＝10-9，λv＝10-7。图8和图9对应的其流体流动规律和敏感性分析过程与图4类似，这里不再赘述。
综上所述，在本发明实施例提出的方法中，确定目标单洞与多重介质复合型储层的初始物性参数；将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。在上述过程中，本发明实施例采用了哑铃状缝洞型储层的数学模型，该模型准确度高，使得解析获得的油井井底的无因次产量拉氏解的准确度高，通过反演后获得的结果可绘制油井产量不稳定分析曲线图版的准确度高，即最后分解出的不同物性参数对油井产量不稳定影响分析曲线图版的准确度高，从而可以更准确地进行不同物性参数对油井产量的敏感性影响的分析。
本发明实施例还提出一种单洞与多重介质复合型储层参数敏感性分析装置，其原理与单洞与多重介质复合型储层参数敏感性分析方法类似，这里不再赘述。
图10为本发明实施例中单洞与多重介质复合型储层参数敏感性分析装置的示意图，该装置包括：
初始物性参数确定模块1001，用于确定目标单洞与多重介质复合型储层的初始物性参数；
解析模块1002，用于将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；
反演模块1003，用于对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；
图版绘制模块1004，用于基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；
敏感性分析模块1005，用于从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。
在一实施例中，单洞与多重介质复合型储层数学模型采用如下公式表示：
其中，kf、km、kv分别为裂缝、基质和溶蚀孔洞的渗透率，μm2；φf、φmφv分别为裂缝、基质和溶蚀孔洞的孔隙度，无因次；φV为大型溶洞的孔隙度，无因次；Ctf、Ctm、Ctv分别为裂缝、基质和溶蚀孔洞的总压缩系数，MPa-1；CtV为大型溶洞的总压缩系数，MPa-1；pf、pm、pv分别为裂缝、基质和溶蚀孔洞在某一时刻某一点的地层压力，MPa；pV为大型溶洞的地层压力，MPa；pi为原始地层压力，MPa；ωf、ωm、ωv分别为裂缝、基质和溶蚀孔洞的弹性储容比，无因次；ωV为大型溶洞的弹性储容比，无因次；λm、λv分别为基质和溶蚀孔洞的窜流系数，无因次；αm、αv分别为基质、溶蚀孔洞的形状因子，无因次；R为大型溶洞的半径，m；r为任一点到井的距离，m；rw为井半径，m；re为泄油半径，m；h为储层厚度，m；μ为原油粘度，mPa·s；B为原油体积系数，无因次；q为油井产量，m3/d；t为生产时间，d；pfD,mD,vD,VD分别为无量纲基质压力、无量纲裂缝压力、无量纲溶蚀孔洞压力和无量纲大型溶洞压力；tD为无量纲生成时间；rD为任一点到井的无量纲距离；RD为无量纲大型溶洞半径；reD为无量纲泄油半径；pwD为无量纲井底流压。
在一实施例中，所述油井产量不稳定分析曲线图版包括无因次产量随无因次时间变化曲线、无因次产量积分随无因次时间变化曲线、无因次产量积分导数随无因次时间变化曲线。
在一实施例中，敏感性分析模块具体用于：
从物性参数中确定待分析的关键物性参数；
从油井产量不稳定分析曲线图版中，分解出关键物性参数对油井产量不稳定影响分析曲线图版；
分析每个关键物性参数对油井产量不稳定影响分析曲线图版，在确定该关键物性参数为敏感性物性参数时，分析该关键物性参数对油井产量的影响。
综上所述，在本发明实施例提出的装置中，确定目标单洞与多重介质复合型储层的初始物性参数；将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响。在上述过程中，本发明实施例采用了哑铃状缝洞型储层的数学模型，该模型准确度高，使得解析获得的油井井底的无因次产量拉氏解的准确度高，通过反演后获得的结果可绘制油井产量不稳定分析曲线图版的准确度高，即最后分解出的不同物性参数对油井产量不稳定影响分析曲线图版的准确度高，从而可以更准确地进行不同物性参数对油井产量的敏感性影响的分析。
本申请的实施例还提供一种计算机设备，图11为本发明实施例中计算机设备的示意图，该计算机设备能够实现上述实施例中的单洞与多重介质复合型储层参数敏感性分析方法中全部步骤，所述电子设备具体包括如下内容：
处理器(processor)1101、存储器(memory)1102、通信接口(CommunicationsInterface)1103和总线1104；
其中，所述处理器1101、存储器1102、通信接口1103通过所述总线1104完成相互间的通信；所述通信接口1103用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输；
所述处理器1101用于调用所述存储器1102中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的单洞与多重介质复合型储层参数敏感性分析方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质，能够实现上述实施例中的单洞与多重介质复合型储层参数敏感性分析方法中全部步骤，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的单洞与多重介质复合型储层参数敏感性分析方法的全部步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种单洞与多重介质复合型储层参数敏感性分析方法，其特征在于，包括：
确定目标单洞与多重介质复合型储层的初始物性参数；
将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；
对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；
基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；
从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响；
所述油井产量不稳定分析曲线图版包括无因次产量随无因次时间变化曲线、无因次产量积分随无因次时间变化曲线、无因次产量积分导数随无因次时间变化曲线；
从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，包括：
从物性参数中确定待分析的关键物性参数；
从油井产量不稳定分析曲线图版中，分解出关键物性参数对油井产量不稳定影响分析曲线图版；
分析每个关键物性参数对油井产量不稳定影响分析曲线图版，在确定该关键物性参数为敏感性物性参数时，分析该关键物性参数对油井产量的影响。
2.一种单洞与多重介质复合型储层参数敏感性分析装置，其特征在于，包括：
初始物性参数确定模块，用于确定目标单洞与多重介质复合型储层的初始物性参数；
解析模块，用于将所述初始物性参数输入至单洞与多重介质复合型储层数学模型中，解析获得目标单洞与多重介质复合型储层的油井井底的无因次产量拉氏解；
反演模块，用于对所述无因次产量拉氏解进行数值反演，获得无因次产量、无因次产量积分、无因次产量积分导数、无因次时间；
图版绘制模块，用于基于无因次产量、无因次产量积分、无因次产量积分导数、无因次时间，绘制油井产量不稳定分析曲线图版；
敏感性分析模块，用于从油井产量不稳定分析曲线图版中，分解出不同物性参数对油井产量不稳定影响分析曲线图版，所述不同物性参数对油井产量不稳定影响分析曲线图版用于分析不同物性参数对油井产量的敏感性影响；
所述油井产量不稳定分析曲线图版包括无因次产量随无因次时间变化曲线、无因次产量积分随无因次时间变化曲线、无因次产量积分导数随无因次时间变化曲线；
敏感性分析模块具体用于：
从物性参数中确定待分析的关键物性参数；
从油井产量不稳定分析曲线图版中，分解出关键物性参数对油井产量不稳定影响分析曲线图版；
分析每个关键物性参数对油井产量不稳定影响分析曲线图版，在确定该关键物性参数为敏感性物性参数时，分析该关键物性参数对油井产量的影响。
3.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1所述方法。
4.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011380695.2A CN112377185B (zh) | 2020-11-30 | 2020-11-30 | 单洞与多重介质复合型储层参数敏感性分析方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011380695.2A CN112377185B (zh) | 2020-11-30 | 2020-11-30 | 单洞与多重介质复合型储层参数敏感性分析方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112377185A CN112377185A (zh) | 2021-02-19 |
| CN112377185B true CN112377185B (zh) | 2023-08-22 |
Family
ID=74590257
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011380695.2A Active CN112377185B (zh) | 2020-11-30 | 2020-11-30 | 单洞与多重介质复合型储层参数敏感性分析方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112377185B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN110532592A (zh) * | 2019-07-15 | 2019-12-03 | 西南石油大学 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
-
2020
- 2020-11-30 CN CN202011380695.2A patent/CN112377185B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN110532592A (zh) * | 2019-07-15 | 2019-12-03 | 西南石油大学 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112377185A (zh) | 2021-02-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3298232A1 (en) | Parallel solution or fully-coupled fully-implicit wellbore modeling in reservoir simulation | |
| CN105178952B (zh) | 确定水平井人工裂缝间距的方法及装置 | |
| CN114169263A (zh) | 基于压裂液返排数据的页岩油气藏裂缝参数反演方法 | |
| CN112069654B (zh) | 一种碳酸盐岩酸化数值模拟方法 | |
| Liu et al. | A simplified and efficient method for water flooding production index calculations in low permeable fractured reservoir | |
| CN112377182B (zh) | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
| CN107451671B (zh) | 用于预测页岩地层压裂后初始产能的方法及系统 | |
| CN114547953A (zh) | 一种基于优化设计图版的压裂施工参数优化方法及系统 | |
| 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 | |
| CN112377185B (zh) | 单洞与多重介质复合型储层参数敏感性分析方法及装置 | |
| EP2678718A2 (en) | Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region | |
| CN111305805B (zh) | 储层裂缝渗吸质量预测方法及系统 | |
| CN110486008B (zh) | 一种径向复合油藏的参数解释方法及系统 | |
| CN112377184B (zh) | 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置 | |
| CN112377183B (zh) | 多洞多缝型碳酸盐岩储层产量分析模型构建方法及装置 | |
| CN115310379A (zh) | 井间干扰条件下压裂水平井的生产动态分析方法及设备 | |
| CN112682033B (zh) | 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
| CN113887110A (zh) | 油井产能预测方法、装置、电子设备及计算机存储介质 | |
| CN104182558A (zh) | 一种缝洞野外露头水驱油数值模拟方法 | |
| CN112377186B (zh) | 哑铃状缝洞型碳酸盐岩储层产量分析模型构建方法及装置 | |
| CN114492219A (zh) | 用于获取缝洞油藏流动特征表征信息的方法及系统 | |
| CN114086933B (zh) | 一种试井解释方法、装置、存储介质以及计算机设备 | |
| CN112419493A (zh) | 页岩储层三维属性模型建立方法及装置 | |
| CN114086934B (zh) | 一种试井解释方法、装置、存储介质以及计算机设备 | |
| CN112377181B (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 |