第九章矿井定向的精度分析第一节用垂球线投点和投向的误差.ppt

上传人:本田雅阁 文档编号:2520471 上传时间:2019-04-05 格式:PPT 页数:92 大小:1.02MB
返回 下载 相关 举报
第九章矿井定向的精度分析第一节用垂球线投点和投向的误差.ppt_第1页
第1页 / 共92页
第九章矿井定向的精度分析第一节用垂球线投点和投向的误差.ppt_第2页
第2页 / 共92页
第九章矿井定向的精度分析第一节用垂球线投点和投向的误差.ppt_第3页
第3页 / 共92页
亲,该文档总共92页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

《第九章矿井定向的精度分析第一节用垂球线投点和投向的误差.ppt》由会员分享,可在线阅读,更多相关《第九章矿井定向的精度分析第一节用垂球线投点和投向的误差.ppt(92页珍藏版)》请在三一文库上搜索。

1、第九章 矿井定向的精度分析 第一节 用垂球线投点和投向的误差,根据煤矿测量规程要求,一井两次独立定向所算得的井下定向边的方位角之差,不应超过2,则一次定向的允许误差是 。若采用两倍中误差作为允许误差,则一次定向的中误差为: 此误差由三部分组成: 井上的连接误差m上; 投向误差; 井下的连接误差m下。,故M0=,一般情况下,一井定向的投向误差和连接误差大致相等。即m2上+m2下2,则投向误差不应大于下列数值: ,若井上与井下的连接误差相等时,则 m上=m下 下面根据上述精度要求,对用垂球线投点的投点误差、投向误差、一井定向和两井定向的误差加以分析。,一、 用垂球线投点的误差来源及估算方法 在井筒

2、中用垂球线投点的误差的主要来源: (1) 气流对垂球线和垂球的作用; (2) 滴水对垂球线的影响; (3) 钢丝的弹性作用; (4) 垂球线的摆动面和标尺面不平行; (5) 垂球线的附生摆动。 下面分别就上述各因素加以讨论。,(一) 气流对垂球线和垂球的作用 井筒中气流对垂球线的影响是十分复杂的,但又是一个很重要的问题。国内外一些矿山测量人员用试验观测的方法进行了不少研究工作。综合分析观测结果可得出如下结论: (1) 井筒中气流所引起的垂球线偏斜是投点误差的最主要来源,也是一井定向的最主要误差来源。 (2) 井筒中气流对垂球线的作用主要发生在马头门处(见图9-1),如对垂球线加防风套筒,可大大

3、减少风流的影响。,(3) 当井深为300-600m时,投点误差不超过1.5-2mm。 投点误差e可用下式进行计算: (9-4) 式中 p钢丝单位长度所受的侧压力; h马头门的高度; H井深; Q垂球的重量。 垂球线因受气流的影响所产生的偏斜值与垂球重量成反比,而与井深成正比。,(二) 井筒内滴水对垂球线的影响 井筒内的滴水、涌水或水管的漏水,将打击垂球线和垂球,破坏其均匀摆动的状态,但这些现象不可能用数学公式来表达。因此,在选择垂球线的悬挂位置时,应注意滴水的影响,并将垂球放入大水桶中稳定。,(三) 钢丝的弹性作用 钢丝弹性的影响表现在两个方面。一方面当缠在绞筒上的钢丝放入井内时,钢丝仍在企图

4、保持原来的环状。这样就使钢丝上各点偏离了其中心位置。为此,应采用直经大于250mm 的绞车、细的钢丝和适当的垂重,以减少其影响。,另一方面是当钢丝自滑轮经定点板放入井筒时,因定点板的中心不是恰好与滑轮槽位于同一铅直线上,故定点板与滑轮间这一段钢丝将成倾斜状态。由于钢丝的弹性,当经过定点板后,钢丝仍将有一小段斜向的位置,往下才逐渐被垂重拉直。为避免这种影响,应在定点板下方钢丝已完全铅直的部分进行地面连接测量。在布置滑轮与定点板时,应使两者间的一段倾斜线与铅垂线的交角尽可能小,同时两定点板应尽可能布置在两垂球线的连线方向上,以减少它对投向的影响。,(四) 垂球线的摆动面与标尺面不平行 当从经纬仪C

