发明名称 基于偏好多目标优化的油田机采过程决策参数优化方法
摘要 本发明提供的基于偏好多目标优化的油田机采过程决策参数优化方法,包括:确定油田机采油过程中的效率影响因素和性能变量;对样本中的载荷变量进行降维处理构建新样本,并归一化新样本;基于归一化后的新样本构建神经网络模型;利用ST‑UPFNN算法估计神经网络模型中由权值阈值所组成的状态变量的最优状态;并利用最优状态变量重构更新后的神经网络模型获得油田机采油过程模型;构建实际产液量的偏好函数;利用多目标进化算法对决策参量各自的上下限进行优化;将优化后的决策变量,带入油田机采油过程模型,计算优化后的决策变量的系统性能的平均值,与实际样本的系统性能的平均值进行比较。利用本发明可以提高油田机采油的生产效率,降低能耗。
申请公布号 CN106502096A 申请公布日期 2017.03.15
申请号 CN201611001501.7 申请日期 2016.11.14
申请人 重庆科技学院 发明人 唐海红;辜小花;杨利平;张堃;李太福;裴仰军;聂玲;王坎
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 重庆蕴博君晟知识产权代理事务所(普通合伙) 50223 代理人 王玉芝;杨明
主权项 一种基于偏好多目标优化的油田机采过程决策参数优化方法,包括:步骤S1:确定油田机采油过程中的效率影响因素,构成效率观测变量集合{x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,L x<sub>n</sub>};以及,选取油田机工艺系统的性能变量,构成性能观测变量集合{y<sub>1</sub>,y<sub>2</sub>};其中,x<sub>1</sub>为冲次决策变量,x<sub>2</sub>为有效冲程决策变量,x<sub>3</sub>~x<sub>5</sub>分别为计算泵效环境变量、含水率环境变量,平均功率因数环境变量,x<sub>6</sub>~x<sub>n</sub>均为载荷环境变量;性能观测变量的个数l=2,y<sub>1</sub>为日产液量,y<sub>2</sub>为日耗电量;步骤S2:根据效率观测变量集合{x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,L x<sub>n</sub>}和性能观测变量集合{y<sub>1</sub>,y<sub>2</sub>},采集通过ST‑UPFNN算法构建神经网络模型的观测变量的样本值矩阵[x<sub>1</sub>,x<sub>2</sub>L x<sub>n</sub>,y<sub>1</sub>,y<sub>2</sub>];其中,设定采样周期为T,在采集观测变量的过程中,如果采样周期小于T,对T周期内的样本求取平均值以作为该T周期的样本[I,Y];如果采样周期大于T,剔除采集到的观测变量;其中,将样本中的I作为输入样本,将样本中的Y作为输出样本;步骤S3:利用主元分析算法对载荷环境变量进行降维,构建新的载荷主元变量{L<sub>z1</sub>,L<sub>z2</sub>,...,L<sub>zd</sub>};其中,构建新的载荷主元变量{L<sub>z1</sub>,L<sub>z2</sub>,...,L<sub>zd</sub>}为d个载荷主元分量,每个载荷主元分量的维度与所述样本[I,Y]的数量相同;步骤S4:重新组合非载荷变量与d个载荷主元分量,构建新的输入样本I<sub>1</sub>,并对新的输入样本I<sub>1</sub>和输出样本Y进行归一化,获得归一化后的样本<img file="FDA0001152436000000011.GIF" wi="146" he="79" />其属于[‑1,1];其中,非载荷变量包括冲次决策变量x<sub>1</sub>、有效冲程决策变量x<sub>2</sub>、计算泵效环境变量x<sub>3</sub>、含水率环境变量x<sub>4</sub>、平均功率因数环境变量x<sub>5</sub>;步骤S5:基于所述归一化后的样本<img file="FDA0001152436000000012.GIF" wi="128" he="86" />构建神经网络模型和所述神经网络模型的初始状态变量X,以及,将所述归一化后的样本<img file="FDA0001152436000000013.GIF" wi="123" he="85" />中的<img file="FDA0001152436000000014.GIF" wi="35" he="63" />作为所述神经网络模型的输入,将所述归一化后的样本<img file="FDA0001152436000000015.GIF" wi="123" he="85" />中的<img file="FDA0001152436000000016.GIF" wi="40" he="63" />作为所述神经网络模型的输出;其中,所述神经网络模型为:<maths num="0001"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>s</mi><mn>2</mn></msub></munderover><mrow><mo>(</mo><mi>f</mi><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>S</mi><mn>1</mn></msub></munderover><msubsup><mi>w</mi><mrow><mi>i</mi><mi>k</mi></mrow><mn>1</mn></msubsup><msub><mi>I</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>b</mi><mi>i</mi><mn>1</mn></msubsup></mrow><mo>)</mo><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>w</mi><mrow><mi>k</mi><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>b</mi><mi>j</mi><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152436000000021.GIF" wi="1398" he="166" /></maths>其中,I<sub>k</sub>为所述训练样本的矢量样本值,并作为所述神经网络模型的输入,<img file="FDA0001152436000000022.GIF" wi="78" he="83" />为网络输入层到隐含层的神经元的连接权值,<img file="FDA0001152436000000023.GIF" wi="54" he="83" />为网络输入层到所述隐含层的神经元的阈值,<img file="FDA0001152436000000024.GIF" wi="77" he="87" />为所述隐含层到网络输出层的神经元的连接权值,<img file="FDA0001152436000000025.GIF" wi="54" he="79" />为所述隐含层到所述网络输出层的神经元的阈值,其中,i=1,2…S<sub>0</sub>;j=1,2…S<sub>1</sub>;k=1,2…S<sub>2</sub>;S<sub>0</sub>为所述网络输入层的神经元的数量,S<sub>1</sub>为所述网络隐含层的神经元的数量,S<sub>2</sub>为所述网络输出层的神经元的数量;所述初始状态变量X为:<maths num="0002"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>w</mi><mn>11</mn><mn>1</mn></msubsup><mi>L</mi></mrow></mtd><mtd><msubsup><mi>w</mi><mrow><msub><mi>s</mi><mn>0</mn></msub><msub><mi>s</mi><mn>1</mn></msub></mrow><mn>1</mn></msubsup></mtd><mtd><mrow><msubsup><mi>b</mi><mn>1</mn><mn>1</mn></msubsup><mi>L</mi></mrow></mtd><mtd><msubsup><mi>b</mi><msub><mi>s</mi><mn>1</mn></msub><mn>1</mn></msubsup></mtd><mtd><mrow><msubsup><mi>w</mi><mn>11</mn><mn>2</mn></msubsup><mi>L</mi></mrow></mtd><mtd><msubsup><mi>w</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><msub><mi>s</mi><mn>2</mn></msub></mrow><mn>2</mn></msubsup></mtd><mtd><mrow><msubsup><mi>b</mi><mn>1</mn><mn>2</mn></msubsup><mi>L</mi></mrow></mtd><mtd><msubsup><mi>b</mi><msub><mi>s</mi><mn>2</mn></msub><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152436000000026.GIF" wi="1356" he="87" /></maths>步骤S6:利用ST‑UPFNN算法估计所述神经网络模型的最优状态变量;步骤S7:将所述最优状态变量作为所述神经网络模型的<img file="FDA0001152436000000027.GIF" wi="344" he="86" />和<img file="FDA0001152436000000028.GIF" wi="86" he="79" />重构神经网络表达式,获得油田机采油过程模型;步骤S8:构建日产液量y<sub>1</sub>的偏好函数perf<sub>c</sub>(y);步骤S9:利用NSGA‑Ⅱ算法对日耗电量y<sub>2</sub>的决策参量各自的上下限进行优化;步骤S10:将优化后的决策变量结合环境变量,带入由ST‑UPFNN算法所建立的油田机采油过程模型,在所述油田机采油过程模型的基础之上进行优化,计算优化后的决策变量的系统性能的平均值,与实际样本的系统性能的平均值进行比较,如果优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值,利用优化后的决策变量对实际生产进行指导;否则重复上述步骤S1‑S9,直至优化后的决策变量的系统性能的平均值大于实际样本的系统性能的平均值为止。
地址 401331 重庆市沙坪坝区大学城东路20号