瀑布沟高心墙土石坝渗流分析

第５期　２０１３年１０月　水利水运工程同　学报　Ｎｏ．５　０ｅｔ．２０１３　ＨＹＤＲ０－Ｓａ匝ＮＣＥ　ＡＮＤＮ团匝ｔＩＮＧ　瀑布沟高心墙土石坝渗流分析　涂扬举　，王文涛　，薛新华。　（１．国电新疆电力有限公司，新疆乌鲁木齐８３００６３；２．国电大渡河流域水电开发有限公司，四川成都　６１００４１；３．四川大学水利水电学院，四ｊｌｌ成都６１００６５）　摘要：在深覆盖层地基上修建高土石坝，其防渗体系的可靠性是一项关键技术问题．防渗墙与土质防渗体连接　处是抵御渗透破坏的关键部位．根据瀑布沟土石坝防渗体系的结构特点，利用有限元方法对瀑布沟土石坝进行　了渗流分析．结果表明：坝体渗流与应力变形计算时，副防渗墙按４０％承担水头较为合适；连接部位的渗透坡降　是非均匀变化的，混凝土结构顶部的渗透坡降较大，心墙底部出口处的渗透坡降较小；坝体与两岸相接部位心　墙底部渗流出口处的坡降最大．研究结论可以为类似工程提供参考和借鉴．　关键词：防渗结构；渗流；渗透坡降；瀑布沟水电站　中图分类号：ＴＵ　４４３　文献标志码：Ａ　文章编号：１００９—６４０Ｘ（２０１３）０５－００７７－０６　对于修建在深厚覆盖层的土质心墙坝而言，主要是通过坝基防渗墙和土质心墙防渗体的连接构成完整　的防渗体系．由于防渗墙的刚度较之覆盖层要大许多，覆盖层的压缩沉降要比防渗墙大得多，防渗墙与覆　盖层在垂直方向上的变形差别必将在防渗墙侧面产生较大的摩擦力．同时由于防渗墙一廊道对心墙的顶托，　使廊道附近的心墙土体内出现应力场和位移场的畸变，产生应力集中，从而出现剪切破坏，并导致心墙土　料出现裂缝，从而严重降低心墙的抗渗性能．工程实践表明，防渗墙与坝体土质防渗体的连接是大坝一地基　整个防渗系统中最薄弱的环节…．因此，防渗墙的上部，尤其是与大坝土质防渗体连接处，是渗流控制的关　键部位．上游防渗墙渗流控制的失效，往往造成大坝防渗体的渗透破坏一ｚ　．　防渗墙与大坝土质心墙的连接目前主要有两种形式：一种是墙顶设廊道的形式；另一种是防渗墙直接　插入土质防渗体一定深度的形式．墙顶设置廊道，墙中问预留灌浆孔，可以在墙体局部失效时进行补强灌浆．　插入式连接由于插入段墙外侧向土压力的作用，墙与填土之间能够紧密结合，不容易发生渗流接触冲刷，因　而是一种结构简单而且效果较好的连接方式．由于坝高的发展和防渗墙允许渗流梯度的限制，一些深厚覆盖　层地基上的高水头土石坝采用了两道防渗墙，如瀑布沟砾石土心墙土石坝等．虽然目前土石坝渗流方面的相　关文献［３—９］已有一些发表，但对于设置两道防渗墙与土质心墙底部的渗流性状进行分析的文章较为鲜见．　因此，以瀑布沟砾石土心墙土石坝为例，分析了瀑布沟水电站心墙土石坝两道防渗墙与土质心墙底部的渗流　性状，可以为今后类似工程两道防渗墙与心墙底部连接处防渗结构设计，以及坝体和坝基防渗设计优化提供　参考，具有十分重要的理论价值和工程意义．　１　工程概况　瀑布沟水电站坝址位于大渡河中游尼日河汇口上游觉托附近，地跨四川省西部汉源县和甘洛两县境内，　下距乌斯河镇７　ｋｍ，距成昆铁路汉源车站９　ｋｍ，上距汉源县城２８　ｋｍ．电站采用坝式开发，是一座以发电为　主，兼有防洪、拦沙等综合利用效益的大型水电站．瀑布沟水电站采用心墙土石坝坝型，最低高程６７０　ｍ，坝　收稿日期：２０１３—０４—０８　作者简介：涂扬举（１９６４一），男，四川邛崃人，教授级高级工程师，主要从事水电工程建设管理工作　通信作者：薛新华（Ｅ—ｍａｉｌ：ｓｃｕｘｘｈ＠１６３．ｔｏｍ）　７８　水　利　水　运　工　程　学　报　２０１３年１０月　高１８６　ｍ，坝顶轴线长５７３　ｍ，宽１４　ｍ，坝体最大底宽约７８０　ｍ．坝顶上游侧设混凝土防浪墙，墙顶高程　８５７．２０　ｍ，大坝上游７９５　ｍ高程设５　ｍ宽马道，马道以上坝坡１：２．０，以下坝坡１：ｌ２．２５．大坝下游坝坡布置　有“之”字形道路，道路布置在１：１．８坡比之外．大坝心墙部分主要是由砾石土防渗体和人工反滤体系组　成，心墙砾石土防渗体的顶高程为８５４　ｍ，顶宽４　ｍ，底高程６７０　ｍ，上、下游坡均为１：０．２５．土料以宽级配砾　石土为主，心墙底部、心墙与岸坡接触带、防渗墙顶和混凝土廊道周围设高塑性黏土．心墙上、下游侧各设两　层反滤，上游两层厚均为４　ｍ，下游两层厚均为６　ｍ，反滤层以外为过渡料和坝壳料．心墙下部河床覆盖层采　用两道各厚１．２　ｍ、间隔１４　ｍ的混凝土防渗墙防渗，其中一道为主墙，一道为副墙，心墙底面以下最大墙深　约７８　ｍ．基础防渗墙与心墙间采用插入式连接，防渗墙插入心墙深度为１５　ｍ，墙间设混凝土廊道，瀑布沟水　电站心墙土石坝剖面见图１．　图１　瀑布沟心墙土石坝剖面　Ｆｉｇ．１　Ｐｒｏｆｉｌｅ　ｏｆ　Ｐｕｂｕｇｏｕ　ｅａｒｔｈ—ｒｏｃｌｄｉｌｌ　ｄａｍ　２　坝基渗流计算与防渗墙的水头分担比例确定　技施阶段委托清华大学和中国科学院力学研究所对瀑布沟心墙堆石坝进行了渗流分析　．利用自编有　限元程序对瀑布沟砾石土心墙坝进行了渗流计算．计算主要分为两个阶段，第一个阶段主要是为了满足心墙　底部局部渗流分析与应力变形分析要求，确定两道防渗墙的水头分担比例而取局部区域进行的渗流计算；第　二个阶段是完成优化论证以后，对整个坝体与坝基进行全面渗流分析．　２．１计算模型和计算网格　坝基渗流计算的主要目的是确定两道防渗墙的水头分担比例．由于坝肩上部的渗流对覆盖层内的渗流　性状影响很小，渗流的上边界只取到高程７５０　ｍ．从心墙内部进入防渗墙之间的渗流量远小于坝基岩体绕渗　的渗流量，因而坝体的渗流对坝基的流场影响可以忽略，在坝基　渗流计算中可以不包含坝体部分，以减小计算的规模，提高模型　的精细程度．　图２为坝基渗流计算的有限元网格，计算的范围取为４５０　至７５０　ｍ高程，顺河向取坝轴线上下游各２００　ｍ，坝轴向取约　９００　ｍ．网格共计１２　９７６个２０节点六面体单元，５６　４９２个节点．　计算所选取的边界条件为：上游覆盖层表面节点、两岸山体表面　和模型上游断面的覆盖层表面节点作为上游已知水头边界；下　图２坝基渗流计算有限元网格　Ｆｉｇ．２　３Ｄ　ｍｅｓｈｅｓ　ｏｆ　ａ　ｍｏｄｅｌ　第５期　涂扬举，等：瀑布沟高心墙土石坝渗流分析　７９　游覆盖层表面节点和模型下游断面的覆盖层表面节点作为下游溢出面边界；其他边界均为不透水边界．计算　上游水位为８５０　ｍ，下游水位６７０　ｍ时坝基的稳定渗流场．　２．２计算参数和计算条件　图３为防渗墙剖面基岩渗流参数与帷幕布置，副防渗墙的帷幕灌浆范围在０＋２３０．００以左取最大深度往　下２０　ｍ，以右取最大深度往下１０　ｍ．主防渗墙帷幕灌浆范围依照工程布置图中大坝剖面图的灌浆范围确定．　１ｎ　Ｔ】ｌ　１ｏ￣１ｍ　Ｔｌｌ１　１　１ｏｏ　Ｔｌｌ　１３　ｌｎ　Ｔ】ｌ　１ｎＩ】ｌ　１　１ｍ　Ｔ　Ｊ】　１０￣１ｍＩ】ｌ　３　１０　Ｔ】ｌ　图３渗流参数与帷幕布置　Ｆｉｇ．３　Ｐａｒａｍｅｔｅｒｓ　ｏｆ　ｓｅｅｐａｇｅ　ａｎｄ　ｃｕｒｔａｉｎ　ａｒｒａｎｇｅｍｅｎｔ　河床覆盖层各层均属强透水层，渗透系数差异不大，在０．０２３—０．１０４　ｃｍ／ｓ范围内，各层中均有局部架　空层分布，架空层渗透系数　值为０．１１６～０．５８０　ｃｍ／ｓ．覆盖层材料按同一参数取值．灌浆帷幕厚度取为　２．６　ｍ．渗流计算采用３组参数，见表１．第１组参数为基本参数，第２组参数防渗墙的渗透系数降低，用于分　析防渗墙参数对水头分担结果的敏感性，第３组参数岩石的渗透系数略为放大，分析岩石参数取值偏差的敏　感性．　表１各材料渗透系数　Ｔａｂ．１　Ｓｅｅｐａｇｅ　ｐａｒａｍｅｔｅｒｓ　ｏｆ　ｄｉｆｆｅｒｅｎｔ　ｍａｔｅｒｉａｌｓ　（ｏｎ・ｓ　）　２．３计算结果分析　３组参数计算的各部位水头消耗比例分别为：第１组，副防渗墙４１．５％，主防渗墙５８．３％，覆盖层０．２％；　第２组，副防渗墙３８．６％，主防渗墙６１．４％，覆盖层０．２％；第３组，副防渗墙４０．１％，主防渗墙５９．５％，覆盖　层０．４％．可见，副墙上消耗的水头比例为４０％左右，主墙上为６０％左右，河床覆盖层上消耗的水头比例可　以忽略不计；在岩石和灌浆帷幕的渗透系数不变的情况下，防渗墙的渗透系数降低为原来的９０％，主防渗墙　承担的水头比例有３％的增长，在防渗墙和帷幕灌浆的渗透系数不变的情况下，岩石的渗透系数取偏大值，　主防渗墙承担的水头比例有１％的增长，故主副防渗墙的水头分担比例基本上为６：４．　综上所述，坝体渗流与应力变形计算时副防渗墙按４０％承担水头来计算防渗墙的应力变形是合适的．　３　坝体与坝基整体三维渗流计算　３．１计算模型与计算网格　如图４所示，整体网格的计算范围左岸截断处断面　距离河谷中心线（坝轴线处桩号０＋２８６．００）８７０　ｍ左右；　右岸截断处距离河谷中心线６００　ｍ左右．上下游断面均　图４三维有限元计算网格　垂直于河流方向，上游断面河谷中心处距离坝轴线　Ｆｉｇ．４　３Ｄ　ｍｅｓｈｅｓ　ｏｆ　ｔｈｅ　ｍｏｄｅｌ　８０　水　利　水　运　工　程　学　报　２０１３年１０月　８５０　ｍ，下游断面河谷中心处距离坝轴线８７０　ｍ．底部断面为高程４００　ｍ的平面．共剖分６２　９０４个２０节点六　面体单元，２６６　１１８个节点，坐标系沿坝轴线桩号增加方向为　轴正方向，垂直坝轴线平面指向下游为Ｙ轴正　方向，沿高程增加方向为ｚ轴正方向　加］．　３．２河谷两侧地下水位反演　根据水文地质资料，左岸ＺＫ４６孑Ｌ年水位动态变化８７３．６５～８９４．０３　ｍ，右岸ＺＨ４７孔年水位动态变化范　围为８９３．４５～８９９．５７　ｍ，左岸平均水力坡度为１４。一１７。，　右岸平均水力坡度为２２。～２５。，两岸地下水向河谷中间　排泄，由于河道在折弯段前后基本顺直，可近似认为天然　地下水的渗流状况为垂直河流方向的二维渗流．通过坝　轴线河谷横剖面的渗流场计算来对三维模型两岸截断处　的地下水位进行反演．通过改变两侧边界水头进行试算，　最终确定的左侧边界水头为９００　ｍ，右侧边界水头为　１　０３０　ｍ，得到的地下水等孔压线如下图５所示．　３．３计算参数和计算条件　／ｍ　图５边界水头反演渗流场　Ｆｉｇ．５　Ｂｏｕｎｄａｒｙ　ｈｅａｄ　ｉｎｖｅｒｓｉｏｎ　ｓｅｅｐａｇｅ　ｆｉｅｌｄ　图６为防渗墙剖面的材料分区与帷幕布置，副防渗墙的帷幕灌浆范围在０＋２３０．００　ｍ以左取最大深度往　下２０　ｍ，以右取最大深度往下１０　ｍ．主防渗墙帷幕灌浆范围依照工程布置图中大坝剖面图的灌浆范围确定．　（ａ）主防渗墙剖面　（ｂ）副防渗墙剖面　图６材料分区与帷幕布置　Ｆｉｇ．６　Ｚｏｎｉｎｇ　ｏｆ　ｍａｔｅｒｉａｌｓ　ａｎｄ　ｃｕｒｔａｉｎ　ａｒｒａｎｇｅｍｅｎｔ　整体渗流计算网格的材料参数见表２，各河床覆盖层均属强透水层，渗透系数差异不大，在０．０２３～　０．１０４　ｃｍ／ｓ范围内，故按同一参数取值．计算工况为稳定渗流工况，上游水位８５０　ｍ，下游水位６７０　ｍ．　表２坝基渗流计算参数　Ｔａｂ．２　Ｓｅｅｐａｇｅ　ｐａｒａｍｅｔｅｒｓ　ｏｆ　ｄａｍ　ｆｏｕｎｄａｔｉｏｎ　材料编号　①　材料种类　１０～１００　Ｌｕ　渗透系数　２５×１０一　．（ｃｍ・ｓ　）　渗透系数　２０ｘ１０—５　．材料编号　⑧　材料种类　固结灌浆　②　③　３～１０　Ｌｕ　１　３　Ｌｕ　５　０ｘ１０一　２　０ｘ１０—５　⑨　⑩　⑥　⑩　帷幕灌浆　防渗墙　１０ｘ１０—５　．１０ｘｌ０—６　．④　⑤　＜１　Ｌｕ　崩积堆积１０ｘ１０—５　．心墙料反滤层　堆石　１ｏ×１０—５　．８０ｘ１０—２　．５０ｘｌ０—３　．⑥　⑦　Ｆ２断层　覆盖层１０ｘ１０—２　８０ｘ１０—２　．．１０ｘ１０　．３．４计算结果分析　通过对瀑布沟工程区水位地质资料和钻孔观测水位等资料的分析研究，建立了瀑布沟工程区天然渗流　场计算分析模型，并通过多组切取边界水位的模拟计算分析，拟合出一组与“观测渗流场”较为相近的三维　天然渗流场及其相应切取边界的地下水位边界条件，为三维渗流控制分析研究提供了基础．计算分析表明：　第５期　涂扬举，等：瀑布沟高心墙土石坝渗流分析　８１　（１）在现设计渗控方案及相应参数（防渗帷幕渗透系数１．０ｘ１０　ｃｍ／ｓ）条件下，左岸岩石外表面混凝土　板下游溢出点高程６８１．８　ｍ，右岸岩石外表面混凝土板下游溢出点高程６８１．２　ｍ，略高于河谷中心水面线，在　河床部位，坝体最大断面心墙下游反滤层后堆石体内水面线高程为６８０．８　ｍ．在考虑边界水头的情况下，防　渗墙处的水头降落要略小于无边界水头的情况．副墙上的水头降落约占总水头的３８．５％，主墙上的水头降　落约占５４．６％，两道防渗墙水头分担比例约为４：６，整体的计算结果与局部计算结果基本一致．在Ｆ２断层　以及左岸坝肩的崩积岩体处由于渗透系数较大，水头有明显的下降，覆盖层内的水头降落较小．　（２）在设计初拟渗控方案及相应参数（防渗帷幕渗透系数１．０×１０　ｃｍ／ｓ）条件库水通过防渗平面向下　游的渗透量为１３　００９　ｍ　／ｄ，其中透过坝体渗流量为８６２　ｍ　／ｄ，０＋２８６　ｍ以左的渗流量为１０　１７６　ｍ　／ｄ，以右　的渗流量为１　９７１　ｍ　／ｄ．　（３）心墙上部的渗透坡降较小，下部较大，心墙７１０　ｍ高程以上的渗透坡降一般不大于２．５．心墙内部高　塑性黏土范围以外一般渗透系数小于４，而廊道和副墙顶部的高塑性黏土渗透坡降较大，最大值达到１６．廊　道与副墙之间高塑性黏土下部的渗透坡降小于１．河谷中心部分心墙下游侧底部与覆盖层相接部位渗透坡　降小于４．５，而与岸坡相接处渗透坡降比较大，最大值达到１０．７．心墙下游面的下部渗透坡降介于２～２．５之　间，下游面出口渗透坡降是很安全的．心墙底部与固结灌浆后的覆盖层接触部位渗透坡降由廊道向下游逐步　减小，渗透坡降在河谷中部也比较小，最大值高塑性黏土内小于５，砾质土内小于３．５，在坡脚比较大．由于底　部设置了０．５　ｍ厚的水泥土，接触面又处于受压状态，土颗粒进入灌浆覆盖层和接触冲刷的可能性不大．　（４）覆盖层内的渗透坡降较小，约为０．０１，小于覆盖层的准许渗透坡降（０．０７—０．４）．最大坡降位于主防　渗墙底部，约为６．７，副防渗墙底部最大坡降约为５．５，主防渗墙顶部与坝体连接处的坡降也较大，约为４．５．　（５）左岸侧由于１０～１００　Ｌｕ岩体范围较大，且有　断层发育，渗流速度较右岸侧大．最大流速位于防渗　帷幕范围以外的断层内部，达到１４　ｍ／ｄ，此处是渗漏的主要通道．因此，左岸侧０＋１２６～０＋２８６　ｍ范围内的防　渗帷幕适当加深，以能覆盖住断层的影响范围为宜．　，　４　结　语　深厚覆盖层中防渗墙与土石坝防渗体的连接往往是大坝防渗体系的关键因素，也是防渗体系的薄弱环　节，对防渗墙的应力状态和安全具有决定性影响，因此，对于设置两道防渗墙与土质心墙底部的渗流性状进　行分析计算是十分必要的．本文通过对瀑布沟砾石土心墙土石坝防渗墙与坝基渗流等进行了分析，可以得出　以下结论：　（１）坝体渗流与应力变形计算时，副防渗墙按４０％承担水头来计算防渗墙的应力变形是比较合适的．　（２）连接部位的渗透坡降是非均匀变化的，防渗墙与廊道顶部的渗透坡降大，渗流出口的渗透坡降小．　（３）设两道防渗墙时，混凝土结构插入土质防渗体的深度可以通过有限元计算进行优化设计．　（４）连接处土体和混凝土结构接触面一般处于紧压闭合状态，高塑性黏土仅设于混凝土结构顶部即可　满足防渗要求．　参考文献：　［１］吴梦喜，余学明，叶发明．高心墙堆石坝坝基防渗墙与心墙连接方案研究［Ｊ］．长江科学院院报，２０１０，２７（９）：５９－６３．　（ｗＵ　Ｍｅｎｇ－ｘｉ，ＹＵ　Ｘｕｅ・ｍｉｎｇ，ＹＥ　Ｆａ—ｍｉｎｇ．Ｒｅｓｅａｒｃｈ　ｏｎ　ｃｏｎｊｏｉｎｉｎｇ　ｓｃｈｅｍｅ　ｆｏｒ　ａｎｔｉ—ｓｅｅｐａｇｅ　ｗａｌｌ　ｏｆ　ｈｉｇｈ　ｒｏｃｋｆｉｌｌ　ｄａｍ　ｗｉｔｈ　ｅａｒｔｈ　ｃｏｒｅ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆＹａｎｇｔｚｅ　Ｒｉｖｅｒ　Ｓｃｉｅｎｔｉｉｆｃ　Ｒｅｓｅａｒｃｈ　Ｉｎｓｔｉｔｕｔｅ，２０１０，２７（９）：５９—６３．（ｉｎ　Ｃｈｉｎｅｓｅ））　［２］刘杰．混凝土防渗墙渗流控制几个问题实例分析［Ｊ］．大坝观测与土工测试，１９９２，１６（５）：３９－４３．（ＬＩＵ　Ｊｉｅ．Ｃａｓｅ　ｓｔｕｄｙ　ｆｏ　ｃｏｎｃｒｅｔｅ　ｄｉａｐｈｒａｇｍ　ｗａｌｌｓ　ｓｅｅｐａｇｅ　ｃｏｎｔｒｏｌ［Ｊ］．Ｄａｍ　Ｏｂｓｅｒｖａｔｉｏｎ　ａｎｄ　Ｇｅｏｔｅｃｈｎｉｃａｌ　Ｔｅｓｔｓ，１９９２，１６（５）：３９－４３．（ｉｎ　Ｃｈｉｎｅｓｅ））　［３］郭成谦．论双排防渗墙的设计［Ｊ］．水力发电，１９９２（９）：１８．２２．（ＧＵＯ　Ｃｈｅｎｇ—ｑｉａｎ．Ｏｎ　ｔｈｅ　ｄｅｓｉｇｎ　ｏｆｔｗｏ　ａｎｔｉ．ｓｅｅｐａｇｅ　ｗａｌｌｓ　［Ｊ］．Ｗａｔｅｒ　Ｐｏｗｅｒ，１９９２（９）：１８—２２．（ｉｎ　Ｃｈｉｎｅｓｅ））　［４］肖白云．混凝土防渗墙墙体材料及接头型式的研究［Ｊ］．水力发电，１９９８（３）：２９—３１．（ＸＩＡＯ　Ｂａｉ—ｙｕｎ．Ｒｅｓｅａｒｃｈ　ｏｎ　８２　水　利　水　运　工　程　学　报　２０１３年ｌ０月　ｍａｔｅｒｉａｌ　ａｎｄ　ｊｏｉｎｔ　ｔｙｐｅ　ｏｆ　ａｎｔｉ－ｓｅｅｐａｇｅ　ｗａｌｌ［Ｊ］．Ｗａｔｅｒ　Ｐｏｗｅｒ，１９９８（３）：２９－３１．（ｉｎ　Ｃｈｉｎｅｓｅ））　［５］陈刚，马光文，付兴友，等．瀑布沟大坝基础防渗墙廊道连接型式研究［Ｊ］．四川大学学报：工程科学版，２００５，３７（３）：　３２－３６．（ＣＨＥＮ　Ｇａｎｇ，ＭＡ　Ｇｕａｎｇ—ｗｅｎ，ＦＵ　Ｘｉｎｇ—ｙｏｕ，ｅｔ　ａ１．Ｒｅｓｅａｒｃｈ　ｆｏｒ　ｔｈｅ　ｊｏｉｎｔ　ｔｙｐｅ　ｂｙ　ｇａｌｌｅｒｙ　ｂｅｔｗｅｅｎ　ｄａｍ　ｉｍｐｅｉｒｏｕｓ　ｗａｌｌ　ａｎｄ　ｃｏｒｅ　ｗａｌｌ　ｏｆ　Ｐｕｂｕｇｏｕ　ｐｒｏｊｅｃｔ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｓｉｃｈｕａｎ　Ｕｎｉｖｅｒｓｉｔｙ（Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ　Ｅｄｉｔｉｏｎ），２００５，３７（３）：３２—３６．（ｉｎ　Ｃｈｉｎｅｓｅ））　［６］傅少君，陈胜宏．瀑布沟堆石坝防渗体自适应有限元分析［Ｊ］．岩土力学，２００６，２７（３）：４９９－５００．（ＦＵ　Ｓｈａｏ—ｊｕｎ，ＣＨＥＮ　Ｓｈｅｎｇ—ｈｏｎｇ．Ａｎａｌｙｓｉｓ　ｏｆ　ｓｅｅｐａｇｅ　ｐｒｅｖｅｎｔｉｏｎ　ｓｙｓｔｅｍ　ｏｆ　Ｐｕｂｕｇｏｕ　ｒｏｃｋｆｉｌｌ　ｄａｍ　ｂｙ　ａｄａｐｔｉｖｅ　ＦＥＭ［Ｊ］．Ｒｏｃｋ　ａｎｄ　Ｓｏｉｌ　Ｍｅｃｈａｎｉｃｓ，　２００６，２７（３）：４９９－５０４．（ｉｎ　Ｃｈｉｎｅｓｅ））　［７］吴秋军，傅少君．子模型方法研究瀑布沟土石坝防渗结构［Ｊ］．武汉大学学报（工学版），２００６，３９（３）：５５－５９．（ＷＵ　Ｑｉｕ—　ｊｕｎ，ＦＵ　Ｓｈａｏ—ｊｎｎ．Ｓｔｕｄｙ　ｏｎ　ｓｅｅｐａｇｅ　ｐｒｅｖｅｎｔｉｏｎ　ｓｔｒｕｃｔｕｒｅ　ｆｏｒ　Ｐｕｂｕｇｏｕ　ｅａｒｔｈ—ｒｏｃｋ　ｄａｍ　ｂｙ　ｓｕｂ—ｍｏｄｅｌ　ｍｅｔｈｏｄ［Ｊ］．Ｅｎｇｉｎｅｅｒｉｎｇ　Ｊｏｕｒｎａｌ　ｏｆ　Ｗｕｈａｎ　Ｕｎｉｖｅｒｓｉｔｙ，２００６，３９（３）：５５—５９．（ｉｎ　Ｃｈｉｎｅｓｅ））　［８］邢玉玲，束一鸣，花加凤，等．膜土联合防渗系统对高土石坝心墙拱效应的影响［Ｊ］．河海大学学报：自然科学版，　２００７，３５（４）：４５６４５９．（ＸＩＮＧ　Ｙｕ—ｌｉｎｇ，ＳＨＵ　Ｙｉ—ｍｉｎｇ，ＨＵＡ　Ｊｉａ—ｆｅｎｇ，ｅｔ　ａ１．Ｅｆｆｅｃｔ　ｏｆ　ｇｅｏｍｅｍｂｒａｎｃｅ　ａｎｄ　ｃｏｒｅ—ｗａｌｌ　ｕｎｉｔｅｄ　ｉｍｐｅｒｉｏｕｓ　ｓｙｓｔｅｍ　ｏｎ　ａｒｃｈ　ｅｆｆｅｃｔ　ｏｆ　ｃｏｒｅ—ｗａｌｌ　ｉｎ　ｈｉｇｈ　ｅａｒｔｈ—ｒｏｃｋ　ｄａｍｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｈｏｈａｉ　Ｕｎｉｖｅｒｉｓｉｔｙ（Ｎａｔｕｒａｌ　Ｓｃｉｅｎｃｅ），２００７，　３５（４）：４５６４５９．（ｉｎ　Ｃｈｉｎｅｓｅ））　［９］段祥宝，李祖贻．瀑布沟水电站大坝三维渗流数值模拟研究［Ｊ］．水电站设计，１９９７，１３（１）：２９—３８．（ＤＵＡＮ　Ｘｉａｎｇ—ｂａｏ，　ＬＩ　Ｚｕ－ｙｉ．Ｎｕｍｅｒｉｃａｌ　ｓｉｍｕｌａｔｉｏｎ　ｒｅｓｅａｒｃｈ　ｏｎ　ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌ　ｓｅｅｐａｇｅ　ｆｌｏｗ　ｆｏｒ　Ｐｕｂｕｇｏｕ　ｅａｒｔｈ　ｒｏｃｋｆｉｌｌ　ｄａｍ［Ｊ］．Ｄｅｓｉｇｎ　ｏｆ　Ｈｙｄｒｏｅｌｅｃｔｒｉｃ　Ｐｏｗｅｒ　Ｓｔａｔｉｏｎ，１９９７，１３（１）：２９—３８．（ｉｎ　Ｃｈｉｎｅｓｅ））　［１０］成都勘测设计研究院．四川省大渡河瀑布沟水电站枢纽工程蓄水安全鉴定设计自检报告［Ｒ］．成都：成都勘测设计研　究院，２００９．（Ｃｈｅｎｇｄｕ　Ｈｙｄｒｏｅｌｅｃｔｒｉｃ　Ｉｎｖｅｓｔｉｇａｔｉｏｎ＆Ｄｅｓｉｇｎ　Ｉｎｓｔｉｔｕｔｅ．Ｋｅｙ　ｐｒｏｊｅｅｔ　ｓｔｏｒａｇｅ　ｓａｆｅｔｙ　ａｐｐｒａｉｓａｌ　ｄｅｓｉｇｎ　ｓｅｌ￣ｅｈｅｅｋ　ｒｅｐｏｒｔ　ｏｆ　Ｐｕｂｕｇｏｕ［Ｒ］．Ｃｈｅｎｇｄｕ：．Ｃｈｅｎｇｄｕ　Ｈｙｄｒｏｅｌｅｃｔｒｉｃ　Ｉｎｖｅｓｔｉｇａｔｉｏｎ＆Ｄｅｓｉｇｎ　Ｉｎｓｔｉｔｕｔｅ，２００９．（ｉｎ　Ｃｈｉｎｅｓｅ））　Ｓｅｅｐａｇｅ　ａｎａｌｙｓｉｓ　ｏｆ　Ｐｕｂｕｇｏｕ　ｈｉｇｈ　ｅａｒｃｈ－ｒｏｃｋｆｉｌｌ　ｄａｍ　ＴＵ　Ｙａｎｇ．ｊｕ　，ＷＡＮＧ　Ｗｅｎ．ｔａｏ　，ＸＵＥ　Ｘｉｎ．ｈｕａ　（１．Ｇｕｏｄｉａｎ　Ｘｉｎｉｆａｎｇ　Ｐｏｗｅｒ　Ｃｏ．，Ｌｔｄ．，Ｗｕｌｕｍｕｑｉ　８３００６３，Ｃｈｉｎａ；２．Ｄａｄｕ　Ｒｉｖｅｒ　Ｈｙｄｒｏｐｏｗｅｒ　Ｄｅｖｅｌｏｐｍｅｎｔ　Ｃｏ．，Ｌｔｄ．，Ｃｈｅｎｇｄｕ　６１００４１，Ｃｈｉｎａ；３．Ｃｏｌｌｅｇｅ　ｏｆ　Ｗａｔｅｒ　Ｒｅｓｏｕｒｃｅ　ａｎｄ　Ｈｙｄｒｏｐｏｗｅｒ，Ｓｉｃｈｕａｎ　Ｕｎｉｖｅｒｓｉｔｙ，　Ｃｈｅｎｇｄｕ　６１００６５，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｒｅｌｉａｂｉｌｉｔｙ　ｏｆ　ｓｅｅｐａｇｅ　ｐｒｅｖｅｎｔｉｏｎ　ｓｙｓｔｅｍ　ｉｓ　ａ　ｋｅｙ　ｐｒｏｂｌｅｍ　ｆｏｒ　ｈｉｇｈ　ｅａｒｔｈ—ｒｏｃｋｆｉｌｌ　ｄａｍ　ｌｏｃａｔｅｄ　ｉｎ　ａ　ｆｏｕｎｄａｔｉｏｎ　ｏｆ　ｌａｒｇｅ　ｏｖｅｒｂｕｒｄｅｎ　ｄｅｐｔｈ．Ｔｈｅ　ｃｏｎｊｏｉｎｉｎｇ　ａｒｅａ　ｏｆ　ｔｈｅ　ａｎｔｉ—ｓｅｅｐａｇｅ　ｓｏｉｌ　ｂｏｄｙ　ａｎｄ　ｔｈｅ　ｃｏｎｃｒｅｔｅ　ｄｉａｐｈｒａｇｍ　ｗａｌｌ　ｉｓ　ａ　ｋｅｙ　ｐａｒｔ　ａｇａｉｎｓｔ　ｓｅｅｐａｇｅ　ｉｎｄｕｃｅｄ　ｆａｉｌｕｒｅ．Ａｃｃｏｒｄｉｎｇ　ｔｏ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｔｈｅ　ｓｅｅｐａｇｅ　ｐｒｅｖｅｎｔｉｏｎ　ｓｙｓｔｅｍ　ｏｆ　ｔｈｅ　Ｐｕｂｕｇｏｕ　ｅａｒｔｈ－ｒｏｃｋｆｉｌｌ　ｄａｍ，ａ　ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ　ｍｅｔｈｏｄ（ＦＥＭ）ｉｓ　ａｐｐｌｉｅｄ　ｔｏ　ａｎａｌｙｚｅ　ｔｈｅ　ｓｅｅｐａｇｅ　ｆｌｏｗ　ｏｆ　ｔｈｅ　ｗｈｏｌｅ　ｄａｍ．Ｔｈｅ　ａｎａｌｙｓｉｓ　ｒｅｓｕｈｓ　ｓｈｏｗ　ｔｈａｔ　ｔｈｅ　ａｐｐｕｒｔｅｎａｎｔ　ｉｍｐｅｒｖｉｏｕｓ　ｗａｌｌ　ｔａｋｉｎｇ　４０％ｗａｔｅｒ　ｈｅａｄ　ｉｓ　ｓｕｉｔａｂｌｅ，ａｎｄ　ｔｈａｔ　ｔｈｅ　ｓｅｅｐａｇｅ　ｇｒａｄｉｅｎｔ　ｉｎ　ｔｈｅ　ｃｏｎｊｏｉｎｅｄ　ａｒｅａ　ｉｓ　ｉｎｈｏｍｏｇｅｎｅｏｕｓ　ａｒｏｕｎｄ　ｔｈｅ　ｃｏｎｃｒｅｔｅ　ｓｔｒｕｃｔｕｒｅ．Ｔｈｅ　ｓｅｅｐａｇｅ　ｇｒａｄｉｅｎｔ　ａｔ　ｔｈｅ　ｔｏｐ　ｏｆ　ｔｈｅ　ｓｔｒｕｃｔｕｒｅ　ｉｓ　ｍｕｃｈ　ｂｉｇｇｅｒ　ｔｈａｎ　ｔｈａｔ　ａｔ　ｔｈｅ　ｏｕｔｌｅｔ　ｏｆ　ｔｈｅ　ｅａｒｔｈ　ｃｏｒｅ　ｗａｌｌ　ｂａｓｅ．Ｔｈｅ　ａｎａｌｙｓｉｓ　ｒｅｓｕｌｔｓ　ｃａｎ　ｂｅ　ｔａｋｅｎ　ａｓ　ｔｈｅ　ｒｅｆｅ　ｒｅｎｅｅｓ　ｆｏｒ　ｓｉｍｉｌａｒ　ｐｒｏｊｅｃｔｓ（ｗｏｒｋｓ）ｉｎ　ｔｈｅ　ｆｕｔｕｒｅ．　Ｋｅｙ　ｗｏｒｄｓ：ｓｅｅｐａｇｅ　ｐｒｅｖｅｎｔｉｏｎ　ｓｔｕｒｃｔｕｒｅ；ｓｅｅｐａｇｅ；ｓｅｅｐａｇｅ　ｇｒａｄｉｅｎｔ；Ｐｕｂｕｇｏｕ　ｈｙｄｒｏｅｌｅｃｔｉｒｃ　ｐｒｏｊｅｃｔ