DB3206T1068-2024近岸池塘养殖海洋灾害风险隐患评估技术规程.docx

ICS 07.060CCS A 453206南通市地方标准DB 3206/T 10682024近岸池塘养殖海洋灾害风险隐患评估技术规程Technical code of practice for hydrodynamic risk assessment of nearshor pond aquaculture2024-01-30 发布2024-01-30 实施南通市市场监督管理局发 布DB 3206/T 10682024目 次前言II1 范围12 规范性引用文件13 术语和定义14 基本要求25 评估步骤26 评估成果4附录 A（规范性） 海水养殖池塘土堤坝筑堤材料及堤身边坡5附录 B（规范性） 气压场和风场计算6附录 C（规范性） 风浪要素计算7附录 D（规范性） 风暴潮数值模型8附录 E（规范性） 浪爬高计算10附录 F（规范性） 池塘内风浪计算12附录 G（规范性） 池塘内壅水面高度13附录 H（规范性） 近岸池塘养殖海洋灾害风险隐患评估技术报告格式要求14参考文献15II前 言本文件按照GB/T 1.12020标准化工作导则 第1部分：标准化文件的结构和起草规则的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。 本文件由南通市自然资源与规划局提出并归口。本文件起草单位：河海大学、南通市自然资源和规划局、国家海洋局烟台海洋环境监测中心站、南通沿海开发集团有限公司、中交上海航道勘察设计研究院有限公司。本文件主要起草人：王岗、陶爱峰、罗锋、刘子铭、严锐、连璐璐、徐辉、杨宏兵、陈立华、纪灵、洪新、李庆杰、侯查伟、武双、付睿丽、范骏、熊焱、李鹏宇、曹力玮、周光准、曾靖伟、田 鹏、石斌、季毅、顾恒瑞。DB 3206/T 10682024近岸池塘养殖海洋灾害风险隐患评估技术规程1 范围本规程规定了近岸池塘养殖海洋灾害风险隐患评估的基本要求、评估步骤、评估成果。 本规程适用于南通市近岸池塘养殖的水动力灾害风险隐患评估。2 规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T 12763.2-2007 海洋调查规范第2部分：海洋水文观测GB/T 12763.3-2020 海洋调查规范 第3部分:海洋气象观测GB/T 15920-2010 海洋学术语 物理海洋学3 术语和定义下列术语和定义适用于本文件。3.1池塘养殖 pond aquaculture利用人工开挖或天然池塘进行水生物经济动物养殖的生产方式。3.2海浪 ocean wave由风引起的海面波动现象，主要包括风浪和涌浪。风浪和涌浪的定义见GB/T 15920-2010。3.3波浪爬高 wave run-up建筑物上波浪上爬的最高点相对静止水面的高度。3.4漫滩 floodplain风暴潮过程淹没海岸线以上陆地的现象。3.5漫堤 dyke overtopping风暴潮过程引起的沿岸涨水越过海堤顶堤顶的现象。3.6风暴潮 storm surge有热带气旋、温带天气系统、海上飑线等风暴过境所伴随的强风和气压骤变而引起局部海面振荡或非周期性异常升高（降低）现象。64 基本要求4.1 总则基于近岸池塘养殖受外海风暴潮与波浪爬高影响下的漫滩影响所导致的灾害影响，开展近岸池塘养殖海洋灾害风险隐患评估，为近岸池塘养殖灾前安全评价、灾中风险预警、灾后受损评估提供依 据，为设施渔业保险服务提供技术支撑。4.2 科学性基于海水池塘养殖设施的受灾机制，综合考虑外海及池塘内部水动力影响，开展海水池塘养殖的风险评估。4.3 评估范围本规程限于海岸线向陆一侧15公里内或地表高程小于平均高潮位2米的海水池塘养殖区。5 评估步骤5.1 资料收集根据风险评估对象，确定评估区域范围，收集和整理研究区域基础地理信息、养殖情况、海洋灾害、气象要素及海防工程等。5.1.1 地理信息资料包括研究区入海河流、地表高程、水深及岸线数据，比例尺不小于1:10万。5.1.2 养殖池塘基本情况包括养殖池塘的分布图、各池塘平面尺寸、池塘水深、堤坝的堤身参数及筑堤材料。5.1.3 海防工程包括海堤的堤线布置、堤身设计参数及设计标准。5.1.4 水文、气象资料包括研究区附近的潮（水）位站、海洋观测站、浮标船舶报、卫星高度计及国家海洋局发布的预报结果等风暴潮过程的潮位、潮流、海浪及气象观测资料。水文、气象资料收集应遵照GB/T 12763.2- 2007和GB/T 12763.3-2020的有关规定执行。5.2 评估对象分析明确外海防潮堤的防护标准、养殖池塘平面布置及堤坝的堤身参数与筑堤材料力学特性。5.2.1 海水池塘a) 通常研究区内的海水养殖池塘面积及形状不完全相同，将池塘统一成长宽比为 3:1 的矩形， 单池面积分别为 30 亩、50、100 亩和 300 亩；b) 池塘水深以 d = 1.5m 和 d = 2.0m 两种情况为例；c) 池塘间堤坝为堤顶高程 2.5m 的梯形，边坡根据筑堤材料的土质而定，具体参考附录 A。5.2.2 外海堤坝的堤顶高程a) 有明确设计文档的外海堤防，其堤顶高程为设计水位+堤顶超高（壅水面+波浪爬高+安全超 高）;b) 未有详实设计材料的，可直接量测。5.3 受灾风险等级确定基于南通市贝类养殖调研数据，根据海水池塘所遭受不同水动力威胁，确定风险等级。海水池塘养殖受灾风险依据其受外海及池塘内水体的水动力影响，确定为高风险区(I级)、中风险区(II级)和低风险区(III级)三级（表 1）。表 1 海水池塘养殖风险等级表风险等级判别标准受灾情况I外海壅水高程+波浪爬高超过外海堤坝设计标准。海水越过外海防护堤坝几乎淹没评估区全部海养殖池塘；养殖物及养殖设施几乎全部损毁；滞留于此的养殖人员面临生命安全危险。II30 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.5m 或 50 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.6m 或 100 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.8m 或300 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 1.2m。池塘中有大量海水溢出；堤坝上有剧烈的水流，对置于其上的养殖设施造成严重损毁；严重威胁滞留的养殖人员生命安全。III30 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.2m 或 50 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.3m 或 100 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.5m 或300 亩的海水池塘堤坝壅水面+波浪爬高超过堤顶高程 0.7m。池塘中有明显海水溢上堤顶；对养殖设施造成一定损坏；威胁滞留的养殖人员人身安全有。5.4 水动力要素计算5.4.1 概述计算海水养殖池塘内风浪、壅水面高程及波浪爬高。5.4.2 外海水动力计算海洋预报部门发布的评估区风暴潮预警级别小于其外海堤防设计水位时池塘养殖不受外海水动力过程影响；若评估区风暴潮预警级别达到或超过其外海堤防设计水位时，需计算评估区外海的波浪爬高与风暴潮壅水。台风气压及风场计算参考附录B，风浪要素计算参考附录C，风暴潮数值计算参考附录D，波浪爬高计算参考附录E。若外海壅水高程+波浪爬高超过外海堤坝设计标准，则确定海水池塘养殖风险为I级。5.4.3 养殖池塘内水动力计算灾害期间池塘内海水在风的作用下产生的波浪爬高和壅水是海水池塘养殖受灾风险的主要动力因素。池塘内风浪计算参考附录 F，池塘内风浪在池塘堤坝的爬高由计算所得风浪结合附录 E 计算得到，池塘内壅水面高度参考附录 G。由计算 30 亩、50 亩和 100 亩海水池塘堤坝壅水面+波浪爬高结果，结合表 1 确定海水池塘养殖风险等级。5.5 风险评估依据所确定的海水池塘养殖受灾风险等级及水动力要素，评估研究区养殖风险。海水池塘养殖受灾风险区划图用红、橙、黄三色标识出评估区各面积单元等级大小，与表1对应，红色对应I级,橙色对应II级，黄色对应III级。6 评估成果6.1 资料汇总汇总资料包括研究区域自然地理环境概述，海水养殖概况，海水池塘分布，典型池塘平面布置， 堤身几何参数，筑堤材料，外海风暴潮壅水、最大波高、波浪爬高，池塘内风浪波高与周期、壅水高程及波浪爬高，受灾风险等级等。6.2 报告编制根据风险分析和风险区划结果绘制相应成果图件并按技术报告格式编制近岸池塘养殖海洋灾害风险评估报告，主要图表及格式要求见附录H。附 录 A（规范性）海水养殖池塘土堤坝筑堤材料及堤身边坡A.1 养殖池塘土堤坝修筑工程主要为细颗粒土，工程土的分类指标应符合表A.1的规定。表 A.1 工程土的分类表土的分类土名分类标准粉土类粉土d500.075mm的颗粒含量小于总质量的50% IP10黏性土类粉质黏土1017淤泥质土类淤泥质粉质黏土36%w55% 1.0e1.5 10IP17淤泥质黏土36%w55% 1.017淤泥类淤泥55%w85% 1.517流泥类流泥85%17表中：d50土的中值颗粒粒径，单位为毫米（mm）； IP塑性指数；w含水率；e孔隙比。A.2 黏性土类的天然状态可根据标准贯入试验击数按表A2确定。表 A.2 根据标贯击数确定黏性土的天然状态击数NN22N44N88N15N15天然状态很软软中等硬坚硬A.3 粉土类和中等、硬、坚硬黏性土类土堤坝放坡比宜为1:11:1.5，其余土类土堤坝放坡比宜为1:1.51:3。A.4 淤泥类和流泥类土堤坝宜采用抛石挤淤、插打木桩或浅层搅拌水泥固化作地基处理。A.5 土堤坝坡面可采用膜袋装土压实或格宾石笼抗冲刷。DB 3206/T 10682024附 录 B（规范性）气压场和风场计算采用梯度风原理计算气压场和风场。在气旋风外围背景风场可以采用CCMP、NCEP、ERA等资料或气象数值预报产品提供风场。热带气旋气压场可以采用气压场模型计算，常用气压场理论模型可参考藤田模型：7) r2 1/ 2P = P - ( P - P01 + 2 R (B.1)式中：P计算点气压，单位为帕(Pa)； P0气旋中心气压，单位为帕(Pa)； P气旋外围气压，单位为帕(Pa)；r为计算点至气旋中心距离，单位为米(m)； 0 R0表征气旋系统特征参数，可根据最大风速半径调整，也可根据台站观测气压调整，单位为 米(m)。气旋模型风场计算可按：r = c W -sin(q+b) +r(B.2)W1 1 cos(q+b) c2W2W1 =- fr f 2 r 2 + r P4ra r2(B.3)W = Vx exp -pr(B.4)2V 500000 式中： y 海面上10m高度处的气旋风速，单位为米每秒(m/s)；W1梯度风速，单位为米每秒(m/s)；W2气旋移行产生的风场，单位为米每秒(m/s)；计算点(x,y)与气旋中心的(xc,yc)连线与x轴夹角，单位为弧度(rad)； 梯度风和海面风的夹角，单位为弧度(rad)；c1，c2为修正系数；f为科氏系数，f= 2sin，为地转角速度，单位为弧度每秒(rad)，为纬度，单位为弧度(rad)；a空气密度，单位为千克每立方米(Kg/m3)；Vx、Vy气旋中心在x和y方向上移动速度，单位为米每秒(m/s)。气旋过程中的中心附近风场可以采用上述梯度风原理计算风场，在气旋风外围风场可以采用CCMP、NCEP、ERA等资料或者气象数值预报产品提供风场，如WRF。DB 3206/T 10682024附 录 C（规范性） 风浪要素计算采用谱浪模型并考虑风能驱动计算波浪场。该模型基于能量守恒原理的波作用量平衡方程数值求解。考虑了水深变化、背景流和障碍物等对波浪传播的影响，同时考虑了包括风摄入波动能量、白浪耗散、水深变化导致的波浪破碎、底摩擦耗散、波-波非线性相互作用等物理过程。荷兰Delft理工大学开发的SWAN模型和丹麦DHI研究所开发的Mike-SW模型等，均是目前应用较广的第3代风浪谱模型， 适用于海岸、湖泊和河口地区的风浪要素计算。笛卡尔坐标系下的波作用量平衡方程为：8 N + c N + c N + c N + c N = S(C.1)式中：tx xy ys sq qsN为波作用量密度随时间的变化率,Cx、Cy、Cs、Cq分别代表在x、y、s、q方向的波浪传播速 度。C = 1 2kdskxx2 1+ sinh(2kd ) k 2+Ux(C.2)C = 1 2kdskyy2 1+ sinh(2kd) k 2+Uy(C.3)sdurr urC =+U d - C k Usd tgs(C.4)C = -1 s d ur + rUk 式中：sk d mm (C.5)kxyUxyr = (k , k) 为波数，d为水深， ur = (U ,U) 为流速，s是沿波向q方向的空间坐标，m为垂直于s的坐标，相对频率=kU+，为波浪的固有频率。右边 S 代表以谱密度表示的源汇项，包括风能输入、波与波之间非线性相互作用和由于底摩擦、白帽、破碎等引起的能量损耗：S = Swind (s,q)+ Snl + Sds + Sbot + Ssurf(C.6)式中：Swind指风输入的能量；Snl 指波与波之间的非线性作用引起的能量耗散；Sds指由白帽引起的能量耗散；Sbot指由底摩阻引起的能量耗散；Ssurf指由于水深变化引起的波浪破碎产生的能量耗散。DB 3206/T 10682024附 录 D（规范性） 风暴潮数值模型风暴潮数值计算的气压场、风场和波浪场由可由附录B和附录C计算确定。 潮流运动可按下列控制方程：连续方程z + (h +z)u + (h +z)v = 0(D.1)13x 方向动量方程txyu + u u + v u - fv = -g z -tbx-1(Sxx + Sxy ) txyxr(h +z)r(h +z) xy-1r(h +z)( Rxx +xRxyy) + Nx (2ux2+ 2uy2) + twx(D.2)r(h +z)y 方向动量方程v + u v + v v + fu = -g z-tby-1(Syx + Syy ) txy-1RyxyRyyr(h +z)2vr(h +z)2vxtwyyr(h +z) ( x+ y ) + Ny (x2 + y2 ) + r(h +z)(D.3)式中：t时间，单位为秒(s)；x、y原点 o 置于某一水平基面的直角坐标系坐标；ru、流速矢量分别为V 沿 x，y 方向的分量，单位为米每秒(m/s)；为相对 xoy 平面的水面水位，单位为米，(m)；h为静水深；g重力加速度，单位为米每立方秒(m/s2)； 为海水密度，单位为千克每立方米(kg/m3)； Nx、Ny分别为 x、y 方向的水平涡动粘性系数；f为科氏系数，f= 2sin，为地转角速度，单位为弧度每秒(rad)，为纬度，单位为弧度(rad)；Sxx、Syy、Syx、Syy波浪辐射应力张量的四个分量，单位为牛顿每米(N/m)。考虑波浪破碎的波浪辐射应力张量的四个分量可按下列表达式确定：rgH 2 Cg1 CgrM 2S=w 2- - sin 2 q -Wx(D.4)xx8C2 Ch + zrgH 2 CgrMWx MwyS = S=wsinqcosq-(D.5)xyyx8Ch +zrgH 2 C1 C rM 2S=w g - +g sin 2 q -Wy(D.6)式中：yy8 C2 Ch + zCg、H、Hw、分别为波浪群速、波浪相速、波高和波向；wMWx 、 MW 由 卷 破 波 和 波 动 Stokes 漂 移 造 成 的 向 岸 质 量 通 量 ， MWx=(E+2Er)/(C)cos,MWx=(E+2Er)/(C) sin； 其中，单位面积波能 E=gH 2/8； 单位面积卷波能量Er,Er=4rC/(gT)PE；r0.9；T 为波浪周期；P=H/(h+)4，=0.5+(H0+L0)1/2，H0 为深水波 高，L0 为深水波长；Rxx、Rxy、Ryx、Ryy 破波波卷产生的切应力， 单位为牛顿每米(N/m)， Rxx=2 Ercos2 ， Ryy=2 Ersin2，Rxy= Ryx= Ersin2；wx 和wy风应力矢量w 沿 x、y 方向的分量，单位为牛顿每平方米(N/m2)；bx 和by波浪、潮流共同作用下的底部剪切应力矢量沿 x、y 方向的分量，单位为牛顿每平方米(N/m2)，可按下列公式确定：t = g rUr up f Uru + Br1c22 gfw 2 Ur u(D.7)bxc2+ 8 r ww wpwt = g rUr vp f Urv + Br1c22 gfw 2 Ur v (D.8)式中：byc2+ 8 r ww wpw水流速度矢量，单位为米每秒(m/s)；波浪底部水质点最大速度矢量，单位为米每秒(m/s)；uw、vw分别为 在 x 和 y 方向的分量，单位为米每秒(m/s)；fw波浪摩阻系数，通常取 0.010.02；c谢才系数, c=(h+)1/6/n； n曼宁糙率系数；B波流相互作用系数：当波流同向时 B=0.917，当波流相互垂直时 B=-0.1983，当波流夹角不定时 B=0.359。定解条件包括边界条件和初始条件，对于边界条件，开边界取流速或潮位的实测或分析值，固边界取流速的法向分量为零；对于初始条件，流速一律取零值，潮位取初始时刻的值。附 录 E（规范性） 浪爬高计算在风的直接作用下，正向来波在单一斜坡上的波浪爬高可按下列要求确定:1) 当斜坡坡率 m=1.55.0、H/L0.025时，可按下列公式计算：pR = KD KV KPHL(E.1)1+ m2m= cota(E.2)式中：RP 累计频率为 P 的波浪爬高，单位为米(m)；K 斜坡的糙率及渗透系数，根据湖面类型按表 E.1 确定；KV 经验系数，可根据风速 U，单位为米每秒(m/s)、堤前水深 d，单位为米(m)、重力加速度g，单位为米每平方秒(m/s2)组成的无维量 U/(gd)1/2，按表 E.2 确定；KP 表示 RP 和平均爬高R 比值 RP/R 的爬高累计频率换算系数，可按 E.3 确定。对不允许越浪的堤防，爬高累计频率宜取 2%；对允许越浪的堤防，应根据越浪量大小，采取相应的防护措施；m 斜坡坡率； 斜坡坡脚，单位为度()；H 堤前波浪的平均波高，单位为米(m)； L 堤前波浪的平均波长，单位为米(m)。2) 当 m1.0、H/L0.025 时，可按下式计算：R = K K K R H(E.3)式中：pD V P 0R0 无风情况下，光滑不透水护面(K=1)、H=1m 时的爬高值，单位为米(m)，可按表 E.4 确定。3) 1.0m1.5 时，可由 m=1.0 和 m=1.5 的计算值按内插法确定。表 E.1 斜坡的糙率及渗透系数 K护面类型K光滑不透水护面（沥青混凝土、混凝土）1.0混凝土板0.95草皮0.90砌石0.80抛填两层块石（不透水堤心）0.600.65抛填两层石块（透水堤心）0.500.55注 1：m1.0，砌石护面取 K=1.0。表 E.2 经验系数 KVU/(gd)1/211.522.533.545KV11.021.081.161.221.251.281.30表 E.3 爬高累积频率换算系数 KPH/dP(%)0.32.131.861.761.701.651.611.481.401.310.99表 E.4 R0 值m=cot00.51.0R01.241.452.20带有平台的复式斜坡堤的波浪爬高，可先确定该断面的这算坡率 me，再按坡率为 me 的单坡断面确定其爬高。折算坡率 me 可下列公式计算：1) m=(m 下 m 上)=0 时：m = m （1-4.0 dw ）K (E.4)e上Lb2) 当m0 时：K = 1+ 3 B(E.5)bLm = (m + 0.3D m- 0.1D m2 )(1- 4.5 dw )K (E.6)e上Lb3) 当m0 时：m = (m + 0.5D m+ 0.08D m2 )(1+ 3.0 dw )K (E.7)e上Lb式中：m 上 平台以上的斜坡坡率；m 下 平台以下的斜坡坡率；dw 平台水深，当平台在静水位以下时取正值；平台在静水位以上时取负值；| dw |表示取绝对值，单位为米(m)；B平台宽度，单位为米(m)； L波长，单位为米(m)。折算坡率法适用条件：m 上=1.04.0，m 下=1.53.0，dw/L=-0.025+0.025，0.05B/L0.25。E.3 当来波波向线与堤轴线的法线成角时，波浪爬高应乘以系数 K，当堤坡坡率 m1 时，K可按表 E.5 确定。表 E.5 系数 K()15203040506090K10.960.920.870.820.760.6附 录 F（规范性） 池塘内风浪计算海水养殖池塘内风浪的有效波高、有效波周期按式(F.1)和(F.2)确定。 gd 0.8 gH gF 0.35 U 2 s = 0.0055tanh 30 (F.1)U 2 U 2 gF 0.35 U 2 gd 0.8 gT gF 0.233 U 2 s = 0.55tanh2 / 3 30 (F.2)U U 2 gF 0.35 式中：Hs有效波高，单位为米（m）；g 重力加速度，单位为米每平方秒（m/s2）；U 2U海面上10m高度处的平均风速，单位为米每秒（m/s）； F风区长度，单位为米（m），取池塘长和宽的大值； d池塘平均水深，单位为米（m）；Ts有效波周期，单位为秒（s）。DB 3206/T 10682024附 录 G（规范性）池塘内壅水面高度14海水养殖池塘内的壅水面高度按式(G.1)确定e =KU 2 F 2gd(G.1)式中：e风壅水面高度，单位为米（m）； K综合摩擦系数，可取3.610-6；g 重力加速度，单位为米每平方秒（m/s2）；U海面上10m高度处的平均风速，单位为米每秒（m/s）； F池塘长，单位为米（m）；d池塘平均水深，单位为米（m）。DB 3206/T 10682024附 录 H（规范性）近岸池塘养殖海洋灾害风险隐患评估技术报告格式要求H.1 封面封面书写内容包括：报告名称；委托单位名称；承担单位名称(盖章)；报告编制日期。H.2 封二封二应书写内容包指:承担单位负责人；任务责任人；技术负责人；报告编写人员；主要参与人员；审核人员。H.3 目录报告应有包含各级标题的目录页。H.4 前言前言包括工作来源、背景、内容和主要成果等。H.5 正文技术报告正文编写内容大纲如下:第1章“评估区域概况”，内容包括自然地理环境概述、海水养殖概况、海洋灾害概况等。第2章“海水池塘概况”，内容包括主要海水池塘分布范围、海水池塘平面形状、堤坝的堤身及筑堤材料情况、养殖活动、外海防潮堤堤线布置、堤身设计参数及设计标准。第3章“水动力要素计算”，内容包括天文潮最大潮位、最大流速、壅水高程、最大波高、波浪爬高及漫堤、漫滩范围，池塘内风浪波高、周期、壅水高程及波浪爬高；以图表形式给出最大潮 位、最大流速、最大波高、波浪爬高及漫堤、漫滩范围。第4章“海上池塘养殖受灾风险评估”，内容包括研究区养殖受灾风险评估方法及评估结果； 以图表形式给出研究区受灾风险等级。第5章“对策与建议”，内容包括根据研究区域养殖受灾风险评估等级结果阐述其具体影响， 提出有针对性的对策建议。H.6 封底印刷版报告宜有封底。封底可放置任务承担单位的名称和地址或其他相关信息，也可为空白页。H.7 报告格式报告文本外形尺寸为A4（210mm297mm）。17参 考 文 献1GB/T14914.2-2019 海洋观测规范 第2部分:海滨观测2GB/T19721.1-2017 海洋预报和警报发布第1部分：风暴潮预报和警报发布3GB/T20014.14-2013 良好农业规范 第14部分：水产池塘养殖基础控制点与符合性规范4GB/T28923.2-2012 自然灾害遥感专题图产品制作要求第2部分：监测专题图产品5GB/T35227-2017地面气象观测规范 风向和风速6GB/T39616-2020卫星导航定位基准站网络实时动态测量（RTK）规范7GB/T50095-2014水文基本术语和符号标准8GB/T50159-2015河流悬移质泥沙测验规范9GB/T50228-2011工程测量基本术语标准10 GB/T 50286-2013 堤防工程设计规范11 GB/T 510152014 海堤工程设计规范12JTS132-2015水运工程水文观测规范13JTS145-2015港口与航道水文规范14JTS147-2017水运工程地基设计规范