CN108256688A - 改进萤火虫算法优化grnn的蒸汽驱开发效果预测方法 - Google Patents
改进萤火虫算法优化grnn的蒸汽驱开发效果预测方法 Download PDFInfo
- Publication number
- CN108256688A CN108256688A CN201810108905.9A CN201810108905A CN108256688A CN 108256688 A CN108256688 A CN 108256688A CN 201810108905 A CN201810108905 A CN 201810108905A CN 108256688 A CN108256688 A CN 108256688A
- Authority
- CN
- China
- Prior art keywords
- firefly
- grnn
- formula
- swarm algorithm
- development effect
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011161 development Methods 0.000 title claims abstract description 63
- 230000000694 effects Effects 0.000 title claims abstract description 51
- 238000010795 Steam Flooding Methods 0.000 title claims abstract description 50
- 241000131894 Lampyris noctiluca Species 0.000 title claims abstract description 43
- 238000005457 optimization Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 30
- 241000254158 Lampyridae Species 0.000 claims abstract description 75
- 238000009499 grossing Methods 0.000 claims abstract description 12
- 230000006870 function Effects 0.000 claims description 16
- 238000010793 Steam injection (oil industry) Methods 0.000 claims description 12
- 239000003921 oil Substances 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 10
- 230000006872 improvement Effects 0.000 claims description 9
- 238000010521 absorption reaction Methods 0.000 claims description 7
- 230000035772 mutation Effects 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 5
- 230000035699 permeability Effects 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 239000010779 crude oil Substances 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 abstract description 3
- 238000013528 artificial neural network Methods 0.000 description 9
- 238000010025 steaming Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0418—Architecture, e.g. interconnection topology using chaos or fractal principles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- General Business, Economics & Management (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Marketing (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Tourism & Hospitality (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及的是改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,它依据现场收集的蒸汽驱开发效果数据,通过在标准萤火虫算法中引入混沌优化算子初始化种群,同时引入差分变异算子完成萤火虫个体自适应变异,增加了种群多样性,加快了收敛速度,提高了种群的可进化能力,利用改进萤火虫算法优化GRNN网络的光滑因子σ,优化目标为拟合误差的均方根误差,建立最优的蒸汽驱开发效果的预测模型,实现蒸汽驱开发效果预测。本发明通过采用改进萤火虫算法优化GRNN对蒸汽驱进行开发效果的预测,提高了预测的速度和精确,可以指导现场蒸汽驱开发方案的编制,实现蒸汽驱开发方案的优化和调整,最大限度的提高经济效益。
Description
技术领域
本发明涉及一种蒸汽驱开发效果预测方法,具体涉及改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法。
背景技术
蒸汽驱作为一种高效开发方式,是稠油油藏开采的主要发展方向,而实现蒸汽驱开发效果预测是蒸汽驱成功的关键。目前蒸汽驱开发效果预测方法有基于矿场资料的经验方法、室内实验模拟方法和数值模拟方法。经验方法操作简单,但是预测精确度低;室内实验方法与现场吻合度低,预测结果与实际情况相差较大。当前,有一些数值模拟软件可用于蒸汽驱的开发效果预测,但是这些软件往往本身很昂贵,而且运行需要应用大型计算机,还会花费大量的运算时间。
随着人工神经网络理论和优化技术的发展,各种新技术和新应用不断涌现。目前已有成功将人工神经网络理论和优化技术相融合用于解决很多工程预测问题。传统的BP神经网络存在学习速度慢、易陷入局部极小值,预测结果精度不高的缺点,广义回归神经网络(GRNN)具有很强的非线性逼近能力和自学习、自适应等特性,能自动从历史数据中提取有关预测的知识,因而非常适合解决预测领域中的一些复杂问题。
发明内容
本发明的目的是提供改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,这种改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法用于解决现有蒸汽驱开发效果预测方法准确度低,成本高的问题。
本发明解决其技术问题所采用的技术方案是:这种改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法:
步骤一:收集现场蒸汽驱开发效果数据作为样本数据,包括原始含油饱和度、原油粘度(mPa·s)、油层深度(m)、渗透率(mD)、油层有效厚度(m)、原始温度(℃)、原始压力(MPa)、井网面积(m2)、注汽速率(m3/d)、注汽干度、注汽压力(MPa),选取累积油汽比作为指标来衡量蒸汽驱开发效果的好坏程度;
步骤二:将所步骤1收集的样本数据根据公式(1)映射到[-1,1]之间,对样本数据作归一化处理得到归一化数据Xnormal:
式中,Xmin、Xmax分别代表蒸汽驱开发效果预测数据集中最小和最大值,X为原始预测数据集中的数据;
步骤三:采用改进萤火虫算法搜寻最优的GRNN中的光滑因子σ,从而获得最优GRNN;改进萤火虫算法是在标准的萤火虫算法中引入混沌优化算子进行种群初始化,同时引入差分变异算子实现萤火虫个体的自适应变异,增加了种群多种性,加快了收敛速度,提高可进化能力;
步骤四:将归一化后的数据作为训练样本,用来训练最优GRNN,从而获得最优的蒸汽驱开发效果预测模型,实现蒸汽驱开发效果预测。
上述方案中改进萤火虫算法中引入混沌优化算子进行种群初始化的具体方法:
步骤1:随机产生D维向量Y0,Y0=(y01,y02,…,y0D),其中D为优化问题的维数,Y0中不能包含0.25、0.5和0.75这3个数值;
步骤2:取混沌序列:根据Logistic映射函数Yn+1=μ·Yn·(1-Yn),其中μ是混沌吸引子,取值为4,得Y1,Y2,Y3,…YN,获得混沌序列Y=(Y1,Y2,Y3…,YN);
步骤3:根据公式(2)将混沌序列Yi的第j个分量载波到优化变量的取值范围(xmin,j,xmax,j)内,xmin,j,xmax,j分别为优化变量的第j个分量的下限值和上限值,得到X1,X2,X3,…,XN,即N个萤火虫个体。
xi,j=xmin,j+yi,j·(xmax,j-xmin,j)i=1,2,…,N;j=1,2,…,D (2)
上述方案中引入差分变异算子实现萤火虫个体的自适应变异,差分变异算子的具体实现方法为:
步骤1:根据公式(3)计算萤火虫个体i的自适应变异步长ωi;
式中:N为萤火虫数目,j为萤火虫个体i按照亮度由小到大的排名;排名靠后者质量相对较高,变异尺度相对较小,排名靠前者质量相对较低,变异尺度相对较大;
步骤2:随机选择2个萤火虫m和k,它们所对应的位置为Xm和Xk,利用差分算法思想,根据公式(4)实现萤火虫个体i的自适应变异,应用式(4)计算Xi变异后萤火虫个体i的位置Xi';最后比较Xi'与Xi,保留较好的,
Xi'=Xi+ωi·(Xm-Xk) (4)
由公式(4)看出,变异步长ωi越大,变异后的萤火虫新位置Xi'与原位置Xi变化的程度也就较大,帮助跳出局部最优;反之,变异步长ωi越小,变异后的萤火虫新位置Xi'与原位置Xi变化的程度也就较小,从而实现算法的局部搜索,可以提高全局搜索能力。
上述方案中改进萤火虫算法的具体实施步骤如下:
步骤1:初始化:
1)设置萤火虫数目N,最大吸引度β0,光强吸收系数γ,步长因子α,最大迭代次数Max_Iter,搜索精度ε;
2)混沌优化算子初始化萤火虫的位置;
2.1)随机产生D维向量Y0,Y0=(y01,y02,…,y0D),其中D为优化问题的维数,Y0中不能包含0.25、0.5和0.75这3个数值;
2.2)取混沌序列:根据Logistic映射函数Yn+1=μ·Yn·(1-Yn),其中μ是混沌吸引子,取值为4,得Y1,Y2,Y3,…YN,由此获得混沌序列Y=(Y1,Y2,Y3…,YN);
2.3)根据公式(2)将混沌序列Yi的第j个分量载波到优化变量的取值范围(xmin,j,xmax,j)内,xmin,j,xmax,j分别为优化变量的第j个分量的下限值和上限值,得到X1,X2,X3,…,XN,即N个萤火虫个体。
xi,j=xmin,j+yi,j·(xmax,j-xmin,j)i=1,2,…,N;j=1,2,…,D (2)
步骤2:迭代寻优:
1)根据目标函数即公式(5),计算萤火虫i的相对亮度Ii;
Ii=f(Xi)(i=1,2,…,N) (5)
式中:f(Xi)表示萤火虫i所在位置Xi的目标函数值;
2)根据公式(6)计算群体中萤火虫i和j的距离rij;
式中:Xi和Xj为萤火虫i和j的位置,D为求解问题的空间维度,xik为第i个萤火虫在k维的坐标;
3)根据式(7)计算群体中萤火虫的吸引度β;
式中:γ为光强吸收系数;rij为萤火虫i和j的距离;β0为最大吸引度;
4)依据相对亮度,根据公式(8)更新萤火虫的空间位置;
Xi_new=Xi+β·(Xj-Xi)+α·(rand-0.5) (8)
式中:Xi和Xj为萤火虫i和j的位置,β为吸引度,α为步长因子,rand是产生[0,1]之间任意数的随机函数;
5)萤火虫种群按照亮度进行排序,执行差分变异算子;
6)进入迭代寻优,重复执行步骤1)-5);
步骤3:终止判断,直至当前迭代次数等于最大迭代次数Max_Iter,或已达到搜索精度ε要求,寻优结束,输出寻优结果。
本发明具有以下有益效果:
1、本发明采用改进萤火虫算法优化GRNN的光滑因子σ,从而获得最优GRNN,建立最优的蒸汽驱开发效果预测模型,实现蒸汽驱开发效果的快速准确的预测,进一步指导现场蒸汽驱开发方案的编制,实现蒸汽驱开发方案的优化和调整,最大限度的提高经济效益。
2、本发明采用GRNN相比于传统的BP神经网络预测,具有泛化能力好、训练速度快、参数调节方便等优点。
3、本发明改进萤火虫算法引入混沌优化算子初始化种群,同时引入差分变异算子实现萤火虫个体的自适应变异,增加了种群的多样性,更容易跳出局部极值,快速找到全局最优解。
4、本发明采用的改进萤火虫优化GRNN的光滑因子σ,建立蒸汽驱开发效果预测模型,预测效果稳定,预测精度高。
附图说明
图1为基于改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法流程图;
图2为改进萤火虫算法(IFA)的流程图;
图3为两种模型预测(FA-GRNN和IFA-GRNN)的均方根误差;
图4为标准萤火虫算法优化GRNN的蒸汽驱开发效果预测模型(FA-GRNN模型)的预测结果;
图5为改进萤火虫算法优化GRNN的蒸汽驱开发效果预测模型(IFA-GRNN模型)的预测结果。
具体实施方式
下面结合附图对本发明作进一步的说明:
参阅图1,这种改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,具体步骤如下:
步骤一:收集现场蒸汽驱开发效果数据作为样本数据,包括原始含油饱和度、原油粘度(mPa·s)、油层深度(m)、渗透率(mD)、油层有效厚度(m)、原始温度(℃)、原始压力(MPa)、井网面积(m2)、注汽速率(m3/d)、注汽干度、注汽压力(MPa)共11个主要因素,选取累积油汽比作为指标来衡量蒸汽驱开发效果的好坏程度;
步骤二:由于收集的数据往往不是在同一个数量级,为了得到较好的预测效果,需要将所收集的数据映射到[-1,1]之间,根据公式(1)对样本数据作如下归一化处理:
式中,Xmin、Xmax分别代表蒸汽驱开发效果预测数据集中最小和最大值,X为原始预测数据集中的数据,Xnormal表示归一化数据。
步骤三:采用改进萤火虫算法搜寻最优的GRNN中的光滑因子σ,从而获得最优GRNN;
改进萤火虫算法是在标准的萤火虫算法中引入混沌优化算子进行种群初始化,同时引入差分变异算子实现萤火虫个体的自适应变异,增加了种群多种性,加快了收敛速度,提高了可进化能力。
所述的改进萤火虫算法的具体实施步骤如下:
步骤1:初始化:
1)设置萤火虫数目N,最大吸引度β0,光强吸收系数γ,步长因子α,最大迭代次数Max_Iter,搜索精度ε;
2)混沌优化算子初始化萤火虫的位置;
2.1)随机产生D维向量Y0,Y0=(y01,y02,…,y0D),其中D为优化问题的维数,Y0中不能包含0.25、0.5和0.75这3个数值;
2.2)取混沌序列:根据Logistic映射函数Yn+1=μ·Yn·(1-Yn),其中μ是混沌吸引子,取值为4,得Y1,Y2,Y3,…YN,由此获得混沌序列Y=(Y1,Y2,Y3…,YN);
2.3)根据公式(2)将混沌序列Yi的第j个分量载波到优化变量的取值范围(xmin,j,xmax,j)内,xmin,j,xmax,j分别为优化变量的第j个分量的下限值和上限值,得到X1,X2,X3,…,XN,即N个萤火虫个体。
xi,j=xmin,j+yi,j·(xmax,j-xmin,j)i=1,2,…,N;j=1,2,…,D (2)
步骤2:迭代寻优:
1)根据目标函数即公式(5),计算萤火虫i的相对亮度Ii;
Ii=f(Xi)(i=1,2,…,N) (5)
式中:f(Xi)表示萤火虫i所在位置Xi的目标函数值;
2)根据公式(6)计算群体中萤火虫i和j的距离rij;
式中:Xi和Xj为萤火虫i和j的位置,D为求解问题的空间维度,xik为第i个萤火虫在k维的坐标;
3)根据式(7)计算群体中萤火虫的吸引度β;
式中:γ为光强吸收系数;rij为萤火虫i和j的距离;β0为最大吸引度;
4)依据相对亮度,根据公式(8)更新萤火虫的空间位置;
Xi_new=Xi+β·(Xj-Xi)+α·(rand-0.5) (8)
式中:Xi和Xj为萤火虫i和j的位置,β为吸引度,α为步长因子,rand是产生[0,1]之间任意数的随机函数;
5)萤火虫种群按照亮度进行排序,执行差分变异算子;
差分变异算子更进一步的具体实现为:
5.1)根据公式(3)计算萤火虫个体i的自适应变异步长ωi;
式中:N为萤火虫数目,j为萤火虫个体i按照亮度由小到大的排名;排名靠后者质量相对较高,变异尺度相对较小,排名靠前者质量相对较低,变异尺度相对较大;
5.2)随机选择2个萤火虫m和k,它们所对应的位置为Xm和Xk,利用差分算法思想,根据公式(4)实现萤火虫个体i的自适应变异,应用式(4)计算Xi变异后萤火虫个体i的位置X'i;最后比较X'i与Xi,保留较好的,
X'i=Xi+ωi·(Xm-Xk) (4)
由公式(4)看出,变异步长ωi越大,变异后的萤火虫新位置Xi'与原位置Xi变化的程度也就较大,帮助跳出局部最优;反之,变异步长ωi越小,变异后的萤火虫新位置Xi'与原位置Xi变化的程度也就较小,从而实现算法的局部搜索,可以提高全局搜索能力。
6)进入迭代寻优,重复执行步骤1)-5);
步骤3:终止判断,直至当前迭代次数等于最大迭代次数Max_Iter,或已达到搜索精度ε要求,寻优结束,输出寻优结果。
步骤四:将归一化后的数据作为训练样本,用来训练最优GRNN,从而获得最优的蒸汽驱开发效果预测模型,实现蒸汽驱开发效果预测。
本发明的原理:人工神经网络(Artificial Neural Network ANN)具有良好的非线性、鲁棒性、容错性、自组织、自适应能力等特点,已成为解决预测问题的有力工具。广义回归神经网络(GRNN)是基于径向基函数神经网络的一种改进。它具有很强的非线性映射能力和柔性网络结构,以及高度的容错性和鲁棒性。与其他神经网络相比,训练过程更为简单,只需确定训练样本和光滑因子σ即可完成网络构建和预测。它具有较高的学习速度、拟合能力强。因而非常适合解决预测领域中的一些复杂问题。
GRNN的光滑因子σ对预测的效果影响较大,需要不断尝试才能获得最佳值。为了提高蒸汽驱开发效果预测的速度和精度,本发明采用改进萤火虫算法优化GRNN的光滑因子σ,从而获得最优GRNN,建立最优的蒸汽驱开发效果预测模型,实现蒸汽驱开发效果的快速准确的预测。
本发明进行的保密性实验:
以辽河油田齐40块为研究对象,通过现场专家分析和研究,确定11个参数即原始含油饱和度、原油粘度(mPa·s)、油层深度(m)、渗透率(mD)、油层有效厚度(m)、原始温度(℃)、原始压力(MPa)、井网面积(m2)、注汽速率(m3/d)、注汽干度、注汽压力(MPa)作为模型的主要影响因素,选取累积油汽比来衡量预测效果的好坏程度,收集100个井组数据,部分样本数据如表1所示。选取90个井组数据作为训练样本和10个井组数据作为测试样本,建立预测模型数据集。
表1部分样本的原始数据
为了验证所提出的改进萤火虫算法(IFA)优化GRNN蒸汽驱开发效果预测模型的性能,采用Matlab2014a开发环境,利用Matlab神经网络GRNN工具箱中的newgrnn函数,建立了两种预测模型,分别是:标准萤火虫算法优化GRNN的蒸汽驱开发效果预测模型(FA-GRNN模型)和改进萤火虫算法优化GRNN的蒸汽驱开发效果预测模型(IFA-GRNN模型)。两种模型采用相同的参数设置:萤火虫数目N=100,最大吸引度β0=1,光强吸收系数γ=1,步长因子α=0.02,最大迭代次数Max_Iter=200,搜索精度ε=1e-08。
从图3中可以看出:FA-GRNN模型预测的均方根误差RMSE=9.71086E-09,迭代次数为169,耗时501.70S;而IFA-GRNN模型预测的均方根误差RMSE=6.13E-09,迭代次数为97,耗时287.96S。因此,IFA-GRNN模型相对于FA-GRNN模型,预测速度更快且预测精度更高。
将标准萤火虫算法FA找到的最佳光滑因子σ值代入到GRNN模型中,进行剩下10组样本的蒸汽驱开发效果预测,结果如图4所示。由图4中可知,大部分样本FA-GRNN模型的预测值和期望值之间比较接近,但有个别样本与期望值偏差较大些,预测效果不是很理想。
将改进萤火虫算法IFA找到的最佳光滑因子σ值代入到GRNN模型中,进行剩下10组样本的蒸汽驱开发效果预测,结果如图5所示。由图5中可知,IFA-GRNN模型的预测值和期望值之间最为接近,预测效果非常好;FA-GRNN模型比IFA-GRNN模型差一些。
从实验结果可以看出,依据现场采集的蒸汽驱开发效果数据,采用改进萤火虫算法寻找GRNN的最佳光滑因子σ,建立的基于改进萤火虫算法优化GRNN的蒸汽驱开发效果预测模型,改善了GRNN泛化能力,在蒸汽驱开发效果预测方面具有较好的预测能力,提高了预测的准确性,为蒸汽驱开发效果预测提供了一种科学、有效的方法。
Claims (4)
1.一种改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,其特征在于:这种改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法:
步骤一:收集现场蒸汽驱开发效果数据作为样本数据,包括原始含油饱和度、原油粘度、油层深度、渗透率、油层有效厚度、原始温度、原始压力、井网面积、注汽速率、注汽干度、注汽压力,选取累积油汽比作为指标来衡量蒸汽驱开发效果的好坏程度;
步骤二:将所步骤1收集的样本数据,根据公式(1)映射到[-1,1]之间,对样本数据作归一化处理得到归一化数据Xnormal:
式中,Xmin、Xmax分别代表蒸汽驱开发效果预测数据集中最小和最大值,X为原始预测数据集中的数据;
步骤三:采用改进萤火虫算法搜寻最优的GRNN中的光滑因子σ,从而获得最优GRNN;改进萤火虫算法是在标准的萤火虫算法中引入混沌优化算子进行种群初始化,同时引入差分变异算子实现萤火虫个体的自适应变异,增加了种群多种性,加快了收敛速度,提高可进化能力;
步骤四:将归一化后的数据作为训练样本,用来训练最优GRNN,从而获得最优的蒸汽驱开发效果预测模型,实现蒸汽驱开发效果预测。
2.根据权利要求1所述的改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,其特征在于:所述的改进萤火虫算法中引入混沌优化算子进行种群初始化的具体方法:
步骤1:随机产生D维向量Y0,Y0=(y01,y02,…,y0D),其中D为优化问题的维数,Y0中不能包含0.25、0.5和0.75这3个数值;
步骤2:取混沌序列:根据Logistic映射函数Yn+1=μ·Yn·(1-Yn),其中μ是混沌吸引子,取值为4,得Y1,Y2,Y3,…YN,获得混沌序列Y=(Y1,Y2,Y3…,YN);
步骤3:根据公式(2)将混沌序列Yi的第j个分量载波到优化变量的取值范围(xmin,j,xmax,j)内,xmin,j,xmax,j分别为优化变量的第j个分量的下限值和上限值,得到X1,X2,X3,…,XN,即N个萤火虫个体。
xi,j=xmin,j+yi,j·(xmax,j-xmin,j)i=1,2,…,N;j=1,2,…,D (2)。
3.根据权利要求1所述的改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,其特征在于:所述的差分变异算子实现萤火虫个体的自适应变异,差分变异算子的具体方法:
步骤1:根据公式(3)计算萤火虫个体i的自适应变异步长ωi;
式中:N为萤火虫数目,j为萤火虫个体i按照亮度由小到大的排名;排名靠后者质量相对较高,变异尺度相对较小,排名靠前者质量相对较低,变异尺度相对较大;
步骤2:随机选择2个萤火虫m和k,它们所对应的位置为Xm和Xk,利用差分算法思想,根据公式(4)实现萤火虫个体i的自适应变异,应用式(4)计算Xi变异后萤火虫个体i的位置X′i;最后比较X′i与Xi,保留较好的,
X′i=Xi+ωi·(Xm-Xk) (4)
由公式(4)看出,变异步长ωi越大,变异后的萤火虫新位置X′i与原位置Xi变化的程度也就较大,帮助跳出局部最优;反之,变异步长ωi越小,变异后的萤火虫新位置X′i与原位置Xi变化的程度也就较小,从而实现算法的局部搜索,可以提高全局搜索能力。
4.根据权利要求3所述的改进萤火虫算法优化GRNN的蒸汽驱开发效果预测方法,其特征在于:所述的改进萤火虫算法的具体实施步骤如下:
步骤1:初始化:
1)设置萤火虫数目N,最大吸引度β0,光强吸收系数γ,步长因子α,最大迭代次数Max_Iter,搜索精度ε;
2)混沌优化算子初始化萤火虫的位置;
2.1)随机产生D维向量Y0,Y0=(y01,y02,…,y0D),其中D为优化问题的维数,Y0中不能包含0.25、0.5和0.75这3个数值;
2.2)取混沌序列:根据Logistic映射函数Yn+1=μ·Yn·(1-Yn),其中μ是混沌吸引子,取值为4,得Y1,Y2,Y3,…YN,由此获得混沌序列Y=(Y1,Y2,Y3…,YN);
2.3)根据公式(2)将混沌序列Yi的第j个分量载波到优化变量的取值范围(xmin,j,xmax,j)内,xmin,j,xmax,j分别为优化变量的第j个分量的下限值和上限值,得到X1,X2,X3,…,XN,即N个萤火虫个体。
xi,j=xmin,j+yi,j·(xmax,j-xmin,j)i=1,2,…,N;j=1,2,…,D (2)
步骤2:迭代寻优:
1)根据目标函数即公式(5),计算萤火虫i的相对亮度Ii;
Ii=f(Xi) (i=1,2,…,N) (5)
式中:f(Xi)表示萤火虫i所在位置Xi的目标函数值;
2)根据公式(6)计算群体中萤火虫i和j的距离rij;
式中:Xi和Xj为萤火虫i和j的位置,D为求解问题的空间维度,xik为第i个萤火虫在k维的坐标;
3)根据式(7)计算群体中萤火虫的吸引度β;
式中:γ为光强吸收系数;rij为萤火虫i和j的距离;β0为最大吸引度;
4)依据相对亮度,根据公式(8)更新萤火虫的空间位置;
Xi_new=Xi+β·(Xj-Xi)+α·(rand-0.5) (8)
式中:Xi和Xj为萤火虫i和j的位置,β为吸引度,α为步长因子,rand是产生[0,1]之间任意数的随机函数;
5)萤火虫种群按照亮度进行排序,执行差分变异算子;
6)进入迭代寻优,重复执行步骤1)-5);
步骤3:终止判断,直至当前迭代次数等于最大迭代次数Max_Iter,或已达到搜索精度ε要求,寻优结束,输出寻优结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810108905.9A CN108256688B (zh) | 2018-02-05 | 2018-02-05 | 改进萤火虫算法优化grnn的蒸汽驱开发效果预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810108905.9A CN108256688B (zh) | 2018-02-05 | 2018-02-05 | 改进萤火虫算法优化grnn的蒸汽驱开发效果预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108256688A true CN108256688A (zh) | 2018-07-06 |
CN108256688B CN108256688B (zh) | 2022-02-25 |
Family
ID=62743673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810108905.9A Expired - Fee Related CN108256688B (zh) | 2018-02-05 | 2018-02-05 | 改进萤火虫算法优化grnn的蒸汽驱开发效果预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108256688B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109146053A (zh) * | 2018-08-30 | 2019-01-04 | 广西民族大学 | 一种模拟生物理想自由分布模型的萤火虫方法 |
CN109471362A (zh) * | 2018-12-18 | 2019-03-15 | 东北大学 | 一种冷热电联产优化系统及方法 |
CN114565237A (zh) * | 2022-02-15 | 2022-05-31 | 石河子大学 | 用于多能组合发电站的综合能源低碳调度方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105608271A (zh) * | 2015-12-21 | 2016-05-25 | 武汉大学 | 一种基于分解与优化的短期风速时序预测方法 |
US20160157146A1 (en) * | 2008-07-07 | 2016-06-02 | Odyssey Wireless, Inc. | Systems/methods of spatial multiplexing |
CN106920241A (zh) * | 2017-03-13 | 2017-07-04 | 华北电力大学 | 一种汽包水位监控方法、装置及系统 |
-
2018
- 2018-02-05 CN CN201810108905.9A patent/CN108256688B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160157146A1 (en) * | 2008-07-07 | 2016-06-02 | Odyssey Wireless, Inc. | Systems/methods of spatial multiplexing |
CN105608271A (zh) * | 2015-12-21 | 2016-05-25 | 武汉大学 | 一种基于分解与优化的短期风速时序预测方法 |
CN106920241A (zh) * | 2017-03-13 | 2017-07-04 | 华北电力大学 | 一种汽包水位监控方法、装置及系统 |
Non-Patent Citations (6)
Title |
---|
倪红梅等: "基于HABC-RBF神经网络的蒸汽驱预测方法", 《吉林大学学报(信息科学版)》 * |
张兆祥等: "稠油油藏蒸汽驱评价新方法", 《稠油油藏蒸汽驱评价新方法》 * |
张军丽等: "《人工萤火虫与差分进化混合优化算法》", 《信息与控制》 * |
李平科等: "蒸汽驱开发采收率预测新方法", 《蒸汽驱开发采收率预测新方法》 * |
王昕等: "《基于萤火虫算法—广义回归神经网络的光伏发电功率组合预测》", 《电网技术》 * |
高伟明: "《萤火虫算法的研究与应用》", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109146053A (zh) * | 2018-08-30 | 2019-01-04 | 广西民族大学 | 一种模拟生物理想自由分布模型的萤火虫方法 |
CN109471362A (zh) * | 2018-12-18 | 2019-03-15 | 东北大学 | 一种冷热电联产优化系统及方法 |
CN114565237A (zh) * | 2022-02-15 | 2022-05-31 | 石河子大学 | 用于多能组合发电站的综合能源低碳调度方法及系统 |
CN114565237B (zh) * | 2022-02-15 | 2023-09-12 | 石河子大学 | 用于多能组合发电站的综合能源低碳调度方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108256688B (zh) | 2022-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106951983A (zh) | 基于采用多父代遗传算法的人工神经网络的喷射器性能预测方法 | |
CN108256688A (zh) | 改进萤火虫算法优化grnn的蒸汽驱开发效果预测方法 | |
CN107829718A (zh) | 基于均衡水驱理念的油藏井网及注采方案优化设计方法 | |
CN107657243B (zh) | 基于遗传算法优化的神经网络雷达一维距离像目标识别方法 | |
CN110405343A (zh) | 一种基于Bagging集成的预测模型与粒子群优化算法的激光焊接工艺参数优化方法 | |
CN110363344A (zh) | 基于miv-gp算法优化bp神经网络的概率积分参数预测方法 | |
CN106779211A (zh) | 一种用于油田注采井网驱替开发的产能预测方法 | |
CN110348075A (zh) | 一种基于改进支持向量机算法的磨削表面粗糙度预测方法 | |
CN105045941A (zh) | 基于无迹卡尔曼滤波的抽油机参数优化方法 | |
CN105046326B (zh) | 基于功图主元分析的抽油机参数优化方法 | |
CN108875159A (zh) | 一种压电驱动器的率相关迟滞建模方法 | |
CN111259517A (zh) | 一种隧道爆破设计方法、装置及设备 | |
CN110211173A (zh) | 一种基于深度学习的古生物化石定位和识别方法 | |
CN105808311A (zh) | 一种基于降维策略的油藏模拟快速拟合方法 | |
CN104680025B (zh) | 基于遗传算法极限学习机的抽油机参数优化方法 | |
CN108843312A (zh) | 油田储层层内非均质性综合评价方法 | |
CN107341521A (zh) | 一种基于煤炭光谱数据对煤炭分类的方法 | |
US20140257577A1 (en) | Optimization of a multi-period model for valuation applied to flow control valves | |
CN112766608A (zh) | 一种基于极限学习机的水泥磨系统电耗指标预测方法 | |
Zhang et al. | Smart well pattern optimization using gradient algorithm | |
Plaksina et al. | Practical handling of multiple objectives using evolutionary strategy for optimal placement of hydraulic fracture stages in unconventional gas reservoirs | |
Anderson et al. | EMPOWER-1.0: an efficient model of planktonic ecosystems written in R | |
CN104680023B (zh) | 基于多目标决策的抽油机参数优化方法 | |
CN107729988A (zh) | 基于动态深度置信网络的蓝藻水华预测方法 | |
CN108843296A (zh) | 一种基于多因素影响下的单井重复压裂效果预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220225 |