中的每一个位置X处，存在四个电阻率rl、r2、r3以 及r4。直方图的横坐标轴表示电阻率，而直方图的纵坐标轴表示频率。标绘了 4个层中的 每一个层的电阻率的分布。还提供了中值、标准偏差值以及众数（mode)值。例如，第一层 (层1)具有与电阻率rl相对应的、以70Ωπι为中心的分布，第二层（层2)具有与电阻率 r2相对应的、以5 Ωπι为中心的分布，第三层（层3)具有与电阻率r3相对应的、以40 Ωπι为 中心的分布，而第四层具有与电阻率r4相对应的、以95 Ω m为中心的分布。
[0034] 图2B示出了图1所示地球的一个位置x处的地球参数（层厚度h)的后验分布 Pr(p |d，I)的示例。横坐标轴上标绘了层厚度，而直方图的纵坐标轴表示频率。标绘了 4 个层中的每一个层的厚度的分布。例如，第一层（层1)具有与厚度hi相对应的、以300为 中心的分布，第二层（层2)具有与厚度h2相对应的、以1400为中心的分布，第三层（层3) 具有与厚度h3相对应的、以500为中心的分布。
[0035] 该后验分布Pr(P |d，I)可以被用于通过根据地球模型参数分布随机抽取来生成 地球的任何数量（N)实现，在图2A和2B例示。
[0036] 图3是根据本发明一实施例的、用于从多个（M个）可能的井位置中确定最佳井位 置的方法的流程图。该方法在S10,通过根据后验分布Pr (P I d，I)抽取N个地球模型而开 始，其中，该后验分布通过现有数据的随机反演来生成。该后验分布Pr(P |d，I)可以通过 任何随机反演技术来生成。例如，在一个实施例中，可以使用马尔可夫链蒙特卡罗（MCMC) 技术，其评估与给定某一观察数据d的地球参数P的矢量相关联的任何地球模型的似 然。随机反演的描述可以在 Trainor-Guitton&Hoversten (2011) "Stochastic inversion for electromagnetic geophysics!Practical challenges and improving convergence efficiency〃Geophysics 76，F373_F386中找到^概括地讲，随机反演对应于似然涵数 f(d| Pc^其描述了给定所观察数据d的候选模型PC出现的概率。在一个实施例中，似然 函数f(d| PcO可以用下面的方程（1)来表达。如果进行有关错误模型的性质的不同假定， 则可以获得似然函数的其它定义。
【主权项】
1. 一种用于从多个可能的井位置中确定最佳井位置的计算机实现的方法，所述方法包 括： 利用计算机，根据后验分布抽取多个地球模型，其中，所述后验分布通过现有数据的随 机反演来生成； 通过所述计算机，利用所述多个地球模型中的一地球模型内的多个提议井位置处的井 生产量与地球参数之间的关系，计算所述井生产量； 通过所述计算机，利用井成本与所述地球参数之间的关系，根据所述多个地球模型计 算成本分布；以及 通过所述计算机，利用根据位置相关随机反演的概率作为权重，计算所述提议井位置 的概率加权值。
2. 根据权利要求1所述的方法，还包括：在计算所述多个提议井位置处的井生产量之 前，在给定与生产量数据相关联的多个井的情况下，导出所述井生产量与所述地球参数之 间的关系。
3. 根据权利要求2所述的方法，其中，导出所述井生产量与地球参数之间的关系包括： 根据现有的同位置处的生产量数据和所述地球参数来导出所述关系。
4. 根据权利要求1所述的方法，还包括：在计算成本分布之前，导出所述井成本与所述 地球参数之间的关系。
5. 根据权利要求4所述的方法，其中，导出所述井成本与所述地球参数之间的关系包 括：输入一关系，在所述关系中，所述井成本本质上取决于深度参数、取决于空间可变性、岩 石硬度、或其任何组合。
6. 根据权利要求1所述的方法，其中，利用随机反演生成所述后验分布包括：利用马尔 可夫链蒙特卡罗技术来生成所述后验分布。
7. 根据权利要求1所述的方法，其中，利用随机反演生成所述后验分布包括：确定给定 所述现有数据的情况下的地球模型出现的概率。
8. 根据权利要求1所述的方法，其中，所述地球参数从由电阻率、渗透率、密度、孔隙 度、以及深度组成的组中选择。
9. 根据权利要求1所述的方法，其中，计算所述提议井位置的所述概率加权值包括：计 算给定井位置处的生产量一成本值，将所述生产量一成本值乘以根据所述随机反演获取的 所述后验分布以获取加权生产量一成本值，以及在所述多个后验地球模型上对所述加权生 产量一成本值求和。
10. 根据权利要求9所述的方法，其中，将所述概率加权值相互比较，以确定与井的最 有价值的位置相对应的最大加权值。
11. 根据权利要求9所述的方法，其中，选择多个地球模型，以使所述概率加权值的均 值和方差随着所述地球模型的数量增加而渐近恒定值。
12. -种用于从多个可能的井位置中确定最佳井位置的计算机系统，所述计算机系统 包括： 计算机可读存储器，被设置成存储井位置和现有井数据；和 处理器，所述处理器与所述计算机可读存储器通信，所述处理器被设置成： 根据后验分布来抽取多个地球模型，其中，所述后验分布通过现有数据的随机反演来 生成； 利用所述多个地球模型中的一地球模型内的多个提议井位置处的井生产量与地球参 数之间的关系，计算所述井生产量； 利用井成本与所述地球参数之间的关系，根据所述多个地球模型计算成本分布；以及 利用根据位置相关随机反演的概率作为权重，计算所述提议井位置的概率加权值。
13. 根据权利要求12所述的计算机系统，其中，所述处理器还被设置成，在给定与生产 量数据相关联的多个井的情况下，导出所述井生产量与所述地球参数之间的关系。
14. 根据权利要求13所述的计算机系统，其中，所述处理器还被设置成根据现有的同 位置处的生产量数据和所述地球参数来导出所述关系。
15. 根据权利要求12所述的计算机系统，其中，所述处理器还被设置成导出所述井成 本与所述地球参数之间的关系。
16. 根据权利要求15所述的计算机系统，其中，所述处理器还被设置成导出所述井成 本与所述地球参数之间的关系，其中，所述井成本本质上取决于深度参数、取决于空间可变 性、岩石硬度、或其任何组合。
17. 根据权利要求12所述的计算机系统，其中，所述处理器还被设置成利用马尔可夫 链蒙特卡罗技术生成所述后验分布。
18. 根据权利要求12所述的计算机系统，其中，所述处理器还被设置成确定给定所述 现有数据的情况下的地球模型出现的概率。
19. 根据权利要求12所述的计算机系统，其中，所述处理器还被设置成从由电阻率、渗 透率、密度、孔隙度、以及深度组成的组中选择所述地球参数。
20. 根据权利要求12所述的计算机系统，其中，所述处理器还被设置成如下计算所述 提议井位置的所述概率加权值：计算给定井位置处的生产量一成本值，将所述生产量一成 本值乘以根据所述随机反演获取的所述后验分布以获取加权生产量一成本值，以及在所述 多个后验地球模型上对所述加权生产量一成本值求和。
21. 根据权利要求20所述的计算机系统，其中，将所述概率加权值相互比较，以确定与 井的最有价值的位置相对应的最大加权值。
22. 根据权利要求12所述的计算机系统，其中，选择多个地球模型，以使所述概率加权 值的均值和方差随着所述地球模型的数量增加而渐近恒定值。
【专利摘要】在此描述了一种用于从多个可能的井位置中确定最佳井位置的系统和计算机实现的方法。该方法包括以下步骤：根据后验分布来抽取多个地球模型，其中，该后验分布通过现有数据的随机反演来生成；利用所述多个地球模型中的一地球模型内的多个提议井位置处的井生产量与地球参数之间的关系，计算该井生产量；利用井成本与该地球参数之间的关系，根据所述多个地球模型计算成本分布；以及利用根据位置相关随机反演的概率作为权重来计算该提议井位置的概率加权值。
【IPC分类】G06N7-00
【公开号】CN104838393
【申请号】CN201380056293
【发明人】G·M·霍弗斯滕, W·特拉尹诺-古顿
【申请人】雪佛龙美国公司, 劳伦斯利弗摩尔国际安全有限责任公司
【公开日】2015年8月12日
【申请日】2013年8月29日
【公告号】CA2883288A1, EP2891114A2, US20140067322, WO2014036300A2, WO2014036300A3