天津市人大常委会


天津市蓄滞洪区管理条例
天津市人大常委会


（1996年10月16日天津市第十二届人民代表大会常务委员会第二十七次会议通过 1996年10月16日公布施行）


第一条 为加强本市蓄滞洪区的安全、建设管理，发挥蓄滞洪区的功能，减少分滞洪水淹没损失，保障经济建设和人民生命财产安全，根据《中华人民共和国水法》及其他有关法律、法规，结合本市实际情况，制定本条例。
第二条 本条例所称蓄滞洪区，是指本市行政区域内临时贮存和调蓄洪水的青甸洼、盛庄洼、黄庄洼、大黄堡洼、永定河泛区、三角淀、七里海、东淀、文安洼、贾口洼、团泊洼、大港行洪道、淀北等十三处低洼地区。
前款规定的蓄滞洪区的运用、安全、建设管理，适用本条例。
第三条 蓄滞洪区的蓄滞洪运用，按照国家防汛指挥机关批准的分滞洪方案实施。分滞洪命令，由市防汛指挥机关发布。
分滞洪命令一经发布，任何单位和个人必须执行。
第四条 蓄滞洪区的安全、建设管理，实行市、区县、乡镇人民政府行政首长负责制。市和区、县水行政主管部门负责日常工作。
市和区、县的计划、城建、规划、土地、农业、乡镇企业、计划生育、公安、交通、环保、卫生等行政主管部门，按照各自的职责，会同同级水行政主管部门对蓄滞洪区的安全、建设实施管理。
第五条 蓄滞洪区的安全建设规划，由区、县水行政主管部门会同同级有关部门编制，经区、县人民政府同意，市水行政主管部门审核，报市人民政府批准。
蓄滞洪区安全建设规划的变更，必须经市人民政府批准。
第六条 蓄滞洪区的经济建设和社会发展规划，应当与蓄滞洪区的安全建设规划相协调。
制定蓄滞洪区的经济建设和社会发展规划，应当征询市水行政主管部门的意见。
第七条 蓄滞洪区的建设与发展必须严格按照规划进行。
第八条 蓄滞洪区的土地利用、开发和各项建设必须符合防洪的要求。在指定的分洪口门附近和洪水主流区域内，禁止修建阻碍行洪的各类建筑物。
第九条 在蓄滞洪区内新建永久性建筑物，应当避开洪水流路，选择较高地形，采取平顶、能避洪救人的结构形式。
在蓄滞洪区内新建、改建、扩建建设项目，建设单位应当按照国家有关规定经水行政主管部门审查同意后，方可按照基本建设程序履行报批手续。
蓄滞洪区内现有的建筑物未达到防洪标准的，应当进行加固或者采取其他补救措施。
第十条 蓄滞洪区内严禁新建、改建、扩建生产或者储存有毒、易爆等严重污染品和危险品的建设项目，避免次生灾害发生。
蓄滞洪区内现有的生产或者储存有毒、易爆等严重污染品和危险品的设施，必须在本条例施行之日起三年内迁出蓄滞洪区或者转产经营，迁出前必须自行采取安全保障措施。
本条例施行前已获批准建设属于本条第一款规定禁止建设的项目，在建未竣工的，应当停止建设；尚未开工的，不得开工建设。
第十一条 各级人民政府，应当严格控制蓄滞洪区内人口的增长，限制向蓄滞洪区内安置移民。
第十二条 蓄滞洪区所在地的区、县人民政府，应当组织水行政主管部门等编制蓄滞洪区居民安全转移预案。
安全楼和其他永久性建筑物上应当标明当地历史最高洪水位线。
第十三条 蓄滞洪区所在地的区、县人民政府，应当根据避洪撤离的需要，结合城乡建设，有计划地修建公路，指定居民撤离线路，落实居民临时安置地点。
第十四条 蓄滞洪区应当设置有线通讯和专用无线通讯两套通讯系统。有线通讯建设要纳入城乡电信网建设计划，优先安排实施；防汛专用无线通讯建设，由各级防汛机构负责组织实施。防汛报警设施可以根据当地条件和群众习惯设置，并将报警方式公布于众，确保报警及时有效。
防汛通讯必须畅通无阻，任何单位和个人不得干扰。
第十五条 禁止损毁蓄滞洪区内的套堤、导流堤。蓄滞洪区内原有的高地、旧堤，未经区、县水行政主管部门批准，不得拆毁。
第十六条 蓄滞洪区围堤的管理范围是从围堤内、外坡脚各向外延伸三十米；保护范围是从管理范围外沿各向外延伸三十米。
围堤管理范围内禁止采砂、取土、弃置砂石、爆破、打井、钻探、挖筑鱼塘、存放物料、建房、开采地下资源、考古发掘。
围堤保护范围内禁止采砂、取土、爆破、打井、钻探、挖筑鱼塘。
第十七条 禁止损毁蓄滞洪区安全设施。
蓄滞洪区所在地的区、县水行政主管部门应当建立健全安全设施档案。
安全设施的产权单位和管理单位，应当加强安全设施的养护管理，保证应急安全使用。
第十八条 市、区县、乡镇人民政府应当将蓄滞洪区安全建设和无线通讯、预警所需资金，列入年度财政预算。
蓄滞洪区的农业税原则返还，用于安全与建设。
第十九条 鼓励和支持在蓄滞洪区发展洪水保险事业。
蓄滞洪区内的单位和个人，根据自愿原则参加洪水保险。
第二十条 建立防洪基金制度。防洪基金主要用于蓄滞洪区的安全设施建设、维护和因分滞洪淹泡造成损失的补偿。
防洪基金的筹集和使用管理具体办法，由市人民政府另行制定。
第二十一条 在蓄滞洪区安全建设和管理工作中做出显著成绩的单位和个人，由各级人民政府和有关部门给予表彰和奖励。
第二十二条 违反本条例第三条第二款规定的，由其上级主管部门对直接责任人，给予行政处分；构成犯罪的，依法追究刑事责任。
第二十三条 违反本条例第八条、第十条、第十五条、第十六条、第十七条第一款规定的，由区、县以上水行政主管部门区别情况责令停产、停业，停止违法行为，限期清除障碍，采取补救措施，并可处以十万元以下罚款；有违法所得的，没收违法所得；造成损失的，责令赔偿损失。


