基于改进人工鱼群算法的互感器Jiles-Atherton模型参数辨识

1、基于改进人工鱼群算法的互感器Jiles-Atherton模型参数辨识*林国营1，宋强1，潘峰1，肖厦颖2，李开成2，王凌云2（1. 广东电网有限责任公司电力科学研究院，广州 510000； 2. 华中科技大学 电气与电子工程学院 强电磁工程与新技术国家重点实验室，武汉 430074）摘要：为开展仿真平台下的电磁式互感器特性研究，需要对试验互感器建立精确可靠的磁滞模型。Jiles-Atherton(J-A)模型广泛应用于互感器铁芯磁滞模型与仿真实验，其5个关键参数的辨识准确度直接影响J-A模型与试验互感器的拟合程度。本文提出一种基于人工鱼群和模拟退火算法的混合智能寻优算法对J-A模型进行参数辨识

2、。改进算法初期使用变步长人工鱼群算法将搜索域快速锁定在全局最优解的附近。当J-A模型拟合达到一定精度后，转而使用并行模拟退火算法继续进行局部的精确搜索。通过Matlab仿真证实：改进混合算法同时解决了鱼群算法后期寻优效率较低以及退火算法难以大范围搜索的问题，且算法稳定性较高，能有效提高电磁式互感器J-A模型参数辨识的时效性与精确度。关键词：J-A模型；参数辨识；人工鱼群算法；模拟退火算法中图分类号：TM45 文献标识码：A 文章编号：1001-1390(2018)00-0000-00Parameter Identification identification of Jiles-Atherto

3、n Model model Based based on an Improved improved Artificial artificial Fish fish Swarm swarm AlgorithmalgorithmLing Guoying1, Song Qiang1, Pan Feng1, Xiao Xiaying2, Li Kaicheng2, Wang Lingyun2(1.Electric Power Research Institute of Guangdong Power Grid Company，, Guangzhou, 510000, China. 2.State Ke

4、y Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China）)Abstract：: In order to research on the characteristics of electromagnetic transformer in the simulation platform,

5、it is necessary to establish a precise and reliable magnetic hysteresis model for test transformer. Jiles-Atherton (J-A) model is widely used in the field of hysteresis modeling and simulation experiments of ferromagnetic materials. The accuracy of the five key parameters of the J-A model directly a

6、ffects the fitting degree of the J-A model and tested sensor. This paper proposes an improved hybrid intelligent optimization algorithm, combining the artificial fish swarm algorithm (AFSA) and simulated annealing algorithm (SAA), which can identify the parameters of J-A model. In the improved hybri

7、d algorithm, the variable-step AFSA is employed to quickly locate the search domain in the vicinity of the global optimal result. When the fitting degree reaches a certain accuracy, the parallel SAA is used to search the global optimal result precisely in a small range. Simulation results show that

8、the improved hybrid algorithm are is able to solve the problem that the AFSA is not efficient enough, and it can effectively improve the time-validity and the accuracy of J-A model parameter identification with higher stability.Keywords：Jiles-Atherton model, parameters identification, AFSA, SAA0引 言

9、*基金项目：国家自然科学基金资助项目();广东电网有限责任公司电力科研究院项目(GDKJXM)电磁式互感器多使用非线性铁磁材料制作。为了为开展仿真平台下的电磁式互感器特性研究，使仿真效果尽可能地接近实际情况，需要对试验互感器搭建精确的铁芯磁滞模型。描述磁滞现象的理论模型有Lucas、JilesAtherton模型（下文均简称为J-A模型）等。其中J-A模型因其控制参数较少而被广泛应用在铁磁材料的磁滞建模与仿真中。J-A模型各方程与参量均具有明确的物理意义，可以准确描述铁芯磁滞回线的非线性关系。通过求解J-A模型的微分方程组便能够得到与真实情况拟合度较高的B-H曲线，从而实现互感器的精确建模。在

