基于达西定律对水量计算问题的研究

2010高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：B我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：2012年8月24日赛区评阅编号（由赛区组委会评阅前进行编号）：2010高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于达西定律对水量计算问题的研究1摘要本文通过设计构造辐射井的地下水降落曲线的数学公式，来建立辐射井水量的计算模型。首先我们就潜水井的出水量，降落漏斗的最大半径，井位降深，流线等一些专有名词进行了解释，给出了相关公式的推导过程，以及该公式的使用范围。在合理的假设下，建立最合理的数学模型，经过推导得到的公式，用于求解：出水量，降深，降落曲线，漏斗半径等。然后我们把模型计算的结果进行对比及检验，并证实了我们所建立的模型合理性、有效性，精准性。在模型评价时，我们正视该模型的缺点，并对其进行改进。针对黄土地区辐射井的地下水降落曲线问题，本文基于对地下水降落曲线高度xT与距离x之间的距离以及集水井中水面的高度pT的变化的分析，建立地下水降落的)(00)(rRppxeTTTT模型。同时，在达西渗透流速定律的基础上，建立集水井的水位高度PT与时间t的关系模型*pkITHutHte，进而根据影响半径21211222112lg2lg2lgSSSSHrSHSrSHSR对其水位降深进行分析。并将观测井与实际值进行分析验证，证明模型的可行性。针对黄土地区辐射井水量计算的问题，首先分析辐射管在含水层对地下水运动状态的影响情况，得出辐射井水量只需要考虑垂直方向上的水流量。而后利用问题一中所求的降落曲线，采用定积分的方法，求出沿管全程的平均高度T,及与之相关的水平距R、剖面矩形宽b、局部抗阻系数，再结合经验公式单位长度流量xpxxTTqk，求出辐射井的水量计算公式：PnTTrRnkQ0。最后用4.19.21和4.21.17这两组数据进行残差检验，误差分别为1.487%和4.52%，因此该数学模型的可行性比较高，可以用来计算黄土地区的水量。本文的特色在于巧妙的运用了达西定律与已知条件相结合，巧妙地分析了集水井的的水位高度与时间的变化关系。同时分别从距井的距离和集水井的水面高度与时间的关系来对辐射井的地下水降落曲线进行分析，考虑较为全面。关键字：水位高度达西定律水位降深渗透系数局部抗阻系数2问题的叙述2.1问题的重述辐射井是由大口径钢筋混凝土管竖井和竖井内沿含水层水平方向布设的多根辐射管组成，水平管呈辐射状分布，垂直竖井主要为集水井,汇集由辐射管流入的地下水,同时作为辐射孔（管）施工场地和提水机具安装场所。一般情况竖井的井壁是封死的。辐射管主要通过水平滤水管进水。辐射井以其出水量大、寿命长、管理费用低、维修方便、便于集中管理等优点，在水力资源开发利用上被广泛应用。从20世纪60年代以来，辐射井技术已在我国推广应用。如辐射井在华北深基坑工程降水中，取得了较好的效果。随着北京东直门地铁站采用辐射井降水技术取得成功，目前辐射井技术在地铁施工又发挥着重要作用。辐射井的出水量，是设计和布置辐射井工作中所需要解决的问题之一。同管井和大口井的计算理论相比,在题中的假设条件下，结合图片和数据：1.设计构造黄土地区辐射井的地下水降落曲线的数学模型；2.分析黄土地区地下水位降深与出水量的关系，建立造黄土地区辐射井的水量计算公式。2.2问题的分析本题研究的是辐射井地下水降落曲线、出水量的数学模型。首先要明确影响辐射井水量计算的可能存在因素：辐射管数量、长度和分布、水位下降时间、渗透系数、孔隙比、井水的种类、地下水的流态等等，其中辐射管数量、长度和分布、渗透系数、孔隙比、井水种类等因素题中已经给出了确切的数据和限定，使我们在构造辐射井地下水降落曲线数学公式和计算水量的过程中可以直接应用。针对问题二，我们对已知数据进行分析，可以看出水降落曲线高度xT与距离x之间近似呈自然对数的函数关系。同时在假设条件成立的情况下，利用达西渗透公式对集水井的水位高度pT随时间t的变化进行分析，并需分析其影响半径对其出水量的影响。针对问题二，由于辐射井的水量主要由各个辐射管的水量汇集而来，而辐射管的水量则是由无数个小的单位长度的水流量累加而成。因此通过根据已有的经验公式和数据，建立单管的水量计算公式，进而求出8个管的辐射井水量计算公式。2.3名词解释层流：laminarflow水流流速彼此不相混杂、运动迹线呈近似平行的流动。紊流：turbulentflow水流流束相互混杂、运动迹线呈不规则的流动。孔隙比：指材料中空隙体积与材料中颗粒体积之比，是反应材料密实程度的重要物理性质指标。一般用e表示e越大材料越疏松。水头损失：液体和固体相互作用而产生的结果。水力坡度：河流水面单位距离的落差。在水力学中，表明了实际液体沿元流单位流程上的水头损失，即总水头线坡度。它是单位重量液体沿流程单位长度上的机械能损失。渗透系数：在各项同质介质中，单位水力梯度下的单位流量，是岩石透水性强弱的数量指标。3模型假设辐射井的出水量，是设计和布置辐射井工作中所需要解决的问题之一。为了更方便更现实地对辐射井理论进行研究，我们做出下列条件假设：（1）假设潜水含水层均质，隔水底板水平，在平面上无限分布，不考虑水和介质骨架的压缩性。（2）假设所研究辐射井为潜水完整井，无越流补给，也无入渗或蒸发。（3）假设竖井的直径3.5米。水平辐射管的长度为120米，8根辐射管均匀的分布在距隔水底板1.2米平面上，辐射管的直径为0.12米。（4）假设黄土含水层的渗透系数k在0.0554~0.1607(米/小时),孔隙比（含水层中空隙的体积与固体颗粒体积之比）为0.75%。（5）假设由实验已知，辐射井的地下水降落曲线在水平集水管（辐射管）延伸范围内，呈凹形的抛物线，近井处水力坡度平缓，远处陡峭。在辐射管的端点，水力坡度陡峭曲线出现凹凸拐点。在辐射管延伸范围以外，降落曲线改变成凸形的抛物线，水力坡度由陡变缓。4符号说明pT集水井水面高度xT降落曲线高度0R水平辐射管的长度（米）。iR水平上距离辐射管的距iT水平上距离为iR的水位高度r大口集水井的半径（米）；rR辐射管的长度（米）；H隔水底板与地面的高度差h△地下水的水头差（米）u——地下水的实际渗透速度；t——水位下降所用时间。u——地下水的平均渗透速度（米/小时）e——孔隙比I——水力坡度R影响半径wh——地下水某时的最低水位（米）；iS第i号观测孔的水位降深ir观测孔距主井的距离Q为单井出水量5模型的建立与模型求解5.1地下水降落曲线的模型5.1.1对降落曲线的分析结合题中相关数据，得出地下水降落曲线高度xT与距离iR之间近似呈自然对数的函数关系，并且随集水井中水面的高度pT的变化而变化，因此，要构建地下水降落的模型，就需要分别从距离集水井的距离x和水位PT两个方面进行分析。流线方向辐射管辐射井的地下水的降落曲线筒井的地下水降落曲线集水竖井t1t3t2rR0地面时间OXZ不透水层∆H图1辐射井工作状况纵剖面图5.1.2地下水随X变化而变化的模型通过查找资料获得地下水在水平距离X方向的变化有如下的经验关系：)(00)(rRppxeTTTT为表示降落曲线弯曲程度的一个经验数值ppTTTTRR1212ln1将0Rx代入（1）式，得到：0TTx，将rx代入（1）式，得到：)(00)(rRpprxeTTTTT现以4月19日21时为例计算地下水的降落结果如下表：井的距离r=1.7515.0030.005060.0080.00110.00121.75)(0xR2.64442.46262.02181.58111.36070.92000.47930.00)(0xRe14.075011.7357.55194.86033.89892.50931.61491．00xT(计算）77.8577.8978.0578.3678.7279.5380.4581.23表1计算值与真实值之间的对比用matlab对计算值进行拟合（程序见附录1）图2计算值的拟合通过对曲线的拟合可以看出，模型比较准确，上图呈现了下水降落曲线高度xT与距离x之间的函数关系。因此可以确定该模型的准确性。用matlab对主要观测孔随时间的变化关系进行编程（附录2），从而得到下图：xT(实测）———78.43——80.4181,160102030405060708090100727476788082848688抽水时间第一次降深集水井观测孔N1观测孔N2观测孔N3'观测孔N5图3主要观测孔随时间的变化关系通过上式，可以看出集水井的水位高度与时间有着明显的关系，从而也会引起降落曲线与距集水井X之间的函数关系发生变化，因而下面就对集水井的水面高度与时间的变化关系进行分析。（1）00~R段（辐射管延伸范围内）：辐射井的水位降落之前pTH（1）水位下降高度Tut（2）达西渗透流速公式/TutuueukI通过查资料得知，地下水的水力坡度与气流动状态相关，即当地下水为层流运动时，'II；当地下水为紊流运动时，'II；当地下水为混合流运动时，'mII。*pkITHutHte然而，水体在流动的过程中，并不是一成不变的。在运动过程中单位质量液体的机械能必将发生损失，因而引起水头损失。产生水头损失的原因有内因和外因两种，外界对水流的阻力是产生水头损失的主要外因，液体的粘滞性是产生水头损失的主要内因，也是根本原因。而水力坡度/IhR△（8）沿程水头损失：22fRvhdg局部水头损失：22mvhg总水头损失：222()222wfmRvvRvhhhdggdg用压强表达：2222fmRpvpdpvp总水头损失：222()222wfmfRpvpvRpvppppdd为沿程阻力系数；为局部阻力系数；d管道直径；v辐射管中水的平均流速g重力加速度通过上式可以计算出水头损失值：22Rvhdg因而求出相应的水力坡度为：22RvIdRg2*2pkvTHutHtgdRe综上，可以得出在00~R段集水井的水位高度Tp与时间t的关系。然而通过分析可以看出：经过一段时间之后，水位下降之后还会升高，这是由水体间的压力所致。（2）潜水井的Dumpuit模型（0~RR段）对上述潜水井应用Dupuit假设，认为径向井的潜水流是近似水平的，因而等水头面仍是共轴的圆柱面，井和过水断面一致，这一假设，在距抽水井r1.5H的区域是足够准确的。同时认为，通过不同过水断面的流量处处相等，并等于水的流量。这时漏斗区潜水流的水头分布满足0HHhhxxyy式。以潜水含水层的底板作基准面，h=H，并用柱坐标形式表示，则方程简化为20hhxr。其边界条件和承压水井相似，为h=H，当wrr时，h=H0,当r=R时，对上式积分，得2hrCr因各断面流量相等，根据通过任意断面的流量220hdh