CN105474046A - 产生虚拟生产测井工具剖面以用于改进历史拟合 - Google Patents

Abstract

用于产生虚拟生产测井工具剖面以用于改进历史拟合并主动控制智能井的系统和方法。在优选实施方案中，建立储层性质的概率密度函数(pdf)，从所述pdf中选择所述储层的样本，使用所述样本执行正态分数变换，通过重新网格化来确定所述储层的概率分布，运行模型以产生含水率和液体速率剖面(Qliq)，以及对每一流入控制装置求解液体速率剖面。公开了用以确定储层模型的液体速率剖面的较快和较准确的系统和方法。

Description

产生虚拟生产测井工具剖面以用于改进历史拟合

相关申请的交叉参考

此处要求2013年9月9日提交的美国临时专利申请号61/875,591的优先权，且其说明书以引用的方式并入本文中。

关于联邦资助研究的声明

不适用。

技术领域

本公开一般来说涉及用于产生虚拟生产测井工具剖面以用于改进历史拟合的系统和方法。更具体地说，本公开涉及产生虚拟生产测井工具剖面以用于改进历史拟合并主动控制智能井。

背景技术

用于在智能完井的过程中使用流入控制装置(ICD)或流入控制阀(ICV)来估计液体速率剖面(Qliq)的常规方法需要连续历史拟合的多个步骤。在此方法的一个实例中，将每一ICD段的生产和注入数据(例如，井生产/注入速率、水饱和度(S_w))和井口处的地面数据(例如，压力(p)、温度(T)、液体速率剖面(Qliq)和含水率(watercut))用作输入。ICD段对应于由给定ICD控制的完井的长度。地面数据用以更新井模型，接着运行井模型以计算更新的操作点(p_u,T_u)。使用本领域中众所周知的标准失配最小化技术来执行局部历史拟合，且确定新的操作点(p_n,T_n)，其对应于地面数据与井模型数据之间的最小化失配。新的操作点(p_n,T_n)用以初始化和运行水力模型，其计算生产和压力测井剖面。使用本领域中众所周知的标准失配最小化技术来执行生产测井工具(PLT)数据的历史拟合以计算新的生产和压力测井剖面。新的生产和压力测井剖面由储存模型使用以使用本领域中众所周知的标准失配最小化技术来对含水率剖面和气油比进行历史拟合。

在上述实例中，过程是耗时的，这是因为它需要标准历史拟合的三个连续步骤。此外，过程在每一ICD段的液体速率剖面(Qliq)方面提供次优结果，这是因为它不考虑与井靠得很近的储层参数(例如，网格单元渗透率)的分布的不确定性。而且，过程在每一ICD段的液体速率剖面(Qliq)方面提供次优结果，这是因为它不考虑储层模型中的储层参数(例如，网格单元渗透率)的最佳解答。

附图说明

下文参看附图描述本公开，附图中相同元件用相同参考数字参考，且其中：

图1是说明用于实现本公开的方法的一个实施方案的流程图。

图2是说明用于执行图1中的步骤104的方法的一个实施方案的流程图。

图3A是说明因为图2中的步骤206产生的作为沿着水平井轨道的储层参数的渗透率的概率密度函数的显示。

图3B是说明因为图1中的步骤110的一次或多次迭代产生的图3A中的沿着水平井轨道的渗透率概率密度函数的分布的显示。

图4是说明用于实现本公开的计算机系统的一个实施方案的框图。

具体实施方式

因此，本公开通过提供用于产生虚拟生产测井工具剖面以用于改进历史拟合并主动控制智能井的系统和方法来克服现有技术中的一个或多个缺陷。

在一个实施方案中，本公开包括一种用于产生虚拟生产测井工具剖面的方法，其包括：a)使用计算机系统建立与一个或多个地质单元模型和相关数据相关联的储层性质的概率密度函数(pdf)；b)从pdf中选择储层性质的样本；c)使用选定的样本对储层性质的pdf执行正态分数变换；d)通过使用对储层性质的pdf执行的正态分数变换进行局部重新网格化来确定储层性质的概率分布；e)使用概率分布运行储层模型以产生含水率和液体速率剖面(Qliq)；f)使用含水率和(Qliq)来定义目标函数；g)定义一个或多个优化约束；以及h)使用目标函数和一个或多个优化约束来对每一流入控制装置段求解(Qliq)，其表示虚拟生产测井工具剖面。

