模型的特色及改进

各种模型的特点与区别~

栅格就是一个规则的阵列（matrix），其中各个像元（pix）互不影响；而矢量图是由一些个坐标和由这些坐标组成的线、面、体，他们之间有着密切的关系。
像.bmp图像就是最典型的栅格图形，.jpeg等也属于栅格图形。CAD图形就是矢量图。

层次模型特点：用一颗“有向树”的数据结构来表示表示各类实体以及实体间的联系，树中每一个节点代表一个记录类型，树状结构表示实体型之间的联系。
网状模型特点：用网络结构表示实体类型及其实体之间联系的模型。
层次模型优缺点：结构清晰，结点间联系简单，只要知道每个结点的双亲结点，就可以知道整个模型结构。但是查询子女结点必须通过双亲结点，因为层次模型对任一结点的所有子树都规定了先后次序，这一限制隐含了对数据库存取路径的控制。
网状模型优缺点：能够更为直接地描述现实客观世界。可表示实体间的多种复杂联系。但是结构比较复杂，其数据定义语言(DDL)、数据操作语言(DML)复杂，用户不容易使用。而且应用环境越大，数据库的结构就变得越复杂，不利于最终用户掌握。


扩展资料
网状模型取消了层次模型的不能表示非数状结构的限制，两个或两个以上的结点都可以有多个双亲结点，则此时有向树变成了有向图，该有向图描述了网状模型。
层次模型的数据操纵主要有查询、插入、删除和更新。进行插入、删除、更新操作时要满足层次模型的完整性约束条件。进行插入数据时，如果没有相应的双亲结点值就不能插入它的子结点值。进行删除数据时，如果删除双亲结点值，则相应的子结点值也被同时删除。
参考资料来源：百度百科-层次模型
参考资料来源：百度百科-网状模型

一、混合井孔的模拟

为了简化成井工艺，增加出水量和降低成本，本区地下水开采多采用混合井（图5-2）。

混合井在地下水流模型中的刻画，目前大都采用由美国地调局1988年推出，并在国际上广泛应用的三维有限差分地下水流模型MODFLOW软件处理。

图5-2 混合抽水井

MODFLOW建议：多层井的流量必须以某种方式人为地分配给每一单层，……把井流量按每层的导水系数大小分配，即Q_i/Q_w=T_i/∑T（1988，1996，2000年）。其中Q_i和Q_w分别为第i层流量和总（井口）流量，T_i和∑T分别为第i层导水系数和总导水系数。

为了便于讨论，又不失其一般，我们以贯穿两个含水层的混合井孔（图5-2）为例进行讨论。如此，上式可表述为Q₁/Q₂=T₁/T₂。

MODFLOW对混合井的这种处理方法没有给出理论上任何的分析、说明，缺乏理论依据，应用中也与实际不符（陈崇希等，1998；Chen Chong-xi和Jiao J J，1999）。这是因为：

1）含水层导水系数对井孔流量的影响，不会如此简单。例如，混合井附近岩性（渗透系数）发生变化，甚至混合井打在岩性透镜体上（图5-3），怎样影响各层流量的分配?含水层厚度发生变化又如何改变流量的分配?MODFLOW无法解答这些经常遇到的实际问题。

图5-3 混合井管贯穿岩性透镜体

2）含水层的参数影响混合井流量的分配，导水系数只是其中一个因素，含水层的弹性给水度（储水系数）也应该起作用。

3）井管的流量分配不仅与含水层的水文地质参数分布有关，还与外边界条件、井径（有效井径）、水泵吸水管的位置及其他抽水井的干扰等因素有关，MODFLOW仅依导水系数对混合井的流量进行预分配，未能考虑其他因素的作用。

4）含水层的导水系数一般不随时间变化，即导水系数比T₁/T₂不会变化，但流量比Q₁/Q₂却不是一个常量。

5）按MODFLOW的方法，一个混合抽水井的流量Q₁/Q₂始终是个常量。然而在模拟过程中，特别是预测时，周围随时可加入或关闭混合的或非混合的抽水井，在这种井群干扰下，原混合井的流量比还会保持常量吗?显然是不可能的。

