• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 20
    • 下载费用:30 金币  

    重庆时时彩的图表工具: 使用松弛变量确定网络中的有效约束的系统和方法.pdf

    关 键 词:
    使用 松弛 变量 确定 网络 中的 有效 约束 系统 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN200880118396.2

    申请日:

    2008.11.25

    公开号:

    CN102160050A

    公开日:

    2011.08.17

    当前法律状态:

    终止

    有效性:

    无权

    法律详情: 未缴年费专利权终止IPC(主分类):G06F 17/11申请日:20081125授权公告日:20131120终止日期:20171125|||授权|||实质审查的生效IPC(主分类):G06F 17/11申请日:20081125|||公开
    IPC分类号: G06F17/11 主分类号: G06F17/11
    申请人: 兰德马克绘图国际公司,哈里伯顿公司
    发明人: 詹姆斯·威廉·Ⅲ·沃茨; 格雷汉姆·弗莱明; 琴·陆
    地址: 美国德克萨斯州
    优先权: 2007.11.28 US 60/990,823
    专利代理机构: 北京路浩知识产权代理有限公司 11002 代理人: 张晶
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN200880118396.2

    授权公告号:

    |||102160050B||||||

    法律状态公告日:

    2018.11.09|||2013.11.20|||2011.09.28|||2011.08.17

    法律状态类型:

    专利权的终止|||授权|||实质审查的生效|||公开

    摘要

    用于确定地面设施网络中的有效约束的系统和方法,其包括在系统方程中使用松弛变量和乘数,以消去无关(无效)约束。

    权利要求书

    1: 一种用于确定网络中有效约束方程的方法, 包括 : 选择所述网络中连接的至少一个约束方程, 每一个约束方程包括各自的松弛变量和各 自的松弛变量乘数 ; 选择所述网络中每一个连接的基础方程, 具有所述至少一个约束方程的所述网络连接 的所述基础方程包括每一个约束方程的所述各自的松弛变量和另一个各自的松弛变量乘 数; 求解每一个各自的松弛变量 ; 根据每一个各自的松弛变量的解, 对于每一个各自的松弛变量乘数和每一个各自的另 一个松弛变量乘数, 在两个预定值之间调整各自的参数 ; 重复求解每一个各自的松弛变量以及对于每一个各自的松弛变量乘数和每一个各自 的另一个松弛变量乘数调整各自的参数的步骤, 直至每一个各自的参数等于所述两个预定 值中的一个 ; 以及 根据每一个各自的参数的所述两个预定值中的一个, 确定每一个约束方程是否有效。
    2: 如权利要求 1 所述的方法, 还包括选择所述网络中每一个节点的质量平衡方程。
    3: 如权利要求 1 所述的方法, 其中, 具有所述至少一个约束方程的所述网络连接的所 述基础方程是液压方程, 并且所述至少一个约束方程是速率约束方程。
    4: 如权利要求 1 所述的方法, 其中, 所述两个预定值是 0 和 1。
    5: 如权利要求 1 所述的方法, 还包括选择所述网络中另一个连接的另一个约束方程。
    6: 如权利要求 4 所述的方法, 其中, 将每一个各自的参数初始设置为 0.5, 并且如果各 自的松弛变量的所述解小于或等于零, 将所述每一个各自的参数调整接近 0, 并且如果所述 各自的松弛变量的解大于零, 将所述每一个各自的参数调整接近 1。
    7: 如权利要求 4 所述的方法, 其中, 如果所述每一个约束方程的各自的参数的值等于 1, 所述网络中每一个约束方程是有效的。
    8: 如权利要求 6 所述的方法, 其中, 通过将一减去所述参数乘以减小因子, 将每一个各 自的参数调整接近 1, 并且通过将所述参数乘以另一个减小因子, 将每一个各自的参数调整 接近 0。
    9: 如权利要求 7 所述的方法, 还包括 : 如果所述每一个约束方程的各自的参数等于 1, 将所述每一个约束方程的各自的松弛 变量加至具有所述至少一个约束方程的所述网络连接的所述基础方程。
    10: 如权利要求 2 所述的方法, 还包括 : 求解所述网络中每一个有效约束方程 ; 求解所述网络中每一个质量平衡方程 ; 以及 对于没有约束方程的每一个网络连接求解所述网络中每一个基础方程。
    11: 一种用于确定网络中有效约束方程的具有计算机可执行指令的计算机可读介质, 包括 : 选择所述网络中连接的至少一个约束方程, 每一个约束方程包括各自的松弛变量和各 自的松弛变量乘数 ; 选择所述网络中每一个连接的基础方程, 具有所述至少一个约束方程的所述网络连接 的所述基础方程包括每一个约束方程的所述各自的松弛变量和另一个各自的松弛变量乘 2 数; 求解每一个各自的松弛变量 ; 根据每一个各自的松弛变量的解, 对于每一个各自的松弛变量乘数和每一个各自的另 一个松弛变量乘数, 在两个预定值之间调整各自的参数 ; 重复求解每一个各自的松弛变量以及对于每一个各自的松弛变量乘数和每一个各自 的另一个松弛变量乘数调整各自的参数的步骤, 直至每一个各自的参数等于所述两个预定 值中的一个 ; 以及 根据每一个各自的参数的所述两个预定值中的一个, 确定每一个约束方程是否有效。
    12: 如权利要求 11 所述的计算机可读介质, 还包括选择所述网络中每一个节点的质量 平衡方程。
    13: 如权利要求 11 所述的计算机可读介质, 其中, 具有所述至少一个约束方程的所述 网络连接的所述基础方程是液压方程, 并且所述至少一个约束方程是速率约束方程。
    14: 如权利要求 11 所述的计算机可读介质, 其中, 所述两个预定值是 0 和 1。
    15: 如权利要求 11 所述的计算机可读介质, 还包括选择所述网络中另一个连接的另一 个约束方程。
    16: 如权利要求 14 所述的计算机可读介质, 其中, 将每一个各自的参数初始设置为 0.5, 并且如果各自的松弛变量的所述解小于或等于零, 将所述每一个各自的参数调整接近 0, 并且如果所述各自的松弛变量的解大于零, 将所述每一个各自的参数调整接近 1。
    17: 如权利要求 14 所述的计算机可读介质, 其中, 如果每一个约束方程的各自的参数 的值等于 1, 所述网络中每一个约束方程是有效的。
    18: 如权利要求 16 所述的计算机可读介质, 其中, 通过将一减去所述参数乘以减小因 子, 将每一个各自的参数调整接近 1, 并且通过将所述参数乘以另一个减小因子, 将每一个 各自的参数调整接近 0。
    19: 如权利要求 17 所述的计算机可读介质, 还包括 : 如果每一个约束方程的各自的参数等于 1, 将每一个约束方程的松弛变量加至具有所 述至少一个约束方程的所述网络连接的所述基础方程。
    20: 如权利要求 12 所述的计算机可读介质, 还包括 : 求解所述网络中每一个有效约束方程 ; 求解所述网络中每一个质量平衡方程 ; 以及 对于没有约束方程的每一个网络连接求解所述网络中每一个基础方程。

    说明书


    使用松弛变量确定网络中的有效约束的系统和方法

        相关申请的交叉引用
         在此要求于 2007 年 11 月 28 日提交的、 第 60/990,823 号美国临时专利申请的优 先权, 并将其说明书以参考的方式并入本申请。
         关于联邦政府资助的研究声明
         不适用。
         技术领域 本发明通常涉及用于确定地面设施网络中的有效约束的系统和方法。 本发明尤其 涉及在系统方程中使用松弛变量和乘数值以消去无关 ( 无效 ) 约束。
         背景技术 在石油和天然气生产领域, 可以将地面设施网络系统表示为节点、 节点之间的连 接、 节点和油藏网格单元 ( 射孔 ) 之间的连接, 以及节点和油库 ( 开采 ) 或源 ( 喷射 ) 之间 的连接的组。图 1 说明了具有约束的示例性的地面设施网络。
         流程网络计算的主要难点是确定哪个可能的约束是有效的。例如, 可能有数以百 计或可能数以千计的以各种方式相互作用的约束。 这最可能使制作简单的试验和误差策略 过分昂贵, 因为每个试验需要分解设施矩阵。此外, 特定的约束可以互不兼容。例如, 在单 一的出油管线内不可以有两个有效的速率约束, 因为这造成奇异矩阵。
         可以将不同种类的速率约束应用到连接, 例如, 地面条件的最大油相速率约束 QOSMAX。同样, 可以将不同种类的最小 / 最大压力约束 (Pmin/Pmax) 应用于节点。速率约束通 常由在连接处的阀施加。 节点处的 Pmin 约束由在节点的下游连接处的阀施加, 并且 Pmax 约束 由在节点的上游连接处的阀施加。 如果没有规定约束 ( 或约束是无效的 ), 可以将液压方程 应用于连接处, 以将流体成分速率和末端节点压力联系起来。
         必须建立方程组, 以求解连接中的流体成分速率和节点处的压力。方程通?;?网络的性质和给定约束, 并可以包括 :
         (a) 节点处的流体成分质量平衡 (mass balance) 方程。
         (b) 连接处的连接方程, 其可以是 (i) 液压方程, (ii) 速率约束方程, (iii) 节点 下游连接处的具有 Pmin 约束的, 或节点上游连接处的具有 Pmax 约束的压力约束方程, 或 (iv) 其它种类的约束, 诸如分离 / 分流约束方程。
         (c) 射孔处的射孔方程。
         将方程 (a) 和 (c) 固定, 同时需要一种算法以为每一个连接选择仅一个可能的连 接方程, 每一个连接有不只一个可能的连接方程。该算法必须决定为给定的连接选择哪一 个可能的连接方程, 需记住不能求解过约束网络系统且每一个连接仅可以有一个连接方 程。 如果网络系统是过约束的, 那么当使用牛顿法求解网络的非线性方程组时, 建立不能在 牛顿步通过线性解法 (linear solver) 求解的系统方程的奇异矩阵。牛顿法是求解非线性 方程组的标准方法中的一种 ; 其使用线性方程组的解, 该解根据非线性方程的泰勒级数展
         开而得到以估计非线性方程的解。该线性化过程可以在多次迭代中重复, 直至解收敛。将 每一次迭代称作牛顿步。
         图 2 说明了过约束的地面 - 设施网络的示例性部分。在该示例中, 假设流体只有 一种成分, 因此系统的未知数是分别为连接 1 处的速率、 连接 2 处的速率和节点 2 处的压力 的 q1、 q2 和 p2 ; 而 p1( 节点 1 处的压力 ) 和 p3( 节点 3 处的压力 ) 是常数 (p1 = 5000, p3 = 150)。对于该示例还假设, 节点 1 和节点 3 之间的压降足够大, 使得当连接 1 和连接 2 处的 阀完全打开时, 这两个速率约束被破坏。被选作连接 1 和连接 2 处的连接方程的两个速率 约束方程和在节点 2 处的质量平衡方程分别为 :
         q1 = 1000,……………………………………………………… (1)
         q2 = 2000,……………………………………………………… (2)
         q1-q2 = 0, ………………………………………………………(3) 很清楚, 这三个方 程是矛盾的, 所以不能都被满足, 并且节点 2 处的压力是不确定的, 因为该压力没有出现在 任何方程中。如果将未知数排序为 q1、 q2 和 p2, 那么系统方程矩阵变为 :
         该矩阵为奇异矩阵所以不能被求解。
         可以将连接 2 处的 QMAX 约束除去, 因为其值大于连接 1 处的 QMAX 约束?;谎灾?, 速率必须等于或小于两个 QMAX 约束中较小的, 该较小的是连接 1 处的 QMAX = 1000。连接 2 处的 QMAX 约束 ( 方程 (2)) 可以由液压方程替代 :
         -p2+p3+b2q2+c2 = 0,………………………………………… (5) 其中, b2 和 c2 是液 压相关性 (hydraulics correlation) 确定的常数。现在可以求解网络系统, 因为方程系统 的矩阵 :
         不是奇异的。由兰德马克图形公司销售的商业 - 油藏 - 模拟应用 Nexus , 通常以图 3 和图 4 说明的方式通过启发式地 (heuristically) 消去无关约束, 做出完全连结的方式的油藏和 地面设施网络系统 ( 包括油井 ) 的模型, 因此, Nexus 只是与其它传统的确定有效约束的 方法一样昂贵。图 3 说明了 Nexus 时间 - 步模拟的全部工作流程 300。工作流程 300 的 开始于第一个牛顿 - 也称作时间步的开端。在步骤 304 中, 参考图 4 中更为详细的描述的 方式, 使用固定的油藏条件求解独立地面网络模型。在步骤 306 中, 工作流程 300 确定独 立 - 地面 - 网络模型牛顿步是否收敛。如果没有收敛, 则工作流程 300 在步骤 304 重复下 一个牛顿步。如果收敛, 则工作流程 300 在步骤 308 求解连结的油藏 / 地面网络模型, 并在 步骤 310 确定独立 - 地面 - 网络模型牛顿步是否收敛。在步骤 312, 工作流程 300 继续至下 一个时间步。图 4 中的方法 400 说明了求解图 3 中步骤 304 的独立 - 地面 - 网络模型的传统方 法。方法 400 开始于独立 - 网络求解的第一个牛顿, 并通过消去无关的约束启发式地确定 有效约束。在步骤 404 中, 在当前时间步对用户规定的网络配置和约束的任何变化进行检 查。如果没有变化, 那么方法 400 进行至步骤 406 ; 否则, 网络在步骤 404(b) 重新初始化并 且随后方法 400 继续至步骤 406。
         在步骤 406 中, 将每一个当前牛顿步的流体速率和节点压力的值与约束值作比 较, 从而确定每一个约束的破坏因子?;谎灾?, 分析每一个连接和节点, 以确定是否有任何 破坏的约束以及其被破坏的程度。选择每一个连接和节点的具有最大破坏因子的约束, 作 为有效约束的候选 (candidate)。
         在步骤 408 中, 将网络中所有的无关约束移除作为有效约束候选。所有保留的约 束为有效约束。 例如, 在图 1 中, 连接 P-1 和 P-2 的 QOSMAX 约束的总和小于集油 1 的 QOSMAX 约束。如果必须满足的 P-1 和 P-2 的约束被满足, 集油 1 的约束也自动被满足。因此, 集油 1 的约束是无关的, 因而应该被移除。
         然而, 当连接集油 1 被由数以百计或者可能数以千计的连接 ( 其可能本身具有多 种约束 ) 组成的复杂网络从连接 P-1 和连接 P-2 分开, 这种检查逻辑的方式会过分昂贵并 且不可靠。例如, 可能遗漏一些无关的约束, 这导致奇异矩阵, 或者可能移除一些应该被保 留的约束, 这导致不准确的解。 在步骤 410 中, 将约束方程作为具有有效约束的每一个连接的连接方程。否则, 将 液压方程作为连接方程。在步骤 412 中, 使用已知的线性解法求解线性化的网络系统方程。 在步骤 414 中, 线性解法确定网络系统方程是否产生奇异矩阵。如果矩阵是奇异的, 方法 继续至步骤 420 并确定网络是否已被重新初始化。如果网络已经被重新初始化, 那么方法 400 在步骤 422 提前终止。如果网络还没有被重新初始化, 那么方法 400 在步骤 404(b) 重 复并重新初始化网络。然而, 如果矩阵不是奇异的, 那么方法 400 确定在步骤 416 是否发生 收敛。如果没有发生收敛, 那么方法 400 在步骤 406 重复, 这是下一个牛顿步。如果获得了 收敛, 那么方法 400 继续至步骤 418 的求解连结的油藏 / 地面网络模型求解。
         当使用方法 400 时, 使用通过步骤 304 确定的有效约束组 ( 具有固定的油藏条 件 ), 以建立网络方程。 然而, 如果油藏条件改变, 这会导致有效约束组与步骤 304 的约束组 不同, 这意味着网络方程不正确。这迫使进行一个或更多的牛顿步。
         用于确定地面设施网络的有效约束的其它已知和公认的方法可能会一样昂贵, 不 可靠和 / 或效率低, 例如, 由 R.Fletcher 发表的 “Practical Methods of Optimization” , 第二版, Linear Programming 中所描述的优化方法, 当应用约束时, 其需要待计算的一定目 标函数的最大值。
         因此, 需要比传统手段更加有效并可靠的系统化方法, 以确定地面 - 设施网络中 的有效约束。
         发明内容 本发明通过提供一种用于确定网络中的有效约束的系统和方法, 克服了现有技术 中的一个或更多的缺点。
         在一个实施例中, 本发明包括一种用于确定网络中的有效约束方程的方法, 其包
         括: i) 选择至少一个所述网络中的连接的约束方程, 每一个约束方程包括各自的松弛变量 和各自的松弛变量乘数 ; ii) 选择所述网络中的每一个连接的基础方程, 具有所述至少一 个约束方程的所述网络连接的所述基础方程包括所述每个约束方程的各自的松弛变量和 另一个各自的松弛变量乘数 ; iii) 求解每一个各自的松弛变量 ; iv) 根据每个各自的松弛 变量的解, 对于每个各自的松弛变量乘数和每个各自的另一个松弛变量乘数, 在两个预定 值之间调整各自的参数 ; v) 重复求解每个各自的松弛变量以及对于每个各自的松弛变量 乘数和每个各自的另一个松弛变量乘数调整各自的参数的步骤, 直至每个各自的参数等于 所述两个预定值中的一个 ; 以及 vi) 根据每个各自的参数的所述两个预定值中的一个, 确 定每个约束方程是否有效。
         在另一个实施例中, 本发明包括一种用于确定网络中的有效约束方程的具有计算 机可执行指令的计算机可读介质。指令是可执行的, 以实现 i) 选择至少一个所述网络的连 接的方程, 每一个约束方程包括各自的松弛变量和各自的松弛变量乘数 ; ii) 选择所述网 络的每一个连接的基础方程, 具有所述至少一个约束方程的所述网络连接的基础方程包括 所述每个约束方程的各自的松弛变量和另一个各自的松弛变量乘数 ; iii) 求解每一个各 自的松弛变量 ; iv) 根据每个各自的松弛变量的解, 对于每个各自的松弛变量乘数和每个 各自的另一个松弛变量乘数, 在两个预定值之间调整各自的参数 ; v) 重复求解每个各自的 松弛变量以及对于每个各自的松弛变量乘数和每个各自的另一个松弛变量乘数调整各自 参数的步骤, 直至每个各自的参数等于所述两个预定值中的一个 ; 以及 vi) 根据每个各自 的参数的所述两个预定值中的一个, 确定每个约束方程是否有效。 从下面的各种实施例和相关附图的描述, 本发明的其它方面、 优点和实施例对本 领域的技术人员将会变得显而易见。
         附图说明
         下面将参考附图描述本发明, 其中, 相同的部分使用相同的参考标号, 其中 : 图 1 说明了示例性的具有约束的地面设施网络。 图 2 说明了过约束地面 - 设施网络的示例性部分。 图 3 为说明 Nexus 工作流程的流程图。 图 4 为说明传统的用于确定地面 - 设施网络中的有效约束的方法的流程图。 图 5 为说明用于实现本发明并确定地面 - 设施网络中的有效约束的方法的流程 图 6 为用于说明根据本发明的有效约束的确定的示例性的地面 - 设施网络。 图 7 为说明用于实施本发明的系统的方框图。图。
         具体实施方式
         下面具体地描述了本发明内容的主题, 但描述本身并不限制本发明的范围。 因而, 结合现有或未来的技术, 本发明的主题也可以以其它方式实施, 以包括不同的步骤或类似 于在此描述的步骤的结合。此外, 虽然在此可以使用术语 “步骤” 以描述采用的方法的不同 组成部分, 但不应该将术语解释为暗示在此公开的各种步骤之间的任何特定的顺序, 除非 通过特定顺序的描述明确限制。松弛变量提供一种方式确保约束方程总是被满足。将松弛变量加入约束方程。然 后, 如果约束是无效的, 松弛变量继续松弛 - 也就是, 其变为不能用其满足原始写入的无效 约束方程的量。
         通常在一组不等式约束方程中使用松弛变量, 其中一些是有效的。 然而, 在本发明 中方程典型地成对出现。 存在基础方程, 如果约束无效, 应用该基础方程, 还存在约束方程, 如果约束是有效的, 应用该约束方程。这两个方程中有一个并且只有一个是有效的。因此 将松弛变量加至无效的方程而不是有效的方程。为此, 定义了两个乘数。将第一个乘数乘 以松弛变量加至基础方程, 并且第二个乘数乘以松弛变量加至约束方程。下面描述的过程 用于确定两个乘数。逐渐地, 过程导致一个乘数为零而另一个为一。乘数为零的方程不包 含松弛变量, 因此是有效的, 乘数为一的方程包含松弛变量, 因此是无效的。乘数作为参数 (w) 的函数计算。设置第一个乘数等于 (w), 并且设置第二个乘数等于 (1-w)??梢允褂闷?它函数的形式, 例如, 以及 计算乘数, 然而, 在此使用 (w) 和 (1-w) 以说明各种示例。那么, 如果 w = 1, 将松弛变量加入基础方程, 约束有效。如果 w = 0, 则将松弛变 量加入约束方程, 约束为无效的。下面的方法论更详细地描述了该过程。
         速率约束
         对于两种成分, 线性化的液压方程为 :
         -pi+po+a1q1+a2q2 = d …………………………………………… (7)
         并且线性化的速率约束为 :
         b1q1+b2q2 = qmax ………………………………………………… (8)
         在给定的连接处, 必须应用这些方程中的一个或另一个。将乘以乘数的松弛变量 加入液压方程。方程 (7) 变为 :
         -pi+po+a1q1+a2q2+wcsc = d……………………………………… (9)
         其中 wc 是乘数, sc 是松弛变量, 并且 c 表示特定的约束。 将乘以另一个乘数 (1-wc) 的相同的松弛变量加入约束方程。方程 (8) 变为 :
         -b1q1-b2q2+(1-wc)sc = -qmax……………………………………… (10)
         方程中的符号发生变化, 从而能够将松弛变量的贡献加入而不是将松弛变量的贡 献减去, 同时也构建了方程, 以这样的方式, 使负松弛变量是有效的而正松弛变量是无效 的。 需要 0 ≤ wc ≤ 1 并且对于解而言, wc 等于零或者一。松弛变量出现在无效方程中 而不出现在有效方程中。因此 wc = 1 表明速率约束是有效的, 而 wc = 0 表明液压方程是有 效的。
         就液压方程而言, 计算的压降可以大于没有松弛变量计算的压降, 但不能小于。 因 此, sc 的正值在该方程中有效, 但负值无效。所以, 如果 sc ≤ 0, 优选 wc = 0, 使液压方程有 效而速率约束无效。现在考虑速率约束, 计算的速率可以小于但不能大于规定。sc 的负值 在该方程中可以有效, 但正值无效。如果 sc > 0, 优选 wc = 1, 使速率约束方程有效而液压 方程无效。
         如前面所述, 可以将每一个约束的 wc 按照如下确定???, 设置所有的 wC = 0.5, 并求解网络。随后调整 wc。如果 sc ≤ 0, 通过减小因子 (R) 减小 wc :
         如果 sc > 0, 那么增加 wC, 通过同样的减小因子 (1-R) 使数量 (1-wC) 减小 :可以使用其它减小因子, 例如, 可以使用随着 wC 接近零或一大小变化的其它减小 因子, 但是在此使用常数 R, 以说明各种示例。 同样, 当 s 大于其小于或等于零时所使用的较 小因子时, 可以使用不同的减小因子, 但在此使用相同的 R。 该过程继续至如果 sc ≤ 0, wCk+1 小于公差, 或如果 sc > 0, 小于公差。计算
         一般使用 R = 0.01 以及公差= 1.0e-8。采用一个另外的步骤??梢苑⑸氖?: 当将其 w 设 置为一后, 松弛变量变为负的, 当将其参数 w 设置为零后, 松弛变量变为正的。当此发生时, 将参数 w 重新设置为 0.5, 并继续进行迭代。
         优选地, 需要包括 sc 而不是 wC 的方程。为此, 将方程 (10) 加入方程 (9) 以获得结 合的速率约束 - 液压方程 :
         -pi+po+(a1-b1)q1+(a2-b2)q2+sc = d-qmax…………………… (13)
         可以通过仅将方程 (9) 和方程 (10) 相加就可以获得不包含 wC 的方程, 因为两个 方程中的松弛变量乘数分别是 wC 和 (1-wC), 并且这两个乘数的和是一。本领域的技术人员 将认识到, 如果这两个乘数的总和不是一, 可以以其它的方式结合这两个方程, 以获得相同 的结果。
         压力约束
         压力约束方程是 :
         pi = pmin ………………………………………………………… (14)
         使用类似上面描述的用于速率约束的过程, 建议应该将该约束从液压方程 ( 方程 (7)) 中减去。但是, pi 将被消去, 这可能会在后面的解中产生问题。因此, 在引入松弛变量 之前将该方程分为两个。结果是 :
         将公式 (9) 加入公式 (3) 产生包含 sc 但不包含 wC 的方程 :分离器流出
         分离器至少有两个流出流。它们中的一个具有压力连续性, 这会导致下面的方程
         pi-po, j = 0……………………………………………………… (17)
         在此, i 表示分离器的入口节点, 而 o, j 表示第 j 个分离器出口节点。将第一个 (j = 1) 出口方程作为基础方程, 并且另一个出口方程作为约束。如果有三个出口, 三个方程 如下 :
         -pi-po, 1+w2s2+w3s3 = 0…………………………………… (18)
         pi-po, 2+(1-w2)s2 = 0……………………………………… (19)
         pi-po, 3+(1-w3)s3 = 0……………………………………… (20)
         在此, s 和 w 的下标 2 和 3 分别表示分别与流出 2 和 3 相关的松弛变量。将方程
         (18) 至 (20) 相加得到 :
         pi+po, 1-po, 2-po, 3+s2+s3 = 0 ………………………… (21)
         现在将处理导致具有 s 但没有 w 的方程的操作的阶段。 一旦获得因子, 分解网络矩 阵比对其求解成本要高得多。如果每次 w 变化, 必须重构网络矩阵, 这发生在每次 w 迭代, 则成本过高。该成本的大部分可以通过产生包括 s 但不包括 w 的方程 ( 方程 (13)、 (16) 和 (21)) 而避免。将未知数排序, 使 s 在最后, 并将方程排序使约束方程 ( 方程 (10)、 (15)、 (19) 和 (20)) 在最后。结果是如下形式的分块矩阵 :
         未知向量 x 包括连接速率和节点压力 ; 除了松弛变量这些是所有的未知数。未知 向量 s 包括所有松弛变量。矩阵 Axx 和 Axs 定义节点质量平衡方程、 液压方程、 和结合的约 束 - 液压方程, 诸如方程 (13)。矩阵 Asx 和 Ass(w) 定义约束方程。矩阵 Axx, Axs, 和 Asx 不包括 任何乘数, 因此不依赖于 w。只有 Ass(w) 依赖 w。这使得可能分解 Axx, 随后整个矩阵如下 :
         进一步求解的第一步是计算 : 然后更新 rs : 留下具有 s 为未知数的矩阵方程为 :然后对矩阵方程求解 s。必要时, 使用方程 (11) 和 (12) 更新 w, 重新计算 Ass(w), 并重新计算方程 (26) 左侧的矩阵。再次求解 s。该过程迭代重复直至所有的 w 为零或一。 当 w 已经收敛时, 返回求解 x?;揪卣蠓匠涛?:
         将其简化为 :建立上面显示的方程, 使得 wc = 1 表明约束 c( 可以是分离器流出流 ) 是有效的, 使得 wc = 0 表明约束 c 是无效的, 使得 sc ≤ 0 表明 wc 应该朝零移动, 并且使得 sc > 0 表明 wc 应该朝一移动。这使得可能 : 在抽象的形式中处理松弛变量、 乘数和其相应的参数 ( 即, w)。 也就是说, 在不需要知晓其相关的哪些约束或约束的种类的情况下, 可以求解约束变量
         并调整参数。
         因此, 本发明的方法使得能够确定有效约束而不必且接整个网络。
         可替代地, 如果约束的状态是明显的, 可以仅将其设置为有效状态, 而不是对其用 松弛变量处理。
         方法描述
         现在参考图 2, 可以将本发明应用于示例性的使用松弛变量的过约束地面 - 设施 网络, 以确保总是满足正确的约束方程。将乘以乘数的松弛变量加入约束方程。那么如果 约束是无效的, 乘数为一, 并且乘以一的松弛变量继续松弛 - 也就是, 其变为不能满足的无 效约束的量。另一方面, 如果约束是有效的, 那么乘数为零, 并且将松弛变量乘以零。结果, 满足了最初写入的约束方程。将松弛变量 s1 和 s2 分别加入液压方程以及连接 1 和连接 2 的约束方程。松弛变量 (s1 和 s2) 分别与参数 w1 和 w2 相关。把这些原理应用于图 2 中的地 面 - 设施网络, 系统方程变为 :
         连接 1 : -p1+p2+b1q1+c1+w1s1 = 0,………………………… (29)
         -q1+(1-w1)s1 = -1000,……………………………………… (30)
         连接 2 : -p2+p3+b2q2+c2+w2s2 = 0,………………………… (31) -q2+(1-w2)s2 = -2000,……………………………………… (32)
         节点 2 : q1-q2 = 0,…………………………………………… (33)
         其中, 方程 (29) 和 (30) 分别是在连接 1 的液压方程 ( 被称为基础方程 ) 和速率 约束方程 ( 二者都具有松弛变量 s1)。方程 (31) 和 (32) 分别是在连接 2 的液压方程和速 率约束方程 ( 二者都具有松弛变量 s2)。方程 (33) 是在节点 2 的质量平衡方程。未知数组 为 q1、 q2、 p2、 s1、 和 s2。
         w1 和 w2 的初始值是 0 和 1 之间的某值, 优选 0.5。一旦将这些参数是规定的, 计算 方程 (29)-(33) 中的系数并且求解所产生的方程组。如果松弛变量的解是负数, 那么使用 方程 (11) 将其相关参数 (w) 调整为接近 0。如果松弛变量的解是正数, 那么使用方程 (12) 将其相关参数 (w) 调整为接近 1。 可以将方程组的系数重新计算, 然后可以对产生的方程求 解, 并且可以再次调整参数??梢越贸绦蛑馗?, 直至所有参数 (w1 和 w2) 为 0 或 1。该条 件等于是说所有的乘数为 0 或 1。将方程 (29) 和 (30) 以及方程 (31) 和 (32) 相加, 产生了 不包含 w 的结合的约束 - 液压方程。那么, 根据方程 (22) 到 (28) 所描述的, 在过程开始时 不能消去所有网络的未知数 ( 速率和压力 ), 这产生了松弛变量的子系统 ; 仅松弛变量的子 系统需要以这种方式重复被求解, 因为仅该子系统的方程受参数的影响。所有的 w 为零或 一后, 松弛变量为已知, 并且网络的未知数可以通过回代法求解。当 w1 或 w2 的值为 0, 相应 的连接的约束是无效的并且液压方程完全被满足。否则, w 为一并且约束是有效的。通过 在两个预定值 ( 例如 0 到 1) 之间向上或向下逐渐地调整参数, 直至每一个等于两个值中的 一个, 因此辨别出了有效的 ( 即, 正确的 ) 约束。
         使用松弛变量使用户能够将未知数与相同数量的方程匹配, 并且不需要计算传统 优化方法所需的目标函数的最大值。
         现在参考图 5, 说明了实施本发明并使用松弛变量确定在地面 - 设施网络中的有 效约束的方法 500。因此, 方法 500 可以用于以比图 4 中说明的传统方法 400 效率更高并更 准确的方式, 在图 3 的步骤 304 求解独立 - 地面 - 网络模型。方法 500 开始于独立网络求
         解的第一个牛顿。在步骤 504 中, 对于在当前时间步用户规定的网络配置和约束的任何变 化进行检查。如果没有变化, 那么方法 500 继续到步骤 506 ; 否则, 网络在步骤 504(b) 重新 初始化并且随后方法 500 继续至步骤 506。
         在步骤 506, 通过写入每一个连接的液压方程建立网络系统方程。 在每一个可能被 约束的连接处, 有基础方程。通过写入可能被约束的连接的所有约束方程完成网络系统方 程建立。每一个约束与松弛变量和乘数相关, 同时乘数根据大于零并小于 1 的预定的参数 (w) 值计算, 参数 (w) 值优选 0.5。在步骤 508 中, 网络系统方程根据给定的参数 (w) 值求 解。例如, 如果松弛变量的解是负数, 将其相关的 w 调整接近 0。但是, 如果松弛变量的解是 正数, 将其相关的 w 调整接近 1。在步骤 510, 再次检查所有 w, 以确定其是否为 0 或 1。如 果任何 w 不是 0 也不是 1, 方法 500 继续在步骤 508 重复。但是, 如果所有 w 是 0 或 1, 那么 方法 500 继续至确定收敛的步骤 512。如果没有获得收敛, 方法 500 在步骤 506 重复。如果 获得收敛, 方法 500 继续至步骤 514, 在步骤 514, 求解连结的油藏 / 地面网络模型。
         现在参考图 6, 说明了根据方法 500 确定有效约束的示例性地面设施网络。在图 6 中, 三个生产井 ( 连接 1, 连接 2 和连接 3) 连接至集油中心 ( 节点 4), 该集油中心最终连接 至油库。此外, 在连接 1、 连接 2、 连接 3 和连接 4 有四个 QOSMAX( 地面条件的最大油相速率 约束 ), 以及在节点 3 和节点 5 有两个 Pmin( 最小压力 ) 约束。总共有六个约束。但是, 因为 在节点 5 的 Pmin 约束作为边界条件工作, 并且其必须被满足, 所以余下的五个约束只需要五 个松弛变量, 其每一个与 w 相关。
         网络中的流体流动用黑油模型制作模型, 黑油包括石油、 天然气和水的成分。 在网 络求解过程中, 将油藏条件假设为常数。
         表 1 示出了每一个约束的松弛变量索引。
         松弛变量索引 1 2 3 4 5
         约束 连接 1 的 QOSMAX 约束 连接 2 的 QOSMAX 约束 连接 3 的 QOSMAX 约束 节点 3 的 Pmin 约束 节点 4 的 QOSMAX 约束表 1. 每一个约束的松弛变量索引 表 2 示出了图 5 的步骤 508 和 510 说明的所有迭代的 w 值和相应的松弛变量 (s)的解。
         表 2. 网络牛顿步的所有线性解法迭代的
         参数 (w) 值和松弛变量 (s) 的解
         如图 5 中的方法 500 说明的, 每一个 w 的初始值是 0.5, 并且在每一次迭代之后, 根 据每一个相应的松弛变量的解的值调整每一个 w。在图 5( 步骤 510) 中, 在每一次迭代的 结束检查 w 值, 并且如果 w 不是 0 也不是 1, 方法 500 在步骤 508 重复。在最后一次 ( 第 7 次 ) 迭代, 所有 w 值为 0 或 1。求解松弛变量以及网络未知数, 并且完成网络求解的该牛顿 步。随后在步骤 512 进行收敛检查。如果解仍没有收敛, 继续下一个牛顿步并且方法在步 骤 506 重复直至获得收敛。否则, 方法继续至图 5 的步骤 514 的连结油藏 / 网络求解。根据 本发明的松弛变量的使用也可以用于在图 5 的步骤 514 求解连结的油藏 / 地面网络模型。
         根据解, 在连接 1、 连接 2、 和连接 4 的约束的 w 的值为 1。因此, 这些约束是有效 的。在连接 3 和节点 3 的约束的 w 的值为 0, 这意味着连接 3 和连接 5 由其液压方程控制。 如表 3 部分列出的, 检查网络的相位速率和节点压力的解, 确认已经恰当地满足了 ( 没有被 破坏 ) 所有约束。
         连接 1 的油相速率 连接 2 的的油相速率 连接 3 的油相速率 节点 3 的油相速率 连接 4 的油相速率
         2000STB/ 天 1000STB / 天 3000STB/ 天 2204.203psi 6000STB/ 天 w=1 w=1 w=0 w=0 w=1 约束起作用 约束起作用 约束不起作用 约束不起作用 约束起作用表 3. 网络的相位速率和节点压力的解和相应的约束状态 因此, 图 6 中的示例明确了松弛变量方法为确定地面 - 设施网络中的有效约束提供了较传统方法的各种优点, 其包括但不限于更高的可靠性、 有效性高、 准确性和仿真速 度。
         系统描述
         本发明可以通过计算机可执行程序指令实现, 诸如程序???, 一般称作由计算机 执行的软件应用或应用程序。 软件可以包括, 例如, 执行特定任务或实现特定抽象数据类型 的例行程序、 程序、 对象、 组件和数据结构。软件形成接口使计算机根据输入源而反应。兰 德马克图形公司 (Landmark Graphics Corporation) 销售的商业软件 Nexus 可以作为接 口应用, 以实现本发明。软件结合接收数据的源, 也可以与其它代码段合作, 以开始响应于 所接收的数据的多种任务??梢越砑娲⒃谌我恢执娲⒚教迳?, 诸如 CD-ROM、 磁盘、 磁泡 存储器和半导体存储器 ( 例如, 各种类型的 RAM 或 ROM) 上。而且, 可以通过诸如光纤、 金属 线、 自由空间的多种载体媒体和 / 或诸如互联网的多种网络传输软件及其结果。
         此外, 本领域的技术人员将认识到本发明可以用多种计算机 - 系统配置实施, 包 括手持装置、 多处理器系统、 基于微处理器的或可编程消费类电子产品、 微型计算机、 大型 计算机等任意数量的计算机 - 系统和计算机网络都可以使用本发明??梢栽诜植际郊扑慊?境中实施本发明, 在分布式计算环境中, 任务可以由通过通信网络链接的远程处理装置执 行。在分布式计算环境中, 程序??榭梢晕挥诒镜睾驮冻碳扑慊?- 存储媒体中, 包括记忆存 储装置。 因此, 可以连同各种硬件、 软件或其结合在计算机系统或其它处理系统中实施本发 明。
         现在参考图 7, 说明了在计算机上实施本发明的系统的方框图。 该系统包括计算单 元, 该计算单元有时称作计算系统, 包括存储器, 应用程序, 用户接口和处理单元。 计算单元 只是一个适合的计算环境的示例并不限制本发明的用途或功能的范围。
         存储器主要存储由计算单元执行以实现在此描述的并在图 5 中说明的本发明的 应用程序, 也可描述为包括计算机可执行指令的程序???。因此, 存储器包括实现图 5 中说 明的方法的油藏模拟??楹?Nexus 。
         虽然所示计算单元具有广义存储器, 该计算单元一般包括多种计算机可读媒体。 通过举例的方式但不限制, 计算机可读媒体可以包括计算机存储媒体和通信媒体。计算系 统存储器可以包括诸如只读存储器 (ROM) 和随机存储器 (RAM) 的形式为易失性存储器和 / 或非易失性存储器的计算机存储媒体?;臼淙?/ 输出系统 (BIOS) 一般存储在 ROM 中, BIOS 包括帮助在计算单元内的组成部分 (element) 之间传送信息的基本例程, 诸如在启动 过程中。RAM 一般包括数据和 / 或程序???, 其可以由处理单元立即访问和 / 或即将开动。 通过举例的方式但不限制, 计算单元包括操作系统、 应用程序、 其它程序??楹统绦蚴?。
         存储器中示出的组件也可以包括在其它可擦除 / 不可擦除、 易失性 / 非易失性计 算机存储媒体中。仅例如, 硬盘驱动器可以读不可擦除、 非易失磁性媒体, 或写不可擦除、 非易失磁性媒体, 磁盘驱动器可以读可擦除、 非易失磁盘中, 或写可擦除、 非易失磁盘, 以及 光盘驱动器可以读诸如 CD ROM 或其它光学媒体的可擦除、 非易失光盘, 或写可擦除、 非易失 光盘??梢栽谑纠圆僮骰肪持惺褂玫钠渌刹脸?/ 不可擦除, 易失性 / 非易失性计算机 存储媒体可以包括但不限于盒式磁带、 闪存卡、 数字多用盘、 数字视频磁带、 固态 RAM、 固态 ROM 等等。上面讨论的驱动器和其相关的计算机存储媒体提供了计算单元的计算机可读指 令、 数据结构、 程序??楹推渌莸拇娲?。用户可以通过用户接口将命令和信息输入计算单元, 该用户接口可以是输入装 置, 诸如键盘和指点装置, 其通常称作鼠标、 轨迹球或触控板。 输入装置可以包括麦克风、 控 制杆、 卫星天线、 扫描仪等等。
         通常这些和其它输入装置通过用户接口连接至处理单元, 用户接口连结至系统总 线, 但是这些和其它输入装置也可以通过诸如并行端口或通用串行总线 (USB) 的其它接口 和总线结构连接。监视器或其它类型的显示装置可以通过诸如视频接口的接口连接至系 统。除了监视器以外, 计算机也可以包括其它外围输出装置, 诸如扬声器和打印机, 其可以 通过输出外围接口连接。
         虽然没有示出计算单元的很多其它内部组件, 本领域的普通技术人员应认识到这 些组件和互连是众所周知的。
         松弛变量方法是系统的和自动的 - 意味着, 一旦系统方程已经建立, 对于网络中 的任何结构, 系统可以被求解。相反, 传统方法的实施依赖于逻辑检测网络中不同结构和 约束的不同组合的程度如何, 这可能遗漏或移除一些无关约束并造成仿真的提前终止。此 外, 如果增加了新类型的约束, 松弛变量方法应该工作, 而不需要新的代码以检查特定的约 束是否是无关的。 因此, 本发明 (1) 通过消去处理每个可能的条件 ( 网络配置和约束类型 ) 的需要, 克服了未被辨别的情况, (2) 避免了修订情况的系统描述的需要, 以及 (3) 如果连 结的油藏 / 网络求解中的网络方程也使用松弛变量方法, 提供了更快的求解方法, 因为传 统方法可以确定独立的网络求解和连结的油藏 / 网络求解之间的不同组有效约束, 这造成 该连结的油藏 / 网络求解的至少又一次牛顿步。 虽然已经结合优选实施例描述了本发明, 本领域的技术人员将理解本发明不限于 这些实施例。 因此, 在不脱离所附权利要求及其等同的精神和范围的情况下, 可以对公开的 实施例作出各种替代实施例和修改。例如, 可以将本发明用于石油和天然气工业之外的应 用, 以确定其它类型网络中的有效约束。
        

    关于本文
    本文标题:使用松弛变量确定网络中的有效约束的系统和方法.pdf
    链接地址://www.4mum.com.cn/p-5867899.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 广东十一选五 山东十一选五吧 幸运七星彩计划软件 365投注抵用金 小霸王越野机车app 11选五中奖规则 太子中心八肖 七星彩基本走势图带连 国彩网官方网站 VR三分彩是官网开奖吗 北京pk10走势图软件 全民奇迹赚钱技巧集锦 双色球开奖结果双今天 时时彩老手心得 江苏快三计划软件免费下载 11选5软件