在另一实施方案中，本公开包括一种非暂时性程序载体装置，其有形地携载用于产生虚拟生产测井工具剖面的计算机可执行指令，指令可被执行以实现以下操作：a)建立与一个或多个地质单元模型和相关数据相关联的储层性质的概率密度函数(pdf)；b)从pdf中选择储层性质的样本；c)使用选定的样本对储层性质的pdf执行正态分数变换；d)通过使用对储层性质的pdf执行的正态分数变换进行局部重新网格化来确定储层性质的概率分布；e)使用概率分布运行储层模型以产生含水率和液体速率剖面(Qliq)；f)使用含水率和(Qliq)来定义目标函数；g)定义一个或多个优化约束；以及h)使用目标函数和一个或多个优化约束来对每一流入控制装置段求解(Qliq)，其表示虚拟生产测井工具剖面。

在又一实施方案中，本发明包括一种非暂时性程序载体装置，其有形地携载用于产生虚拟生产测井工具剖面的计算机可执行指令，指令可被执行以实现以下操作：a)建立与一个或多个地质单元模型相关联的储层性质的概率密度函数(pdf)；b)从pdf中选择储层性质的样本；c)使用选定的样本对储层性质的pdf执行正态分数变换；d)通过使用对储层性质的pdf执行的正态分数变换进行局部重新网格化来确定储层性质的概率分布；e)使用概率分布运行储层模型以产生含水率和液体速率剖面(Qliq)；f)使用基于含水率和(Qliq)的目标函数和一个或多个优化约束来对每一流入控制装置段求解(Qliq)，其表示虚拟生产测井工具剖面；以及g)重复步骤b)至f)直到虚拟生产测井工具剖面被优化为止。

具体地描述了本公开的主题，然而，描述本身无意限制本公开的范围。因此，标的也可结合其它技术用其它方式体现以包括不同步骤或与本文中描述的步骤类似的步骤的组合。此外，尽管术语“步骤”在本文中可用以描述使用的方法的不同元素，但术语不应解释为暗示本文中公开的各步骤当中或之间的任何特定次序，除非另外明确地受特定次序的描述所限制。尽管以下描述指石油和天然气行业，但本公开的系统和方法不限于此且也可应用于其它行业以实现类似结果。

方法描述

现在参看图1，说明了用于实现本公开的方法100的一个实施方案的流程图。

在步骤102中，输入因为使用本领域中众所周知的技术运行井模型而产生的表示压力(p)和温度(T)的操作点以及因为使用本领域中众所周知的技术运行水力模型而产生的压力(例如，井底压力)、生产(例如，气油比)剖面和注入速率以用于储层模型。

在步骤104中，使用操作点(p,T)、压力、生产剖面和本领域中众所周知的技术来运行储层模型以产生含水率、液体速率剖面(Qliq)和饱和度剖面(Sw)。