10、难以直接获取方程参数的场合中经常需要使用数学方法对参数进行辨识，即通过求解最优拟合参数来逼近理论参数。若求解结果的误差足够小，则认为拟合参数已经精确到可以取代理论参数1。J-A磁滞模型的参数虽具有实际意义，却难以通过实测直接得到，所以要对其进行准确可靠的参数辨识。目前应用于J-A模型参数辨识的常用算法主要是一些基于仿生学的智能优化算法，如遗传算法、粒子群算法、模拟退火算法等，但是对于处理复杂的优化问题，单一的智能寻优方法往往存在显著缺陷，不利于获取精确的全局最优解，如模拟退火算法对初始温度、温度衰减系数等参数设置较为敏感，且搜索后期收敛速度较慢。近年来，大多数文献采用多种算法相结合的方式来寻找

11、J-A模型的最佳参数组合。文献2使用了粒子群算法（PSO）结合Simulink动态仿真的方法来拟合J-A磁滞回线模型，并获取相关参数。该方法虽表现出较好的鲁棒性和收敛性，但是算法复杂度较高，运行时间较长。文献34采用神经网络和遗传算法相结合的方法，先使用神经网络算法训练寻优函数，再通过遗传算法进行极值寻优，相比单优化算法达到了较高的拟合精度，但也导致了算法时效性不足。文献5提出了一种结合混沌优化、遗传算法和模拟退火算法的混合优化方法，牺牲算法简便性以获取较快的收敛速度。文献6采用模糊控制方法获取动态反馈系数进行蛙跳算法步长修正，减小了原本蛙跳算法的随机性且具有一定的自适应能力，但对整体收敛品质

12、提升不高。本文使用一种基于变步长人工鱼群和并行模拟退火算法的改进混合算法进行电磁式电流互感器铁芯J-A模型的关键参数辨识。人工鱼群算法是一种通过模拟鱼群行为进行随机搜索的智能优化算法，能够快速将搜索域缩小至全局最优解的附近，但是该算法盲目性较大，导致后期寻优难以获得更为精确的结果，同时减慢了寻优速度。而基于Meteopolis准则的模拟退火算法在搜索初期温度较高时以较高的概率接受差解，从而不易陷入局部最优解；在温度较低时，允许差解的概率变低从而使搜索更为精确。本文提出的改进混合算法，在寻优前期利用人工鱼群算法快速接近全局最优解，并引入变步长因子减小算法的随机性；再使用模拟退火算法在较小范围内进

13、行局部搜索，同时利用并行控制策略进行多线程运算。改进混合算法解决了鱼群算法效率较低以及退火算法难以大范围搜索的问题，能有效提高J-A模型参数辨识的时效性与精确度。1互感器JilesAtherton模型互感器等效模型的准确性和有效性，关键在于铁芯磁滞现象的物理描述与数学表达7。目前，描述磁滞现象较为经典的模型有Lucas模型和Jiles-Atherton(J-A)模型等8-10。其中J-A模型所需辨识的参数较少，实现起来较为方便，PSCAD等多款电力仿真软件中都嵌入了基于J-A理论的互感器等值模型。该模型适用性广且便于扩展研究。在J-A理论模型的基础上又出现了多种改进策略，不同文献对于J-A理论

14、的引用和转述也不尽相同，本文采用文献7提出的改进J-A模型作为互感器建模仿真的基础。磁化强度M、磁场强度H与磁通密度B之间的关系如式（1）： （1）式中，为真空磁导率。根据能量守恒原理得到能量方程式： （2）式中He表示有效磁场强度；为表征磁畴内部耦合的平均场参数；Man为无磁滞磁化强度，即理想材料的磁化强度；Mirr表示M中的不可逆磁化分量；Mrev表示M中的可逆磁化分量。He、Man、Mirr和Mrev的表达式如下： （3） （4） （5） （6）再由式（5）和式（6）得到： （7）将式（3）式（7）带入能量守恒式（2），可得： （8）式中是一个符号函数，当时， 时，。为了消除J-A的非物

15、理解，对式（8）进行修正，修正后的J-A模型方程为： （9） （10） （11）上述公式中修正系数通常取0.96。Ms、a、c、k为J-A模型中待求的5个未知参数，都具有明确的物理意义:Ms代表饱和磁化强度，a为朗之万函数参数；为磁畴内部耦合的平均场参数；c为可逆磁化系数；k为损耗系数。因此,对于每种铁磁材料，其参数是确定的，且能准确反应该材料的磁化特性。2基于改进混合算法的智能寻优算法由上一节描述可知，互感器J-A模型参数辨识问题可以看作一个计算量非常庞大的非线性方程组求解问题，群体智能优化算法因其在大数据环境下良好的计算能力而常被应用于解决此类问题。2.1 人工鱼群算法人工鱼群算法（Art