6）从另一角度分析，混合观测孔是混合抽水井的特殊情况（Q_w=0），对于两层混合的观测孔，其孔中水位（混合水位）必介于两含水层水位之间，即混合观测孔对于其中一含水层（例如1含水层）起抽水作用（Q₁＞0），对于另一含水层（2含水层）起注水作用（Q₂＜0）。如此，Q₁/Q₂＜0。而两含水层的导水系数的比值肯定是正值，即T₁/T₂＞0。如此，两个比值怎能相等?有人认为：MODFLOW并没有说上式可用于混合观测孔。对此，很容易论证：当混合抽水井的抽水量足够小，足以保持混合井中的水位介于上下两含水层的水头之间，上述论证一样成立。

7）按各层导水系数T的比例来预先人为划分各层的流量。这是缺乏理论依据的，因为这种做法要求各分层的有效井径相等和井壁处的水力坡度上下处处相等，这两个条件不可能人为控制，预先也不得而知。

显然，用导水系数的比值预先给定各分层的流量是不妥的。实质上，这种方法不是模拟，而是“处理”，一种未考虑机理的“处理”。

“防止地下水模拟失真，提高仿真性”是水文地质模拟工作者的核心任务，而地下水流系统中普遍存在的混合井，又是当前国内外模拟失真的主要问题之一，应当引起我们的足够重视。

研究区的混合井，既有抽水井，也有观测孔。三维流场中，即使属于均质含水层（无需多层含水系统），常规（理论上非点状滤管）抽水井和观测孔都属于混合井孔，因为滤管中不同深度处的水头是不相等的，因此滤管中的水要发生垂直流动，即使在不抽水的观测孔中也一样。这就是地下水流对混合井孔响应的本质所在。

本研究采用“渗流-管流耦合模型”（陈崇希等，1992，1996）来刻画混合井孔，以“渗流”刻画地下水的运动，“管流”刻画井孔中的水流，以解决混合井孔的模拟问题，大大地提高了模型的仿真性。“渗流-管流耦合模型”的基本思路如下（图5-4）：

图5-4 计算框图

1）将产生管流的混合井视为渗透系数很大的“圆柱形透镜体”，含水层通过这个具有很大导水性的圆柱体（井筒）强烈地交换水量，那么混合抽水问题就可以视为一个特殊的“越流系统”。

2）求越流系统“圆柱形透镜体”的“渗透系数”。需要强调的是，“圆柱形透镜体”的“渗透系数”不是常量，而是随流态的变化而变化。

具体的计算方法如下：

由流体力学知识可知，当管流呈层流状态时，其水头损失可依Darcy-Weisbach方程计算：

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：ΔH为水头损失；f为摩擦系数；l为管长，m；d为管内直径，m；u为管内平均流速，m/d；g为重力加速度，m/d²。

当管流为层流时

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：Re为雷诺数。

河西走廊疏勒河流域地下水资源合理开发利用调查评价

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：ν为流体的运动黏度；μ为流体的动力黏度；ρ为流体密度。

将式（5-2）至（5-4）代入式（5-1），得

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：J为水力坡度；γ为流体重度。

将上式写成渗透流速的形式，对于管流，其空隙率n=1，则渗透流速v=nu=u，故有

河西走廊疏勒河流域地下水资源合理开发利用调查评价

将此式与达西定律v=KJ对比，可得出层流状态下管流的等效渗透系数K_L的表达式

河西走廊疏勒河流域地下水资源合理开发利用调查评价

此式由陈崇希（1966）和J.Bear（1972）分别获得。

根据流体力学知识可知，当雷诺数大于100000时，摩擦系数f与雷诺数Re无关，而取决于管内壁的相对粗糙度（），这时的水头损失ΔH与流速的平方u²成正比。在3000＜Re＜100000的范围内，其中存在f=，即ΔH∝u^1.7⁵的区段。在紊流条件下，除了上述两个区段外，还存在两个过渡区，共分为4个区段。可见，紊流是个比较复杂的问题。为了解决此问题，这里提出等效渗透系数K_N的概念。