在步骤106中，使用含水率、液体速率剖面(Qliq)和本领域中众所周知的技术来定义目标函数。将目标函数定义为建模的与测量的参数之间的失配，此处为含水率(wc)＝f(k,Qliq)，其中f表示非线性关系。因此，优化是优选的。控制变量对应于动态数据(总Qliq)和静态数据(例如，储层参数，像渗透率(k)。目标是使储层模型与动态数据一致以同时将井含水率减到最小。

可将具有相关联的统计模式的控制变量、平均值(μ)和标准偏差(σ)表示为：

·渗透率k(μ_k,σ_k)

·Qliq(μ_Q,σ_Q)

因此，可将定义的目标函数表示为：

$\begin{array}{c}OF=\underset{m=1}{\overset{M}{\Σ}}\[\frac{w_{1m}}{s_{1m}}{({\mathrm{\μ}}_{ym}-T_m)}^2+\frac{w_{2m}}{s_{2m}}{\mathrm{\σ}}_{ym}^2\]\\ =\[(\frac{w_{1k}}{s_{1k}}{\left({\mathrm{\μ}}_k-T_k\right)}^2+\frac{w_{2k}}{s_{2k}}{\mathrm{\σ}}_k^2)+(\frac{w_{1Q}}{s_{1Q}}{\left({\mathrm{\μ}}_Q-T_Q\right)}^2+\frac{w_{2Q}}{s_{2Q}}{\mathrm{\σ}}_Q^2)\]\end{array}---\left(1\right)$

其中(μ)对应于测量的平均值且(T)对应于建模的(目标)平均值。权数(w₁)和(w₂)和比例因子(s₁)和(s₂)对应于“目标平均值”和“最小化变化”目标分量。在应用的第一近似中(s_1k＝s_2k＝s_1Q＝s_1Q＝1(归一化为统一的)且w_1k＝w_2k＝w_1Q＝w_1Q＝f(d))，其中(d)对应于来自步骤122中的局部重新网格化的网格单元分割器。重要的是应注意，当在正态分数变换域中操作时，μ_k＝0且σ_k＝1。

在步骤108中，使用注入速率、操作点(p,T)和本领域中众所周知的技术来定义一个或多个优化约束。变化约束可表示为：p+nσ_p≤p_i≤p-nσ_p和T+nσ_T≤T_i≤T–nσ_T，其中(p_i)和(T_i)对应于测量的压力或温度的第i个样本。平均值约束可表示为：g(μ_p,μ_T)≥0。

在步骤110中，使用众所周知的随机优化技术中的目标函数和优化约束来对每一ICD段的液体速率剖面(Qliq)求解。众所周知的随机优化技术的实例包括马尔科夫链蒙特卡罗或自适应模拟退火，其中目标函数可用高斯形式定义。结果表示虚拟生产测井工具剖面，下文中称作VPLT剖面。VPLT剖面保证每一井的总含水率分布的最小化失配，其可表示为：Qliq_{|每一ICD优化的}＝Qliq_{|wc_失配＝最小值}。在图3B中，举例来说，将图3A中的沿着水平井轨道的渗透率概率密度函数的分布的显示300B说明为此步骤的一次或多次迭代的结果。分布308对应于在步骤114至122之后的储层性质的分布，其受井轨道的由个别ICD控制的部分所约束。因此，通过沿着井轨道累积每一相应ICD段306内的储层性质值来建立每一分布308。图3B中的每一ICD段306分隔分布308且将图3A中的网格单元302组合为由ICD段306表示的较大伪单元。通过参看步骤122描述的网格单元分割器来确定组合为较大伪单元的网格单元302的数目。

在步骤112中，方法100通过将来自步骤110的上一VPLT剖面与表示最小含水率和气油比的预定VPLT剖面进行比较来确定VPLT剖面是否经过优化。如果来自步骤110的上一VPLT剖面小于或等于预定VPLT剖面，那么VPLT剖面经过优化。如果VPLT剖面经过优化，那么方法100结束。如果VPLT剖面未经过优化，那么方法100进行至步骤118。

在步骤114中，输入一个或多个地质单元模型(n)和相关数据。地质单元模型表示组合用储层性质(诸如，例如，渗透率(k))的(i,j,k)地质单元网格的静态地质单元模型。因此，每一地质单元模型包括一个或多个网格单元(c)。相关数据可包括(例如)一个或多个井轨道(t)作为动态/模拟模型的一部分，动态/模拟模型描绘研究中的井的(i,j,k)路径。

在步骤116中，建立或更新与步骤114中输入的数据相关联的储层性质概率密度函数(pdf)。参看图2进一步描述用于执行此步骤的方法的一个实施方案。

在步骤118中，从步骤116中建立或更新的pdf中选择储层性质的样本。

在步骤120中，使用在步骤118中选择的样本和本领域中众所周知的技术对在步骤116中建立或更新的储层性质的pdf执行正态分数变换。对于所有值i＝(1,…,N)以升序对的值进行排序。接着使用来计算用于秩(m)的观察的累积频率或p_m分位数，其中(w)是在步骤118中选择的具有秩(m)的样本的权数。如果样本的权数(w)不可获得，那么使用默认权数(w_m＝1/N)。具有秩(m)的样本的正态分数变换是标准正态分布的p_m分位数，其表示为：其中G(.)是累积标准正态分布。为了简单起见，可参考正态分数变换而不指定秩(m)(即，)。

在步骤122中，通过使用来自步骤120的正态分数变换进行局部重新网格化来确定储层性质的概率分布。在网格模型区域中执行局部重新网格化，其中保证网格块的更新的密度以便呈现与ICD网格分段一致的储层性质的解且呈现相关估计量的较高准确度和改进的取样统计。相应地更改在装备有ICD的井附近的网格块解以对井筒附近的压力或流体流动的显著变化建模。首先选择体现所有ICD子段的关注体积。定义网格单元分割器(d)(对于每一I,J,K网格单元位置或网格单元的子集)。子集的大小对应于ICD段的大小。对于与I,J,K网格单元位置相关联的每一正态分数变换值应用对应位置的网格单元分割器。这产生网格单元与相关联的正态分数变换值的局部相乘。局部重新网格化保证高斯域中的无偏估计量。局部重新网格化的结果是通过与ICD段相关联的分割器加权的新分布。因此，正态分数变换变成其用作另一输入以运行步骤104中的储层模型。步骤104的下一迭代将运行储层模型以产生另一含水率、液体速率剖面(Qliq)和将含水率的失配减到最小的饱和度剖面(S_w)。

方法100产生VPLT剖面以用于主动控制智能井。每一VPLT剖面表示与完井剖面(也称作流入控制阀ICD段)拟合的优化的液体速率剖面(Qliq)。因此，方法100实时地优化井下阀设置以通过减小含水率和/或气油比来最大化石油采收率。通过计算每一ICD段的水流速，可促进任何类型的注水优化。

现在参看图2，说明了用于执行图1中的步骤116的方法200的一个实施方案的流程图。

在步骤201中，从步骤114中输入的地质单元模型的总数中自动地选择或者可使用参看图4进一步描述的客户端接口和/或视频接口来选择地质单元模型(n)。

在步骤202中，从步骤114中输入的井轨道的总数中自动地选择或者可使用参看图4进一步描述的客户端接口和/或视频接口来选择选定的地质单元模型的井轨道(t)。

在步骤203中，从步骤114中输入的网格单元的总数中自动地选择或者可使用参看图4进一步描述的客户端接口和/或视频接口来选择选定的储层模型的网格单元(c)。

在步骤204中，使用本领域中众所周知的技术识别选定的地质单元模型(n)、井轨道(t)和网格单元(c)的储层性质(k)。

在步骤205中，使用步骤114中输入的数据和本领域中众所周知的技术来产生储层性质(k)的一个或多个直方图(h_k ^t(c))。

在步骤206中，使用直方图(h_k ^t(c))、以下等式：p_k ^t(c)＝h_k ^t(c)/c和本领域中众所周知的技术来建立或更新储层性质(k)的概率密度函数(pdf)。在图3A中，举例来说，将沿着水平井轨道的作为储层性质(k)的渗透率的概率密度函数的显示说明为此步骤的结果。分布304等效于在步骤201中选择的地质单元模型(n)的网格解。因此，每一分布304的大小对应于来自地质单元模型(n)的相应网格单元302。因为图3A中的每一网格单元302小于图3B中的ICD段306表示的伪单元，所以图3A中的每一分布304的实现的统计小于图3B中的每一分布308的实现的统计。因此，将分布从分布304改进至分布308的图1的步骤122中的局部重新网格化减少相关联的不确定性。

在步骤208中，方法200确定选定的储层模型的剩余网格单元的总数中是否有另一网格单元(c)供选择。如果有另一网格单元(c)供选择，那么方法200返回至步骤203。如果没有另一网格单元(c)供选择，那么方法200进行至步骤209。

在步骤209中，方法200确定选定的储层模型的剩余井轨道的总数中是否有另一井轨道(t)供选择。如果有另一井轨道(t)供选择，那么方法200返回至步骤202。如果没有另一井轨道(t)供选择，那么方法200进行至步骤210。

在步骤210中，方法200确定剩余储层模型的总数中是否有另一储层模型(n)供选择。如果有另一储层模型(n)供选择，那么方法200返回至步骤201。如果没有另一储层模型(n)供选择，那么方法200将储层性质(k)的更新的pdf传回至步骤120。

系统描述

本公开可通过计算机可执行指令程序实现，诸如一般来说称作由计算机执行的软件应用或应用程序的程序模块。软件可包括(例如)执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件和数据结构。软件形成允许计算机根据输入源进行反应的接口。用于生产的是由LandmarkGraphicsCorporation销售的商业软件应用，可用作实现本公开的接口应用。软件还可与其它代码段合作以响应于结合接收数据的源接收的数据而起始多个任务。软件可存储和/或携载在任何多个存储器上，诸如CD-ROM、磁盘、磁泡存储器和半导体存储器(例如，各种类型的RAM或ROM)。此外，软件和其结果可经由多种携载媒体(诸如光纤、金属线)和/或通过多种网络中的任一者(诸如，因特网)传输。

此外，本领域技术人员将了解，本公开可借助多种计算机系统配置实践，包括手持式装置、多处理器系统、基于微处理器的或可编程消费者电子装置、小型计算机、大型计算机，和类似物。可接受任何数目的计算机系统和计算机网络与本公开一起使用。本公开可在分布式计算机环境中实践，在所述环境中通过经由通信网络链接的远程处理装置来执行任务。在分布式计算环境中，程序模块可位于远程和本地计算机存储媒体(包括存储器存储装置)中。因此，本公开可结合各种硬件、软件或其组合在计算机系统或其它处理系统中实现。

现在参看图4，框图说明了用于在计算机上实现本公开的系统的一个实施方案。系统包括计算单元，有时称作计算系统，其含有存储器、应用程序、客户端接口、视频接口和处理单元。计算单元仅仅是合适的计算环境的一个实例，且无意提出关于本公开的使用范围或功能性的任何限制。

存储器主要存储应用程序，其也可描述为含有计算机可执行指令的程序模块，计算机可执行指令由计算单元执行以用于实现本文中描述的且在图1至图3中说明的本公开。因此，存储器包括VPLT剖面模块，其使得能够执行参看图1中的步骤106至122描述的方法。上述模块和应用可集成图4中说明的剩余应用程序的功能性。具体地说，用于生产的可用作接口应用以执行图1中的步骤102至104。尽管用于生产的可用作接口应用，但也可替代地使用其它接口应用，或者VPLT剖面模块可用作独立应用。

尽管计算单元示出为具有一般化存储器，但计算单元通常包括多种计算机可读媒体。举例来说(且不限于此)，计算机可读媒体可包括计算机存储媒体和通信媒体。计算系统存储器可包括计算机存储媒体，其呈易失性和/或非易失性存储器的形式，诸如只读存储器(ROM)和随机存取存储器(RAM)。基本输入/输出系统(BIOS)通常存储在ROM中，BIOS含有帮助(诸如在启动期间)在计算单元内的元件之间传送信息的基本例程。RAM通常含有处理单元可立即存取和/或目前正在处理单元上操作的数据和/或程序模块。举例来说(且不限于此)，计算单元包括操作系统、应用程序、其它程序模块，和程序数据。

存储器中所示的组件也可包括在其它可移动/不可移动、易失性/非易失性计算机存储媒体中，或其可在计算单元中通过应用程序接口(“API”)或云计算实现，应用程序接口(“API”)或云计算可驻留在通过计算机系统或网络连接的单独的计算单元上。仅举例来说，硬盘驱动器可从不可移动的非易失性磁性媒体读取或写入至非易失性磁性媒体，磁盘驱动器可从可移动的非易失性磁盘读取或写入至可移动的非易失性磁盘，且光盘驱动器可从可移动的非易失性光盘(诸如CDROM)或其它光学媒体读取或写入至可移动的非易失性光盘或其它光学媒体。可用于示例性操作环境中的其它可移动/不可移动、易失性/非易失性计算机存储媒体可包括(但不限于)磁带盒、快闪存储卡、数字通用光盘、数字录像带、固态RAM、固态ROM，和类似物。上文论述的驱动器和其相关联的计算机存储媒体提供对计算机可读指令、数据结构、程序模块和用于计算单元的其它数据的存储。

客户端可通过客户端接口将命令和信息键入至计算单元中，客户端接口可以是输入装置，诸如键盘和指向装置，其通常称作鼠标、轨迹球或触控板。输入装置可包括麦克风、操纵杆、碟形卫星天线、扫描仪，或类似物。这些和其它输入装置经常通过耦合至系统总线的客户端接口连接至处理单元，但可通过其它接口和总线结构(诸如并行端口或通用串行总线(USB))连接。

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

尽管未示出计算单元的许多其它内部组件，但本领域技术人员将了解，此类组件和其互连是众所周知的。

尽管已结合目前优选的实施方案描述了本公开，但本领域技术人员应理解，其无意将本公开限于那些实施方案。因此，预期可对公开的实施方案进行各种替代实施方案和修改，而不脱离由所附权利要求书和其等效物界定的本公开的精神和范围。

Claims (20)

1.一种用于产生虚拟生产测井工具剖面的方法，其包括：

a)使用计算机系统建立与一个或多个地质单元模型和相关数据相关联的储层性质的概率密度函数(pdf)；