16、ificial Fish Swarm Algorithm， AFSA）思路来自于自然界鱼群总是朝向食物浓度高的地方移动这一特征8，将计算模型的目标函数作为搜索目标，模拟鱼群的行为方式，在一定范围的搜索域中寻找最佳目标及其对应的控制变量参数。相比于其他的智能寻优算法，AFSA只将寻优目标值作为唯一辨识条件，对约束范围的依赖性小，且对初值与参数设定不敏感，适应性广9；引入变步长因子的AFSA，其移动步长根据条件灵活应变，具有自适应能力，整体收敛速度较快。将计算模型的关键参数作为人工鱼的位置坐标，若存在k个模型参数，那么在共有N条人工鱼的鱼群中，第i条人工鱼位置Xi=(x1i,x2i,xki)，其中

17、xki 为第k个参数的取值。将这些参数带入计算模型和目标函数中得到因变量Yi作为该位置的食物量（即寻优目标），每次迭代通过模拟各种鱼群行为来更新人工鱼的位置，使得鱼群朝向食物量最大的区域聚集10。以求极大值问题为例，AFSA三种典型的人工鱼行为描述如下：（1）觅食行为。人工鱼当下位置为Xi，在其感知范围内随机另选位置Xj，比较这两个位置的食物量Yi和Yj，若YjYi，则向Xj方位前进一步。否则重新选择位置Xj，若反复试验一定次数后Yj仍不能超出Yi，则判定其不满足前进条件，人工鱼随机移动一步。（2）聚群行为。设第i条人工鱼Xi当前可视范围内的伙伴数量及中心位置分别为nf和Xc，当Yc/nfYi

18、(为拥挤度)时，表示伙伴中心的食物密度较高，则Xi朝该方位移动一步，否则执行觅食行为。（3）追尾行为。类似于聚群行为，记当前可视范围的伙伴中Xj的食物量Yj最大，如果Yj/nfYi，表明伙伴Xj所处位置具有较多食物，则Xi朝该方位移动一步，否则执行觅食行为。虽然AFSA前期能较快收敛，缩小搜索域的约束范围，但其在小范围内的精确检索却较为缓慢，无法快速获取高精度要求下的最优目标值 11。这一缺陷在处理复杂优化问题时表现的尤为明显。该算法一旦进入到目标值更迭不明显的搜索域时，AFSA的收敛速度大幅减慢，难以逼近全局最优解。总的来说，AFSA后期较差的搜索性能抑制了算法的总体寻优能力，搜索至一定精度

19、后难以得到进一步优化，极大地减慢了运算速度，严重影响其搜索质量和效率12-13。2.2 改进混合智能寻优算法模拟退火（Simulated Annealing, SA）算法最早是由Metropolis等提出，后来逐渐发展成一种迭代自适应启发式概率性搜索算法14。该算法具有较高的鲁棒性且局部搜索能力强。为了使J-A模型参数辨识更加精确，尽可能地缩小全局误差，本文在人工鱼群算法的基础上引入模拟退火算法，使混合算法在前期使用AFSA进行快速收敛，后期利用SA进行局部搜索，进一步精确寻优。SA在Metropolis准则的基础上判定是否接受搜寻优过程中的产生的新解。假设计算模型的目标函数为f(x)，自变量

20、为x。当温度为T时，若当前解为x1，不同于x1的一个新解为x2，其目标差df=f(x1)-f(x2)，则Metropolis准则为： （12）如果dfYi(Xi)，则改变当前人工鱼的位置坐标Xi=Xi+rand*step*| X-Xi |，其中rand是小于1的随机数，目的是对定步长step进行自适应修正。取两种行为获取的较大目标函数值更新Yi，人工鱼随之移动到对应位置Xi；（c）获取整个鱼群中最高的目标函数Ybest=maxYi，以及Ybest对应的人工鱼位置Xbest，即得到搜索至此的最佳J-A模型参数组合；（d）若此时Ybest小于前期阈值Ymax1，则进入下一次迭代，迭代次数gen=g

