CN105283867A - 用于基于平均有效裂缝长度的分布优化现有井和设计新井的系统和方法 - Google Patents

Abstract

本文公开了用于基于对于相对于不同储层性质的相应每个压裂阶段而言的每个平均有效裂缝长度的分布来优化现有井和设计新井的系统和方法。

Description

用于基于平均有效裂缝长度的分布优化现有井和设计新井的系统和方法

相关申请的交叉引用

不适用。

关于联邦资助研究的声明

不适用。

发明领域

本发明整体涉及用于基于平均有效裂缝长度的分布优化现有井和设计新井的系统和方法。更具体地讲，本发明涉及基于对于相对于不同储层性质的相应的压裂阶段而言的每个平均有效裂缝长度的分布来优化现有井和设计新井。

发明背景

井生产剖面的历史拟合代表油田开发规划中的重要组成部分。在非常规的场合中，必要的是以储层井为单位使用单井模拟器来历史拟合流体的产量衰减曲线并且准确地估计与由裂缝系统诱发的改造储层(引流)体积(SRV)相关联的改善的渗透率(k_imp)。重要的是应注意到，改善的渗透率不仅与渗透基质相关联，而且还对应于裂缝系统的增强的流体流动性质。

传统上，当微震监测数据无法获得时，历史拟合过程假设诱发的裂缝系统的过分简化的模型，该模型由具有相同(x_eff)、相同(SRV)的若干阶段的双翼裂缝构成，并且每个阶段仅一个双翼裂缝，如由图5A中的简单的示意性模型所示。在图5A中，双翼裂缝为大体上垂直于井轴线延伸的细长裂缝。在裂缝平面中，每个双翼裂缝实际上在两个方向上延伸相同的长度。双翼裂缝通常用两个主要参数建模：裂缝长度，也称为有效裂缝长度；和裂缝宽度。这两个参数通常与由裂缝接触的渗透基质的改善的渗透率有关。改善的渗透率和有效裂缝长度之间的关系因此通过式1来描述：

现在参见图1，示出了常规方法100的流程图，该方法用于使用单井储层模拟器来历史拟合生产剖面。

在步骤102中，将标准储层性质(例如，地层厚度、BHP、基质孔隙率和渗透率、岩石类型、标准裂缝设计性质(例如，简单双翼裂缝的有效裂缝长度和裂缝宽度)和生产数据分布(例如，天然气/石油/水速率和BHP)输入单井储层模拟器中。

在步骤104中，由单井储层模拟器使用本领域熟知的用于历史拟合的技术和从步骤102输入的数据来执行历史拟合。

在步骤106中，显示作为在步骤104中执行的历史拟合的结果的SRV的改善的渗透率(k_imp)。用于历史拟合的该常规方法100确定改善的渗透率的标准估计，但常常给出井生产性能的次优预测。因此，挑战在于更准确地估计有效裂缝长度，该长度表示比包括双翼裂缝的若干阶段的模型更现实的裂缝系统，这些阶段具有相同(x_eff)、相同(SRV)，并且每个阶段仅一个双翼裂缝。

附图简述

下文参照附图来描述本发明，其中类似的元件用类似的附图标记参考，并且在附图中：

图1是流程图，示出了用于使用单井储层模拟器来历史拟合生产剖面的常规方法。

图2是流程图，示出了用于实施本发明的方法的一个实施方案。

图3是流程图，示出了用于执行图2中的步骤204的方法的一个实施方案。

图4A是示出与裂缝簇(fracturecluster)相关联的一系列微震成像事件的显示画面。

图4B是示出基于图4A中的微震成像事件的时间相关性的3D裂缝平面的显示画面。

图5A是诱发的裂缝系统的简单示意性模型，示出了具有相同(x_eff)、相同(SRV)和每个阶段仅一个裂缝的双翼裂缝。

图5B是诱发的裂缝系统的复杂示意性模型，示出了多个复杂的裂缝网络，每个网络具有不同(x_eff)、不同(SRV)和每个阶段多个裂缝。

图6是框图，示出了用于实施本发明的计算机系统的一个实施方案。

优选实施方案的详细描述

因此，通过提供用于基于对于相对于不同储层性质的相应的压裂阶段而言的每个平均有效裂缝长度的分布来优化现有井和设计新井的系统和方法，本发明克服了现有技术中的一个或多个缺点。

在一个实施方案中，本发明包括一种用于优化按改造储层体积计的井产量的方法，该方法包括：i)输入一个或多个复杂储层性质和一个或多个复杂裂缝网络性质，复杂裂缝网络性质包括对应于复杂裂缝网络模型中的簇的数据；ii)基于复杂储层性质和复杂裂缝网络性质确定平均有效裂缝长度的分布；iii)使用计算机处理器从平均有效裂缝长度的分布对平均有效裂缝长度抽样；以及iv)使用平均有效裂缝长度的分布和抽样的平均有效裂缝长度通过历史拟合来优化井产量，以改善改造储层体积的渗透率。