当管流呈紊流状态时，式（5-1）可改写为

河西走廊疏勒河流域地下水资源合理开发利用调查评价

如果定义紊流状态下管流的等效渗透系数

河西走廊疏勒河流域地下水资源合理开发利用调查评价

则

河西走廊疏勒河流域地下水资源合理开发利用调查评价

该式也具有达西定律的形式。

本来，紊流的运动规律有别于达西定律，而且不同流态区具有不同的形式。引入等效渗透系数后，将5个流态分区（1个层流区和4个紊流区）的运动规律统一为达西定律形式，且与地下水渗流定律一致。即

河西走廊疏勒河流域地下水资源合理开发利用调查评价

其中K_e为渗流-管流耦合模型的等效渗透系数。

河西走廊疏勒河流域地下水资源合理开发利用调查评价

需注意的是，紊流条件下的等效渗透系数K_N与流速v及摩擦系数f有关，它是个随雷诺数而变化的量。另外，雷诺数Re的确定与流速v有关，而流速v又依赖于摩擦系数f，摩擦系数f的确定反过来又取决于雷诺数Re，因此必须采用迭代法来确定三者，进而确定等效渗透系数K_e。具体计算流程如图5-4。

二、混合观测孔水位的模拟

由于种种原因，几乎没有无混合观测孔的地区。什么是混合观测孔，混合观测孔中的水位（混合水位）如何形成，如何根据混合水位来确定初始水头的分布，如何利用混合水位求取各分层水文地质参数，等等，是“防止模拟失真，提高仿真性”（陈崇希，2003）的重要问题之一。然而MODFLOW（McDonald等，1988）对混合观测孔没有做任何分析、讨论，也没有混合观测孔的模块。简单地放弃混合观测孔的信息是不对的，因为“混合观测孔”只是孔口的流量为零，混合观测孔内却存在“抽水”与“注水”。

D.Sokol（1963）采用Thiem“影响半径”稳定井流模型证明混合观测孔的水位等于以导水系数为权重的各层水头的代数平均值。该论文存在如下几个主要问题：①Thiem“影响半径”模型是不能形成稳定流的（陈崇希，1966，1975，1983）；②Sokol隐含着假定，尽管两含水层的导水系数不等，流量也不同，但取其“影响半径”相等；③Sokol还隐含着假定，井管的阻力可以忽略不计，即井管处处水头相等。在这些假定条件下得出的结论是不可信的。

Hantush（1961）和Бочевер等（前苏联）（1961）在对承压含水层非完整井流的研究中，都认为观测孔中的水头降深s（r，l′，d′，t）反映该孔滤水管中各点降深s_p（r，z，t）的平均降深，即

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：l′，d′是观测孔滤管顶点、底点的标高。本书称之为Hantush-Бочевер方程。

对于潜水三维流的研究，Neuman（1972）也认为，观测孔中的水头降深可视为滤管内各点降深的平均值，即

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：r是观测孔至抽水井的距离；z₁、z₂是观测孔滤管顶点、底点的标高。

上述诸学者的见解与确定观测孔中水位的方法在水文地质界一直沿用至今，40多年来未见到异议。然而，上述观测孔中水位的确定方法，缺乏对形成机理的基本分析。简单分析之，三维流中的垂直观测孔，由于观测孔滤管中的水头不等，会发生垂向流动，而水的流动则会导致水头的变化，并引起井孔周围地下水的运动与水头再分布。因此，与测压计式观测孔不同，一般观测孔实质上是一部分滤管进水（抽水）而另一部分滤管出水（注水），即观测孔并不简单地反映含水层地下水的水头，而是它兼有“抽水”与“注水”作用的井孔，只是孔口的流量为零，井管内的“抽水”量与“注水”量绝对值相等而已（如果忽略微小的井筒储存量的变化）。观测孔非“观测”孔也！