5、对垂球线的摆动极边位置L和R进行多次观测,在标尺MN上读取一系列读数l和r,然后取其平均值求得标尺上的读数时,则垂球线的摆动方向LR与标尺面MN平行和不平行时引起的差距aa0为: 式中 垂球线的摆幅(即LR); S经纬仪至标尺的距离; 垂球线摆动方向与标尺间的夹角。,(五) 垂球线的附生摆动 在理想的条件下,井筒内垂球线的摆动,应像钟摆一样具有均匀而逐渐衷减的摆动。但由大量的实际观测资料发现,垂球线各相邻摆幅的平均中点位置的连线,并没有成为一条直线,而是向左右偏移的曲线。 当垂球有了附生摆动后,按标尺读数所求得的平均位置,就不等于其真正的稳定位置,从而产生了投点误差。,产生附生摆动的原因: (

6、1)井筒内气流变化的影响; (2)滴水的打击; (3)气流对钢丝的摩擦作用; (4)地面垂球线固定点的振动; (5)钢丝的弹性。 减小措施:将垂球浸入稳定液中。,二、减少投点误差的措施 1.增大C,位置合理; 2.尽量减小马头门处气流的影响; 3.小直径、高强度钢丝,加大锤重Q,浸入稳定液; 4.摆动观测时, 垂球线摆动的方向尽量和标尺面平行; 5.减小滴水对钢丝及垂球的影响,加桶盖。,三、用重球线投向的误差 钢丝投点产生投点误差e, 一井定向时投二根钢丝, 产生投向误差,投点误差:风流、滴水等影响,钢丝地面井下投影不重合,线量偏差 投向误差:由投点误差所引起的垂球线连线的方向误差,a,A,B

7、,b,a,A,b,B,B,A,a,b,c,c,c,e,2A=(eA/c) 2/2 2B=(eB/c) 2/2 = e/ c,A0,B0,Bi,eB,i,i,c,第二节 三角形连接法的误差和有利形状,C,B,A,C,D,E,D,E,a,b,c,a,b,c,一、连接三角形中垂球线处角度的误差及三角形最有利形状 (1) 连接三角形最有利的形状为锐角不大于2的延伸三角形。 (2) 计算角和的误差,随测量角的误差(m只含测角方法误差)增大而增大,随比值ac的减小而减小。故在连接测量时,应使连接点C和C尽可能靠近最近的垂球线,并精确地测量角度。,(3) 两垂球线间的距离c越大,则计算角的误差越小。 (4)

8、 在延伸三角形时,量边误差对定向精度的影响较小。,二、连接角的误差对连接精度的影响 首先,讨论经纬仪在连接点C上的对中误差对连接精度的影响。 假设经纬仪在连接点C上的对中有线量误差eT而对中在C1点上,则连接边就成了C1 D。,因为在定向时,连接三角形的各测量元素(角和a、b、c边)都是根据经纬仪中心来测得的,所以仪器在C点的对中误差对连接三角形的解算没有影响,而只是对垂球线的方位角AB的确定有影响。当经纬仪对中无误差时,则 AB=DC+-2180,当经纬仪有对中误差时,则 由此而引起的确定方位角AB的误差为: 由图9-7可知: 故经纬仪对中不正确对 的影响为 。,故经纬仪对中不正确对 的影响

9、为 。由第七章式(7-21)可得中误差 (9-15) 由上式可知,连接边d越长,则此项误差就越小,它与CA的长短无关。,其次,在连接测量时,还要考虑到D点上的觇标对中误差meD,即 因此,在c点测连接角的误差,对连接精度的影响m为 式中 mi 测量方法误差。,当eT=eD=e1 时,则 (9-17) 由此可知,欲减少测量连接角的误差影响,主要应使连接边d尽可能长些,并提高仪器及觇标的对中精度。煤矿测量规程要求d尽量大于20m。 上述公式对估算井下连接测量的误差也同样适用。,三、 三角形连接法连接时一井定向的总误差 根据式(9-7)得定向总误差为: 式中各项误差的计算方法汇集如下: m2和m2一

