国家安全生产监督管理局 国家煤矿安全监察局


关于进一步加强煤矿安全生产工作的紧急通知


各产煤省、自治区、直辖市安全生产委员会和经贸、煤炭管理部门及安全生产监督管理部门，各省级煤矿安全监察局及北京、新疆生产建设兵团煤矿安全监察办事处：

进入二季度以来，各地区、各有关部门认真贯彻落实4月8日黄菊副总理在全国安全生产电视电话会议上的重要讲话和国务院办公厅《关于进一步加强安全生产工作紧急通知》精神，在组织开展煤矿安全生产大检查、深化煤矿安全整治、防范事故方面做了大量工作，取得了一定效果。但是，各地工作进展不平衡，煤矿重特大事故仍时有发生，全国煤矿安全生产形势依然严峻。特别是5月13日至5月24日的12天中，接连发生3起煤矿特大事故和2起特大事故险情：5月13日，安徽省淮北矿业集团公司芦岭煤矿（国有重点矿）发生瓦斯爆炸事故，死亡86人；5月19日，山西省晋中市寿阳县富东煤业有限公司（国有地方矿）发生透水事故，33人被困，后经全力抢救全部脱险；5月20日，山西省临汾地区安泽县永泰煤矿（无证非法矿）发生瓦斯爆炸事故，死亡25人；5月21日，云南省丽江地区华坪县永兴煤炭有限责任公司基佐煤矿（国有地方矿）发生瓦斯爆炸事故，死亡24人；5月24日，河南省安阳市龙安区安里煤矿（乡镇有证矿）发生透水事故，15人被困，目前正在全力抢救。特大事故接连发生，反映出一些地方人民政府安全生产监督管理工作不严，经贸、煤炭管理部门督促检查不力，一些煤炭企业未能正确处理安全生产与效益的关系，重生产轻安全问题比较突出，安全管理工作松懈，安全措施不落实。同时也反映出对煤矿安全执法力度不够、对事故查处不严的问题，在一些地方依然存在。党中央、国务院领导同志对此高度重视，相继作出重要批示，要求“对这一苗头和倾向，各地特别是重点产煤省区必须本着对人民生命负责的精神，采取一些强硬措施尽快加以制止；各级煤矿安全监察部门也要通过明察暗访、严格执法、严厉追究责任等措施，真正把监管责任落实下去。”为贯彻落实党中央、国务院领导同志的重要批示精神，遏制煤矿重、特大事故上升的势头，特提出以下要求：

一、认真贯彻落实党中央、国务院关于“坚持一手抓防治非典型肺炎这件大事，一手抓经济建设这个中心不动摇”的重大战略决策，在全力抗击“非典”的同时，继续扎实有效地抓好安全生产工作。各级安全生产监管部门、煤炭工业管理部门和煤矿安全监察机构，一方面，在当地政府领导下把防控“非典”作为一项重大任务，切实抓紧抓好；另一方面，决不能放松安全监督管理和煤矿安全监察，必须采取切实有效的措施，在解决实际问题上狠下功夫，切实把安全监管工作做到位，把煤矿安全监察工作做实，为夺取抗击“非典”和经济发展的双胜利，创造安全稳定的环境。