显然，上述Hantush、Бочевер等以及Neuman提出的“积分平均水位”都未涉及水流的机理，是缺乏物理基础的纯数学方法。作者曾用“渗流-管流耦合模型”来模拟混合观测孔中的水位（混合水位）。算例的条件是：

一个均质、等厚、水平无限延伸的承压含水层，于顶面处有一微小的半球状井定流量抽水，抽水延续时间t=10.02d，在不同径距r处有不同孔径的完整观测孔（此即混合观测孔），模拟这些观测孔中的水头分布和流量变化。含水层的厚度M=100m，渗透系数K=100m/d，单位弹性给水度μ_s=0.00001m^-1，抽水流量Q=36000m³/d。模拟结果示于图5-5和图5-6。由图5-5可见，孔径≤0.1m的观测孔中，不同深度处的水头值明显存在差异，与Hantush-Бочевер的积分平均相差也很大。

图5-5 不同径距r观测孔中水头降深分布及Hantush-Бочевер公式计算值s_ga

图5-6 不同径距观测孔流量分布

由于抽水井（点汇）位于承压含水层顶面处，因此所有径距r处的顶部处的水头降深s大于底部，即底部的水头H大于顶部的水头，于是观测孔中的水从下而上流动。基于水流连续性原理，底部观测孔要吸水（抽水）而顶部观测孔要排水（注水），图5-6表示观测孔中的流量分布。图5-6a表明：径距21.19m，直径0.20m的观测孔中，大约在z=35m以下为“抽水”，35m以上为“注水”；观测孔中垂向最大流量超过500m³/d。这就是我们所说的，观测孔非“观测”孔也！

基于上述分析，本项研究采用“渗流-管流耦合模型”来模拟混合观测孔中的水位，使混合水位成为有用的信息，即用于初始水头分布和拟合求参。

三、泉流量动态的模拟

泉是地下水转化为地表水的主要形式之一，不同类型泉的出现及其流量动态，是水文地质工作者要重点分析的问题，因为这些数据为认识水文地质条件提供了极为重要的信息。泉是数值模型中主要的模拟要素；泉流量的预测是地下水资源评价、管理不可缺省的内容。

美国地质调查局（2001）建议采用MODFLOW软件的Drain模块迭代计算泉流量，实质上将泉作为第二类边界来刻画，即先不考虑泉的存在求解水头分布，再用泉所在格点的地下水水头与泉口标高之差乘以某比例系数（该系数缺乏物理意义）计算出泉流量，再以此流量置于该格点（第二类边界）重新求解水头分布，再次重新计算出新的泉流量，如此迭代直至收敛为止。当地下水水位低于泉口标高时，流量为零；当地下水水位高于泉口标高时，泉流量与水头差成正比，其比例系数必须通过流量的拟合来反算。

我国以往数值模型中关于泉流的刻画，采用的像抽水井一样输入（给定）流量。如此处理，预测怎么办?泉流量是未知的。近年来，我国也有类似上述MODFLOW建议的迭代算法。例如有的研究报告将鄂尔多斯大向斜的东侧山西省柳林泉（实为单斜自流斜地上升泉）刻画为平面二维流模型，进而用承压完整井的公式计算泉的流量，即Q=qs。这里存在几个问题：①柳林泉是上升泉，泉附近的钻孔愈深其地下水头也愈高，这是单斜自流斜地上升泉的特点，这种条件具有明显的三维流特征；②用承压完整井的公式计算泉的流量，特别是柳林泉的流量，不妥；③比例系数q在不同时间（不同条件）并不是常量。

“泉，必需将其放在水文地质体中去研究其成因类型，才能正确模拟仿真”。本模型提出的方法不同于MODFLOW等，是将泉口标高作为第一类边界条件（当泉流存在时），并将泉口格点与同层周边格点及下格点分别依达西定律（如果地下水流属于线性流，例如陈崇希等（1995），Chen Chongxi等（1996），陈崇希等（2003）；或采用非达西定律，例如陈崇希（1995），成建梅等（1998），ChengJianmei等（1998））和水均衡原理建立关系，整个计算方法建立在流动机理之上，因此由模型运转直接产生泉流量，无需通过迭代求解（且有的情况下迭代并不收敛）。正因为此，在模型识别阶段才有可能将泉流量作为拟合对象。这一点十分重要。

