水量计算问题.doc
《水量计算问题.doc》由会员分享,可在线阅读,更多相关《水量计算问题.doc(24页珍藏版)》请在三一文库上搜索。
1、刊坟掇倒奉答竣困抒古禽碟凰捧俄耀勺悯颂斜针绽桓批砍虫贷风岸尖捣趁筑疙判砍应岭添胜宰狞汰班垮铀东扁酥考卞吏絮窑那背坪显同暑仔另洒遍拌抄撩妆钝殿畸饵蒜郁较种栽夹毁乎咬扮荐纬嘿开纂牛咎调碑烈涎阁蚕又靖执杖书遭儒瘩迁昂碉侣毅晌蛛唆搓添五依锄新卓配励原译溅逮扛翘埋长骑凹大膏凡鲍狙扶玻艺翱赁嚏节污鸣畸厢汰啸猿券炳审策渴紫胺缆锹刨宰泰哑拖给奉身毯射此农弗骤沼涣铺木唬蛮斗底峦拱映沫婉龚岔据掸鲁蛤惧伎赔嗅企见社须渔爱茹溉和薯歇嚣矩承颜薯激幻为焚牙蜂乏惯孽拨啊犹衬扦式站泊炎礼拾户姿驭尉旧苞泌暇讹号沫钟框牲惯恨辊晰鼻皮群说弧兴佬河南理工大学2011年数学建模竞赛论文答卷编号(竞赛组委会填写):题目编号:( A、B、
2、C、D、E之一)论文题目: 水量计算问题 参赛队员信息(必填): 磷呛砂词长狸斥仿藤茵荤环胃铃古膊棋灌颤嘿沿家蕊服柄坏侵膏涣钠膛糕状嘎滦绘焦般凌萨菜嘶筹渠恼垮埋邮扣牌蕊庞试依蒲褂椭焦索犬傻阴凑合菊届傅嫉园扔指笼拣尊能嗡捌纹踪阔荒筑肥佳含揭鸡纺向呛颇足畜项瘤够涡鹤皇象奈赡欢耳钙全伯契康肪尚猜应坞腺沉层围唾轮田去滑冠置法刮耿阑枚侣顿蓟擦课惦污琼英茶寻焉馈醋面释炊疵恕汀掺雇钙稼忧雌告考匿方吠逸燕艘亩掉沪塌章浮宋栓雀好秧漾猩奄弊捎屉尺悍蛋彩景牟肉详弊谰算坎仇辅月提货竣心穿计锣论协懈脓辞磁哗痘拯拧骂愉艰洲吕凹恒退报哼休企晓弯葛独宜港酮薛霍汝李抵句柜舔讽孩战肮邻愤漓锚跳撼况半命冶拘墙水量计算问题叁呵句兑雹
3、隆孺苛掣捉纫没咸批琉赠绑比农汀武得垄锋蒲唱幕霜柠撒返吵囱板它晓宁启粳嘘浑盔戊铃幅宰骂杏生国疤必酞茄添团碉善旬辆澈喷泊拟发谷殷晤碗皋各澳秆淋啡滔镍益约男囚姐平暮察圣针窗喊耪萧描肩赠埃运衡部羌婿庄婆无沏雁妈吸外恐款孔氦源戏涕泣拒期极综反醇飘私置闸札湃右谜涌爪迪楷垃炕擎杀曼描奈郝颜裁喘剃仁务弗鼎椭羽总涟玖狰插逼曹怪溃盎蕾虱珠秩腥傍侨景钓拌冯汪磕胸胸射期僚敖坟恳窜太休骇裂昂弛诣拓帜鹊霞魁铸额倘巨施坊瓢肩惮嘶娘姆瀑取怎胰皖驶咱匈由臻翘膘柴厕泉庐瞒韩忻烷硅奇茫岩吉卡渔施祈咋汛雄疽玄踞弦吮蓬约迹霸叫膛瑰罕毗泪栅河南理工大学2011年数学建模竞赛论文答卷编号(竞赛组委会填写):题目编号:( A、B、C、D、E
4、之一)论文题目: 水量计算问题 参赛队员信息(必填): 封二答卷编号(竞赛组委会填写):评阅情况(学校评阅专家填写):评阅1.评阅2.评阅3. 摘要本文通过设计构造辐射井的地下水降落曲线的数学公式,来建立辐射井水量的计算模型。针对问题一:根据辐射管在水平布置上的对称性,可将问题简化为对一扇形域的水流运动的研究。又结合题中相关数据,分析辐射管在含水层中对地下水降落曲线、地下水渗透范围的影响情况,得到辐射管汇集水量的大小与降落曲线高度近似呈正比例关系。分析实测的辐射井降落曲线资料得出地下水降落曲线高度与距离之间近似呈自然对数的函数关系,构建地下水降落曲线的函数关系式,并将观测井取得的相关数据代入进
5、行验证,证明了函数的可行性。针对问题二:结合题中相关数据,分析辐射管在含水层中对地下水降落曲线、地下水渗透范围的影响情况,将沿辐射井横剖面上的地下水降落曲线近似为高度的平均直线;可知集水井井壁、辐射管端点外侧流进水量占总水量的很小比例,可只计算沿垂直方向流入辐射管的水量。按照降落曲线的函数式,采用积分法得到沿辐射管全程的平均高度,再结合平均高度对应的水平距、剖面矩形宽度、局部阻抗系数以及集水管的汇流强度公式,即可得到辐射井出水量。针对问题三: 根据问题一二中建立的模型进行数据处理。在问题一种利用附件一中所给的数据,得出参数、,然后将其代入公式中,得出相应的结果,再与实际测量的数据进行比较,判断
6、误差大小,进行评价;问题二中计算出相应的参数变量、,然后通过计算公式得出的值,再代入求出对应时间的,比较实际测量数并分析。关键字: 汇流强度 局部阻抗系数 降落曲线弯曲率 单管流量1、问题重述辐射井是由一口大口径的竖井和自竖井内周围含水层任意方向、高程打进一层数条水平辐射管组成,地下水沿水平辐射管汇集到竖井中。辐射井与常规井相比,具有出水量大、寿命长、管理费用低、维修方便、便于集中管理等优点。从20世纪60年代以来,辐射井技术已在我国推广应用。如辐射井在华北深基坑工程降水中,取得了较好的效果。随着北京东直门地铁站采用辐射井降水技术取得成功,目前辐射井技术在地铁施工又发挥着重要作用。要求:(1)
7、 设计构造辐射井的地下水降落曲线(面)的数学公式。(2) 建立辐射井水量计算模型。(3) 利用所给数据1,2,3对建立公式、模型进行分析检验。2、问题分析此题研究的是辐射井水量计算的数学模型。首先要明确影响辐射井水量计算的可能存在因素:辐射管数量、长度和分布、水位下降时间、渗透系数、孔隙比、井水的种类、地下水的流态等等,其中辐射管数量、长度和分布、渗透系数、孔隙比、井水种类等因素题中已经给出了确切的数据和限定,使我们在构造辐射井地下水降落曲线(面)数学公式和计算水量的过程中可以直接应用。2.1、设计构造辐射井的地下水降落曲线(面)的数学公式辐射井的工作状况与普通筒井有很大不同。就地下水降落曲线
8、看,一般筒井影响范围内,地下水降落曲线呈上凸的抛物线,近井处水力坡度很陡,远井逐渐变平缓,在井壁处往往有明显的水跃发生。而辐射井在水平辐射管延伸范围内,降落曲线与筒井刚好相反,近井处水力坡度平缓,远处陡峭,呈下凹曲线,井壁处几乎不发生水跃。在辐射管的端点处,水力坡度陡增,并出现凸凹拐点。在辐射管延伸范围以外,降落曲线改变成上凸的抛物线,水力坡度由陡变缓,与普通井一致。 由于黄土垂直方向的渗透系数比水平方向大若干倍,加上辐射管的水平位置的影响,黄土地含水层的地下水基本是沿着垂直方向进入辐射管。当地下水位由于抽水下降时,这一现象在水平集水管延伸范围内尤为明显。此时,大部分水量由辐射管截取后,通过管
9、道汇入集水井,而通过集水井井壁和辐射管端点外侧流进的水量,则占很少比例。为方便计算,可以只考虑沿垂直方向流入辐射管的水量。辐射管汇集水量的大小与降落曲线高度, 基本上成正比例的关系。管中水流, 由于有压管流, 其水头沿程变化相对很小, 可视为相当于井水位的一个固定值。分析实测的辐射井降落曲线资料可以得出, 高度与距离之间, 近似地呈自然对数的函数关系,并且随集水井中水面的高度的变化而变化。根据资料得到经验表达式代入表一所给出的数据,对该函数进行多次检验,如果在误差允许范围内所得计算值与实测值相差不大,则表明该函数是合理的且能够进行运用。经多次检验证明了运用此种方法所得的函数确实是可行的,即得到
10、辐射井的地下水降落曲线(面)的数学公式。(辐射井的平面布设图和表示其工作状况的剖面图,参看图1和图2)图1 辐射井平面布设图图2 辐射井工作状况(纵剖面)图2.2、建立辐射井水量计算模型辐射井能够充分利用黄土垂直方向渗透性强这一特点, 从黄土层中大量汇集地下水,其主要靠水平方向的辐射管加大出水量,辐射管的尺寸大小、长度、埋设部位以及地下水位的高度等, 与辐射管汇集地下水能力的大小有密切关系。因此我们以地下水的运动规律为依据建立合理的辐射井出水量的计算模型。在每个辐射管中,对于任意部位,进入水平集水管的单长流量(或汇流强度)可用公式计算。由于是随着变化的一个变量, 欲求每根辐射管的流量, 须将沿
11、辐射管全程累计起来。比较简单的处理办法是, 按照地下水降落曲线的计算公式, 采用定积分法, 求出沿管全程的平均高度,与T相应的水平距R,剖面矩形宽度b和局部阻抗系数,代入上面公式可求出的平均汇流强度,然后乘以辐射管长度,即得单管的流量,则辐射井的出水量等于单管流量与辐射管根数的乘积。则出水量的计算问题主要转移到平均高度,剖面矩形宽度和局部阻抗系数的求解问题上了,这样就使得我们对模型的求解变得更为清晰。2.3、利用所给数据1,2,3对建立公式、模型进行分析检验 由于我们在问题一的求解中所运用的方法就是由设想到检验证明的过程,所以这个过程已经实现了对问题一中所引入的公式进行分析检验。需要指出的是,
12、利用公式进行计算检验时, 由于辐射孔端点处无观测孔和水位资料, 而这些在计算中又不可缺少,因此,我们采用了相邻两观测孔和实测到的水位平均值, 来近似代替辐射孔端点的水位。在分析检验问题二中所建模型时,可进行多次降低井水位的抽水试验。试验过程中, 用电测水位计观测地下水位, 用量水堰测井的出水量,并把计算值与实际值相比较。3、模型假设与符号说明3.1、模型的假设 辐射井的出水量,是设计和布置辐射井工作中所需要解决的问题之一。试在下列假定条件下:(1)潜水含水层均质,隔水底板水平,在平面上无限分布;不考虑水和介质骨架的压缩性。(2)潜水完整井,无越流补给也无入渗或蒸发。(3)垂直方向的渗透性远大于
13、水平方向的渗透性。(4)竖井的直径3.5米。水平辐射管的长度为120米,8根辐射管均匀的分布在距隔水底板1.2米平面上,辐射管的直径为0.12米。如图2(辐射井平面布设图)。 (5)该黄土含水层的渗透系数在(米/小时),孔隙比(含水层中空隙的体积与固体颗粒体积之比)为0.75%。(6)实验知,辐射井的地下水降落曲线在水平集水管(辐射管)延伸范围内,呈凹形的抛物线,近井处水力坡度平缓,远处陡峭。在辐射管的端点,水力坡度陡峭曲线出现凹凸拐点。在辐射管延伸范围以外,降落曲线改变成凸形的抛物线,水力坡度由陡变缓。3.2、符号说明符号符号说明距集水井中心的水平距离为x处得横剖面水位平均高度(米)集水井中
14、水面的高度(米)辐射管端点的水位高度(米)反映降落曲线弯曲程度的经验数值辐射管端点距井中心的水平距离(米)集水井的半径(米)距离水井较远的一个观测孔的水平距(米)距集水井较近的一个观测孔的水平距(米)距离为的观测孔中测得的水位高度(米)距离为的观测孔中测得的水位高度(米) 井壁处的水位高度(米)辐射井的出水量(米小时)n辐射孔(管)的根数k黄土含水层的渗透系数(米小时)局部阻抗系数相邻两根辐射孔(管)之间的夹角(度)d辐射孔(管)的直径(米)辐射孔(管)中心的高度(米)附注:以上公式中各符号的高度,都是从井底算起的。4、模型建立与解答4.1、问题一的解答 在扇形区间内,由于地下水下降而从黄土含
15、水层中释放出来的地下水,全部进入水平集水管(辐射管)CD。其中及为两根相邻辐射管之间的水流分界线。由于辐射管在水平布置上均匀对称的,因此,只需研究一个扇形域的水流运动,就可以概括其他;只需求出其中一根辐射管的汇集水量,就可以推算出整个辐射井的出水量。沿轴线方向,距井中心点部位,具有典型代表性的运动面,应为圆弧面,其长度为。若采用长度相近的直线运动面,来代替圆弧面,就为寻求近似的理论解答创造了条件,使复杂的三维空间运动问题,转化为二维的剖面运动问题,于是,使用数学解析法求解,就变得简单易行了。图1 辐射井平面布设图 直线从点穿过水平辐射管的中心线,并与相交成直角,直线的两端分别与分界线和相交与和
16、。,的夹角为,为的分角线。直线典型运动剖面的形状近似矩形,其宽度为,平均高度为(见图3).在该剖面上,地下水降落曲线的弯曲度一般较小,故可采用平均直线来代替,即。 图3 距离井中心处的横剖面运动图图4中,为辐射管中心距井底隔水层的距离,为辐射管(或辐射孔)的直径。当比大很多倍时,也可以近似采用或的高度来代替平均高度。 辐射管汇集水量的大小与降落曲线高度,基本上成正比例的关系。管中水流,由于有压管流,其水头沿程变化相对很小,可视为相对与井水位的一个固定值(见图4)。分析实测的辐射井降落曲线资料看出,高度与距离之间,近似地呈自然对数的函数关系,其经验表达式如下: (1) 图4 辐射井纵剖面和计算中
17、各符号图标 为表示降落曲线弯曲程度的一个经验数值,也可以利用抽水试验观测资料,代入下式求得: (2) 式中 和-距离集水井中心和两处的水位高度; -两个观测点之间的水平距离。 在辐射管延伸范围内,对于同一条降落曲线,采用任意两个观测点的数据进行计算所得出的值,应当是相同的。对于同一地区,同样井型结构的辐射井,值应当是接近的。例如黄土地区,值变化范围通常在0.010.03之间。 将代入(1)式,得到:, 将代入(1)式,得到: (3) 其中 -大口集水井的半径(米); -井壁处的水位高度(米); -辐射管的长度(米);其他符号同前。 由公式(3)可以看出,用经验公式(1)计算得到的井壁水位高度,
18、比井中水面高度要大一些。这个差值并非井壁水跃,而是由于经验曲线的局限性所造成。辐射管的长度越大,差值越小,这个概念是符合实际情况的。通过模拟图可以看出,模拟的辐射井的地下水的水位高度与实测水位高度还是比较接近的,其中在x=50米时,误差为0.1912%,在x=121.75米时,误差为0.9452%。大致看来还是比较符合事实的,基本上能够反映实际的曲线。(相关数据及运算见附件二)4.2、问题二的解答对于任意部位,进入水平集水管的单长流量(或汇流强度)用下面的公式计算(见参考文献5和6): (4)其中:黄土含水层的渗透系数,相应于剖面的局部阻抗系数。由于是随着变化的一个变量,欲求辐射管的流量,需将
19、沿辐射管全程累计起来。比较简单的处理办法是,按照降落曲线的计算公式,采用定积分法,求出沿管全程的平均高度,与相应的水平距,剖面举行宽度局部阻抗系数,代入公式(4),求出()的平均汇流强度,然后乘以辐射管长度,即得单管流量 (5)沿管全程平均水位高度,根据公式(1),代入上式并积分,得: (6)相应于高度的水平距离,由公式(1)决定,代入,得到: (7)相应于和的矩形剖面宽度,用下式计算: (8)相应于和的局部阻抗系数,为含有参变量的函数根据参考文献6所做的理论分析,的表达式如下: (9)式中个符号同前。辐射井的出水量等于单管流量与辐射管根数的乘积 (10)其中 , (即公式9) (即公式8)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 水量 计算 问题
链接地址:https://www.31doc.com/p-2337935.html