b)从所述pdf中选择所述储层性质的样本；

c)使用所述选定的样本对所述储层性质的所述pdf执行正态分数变换；

d)通过使用对所述储层性质的所述pdf执行的所述正态分数变换进行局部重新网格化来确定所述储层性质的概率分布；

e)使用所述概率分布运行储层模型以产生含水率和液体速率剖面(Qliq)；

f)使用所述含水率和(Qliq)来定义目标函数；

g)定义一个或多个优化约束；以及

h)使用所述目标函数和所述一个或多个优化约束来对每一流入控制装置段求解(Qliq)，其表示所述虚拟生产测井工具剖面。

2.如权利要求1所述的方法，其进一步包括重复步骤b)至h)直到所述虚拟生产测井工具剖面被优化为止。

3.如权利要求1所述的方法，其进一步包括使用操作点、压力剖面和生产剖面来运行所述储层模型。

4.如权利要求3所述的方法，其中使用一个或多个注入速率和所述操作点来定义所述一个或多个优化约束。

5.如权利要求2所述的方法，其中当上一虚拟生产测井工具剖面小于或等于预定虚拟生产测井工具剖面时，所述虚拟生产测井工具经过优化。

6.如权利要求5所述的方法，其中所述预定虚拟生产测井工具剖面表示最小含水率和气油比。

