<p>&nbsp; &nbsp; &nbsp; &nbsp;GEO5土坡模块可以分析考虑水位骤降下的边坡稳定性，但在实际的使用中，有工程师反映不知道该怎么使用，也有人说输入了地下水之后边坡安全系数并没有发生改变，十分困惑，所以本文将对GEO5中分析水位骤降的方法进行详细说明。</p><p><strong>1、注意事项</strong></p><p>&nbsp; &nbsp; &nbsp; &nbsp;无论是分析水位骤降下边坡稳定性，还是一般情况下分析有地下水位的边坡稳定性，都需要注意的是在岩土材料输入的时候选择有效应力法进行计算，只有选择了有效应力法，软件才会考虑孔隙水压力对条块的作用。如果选择了总应力法或者总应力c_cu,φ_cu，软件都不会考虑坡内地下水位对边坡的影响，但坡外水位的有利作用软件还是会考虑。有效应力法和总应力法不同选择的区别可以查看<a href="http://www.wen.kulunsoft.com/a ... BGEO5中有效应力法、总应力法，水土分算、水土合算的说明</a>。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B alt="blob.png" width="435" height="415" style="width: 435px; height: 415px;"/></p><p><strong>2、传统分析方法</strong></p><p>&nbsp; &nbsp; &nbsp; &nbsp;传统分析水位骤降的方法是通过设置初始地下水位和水位骤降后的地下水位面来分析，最简单的做法是认为坡内的水来不及排出，那么水位骤降后坡内的水位保持不变，只改变坡外的静水面，随着水位的下降，边坡安全系数将逐渐降低。</p><p style="text-align: left;">&nbsp; &nbsp; &nbsp; &nbsp;在GEO5土坡模块中，选择【地下水】中的类型为“水位骤降”，可以直接定义边坡的初始地下水位和骤降后的地下水位：</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1584080821926652.png" alt="image.png"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;定义完成后，和一般的边坡计算一样直接进行分析即可。下图展示的是相同的初始地下水位，不同水位骤降情况的边坡安全系数。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1584087011538512.png" alt="image.png"/></p><p><strong>3、结合GEO5中的初始孔压折减系数分析</strong></p><p>&nbsp; &nbsp; &nbsp; &nbsp;传统的考虑坡内水来不及排出的方法实际上是一种偏保守的方法，因为水位骤降其实也是有一个过程的，那么坡内的水或多或少都会渗出坡外，如果是对于渗透性较好的土体，那么坡内的水位还会有明显的下降，但是针对这个问题，再去使用非稳定流分析浸润线就会显得有点麻烦。所以，GEO5通过巧妙地设置初始孔压折减系数X这样一个值，使得我们可以去考虑有水排出的情况。</p><p>&nbsp; &nbsp; &nbsp; &nbsp;当我们在【地下水】中选择的地下水类型为“水位骤降”时，需要在【岩土材料】中输入初始孔压折减系数的值：</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1584083072289482.png" alt="image.png" width="368" height="63" style="width: 368px; height: 63px;"/></p><p>这里X的取值范围为[0,1]，当土体完全透水时X=1，完全不透水时X=0，其他情况介于0和1之间，X值的作用原理可查看GEO5的帮助文档，或者直接点击<a href="http://www.wen.kulunsoft.com/a ... BGEO5土坡模块中地下水类型</a>。</p><p>&nbsp; &nbsp; &nbsp; &nbsp;这里需要对三种情况的取值进一步说明：</p><p>（1）X=1</p><p>&nbsp; &nbsp; &nbsp; &nbsp;X=1意味着土体完全透水，它的实际意义是：不考虑骤降后的水位与初始水位之间土体的孔隙水压力，所以X=1时，坡内不同的地下水位面会得到不同的结果。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1584085542948660.png" alt="image.png"/></p><p>（2）X=0</p><p>&nbsp; &nbsp; &nbsp; &nbsp;X=0意味着土体完全不透水，它的实际意义是：认为骤降后的水位与初始水位之间土体仍然处于饱和状态，所以X=0时，坡内不同的地下水位面会得到相同的结果。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1584085619188002.png" alt="image.png"/></p><p>（3）0&lt;X&lt;1</p><p>&nbsp; &nbsp; &nbsp; &nbsp;0&lt;X&lt;1其实模拟的是真实的情况，即水位骤降后考虑部分水的排出，既不是完全透水也不是完全不透水，在相同水位条件下，边坡安全系数将位于X=0和X=1之间。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1584086262289515.png" alt="image.png" width="346" height="211" style="width: 346px; height: 211px;"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;至于X如何取值，则需要根据实际岩土材料的渗透性以及水位骤降的速度和阶段综合选取。另外，通过以上分析，我们也不难发现，如果采用传统的通过控制坡内水位面不变化的方法来分析，那么X值无论设置为多少，对最终结果都没有影响。</p><p><br/></p>