10、样可用式(9-17)计算,即 投向误差可按式(9-6)计算,即,m(或m )在2,178的延伸三角形中可用式(9-14)计算,即 由于连接边的方位角 是由地面近井点设导线测出的,故 可按支导线的误差累积公式计算,即 式中 m地面近井导线的测角中误差; n 近井导线的角数。,四、 按正弦公式解算三角形时所用检查方法的可靠性 按正弦公式解算三角形时,曾用两种方法检查测量和计算的正确性。其一是对比两垂球线间距离的丈量值和计算值;其二是用三角形内角和是否等于180来检查。下面就分别讨论这两种检查方法的可靠程度。,(一) 两垂球线间距离检查的可靠性 若两垂球线间距离的丈量值为c,而计算值为c,则其差数d

11、=c-c的误差为: (9-18) 因,按前式取各偏导数,并令c=c后代入上式得 当用正弦公式解延伸三角形时,cos1,cos-1。将上式代入式(9-18)得 (9-20),上式等号右边三项为量边误差对差数d的影响,而最后一项为测角误差的影响。因在延伸三角形中,sin0,所以测角误差的影响反映不出来。因此,这种检查方法只能检查量边的正确性,而不能检查测角的正确性。,当三角形用正弦公式解算时,式(9-20)可近似写成为: 若 则 (9-21) 当ml=0.5mm时,md=0.531.0mm。取允许误差为中误差的二倍,则,煤矿测量规程规定,两垂球线间距离的丈量值与计算值之差,井上不应超过2mm。 考

12、虑到井下的工作条件较困难,故对井下放宽到不超过4 mm.,(二) 内角和检查的可靠性 三角形中三内角和数公式为 式中角度是实际测的,而及是按下式计算的: sin=a/csin sin=b/csin 因此,和数S是角度及边长a、b、c的实测值的函数。当测角量边均有误差时,则和数S的误差mS为:,将上列各偏导数值代入,则得 (9-23),上式等号右边第一项为量边误差对三内角和的影响,而第二项则为测角误差的影响。在延伸三角形中,sin0,tan及tan都近似等于零。所以三内角和不能检查量边的正确性,也不能检查测角的正确性。为此,现行规程建议在C点上测量、及三个角度(见图9-6),以资检查,三内角和可

