[0096] S231、基于一维圣维南方程组，建立潮波传播的运动方程如下：
[0098]同时，建立潮波传播的连续方程如下：
[0100] 其中，v表示X方向的沿程流速，g表示重力加速度，h表示水深，Ib表示河床断面的 床面坡降，η为曼宁糙率，S表示河床断面的断面面积，Q表示河床断面的断面流量；
[0101 ] S232、结合拟合获得的沿程指数将潮波传播的连续方程简化如下：
[0103] 其中，β表示拟合获得的沿程指数；
[0104] 简化原理：因
，且断面流量等于断面面积与沿程流速的乘 积，则Q = Sv，因此连续方
[0106] S233、将步骤S231的欧拉形式的运动方程转换为如下的拉格朗日形式的运动方 程：
[0108] 上式中，c表示波速，a表示振幅，ε表示相位，Η表示潮差；
[0109] S234、根据步骤S232中的公式对S233的公式进行简化，获得高潮、低潮时刻潮波传 播的运动方程如下：
[0111] 上式中，?表示实测的平均水深数据，f为自定义的过渡参数，且
[0112] S235、求解步骤S234的公式，获得河口潮差预报模型如下：
[0114] 其中，y表示一无量纲参数，α表示河口湾的参数，γ表示河口地形沿程衰减的尺度 参数，Η0表示河口起点处的实测潮差值。
[0115] 另外，曼宁糙率η、波速c，振幅a、相位ε均为经验值，通过预设设定。
[0116] S23、建立潮波传播的一维圣维南方程组，结合河床断面与沿程距离的指数曲线对 其进行简化后，构建河口潮差预报模型。
[0117] 更详细的，步骤S234中，高潮、低潮时刻潮波传播的运动方程的简化推算过程如 下：
[0118]用a表示潮波的振幅，c表示波速，Η表示潮差，则有：
[0120]将这些内容代入式步骤S232的简化公式，即为
[0128] 由于高潮（HW )、低潮（LW )时刻是潮汐涨落过程中比较特殊的时刻，
;而且潮差与高低潮水位可以建立起简单的关系式，即H = hHW-
为平均水面坡降。其中下标HW表示高潮，LW表示低潮。涨憩(HWS)或落憩(LWS)滞后于高潮、 低潮时刻，用ε表示高潮与涨憩的相位差，则vhw=a s i η ε，vlw=_as i η ε。
[0129] 因此可获得高潮时刻的控制方程为：
[0131]低潮时刻的控制方程为：
[0133]由于huKhm，波速CLW〈CHW，基本上，高潮时波速与水深的比值与低潮时相同，假定
c +c ?& 高潮时刻的彳空制方程减去彳氐潮时·刻的彳空制方 程，可获得：
[0139]等式右边第一项表示河口宽度的沿程距离所引起的潮汐增强或减弱，等式右边第 二项表示潮汐运动为克服平均水面坡降所产生的变化，右边第三项表示河道摩擦对潮汐变 化的影响。根据量纲分析，大部分河口平均水深的沿程坡降（等式右边第二项）相比起右边 第一项来说小很多，基本上可忽略不计，上式简化为
[0141] 上式即是描述河口湾潮汐沿程运动的简单控制方程，即河口潮差预报模型。
[0142] 实施例二
[0143] 以珠江河口为示例，采用实施例一的方法，开展城市河涌防洪排涝动态预报控制。
[0144] 步骤1、河口起点处的实时的潮位数据和平均水深数据，采集装置包括水位计，可 将采集得到的实时雨量数据存储于数据存储器中，采集珠江河口伶仃洋、磨刀门、黄茅海起 点处的潮位、平均水深等数据，采集的河口起点处的实测潮差值Ho和平均水深数据E如下表 1所示：
[0145] 表1
[0146]
[0147] 步骤2、数据传输与接收
[0148] 采用无线通讯方式，通过GPRS通信方式将采集的数据发送到控制终端的数据接收 器，该数据接收器可包括接收天线、电台、计算机和稳压电源等组成。
[0149] 步骤3、结合所采集的数据，拟合河床断面与河口沿程距离之间的关系曲线
[0150] 例如，剖分珠江河口各河床断面，统计河口各断面的断面面积、河口宽度两个参数 的沿程距离，绘制河床断面与沿程距离的指数曲线，分别如图2、图3、图4所示，然后根据下 式建立河床断面与沿程距离的指数曲线，进而进行曲线拟合沿程指数β:
[0152]最后获得的河床断面的断面面积的沿程距离曲线、河床断面的河口宽度的沿程距 离曲线分别如下：
[0153] 伶仃洋:断面面积沿程距离:S = 348225e-Q·Q271x R2 = 0.9049
[0154] 伶仃洋:河口宽度沿程距离:B = 49728e-Q'Q274x R2 = 0.8555
[0155] 黄茅海：断面面积沿程距离：3 = 705946-°·'931 R2 = 0.9407
[0156] 黄茅海:河口宽度沿程距离:8 = 297486-°··831 R2 = 0.9128
[0157] 磨刀门：断面面积沿程距离:S = 10835e-Q'Q165x R2 = 0.758
[0158] 磨刀门：河口宽度沿程距离:B = 49728e-Q'Q36x R2 = 0.8119 [0159 ] R2表示拟合度，R2越接近1，表示拟合度越好。
[0160]步骤4、河口潮差预报模型确定
[0161 ]将拟合获得的沿程指数代入下式的河口潮差预报模型的公式，得到河口沿程的潮 差：
[0163] 本实施例中，预设设定各经验值的取值范围如下：曼宁糙率η为0.015~0.045、波 速c为0.8~1.2，振幅a为0.2~0.8、相位ε为2小时。更具体的，珠江口各河口湾的相关参数 如下表2所示：
[0164] 表2
?〇166]^步骤5、采用步骤4的河口潮差预报模型进行潮差预测，所获得的预测曲线如图5所 示，图5中显示的是预测计算得到的潮差值与河口起点处的实测潮差值的比值，同时以沿程 站位的平均潮差作为验证数据，在图中标出，因此可以看出模型预测结果与沿程观测数据 吻合。图5中三种不同标记处的直线分别代表对应的河口的预测曲线。
[0167] 因此，本发明的潮差快速预报方法及系统能够有效、准确地实时动态预测河口沿 程的潮差，从而预报河口潮差沿程距离情况。而且本发明采用解析模型进行求解，与现有技 术中采用数学模型的预报方法相比，其计算复杂度较低，计算效率高。
[0168] 以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施 例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替 换，这些等同的变型或替换均包含在本申请权利要求所限定的范围内。
【主权项】
1. 一种河口潮差快速预报方法，其特征在于，包括步骤： 51、 采集河口起点处的实时的潮位数据和平均水深数据； 52、 结合所采集的数据，拟合河床断面与河口沿程距离之间的关系曲线，并构建河口潮 差预报模型； 53、 将采集的实时的潮位数据和平均水深数据代入河口潮差预报模型后，计算获得河 口沿程各断面的潮差预测值。2. 根据权利要求1所述的一种河口潮差快速预报方法，其特征在于，所述步骤S2，包括： 521、 收集河口河道地形图，W河宽的两倍距离作为间隔，剖分河床断面； 522、 建立河床断面与沿程距离的指数曲线，并拟合获得沿程指数； 523、 建立潮波传播的一维圣维南方程组，结合拟合获得的沿程指数对其进行简化后， 构建河口潮差预报模型。3. 根据权利要求2所述的一种河口潮差快速预报方法，其特征在于，所述步骤S21，其具 体为： 收集河口河道地形图，W河宽的两倍距离作为间隔，剖分河床断面，并记录河床断面的 断面间距、离岸距离、断面面积W及河口宽度作为河床断面信息。4. 根据权利要求2所述的一种河口潮差快速预报方法，其特征在于，所述步骤S22，其具 体为： 根据下式建立河床断面与沿程距离的指数曲线，进而根据剖分的河床断面的信息，拟 合获得沿程指数β:上式中，S表示河床断面的断面面积，Β表示河床断面的河口宽度，So表示河口起点处河 床断面的断面面积，Bo表示河口起点处河床断面的河口宽度，X表示沿程距离值。5. 根据权利要求2所述的一种河口潮差快速预报方法，其特征在于，所述步骤S23包括： 5231、 基于一维圣维南方程组，建立潮波传播的运动方程如下：同时，建立潮波传播的连续方程如下：其中，V表示X方向的沿程流速，g表示重力加速度，h表示水深，Ib表示河床断面的床面坡 降，η为曼宁糖率，S表示河床断面的断面面积，Q表示河床断面的断面流量； 5232、 结合拟合获得的沿程指数将潮波传播的连续方程简化如下：其中，β表示拟合获得的沿程指数； 5233、 将步骤S231的欧拉形式的运动方程转换为如下的拉格朗日形式的运动方程：上式中，c表示波速，a表示振幅，ε表示相位，Η表示潮差； 5234、 根据步骤S232中的公式对S233的公式进行简化，获得高潮、低潮时刻潮波传播的 运动方程如下：上式中，拓表示实测的平均水深数据，且5235、 求解步骤S234的公式，获得河口潮差预报模型如下：其中，y表示一无量纲参数，α表示河口湾的参数，丫表示河口地形沿程衰减的尺度参 数，化表示河口起点处的实测潮差值。 S23、建立潮波传播的一维圣维南方程组，结合河床断面与沿程距离的指数曲线对其进 行简化后，构建河口潮差预报模型。6.根据权利要求2所述的一种河口潮差快速预报方法，其特征在于，所述步骤S3，其具 体为： 将拟合获得的沿程指数和采集的实时的潮位数据和平均水深数据代入河口潮差预报 模型的公式中进行求解，从而计算获得河口沿程各断面的潮差预测值。
【专利摘要】本发明公开了一种河口潮差快速预报方法，包括步骤：S1、采集河口起点处的实时的潮位数据和平均水深数据；S2、结合所采集的数据，拟合河床断面与河口沿程距离之间的关系曲线，并构建河口潮差预报模型；S3、将采集的实时的潮位数据和平均水深数据代入河口潮差预报模型后，计算获得河口沿程各断面的潮差预测值。本发明计算复杂度较低、计算效率高、预测准确度高，操作简单，可及时地进行潮差预测，可广泛应用于水利行业中。
【IPC分类】G06Q10/04, G06F17/50, G06F17/11
【公开号】CN105701066
【申请号】CN201610131547
【发明人】谭超, 黄本胜, 刘达, 邱静, 黄广灵
【申请人】广东省水利水电科学研究院
【公开日】2016年6月22日
【申请日】2016年3月8日