描述孔隙含水层的基本水文地质参数是（　）。

描述孔隙含水层的基本水文地质参数是（ ）。

A、孔隙度、给水度、渗透系数

B、空隙率、给水度、弹性系数

C、有效孔隙度、渗透系数、实际流速

D、有效孔隙度、弹性系数、水力传导系数

参考答案

【正确答案：A】

描述孔隙含水层的基本水文地质参数有：①孔隙度，指某一体积岩石（包括孔隙在内）中孔隙体积所占的比例；

②给水度，即把地下水位下降一个单位深度，从地下水位延伸到地表面的单位水平面积岩石柱体，在重力作用下释出的水的体积；

③渗透系数，又称水力传导系数。在各向同性介质中，它定义为单位水力梯度下的单位流量，表示流体通过孔隙骨架的难易程度；在各向异性介质中，渗透系数以张量形式表示。

地下水资源模型计算

数值模型模拟计算方法适用于非均质性、各向异性的复杂地下水系统，包括存在越流和具有不规则形状各类边界条件等情况。但是该方法对资料的要求比较严格，要求研究程度较高和资料较丰富。应用数值模型方法的一般程式为：①气象、水文、水文地质资料分析→②水文地质概念模型概化→数学模型建立（水动力方程和定解条件）→③选择计算程序→④模型设计→⑤模型识别和检验→⑥计算模拟。

一、水文地质概念模型建立

在对黑河流域地下水系统做了全面、深入分析的基础上，根据研究目的，对地下水系统的组成要素和相互关系作出合理的简化和假设，并且用文字、框图、平面图、剖面图等形式把系统再现出来，即为地下水系统概念模型。

（一）地下水系统空间结构概化与边界确定

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 黑河流域模型采用的有关水文地质参数

常用水文地质参数的类型

1.渗透系数K

根据达西定律,渗透系数是水力坡度等于1时的渗透流速。对于具体工程,土层的渗透系数关系到降水设计方案的选择、水位降深的大小及基坑涌水量的大小,影响到降水时间的长短及工期。渗透系数选取正确与否直接关系到降水的成败,该参数是基坑降水设计中最重要的水文地质参数之一。土层的渗透系数可由岩土工程勘察报告提供。对于勘察报告中没有提供该参数或提供的参数未经试验取得,对一些中小工程,可采用经验值对于一些重大工程,应进行水文地质补充勘察、试验,来确定水文地质参数。

影响渗透系数主要因素为渗透流体和土的颗粒大小、形状、级配以及密度。渗透流体的影响主要是粘滞度,而粘滞度又受温度影响。温度越高,粘滞度越低,渗流速度越大。

土颗粒的影响是颗粒越细,渗透性越低级配良好的土,因细小颗粒充填在大颗粒的孔隙中,减小孔隙了尺寸,从而降低渗透性。土的密度增加,孔隙减小,渗透性也会降低。

影响粘性土渗透性的主要因素为颗粒的矿物成分、形状和结构(孔隙大小和分布)。粘土颗粒的形状为扁平的,有定向排列作用,因此渗透性具有显著的各向异性性质。层状粘土水平方向的渗透性往往远大于垂直方向而黄土和黄土状土,由于垂直大孔隙发育,其中的垂直方向的渗透性大于水平方向。

2.降水影响半径R

根据裘布依理论,井点系统开始抽水后,地下水位围绕抽水井形成了降落漏斗。随着抽水时间的延长,地下水流出现相对稳定状态,降落漏斗的曲线逐渐向外扩大直至达到稳定。在距离降水井距离为R的地方,观测不到地下水位的变化,该稳定的降落漏斗的半径即为降水影响半径R。

当要求计算精度不高时,可采用经验值或经验公式计算。对计算精度要求较高的工程应采用现场抽水试验的方法确定降水影响半径。

3.给水度μ

给水度表示潜水含水层的释水能力,它表示单位面积的含水层当潜水面下降一个单位长度时,在重力作用下所能释放出的水量。给水度大,说明含水层能够释放的水量大,反之则小。

给水度大小与含水层岩性有关。松散沉积物含水层的颗粒粗、大小均一,则给水度大反之,颗粒细、大小不均,则给水度小。

在基坑降水设计计算中,给水度可采用经验值。对重要工程可采用室内实验、室外抽水试验来确定该值。

4.贮水系数S

贮水系数S(或弹性给水度μ*)是指承压含水层的测压水位下降或上升1个单位时,单位水平面积的含水层(厚度为M)释出或存储的水的体积称之为贮水系数。无量纲。

5.导水系数T

导水系数是表示含水层导水能力的大小的参数,它是渗透系数与含水层厚度的乘积。

6.导压系数α

压力传导系数是表示水压力向四周扩散、传递的速率,为导水系数与贮水系数的比值。贮水系数、导水系数可由现场抽水试验确定。

渗透系数和降水影响半径是进行稳定井流计算的主要水文参数,进行非稳定流计算则需用到贮水系数、给水度、导水系数。

水文参数是什么

水文参数是表征与岩石性质、水文气象等因素的数量指标，主要包括：

1、降水入渗系数

2、潜水蒸发强度

3、灌溉水回渗补给系数等

狭义的水文地质参数是表征含水介质水文地质性能的数量指标，主要包括：

1、含水层的渗透系数和导水系数

2、层压含水层的储水系数

3、潜水含水层的给水度

4、弱透水层的越流系数

5、含水介质的水动力弥散系数等

水文参数和狭义的水文地质参数统称为水文地质参数。