在另一个实施方案中，本发明包括一种非暂时性程序载体装置，其有形地承载计算机可执行的指令以用于优化按改造储层体积计的井产量，该装置包括：i)输入一个或多个复杂储层性质和一个或多个复杂裂缝网络性质，复杂裂缝网络性质包括对应于复杂裂缝网络模型中的簇的数据；ii)基于复杂储层性质和复杂裂缝网络性质确定平均有效裂缝长度的分布；iii)从平均有效裂缝长度的分布对平均有效裂缝长度抽样；以及iv)使用平均有效裂缝长度的分布和抽样的平均有效裂缝长度通过历史拟合来优化井产量，以改善改造储层体积的渗透率。

在又一个实施方案中，本发明包括非暂时性程序载体装置，其有形地承载计算机可执行的指令以用于优化按改造储层体积计的井产量，该装置包括：i)输入一个或多个复杂储层性质和一个或多个复杂裂缝网络性质，复杂裂缝网络性质包括对应于复杂裂缝网络模型中的簇的数据；ii)通过以下方式确定平均有效裂缝长度的分布：a)读取针对每个井的每个压裂阶段中的每个裂缝平面的有效裂缝长度；b)使用用于相应的压裂阶段的每个有效裂缝长度来计算每个压裂阶段的平均有效裂缝长度；和c)通过将每个相应的压裂阶段的平均有效裂缝长度与每个储层或井测井性质关联来构建平均有效裂缝长度的分布；iii)从平均有效裂缝长度的分布对平均有效裂缝长度抽样；以及iv)使用平均有效裂缝长度的分布和抽样的平均有效裂缝长度通过历史拟合来优化井产量，以改善改造储层体积的渗透率。

虽然已详细地描述了本发明的主题，但描述本身并非意图限制本发明的范围。因此，该主题也可能结合其它技术以其它方式具体化，以包括类似于本文所述那些的不同步骤或步骤的组合。此外，虽然术语“步骤”在本文中可用来描述所采用的方法的不同要素，但该术语不应解释为暗示在本文所公开的各个步骤之间的任何特定的次序，除非由说明书明确地限制到特定的次序。虽然以下描述涉及石油和天然气行业，但本发明的系统和方法不限于此，并且也可以应用于其它行业以达到类似的结果。

方法描述

现在参见图2，示出了用于实施本发明的方法200的一个实施方案的流程图。方法200使用单井储层模拟器优化历史拟合生产剖面。

在步骤202中，使用参照图6进一步描述的客户端接口和/或视频接口将标准储层性质(例如，地层厚度、BHP、基质孔隙率和渗透率、岩石类型)、来自使用空间上分布在储层上且用井数据约束的映射性质(例如，TOC、孔隙率和脆性)的先进岩石物理测井解释的复杂储层性质(例如，岩石物理性质(例如，HC含量、粘土含量))、复杂裂缝网络(“CFN”)性质(例如，对应于CFN模型中的簇的数据)以及生产数据分布(例如，天然气/石油/水速率和BHP)输入单井储层模拟器中。簇提供裂缝系统的准确得多的表示，因为压裂不仅产生细长的双翼裂缝，而且产生较小的复杂裂缝的网络，这些复杂裂缝优选地全部互连并且在彼此之间连通，形成CFN。每个CFN受其它岩石性质诸如上文提及的标准储层性质和标测性质影响。

