2012数学建模B题解答.docVIP

一 问题重述 辐射井是由一大口径的竖井和自竖井内周围含水层任意方向、高程打进一层数条水平辐射管组成，地下水沿水平辐射管汇集到竖井中。辐射井出水量较普通管井、筒井及大口井的高出数至十倍以上，具有较强的取水能力。由于其辐射管水平伸入含水层达数十至上百米以上，其水位下降范围大，降速快，具有良好的排水效益，被广泛应用于地基施工降水、尾矿坝排水等工程领域。经过多年的应用和发展，辐射井在成井技术方面已逐渐成熟。然而，长期以来，对辐射井井流计算理论却滞后于辐射井成井技术的发展，使该项技术的进一步应用受到极大的阻碍。由于辐射井结构的特殊性，地下水运动系三维流，因而其渗流计算不能随意套用其它水井的计算公式。在实际工程中，采用辐射井作为基础工程降水施工、农田地下水人工调控，均需考虑辐射井水位降深计算问题。然而，目前多按抽水试验建立一些经验公式来近似计算辐射井的出水量，难以求解辐射井的降落曲线。这一重大理论和实际问题，因其研究不够，长期未获得满意解决，至今尚未建立系统完整的辐射井取水和降水的水量计算理论。 在下列假定条件下：潜水含水层均质各向同性，隔水底板水平，在平面上无限分布；不考虑水和介质骨架的压缩性；潜水完整井，无越流补给也无入渗或蒸发。要求依据黄土地区测得的实验数据（见后图1、2、3及表1）： 设计构造黄土地区辐射井的地下水降落曲线的数学模型； 分析黄土地区地下水位降深与出水量的关系，建立造黄土地区辐射井的水量计算公式。 二 问题分析 2.辐射井的地下水降落曲线即辐射流的流动特征?? 地下水在天然流动过程中，以层流和缓变流的形式运动，无争议。但是地下水向集水井辐射流动过程中，由于渗透流场内的水力坡度和渗透速度增大，地下水辐射流是否还保持层流和缓变流的运动形式，迄今仍有较大的分歧，根据孔隙含水层和裂隙含水层中大量抽水试验资料的分析研究，初步证实，在抽水井作用酌汇点辐射流场内，除井壁附近地下水可能有层流转化为混合流或紊流以及出现明显的空间留动特征外，在广大降落漏斗的范围内仍基本保持层流宇缓变流的运动形式。?? 3.辐射井出水量的计算?? 依据潜水井出水量方程推导涌水量计算公式，并且参考实际水利工作中总结得到的经验公式加以修正，最终给出出水量计算。 条件假设与符号 3.1?条件的假设 （1）潜水含水层均质，隔水底板水平，在平面上无限不考虑水和介质骨架的压缩性； （2）潜水完整井，无越流补给，也无入渗或蒸发； （3）垂直方向的渗透性远大于水平方向的渗透性； （4）地下水处于稳定流动状态； （5）地下水成层流运动，遵循达尔西直线渗透定律； （6）地下水为缓变流，可将空间流简化为平面流； （7）假定静水位是水平的，降落漏斗的供水边界是圆形的； .2符号的 符号 符号说明 距集水井中心的水平距离为x处得横剖面水位平均高度（米） 集水井中水面的高度（米） 辐射管端点的水位高度（米） 反映降落曲线弯曲程度的经验数值 辐射管端点距井中心的水平距离（米） 集水井的半径（米） 距离水井较远的一个观测孔的水平距（米） 距集水井较近的一个观测孔的水平距（米） 距离为的观测孔中测得的水位高度（米） 距离为的观测孔中测得的水位高度（米） 井壁处的水位高度（米） 辐射井的出水量（米小时） n 辐射孔（管）的根数 k 黄土含水层的渗透系数（米小时） 局部阻抗系数 相邻两根辐射孔（管）之间的夹角（度） d 辐射孔（管）的直径（米） 辐射孔（管）中心的高度（米） ?????? ????????????????? 四??辐射井 4.1 对辐射井工作状况纵剖面图和测得的实验数据的分析 据图形可以很明显的看出：辐射井的工作状况与普通筒井有很大不同。就地下水降落曲线看，一般筒井影响范围内,地下水降落曲线呈上凸的抛物线，近井处水力坡度很陡，远井逐渐变平缓，在井壁处往往有明显的水跃发生。而辐射井在水平辐射管延伸范围内，降落曲线与筒井刚好相反，近井处水力坡度平缓，远处陡峭，呈下凹曲线，井壁处几乎不发生水跃。在辐射管的端点处，水力坡度陡增，并出现凸凹拐点。在辐射管延伸范围以外，降落曲线改变成上凸的抛物线，水力坡度由陡变缓，与普通井一致。 由于黄土垂直方向的渗透系数比水平方向大若干倍，加上辐射管的水平位置的影响，黄土地含水层的地下水基本是沿着垂直方向进入辐射管。当地下水位由于抽水下降时，这一现象在水平集水管延伸范围内尤为明显。此时，大部分水量由辐射管截取后，通过管道汇入集水井，而通过集水井井壁和辐射管端点外侧流进的水量，则占很少比例。为方便计算，可以只考虑沿垂直方向流入辐射管的水量。辐射管汇集水量的大小与降落曲线高度, 基本上成正比例的关系。管中水流, 由于有压管流, 其水头沿程变化相对很小, 可视为相当于井水位的一个固定值。