试算法水文
Ⅰ 地下水资源模型计算
数值模型模拟计算方法适用于非均质性、各向异性的复杂地下水系统，包括存在越流和具有不规则形状各类边界条件等情况。但是该方法对资料的要求比较严格，要求研究程度较高和资料较丰富。应用数值模型方法的一般程式为：①气象、水文、水文地质资料分析→②水文地质概念模型概化→数学模型建立（水动力方程和定解条件）→③选择计算程序→④模型设计→⑤模型识别和检验→⑥计算模拟。
一、水文地质概念模型建立
在对黑河流域地下水系统做了全面、深入分析的基础上，根据研究目的，对地下水系统的组成要素和相互关系作出合理的简化和假设，并且用文字、框图、平面图、剖面图等形式把系统再现出来，即为地下水系统概念模型。
（一）地下水系统空间结构概化与边界确定
1.图式表示地下水系统空间结构
根据黑河流域水文地质图和水文地质剖面图，梳理和划分主要含水层、隔水层与弱透水层，阐明它们的产状、分布范围和厚度等，确定透水、阻水等断层属性。分析地下水系统的各类等值线图，包括第四系基底埋深等值线图、地下水水头等值线图、含水层顶底板高程等值线图、含水层和隔水层的厚度等值线图等。
2.确定地下水系统边界
地下水系统的边界，包括自然边界（固定边界）和水力边界（可移动边界）。自然边界包括不透水岩层、不透水断层或断裂带、较大的地表水体等；水力边界包括地下水分水岭和地下水流线等。
数值模型模拟研究，其对象的底界一般为不透水岩层。侧向边界可以是自然边界，也可以是水力边界或无穷远边界（边界水头或流量不受输入条件的影响）。模拟顶界，对于承压水系统而言，一般为不透水边界或越流边界，对于潜水系统一般采用大气边界（蒸发和入渗）。地下水系统内部边界包括零流量边界（不透水岩体）和流量边界（河流、湖泊或水库的渗流带）等。
3.水文地质参数
水文地质参数是数值模型模拟研究的灵魂，一般包括含水层组的渗透系数、导水系数、给水度、储水率、储水系数、孔隙度、垂向渗透系数和越流系数，以及包气带的降水入渗系数、河道渗漏系数、井灌回归系数、田间与渠道渗漏系数、潜水蒸发系数和陆面蒸发系数等。
确定降水入渗补给系数、灌溉渗漏系数、蒸发系数等方法，有水文分析法（降水量、河流径流量曲线、地下水水头动态曲线等）、直接试验法（地渗仪、张力计、同位素示踪等）、计算法（氯质量平衡法、非饱和模型法等）、经验公式法和ZFP零通量面实测法等。
（二）地下水流系统概化
对地下水流系统进行概化，包括确定地下水的基本流向、地下水补给要素组成、排泄模式、地下水与地表水之间转化关系、不同层位含水层之间水力关系等。主要依据有，地下水水头等值线图、水化学信息、同位素信息、地下水温度信息和水位动态曲线等。
根据地下水流状态及其特征，确定所研究的地下水流系统具体属性，例如稳定流或非稳定流，一维流、二维流、准三维流或三维流等。
（三）模型输入量计算
降水入渗、地表水入渗（河渠）、地下水侧向流入、灌溉入渗、蒸发蒸腾、泉水排泄、基流排泄、地下水侧向流出、开采等。
二、建立数学模型
根据建立的水文地质概念模型，选择适宜数学模型。一般由描述地下水运动规律的偏微分方程和反映地下水系统边界条件及初始条件的定解条件组成。
非均质承压水三维非稳定流偏微分方程为
西北内陆黑河流域水循环与地下水形成演化模式
非均质无压水三维非稳定流偏微分方程有下列几种情况：第一类边界条件（狄利克雷边界）为
西北内陆黑河流域水循环与地下水形成演化模式
第二类边界条件（纽曼边界）为
西北内陆黑河流域水循环与地下水形成演化模式
初始条件为
西北内陆黑河流域水循环与地下水形成演化模式
三、计算程序、模型设计与识别
（一）计算程序与模型设计
计算程序分为一维流、二维流、准三维流或三维流模型，以及对均质、非均质、各向同性或各向异性和对不同输入项的处理能力。目前，可供软件有MODFLOW、FEWFLOW、PM、GMS、GWVISTA、MODME、PM等，它们多为有限差分法和有限元法。模型设计，包括网格剖分（规则剖分或不规则剖分、三角剖分或矩形剖分）、选择时间步长（试算法）、设置模型边界、设置初始条件、数据输入（降水入渗速率、田间灌溉入渗速率、蒸发速率、水井位置及开采或回灌强度、地下水与地表水相互作用的时空分布、泉的时空分布、边界水位或边界流量、观测井位置及观测水位等）。
（二）模型识别与检验
1.模型识别
模型识别亦称反演问题，即利用实测地下水动态资料和抽水试验资料，反求水文地质参数或源汇项和定解条件的过程。模型识别是为了解决选用的偏微分方程是否合适问题，确定模型中的水文地质参数和源汇项及定解条件，从而建立一个能再现地下水系统实际功能（水头或浓度）的模拟模型。模型识别一般采用试估-校正法。就是选择一合适的时段，根据水文地质条件和经验数据估算一组水文地质参数输入模型，利用所选时段的输入输出数据，求解模型。然后以模型计算结果与实测结果比较，如果拟和结果不符合精度要求，适当调整参数，重复上述过程，直到符合精度要求为止。也可以采用试估-校正法与最优化方法相结合的方法。首先用试估-校正法粗调，然后用最优化方法细调，即用最优化方法求得一组最佳的参数值，使得计算水头值与观测值之间的差值在给定的约束条件下，达到极小。
模型识别的结果具有多解性。要识别的参数数目应少于总数据数目。也就是说必须要有已知量。已知量愈多，反求的参数愈精确，由此建立的模型的适用性就愈好。正因为模型识别结果的多解性，所以对于同一个问题，不同的人所求得的参数组合不同，甚至同一个人在不同的时间所求得的参数也不同。显然，模型识别的参数不一定是含水层所固有的参数。因此，有人称模型识别的参数为“模型参数”，以示区别。尽管模型参数不能完全反映实际系统的参数，但是模型参数有其特殊作用，它能够使得数学模型在行为和功能上代替实际的地下水系统，成为地下水系统的“复制品”。
2.模型检验
为了检验识别后的模型的可靠性，需要采用同一系统的另一时段的数据资料输入模型进行检验。如果计算结果符合实际资料，则可以说明模型能真实反映实际系统。需要指出的是，在模型识别和模型检验阶段所用的两组数据资料，必须是相对独立的不同时间段的资料。
模型灵敏度分析的目的是了解参数变化对计算结果的影响，同时识别重要参数。灵敏度分析一般在模型识别之前进行，也可以在模型识别之后进行。
选取要分析的一个参数（θ），然后固定其余参数，改变θ的数值分析计算结果。这时计算水头（g）就是θ的函数，即g=f（θ）。则有如下定义：在θ=θ0附近，水头变量g（θ）相对于原值g*（θ）的变化率和参数θ相对于θ0的变化率之比称为水头对参数θ的灵敏度，以下式表示：
西北内陆黑河流域水循环与地下水形成演化模式
四、黑河流域模拟区水文地质条件概化
地下水数量转化研究的数值模型模拟区，选择了张掖盆地和酒泉东盆地，包括张掖、临泽、高台的所有灌区和民乐及山丹的个别灌区，还有肃南县明花区，面积近9000 km2。
数值模拟区是只有侧向流入而没有侧向流出的山间断陷盆地，其间充填了巨厚的松散沉积物，构成赋存地下水的天然场所，为连续和统一的第四纪含水岩系综合体，周边山体为天然的地质边界。在张掖盆地，地下水自南东向北西运动，泄于黑河干流而流出区外。西部酒泉东盆地，地下水由南西向北东运动，榆木山至高台县城一线为两盆地天然汇水线。
数值模拟区地下水的主要补给来源是河水（含雨洪水）、渠系引水和田间灌溉水的垂直入渗，而泉水溢出、蒸发和人工开采是主要排泄方式。
据均衡计算结果，1999年区内补给量为11.94×108 m3，排泄量为14.09×108 m3，均衡差为-2.15×108 m3，数值模拟区处于负均衡状态，地下水水位呈下降态势。
数值模拟区周边皆为二类流量边界。山区边界沿山前大断裂分布，流入量主要为基岩裂隙水侧向流入和沟谷潜流。东部民乐、山丹断面和西部明花区断面为区外侧向流入量，利用断面法求得。南部新坝-红崖子隐伏断层使地下水流不连续，作为该段边界，概化的水文地质模型如图5-1。
五、数学模型概化
数值模拟区南半部为潜水、北半部为承压水，适宜采用潜水-承压水数学模型。但是各灌区开采地下水的程度不同，一些地带已将潜水与承压水连通，承压水头与潜水水位动态变化具有一致性。因此，将模型概化为非均质各向同性二维流潜水模型。鉴于区域面积大，地下水水位年变幅小，与含水层厚度相比可忽略，所以用导水系数（T）近似代替渗透系数（K）与含水层厚度（H）之积。
数学模型及定解条件如下：
图5-1 黑河流域模拟区水文地质模型概化图
西北内陆黑河流域水循环与地下水形成演化模式
式中：T——含水层导水系数（m2/d）；
μ——含水层给水度（无量纲）；
Wb——各项补给项强度之和（m3/km2·d）；
Wp——各项排泄项强度之和（m3/km2·d）；
q——流量边界单宽流量（m3/km2·d）；
Γ2——流量边界代号；
n——边界上的内法线方向。
采用线性插值，伽辽金有限元法解上述方程组，见程序框图（图5-2）。
图5-2 黑河流域数值模型模拟程序求解流程
六、定 解 条 件
（一）初始条件
以1999年水位统测结果为基础，结合地下水动态长观资料，绘制1月份等水位线图为初始流场。采用三角剖分法将计算区剖分成1421个单元，799个结点。其中内结点624个，边界点175个。水位观测点33个，均分布于结点上（图5-3）。同时尽量把结点布置在概化的灌区边界上。
（二）计算时段
以1999年元月初至12月末每个自然月实际天数为时段长度，全年共分12个时段。
（三）水文地质参数
根据黑河勘察报告研究成果，数值模拟区参数取值范围T值为100～6500 m2/d，μ值为0.1～0.25之间。参数分区以灌区为基础，按不同埋深划分。
（四）源汇项
计算区地下水主要靠河水、渠系引水、灌溉水、降水凝结水入渗及边界流入补给。消耗于蒸发蒸腾、泉水溢出和人工开采。有关参数的选取，主要依据黑河报告和各县水利部门研究成果，补给量与排泄量通过水量均衡方法计算求得。
由于数值模拟区范围较大，而且区内农业发达、干支渠密布，沿主要河流（黑河）引水口众多，所能收集到的水文和水利资料有限，所以剖分不宜过细，可将河水（含雨洪）、渠系水、灌溉水和降凝水入渗及人工开采处理为面状量，把各灌区不同埋深均衡计算结果以单位面状量进入模型，补给项为正，排泄项为负。非灌溉期（1～3月，10～12月）的渠系水和灌溉水入渗及人工开采量强度为0，灌溉期（4～9月）摊分全年入渗量。
图5-3 黑河流域数值计算区剖分图
1999年河水入渗量占当年黑河（莺落峡）径流量的32%，每月径流量占全年径流量的比例分配到12个时段。降水、蒸发强度按各月份所占全年比值分配到12个时段。1～3月和10～12月的降水为0，4～6月降水占30%，7～9月降水占70%。按地下水水位不同埋深，计算蒸发量，其中1～3月占13%，4～6月占41%，7～9月占35%，10～12月占11%。
泉水溢出带均分布于细土平原、地下水水位埋深小于3.5m的地带，各泉沟及黑河河床地下水水位高于河床标高，实际为线状量。但是因剖分单元较大，无法准确描述，所以将线状量处理成面状量，假设地下水水位埋深小于3.5m带为泉水溢出带，具体做法将所有结点地面高程减去3.5m，于是该区地下水水位埋深值为负。将1999年泉水溢出量除以该区面积，再除以平均水头差1.5m，获得单位水头差条件下泉水溢出强度，引入模型。然后根据各时段水头变化，获得不同时段的泉水溢出量。
数值模拟区边界为透水边界或弱透边界，均给出单宽流量，全年一致，不再按时段划分。
七、数 模 结 果
按上述补给与排泄要素及其参数，采用观测点的地下水水位拟合，对1999年实施模型进行识别。
区内共有观测点33个，集中在张掖、临泽、高台的细土平原带。在调参过程中，不断缩小拟合点误差，兼顾初始流场与计算流场形态一致，并且每个节点水位偏差不宜过大。调参结果，数值模拟区共有60个参数分区，如图5-4和表5-2所示。观测点拟合结果如图5-5和图5-6所示，地下水流场拟合情况如图5-7所示。
图5-4 黑河流域数值模拟参数分区图
表5-2 黑河流域模型采用的有关水文地质参数
Ⅱ 抽水试验求水文地质参数
2.4.6.1 抽水试验方法选择
抽水试验是地下水试验与求参数的常用方法。在以往的水文地质区域调查中，普遍使用的是稳定流抽水。稳定流抽水施工所需时间较短，操作简单。然而随着地下水资源研究程度的提高，稳定流已不能满足地下水资源研究的需求。这主要是因为稳定流抽水试验只能求取含水层水平渗透系数和导水系数。稳定流试验在抽水孔中进行，由于施工不当，或因抽水井水位波动大，甚至水花的飞溅等都会影响数据的准确性。而且稳定流计算结果是不能用来预测地下水资源动态变化的，而非稳定流抽水必须用一个孔组，数据在观测孔中测试。根据含水层特点，抽水试验资料选择不同的模型整理，不但可以求K、T，而且可以求给水度μ、垂向渗透系数Kz、弱透水层越流系数K'/m'、承压含水层弹性释放系数s、压力传导磨衡系数a等。因此获取的信息量比稳定流试验要多的多。
因此要求:
(1)偏远地区，施工比较困难，地下水开采程度低，地下水评价精度要求低的地区，可选择稳定流抽水求参。
(2)对于地下水资源评价精度要求比较高的地区，原则上都要选择非稳定抽水试验来求参。
2.4.6.2 稳定流抽水求参
2.4.6.2.1 抽水设计要符合裘布依公式
稳定流抽水试验主要是求渗透系数K，其准确程度取决于钻孔施工质量、选用计算公式、抽水引起的地下水运动规律、边界条件与裘布依公式的基本假设条件是否相符等。
裘布依(A.Dupuit)公式的基本假定为:
(1)含水层均质、水平;
(2)承压水顶底板是隔水的;潜水井边水力坡度小于1/4，底板隔水，抽水前地下水是静止的，即天然水力坡度等于零;
(3)半径R的圆柱面上保持常水头，抽水井内水头上下一致。
抽水过程中可能出现的问题是:大降深抽水出水量足够大时，井壁和周围含水层容易产生三维流，井周产生紊流，井壁附近潜水水力坡度增大，I＞1/4使裘布依假定失效等等。滤水管长度小于含水层厚度，井壁边界无法保持相等水头。在抽水后，形成下降漏斗，大部分含颂信水层不存在圆柱形常水头边界，距主孔很近的范围内(r≤0.178R)水位属对数关系。当观测孔距主孔距离r＞0.178R后，水位就变成贝塞尔函数关系，贝塞尔函数的斜率比对数函数小，因此观测孔越远，计算出的K值越大。当含水层具有越流渗透补给时，通过不同半径圆柱面的流量不等，离主井越近，流量越大，动水位与半径的贝塞尔函数成正比，所以有越流补给时，只有r≤0.178R时，裘布依才是适用的。在天然径流条件下，等水位线不是一个同心圆，一般是下游半径较长的椭圆形。观测孔取得的降深是角度θ的函数，即上游偏小，下游偏大，只有在垂直地下水水流方向上的降深值无变化，因此观测孔的布置方向也是影响K值的因素之一。
在实际工作中，建议使用的抽水设计方法是:
(1)采用较小降深抽水;
(2)观测孔距主井适宜的范围是:1.6M≤r≤0.178R(其中:R为引用半径，M为含水层厚度);
(3)每个抽水试验一般要做3个降深，抽水试验最好安排在地下水非开采期，并将抽出的水引出试验区外，以免干扰水位下降。
2.4.6.2.2 稳定流常用计算公式
(1)承压含水层完整井单孔:
地下水资源调查评价技术方法汇编
(2)承压含水层完整井单孔二次以上降深:
地下水资源调查评价技术方法汇编
其中:二次降深，
三次降深， (Qi为三次降深的三个流量，Swi为三次降深的抽水井水位降深)。
式中:Q———抽水井出水量(m3/d);
K———渗透系数(指水平渗透系数)(m/d);
R———影响半径(m);
rw———抽水井半径(m);
Sw———抽水井水位降深(m);
S1、S2———观测孔水位降深(m);
M———含水层厚度(m);
h———动水位至含水层底板深度(m)。
(3)承压含水层完整井有一个观测孔:
地下水资源调查评价技术方法汇编
(4)承压含水层完整井有二个观测孔:
地下水资源调查评价技术方法汇编
式中:h1、h2———含水层底板至观测孔水位降深高度;
r1、r2———抽水孔至观测孔距离，其他同上。
(5)承压含水层完整井岸边抽水(单孔，b＜0.5R):
地下水资源调查评价技术方法汇编
(6)承压含水层完整井岸边抽水(有一个观测孔，位于近河一边):
地下水资源调查评价技术方法汇编
式中:b———抽水孔距河岸距离，其他同上。
(7)承野游轮压含水层非完整井(单孔，井壁进水):
地下水资源调查评价技术方法汇编
式中:l———观测孔底至含水层顶板距离。
(8)承压含水层非完整井(一个观测孔):
地下水资源调查评价技术方法汇编
式中:l———观测孔底至含水层顶板距离，等于过滤管有效进水长度。
(9)承压含水层非完整井(单井、井壁井底进水):
地下水资源调查评价技术方法汇编
(10)潜水-承压水完整井(单井):
地下水资源调查评价技术方法汇编
(11)潜水完整井(单孔):
地下水资源调查评价技术方法汇编
式中:H———含水层厚度。
(12)潜水完整井(一个观测孔):
地下水资源调查评价技术方法汇编
(13)潜水非完整井(单井):
地下水资源调查评价技术方法汇编
含水层厚度很大时，应计算有效带厚度代替含水层厚度。
2.4.6.3 非稳定流抽水求参
2.4.6.3.1 非稳定抽水试验的设计
地下水非稳定流理论对含水层抽水过程的认识与稳定流理论的不同之处主要在于，非稳定流理论将含水层看作弹性体，在无限边界含水层中抽水时，整个流场的各运动要素是随时间而变化的，即流向钻孔的地下水是非稳定的流动。经过一定时间后地下水流才趋于稳定流动。非稳定流理论的基本公式———泰斯(C.V.Theis)公式的基本假设条件是:
(1)含水层均质、等厚、水平埋藏。
(2)没有垂向和水平补给。
(3)地下水初期水力坡度为零。
(4)地下水是平面流。
(5)含水层在平面上是无限边界。
泰斯公式与裘布依公式比较，其优点在于反映了地下水运移普遍存在的非稳定过程，公式中考虑了时间因素，因此在一定条件下可以预测含水层中任一点的水位降深及降落漏斗展布的范围。有利于求取除K、T以外的其他参数，如弹性释水系数sa(潜水为给水度μ)、压力传导系数a等。根据泰斯公式发展的其他模型和计算公式，还可计算弱透水层越流系数K'/M'、垂向渗透系数Kz等。
抽水试验设计须考虑的主要方面有:
(1)抽水前要进行试抽，了解抽水孔的出水量、水位降深和观测孔水位降深情况，选择一个较小的适当流量，以免抽水时掉泵和形成大降深。在1.6M≤r≤0.178R处设置观测孔，以避免三维流、紊流和远处计算K值偏大等问题的干扰。
(2)观测孔设置在垂直于地下水流动的方向上。
(3)抽水试验选择时间段内周边地区无地下水开采，抽水井抽出水量引出区外，避免引起对水位降深的干扰。
(4)抽水流量必须保持基本稳定，最大流量与最小流量之比不应大于1.05。
(5)抽水时间的长短，要根据抽水过程中所绘制的水位降深(S)与时间(t)的双对数曲线所显示的抽水阶段来决定。当曲线平稳的第二阶段末期出现曲线上翘，显示达到第三阶段后，再略延长一段时间抽水试验就可结束。所需抽水时间的长短与含水层岩性有关。
2.4.6.3.2 承压完整井非稳定流抽水求参
非稳定承压完整井计算公式:以固定流量Q抽水时，距抽水井距离r处任一时间t的水位降深，可简化为:
地下水资源调查评价技术方法汇编
(1)试算法。
压力传导系数a，导水系数T，渗透系数K，弹性释水系数s，t1、t2时刻测得抽水孔水位降S2，观测孔水位降S1
地下水资源调查评价技术方法汇编
用此公式通过试算法求a。
设 为纵坐标，a为横坐标。用已知观测时间t1、t2和任意给定的a1、a2、…、an代入上式，求相应的β1、β2、…、βn值，绘制β=f(a)关系曲线。根据抽水孔、观测孔实测所获得的S1、S2，得实测
β=f(a)关系曲线上得到实际a值。将所计算的a值代入上述S1或S2计算公式中求得导水系数T，渗透系数 弹性释水系数
为避免作图的不方便，注意时间t，采取抽水2h后观测，且t1、t2间隔不小于4～5h(图2.4.3)。
(2)降深-时间双对数量板法:
非稳定流计算公式:
地下水资源调查评价技术方法汇编
图2.4.3 试算法关系曲线
地下水资源调查评价技术方法汇编
(2.4.35)式至(2.4.38)式式中:
——井函数自变量;
S(r，t)———距抽水孔r处，任一时间(t)的水位降深;
T=K·M———导水系数;
——压力传导系数;
r———观测孔距抽水孔距离;
sa———弹性释水系数;
K———渗透系数;
W(u)———井函数，可查表。
配线的做法是:
(1)将观测孔不同时间测得的水位降深值，点绘在透明的双对数纸上。然后将对数纸重叠在理论标准曲线(即量板)上。使实测点完全重合在理论标准曲线上(注意:对数纸与量板要采用同一模数，且纵、横坐标必须平行)。
(2)读出相应的W(u)、S和1/u，t值代入S(r，t) 式中求得T、a。随之又可求出K、S。此方法主要用于一个观测孔。
(3)降深-距离双对数量板法。
与降深-时间曲线法一样，点绘同一时间各观测孔S-r2关系曲线，重叠在W(u)-u理论曲线上(注意纵横坐标平行)，求a、T以及K、S。
本方法主要用于有数个观测孔的条件下。
(4)直线解析法(图2.4.4)。
设在t1时间测定降深S1，t2时间测定降深S2，有S2-S1=ΔS
地下水资源调查评价技术方法汇编
图2.4.4 S-lgt曲线
当ΔS=0时，t1=t0有:
地下水资源调查评价技术方法汇编
同样 求出渗透系数和弹性释水系数。
采用直线解析法常因人为误差导致直线斜率和截距的不准确，而影响计算结果。实际工作中可用最小二乘法推求直线方程斜率和截距后，再用上述方法求参。
(5)水位恢复法。
此方法优点是排除了抽水过程中的一些干扰因素，是常被采用的方法。计算公式是:
地下水资源调查评价技术方法汇编
得T、a后，同样也可求出K、S。
2.4.6.3.3 承压非完整井非稳定流抽水求参
非完整井抽水时，水流越接近井孔，流线越弯曲集中，其运动状态不符合泰斯公式平面流的假设条件。但当观测孔布置在距抽水孔r≥1.6M时，地下水流线趋于平行，因此在r≥1.6M距离处的观测孔内取得的不同抽水时间t和相应水位降S值，同样可以利用泰斯公式计算T、a值。
根据抽水资料绘制S=f(lgt)曲线(图2.4.5)，在曲线上任意两点P1、P2，解得该曲线P1、P2两点斜率(m1、m2):
地下水资源调查评价技术方法汇编
图2.4.5 S-lgt曲线
地下水资源调查评价技术方法汇编
式中:m1、m2———S=f(lgt)曲线上相应lgt1，lgt2点的斜率;
t1、t2———测得观测孔水位降深S1、S2时的时间。
2.4.6.3.4 潜水完整井非稳定流求参
潜水抽水时，由于孔隙水具有延迟重力排水作用，所以瞬时释放水量的假定是不适宜的。在抽水开始很短的早期，降深很小时，可以认为存在弹性释放水量。随着抽水时间的延长，含水层出现延迟释放水量的情况，我国大部分孔隙含水层中已被证实大都属于这种类型，因此不考虑延迟释水的计算方法常常使计算结果不合理。
这里推荐较符合大部分平原(盆地)的冲洪积、冲湖积沉积的孔隙含水层条件，在实践中反映比较有效的、考虑延迟给水的布尔顿、纽曼和二元结构模型，以供参考。
(1)潜水布尔顿(S.N.Boulton)公式。
含水层均质、等厚，底板水平埋藏，考虑含水层滞后重力释水。
布尔顿模型的计算公式为: 为潜水完整井布尔顿井函数。
抽水前期
地下水资源调查评价技术方法汇编
抽水后期
地下水资源调查评价技术方法汇编
(2)纽曼(S.P.Neuman)公式。
含水层不厚，各向异性，潜水面无垂向补给，水位降远远小于含水层厚度，考虑了抽水时含水体内垂直方向水力梯度变化。计算公式为:
(ts.y，β);Sd(ts.y，β)为潜水完整井纽曼模型井函数。
前期
地下水资源调查评价技术方法汇编
后期
地下水资源调查评价技术方法汇编
因此，纽曼模型还可以计算垂向渗透系数Kz。
式中:Kr———水平渗透系数;
Kz———垂向渗透系数;
ss———比弹性释水系数，Ss=Sa/M，M为含水层厚度;
sa———抽水前期弹性释水系数;
sy———抽水后期水位变动带延迟释水率(相当于μ);
r———观测孔与抽水孔距离;
S———观测孔水位降深;
Q———抽水孔抽水量。
(3)二元结构计算公式。
潜水-微承压水含水层分为上下两个部分，上部为弱透水层潜水，有自由水面，垂向渗透系数Kz，水位变动带释水率sy，弱透水层厚度M'，水位降深S'下部为微承压含水层，其厚度M，弹性释水系数sa，导水系数T，水头略高于弱透水层自由水面。抽水时，下部弱承压含水层有汇点径向流，水头迅速下降，与自由水面逐渐合成一体。上部弱透水层向下释水补给下部微承压含水层。我国平原中许多地区存在这种上细下粗的二元含水层结构和水动力特征。
下部微承压含水层水位降深的计算公式为:
前期
地下水资源调查评价技术方法汇编
后期
地下水资源调查评价技术方法汇编
式中:前期
地下水资源调查评价技术方法汇编
后期
地下水资源调查评价技术方法汇编
用S-lgt双对数量板法，采用S·N·布尔顿、S·P·纽曼和二元结构计算公式求参，都可以得到较满意的结果。这里以布尔顿公式为例，简述其方法和注意事项。
主要步骤(图2.4.6):
(1)将抽水资料用双对数纸点绘lgS=f(lgt)曲线，并绘在标准曲线A上。注意纵横坐标保持平行，尽可能将初期曲线与标准曲线A重合。
(2)记下重合曲线上 值，任选一点并在标准曲线上读出S、1/ua、 及t坐标值，求出T、sa。
(3)将资料曲线沿水平方向移动，尽可能使资料后期曲线与标准曲线Y重合(注意曲线前段r/D值与后段r/D值一致)，同样读出 、t值，求出T、sy。
图2.4.6 非稳定流潜水标准曲线图
以上步骤同样可以应用到纽曼公式和二元结构公式中，只要采用相应的井函数。前期与后期水位降公式以及各自标准曲线特征值 即可。同样要注意前期曲线与后期曲线配线时要在同一特征值的标准曲线上。只要认真按上述步骤操作，一般双对数量板法计算结果较为满意。
2.4.6.3.5越流含水层求参
(1)承压含水层受上部弱透水层补给，弱透水层储水系数忽略不计。有一个抽水孔，一个观测孔(必须打入越补含水层中)任一点水位降的解为:
地下水资源调查评价技术方法汇编
地下水资源调查评价技术方法汇编
(2)考虑弱透水层释水，越流供给层为弱透水层，可位于越流层之上或之下。任一点水位降的解:
地下水资源调查评价技术方法汇编
式中: 为井函数自变量;
m'———弱透水层厚度;
K———越流含水层渗透系数;
K'———弱透水层渗透系数;
S———任一点水位降深;
M———越补层厚度;
T———导水系数;
a———导压系数;
s———越补层释水系数;
s'———弱透水层释水系数;
rm———抽水孔的半径;
r———计算点与抽水孔轴心的距离;
k'/m'———越流系数。
Ⅲ 水文地质参数变化
一、太原盆地水文地质参数计算
水文地质参数的选取直接影响着地下水资源计算量的大小和可信度，研究水文地质参数具有十分重要的意义。本次相关的水文地质参数主要有降水入渗补给地下水系数（α）、潜水蒸发极限深度（L）、蒸发强度（ε）、灌溉回渗地下水系数（β）、疏干给水度（μ）、导水系数（T）、弹性储水系数（s）、渗透系数（K）、河流渗漏补给系数、渠系渗漏补给系数等。
（一）降水入渗补给地下水系数（α）
影响降水对地下水的补给量的因素很多，主要有地形、包气带岩性及结构、地下水位埋深、降水特征及土壤前期含水量等。
降水入渗补给系数为降水入渗补给地下水量与降水量之比值。年降水入渗补给系数为年内所有场次降水对地下水入渗补给量总和与年降水总量的比值，其表达式为:
山西六大盆地地下水资源及其环境问题调查评价
式中:α年是年降水入渗补给系数；pri是场次降水入渗补给量，mm；P是年降水量，mm；n是年降水场次数。
用长期动态观测孔求取年降水入渗系数的计算方法:
山西六大盆地地下水资源及其环境问题调查评价
式中:μ∑Δh次是年内各次降水入渗补给地下水量之和；P年是年降水量；Δh次是某次降水引起的地下水位升幅值。
根据动态资料分析计算，在前人试验的基础上，综合考虑各方面的因素，给出盆地区降水入渗补给地下水系数（详见第四章）。
（二）地下水蒸发极限深度（L）、蒸发强度（ε）
蒸发极限深度就是指浅层水停止蒸发或蒸发量相当微弱时，浅层水位埋深值。蒸发强度就是在极限蒸发深度以上，单位时间浅层水的蒸发量。
影响地下水蒸发的主要因素是地下水位埋深、包气带岩性和水面蒸发强度等。
理论上，当水位埋深处于蒸发极限深度时，地下水在无补给、无开采的条件下，动态曲线近于平直。
地下水蒸发极限深度（L）
蒸发极限深度通常采用迭代法、试算法和经验公式计算（L），公式如下:
迭代法:
试算法:
经验公式法:
式中:ΔT1、ΔT2为计算时段，d；H1、H2、H3为时段内水位埋深，m；Z1、Z2为时段内水面蒸发强度，m/d；
经计算，太原盆地孔隙水区不同岩性的蒸发极限深度依包气带岩性不同分别为:亚砂、亚粘土互层为3.5m，亚砂土为4.0m，粉细砂、亚砂土互层为4.5m。
地下水蒸发强度
计算公式:
式中:Z0是液面蒸发强度，mm/d；ΔH是浅层水降落间段的平均水位埋深，mm；Z是蒸发强度，mm/d。
由本区浅层水水位埋深图（详见第四章）可看出，水位埋深小于4m的区域在北部太原市和南部平遥、介休一带，根据上式计算太原、平遥、介休等地的地下水蒸发强度见表3-1。
表3-1 太原盆地孔隙水区地下水蒸发强度
（三）灌溉回渗地下水系数（β）
是指田间灌溉补给地下水的量与灌溉总量的比值。影响灌溉回渗系数和因素主要有岩性、水位埋深、土壤含水率、灌溉定额等多种。
计算公式:
式中:μ是给水度；Δh是由灌溉引起的地下水位平均升高值，m；Q是灌溉水量，m3；F是面积，m2。
本次工作在盆地太原市小店区郜村、汾阳市贾家庄镇东马寨村和榆次市杨盘等3个地方布置了3组灌溉入渗试验，地表岩性郜村为粉质粘土、东马寨上部为粉质粘土，下部为粉土，杨盘为粉土，化验室给水度试验结果分别为0.195、0.11、0.143。郜村在37m×37m的面积上布置10眼观测孔，水位埋深1.2～1.3m，累计灌溉水量160m3，10个孔平均水位上升值为0.1912m，根据上式计算得灌溉入渗地下水系数为0.32；东马寨村水位埋深1.95～2.44m，在26m×26m的面积上布置10眼观测孔，灌溉水量60m3，观测孔平均水位上升值为0.465m，计算得灌溉入渗地下水系数为0.58；杨盘布3个观测孔，水位埋深5.76～6.01m，灌溉面积100m2，灌溉水量100m3，平均水位上升高度为0.27m，计算得灌溉入渗系数为0.039。
从以上试验数据可以看出，不同水位埋深、不同岩性地区灌溉入渗系数有很大区别。综合考虑各种因素，灌溉回渗地下水系数选用值见表3-2。
表3-2 灌溉回渗地下水系数
（四）弹性贮水系数S、导水系数T、给水度μ、渗透系数K
盆地区大部分地区都进行过1∶5万比例尺的农田供水水文地质勘查，做过大量单孔和多孔抽水试验，本次在文水文倚、汾阳等5地分别作了5组抽水试验，用非稳定流公式，降深-时间半对数法计算结果如下:文倚导水系数T＝1983.59～2181.95m2/d，渗透系数K＝32.19～35.4m/d，弹性贮水系数S＝1.79×10-3；汾阳县贾家庄镇东马寨村抽水试验求得导水系数T＝325.84～376.5m2/d，渗透系数K＝5.65～6.53m/d。结合以往本区的工作成果，给出太原盆地浅层孔隙潜水和中深层孔隙承压水水文地质参数，详见参数分区图3-13和参数分区表3-3。
表3-3 太原盆地中深层孔隙承压水及浅层孔隙潜水参数分区
图3-13 太原盆地参数计算分区图
二、大同盆地水文地质参数计算
由本区浅层水2004年水位埋深图可看出，水位埋深小于4m的区域主要分布于盆地中部冲积平原区，盆地南部怀仁、山阴、应县、朔州分布面积较大。根据计算和以往试验资料，本区蒸发强度确定值见下表（表3-4）。
表3-4 大同盆地孔隙水区地下水蒸发强度
据“山西省雁同小经济区水资源评价、供需平衡研究报告”中搜集的本区灌溉回渗试验数据取得不同水位埋深、不同岩性、不同灌溉定额的灌溉回渗系数，灌溉回渗系数选定值见表3-5。
盆地区大部分地区都进行过1/5万比例尺的农田供水水文地质勘查，做过大量单孔和多孔抽水试验。本次工作搜集本区以往抽水试验孔117个，本次在大同县党留庄乡、怀仁县金沙滩镇、怀仁县新发村、怀仁县榆林村、山阴县张庄乡、朔州市城区沙塄乡等6地分别作了6组抽水试验，采用AquiferTest计算程序，非稳定流方法计算，本次抽水孔具体情况和计算结果见表3-6和表3-7 。
表3-5 灌溉回渗地下水系数
表3-6 大同盆地本次抽水试验数据统计
表3-7 大同盆地本次抽水试验计算成果表
结合以往本区的工作成果，给出大同盆地浅层孔隙潜水和中深层孔隙承压水水文地质参数，详见参数分区图3-14、图3-15和参数分区表3-8、表3-9 。
图3-14 大同盆地降水入渗系数分区图
图3-15 大同盆地浅层、中深层孔隙水参数分区图
表3-8 大同盆地浅层孔隙潜水参数分区表
续表
表3-9 大同盆地中深层孔隙承压水参数分区
三、忻州盆地
忻州盆地地下水资源较为丰富，开采条件优越，20世纪70年代之前地下水开采规模较小；70年代初至80年代末随着农业灌溉的普及，工业生产的发展和城市规模的扩大，地下水开采量迅速增加。开采对象以浅层水为主，造成浅层水水位普遍有所下降（但下降幅度不大）。从20世纪90年代至今，虽然地下水开采量具有逐年增大的趋势，但增加幅度较小，且中层井数量逐渐增多，形成了浅层水、中层水混合开采的新模式，地下水位总体处于动态平衡状态。受地下水人工开采的影响，降水入渗系数及导水系数等水文地质参数发生了一定程度的变化。
区内降水入渗系数的变化除了与年降水量及降水特征有关外，主要与浅层地下水位埋深关系较为密切。已有资料表明，在山前倾斜平原区，浅层水位埋深一般大于7m，因水位下降使降水入渗系数发生了不同程度的减小。在冲积平原区浅层水位埋深一般小于7m，水位下降的结果引起了降水入渗系数有所增大。不同地貌单元降水入渗系数的变化见第五章。
从20世纪70年代以来，区内含水层的导水系数发生了较为明显的减小，主要体现在因浅层地下水位下降，使浅层含水层上部处于疏干状态，含水层厚度减小，直接导到导水系数减小。因浅层水水位下降幅度不同，导水系数减小的程度也存在差异，从本次地下水侧向补给量计算断面附近的井孔资料分析，含水层厚度一般减小了3～6m，导水系数由70年代中期的60～250m2/d，减少到目前的50～200m2/d左右。
忻州盆地给水度根据不同地貌单元含水层岩性、分选性及富水性综合确定见表3-10及图3-16 。
表3-10 忻州盆地浅层含水层给水度分区
图3-16 忻州盆地给水度分区图
四、临汾盆地
经过搜集以往资料，调查和计算确定临汾盆地降水入渗系数见表3-11。临汾盆地渗透系数及给水度分区见图3-17，表3-12。
表3-11 临汾盆地平原区降水入渗系数统计
图3-17 研究区渗透系数及给水度分区图
表3-12 临汾盆地参数分区表
五、运城盆地
运城盆地地下水长观网建站年代较远，积累了大量的地下水位监测资料，且经过多次的地质、水文地质勘察、地下水资源评价工作，取得了大量的降水入渗值，参考前人综合成果，结合目前包气带岩性、地下水位埋深，给出运城盆地降水入渗补给系数，见表3-13。
表3-13 运城盆地平原区降水入渗系数统计
渠系有效利用系数除受岩性、地下水埋深影响外，还与渠道衬砌程度有关。修正系数r为实际入渗补给地下水量与渠系损失水量Q损的比值，是反映渠道在输水过程中消耗于湿润土壤和侵润带蒸散损失量的一个参数，它受渠道输水时间、渠床土质及有无衬砌、地下水埋深等因素的影响。一般通过渠道放水试验获得。本次评价主要参考运城市水利局相关试验成果，见表3-14。
表3-14 运城盆地万亩以上灌区η、r、m值统计
灌溉回归补给系数β值与岩性、植被、地下水埋深及灌溉定额有关，一般通过灌溉入渗试验求得，本次评价主要参照运城市水利部门资料综合确定，详见表3-15。
表3-15 运城盆地灌溉回归系数β取值
河道渗漏补给系数是河道渗漏补给地下水量与河道来水量的比值。其值大小与河床下垫面岩性、流量、地下水位埋深及渗漏段长度有关。运城盆地沿中条山前发育数条季节性河流，河床下垫面主要为砂卵砾石，当洪雨季节，地表河床水位远高于地下水位，为地表水的入渗造就了十分便利的条件。根据河道渗漏资料，可建立如下数学模型:
山西六大盆地地下水资源及其环境问题调查评价
式中:m河是河道渗漏补给系数；A是计算系数，A＝（1-λ）×（1-φ）L，φ是单位千米损失率；L是河道渗漏长，km，Q径是河道来水量，m3/s。
据运城市水利部门研究成果，A值约为0.090。
含水层的渗透系数主要由野外抽水试验通过稳定流及非稳定流计算公式求得，各勘探部门在运城盆地先后进行过各种勘察，进行了大量的抽水试验工作，积累了丰富的资料，参考本次抽水试验成果对以往参数进行了修正，取值结果见表3-16 。
表3-16 运城盆地松散岩类K值选定表
降雨入渗补给系数在同岩性、同降雨量情况下，随地下水位埋深的增大，降雨入渗补给系数会达到一个最大值之后趋于减少或变为常数。运城盆地北部的峨嵋台塬及闻喜北塬，其地下水位埋藏深，地表主要以黄土类为主，降水入渗主要依靠黄土垂直节理裂隙及“流海缝”以“活塞式”注入地下，多年来其降水入渗系数基本为常量，经用动态分析法计算其降水入渗系数在0.108～0.11间；在盆地中部的冲湖积平原区，其地表岩性主要以Qp3＋Qh冲湖积相的亚砂土、亚粘土、粉细砂为主，由于开采强烈，区域水位严重下降，地表数米至几十米内均为饱气带，为降水入渗准备了调蓄空间，加强了降水向地下水的转化。根据盆地地下水长观孔资料及次降雨资料，计算出盆地冲湖积平原地带，降水入渗系数在0.1～0.162之间，总体上上游大于下游。而在东部及南部的山前倾斜平原区，地下水位埋深一般大于5m、乃至几十米，地表岩性大多为亚砂土及亚粘土，尤其是在一些沟口附近，从地表往下几十米范围内为干砂卵砾石，一般降雨基本上不产生地表径流，这无疑加大了降水的转化。据相关资料计算，降水入渗系数高达0.21～0.30。因过去所做的工作不系统，没有对降雨入渗系数进行系统分类，不便比较，但根据运城盆地饱气带岩性、地下水变动情况，除峨嵋台塬及黄土丘陵区变化不大外，其他地区降雨入渗系数无疑有增大趋势。
盆地内抽水井的含水层，大多为数个含水层混合开采。现根据本次抽水计算值，对历次研究成果中的K值加以修正，得出运城盆地各个地貌单元的渗透系数。总体来说，黄河岸边低阶地区K值最大为11.3～14.6m/d，中条山山前倾斜平原次之，为5.45～6.12m/d，最次为闻喜北垣K＝1.10m/d左右。
根据地貌单元、含水层岩性、地下水水力特征及各参数特征，将运城盆地划分为10个参数分区，见表3-17及图3-18。
表3-17 运城盆地水文地质参数分区
六、长治盆地
根据水文地质条件，长治盆地参数分区见图3-19，表3-18 。
图3-18 运城盆地水文地质参数分区表
图3-19 长治盆地参数分区图
表3-18 长治盆地浅层孔隙潜水参数分区
（一）降水入渗补给系数变化
根据《太原市地下水资源评价报告》研究成果，盆地区亚砂土、极细砂、细砂的降水入渗系数随着地下水位埋深的增大而增大，当水位埋深超过一定值以后，降水入渗系数开始趋于稳定；降水量越大，降水入渗系数在相同的岩性和地下水位埋深条件下也越大。对于亚砂土、极细砂、细砂在相同水位埋深和降水情况下，细砂的降水入渗系数>极细砂的>亚砂土的。总体来说，颗粒越粗，降水入渗系数也越大。
α随降水量的变化，非饱和带在降水入渗补给地下水过程中起调节作用，降水入渗补给过程要滞后于降水过程，其滞后时间的长短、特征与非饱和带的重力水蓄水库容关系密切，地下水埋深越大，其蓄水库容也越大，调节能力也越强，滞后现象也越明显。
在亚砂土、极细砂和细砂3种岩性中，降水量相等时，降水入渗系数从大到小的顺序为细砂、极细砂、亚砂土。场次降水量的影响表现为α次先是随着降水量的增大而变大，当降水量超过一定数值后，α次反而呈减少趋势，这个降水量即是最佳降水量。α年与α次有相同的规律性，从入渗机制分析，α年也存在最佳年降水量。
当地下水埋深为零时，降水入渗补给系数亦为零，然后随埋深的增加由小变大；当地下水埋深到达某一定值时，降水入渗补给系数达到最大值即最佳降水入渗补给系数，并由此随埋深的增加由大到小，到达一定的埋深时，趋于定值。地下水埋深对降水入渗补给系数的影响，可从3方面来说明。
埋深反映了蓄水库容的大小。当埋深为零时，即蓄水库容为零，这时无论降水量多大，均无入渗补给的可能。当埋深增加时，地下水库得到了降水入渗补给量，此时降水入渗补给系数大于零，降水入渗补给系数随埋深的增加而增大。当地下水达到最佳埋深时，其对应的降水入渗补给系数为最佳降水入渗补给系数，原因是由于条件一致的地区中的依次降水，其入渗补给量随地下水埋深的变化必存在一个最大值。当地下水埋深较小时，由于地下水蓄水库容较小，形成蓄满产流，不能使降水全部入渗；当地下水埋深再增大时，则损失较最佳埋深为大，故降水入渗补给系数随埋深的增加而减小。对于不同级别的降水量，α最大值出现的地下水位埋深区域也不同。最佳埋深与岩性和降水量有关。
地下水埋深在某种程度上反映了土壤水分的多少。土壤水垂直分布大体可概化为3种状况。第1种情况是地下水埋深较小，毛管上升水总能到达地表；第2种情况是地下水埋深较大时，毛管上升水无法到达地表；第3种情况是地下水埋深介于两者之间，在此埋深内，由于地下水位是升降变化，毛管上升水有时达到地表，有时达不到地表。这3种情况将对降水入渗补给量有不同的影响。第1种情况，降水一开始，水即可通过毛管在重力作用下迅速向下移动，地下水位在降水开始后很快上升。第2种情况，降水首先应满足土壤缺水的需要，而后在重力作用下通过空隙下渗补给地下水。其渗漏途径较第1种情况长，入渗方式也有差异。
图3-20 渗透系数与深度关系图
不同地下水位埋深条件对降水入渗补给系数取值的影响。盆地太谷均衡实验场的水分势能实验最大深度为8.2m，有观测点41个。多年资料的分析结果表明，土壤水分势能变化从地面往下可分为3个变化带———剧烈变化带、交替变化带和稳定带，剧烈变化带埋深为0～1.1m，土壤水分势能变幅大于200×133Pa；交替变化带埋深1.1～3.6m，土壤水分势能变幅大于（100～200）×133Pa之间；埋深3.6m以下为稳定带，其土壤水分势能变幅小于100×133Pa，其中埋深在4.5～5.0m以下的稳定特性更为明显，其土壤水分势能的变幅一般不超过50×133Pa，其土壤水分全年为下渗状态。表明埋深在5.0m以下为稳定入渗补给，反映在降水入渗补给系数上随埋深增加，α年将趋于稳定，故当埋深大于5.0m时，α年值可取定值，不再随埋深而变化。原因是地下水埋深已到达或超过地下水极限埋深，损失趋于定值，水分不向上运动，必然向下运动，故形成了降水入渗补给系数随地下水埋深变化的稳定值。
（二）渗透系数变化
孔隙含水介质的渗透能力不仅取决于粒径大小、颗粒级配、胶结程度，还与其埋深有关。同一岩性的孔隙含水介质，随着深度的增加，介质被压密，渗透系数会减小。
根据河北平原山前冲洪积扇扇顶区数百个钻孔资料的统计，各种含水介质的渗透系数随埋深增加呈指数衰减，部分深层不同岩性渗透系数随埋深的变化规律参考下述经验公式:
岩性为卵砾石时，渗透系数与埋深关系式:
K＝K0e-0.0131h R＝0.877
岩性为砂砾石时，渗透系数与埋深关系式:
K＝K0e-0.0116h R＝0.869
岩性为中粗砂时，渗透系数与埋深关系式:
K＝K0e-0.0057h R＝0.896
K为埋深处的渗透系数；K0为地表浅层的渗透系数；h为埋深；R为相关系数。
因此，对于同一种岩性，其渗透系数大小与深度有关（图3-20）。
Ⅳ 防洪评价计算跨河桥梁工程
1背景
防洪影响分析主要包括以下的工作内容：通过桥梁项目建设区自然地理、社会经济及工程现状的认真分析，深入分析并确定流域河道堤防的水利规划。根据建设项目的情况、可能带来的影响、所在河段的水文泥沙特性、防洪评价的主要任务，选取有代表性的水文系列进行工程实施后的冲刷与淤积计算。计算水文系列的选取要能反映冲刷和淤积的不利水、沙条件组合。根据建设项目的基本情况、所在河段的防洪任务与防洪要求、防洪工程与河道整治工程布局及其它国民经济设施的分布情况等，以及河道演变分析成果、防洪评价计算，对建设项目的防洪影响进行综合评价。
根据建设项目雍高和冲刷与淤积分析计算，分析工程对河道蔽慧含行洪安全的影响范围和程度。对施工方案占用河道过水断面的建设项目，还需根据施工设计方案及工期的安排，分析工程施工对河道泄洪能力的影响。结合工程的施工进度安排，合理分析不同时段的防洪要求。结合河道冲淤变化的计算，评价项目建设是否会影响河势的稳定；对工程施工临时建筑物可能影响河势稳定的建设项目，应根据有关计算，分析工程施工其对河势稳定的影响。根据调查及有关部门提供的资料，分析项目建设对第三人合法水事权益的影响。
2防洪评价计算
2.1资料选用和设计水位的计算方法
根据《河道管理范围内建设项目防洪评价报告编制导则（试行）》等规范的有关规定，在进行防洪评价编制时，所在河段有规划水文成果的，应依据流域规划。路线纵坡的确定、桥涵方案的选择必须要设计水位。计算的过程一般采用试算法或者绘制水位—流量关系曲线求得。当桥位断面与水文计算断面间的形态断面较规整，突变点不多，河段顺直，河底纵坡均一时，宜按水文断面的水位~流量关系曲线，以设计流量确定设计水位后，利用水面比降推算出桥位计算断面的设计水位。当桥位计算断面和水文计算断面上、下游有卡口、人工建筑物或断面形状和面积相差较大、河底纵坡有明显曲折时，宜采用水力学方法的试算法求算设计流量时的水面线，推求桥位断面的设计水位。
2.2壅水高度计算
大桥建成后，桥梁桥墩布置于渠槽内，阻水明显，使桥位附近渠段水位、流速发生变化，断面过水面积将会减少，从而造成桥梁上游洪水位一定的雍高。
2.3壅水曲线长度计算
壅水曲线长度近似估算式为：L=2△Z/I式中：I———计算河段天然河底比降△Z———桥前最大壅水高（m）
2.4冲刷与淤积分析
计算天然状态下，由于流域的来水、来沙及河床边界条件的不断变化，河床形态总是处在不断的冲淤变化过程中。但在相当长的一个时段内，冲淤量可以相互补偿，渠道处于一个相对的动力平衡状态。渠道上建桥后，破坏了原有的这种平衡状态，由于压缩水流，致使桥下流速增大，水流挟沙能力增强，在桥下产生冲刷。随着冲刷的发展，桥下河床加深，过水面积加大，流速逐渐下降，待桥下流速降低到河床质的容许不冲刷流速时，河道内达到新的输沙平衡状态，冲刷停止。桥墩墩台附近河床床面总的冲刷深度，应为河床演变、一般冲刷和局部冲刷深度的综合。
2.4.1一般冲刷计算根据地质资料，冲刷计算采用采用《公路桥位勘测设计规范》推荐的计算公式，冲止流速公式一般用于桥梁的一般冲刷计算，其基本假定是桥下河槽断面内任一垂线平均流速达到冲止流速时，此垂线处冲刷停止，至此，一般冲刷深度达到最大，可近似计算河道冲刷，公式形式如下：
2.4.2局部冲刷计算流向桥墩的水流受到桥墩的阻挡，水流的绕流使流线急剧弯曲，形成局部冲刷坑。冲刷过程中，坑底流速逐渐降低，坑底泥沙逐渐粗化，使水流的冲刷作用与床沙的抗冲作用趋于平衡，局部冲刷坑达到最深。计算采用《公路桥位勘测设计规范》推荐的河床桥墩局部冲刷公式。
3结语
防洪评价管理是一项复杂的系统工程，涉及多部门、多专业碧手、多学科，情况复杂，许多问题仍需在今后的工作中不断总结和完善。关于分析计算是否需要采用模型分析，一些计算经验公式与实际的误差，数学模型与物理模型之间的误差等，以及哪种方法更实用和更符合实际，仍需大量的实践以进一步检验，需要在今后工作和工程实践中不断总结和分析。今后如有可能，在对大量工程实际观测资料分析的基础上，应对公宏笑式和模型实验结果进行验证分析，不断对计算方法和公式进行修正和改进，使之更符合实际要求。随着工程实践的不断深入，防洪评价计算方法和理论会更加成熟和更切合实际，从而使工程建设对防洪和水利管理的负面影响得到进一步减轻和降低，补救措施更加合理和符合实际，为我国的经济建设提供强大的技术支撑。
更多关于工程/服务/采购类的标书代写制作，提升中标率，您可以点击底部官网客服免费咨询：https://bid.lcyff.com/#/?source=bdzd
Ⅳ 地下水流基本问题的计算原理
一、地下水向河渠的运动
（一）河渠间地下水的稳定运动
1.潜水的稳定运动
河渠间潜水的运动由于受大气降水入渗补给或蒸发消耗的影响，应该属于非稳定运动。但为了简化计算，当入渗在时间和空间分布上都较均匀时，可以把潜水运动看作稳定运动。
假设条件：
（1）含水层为均质各向同性，隔水底板水平；
（2）河渠间潜水有垂向均匀入渗补给或蒸发消耗，设其强度W为常数；
（3）河渠基本上平行展布，潜水流可视为一维、渐变流并趋于稳定。
基于上述假设条件，取垂直于河渠的单位宽度进行研究并按图4-24取坐标，根据式（4-62），可写出上述问题的数学模型为：
图4-24 河渠间潜水的运动
续表