Fluent模型在河道堤防渗透数值计算中的适用性和精度对比分析

第０４期（总第４１９期）　【文章编号】１００９—２８４６（２０１７）０４—００２３－０５　吉　林　水利　２０１７年０４月　Ｆｌｕｅｎｔ模型在河道堤防渗透数值计算中的　适用性和精度对比分析　刘亿军。曲春阳　（辽宁天阳工程技术咨询服务有限公司，辽宁　沈阳　１１００００）　［摘要】本文结合渗流物理模型试验，分析了Ｆｌｕｅｎｔ模型在河道堤防渗透计算的适用性，并对比分析了Ｆｌｕｅｎｔ　模型和传统渗流计算模型在河道堤防渗透计算的精度。研究结果表明：Ｆｌｕｅｎｔ模型在河道堤防渗透计算具有较　好的适用性．计算的渗透系数和渗流试验分析的堤防渗透系数之间的相对误差值小于１５％；Ｆｌｕｅｎｔ模型由于采　用有限体积方法对非稳定渗透流方程进行离散，计算精度好于传统渗流计算方程，计算误差减少１２．２％，相关　系数提高０．２８。研究成果对于设计河道堤防的渗透合理计算提供方法参考。　【关键词】渗流物理模型试验；Ｆｌｕｅｎｔ模型；传统渗流方程；非稳定渗透流；模型适用性和精度对比分析　［中图分类号］Ｔｖ　２１４　［文献标识码］Ｂ　２　Ｆｌｕｅｎｔ模型计算原理　河道堤防设计需要计算堤防的渗透系数，从　而对设计河道堤防的稳定性进行评价。河道堤防　渗透系数的准确计算对于河道堤防设计的合理性　至关重要．因为河道堤防渗透系数和河道堤防的　Ｆｌｕｅｎｔ模型将区域设定为多孔传输介质的水　流流体，结合达西定律计算河道堤防的渗透孔隙　率以及渗透水流的阻力系数，其渗透计算方程为：　：ＫＪ＝一Ｋ　Ⅱｓ　一　稳定性密切相关，为此许多学者开展对河道堤防　渗透计算的研究，并取得一定的研究成果【ｔ－５］，但是　ｃＺｓ　（１）　Ｆｌｕｅｎｔ模型对传统达西计算定律进行了简化　计算，计算方程为：　如　一　传统渗流模型大都采用稳定流方程进行计算，而　河道堤防大都为非稳定渗透流．因此需要进行非　（２）　稳定渗透流计算，从而提高河道堤防设计渗透系　数计算的合理性。当前，Ｆｌｕｅｎｔ模型采用有限体积　方法对非稳定水流方程进行离散．在非稳定水动　方程（１）和（２）中，Ｋ为Ｆｌｕｅｎｔ模型计算的渗　透系数值（ｃｍ／ｓ）；Ｉ，为渗透水流的水力学比降（％）；　ｈ为Ｆｌｕｅｎｔ模型计算的渗透流水头（ｃｍ）；Ｓ为渗透　力学计算中得到具体应用［６－　Ｏｌ，但在河道堤防非稳　定渗透流数值计算中的应用较少，为此本文结合　渗流物理模型试验。通过对比观测试验分析Ｆｌｕｅｎｔ　模型在河道堤防渗透数值计算中的适用性．并对　流水力半径（ｃｍ）；ｒ为渗透流的容重系数（ｋｇ／ｍ　）。　结合方程（１）和（２）可以推算渗透水流的孔隙　率，计算方程为：　＝比Ｆｌｕｅｎｔ模型和传统渗流模型在河道堤防渗透的　计算精度。研究成果对于河道堤防非稳定渗透流　ＫＪ＝－Ｋ单一坠孕：一　孕　ａｓ　ｙ　ａｓ　ａｓ　（３）　（４）　计算提供方法参考。　由方程：　＝一　则Ｆｌｕｅｎｔ模型的水流阻力系数计算方程计算　的表达式：　［收稿日期】２０１７一Ｏ１—２３　［作者简介］刘亿军（１９８６一），男，山西清徐人，助工，现从事水利工程设计工作。　２３—　一吉林水利　Ｆｌｕｅｎｔ模型在河道堤防渗透数值计算中的适用性和精度对比分析刘亿军等２０１７年０４月　一：堕：　：　ｖ　Ｑ　Ｑ　（５）　、　在方程（６）一（８）Ｋ和　的变量含义同方程　９８００　（２）；ｚ为渗透水流在垂向上的方向参数。　方程（５）中，　为计算的渗透系数值（ｃｎｄｓ）；ｏｒ　和　为河道堤防的地质参数。　３模型适用性及精度对比分析　３．１　渗流物理模型试验概况　Ｆｌｕｅｎｔ模型结合有限体积方法对非稳定渗透　流进行求解计算，求解计算方程为：　Ｏ（Ｋｘｘ　Ｏｘ）＋　（＝Ｉｖ　（　ｄｖ　）＋　Ｏ（　Ｏｚ渗透物理模型试验河道长度为４５　ｍ，由三个　部分组成．分别为过度连接段、动床河段以及稳水　（６）　ｚ　河段，其中动水段试验河长度为３０ｍ，在试验河槽　模型假定河道堤防各渗透方向的参数为常数　的上、下游各没置５ｍ的过渡连接河段，过渡连接　值，则方程（６）可转化为以下方程：　河段与动水河床段之间采用长度为２０ｃｍ的粗砂　Ｋｘ　＋研　＋Ｋｚ　：０　（７）　段进行连接，整个试验河槽平面布置图见图１．渗　ａ　‘　ａｙ　ａ　‘　透物理模型河道堤防布设场景图见图２。渗透物理　若模型各渗透方向为同向转换时．其计算方　模型试验布设５个试验断面，结合放水试验．进行　程变化为：　５个试验断面渗透流试验，测定各断面的渗透系　数　＋　＋　：０　（８）　ａ　ａｖ　５咖…十　…５ｏ０ｏ……一ｒ…　一………　ｒ　、一　，Ｉ：２，　：　。　【　一　Ｉ　２　Ｉ＝　｝１Ｉ【ｍ０　。　Ｉ一　…　一～　…　唧　…一　…　…一…　～．．３　三二　１、。　～＿　５　。　：２　。　０　：００　一　¨ｌ　冀。。　　●■■■●＿■■＿＿●＿■■■＿■■■●●■＿＿＿＿＿●■＿＿＿＿■■■＿＿＿＿＿＿＿■一　图Ｉ　模式试验河段平面布置图　表ｌ　各试验断面物理属性参数表　图２渗透试验河道堤防场景图　３．２模型的适用性分析　为对Ｆｌｕｅｎｔ模型进行河道堤防渗透计算适Ｈｊ　性分析，结合渗透物理模型试验对模型进行试验　分析，各试验断面物理参数见表１。进行了３组放　一２４一　吉林水利　表２　Ｆｌｕｅｎｔ模型在河道堤防渗透数值计算中的适用性和精度对比分析刘亿军等基于Ｆｌｕｅｎｔ模型的渗透系数　２０１７年ｏ４月　价结果见表２和图３。　计算精度分析结果　组序　．　，　表１为各试验断面地质参数，从表１中可以　看出。渗透物理模型试验各试验断面的粘聚力在　（％）（１０－Ｓｅｒｅ／ｓ）　试验测定渗透系数计算渗透系数值相对误差绝对误差　（１Ｏ　ＴＩｌ，ｓ）　（１０－ＳｃＩＩｌ，ｓ）　２．６９－－３９．５ｋＰａ之间变化，其内摩擦角在８．２。一　２３．１ｏ之间变化，其物理参数的变化主要受气土质　和土层厚度有关。表２为采用Ｆｌｕｅｎｔ模型计算的　各试验断面在不同放水试验下的渗透系数和试验　测定的渗透系数之间的精度分析结果，从表２中　可以看出，采用Ｆｌｕｅｎｔ模型计算的渗透系数和试　验测定的渗透系数之间相对误差在一３．８％一１４．５％　之间，均小于１５％，在水利工程设计规范中要求计　算规范小于２０％，且两种渗透系数之间的绝对误　差值在０．１ｘｌ０　一１．２ｘ１０　ｅｍ／ｓ之间，符合水利工程　设计规范中的精度范围要求。图３为Ｆｌｕｅｎｔ模型　计算的渗透系数和试验测定的渗透系数过程拟合　度和相关性分析结果，从图中可以看出，在３组放　水试验中，Ｆｌｕｅｎｔ模型计算各试验断面的渗透系数　和试验测定的渗透系数过程吻合度较高，且相关　水试验，并测定放水试验５个试验断面１５组渗透　系数达到０．７６５１，相关性较高。综上，Ｆｌｕｅｎｔ模型在　河道堤防的渗透计算具有较好的精度。适用于河　道堤防的渗透计算。　系数值。结合对比试验测定的渗透系数分析Ｆｌｕｅｎｔ　模型在河道堤防渗透系数计算的适用性，精度评　（ａ）过程拟合分析　（ｂ）相关性分析　图３　Ｆｌｕｅｎｔ模型计算的渗透系数和试验测定的渗透系数过程拟合和相关性分析　３．３不同模型精度对比分析　果，从表中可以看出，Ｆｌｕｅｎｔ模型计算的渗透系数　和试验测定的渗透系数结果好于传统渗透模型。　相比于传统渗透模型．Ｆｌｕｅｎｔ模型计算的渗透系数　和试验测定的渗透系数之间的相对误差均值减少　１２．２％。图４为不同模型计算的渗透系数和试验测　定的渗透系数相关分析结果，从图中可以看出，　在Ｆｌｕｅｎｔ模型渗透数值计算适用性分析的基　础上，结合２组渗透放水试验，基于放水试验５个　试验断面测定的１０组渗透系数数据，对比Ｆｌｕｅｎｔ　模型和传统渗透计算模型的渗透系数计算精度，　对比分析结果见表３和图４。　表３为不同模型的渗透系数计算精度对比结　Ｆｌｕｅｎｔ模型计算的渗透系数和试验测定的渗透系　一２５—　吉林水利　Ｆｌｕｅｎｔ模型在河道堤防渗透数值计算中的适用性和精度对比分析刘亿军等２０１７年０４月　数之间的相关系数为０．７８５２，而传统渗透模型计　算的渗透系数和试验测定的渗透系数之间的相关　系数为０．５０４２，相比于传统渗透模型，Ｆｌｕｅｎｔ模型　采用有限体积方法对非稳定渗流方程进行求解。　而传统渗透模型假定渗透水流为稳定的水流。对　方程进行的简化计算，因此计算精度要低于Ｆｌｕｅｎｔ　模型。　在相关系数上提高０．２８。这主要因为Ｆｌｕｅｎｔ模型　表３　不同模型渗透系数计算精度对比结果　（ａ）Ｆｌｕｅｎｔ模型　（ｂ）传统渗透模型　圈４不同模型计算的渗透系数和试验测定的渗透系数相关分析结果　散求解问题，河道堤防渗透计算精度好于传统渗　４　结论　本文结合渗透物理模型试验，采用对比观测　方式分析Ｆｌｕｅｎｔ模型在河道堤防渗透计算的适用　性．并对比分析Ｆｌｕｅｎｔ模型和传统渗透模型在河　透模型，可进一步推广应用于河道堤防设计中的　渗透计算。口　参考文献：　［１】范朝晖，宋建忠．辽西河道堤防护坡方法分析［Ｊ】．吉林水利，２００７，　（０６）：１５—１７．　【２】张德荣，陈小平．土堤坝渗透破坏及抢险措施［Ｊ】．吉林水利，２００３，　（０６）：３１－３２．　道堤防渗透计算精度，研究取得以下结论：　１）Ｆｌｕｅｎｔ模型可用于河道堤防的渗透计算，模　型计算值和试验测定值之间的相对误差小于　［３］韩璐．饮马河堤防防渗方案优选设计【Ｊ】．吉林水利，２０１６，（１２）：　２９—３０＋３３．　１５％，符合水利工程设计规范的精度要求；　２）Ｆｌｕｅｎｔ模型可解决非稳定渗透流的方程离　１４］杨扬，徐伟，李再兴，孙健．基于ＳＥＥＰ／Ｗ程序的堤防截渗墙防渗　性能研究『Ｊ］．黄河水利职业技术学院学报，２０１６，（Ｏ１）：１－５．　【５】冷爱国．垂直防渗在堤防防渗加固中的比较应用【Ｄ］．合肥工业大　一２６一　吉林水利　Ｆ１ｕｅｎｔ模型在河道堤防渗透数值计算中的适用性和精度对比分析刘亿军等２０１７年０４月　学。２００５．　大学．２０１１．　【６Ｊ刘宝星．圃中环沉沙排沙过滤池水工模型试验研究与数值分析　［９１￣ｌ鑫煜．南排河中心渔港工程潮流泥沙数学模型【Ｄ】．天津大学，　［Ｄ】．新疆农业大学，２０１３．　２０ｏ８．　［７】张瑜．基于Ｆｌｕｅｎｔ的堤防工程渗流场数值分析【Ｄ】．太原理工大　【１ｏ］刘宝星．圆中环沉沙排沙过滤池水工模型试验研究与数值分析　学，２０１３．　［Ｄ】．新疆农业大学，２０１３．　【８】林长强．基于ＦＬＵＥＮＴ的土石坝逐渐溃坝水流模拟［Ｄ］．华中科技　Ａｐｐｌｉｃａｂｉｌｉｔｙ　ａｎｄ　Ｃｏｎｔｒａｓｔ　Ａｎａｌｙｓｉｓ　ｏｆ　Ｆｌｕｅｎｔ　Ｍｏｄｅｌ　ｉｎ　Ｎｕｍｅｒｉｃａｌ　Ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　Ｒｉｖｅｒ　Ｄｙｋｅ　Ｉｎｆｉｌｔｒａｔｉｏｎ　Ｌｉｕ　Ｙｉｊｕｎ，Ｑｕ　Ｃｈｕｎｙａｎｇ　Ａｂｓｔｒａｃｔ：Ｉｎ　ｔｈｉｓ　ｐａｐｅｒ．ｔｈｅ　ａｐｐｌｉｃａｂｉｌｉｔｙ　ｏｆ　Ｆｌｕｅｎｔ　ｍｏｄｅｌ　ｉｎ　ｓｅｅｐａｇｅ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｆｉｖｅｒ　ｅｍｂａｎｋｍｅｎｔ　ｉｓ　ａｎａ’　ｌｙｚｅｄ　ｂｙ　ｔｈｅ　ｐｈｙｓｉｃａｌ　ｍｏｄｅｌ　ｔｅｓｔ，ａｎｄ　ｔｈｅ　ａｃｃｕｒａｃｙ　ｏｆ　Ｆｌｕｅｎｔ　ｍｏｄｅｌ　ａｎｄ　ｔｒａｄｉｔｉｏｎａｌ　ｓｅｅｐａｇｅ　ｃａｌｃｕｌａｔｉｏｎ　ｍｏｄｅｌ　ｉｎ　ｔｈｅ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｉｎｆｉｌｔｒａｔｉｏｎ　ｏｆ　ｒｉｖｅｒ　ｅｍｂａｎｋｍｅｎｔ　ａｒｅ　ａｎａｌｙｚｅｄ　ａｎｄ　ｃｏｍｐａｒｅｄ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　Ｆｌｕｅｎｔ　ｍｏｄｅｌ　ｈａｓ　ｇｏｏｄ　ａｐｐｌｉｃａｂｉｌｉｔｙ　ｉｎ　ｓｅｅｐａｇｅ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｆｉｖｅｒ　ｅｍｂａｎｋｍｅｎｔ，ａｎｄ　ｔｈｅ　ｒｅｌａｔｉｖｅ　ｅｒｒｏｒ　ｂｅｔｗｅｅｎ　ｔｈｅ　ｃａｌｃｕｌａｔｅｄ　ｐｅｒｍｅａｂｉｌｉｔｙ　ｃｏｅｆｆｉｃｉｅｎｔ　ａｎｄ　ｓｅｅｐａｇｅ　ｃｏｅｆｉｆｃｉｅｎｔ　ｉｓ　ｌｅｓｓ　ｔｈａｎ　１５％．Ｆｌｕｅｎｔ　ｍｏｄｅｌ　ｉｓ　ｍｏｒｅ　ｓｕｉｔｂａｌｅ　ｆｏｒ　ｕｎｓｔｅａｄｙ　ｆｌｏｗ　ｄｕｅ　ｔｏ　ｔｈｅ　ｆｉｎｉｔｅ　ｖｏｌｕｍｅ　ｍｅｔｈｏｄ　Ｔｈｅ　ｃａｌｃｕｌａｔｉｏｎ　ａｃｃｕｒａｃｙ　ｉｓ　ｂｅｔｔｅｒ　ｔｈａｎ　ｔｈａｔ　ｏｆ　ｔｈｅ　ｔｒａｄｉｔｉｏｎａｌ　ｓｅｅｐａｇｅ　ｃａｌｃｕｌａｔｉｏｎ　ｅｑｕａｔｉｏｎ．Ｔｈｅ　ｃａｌｃｕｌａｔｉｏｎ　ｅｒｒｏｒ　ｉｓ　ｒｅｄｕｃｅｄ　ｂｙ　１２．２％．ａｎｄ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｃｏｅｆｆｉｃｉｅｎｔ　ｉｓ　ｉｎ－　ｃｒｅａｓｅｄ　ｂｖ　０．２８．ｈＴｅ　ｒｅｓｅａｒｃｈ　ｒｅｓｕｌｔｓ　ｃａｎ　ｐｒｏｖｉｄｅ　ａ　ｒｅｆｅｒｅｎｃｅ　ｆｏｒ　ｔｈｅ　ｄｅｓｉｇｎ　ｏｆ　ｉｒｖｅｒ　ｄｉｋｅ　ｉｎｆｉｌｔｒａｔｉｏｎ　ｒｅａｓｏｎａｂｌｅ　ｃａｌｃｕｌａｔｉｏｎ．　Ｋｅｙｗｏｒｄｓ：ｓｅｅｐａｇｅ　ｐｈ）ｒｓｉｃａｌ　ｍｏｄｅｌ　ｔｅｓｔ；Ｆｌｕｅｎｔ　ｍｏｄｅｌ；ｔｒａｄｉｔｉｏｎａｌ　ｓｅｅｐａｇｅ　ｅｑｕａｔｉｏｎ；ｕｎｓｔｅａｄｙ　ｓｅｅｐａｇｅ　ｆｌｏｗ；　ｍｏｄｅｌ　ａｐｐｌｉｃａｂｉｌｉｔｙ　ａｎｄ　ａｃｃｕｒａｃｙ　ｃｏｍｐａｒｉｓｏｎ　ａｎａｌｙｓｉｓ　（上接第２２页）　究ｆＪ１．中国农村水利水电，２０１６，０５：５８—６０＋６４．　参考文献：　【６】宋楠．基于ＰＰＲ模型在辽宁中部区域灌溉用水预测【Ｊ】．水科学与　【１】刘艳梅，宋品玉。李强．阜蒙县节水农业存在问题及对策［Ｊ】．水土　工程技术，２０１６，０３：３９－－４１．　保持应用技术，２０ｌ２，ｏ４：４６—４７．　【７】徐世民，景淑娟．基于混沌蚁群算法的ＳＶＭ模型在水库枯水期径　【２】张兵，黄文生，周瑛，袁威．优化灌溉决策模型研究综述［Ｊ】．水土保　流预测中的应用们．水利技术监督，２０１６，０６：９９—１０２．　持应用技术，２０ｏ９。０５：２４—２６．　【８】雷智昌．水库优化调度混沌蚁群算法数学模型实例研究［Ｊ］．陕西　ｆ３】曲春阳，刘亿军．基于改进ＢＰ神经网络模型的区域农业灌溉用　水利，２０１０，Ｏｌ：９９—１０１．　水预测研究叨．吉林水利，２０１６，１０：１－３＋６．　［９】王志良，潘文学．混沌蚁群算法在水库优化调度中的应用叨．水利　［４】齐英．灌溉用水预测模型应用与分析叫．水利规划与设计，　与建筑工程学报，２００７。ｏ４：３１—３４．　２０１６，０３：１－２＋２１．　【ｌ０１曲晓宁．电力市场环境下梯级水库发电调度风险顸警系统研究　［５】王鑫东．加权马尔种夫链模型在农业灌溉用水预测中的应用研　『Ｄ】．郑州大学，２０１３．　Ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ａｇｒｉｃｕｌｔｕｒａｌ　ｉｒｒｉｇａｔｉｏｎ　ｗａｔｅｒ　ｑｕａｎｔｉｔｙ　ｉｎ　ｃｅｎｔｒａｌ　Ｌｉａｏｎｉｎｇ　ｐｒｏｖｉｎｃｅ　ｂａｓｅｄ　ｏｎ　ｃｈａｏｔｉｃ　ａｎｔ　ｃｏｌｏｎｙ　ａｌｇｏｒｉｔｈｍ　Ｓｏｎｇ　Ｎａｎ　Ａｂｓｔｒａｃｔ：Ｉｎ　ｔｈｉｓ　ｐａｐｅｒ，ｔｈｅ　ｃｈａｏｔｉｃ　ａｎｔ　ｃｏｌｏｎｙ　ａｌｇｏｒｉｔｈｍ　ｉｓ　ｕｓｅｄ　ｔｏ　ｏｐｔｉｍｉｚｅ　ｔｈｅ　ｃｏｍｐｕｔａｔｉｏｎａｌ　ｓｔｒｕｃｔｕｒｅ　ｏｆ　ＧＭ　（１，１）ｍｏｄｅｌ　ｔｏ　ｉｍｐｒｏｖｅ　ｔｈｅ　ｃｏｍｐｕｔａｔｉｏｎａｌ　ｅｆｉｆｃｉｅｎｃｙ　ａｎｄ　ｐｒｅｃｉｓｉｏｎ　ｏｆ　ＧＭ（１，１）ｍｏｄｅ１．Ｔｈｅ　ＧＭ（１，１）ｍｏｄｅｌ　ｂａｓｅｄ　ｏｎ　ｔｈｉｓ　ａｌｇｏｒｉｔｈｍ　ｉｓ　ｕｓｅｄ　ｉｎ　ｔｈｅ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ａｇｒｉｃｕｌｔｕｒｌａ　ｉｒｒｉｇａｔｉｏｎ　ｗａｔｅｒ　ｉｎ　ｃｅｎｔｒａｌ　Ｌｉａｏｎｉｎｇ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｔｈｅ　ｅｒｒｏｒ　ｏｆ　ａｎｎｕａｌ　ｉｒｒｉｇａｔｉｏｎ　ｗａｔｅｒ　ｉｓ　ｌｅｓｓ　ｔｈａｎ　２０％ａｎｄ　ｔｈｅ　ｃａｌｃｕｌａｔｉｏｎ　ｅｒｒｏｒ　ｉｓ　ｌｅｓｓ　ｔｈａｎ　２５％ｉｎ　ｄｉｆｅｒｅｎｔ　ｍｏｎｔｈｓ．ｈＴｅ　ＧＭ（１，１）ｍｏｄｅｌ　ｂａｓｅｄ　ｏｎ　ｃｈａｏｔｉｃ　ａｎｔ　ｃｏｌｏｎｙ　ｏｐｔｉｍｉｚａｔｉｏｎ　ａｌｇｏｒｉｔｈｍ　ｉｓ　ｂｅｔｔｅｒ　ｔｌｌａｎ　ｔｈｅ　ａ—　ｇｒｉｃｕｈｕｒｌａ　ｉｒｒｉｇａｔｉｏｎ　ｗａｔｅｒ　ｑｕａｎｔｉｔｙ　Ｐｒｅｄｉｃｔｉｏｎ　ａｃｃｕｒａｃｙ　ａｎｄ　ａｐｐｌｉｃａｂｉｌｉｔｙ．ｈＴｅ　ｒｅｓｅａｒｃｈ　ｒｅｓｕｌｔｓ　ｃａｎ　ｐｒｏｖｉｄｅ　ｒｅｆ－　ｅｒｅｎｃｅ　ｖａｌｕｅ　ｆｏｒ　ｗａｔｅｒ－ｓａｖｉｎｇ　ａｇｉｃｕｌｔｕｒｌａ　ｗａｔｅｒ　ｕｓｅ　ｐｌａｎｎｉｎｇ　ｉｎ　ｃｅｎｔｒｌａ　Ｌｉａｏｎｉｎｇ．　Ｋｅｙｗｏｒｄｓ：ｃｈａｏｔｉｃ　ａｎｔ　ｃｏｌｏｎｙ　ａｌｇｏｒｉｔｈｍ；ＧＭ（１，１）ｍｏｄｅｌ；ａｇｒｉｃｕｌｔｕｒｌａ　ｉｒｒｉｇａｔｉｏｎ　ｗａｔｅｒ　ｆｏｒｅｃａｓｔ；Ｌｉａｏｎｉｎｇ　ｃｅｎ－　ｔｒａｌ　ｒｅｇｉｏｎ　一２７—