<p>&nbsp; &nbsp; &nbsp; &nbsp;GEO5有限元渗流分析得到的浸润面可以直接导入到GEO5土坡模块中使用，这对于计算有地下水位的边坡稳定性十分方便。本文将简述操作方法及注意事项。</p><p>&nbsp; &nbsp; &nbsp; &nbsp;首先，将我们绘制的DXF文件以多段线形式导入到土坡模块中建立边坡模型，编辑好模型尺寸和材料参数后，复制模型数据。</p><p style="text-align:center"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583413604240387.png" alt="image.png"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;然后，在GEO5有限元模块中粘贴数据，建立和土坡模块相同的模型（尺寸相同、坐标不偏移），并在【分析设置】中选择分析类型为“稳定流”或“非稳定流”。<br/></p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583413664681364.png" alt="image.png"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;输入岩土材料的渗流参数，并生成网格。然后，在工况1当中定义线渗流边界条件，不同的线渗流边界的概念可查看<a href="http://www.wen.kulunsoft.com/d ... gt%3B。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583413924201413.png" alt="image.png"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;下一步直接进行渗流分析，得到如下图所示的浸润面，然后点击界面右侧“GEO剪贴板”中的复制计算地下水位。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583413962314702.png" alt="image.png"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;这样，浸润面就已经复制到了剪贴板当中。此时回到最初建好的土坡模块当中，在【地下水】中选择地下水类型为“地下水位”，并在右侧“GEO剪贴板”中粘贴地下水位。这样，有限元渗流分析得到的浸润面就直接导入到了土坡模块当中，接下来就可以进行有地下水位面的边坡稳定性分析。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583414017825248.png" alt="image.png"/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;需要注意的是，我们在有限元当中生成浸润面的时候，可能会出现下面这种奇怪的浸润面形态：</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583414101693454.png" alt="image.png" width="416" height="196" style="width: 416px; height: 196px;"/><br/></p><p>&nbsp; &nbsp; &nbsp; &nbsp;出现这种情况是因为下游水头高于了地形面，而整个坡面设置的线边界条件又都是溢出边界。由于溢出边界意味着该位置的孔隙水压力为0，所以在两个边界条件交接的位置会出现自相矛盾的情况。这个时候只需要根据下游的实际水位更改对应坡面的渗流边界条件即可。</p><p>&nbsp; &nbsp; &nbsp; &nbsp;解决方法是，在坡面对应位置添加一个自由点：</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583414221457587.png" alt="image.png"/></p><p>然后重新生成网格，并将原来下部的溢出边界改为孔隙水压力边界：</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583414290646284.png" alt="image.png"/></p><p>最后就可以得到正常的浸润面，如下所示：</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1583414337737486.png" alt="image.png" width="412" height="210" style="width: 412px; height: 210px;"/></p><p><br/></p>