在步骤204中，确定平均有效裂缝长度的分布。参照图3进一步描述用于执行此步骤的方法的一个实施方案。

在步骤205中，从在步骤204中确定的平均有效裂缝长度的分布(离散的或连续的)对平均有效裂缝长度抽样。任何熟知的标准概率抽样技术(例如，随机抽样器)均可以用于抽样。这样，可以生成具有较低中值和较高概率场景(例如，P10、P50和P90模型)的估计的改善的渗透率(k_imp)的不确定性图。

在步骤206中，由SRV使用从步骤202输入的标准储层性质和生产数据分布、来自步骤204的平均有效裂缝长度的分布、来自步骤205的抽样的平均有效裂缝长度和本领域熟知的用于历史拟合的技术来执行历史拟合。

在步骤208中，使用参照图6进一步描述的视频接口显示作为在步骤206中执行的历史拟合结果的优化的改善的渗透率(k_imp)。方法200将提供对井生产动态的更准确的预测，所述井生产动态可以用来优化现有井和设计新井，因为它基于并且包括复杂储层性质和CFN性质，这优化了平均有效裂缝长度的分布。换句话讲，CFN不再与有效裂缝长度/宽度和改善的渗透率相关，而是与SRV相关。因此，目标是生成使SRV最大化的CFN并开发更准确地表示实际SRV的模型。

现在参见图3，示出了用于执行图2中的步骤204的方法300的一个实施方案的流程图。

在步骤301中，从在步骤202中输入的井的总数(W)中自动地选择井(w)，或者备选地可以使用参照图6进一步描述的客户端接口和/或视频接口来选择。

在步骤302中，从在步骤202中输入的每个井(w)的压裂阶段的总数(S)中自动地选择压裂阶段(s)，或者备选地可以使用参照图6进一步描述的客户端接口和/或视频接口来选择。

在步骤303中，从在步骤202中输入的每个压裂阶段(s)的裂缝平面的总数(F)中自动地选择裂缝平面(f)，或者备选地可以使用参照图6进一步描述的客户端接口和/或视频接口来选择。假设每个压裂阶段(s)内的裂缝平面(f)分布为簇并且不是简化的单个双翼裂缝。

在步骤304中，从对应于在步骤202中输入的CFN模型的数据读取用于所选裂缝平面(f)、压裂阶段(s)和井(w)的有效裂缝长度对应于CFN模型的数据可包括例如每个压裂阶段的簇的3D裂缝平面的数目。3D裂缝平面基于微震成像事件的时间分析来构造。在图4A中，示出了与裂缝簇相关联的一系列解释的微震成像事件的显示画面400a。在图4B中，示出了基于图4A中的微震成像事件的时间相关性的3D裂缝平面的显示画面400b。显示画面400b中的3D裂缝平面被井轨迹穿过，以示出压裂过程的解释结果。基于从步骤202输入的此数据，压裂阶段(s)和井(w)的所选裂缝平面(f)的最长轴线的尺寸或长度可被读取并指定为该所选裂缝平面(f)的有效裂缝长度。在图5B中，诱发的裂缝系统的复杂示意性模型示出了多个复杂的裂缝网络，每个网络具有不同(x_eff)、不同(SRV)和每个压裂阶段的多个裂缝。相比基于图5A中所示双翼裂缝的诱发的裂缝系统的简化模型，就裂缝系统的准确得多的表示而言，图5B中的更复杂模型的优点是显而易见的。

在步骤305中，使用在步骤304中读取的每个有效裂缝长度和公式2来计算压裂阶段(s)的平均有效裂缝长度