21、en+1，继续执行（2）（3）的迭代过程，直至YbestYmax1，前期搜索结束，转而进入并行模拟退火算法；（4）进入并行模拟退火算法，流程图如图3所示。（a）参数初始化：初温T0，终止温度Tend，链长L，降温速率q，最大允许误差per_erro；（b）主线程和分线程同时输入初解作为模拟退火当下温度的初始解，初解为前期搜索获取的最优解，有：S1=Xbest；（c）主线程负责记录初始解，分线程进入温度迭代过程，根据Metropolis法则，进行1次解更新，更新后的解和初始解一起比较选出当前最优解ybest，以及对应的bestx和当前误差erro；图2 变步长人工鱼群算法流程图Fig.2 Flo

22、w chart of the variable-step artificial fish swarm algorithm（d）若erro依然大于计算允许的最大误差per_erro，则将当前局部最优解作为主线程和分线程下一退火温度的初始解，令T0=q*T0，继续迭代步骤（c），直至erroYmax2，结束搜索，对应的参数组合Ms，a，c，k即为J-A模型参数辨识的结果。图3 并行模拟退火算法流程图Fig.3 Flow chart of parallel simulated annealing algorithm3改进混合算法的参数辨识及误差对比分析为了验证上述改进混合算法的时效性和精确性，利用M

23、atlab工具进行该算法下的J-A模型参数辨识。仿真计算程序由Matlab2016编译运行，计算机CPU型号为i7-5930K、3.5 GHz,内存为32 GB。本文使用的实测数据来自文献6，J-A模型的五个参数的取值范围及理论值如表1。实测数据为1 000组B-H值，将H带入式（1）、式（4）和式（6），计算对应的B值，使得目标函数fitness达到全局最优解，并得到最优解所对应的五个参数a,c,k,和Ms的值，即为辨识结果。本节分别用粒子群算法、人工鱼群算法、以及改进混合算法对J-A模型进行参数辨识。其中鱼群算法及其改进算法均设定鱼群数为20只，最大迭代次数为20次。参数辨识结果对比如表2

24、。表1 J-A模型参数范围及理论值Tab.1 The range and theoretical values of J-A parameters参数取值范围理论值MS1*1062*1061.427*106a110041.9k5015097.45*10-63*10-51.366*10-5c0.010.050.0284表2 参数辨识结果及计算时间对比Tab.2 Comparison of the identified J-A parameters value and computation time 辨识方法粒子群人工鱼群改进混合算法MS/(A/m)1.450*1061.437*1061.425

25、*106a/(A/m)50.239647.0137742.3358k/(A/m)95.86793.7475995.13122.098*10-52.0*10-51.226*10-5c0.03820.044150.0195Fitness1.47001.47091.4744相对误差(%)0.4540.3930.189耗时(s)1101.2658.2400.5改进混合算法和鱼群算法的相对理论曲线的拟合情况对比如图4所示，可发现两种方法的曲线拟合效果均较好，几乎无可视误差。图4 B-H曲线辨识结果对比Fig.4 Comparison of identified B-H curve由表1可以看出，改进混合

26、算法的辨识精度最高，单一的人工鱼群算法精度次之，而粒子群算法精度最低，算法耗时排序则正好相反，综合来看，改进混合算法可以在最短时间内辨识出更加精确的全局最优解，参数辨识时效性和精准度均优于其他两种方法，尤其是与人工鱼群算法相比，改进混合算法的最优解相对误差减少了0.2个百分点，仅为鱼群算法相对误差的一半。图5 优化值进化曲线对比Fig.5 Comparison of fitness value curves人工鱼群算法及改进混合算法收敛速度对比如图4。由图4明显可以看到，鱼群算法在前5次迭代时目标值优化速度较快，但从第6次迭代开始目标值停留在1.48附近不再继续增加，直至达到最大迭代次数，相对