<p><strong>项目名称</strong>：某降水分析项目</p><p><strong>使用软件：</strong>GEO5岩土工程有限元分析</p><p><strong>设计方案：</strong>放置两个降水井，岩土材料为粉土。</p><p><strong>项目背景：</strong></p><p>&nbsp;<img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1505181133715953.png" alt="blob.png"/></p><p><strong>项目特点：</strong>此项目是位于国外，但是由国内某著名设计院设计的。</p><p><strong>软件优势：</strong>GEO5有限元渗流分析能简单快速的计算出渗流结果，软件同时支持18种语言与计算书，可直接中文界面下设计，计算书可选择合适的语言，无需单独翻译。</p><p><strong>过程</strong><strong>：</strong></p><p>&nbsp;<img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1505181142516064.png" alt="blob.png"/></p><p><strong>接触面</strong></p><table><tbody><tr class="firstRow"><td style="border-width: 1px; border-style: solid;"><p><strong>编号</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>位置</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>渗透性</strong></p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>1</p></td><td style="border-width: 1px; border-style: solid;"><p>编号5网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>dn&nbsp;= 700.0 mm, kn&nbsp;= 5.00E+04 m/天, ks&nbsp;= 5.00E+04 m/天</p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>2</p></td><td style="border-width: 1px; border-style: solid;"><p>编号6网格线</p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;"><p>dn&nbsp;= 700.0 mm, kn&nbsp;= 5.00E+04 m/天, ks&nbsp;= 5.00E+04 m/天</p></td></tr></tbody></table><p><strong>点渗流边界</strong></p><table><tbody><tr class="firstRow"><td style="border-width: 1px; border-style: solid;"><p><strong>编号</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>位置</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>渗流边界类型</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>参数</strong></p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>1</p></td><td style="border-width: 1px; border-style: solid;"><p>编号6网格点</p></td><td style="border-width: 1px; border-style: solid;"><p>孔隙水压力 - 水位坐标</p></td><td style="border-width: 1px; border-style: solid;"><p>z水位&nbsp;= 81.60 m</p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>2</p></td><td style="border-width: 1px; border-style: solid;"><p>编号8网格点</p></td><td style="border-width: 1px; border-style: solid;"><p>孔隙水压力 - 水位坐标</p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;"><p>z水位&nbsp;= 81.60 m</p></td></tr></tbody></table><p><strong>线渗流边界条件</strong></p><table data-sort="sortDisabled"><tbody><tr class="firstRow"><td rowspan="2" colspan="1" style="border-width: 1px; border-style: solid;"><p><strong>编号</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>线渗流边界条件</strong></p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;"><p><strong>位置</strong></p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;" rowspan="2" colspan="1"><p><b>位置</b></p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;" rowspan="2" colspan="1"><p><b>边界条件类型</b></p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;" rowspan="2" colspan="1"><strong>参数</strong></td></tr><tr><td style="word-break: break-all; border-width: 1px; border-style: solid;"><span style="line-height: 22.5px;"><b>新</b></span></td><td style="word-break: break-all; border-width: 1px; border-style: solid;"><strong>修改</strong></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>1</p></td><td style="border-width: 1px; border-style: solid;"><p>是</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>编号1网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>孔隙水压力边界</p></td><td style="border-width: 1px; border-style: solid;"><p>z水位&nbsp;= 119.00 m</p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>2</p></td><td style="border-width: 1px; border-style: solid;"><p>是</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>编号2网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>不透水边界</p></td><td style="border-width: 1px; border-style: solid;"><br/></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>3</p></td><td style="border-width: 1px; border-style: solid;"><p>是</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>编号3网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>孔隙水压力边界</p></td><td style="border-width: 1px; border-style: solid;"><p>z水位&nbsp;= 119.00 m</p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>4</p></td><td style="border-width: 1px; border-style: solid;"><p>是</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>编号4网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>不透水边界</p></td><td style="border-width: 1px; border-style: solid;"><br/></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>5</p></td><td style="border-width: 1px; border-style: solid;"><p>是</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>编号7网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>不透水边界</p></td><td style="border-width: 1px; border-style: solid;"><br/></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>6</p></td><td style="border-width: 1px; border-style: solid;"><p>是</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>编号8网格线</p></td><td style="border-width: 1px; border-style: solid;"><p>不透水边界</p></td><td style="word-break: break-all; border-width: 1px; border-style: solid;"><br/></td></tr></tbody></table><p><strong>结果：</strong></p><table data-sort="sortDisabled"><tbody><tr class="firstRow"><td><p><strong>名称 : 分析</strong></p></td><td><p><strong>工况阶段 : 1</strong></p></td></tr><tr><td style="word-break: break-all;" height="65" rowspan="1" colspan="2"><p>结果 : 全量; 变量 : 孔隙水压力 u&nbsp;渗流; 范围 : &lt;0.00; 1994.00&gt; kPa</p><p>∑Q [m3/天/m]</p><p><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1505181211306436.png" alt="blob.png" width="486" height="239" style="width: 486px; height: 239px;"/><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1505181217365438.png" alt="blob.png"/><br/></p></td></tr></tbody></table><p><strong>计算总的流出量 /流入量</strong></p><table data-sort="sortDisabled"><tbody><tr class="firstRow"><td style="border-width: 1px; border-style: solid;" rowspan="2" colspan="1"><p><strong>位置</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>流入</strong></p></td><td style="border-width: 1px; border-style: solid;"><p><strong>流出边界</strong></p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p><strong>[m</strong><strong>3</strong><strong>/天/m]</strong></p></td><td style="border-width: 1px; border-style: solid; word-break: break-all;"><strong>[m</strong><strong>3</strong><strong>/天/m]</strong></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>点渗流边界条件编号1</p></td><td style="border-width: 1px; border-style: solid;"><p>476.927</p></td><td style="border-width: 1px; border-style: solid;"><br/></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>点渗流边界条件编号2</p></td><td style="border-width: 1px; border-style: solid;"><p>476.893</p></td><td style="border-width: 1px; border-style: solid;"><br/></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>线渗流边界条件编号1</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>-477.177</p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>线渗流边界条件编号3</p></td><td style="border-width: 1px; border-style: solid;"><br/></td><td style="border-width: 1px; border-style: solid;"><p>-476.644</p></td></tr><tr><td style="border-width: 1px; border-style: solid;"><p>总数</p></td><td style="border-width: 1px; border-style: solid;"><p>953.820</p></td><td style="border-width: 1px; border-style: solid;"><p>-953.820</p></td></tr></tbody></table><p><br/></p>

