基于BP神经网络的土坝渗流场反演研究

2010年西安理工大学研究生学术年会-1-基于BP神经网络的土坝渗流场反演研究任喜平，钟声（西安理工大学，陕西西安，710048）摘要：土坝坝体、坝基土壤或岩体、两岸土体或岩体的渗透系数则是分析土坝渗流场的一个重要参数。它关系到人们是否能够真正了解大坝的实际安全状况。本文将现场监测和反演分析结合起来反演渗透系数，对于材料比较复杂的地基尤为如此，地基渗透系数牵扯到地基沉降、渗漏、滑体等众多安全因素。从该角度出发，针对监测资料中的测压管水头与渗透系数之间的非线性关系，引入BP神经网络法、遗传算法，用于反演土石坝渗透系数，对比参照实际设计渗透系数。关键词：BP神经网络法；遗传算法；渗透系数；非线性中图分类号：文献标识码：文章编号：InfiltratestheflowfieldinversionresearchbasedontheBPneuralnetwork'searthstonedamRenXi-ping，Zhongsheng（Xi’anUniversityofTechnology,Xi’an,Shaanxi710048,China）Abstract:Theearthdam,thebaseofdamsoileithertherockmass,bothbankssoilbodyortherockmasspenetrationcoefficientisanalyzestheearthdamtoinfiltratetheflowfieldakeyparameter.Itrelatesthepeoplewhethercanunderstandthedamtrulytheactualsafecondition.Thisarticlefromthereforemonitorandtheinversionanalysisonthespotunifiestheinversionpenetrationcoefficient.Thisarticlefromthisangleembarking,inviewofbetweenmonitoringdata'sinmeasuringtubefloodpeakandpenetrationcoefficient'smisalignmentrelations,introducestheBPnervenetworkmethod,thegeneticalgorithm,usesintheinversionearthstonedampenetrationcoefficient,thecontrastreferenceactualdesignpenetrationcoefficient.Keywords:BPneuralnetwork;Geneticalgorithm;Penetrationcoefficient;Misalignment1引言人工神经网络的发展，经历了一个曲折艰难的过程。早在1943年，心理学家WarrenMcCulloch和数理逻辑学家WalterPitts从信息处理的角度出发，采用数理逻辑的方法，对生物神经细胞的动作进行研究，提出了形式神经元的模型，称为MP模型。通过MP模型提出了神经元的形式化数学描述和网络结构方法，证明了单个神经元能执行逻辑功能，从而开创了人工神经网络研究的新纪元。1949年心理学家DonaldHebb明确提出了突触联系强度可变假设，认为突触联系强度是随前后神经元的活动而变化，从而说明了神经元连接强度是可以修正的，为神经网络的产生奠定了基础。随后，感知器（Perception）、自适应线形神经元（AdaptiveLinearNeuron）等网络与运行规则相继被提出，早期主要应用子模式识别、天气预报和自适应控制。然而，简单的感知器只能局限于处理线形问题，加之当时VonNeumann数字计算机和人工智能正处于发展的黄金时期，暂时掩盖了发展新型计算机和人工智能新途径的必要性和迫切性，使人工收稿日期：作者简介：任喜平（1984-），男（汉族），会宁人，西安理工大学研究生，硕士，主要研究方向为水工结构、水工建筑物应力仿真分析及监测控制。Email：[e-mail].2010年西安理工大学研究生学术年会-2-神经网络处于低潮期。直到1982年，美国加州工学院物理学家JohnHopfied提出了Hopfied神经网络模型，引入了“能量模型”概念，给出了网络稳定的依据，并成功解决了著名的“旅行商问题”，成为了人工神经网络走向成熟的里程碑。1974年P.Werbos在其博士论文中提出了第一个适合多层网络的学习算法，但该算法并未受到足够的重视和广泛的应用，直到20世纪80年代中期，美国加利福尼亚的Rumelhart、Mcclelland和他的同事们成立了一个PDP（paralleldistributedprocession）小组，研究并行分布式信息处理方法，探索人类认知的微结构，并于1986年发表了ParallelDistributedProcessing一书，将该算法应用于神经网络的研究，才使之成为迄今为止最著名的多层网络学习算法—BP算法，由此算法训练的神经网络，称之为BP神经网络。在人工神经网络的实际运用中，BP神经网络广泛应用于函数逼近、模糊识别/分类、数据压缩等，80%~90%的人工神经网络模型是采用BP网络或它的变化形式，它也是前馈网络的核心部分，体现了人工神经网络最精华的部分。BP算法采用梯度下降方法在权矢量空间总求取误差函数的极小值，能够使误差函数极小化的权值组合即为学习问题的解答。2BP神经网络的基本思想BP网络的基本原理是，学习过程是由信号的正向传播与误差的反向传播两个过程组成的。正向传播时，输入样本从输入层输入，经隐层逐层处理后，传出输出层。若输出层的实际输出与期望的输出不符，则转入误差的反向传播阶段。误差反传是将输出误差以某种形式通过隐层向输入层反传，并将误差分摊给各个层的所有单元，从而获得各层单元的误差信号，此信号作为修正各单元权值的依据。这种信号正向传播与误差反向传播的各层权值调整过程，也就是网络不断学习的过程。如此反复进行，直至误差满足设定要求。BP神经元与其他神经元类似，不同的是BP神经元的传输函数为非线性函数，最常用的是logsig和tansig函数，有的输出层也采用线性函数（purelin）其输出为log()asigWpb(2.1)BP网络一般为多层神经网络。由BP神经元构成的图2-1，BP网络的信息从输入层流向输出层，因此是一种多层前馈神经网络。图2-1BP神经网络的一般结构Figure2-1BPneuralnetworkgeneralstructure输入层X隐含层A输出层Y2010年西安理工大学研究生学术年会-3-如果多层BP网络的输出层采用S形传输函数（如logsig），其输出值将会限制在一个较小的范围内（0,1）；而采用线形传输函数则可以取任意值。由于误差反传网络增加了隐层，并有相应的学习规律可循，使其具有对非线性模式的识别能力，特别是数学意义明确，步骤分明的学习算法。尽管误差反向算法存在收敛速度慢、局部极值等缺点。但可以通过各种改进措施来提高它的收敛速度、克服局部极值现象，而且具有简单、易行、计算量小，并行性强等特点。目前是多层前向网络的首选算法。3工程实例黄家畔水库工程位于无定河一级支流大理河上游，枢纽工程位于靖边县小河乡黄家畔村，水库枢纽至靖边县有45km便道相通，S204省道从靖边县城通过，交通便利。水库枢纽由大坝、输水建筑物组成。大坝为水坠均质土坝，最大坝高32.8m，坝顶长120.0m，坝顶宽6.0m，总库容167万m3，兴利库容47万m3，灌溉面积3000亩。是一座以防汛、拦泥为主，兼有生态补水、灌溉功能的小（1）型水库。枢纽工程属Ⅳ等，主要建筑物等级4级。该水库1973年9月动工修建，1977年10月建成投入运行，水库运行30多年来，在防汛、拦沙等方面发挥了重要的作用，取得了较大的社会效益和较好的经济效益。黄家畔水库枢纽位于靖边县小河乡黄家畔村大理河的上游，大理河自北西流向南东，河谷呈“U”型发育，河床高程950m，河谷宽60m～80m，河谷覆盖第四系全新统含砾粉质壤土，厚度1.0～2.0m。下伏第三系（N）棕红色泥质细砂岩、中厚层状、泥钙质胶结，易风化。3.1黄家畔监测资料3-1黄家畔横剖面图Figure-1Huangjiapanaxlesectionaldrawing图3-1为黄家畔横剖面图，由于一般均质土石坝坝体内土质比较单一所以本算例反演相对材料比较复杂的地基渗透系数，以下是监测资料给出的不同水位下的10组测压管水头数据，每组数据取5个相同位置下的测压管水头，具体数值下表列出：表3-1测压管水头数据（单位:m）Table3-1measuringtubefloodpeakdata(unit:m)测压管水头12345T3wT3w42alQ42alQ4Qs2010年西安理工大学研究生学术年会-4-118.4315.7611.458.876.98216.4414.9410.698.186.38315.3814.219.767.335.64414.4513.289.026.735.22513.9312.828.606.364.98613.0311.937.755.524.21712.6211.527.395.173.90812.2211.026.944.713.52911.5710.366.344.123.021010.859.705.803.582.56以及相应的地基渗透系数期望值，即实际地基渗透系数参数值。表3-2实际渗透系数参数值(单位:m/s)Table4-2actualpenetrationcoefficientparametervalue(unit:m/s)1234572107510611061.510621067891063.21064.51064.781065.54106610根据以上数据，本文将分别利用BP神经网络与基于遗传算法的BP神经网络法反演黄家畔渗透系数，即预测值，并将其与目标期望值进行比较。从而应证两种算法的可行性。最后两者互相比较，选出其中较为优越的方法，并对两种方法做出评论。3.2黄家畔BP神经网络成果分析运用第2章BP网络算法为基础，将上述两种数据代入自编的BP神经网络程序，并将1-6组以及第9组数据作为网络训练值，将第7，8，10组数据做为网络测试值，将他们归一化后进行运算，得到相应的归一化后的渗透系数预测值，而后将其反归一化后，与原目标期望值进行比较。本程序利用traingdx方法进行训练，迭代次数为5000次，目标为0.01，程序在第118次收敛，说明了该训练方法可行性，以下是将上述两种数据代入该程序后的收敛图形2010年西安理工大学研究生学术年会-5-图3-2BP网络收敛图形Figure3-2BPnetworkrestraininggraph经过网络训练，得到相应的1-6组以及第9组数据的训练预测输出为下表：表3-3渗透系数训练输出值(单位:m/s)Table3-3penetrationcoefficienttrainingvalueofexports(unit:m/s)123456972.391075.401061.041061.571062.081063.281065.54310经过网络测试运算，得到相应的第7，8，10组数据的测试输出为下表：表3-4渗透系数测试输出值(单位:m/s)Table3-4penetrationcoefficienttestvalueofexports(unit:m/s)781063.951064.871065.7810将经由程序得出的这10组预测输出值与原目标期望值进行比较，发现最大的误差为19.5%。下表为BP神经网络反演得出的结果与原参数比较误差值：表