基于小波变换的小电流接地故障选线仿真研究

第２７卷第２期　电力科学与工程　Ｖｏ１．２７．Ｎｏ．２　４１　２０１１年２月　Ｅｌｅｃｔｒｉｃ　Ｐｏｗｅｒ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｆｅｂ．，２０１　１　基于小波变换的小电流接地故障选线仿真研究　窦新宇　，李春明　（１．唐山学院，河北唐山０６３０００；２．内蒙古工业大学，内蒙古呼和浩特０１００５１）　摘要：分析了小电流接地系统发生单相接地故障时的运行特点，提取故障暂态电流对其进行小波变换，　选出重构后小波模值极大且极性与其他线路相反的线路作为故障线路。在Ｍａｔｌａｂ　６．５环境下建立了相关　模型并进行了仿真验证；仿真结果表明：该方法不受故障条件的影响，可准确、可靠地实现配电网单相　接地故障选线，适用于中性点经消弧线圈接地或中性点不接地系统。　关键词：小电流接地；故障选线；小波变换；仿真　中圈分类号：ＴＭ７１１　文献标识码：Ａ　Ｏ　引　言　单相接地故障是小电流接地系统中最常见的故　　ｌ上上上　障，由于故障电流很小，因此，其故障选线的准确　工ＴｌＩｌａ　占　性一直是电力系统亟待解决的难题之一。随着对此　Ｅｂ　Ｅ　问题的深入研究，出现了各种故障选线方法，如零　Ｉ　序电流原理、首半波原理和五次谐波分量原理等。　上．．Ｌ．　Ｔ　Ｉｌ１－ｉ（　以上几种方法分别适用于不同的选线情况。　本文从发生单相接地故障时的暂态情况　出　一　Ｉ　发，在研究系统各线路零序电流的基础上，利用　｛　ｌ　ｃ　小波变换进行相关信息的提取；再利用小波对信　ｆ　号奇异性的检测准确定位故障时刻；最后利用小　ｊ＿上上　Ｔ　Ｔｑ－ｌ。　波模极大值原则实现故障选线。利用Ｍａｔｌａｂ软件　Ｉ　进行了大量仿真，表明这种方法是可行的。　ｌ小电流接地系统单相接地故障分析　１．１小电流接地系统模型　图１小电流接地系统单相接地　Ｆｉｇ．１　Ｓｉｎｇｌｅ－ｐｈａｓｅ—ｔｏ－ｇｒｏｕｎｄ　ｆａｕｌｔ　ｉｎ　ｎｅｕｔｒａｌ　图１给出了一个５条出线的小电流接地系统　ｎｏｎｅｆｆｅｃｔｉｖｅｌｙ　ｇｒｏｕｎｄｅｄ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｓｙｓｔｅｍ　单相接地故障示意图。按照一般方法，忽略其线　路电感、电阻以及对地电导。当开关Ｋ闭合时为　故障时，故障相对地电压降低，非故障相对地电　中性点经消弧线圈接地系统，打开时为中性点不　压升高。流过故障点的暂态零序电流由暂态电容　接地系统。图中Ｅ分别为等效的三相电源电动势；　电流和暂态电感电流叠加而成。暂态电感电流衰　Ｌ，Ｃ分别为消弧线圈的等效电感和各线路的单相　减很快，暂态电容电流衰减较慢。两者频率相差　对地分布电容。　悬殊，不能相互抵消，故在暂态初始过程，消弧　１．２暂态故障特征分析　线圈的补偿作用可以忽略不计，中性点不接地或　当中性点经消弧线圈接地系统发生单相接地　经消弧线圈接地系统的暂态零序电流分布相同　Ｊ。　收稿日期：２０１０—０６—０５。　作者简介：窦新字（１９８４一），男，助教，研究方向为电力系统的故障诊断，智能控制系统，Ｅ—ｍａｉｌ：ｄｏｕｘｉｎｙｕ＠　ｙａｈｏｏ．ｅｎ。　４２　电力科学与工程　此时，非故障线路母线处暂态零序电流为自身对　波变换模极大值的极性和大小判别出故障线路，　　地暂态电容电流，方向由母线流向线路；故障线　进而归纳选线判据如下　：（１）通过零序电压的突变来确定故障时刻。　路母线处暂态零序电流为所有非故障线路对地暂　态零序电流之和，方向由线路流向母线。　２基于小波的选线方法　２．１小波变换基本原理　通过对小电流接地系统单相接地故障时的零　序电压、零序电流进行奇异性检测，可以确定出　它们在故障后突变部分的极性和大小，比较其在　各条出线上的不同情况，由此识别出故障线路。　数学上将无限次可导的函数称为没有奇异性，　若函数在某处有间断点或某阶导数不连续，则称　其在此处有奇异点。奇异性检测就是要将信号的　奇异点识别出来并判断它的奇异性。通常用Ｌｉｐｓ—　ｃｈｉｔｚ指数来表达信号的奇异性　。由于小波变换　极大值在多尺度上的表现与Ｌｉｐｓｃｈｉｔｚ指数之问存　在对应关系，这为通过小波变换检测信号奇异点　提供了依据。即小波变换后的模极大值能够反映　接地故障的某些特性，故本文利用此理论实现故　障选线。　２．２小波函数和分解层数的选取　为了减小频谱的泄漏和混叠，要求小波函数　具有较好的频域特性。其中ｄｂＮ小波是工程上应　用较多的小波函数，该小波系特点为随Ⅳ的增大，　时域支集变长，时间局部性变差；同时，正则性　增加，频域局部性变好。综合考虑在时域与频域　的需要，通过采用多种小波进行仿真，本文使用　ｄｂｌ０小波可以得到较为理想的结果。　对于分解层数的选择，层数过多则频带宽度　过小，频带对应的采样点数过少，使得判据灵敏　度降低；反之层数过少则频带宽度过大，信息量　增加，可能引入更多的干扰成分，降低可靠性。　因此，一般选择分解到６层　’。　。　２．３故障选线原理及判据　小电流接地系统发生单相接地故障时产生一　个零序突变，所有非故障线路零序电流突变的极　性相同，故障线路的极性与非故障线路的极性相　反，且故障线路零序电流突变的幅值等于所有非　故障线路零序电流突变幅值之和。利用小波变换　对奇异值检测的特性，通过比较暂态零序电流小　（２）在故障时刻附近，若每一条线路小波变　换模极大值的极性与其他线路相反，且幅值最大，　则该线路为故障线路。　（３）在故障时刻附近，若所有线路的小波变　换模极大值的极性相同，则为母线故障。　３　系统仿真　３．１系统仿真模型　该系统为５条出线的１０　ｋＶ系统，包括三相　电源、升压变压器、输电线（架空）、降压变压　器、用户负载。升压变压器低压侧中性点通过开　关和消弧线圈相连。开关打开时为不接地系统，　闭合时为经消弧线圈接地系统。　（１）线路参数：电力输电线的参数严格来讲　是均匀分布的。由于配电网中存在电力负荷的频　繁切换使三相参数不可能平衡，采用集中参数模　型存在着较大误差。因此采用分布参数模型。对　架空线路进行仿真，其参数选定如表ｌ。　表１架空线参数　Ｔａｂ．１　Ｔｒａｎｓｍｉｓｓｉｏｎ　ｌｌｎｅ　ｐａｒａｍｅｔｅｒｓ　（２）变压器参数：升压变压器为１１０／１０　ｋＶ，　ＹＹ联接，降压变压器为１０／６　ｋＶ，Ｄ１１Ｙ联接。　（３）线路长度：Ｌ　＝５　ｋｍ，　＝１０　ｋｍ，　＝　１４　ｋｍ，Ｌ４＝１６　ｋｍ，Ｌ５＝ｌ８　ｋｍ。　（４）消弧线圈参数：在仿真消弧线圈接地系　统时，系统按过补偿１０％整定。根据线路参数及　长度可以计算出Ｒ　．＝６　ｎ，Ｌ＝０．７１９　Ｈ。　（５）仿真时间为０．１　ｓ，微分方程解算器选择　变步长ｏｄｅｌ５，在含架空线的暂态电容电流自由振　荡的频率一般集中在３００—３　０００　Ｈｚ，本课题选取　的采样频率为６　０００　Ｈｚ。　３．２仿真算例与结果分析　假设线路　的Ａ相发生接地故障，故障时刻　为Ｔ＝０．０２５　Ｓ，距离母线８　ｋｍ，接地电阻为１００　Ｑ。　图２为系统零序电压波形和各个出线零序电流波　形，图３是系统三相电压，图４为用ｄｂｌ０小波对　第２期　窦新字，等基于小波变换的小电流接地故障选线仿真研究　４３　各线路零序电流进行第６层重构后的局部波形。　从图２中可以看出，系统零序电压在０．０２５　Ｓ　附近发生突变，说明系统发生单相接地故障；通过　以上分析，从图３中可以知道故障发生在Ａ相；最　后通过图４可知线路　上的模值在Ｔ＝０．０２５　Ｓ附　近均大于其他线路的值，且方向相反，按照前面介　蚰如加ｍ　０　表２仿真结果　加如们　Ｏ　０　Ｔａｂ．２　Ｓｉｍｕｌａｔｉｏｎ　ｒｅｓｕｌｔｓ　绍的选线方法可知，线路厶的Ａ相发生故障。　表２是仿真的几种情况的结果概括。　＞　寸ｅ　×　Ｈ皇　恤　恃　建　性　盼　图２系统零序电压波形和各个出线零序电流　Ｆｉｇ．２　Ｖｏｌｔａｇｅ　ａｎｄ　ｃｕｒｒｅｎｔ　ｏｆ　ｚｅｒｏ　ｓｅｑｕｅｎｃｅ　＞　可ｅ　×　０　ｆ／ｓ　图３系统三相电压　Ｆｉｇ．３　Ｔｈｒｅｅ－ｐｈａｓｅ　ｖｏｌｔａｇｅ　ｏｆ　ｓｙｓｔｅｍ　２　４　２　４５　２　５　２　５５　采样点　图４小波重构　Ｆｉｇ．４　Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ　ｆｏｒ　ｗａｖｅｌｅｔ　４结论　利用小波变换识别暂态故障线路的方法是可　行的。仿真结果表明，在经高阻接地故障时，仍　能得出正确的选线结果。因而，利用小波变换进　行单相接地故障选线具有很强的抗过渡电阻能力　且不受故障初相角的影响，该方法适用于中性点　经消弧线圈接地或中性点不接地系统。　参考文献：　［１］薛永端，冯祖仁，徐丙垠．中性点非直接接地电网单　相接地故障暂态特征分析［Ｊ］．西安交通大学学报，　２００４，３８（２）：１９５—１９９．　Ｘｕｅ　Ｙｏｎｇｄｕａｎ，Ｆｅｎｇ　Ｚｕｒｅｎ，Ｘｕ　Ｂｉｎｇｙｉｎ．Ａｎａｌｙｓｉｓ　ｏｆ　ｔｒａｎｓｉｅｎｔ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｓｉｇｎａｌｓ　ｐｈａｓｅ　ｅａｃｈ　ｆａｕｌｔ　ｉｎ　ｎｏｎ—ｓｏｌｉｄｌｙ　ｅａｒｔｈｅｄ　ｎｅｔｗｏｒｋ『＿Ｉ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｘｉ’ａｎ　Ｊｉａｏ—　ｔｏｎｇ　Ｕｎｉｖｅｒｓｉｔｙ，２００４，３８（２）：１９５—１９９．　［２］Ｚｈａｎｇ　Ｘｉｎｈｕｉ，Ｘｕｅ　Ｙｏｎｇｄｕａｎ．Ｓｉｍｕｌａｔｉｏｎ　ａｎａｌｙｓｉｓ　ｏｆ　ｚｅｒｏ　ｍｏｄａｌ　ｔｒａｎｓｉｅｎｔ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｓｉｎｇｌｅ—ｐｈａｓｅ　ｅａｒｔｈ　ｆａｕｌｔ　［Ａ］．２００６　Ｃｈｉｎａ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｃｏｎｆｅｒｅｎｃｅ　ｏｎ　Ｅｌｅｃｔｉｒｃｉｔｙ　Ｄｉｓｔｉｒｂｕｔｉｏｎ［ｃ］．ＩＥＴ　Ｃｏｎｆｅｒｅｎｃｅ　Ｐｕｂｌｉｃａｔｉｏｎｓ，２００６．　２１８—２２３．　［３］唐捷．小电流接地系统单相接地故障选线新方法　［Ｊ］．高电压技术，２００７，３３（１２）：１４３—１４７．　Ｔａｎｇ　Ｊｉｅ．Ｎｅｗ　ｍｅｔｈｏｄ　ｆｏｒ　ｌｉｎｅ　ｆａｕｌｔ．ｄｅｔｅｃｔｉｏｎ　ｕｎｄｅｒ　ｓｉｎ．　ｇｌｅ－－ｐｈａｓｅ・－ｔｏ－－ｅａｒｔｈ　ｆａｕｌｔ　ｏｆ　ｌｏｗ　ｃｕｒｒｅｎｔ　ｎｅｕｔｒａｌ　ｇｒｏｕｎｄｉｎｇ　ｐｏｗｅｒ　ｓｙｓｔｅｍ［Ｊ］．Ｈｊｇｈ　Ｖｏｌｔａｇｅ　Ｅｎｇｉｎｅｅｒｉｎｇ，２００７，３３　（１２）：１４３—１４７．　［４］刘明才．小波分析及其应用［Ｍ］．北京：清华大学　出版社，２００５．　［５］Ｃｈｕｌ　Ｈｗａｎ　Ｋｉｍ，Ａｇｇａｒｗａｌ　Ｒ．Ｗａｖｅｌｅｔ　ｔｒａｎｓｆｏｒｍｓ　ｉｎ　ｐｏｗ—　ｅｒ　ｓｙｓｔｅｍｓ．ｐａｒｔ　２：Ｅｘａｍｐｌｅｓ　ｏｆ　ａｐｐｌｉｃａｔｉｏｎ　ｔｏ　ａｃｔｕａｌ　ｐｏｗ—　ｅｒ　ｓｙｓｔｅｍ　ｔｒａｎｓｉｅｎｔｓ［Ｊ］．Ｐｏｗｅｒ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｊｏｕｒｎａｌ，　２００１，１５（４）：１９３—２０２．　（下转第７１页）　第２期　李义成超临界锅炉受热面氧化皮的检测分析与对策　Ｄｅｔｅｃｔｉｏｎ　Ａｎａｌｙｓｉｓ　ａｎｄ　Ｃｏｕｎｔｅｒｍｅａｓｕｒｅｓ　ｏｎ　Ｓｕｐｅｒｃｒｉｔｉｃａｌ　Ｂｏｉｌｅｒ　Ｈｅａｔｉｎｇ　Ｓｕｒｆａｃｅ　Ｏｘｉｄｅ　Ｓｋｉｎ　Ｌｉ　Ｙｉｃｈｅｎｇ　ｆ　Ａｎｈｕｉ　Ｈｕａｉｎａｎ　Ｐｉｎｇｗｅｉ　Ｐｏｗｅｒ　Ｍａｉｎｔｅｎａｎｃｅ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｃｏ．，Ｌｔｄ．，Ｈｕａｉｎａｎ　２３２０８９，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：ｉｎ　ｔｈｉｓ　ｐａｐｅｒ，Ｐｉｎｇ　Ｐｏｗｅｒ　Ｐｌａｎｔ　Ｎｏ．３　ｂｏｉｌｅｒ　ｈｅａｔｉｎｇ　ｓｕｒｆａｃｅ　ｏｘｉｄｅ　ｄｅｔｅｃｔｉｏｎ　ｍｅｔｈｏｄｓ　ｗｅｒｅ　ｉｎｔｒｏｄｕｃｅｄ，　ｆｏｃｕｓｅｄ　ｏｎ　ａｎａｌｙｚｉｎｇ　ｔｈｅ　ｃａｕｓｅｓ　ｏｆ　ｏｘｉｄａｔｉｏｎ　ｏｆ　ｓｋｉｎ，ａｎｄ　ｐｕｔ　ｆｏｒｗａｒｄ　ｓｐｅｃｉｆｉｃ　ｃｏｎｔｒｏｌ　ｍｅａｓｕｒｅｓ，ｉｎ　ｏｒｄｅｒ　ｔｏ　ｓｏｌｖｅ　ｔｈｅ　ｓｕｐｅｒ—ｃｒｉｔｉｃａｌ　ｕｎｉｔ　ｆｕｒｎａｃｅ　ｔｕｂｅ　ｏｘｉｄａｔｉｏｎ　ｏｆ　ｓｋｉｎ　ｐｅｅｌｉｎｇ　ｏｆｆ　ｃａｕｓｅｄ　ｂｙ　ｔｈｅ　ｂｏｉｌｅｒ，ｔｈｉｓ　ｐａｐｅｒ　ｐｒｏｖｉｄｅｓ　ａ　ｇｏｏｄ　ｗａｙ．　Ｋｅｙ　ｗｏｒｄｓ：ｏｘｉｄｅ　ｓｋｉｎ；ｄｅｔｅｃｔｉｏｎ；ａｎａｌｙｓｉｓ；ｃｏｕｎｔｅｒｍｅａｓｕｒｅｓ　（上接第４３页）　［６］Ｃｈａｎｄａ　Ｄ，Ｋｉｓｈｏｒｅ　Ｎ　Ｋ，Ｓｉｎｈａ　Ａ　Ｋ．Ａ　ｗａｖｅｌｅｔ　ｍｕｌｔｉ—　ｇｒｏｕｎｄ　ｆａｕｌｔ　ｆｅｅｄｅｒ　ｉｎ　ｎｅｕｔｒａｌ　ｎｏｎｅｆｆｅｃｔｉｖｅｌｙ　ｇｒｏｕｎｄｅｄ　ｄｉｓ—　ｒｅｓｏｌｕｔｉｏｎ　ａｎａｌｙｓｉｓ　ｆｏｒ　ｌｏｃａｔｉｏｎ　ｏｆ　ｆａｕｌｔｓ　ｏｎ　ｔｒａｎｓｍｉｓｓｉｏｎ　ｔｒｉｂｕｔｉｏｎ　ｓｙｓｔｅｍ　ｕｓｉｎｇ　ｗａｖｅｌｅｔ　ｔｒａｎｓｆｏｒｍ［Ｊ］．ＩＥＥＥ　ｌｉｎｅｓ［Ｊ］．Ｅｌｅｃｔｒｉｃｌａ　Ｐｏｗｅｒ　ａｎｄ　Ｅｎｅｒｇｙ　Ｓｙｓｔｅｍｓ，２００３，　Ｔｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｄｅｌｉｖｅｒｙ，２００８，２３（４）：１８２９—　（２５）：５９—６９．　１８３７．　［７］Ｄｏｎｇ　Ｘｉｎｚｈｏｕ，Ｓｈｉ　Ｓｈｅｎｘｉｎｇ．Ｉｄｅｎｔｉｆｙｉｎｇ　ｓｉｎｇｌｅ—ｐｈａｓｅ—ｔｏ—　Ｓｉｍｕｌａｔｉｏｎ　Ｒｅｓｅａｒｃｈ　ｏｎ　Ｓｉｎｇｌｅ－ｐｈａｓｅ－ｔｏ—ｇｒｏｕｎｄ　Ｆａｕｌｔ　Ｆｅｅｄｅｒ　ｉｎ　Ｎｅｕｔｒａｌ　Ｎｏｎｅｆｆｅｃｔｉｖｅｌｙ　Ｇｒｏｕｎｄｅｄ　Ｄｉｓｔｒｉｂｕｔｉｏｎ　Ｓｙｓｔｅｍ　Ｕｓｉｎｇ　Ｗａｖｅｌｅｔ　Ｔｒａｎｓｆｏｒｍ　Ｄｏｕ　Ｘｉｎｙｕ　，Ｌｉ　Ｃｈｕｎｍｉｎｇ　（１．Ｔａｎｇｓｈａｎ　Ｃｏｌｌｅｇｅ，Ｔａｎｇｓｈａｎ　０６３０００，Ｃｈｉｎａ；　２．Ｉｎｎｅｒ　Ｍｏｎｇｏｌｉａ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｔｅｃｈｎｏｌｏｇｙ，Ｈｏｈｈｏｔ　０１００５　１，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ａ　ｓｃｈｅｍｅ　ｏｆ　ｓｉｎｇｌｅ—ｐｈａｓｅ・ｔｏ—ｇｒｏｕｎｄ　ｆａｕｌｔ　ｆｅｅｄｅｒ　ｉｄｅｎｔｉｆｉｃａｔｉｏｎ　ｉｎ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｎｅｔｗｏｒｋｓ　ｗｉｔｈ　ｔｈｅ　ａｐｐｌｉｃａ—　ｔｉｏｎ　ｏｆ　ａ　ｗａｖｅｌｅｔ　ｔｒａｎｓｆｏｒｍ　ｔｅｃｈｎｉｑｕｅ　ｉｓ　ｐｒｅｓｅｎｔｅｄ　ｉｎ　ｔｈｉｓ　ｐａｐｅｒ．Ｔｈｅ　ｌｏｃａｌ　ｍｏｄｕｌｕｓ　ｍａｘｉｍａ　ｏｆ　ｔｈｅ　ｗａｖｅｌｅｔ　ｔｒａｎｓｆｏｒｍ　ａｒｅ　ｅｘｔｒａｃｔｅｄ　ｔｏ　ｄｅｔｅｒｍｉｎｅ　ｔｈｅ　ｔｉｍｅ　ｏｆ　ｔｈｅ　ｓｉｎｇｌｅ—ｐｈａｓｅ—ｔｏ－ｇｒｏｕｎｄ　ｆａｕｌｔ．Ｔｈｅ　ｗａｖｅｌｅｔ　ｔｒａｎｓｆｏｒｍｓ　ｏｎ　ａｌｌ　ｆｅｅｄｅｒｓ　ａｔ　ｔｈｅ　ｔｉｍｅ　ａｒｅ　ｃｏｍｐａｒｅｄ　ｉｎ　ｍａｇｎｉｔｕｄｅ　ａｎｄ　ｐｏｌａｒｉｔｙ　ｗｉｔｈ　ｅａｃｈ　ｏｔｈｅｒ　ｔｏ　ｉｄｅｎｔｉｆｙ　ｔｈｅ　ｆａｕｌｔｅｄ　ｆｅｅｄｅｒ．Ｒｅｌａｔｉｖｅ　ｓｙｓｔｅｍ　ｍｏｄｅｌｓ　ｃ０ｕｌｄ　ｂｅ　ｂｕｉｌｔ　ａｎｄ　ｓｉｍｕｌａｔｅｄ　ｂａｓｅｄ　ｏｎ　ｔｈｅ　ｅｎｖｉｒｏｎｍｅｎｔ　ｏｆ　ＰＳＢ　ｉｎ　Ｍａｔｌａｂ６．５．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｔｈｅ　ｍｅｔｈｏｄ　ｃａｎ　ｃｈｏｏｓｅ　ｔｈｅ　ｆａｕｌｔｙ　ｌｉｎｅ　ｅｘａｃｔｌｙ　ｕｎｄｅｒ　ａｌｌ　ｋｉｎｄｓ　ｏｆ　ｆａｕｌｔ　ｃｏｎｄｉｔｉｏｎｓ．　Ｋｅｙ　ｗｏｒｄｓ：ｌｏｗ—ｃｕｒｒｅｎｔ　ｇｒｏｕｎｄ　ｓｙｓｔｅｍ；ｆａｕｌｔｙ　ｌｉｎｅ　ｄｅｔｅｃｔｉｏｎ；ｗａｖｅｌｅｔ　ｔｒａｎｓｆｏｒｍ；ｓｉｍｕｌａｔｉｏｎ