第二十四条 违反本条例第十四条第二款规定的，由有关行政主管部门按照有关规定予以处理。
第二十五条 水行政主管部门和其他有关主管部门的工作人员玩忽职守、滥用职权、徇私舞弊的，由其所在单位或者上级主管部门给予行政处分；构成犯罪的，依法追究刑事责任。
第二十六条 违反本条例规定，应当给予治安处罚的，由公安机关依照《中华人民共和国治安管理处罚条例》的规定处罚。
第二十七条 当事人对行政处罚决定不服的，可以在接到处罚决定书之日起十五日内，向作出处罚决定的上一级机关申请复议；对复议决定不服的，可以在接到复议决定书之日起十五日内，向人民法院起诉。当事人也可以在接到处罚决定书之日起十五日内，直接向人民法院起诉。当事
人逾期不申请复议或者不向人民法院起诉又不履行处罚决定的，由作出处罚决定的机关申请人民法院强制执行。
第二十八条 本条例所称安全设施，是指为保障蓄滞洪区人民生命安全而修建的围村埝、安全台、安全楼、安全房、撤离路和购置的救生船，以及通讯、预警设备等避洪救生设施。
第二十九条 本条例自公布之日起施行。



1996年10月16日
国土资源部


关于统一图幅理论面积与图斑椭球面积计算要求的通知

（国土调查办发〔2008〕32号）


各省、自治区、直辖市第二次土地调查领导小组办公室，国土资源厅（国土环境资源厅、国土资源局、国土资源和房屋管理局、房屋土地资源管理局），解放军土地管理局、新疆生产建设兵团国土资源局：


面积计算是第二次土地调查的一项重要内容，国务院第二次全国土地调查领导小组办公室组织有关专家，依据《第二次全国土地调查技术规程》，对图幅理论面积与图斑椭球面积计算公式进行了细化，明确了面积计算方法，统一了公式中的有关参数，现将《图幅理论面积与图斑椭球面积计算公式及要求》予以印发，请各地严格遵照执行。



