绍兴市人民政府办公室关于转发市安监局绍兴市重大安全生产事故隐患整治监督管理办法的通知
浙江省绍兴市人民政府办公室
绍兴市人民政府办公室关于转发市安监局绍兴市重大安全生产事故隐患整治监督管理办法的通知
绍政办发〔2005〕140号
各县(市、区)人民政府,市政府各部门:
市安监局制订的《绍兴市重大安全生产事故隐患整治监督管理办法》已经市政府同意,现转发给你们,请认真贯彻执行。
绍兴市人民政府办公室
二○○五年十一月八日
绍兴市重大安全生产事故隐患整治监督管理办法
绍兴市安全生产监督管理局
(2005年9月14日)
第一条 为实施对重大安全生产事故隐患的有效监控,督促生产经营单位进行整治,防止重大事故的发生,根据国家有关法律、法规,制定本办法。
第二条 本办法适用于本市行政区域内的各类生产经营单位。
第三条 重大安全生产事故隐患是指生产经营作业场所、设备及设施的不安全状态,人的不安全行为和管理上的缺陷,可导致重大事故发生的隐患,即可导致群死群伤或重大经济损失的隐患。
第四条 重大安全生产事故隐患整治监督管理原则:
(一)生产经营单位全面负责;
(二)各级政府属地管理,分级监管、督办;
(三)确定重点隐患,限期完成整治。
第五条 重大安全生产事故隐患整治监督管理的职责划分:
(一)各生产经营单位作为安全生产责任主体,对本单位存在的重大安全生产事故隐患排查、监控和整治全面负责,发现重大安全生产事故隐患应当及时报告主管部门,单位主要负责人是隐患整治的第一责任人。
(二)各级人民政府对本行政区域内重大安全生产事故隐患整治监督管理全面负责。
各级人民政府安全生产综合监督管理部门对行政区域内的重大安全生产事故隐患整治实施综合监督管理。
各级人民政府负有行业安全生产监督管理职责的部门(以下简称行业主管部门),对本行业生产经营单位的重大安全生产事故隐患整治监督管理负责。
(三)社区居民委员会、村民委员会发现其辖区内生产经营单位存在事故隐患或有安全生产违法行为时,应当及时向当地人民政府(街道办事处)或有关部门报告,并配合当地人民政府(街道办事处)或有关部门对辖区内存在的安全生产事故隐患进行整治。
第六条 各级人民政府应当加强对行政区域内重大安全生产事故隐患整治监督管理工作的领导,组织有关生产经营单位对可能发生重大安全生产事故的隐患进行排查、登记、建档、整治,并对下一级人民政府重点监督整治的重大事故隐患整治工作进行监督检查。
第七条 市、县(市、区)人民政府每半年分别组织有关部门对本级重点监督管理的重大安全生产事故隐患整治情况进行监督检查;乡镇人民政府街道办事处每季度对本级重点监督管理的重大安全生产事故隐患整治情况进行监督检查。各级开展监督检查情况应当及时书面报上一级人民政府及安全生产综合监督管理部门。
各级人民政府确定的重点监督整治隐患项目情况、隐患整治情况及验收合格情况,于每年1月底和7月底前报上一级人民政府及安全生产综合监督管理部门。
各级人民政府对重点监控整治的重大安全生产隐患,要明确事故隐患整治内容、整治责任单位和责任人、监督单位和责任人及完成整治时间。
第八条 市、县市、区人民政府有关行业主管部门应当对本行业内重点监控整治的重大安全生产事故隐患下达整治通知书,并对隐患整治进行监督检查,协调处理隐患整治中的重大问题。重大安全生产事故隐患排除前或排除过程中,无法保证安全的,要依法责令暂时停产停业整顿或者停止使用。
有关行业主管部门下达重大安全生产事故隐患通知书,应当抄送本级安全生产综合监督管理部门。
重大安全生产事故隐患整治结束后,事故隐患整治责任单位应当及时向履行该隐患整治监督管理职责的行业主管部门申请验收,有关行业主管部门应当会同本级安全生产综合监督管理部门及时组织有关专家,对事故隐患整治完成情况进行验收,验收合格的,应当及时予以书面确认并注销;验收不合格的,责令重新制定整治方案,继续整治。
第九条 重大安全生产事故隐患整治通知书应当包括下列内容:
(一)隐患的基本情况;
(二)隐患的类别;
(三)隐患整治要求和整治期限;
(四)隐患整治责任单位及责任人;
(五)隐患整治监管、督办单位及责任人。
第十条 确定列入重点监控整治重大安全生产事故隐患的单位,必须按照隐患整治通知书的要求,制定隐患整治实施方案,报经履行该隐患监控职责的行业主管部门及本级安全生产综合监督管理部门审查同意后组织实施。
重大安全生产事故隐患整治方案应包括以下内容:
(一)事故隐患类别;
(二)影响范围和危害程度;
(三)整治内容、措施和目标;
(四)整治资金来源及其保障措施;
(五)实施整治方案的时间安排及人员组成。
第十一条 存在重大安全生产事故隐患的单位在组织实施隐患整治期间,应成立由单位主要负责人负责的隐患整治管理小组,并履行以下职责:
(一)掌握本单位重大安全生产事故隐患的分布情况及发生事故的可能性及其程度,负责重大安全事故隐患的现场管理;
(二)制订应急预案,并报当地人民政府及有关部门备案;
(三)进行安全教育,提高职工的安全防范和自我保护能力;
(四)随时掌握重大事故隐患的动态变化;
(五)定期组织救援设施、设备调配和人员疏散演习,确保救援装备、器材完好有效。
第十二条 各级人民政府和有关部门要将重点监控整治的重大安全生产事故隐患向社会公布,并定期发布事故隐患整治信息,接受社会监督。
履行重大安全生产事故隐患整治监督管理情况,列入各级人民政府及安全生产综合监督管理部门和有关行业主管部门年度安全生产工作目标责任制考核内容。
第十三条 县市、区人民政府和有关部门对本行政区域内存在的重大安全生产事故隐患,超出其管辖或者职责范围的,应当立即向有管辖权或负有职责的上级人民政府或部门报告;情况紧急的,可以依法采取包括责令暂时停产停业或者停止使用等措施;上级人民政府或有关部门接到报告后,应当立即组织查处。
第十四条 各级安全生产综合监督管理部门应采取措施,确保有效实施本办法各项规定,并定期向本级政府报告本办法实施情况。
第十五条 本办法由市安全生产监督管理局负责解释。本办法自发布之日起施行。
关于统一图幅理论面积与图斑椭球面积计算要求的通知
国土资源部
关于统一图幅理论面积与图斑椭球面积计算要求的通知
(国土调查办发〔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