13、以检查计算的正确性。,第三节 两井定向的误差 两井定向也和一井定向一样,是由投点、井上连接和井下连接三个部分组成的。因此,井下连接导线某一边方位角的总误差为: (9-24) 式中为投向误差,同样可按式(9-6)计算。但此时因两垂球线间的距离c加大,投向误差对定向精度的影响就不像一井定向那样起主要作用了。,煤矿测量规程规定,两井两次独立定向所算得的井下定向边的方位角之差,不应超过1。则一次定向的中误差为 若忽略投向误差,认为井上、下连接误差大致相同,则 下面分别研究井上、下连接误差m上和m下的估算方法。,一、 地面连接误差 两井定向时,井下连接导线某一边的方位角是按下式计算的(参阅图3-13及式

14、(3-14): (9-25) 式中 两垂球线的连线在地面坐标系统中的方位角; 两垂球线的连线在井下假定坐标系统中的方位角; 该边在假定坐标系统中的假定方位角。,式(9-25)中仅方位角 与地面连接有关,故地面连接误差 。 由第三章可知,两井定向的地面连接,根据两井距离的远近,可以采取两种不同的方案,现分述其连接误差如下。,(一) 由一个近井点向两垂球线敷设连接方案的误差 地面连接误差包括由近井点T到结点和由结点到两垂球线A、B所设两部分导线的误差。为了研究方便起见,假定一坐标系统:AB为y轴,垂直于AB的方向线为x轴。则 (9-26),Ry1A,Ry2A,Ry2B,A,B,1,2,终点,D,近

15、井点T,y,x,后视方向,式中 c两垂球线间的距离; mxA由结点到垂球线A间所测设的支导线误差所引起的A点在x轴方向上的位置误差 mxB由结点到垂球线B间所测设的支导线误差所引起的B点在x轴方向上的位置误差; n 由近井点到结点间的导线测角数; m由近井点到结点间导线的测角误差。,其中 上式中,式中 RyA由结点到垂球线A间的导线上各点到A的距离在AB线上的投影; RyB由结点到垂球线B间的导线上各点到B的距离在AB线上的投影; 导线各边与AB连线间的夹角。 在这种情况下,量边的系统误差对方位角 没有影响。,故量边误差对A、B点位的影响可用下式计算: 式中 a 量边的偶然误差影响系数; l

16、导线边长。,(二) 分别由两个近井点向相应的两垂球线连接方案 的误差 如图9-9所示,同样假定AB为y轴,垂直于AB的方向为x轴。则方位角AB的误差用下式计算: (9-27) 其中,近井点S,近井点T,后视方向,后视方向,A,B,x,y,式中 , 近井点S和T处的起始方位角中误差所引起的A、B垂球线在x轴上的误差; , 近井点S和T的x坐标误差,可按相对点位误差椭圆来求算。,二、井下连接误差 图9-10 为井下连接导线图,共测了n-1 个角和n条边。井下连接误差是由井下导线的测角误差m和量边误差ml所引起的,即 (9-28) 式中 , 测角和量边误差所引起的井下导线某边的方位角误差。,A,B,

17、x,x,x,1,2,n-1,l1,l2,ln,R,R,x,x,(一) 由井下导线测角误差所引起的连接误差 (9-29) 由式(9-25)对井下导线的角度取偏导数,得 因为方位角 是由地面连接测量算得的,与井下测量无关,故 。因此,上式可写为: (9-30) 由于井下导线各边的假定方位角 是由不同的角度算得的,因此对不同的边来说,其 之值也不同。,将 及 对的偏导数值代入上式,然后再代入式(9-29),即可求得不同边的方位角误差。经简化,即可得出由井下导线测角误差所引起的不同边的连接误差计算公式:,1,B,2,A,3,i-1,i+1,i,n-2,n-1,R1A,R2A,R3A,R(i-1)A,R

18、(n-1)B,R(n-2)B,R(i+1)B,RiB,上式中RA(见图9-11)为由导线点1、2、3、(i-1)到垂球线A的距离在AB连线上的投影;而RB则为由导线点i、i+1、(n-1)到垂球线B的距离在AB连接上的投影。,(二) 由井下导线量边误差所引起的连接误差 (9-32) 因 则 由于 及 均与井下量边无关,因此 求算偏导数,并将各偏导数代入式(9-32)中,得,考虑到量边中包括系统误差和偶然误差的影响,而量边的系统误差对方位角没有影响,因此,用钢尺量边时,上式可写成: (9-33) 上式即为计算井下导线量边误差而引起的任一边方位角的误差公式。式中i为井下导线各边与AB连线的夹角。,

19、(三) 由井下导线测角量边误差所引起的各边的连接总误差 (1) 第二边的井下连接误差为: (9-34) (2) 其他各边可类推。第i边则为 (9-35),三、 井上下两垂球线间距离的容许差值 在两井定向中,两垂球线之间的距离是由坐标反算得来的。据地面连接所算得的距离c同井下连接按假定坐标系统所算得的距离c加上改正数Hc/R后,在理论上应该相等。 但由于投点误差和井上下连接误差的影响,两者不可能相等,其差值为 但考虑到投点误差的影响很小,可忽略不计,故可把fc看做是井上、下连接误差所引起的。将连接导线看做始点为A、终点为B的支导线,根据第七章第五节的分析,并按煤矿测量规程要求,取二倍中误差作为容

20、许误差,,则得 (9-36) 式中 为导线测角中误差;Rxi为井下、地面(不包括近井点到结点)的连接导线各点到AB连线的垂直距离;mli、i分别为井下、地面(不包括近井点到结点)的连接导线各边的量边误差及各边与AB边线的夹角;H为井筒深度;R为地球平均曲率半经。 关于两井定向的平差,即差值fc的分配问题,通常用近似平差法解决。,第四节 陀螺经纬仪定向的误差分析及导线平差 一、 陀螺经纬仪定向的精度评定 陀螺经纬仪的定向精度主要以陀螺方位角一次测定中误差mT和一次定向中误差m表示。 (一) 陀螺方位角一次测定中误差 在待定边进行陀螺定向前,陀螺仪需在地面已知坐标方位角边上测定仪器常数。按煤矿测量

21、规程的规定,前后共需测46次,这样就可按白塞尔公式计算陀螺方位角一次测定中误差,即仪器常数一次测定中误差为:,m2=m2T =v2/(n-1) 式中,vi =i 平 n -测定仪器常数的次数 则测定仪器常数平均值的中误差为: m2平=m2T平 =m2T/ n (二) 一次定向中误差 井下陀螺定向边的坐标方位角为: 所以一次定向中误差为:,式中 - 仪器常数平均值中误差; - 待定边陀螺方位角平均值中误差; - 确定子午线收敛角的中误差。 因 0, 故 因井下定向边少, 观测次数少,需近似计算。一般井上下由同一个观测者用同一台仪器相同方法观测,则可认为 ,此时(3-2-3)程序:,当井下有多条陀

22、螺定向边时,或用同一台仪器在不同矿的井下进行过多条边的定向时,则可按双次观测列求算出井下 : 然后用下式求一次定向中误差;,二、 陀螺经纬仪一次测定方位角的中误差分析 如前所述,陀螺经纬仪的测量精度,以陀螺方位角一次测定中误差表示。不同的定向方法,其误差来源也有差异。这里仅以目前国内外最常用的跟踪逆转点法和中天法为例作一分析。其中所用的一些数据,是根据具体的仪器试验分析所得的,有一定的局限性,只能作为参考。但对掌握误差分析方法而言,这一点却是无关紧要的。,(一) 跟踪逆转点法定向时的误差分析 以JT15型陀螺经纬仪为例来进行探讨。按跟踪逆转点法进行陀螺定向时,主要误差来源有: 经纬仪测定方向的

23、误差; 上架式陀螺仪与经纬仪的连接误差; 悬挂带零位变动误差; 灵敏部摆动平衡位置的变动误差; 外界条件,如风流、气温及震动等因素的影响。 根据对JT15陀螺经纬仪的测试结果,对上述因素作如下分析。,1. 经纬仪测定方向的误差 一条测线一次观测的程序为:仪器在测站对中整平;测前以一测回测定测线方向值;以5个连续跟踪逆转点在度盘上的读数确定陀螺北方向值;测后以一测回测定测线方向值。这样,此项误差包括: (1) 对中误差 一般陀螺定向边都较长,当测线边长d=50m时,取eT=ec=0.8mm,则觇标对中误差 和仪器对中误差 为:,2. 上架式陀螺仪与经纬仪的连接误差 陀螺仪与经纬仪靠固定在照准部上

24、的过渡支架来连接。每次定向都要把陀螺仪安置在经纬仪支架上,这样由于每次拆装连接而造成的方向误差,根据用WILDT3经纬仪对三台仪器多次的实际测试,求得其连接中误差mE2,取mE=2。,3. 悬挂带零位变动误差 悬挂带对陀螺摆动系统的指向起阻碍作用,在实际观测时采用跟踪的方法可以消除悬挂带扭力的大部分影响。悬挂带材料的力学性质的优劣、陀螺运转造成的温升、外界气候的变化以及摆动系统的机械锁紧和释放等因素的影响,均会引起零位变位。根据对三台陀螺经纬仪的167次测试结果,求得悬挂带零位变动中误差ma=4。,4. 灵敏部摆动平衡位置的变动误差 影响摆动平衡位置变动的主要因素是:电源电压频率的变化引起角动

25、量的变化,灵敏部内部温度的变化引起重心位移以及由于温升造成悬挂带和导流丝的形变等因素,都会造成平衡位置的变动。由此而造成的误差多呈系统性,按JT15陀螺经纬仪灵敏部结构形式进行的98次试验,摆动平衡位置的最大离散度为1216,中误差mb=6,5. 外界条件,如风流、气温及震动等影响 这些条件的影响程度较为复杂,无法精确地一一测试,可取m外=5。 所以,测线陀螺方位角一次测定中误差为: 误差分析的结果说明,JT15陀螺经纬仪的设计精度是合理可行的。,(二) 中天法定向时的误差分析 以WILD GAK1-T2型陀螺经纬仪为例进行分析探讨。用中天法定向时的主要误差来源有: 经纬仪测定测线方向的误差;

26、 陀螺仪与经纬仪的连接误差: 悬挂带零位变动误差和摆动平衡位置的变动误差; 中天时间的测定误差和摆幅的读数误差; 外界条件的影响。 上述误差中与逆转点法定向时相同的部分,这里就不再重复分析。,1. 经纬仪测定方向的误差 (1) 对中误差 当定向边边长d=50m时,仪器及觇标的对中误差由上面分析得: (2) 测线前后两测回的平均值中误差 则,2. 陀螺仪与经纬仪的连接误差 取mE=2 3.悬挂带零位变动误差和摆动平衡位置的变动误差 悬挂带零位变动误差将综合考虑于零位改正中,摆动平衡位置的变动误差mb=6。,4. 中天时间的测定误差和摆幅的读数误差 由式(3-47) 知 如考虑到零位改正(见式(3

27、-42),则 (9-44) 因式(9-44) 中N和的测量误差影响较小,一般可忽略不计。考虑到N=cat,则 (9-45),故 根据室内在4 h内测定彼此不相关的不跟踪摆动周期共45个所求得的平均值T2平=38834 s,用白塞尔公式求算得mT2=0.137s; 由图3-28可知T2=ti+1-ti 故 则,(7) NT的测定误差mNT随N 的增加而增大,根据实际资料分析:近似陀螺北偏离陀北的H值小于5为宜。取H=5。 将上述有关数值代入式(9-45),得 则 =8.86 所以,测线陀螺方位角一次测定中误差为:,三、 陀螺仪定向导线的平差 由于目前陀螺经纬仪的定向精度在1560之间,所以陀螺定

28、向边不能完全作为坚强边来控制7和15基本导线(15陀螺可控制15导线),因而陀螺定向边应和导线边一起作联合平差。下面介绍两种类型的陀螺仪定向导线的平差方法。,1. 具有两条陀螺定向边导线的平差 图9-12 具有两条陀螺定向边导线的平差示意图 图9-12中的AB及CD边为陀螺定向边,其坐标方位角分别为1与2 ,平差步骤如下: (1) 求算陀螺定向边AB与CD的定向中误差m1与m2及导线测角中误差mm1与m2可按式(9-40)、(9-41)求算,m或按导线的实际情况来求,或按闭合导线的闭合要求,或按双次观测列求得。,(2) 按条件观测平差,列出角改正数条件方程式 如图9-12所示,导线的角闭合差为

29、: 改正数条件方程式为: 式中 分别为陀螺定向边坐标方位角 的改正数; 导线中角度 的改正数; n 导线中角度个数。,(3) 确定定向边方位角和角度的权 当等精度观测时,取导线的测角中误差m为单位权中误差。 即p=1 因为p=2/m2d,则定向边方位角的权为 (9-46) 权倒数为: (9-47),(4) 组成法方程式 NK+W=0 其中 (9-48) 解法方程式得 (9-49),(5) 计算各改正数 导线各角度的改正数为: (9-50) 定向边AB的方位角 的改正数为: (9-51) 定向边CB的方位角 的改正数为: (9-52) 将各观测值加入所求得的相应的改正数v,就可得到各方位角和导线

30、角的最或是值。,2.具有三条陀螺定向边导线的平差 如图9-13所示,AB、CB、EF为陀螺定向边,其相应的坐标方位角为1、2、3,这时可将整个导线分为两部分,即导线和导线。平差步骤如下: 图9-13 具有三条陀螺定向边导线的平差图,(1) 求定向边定向中误差及测角中误差 陀螺定向边AB、CD、EF的定向中误差 及导线测角中误差 (等精度观测时)的计算方法同前。,(2) 按条件观测平差列出角改正数条件方程式 导线、的角闭合差为 改正数条件方程式为: 式中 分别为陀螺定向边坐标方位角 的改正数; 导线中角度 的改正数; 导线中角度 的改正数; n1、n2 分别为导线、中角度的个数。,(3) 确定定向边方位角和角度的权 当导线等精度观测时,取导线的测角中误差 为单位权中误差,即 ,则定向边坐标方位角的权为: (9-53) 权倒数为: (9-54),(4) 组成法方程式 式中 (9-55) 解法方程式,求得 (9-56),(5) 计算各改正数 导线各角度的改正数为: (9-57) 导线各角度的改正数为: (9-58) 陀螺定向边方位角的改正数为: (9-59),(6) 计算各观测值的最或是值 设 为定向边方位角的最或是值, 分别为导线、各角度的最或是值 为导线、角度最或是值之和,则有 具有多个陀螺定向边导线的平差仍可按上述原则进行,只不过多列几个法方程式。,

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 其他


经营许可证编号:宁ICP备18001539号-1