7.如权利要求1所述的方法，其中通过以下操作建立所述pdf：

a)读取所述一个或多个地质单元模型中的一者、一个或多个井轨道中的一者和一个或多个网格单元中的一者的所述储层性质；

b)产生所述储层性质的一个或多个直方图；

c)使用所述一个或多个直方图建立所述储层性质的所述pdf；以及

d)针对所述一个或多个网格单元中的每一者、所述一个或多个井轨道中的每一者和所述一个或多个地质单元模型中的每一者重复步骤a)至c)，其中每一迭代更新所建立的所述pdf。

8.一种非暂时性程序载体装置，其有形地携载用于产生虚拟生产测井工具剖面的计算机可执行指令，所述指令可被执行以实现以下操作：

a)建立与一个或多个地质单元模型和相关数据相关联的储层性质的概率密度函数(pdf)；

b)从所述pdf中选择所述储层性质的样本；

c)使用所述选定的样本对所述储层性质的所述pdf执行正态分数变换；

d)通过使用对所述储层性质的所述pdf执行的所述正态分数变换进行局部重新网格化来确定所述储层性质的概率分布；

e)使用所述概率分布运行储层模型以产生含水率和液体速率剖面(Qliq)；

f)使用所述含水率和(Qliq)来定义目标函数；