27、而言整体收敛速度较慢。而采用改进混合算法时，在第6次迭代后跳出鱼群算法的约束，转用并行模拟退火算法进行局部寻优，使得后期也达到了较高的收敛速度，快速收敛至第10次迭代，之后目标值稳定在1.475直至收敛结束，具有较优的搜索质量。多次使用改进混合算法进行参数辨识仿真实验可以发现，寻优迭代曲线基本一致，算法的随机性降低，收敛速度较为稳定，平均4.8次迭代即可达到前期阈值，平均11.5次可达到全局最优解。故而可设定最大迭代次数在15次以内，改进算法在300s以内即可得出互感器J-A模型参数辨识结果，准确度高且搜索速度快。4结束语文中提出了一种基于变步长人工鱼群算法和并行模拟退火算法的改进混合算法，用

28、于互感器J-A模型的5个关键参数辨识问题中。改进算法在搜索前期利用变步长人工鱼群算法接近全局最优解，搜索后期利用并行模拟退火算法进行进一步的局部搜索，以获取更为准确的全局最优解及其对应的参数辨识结果。建模仿真实验也验证了改进算法在较少的收敛次数下能获取更加精准的全局最优解，且计算速度不受传统算法随机性制约，稳定性较高。相比于传统的粒子群算法和人工鱼群算法，在耗时较短的情况下即可将精度提高近一倍之多，有效提升了互感器J-A模型参数辨识的效率和精度。该算法可应用于互感器特性研究的仿真实验中，能有效构建更为准确可靠的互感器模型，具有较高的实际应用价值，适用性较广。参 考 文 献1 李甘. 电流互感器

29、仿真技术研究及其分析评估平台开发D. 华中科技大学, 2014.Li Gan. Study on Current Transformer Simulation Technology and Development of Its Analysis & Evaluation PlatformD. Huazhong University of Science and Technology, 2004.2 郝晓亮. 基于智能优化算法的Jiles-Atherton磁滞模型参数计算研究D. 浙江师范大学, 2015.Hao Xiaoliang. Parameter Calculation Research

30、 of the Jiles-Atherton Hysteresis Model Using Intelligent Optimization Algorithms D. Zhejiang Normal University, 2015.3 李慧奇, 杨延菊, 邓聘, 等. 基于神经网络结合遗传算法的Jiles-Atherton磁滞模型参数计算J. 电网与清洁能源, 2012, (4): 19-22.Li Huiqi, Yang Yanju, Deng Pin, et al. Parameters Calculation of Jiles-Atherton Hysteresis Model Ba

31、sed on Combination of Genetic Algorithms and Neural NetworksJ. Power System and Clean Energy, 2012, (4): 19-22.4 杨延菊. 基于神经网络的Jiles-Atherton磁滞模型的实现D. 华北电力大学, 2012.Yang Yanju. Implement of Jiles-Atherton Hysteresis Model Based on Neural NetworksD. North China Electric Power University, 2012.5 李青峰. 基于混

32、合算法的J-A磁滞模型的实现D. 华北电力大学（河北）, 2010. Li Qingfeng. Achievement of J-A hysteresis model based on a Hybrid AlgorithmD. North China Electric Power University, 2010.6 王洋, 刘志珍. 基于蛙跳模糊算法的Jiles Atherton铁心磁滞模型参数确定J. 电工技术学报, 2017, (4): 154-161.Wang Yang, Liu Zhizhen. Determination of Jiles Atherton Core Hystere

33、sis Model Parameters Based on Fuzzy-Shuffled Frog Leaping AlgorithmJ. Transactions of China Electrotechnical Society, 2017, (4): 154-161.7 李贞, 李庆民, 李长云, 等. J-A磁化建模理论的质疑与修正方法研究J. 中国电机工程学报, 2011, (3): 124-131.Li Zhen, Li Qingmin, Li Changyun, et al. Queries on the J-A modelling theory of the magnetiza

34、tion processin ferromagnets and proposed correction methodJ. Proceedings of the CSEE, 2011, 31(3): 124-131. 8 J.R.Lucas, P.G.McLaren, and R.P.Jayasinghe. Improved simulation models for current and voltage transformers in relay studiesJ. IEEE Trans. On Power Delivery, 1992, 7(1): 152-159.9 Jiles D C, Atherton D L. Theory of ferromagnetic hysteresisJ. Journal of Applied Physics, 1984, 55(6): 2115-2120.10 Jiles D C, Atherton D L. Theory of ferromagnetic hysteresisJ. Journal of Magnetism and Magnetic Materials, 1986, 61(1