二、继续深入开展煤矿安全生产大检查。各省（区、市）安委会要组织经贸、煤炭管理部门、安全生产监管部门和煤矿安全监察机构，认真分析本地煤矿安全生产形势，查找存在的突出问题和薄弱环节，继续深入贯彻落实4月8日全国安全生产电视电话会议和国务院办公厅《关于进一步加强安全生产工作紧急通知》精神，按照5月21日召开的全国安全生产电话会议要求，继续搞好煤矿安全生产大检查。一是要求国有煤炭企业对本企业安全生产情况进行一次全面检查。特别要切实抓好“一通三防”和防治水工作，不折不扣地落实“先抽后采，监测监控，以风定产”十二字方针，坚持做到“六不准”：不准超通风能力生产；瓦斯超限不准作业；无风电、瓦斯电闭锁的不准掘进；高瓦斯和突出矿井无瓦斯抽放系统不准生产；采区接替不准违反有关规定，随意并巷减尺；从业人员未经培训合格、未佩带自救器的不准下井。二是要求各市（地）、县（市）人民政府对所属煤矿安全生产情况进行抽查，抽查覆盖面不低于50%。三是各省（区、市）安委会要组织有关部门组成煤矿安全生产督查组，对重点产煤市（地）、县（市）和重点煤炭企业进行安全生产督促检查，抽查覆盖面不低于30%。重点督查地方各级政府对煤矿安全监管的责任是否落实，措施是否到位，煤矿安全整治是否深入；是否存在已关闭非法煤矿死灰复燃问题；是否存在国有大矿超通风能力突击生产问题。发现问题，要立即采取强硬措施予以纠正；对非法生产的小煤矿，坚决依法取缔、并强制予以关闭。

三、对发生特大事故的市（地）坚决实施停产整顿。各省（区、市）安委会，要组织有关部门认真研究对煤矿安全生产问题突出、事故隐患严重的市（地）或县（市），采取停产整顿措施。特别是对今年4月1日以来发生10人以上特大事故（含3起特大事故险情）的市（地），要责令其辖区内所有煤矿停产整顿，认真整改事故隐患，经市（地）人民政府验收通过后，方可恢复生产。这些市（地）是：贵州安顺地区、贵州六盘水市、河北邯郸市、湖南娄底市、山西临汾地区、山西晋中市（发生了特大事故险情）、安徽淮北市、云南丽江地区、河南洛阳市（发生了特大事故险情）、河南安阳市（发生了特大事故险情）。国务院安全生产检查组将对这些市（地）停产整顿情况进行督促检查。

四、加大煤矿安全监察执法力度。各省级煤矿安全监察局除留1名负责人在机关主持日常工作外，主要负责人和领导班子其他成员一律带队深入下去进行煤矿安全监察，局机关公务员至少抽调一半深入一线进行煤矿安全监察；各煤矿安全监察办事处除留值班人员外，一律分片包干，深入煤矿进行执法监察。要重点监察“一通三防”、防治水和以风定产措施的落实情况，打击小煤矿非法生产措施落实情况。要通过明察暗访等形式，彻底排查事故隐患，特别是对无证非法开采的小煤矿，要依法予以取缔，加大经济处罚力度，没收全部非法所得，并处以1倍以上5倍以下罚款。对于监察所发现的问题，要依据《煤矿安全监察条例》及时下达监察执法文书，并通报当地人民政府予以落实，如果落实不下去，要逐级向上级地方政府通报，并依此作为追究地方政府和部门负责人行政责任的重要依据之一。

五、认真执行国务院第302号令，加大煤矿重、特大事故查处力度。对今年以来发生的重大、特大煤矿事故，凡已经结案的，要加快落实处理决定；尚未结案的，要严格按照“四不放过”原则，抓紧调查处理和结案工作，对于特大事故，要认真执行《国务院关于特大安全事故行政责任追究的规定》，严肃追究有关地方政府和部门负责人的行政责任。对经过两年煤矿安全专项整治，仍存在无证非法开采小煤矿的市（地）、县（市），也要严肃追究当地政府有关责任人的行政责任。

六、在煤矿安全生产大检查工作中，要切实落实有关防控“非典”的措施，保护安全监管监察人员的身体健康，确保各方面有关人员不被感染。

在国务院安全生产检查组赴各地检查之前，国家安全生产监督管理局（国家煤矿安全监察局）近日将组织煤矿安全督查组，由国家局负责人带队，分赴部分重点产煤省（区、市）进行煤矿安全专项督促检查。

请各地将贯彻落实此通知和“5.21”全国安全生产电话会议的工作情况，于6月15日前书面报告国家局。



国家安全生产监督管理局(国家煤矿安全监察局)
二ＯＯ三年五月二十五日



国土资源部


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

（国土调查办发〔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