${\hat{x}}_{eff,s}^w=\frac{1}{F}\underset{f=1}{\overset{F}{\mathrm{\Σ}}}{x}_{eff,s,f}^w---\left(2\right)$

其中对应于所选压裂阶段(s)内的所选裂缝平面(f)的有效裂缝长度。

在步骤306中，方法300确定是否存在要从裂缝平面的总数(F)中选择的另一个裂缝平面(f)。如果存在要选择的另一个裂缝平面(f)，那么方法300返回到步骤303以从裂缝平面的总数(F)中选择另一个裂缝平面(f)。如果不存在要选择的另一个裂缝平面(f)，那么方法300前进至步骤307。

在步骤307中，方法300确定是否存在要从压裂阶段的总数(S)中选择的另一个压裂阶段(s)。如果存在要选择的另一个压裂阶段(s)，那么方法300返回到步骤302以从压裂阶段的总数(S)中选择另一个压裂阶段(s)。如果不存在要选择的另一个压裂阶段(s)，那么方法300前进至步骤308。

在步骤308中，方法300确定是否存在要从井的总数(W)中选择的另一个井(w)。如果存在要选择的另一个井(w)，那么方法300返回到步骤301以从井的总数(W)中选择另一个井(w)。如果不存在要选择的另一个井(w)，那么方法300前进至步骤309。

在步骤309中，从在步骤202中输入的复杂储层性质的总数(P)中自动地选择储层或井测井性质(p)，或者备选地可以使用参照图6进一步描述的客户端接口和/或视频接口来选择。

在步骤310中，将在步骤305中计算的每个相应的压裂阶段(s)的平均有效裂缝长度与在步骤309中选择的储层或井测井性质(p)关联，以构建平均有效裂缝长度的分布(离散的或连续的)离散条件分布(柱状图)可使用公式3构建：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{\mathrm{Pr}ob(P=p\∩X_{eff}={\hat{x}}_{eff,s}^w)}{\mathrm{Pr}ob(P=p)}---\left(3\right)$

其中“Prob”表示“概率”，(x_eff)限定作为概率因变量的平均有效裂缝长度的总抽样域，并且(P)限定作为概率自变量的复杂储层性质的总抽样域。

备选地，可以使用公式4构建连续条件分布(pdf)：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{{\mathrm{Prob}}_{P,X_{eff}}(p,{\hat{x}}_{eff,s}^w)}{{\mathrm{Prob}}_P\left(p\right)}---\left(4\right)$

其中限定(P)和(x_eff)的联合密度(pdf)，而(Prob_P(p))限定(P)的边际密度。出于pdf归一化目的，有必要保持Prob_P(p)>0。

在步骤312中，方法300确定是否存在要从复杂储层性质的总数(P)中选择的另一个储层或井测井性质(p)。如果存在要选择的另一个储层或井测井性质(p)，那么方法300返回到步骤309以从复杂储层性质的总数(P)中选择另一个储层或井测井性质(p)。如果不存在要选择的另一个储层或井测井性质(p)，那么方法300将平均有效裂缝长度的分布返回至步骤204。

系统描述

本发明可以通过诸如程序模块的计算机可执行指令程序来实施，该程序通常称为由计算机执行的软件应用程序或应用程序。软件可包括例如执行特定任务或实施特定抽象数据类型的例程、程序、对象、组件和数据结构。软件形成接口以允许计算机根据输入的源做出反应。Desktop(EarthModeling)可被用作用来实施本发明的接口应用程序，其为由LandmarkGraphicsCorporation销售的商用软件应用程序。该软件也可与其它代码段配合，以响应于结合所接收数据的源而接收的数据来引发多种任务。该软件可以被存储和/或承载在任意多种存储器上，诸如CD-ROM、磁盘、磁泡存储器和半导体存储器(例如，各种类型的RAM或ROM)。此外，该软件及其结果可以通过诸如光纤、金属线的多种载体媒介和/或通过诸如因特网的多种网络中的任一者来传输。

