技术领域本发明属于油藏生产过程优化领域，特别是涉及一种基于混合遗传算法的油藏生产过程设定点优化方法。
背景技术：
包括油藏生产过程设定点优化在内的复杂工业生产过程优化与控制方法研究一直是一个热门领域，蕴藏着极大的经济效益，具有广阔的应用前景，其发展对我国国民经济意义重大。注水是油藏采油过程中最经济地提高采收率的开采方式之一，同时也是我国大部分油藏主要采油方式之一。在油藏注水采油生产过程中，一般在油藏地面的适当位置挖掘有若干个注水井和生产井，每个注水井上设有一套注水控制系统，每个生产井上也设有一套生产控制系统，以保证向各个注水井以适当速率注水并从各个生产井以适当速率抽取油水混合物，从而将油藏中的石油开采出来。过程优化控制理论与方法一直是一个热门领域，如何寻求工艺变量的最优设定值并实时地对系统进行设定以使生产过程最优地运行，是当今国际复杂工业生产过程优化与控制领域的研究热点之一。油藏生产过程长期以来积累了大量的注水井的注入率和生产井的产出率等历史数据，通过这些数据的研究与分析，可以得到有关各注水井和生产井之间的连通信息，实时优化各注水井的注入率和各生产井的产出率，并及时调整各注水井和各生产井对应控制系统的设定点，可使石油开采过程保持在最优工况下进行。然而，很少有文献报道对油藏生产过程中各注入率和产出率的设定点优化，传统的基于梯度信息的优化方法难以对油藏生产过程设定点优化问题进行有效的求解。
技术实现要素：
有鉴于现有技术的上述缺陷，本发明所要解决的技术问题是提供一种更加有效的油藏生产过程设定点优化方法。为实现上述目的，本发明提供了一种基于混合遗传算法的油藏生产过程设定点优化方法，按以下步骤进行：步骤一、设定在油藏地面挖掘m个注水井和n个生产井，设定生产井j的总产出率为qj，m、n均为正整数；建立油藏容量模型：其中ik是注水井k的注水率，τj为生产井j的时间常数，λkj为注水井k与生产井j之间的权重，设定生产井j在t时刻的石油产出率为qcj，建立油藏生产率模型求取aj和bj；步骤二、建立油藏生产过程设定点优化模型：p0当前石油销售价格，pwj为从油水混合液中分离水的单位费用，qwj为生产井j在t时刻的水产出率；步骤三、以油藏容量模型、油藏生产率模型为约束条件，以所述R为目标函数，得到油藏生产过程设定点优化问题：z(x)=Rmaxs.t.qj=Σk=1mλkjik-τjdqjdtqcj=qj/[1+aj(Σk=1mλkj∫t0tfik(s)ds)bj]lk≤ik≤ukj=1,2,...,n,k=1,2,...,m.;]]>步骤四、对于油藏生产过程设定点优化问题，将其离散化为一般的非线性约束优化问题：z′(x)=minf(x→)x→=(x1,x2,...,xd)s.t.hj(x→)=0,j=1,2,...pgj(x→)≥0,j=p+1,p+2,...,mli≤xiui,i=1,2,...,d;]]>步骤五、对z′(x)，设定Lagrange乘子向量λk、和罚参数向量σk、增广Lagrange乘子法第k步所求得子问题为：其中是修改增广Lagrange函数；x为决策变量，λ为拉格朗日乘子参数；σ为罚函数参数；其中，步骤六、由遗传算法对步骤五所述子问题进行全局寻优以得到下一迭代点，然后修正乘子向量和罚参数向量，检查收敛准则是否满足，当收敛准则不满足时重新构造界约束子问题，收敛准则满足时终止算法。较佳的，步骤六所述遗传算法在进化过程中随机组合一个交叉算子和一个变异算子产生若干个新的个体。本发明的有益效果是：本发明能有效地处理复杂的非线性约束优化问题，能够对油藏生产过程设定点更加有效的进行优化。附图说明图1是本发明一具体实施方式的流程示意图。图2是组合交叉算子和变异算子示意图。图3是第一种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。图4是第二种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。图5是第三种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。图6是第四种生产井的实际石油产出率与石油生产率模型计算所得到的石油产出率的比较示意图。图7是第一种生产井的实际总产出率和容量模型计算值的比较示意图。图8是第二种生产井的实际总产出率和容量模型计算值的比较示意图。图9是第三种生产井的实际总产出率和容量模型计算值的比较示意图。图10是第四种生产井的实际总产出率和容量模型计算值的比较示意图。图11是p0＝60，pw＝1，r＝0.0时的注入率最优设定值示意图。图12是p0＝60，pw＝1，r＝0.2时的注入率最优设定值示意图。具体实施方式下面结合附图和实施例对本发明作进一步说明：如图1至图12所示，一种基于混合遗传算法的油藏生产过程设定点优化方法，由两层循环形式组成，内层循环中由改进遗传算法对子问题进行全局寻优以得到下一迭代点，外层循环则修正乘子向量和罚参数向量、检查收敛准则是否满足、重新构造界约束子问题或在收敛准则满足时终止算法，按以下步骤进行：步骤一、设定在油藏地面挖掘m个注水井和n个生产井，设定生产井j的总产出率为qj，m、n均为正整数；建立油藏容量模型：其中ik是注水井k的注水率，τj为生产井j的时间常数，λkj为注水井k与生产井j之间的权重，设定生产井j在t时刻的石油产出率为qcj，建立油藏生产率模型求取aj和bj；步骤二、建立油藏生产过程设定点优化模型：p0当前石油销售价格，pwj为从油水混合液中分离水的单位费用，qwj为生产井j在t时刻的水产出率；步骤三、以油藏容量模型、油藏生产率模型为约束条件，以所述R为目标函数，得到油藏生产过程设定点优化问题：z(x)=Rmaxs.t.qj=Σk=1mλkjik-τjdqjdtqcj=qj/[1+aj(Σk=1mλkj∫t0tfik(s)ds)bj]lk≤ik≤ukj=1,2,...,n,k=1,2,...,m.;]]>步骤四、对于油藏生产过程设定点优化问题，将其离散化为一般的非线性约束优化问题：z′(x)=minf(x→)x→=(x1,x2,...,xd)s.t.hj(x→)=0,j=1,2,...pgj(x→)≥0,j=p+1,p+2,...,mli≤xiui,i=1,2,...,d;]]>步骤五、对z′(x)，设定Lagrange乘子向量λk、和罚参数向量σk、增广Lagrange乘子法第k步所求得子问题为：其中是修改增广Lagrange函数；其中，步骤六、由遗传算法对步骤五所述子问题进行全局寻优以得到下一迭代点，然后修正乘子向量和罚参数向量，检查收敛准则是否满足，当收敛准则不满足时重新构造界约束子问题，收敛准则满足时终止算法。本实施例中，步骤六所述遗传算法在进化过程中随机组合一个交叉算子和一个变异算子产生若干个新的个体。本发明由软件Eclipse对非均质油藏Synfield进行模拟所得到的数据，用来验证本文所提出的油藏生产过程设定点优化模型以及所提出的混合优化方法的有效性。Synfield油藏是一个有5个注水井(分别用I01，I02，I03，I04和I05表示)和4个生产井(分别用P01，P02，P03和P04表示)的多层各向异性油藏，油藏中水平方向的渗透率是竖直方向的3倍。用Eclipse对Synfield模拟时，时间长度为100个月，即用5个注水井的100个月注入率模拟产生4个生产井相应的100个月的产出率数据.将1-50个月的数据获得设定点优化模型中的参数，然后用混合优化方法求解51-100个月时段的油藏生产过程设定点优化模型，得到5个注水井在该时段内每个月的最优注入率设定值，再以这些最优设定值用Eclipse对Synfield进行模拟产生4个生产井对应的新50组产出率数据，最后将新数据与对比数据进行比较，分析油藏生产过程设定值优化前后的经济效益.由油藏容量模型和1-50个月的数据得到了油藏中各个井之间的连通性权重及每个生产井的时间常数，如表1所示.表1由模型和数据所得权重及时间常数从表1可知，对于非均质油藏，对称注水井和生产井之间的权重是不对称的，权重的大小恰当地量化了各注水井与生产井之间连通性大小。由图7至图10可知，容量模型能较好地描述注水油藏中水的衰减和滞后特性，以及各注水井注入率和各生产井产出率之间的相互关系.生产井P01，P02，P03和P04的总产出率和容量模型计算量之间拟合的确定性系数R2分别为0.9959、0.9970、0.9973和0.9991。基于前面所得到的权重λ和模型数据，对油藏生产率模型(3)中参数aj和bj的估计结果如表2.表2由模型数据所得参数aj和bj由图3至图6可知，4个生产井的石油产出率都拟合的很好，由石油生产率模型计算所得到的石油产出率与其实际石油产出率拟合的确定性系数R2分别为0.9629、0.9590、0.9608和0.9698，说明所建立的石油生产率模型较准确地描述了它们的石油产出率，可对它们的石油产出率较准确预测。利用油藏生产过程设定点优化模型求各注水井的注入率和各生产井的产出率的最优设定值，对每个注水井，本实施例设其注入率下限为零，上限分别为其100个月数据中的最大值(RB/天)，即lk＝0，k＝1，2，3，4，5，u1＝3603.739，u2＝2000.666，u3＝1464.216，u4＝1351.601，u5＝2220.153.本实施例优化时间段从第51个月开始到第100个月，以月为步长离散化问题z(x)得到下面非线性约束优化问题：maxΣk=150{poΣj=14qcj(k+50)-pwΣj=14qwj(k+50)