发明名称 基于粒子群与线性加法器结合的地震数据反演方法
摘要 基于粒子群与线性加法器结合的地震数据反演方法,包括步骤1.设定粒子空间及其初始参数;步骤2.产生粒子空间中各个粒子的个体极值及全局极值;步骤3.拟合计算出每个粒子的粒子适用值;步骤4.比较并优化粒子空间中各个粒子的个体极值及全局极值;步骤5.将步骤4中得到的新的个体极值与全局极值带回步骤2,重复步骤2-4,完成一次迭代,重复多次迭代,直到满足结束条件。采用本发明所述的基于粒子群与线性加法器结合的地震数据反演方法,使用粒子群算法仅需要调整两个参数即可实现对非线性问题的求解,得到反演结果并应用于地震数据资料解释中,以达到准确预测地层信息的目的。
申请公布号 CN103454677B 申请公布日期 2016.03.02
申请号 CN201310425517.0 申请日期 2013.09.18
申请人 中海石油(中国)有限公司湛江分公司;成都晶石石油科技有限公司 发明人 谢玉洪;刘力辉;周家雄;陈殿远;刘兵;马光克;刘仕友
分类号 G01V1/28(2006.01)I 主分类号 G01V1/28(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 梁田
主权项 基于粒子群与线性加法器结合的地震数据反演方法,包括测量一一对应的多个数据点(Xi,Yi),其中Xi,Yi均为向量,i=1、2、…m,m为数据量个数,其特征在于,还包括如下步骤:步骤1.设定解空间的惯性权重参数ω,加速系数c1、c2,设定粒子群中各个粒子的初始个体极值和粒子群的初始全局极值,设定各个粒子的初始速度;步骤2.以Xi为基础,利用①式和②式,随机产生各个粒子的初始位置<img file="920823dest_path_image001.GIF" wi="23" he="24" />与初始速度<img file="dest_path_image002.GIF" wi="21" he="24" />,每个粒子的初始个体极值作为其初始位置, 全部个体极值中的最佳值为初始全局极值;<img file="650355dest_path_image003.GIF" wi="372" he="29" />‑‑‑‑‑①<img file="dest_path_image004.GIF" wi="99" he="24" />‑‑‑‑‑‑②其中<img file="784664dest_path_image005.GIF" wi="23" he="24" />为Xi在第d维的分量,<img file="dest_path_image006.GIF" wi="21" he="24" />为Xi对应的个体极值在第d维的分量,<img file="601311dest_path_image007.GIF" wi="24" he="27" />为全局极值在第d维的分量,rand()是在[0,1]范围内变化的随机数;上标n或n+1为迭代次数,n=0时,表示相应参数的初始值;步骤3.根据步骤2中的Xi拟合出<img file="dest_path_image008.GIF" wi="18" he="24" />,将<img file="649907dest_path_image008.GIF" wi="18" he="24" />与Yi比较,计算出每个粒子的粒子适用值;步骤4.比较粒子适用值与它的个体极值,若前者优于后者,则以粒子适用值作为该粒子新的个体极值;比较各个粒子的粒子适用值与全局极值,若前者中有优于后者的粒子,则以该粒子的粒子适用值作为新的全局极值;步骤5.将步骤4中得到的新的个体极值与全局极值带回步骤2,重复步骤2‑4,完成一次迭代,重复多次迭代,直到满足结束条件;步骤3中,计算Xi粒子的粒子适用值<img file="173292dest_path_image009.GIF" wi="51" he="23" />为利用下式;<img file="dest_path_image010.GIF" wi="174" he="42" />其中<img file="917651dest_path_image011.GIF" wi="18" he="23" />为根据全部Xi拟合出的训练值,i=1、2、…m;所述步骤3中利用多层线性计算器对<img file="769063dest_path_image011.GIF" wi="18" he="23" />进行拟合,具体为步骤301.利用<img file="dest_path_image012.GIF" wi="99" he="24" />,计算出<img file="621350dest_path_image013.GIF" wi="23" he="24" />,其中<img file="dest_path_image014.GIF" wi="99" he="24" />,计算出每一维度<img file="327138dest_path_image013.GIF" wi="23" he="24" />的分量,从而计算出<img file="475354dest_path_image015.GIF" wi="17" he="23" />;其中b为预设常数,下标i表示第i个粒子,d表示粒子的维度,上标n表示迭代次数;步骤302.利用<img file="dest_path_image016.GIF" wi="106" he="64" />计算出<img file="266592dest_path_image008.GIF" wi="18" he="24" />,其中k为输入个数。
地址 524000 广东省湛江市坡头区南调路