油藏数值模拟自动历史拟合方法研究及应用.pdfVIP

石 油 勘 探 与 开 发 82 2003年 8月 PETR0IEUM EXPL0RAT10N AND DEVEL0PMENT Vo1．30 No．4 文章编号 ：1000—0747(2003)04—0082—03 油藏数值模拟 自动历史拟合方法研究及应用 陈兆芳 ，张建荣 ，陈月明 ，李卉朋。 (1．中国石化胜利油 田现河采油厂 ；2．石油大学(华东)；3．中国石油大港油田第八采油作业 区) 摘要 ：油藏数值模拟 的历史拟合过程 中模型参数多、可调整 自由度大，目前普遍采用实验校正法进行拟合 ，主观随意性大 且费时多。引入最优化理论建立 目标函数，以渗流方程和油藏 的边界条件 、初始条件为约束条件 ，将历史拟合转化为非线 性规划中求解定义在 维空间的函数极值 问题 ，用改进 的共轭梯度法求解 ，基本解决 了收敛 问题。拟合过程 中调整参数 的原则为 ：不确定优先，敏感优先，先全局后个别 。为减少调整参数，用解析函数表达相对渗透率与饱和度 的关系，将模拟 区域按油藏地质和井网条件划分为若干小区，分区拟合。对某封 闭小断块油藏的区块与单井综合含水和地层压力进行 自 动历史拟合 ，调整 的主要参数是相对渗透率和渗透率 ，计算了 126次，拟合结果与实际动态相符 。图2参8 关键词 ：油藏数值模拟 ；自动历史拟合 ；最优化 ；目标函数 ；约束条件 中图分类号 ：TE319 文献标识码 ：A 在油藏数值模拟 中，目前 国内外仍普遍采用实验 校正法进行历史拟合 ，即靠人工反复试算 ，根据计算者 )孥一V(b(x)Vp)一 z ∽ 经验来修改油层参数 。 自动历史拟合方法从 20世纪 6(z)ap — r(S)(户一 Po) 60年代已开始进行研究[I-73，但一直处于探索阶段。 Ull p(x，0)= Pi(z) 1自动历史拟合模型及其计算方法 S (z，0)： S (z) 地层识别 中的最佳拟合包括 3方面的识别 ：参数 其中 一 ； (如渗透率 、厚度 、孔隙度等)，边界 (如边界位置、边界 物性 、边水大小及能量等等)，过程特性 (如相对渗透率 口一 C S b— Kh／p 与饱和度的变化关系等)。这是 自动历史拟合要解决 这样，拟合就转化成非线性规划 中求一个定义在多 的问题 ，即找出实际油藏动态变化与同一拟合期 中模 维空间的函数极值问题 ，自变量应满足用偏微分方程表 型的油藏动态变化 的差异 ，以降低此差异为 目标 ，利用 达的渗流方程、边界条件和初始条件等约束条件 。 最优化方法修正模型参数及结构n]。 1．2问题的求解 1．1模型的建立 1．2．1一维线搜索方法 假设需拟合 的参数为压力和综合含水率 ，可 以建 在 自动历史拟合 中，线搜索就是在其它参数不变 立如下 目标函数 ： 的条件下，寻找某一参数的最佳值 ，即求解或者非精确 地求解多维函数 ，(z)在一个一维子空间的极小值 ： J一∑W~(poh一户)+∑w (瑞，一舟，) minf(x+ad) l一 1 J一 1 其中 户 一 ，l(z1，z2，…，X) 其 中 z∈R d∈R dV 厂(z) 0 儡 ．一，2(z1，X2，…，X