附：图幅理论面积与图斑椭球面积计算公式及要求



二〇〇八年三月二十七日

图幅理论面积与图斑椭球面积计算公式及要求

一、 图幅理论面积计算公式
（1）
式中：
a—椭球长半轴(单位：米)，α—椭球扁率，b—椭球短半轴(单位：米)。
е²﹦(a²﹣b²)/a²。
A﹦1﹢（3/6）е²﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8。
B﹦ （1/6）е²﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8。
C﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8。
D﹦ （1/112）е6﹢ （45/2304）е8。
E﹦ （5/2304）е8。
ΔL—图幅东西图廓的经差（单位：弧度）。
(B2﹣B1)—图幅南北图廓的纬差（单位：弧度），Bm﹦（B1﹢B2）/2。

二、椭球面上任意梯形面积计算公式
（2）
其中：A,B,C,D,E 为常数，按下式计算:
е²﹦(a²﹣b²)/a²
A﹦1﹢（3/6）е²﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8
B﹦ （1/6）е²﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8
C﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8
D﹦ （1/112）е6﹢（45/2304）е8
E﹦ （5/2304）е8
式中：a—椭球长半轴(单位：米)，b—椭球短半轴(单位：米)；
ΔL—图块经差(单位：弧度)； (B2﹣B1)—图块纬差(单位：弧度)
Bm﹦（B1﹢B2）/2。

三、高斯投影反解变换（ ）模型
（若坐标不带带号，则不需减去带号×1000000；）



+中央子午线经度值（孤度） （3）
式中：

公式说明：若坐标为没有带号前缀格式，则不需减去带号×1000000；若坐标为有带号前缀格式，则需减去带号×1000000。

四、计算用到的常数、椭球参数
在计算图幅理论面积与任意图斑椭球面积时，有关常数及保留的位数按给定数值计算。
常数：
π﹦3.14159265358979
206264.8062471
80椭球常数：
= 6378140 = 1/ 298.257
= 6356755.29
= 6.69438499958795E-03
= 6.73950181947292E-03
= 6399596.65198801
相关常数:
k0 = 1.57048687472752E-07
k1 = 5.05250559291393E-03
k2 = 2.98473350966158E-05
k3 = 2.41627215981336E-07
k4 = 2.22241909461273E-09

五、计算中的取位及要求
① 高斯投影反解变换后的B,L以秒为单位，保留到小数点后6位，四舍五入。
② 采用计算机计算时，所有变量数据类型均要定义为双精度。
③ 面积计算结果以平方米为单位，保留一位小数，四舍五入。
④ 各种比例尺标准分幅图经差、纬差见表1。
⑤ 在用大地坐标生成标准分幅图框时，要求在每条边框线的整秒处插入加密点。
表1 各种比例尺标准分幅图经差、纬差表
比例尺 1:100万 1:50万 1:25万 1:10万 1:5万 1:2.5万 1:1万 1:5千
经差 6º 3º 1º30′ 30′ 15′ 7′30″ 3′45″ 1′52.5″
纬差 4º 2º 1º 20′ 10′ 5′ 2′30″ 1′15″

