不同蓄水方式对堆石坝应力变形的影响分析——以冶勒水电站工程为例

第４２卷第１１期　２　０　１　１年６月　文章编号：１００１—４１７９（２０１１）１１—００３５—０４　人　民　长　江　Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　Ｖｏ１．４２．Ｎｏ　１１　Ｊｕｎｅ，　２０１１　不同蓄水方式对堆石坝应力变形的影响分析　——以冶勒水电站工程为例　胡升伟　，张　权　，龚　军　，陈建康　（１．四川大学水利水电学院，四川成都６１００６５；　２．黄河勘测规划设计有限公司，河南郑州４５０００３；　３．中　国水电顾问集团成都勘测设计研究院，四川成都６１００７２）　摘要：水库先期蓄水和后期蓄水对大坝的应力变形有不同的影响。采用ｘＦ肯一张Ｅ—Ｂ模型对冶勒水电站堆　石坝进行三维有限元建模，结合不同坝体材料的计算参数，在先期蓄水和后期蓄水条件下，拟定两种计算方　案，对沥青混凝土堆石坝在竣工工况和正常蓄水位工况下的应力变形进行了三维有限元分析，并从应力变形　角度对其安全性进行了初步评价。通过计算分析和比较，这两种蓄水方式下大坝的应力变形均符合堆石坝的　一般分布规律，且量值在合理的范围内。　关　键　词：有限元；蓄水方式；应力变形；邓肯一张Ｅ—Ｂ模型　中图法分类号：ＴＶ６４１．４１　文献标志码：Ａ　冶勒水电站位于长江的支流楠亚河上游，是楠亚　下对大坝的应力变形进行深入的探究分析，并从应力　河的龙头水库，电站拦河坝为沥青混凝土心墙堆石坝，　坝顶高程２　６５４．５０　ｍ，最大坝高１２５．００　ｍ，水库正常　蓄水位２　６５０．００　ｍ，死水位２　６００．００　ｍ。由于冶勒水　电站的重要性和复杂性，对大坝的施工逐级加荷过程　变形的分析结果对大坝的安全性进行初步评价。　１．１　有限元模拟范围及网格划分　使用ＭＩＤＡＳ／ＧＴＳ岩土与隧道分析软件中的邓肯　一张Ｅ—Ｂ模型对冶勒沥青混凝土心墙堆石坝进行三　进行三维有限元模拟是很有必要的，特别是在两种蓄　水方式下，拟定出两种计算方案，根据所给计算参数，　对大坝应力变形进行计算分析和比较，从而了解大坝　的工作性态和内部结构的变化，并及时根据分析计算　成果提出合理化的意见和建议，为大坝的安全和稳定　提供强有力的理论支持。　维有限元建模。模型范围为：竖向５３７．６　ｍ，坝轴向　９２７　ｍ，顺河向９８９　ｍ。考虑大坝受力和变形的实际情　况，采用位移约束条件。将地基底面设置为全约束，左　右岸、顺河向方向施加法向约束，坝体部分均不设约　束。考虑到模型中仪器埋设位置应与实际相同，整个　模型自由剖分，采用四面体单元，网格节点总数为　２６　４３７个，单元总数为１３９　６２８个，其中地基和周围山　体４２　２５４个，坝体９３　７３４个。　邓肯一张模型属于非线性弹模型中的变模量模　型　，它是以广义虎克定律为基础，利用土体常规三　轴试验得到的应力一应变曲线建立模型参数关系的。　该模型假定弹模Ｅ、泊松比　、体变模量Ｋ和剪切模量　Ｇ是应力状态的函数，与应力路径无关。　１．２　计算参数　本文以冶勒沥青混凝土心墙堆石坝为工程实例，　通过对坝体、坝基的设计参数和反演参数进行计算和　１　不同蓄水方式的应力应变分析条件　采用邓肯一张Ｅ—Ｂ模型借助三维有限元建模，　根据不同的蓄水方式，拟定两种计算方案，在计算参数　收稿日期：２０１１—０２—１８　调整，得出的邓肯一张模型的计算参数见表１。混凝　土为线性材料，弹性模量Ｅ＝３０　ＧＰａ，泊松比取０．１６７，　密度为２４．０　ｋＮ／ｍ　，左岸基岩为石英闪长岩，为线性　通讯作者：胡升伟，男，硕士研究生，主要从事水工结构工程及基础工程工作。Ｅ—ｍａｉｌ：ｈｕｓｈｅｎｇｗｅ＠１６３．ｔｏｍ　３６　人　民　长　江　材料，弹性模量Ｅ＝１９　ＧＰａ，泊松比取０．２，密度为　第ｌ２～ｌ９级：坝体填筑到坝顶２　６５４．５０　ｎｌ高程；　２９．３　ｋＮ／ｍ　。　表１邓肯一张模型材料计算参数　（部位？）（　密ｋＮ・ｍ　）度　。（。）（。）ｋＰ：　ａ　，　Ｄ　注：　为土体的摩擦角；ｃ为黏聚力；　为破坏比；Ｆ、Ｇ、Ｄ为反映弹模与　泊松比的试验参数，Ｋ　Ｋ、　、ｍ为反映弹模与泊松比的试验参数。　１．３　加载程序　冶勒沥青混凝土心墙堆石坝的计算分析分两种方　案进行，两方案的计算参数相同。　１．３．１　先期蓄水　方案１在有限元计算时，认为冶勒沥青混凝土心　墙与坝体堆石同步填筑上升，荷载采用逐级施加的方　式。计算分２１级加荷，对施工逐级加荷过程进行模　拟。在加荷的过程中，模拟第１次蓄水，坝体施工到　２　６０５．００　ｒｎ高程时，分２级蓄水到２　６００．００　ｍ高程，然　后停止蓄水，继续填筑坝体直到坝顶。坝体填筑完成　后，分２级蓄水到正常蓄水位２　６５０．００／／１。载荷施加　的具体方法如下：　第１级：施加左岸基岩以及右岸、河床覆盖层自重　载荷；　第２—９级：坝体填筑到２　６０５．００　ｍ高程；　第１０～Ｉ１级：分２级施加水荷载到２　６００．００　ｍ高　程；　第２０—２ｌ级：分２级施加水荷载到正常蓄水位　２　６５０．Ｏ０　ｎｌ高程。　坝体填筑程序见图１。　１．３．２　后期蓄水　方案２在有限元计算时，先施工、后蓄水，不模拟　第１次蓄水，在坝体填筑到坝顶高程以后，分４级蓄水　到正常蓄水位。载荷施加的具体方法如下：　第１级：施加左岸基岩以及右岸、河床覆盖层自重　载荷；　第２～１７级：坝体填筑到坝顶２　６５４．５０　ｍ高程；　第１８～２１级：分４级施加水荷载到正常蓄水位　２　６５０．００　ｍ。　坝体填筑程序见图２。　２计算成果与分析　基于上述计算条件，对冶勒沥青混凝土心墙堆石　坝进行三维有限元静力计算，计算成果包括：①最大　横剖面（０＋２２０．００　Ｉｎ断面）水平位移、竖直沉降、主应　力分布；②沥青混凝土心墙坝轴向水平位移、竖直沉　降、主应力分布；③心墙上游面单元中主应力与上游　水压力之比等值线。工况为竣工工况和正常蓄水位工　况。符号规定为：压应力为正，拉应力为负，单位为　ＭＰａ。位移符号规定为：沉降向上为正，顺河向水平位　移向下游为正，坝轴向水平位移向左岸为正，单位为　ｃｍ。　２．１　坝体的应力与变形分析　方案１考虑了第１次蓄水，竣工工况下坝体顺河　向水平位移在水压力作用下不对称，总体偏向下游；两　种方案在正常蓄水位工况下，由于水压力较大，顺河向　位移均不对称，总体偏向下游且增幅明显，最大增幅为　方案１。各方案顺河向水平位移最大值发生在下游坝　壳约１／３坝高处，顺河向水平位移最大值结果见表２。　对于最大横剖面，两种方案各工况下的沉降均为　８０　ｃｍ左右，占最大坝高的比例均小于１％，沉降最大　值均发生在１／２坝高处，且正常蓄水位工况下的沉降　较竣工工况大，沉降最大值结果见表２。水平位置均　发生在心墙中轴线偏向上游，主要原因是堆石体Ｉ的　材料切线弹性模量参数较过渡层以及下游堆石体Ⅱ的　小，且堆石Ｉ区在坝体材料上游分区中占据大部分。　但是，方案１最大沉降发生的水平位置距心墙中轴线　的距离较方案２远，主要原因是方案１考虑了第１次　蓄水，坝体后期填筑过程中，上部和下部所受水压力不　均所致。　第１１期　胡升伟，等：不同蓄水方式对堆石坝应力变形的影响分析　３７　图１堆石坝坝体分级施工填筑（方案１）　上　图２堆石坝坝体分级施工填筑（方案２）　表２最大横剖面（０＋２２０．ｏｏ　ｍ断面）位移最大值计算结果　下游堆石区和心墙下游侧的小主应力均呈增加趋势。　２．２　沥青混凝土心墙应力与变形分析　对于坝轴向水平位移，坝体向右岸位移比向左岸　位移稍大，主要原因是坝体左岸为倾斜的基岩岸坡，河　谷下部以及右岸为深厚覆盖层，坝体填筑荷载必然使　坝体本身及覆盖层有向右位移的趋势。蓄水后，两种　方案坝轴向水平位移均有小幅增加，坝轴向水平位移　最大值见表３。　表３沥青混凝土心墙轴向水平位移最大值计算结果　坝体大小主应力极值均发生在沥青混凝土心墙底　部附近。对比分析主应力等值线可知，对于最大横剖　面，方案１竣工工况的小主应力在心墙２　６００．Ｏ０　ｍ高　程以下出现应力集中，主要原因是方案１考虑第１次　蓄水，初次蓄水高程为２　６００．Ｏ０　ｍ，在此高程以下心墙　受到水压力作用，而小主应力方向近乎顺河流方向。　坝体大主应力等值线显示，沥青混凝土心墙的应力并　不比同高程处的过渡层应力有明显降低，拱效应不明　显。正常蓄水位工况下，由于水压力的作用，两种方案　从竖直沉降等值线分析可知，对于沥青混凝土心　计算的上游堆石区的大主应力均低于下游堆石区的大　墙，两种方案各工况下的沉降最大值均发生在约１／２　主应力。同时由于心墙在水压力作用下向下游变形，　坝高处，量大值在７８　ｃｍ左右，占最大坝高的比例均小　上游堆石区和心墙上游侧的小主应力均呈减小趋势，　于１％，心墙沉降最大值结果见表３。从表中可以看　３８　人　民　长　江　出，不考虑第１次蓄水，心墙沉降较大，主要原因是不　考虑第１次蓄水，坝体最大沉降发生的水平位置距心　墙中轴线较近。对比各方案竣工工况和正常蓄水工况　３结论　本文采用邓肯一张Ｅ—Ｂ模型，以是否考虑第１　下的沉降，可以看出，蓄水后心墙沉降有所减小，但减　次蓄水以及采用不同坝体材料计算参数为标准，拟定　在先期蓄水和后期蓄水方式下的两种计算方案，通过　计算成果的比较分析，可以得出以下结论。　小幅度不大，主要原因是蓄水后坝体最大沉降偏离心　墙中轴线向上游。　沥青混凝土心墙大小主应力最大值均发生在心墙　（１）不同蓄水方式条件下，坝体的应力和变形也　不同。顺河向水平位移总体偏向下游，且先期蓄水条　底部，且量值沿高程从坝顶向下逐渐增大。对比分析　各方案主应力等值线，对于沥青混凝土心墙，各工况　件下增幅较大；两种方案各工况下的沉降均为８０　ｃｍ　下，两种方案主应力相差不大，最大小主应力为１．９　ＭＰａ左右，最小小主应力为一０．０９　ＭＰａ左右，最大大　主应力为３．９　ＭＰａ左右，最小大主应力为０．Ｏ１　ＭＰａ左　右，蓄水后，心墙内大小主应力极值均有小幅增加，小　主应力拉力区有所减小。　对于冶勒沥青混凝土心墙堆石坝，由于左右岸严　重不对称，坝体及心墙有整体向右位移的趋势，从坝轴　向水平位移等值线分析出左岸向右的位移稍大，这样　容易导致左岸坝肩处出现拉应力，甚至拉裂缝。从两　种方案各工况下的小主应力等值线可以看出，左岸坝　肩出现拉应力，竣工工况下最大拉应力为一０．１１１　ＭＰａ，正常蓄水位工况下最大拉应力为一０．０５２　ＭＰａ，　拉应力区范围不大。如果按照ｃ／ｔａｎ￣ｐ：０．２／ｔａｎ２０”＝　０．３９　ＭＰａ作为沥青混凝土的抗拉强度，心墙不会开　裂。　２．３　沥青混凝土心墙抗水力劈裂能力分析　对冶勒沥青混凝土心墙堆石坝，本文采用心墙上　游面水压力大于同高程心墙上游面单元心墙中主应力　作为判别发生水力劈裂的标准。图３分别给出了两种　方案正常蓄水位工况下沥青混凝土心墙上游面单元中　主应力与上游水压力之比等值线。可以看出，两种方　案的该比值均大于１．０，说明这两种方案的心墙均能　承受高水压的作用，不会发生水力劈裂现象。　（ｂ）方案２　图３心墙上游水压力与心墙中主应力之比等值线　左右，占最大坝高的比例均小于１％，且正常蓄水位工　况下的沉降较竣工工况大；坝体的大小主应力最大值　均发生在沥青混凝土心墙底部附近，其中在先期蓄水　的竣工工况下，小主应力在高程２　６００　ｍ以下会出现　应力集中。　（２）不同蓄水方式对沥青混凝土心墙应力变形的　影响显著。蓄水后，两种方案坝轴向水平位移均有小　幅增加；通过对竖直等值线分析，两种方案在各工况下　的沉降最大值均发生在约１／２坝高处，量值在７８　ｃｍ　左右，占最大坝高的比例均小于１％，后期蓄水时，心　墙沉降较大；沥青混凝土心墙大小主应力极值均发生　在心墙底部，且量值沿高程从坝顶向下逐渐增大，通过　计算比较分析，心墙不会开裂。从计算结果看，坝体与　心墙存在沉降差，但其差值在合理的范围内，不会发生　断裂现象。　（３）不同蓄水方式下，采用心墙上游面水压力大　于同高程心墙上游面单元主应力作为发生水力劈裂的　判别标准，计算结果表明，该沥青混凝土心墙堆石坝不　会发生水力劈裂现象。　参考文献：　［１］Ｄｕｎｃａｎ　Ｊ　Ｍ，Ｃｈａｎｇ　ｃ　Ｙ．Ｎｏｎｌｉｎｅａｒ　ａｎａｌｙｓｉｓ　ｏｆ　ｓｔｒｅｓｓ　ａｎｄ　ｓｔｒａｉｎ　ｉｎ　ｓｏｉｌｓ，Ｊ．Ｓ．Ｍ．Ｆ．Ｄ．ＡＳＣＥ，１９７０，９６（ＳＭ５）．　［２］　Ｋｕｌｈａｗｙ　Ｆ　Ｈ，Ｄｕｎｃａｎ　Ｊ　Ｍ．Ｆｉｎｉｔｅ　Ｅｌｅｍｅｎｔ　Ａｎａｌｙｓｉｓ　ｏｆ　Ｓｔｒｅｓｓ　ａｎｄ　Ｍｏｖｅｍｅｎｔ　ｉｎ　Ｅｍｂａｎｋｍｅｎｔ　ｄｕｒｉｎｇ　Ｃｏｎｓｔｒｕｃｔｉｏｎ，Ｊ．Ｓ．Ｍ．Ｆ．Ｄ．，　ＡＳＣＥ，１９７２，９８（ＳＭ７）．　［３］　Ｋｕｌｈａｗｙ　Ｆ　Ｈ，Ｄｕｎｃａｎ　Ｊ　Ｍ．Ｓｔｒｅｓｓｅｓ　ａｎｄ　Ｍｏｖｅｍｅｎｔｓ　ｉｎ　Ｏｒｏｖｉｌｌ　Ｄａｍ，　ＡＳＣＥ，１９７２，９８（ＳＭ７）．　［４］　钱家欢，殷宗泽．土工原理与计算（第二版）［Ｍ］．北京：中国水利　水电出版社，１９９６．　（编辑：徐诗银）　（下转第８２页ｊ　８２　人　民　长　江　２０１１生　［Ｊ］．中华昆虫．１９９７，１７（３）：１５２—１６２．　　［３］　姜建国，沈韫芬．用于评价水污染的生物指数［Ｊ］．云南环境科学，２０００。１９（Ｓ）：２５１—２５３．　［９］Ｎｏｖａｋ　Ｍ，Ｂｏｄｅ　Ｒ．Ｐｅｒｃｅｎｔ　ｍｏｄｅｌ　ａｆｆｉｎｉｔｙ：ａ　ｎｅｗ　ｍｅａｓｕｒｅ　ｏｆｍａｅｒｏｉｎ。　ｖｅｒｔｅｂｒａｔｅ　ｃｏｍｍｕｎｉｔｙ　ｃｏｍｐｏｓｉｔｉｏｎ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　ｔｈｅ　Ｎｏｒｔｈ　Ａｍｅｒｉｃａｎ　Ｂｅｎｔｈｏｌｏｇｉｅａｌ　Ｓｏｃｉｅｔｙ．１９９２，１１（１）：８０—８５．　［４］　尤平，任辉．底栖动物及其在水质评价和监测上的应用［Ｊ］．淮北　煤师院学报：自然科学版，２００１，２２（４４）：４４—４８．　［５］　周长发，归鸿，周开亚．中国蜉蝣目稚虫科检索表（昆虫纲）［Ｊ］．　［１Ｏ１　ＤＲ　Ｂ．Ｔｈｅ　ｕｓｅ　ｏｆ　ｐｅｒｃｅｎｔ　ｍｏｄｅｌ　ａｆｉｆｎｉｔｙ　ｔｏ　ａｓｓｅｓｓ　ｔｈｅ　ｅｆｆｅｃｔｓ　ｏｆ　ａｇｒｉ・　ｃｕｌｔｕｒｅ　ｏｎ　ｂｅｎｔｈｉｃ　ｉｎｖｅｒｔｅｂｒａｔｅ　ｃｏｍｍｕｎｉｔｉｅｓ　ｉｎ　ｈｅａｄｗａｔｅｒ　ｓｔｒｅａｍｓ　ｏｆ　ｓｏｕｔｈｅｒｎ　Ｏｎｔａｒｉｏ［Ｊ］，Ｃａｎａｄａ．Ｆｒｅｓｈｗａｔｅｒ　Ｂｉｏｌｏｇｙ．１９９６，３６（２）：　３９７—４１０．　南京师范大学学报：自然科学版，２００３，２６（２）：６５—６８．　【６］　王洪铸．中国小蚓类研究一附中国南极长城站附近地区两新种　（Ｉ］．北京：高等教育出版社，２００２．　Ｔｉａｎ　Ｌ．Ａｑｕａｔｉｃ　ｉｎｓｅｃｔｓ　ｏｆ　Ｃｈｉｎａ　ｕｓｅｆｕｌ　ｆｏｒ　ｍｏｎｉｔｏ。　［７］　Ｍｏｍｃ　Ｊ，Ｙａｎｇ　Ｌ，［１　１］　Ｂａｒｔｏｎ　Ｄ，Ｆａｒｍｅｒ　Ｍ．Ｔｈｅ　ｅｆｆｅｃｔｓ　ｏｆ　ｃｏｎｓｅｒｖａｔｉｏｎ　ｔｉｌｌａｇｅ　ｐｒａｃｔｉｃｅｓ　ｏｎ　ｂｅｎｔｈｉｃ　ｉｎｖｅｒｔｅｂｒａｔｅ　ｃｏｍｍｕｎｉｔｉｅｓ　ｉｎ　ｈｅａｄｗａｔｅｒ　ｓｔｒｅａｍｓ　ｉｎ　ｓｏｎｔｈｗｅａ１．　ｅｒｎ　Ｏｎｔａｒｉｏ［Ｊ］，Ｃａｎａｄａ．Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｐｏｌｌｕｔｉｏｎ．１９９７，９６（２）：　２０７—２１５．　ｒｉｎｇ　ｗａｔｅｒ　ｑｕａｌｉｔｙ［Ｍ］．Ｎａｎｊｉｎｇ：Ｈｏｈａｉ　Ｕｎｉｖｅｒｓｉｔｙ　Ｐｒｅｓｓ，１９９４．　　［８］　刘月荚．中国经济动物志：淡水软体动物［Ｍ］．北京：科学出版社，１９７９．　（编辑：常汉生）　Ａｓｓｅｓｓｍｅｎｔ　ｏｎ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｏｆ　Ｗｕｊｉａｎｇ　Ｒｉｖｅｒ　ｂｙ　ｐｅｒｃｅｎｔ　ｍｏｄｅｌ　ａｆｉｎｉｆｔｙ　ＣＨＩ　Ｓｈｉｙｕｎ，ＨＵ　Ｊｕｘｉａｎｇ，ＺＨＡＯ　Ｘｉａｎｆｕ，ＭＡ　Ｐｅｉｍｉｎｇ，ＬＩ　Ｍｉｎｇ　（１ｎｓｉｔｉｔｕｔｅ　ｏｆ　Ｈｙｄｒｏｅｃｏｌｏｇｙ，Ｍｉｎｉｓｔｒｙ　ｏｆ　Ｗａｔｅｒ　Ｒｅｓｏｕｒｃｅｓ　ａｎｄ　Ｃｈｉｎｅｓｅ　Ａｃａｄｅｍｙ　ｏｆ　Ｓｃｉｅｎｃｅｓ，Ｗｕｈａｎ　４３００７９，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：　Ｉｎ　ｏｒｄｅｒ　ｔｏ　ａＳＳｅＳＳ　ｔｈｅ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｂｙ　ａｑｕａｔｉｃ　ｏｒｇａｎｉｓｍｓ，１４　ｓａｍｐｌｉｎｇ　ｓｅｃｔｉｏｎｓ　ａｒｅ　ａｒｒａｎｇｅｄ　ｉｎ　ｔｈｅ　ｍａｉｎ　ｓｔｅｍ　ａｎｄ　ｂｒａｎｃｈｅｓ　ｏｆ　Ｗｕｊｉａｎｇ　Ｒｉｖｅｒ，ａｎｄ　ｂａｓｅｄ　ｏｎ　ｔｈｅ　ｓａｍｐｌｅｄ　ｄａｔａ，ｔｈｅ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｉｎ　ｔｈｅ　ｄｏｗｎｓｔｒｅａｍ　ｏｆ　Ｗｕｊｉａｎｇ　Ｒｉｖｅｒ　ｄｕｒｉｎｇ　ｔｈｅ　ｈｙ－　ｄｒｏｐｏｗｅｒ　ｃｏｎｓｔｒｕｃｔｉｏｎ　ｐｅｒｉｏｄ　ｉｓ　ａｓｓｅｓｓｅｄ　ｗｉｔｈ　Ｐｅｒｃｅｎｔ　Ｍｏｄｅｌ　Ａｆｆｉｎｉｔｙ（ＰＭＡ）ｅｓｔａｂｌｉｓｈｅｄ　ｂｙ　Ｎｏｖａｋ，ｅｔｃ．Ｔｈｅ　ａｓｓｅｓｓｅｄ　ｒｅｓｕｌｔｓ　ａｒｅ　ｃｏｍｐａｒｅｄ　ｗｉｔｈ　ｔｈｏｓｅ　ｂｙ　ｓｐｅｃｉｅｓ　ａｂｕｎｄａｎｃｅ，ＥＰＴ％ａｎｄ　ｃｏｒｒｅｌａｔｉｖｅ　ｉｎｄｅｘ　ｏｆ　ｓｐｅｃｉｅｓ　ｄｉｖｅｒｓｉｔｙ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｆｏｒ　ｔｈｅ　ｂｒａｎ—　ｅｈｅｓ，ｅｘｃｅｐｔ　ｔｈｅ　Ｙｕｊｉａｎｇ　Ｒｉｖｅｒ　ｉｓ　ｍｏｄｅｒａｔｅ，ｔｈｅ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｏｆ　ｏｔｈｅｒ　ｓｅｃｔｉｏｎｓ　ｉｓ　ｉｎ　ａ　ｇｏｏｄ　ｓｔａｔｅ；ｔｈｅ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｏｆ　ｔｈｅ　ｍａｉｎ－　ｓ￣ｅａｍ　ｉｓ　ｍｏｄｅｒａｔｅ　ｅｘｃｅｐｔ　ｔｈａｔ　ｔｈｅ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｏｆ　Ｐｅｎｇｓｈｕｉ　ａｎｄ　Ｈｏｎｇｄｕ　ｓｔａｔｉｏｎｓ　ｉｓ　ｇｏｏｄ．Ｏｖｅｒａｌｌ，ｔｈｅ　ｗａｔｅｒ　ｑｕａｌｉｔｙ　ｏｆ　ｔｒｉｂｕｔａｒｉｅｓ　ｉｓ　ｂｅｔｔｅｒ　ｔｈａｎ　ｔｈａｔ　ｏｆ　ｔｈｅ　ｍａｉｎｓｔｒｅａｍ．Ｉｔ　ｉｓ　ｒｅｖｅａｌｅｄ　ｔｈａｔ　ｏｒｇａｎｉｃ　ｐｏｌｌｕｔｉｏｎ　ｉｓ　ｎｏｔ　ｔｈｅ　ｍａｉｎ　ｐｏｌｌｕｔｉｏｎ　ｓｏｕｒｃｅ，ｗｈｏｓｅ　ｒｅａｓｏｎ　ｍａｙ　ｂｅ　ｔｈｅ　ｈｙｄｒｏｐｏｗｅｒ　ｃｏｎｓｔｒｕｃｔｉｏｎ　ａｎｄ　ｆｒｅｑｕｅｎｔ　ｓｈｉｐｐｉｎｇ　ａｃｔｉｖｉｔｉｅｓ．Ｔｈｅ　ａｎａｌｙｓｉｓ　ｓｈｏｗｓ　ｔｈａｔ　ｔｈｅ　ｉｎｄｅｘ　ｏｆ　ＰＭＡ　ｈａｓ　ｓｉｇｎｉｉｃａｎｔｆｌｙ　ｃｏｒｒｅｌａ－　ｌｉｏｎ　ｗｉｔｈ　ｓｐｅｃｉｅｓ　ａｂｕｎｄａｎｃｅ，ＥＰＴ％ａｎｄ　ｃｏｒｒｅｌａｔｉｖｅ　ｉｎｄｅｘ　ｏｆ　ｓｐｅｃｉｅｓ　ｄｉｖｅｒｓｉｔｙ．　Ｋｅｙ　ｗｏｒｄｓ：Ｐｅｒｃｅｎｔ　Ｍｏｄｅｌ　Ａｆｆｉｎｉｔｙ；ｚｏｏｂｅｎｔｈｏｓ；ｓｐｅｃｉｅｓ　ｄｉｖｅｒｓｉｔｙ　ｉｎｄｅｘ；ｗａｔｅｒ　ｑｕａｌｉｔｙ；Ｗｕｊｉａｎｇ　Ｒｉｖｅｒ　，…，…，一　，…，　ｏ●】Ｈ＾ｌ　（上接第３８页）　Ａｎａｌｙｓｉｓ　ｏｎ　ｉｎｆｌｕｅｎｃｅ　ｏｆ　ｓｔｒｅｓｓ—－ｄｅｆｏｒｍａｔｉｏｎ　ｏｆ　ａｓｐｈａｌｔ　ｃｏｎｃｒｅｔｅ　ｃｏｒｅ　ｒｏｃｋｆｉｌｌ　ｄａｍ　ｗｉｔｈ　ｄｉｆｅｒｅｎｔ　ｉｍｐｏｕｎｄｍｅｎｔ　ｗａｙｓ：ｃａｓｅ　ｏｆ　Ｙｅｌｅ　Ｈｙｄｒｏｐｏｗｅｒ　Ｓｔａｔｉｏｎ　ＨＵ　Ｓｈｅｎｇｗｅｉ　，ＺＨＡＮＧ　Ｑｕａｎ　，ＧＯＮＧ　Ｊｕｎ　，ＣＨＥＮ　Ｊｉａｎｋａｎｇ　（１．Ｃｏｌｌｅｇｅ　ｆ　ｏＷａｔｅｒ　Ｒｅｓｏｕｒｃｅ　ａｎｄ　Ｈｙｄｒｏｐｏｗｅｒ，Ｓｉｃｈｕａｎ　Ｕｎｉｖｅｒｓｉｔｙ，Ｃｈｅｎｇｄｕ　６１００６５，Ｃｈｉｎａ；２．Ｙｅｌｌｏｗ　Ｒｉｖｅｒ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｃｏｎ—　ｓｕｌｔｉｎｇ　Ｃｏ．，Ｌｔｄ．，Ｚｈｅｎｇｚｈｏｕ　４５０００３，Ｃｈｉｎａ；３．Ｈｙｄｒｏｃｈｉｎａ　Ｃｈｅｎｇｄｕ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｃｏｒｐｏｒａｔｉｏｎ，Ｃｈｅｎｇｄｕ　６１００７２，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：　Ｅａｒｌｙ　ａｎｄ　ｌａｔｅ　ｉｍｐｏｕｎｄｍｅｎｔ　ｗａｙｓ　ｈａｖｅ　ｄｉｆｆｅｒｅｎｔ　ｅｆｆｅｃｔｓ　ｏｎ　ｓｔｒｅｓｓ—ｄｅｆｏｒｍａｔｉｏｎ　ｏｆ　ｒｅｓｅｒｖｏｉｒ　ｄａｍ．Ｂｙ　ａｄｏｐｔｉｎｇ　Ｄｕｎ－　ｃａｎ—Ｃｈａｎｇ　Ｅ—Ｂ　ｍｏｄｅｌ，ｗｅ　ｅｓｔａｂｌｉｓｈｅｄ　３　Ｄ　ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ　ｍｏｄｅｌ　ｆｏｒ　ａｓｐｈａｌｔ　ｃｏｎｃｒｅｔｅ　ｃｏｒｅ　ｒｏｅｋｆｉｌｌ　ｄａｍ　ｉｎ　Ｙｅｌｅ　Ｈｙｄｒｏｐｏｗｅｒ　Ｓｔａ・　ｔｉｏｎ．Ｃｏｍｂｉｎｉｎｇ　ｃａｌｃｕｌａｔｉｏｎ　ｐａｒａｍｅｔｅｒｓ　ｏｆ　ｄｉｆｆｅｒｅｎｔ　ｍａｔｅｒｉａｌｓ　ｏｆ　ｄａｍ，ｗｅ　ｗｏｒｋｅｄ　ｏｕｔ　ｔｗｏ　ｎｕｍｅｒｉｃａｌ　ｐｒｏｃｅｄｕｒｅｓ　ｕｎｄｅｒ　ｔｈｅ　ｅｏｎｄｉ—　ｔｉｏｎｓ　ｏｆ　ｅａｒｌｙ　ａｎｄ　ｌａｔｅ　ｗａｔｅｒ　ｓｔｏｒａｇｅ　ｗａｙｓ．Ｆｕｒｔｈｅｒｍｏｒｅ，ｗｅ　ｃａｒｒｉｅｄ　ｏｕｔ　３　Ｄ　ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ　ａｎａｌｙｓｉｓ　ｏｎ　ｓｔｒｅｓｓ—ｄｅｆｏｒｍａｔｉｏｎ　ｏｆ　ｔｈｅ　ａｓｐｈａｌｔ　ｃｏｎｃｒｅｔｅ　ｃｏｒｅ　ｒｏｃｋｆｉｌｌ　ｄａｍ　ｉｎ　ｃｏｍｐｌｅｔｉｏｎ　ｕｎｄｅｒ　ｔｈｅ　ｃｏｎｄｉｔｉｏｎｓ　ｏｆ　ｎｏｒｍａｌ　ｗａｔｅｒ　ｓｔｏｒａｇｅ　ｌｅｖｅｌ，ａｎｄ　ｐｒｅｌｉｍｉｎａｒｉｌｙ　ｅｖａｌｕａｔｅｄ　ｉｔｓ　ｓａｆｅｔｙ　ｆｒｏｍ　ｔｈｅ　ｐｅｒｓｐｅｃｔｉｖｅ　ｏｆ　ｓｔｒｅｓｓ—ｄｅｆｏｒｍａｔｉｏｎ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｔｈｒｏｕｇｈ　ｔｈｅ　ａｎａｌｙｓｉｓ　ａｎｄ　ｃｏｍｐａｒｉｓｏｎ，ｓｔｒｅｓｓ—ｄｅｆｏｒｍａ－　ｔｉｏｎ　ｏｆ　ｄａｍ　ｕｎｄｅｒ　ｔｈｅ　ｔｗｏ　ｗａｔｅｒ　ｉｍｐｏｕｎｄｍｅｎｔ　ｗａｙｓ　ａｒｅ　ｉｎ　ａｃｃｏｒｄａｎｃｅ　ｗｉｔｈ　ｔｈｅ　ｇｅｎｅｒａｌ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｒｏｃｋｆｉｌｌ　ｄａｍ，ａｎｄ　ｔｈｅ　ｖａｌｕｅｓ　ａｒｅ　ｗｉｔｈｉｎ　ａ　ｒｅａｓｏｎａｂｌｅ　ｓｃｏｐｅ．　Ｋｅｙ　ｗｏｒｄｓ：　ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ；ｉｍｐｏｕｎｄｍｅｎｔ　ｗａｙｓ；ｓｔｒｅｓｓ—ｄｅｆｏｒｍａｔｉｏｎ；Ｄｕｎｃａｎ—Ｃｈａｎｇ　Ｅ—Ｂ　ｍｏｄｅｌ