具体地说（图5-7），根据达西定律，泉流量Q为水平向的流量与垂向上的流量之和，如果以水平向的流量为主，则为下降泉，否则为上升泉。具体可通过下式计算：

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：Q为泉流量；K_zsp为泉流区域的垂直渗透系数；H_b为泉口下层结点的水头值；Z_sp为泉出露标高；D_zsp为泉流过程中的渗流长度（含水层厚度）；S为泉口结点控制面积；C_A为泉颈修正系数；C_s为泉径修正系数；T_ij，T_ik为流段ij和ik的平均导水系数；，，，为线段的长度；h_i，h_j，h_k为结点i，j，k的水头值；e为结点i周围的均衡单元，如ipbq四边形。

图5-7 泉的流量计算示意图

在以往的地下水资源评价数值模拟中，没有将泉的流量作为拟合要素。实际上，象MODFLOW等的做法也不可能将泉的流量作为拟合要素，因为他们在模型识别其间要求解一个缺乏物理意义的系数。

本模型选择了3个典型点泉和22条线泉（泉沟），模拟其流量动态，并将其实测流量作为一个重要的拟合对象用以求取含水系统的水文地质参数和预测泉流量的动态。如此将会提高模型的仿真性。

四、大气降雨、地表水等入渗补给潜水滞后性的刻画

降水、地表水、渠系等入渗补给潜水是地下水的重要补给来源，若处理不当，则会导致模型失真。因此，正确的刻画入渗补给的滞后性对于模型仿真性的提高具有重要的意义。

目前地下水资源评价中降水、地表水、渠系等入渗补给量的计算主要采用下列两种方法：

1）未考虑滞后的入渗系数法。该法对于小埋深区域大体上可以用来刻画入渗对潜水的补给，但有些地区潜水埋深较大，若不考虑降雨、地表水等入渗补给的滞后特征，则会导致模型失真。

2）平移滞后入渗系数法。有的研究报告把潜水位出现高峰时间与降雨高峰出现时间的时间差，作为统一的滞后时间。例如若该时间差为3个月，则1月份的降雨量都在4月份入渗补给，2月份的降雨量则在5月份入渗补给；若3月份无降雨，则6月份无降雨入渗补给。如此处理滞后补给也存在明显的不合理性。因为潜水位高峰出现的时间，是所有影响地下水位动态诸因素作用的总和，除降雨补给外，还有地下水蒸发，灌溉入渗，越流和地下水开采等因素作用下产生地下水不平衡流动的结果，特别是地下水开采这个人为因素在本区已经成为地下水动态的主要控制因素，它已经完全掩盖了地下水的天然动态，因此用该法研究降雨入渗补给的滞后性是不恰当的。即使这一地区完全处于未开采状态，甚至忽略除降雨以外的其他因素对地下水动态的影响，仅就降雨入渗补给单一因素对潜水位动态的影响而言，该法也存在明显的问题。降雨过程是间断的，不连续的，但它对潜水的补给则是连续的，这是因为地下水流动，特别是非饱和流动，其介质具明显的储存性而有调节功能，从而使补给具有滞后性及连续性。一年内虽然降雨仅有几十天，其余时间均无降雨，但在潜水位深埋区，补给却天天作用着（浅埋区由于存在强蒸发作用而使情况复杂化了），这种连续补给作用是不可能用一个时间差的平移来刻画其滞后性的。

上述“入渗系数法”在潜水小埋深区有一定的适用性。但由于本区部分地段潜水埋深较大，如昌马洪积扇顶部潜水水位最大埋深达290m上下。一个月的降雨量乘以入渗系数的水量不可能在当月内全部补给其下的潜水含水层。因此，不考虑降雨补给的滞后特征，则会导致模型的失真。

本研究采用“降雨补给滞后权系数法”（陈崇希等，1991，1998）来刻画降雨入渗补给。该方法既能基本反映客观实际情况，又很实用。

