本发明涉及油气田开发领域，尤其涉及一种储层产能的预测方法及装置。
背景技术：
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
产能评价是油气田开发方案设计和油气井合理配产的重要基础工作。目前通用的产能评价方法是依靠现场获取的动静态资料紧密结合来进行油气层产能的预测。其中，静态资料主要是地震、岩心和测井等资料，地震资料注重区域层面的考虑，从开发的角度来看，其对单井的描述往往显得分辨率过低，而岩心、测井资料只能反映井壁附近很小区域内的地层特征，特别是对于缝洞发育的储层更是难以代表实际情况，这些静态资料所能提供的信息只占储层整体信息量的数亿分之一；动态资料主要是试井资料，试井资料综合反映了流动所波及到的足够大地层范围内的流动特征，但是试井解释面临着多解性强和对压力的精度要求过高的限制，而这一问题是高度复杂的碳酸盐岩储层中更是不容忽视。
现有储层(例如，碳酸盐岩气藏储层)产能评价是在储层特征认识的基础上，围绕缝洞型储层在试井资料基础上，采用多重介质模型进行气井产能评价。这种产能评价方法在进行准确产能预测之前需要明确两方面的资料，一方面需要弄清储层发育情况(包括储层类型、储集空间及其连通性等)，另一方面需要进行不稳定试井，从而建立准确的多重介质或者是单重介质模型进行产能评价。由此可以看出，目前的储层产能评价方法对于资料要求高，而且必须在钻井完成后，通过测井和试气工作，以及测井解释和模型建立，才能获得准确的储层产能数据。
技术实现要素：
本发明实施例提供一种储层产能的预测方法，用以解决现有储层产能评估方法对资料要求高，且无法对开发建产阶段的储层产能进行评估的技术问题，该方法包括：获取储层的上覆灰岩厚度，其中，上覆灰岩厚度为储层在钻井过程中实时获取到的上覆灰岩厚度；基于上覆灰岩厚度与无阻流量的关系模型，根据储层的上覆灰岩厚度确定储层的无阻流量，其中，储层的无阻流量用于表征储层的产能，上覆灰岩厚度与无阻流量的关系模型是预先根据已开发储层的上覆灰岩厚度与无阻流量建立的关系模型。
本发明实施例还提供一种储层产能的预测装置，用以解决现有储层产能评估方法对资料要求高，且无法对开发建产阶段的储层产能进行评估的技术问题，该装置包括：待开发储层上覆灰岩厚度获取模块，用于获取储层的上覆灰岩厚度，其中，上覆灰岩厚度为储层在钻井过程中实时获取到的上覆灰岩厚度；无阻流量预测模块，用于基于上覆灰岩厚度与无阻流量的关系模型，根据储层的上覆灰岩厚度确定储层的无阻流量，其中，储层的无阻流量用于表征储层的产能，上覆灰岩厚度与无阻流量的关系模型是预先根据已开发储层的上覆灰岩厚度与无阻流量建立的关系模型。
本发明实施例还提供一种计算机设备，用以解决现有储层产能评估方法对资料要求高，且无法对开发建产阶段的储层产能进行评估的技术问题，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述的储层产能的预测方法。
本发明实施例还提供一种计算机可读存储介质，用以解决现有储层产能评估方法对资料要求高，且无法对开发建产阶段的储层产能进行评估的技术问题，所述计算机可读存储介质存储有执行上述储层产能的预测方法的计算机程序。
本发明实施例中，通过预先构建的上覆灰岩厚度与无阻流量的关系模型，可以根据钻井过程中获取到的待开发储层的上覆灰岩厚度，确定待开发储层的无阻流量，由于无阻流量可以表征待开发储层的产能大小，因而，通过本发明实施例，提供了一种无需钻井完成便可预测储层产能的方案，实现了在开发建产阶段对储层产能进行快速预测的技术效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中提供的一种储层产能的预测方法流程图；
图2为本发明实施例中提供的一种上浮灰岩识别结果示意图；
图3为本发明实施例中提供的一种上浮灰岩厚度与无阻流量的关系示意图；
图4为本发明实施例中提供的一种待开发储层无阻流量预测结果示意图；
图5为本发明实施例中提供的一种储层产能的预测装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
本发明实施例中还提供了一种储层产能的预测方法，图1为本发明实施例中提供的一种储层产能的预测方法流程图，如图1所示，该方法包括如下步骤：
s101，获取待开发储层的上覆灰岩厚度，其中，上覆灰岩厚度为待开发储层在钻井过程中实时获取到的上覆灰岩厚度；
s102，基于上覆灰岩厚度与无阻流量的关系模型，根据待开发储层的上覆灰岩厚度确定待开发储层的无阻流量，其中，待开发储层的无阻流量用于表征待开发储层的产能，上覆灰岩厚度与无阻流量的关系模型是预先根据已开发储层的上覆灰岩厚度与无阻流量建立的关系模型。
需要说明的是，上述待开发储层可以是任意一种具有连通孔隙、允许油气在其中储存和渗滤的岩层。在基于上覆灰岩厚度与无阻流量的关系模型，根据待开发储层的上覆灰岩厚度确定待开发储层的无阻流量之前，需要先建立上覆灰岩厚度与无阻流量建立的关系模型，建立关系模型的具体步骤为：获取已开发储层的上覆灰岩厚度；获取已开发储层的无阻流量；根据已开发储层的上覆灰岩厚度和无阻流量，构建上覆灰岩厚度与无阻流量的关系模型。
由上可知，在本发明上述实施例中，通过预先构建的上覆灰岩厚度与无阻流量的关系模型，在获取到待开发储层在钻井过程中的上覆灰岩厚度后，根据待开发储层的上覆灰岩厚度确定待开发储层的无阻流量，由于无阻流量可以表征待开发储层的产能大小，因而，通过本发明实施例，提供了一种无需钻井完成便可预测储层产能的方案，实现了在开发建产阶段对储层产能进行快速预测的技术效果。
作为一种可选的实施例，本发明实施例以叠合岩溶储层为例来进行说明。叠合岩溶储层是指由于区域构造抬升、局部构造抬升、海平面下降以及深部热液流体窜入导致碳酸盐岩地层遭受持续时间不同、影响范围不一多期大气淡水林滤或者是深部热液溶蚀，形成孔缝洞发育尺度、缝洞体规模等存在明显叠合特征的岩溶型储层。叠合岩溶包括构造不整合岩溶(表生期岩溶)、层内准同生岩溶及同生岩溶。
叠合岩溶受两期构造不整合影响，在平面上不同的区域两期风化壳溶蚀程度存在差异，在局部地区存在两期不整合残留地层。本发明实施例针对两期构造不整合风化壳形成的叠合岩溶，可以提供一种预测叠合岩溶储层气井产能的方法，通过识别气藏早期(即气田开发方案建立前的阶段，也即在主要依靠探井和评价井论证气田地质特征和开发特征的阶段)探井及评价井叠合岩溶风化壳储层上覆灰岩，建立上覆灰岩厚度与早期探井和评价井无阻流量的关系图版，依据新钻井上覆灰岩厚度可以快速判别该井无阻流量。受叠合岩溶特征影响，不同区域发育溶蚀孔隙、溶蚀扩大缝、高角度缝等储集空间特征存在较大差异，导致气藏产层中缝洞型、孔洞型及孔隙型储层构型厚度及各类储层构型所占比例不同，从而影响气井产能。因此，上覆灰岩厚度可以间接反映下伏地层有效储层发育情况，从而通过对上覆灰岩厚度与完钻井无阻流量的关系研究，建立气藏开发区上覆灰岩厚度与气井产能的关系图版(用于反映即上覆灰岩厚度与无阻流量的关系)。
气藏开发建产过程中，依据完钻井高分辨率地震资料可以进行开发建产区筛选、井位优选和轨迹优化，之后可以进行开发井钻井。在开发井钻遇风化壳储层之后，可以依据随钻测井及岩屑分析可以快速分析出风化壳储层上覆灰岩厚度，结合已经建立的上覆灰岩厚度与无阻流量关系图版可以快速判别该井无阻流量。
上述s101可以根据待开发储层的钻井、随钻测井等资料，可以识别上覆灰岩地层，由于灰岩地层的岩屑、测井等特征比较明显，依靠测井、岩屑录井等技术手段准确识别灰岩，是已经十分成熟的技术。因而，本发明实施例对于灰岩的识别不做重要阐述。
本发明实施例选取四川盆地安岳气田震旦系磨溪区块叠合岩溶为例来对本发明进行详细说明，磨溪区块震旦系气藏目前处于产能建设阶段，目前25口完钻井(直井和大斜度井)完成测试，无阻流量介于2～217.58×104m3/d，平均40×104m3/d，气井测试产量差异较大。气井产能控制因素的研究相对复杂，由于储集空间相对复杂，气井产能的评价也存在较大的难度，同时评价结果也存在误差。由于气藏开发建产过程中要依据气井产能来进行地面官网建设，也就是地面管网的建设很大程度上依据不同区域气井产能的大小来进行规划和设计。因此进行快速气井产能预测对于该类气藏建产节奏的加快以及地面建设的设计具有非常重要的意义。
通过本发明实施例提供的储层产能的预测方法，首先识别早期探井及评价井叠合岩溶风化壳储层上覆灰岩，上覆灰岩识别结果如图2所示，上扬子地区在晚震旦世灯影组沉积期及早寒武世麦地坪组沉积期发生的3幕地壳升降运动，每幕运动均导致地层的抬升、剥蚀，形成风化壳不整合面。桐湾运动ⅰ幕发生在灯影组灯二段沉积末期，表现为灯三段区域性碎屑岩假整合于灯二段白云岩之上。桐湾运动ⅱ幕发生在灯影组沉积期末，表现为灯影组与下寒武统麦地坪组假整合接触。桐湾运动ⅲ幕发生在早寒武世麦地坪组沉积期末，表现为下寒武统麦地坪组与筇竹寺组假整合接触。三幕桐湾运动主要表现为升降运动，形成了3个区域性不整合面，分别位于震旦系灯影组灯三段与灯二段之间、灯影组与下寒武统麦地坪组之间、麦地坪组与筇竹寺组之间。桐湾运动ⅰ幕导致灯二段岩溶储层的发育，桐湾运动ⅱ幕和ⅲ幕在区域上大部分具有叠合特征，导致灯四段岩溶储层的发育。两幕不整合面之间为残存麦地坪组灰岩段，依据灰岩、云岩以及泥岩之间的测井特征可以有效识别灰岩厚度。
由于灰岩和云岩均具有较低的gr值(一种测井值，也叫伽马曲线，可以反映地层泥质含量的大小)，但是灰岩pe值(一种测井曲线值，叫做光电吸收截面值，可以用来判别灰岩和云岩)要高于白云岩，即可据此有效判别上覆灰岩。
在灰岩厚度识别基础上，可以依据完钻井试气资料，对已完钻井上覆灰岩厚度与其计算出的无阻流量值进行拟合，可以得到上覆灰岩厚度与测试井无阻流量之间的关系图版，如图3所示，灰岩厚度与无阻流量之间的关系图版具有明显的分段特征：
ⅰ段(灰岩厚度0米)：
ⅱ段(灰岩厚度介于0～5米)，拟合曲线为：
y＝1417.6e-1.8811x，拟合曲线的相关系数r2＝0.786；
ⅲ段(灰岩厚度5～20米)，拟合曲线为：
y＝1980.5e-0.418x，拟合曲线的相关系数r2＝0.8991
ⅳ段(灰岩厚度大于20米)，拟合曲线为：
y＝5269.7e-0.17x，拟合曲线的相关系数r2＝0.9846；
其中，y表示储层的无阻流量，单位为104m3/d；x表示储层的上覆灰岩厚度，单位为m。
由上可以看出，各段拟合方程相关系数较高(拟合曲线的相关系数，r2＝1代表最好，r2＝0代表最差，一般大于0.7即为很好)，具有较好的示范意义。因而，作为一种可选的上述s102可以包括如下至少一种情形：
第一种情形，如果待开发储层的上覆灰岩厚度在0～5米的范围内，则通过如下关系模型计算待开发储层的无阻流量：y＝1417.6e-1.8811x；
第二种情形，如果待开发储层的上覆灰岩厚度在5～20米的范围内，则通过如下关系模型计算待开发储层的无阻流量：y＝1980.5e-0.418x；
第三种情形，如果待开发储层的上覆灰岩厚度大于20米的范围，则通过如下关系模型计算待开发储层的无阻流量：y＝5269.7e-0.17x。
对于气藏开发井，在钻井的过程中可以依据钻井、测井资料快速判别上覆灰岩厚度，将该厚度值带入所建立模板，即可快速得到该井无阻流量。图4为本发明实施例中提供的一种待开发储层无阻流量预测结果示意图，如图4所示，ii类是大于0小于等于5米灰岩井，iii类是大于5米小于等于20米灰岩井，ⅳ类是大于20米灰岩井。需要指出的是，目前该方法对于上覆灰岩为0的井(ⅰ段)因不能建立相对应的拟合方程，尚不能快速预测其无阻流量，只能钻完井之后，依靠传统方法才能进行产能评价。
以新钻井m井为例，该井在钻达目的层之后，依据岩屑及随钻测井资料判别该井上覆地层灰岩厚度为7.2米，依据所建立关系图版发现该井处于ⅲ段，将该值带入ⅲ段的关系模型，可得气井无阻流量为97.7×104m3/d。m井完钻之后，获取较准确测井和试气资料基础上，经二项式产能方程得到该井无阻流量为102.5×104m3/d，误差为4.7％。由此可见，依靠该方法相对于传统方法误差在可接受范围之内。
随着完钻井资料的增大，通过对关系图版进行优化，可以不断提高预测产能的准确率。因而，为了提高本发明实施例储层长能预测的准确率，一种可选的实施例中，上述方法还可以包括：根据待开发储层完成钻井后的上覆灰岩厚度与无阻流量，更新上覆灰岩厚度与无阻流量的关系模型。
与现有的储层产能评价方法相比，本发明实施例提供的储层产能的预测方案，只需要上覆灰岩厚度这一静态资料，对于资料要求比较少。另外，在气井尚未完钻、也无需进行测井和试气工作，即可在早期钻遇目的层后进行快速产能评价，工作量小，简单易行。
本发明实施例中还提供了一种储层产能的预测装置，如下面的实施例所述。由于这些模块解决问题的原理与储层产能的预测方法相似，因此这些模块的实施可以参见方法的实施，重复之处不再赘述。
图5为本发明实施例中提供的一种储层产能的预测装置示意图，如图5所示，该装置包括：
待开发储层上覆灰岩厚度获取模块51，用于获取待开发储层的上覆灰岩厚度，其中，上覆灰岩厚度为待开发储层在钻井过程中实时获取到的上覆灰岩厚度；
待开发储层无阻流量预测模块52，用于基于上覆灰岩厚度与无阻流量的关系模型，根据待开发储层的上覆灰岩厚度确定待开发储层的无阻流量，其中，待开发储层的无阻流量用于表征待开发储层的产能，上覆灰岩厚度与无阻流量的关系模型是预先根据已开发储层的上覆灰岩厚度与无阻流量建立的关系模型。
由上可知，在本发明上述实施例中，通过预先构建的上覆灰岩厚度与无阻流量的关系模型，在通过待开发储层上覆灰岩厚度获取模块51获取到待开发储层在钻井过程中的上覆灰岩厚度后，通过待开发储层无阻流量预测模块52基于预先构建的上覆灰岩厚度与无阻流量的关系模型，根据待开发储层的上覆灰岩厚度确定待开发储层的无阻流量，由于无阻流量可以表征待开发储层的产能大小，因而，通过本发明实施例，提供了一种无需钻井完成便可预测储层产能的方案，实现了在开发建产阶段对储层产能进行快速预测的技术效果。
在一种可选的实施例中，上述待开发储层无阻流量预测模块52包括：第一预测子块521，用于如果待开发储层的上覆灰岩厚度在0～5米的范围内，则通过如下关系模型计算待开发储层的无阻流量：y＝1417.6e-1.8811x；第二预测模块522，用于如果待开发储层的上覆灰岩厚度在5～20米的范围内，则通过如下关系模型计算待开发储层的无阻流量：y＝1980.5e-0.418x；第三预测模块523，用于如果待开发储层的上覆灰岩厚度大于20米的范围，则通过如下关系模型计算待开发储层的无阻流量：y＝5269.7e-0.17x；其中，y表示储层的无阻流量；x表示储层的上覆灰岩厚度。
在一种可选的实施例中，上述装置还可以包括：关系模型更新模块53，用于根据待开发储层完成钻井后的上覆灰岩厚度与无阻流量，更新上覆灰岩厚度与无阻流量的关系模型。
在一种可选的实施例中，上述装置还可以包括：已开发储层上覆灰岩厚度获取模块54，用于获取已开发储层的上覆灰岩厚度；已开发储层无阻流量获取模块55，用于获取已开发储层的无阻流量；关系模型构建模块56，用于根据已开发储层的上覆灰岩厚度和无阻流量，构建上覆灰岩厚度与无阻流量的关系模型。
基于上述任意一种可选的或优选的装置实施例，作为一种可选的实施例，上述装置实施例中的待开发储层可以为叠合岩溶储层。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法实施例中任意一种可选的或优选的储层产能的预测方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法实施例中任意一种可选的或优选的储层产能的预测方法的计算机程序。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
1.一种储层产能的预测方法，其特征在于，包括：
获取待开发储层的上覆灰岩厚度，其中，所述上覆灰岩厚度为所述待开发储层在钻井过程中实时获取到的上覆灰岩厚度；
基于上覆灰岩厚度与无阻流量的关系模型，根据所述待开发储层的上覆灰岩厚度确定所述待开发储层的无阻流量，其中，所述待开发储层的无阻流量用于表征所述待开发储层的产能，所述上覆灰岩厚度与无阻流量的关系模型是预先根据已开发储层的上覆灰岩厚度与无阻流量建立的关系模型。
2.如权利要求1所述的方法，其特征在于，基于上覆灰岩厚度与无阻流量的关系模型，根据所述待开发储层的上覆灰岩厚度确定所述待开发储层的无阻流量，包括：
如果所述待开发储层的上覆灰岩厚度在0～5米的范围内，则通过如下关系模型计算所述待开发储层的无阻流量：y＝1417.6e-1.8811x；
如果所述待开发储层的上覆灰岩厚度在5～20米的范围内，则通过如下关系模型计算所述待开发储层的无阻流量：y＝1980.5e-0.418x；
如果所述待开发储层的上覆灰岩厚度大于20米的范围，则通过如下关系模型计算所述待开发储层的无阻流量：y＝5269.7e-0.17x；
其中，y表示储层的无阻流量；x表示储层的上覆灰岩厚度。
3.如权利要求1所述的方法，其特征在于，所述方法还包括：
根据所述待开发储层完成钻井后的上覆灰岩厚度与无阻流量，更新所述上覆灰岩厚度与无阻流量的关系模型。
4.如权利要求1所述的方法，其特征在于，在基于上覆灰岩厚度与无阻流量的关系模型，根据所述待开发储层的上覆灰岩厚度确定所述待开发储层的无阻流量之前，所述方法还包括：
获取已开发储层的上覆灰岩厚度；
获取所述已开发储层的无阻流量；
根据所述已开发储层的上覆灰岩厚度和无阻流量，构建所述上覆灰岩厚度与无阻流量的关系模型。
5.如权利要求1至4中任一所述的方法，其特征在于，所述待开发储层为叠合岩溶储层。
6.一种储层产能的预测装置，其特征在于，包括：
待开发储层上覆灰岩厚度获取模块，用于获取待开发储层的上覆灰岩厚度，其中，所述上覆灰岩厚度为所述待开发储层在钻井过程中实时获取到的上覆灰岩厚度；
待开发储层无阻流量预测模块，用于基于上覆灰岩厚度与无阻流量的关系模型，根据所述待开发储层的上覆灰岩厚度确定所述待开发储层的无阻流量，其中，所述待开发储层的无阻流量用于表征所述待开发储层的产能，所述上覆灰岩厚度与无阻流量的关系模型是预先根据已开发储层的上覆灰岩厚度与无阻流量建立的关系模型。
7.如权利要求6所述的装置，其特征在于，所述待开发储层无阻流量预测模块包括：
第一预测子块，用于如果所述待开发储层的上覆灰岩厚度在0～5米的范围内，则通过如下关系模型计算所述待开发储层的无阻流量：y＝1417.6e-1.8811x；
第二预测模块，用于如果所述待开发储层的上覆灰岩厚度在5～20米的范围内，则通过如下关系模型计算所述待开发储层的无阻流量：y＝1980.5e-0.418x；
第三预测模块，用于如果所述待开发储层的上覆灰岩厚度大于20米的范围，则通过如下关系模型计算所述待开发储层的无阻流量：y＝5269.7e-0.17x；
其中，y表示储层的无阻流量；x表示储层的上覆灰岩厚度。
8.如权利要求6所述的装置，其特征在于，所述装置还包括：关系模型更新模块，用于根据所述待开发储层完成钻井后的上覆灰岩厚度与无阻流量，更新所述上覆灰岩厚度与无阻流量的关系模型。
9.如权利要求6所述的装置，其特征在于，所述装置还包括：
已开发储层上覆灰岩厚度获取模块，用于获取已开发储层的上覆灰岩厚度；
已开发储层无阻流量获取模块，用于获取所述已开发储层的无阻流量；
关系模型构建模块，用于根据所述已开发储层的上覆灰岩厚度和无阻流量，构建所述上覆灰岩厚度与无阻流量的关系模型。
10.如权利要求6至9任一所述的装置，其特征在于，所述待开发储层为叠合岩溶储层。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述的储层产能的预测方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述储层产能的预测方法的计算机程序。
技术总结