此外，本领域的技术人员将了解，本发明可以用多种计算机系统配置实践，包括手持装置、多处理器系统、基于微处理器的或可编程的消费电子产品、微型计算机、大型计算机等。许多计算机系统和计算机网络对于用于本发明而言是可以接受的。本发明可以在分布式计算环境中实施，在该环境中，任务由通过通信网络连接的远程处理装置来执行。在分布式计算环境中，程序模块可位于包括存储器存储装置的本地计算机存储介质和远程计算机存储介质两者中。因此，本发明可以结合各种硬件、软件或它们的组合在计算机系统或其它处理系统中实施。

现在参见图6，框图示出了用于在计算机上实施本发明的系统的一个实施方案。该系统包括有时称为计算系统的计算单元，其包含存储器、应用程序、客户端接口、视频接口和处理单元。计算单元仅是合适的计算环境的一个实例，而并非意图暗示对本发明的使用范围或功能的任何限制。

存储器主要存储应用程序，应用程序也可描述为包含计算机可执行指令的程序模块，计算机可执行指令由计算单元执行以用于实施本发明，如本文所述和图2-3中所示。因此，存储器包括井优化模块，其能够实现结合图2中的步骤204-205描述的方法。存储器也包括单井模拟器，其能够执行图2中的步骤206。QuickLook_TM和Knoesis/Slate_sm是可以使用的由HalliburtonCompany销售的单井模拟器的实例。上述模块和应用程序可以整合来自图6中所示其余应用程序的功能。特别地，Desktop(EarthModeling)可用作接口应用程序以执行图2中的步骤202和208。ASCII文件也包括在存储器中以用于存储从图2中的步骤202输入的数据。虽然Desktop(EarthModeling)和单井模拟器可以用作接口应用程序，但也可以使用其它接口应用程序，或者井优化模块可以用作独立式应用程序。

虽然计算单元示为具有一般化存储器，但计算单元通常包括多种计算机可读介质。作为实例，而非限制，计算机可读介质可包括计算机存储介质和通信介质。计算系统存储器可包括诸如只读存储器(ROM)和随机存取存储器(RAM)的易失性和/或非易失性存储器形式的计算机存储介质。包含在诸如启动期间帮助在计算单元内的元件之间传输信息的基本例程的基本输入/输出系统(BIOS)通常存储在ROM中。RAM通常包含可由处理单元即时存取和/或目前正处理的数据和/或程序模块。作为实例，而非限制，计算单元包括操作系统、应用程序、其它程序模块和程序数据。

存储器中所示的部件也可包括在其它可移动/不可移动、易失性/非易失性计算机存储介质中，或者它们可以通过应用程序接口(“API”)或云计算实施在计算单元中，应用程序接口(“API”)或云计算可以驻留在通过计算机系统或网络连接的单独的计算单元上。仅仅是举例，硬盘驱动器可以从不可移动、非易失性磁介质读取或向其写入，磁盘驱动器可以从不可移动、非易失性磁盘读取或向其写入，并且光盘驱动器可以从诸如CDROM或其它光学介质的不可移动、非易失性光盘读取或向其写入。其它能在示例性操作环境中使用的可移动/不可移动、易失性/非易失性计算机存储介质可包括但不限于磁带盒、闪存卡、数字多用盘、数字录像带、固态RAM、固态ROM等。以上讨论的驱动器及其相关的计算机存储介质为计算单元提供计算机可读指令、数据结构、程序模块和其它数据的存储。

客户可以将命令和信息通过客户端接口输入计算单元中，客户端接口可以是诸如键盘和定点装置(通常被称为鼠标、轨迹球或触控板)的输入装置。输入装置可包括麦克风、操纵杆、卫星天线、扫描仪等。这些和其它输入装置通常通过联接到系统总线的客户端接口连接到处理单元，但也可通过诸如并行端口或通用串行总线(USB)的其它接口和总线结构来连接。