g)定义一个或多个优化约束；以及

h)使用所述目标函数和所述一个或多个优化约束来对每一流入控制装置段求解(Qliq)，其表示所述虚拟生产测井工具剖面。

9.如权利要求8所述的程序载体装置，其进一步包括重复步骤b)至h)直到所述虚拟生产测井工具剖面被优化为止。

10.如权利要求8所述的程序载体装置，其进一步包括使用操作点、压力剖面和生产剖面来运行所述储层模型。

11.如权利要求10所述的程序载体装置，其中使用一个或多个注入速率和所述操作点来定义所述一个或多个优化约束。

12.如权利要求9所述的程序载体装置，其中当上一虚拟生产测井工具剖面小于或等于预定虚拟生产测井工具剖面时，所述虚拟生产测井工具经过优化。

13.如权利要求12所述的程序载体装置，其中所述预定虚拟生产测井工具剖面表示最小含水率和气油比。

14.如权利要求8所述的程序载体装置，其中通过以下操作建立所述pdf：

a)读取所述一个或多个地质单元模型中的一者、一个或多个井轨道中的一者和一个或多个网格单元中的一者的所述储层性质；

b)产生所述储层性质的一个或多个直方图；

c)使用所述一个或多个直方图建立所述储层性质的所述pdf；以及