六、任意图斑椭球面积计算方法
任意封闭图斑椭球面积计算的原理：将任意封闭图斑高斯平面坐标利用高斯投影反解变换模型，将高斯平面坐标换算为相应椭球的大地坐标，再利用椭球面上任意梯形图块面积计算模型计算其椭球面积，从而得到任意封闭图斑的椭球面积。
1、计算方法：
任意封闭区域总是可以分割成有限个任意小的梯形图块，因此，任意封闭区域的面积 ，式中Si为分割的任意小的梯形图块面积（i=1,2,…n）用公式(2)计算。
求封闭区域（多边形如图1）ABCD的面积 ，其具体方法为：
（1）对封闭区域（多边形）的界址点连续编号（顺时针或逆时针）ABCD，提取各界址点的高斯平面坐标A(X1,Y1)，B(X2,Y2)，C(X3,Y3)，D(X4,Y4)；
（2）利用高斯投影反解变换模型公式（3），将高斯平面坐标换算为相应椭球的大地坐标A(B1,L1)，B(B2,L2)，C(B3,L3)，D(B4,L4)；
（3）任意给定一经线L0（如L0＝60°），这样多边形ABCD的各边AB、BC、CD、DA与L0就围成了4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）；
（4）由于在椭球面上同一经差随着纬度升高，梯形图块的面积逐渐减小，而同一纬差上经差梯形图块的面积相等，所以，将梯形图块ABB1A1按纬差分割成许多个小梯形图块AEiFiA1，用公式（2）计算出各小梯形图块AEiFiA1的面积Si，然后累加Si就得到梯形图块ABB1A1的面积，同理，依次计算出梯形图块BCC1B1、CDD1C1、DAA1D1的面积（注：用公式（2）计算面积时，B1、B2分别取沿界址点编号方向的前一个、后一个界址点的大地纬度，ΔL为沿界址点编号方向的前一个、后一个界址点的大地经度的平均值与L0的差）；
（5）多边形ABCD的面积就等于4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）面积的代数和。

图1 椭球面上任意多边形计算面积
则任意多边形ABCD的面积P为：
P=ABCD= BCC1B1+ CDD1C1+ DAA1D1- ABB1A1
2、计算要求
① 利用图形坐标点将高斯坐标系下的几何图形反算投影到大地坐标系，进行投影变换。
② 任意指定一条经线L0，从选定多边形几何形状的起始点开始，沿顺时针方向依次计算相邻两点构成的线段，以及两点到指定经线的平行线构成的梯形面积。将该梯形沿纬度变化方向（Y轴）进行切割，至少需切割为2个部分。
③ 计算过程中应顺同一方向依坐标点逐个计算相邻两点连线与任意经线构成的梯形面积，坐标点不得有遗漏。若多边形包含内多边形（洞），则该多边形面积为外多边形面积减去所有内多边形面积之和。
④ 计算所有梯形面积的代数和即为该多边形的面积。
七、算法伪代码描述
为了确保编程使用的参数、算法一致，保证不同软件计算的椭球面积一致，我们用算法伪代码描述的方法对编程进行统一，在利用计算机编制椭球面积计算软件时，计算参数与计算顺序应严格按照以下代码执行。
1、参数说明
双精度类型：
圆周率值：PI = 3.14159265358979
中央经线:CenterL
RHO = 206264.8062471

A:ParamA
B:ParamB
C:ParamC
D:ParamD
E:ParamE

Const ZERO As Double = 0.000000000001

80椭球常数

椭球长半轴：aRadius = 6378140
椭球短半轴：bRadius = 6356755.29
椭球扁率：ParaAF = 1/ 298.257
椭球第一偏心率：ParaE1 = 6.69438499958795E-03
椭球第二偏心率：ParaE2 = 6.73950181947292E-03
极点子午圈曲率半径：ParaC = 6399596.65198801

k0:Parak0 = 1.57048687472752E-07
k1:Parak1 = 5.05250559291393E-03
k2:Parak2 = 2.98473350966158E-05
k3:Parak3 = 2.41627215981336E-07
k4:Parak4 = 2.22241909461273E-09

2、算法描述

初始化参数

Double e;
Double a;

e = ParaE2;
ParaC = aRadius / (1 - ParaAF);

ParamA = 1 + (3 / 6) * e + (30 / 80) * Power(e, 2) + (35 / 112) * Power(e, 3) + (630 / 2304) * Power(e, 4);

ParamB = (1 / 6) * e + (15 / 80) * Power(e, 2) + (21 / 112) * Power(e, 3) + (420 / 2304) * Power(e, 4);

ParamC = (3 / 80) * Power(e, 2) + (7 / 112) * Power(e, 3) + (180 / 2304) * Power(e, 4);
ParamD = (1 / 112) * Power(e, 3) + (45 / 2304) * Power(e, 4);

ParamE = (5 / 2304) * Power(e, 4);

