三参数非线性方程自动拟合方法的研究及应用.pdfVIP

石油天然气学报(江汉石油学院学报) 2009年4月第31卷第2期 of andGas V01．31No．2 JournalOil Technology(J．JPI)Apr．2009 三参数非线性方程自动拟合方法的 研究及应用， 邹存友。常毓文，唐玮，王东辉，胡丹丹 (中围石油勘探开发研究院，北京1ooos3) [摘要]针对油藏工程中广泛存在的三参数非线性方程的求解问题．提出了一种简单有效的自动拟合方 法。其基本思路是：首先，将三参数非线性方程进行归一化处理，并确定待定参数的初始分布区闻；其 次，以待定参数最优解使非线性方程的线性化程度最高为准则，在待定参数分布区间内，通过等步长进 行迭代、对比以不断缩小参数分布范围．在计算精度的控制下，最终获得最优解。将上述算法开发成计 算机分析软件，并应用软件对我国和世界范围内的大量油田资料进行了分析。经实际验证，该算法及应 用软件具有计算速度快、精度高、易操作的优点。而且不存在收敛问题。据初步统计，将待定参数分布 区间分成100等份．迭代2～7次即可完全满足精度要求。最后，通过三个油田实例说明了这种方法的应 用效果。 [关键词]三参数；非线性方程；分布区问；迭代；自动拟合 [中图分类号]TE328 [文献标识码]A 线性回归是油藏工程中广泛使用的方法，如用于求得油田可采储量的水驱曲线系列，以甲型水驱曲 线法r-3为例，须对(N，，lgW，)进行线性回归，当得到直线的截距和斜率后，再利用截距和斜率值进一步 求得油田(或区块)的可采储量以及利用童宪章或陈元千公式求得油田的动态地质储量等。再如，Arps[2’3] 的指数递减分析需利用(￡，IgQ)进行线性回归求得递减率。这样的例子不胜枚举，其特征是用于回归的 数据点是确定的，方程中的待定参数只有2个，因而只需用实际点线性回归就可以达到正确分析的目 的，这里不再赘述。 然而，在众多的油藏工程方法中，存在这样一种类型，即控制方程中有3个待定参数，比如绝大部 性方程。对于此类问题，目前广泛使用的是所谓的“试凑法”，这样做无疑加大了工作量，而且也无法 满足精度要求。笔者就此类问题的自动拟合算法进行了研究，以最小二乘法为基础，通过逐次叠代缩小 步长逼近真实解的过程，很好地解决了此类问题的快速求解。文中以广义翁氏模型、HC模型和Arps 的双曲递减分析为例，详细地阐述了该方法的实现原理和过程。最后，以我国大庆油区的喇萨杏油田、 俄罗斯著名的巴夫雷油田和我国的M油田为例，说明了此方法的应用效果。 1问题提出 在油藏工程方法中广泛存在的三参数非线性方程的求解，目前大多数采用试凑的方法实现。这无疑 加大了操作时间和精力，并且不能保证计算精度。下面以广义翁氏模型、HC模型和Arps双曲递减分 析为例进行说明。 [收稿日期]2008一12—16 [作者筒介]邹存友(1978一)，男，2002年大学毕业，博士(后)。现主要从事油气田开发规划、试井分析和数值模拟等方面的研 究工作。 石油天然气学报(江汉石油学院学报) 2009年4月 理论推导，该预测模型的关系式如下： Q—at6e一‘‘7‘’ (1) Qm。=a(bc／2．718)6 (2) t。一bc (3) NR—acb-l-1r(6+1) (4) 发生的时间，a；F为伽玛函数；a、6、c为模型常数。 HC模型嘲是胡建国和陈元千根据大量实际开发资料推导建立的产量预测模型。其基本关系式为： NP一寿 (5) Q一赫 (6) Qm。一盟型等