监视器或其它类型的显示装置可以经由诸如视频接口的接口连接到系统总线。图形用户接口(“GUI”)也可与视频接口一起使用以从客户端接口接收指令并将指令传输到处理单元。除了监视器之外，计算机还可包括诸如扬声器和打印机的其它外围输出装置，这些装置可通过输出外围接口连接。

虽然未示出计算单元的许多其它内部部件，但本领域的普通技术人员将了解，这样的部件和其互连是熟知的。

虽然已结合目前优选的实施方案描述了本发明，但本领域的技术人员将理解，并非意图将本发明限制到这些实施方案。因此，可以构思到，在不脱离由所附权利要求及其等同形式限定的本发明的精神和范围的情况下，可以对所公开的实施方案提出各种备选的实施方案和修改。

Claims (20)

1.一种用于优化按改造储层体积计的井产量的方法，其包括：

输入一个或多个复杂储层性质和一个或多个复杂裂缝网络性质，所述复杂裂缝网络性质包括对应于复杂裂缝网络模型中的簇的数据；

基于所述复杂储层性质和所述复杂裂缝网络性质来确定平均有效裂缝长度的分布；

使用计算机处理器从平均有效裂缝长度的所述分布对平均有效裂缝长度抽样；以及

使用平均有效裂缝长度的所述分布和所述抽样的平均有效裂缝长度通过历史拟合来优化井产量，以改善所述改造储层体积的渗透率。

2.根据权利要求1所述的方法，其中所述历史拟合由单井储层模拟器执行。

3.根据权利要求1所述的方法，其中平均有效裂缝长度的所述分布通过以下方式确定：

读取针对每个井的每个压裂阶段中的每个裂缝平面的有效裂缝长度；

使用用于相应的压裂阶段的每个有效裂缝长度来计算每个压裂阶段的所述平均有效裂缝长度；以及

通过将每个相应的压裂阶段的所述平均有效裂缝长度与每个储层或井测井性质关联来构建平均有效裂缝长度的所述分布。

4.根据权利要求3所述的方法，其中平均有效裂缝长度的所述分布是通过下式构建的离散条件分布：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{\mathrm{Pr}ob(P=p\∩X_{eff}={\hat{x}}_{eff,s}^w)}{\mathrm{Pr}ob(P=p)}$

或通过下式构建的连续条件分布：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{{\mathrm{Prob}}_{P,X_{eff}}(p,{\hat{x}}_{eff,s}^w)}{{\mathrm{Prob}}_P\left(p\right)}$

5.根据权利要求3所述的方法，其中每个裂缝平面的最长轴线被读取并指定为每个相应的裂缝平面的所述有效裂缝长度。

6.根据权利要求3所述的方法，其中每个压裂阶段的所述平均有效裂缝长度通过下式计算：

${\hat{x}}_{eff,s}^w=\frac{1}{F}\underset{f=1}{\overset{F}{\Σ}}{x}_{eff,s,f}^w$

7.根据权利要求3所述的方法，其中每个储层或井测井性质为复杂储层性质。

8.根据权利要求3所述的方法，其中每个井的每个压裂阶段包括多个裂缝平面，在相应的压裂阶段内的每个裂缝平面具有不同的有效裂缝长度。

9.一种非暂时性程序载体装置，所述非暂时性程序载体装置有形地承载计算机可执行的指令以用于优化按改造储层体积计的井产量，所述指令为可执行的以实施：

输入一个或多个复杂储层性质和一个或多个复杂裂缝网络性质，所述复杂裂缝网络性质包括对应于复杂裂缝网络模型中的簇的数据；

基于所述复杂储层性质和所述复杂裂缝网络性质来确定平均有效裂缝长度的分布；

从平均有效裂缝长度的所述分布对平均有效裂缝长度抽样；以及

使用平均有效裂缝长度的所述分布和所述抽样的平均有效裂缝长度通过历史拟合来优化井产量，以改善所述改造储层体积的渗透率。

