手持GPS坐标系转换方法.pdf

第2 5 卷第 1 期 镑 他鹭 GEOLOGY OF S HAANXI 2 o o 7 年 6 月 文章编号1 0 0 1 --6 9 9 6 2 0 0 7 O 1 一O 0 9 6 一O 7 手持 GP S坐标系转换方法 杜大彬，张宽房，张开盾，李明贵 陕西省地质调查 院，西安7 1 0 0 5 8 摘要 导航型手持 GP S目前在 中小比例地质调查等领域得到广 泛应用，由于坐 标系之I．-1 存在差异，在实际应用过程中，必须将手持机 的 WG S 8 4坐标 系转换为我 国应用的 B J 5 4或西安 8 O坐标系。坐标转换 的准确与否，直接影响到工程测量定位 的精度 ，传统的坐标转换计算所需要的起算资料不易收集，计算过程过于繁琐 ，非 专业人 员难 以掌握。本文根据收集的三角点 B J 5 4坐标 或西安 8 O坐标 ，和现场 测定的过渡坐标，求出各参数在本工作地 区的变化率，建立参数方程 ，反 向求出适 合于当地 的各项改正参数，方法简便易行，为手持 GP S定位的坐标转换方法提出 一 种 新 的思路 。 关 键 词 坐标 转换 ；WGS 8 4坐标 系 ；B J 5 4坐标 系；过渡 坐标 ；变 化率 中图分 类号 P 2 2 8 ． 4 文献 标识 码 B 随着技术的不断完善，导航型 GP S的定位精度及功能较之以前有很大提高。它以其全 天候工作 、携带方便 、数据记录及 回放快捷等功能 ，倍受使用者青睐。经过参数校正后 的 GP S ，其平面精度完全可以取代地形图定点 ，因而在 中小 比例尺地质矿产调查数字填图、地 球物理、地球化学勘探野外作业的点位测量 中有着广泛的应用前景。 坐标 系转换 问题提出 由于 GP S卫星星历是以 WG S 8 4坐标系 经纬度坐标为依据而建立的 ，我 国目前应 用的地形图一般采用 1 9 5 4年北京坐标 以下简称 B J 5 4坐标系或西安 8 O大地坐标系，不 同的坐标系之间存在平移和旋转关系，在不同地区，同一点位的 WGS 8 4坐标值与我国应用 的坐标系的坐标值 ，有约 6 O ～1 5 0 I n的差值 。在实际应用 中，不 同的坐标 系必须进行坐标 转换 。由于手持机测量通常是短时间近似测量 ，采用单次测量或多次测量值取平均值，一般 不作差分处理 ，从某种意义上讲，手持机的相对定位精度受其接收信号强度影响，坐标转换 参数的准确与否，直接影响其绝对定位精度。 坐标转换的关键是求出不 同坐标系之间的坐标转换参数 ，在实际工作过程中，坐标系统 收稿 日期 2 O O 7 一O 5 一O 8 作者简介 杜大彬 ，男 ，3 7岁，工程测绘工程师 ，主要从事物化探及地质测量工作。 维普资讯 第 1 期 杜 大彬等 手持 GP S坐标系转换方法 的转换通常采用方法是在应用区域内 G P S“ B ”级 网内，收集三个 以上网点 的 WG S 8 4坐标 系 B、L、H值及我国坐标系 B J 5 4或西安 8 0 B 、L、h 、x 高程异常 ，按其参考球体的 投影方式，计算各参数的差值。由于各地 GP S建网及重力研究工作程度不 同，通常在某些 地区 ，常用参数尤其是高程异常，一般不易收集 ，并且其计算过程较为繁琐。 为了寻求一种快捷 、方便、精度满足工作要求的 G P S坐标转换方法 ，作者经反复试验 ， 总结出坐标转换的一些规律 。以台湾 GAR MI N仪器公司的 E T RE X VI S TA 展望机型使 用为例 ，这里给出一种只用一个三角点，推算其 B J 5 4 西安 8 0 坐标改正参数的方法 。 2 参数变化在坐标系转换的规律 作者曾在陕南某地从事物探 电法工作 时，特 意在一 已知三角点作 GP S参数变化试验 ， 、 该三角点的 B J 5 4坐标值 为XXX X 0 4 3 3 ． 2 1 7 ；YXX X 6 7 6 0 5 ． 1 1 0 ，三角 点位于汉 江南 岸，视野开阔，有利 于 G P S观测。在观测 时设 置 当地 中央经 线、DA、D F等参 数，DX、 DY、DZ均为 0 ，在星况稳定且仪器显示估计误差为 5 m 时，在 已知点上读取若干组数据 ， 取得其平均值为 xX X X0 4 4 5 ；YX X X 6 4 4 。此值作为 WGS 8 4与 B J 5 4坐标系之间转换 的 过渡 坐标 。 如下表所示，利用 GP S自带软件‘ MAP S OUR C E’ ，分别变化不 同的参数 ，研究 D X、 D Y、D Z变化时，过渡坐标 X 0 、Y 0的变化规律。 表 1 同点位变化不同坐标参数过渡坐标 变化规律 统计 表 Ta b． 1 Cha n ge s i n t r a n s i t i o na l c o o r di n a t e s o f di f f e r e nt c o o r di n a t e p a r a me t e r s i n d i f f e r e n t po i n t s X0 Y0 DX DY DZ X1 Y1 △ X △ Y XXX0 4 4 5 XXX64 4 1 O O O O O XXXO2 5 O XXX5 8 4 1 9 5 94 0 XXX0 4 4 5 XXX6 4 4 1 0 O 0 O O XXX06 4 0 XXX7 0 4 1 9 5 9 4 0 XXX0 4 4 5 XXX6 4 4 O 1 0 0 0 O XXXO9 5 2 XXX9 8 6 5 O7 3 4 2 XXX0 4 4 5 XXX6 4 4 O 一 1 0 O 0 O XXX9 9 3 8 XXX3 0 2 5 O7 3 4 2 XXX0 4 4 5 XXX6 4 4 O O ，1 0 O 0 XXX9 6 0 5 XXX6 3 3 8 4 0 1 1 XXX0 4 4 5 XXX6 4 4 O O 一 1 O O O XXX1 3 8 5 XXX6 5 5 8 4 0 1 1 为便于理解，这里将 D X、D Y、D Z分 别设 置为 0时现场测定 的坐标定义 为 X 0 、Y0 ； X 1 、Y1为相应 DX、DY、D Z变化后对应于过渡坐标 X 0 、Y0变化后 的值 ，△X X 1 一X0 ， △YYI Y 0 ；X、Y的变化率为 △X ／ DX*1 0 0 、△Y／ DX*1 0 0 、△X ／ DY*1 0 0 、 △ Y ／ DY*1 0 0 、△ X ／ D Z*1 0 0 、△ Y ／ D Z*1 0 0 。 一 通过上表数据，我们可以看出 D X、D Y、D Z在各方 向上的变化，相应影响过渡坐标 的变化均呈线性，并有以下规律 1 当 DY、D Z不变时，D X十_ X1 ，Y1十；DX _ X1 十，Y1 ；x、Y的变化 率分 别为 一1 9 ． 5 、 9 4 ． 0 维普资讯 9 8 陕西地质 第 2 5 卷 2 当 DX、D Z不变时，D Y十 - - - Xl十，Y1 十；DY ---- -X1 J ， ，Y1 J ， ；x、Y的变化率 分别为 5 0 ． 7 、3 4 ． 2 3 D X、D Y不变时，DZ十---- -X 1 J ， ，Y 1 一 ；DZ J ， ---- X l十，Y 1 一 ；X、Y的变化率分 别为--8 4 ． O 、一1 ． 1 以上 十表示增大 ，J ， 表示减小，一表示基本不变化。 3 快速确定转换参数的数学基础 由于上述各参数的改变 ，不影响高程测量，因此，根据 以上规律 ，可列以下方程 r 一 0 ． 1 9 5 DX 0 ．5 0 7 DY 一 0 ． 84 DZ 一 1 1 ． 8 1 0 ． 9 4 DX 0 ． 3 4 2 DY 一 0 ． 0 1 1 DZ 一一 3 8 ． 9 X5 4 一 X过 渡 Y 5 4 一 y 过 菠 由此可见，D X、DY、D Z三参数是一个多解方程，在本例中任意求得几组参数 ，例如 ① DX一7 4 ． 9、D Y9 5 、DZ 8 8 ． 8 ；② D X--2 8 ． 9、D Y一 ～3 4 ． 4 、DZ 0 ；③ D X 一 1 7 ． 9、D Y---- - 6 5 ． 3 、D Z --2 1 ． 2 ；用三台手持机分别置人以上三组参数 ，在同一地区多个 已知点上 测试 ，测量结果 全部正 确 。 在实际工作中，根据以上规律 ，我们可以通过计算过渡坐标与 B J 5 4 或西安 8 0 坐标 的差值，除以各方向的变化率 ，按以上规律 ，列出求解方程 ，通过求解 ，可快速得 出以上各 参 数值 。 由以上参数变化规律，也可 以用趋近法快速求出各 自参数。因 DZ的变化对于 Y值几乎 不构成影响，可以先根据△x，确定 DX大致值，然后根据△Y的值确定 DY的值 ，最后再 根据△x精确确定 D z的值。 4 应用实例 作者近年一直从事地球物理勘探工作 ，曾在不同地区多个项 目中使用手持 G P S进行测 量定位 ，物探测量作业方法，要求用手持 GP S按设计的测量线路 ，在现场对每一个物理测 点进行实地定位。为保持与工作设计一致，必须对工作中使用 的 G P S进行坐标转换 ，使之 与工作用图相统一 。以上方法在不同地 区经实验均得以验证。以三个项 目 G P S坐标转换为 例，特将此方法在工作中的应用过程作具体介绍。 4 ． 1 实例一 内蒙某工作区手持 G P S坐标转换参数求取 该工作区位于内蒙古 自治区中部，共 4 个 1 5万图幅，处于 1 9 度带，工作用地形图和 收集三角点成果均为 B J 5 4坐标系，测网布设为规则测网方式，测网网度 5 0 0 mX 1 0 0 m，磁 测定位所用 G P S型号为E T R E X VI S T A，标本采集用 G P S型号为 GP S 7 2 ，收集三角点资 料 为 2 0 0 6年 5月 1 0日上午 ，校 正三 角点 时，在 三角 点‘ X 2 ’上 测得 过 渡坐标 ，X 0 维普资讯 第 1 期 杜 大彬 等 手持 GP S坐标 系转换 方法 9 9 4 6 2 7 5 4 0 、YO 一4 0 5 3 5 9 。求 出各 自变 化率后 ， 由以上方 法列 方程 表 2 内蒙某工作 区四幅 1 5万 地质 矿产调查 工作区收 集三 角点坐标 Ta b ． 2 Th e t r i g o no me t r i c c o or d i na t e s a c q ui s i t e d f r o m g e o l o gi c a l s ur v e y i n g a t 15 0， 0 00 I nn e r M o go l i a 点 名 等 级 X Y H X1 I V XXX6 9 6 7 ． 5 9 XXx8 5 3 7 O ． 7 7 1 7 O 2 ． 4 X2 Ⅲ XXX7 4 9 9 ． 7 9 XXX0 5 3 1 8 ． 9 9 1 6 0 9． 6 X3 Ⅲ XXX0 6 1 3 ． 8 2 XXX1 2 9 21 ． 7 2 1 5 7 7 ． O ， 一 0 ． 2 3 9 DX 4 - 0 6 2 2 DY一 0 ． 7 4 6 DZ 一一 4 0 ． 2 Xs 4一 X 1 0 ． 9 3 7 D X 4-0 ． 3 4 8 D Y一0 ． 0 1 0 D Z一一4 0 ． 0 y 5 4 一y 求得① D X一5 1 ，D Y 2 4 ． 9 ，D Z 9 1 ；② D X一4 2 ，DY 0 ，D Z 6 7 ． 3 ；③ D X 一0 ，DY一 一 l 1 6 ，D Z一 一 4 2 ． 9分 别 置 人 各 组 参 数 后 ，现 场 测 得 X X XX 7 5 0 0 ，Y XXX31 9。 按照‘ X1 ’三角 点的坐标 导 航 ，在 坐标 处 地 表见 一 约 近 1 m 深 坑 ，中心 位置 与 该 点坐 标相一致 ，深挖后未见任何标记，在坑东旁约 1 0 m处有一敖包，仔细查找后 ，发现该点标 石被 当地牧 民挖 出后 垒在 敖包 中 。 当日下 午，利 用 x3三 角 点 的 坐 标 导 航 ，找 出‘ x 3’后 ，测 得 该 点 坐 标 为 X X X X 0 6 1 5 ，YX X X9 2 2，和 收集 资料 基本 一致 。 4 ． 2实例题 新疆 A 工作 区手持 GP S坐标 转换 参数 求取 该工区位于新疆西部，东西向两个 1 5万图幅，处于 1 3 度带，所用工作布置图及收集 三角点成果均为西安 8 O坐标系。测 网布设 为半 自由网方式 。磁测定位和标本采集采用 GP S 型号均为 E TR E X VI S T A。收集 三角点 成果 如下 表 表 3 新疆 A 工区 1 5万地质 矿产 调查 收集三角点坐 标 Ta b ．3 Th e t r i g on o m e t r i c c o or d i na t e s a c q ui s i t e d f r o m ge o l og i c a l s ur v e y i n g a t 15 0， 0 0 0 i n r e g i o n A ， Si n ki a ng 点名 等级 x I Y H I x 1 Ⅲ X X X 2 3 9 8 ． 9 3 l X X X 9 2 1 2 8 ． 0 6 3 6 3 7 ． 1 O I X 2 Ⅱ X X X 6 6 6 2 7 ． 7 0 I X X X 1 2 3 3 2 ． 4 4 2 5 2 7 ． 1 6 2 0 0 6 年 6月 2 7日上午 ，参数 改正前 ，分两 个作业组 在三角点‘ x 2 ，和‘ x1 ，上作 G P S观测。‘ x 2 ’在测得该点的过渡坐标 X O 4 5 1 6 6 7 8 、Y0 7 1 2 4 1 0 。由以上方法求各项 参数的变化率，列出以下方程 r 0 ． 1 6 9 DX 0 。 6 3 2 DY 一 0 ． 7 5 8 DZ 一 5 0 。 3 1 0 ． 9 7 2 DX 0 ． 2 3 5 DY 0 ． 0 2 2 DZ 一 7 7 ． 6 X8 o x过 渡 y8 o Y过 渡 维普资讯 i 0 0 陕西地质 第 2 5卷 求得① DX一一1 0 7 ． 3 ，D Y一 一1 1 9 ，D Z ----一5 6 ． 8 ；② D X一一8 O ． 9 ，D Y一 0 ，D Z一 4 8 ． 3 ；③D XO ，D Y----3 6 4 ． 8 ，D Z ----3 7 0 ． 7 分别置入各组参数后，现场测得 XX XX 6 6 6 2 8 ， y_--X X X3 3 2 。另一组测 得 ‘ X1 ’三角 点的测 量过 渡坐标 为 X_--X X X 2 4 5 0 ，y_--X X X 2 0 6 ，经 软件改正计算后 ，XX X X 2 3 9 9 、YXX Xl 1 2 7 。与收集资料吻合。 4 ． 3实例 3 新疆 B工 区手持 G P S坐标 转换 参数 求取 本工区位于 A工区以北并南北相接，为同地区跨年度不同项 目，如上同样为东 西向共 两个 1 5万图幅，处于 1 3度带 ，所用工作布置图以及所 收集 的三角点成果均为 B J 5 4坐标 系 ，测网布设为半 自由网方式。磁测定位和标本采集采用 GP S型号均为 E T RE X VI S T A。 收集三角点成果为 表 4新疆 B工 区 1 5万地 质矿产调查收集三角点坐标 Ta b ． 4 The t r i go n o m e t r i c c o o r di na t e s a c q ui s i t e d f r o m g e o l o gi c a l s ur ve y i ng a t 15 0， 0 0 0 i n r e g i o n B，Si n ki a ng 点 名 等 级 X Y H X1 Ⅲ X Xx7 8 4 3 ． 1 6 XX X2 3 0 3 1 ． 2 4 3 1 O 7 ． 0 X 2 Ⅱ XXX 6 9 7 4 ． 1 3 X XX O 9 1 6 4 ． 6 9 3 4 7 6 ． 5 X 3 Ⅱ XXX 8 1 2 6 ． 6 5 X XX 2 5 2 7 9 ． 4 4 3 2 1 0 ． 3 2 0 0 6 年 6月 2 5日，分三个小组寻找 ‘ X 1 ’ 、‘ X 2 ’及 ‘ X 3 ’三个三角点 ，并作 G P S观 测 ，第一组在 ‘ X1 ’测得过渡坐标 ，X 0 一X XX 7 8 8 0 、YO XX X 9 7 7 。由以上方法列方程 r O ．1 69 DX 0 ．6 3 3 DY 一 0 ． 7 57 DZ 一 一 3 6 ．8 1 0 ． 9 7 3 DX 一 0 ． 2 3 3 Dy 0 ． 0 2 3 DZ 一 5 4 ． 2 X5 4一 X过渡 y5 4一 y过 渡 求 得① DX 4 3 ． 2 ，D Y 一5 0 ． 7 ，D Z 1 5 ． 9 ；② D X3 9 ． 3 ，DY---- 一6 8 ． 6 ，DZ 0 ； ③ D X5 4 ． 3 ，DY----O ，D Z 6 0 ． 7分别置人各组参数后 ，经转换计算后 X XX X 7 8 4 3 ，Y XX X 0 3 1 。另一组测得 ‘ 库尔嘎克 ’三角点的测量过渡坐标为 XXX X 7 0 1 1 ，YX X X1 1 1 ， 经软件改正计算后 ，XXX X 6 9 7 4 、YX XX 1 6 5 。与收集资料一致 。‘ X3 ’被破坏。 5 注意事项 以上坐标转换方法是在 1 5万工作 比例尺，6度带投影方式下进行 的，但此方法同样 适用于 3度带投影的坐标转换 。此方法的关键是求出各参数变化率。 坐标转换的前提是必须正确设置位置格式参数和地 图基准参数 。手持机一般默认为经纬 度坐标显示格式，在位置显示格式菜单 中，选择“ US E R UTM GR I D” ，在输人参数 页面 中，分别输入相关参数。其中中央经线为工作 区域坐标分带的中央子午线经度 ，投影 比例为 维普资讯 第 1期 杜 大彬 等手持 GP S坐标系 转换方法 1 O 1 1 ，东西偏差为 5 0 0 0 0 0 m，南北偏差为 0 m。 在地图基准菜单 中，选择用户定义模式 “ US E R” ，在输入参数页面中，共有 D X、DY、 D Z 、DA、DF等 5项参数 。转换 目标坐标为 B J 5 4坐标 系时， DA一一1 0 8 ，DF 0 ． 0 0 0 0 0 0 5 为 固定常量 。转换 目标坐 标为 西安 8 O坐标 系时 ，DA 一3 ，D F0 ． 0 0 0 0 0 0 0 0 3 。 按照部分 G P S使用说明书介绍，在同一地区完成 GP S参数校正后 ，选择 5个以上 的三 角点验证 ，当最大误差不 大于 1 5 m 时，即可使 用。笔者认为这是一种不严密说法 ，由于 WG S 8 4系统所采用的球体模 型与我国使用的 B J 5 4及西安 8 O系统存 在差异 ，在不同地 区甚 至同一地区跨度大时均需作参数 改正。在 同一地 区，当南北两三角点跨度 大于 2 0 0 k m 时， 由两已知点坐标换算所影响的计算误差就有可能超过 1 0 m，根据经验 ，在不同方向距离每 变化 2 0 k m，GP S坐标换算本身大约影响 1 m。 实际上 G P S定位的精度由相对误差和绝对误差构成 ，相对误差的精度取决于测量时刻 卫星的星历以及接收机所处位置的 G P S信号 的接收情况影 响。绝对误差 与测量过程无关 ， 转换参数的正确与否，直接影响到我们所使用的测量成果 。 由于 GP S的工作原 理 ，手持机 的定 位模 式 是 一种 未 经 差分 的测 量模 式 ，所有 测 量 工作 应在星况稳定时进行 。所以此方法适应于同一分带区东西、南北方向跨度不大时使用 ，一般 用一个三角点就可以完成校正工作，不需在其他三角点上验证，但在实际工作中，仍然提倡 在已知三角点或明显标志点上进行验证 。虽然以上方程参数存在多解 ，在实际应用中，尽可 能使用参数绝对值小的参数改正 ，并且校正参数时，所选择的三角点位置最好位于工区的中 部 ，如果同一工作区跨度超过一定范围时，在工区两侧地带采用不同的校正参数 ，可以提高 测量结果 的绝 对精 度 。 以上为作者在工作过程 中总结出的一些经验，由于时间及水平所限，错误或疏漏在所难 免，诚请广大同行批评指正 。 [ 参考文献] 略 CooRDI NATE CoNVERS I oN BY PALM ToP GPS Du Da b i n， Zha ng Ku a n f a ng， Zh a ng Ka i - d un a nd Li M i ng g ui Ge olo g i c a l S u r v e y o f S h a a n x i Pr o v i n c e ，Xi ’ a n 7 1 0 0 5 8 Ab s t r a c t Na v i ga t i o n t y pe p a l mt o p GPS h a s be e n wi d e l y a ppl i e d i n g e o l o gi c a l s u r v e y i n g a t s ma l l s c a l e ． W he n i t i s us e d i n g e o l o g i c a l s ur v e y i n g，t he W GS 8 4 s ys t e m o f pa l mt o p GPS ha s t o b e c o n v e r s e d i n t o B J 5 4 o r Xi ’a n 8 O s y s t e ms t h a t a r e u s e d i n Ch i n a ，s i n c e t h e d i s c r e p a n 维普资讯 1 O 2 陕西地质 第 2 5卷 c i e s i n di f f e r e nt c o o r d i n a t e s y s t e ms ． The p r e c i s i on o f c o or d i n a t e c o nv e r s i o n wo ul d d i r e c t l y i mpa c t t h e l o c a t i o n o f e ng i ne e r i n g s ur v e y ．Ac q u i s i t i o n o f b e g i ns r e c k on i n g ma t e r i a l s r e q ui r e d f o r c o nv e r s i o n c a l c u l a t i on of t r a d i t i o n a l c o o r di n a t e s i s no t e a s y a n d t h e pr o c e s s of c a l c ul a t i o n i s c o mp l i c a t e d a n d h a r d t o b e u s e d b y n o n p r o f e s s i o n a l s ．Th e t r i g o n o me t r i c B J 4 5 c o o r d i n a t e s o r t h e Xi ’ a n c o o r d i n a t e s a n d me a s u r e d t r a n s i t i o n a l c o o r d i n a t e s c a n b e u s e d t o r e c k o n t h e v a r i a t i o n s o f e a c h p a r a me t e r i n t hi s a r e a，t o c r e a t e a p a r a me t e r e qu a t i o n a nd r e v e r s e l y c a l c u l a t e t h e c or e c t e d p a r a me t e r s i n l o c al i t y ． The me t ho d i s s i mpl e a nd c o nv e ni e nt ，wh i c h i s ne w i n c on ve r s i on o f me n s u r a t e d c oo r d i n a t e s by p a l mt o p GPS ． Ke y wo r d s c o n v e r s i o n o f c o o r d i n a t e s ；B J 5 4 c o o r d i n a t e s y s t e m；t r a n s i t i o n a l c o o r d i n a t e s ； v a r i a t i o n 上接 第 8 4页 AS S ES S M ENT oF W ATER FI LLI NG FACToRS I N HENGGoU CoAL M I NE oF THE W UBU CoAL FI ELDS Ya n g Z h i g u o a n d L i a n g X i u j u a n 。 1 ． Th e Ge o l o g i c a l Te a m No ． 3 S h a a n x i Bu r e a u o f Ge o l o g y a n d Mi e n r a l Re s o u r c e s ，B a o j i 7 2 1 3 0 0； 2 ． I n v e s t i g a t i o n a n d S u r v e y i n g I n s t i t u t e o f Ba o j i Ci t y ，B a o j i 7 2 1 0 0 0 Ab s t r a c t The W u bu c o a l f i e l ds a r e a bun d a nt i n c oa l r e s ou r c e s ． M i ni n g c ond i t i o ns i n t he f i e l d s a r e c o mp l i c a t e d，i n wh i c h c o mp l e x c o n d i t i o n s o f g e o h y d r o l o g y i s t h e ma j o r f a c t o r t h a t i m p a c t c o a l mi ni ng i n t h i s c o a l f i e l d s ． Ac c o r di ng t o e x pl or a t i o n i n f o r ma t i o n o f t he H e ng g ou c o a l m i ne，t h i s p a pe r h a s ma de a l l a r o und a na l y s i s o f t h e g e oh y dr o l o gi c a l c o ndi t i o ns i n t he c o a l f i e l ds ，a s s e s s e d t he wa t e r f i l l i ng f a c t o r s i n c o a l m i ne s a nd r e a c he d t h e f o l l owi ng c o n c i d e r a t i o nswa t e r f r o m t he Ye l l o w Ri ve r i s no i m p a c t s t o c o a l mi ni ng i n t h e He ng go u c o a l mi ne s； t he Aoh ui c o nf i ne d wa t e r ma y no t c a u s e d a ma ge s t o m i ni ng t h e c o a l b e ds o f t h e Sh a nx i Fo r ma t i on i n no r ma l s e c t i on s，bu t wo ul d b e a s o ur c e o f da ng e r t o c o a l mi ni n g i n t he Ta i yu a n Fo r ma t i on ． I t i s r e q ui r e d t O d o f u r t he r de t a i l g e o l o gi c a l e x pl o r a t i o n f o r ma k i ng s ur e t he p os s i bi l i t y o f c oa l mi n i ng i n t he Ta i y u a n For m a t i o n． Ke y wo r d s W u bu；t h e He n gg o u c o a l f i e l d s；c o a l mi n e；wa t e r - f i l l i ng f a c t o r ；a s s e s s me n t 维普资讯