参数初始化结束

中央经线转换为弧度
CenterL = TransDegreeToArc(CenterL)

选定本初子午线为参考经线
StandardLat = 0

For 起始点 To 倒数第二点

由高斯坐标反解计算经纬度值
ComputeXYGeo (PntColl.Point(i).y, PntColl.Point(i).x, B, L, CenterL)
ComputeXYGeo (PntColl.Point(i + 1).y, PntColl.Point(i + 1).x, B1, L1, CenterL)
将经纬度转换为弧度值
B = B / RHO
L = L / RHO
B1 = B1 / RHO
L1 = L1 / RHO

计算梯形面积
Double AreaVal;//梯形面积值
Double lDiference ;//经差
Double bDiference; //纬差
Double bSum;//纬度和
Double ItemValue(5);//计算变量

bDiference = (B1 - B0)；

bSum = (B1 + B0) / 2；

lDiference = (L1 + L) / 2;

ItemValue(0) = ParamA * Sin(bDiference / 2) * Cos(bSum);
ItemValue(1) = ParamB * Sin(3 * bDiference / 2) * Cos(3 * bSum);
ItemValue(2) = ParamC * Sin(5 * bDiference / 2) * Cos(5 * bSum);
ItemValue(3) = ParamD * Sin(7 * bDiference / 2) * Cos(7 * bSum);
ItemValue(4) = ParamE * Sin(9 * bDiference / 2) * Cos(9 * bSum);
AreaVal = 2 * bRadius * lDiference * bRadius * (ItemValue(0) - ItemValue(1) + ItemValue(2) - ItemValue(3) + ItemValue(4));

areaSum = areaSum + AreaVal;
Next

End Sub

3、高斯坐标反解算法

Public Sub ComputeXYGeo(x As Double, y As Double, B As Double, L As Double, center As Double)

Dim y1 As Double
Dim bf As Double

y1 = y - 500000

Dim e As Double

e = Parak0 * x

Dim se As Double

se = Sin(e)
bf = e + Cos(e) * (Parak1 * se - Parak2 * Power(se, 3) + Parak3 * Power(se, 5) - Parak4 * Power(se, 7))

Dim v As Double
Dim t As Double
Dim N As Double
Dim nl As Double
Dim vt As Double
Dim yn As Double
Dim t2 As Double
Dim g As Double

g = 1

t = Tan(bf)
nl = ParaE1 * Power(Cos(bf), 2)
v = Sqr(1 + nl)
N = ParaC / v
yn = y1 / N
vt = Power(v, 2) * t
t2 = Power(t, 2)
B = bf - vt * Power(yn, 2) / 2 + (5 + 3 * t2 + nl - 9 * nl * t2) * vt * Power(yn, 4) / 24 - (61 + 90 * t2 + 45 * Power(t2, 2)) * vt * Power(yn, 6) / 720

B = TransArcToDegree(B)

Dim cbf As Double

cbf = 1 / Cos(bf)
L = cbf * yn - (1 + 2 * t2 + nl) * cbf * Power(yn, 3) / 6 + (5 + 28 * t2 + 24 * Power(t2, 2) + 6 * nl + 8 * nl * t2) * cbf * Power(yn, 5) / 120 + center
L = TransArcToDegree(L)
End Sub

弧度转换为度
Public Function TransArcToDegree(arc As Double) As Double
Dim degree As Double
Dim min As Double
Dim sec As Double
Dim ret As Double
Dim tmp As Double
ret = arc * 180 / PI
degree = FormatValue(ret, 100, 100)
tmp = (ret - degree) * 60
min = FormatValue(tmp, 100, 100)
sec = (tmp - min) * 60
//秒保留到小数点后6位，四舍五入
sec = Format(sec, "####.000000") 'FormatValue(sec, 10000000, 100)
TransArcToDegree = degree * 3600 + min * 60 + sec
End Function

Private Function FormatValue(inputVal As Double, precsion As Long, scaleNum As Long) As Double
FormatValue = (Int(inputVal * precsion) - Int(inputVal * precsion) Mod scaleNum) / precsion
End Function