surfer河道演变分析.docVIP

Surfer在河道演变分析中的应用 1.2绘制数字高程模型图 经过前期数据处理后，就可以绘制数字高程模型图了。具体步骤如下： 步骤一，把数据文件转换成grd文件：①打开菜单“网格|数据”在open对话框中选择数据文件；②打开“网格|数据”对话框.在“数据列”中选择要进行grid的网格数据（X和Y坐标）以及格点上的值（Z列）（不用选择，因只有3列数据且它们的排列顺序已经是X，Y，Z了，如果是多列数据，则可在下拉菜单中选择所需要的列数据）。选择好X，Y，Z值后，在“插值模式”中选择一种插值方法（如需要比原始数据的网格X和Y更密的Z数据，或网格为非均匀），则在grid的过程中，Surfer会自动插值计算，生成更密网格的数据。如果只是想绘制原始数据的图，不想插值，则最好选择反距离加权插值法（Inverse Distance To A Power）或克里金法（Kriging Method）。因为这两种方法在插值点与取样点重合时，插值点的值就是样本点的值，而其他方法不能保证如此。在Output Grid File中输入将输出的文件命名，然后在“网格点几何分布”中设置网格点数，确认，画图所需要的grd文件就生成了。不过，为了便于后面对各年地形进行比较分析或冲淤分析，尽量使每个grd文件的几何分布一直，即同样的XY坐标范围和插值的网格密度。 步骤二，将河道边界白化。在Surfer中默认的插值区域为数据文件中离散点坐标x，y的最小值和最大值所围成的矩形，经过插值生成的图形边界为矩形，但在实际情况下，河道边界可能是不规则的，或者需要显示某些特定区域的形态（如潜洲）、添加图签等，这时就用到Surfer的白化（Grid Blank）功能。 白化文件[.bln]格式 [.bln]文件是以ASCII文件格式存储的用来描述白化边界及白化信息的文件，其格式如下： length，flag″Pname 1″ x1，y1 x2，y2 ... xn，yn x1，y1 x1，y1 x2，y2 ... xn，yn x1，y1 图1. 长江宜都河段三维数字高程图 图2.长江宜都河段等值线图 1.3提取等值线 数据处理基本方法同上，只是在步骤三中要选择“地图│等值线图”命令（map│contour map），导入白化后的文件，即可生成一个等值线图，在等值线图属性中“常规”选项卡中，选择“平滑等值线”，在“等级”选项卡中可更改默认等级，设置要在图中显示出的某一特定等级的等值线（如航行0m或3m等深线等）并选择标注，这样该等值线就生成了，如图2中标注为35.79的即为当地的航行基准面，而黑色较粗标注为32.79的为当地的航行基准面以下3m水深线。如果需要在AutoCAD河势图中显示该等值线，则需要将该等值线图在Surfer中输出为DXF文件，用AutoCAD打开，将要显示的等值线选中并原坐标复制到河势图中即可。 2.绘制河道冲淤分布图和计算冲淤量 2.1绘制河道冲淤分布图 按照1.2中步骤一和步骤二的方法，可以生成白化后的grd文件。选取“网格│数学”（grid│math）命令可以进行grd文件的相减操作，先在对话框中选中一个某时期河道白化后的grd文件 （A文件），确认，然后另选一个要与之对比的另一时期的河道白化文件（B文件），再将第三个文件（C文件）命名为“A-B.grd”，最后将运算公式写为“C=A-B”。由于文件运算的是Z值，那么“A-B.grd”文件即反映为B时期到A时期该分析河段各坐标点的Z值变化量。注意：必须是同样的网格信息，即平面区域和间距均相等，否则无法运算。然后选取“地图│影像图” 命令（map│image）即可绘制成河道冲淤分布图，双击该图，在属性各个选项中进行调整，使其具有较好的视觉效果，如图3即反映的是长江宜都河道在三峡蓄水后某年度的冲淤分布。 图3.长江宜都河段某年度冲淤分布图 图4. 断面bln文件的格式 2.2冲淤量的计算 Surfer可以对三维地形进行体积和面积计算，只需指定上表面和下表面即可。由于在2.1中相减生成的grd文件反映的是冲淤变化，那么利用数据计算功能“体积”还可以计算出冲淤量。具体操作如下：选取“网格│体积”（grid│volume）命令，指定上表面为该网格文件，下表面为高程为0的水平面，确认之后获得的报告文件中含有“cut volumes（正体积）”、“fill volumes（负体积）”、“net volumes（净体积）” “planar areas（平面积）”等数据指标，其中“cut volumes（正体积）”表示上表面在下表面之上的那部分实体的净体积，在这里即为泥沙的淤积量；“fill volumes（负体积）” 表示上表面在下表面之下的那部分实体的净体积，即为泥沙