马奎特法反演matlab,基于演化算法的复电阻率频谱参数反演

1引言复电阻率法(简称CR法)作为频率域激电法中的一种,是最近十几年发展起来的。它在相当宽的频段(2-8~210Hz)范围内测量复电阻率幅值和相位离散频谱,克服了一般激电法信息量少以及不能满足电磁耦合校正和评价激电异常要求的缺点,能较好地分离、提取激电效应和电磁效应。复电阻率法作为直接检测构造的含油气性方法,其勘探成本低,周期短,见效快,在石油勘探和油气预测中的应用被逐渐推广。该方法的基本原理、野外资料采集仪器及施工技术已日益成熟,准确地提取谱参数、进而分别提取激电效应和电磁效应仍是当前待研究的一个重要内容[1]。国内外学者在频谱参数的反演方面做了大量的工作,并且取得了极大的成果。1988年Jaggar等人提出了基于马奎特算法的反演[2]。2001年Xiang等人提出一种新的直接反演算法,该方法不需要对参数初始值进行估计[3]。2002年阮百尧提出一种新的基于一个ColeCole模型的递推反演算法,这种方法不需要任何中间变量的介入,也不要对参数初始值进行估计[4]。前面几种方法都只对一个ColeCole模型反演有一定的效果,本文提出基于郭涛算法[5](演化计算的一种)的频谱参数反演方法,这种方法不需要任何中间变量的介入,也不要对参数初始值进行估计,只要给定参数的大概范围,就能保证解的精确性,稳定性和较快的收敛性,该方法不仅对一个ColeCole模型的频谱参数反演效果好,同时还对用于分离电磁耦合效应的两个ColeCole模型的频谱参数反演效果也很理想,并且能允许观测数据存在一定的误差。2方法原理21郭涛算法演化计算(EvolutionaryComputation,EC)是模拟自然界的生物进化过程而发展起来的一种计算方法,1999年郭涛提出了一类基于子空间搜索(多父体重组)和群体爬山法结合的群体随机搜索算法用来求解带不等式约束条件的函数优化问题[5]。该算法的基本原理:随机产生N个初始模型(通常要求它们均匀分布在可行解空间中),计算出每个模型的适应值,并找出其中最好的和最差的模型;然后,在随机的N个初始模型中选出m个模型,进行多父杂交,产生一个子体,让它与初始模型中最好的和最差的比较,淘汰最差的;一次次循环下去,最终得到的模型就是最好的模型,也就是我们所需的解。其具体算法见参考文献[5~7]。N的选取,可根据问题的维数n和f(x)场景的复杂性而定,当n较大且场景复杂时,N可取大些,反之,则取小些,一般20N150。m的选取,根据经验m=7、8、9或10较合适。22郭涛算法的特点[7]1)算法采用了演化计算中的群体搜索策略,保证了搜索空间的全局性。算法采用了随机子空间的随机搜索策略,特别是子空间中随机搜索的非凸性,即X=mi=1aiXimi=1ai=1,-0.5ai1.5使算法的子空间可覆盖多父体的凸组合空间,保证了随机搜索的遍历性,即解空间中不存在搜索不到的死角。2)算法采用了劣汰策略,每次只把群体中适应最差(目标函数值最大)的个体淘汰出局。淘汰压力最小,既保证群体的多样性,也保证了适应值最好(目标函数值最小)的个体可以永远生存下去。这种群体爬山策略,保证了整个群体最后集体达到最深的谷底。当最优解不唯一时,算法可以一次同时找到多个最优解。3复电阻率频谱参数的反演实现31单个柯尔模型[8]W.H.Pelton等基于对大量岩、矿石标本及露头测定结果认为,对岩、矿石由激电效应引起的复电阻率随频率的变化可以用柯尔柯尔模型近似描述:(i)=(0)1-m1-11+(i)c(1)式中:(i)供交流电时,可极化岩石的复电阻率;(0)交流电频率为零赫兹的电阻率值;m充