水位数据处理计算程序开发和编制.doc

水位数据处理计算程序开发和编制 上海达华测绘公司 寿顺宝 摘要：针对在水运工程水深测量中遇到水位控制、水位改正的实际问题，提出了水位加权系数平均值计算方法，并应用vissual Basic 6.0编写计算程序。该方法的应用是对传统作业方法的一大改进，免去了人工绘制联合潮位曲线，提高读取水位值的准确度和分辨率。 关键词：水位 数据处理 程序 开发 编制 前言 水位改正是水深测量工作中的重要部分。为了正确表示水底形态，需要将瞬时水面测得的深度化算到深度基准面起算的深度。 水位改正的方法较多，在水运工程水深测量中，大部分采用分带法。因为分带法适用潮差线和相同潮时线基本上与岸线垂直的测区，同时要求相邻水位站之间距离要满足最大潮差不大于1m、最大潮时差不大于2m，否则须增设定点水位站或采用其他解析法。 分带方法一般可分三种：单站分带水位改正；双站分带水位改正；三角（三站）分带改正。单站水位改正、双站水位改正往往是绘制测区邻近水位站的水位曲线图进行分区、分带内插水位改正，而三角分带改正是采用模拟法，一般将测区按一定的比例缩小，使三个水位站和测区能绘制在同一幅图上，用模拟杆分别竖插在三个水位站上，用某时刻的水位读数或潮高值模拟瞬时该时刻海平面进行分区分带内插水位改正，直接用几何法量取测深点在这模拟瞬时海平面上的水位值或潮高，这就是所谓的“三角分带法”。在实际操作中这个方法由 由于受客观条件的限制和实际需要，水位站的布设密度一般是不足的，对测区面积较大的测区而言，当潮汐传播存在纵横向比降时，使得“模拟法”的三角分带有许多麻烦和操作误差难以达到精度要求，故“三角分带法”在生产过程控制中难以应用。 近年来，引进了HYPACK MAX 水道测量软件，虽然能在水位改正的自动化处理取得了比较满意的效果，但是，在进行水位数据处理时，只能采用建立单站水位曲线的方法进行水位改正，对于在同一测区必须使用二个或二个以上水位站进行分带改正的水位数据无法处理。为此，在这种情况下我们仍采用传统方法，传统的做法是：根据整点或半点水位观测值，先在厘米方格纸上人工绘制组合（一条或多条）水位曲线图，要求非观测时刻的水位改正值，在组合水位曲线图上摘取瞬时潮高（或潮位）值，然后人工将潮高（或潮位）数据输入计算机进行数据处理，这种方法，对于大潮汛、潮差大的时段分带改正，免不了会出现一些“跳带”现象，造成局部区域潮位修正失真。 于是，开发和编制加权系数平均值计算法确定测深点的水位改正值的方法。 加权系数平均值的计算理论依据及条件和数学模型 在感潮河口段往往由于受水下地貌形态和河岸状态的影响，水位的推移速度要比外海快的多，除某点正处在潮汐的高潮或低潮时间外，该点的瞬时水位与其测区内潮波同相的各水位站的水位是呈相关关系，即各站间潮时和潮高变化与距离成正比。此外，在有限的感潮河口范围内，视同相潮波点的集合为一个空间平面。故将瞬时海面视作平面，测深定位点到任何一个水位站沿线瞬时水位变化是呈线性关系的。因此，测深定位点的水位值可以用具有相关关系的水位站的瞬时水位加权平均值来表示，于是，建立数学模型如式（1）所示。 Hi（x，y，t）= Ki * hi （x，y，t）·····（1） 式中：Hi（x，y，t）--为测深定位点的瞬时潮位值; hi（x，y，t）--为第i个潮位站的瞬时潮位值; Ki--为第i个潮位站的加权系数。 i与该站到测深定位点的距离di成反比，即Ki = 1/di . 测深定位点离水位站越近，加权系数Ki越大。 2.hi是从深度基准面起算的同相潮高 3．加权系数应为不大于1的正实数。 4．加权系数和等于1。 5．当测深定位点无限接近某一水位站时，该站的加权系趋于1，而其余水位站的加权系数趋于零。 根据《水运工程测量规范》，为满足计算水位的精度要求，水位站的布设可根据港口的实际情况，尽量分布合理，例如：长江口测区从吴淞口到长江口灯船，布设有吴淞、长兴、横沙、北槽中、中浚、牛皮礁、绿华山等水位站，对于某一测深定位点，水位改正值可根据实际情况采用一个或二个、或三个水位站的瞬时水位加权平均计算。表（一）列出了二个、三个水位站加权系数的计算公式。 二个、三个水位站加权系数的计算公式表 表（一） 站数 二个站 三个站 满足条件 0≤Ki≥1 K1:k2=1/d1:1/d2 K1+k2=1 0≤Ki≥1 K1:K2:K3=1/d1:1/d2:1/d3 K1+k2+k3=1 加权系数 D=d1+d2 K1=d2/D K2=d1/D D=(d1*d2)+(d2*d3)+(d3*d1) K1=(d2*d3)/D K2=(d3*d1)/D K3=(d1*d2)/D 加权潮位平均值 H1=K1*Ha H2=k2*Hb Hpoin