d)针对所述一个或多个网格单元中的每一者、所述一个或多个井轨道中的每一者和所述一个或多个地质单元模型中的每一者重复步骤a)至c)，其中每一迭代更新所建立的所述pdf。

15.一种非暂时性程序载体装置，其有形地携载用于产生虚拟生产测井工具剖面的计算机可执行指令，所述指令可被执行以实现以下操作：

a)建立与一个或多个地质单元模型相关联的储层性质的概率密度函数(pdf)；

b)从所述pdf中选择所述储层性质的样本；

c)使用所述选定的样本对所述储层性质的所述pdf执行正态分数变换；

d)通过使用对所述储层性质的所述pdf执行的所述正态分数变换进行局部重新网格化来确定所述储层性质的概率分布；

e)使用所述概率分布运行储层模型以产生含水率和液体速率剖面(Qliq)；

f)使用基于所述含水率和(Qliq)的目标函数和一个或多个优化约束来对每一流入控制装置段求解(Qliq)，其表示所述虚拟生产测井工具剖面；以及

g)重复步骤b)至f)直到所述虚拟生产测井工具剖面被优化为止。

16.如权利要求15所述的程序载体装置，其进一步包括使用操作点、压力剖面和生产剖面来运行所述储层模型。

17.如权利要求16所述的程序载体装置，其中使用一个或多个注入速率和所述操作点来定义所述一个或多个优化约束。

18.如权利要求15所述的程序载体装置，其中当上一虚拟生产测井工具剖面小于或等于预定虚拟生产测井工具剖面时，所述虚拟生产测井工具经过优化。

19.如权利要求18所述的程序载体装置，其中所述预定虚拟生产测井工具剖面表示最小含水率和气油比。

20.如权利要求15所述的程序载体装置，其中通过以下操作建立所述pdf：

a)读取所述一个或多个地质单元模型中的一者、一个或多个井轨道中的一者和一个或多个网格单元中的一者的所述储层性质；

b)产生所述储层性质的一个或多个直方图；

c)使用所述一个或多个直方图建立所述储层性质的所述pdf；以及

d)针对所述一个或多个网格单元中的每一者、所述一个或多个井轨道中的每一者和所述一个或多个地质单元模型中的每一者重复步骤a)至c)，其中每一迭代更新所建立的所述pdf。