根据地下水非饱和流动理论，一次降雨过程对潜水补给量的分布曲线如图5-8所示。它是一单峰曲线。曲线的形态取决于潜水位的埋深和包气带的岩性，当然还与包气带的初始含水率的分布有关。我们将后面的因素固定起来，专门讨论前面因素的作用。

图5-8 一次降雨对补给量的分布曲线图

我们以包气带地下水非饱和流动理论为基础，考察不同水位埋深和岩性条件下，某时段（例如月、旬等）降雨对该时段及其后各时段补给强度的分布特征。为了使问题的讨论能与地下水资源评价数值法相匹配，将连续的特征曲线改为离散点的形式（因为数值法在时间上是离散的）。不失一般性，我们将时段取为1个月。

对于潜水位埋深较大和（或）包气带渗透系数较小的条件，入渗强度分布是一单峰函数，当潜水位埋深很小时，可能是单调减函数（图5-9）。若在0月有某降雨量RA（mm）作用，则该月及后各月降雨入渗补给量之和∑RE_k（mm）就是该时段降雨量形成的全部补给量。该值∑RE_k与0月的降雨量RA之比就是（总）入渗补给系数α。各月份的入渗量与0月降雨量RA之比，则是0月降雨在各月的入渗补给系数α_k。即

河西走廊疏勒河流域地下水资源合理开发利用调查评价

图5-9 降雨对不同埋深（岩性）条件的入渗补给历时分布曲线图

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：RA为0月的降雨量，mm；RE_k为0月的降雨在k月的入渗补给量，mm；RE为（总）补给量，mm；α_k为0月降雨在k月的入渗补给系数；α为（总）入渗补给系数。

如果令图5-9上0月的降雨量为100mm/月，则该图纵坐标转变为以百分数表示的月入渗系数α_k（k=0，1，2，3，…）。

月入渗补给系数组已能很好地反映降雨入渗补给滞后性的特征，为了使刻画滞后性的特征数组更通用，我们定义：

河西走廊疏勒河流域地下水资源合理开发利用调查评价

因此

河西走廊疏勒河流域地下水资源合理开发利用调查评价

该式表明ω_k满足权的定义，所以我们称ω_k（k=0，1，2，…）为滞后补给权系数（组）。

引入滞后补给权系数（组）的优点是，不同入渗系数的诸分区（不同区，同月的入渗补给量不等），有可能取用同一滞后补给权系数来计算各月的补给量，也就是说，滞后补给权系数更具通用性。于是我们的任务转到寻找一个能满足上述条件的特征函数，即：

1）它们是单峰函数或单调减函数，并随峰值（或最大值）出现时间的后延，其峰值（或最大值）减小；

2）各权系数之和为1；

3）适应性较强，能较好地与不同潜水位埋深，不同岩性条件下滞后入渗补给量相拟合。

经反复分析研究，我们确定采用下面的离散函数来表征滞后入渗补给权系数的分布。

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：ω（j-i）为i月降雨，j月补给的权系数；h为潜水位埋深，m；B为与包气带岩性等因素有关的经验系数，m；可称为埋深区间系数。该系数是包气带渗透系数的增函数。

该式中h是已知的，仅有一个与岩性等因素有关的经验系数B是待求的。这个系数可以通过模型识别来确定，这样就可以利用上式计算出滞后入渗补给权系数。各区（总）入渗系数也是通过模型识别来确定的，由此可计算出在某月i降雨量RA（i）作用下，该月及其后各月j的补给量RE（j-i），即

河西走廊疏勒河流域地下水资源合理开发利用调查评价

至此，我们分析了某月降雨量对该月及其后各月入渗补给的贡献。对于实际问题——多月降雨对某月j的补给量问题，为了简化计算，我们近似用该月及其以前各月i降雨量对潜水补给贡献的代数和来计算该月j获得的总入渗补给量，即

河西走廊疏勒河流域地下水资源合理开发利用调查评价