<p>　　GEO5有限元渗流模块计算不收敛时，应适当简化模型、结合实际经验调整相关参数可增加计算的收敛性。下面以近期某客户发来的项目为例，说明如何简化模型及调整哪些参数。</p><p>　　该项目为土石坝渗流分析，模型如下所示：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694317188428.png" alt="blob.png"/></p><p>　　加密网格后，显示如下：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694337416966.png" alt="blob.png"/></p><p>　　由于模型过于复杂，划分网格后现错误提示如下：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694355154217.png" alt="blob.png"/></p><p>　　在GEO5有限元中，当网格划分出现过多错误提示时，若不修改，可能会导致计算结果不收敛。因此需要适当的简化模型，针对本项目，简化之处有：</p><p>　　1、模型中有些地方较复杂，建议简化，以提高网格划分质量，如下图：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694378764216.png" alt="blob.png"/></p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694388109333.png" alt="blob.png"/></p><p>　　简化后，如下：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694408739290.png" alt="blob.png"/></p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694418135747.png" alt="blob.png"/></p><p>　　此时再划分网格，软件不再显示网格质量差的提示。</p><p>　　2、大坝最上方为混凝土墙，可以直接用不排水边界条件代替即可，简化模型。</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694437110412.png" alt="blob.png"/></p><p>　　重新划分完网格之后，此时软件不再有错误提示，显示如下：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694454378189.png" alt="blob.png"/></p><blockquote><p>注：针对复杂模型，建议地层线用dxf多段线导入，其他内部点线用dxf模板导入后用有限元中的自由点和自由线定义，这样方便后期修改，本模型中内部点线就是这样建成的。</p><p>GEO5多段线建模，<a href="/dochelp/19" target="_blank">点击这里</a>查看。</p></blockquote><p>　　修改完模型之后，计算结果还是不收敛，最后查明原因为相关渗透参数输入不正确。将渗透系数修改完之后计算收敛，结果如下图所示：</p><p style="text-align: center;"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1495694477664285.png" alt="blob.png"/></p><blockquote><p>注：关于渗透系数如何取值，<a href="/dochelp/912" target="_blank">点击这里</a>查看软件自带帮助中的相关说明。</p></blockquote><p>　　至此，关于有限元渗流模块中的注意事项至此结束，如有更好的想法欢迎在下方留言讨论。</p><p>案例源文件：<img src="http://www.wen.kulunsoft.com/s ... ot%3B style="line-height: 16px; vertical-align: middle; margin-right: 2px;"/><a href="http://wen.kulunsoft.com/uploa ... ot%3B title="有限元渗流稳定分析-简化- modified van Genuchten.rar" style="line-height: 16px; font-size: 12px; color: rgb(0, 102, 204);">有限元渗流稳定分析-简化- modified van Genuchten.rar</a></p>

<p>　　本工程案例为碎石土心墙堆石坝，坝高136m，坝顶宽12m，坝底宽71m。图1为坝体剖面图。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301658445588.png" alt="image.png"/></p><p style="text-align: center;">图1 坝体剖面图</p><p>　　首先在工况阶段[建模]中进行分析设置、添加岩土材料、导入几何模型及划分网格等。GEO5支持直接导入在CAD等软件中已经建好的DXF模型文件。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301672475580.png" alt="image.png"/></p><p style="text-align: center;">图2 导入坝体模型</p><p>　　在[分析设置]中选择分析类型为“稳定流”，勾选“详细结果”复选框。</p><p style="text-align: center;">&nbsp;<img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301684988915.png" alt="image.png"/></p><p style="text-align: center;">图3 分析设置</p><p>　　添加完所有的坝体材料以后，在[指定材料]界面中将材料指定给各自对应的部位。</p><p style="text-align: center;">&nbsp;<img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301694775477.png" alt="image.png"/></p><p style="text-align: center;">图4 指定坝体材料</p><p>　　接着通过创建自由线来确定混凝土防渗墙在地基中的位置，最后对模型进行网格划分操作。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301706437722.png" alt="image.png"/></p><p style="text-align: center;">图5 生成网格</p><p>　　点击进入工况阶段[1]，在[梁]设置界面中添加防渗墙。接着通过[线渗透边界]设置渗流边界条件，这里将上游水位高程设置为62.0m，下游边界类型设置为“溢出边界”。</p><p style="text-align: center;">&nbsp;<img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301718861785.png" alt="image.png"/></p><p style="text-align: center;">图6 边界条件设置</p><p>　　设置好边界条件以后，在[分析]界面中点击[开始分析按钮]，通过短暂的求解即可得到坝体的渗流特征计算结果。</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301730747915.png" alt="image.png"/></p><p style="text-align: center;">图7 孔隙水压力分布云图</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301742458576.png" alt="image.png"/></p><p style="text-align: center;">图8 总水头分布云图</p><p style="text-align: center;"><img src="http://www.wen.kulunsoft.com/u ... ot%3B title="1597301755684981.png" alt="image.png"/></p><p style="text-align: center;">图9 渗流矢量图</p>