10.根据权利要求9所述的程序载体装置，其中所述历史拟合由单井储层模拟器执行。

11.根据权利要求9所述的程序载体装置，其中平均有效裂缝长度的所述分布通过以下方式确定：

读取针对每个井的每个压裂阶段中的每个裂缝平面的有效裂缝长度；

使用用于相应的压裂阶段的每个有效裂缝长度来计算每个压裂阶段的所述平均有效裂缝长度；以及

通过将每个相应的压裂阶段的所述平均有效裂缝长度与每个储层或井测井性质关联来构建平均有效裂缝长度的所述分布。

12.根据权利要求11所述的程序载体装置，其中平均有效裂缝长度的所述分布是通过下式构建的离散条件分布：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{\mathrm{Pr}ob(P=p\∩X_{eff}={\hat{x}}_{eff,s}^w)}{\mathrm{Pr}ob(P=p)}$

或通过下式构建的连续条件分布：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{{\mathrm{Prob}}_{P,X_{eff}}(p,{\hat{x}}_{eff,s}^w)}{{\mathrm{Prob}}_P\left(p\right)}$

13.根据权利要求11所述的程序载体装置，其中每个裂缝平面的最长轴线被读取并指定为每个相应的裂缝平面的所述有效裂缝长度。

14.根据权利要求11所述的程序载体装置，其中每个压裂阶段的所述平均有效裂缝长度通过下式计算：

${\hat{x}}_{eff,s}^w=\frac{1}{F}\underset{f=1}{\overset{F}{\Σ}}{x}_{eff,s,f}^w$

15.根据权利要求11所述的程序载体装置，其中每个储层或井测井性质为复杂储层性质。

16.根据权利要求11所述的程序载体装置，其中每个井的每个压裂阶段包括多个裂缝平面，在相应的压裂阶段内的每个裂缝平面具有不同的有效裂缝长度。

17.一种用于优化按改造储层体积计的井产量的方法，所述方法包括：

输入一个或多个复杂储层性质和一个或多个复杂裂缝网络性质，所述复杂裂缝网络性质包括对应于复杂裂缝网络模型中的簇的数据；

通过以下方式确定平均有效裂缝长度的分布：

读取针对每个井的每个压裂阶段中的每个裂缝平面的有效裂缝长度；

使用用于相应的压裂阶段的每个有效裂缝长度来计算每个压裂阶段的平均有效裂缝长度；以及

通过将每个相应的压裂阶段的所述平均有效裂缝长度与每个储层或井测井性质关联来构建平均有效裂缝长度的所述分布；

从平均有效裂缝长度的所述分布对所述平均有效裂缝长度抽样；以及

使用平均有效裂缝长度的所述分布和所述抽样的平均有效裂缝长度通过历史拟合来优化井产量，以改善所述改造储层体积的渗透率。

18.根据权利要求17所述的方法，其中平均有效裂缝长度的所述分布是通过下式构建的离散条件分布：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{\mathrm{Pr}ob(P=p\∩X_{eff}={\hat{x}}_{eff,s}^w)}{\mathrm{Pr}ob(P=p)}$

或通过下式构建的连续条件分布：

${\tilde{x}}_{eff,s|p}^w=\mathrm{Pr}ob(X_{eff}={\hat{x}}_{eff,s}^w|P=p)=\frac{{\mathrm{Prob}}_{P,X_{eff}}(p,{\hat{x}}_{eff,s}^w)}{{\mathrm{Prob}}_P\left(p\right)}$

19.根据权利要求17所述的方法，其中每个裂缝平面的最长轴线被读取并指定为每个相应的裂缝平面的所述有效裂缝长度。

20.根据权利要求17所述的方法，其中每个压裂阶段的所述平均有效裂缝长度通过下式计算：

${\hat{x}}_{eff,s}^w=\frac{1}{F}\underset{f=1}{\overset{F}{\Σ}}{x}_{eff,s,f}^w$