式中：RES（j）为j月及其以前各月的降雨量在j月的总入渗补给量，mm；RE（j-i）为j月降雨对j（j＞i）月的入渗补给量，mm；α为入渗补给系数；ω（j＞i）为i月降雨对j（j＞i）月的滞后补给权系数；RA（i）为i月的降雨量，mm。

依式（5-17），获得的降雨入渗补给滞后权系数曲线共75条，它对应于75个潜水位埋深区间，每个区间埋深差Bm（即埋深区间系数）即每增加埋深Bm，换一条权系数曲线的序号）。图5-10中仅绘出12条典型的权系数曲线，而各两条曲线之间，尚有5～7条曲线图中并未绘出。

降雨补给滞后权系数法的提出和使用，不仅大大提高了模拟的仿真性，而且使用十分方便。地表水入渗补给滞后性的机理与降雨入渗的类似，也采用上述方法刻画。

图5-10 降雨入渗补给滞后权系数曲线

五、参数的分区

本模型含水介质厚度大，垂向上共分6个模拟层，即便是同一岩性（如亚粘土）而处于不同埋藏深度，由于其固结程度不同而导致其水文地质参数不相等。这种情况下如何给定初始参数估计值的空间分布又是一个难题。若按常规惯用的方法，则同一岩性根据不同层位（不同埋藏深度）要给出不同参数分区，如此的参数分区数过多，给调参带来很大的工作量。

对此问题，本模型采用如下方法，处于不同埋藏深度同一岩性的K_h值作为同一分区，但其参数值随埋深的负指数衰减。这样，每个参数分区只要1个零埋深的水平渗透系数K_h0和1个衰减系数γ刻画，即K_h=K_h0×e^-γh（h为埋深值）。其他参数，如K_z值、μ_d值和μ_s值的空间分布也按同样方法给定。如此处理，就大大地减少了模型参数的数量，相应地也减少了模型调参的工作量。

六、自流井的模拟

自流井普遍地存在于我国西北内陆盆地，冲洪积扇的前缘洼地，大型岩溶上升泉的周围以及平原区的深部承压含水层等区段。我国至今还存在大面积的地下水自流区（尽管由于大量开采地下水而缩小了其自流面积）。因此模型中如何刻画自流井是一个重要的水文地质问题。然而，美国地调局推出的MODFLOW软件尚无模拟自流井流量的模块，这似乎不是疏忽，而是存在一定的难度。常规的地下水流模型已经解决不了这类问题。至今，对于数值模型中的自流井，仍是无可奈何地采用人为地给定其自流量。然而，自流井与抽水井不同，后者人们一般可以人为给定抽水流量，而井中水位是作为其响应；当井中水位降至井底时，人们就不能自由地抽取了（尽管如此，大多含有抽水井的数值模型，未能客观地将“井中水位作为抽水量的约束”来认识）。然而，自流井就不能人为地给定自流量了，因为自流量是自流井对周围环境（气象、水文、地下水的开采等因素）的响应，自流井的流量是属于预测范畴，是输出信息，而不是输入信息。

本项目采用“渗流-管流耦合模型”模拟自流井，只要在自流井的井口处设置水头已知边界，即水头等于井口标高，即可模拟出自流井的流量。该法已在“陕西渭北东部岩溶水资源数值模拟研究”中获得成功的应用（陈崇希等，2003）。

七、初始水头的形成

初始水头的分布是地下水不稳定流数值模拟不可缺少的条件。通常利用一定数量观测孔的水位通过插值获得各结点（格点）的初始水头值。然而，当研究区范围较大，观测孔较少，观测孔所在的层位不一，而且含有混合孔时，难以用插值的办法获得每一模拟层的初始水头分布。

本项目采用“参数-初始水头迭代法”（陈崇希、林敏等，1991；陈崇希、裴顺平，2001），来确定各层的初始水头分布。具体处理办法是：取模拟的初始时间为2003年7月，用水文地质参数初始估计值从2001年7月开始模拟，至2003年7月将模拟所得的该月水头值与实测值拟合，粗略确定模拟时段的初始水头（2003年7月）；再以此初始水头按常规方法求参，如此迭代求解。当拟合达到一定程度之后，必须全面地检验模型是否符合基本规律，否则要再次重新调整，直至满足要求。上述方法虽然工作量很大，但效果良好。