<p style="text-align: justify;">我国许多大都市人口在数百万至数千万以上，城市人口高度集中，现代工业迅猛发展，对城市的环境地质影响作用也愈来愈强烈。我国幅员辽阔，国土南北跨越的纬度近50°，东西跨越的经度有60多度，大陆海岸线长度逾18000km。自然条件，如气候、地形、地貌、地质等差别显著，大多数城市所处的自然条件和地质环境比较复杂，城市发展中遇到的地质灾害各不相同。但几乎所有的地质灾害都或多或少与地下水有关。</p><p style="text-align: justify;">通过对地下水的动态、运动规律、物理化学特征等的研究，可以加深对地质灾害形成机理的认识，从而实现对地质灾害的预测预报和防止治理。</p><p style="text-align: justify;"><span style="color: #FF0000;"><strong>（一）城市地质灾害的类型</strong></span></p><p style="text-align: justify;">根据产生地质灾害的动力作用性质，可将地质灾害分为3类：</p><p style="text-align: justify;">Ⅰ类为由内动力地质作用产生的地质灾害；</p><p style="text-align: justify;">Ⅱ类为有外动力地质作用产生的地质灾害；</p><p style="text-align: justify;">Ⅲ类为由人类活动作用产生的地质灾害（如图1）。</p><p style="text-align: justify;">其中常发生的地质灾害主要有：</p><p style="text-align:center"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1489729412207241.png" alt="1489729412207241.png" width="480" height="308" border="0" vspace="0" style="width: 480px; height: 308px;"/></p><p style="text-align: center;">图1 地质灾害的类型</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>1、地震与地裂缝</strong></span></p><p style="text-align: justify;">我国人口在100万以上的大城市，70%位于地震裂度大于7度的地区内。我国是一个多地震的国家，8级以上的地震平均每10年1次，7级以上的地震平均每年1次，而5级以上的地震平均每年14次之多。我国地震活动强烈的地区，多分布在地壳不稳定的大陆板块和大洋板块接触带及板块断裂破碎带上，从地区分布上看主要是东南部的台湾和福建沿海，华北太行山沿线和京津唐地区，西南青藏高原及其边缘的四川、云南省西部，西北的新疆、甘肃和 宁夏。有资料记载以来，我国最大地震为8.5级，山东、西藏、宁夏各发生一次。1556年陕西华县发生8级大地震，死亡80多万人。1976年7月28日夜间，河北唐山发生7.8级大地震，有24万人丧生。</p><p style="text-align: justify;">20世纪60年代西安市在东南郊一带的小寨、 雁塔路、南沙坡村和秦川机械厂等地相继发现有地裂缝。1976年唐山发生7.8级大地震后，地裂缝活动日趋明显，特别是1977年西安北郊发生2.9级地震后，又陆续发现了西郊劳动公园、北郊八府庄和辛家庙等多处地裂缝，截止2002年经勘察确认的地裂缝已有13条。</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>2、地面沉降与塌陷</strong></span></p><p style="text-align: justify;">我国100万人口以上的大城市大约有30多个城市出现地面沉降区，还有一些大城市位于塌陷性黄土或胀缩性强的膨胀土地区、软土地区和岩溶区。</p><p style="text-align: justify;">由于过度强烈开采地下水、石油和天然气，平原地区的城市常出现地面沉降，位于地下隐伏岩溶发育地区的城市发生岩溶塌陷，如上海市、西安市、天津市和北京市等大都市形成了330~4000km²的沉降区或沉降洼地，其中西安市最大沉降量达2600mm以上，天津市最大沉降量达3090mm。沉降导致地面高程损失、雨后积水、市政设施破坏、河流泄洪能力下降、市区内河成为“地上悬河”、沿海风爆加剧、防汛设施的防御标准降低、土壤盐渍化等灾害。</p><p style="text-align: justify;">根据国土资源部的统计数字显示，40年来，仅上海市因地面沉降而造成的直接经济损失达2900亿元，其中潮损1755亿元、涝损848亿元、安全高程损失189亿元。</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>3、水环境恶化</strong></span></p><p style="text-align: justify;">城市“三废”排放增多，以及农业化肥和农药的大量使用，使区域水环境有恶化的趋势。</p><p style="text-align: justify;">如北京市全年有近500万t的建筑、工业和生活垃圾需要堆放、填埋、消纳，这些垃圾影响或污染着地表水体、地下水质、土壤地球化学背景和周边环境。由于地下水超采和人为活动因素干扰，北京市局部地区地下水硬度升高超标，并发现硝酸盐氮有面型污染趋势。</p><p style="text-align: justify;">西安市全年排放的各类垃圾达150万t，1990年前各类垃圾点有600余处，占地超过380万m²，分布在近郊城乡结合部位，严重污染水体、空气和土壤。 水资源匮乏和水环境恶化已成为制约我国许多大都市经济可持续发展的一个瓶颈。</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>4、滑坡、崩塌、泥石流等边坡失稳</strong></span></p><p style="text-align: justify;">位于丘陵山区或河流、海洋、湖泊附近的城市，滑坡、崩塌、泥石流等边坡失稳问题时有发生。</p><p style="text-align: justify;">在西安崩塌、滑坡、泥石流灾害隐患主要分布在秦岭北麓、骊山周边、黄河台塬边坡地带及横岭黄土丘陵区。</p><p style="text-align: justify;">上海市沿江沿海岸线，包括崇明、长兴、横沙等岛屿总长有460km，均为泥质海岸线。岸坡稳定性是一种缓变性地质灾害，它可由地震、地壳运动等内动力因素作用发生，也可由风化、气象、海洋、地表水、地下水等外动力因素和人为因素诱发。</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>5、流砂、管涌、软土变形</strong></span></p><p style="text-align: justify;">第四系地层中，砂质粉土和粉细沙层往往是隧道施工、深基坑开挖等地下工程最大的事故隐患。在上海城市建设中，流砂、管涌灾害比较普遍，如2003年7月1日凌晨上海轨道交通4号线发生大量流沙，引起地面大幅沉降，楼房倾斜、倒塌，堤防开裂，直接经济损失1.5亿元左右；又如上海市长寿路某大厦深基坑开挖，流砂引发塌方的场面不可收拾，工地停工处理事故长达1年。</p><p style="text-align: justify;">河口—滨海—浅海相成因的淤泥质粉质粘土和淤泥质粘土，一般为高压缩性土，是不良工程地质层，不可作为建筑物的持力层，如果处理不当，会给建筑物的地基质量留下隐患，引发事故。</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>6、其它问题</strong></span></p><p style="text-align: justify;">中国大陆海岸线长度逾18000km，由于沿海地区地下水严重超采，经常引起海水入侵问题。如山东省沿海共有19个县（市、区）发生海水入侵，面积超过1000km²，其中莱州市海水入侵面积达234km²。天津市、河北省秦皇岛市也有类似的海水入侵问题。</p><p style="text-align: justify;"><span style="color: #FF0000;"><strong>（二）地质灾害中的地下水环境效应</strong></span></p><p style="text-align: justify;">可以说无论内动力地质作用、外动力地质作用产生的地质灾害，还是人类活动作用产生的地质灾害，无不与地下水有关（如图2）。</p><p style="text-align:center"><img src="http://wen.kulunsoft.com/uploa ... ot%3B title="1489729508685324.png" alt="blob.png" style="line-height: 22.5px; white-space: normal;"/></p><p style="text-align: center;">图2 地下水与地质灾害的相关性图示</p><p style="text-align: justify;"><span style="color: #00B050;"><strong>1、地下水的作用形式</strong></span></p><p style="text-align: justify;">地下水作为地质灾害最敏感的触发因素之一，在地质灾害的萌发、发育、形成和发展过程中起着至关重要的作用，作用的形式可分为物理作用、化学作用和生物作用等。地下水对地质灾害的物理作用是最经常、最普遍的作用，贯穿于地质灾害的萌发、发育、形成和发展的全过程，但在不同的阶段，其作用的强弱有所区别。</p><p style="text-align: justify;">物理作用主要表现为地质体内的地下水，通过温度、物态变化和渗透、潜蚀作用，改变地质体的物理力学特性和受力状态，促进或影响各种地质灾害的萌发、发育、形成和发展。</p><p style="text-align: justify;"><span style="line-height: 1.5em;">化学作用主要表现为地质体内的地下水作为一种天然溶液，在渗透、潜蚀的同时，与矿物岩石发生各种化学反应，如氧化反应、溶解反应、水化反应、水解反应、碳酸化反应等，从而改变地质体的物理力学特性和受力状态，影响地质灾害的萌发、发育、形成和发展。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">生物作用是指生物在其生命活动中，对地质灾害的萌发、发育、形成和发展所起的作用。生物的生命活动离不开水，生物作用可以是机械的，也可以是化学的。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">2、地下水作用的效应</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">地下水对地质灾害的作用效应分为控制效应、辅助效应、次生效应。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">控制效应是指地下水作为一种控制性因素，在地质灾害的萌发、发育、形成和发展的全过程或某一阶段起控制性作用。控制效应大都是基于地下水的物理作用来实现的，如水库诱发地震，地下水往往是造成地应力集中，岩体破坏的直接原因，地下水几乎在滑坡的萌发、发育、形成和发展整个过程中都起着控制性作用，雨季产生的滑坡占总数的90%以上，以至于形成了“无水不成滑坡”的观点。众所周知，地下水的超采是导致地面沉降、海水入侵的最直接的原因；流砂、管涌、软土变形的发生，也离不开地下水。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">辅助效应是指地下水作为一种辅助性因素，在地质灾害的萌发、发育、形成和发展的全过程或某一阶段起催化剂的作用。辅助效应可以通过地下水的物理作用、化学作用或生物作用来实现的。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">次生效应指地下水作为一种次生因素，在地质灾害的萌发、发育、形成和发展的全过程或某一阶段，以某种物理、化学或生物现象表现出来。次生效应同样可以通过地下水的物理作用、化学作用或生物作用来实现。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">3、地下水的环境响应</strong><strong style="line-height: 1.5em;"></strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">由于发生地质灾害的地质环境、动力作用性质、类型各异，因此各种地质灾害在萌发、发育、形成和发展的全过程或某一阶段引起的地下水的环境响应是各不相同的。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">地震地质灾害中地下水的环境响应最为显著和多样化。大震前，有时天气大旱，但地下水（井水）却猛涨，甚至溢出地表；在多雨的季节里，井水本应逐渐上升，但却猛降，甚至干涸等。产生的原因是因地 震孕育过程中，地应力不断增强，尤其在震中区附近，因地应力的作用，地壳活动随之加强。压性区水位会逐渐抬升；张性区在张应力作用下，水位会逐渐下降。尤其是地壳局部区域在地应力作用下遭受破坏，发生变形或加速位移，以及由岩层破坏而引起上下层间水的贯通，都会使水位产生急剧的升降运动，这就是利用地下水预报地震的依据。地震时，由于地裂缝切过地下含水层，地下水受到挤压，并沿着裂缝夹带着泥砂涌出地表，即形成喷砂、冒水现象。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">城市地面沉降是摩天大楼的重量压塌了地面的结果吗？研究发现，引起地面沉降的最主要原因是大量开采地下水。2002年全国218个进行地下水动态监测城市和地区中，共有108个城市和地区的地下水水位出现不同程度的下降，占统计总数的49.6%，其中北方城市有66个。过量开采地下水形成 了规模不等的降落漏斗，华北平原的降落漏斗以天津、沧州、德州、衡水、唐山、廊坊、石家庄为中心，形成了大面积区域性承压水头下降区。总面积约14万km²的华北平原地区是世界上超采地下水最严重的地区，也是地下水降落漏斗和地面沉降面积最大、类型最复杂的地区。南方最大的降落漏斗是苏、锡、常地区，杭、嘉、湖平原地区以及上海市的区域水位降落漏斗相互迭交，形成了长江三角洲区域水位降落漏斗。目前，上海、嘉兴、宁波等地的沉降速度虽趋缓，但总体沉降范围却在迅速扩展，杭、嘉、湖的沉降正向整个平原蔓延，苏、锡、常地区的沉降速度也在加大，长江三角洲地区的地面沉降在区域上有连成一片的趋势。另外，全国31个省（区、市）都不同程度地存在着与饮用水水质有关的地方病。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">滑坡灾害中地下水的作用日益被人们所重视，以至于形成了“无水不成滑坡”的观点。边坡的充水张裂隙将承受裂隙水的静水压力作用；边坡地下水的渗透将对边坡体产生动水压力，地下水对边坡岩土体产生软化或泥化作用，地下水的溶蚀和潜蚀对边坡产生直接的破坏作用。根据边坡的地貌形态，结合边坡地下水环境响应的特征，如边坡泉的分布、动态、水质和混浊度的变化，可预测和评价边坡的稳定性。</span></p><p style="text-align: justify;"><span style="color: #FF0000;"><strong style="line-height: 1.5em;">（三）几个关键问题</strong></span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">1、地震预测预报</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">由于地震成因的复杂性和发震的突发性，以及人们现时的科学水平所限，地震预报是一个世界性的难题。在震前的一段时间内，震区附近总会出现一些异常变化。如地下水的变化，突然升、降或变味、发浑、发响、冒泡。如何利用地下水动态的宏观和微观观测数据，综合分析，对发震的时间、地点和震级进行预报值得深入研究。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">2、非线性水流方程与沉降方程耦合</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">地下水开采引起土层的固结，宏观上表现为地面沉降。在土层固结过程中，土层被压密，孔隙度和孔隙比减小，因而改变土的渗透系数K和贮水率μs。如果把K和μs视为常量，则属于线性固结问题；若随着固结过程参数K和μs是变量，则是非线性固结问题。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">地下水开采条件下，地下水的流动问题与固结问题是同时发生和发展的。固结问题蕴涵在地下水流动问题之中，同时固结问题又通过对地下水流动介质性质的改变而直接影响地下水流动问题。如何建立非线性水流方程与沉降方程的耦合模型，是正确计算、预测和分析地下水在开采条件下地面沉降量的前提。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">3、多尺度非均质中地下水流与溶质运移</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">地质体是非均质的，各个测度上都存在非均质性。为了掌握地下水及污染物在不同尺度多孔介质中的运移规律，进一步对地下水的污染进行预测、控制和治理，必需要回答一些重要的问题，如用不同方法、不同尺度得到的渗透系数K如何联系起来；弥散系数尺度效应的结构，哪种测量尺度对被污染的地下水流运动和扩散起控制作用，在解决实际问题时，弥散参数如何选取，如何采用随机的方法来模拟、预测地下水污染运移，如何建立溶质运移理论与模拟试验之间的关系等。因此，从理论、试验和实际应用方面研究多尺度非均质介质中渗流和弥散对于解决上述问题有很重要的意义。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">4、非饱和带水动力学</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">降雨入渗在边坡体内形成的所谓暂态饱和区及暂态水压力对边坡稳定性有极大得影响。然而，至今尚无可应用的定量研究成果，已有的工程设计均采用对暂态饱和区及暂态水压力进行假定的方法。近年来，随着非饱和水动力学的发展，为如何恰当估计基质吸力对边坡稳定性的影响提供了新的理论基础及相应的分析方法。非饱和水动力学理论认为，降雨影响边坡稳定性、诱发滑坡的主要作用机理是：降雨入渗使得边坡非饱和带土体的基质吸力降低、产生暂态饱和区，而基质吸力的降低，使得边坡非饱和带土体的抗剪强度下降，进而导致边坡稳定性降低，甚至滑坡。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">5、裂隙介质水动力学</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">研究裂隙介质地下水运动规律，究其实质最终可归结为3个方面的问题：其一是介质；其二是水；其三是水与介质的相互关系。介质研究的核心是其透水性，包括岩体的透水性、结构面的透水性及其空间分布规律。水主要是指重力水，研究的核心是地下水的质、量和力三方面问题。质指的是地下水的水质，量指的是地下水的流量（数量），力指的是地下水的静水压力和渗透力。由于裂隙介质本身的复杂性，使的研究裂隙介质中地下水运动规律难度很大。</span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">近年来，裂隙岩体地下水受到了严重污染的威胁，诸如核废料的地质贮存、垃圾填埋造成污水下渗、海水入侵、输油管道老化而引起的渗漏等。这些与人类生活密切相关的环境问题迫切要求我们对裂隙介质地下水中污染物运移的机制进行研究，以便对其进行预报和控制。此外，随着地热开发、稠油热采以及一些大型水利水电工程中的地下水热量运移问题的提出，也为裂隙介质地下水传热研究开辟了一些具有重要应用背景的新研究领域。</span></p><p style="text-align: justify;"><span style="color: #00B050;"><strong style="line-height: 1.5em;">6、水与介质相互作用</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">水与介质的相互关系重点是研究介质中渗流场与化学场、温度场和应力场之间的耦合作用，几乎所有的城市地质灾害问题都涉及到水与介质相互作用。渗流场、化学场、温度场和应力场之间的耦合机理，是分析解决相关城市地质灾害问题的理论基础，如地下水中的溶质和贮存在地下的热量，随着水流的运动在空间上发生迁移，就是一个典型的渗流场、化学场、温度场三场耦合问题。水库诱发地震、地面沉降都是典型的渗流场和应力场耦合问题。对于这些问题如何建立对应的数学模型，如何计算分析，有待进一步研究。</span></p><p style="text-align: justify;"><span style="color: #FF0000;"><strong style="line-height: 1.5em;">（四）结语</strong></span></p><p style="text-align: justify;"><span style="line-height: 1.5em;">地质灾害大多是内外动力地质作用产生的，不依人的意志转移的，但许多城市地质灾害与人类活动密切相关。地下水作为地质灾害萌发、发展和形成过程的重要因素，加强对其动态、运动规律、物理化学特征等的研究，有助于加深对地质灾害，尤其是人类活动引起的城市地质灾害形成机理的认识，从而实现对地质灾害的预测预报和控制，做到人与自然的和谐共处。</span></p>