实验六 地下水流动系统数值模拟演示(计算机)

一、实验目的

运用计算机模拟技术，通过对不同条件下，地下水流动特征的演示，了解并熟悉地形盆地及流动系统的概念。

二、实验内容

1.模拟一个简单地形盆地地下水流动系统和复杂地形盆地地下水流动系统。

2.按实验中的提问，独立设计一个流动模式，进行流动系统模拟演示。

三、模拟程序简介

实验采用的是荷兰自由大学地球科学院的教学演示程序，程序名为“FLOWNETD”，用于地下水稳定二维流模拟。

四、实验步骤

1.简单地形盆地地下水流动系统演示

(1)均质介质场

模型条件:左、右、下边界模拟隔水层，上边界为单斜的地形坡度，模拟从左到右高程逐渐降低的势边界，中间模拟均质的砂层。渗透系数为Kx=Ky=1m/d，模型(1)。

问题:在上述模拟条件下，盆地的流线是什么样的?请用铅笔在图Ⅰ6-1a中试画。

*** 作:①开机进入“FLOWNETD”程序， *** 作到屏幕显示第3页②键入(1)模式程序名:EXAM1.1→F1→PgDn③屏幕显示第4页，观察上边界水头分布→PgDn④屏幕显示第5页，观察介质场的分布与数值→PgDn屏幕显示第6页，继续→PgDn，计算完毕，屏幕显示最终流网结果(第7页)⑤按键盘右侧的“+”键，观察流线上的等时点的分布与运动情况。

注意:计算完成后，按PgUp键返回第6页，认真阅读屏幕第6页的内容。通过改变对应的描述，可以显示不同的计算结果图，如均衡计算、加密流线或等势线、改变显示比例、显示介质场等。

思考题:流线分布有何特征?等时点分布(疏密)说明什么?

(2)非均质介质场(供选做)

在(1)模型条件中，置入一个渗透系数较高或较低的岩层或透镜体，构成非均质介质。图Ⅰ6-1b、c、d中K1的位置与数值可以任意给定。

问题:图Ⅰ6-1b、c、d模式的流网与图Ⅰ6-1a有何不同，请用铅笔在图Ⅰ6-1b、c、d上示意画出流慎碧线。

*** 作:重复(1)的演示 *** 作团孝锋步骤至④，按照计算机屏幕下方的提示，确定K1的数据，并置入到模型中，形成非均质介质场(见图Ⅰ6-1b、c、d)继续(1)演示中的步骤⑤，直至屏幕显示计算结果。

2.复杂地形盆地地下水流动系统演示

(1)均质介质场模拟演示

模型条件:左、右、下边界为隔水层，上边界为起伏的地形，中间为均质的砂层，如图Ⅰ6-2a所示。

注意:演示前阅读并思考图Ⅰ6-2a所示的流动系统特征。

*** 作:①开机 *** 作至屏幕显塌晌示第3页②键入图Ⅰ6-2a模式程序名:Toth→2.1→F1③重复简单地形盆地地下水流动系统演示(1)中的演示步骤③和步骤⑤。屏幕显示复杂的地下水流动模式，观察流线、等势线与等时点的分布特征。

(2)非均质介质场模拟对比(供选做)

在图Ⅰ6-2a模式中，任意加入一个透镜体，如图Ⅰ6-2c所示，重复演示并观察流动特征。

演示 *** 作:修改参数同简单地形盆地地下水流动系统中非均质介质场中的模式， *** 作相同，其他 *** 作仍与简单地形盆地地下水流动系统中均质介质模拟演示相同。

思考题:对照(1)与(2)的模拟演示结果，如图Ⅰ6-2a与图Ⅰ6-2c所示，指出它们的不同之处。

五、实验成果

1.提交简单地形盆地地下水流动演示结果图系(见图Ⅰ6-1)，并对模拟结果进行简要分析。

2.提交一张复杂地形盆地地下水流动系统演示结果图(可以利用图Ⅰ6-2b或图6-2c，根据演示结果草绘)，并对模拟结果进行简要分析。

图Ⅰ6-1a 简单地形盆地地下水流动系统图(K=C，C为常数)

图Ⅰ6-1b 简单地形盆地地下水流动系统图(K1=5K2)(问题:流线密度如何分布?画出模拟流网)

图Ⅰ6-1c 简单地形盆地地下水流动系统图(K1=5K2)(问题:流线密度如何分布?画出模拟流网)

图Ⅰ6-1d 简单地形盆地地下水流动系统图(K1=5K2)(问题:流线密度如何分布?画出模拟流网)

图Ⅰ6-2a 复杂地形盆地地下水流动系统图(K=C，C为常数)

图Ⅰ6-2b 复杂地形盆地地下水流动系统图(K1=5K2)(问题:流动系统数量，系统级次如何?画出模拟流网)

图Ⅰ6-2c 复杂地形盆地地下水流动系统图(K1=5K2)(问题:流动系统数量，系统级次如何?画出模拟流网)

图Ⅰ6-2d复杂地形盆地地下水流动系统图( (问题:流动系统数量，系统级次如何?画出模拟流网)

不太明白你的意思，你看下面的代码可行吗？

#include <stdio.h>

#include <string.h>

void main(){

char buffer[256]

while (gets(buffer)){

if (0==strcmp("case",buffer) ||

0==strcmp("char",buffer) ||

0==strcmp("const",buffer) ||

0==strcmp("continue",buffer) 迟悉冲档

){

printf("OK\n")

}

else {

printf("NO\n")

}

if (0==strcmp("quit",buffer))break  // 码判乎 退出

}

}

是的，

1.任何DFA都不能蚂大识别e(空)符号串

2.一个DFA,只能包含唯一的开始状誉物锋态3.DFA识别的符号串集合,可以庆晌是有限的

4.一个DFA,所有的映射必须是单值...

---