---

沭阳县13444649822： 关系模型具有哪些特点? ？
勇岚百为： 关系模型:用二维表的形式表示实体和实体间联系的数据模型.关系模型的特点:关系模型与非关系模型不同,它是建立在严格的数学概念基础上的.2.关系模型的概念...

沭阳县13444649822： 分布系数法的优缺点 - ？
勇岚百为： 增长系数法的特点:优点:(1)结构简单,实用的比较多,不需要交通小区之间的距离和时间.(2)可以适用于小时交通量或日交通量等的预测,也可以获得各种交通目的的od交通(3)对于变化较小的od表预测非常有效.(4)预测铁路车站间的...

沭阳县13444649822： 数学建模竞赛论文基本步骤是什么? - ？
勇岚百为： 基本步骤: 答卷的基本内容 1、问题的叙述,背景的分析等 . 2、模型的假设,符号说明(列表). 3、 模型的建立:问题分析,引用的数学命题,公式推导,模型Ⅰ,模型Ⅱ 等.. 4、 模型的求解:计算方法设计或选择,计算步骤(框图),所采用的软件名称等. 5、模型的结果:误差分析,模型检验. 6、模型评价:特色,优缺点,改进方法,推广. 7、 参考文献. 8、 附录:图表、程序等.

沭阳县13444649822： 数学建模 范文 - ？
勇岚百为： 数学建模论文写作: 一、写好数模答卷的重要性 评定参赛队的成绩好坏、高低,获奖级别,数模答卷,是唯一依据.二、答卷的基本内容,需要重视的问题 1.评阅原则假设的合理性,建模的创造性,结果的合理性,表述的清晰程度.2.答卷的...

沭阳县13444649822： 房屋建筑模型主要有哪些作用? ？
勇岚百为： 作用一、更好的表现设计效果 (1)沙盘模型是向观者展示其设计特色的一种很好的表达方式; (2)沙盘模型是设计师与业主之间进行交流的重要手段. 作用二、完善相关设计构思 (1)沙盘模型是进一步完善和优化设计的过程; (2)设计人员亲自动手制作模型,是从二维平面到三维立体形象的体验. (3)通过亲身感受与参与制作,可以进一步激发设计师的灵感,发现设计中存在的问题,并进行改进与优化,使设计方案达到理想的状态.

沭阳县13444649822： 如何写一篇好的数学建模论文 - ？
勇岚百为： 如何撰写数学建模论文 兼谈数学建模竞赛答卷要求 当我们完成一个数学建模的全过程后,就应该把所作的工作进行小结,写成论文.撰写数学建模论文和参加大学生数学建模时完成答卷,在许多方面是类似的.事实上数学建模竞赛也包含了学生写...

沭阳县13444649822： 细胞膜的结构模型中的流动镶嵌模型特性是什么? - ？
勇岚百为： 1、细胞膜的流动镶嵌模型:从流动和镶嵌两个方面学习和记忆. 2、流动:组成细胞膜的分子会运动,表现出细胞膜的结构特点是一定的流动性.磷脂双分子层是细胞膜的基本骨架结构,磷脂分子会运动,蛋白质分子的分布不均匀,蛋白质分子会运动. 3、镶嵌:蛋白质分子在磷脂双分子层中的分布是不均匀的,有的镶,有的嵌,有的横跨,还有贯穿等. 4、因为分子的运动实现了物质的跨膜运输,实现了变形虫的变形运动等细胞的生命现象.

沭阳县13444649822： 数学建模步骤 - ？
勇岚百为： 摘要在整篇论文评阅中占有重要权重,务必认真书写(篇幅不能超过一页).全国评阅时将首先根据摘要和论文整体结构及概貌对论文优劣进行初步筛选.摘要写得不好,论点不明,条理不清,评委不再阅读正文,论文即遭被淘汰. 摘要是全文...