预测致密油待评价区域可采储量的方法及装置的制造方法
预测致密油待评价区域可采储量的方法及装置的制造方法
【技术领域】
[0001]本申请涉及致密油开发技术领域，尤其是涉及一种预测致密油待评价区域可采储 量的方法及装置。
【背景技术】
[0002] 致密油资源潜力巨大，目前且尚处于勘探初期，是未来石油产量的重要增长点。目 前致密油可采储量评价体系仍以常规油气评价为主。目前涉及致密油储量评价的方法主要 技术有:容积法、概率统计法、类比法、井距外推法、产量递减法等。其中：
[0003] ①容积法，基于测井和地震资料的刻画出待评价区地下致密油的空间分布情况， 计算出致密油层的体积，通过油层的孔隙度估算出致密油的储量。
[0004] ②概率统计法，利用生产井最终估算可采储量(Estimated Ultimate Recovery， EUR)，建立生产井EUR概率分布图，确定待评价区的钻井数N以后，从生产井EUR概率分布图 上通过蒙特卡洛随机模拟对应得到单井EUR，选取N次，得到N个单井EUR值，求得平均值然后 输入。重复该过程5000到10000次，得到5000-10000个平均值，根据这些数值做出新的EUR概 率统计分布图，根据不同概率处获得未钻探可采储量。
[0005] ③类比法，利用在产致密油区的地质条件和开发井特征，建立刻度区，通过待评价 区地质条件，类比已有刻度区，确定待评价区开发井特征，进而推算出待开发区的储量。
[0006] ④外推井距法，利用已投产生产井的生产曲线，计算出单井的EUR，即为1P级储量。 对井周围按照一个井距进行外推，外推区域的单井EUR与已钻井EUR相当，即生产井一个井 距以内区域的预测产量为2P储量。继续向外推第二个井距的产量，仍认为与已投产井的产 量相当，即生产井2个井距以内区域的预测产量，为3P储量。
[0007] ⑤递减曲线法，利用具有半年以上生产历史的生产井，通过产量递减模型，计算出 来的单井最终估算可采储量。
[0008] 然而，容积法基于常规油气地下聚集模式来评价致密油储量，对于均质厚层致密 油层有一定适用性，但对于大多数致密油层而言，均具有强非均质性，横向纵向变化都很 大，其地质特征参数几乎无法取准，因而最终评价结果精度极低，该方法仅适用于评价致密 油资源量。
[0009] 概率统计法是基于数学统计的方法来计算致密油储量，其中未考虑地质因素的影 响。但是，地质因素的差异对于任何地区油气藏储量都会有非常大的影响。尤其是可以反映 生油能力的烃源条件和反映储集能力的储集条件。另外，概率统计法也未考虑工程因素的 影响，例如水平井段长度以及压裂后的泄流半径，这些因素也会直接造成单井产量大幅变 化。因而，概率统计法忽视了地质因素，无法获得较为精确的致密油储量。
[0010]类比法需要大量的与被评价区域类似的对比数据，即刻度区数据。在类比过程中， 需要通过给评价区各个地质要素进行打分，来给出类比相似系数，以确定被评价区与刻度 区的相似度。但这个确定类比相似系数的过程受到评价人主观影响较大，结果也因此差异 明显。另外，若刻度区与被评价区不属于同一致密岩层，其类比预测精度本身也会比较低。
[0011] 外推井距法，同样缺乏对地质因素差异性的考虑，其默认致密油生产井周边2个井 距内的地质条件几乎一致，未来新钻生产井的产能也会基本相同；但事实上，因为地质因素 总会发生变化的，在已投产生产井的一侧地质因素可能更好，另一侧可能更差。最重要的 是，在致密油井较少地区，外推井距法会极大的低估了致密油的实际可采储量。对于目前大 部分尚处于勘探开发初期的致密油开采区，该方法评价出的储量误差极大。因此，外推井距 法，并不能很好的评价致密油储量。
[0012] 递减曲线法，受控于生产井的生产历史数据且仅能对单井的可采储量进行预测。 但在致密油新区块的勘探过程中，是不具备开发井及生产历史资料的，所以该方法仅能应 用于对已开发的致密油区进行评价，无法对勘探阶段的致密油区块进行预测。
[0013] 综上所述，由于致密油资源与常规油气资源在地质条件、开发手段等方面具有显 著的差异，现有的技术方案均难以获得较为准确的致密油储量。
【发明内容】
[0014] 本申请实施例的目的在于提供一种预测致密油待评价区域可采储量的方法及装 置
[0015] 为达到上述目的，一方面，本申请实施例提供了一种预测致密油待评价区域可采 储量的方法，包括以下步骤：
[0016] 将致密油区域中已有生产井区域划分为若干区块；
[0017]确定每个区块的地质特征参数；
[0018] 依据每个区块内各生产井的历史生产数据对应确定该区块的各生产井的单井可 采储量，并据此确定每个区块的单井平均可采储量；
[0019] 利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储量正相 关的地质特征影响参数；
[0020] 根据各个区块的单井平均可采储量与其地质特征影响参数建立多元线性模型；
[0021] 将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归分 析，以确定所述多元线性模型的模型系数，从而获得所述致密油区域的单井平均可采储量 预测模型；
[0022] 根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采 储量。
[0023] 本申请实施例的预测致密油待评价区域可采储量的方法，每个区块内的各个生产 井处于同一构造带且每个区块具有至少两口达到设定生产历史年限的生产井。
[0024] 本申请实施例的预测致密油待评价区域可采储量的方法，在所述利用相关性分析 法从每个区块的地质特征参数选出与其单井平均可采储量正相关的地质特征影响参数中：
[0025] 当区块的地质特征参数对该区块的单井平均可采储量的影响呈抛物线型关系时， 根据公式
·将该地质特征参数对该区块的单井平均可采储量的影响变换成 正相关关系，其中，X为变换后的地质特征参数，X〇为原地质特征参数，Xm为地质特征参数的 最优经验值。
[0026] 本申请实施例的预测致密油待评价区域可采储量的方法，所述多元线性模型为：
[0027] +名2 +... + <，其中，EUR为致密油区域的单井平均可采储量，e为 自然数，仏为第i个地质特征影响参数，xi为模型系数。
[0028] 本申请实施例的预测致密油待评价区域可采储量的方法，在所述将各个区块的地 质特征影响参数代入所述多元线性模型，并进行多元线性回归分析之前，还包括：
[0029] 根据公式4=^1将所述将各个区块的地质特征影响参数进行去量纲的标准化处 Am 理，其中Ai为标准化后的第i个地质特征影响参数，Ai〇为第i个地质特征影响参数，Aim为第i 个地质特征影响参数的最优经验值；
[0030] 将所述多元线性模型进行两边取对数处理；
[0031 ]对应的，所述将各个区块的地质特征影响参数代入所述多元线性模型，具体为： [0032]将各个区块的经去量纲的标准化处理后的地质特征影响参数代入经两边取对数 处理后的多元线性模型。
[0033]本申请实施例的预测致密油待评价区域可采储量的方法，所述确定每个区块的地 质特征参数，包括：
[0034]利用钻井资料确定各区块内致密油层系的产层厚度；
[0035]利用测井资料和测试分析资料确定各区块内致密油层系的平均孔隙度、产层泥地 比以及供烃源岩的有机质丰度、有机质成熟度；以及，
[0036]利用地层测试资料确定各区块内地层压力系数或压力梯度、气油比、油水比和原 油重度。
[0037] 本申请实施例的预测致密油待评价区域可采储量的方法，所述根据所述单井平均 可采储量预测模型预测所述致密油区域中待评价区域的可采储量，具体包括：
[0038] 获取所述致密油区域中待评价区域的地质特征参数；
[0039] 将所述致密油区域中待评价区域的地质特征参数代入所述单井平均可采储量预 测模型，获得所述致密油区域中待评价区域的单井平均可采储量；
[0040] 获取所述致密油区域中待评价区域的面积，以及所述致密油区域中已有生产井区 域内生产井平均井控面积；
[0041] 根据所述致密油区域中待评价区域的单井平均可采储量、所述致密油区域中待评 价区域的面积、以及所述致密油区域中已有生产井区域内生产井平均井控面积预测所述致 密油区域中待评价区域的可采储量。
[0042] 本申请实施例的预测致密油待评价区域可采储量的方法，在所述根据所述单井平 均可采储量预测模型预测所述致密油区域中待评价区域的可采储量之前，还包括：
[0043] 根据所述单井平均可采储量预测模型及所述致密油区域中各个区块的地质特征 参数，获取所述致密油区域中各个区块的单井平均可采储量预测值；
[0044]判断所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差是否均小于设定误差阈值；
[0045]如果所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差均小于设定误差阈值，则根据所述单井平均可采储量预测模 型预测所述致密油区域中待评价区域的可采储量。
[0046] 本申请实施例的预测致密油待评价区域可采储量的方法，所述相关性分析法包括 Pearson相关性分析法。
[0047] 本申请实施例的预测致密油待评价区域可采储量的方法，所述设定生产历史年限 包括半年。
[0048] 另一方面，本申请实施例还提供了一种预测致密油待评价区域可采储量的装置， 包括：
[0049]有井区域分块模块，用于将致密油区域中已有生产井区域划分为若干区块，其中 每个区块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年 限的生产井；
[0050] 地质特征参数获取模块，用于确定每个区块的地质特征参数；
[0051] 单井储量获取模块，用于依据每个区块内各生产井的历史生产数据对应确定该区 块的各生产井的单井可采储量，并据此确定每个区块的单井平均可采储量；
[0052] 地质特征参数优选模块，用于利用相关性分析法从每个区块的地质特征参数选出 与其单井平均可采储量正相关的地质特征影响参数；
[0053] 预测模型建立模块，用于根据各个区块的单井平均可采储量与其地质特征影响参 数建立多元线性模型；
[0054] 预测模型建立模块，用于将各个区块的地质特征影响参数代入所述多元线性模 型，并进行多元线性回归分析，以确定所述多元线性模型的模型系数，从而获得所述致密油 区域的单井平均可采储量预测模型；
[0055] 待评价区域预测模块，用于根据所述单井平均可采储量预测模型预测所述致密油 区域中待评价区域的可采储量。
[0056] 本申请实施例利用待评价区邻区已投产生产井的生产情况，同时结合地质因素的 变化对致密油单井可采储量进行合理预测，利用相关性分析优选待评价区主控因素，避免 了以往地质预测模型使用相同地质参数来预测致密油单井EUR的不适用问题;而且，还利用 多元回归的方法来确定EUR预测模型中各个地质参数的系数，即解决了人为主观确定相似 系数的误差大的缺点。因此，本申请实施例的评价结果更为精确可靠，为后期致密油区块经 济评价提供了依据。
【附图说明】
[0057] 此处所说明的附图用来提供对本申请实施例的进一步理解，构成本申请实施例的 一部分，并不构成对本申请实施例的限定。在附图中：
[0058] 图1为本申请一实施例的预测致密油待评价区域可采储量的方法的流程图；
[0059] 图2为本申请一实施例的区块划分示意图；
[0060] 图3为本申请一实施例的递减曲线模型；
[0061] 图4为本申请一实施例的模拟EUR与实际EUR对比曲线示意图。
【具体实施方式】
[0062] 为使本申请实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附 图，对本申请实施例做进一步详细说明。在此，本申请实施例的示意性实施例及其说明用于 解释本申请实施例，但并不作为对本申请实施例的限定。
[0063] 下面结合附图，对本申请实施例的【具体实施方式】作进一步的详细说明。
[0064] 参考图1所示，本申请实施例的预测致密油待评价区域可采储量的方法包括以下 步骤：
[0065]步骤S1，将致密油区域中已有生产井区域划分为若干区块。为了有利于达到更为 准确的评价效果每个区块内的各个生产井处于同一构造带(斜坡带、背斜区、凹陷中心等） 且每个区块具有至少两口达到设定生产历史年限（例如半年）的生产井。此外，实验研究表 明，在满足上述要求的前提下，当每个区块的面积大小越相近，每个区块内的钻井分布越相 对聚集且单井产量越相近时，预测出的产量的偶然误差也就越小。
[0066] 步骤S2,确定每个区块的地质特征参数。比如利用钻井资料确定各区块内致密油 层系的产层厚度等；利用测井资料和测试分析资料确定各区块内致密油层系的平均孔隙 度、产层泥地比以及供烃源岩的有机质丰度、有机质成熟度等；以及，利用地层测试资料确 定各区块内地层压力系数或压力梯度、气油比、油水比和原油重度等。
[0067] 步骤S3,依据每个区块内各生产井的历史生产数据对应确定该区块的各生产井的 单井可采储量，并据此确定每个区块的单井平均可采储量。其中，若生产井为水平井，还需 利用公式= + X 对其压裂水平段长进行标准化，然后再进行小区块内的平均，其 中，L为某水平井压裂水平段长度；ζ为研究区标准压裂水平段长度;EUR°为该井原始EUR。 其中，生产井的历史生产数据可以包括产量地建模型。
[0068] 步骤S4,利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储 量正相关的地质特征影响参数。相关性分析法例如可以采用Pearson相关性分析法。
[0069] 其中，当区块的地质特征参数对该区块的单井平均可采储量的影响呈抛物线型关 系时，根据公式
J将该地质特征参数对该区块的单井平均可采储量的影响 变换成正相关关系，其中，X为变换后的地质特征参数，X〇为原地质特征参数，Xm为地质特征 参数的最优经验值。
[0070] 需要说明的是，本申请实施例中，地质特征影响参数实际为一种地质特征参数，由 于该地质特征参数会影响单井平均可采储量，为了区分这种地质特征参数与其他地质特征 参数，在这里将其称为地质特征影响参数。
[0071] 步骤S5,根据各个区块的单井平均可采储量与其地质特征影响参数建立多元线性 模型。其中，多元线性模型还可以参考，多元线性模型具体可为：
[0072] fY +為 12 +…+ ^^其中^皿为致密油区域的单井平均可采储量…为 自然数，仏为第i个地质特征影响参数，xi为模型系数。
[0073] 步骤S6,将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线 性回归分析，以确定所述多元线性模型的模型系数，即确定系数xO~xi的最优解，从而获得 所述致密油区域的单井平均可采储量预测模型。
[0074] 步骤S7,根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域 的可采储量。具体的：
[0075] 获取所述致密油区域中待评价区域的地质特征参数；
[0076] 将所述致密油区域中待评价区域的地质特征参数代入所述单井平均可采储量预 测模型，获得所述致密油区域中待评价区域的单井平均可采储量R;
[0077] 获取所述致密油区域中待评价区域的面积S1，以及所述致密油区域中已有生产井 区域内生产井平均井控面积S2;
[0078] 根据所述致密油区域中待评价区域的单井平均可采储量、所述致密油区域中待评 价区域的面积、以及所述致密油区域中已有生产井区域内生产井平均井控面积预测所述致 密油区域中待评价区域的可采储量，即根据公式可计算出致密油区域中待评价区域 的可采储量。
[0079] 本申请另一实施例中，在所述将各个区块的地质特征影响参数代入所述多元线性 模型，并进行多元线性回归分析之前，还可以包括：
[0080] 根据公式4 将所述将各个区块的地质特征影响参数进行去量纲的标准化处 im 理，其中Ai为标准化后的第i个地质特征影响参数，Ai〇为第i个地质特征影响参数，Aim为第i 个地质特征影响参数的最优经验值；
[0081 ]将所述多元线性模型进行两边取对数处理；
[0082] 对应的，所述将各个区块的地质特征影响参数代入所述多元线性模型，具体为：
[0083] 将各个区块的经去量纲的标准化处理后的地质特征影响参数代入经两边取对数 处理后的多元线性模型。
[0084] 本申请另一实施例中，在所述根据所述单井平均可采储量预测模型预测所述致密 油区域中待评价区域的可采储量之前，还包括：
[0085] 根据所述单井平均可采储量预测模型及所述致密油区域中各个区块的地质特征 参数，获取所述致密油区域中各个区块的单井平均可采储量预测值；
[0086] 判断所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差是否均小于设定误差阈值；
[0087]如果所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差均小于设定误差阈值，则根据所述单井平均可采储量预测模 型预测所述致密油区域中待评价区域的可采储量。 从而可以确保只有准确性满足要求的预 测模型才可以用于预测所述致密油区域中待评价区域的可采储量。
[0088] 虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解， 这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行（例如使用并行 处理器或多线程环境），比如步骤S2和步骤S3之间的顺序就调换，也可以同时进行。
[0089] 下面以北美某区块致密油地质情况为例，介绍本申请实施例的具体应用。
[0090] 首先，针对北美某区块致密油地质情况，确定出致密油主要为源内成藏，未经历明 显侧向运移;根据研究区井位分布情况，结合地质条件、地理特征等，将区块划分为面积相 近的N个小区块(如图2所示）。
[0091] 其次，利用每个小区块钻井、测井、测试分析等资料，确定每个小区块的地质因素 参数，厚度、孔隙度、渗透率、有机质丰度、有机质成熟度、产层泥地比、地层压力系数等。具 体的，利用钻井资料确定各小区块内致密油层系的产层厚度;利用测井资料和测试分析资 料可确定各小区块内致密油层系平均孔隙度、产层泥地比，和供烃源岩的有机质丰度、有机 质成熟度;利用地层测试资料确定地层压力系数(或压力梯度）、气油比、油水比、原油重度 (API)。
[0092]然后，依据每个区块内各生产井的历史生产数据(例如图3所示的产量递减模型） 对应确定该区块的各生产井的单井可采储量，并算出每个小区块的单井平均可采储量。 [0093]其次，利用Pearson相关性分析法确定影响EUR的地质因素，当概率α〈〇. 〇1时，即认 为两者明显相关(表1)。
[0094] 表1研究区地质因素与EUR相关性分析参数表
[0095]
[0096] 当部分地质因素对EUR的影响可能是呈抛物线状型关系，因此要根据研究区公认 的最优条件值选取参数最优值。例如统计发现，研究区产层段泥地比Μ为0.55时为水平井产 量较高，泥地比Μ过大或过小，产量都较低。因此认为泥地比Μ是与EUR成抛物线相关，需要利 用进行变换，变成正相关关系。其中，Μ为变换后的泥地比，M0是原来的泥地 比，Mm为?.
55，埘％区统计规律得到。
[0097]由此，选取影响EUR的主要相关地质因素:产层总厚度Η、泥地比Μ、压力系数Ρ、成熟 度Ro、下T0C、含油饱和度So、孔隙度Φ (见表2)。
[0098] 表2各区块原始地质特征参数表
[0099]
[0100] 再次，已选取地质特征参数与EUR均具有正相关关系，而多个因素均对EUR值有影 响，因此建立单井可采储量EUR与地质因素参数关系模型：
[0101] EURi = X X Mf3 X pp X H X TOCp X Χφ?其中，Η为产层总厚度， Μ为泥地比，Ρ压力系数，Ro成熟度，T0C为下部T0C值，So为含油饱和度So，Φ为孔隙度，xl~ x8分别为各个系数。
[0102] 在将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归 分析之前，将上述公式两边对数处理，可得
[0103]
[0104] 其中，可以先利用公式ln(A:) = h(^x100 )对模型中所有参数进行去量纲的标准 Am 化处理(见表3)，再做上述对数处理(各区块对数变换后地质特征参数见表4)。
[0105] 表3各区块标准化后地质特征参数表
[0106]
[0107] 表4各区块对数变换后地质特征参数表
[0108]
[0109] 然后，将处理过的各个小区块的地质因素参数代入模型，进行多元线性回归分析， 确定系数XI~X8的最优解(见表5)，从而获得所述致密油区域的单井平均可采储量预测模 型：
[0110] EUR = e3.料3 X: H0.1652 X M:0'0615 X 把.1495 X: TOC0'2挪 X 2848 X φ L.1 棚：
[0111] 表5多元回归后各地质特征参数系数表
[0112]
[0113]
[0114] 最后，根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的 可采储量。其中，致密油区域中待评价区域的地质因素参数表(见表6)，将这些地质因素参 数代入单井平均可采储量预测模型，预测单井平均可采储量，计算出待评价区的预测单井 平均可采储量(模拟EUR)。预测结果模拟EUR为260355bbl。根据致密油区域中待评价区域面 积和已钻井平均井控面积，可计算出区块可钻水平井数为10口，然后将可钻水平井数乘以 模拟EUR，可得致密油区域中待评价区域的致密油预测可采储量为260万桶。
[0115]表6预测区块地质因素参数表
[0116]
[0117]根据后期实际资料，实钻一口生产井，基于半年的生产数据，利用产量递减法计算 出的实际EUR为254056bbl，与预测值相对误差为2.48%，进一步说明了该方法的可行性。
[0118] 如图4所示，致密油区域中各个区块的单井平均可采储量预测值与致密油区域中 各个区块的实际单井平均可采储量的相关系数达到0.85以上，平均误差为7.1%，说明预测 结果较为可靠(见表7)。
[0119] 表7致密油区域中各区块实际EUR与模拟EUR对比表
[0120]
[0122]~与上述预测致密油待评价区域可采储量的方法对应，本申请实施例的预测致密油 待评价区域可采储量的装置包括：
[0123] 有井区域分块模块，用于将致密油区域中已有生产井区域划分为若干区块，其中 每个区块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年 限的生产井；
[0124] 地质特征参数获取模块，用于确定每个区块的地质特征参数；
[0125] 单井储量获取模块，用于依据每个区块内各生产井的历史生产数据对应确定该区 块的各生产井的单井可采储量，并据此确定每个区块的单井平均可采储量；
[0126] 地质特征参数优选模块，用于利用相关性分析法从每个区块的地质特征参数选出 与其单井平均可采储量正相关的地质特征影响参数；
[0127] 预测模型建立模块，用于根据各个区块的单井平均可采储量与其地质特征影响参 数建立多元线性模型；
[0128] 预测模型建立模块，用于将各个区块的地质特征影响参数代入所述多元线性模 型，并进行多元线性回归分析，以确定所述多元线性模型的模型系数，从而获得所述致密油 区域的单井平均可采储量预测模型；
[0129] 待评价区域预测模块，用于根据所述单井平均可采储量预测模型预测所述致密油 区域中待评价区域的可采储量。
[0130] 至于上述模块的具体情况，可参见上述预测致密油待评价区域可采储量的方法， 在此不再赘述。
[0131] 本申请实施例以影响致密油产量的关键地质因素与致密油单井EUR的关系，解决 了现有技术方案中均未综合考虑地质因素与单井EUR关系的缺点；选择以待评价区相邻地 区已钻井数据为基础，形成适合该区预测模型，避免了其它地区预测模型对待评价区不适 用的问题，并且利用相关性分析优选待评价区主控因素，避免了以往地质预测模型使用相 同地质参数来预测致密油单井EUR的不适用问题；而且，还利用多元回归的方法来确定EUR 预测模型中各个地质参数的系数，即解决了人为主观确定相似系数的误差大的缺点。此外， 在于利用待评价区面积确定可钻生产井数，进而获得可采储量，避免了以往预测出的可采 资源量无法直接应用于经济评价中。因此，本申请实施例的评价结果更为精确可靠，为致密 油区块经济评价提供了依据。
[0132] 本领域技术人员还可以了解到本申请实施例列出的各种说明性逻辑块、单元和步 骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定 的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种 方法实现所述的功能，但这种实现不应被理解为超出本申请实施例保护的范围。
[0133] 本申请实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器， 数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门 或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处 理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制 器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器， 多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置 来实现。
[0134] 本申请实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软 件模块、或者这两者的结合。软 件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储 器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介 中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并 可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可 以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用 户终端中的不同的部件中。
[0135] 在一个或多个示例性的设计中，本申请实施例所描述的上述功能可以在硬件、软 件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的 媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以 是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但 不限于RAM、ROM、EEPROM、⑶-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任 何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理 器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例 如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双 绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的 电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、错射盘、光盘、DVD、软盘 和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合 也可以包含在电脑可读媒介中。
[0136] 以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申 请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包 含在本申请的保护范围之内。
【主权项】
1. 一种预测致密油待评价区域可采储量的方法，其特征在于，包括以下步骤： 将致密油区域中已有生产井区域划分为若干区块； 确定每个区块的地质特征参数； 依据每个区块内各生产井的历史生产数据对应确定该区块的各生产井的单井可采储 量，并据此确定每个区块的单井平均可采储量； 利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储量正相关的 地质特征影响参数； 根据各个区块的单井平均可采储量与其地质特征影响参数建立多元线性模型； 将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归分析， 以确定所述多元线性模型的模型系数，从而获得所述致密油区域的单井平均可采储量预测 丰旲型； 根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采储量。2. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，每个区 块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年限的生 产井。3. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，在所述 利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储量正相关的地质 特征影响参数中： 当区块的地质特征参数对该区块的单井平均可采储量的影响呈抛物线型关系时，根据 公式：将该地质特征参数对该区块的单井平均可采储量的影响变换成正相 关关系，其中，X为变换后的地质特征参数，Xo为原地质特征参数，XmS地质特征参数的最优 经验值。4. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述多 元线性模型为：其中，EUR为致密油区域的单井平均可采储量，e为自然 数，A1为第i个地质特征影响参数，xi为模型系数。5. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，在所述 将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归分析之前， 还包括： 根据公式将所述将各个区块的地质特征影响参数进行去量纲的标准化处理， 其中仏为标准化后的第i个地质特征影响参数，Aiq为第i个地质特征影响参数，Aim为第i个地 质特征影响参数的最优经验值； 将所述多元线性模型进行两边取对数处理； 对应的，所述将各个区块的地质特征影响参数代入所述多元线性模型，具体为： 将各个区块的经去量纲的标准化处理后的地质特征影响参数代入经两边取对数处理 后的多元线性模型。6. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述确 定每个区块的地质特征参数，包括： 利用钻井资料确定各区块内致密油层系的产层厚度； 利用测井资料和测试分析资料确定各区块内致密油层系的平均孔隙度、产层泥地比以 及供烃源岩的有机质丰度、有机质成熟度；以及， 利用地层测试资料确定各区块内地层压力系数或压力梯度、气油比、油水比和原油重 度。7. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述根 据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采储量，具体包 括： 获取所述致密油区域中待评价区域的地质特征参数； 将所述致密油区域中待评价区域的地质特征参数代入所述单井平均可采储量预测模 型，获得所述致密油区域中待评价区域的单井平均可采储量； 获取所述致密油区域中待评价区域的面积，以及所述致密油区域中已有生产井区域内 生产井平均井控面积； 根据所述致密油区域中待评价区域的单井平均可采储量、所述致密油区域中待评价区 域的面积、以及所述致密油区域中已有生产井区域内生产井平均井控面积预测所述致密油 区域中待评价区域的可采储量。8. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，在所述 根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采储量之前， 还包括： 根据所述单井平均可采储量预测模型及所述致密油区域中各个区块的地质特征参数， 获取所述致密油区域中各个区块的单井平均可采储量预测值； 判断所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单井平 均可采储量的之间的误差是否均小于设定误差阈值； 如果所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单井平 均可采储量的之间的误差均小于设定误差阈值，则根据所述单井平均可采储量预测模型预 测所述致密油区域中待评价区域的可采储量。9. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述相 关性分析法包括Pearson相关性分析法。10. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述 设定生产历史年限包括半年。11. 一种预测致密油待评价区域可采储量的装置，其特征在于，包括： 有井区域分块模块，用于将致密油区域中已有生产井区域划分为若干区块，其中每个 区块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年限的 生产井； 地质特征参数获取模块，用于确定每个区块的地质特征参数； 单井储量获取模块，用于依据每个区块内各生产井的历史生产数据对应确定该区块的 各生产井的单井可采储量，并据此确定每个区块的单井平均可采储量； 地质特征参数优选模块，用于利用相关性分析法从每个区块的地质特征参数选出与其 单井平均可采储量正相关的地质特征影响参数； 预测模型建立模块，用于根据各个区块的单井平均可采储量与其地质特征影响参数建 立多元线性模型； 预测模型建立模块，用于将各个区块的地质特征影响参数代入所述多元线性模型，并 进行多元线性回归分析，以确定所述多元线性模型的模型系数，从而获得所述致密油区域 的单井平均可采储量预测模型； 待评价区域预测模块，用于根据所述单井平均可采储量预测模型预测所述致密油区域 中待评价区域的可采储量。
【专利摘要】本申请实施例提供了一种预测致密油待评价区域可采储量的方法及装置。本申请实施例利用待评价区邻区已投产生产井的生产情况，同时结合地质因素的变化对致密油单井可采储量进行合理预测，利用相关性分析优选待评价区主控因素，避免了以往地质预测模型使用相同地质参数来预测致密油单井EUR的不适用问题；而且，还利用多元回归的方法来确定EUR预测模型中各个地质参数的系数，即解决了人为主观确定相似系数的误差大的缺点。因此，本申请实施例的评价结果更为精确可靠，为后期致密油区块经济评价提供了依据。
【IPC分类】G06Q10/04, G06Q50/02
【公开号】CN105488583
【申请号】CN201510816016
【发明人】张新顺, 王红军, 马锋, 刘祚冬, 汪永华, 刘淑文
【申请人】中国石油天然气股份有限公司
【公开日】2016年4月13日
【申请日】2015年11月23日
【技术领域】
[0001]本申请涉及致密油开发技术领域，尤其是涉及一种预测致密油待评价区域可采储 量的方法及装置。
【背景技术】
[0002] 致密油资源潜力巨大，目前且尚处于勘探初期，是未来石油产量的重要增长点。目 前致密油可采储量评价体系仍以常规油气评价为主。目前涉及致密油储量评价的方法主要 技术有:容积法、概率统计法、类比法、井距外推法、产量递减法等。其中：
[0003] ①容积法，基于测井和地震资料的刻画出待评价区地下致密油的空间分布情况， 计算出致密油层的体积，通过油层的孔隙度估算出致密油的储量。
[0004] ②概率统计法，利用生产井最终估算可采储量(Estimated Ultimate Recovery， EUR)，建立生产井EUR概率分布图，确定待评价区的钻井数N以后，从生产井EUR概率分布图 上通过蒙特卡洛随机模拟对应得到单井EUR，选取N次，得到N个单井EUR值，求得平均值然后 输入。重复该过程5000到10000次，得到5000-10000个平均值，根据这些数值做出新的EUR概 率统计分布图，根据不同概率处获得未钻探可采储量。
[0005] ③类比法，利用在产致密油区的地质条件和开发井特征，建立刻度区，通过待评价 区地质条件，类比已有刻度区，确定待评价区开发井特征，进而推算出待开发区的储量。
[0006] ④外推井距法，利用已投产生产井的生产曲线，计算出单井的EUR，即为1P级储量。 对井周围按照一个井距进行外推，外推区域的单井EUR与已钻井EUR相当，即生产井一个井 距以内区域的预测产量为2P储量。继续向外推第二个井距的产量，仍认为与已投产井的产 量相当，即生产井2个井距以内区域的预测产量，为3P储量。
[0007] ⑤递减曲线法，利用具有半年以上生产历史的生产井，通过产量递减模型，计算出 来的单井最终估算可采储量。
[0008] 然而，容积法基于常规油气地下聚集模式来评价致密油储量，对于均质厚层致密 油层有一定适用性，但对于大多数致密油层而言，均具有强非均质性，横向纵向变化都很 大，其地质特征参数几乎无法取准，因而最终评价结果精度极低，该方法仅适用于评价致密 油资源量。
[0009] 概率统计法是基于数学统计的方法来计算致密油储量，其中未考虑地质因素的影 响。但是，地质因素的差异对于任何地区油气藏储量都会有非常大的影响。尤其是可以反映 生油能力的烃源条件和反映储集能力的储集条件。另外，概率统计法也未考虑工程因素的 影响，例如水平井段长度以及压裂后的泄流半径，这些因素也会直接造成单井产量大幅变 化。因而，概率统计法忽视了地质因素，无法获得较为精确的致密油储量。
[0010]类比法需要大量的与被评价区域类似的对比数据，即刻度区数据。在类比过程中， 需要通过给评价区各个地质要素进行打分，来给出类比相似系数，以确定被评价区与刻度 区的相似度。但这个确定类比相似系数的过程受到评价人主观影响较大，结果也因此差异 明显。另外，若刻度区与被评价区不属于同一致密岩层，其类比预测精度本身也会比较低。
[0011] 外推井距法，同样缺乏对地质因素差异性的考虑，其默认致密油生产井周边2个井 距内的地质条件几乎一致，未来新钻生产井的产能也会基本相同；但事实上，因为地质因素 总会发生变化的，在已投产生产井的一侧地质因素可能更好，另一侧可能更差。最重要的 是，在致密油井较少地区，外推井距法会极大的低估了致密油的实际可采储量。对于目前大 部分尚处于勘探开发初期的致密油开采区，该方法评价出的储量误差极大。因此，外推井距 法，并不能很好的评价致密油储量。
[0012] 递减曲线法，受控于生产井的生产历史数据且仅能对单井的可采储量进行预测。 但在致密油新区块的勘探过程中，是不具备开发井及生产历史资料的，所以该方法仅能应 用于对已开发的致密油区进行评价，无法对勘探阶段的致密油区块进行预测。
[0013] 综上所述，由于致密油资源与常规油气资源在地质条件、开发手段等方面具有显 著的差异，现有的技术方案均难以获得较为准确的致密油储量。
【发明内容】
[0014] 本申请实施例的目的在于提供一种预测致密油待评价区域可采储量的方法及装 置
[0015] 为达到上述目的，一方面，本申请实施例提供了一种预测致密油待评价区域可采 储量的方法，包括以下步骤：
[0016] 将致密油区域中已有生产井区域划分为若干区块；
[0017]确定每个区块的地质特征参数；
[0018] 依据每个区块内各生产井的历史生产数据对应确定该区块的各生产井的单井可 采储量，并据此确定每个区块的单井平均可采储量；
[0019] 利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储量正相 关的地质特征影响参数；
[0020] 根据各个区块的单井平均可采储量与其地质特征影响参数建立多元线性模型；
[0021] 将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归分 析，以确定所述多元线性模型的模型系数，从而获得所述致密油区域的单井平均可采储量 预测模型；
[0022] 根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采 储量。
[0023] 本申请实施例的预测致密油待评价区域可采储量的方法，每个区块内的各个生产 井处于同一构造带且每个区块具有至少两口达到设定生产历史年限的生产井。
[0024] 本申请实施例的预测致密油待评价区域可采储量的方法，在所述利用相关性分析 法从每个区块的地质特征参数选出与其单井平均可采储量正相关的地质特征影响参数中：
[0025] 当区块的地质特征参数对该区块的单井平均可采储量的影响呈抛物线型关系时， 根据公式
·将该地质特征参数对该区块的单井平均可采储量的影响变换成 正相关关系，其中，X为变换后的地质特征参数，X〇为原地质特征参数，Xm为地质特征参数的 最优经验值。
[0026] 本申请实施例的预测致密油待评价区域可采储量的方法，所述多元线性模型为：
[0027] +名2 +... + <，其中，EUR为致密油区域的单井平均可采储量，e为 自然数，仏为第i个地质特征影响参数，xi为模型系数。
[0028] 本申请实施例的预测致密油待评价区域可采储量的方法，在所述将各个区块的地 质特征影响参数代入所述多元线性模型，并进行多元线性回归分析之前，还包括：
[0029] 根据公式4=^1将所述将各个区块的地质特征影响参数进行去量纲的标准化处 Am 理，其中Ai为标准化后的第i个地质特征影响参数，Ai〇为第i个地质特征影响参数，Aim为第i 个地质特征影响参数的最优经验值；
[0030] 将所述多元线性模型进行两边取对数处理；
[0031 ]对应的，所述将各个区块的地质特征影响参数代入所述多元线性模型，具体为： [0032]将各个区块的经去量纲的标准化处理后的地质特征影响参数代入经两边取对数 处理后的多元线性模型。
[0033]本申请实施例的预测致密油待评价区域可采储量的方法，所述确定每个区块的地 质特征参数，包括：
[0034]利用钻井资料确定各区块内致密油层系的产层厚度；
[0035]利用测井资料和测试分析资料确定各区块内致密油层系的平均孔隙度、产层泥地 比以及供烃源岩的有机质丰度、有机质成熟度；以及，
[0036]利用地层测试资料确定各区块内地层压力系数或压力梯度、气油比、油水比和原 油重度。
[0037] 本申请实施例的预测致密油待评价区域可采储量的方法，所述根据所述单井平均 可采储量预测模型预测所述致密油区域中待评价区域的可采储量，具体包括：
[0038] 获取所述致密油区域中待评价区域的地质特征参数；
[0039] 将所述致密油区域中待评价区域的地质特征参数代入所述单井平均可采储量预 测模型，获得所述致密油区域中待评价区域的单井平均可采储量；
[0040] 获取所述致密油区域中待评价区域的面积，以及所述致密油区域中已有生产井区 域内生产井平均井控面积；
[0041] 根据所述致密油区域中待评价区域的单井平均可采储量、所述致密油区域中待评 价区域的面积、以及所述致密油区域中已有生产井区域内生产井平均井控面积预测所述致 密油区域中待评价区域的可采储量。
[0042] 本申请实施例的预测致密油待评价区域可采储量的方法，在所述根据所述单井平 均可采储量预测模型预测所述致密油区域中待评价区域的可采储量之前，还包括：
[0043] 根据所述单井平均可采储量预测模型及所述致密油区域中各个区块的地质特征 参数，获取所述致密油区域中各个区块的单井平均可采储量预测值；
[0044]判断所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差是否均小于设定误差阈值；
[0045]如果所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差均小于设定误差阈值，则根据所述单井平均可采储量预测模 型预测所述致密油区域中待评价区域的可采储量。
[0046] 本申请实施例的预测致密油待评价区域可采储量的方法，所述相关性分析法包括 Pearson相关性分析法。
[0047] 本申请实施例的预测致密油待评价区域可采储量的方法，所述设定生产历史年限 包括半年。
[0048] 另一方面，本申请实施例还提供了一种预测致密油待评价区域可采储量的装置， 包括：
[0049]有井区域分块模块，用于将致密油区域中已有生产井区域划分为若干区块，其中 每个区块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年 限的生产井；
[0050] 地质特征参数获取模块，用于确定每个区块的地质特征参数；
[0051] 单井储量获取模块，用于依据每个区块内各生产井的历史生产数据对应确定该区 块的各生产井的单井可采储量，并据此确定每个区块的单井平均可采储量；
[0052] 地质特征参数优选模块，用于利用相关性分析法从每个区块的地质特征参数选出 与其单井平均可采储量正相关的地质特征影响参数；
[0053] 预测模型建立模块，用于根据各个区块的单井平均可采储量与其地质特征影响参 数建立多元线性模型；
[0054] 预测模型建立模块，用于将各个区块的地质特征影响参数代入所述多元线性模 型，并进行多元线性回归分析，以确定所述多元线性模型的模型系数，从而获得所述致密油 区域的单井平均可采储量预测模型；
[0055] 待评价区域预测模块，用于根据所述单井平均可采储量预测模型预测所述致密油 区域中待评价区域的可采储量。
[0056] 本申请实施例利用待评价区邻区已投产生产井的生产情况，同时结合地质因素的 变化对致密油单井可采储量进行合理预测，利用相关性分析优选待评价区主控因素，避免 了以往地质预测模型使用相同地质参数来预测致密油单井EUR的不适用问题;而且，还利用 多元回归的方法来确定EUR预测模型中各个地质参数的系数，即解决了人为主观确定相似 系数的误差大的缺点。因此，本申请实施例的评价结果更为精确可靠，为后期致密油区块经 济评价提供了依据。
【附图说明】
[0057] 此处所说明的附图用来提供对本申请实施例的进一步理解，构成本申请实施例的 一部分，并不构成对本申请实施例的限定。在附图中：
[0058] 图1为本申请一实施例的预测致密油待评价区域可采储量的方法的流程图；
[0059] 图2为本申请一实施例的区块划分示意图；
[0060] 图3为本申请一实施例的递减曲线模型；
[0061] 图4为本申请一实施例的模拟EUR与实际EUR对比曲线示意图。
【具体实施方式】
[0062] 为使本申请实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附 图，对本申请实施例做进一步详细说明。在此，本申请实施例的示意性实施例及其说明用于 解释本申请实施例，但并不作为对本申请实施例的限定。
[0063] 下面结合附图，对本申请实施例的【具体实施方式】作进一步的详细说明。
[0064] 参考图1所示，本申请实施例的预测致密油待评价区域可采储量的方法包括以下 步骤：
[0065]步骤S1，将致密油区域中已有生产井区域划分为若干区块。为了有利于达到更为 准确的评价效果每个区块内的各个生产井处于同一构造带(斜坡带、背斜区、凹陷中心等） 且每个区块具有至少两口达到设定生产历史年限（例如半年）的生产井。此外，实验研究表 明，在满足上述要求的前提下，当每个区块的面积大小越相近，每个区块内的钻井分布越相 对聚集且单井产量越相近时，预测出的产量的偶然误差也就越小。
[0066] 步骤S2,确定每个区块的地质特征参数。比如利用钻井资料确定各区块内致密油 层系的产层厚度等；利用测井资料和测试分析资料确定各区块内致密油层系的平均孔隙 度、产层泥地比以及供烃源岩的有机质丰度、有机质成熟度等；以及，利用地层测试资料确 定各区块内地层压力系数或压力梯度、气油比、油水比和原油重度等。
[0067] 步骤S3,依据每个区块内各生产井的历史生产数据对应确定该区块的各生产井的 单井可采储量，并据此确定每个区块的单井平均可采储量。其中，若生产井为水平井，还需 利用公式= + X 对其压裂水平段长进行标准化，然后再进行小区块内的平均，其 中，L为某水平井压裂水平段长度；ζ为研究区标准压裂水平段长度;EUR°为该井原始EUR。 其中，生产井的历史生产数据可以包括产量地建模型。
[0068] 步骤S4,利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储 量正相关的地质特征影响参数。相关性分析法例如可以采用Pearson相关性分析法。
[0069] 其中，当区块的地质特征参数对该区块的单井平均可采储量的影响呈抛物线型关 系时，根据公式
J将该地质特征参数对该区块的单井平均可采储量的影响 变换成正相关关系，其中，X为变换后的地质特征参数，X〇为原地质特征参数，Xm为地质特征 参数的最优经验值。
[0070] 需要说明的是，本申请实施例中，地质特征影响参数实际为一种地质特征参数，由 于该地质特征参数会影响单井平均可采储量，为了区分这种地质特征参数与其他地质特征 参数，在这里将其称为地质特征影响参数。
[0071] 步骤S5,根据各个区块的单井平均可采储量与其地质特征影响参数建立多元线性 模型。其中，多元线性模型还可以参考，多元线性模型具体可为：
[0072] fY +為 12 +…+ ^^其中^皿为致密油区域的单井平均可采储量…为 自然数，仏为第i个地质特征影响参数，xi为模型系数。
[0073] 步骤S6,将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线 性回归分析，以确定所述多元线性模型的模型系数，即确定系数xO~xi的最优解，从而获得 所述致密油区域的单井平均可采储量预测模型。
[0074] 步骤S7,根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域 的可采储量。具体的：
[0075] 获取所述致密油区域中待评价区域的地质特征参数；
[0076] 将所述致密油区域中待评价区域的地质特征参数代入所述单井平均可采储量预 测模型，获得所述致密油区域中待评价区域的单井平均可采储量R;
[0077] 获取所述致密油区域中待评价区域的面积S1，以及所述致密油区域中已有生产井 区域内生产井平均井控面积S2;
[0078] 根据所述致密油区域中待评价区域的单井平均可采储量、所述致密油区域中待评 价区域的面积、以及所述致密油区域中已有生产井区域内生产井平均井控面积预测所述致 密油区域中待评价区域的可采储量，即根据公式可计算出致密油区域中待评价区域 的可采储量。
[0079] 本申请另一实施例中，在所述将各个区块的地质特征影响参数代入所述多元线性 模型，并进行多元线性回归分析之前，还可以包括：
[0080] 根据公式4 将所述将各个区块的地质特征影响参数进行去量纲的标准化处 im 理，其中Ai为标准化后的第i个地质特征影响参数，Ai〇为第i个地质特征影响参数，Aim为第i 个地质特征影响参数的最优经验值；
[0081 ]将所述多元线性模型进行两边取对数处理；
[0082] 对应的，所述将各个区块的地质特征影响参数代入所述多元线性模型，具体为：
[0083] 将各个区块的经去量纲的标准化处理后的地质特征影响参数代入经两边取对数 处理后的多元线性模型。
[0084] 本申请另一实施例中，在所述根据所述单井平均可采储量预测模型预测所述致密 油区域中待评价区域的可采储量之前，还包括：
[0085] 根据所述单井平均可采储量预测模型及所述致密油区域中各个区块的地质特征 参数，获取所述致密油区域中各个区块的单井平均可采储量预测值；
[0086] 判断所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差是否均小于设定误差阈值；
[0087]如果所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单 井平均可采储量的之间的误差均小于设定误差阈值，则根据所述单井平均可采储量预测模 型预测所述致密油区域中待评价区域的可采储量。 从而可以确保只有准确性满足要求的预 测模型才可以用于预测所述致密油区域中待评价区域的可采储量。
[0088] 虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解， 这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行（例如使用并行 处理器或多线程环境），比如步骤S2和步骤S3之间的顺序就调换，也可以同时进行。
[0089] 下面以北美某区块致密油地质情况为例，介绍本申请实施例的具体应用。
[0090] 首先，针对北美某区块致密油地质情况，确定出致密油主要为源内成藏，未经历明 显侧向运移;根据研究区井位分布情况，结合地质条件、地理特征等，将区块划分为面积相 近的N个小区块(如图2所示）。
[0091] 其次，利用每个小区块钻井、测井、测试分析等资料，确定每个小区块的地质因素 参数，厚度、孔隙度、渗透率、有机质丰度、有机质成熟度、产层泥地比、地层压力系数等。具 体的，利用钻井资料确定各小区块内致密油层系的产层厚度;利用测井资料和测试分析资 料可确定各小区块内致密油层系平均孔隙度、产层泥地比，和供烃源岩的有机质丰度、有机 质成熟度;利用地层测试资料确定地层压力系数(或压力梯度）、气油比、油水比、原油重度 (API)。
[0092]然后，依据每个区块内各生产井的历史生产数据(例如图3所示的产量递减模型） 对应确定该区块的各生产井的单井可采储量，并算出每个小区块的单井平均可采储量。 [0093]其次，利用Pearson相关性分析法确定影响EUR的地质因素，当概率α〈〇. 〇1时，即认 为两者明显相关(表1)。
[0094] 表1研究区地质因素与EUR相关性分析参数表
[0095]
[0096] 当部分地质因素对EUR的影响可能是呈抛物线状型关系，因此要根据研究区公认 的最优条件值选取参数最优值。例如统计发现，研究区产层段泥地比Μ为0.55时为水平井产 量较高，泥地比Μ过大或过小，产量都较低。因此认为泥地比Μ是与EUR成抛物线相关，需要利 用进行变换，变成正相关关系。其中，Μ为变换后的泥地比，M0是原来的泥地 比，Mm为?.
55，埘％区统计规律得到。
[0097]由此，选取影响EUR的主要相关地质因素:产层总厚度Η、泥地比Μ、压力系数Ρ、成熟 度Ro、下T0C、含油饱和度So、孔隙度Φ (见表2)。
[0098] 表2各区块原始地质特征参数表
[0099]
[0100] 再次，已选取地质特征参数与EUR均具有正相关关系，而多个因素均对EUR值有影 响，因此建立单井可采储量EUR与地质因素参数关系模型：
[0101] EURi = X X Mf3 X pp X H X TOCp X Χφ?其中，Η为产层总厚度， Μ为泥地比，Ρ压力系数，Ro成熟度，T0C为下部T0C值，So为含油饱和度So，Φ为孔隙度，xl~ x8分别为各个系数。
[0102] 在将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归 分析之前，将上述公式两边对数处理，可得
[0103]
[0104] 其中，可以先利用公式ln(A:) = h(^x100 )对模型中所有参数进行去量纲的标准 Am 化处理(见表3)，再做上述对数处理(各区块对数变换后地质特征参数见表4)。
[0105] 表3各区块标准化后地质特征参数表
[0106]
[0107] 表4各区块对数变换后地质特征参数表
[0108]
[0109] 然后，将处理过的各个小区块的地质因素参数代入模型，进行多元线性回归分析， 确定系数XI~X8的最优解(见表5)，从而获得所述致密油区域的单井平均可采储量预测模 型：
[0110] EUR = e3.料3 X: H0.1652 X M:0'0615 X 把.1495 X: TOC0'2挪 X 2848 X φ L.1 棚：
[0111] 表5多元回归后各地质特征参数系数表
[0112]
[0113]
[0114] 最后，根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的 可采储量。其中，致密油区域中待评价区域的地质因素参数表(见表6)，将这些地质因素参 数代入单井平均可采储量预测模型，预测单井平均可采储量，计算出待评价区的预测单井 平均可采储量(模拟EUR)。预测结果模拟EUR为260355bbl。根据致密油区域中待评价区域面 积和已钻井平均井控面积，可计算出区块可钻水平井数为10口，然后将可钻水平井数乘以 模拟EUR，可得致密油区域中待评价区域的致密油预测可采储量为260万桶。
[0115]表6预测区块地质因素参数表
[0116]
[0117]根据后期实际资料，实钻一口生产井，基于半年的生产数据，利用产量递减法计算 出的实际EUR为254056bbl，与预测值相对误差为2.48%，进一步说明了该方法的可行性。
[0118] 如图4所示，致密油区域中各个区块的单井平均可采储量预测值与致密油区域中 各个区块的实际单井平均可采储量的相关系数达到0.85以上，平均误差为7.1%，说明预测 结果较为可靠(见表7)。
[0119] 表7致密油区域中各区块实际EUR与模拟EUR对比表
[0120]
[0122]~与上述预测致密油待评价区域可采储量的方法对应，本申请实施例的预测致密油 待评价区域可采储量的装置包括：
[0123] 有井区域分块模块，用于将致密油区域中已有生产井区域划分为若干区块，其中 每个区块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年 限的生产井；
[0124] 地质特征参数获取模块，用于确定每个区块的地质特征参数；
[0125] 单井储量获取模块，用于依据每个区块内各生产井的历史生产数据对应确定该区 块的各生产井的单井可采储量，并据此确定每个区块的单井平均可采储量；
[0126] 地质特征参数优选模块，用于利用相关性分析法从每个区块的地质特征参数选出 与其单井平均可采储量正相关的地质特征影响参数；
[0127] 预测模型建立模块，用于根据各个区块的单井平均可采储量与其地质特征影响参 数建立多元线性模型；
[0128] 预测模型建立模块，用于将各个区块的地质特征影响参数代入所述多元线性模 型，并进行多元线性回归分析，以确定所述多元线性模型的模型系数，从而获得所述致密油 区域的单井平均可采储量预测模型；
[0129] 待评价区域预测模块，用于根据所述单井平均可采储量预测模型预测所述致密油 区域中待评价区域的可采储量。
[0130] 至于上述模块的具体情况，可参见上述预测致密油待评价区域可采储量的方法， 在此不再赘述。
[0131] 本申请实施例以影响致密油产量的关键地质因素与致密油单井EUR的关系，解决 了现有技术方案中均未综合考虑地质因素与单井EUR关系的缺点；选择以待评价区相邻地 区已钻井数据为基础，形成适合该区预测模型，避免了其它地区预测模型对待评价区不适 用的问题，并且利用相关性分析优选待评价区主控因素，避免了以往地质预测模型使用相 同地质参数来预测致密油单井EUR的不适用问题；而且，还利用多元回归的方法来确定EUR 预测模型中各个地质参数的系数，即解决了人为主观确定相似系数的误差大的缺点。此外， 在于利用待评价区面积确定可钻生产井数，进而获得可采储量，避免了以往预测出的可采 资源量无法直接应用于经济评价中。因此，本申请实施例的评价结果更为精确可靠，为致密 油区块经济评价提供了依据。
[0132] 本领域技术人员还可以了解到本申请实施例列出的各种说明性逻辑块、单元和步 骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定 的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种 方法实现所述的功能，但这种实现不应被理解为超出本申请实施例保护的范围。
[0133] 本申请实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器， 数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门 或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处 理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制 器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器， 多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置 来实现。
[0134] 本申请实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软 件模块、或者这两者的结合。软 件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储 器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介 中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并 可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可 以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用 户终端中的不同的部件中。
[0135] 在一个或多个示例性的设计中，本申请实施例所描述的上述功能可以在硬件、软 件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的 媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以 是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但 不限于RAM、ROM、EEPROM、⑶-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任 何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理 器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例 如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双 绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的 电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、错射盘、光盘、DVD、软盘 和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合 也可以包含在电脑可读媒介中。
[0136] 以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申 请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包 含在本申请的保护范围之内。
【主权项】
1. 一种预测致密油待评价区域可采储量的方法，其特征在于，包括以下步骤： 将致密油区域中已有生产井区域划分为若干区块； 确定每个区块的地质特征参数； 依据每个区块内各生产井的历史生产数据对应确定该区块的各生产井的单井可采储 量，并据此确定每个区块的单井平均可采储量； 利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储量正相关的 地质特征影响参数； 根据各个区块的单井平均可采储量与其地质特征影响参数建立多元线性模型； 将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归分析， 以确定所述多元线性模型的模型系数，从而获得所述致密油区域的单井平均可采储量预测 丰旲型； 根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采储量。2. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，每个区 块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年限的生 产井。3. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，在所述 利用相关性分析法从每个区块的地质特征参数选出与其单井平均可采储量正相关的地质 特征影响参数中： 当区块的地质特征参数对该区块的单井平均可采储量的影响呈抛物线型关系时，根据 公式：将该地质特征参数对该区块的单井平均可采储量的影响变换成正相 关关系，其中，X为变换后的地质特征参数，Xo为原地质特征参数，XmS地质特征参数的最优 经验值。4. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述多 元线性模型为：其中，EUR为致密油区域的单井平均可采储量，e为自然 数，A1为第i个地质特征影响参数，xi为模型系数。5. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，在所述 将各个区块的地质特征影响参数代入所述多元线性模型，并进行多元线性回归分析之前， 还包括： 根据公式将所述将各个区块的地质特征影响参数进行去量纲的标准化处理， 其中仏为标准化后的第i个地质特征影响参数，Aiq为第i个地质特征影响参数，Aim为第i个地 质特征影响参数的最优经验值； 将所述多元线性模型进行两边取对数处理； 对应的，所述将各个区块的地质特征影响参数代入所述多元线性模型，具体为： 将各个区块的经去量纲的标准化处理后的地质特征影响参数代入经两边取对数处理 后的多元线性模型。6. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述确 定每个区块的地质特征参数，包括： 利用钻井资料确定各区块内致密油层系的产层厚度； 利用测井资料和测试分析资料确定各区块内致密油层系的平均孔隙度、产层泥地比以 及供烃源岩的有机质丰度、有机质成熟度；以及， 利用地层测试资料确定各区块内地层压力系数或压力梯度、气油比、油水比和原油重 度。7. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述根 据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采储量，具体包 括： 获取所述致密油区域中待评价区域的地质特征参数； 将所述致密油区域中待评价区域的地质特征参数代入所述单井平均可采储量预测模 型，获得所述致密油区域中待评价区域的单井平均可采储量； 获取所述致密油区域中待评价区域的面积，以及所述致密油区域中已有生产井区域内 生产井平均井控面积； 根据所述致密油区域中待评价区域的单井平均可采储量、所述致密油区域中待评价区 域的面积、以及所述致密油区域中已有生产井区域内生产井平均井控面积预测所述致密油 区域中待评价区域的可采储量。8. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，在所述 根据所述单井平均可采储量预测模型预测所述致密油区域中待评价区域的可采储量之前， 还包括： 根据所述单井平均可采储量预测模型及所述致密油区域中各个区块的地质特征参数， 获取所述致密油区域中各个区块的单井平均可采储量预测值； 判断所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单井平 均可采储量的之间的误差是否均小于设定误差阈值； 如果所述致密油区域中每个区块的单井平均可采储量预测值与其每个区块的单井平 均可采储量的之间的误差均小于设定误差阈值，则根据所述单井平均可采储量预测模型预 测所述致密油区域中待评价区域的可采储量。9. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述相 关性分析法包括Pearson相关性分析法。10. 根据权利要求1所述的预测致密油待评价区域可采储量的方法，其特征在于，所述 设定生产历史年限包括半年。11. 一种预测致密油待评价区域可采储量的装置，其特征在于，包括： 有井区域分块模块，用于将致密油区域中已有生产井区域划分为若干区块，其中每个 区块内的各个生产井处于同一构造带且每个区块具有至少两口达到设定生产历史年限的 生产井； 地质特征参数获取模块，用于确定每个区块的地质特征参数； 单井储量获取模块，用于依据每个区块内各生产井的历史生产数据对应确定该区块的 各生产井的单井可采储量，并据此确定每个区块的单井平均可采储量； 地质特征参数优选模块，用于利用相关性分析法从每个区块的地质特征参数选出与其 单井平均可采储量正相关的地质特征影响参数； 预测模型建立模块，用于根据各个区块的单井平均可采储量与其地质特征影响参数建 立多元线性模型； 预测模型建立模块，用于将各个区块的地质特征影响参数代入所述多元线性模型，并 进行多元线性回归分析，以确定所述多元线性模型的模型系数，从而获得所述致密油区域 的单井平均可采储量预测模型； 待评价区域预测模块，用于根据所述单井平均可采储量预测模型预测所述致密油区域 中待评价区域的可采储量。
【专利摘要】本申请实施例提供了一种预测致密油待评价区域可采储量的方法及装置。本申请实施例利用待评价区邻区已投产生产井的生产情况，同时结合地质因素的变化对致密油单井可采储量进行合理预测，利用相关性分析优选待评价区主控因素，避免了以往地质预测模型使用相同地质参数来预测致密油单井EUR的不适用问题；而且，还利用多元回归的方法来确定EUR预测模型中各个地质参数的系数，即解决了人为主观确定相似系数的误差大的缺点。因此，本申请实施例的评价结果更为精确可靠，为后期致密油区块经济评价提供了依据。
【IPC分类】G06Q10/04, G06Q50/02
【公开号】CN105488583
【申请号】CN201510816016
【发明人】张新顺, 王红军, 马锋, 刘祚冬, 汪永华, 刘淑文
【申请人】中国石油天然气股份有限公司
【公开日】2016年4月13日
【申请日】2015年11月23日
最新回复(0)